diff options
| -rw-r--r-- | README | 9 | ||||
| -rw-r--r-- | par/misner1.2-025-init.par | 3 | ||||
| -rw-r--r-- | test/Kerr.par | 18 | ||||
| -rw-r--r-- | test/Kerr/BH_diagnostics.ah1.gp | 26 |
4 files changed, 30 insertions, 26 deletions
@@ -61,6 +61,9 @@ This thorn can use either of two different interpolator APIs * CCTK_InterpLocalUniform() interpolator API, which at present is only supplied by the CactusBase/InterpLocal thorn. +This thorn uses various Cactus reduction operations, so you'll need +a (driver-specific) reduction thorn (e.g. CactusPUGH/PUGHReduce). + This thorn is written in C++, so you'll need a C++ compiler -- in fact a fairly modern one -- to compile this thorn. See the "Compilation Notes" section below for details. @@ -151,9 +154,11 @@ IBM Power4 (AIX) Please let me (Jonathan) know of any other compilers which are ok. +It's uncertain whether this thorn will work or not using +* gcc version 2.96 20000731 (Red Hat Linux 7.1 2.96-98) + (This is the system gcc on Red Hat GNU/Linux 7.2 systems) + The code will almost certainly *not* compile using * any "egcs" version of gcc, or more generally, any version of gcc earlier than 2.95 -* gcc version 2.96 20000731 (Red Hat Linux 7.1 2.96-98) - (This is the system gcc on Red Hat GNU/Linux 7.2 systems) * any version of Microsoft Visual C++ up to and including mid-2002 diff --git a/par/misner1.2-025-init.par b/par/misner1.2-025-init.par index 42fe9f9..87b01f9 100644 --- a/par/misner1.2-025-init.par +++ b/par/misner1.2-025-init.par @@ -17,7 +17,7 @@ ActiveThorns = "CartGrid3D PUGH" Driver::global_nx = 83 Driver::global_ny = 83 Driver::global_nz = 115 -PUGH::ghost_size = 2 +PUGH::ghost_size = 20 Grid::domain = "octant" Grid::type = "byspacing" @@ -68,6 +68,7 @@ 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-init.h" diff --git a/test/Kerr.par b/test/Kerr.par index 3f51f3b..ad70284 100644 --- a/test/Kerr.par +++ b/test/Kerr.par @@ -8,10 +8,10 @@ cactus::cctk_itlast = 0 ActiveThorns = "PUGH" -driver::ghost_size = 2 -driver::global_nx = 41 -driver::global_ny = 41 -driver::global_nz = 23 +driver::ghost_size = 3 +driver::global_nx = 31 +driver::global_ny = 31 +driver::global_nz = 19 ActiveThorns = "CartGrid3D" grid::domain = "bitant" @@ -37,24 +37,20 @@ IO::parfile_write = "no" ######################################## -ActiveThorns = "LocalInterp PUGHInterp AHFinderDirect" +ActiveThorns = "LocalInterp PUGHInterp PUGHReduce AHFinderDirect" 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::h_base_file_name = "Kerr/Kerr.h" +AHFinderDirect::h_base_file_name = "Kerr/h" AHFinderDirect::BH_diagnostics_base_file_name = "Kerr/BH_diagnostics" - -AHFinderDirect::Theta_norm_for_convergence = 1.0e-10 -AHFinderDirect::Delta_h_norm_for_convergence = 1.0e-10 +AHFinderDirect::Theta_norm_for_convergence = 1.0e-12 AHFinderDirect::N_horizons = 1 AHFinderDirect::origin_x[1] = 0.5 AHFinderDirect::origin_y[1] = 0.7 AHFinderDirect::origin_z[1] = 0.0 -AHFinderDirect::N_zones_per_right_angle[1] = 18 AHFinderDirect::initial_guess_method[1] = "coordinate sphere" AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -0.2 diff --git a/test/Kerr/BH_diagnostics.ah1.gp b/test/Kerr/BH_diagnostics.ah1.gp index 5cfef85..8702764 100644 --- a/test/Kerr/BH_diagnostics.ah1.gp +++ b/test/Kerr/BH_diagnostics.ah1.gp @@ -1,14 +1,16 @@ # apparent horizon 1/1 # -# column 1 = cctk_iteration -# column 2 = cctk_time -# column 3 = centroid_x -# column 4 = centroid_y -# column 5 = centroid_z -# column 6 = mean radius -# column 7 = xy-plane circumference -# column 8 = xz-plane circumference -# column 9 = yz-plane circumference -# column 10 = area -# column 11 = m_irreducible -0 0.000 0.000039 0.000009 0.000000 1.996522276 12.56580486 10.83709599 11.22743101 45.23550768 0.9486473693 +# column 1 = cctk_iteration +# column 2 = cctk_time +# column 3 = centroid_x +# column 4 = centroid_y +# column 5 = centroid_z +# column 6 = min radius +# column 7 = max radius +# column 8 = mean radius +# column 9 = xy-plane circumference +# column 10 = xz-plane circumference +# column 11 = yz-plane circumference +# column 12 = area +# column 13 = m_irreducible +0 0.000 0.000039 0.000009 0.000000 1.037097639 2.757422636 1.996522276 12.56580486 10.83709599 11.22743101 45.23550768 0.9486473693 |