more tests...

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

6 files changed, 260 insertions, 2 deletions

diff --git a/run/test-ahfinderdirect/Schw/large-ahf-np.par b/run/test-ahfinderdirect/Schw/large-ahf-dx8.par
index 79db2a6..aa02fd8 100644
--- a/run/test-ahfinderdirect/Schw/large-ahf-np.par
+++ b/run/test-ahfinderdirect/Schw/large-ahf-dx8.par
@@ -26,8 +26,9 @@ Exact::exact_model = "Schwarzschild/EF"
 Exact::Schwarzschild_EF__mass = 1.0
 
 AHFinderDirect::method = "Newton solve"
-AHFinderDirect::h_file_name      = "large-ahf-np.h.dat"
-AHFinderDirect::H_of_h_file_name = "large-ahf-np.H.dat"
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "large-ahf-dx8.h.dat"
+AHFinderDirect::H_of_h_file_name = "large-ahf-dx8.H.dat"
 
 AHFinderDirect::origin_x = 0.0
 AHFinderDirect::origin_y = 0.0
@@ -42,6 +43,8 @@ AHFinderDirect::geometry_interpolator_pars = "order=3"
 AHFinderDirect::interpatch_interpolator_name = "generalized polynomial interpolation"
 AHFinderDirect::interpatch_interpolator_pars = "order=3"
 
+AHFinderDirect::Jacobian_perturbation_amplitude = 1.0e-8
+
 AHFinderDirect::initial_guess_method = "ellipsoid"
 AHFinderDirect::initial_guess__ellipsoid__x_radius = 1.8
 AHFinderDirect::initial_guess__ellipsoid__y_radius = 1.8
diff --git a/run/test-ahfinderdirect/Schw/offset-ahf-dx8.par b/run/test-ahfinderdirect/Schw/offset-ahf-dx8.par
new file mode 100644
index 0000000..a63ebca
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw/offset-ahf-dx8.par
@@ -0,0 +1,51 @@
+# 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::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "offset-ahf-dx8.h.dat"
+AHFinderDirect::H_of_h_file_name = "offset-ahf-dx8.H.dat"
+
+AHFinderDirect::Jacobian_perturbation_amplitude = 1.0e-8
+
+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 = 5.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 = 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/offset-ahf-it5.par b/run/test-ahfinderdirect/Schw/offset-ahf-it5.par
new file mode 100644
index 0000000..e8515c7
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw/offset-ahf-it5.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::max_Newton_iterations = 5
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "offset-ahf-it5.h.dat"
+AHFinderDirect::H_of_h_file_name = "offset-ahf-it5.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 = 5.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 = 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/offset-ahf-np6.par b/run/test-ahfinderdirect/Schw/offset-ahf-np6.par
new file mode 100644
index 0000000..d4879c8
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw/offset-ahf-np6.par
@@ -0,0 +1,52 @@
+# 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::output_h_and_H_at_each_Newton_iteration = "true"
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "offset-ahf-np6.h.dat"
+AHFinderDirect::H_of_h_file_name = "offset-ahf-np6.H.dat"
+AHFinderDirect::Jacobian_method = "numerical perturbation"
+AHFinderDirect::Jacobian_perturbation_amplitude = 1.0e-6
+
+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 = 5.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 = 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/offset-ahf-np8.par b/run/test-ahfinderdirect/Schw/offset-ahf-np8.par
new file mode 100644
index 0000000..c197027
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw/offset-ahf-np8.par
@@ -0,0 +1,52 @@
+# 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::output_h_and_H_at_each_Newton_iteration = "true"
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "offset-ahf-np8.h.dat"
+AHFinderDirect::H_of_h_file_name = "offset-ahf-np8.H.dat"
+AHFinderDirect::Jacobian_method = "numerical perturbation"
+AHFinderDirect::Jacobian_perturbation_amplitude = 1.0e-8
+
+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 = 5.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 = 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/offset-ahf.par b/run/test-ahfinderdirect/Schw/offset-ahf.par
new file mode 100644
index 0000000..0dea631
--- /dev/null
+++ b/run/test-ahfinderdirect/Schw/offset-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::output_h_and_H_at_each_Newton_iteration = "true"
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_file_name      = "offset-ahf.h.dat"
+AHFinderDirect::H_of_h_file_name = "offset-ahf.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 = 5.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 = 1.8
+AHFinderDirect::initial_guess__ellipsoid__y_radius = 1.8
+AHFinderDirect::initial_guess__ellipsoid__z_radius = 1.8