Implement Driscoll&Healy integration

Implements a more accurate integration over the sphere using an
algorithm by Driscoll & Healy. This algorithm uses Gaussian
integration weights, leading (almost) to exponential convergence.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Multipole/trunk@83 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843

1 files changed, 3 insertions, 2 deletions

diff --git a/param.ccl b/param.ccl
index 9af3650..b8ed861 100644
--- a/param.ccl
+++ b/param.ccl
@@ -29,8 +29,9 @@ CCTK_STRING coord_system "What is the coord system?"
 
 KEYWORD integration_method "How to do surface integrals" STEERABLE=always
 {
-  "midpoint" :: "Midpoint rule (1st order)"
-  "Simpson"  :: "Simpson's rule (4th order)"
+  "midpoint"      :: "Midpoint rule (1st order)"
+  "Simpson"       :: "Simpson's rule (4th order) [requires even ntheta and nphi]"
+  "DriscollHealy" :: "Driscoll & Healy (exponentially convergent) [requires even ntheta]"
 } "midpoint"
 
 CCTK_INT out_every "How often to output" \