concerning this item:
> Memory usage during regrid:
> * https://trac.einsteintoolkit.org/ticket/2194
> * still waiting for example parfile
i tried to add one such example to the ticket, but was rejected by the spam
filter. so i'm attaching
the parfile to this email.
cheers,
Miguel
--
Enviado a partir do meu dispositivo Android com o K-9 Mail. Peço desculpa pela
brevidade.
#------------------------------------------------------------------------------
ActiveThorns = "
ADMBase
ADMCoupling
ADMMacros
AEILocalInterp
AHFinderDirect
Boundary
Carpet
CarpetInterp
CarpetIOASCII
CarpetIOBasic
CarpetIOHDF5
CarpetIOScalar
CarpetLib
CarpetMask
CarpetReduce
CarpetRegrid2
CarpetSlab
CarpetTracker
CartGrid3D
CoordBase
CoordGauge
Dissipation
# Formaline
Fortran
GenericFD
GSL
HDF5
InitBase
IOUtil
LocalInterp
LoopControl
ML_ADMConstraints
ML_BSSN
ML_BSSN_Helper
MoL
Multipole
NaNChecker
NewRad
# PunctureTracker
QuasiLocalMeasures
ReflectionSymmetry
RotatingSymmetry180
Slab
SpaceMask
SphericalSurface
StaticConformal
SummationByParts
SymBase
SystemStatistics
TerminationTrigger
Time
TimerReport
TmunuBase
TwoPunctures
Vectors
WeylScal4
"
#------------------------------------------------------------------------------
# Grid setup
#------------------------------------------------------------------------------
CartGrid3D::type = "coordbase"
Carpet::domain_from_coordbase = yes
CoordBase::domainsize = "minmax"
# make sure all (xmax - xmin)/dx are integers!
CoordBase::xmin = 0.00
CoordBase::ymin = -120.00
CoordBase::zmin = 0.00
CoordBase::xmax = +120.00
CoordBase::ymax = +120.00
CoordBase::zmax = +120.00
CoordBase::dx = 2.00
CoordBase::dy = 2.00
CoordBase::dz = 2.00
driver::ghost_size = 3
CoordBase::boundary_size_x_lower = 3
CoordBase::boundary_size_y_lower = 3
CoordBase::boundary_size_z_lower = 3
CoordBase::boundary_size_x_upper = 3
CoordBase::boundary_size_y_upper = 3
CoordBase::boundary_size_z_upper = 3
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_y_lower = 0
CoordBase::boundary_shiftout_z_lower = 1
CarpetRegrid2::symmetry_rotating180 = yes
ReflectionSymmetry::reflection_x = no
ReflectionSymmetry::reflection_y = no
ReflectionSymmetry::reflection_z = yes
ReflectionSymmetry::avoid_origin_x = yes
ReflectionSymmetry::avoid_origin_y = yes
ReflectionSymmetry::avoid_origin_z = no
# Mesh refinement
#------------------------------------------------------------------------------
Carpet::max_refinement_levels = 7
CarpetRegrid2::num_centres = 3
CarpetRegrid2::num_levels_1 = 7
CarpetRegrid2::position_x_1 = +1.168642873
CarpetRegrid2::radius_1[ 1] = 64.0
CarpetRegrid2::radius_1[ 2] = 16.0
CarpetRegrid2::radius_1[ 3] = 8.0
CarpetRegrid2::radius_1[ 4] = 4.0
CarpetRegrid2::radius_1[ 5] = 2.0
CarpetRegrid2::radius_1[ 6] = 1.0
CarpetRegrid2::movement_threshold_1 = 0.16
CarpetRegrid2::num_levels_2 = 7
CarpetRegrid2::position_x_2 = -1.168642873
CarpetRegrid2::radius_2[ 1] = 64.0
CarpetRegrid2::radius_2[ 2] = 16.0
CarpetRegrid2::radius_2[ 3] = 8.0
CarpetRegrid2::radius_2[ 4] = 4.0
CarpetRegrid2::radius_2[ 5] = 2.0
CarpetRegrid2::radius_2[ 6] = 1.0
CarpetRegrid2::movement_threshold_2 = 0.16
CarpetRegrid2::num_levels_3 = 7
CarpetRegrid2::active_3 = no
CarpetRegrid2::radius_3[ 1] = 64.0
CarpetRegrid2::radius_3[ 2] = 16.0
CarpetRegrid2::radius_3[ 3] = 8.0
CarpetRegrid2::radius_3[ 4] = 4.0
CarpetRegrid2::radius_3[ 5] = 2.0
CarpetRegrid2::radius_3[ 6] = 1.0
CarpetRegrid2::movement_threshold_3 = 0.16
Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time = 2
CarpetRegrid2::freeze_unaligned_levels = yes
CarpetRegrid2::regrid_every = 128
CarpetRegrid2::verbose = no
Carpet::grid_structure_filename = "carpet-grid-structure"
Carpet::grid_coordinates_filename = "carpet-grid-coordinates"
Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128,
256]"
Time::dtfac = 0.25
# Initial Data
#------------------------------------------------------------------------------
ADMBase::metric_type = "physical"
ADMBase::initial_data = "twopunctures"
ADMBase::initial_lapse = "twopunctures-averaged"
ADMBase::initial_shift = "zero"
ADMBase::initial_dtlapse = "zero"
ADMBase::initial_dtshift = "zero"
ADMBase::lapse_timelevels = 3
ADMBase::shift_timelevels = 3
ADMBase::metric_timelevels = 3
# TwoPunctures::target_M_plus = 0.5
# TwoPunctures::target_M_minus = 0.5
TwoPunctures::par_m_plus = 0.453
TwoPunctures::par_m_minus = 0.453
TwoPunctures::par_b = 1.168642873
TwoPunctures::center_offset[0] = 0
TwoPunctures::par_P_plus[0] = 0.0
TwoPunctures::par_P_plus[1] = 0.3331917498
TwoPunctures::par_P_plus[2] = 0.
TwoPunctures::par_P_minus[0] = 0.0
TwoPunctures::par_P_minus[1] = -0.3331917498
TwoPunctures::par_P_minus[2] = 0.
TwoPunctures::par_S_plus[0] = 0.
TwoPunctures::par_S_plus[1] = 0.
TwoPunctures::par_S_plus[2] = 0.
TwoPunctures::par_S_minus[0] = 0.
TwoPunctures::par_S_minus[1] = 0.
TwoPunctures::par_S_minus[2] = 0.
# TwoPunctures::give_bare_mass = yes
TwoPunctures::TP_epsilon = 1.0e-2
TwoPunctures::TP_Tiny = 1.0e-2
TwoPunctures::npoints_A = 30
TwoPunctures::npoints_B = 30
TwoPunctures::npoints_phi = 16
TwoPunctures::grid_setup_method = "evaluation"
TwoPunctures::keep_u_around = yes
TwoPunctures::verbose = yes
InitBase::initial_data_setup_method = "init_some_levels"
Carpet::init_fill_timelevels = yes
Carpet::init_3_timelevels = no
# Evolution
#------------------------------------------------------------------------------
ADMBase::evolution_method = "ML_BSSN"
ADMBase::lapse_evolution_method = "ML_BSSN"
ADMBase::shift_evolution_method = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"
ML_BSSN::harmonicN = 1 # 1+log
ML_BSSN::harmonicF = 2.0 # 1+log
ML_BSSN::ShiftGammaCoeff = 0.75
ML_BSSN::BetaDriver = 1.0
ML_BSSN::advectLapse = 1
ML_BSSN::advectShift = 1
ML_BSSN::MinimumLapse = 1.0e-8
ML_BSSN::conformalMethod = 1 # 1 for W
# Boundary conditions
#------------------------------------------------------------------------------
ML_BSSN::initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::rhs_boundary_condition = "NewRad"
Boundary::radpower = 2
ML_BSSN::ML_log_confac_bound = "none"
ML_BSSN::ML_metric_bound = "none"
ML_BSSN::ML_Gamma_bound = "none"
ML_BSSN::ML_trace_curv_bound = "none"
ML_BSSN::ML_curv_bound = "none"
ML_BSSN::ML_lapse_bound = "none"
ML_BSSN::ML_dtlapse_bound = "none"
ML_BSSN::ML_shift_bound = "none"
ML_BSSN::ML_dtshift_bound = "none"
# Spatial finite differencing
#------------------------------------------------------------------------------
SummationByParts::order = 4
ML_BSSN::fdOrder = 4
Dissipation::order = 5
Dissipation::vars = "
ML_BSSN::ML_log_confac
ML_BSSN::ML_metric
ML_BSSN::ML_trace_curv
ML_BSSN::ML_curv
ML_BSSN::ML_Gamma
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
ML_BSSN::ML_dtlapse
ML_BSSN::ML_dtshift
"
# Integration method
#------------------------------------------------------------------------------
MoL::ODE_Method = "RK4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1
Carpet::num_integrator_substeps = 4
# Spherical surfaces
#------------------------------------------------------------------------------
SphericalSurface::nsurfaces = 3
SphericalSurface::maxntheta = 66
SphericalSurface::maxnphi = 124
SphericalSurface::verbose = no
# Horizon 1
SphericalSurface::ntheta [0] = 41
SphericalSurface::nphi [0] = 80
SphericalSurface::nghoststheta [0] = 2
SphericalSurface::nghostsphi [0] = 2
CarpetMask::excluded_surface [0] = 0
CarpetMask::excluded_surface_factor [0] = 1.0
# Horizon 2
SphericalSurface::ntheta [1] = 41
SphericalSurface::nphi [1] = 80
SphericalSurface::nghoststheta [1] = 2
SphericalSurface::nghostsphi [1] = 2
CarpetMask::excluded_surface [1] = 1
CarpetMask::excluded_surface_factor [1] = 1.0
# Common horizon
SphericalSurface::ntheta [2] = 41
SphericalSurface::nphi [2] = 80
SphericalSurface::nghoststheta [2] = 2
SphericalSurface::nghostsphi [2] = 2
CarpetMask::excluded_surface [2] = 2
CarpetMask::excluded_surface_factor [2] = 1.0
CarpetMask::verbose = no
# Puncture tracking
#------------------------------------------------------------------------------
CarpetTracker::surface[0] = 0
CarpetTracker::surface[1] = 1
CarpetTracker::surface[2] = 2
# PunctureTracker::track [0] = yes
# PunctureTracker::initial_x [0] = 3
# PunctureTracker::which_surface_to_store_info[0] = 0
# PunctureTracker::track [1] = yes
# PunctureTracker::initial_x [1] = -3
# PunctureTracker::which_surface_to_store_info[1] = 1
# PunctureTracker::verbose = no
# Wave extraction
#------------------------------------------------------------------------------
WeylScal4::fd_order = "4th"
WeylScal4::offset = 1e-8
Multipole::nradii = 4
Multipole::radius[0] = 30
Multipole::radius[1] = 50
Multipole::radius[2] = 80
Multipole::radius[3] = 100
Multipole::ntheta = 120
Multipole::nphi = 240
Multipole::variables = "
WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='Psi4'}
"
# Multipole::integration_method = "Simpson"
Multipole::out_every = 64
Multipole::l_max = 4
Multipole::output_hdf5 = no
Multipole::output_ascii = yes
# Horizons
#------------------------------------------------------------------------------
# AHFinderDirect::verbose_level = "algorithm
highlights"
AHFinderDirect::verbose_level = "physics details"
AHFinderDirect::output_BH_diagnostics = "true"
AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
AHFinderDirect::N_horizons = 3
AHFinderDirect::find_every = 128
AHFinderDirect::output_h_every = 0
AHFinderDirect::max_Newton_iterations__initial = 50
AHFinderDirect::max_Newton_iterations__subsequent = 50
AHFinderDirect::max_allowable_Theta_growth_iterations = 10
AHFinderDirect::max_allowable_Theta_nonshrink_iterations = 10
AHFinderDirect::geometry_interpolator_name = "Lagrange polynomial
interpolation"
AHFinderDirect::geometry_interpolator_pars = "order=4"
AHFinderDirect::surface_interpolator_name = "Lagrange polynomial
interpolation"
AHFinderDirect::surface_interpolator_pars = "order=4"
AHFinderDirect::move_origins = yes
AHFinderDirect::reshape_while_moving = yes
AHFinderDirect::predict_origin_movement = yes
AHFinderDirect::origin_x [1] = +1.168642873
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +1.168642873
AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.25
AHFinderDirect::which_surface_to_store_info [1] = 0
AHFinderDirect::set_mask_for_individual_horizon [1] = no
AHFinderDirect::reset_horizon_after_not_finding [1] = no
# AHFinderDirect::track_origin_from_grid_scalar [1] = yes
# AHFinderDirect::track_origin_source_x [1] =
"PunctureTracker::pt_loc_x[0]"
# AHFinderDirect::track_origin_source_y [1] =
"PunctureTracker::pt_loc_y[0]"
# AHFinderDirect::track_origin_source_z [1] =
"PunctureTracker::pt_loc_z[0]"
AHFinderDirect::max_allowable_horizon_radius [1] = 3
AHFinderDirect::dont_find_after_individual_time [1] = 30.0
AHFinderDirect::origin_x [2] = -1.168642873
AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -1.168642873
AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.25
AHFinderDirect::which_surface_to_store_info [2] = 1
AHFinderDirect::set_mask_for_individual_horizon [2] = no
AHFinderDirect::reset_horizon_after_not_finding [2] = no
# AHFinderDirect::track_origin_from_grid_scalar [2] = yes
# AHFinderDirect::track_origin_source_x [2] =
"PunctureTracker::pt_loc_x[1]"
# AHFinderDirect::track_origin_source_y [2] =
"PunctureTracker::pt_loc_y[1]"
# AHFinderDirect::track_origin_source_z [2] =
"PunctureTracker::pt_loc_z[1]"
AHFinderDirect::max_allowable_horizon_radius [2] = 3
AHFinderDirect::dont_find_after_individual_time [2] = 30.0
AHFinderDirect::origin_x [3] = 0
AHFinderDirect::find_after_individual_time [3] = 15.0
AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0
AHFinderDirect::initial_guess__coord_sphere__radius [3] = 1.0
AHFinderDirect::which_surface_to_store_info [3] = 2
AHFinderDirect::set_mask_for_individual_horizon [3] = no
AHFinderDirect::max_allowable_horizon_radius [3] = 6
# Isolated Horizons
#-------------------------------------------------------------------------------
QuasiLocalMeasures::verbose = yes
QuasiLocalMeasures::veryverbose = no
QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation"
QuasiLocalMeasures::interpolator_options = "order=4"
QuasiLocalMeasures::spatial_order = 4
QuasiLocalMeasures::num_surfaces = 3
QuasiLocalMeasures::surface_index [0] = 0
QuasiLocalMeasures::surface_index [1] = 1
QuasiLocalMeasures::surface_index [2] = 2
# Check for NaNs
#-------------------------------------------------------------------------------
Carpet::poison_new_timelevels = yes
CarpetLib::poison_new_memory = yes
Carpet::check_for_poison = no
NaNChecker::check_every = 512
NanChecker::check_after = 0
NaNChecker::report_max = 10
# NaNChecker::verbose = "all"
NaNChecker::action_if_found = "terminate"
NaNChecker::out_NaNmask = yes
NaNChecker::check_vars = "
ML_BSSN::ML_log_confac
"
# Timers
#-------------------------------------------------------------------------------
Cactus::cctk_timer_output = "full"
TimerReport::out_every = 5120
TimerReport::n_top_timers = 40
TimerReport::output_all_timers_together = yes
TimerReport::output_all_timers_readable = yes
TimerReport::output_schedule_timers = no
# I/O thorns
#-------------------------------------------------------------------------------
Cactus::cctk_run_title = $parfile
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOASCII::one_file_per_group = yes
IOHDF5::use_checksums = no
IOHDF5::one_file_per_group = no
IOBasic::outInfo_every = 4
IOBasic::outInfo_reductions = "minimum maximum"
IOBasic::outInfo_vars = "
Carpet::physical_time_per_hour
ML_BSSN::ML_log_confac
SystemStatistics::maxrss_mb
"
# for scalar reductions of 3D grid functions
IOScalar::outScalar_every = 64
IOScalar::outScalar_reductions = "minimum maximum average norm2"
IOScalar::outScalar_vars = "
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
SystemStatistics::maxrss_mb
SystemStatistics::process_memory_mb
Carpet::memory_procs
"
# output just at one point (0D)
IOASCII::out0D_every = 32
IOASCII::out0D_vars = "
Carpet::timing
# PunctureTracker::pt_loc
QuasiLocalMeasures::qlm_scalars{out_every = 128}
"
IOASCII::output_symmetry_points = no
IOASCII::out3D_ghosts = no
# 1D text output
IOASCII::out1D_every = 128
IOASCII::out1D_d = no
IOASCII::out1D_x = yes
IOASCII::out1D_y = no
IOASCII::out1D_z = no
IOASCII::out1D_vars = "
ADMBase::lapse
ML_BSSN::ML_Ham
ML_BSSN::ML_mom
"
# 1D HDF5 output
#IOHDF5::out1D_every = 256
#IOHDF5::out1D_d = no
#IOHDF5::out1D_x = yes
#IOHDF5::out1D_y = no
#IOHDF5::out1D_z = no
#IOHDF5::out1D_vars = "
# ADMBase::lapse
#"
# # 2D HDF5 output
# IOHDF5::out2D_every = 256
# IOHDF5::out2D_xy = yes
# IOHDF5::out2D_xz = no
# IOHDF5::out2D_yz = no
# IOHDF5::out2D_vars = "
# ADMBase::lapse
# "
# # 3D HDF5 output
# IOHDF5::out_every = 8192
# IOHDF5::out_vars = "
# ADMBase::lapse
# "
Carpet::verbose = no
Carpet::veryverbose = no
Carpet::schedule_barriers = no
Carpet::storage_verbose = no
CarpetLib::output_bboxes = no
Cactus::cctk_full_warnings = yes
Cactus::highlight_warning_messages = no
# Checkpointing and recovery
#-------------------------------------------------------------------------------
CarpetIOHDF5::checkpoint = yes
IO::checkpoint_dir = "../checkpoints"
IO::checkpoint_ID = no
IO::checkpoint_every_walltime_hours = 23
IO::checkpoint_on_terminate = yes
IO::out_proc_every = 2
IO::checkpoint_keep = 1
IO::recover = "autoprobe"
IO::recover_dir = "../checkpoints"
IO::abort_on_io_errors = yes
CarpetIOHDF5::open_one_input_file_at_a_time = yes
CarpetIOHDF5::compression_level = 9
# Run termination
#-------------------------------------------------------------------------------
TerminationTrigger::max_walltime = 72 # hours
TerminationTrigger::on_remaining_walltime = 60 # minutes
TerminationTrigger::output_remtime_every_minutes = 60
TerminationTrigger::termination_from_file = yes
TerminationTrigger::termination_file = "terminate.txt"
TerminationTrigger::create_termination_file = yes
Cactus::terminate = "time"
Cactus::cctk_final_time = 300.0
_______________________________________________
Users mailing list
[email protected]
http://lists.einsteintoolkit.org/mailman/listinfo/users
