add beginnings of support for using new CCTK_InterpGridArrays()

global interpolator & multiprocessor operation

... at present the CCTK_InterpGridArrays() support works fine,
    but only on 1 processor, I seem to get an infinite loop on
    multiple processors


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

1 files changed, 36 insertions, 20 deletions

diff --git a/run/test-ahfinderdirect/misc/Kerr.par b/run/test-ahfinderdirect/misc/Kerr.par
index 9ba6be0..003c5ad 100644
--- a/run/test-ahfinderdirect/misc/Kerr.par
+++ b/run/test-ahfinderdirect/misc/Kerr.par
@@ -1,22 +1,25 @@
-# parameter file for patch system test
-ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling StaticConformal CoordGauge Exact AHFinderDirect"
+# This parameter file sets up Kerr/Kerr-Schild initial data, then
+# finds the apparent horizon in it.  The local coordinate system origin
+# and the initial guess are both deliberately de-centered with respect
+# to the black hole, to make this a non-trivial test for the apparent
+# horizon finder.
 
 # flesh
 cactus::cctk_itlast = 0
 
-# PUGH
-driver::ghost_size = 0
-driver::global_nx = 61
-driver::global_ny = 61
-driver::global_nz = 61
+ActiveThorns = "PUGH"
+driver::ghost_size = 2
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 17
 
-# CartGrid3D
-grid::domain = "full"
-grid::type = "byrange"
-grid::xyzmin = -3.0
-grid::xyzmax =  3.0
+ActiveThorns = "CartGrid3D"
+grid::domain = "bitant"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
 
-# ADMBase
+ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge Exact"
 ADMBase::initial_lapse = "exact"
 ADMBase::initial_shift = "exact"
 ADMBase::initial_data = "exact"
@@ -24,17 +27,30 @@ ADMBase::lapse_evolution_method = "static"
 ADMBase::shift_evolution_method = "static"
 ADMBase::metric_type = "physical"
 
-# Exact
 Exact::exact_model = "Kerr/Kerr-Schild"
 Exact::Kerr_KerrSchild__mass = 1.0
 Exact::Kerr_KerrSchild__spin = 0.6
 
-AHFinderDirect::find_AHs_at_poststep = "false"
+########################################
 
-AHFinderDirect::patch_system_type = "full sphere"
-AHFinderDirect::N_zones_per_right_angle = 18
+ActiveThorns = "LocalInterp PUGHInterp AHFinderDirect"
+
+AHFinderDirect::find_AHs_at_postinitial = "true"
+AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::final_Theta_update_if_Delta_h_converged = "true"
+AHFinderDirect::how_often_to_output_Theta = 1
+AHFinderDirect::h_base_file_name     = "Kerr.h"
+AHFinderDirect::Theta_base_file_name = "Kerr.Theta"
 
 AHFinderDirect::N_horizons = 1
-AHFinderDirect::initial_guess_method = "Kerr/Kerr-Schild"
-AHFinderDirect::initial_guess__Kerr_KerrSchild__mass[1] = 1.0
-AHFinderDirect::initial_guess__Kerr_KerrSchild__spin[1] = 0.6
+AHFinderDirect::origin_x[1] = 0.5
+AHFinderDirect::origin_y[1] = 0.7
+AHFinderDirect::origin_z[1] = 0.0
+
+AHFinderDirect::initial_guess_method = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.2
+AHFinderDirect::initial_guess__coord_sphere__y_center[1] =  0.3
+AHFinderDirect::initial_guess__coord_sphere__z_center[1] =  0.0
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0