From c5938832b234565efde9808b7016d78b3035ef91 Mon Sep 17 00:00:00 2001 From: jthorn Date: Wed, 19 Nov 2003 16:53:29 +0000 Subject: make max_allowable_horizon_radius into an array parameter so it can/must be specified for each horizon --> gives a bit more flexibility for screening out "bad" cases (alas this isn't in the thornguide yet :() git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1205 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- par/misner-max-allowable-horizon-radius=0.75.par | 97 ++++++++++++++++++++++++ par/misner-max-allowable-horizon-radius=0.8.par | 97 ++++++++++++++++++++++++ 2 files changed, 194 insertions(+) create mode 100644 par/misner-max-allowable-horizon-radius=0.75.par create mode 100644 par/misner-max-allowable-horizon-radius=0.8.par (limited to 'par') diff --git a/par/misner-max-allowable-horizon-radius=0.75.par b/par/misner-max-allowable-horizon-radius=0.75.par new file mode 100644 index 0000000..da8e1c3 --- /dev/null +++ b/par/misner-max-allowable-horizon-radius=0.75.par @@ -0,0 +1,97 @@ +# 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 = "CoordBase 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::N_horizons = 2 + +AHFinderDirect::max_allowable_horizon_radius[1] = 1.0 # no problem +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::max_allowable_horizon_radius[2] = 0.75 # should prevent this + # AH from being found +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 diff --git a/par/misner-max-allowable-horizon-radius=0.8.par b/par/misner-max-allowable-horizon-radius=0.8.par new file mode 100644 index 0000000..7075965 --- /dev/null +++ b/par/misner-max-allowable-horizon-radius=0.8.par @@ -0,0 +1,97 @@ +# 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 = "CoordBase 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::N_horizons = 2 + +AHFinderDirect::max_allowable_horizon_radius[1] = 1.0 # no problem +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::max_allowable_horizon_radius[2] = 0.8 # should still allow + # finding this AH +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 -- cgit v1.2.3