tests with geometry hardwired to Schw/EF

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

4 files changed, 211 insertions, 0 deletions

diff --git a/run/test-ahfinderdirect/Schw-hardwire/large-ahf-posn.par b/run/test-ahfinderdirect/Schw-hardwire/large-ahf-posn.par
new file mode 100644
index 0000000..97c7cb2
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw-hardwire/large-ahf-posn.par
@@ -0,0 +1,57 @@
+# 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 = "Schwarzschild/EF"
+Exact::Schwarzschild_EF__mass = 1.0
+
+AHFinderDirect::method = "Newton solve"
+AHFinderDirect::h_file_name      = "large-ahf-posn.h.dat"
+AHFinderDirect::H_of_h_file_name = "large-ahf-posn.H.dat"
+
+AHFinderDirect::origin_x = 0.0
+AHFinderDirect::origin_y = 0.0
+AHFinderDirect::origin_z = 0.0
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 5.0
+
+AHFinderDirect::hardwire_Schwarzschild_EF__x_posn = 0.5
+AHFinderDirect::hardwire_Schwarzschild_EF__y_posn = 0.7
+AHFinderDirect::hardwire_Schwarzschild_EF__z_posn = 0.9
+
+AHFinderDirect::hardwire_Schwarzschild_EF = "true"
+
+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_center = 0.5
+AHFinderDirect::initial_guess__ellipsoid__y_center = 0.7
+AHFinderDirect::initial_guess__ellipsoid__z_center = 0.9
+AHFinderDirect::initial_guess__ellipsoid__x_radius = 1.8
+AHFinderDirect::initial_guess__ellipsoid__y_radius = 1.8
+AHFinderDirect::initial_guess__ellipsoid__z_radius = 1.8
diff --git a/run/test-ahfinderdirect/Schw-hardwire/large-ahf.par b/run/test-ahfinderdirect/Schw-hardwire/large-ahf.par
new file mode 100644
index 0000000..f1f85bf
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw-hardwire/large-ahf.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 = "Schwarzschild/EF"
+Exact::Schwarzschild_EF__mass = 1.0
+
+AHFinderDirect::method = "Newton solve"
+AHFinderDirect::h_file_name      = "large-ahf.h.dat"
+AHFinderDirect::H_of_h_file_name = "large-ahf.H.dat"
+
+AHFinderDirect::origin_x = 0.0
+AHFinderDirect::origin_y = 0.0
+AHFinderDirect::origin_z = 0.0
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 5.0
+
+AHFinderDirect::hardwire_Schwarzschild_EF = "true"
+
+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 = 1.8
+AHFinderDirect::initial_guess__ellipsoid__y_radius = 1.8
+AHFinderDirect::initial_guess__ellipsoid__z_radius = 1.8
diff --git a/run/test-ahfinderdirect/Schw-hardwire/large-posn.par b/run/test-ahfinderdirect/Schw-hardwire/large-posn.par
new file mode 100644
index 0000000..fae757a
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw-hardwire/large-posn.par
@@ -0,0 +1,54 @@
+# 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 = "Schwarzschild/EF"
+Exact::Schwarzschild_EF__mass = 1.0
+
+AHFinderDirect::method = "horizon function"
+AHFinderDirect::h_file_name      = "large-posn.h.dat"
+AHFinderDirect::H_of_h_file_name = "large-posn.H.dat"
+
+AHFinderDirect::origin_x = 0.0
+AHFinderDirect::origin_y = 0.0
+AHFinderDirect::origin_z = 0.0
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 5.0
+
+AHFinderDirect::hardwire_Schwarzschild_EF = "true"
+AHFinderDirect::hardwire_Schwarzschild_EF__x_posn = 0.5
+AHFinderDirect::hardwire_Schwarzschild_EF__y_posn = 0.7
+AHFinderDirect::hardwire_Schwarzschild_EF__z_posn = 0.9
+
+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 = "Kerr/Kerr"
+AHFinderDirect::initial_guess__Kerr__spin = 0.0
+AHFinderDirect::initial_guess__Kerr__x_posn = 0.5
+AHFinderDirect::initial_guess__Kerr__y_posn = 0.7
+AHFinderDirect::initial_guess__Kerr__z_posn = 0.9
diff --git a/run/test-ahfinderdirect/Schw-hardwire/large.par b/run/test-ahfinderdirect/Schw-hardwire/large.par
new file mode 100644
index 0000000..4b8c514
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw-hardwire/large.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 = "Schwarzschild/EF"
+Exact::Schwarzschild_EF__mass = 1.0
+
+AHFinderDirect::method = "horizon function"
+AHFinderDirect::h_file_name      = "large.h.dat"
+AHFinderDirect::H_of_h_file_name = "large.H.dat"
+
+AHFinderDirect::origin_x = 0.0
+AHFinderDirect::origin_y = 0.0
+AHFinderDirect::origin_z = 0.0
+AHFinderDirect::patch_system_type = "full sphere"
+AHFinderDirect::N_ghost_points = 2
+AHFinderDirect::N_overlap_points = 1
+AHFinderDirect::delta_drho_dsigma = 5.0
+
+AHFinderDirect::hardwire_Schwarzschild_EF = "true"
+
+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