new testsuite test

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

1 files changed, 129 insertions, 0 deletions

diff --git a/test/misner1.2-025.par b/test/misner1.2-025.par
new file mode 100644
index 0000000..3708885
--- /dev/null
+++ b/test/misner1.2-025.par
@@ -0,0 +1,129 @@
+# 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.
+
+########################################
+ 
+#
+# ***** grid *****
+#
+
+ActiveThorns = "PUGH"
+
+Driver::global_nx   = 43
+Driver::global_ny   = 43
+Driver::global_nz   = 83
+driver::ghost_size  = 2
+
+ActiveThorns = "CoordBase CartGrid3D"
+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"
+SpaceMask::use_mask                = "true"
+
+ADMBase::evolution_method       = "none"
+ADMBase::initial_lapse          = "Cadez"
+ADMBase::lapse_evolution_method = "static"
+ADMBase::initial_shift          = "zero"
+ADMBase::shift_evolution_method = "static"
+
+########################################
+
+#
+# automagically create output directory
+#
+ActiveThorns = "IOUtil PUGHSlab IOASCII"
+IO::out_dir          = "misner1.2-025"
+IO::parfile_write    = "no"
+IOASCII::out2d_every = 1
+IOASCII::out2d_vars  = "SpaceMask::emask SphericalSurface::sf_radius[1] SphericalSUrface::sf_radius[2]"
+
+########################################
+
+#
+# SphericalSurface
+#
+ActiveThorns = "SphericalSurface"
+SphericalSurface::nsurfaces = 3
+SphericalSurface::maxntheta = 62
+SphericalSurface::maxnphi   = 62
+
+SphericalSurface::nghoststheta[1] = 2
+SphericalSurface::ntheta[1] = 32
+SphericalSurface::nphi  [1] = 32
+
+SphericalSurface::nghostsphi  [2] = 2
+SphericalSurface::ntheta[1] = 62
+SphericalSurface::nphi  [1] = 62
+
+########################################
+
+#
+# AHFinderDirect
+#
+
+ActiveThorns = "PUGHInterp PUGHReduce AEILocalInterp AHFinderDirect"
+
+AHFinderDirect::Theta_norm_for_convergence = 1.0e-12
+
+AHFinderDirect::N_horizons = 2
+
+AHFinderDirect::set_mask_for_all_horizons = "true"
+# only test old-style mask for now,
+# since we don't have a thorn to set up the new-style bit fields yet
+AHFinderDirect::set_old_style_mask = "true"
+AHFinderDirect::set_new_style_mask = "false"
+# have mask directly show horizon position for testing (no buffer zone)
+AHFinderDirect::mask_radius_multiplier[1] = 1.0
+AHFinderDirect::mask_radius_multiplier[2] = 1.0
+AHFInderDirect::mask_radius_offset[1]     = 0.0
+AHFInderDirect::mask_radius_offset[2]     = 0.0
+
+AHFinderDirect::which_surface_to_store_info[1] = 1
+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::which_surface_to_store_info[2] = 2
+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