From cc221a875a640142697df54a2bcff6638f7b1a88 Mon Sep 17 00:00:00 2001 From: jthorn Date: Sun, 16 Feb 2003 22:21:49 +0000 Subject: * general tidying of multiprocessor code * add current iteration number to what's sent to processor #0 for printing at each iteration * add a new verbosity level "algorithm debug" for what used to be #ifdef DEBUG debugging code * add new parameters & logic to set CCTK_Warn() levels for various warnings (eg AH is outside the gid, metric isn't +ve definite, etc) * restore the ability to output the horizon shape to a data file, which I mistakenly deleted a few days ago :( git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@947 f88db872-0e4f-0410-b76b-b9085cfa78c5 --- par/Kerr-tiny.par | 4 ++-- par/Kerr.par | 5 ++--- par/README | 5 +++-- par/misner-run.par | 18 +++++++++++++----- 4 files changed, 20 insertions(+), 12 deletions(-) (limited to 'par') diff --git a/par/Kerr-tiny.par b/par/Kerr-tiny.par index 3c91a90..697f486 100644 --- a/par/Kerr-tiny.par +++ b/par/Kerr-tiny.par @@ -11,10 +11,10 @@ cactus::cctk_itlast = 0 ActiveThorns = "PUGH" -driver::ghost_size = 2 +driver::ghost_size = 3 driver::global_nx = 31 driver::global_ny = 31 -driver::global_nz = 17 +driver::global_nz = 19 ActiveThorns = "CartGrid3D" grid::domain = "bitant" diff --git a/par/Kerr.par b/par/Kerr.par index c9e1077..59f24b7 100644 --- a/par/Kerr.par +++ b/par/Kerr.par @@ -8,10 +8,10 @@ cactus::cctk_itlast = 0 ActiveThorns = "PUGH" -driver::ghost_size = 2 +driver::ghost_size = 3 driver::global_nx = 31 driver::global_ny = 31 -driver::global_nz = 17 +driver::global_nz = 19 ActiveThorns = "CartGrid3D" grid::domain = "bitant" @@ -38,7 +38,6 @@ AHFinderDirect::find_AHs_at_postinitial = "true" AHFinderDirect::find_AHs_at_poststep = "false" AHFinderDirect::print_timing_stats = "true" -AHFinderDirect::final_Theta_update_if_Delta_h_converged = "true" AHFinderDirect::how_often_to_output_Theta = 1 AHFinderDirect::h_base_file_name = "Kerr.h" AHFinderDirect::Theta_base_file_name = "Kerr.Theta" diff --git a/par/README b/par/README index 0f3f353..91a7e27 100644 --- a/par/README +++ b/par/README @@ -7,8 +7,9 @@ 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. +misner1.2-025-init.par sets up Misner initial data (mu=1.2, so there's +a common horizon as well as two separate horizons) 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 diff --git a/par/misner-run.par b/par/misner-run.par index a5c0caf..1ddbf79 100644 --- a/par/misner-run.par +++ b/par/misner-run.par @@ -14,7 +14,7 @@ ######################################## -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" +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" ######################################## @@ -181,19 +181,27 @@ AHFinderDirect::find_AHs_at_postinitial = "true" AHFinderDirect::find_AHs_at_poststep = "true" AHFinderDirect::verbose_level = "algorithm highlights" -AHFinderDirect::final_Theta_update_if_Delta_h_converged = "true" AHFinderDirect::h_base_file_name = "misner.h" AHFinderDirect::how_often_to_output_h = 10 # time steps -AHFinderDirect::N_horizons = 1 +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[1] = "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