#2398: Multipole: Default (midpoint) integration method incorrectly implemented, yielding wrong results... especially with m=odd modes
Reporter:Zach Etienne
Status:open
Milestone:
Version:
Type:bug
Priority:major
Component:EinsteinToolkit thorn

Comment (by Zach Etienne):

Here is my proposed patch for the upcoming release branch, which replaces the default integration with trapezoidal. I checked the test/ directory, and all parfiles select the integration_method explicitly, so no patch should be needed for those.

diff --git a/Multipole/param.ccl b/Multipole/param.ccl
index b2ac27e5..8125e48c 100644
--- a/Multipole/param.ccl
+++ b/Multipole/param.ccl
@@ -25,11 +25,11 @@ CCTK_STRING coord_system "What is the coord system?"

 KEYWORD integration_method "How to do surface integrals" STEERABLE=always
 {
-  "midpoint

--

Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2398/multipole-default-midpoint-integration