Peter Diener's suggestions of checking radius of convergence for initial guess

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

2 files changed, 99 insertions, 0 deletions

diff --git a/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess-2.5.par b/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess-2.5.par
new file mode 100644
index 0000000..bb0fd58
--- /dev/null
+++ b/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess-2.5.par
@@ -0,0 +1,50 @@
+# parameter file for patch system test
+ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling CoordGauge Exact AHFinderDirect"
+
+# PUGH
+driver::ghost_size = 0
+driver::global_nx = 61
+driver::global_ny = 61
+driver::global_nz = 61
+
+# CartGrid3D
+grid::domain = "full"
+grid::type = "byrange"
+grid::xyzmin = -3.0
+grid::xyzmax =  3.0
+
+# ADMBase
+ADMBase::initial_lapse = "exact"
+ADMBase::initial_shift = "exact"
+ADMBase::initial_data = "exact"
+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::method = "Newton solve"
+AHFinderDirect::h_file_name      = "small-ahf-badguess-2.5.h.dat"
+AHFinderDirect::H_of_h_file_name = "small-ahf-badguess-2.5.H.dat"
+
+AHFinderDirect::origin_x = 0.5
+AHFinderDirect::origin_y = 0.7
+AHFinderDirect::origin_z = 0.9
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 9.0
+
+AHFinderDirect::geometry_interpolator_name = "generalized polynomial interpolation"
+AHFinderDirect::geometry_interpolator_pars = "order=3"
+AHFinderDirect::interpatch_interpolator_name = "generalized polynomial interpolation"
+AHFinderDirect::interpatch_interpolator_pars = "order=3"
+
+AHFinderDirect::initial_guess_method = "ellipsoid"
+AHFinderDirect::initial_guess__ellipsoid__x_radius = 2.5
+AHFinderDirect::initial_guess__ellipsoid__y_radius = 2.5
+AHFinderDirect::initial_guess__ellipsoid__z_radius = 2.5
+AHFinderDirect::initial_guess__ellipsoid__x_center = -0.25
diff --git a/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess.par b/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess.par
new file mode 100644
index 0000000..f52c6d2
--- /dev/null
+++ b/run/test-ahfinderdirect/Kerr-offset/small-ahf-badguess.par
@@ -0,0 +1,49 @@
+# parameter file for patch system test
+ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling CoordGauge Exact AHFinderDirect"
+
+# PUGH
+driver::ghost_size = 0
+driver::global_nx = 61
+driver::global_ny = 61
+driver::global_nz = 61
+
+# CartGrid3D
+grid::domain = "full"
+grid::type = "byrange"
+grid::xyzmin = -3.0
+grid::xyzmax =  3.0
+
+# ADMBase
+ADMBase::initial_lapse = "exact"
+ADMBase::initial_shift = "exact"
+ADMBase::initial_data = "exact"
+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::method = "Newton solve"
+AHFinderDirect::h_file_name      = "small-ahf-badguess.h.dat"
+AHFinderDirect::H_of_h_file_name = "small-ahf-badguess.H.dat"
+
+AHFinderDirect::origin_x = 0.5
+AHFinderDirect::origin_y = 0.7
+AHFinderDirect::origin_z = 0.9
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 9.0
+
+AHFinderDirect::geometry_interpolator_name = "generalized polynomial interpolation"
+AHFinderDirect::geometry_interpolator_pars = "order=3"
+AHFinderDirect::interpatch_interpolator_name = "generalized polynomial interpolation"
+AHFinderDirect::interpatch_interpolator_pars = "order=3"
+
+AHFinderDirect::initial_guess_method = "ellipsoid"
+AHFinderDirect::initial_guess__ellipsoid__x_radius = 2.0
+AHFinderDirect::initial_guess__ellipsoid__y_radius = 2.0
+AHFinderDirect::initial_guess__ellipsoid__z_radius = 2.0