Comment (by Roland Haas):

I updated the test proposed (since for eg trapezoidal rule, the 0,0 mode was not integration to 1+/01e-6 but only to 0.996 or so), and added a test for the midpoint rule as well.

Updated pull request is here:

https://bitbucket.org/einsteintoolkit/einsteinanalysis/pull-requests/9/multipole-turn-orthogonality-test-into/diff

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2210/multipole-added-code-to-test-mode