From abc0ce08d5e6314062fb7dbb915ca1b5b4242618 Mon Sep 17 00:00:00 2001 From: jthorn Date: Fri, 25 Jul 2003 12:05:35 +0000 Subject: test case for Michael Koppitz's bug where if one horizon was entirely outside the grid, we would spuriously miss finding all the others git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1144 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- par/bl-bad.par | 93 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 93 insertions(+) create mode 100644 par/bl-bad.par (limited to 'par') diff --git a/par/bl-bad.par b/par/bl-bad.par new file mode 100644 index 0000000..0a1976b --- /dev/null +++ b/par/bl-bad.par @@ -0,0 +1,93 @@ +# This parameter file sets up Brill-Lindquist initial data and finds the +# apparent horizons in it with AHFinderDirect. Notice that the initial +# guess for one of the AHs is entirely outside the grid. AHFinderDirect +# should find the other one ok, though... + +######################################## + +# +# ***** grid ***** +# + +ActiveThorns = "CartGrid3D PUGH" + +Driver::global_nx = 42 +Driver::global_ny = 42 +Driver::global_nz = 82 +PUGH::ghost_size = 3 + +Grid::domain = "full" +Grid::type = "byspacing" +Grid::dxyz = 0.05 +Grid::avoid_origin = "true" + +######################################## + +# +# physics +# + +ActiveThorns = "ADMBase ADMCoupling CoordGauge StaticConformal SpaceMask" +ActiveThorns = "IDAnalyticBH" + +ADMBase::initial_data = "bl_bh" +IDAnalyticBH::bl_nbh = 2 + +IDAnalyticBH::bl_m_1 = 1.0 +IDAnalyticBH::bl_x0_1 = 0.0 +IDAnalyticBH::bl_y0_1 = 0.0 +IDAnalyticBH::bl_z0_1 = 1.0 + +IDAnalyticBH::bl_m_2 = 1.0 +IDAnalyticBH::bl_x0_2 = 0.0 +IDAnalyticBH::bl_y0_2 = 0.0 +IDAnalyticBH::bl_z0_2 = -1.0 + +######################################## + +# +# 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" + +######################################## + +# +# AHFinderDirect +# + +ActiveThorns = "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::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.5 + +AHFinderDirect::origin_x[2] = 0.0 +AHFinderDirect::origin_y[2] = 0.0 +AHFinderDirect::origin_z[2] = -5.0 +AHFinderDirect::initial_guess_method[2] = "coordinate sphere" +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__y_center[2] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__z_center[2] = -5.0 +AHFinderDirect::initial_guess__coord_sphere__radius[2] = 0.5 -- cgit v1.2.3