#215: Driscoll&Healy integration for Multipole ------------------------------------+--------------------------------------- Reporter: eschnett | Owner: eschnett Type: enhancement | Status: closed Priority: major | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: fixed | Keywords: Multipole ------------------------------------+--------------------------------------- Changes (by eschnett):
* status: assigned => closed * resolution: => fixed
Comment:
I found the missing factor pi. This factor is actually present in the (unnumbered!) equations in the Driscoll&Healy paper. In my tests, I left it out because I implicitly assumed that theta ranges [0...1], whereas of course it ranger [0...pi]. This made the results too large by a factor of pi.
There was also an error in the index bounds. Driscoll&Healy want a point on the equator, which corresponds to even ntheta, not odd ntheta as previously in the code.
Applied.