change the default geometry interpolator from order 3 down to order 2

this will make AHFinderDirect faster (typically about 2.5 times faster),
but less accurate (on analytic data at dx=0.1m or so, typical accuracies
of 1e-4m or so, instead of 1e-5m or so for order=3).  In practice, for
most of our runs the AHFinderDirect (in)accuracy is dominated by the
inaccuracy of the numerically computed g_ij and K_ij fed into it, so
the slightly lower AH-finding accuracy is unimportant.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1178 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 4 insertions, 7 deletions

diff --git a/param.ccl b/param.ccl
index 6e6c36c..a64c354 100644
--- a/param.ccl
+++ b/param.ccl
@@ -1018,18 +1018,15 @@ string geometry_interpolator_name \
 } "Hermite polynomial interpolation"
 
 #
-# The defaults here are set up for 3rd order Hermite interpolation,
-# which requires driver::ghost_size >= 3.  If you want to run with
-# driver::ghost_size == 2, you need to lower the interpolation order
-# to 2; see  par/misner-ghost2-mask-noshrink.par  for an example of
-# how to do this.
+# The defaults here are set up for 2nd order Hermite interpolation,
+# which requires driver::ghost_size >= 2.
 #
 string geometry_interpolator_pars \
   "parameters for the geometry interpolator"
 {
 ".*" :: \
   "any string acceptable to Util_TableSetFromString() and to the interpolator"
-} "order=3 \
+} "order=2 \
    boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10 1.0e-10} \
    boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0 0.0 0.0}"
 
@@ -1166,7 +1163,7 @@ string surface_interpolator_pars \
 {
 ".*" :: \
   "any string acceptable to Util_TableSetFromString() and to the interpolator"
-} "order=3 \
+} "order=2 \
    boundary_off_centering_tolerance={1.0e-10 1.0e-10 1.0e-10 1.0e-10} \
    boundary_extrapolation_tolerance={0.0 0.0 0.0 0.0}"