From 9acf611ea96b479aee04f4e6dfb6dd30d440b3c1 Mon Sep 17 00:00:00 2001 From: jthorn Date: Thu, 20 Feb 2003 23:45:00 +0000 Subject: src/driver/Newton.cc * tweak printing for dummy horizons git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@948 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- run/test-ahfinderdirect/Misner/README | 1 + .../Misner/misner1.2-025-N30-guess300/misner.par | 110 +++++++++++++++++++++ .../misner1.2-025-guess300.par | 109 ++++++++++++++++++++ .../Misner/misner1.2-025-guess400/misner.par | 109 ++++++++++++++++++++ 4 files changed, 329 insertions(+) create mode 100644 run/test-ahfinderdirect/Misner/README create mode 100644 run/test-ahfinderdirect/Misner/misner1.2-025-N30-guess300/misner.par create mode 100644 run/test-ahfinderdirect/Misner/misner1.2-025-guess300/misner1.2-025-guess300.par create mode 100644 run/test-ahfinderdirect/Misner/misner1.2-025-guess400/misner.par (limited to 'run') diff --git a/run/test-ahfinderdirect/Misner/README b/run/test-ahfinderdirect/Misner/README new file mode 100644 index 0000000..fec0cfa --- /dev/null +++ b/run/test-ahfinderdirect/Misner/README @@ -0,0 +1 @@ +These are some tests of apparent horizon finding in Misner initial data. diff --git a/run/test-ahfinderdirect/Misner/misner1.2-025-N30-guess300/misner.par b/run/test-ahfinderdirect/Misner/misner1.2-025-N30-guess300/misner.par new file mode 100644 index 0000000..09f9612 --- /dev/null +++ b/run/test-ahfinderdirect/Misner/misner1.2-025-N30-guess300/misner.par @@ -0,0 +1,110 @@ +# 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. + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "CartGrid3D PUGH" + +Driver::global_nx = 99 +Driver::global_ny = 99 +Driver::global_nz = 115 +PUGH::ghost_size = 6 + +Grid::domain = "octant" +Grid::type = "byspacing" +Grid::dxyz = 0.025 +Grid::avoid_origin = "false" + +######################################## + +# +# physics +# + +ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal 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" + +######################################## + +# +# ***** output ***** +# +ActiveThorns = "PUGHSlab IOBasic IOUtil IOASCII" + +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 +# + +ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect" + +AHFinderDirect::find_AHs_at_postinitial = "true" +AHFinderDirect::find_AHs_at_poststep = "false" +AHFinderDirect::print_timing_stats = "true" +##AHFinderDirect::verbose_level = "algorithm details" + +AHFinderDirect::h_base_file_name = "misner.h" + +AHFinderDirect::N_horizons = 2 + +AHFinderDirect::N_zones_per_right_angle[1] = 30 +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.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.3 + +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.8 +AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.8 +AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.4 diff --git a/run/test-ahfinderdirect/Misner/misner1.2-025-guess300/misner1.2-025-guess300.par b/run/test-ahfinderdirect/Misner/misner1.2-025-guess300/misner1.2-025-guess300.par new file mode 100644 index 0000000..54af2ba --- /dev/null +++ b/run/test-ahfinderdirect/Misner/misner1.2-025-guess300/misner1.2-025-guess300.par @@ -0,0 +1,109 @@ +# 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. + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "CartGrid3D PUGH" + +Driver::global_nx = 99 +Driver::global_ny = 99 +Driver::global_nz = 115 +PUGH::ghost_size = 6 + +Grid::domain = "octant" +Grid::type = "byspacing" +Grid::dxyz = 0.025 +Grid::avoid_origin = "false" + +######################################## + +# +# physics +# + +ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal 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" + +######################################## + +# +# ***** output ***** +# +ActiveThorns = "PUGHSlab IOBasic IOUtil IOASCII" + +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 +# + +ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect" + +AHFinderDirect::find_AHs_at_postinitial = "true" +AHFinderDirect::find_AHs_at_poststep = "false" +AHFinderDirect::print_timing_stats = "true" +##AHFinderDirect::verbose_level = "algorithm details" + +AHFinderDirect::h_base_file_name = "misner1.2-025-guess300.h" + +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.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.3 + +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.8 +AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.8 +AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.4 diff --git a/run/test-ahfinderdirect/Misner/misner1.2-025-guess400/misner.par b/run/test-ahfinderdirect/Misner/misner1.2-025-guess400/misner.par new file mode 100644 index 0000000..33a2795 --- /dev/null +++ b/run/test-ahfinderdirect/Misner/misner1.2-025-guess400/misner.par @@ -0,0 +1,109 @@ +# 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. + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "CartGrid3D PUGH" + +Driver::global_nx = 99 +Driver::global_ny = 99 +Driver::global_nz = 115 +PUGH::ghost_size = 6 + +Grid::domain = "octant" +Grid::type = "byspacing" +Grid::dxyz = 0.025 +Grid::avoid_origin = "false" + +######################################## + +# +# physics +# + +ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal 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" + +######################################## + +# +# ***** output ***** +# +ActiveThorns = "PUGHSlab IOBasic IOUtil IOASCII" + +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 +# + +ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect" + +AHFinderDirect::find_AHs_at_postinitial = "true" +AHFinderDirect::find_AHs_at_poststep = "false" +AHFinderDirect::print_timing_stats = "true" +##AHFinderDirect::verbose_level = "algorithm details" + +AHFinderDirect::h_base_file_name = "misner.h" + +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.0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 0.4 + +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.8 +AHFinderDirect::initial_guess__coord_ellipsoid__y_radius[2] = 0.8 +AHFinderDirect::initial_guess__coord_ellipsoid__z_radius[2] = 1.4 -- cgit v1.2.3