some more trial par files

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

2 files changed, 154 insertions, 0 deletions

diff --git a/par/Kerr-umfpack.par b/par/Kerr-umfpack.par
new file mode 100644
index 0000000..4a2ceac
--- /dev/null
+++ b/par/Kerr-umfpack.par
@@ -0,0 +1,59 @@
+# 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
+
+ActiveThorns = "PUGH"
+driver::ghost_size = 2
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 19
+
+ActiveThorns = "CartGrid3D"
+grid::domain = "bitant"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
+
+ActiveThorns = "ADMBase ADMCoupling StaticConformal CoordGauge SpaceMask Exact"
+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_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.6
+
+########################################
+
+ActiveThorns = "IOUtil"
+IOUtil::parfile_write = "no"
+
+########################################
+
+ActiveThorns = "always AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+
+AHFinderDirect::verbose_level = "algorithm details"
+AHFinderDirect::print_timing_stats = "true"
+AHFinderDirect::Jacobian_store_solve_method = "row-oriented sparse matrix/UMFPACK"
+
+AHFinderDirect::output_Theta_every = 1
+AHFinderDirect::h_base_file_name     = "Kerr.h"
+AHFinderDirect::Theta_base_file_name = "Kerr.Theta"
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::origin_x[1] = 0.5
+AHFinderDirect::origin_y[1] = 0.7
+AHFinderDirect::origin_z[1] = 0.0
+
+AHFinderDirect::initial_guess_method[1] = "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
diff --git a/par/misner-umfpack.par b/par/misner-umfpack.par
new file mode 100644
index 0000000..11c8ef1
--- /dev/null
+++ b/par/misner-umfpack.par
@@ -0,0 +1,95 @@
+# This parameter file sets up Misner initial data and finds the
+# apparent horizons in it with AHFinderDirect.
+#
+# Note that the grid is very high-resolution, but small (the outer
+# boundaries are very close in).  This isn't a problem for the purposes
+# of this test.  This run takes about 185 megabytes of memory, and
+# 36 cpu seconds on a 1.07GHz PIII laptop.
+
+########################################
+ 
+#
+# ***** grid *****
+#
+
+ActiveThorns = "CartGrid3D PUGH"
+
+Driver::global_nx   = 83
+Driver::global_ny   = 83
+Driver::global_nz   = 115
+PUGH::ghost_size    = 2
+
+Grid::domain        = "octant"
+Grid::type          = "byspacing"
+Grid::dxyz          = 0.025
+Grid::avoid_origin  = "true"
+
+########################################
+
+#
+# physics
+#
+
+ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask"
+ActiveThorns = "IDAnalyticBH"
+
+ADMBase::initial_data              = "misner_bh"
+IDAnalyticBH::mu                   = 1.2
+
+########################################
+
+#
+# numerical methods
+#
+
+Cactus::cctk_itlast = 0
+ADMBase::metric_type            = "static conformal"
+StaticConformal::conformal_storage = "factor+derivs+2nd derivs"
+ADMBase::evolution_method       = "none"
+ADMBase::initial_lapse          = "Cadez"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::initial_shift          = "zero"
+ADMBase::shift_evolution_method = "static"
+
+########################################
+
+ActiveThorns = "IOUtil"
+IOUtil::parfile_write = "no"
+
+########################################
+
+#
+# AHFinderDirect
+#
+
+ActiveThorns = "always AEILocalInterp PUGHInterp PUGHReduce AHFinderDirect"
+
+AHFinderDirect::print_timing_stats = "true"
+##AHFinderDirect::verbose_level = "algorithm details"
+
+AHFinderDirect::h_base_file_name              = "misner.h"
+AHFinderDirect::BH_diagnostics_base_file_name = "misner.BH_diagnostics"
+
+AHFinderDirect::Jacobian_store_solve_method = "row-oriented sparse matrix/UMFPACK"
+AHFinderDirect::N_horizons = 2
+
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+AHFinderDirect::origin_z[1] = 1.0
+AHFinderDirect::initial_guess_method[1] = "coordinate sphere"
+AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0
+AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 1.1
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.60
+
+AHFinderDirect::N_zones_per_right_angle[2] = 30
+AHFinderDirect::origin_x[2] = 0.0
+AHFinderDirect::origin_y[2] = 0.0
+AHFinderDirect::origin_z[2] = 0.0
+AHFinderDirect::initial_guess_method[2] = "coordinate ellipsoid"
+AHFinderDirect::initial_guess__coord_ellipsoid__x_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_ellipsoid__y_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_ellipsoid__z_center[2] = 0.0
+AHFinderDirect::initial_guess__coord_ellipsoid__x_radius[2] = 0.6
+AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.6
+AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.8