demo par files

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

5 files changed, 420 insertions, 0 deletions

diff --git a/par/Kerr-tiny.par b/par/Kerr-tiny.par
new file mode 100644
index 0000000..8172bed
--- /dev/null
+++ b/par/Kerr-tiny.par
@@ -0,0 +1,54 @@
+# 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.
+#
+# This parameter is "tiny" in the sense that it sets only a small number
+# of AHFinderDirect parameters.
+
+ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling StaticConformal CoordGauge Exact AHFinderDirect"
+
+# flesh
+cactus::cctk_itlast = 0
+
+# PUGH
+driver::ghost_size = 2
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 17
+
+# CartGrid3D
+grid::domain = "bitant"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
+
+# 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::find_AHs_at_poststep = "false"	# no time evolution
+AHFinderDirect::h_base_file_name     = "Kerr-tiny.h"
+
+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 = "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/Kerr.par b/par/Kerr.par
new file mode 100644
index 0000000..f049d8a
--- /dev/null
+++ b/par/Kerr.par
@@ -0,0 +1,57 @@
+# 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.
+ActiveThorns = "CartGrid3D LocalInterp PUGH ADMBase ADMCoupling StaticConformal CoordGauge Exact AHFinderDirect"
+
+# flesh
+cactus::cctk_itlast = 0
+
+# PUGH
+driver::ghost_size = 2
+driver::global_nx = 31
+driver::global_ny = 31
+driver::global_nz = 17
+
+# CartGrid3D
+grid::domain = "bitant"
+grid::avoid_origin = "false"
+grid::type = "byspacing"
+grid::dxyz = 0.2
+
+# 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::find_AHs_at_postinitial = "true"
+AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::verbose_level = "algorithm highlights"
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::how_often_to_output_H_of_h = 1
+AHFinderDirect::h_base_file_name      = "Kerr.h"
+AHFinderDirect::H_of_h_base_file_name = "Kerr.H"
+
+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 = "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/README b/par/README
new file mode 100644
index 0000000..0f3f353
--- /dev/null
+++ b/par/README
@@ -0,0 +1,16 @@
+This directory contains sample parameter files.
+
+Kerr.par is a good one to start with if you want a simple example
+of AHFinderDirect usage.  It's my standard "sanity-check" test case
+for AHFinderDirect, and runs in a < 10 seconds on a laptop.  It sets
+up Kerr/Kerr-Schild initial data and finds the apparent horizon in
+it, with the coordinate origin and initial guess offset to make this
+a nontrivial test of the apparent horizon finder.
+
+misner-init.par is another fairly simple parameter file, which sets up
+Misner initial data and then finds the apparent horizons in it.
+
+misner-run.par is a more complicated parameter file, which sets up the
+same Misner initial data, but now time-evolves it, finding the apparent
+horizons at each time step.  (It also runs AHFinder every 20 time steps
+for comparison.)  It takes a lot of cpu and memory to run.
diff --git a/par/misner-init.par b/par/misner-init.par
new file mode 100644
index 0000000..a4bb26b
--- /dev/null
+++ b/par/misner-init.par
@@ -0,0 +1,93 @@
+# 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 213 megabytes of memory on my laptop.
+
+ActiveThorns = "CartGrid3D LocalInterp IOBasic IOUtil IOASCII   PUGH PUGHSlab   ADMBase ADMCoupling CoordGauge StaticConformal IDAnalyticBH   AHFinderDirect"
+
+########################################
+
+#
+# physics
+#
+
+ADMBase::initial_data              = "misner_bh"
+IDAnalyticBH::mu                   = 2.0
+
+########################################
+ 
+#
+# ***** grid *****
+#
+
+Driver::global_nx   = 105
+Driver::global_ny   = 105
+Driver::global_nz   = 80
+PUGH::ghost_size    = 2
+
+Grid::domain        = "bitant"
+Grid::type          = "byspacing"
+Grid::dxyz          = 0.040
+Grid::avoid_origin  = "false"
+
+########################################
+
+#
+# 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"
+
+########################################
+
+#
+# ***** output *****
+#
+
+IOUtil::parfile_write    = "no"
+IOUtil::checkpoint_ID    = "yes"
+IOUtil::checkpoint_every = -1	# no checkpointing after the initial data
+
+IOBasic::outScalar_every = 1
+IOASCII::out1D_every     = 5
+IOBasic::outInfo_every   = 1
+
+IO::recover              = "no"
+IO::verbose              = "standard"
+
+IOBasic::outScalar_vars	= "ADMAnalysis::grr   ADMAonstraints::momx ADMAonstraints::ham"
+IOBasic::outInfo_vars	= " admconstraints::ham ADMBase::alp "
+IOASCII::out1D_vars	= " ADMAnalysis::grr   ADMBase::gxx ADMBase::gyy ADMBase::gzz   StaticConformal::psi   ADMConstraints::momx ADMConstraints::ham "
+
+########################################
+
+#
+# AHFinderDirect
+#
+
+AHFinderDirect::find_AHs_at_postinitial = "true"
+AHFinderDirect::find_AHs_at_poststep    = "false"
+AHFinderDirect::print_timing_stats = "true"
+
+AHFinderDirect::h_base_file_name       = "misner.h"
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+AHFinderDirect::origin_z[1] = 1.0
+AHFinderDirect::delta_drho_dsigma = 7.5
+
+AHFinderDirect::initial_guess_method = "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.0
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.3
diff --git a/par/misner-run.par b/par/misner-run.par
new file mode 100644
index 0000000..dbaacfb
--- /dev/null
+++ b/par/misner-run.par
@@ -0,0 +1,200 @@
+# This parameter file sets up Misner initial data and time-evolves
+# it with maximal slicing and gamma-driver shift.  It finds apparent
+# horizons with both AHFinderDirect and AHFinder, to allow comparison
+# of the two apparent horizon finders' results.
+#
+# Note that the grid is very high-resolution to ensure that both
+# apparent horizon finders find the horizons properly.  Even with
+# the fairly large number of grid points (the evolution takes around
+# 830 megabytes of memory on an AEI xeon), this means that the outer
+# boundary is very close in, so the boundary conditions are pretty
+# inaccurate.  But this doesn't matter for the purposes of this
+# test
+#
+
+########################################
+
+ActiveThorns = "BAM_Elliptic ADMBase ADMConstraints StaticConformal Ellbase CartGrid3D Time PUGH PUGHreduce PUGHslab IOAscii IOUtil IOBasic ADMCoupling ADMMacros ADMAnalysis SpaceMask CoordGauge ADM_BSSN Boundary PUGHInterp Maximal NaNChecker AHFinder LocalInterp AHFinderDirect"
+
+########################################
+
+#
+# physics
+#
+
+admbase::initial_data                  = "misner"
+bam_elliptic::misner_mu                = 2.0
+
+#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+ 
+#
+# grid
+#
+
+driver::global_nx                       = 83
+driver::global_ny                       = 83
+driver::global_nz                       = 116
+
+time::dtfac                             = 0.25
+
+grid::type                              = "byspacing"
+grid::avoid_origin                      = "true"
+grid::dxyz                              = 0.04
+grid::domain                            = "bitant"
+
+pugh::ghost_size                        = 2
+
+cactus::cctk_itlast                     = 2000	# t=20
+
+########################################
+
+#
+# output
+#
+
+IOBasic::outScalar_every                = 1
+IOASCII::out1D_every                    = 10
+IOBasic::outInfo_every                  = 1
+
+IO::recover                             = "no"
+IO::verbose                             = "standard"
+
+IOUtil::parfile_write = "no"
+IOBasic::outScalar_vars	= "admanalysis::grr admbase::alp admbase::betax admbase::betay admbase::betaz adm_bssn::adm_bs_gyy adm_bssn::adm_bs_k adm_bssn::adm_bs_gx adm_bssn::adm_bs_phi admconstraints::momx admconstraints::ham"
+IOBasic::outInfo_vars	= "admconstraints::ham admbase::alp"
+IOASCII::out1D_vars	= "admanalysis::grr ADMBase::gxx admbase::gyy ADMBase::gzz admbase::alp StaticConformal::psi admconstraints::momx admconstraints::ham"
+
+########################################
+
+NaNChecker::check_every     = 1
+NaNChecker::check_vars      = "admbase::lapse admbase::shift admbase::metric admbase::curv adm_bssn::adm_bssn_phi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_k"
+NaNChecker::action_if_found = "abort"
+
+########################################
+
+#
+# standard ADM_BSSN evolution
+#
+
+admbase::metric_type                   = "static conformal"
+staticconformal::conformal_storage     = "factor+derivs+2nd derivs"
+admbase::evolution_method              = "adm_bssn"
+
+adm_bssn::method                        = "icn"
+adm_bssn::ICN_itnum                     = 3
+adm_bssn::bound                         = "radiative"
+adm_bssn::stencil_size                  = 2
+adm_bssn::advection                     = "upwind2"
+boundary::radpower                      = 2
+
+########################################
+
+#
+# maximal slicing
+#
+
+admbase::initial_lapse                  = "ssapprox"
+admbase::lapse_evolution_method         = "maximal"
+adm_bssn::forceKzero                    = "yes"
+maximal::max_every                      = 5
+maximal::max_solver                     = "bam"
+maximal::max_bound                      = "robin"
+maximal::max_robin_falloff              = 1
+maximal::max_robin_inf                  = 1.0
+maximal::maximal_thresh                 = 1e-08
+bam_elliptic::maximal_maxcycles         = 250
+bam_elliptic::bam_laplacian             = "delginv"
+bam_elliptic::maximal_verbose           = "no"
+bam_elliptic::bam_fullV                 = "no"
+bam_elliptic::maximal_coarsetol         = 1e-10
+bam_elliptic::bam_solve_nmin            = 10
+bam_elliptic::bam_solve_nmax            = 500
+bam_elliptic::bam_solve_normmax         = 1e-10
+bam_elliptic::bam_persist               = "yes"
+
+bam_elliptic::bam_bound = "bamrobin"
+
+########################################
+
+#
+# gamma-driver shift
+#
+
+admbase::initial_shift                 = "zero"
+admbase::shift_evolution_method        = "gamma2"
+
+adm_bssn::ShiftGammaCoeff               = 0.75
+adm_bssn::ShiftAlpPower                 = 1
+adm_bssn::ShiftPsiPower                 = 4
+adm_bssn::ShiftPhiPower			= 0
+
+adm_bssn::BetaDriver                    = 2.8
+adm_bssn::BetaDriverPsi			= 0.0
+adm_bssn::BetaDriverPhi			= 0.0
+
+#########################################
+
+# AHFinder
+#
+
+AHFinder::ahf_active    = "yes"
+AHFinder::ahf_find3     = "no"
+AHFinder::ahf_flow      = "no"		# minimization algorithm
+AHFinder::ahf_verbose   = "no"
+AHFinder::ahf_veryverbose   = "no"
+AHFinder::ahf_guessverbose   = "yes"
+
+AHFinder::ahf_findevery = 20
+AHFinder::ahf_guessold  = "yes"
+
+AHFinder::ahf_2Doutput  = "yes"
+AHFinder::ahf_gaussout  = "yes"
+
+AHFinder::ahf_phi        = "no"		# axisymmetry
+
+AHFinder::ahf_refx       = "yes"	# reflection symmetry x <--> -x
+AHFinder::ahf_refy       = "yes"	# reflection symmetry y <--> -y
+
+AHFinder::ahf_lmax       = 10
+AHFinder::ahf_maxiter    = 15
+
+AHFinder::ahf_offset     = "yes"
+AHFinder::ahf_xc = 0.0
+AHFinder::ahf_yc = 0.0
+AHFinder::ahf_zc = 1.0
+AHFinder::ahf_wander     = "no"		# for easier cmp with AHFinderDirect
+
+# initial guess
+AHFinder::ahf_sloppyguess    = "yes"
+AHFinder::ahf_r0 = 0.3
+
+# surface integrals
+AHFinder::ahf_ntheta     = 50		# as low as 50
+
+########################################
+
+#
+# AHFinderDirect
+#
+
+AHFinderDirect::find_AHs_at_postinitial = "true"
+AHFinderDirect::find_AHs_at_poststep    = "true"
+AHFinderDirect::verbose_level = "algorithm highlights"
+
+AHFinderDirect::final_H_update_if_exit_x_H_small = "true"
+AHFinderDirect::h_base_file_name      = "misner.h"
+AHFinderDirect::how_often_to_output_h = 10	# time steps
+
+AHFinderDirect::N_horizons = 1
+AHFinderDirect::origin_x[1] = 0.0
+AHFinderDirect::origin_y[1] = 0.0
+AHFinderDirect::origin_z[1] = 1.0
+AHFinderDirect::delta_drho_dsigma = 7.5
+
+AHFinderDirect::initial_guess_method = "sphere"
+AHFinderDirect::initial_guess__sphere__x_center[1] = 0.0
+AHFinderDirect::initial_guess__sphere__y_center[1] = 0.0
+AHFinderDirect::initial_guess__sphere__z_center[1] = 1.0
+AHFinderDirect::initial_guess__sphere__radius[1] = 0.3
+
+########################################