From bf34a0ebe6720a208ca3450d140a35d0f7b0997d Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Sun, 27 Dec 2009 23:56:35 +0100 Subject: Update example parameter files: - remove outdated examples - ensure the vacuum example works - add matter example --- par/bl_bssn_1.par | 171 ----------------- par/bl_bssn_2.par | 171 ----------------- par/gw_bssn_1.par | 206 --------------------- par/gw_bssn_2.par | 206 --------------------- par/head_on.par | 368 ------------------------------------- par/ks-mclachlan.par | 400 ++++++++++++++++++++++++++++++++++++++++ par/ks_bssn_0_8th.par | 172 ------------------ par/ks_bssn_1.par | 173 ------------------ par/ks_bssn_1_8th.par | 172 ------------------ par/ks_bssn_2.par | 173 ------------------ par/ks_bssn_2_8th.par | 172 ------------------ par/ks_bssn_3.par | 172 ------------------ par/ks_bssn_evol.par | 327 --------------------------------- par/ks_ccatie_1.par | 152 ---------------- par/m_bssn.par | 74 -------- par/old/bl-mclachlan-ab4.par | 380 ++++++++++++++++++++++++++++++++++++++ par/old/bl_bssn_1.par | 171 +++++++++++++++++ par/old/bl_bssn_2.par | 171 +++++++++++++++++ par/old/gw_bssn_1.par | 206 +++++++++++++++++++++ par/old/gw_bssn_2.par | 206 +++++++++++++++++++++ par/old/head_on.par | 368 +++++++++++++++++++++++++++++++++++++ par/old/ks_bssn_0_8th.par | 172 ++++++++++++++++++ par/old/ks_bssn_1.par | 173 ++++++++++++++++++ par/old/ks_bssn_1_8th.par | 172 ++++++++++++++++++ par/old/ks_bssn_2.par | 173 ++++++++++++++++++ par/old/ks_bssn_2_8th.par | 172 ++++++++++++++++++ par/old/ks_bssn_3.par | 172 ++++++++++++++++++ par/old/ks_bssn_evol.par | 327 +++++++++++++++++++++++++++++++++ par/old/ks_bssn_evol_mp.par | 309 +++++++++++++++++++++++++++++++ par/old/ks_bssn_evol_mp_hi.par | 309 +++++++++++++++++++++++++++++++ par/old/ks_ccatie_1.par | 152 ++++++++++++++++ par/old/m_bssn.par | 74 ++++++++ par/qc0-mclachlan.par | 297 ++++++++---------------------- par/tov-mclachlan.par | 404 +++++++++++++++++++++++++++++++++++++++++ 34 files changed, 4584 insertions(+), 2933 deletions(-) delete mode 100644 par/bl_bssn_1.par delete mode 100644 par/bl_bssn_2.par delete mode 100644 par/gw_bssn_1.par delete mode 100644 par/gw_bssn_2.par delete mode 100644 par/head_on.par create mode 100644 par/ks-mclachlan.par delete mode 100644 par/ks_bssn_0_8th.par delete mode 100644 par/ks_bssn_1.par delete mode 100644 par/ks_bssn_1_8th.par delete mode 100644 par/ks_bssn_2.par delete mode 100644 par/ks_bssn_2_8th.par delete mode 100644 par/ks_bssn_3.par delete mode 100644 par/ks_bssn_evol.par delete mode 100644 par/ks_ccatie_1.par delete mode 100644 par/m_bssn.par create mode 100644 par/old/bl-mclachlan-ab4.par create mode 100644 par/old/bl_bssn_1.par create mode 100644 par/old/bl_bssn_2.par create mode 100644 par/old/gw_bssn_1.par create mode 100644 par/old/gw_bssn_2.par create mode 100644 par/old/head_on.par create mode 100644 par/old/ks_bssn_0_8th.par create mode 100644 par/old/ks_bssn_1.par create mode 100644 par/old/ks_bssn_1_8th.par create mode 100644 par/old/ks_bssn_2.par create mode 100644 par/old/ks_bssn_2_8th.par create mode 100644 par/old/ks_bssn_3.par create mode 100644 par/old/ks_bssn_evol.par create mode 100644 par/old/ks_bssn_evol_mp.par create mode 100644 par/old/ks_bssn_evol_mp_hi.par create mode 100644 par/old/ks_ccatie_1.par create mode 100644 par/old/m_bssn.par create mode 100644 par/tov-mclachlan.par (limited to 'par') diff --git a/par/bl_bssn_1.par b/par/bl_bssn_1.par deleted file mode 100644 index 862e994..0000000 --- a/par/bl_bssn_1.par +++ /dev/null @@ -1,171 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.2 -CoordBase::dy = 0.2 -CoordBase::dz = 0.2 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Schwarzschild/BL" -Exact::Schwarzschild_BL__mass = 1.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -ML_BSSN::log_confac_bound = "none" -ML_BSSN::metric_bound = "none" -ML_BSSN::Gamma_bound = "none" -ML_BSSN::trace_curv_bound = "none" -ML_BSSN::curv_bound = "none" -ML_BSSN::lapse_bound = "none" -ML_BSSN::dtlapse_bound = "none" -ML_BSSN::shift_bound = "none" -ML_BSSN::dtshift_bound = "none" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" diff --git a/par/bl_bssn_2.par b/par/bl_bssn_2.par deleted file mode 100644 index caefed4..0000000 --- a/par/bl_bssn_2.par +++ /dev/null @@ -1,171 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.1 -CoordBase::dy = 0.1 -CoordBase::dz = 0.1 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Schwarzschild/BL" -Exact::Schwarzschild_BL__mass = 1.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -ML_BSSN::log_confac_bound = "none" -ML_BSSN::metric_bound = "none" -ML_BSSN::Gamma_bound = "none" -ML_BSSN::trace_curv_bound = "none" -ML_BSSN::curv_bound = "none" -ML_BSSN::lapse_bound = "none" -ML_BSSN::dtlapse_bound = "none" -ML_BSSN::shift_bound = "none" -ML_BSSN::dtshift_bound = "none" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" diff --git a/par/gw_bssn_1.par b/par/gw_bssn_1.par deleted file mode 100644 index 10d9255..0000000 --- a/par/gw_bssn_1.par +++ /dev/null @@ -1,206 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - Time - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - Periodic - - Dissipation - Exact - Formaline - GSL - SphericalSurface - TmunuBase - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 100 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 1.0 -CoordBase::ymax = 0.0 -CoordBase::zmax = 0.0 -CoordBase::dx = 0.02 -CoordBase::dy = 0.02 -CoordBase::dz = 0.02 - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 1 -CoordBase::boundary_shiftout_z_upper = 1 - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - -Periodic::periodic = yes - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "exact" -ADMBase::initial_dtshift = "exact" - -Exact::exact_model = "Minkowski/gauge wave" -Exact::Minkowski_gauge_wave__what_fn = "expsin" -Exact::Minkowski_gauge_wave__amplitude = 0.01 - -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" - -ADMMacros::spatial_order = 4 - - - -Time::dtfac = 0.5 - -MoL::ODE_Method = "RK4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -Dissipation::epsdis = 0.25 -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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift -" - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - 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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs - ML_BSSN::ML_lapserhs - ML_BSSN::ML_dtlapserhs - ML_BSSN::ML_shiftrhs - ML_BSSN::ML_dtshiftrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - 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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs - ML_BSSN::ML_lapserhs - ML_BSSN::ML_dtlapserhs - ML_BSSN::ML_shiftrhs - ML_BSSN::ML_dtshiftrhs -" diff --git a/par/gw_bssn_2.par b/par/gw_bssn_2.par deleted file mode 100644 index ff09c74..0000000 --- a/par/gw_bssn_2.par +++ /dev/null @@ -1,206 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - Time - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - Periodic - - Dissipation - Exact - Formaline - GSL - SphericalSurface - TmunuBase - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 200 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 1.0 -CoordBase::ymax = 0.0 -CoordBase::zmax = 0.0 -CoordBase::dx = 0.01 -CoordBase::dy = 0.01 -CoordBase::dz = 0.01 - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 1 -CoordBase::boundary_shiftout_z_upper = 1 - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - -Periodic::periodic = yes - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "exact" -ADMBase::initial_dtshift = "exact" - -Exact::exact_model = "Minkowski/gauge wave" -Exact::Minkowski_gauge_wave__what_fn = "expsin" -Exact::Minkowski_gauge_wave__amplitude = 0.01 - -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" - -ADMMacros::spatial_order = 4 - - - -Time::dtfac = 0.5 - -MoL::ODE_Method = "RK4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -Dissipation::epsdis = 0.25 -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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift -" - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - 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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs - ML_BSSN::ML_lapserhs - ML_BSSN::ML_dtlapserhs - ML_BSSN::ML_shiftrhs - ML_BSSN::ML_dtshiftrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - 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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs - ML_BSSN::ML_lapserhs - ML_BSSN::ML_dtlapserhs - ML_BSSN::ML_shiftrhs - ML_BSSN::ML_dtshiftrhs -" diff --git a/par/head_on.par b/par/head_on.par deleted file mode 100644 index e6c326d..0000000 --- a/par/head_on.par +++ /dev/null @@ -1,368 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - time - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - CarpetRegrid2 - CarpetInterp - LoopControl - - ReflectionSymmetry - - TwoPunctures - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN - ML_BSSN_Helper - SphericalSurface - TmunuBase - Dissipation - - carpetiohdf5 - - AHFinderDirect - CarpetTracker - - LocalInterp AEILocalInterp LocalReduce -" - -# output - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 64 -IOScalar::outScalar_vars = " - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 64 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - 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_dtlapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtshift - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" - -### Checkpointing - -CarpetIOHDF5::checkpoint = yes -IO::checkpoint_ID = no -IO::recover = "autoprobe" -IO::checkpoint_every = 512 -IO::out_proc_every = 1 -IO::checkpoint_keep = 1 -IO::checkpoint_dir = $parfile -IO::recover_dir = $parfile -Carpet::regrid_during_recovery = no -CarpetIOHDF5::use_grid_structure_from_checkpoint = yes - -#--- driver - -#Cactus::cctk_itlast = 0 -Cactus::terminate = time -Cactus::cctk_final_time = 80 -Carpet::use_buffer_zones = yes - -CarpetLib::interleave_communications = yes -CarpetLib::vary_tags = yes -CarpetLib::combine_sends = yes -CarpetLib::reduce_mpi_waitall = yes -CarpetLib::print_memstats_every = 1024 -Carpet::regrid_in_level_mode = no -Carpet::output_timers_every = 1024 -CarpetLib::print_timestats_every = 1024 -Carpet::print_timestats_every = 1 -#Carpet::init_each_timelevel = yes - - -carpet::verbose = no -carpet::veryverbose = no - -#--- Carpet prolongation order settings - -Carpet::prolongation_order_space = 5 -Carpet::prolongation_order_time = 2 - - -# MoL time integration - -MoL::ODE_Method = RK4 -Carpet::num_integrator_substeps = 4 -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 -time::dtfac = 0.25 - -#--- grid and symmetries - -driver::ghost_size = 3 -grid::type = "CoordBase" -CoordBase::domainsize = "minmax" -CoordBase::xmax = 65.536 -CoordBase::ymax = 65.536 -CoordBase::zmax = 65.536 -CoordBase::dx = 2.048 -CoordBase::dy = 2.048 -CoordBase::dz = 2.048 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes - -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 - - -carpet::max_refinement_levels = 7 -CarpetRegrid2::num_centres = 2 -CarpetRegrid2::num_levels_1 = 7 -CarpetRegrid2::position_x_1 = +3.0 -CarpetRegrid2::radius_1 [1] = 25.6 # 0.256 -CarpetRegrid2::radius_1 [2] = 12.8 # 0.128 -CarpetRegrid2::radius_1 [3] = 6.4 # 0.256 -CarpetRegrid2::radius_1 [4] = 3.2 # 0.128 -CarpetRegrid2::radius_1 [5] = 1.6 # 0.064 -CarpetRegrid2::radius_1 [6] = 0.8 # 0.032 - -CarpetRegrid2::num_levels_2 = 7 -CarpetRegrid2::position_x_2 = -3.0 -CarpetRegrid2::radius_2 [1] = 25.6 -CarpetRegrid2::radius_2 [2] = 12.8 -CarpetRegrid2::radius_2 [3] = 6.4 -CarpetRegrid2::radius_2 [4] = 3.2 -CarpetRegrid2::radius_2 [5] = 1.6 -CarpetRegrid2::radius_2 [6] = 0.8 - -### Refinement tracking - -CarpetTracker::surface [0] = 0 -CarpetTracker::surface [1] = 1 -CarpetRegrid2::regrid_every = 32 - - -#--- initial data - -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::par_b = 3.0 -TwoPunctures::TP_epsilon = 1e-6 -TwoPunctures::par_m_plus = 0.5 -TwoPunctures::par_m_minus = 0.5 -#TwoPunctures::grid_setup_method = evaluation -TwoPunctures::verbose = yes -TwoPunctures::do_residuum_debug_output = yes -TwoPunctures::do_initial_debug_output = yes -Carpet::init_3_timelevels = yes - - -#--- ML_BSSN - -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_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::LapseAdvectionCoeff = 1.0 -ML_BSSN::ShiftAdvectionCoeff = 1.0 - -ML_BSSN::ML_log_confac_bound = "radiative" -ML_BSSN::ML_metric_bound = "radiative" -ML_BSSN::ML_Gamma_bound = "radiative" -ML_BSSN::ML_trace_curv_bound = "radiative" -ML_BSSN::ML_curv_bound = "radiative" -ML_BSSN::ML_lapse_bound = "radiative" -ML_BSSN::ML_dtlapse_bound = "radiative" -ML_BSSN::ML_shift_bound = "radiative" -ML_BSSN::ML_dtshift_bound = "radiative" - -ADMMacros::spatial_order = 4 - -#--- Dissipation - -Dissipation::vars = " - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv -" -Dissipation::order = 5 - -### Horizons - -AHFinderDirect::N_horizons = 3 -AHFinderDirect::find_every = 6 -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::verbose_level = "physics details" -AHFinderDirect::move_origins = yes - -AHFinderDirect::origin_x [1] = 3.0 -AHFinderDirect::initial_guess__coord_sphere__x_center [1] = 3.0 -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::origin_x [2] =-3.0 -AHFinderDirect::initial_guess__coord_sphere__x_center [2] =-3.0 -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::origin_x [3] = 0 -AHFinderDirect::find_after_individual [3] = 1024 -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 - - -# Horizon surfaces - -SphericalSurface::nsurfaces = 3 -SphericalSurface::maxntheta = 73 -SphericalSurface::maxnphi = 76 - -SphericalSurface::ntheta [0] = 37 -SphericalSurface::nphi [0] = 76 -SphericalSurface::symmetric_z [0] = yes -SphericalSurface::nghoststheta [0] = 2 -SphericalSurface::nghostsphi [0] = 2 - -SphericalSurface::ntheta [1] = 37 -SphericalSurface::nphi [1] = 76 -SphericalSurface::symmetric_z [1] = yes -SphericalSurface::nghoststheta [1] = 2 -SphericalSurface::nghostsphi [1] = 2 - -SphericalSurface::ntheta [2] = 37 -SphericalSurface::nphi [2] = 76 -SphericalSurface::symmetric_z [2] = yes -SphericalSurface::nghoststheta [2] = 2 -SphericalSurface::nghostsphi [2] = 2 - - -#--- Analysis - -admconstraints::constraints_timelevels = 3 -admconstraints::constraints_persist = yes -admconstraints::bound = "static" - -nanchecker::check_every = 1 -nanchecker::check_vars = " -ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv -" -nanchecker::action_if_found = "terminate" diff --git a/par/ks-mclachlan.par b/par/ks-mclachlan.par new file mode 100644 index 0000000..7fc652f --- /dev/null +++ b/par/ks-mclachlan.par @@ -0,0 +1,400 @@ +Cactus::cctk_run_title = "Kerr-Schild" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 100.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "Fortran" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "Slab" + + + +ActiveThorns = "InitBase" + + + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" + +Carpet::verbose = no +Carpet::veryverbose = no +Carpet::schedule_barriers = no +Carpet::storage_verbose = no +Carpet::timers_verbose = no +CarpetLib::output_bboxes = no + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 10 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::init_fill_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 5120 +CarpetLib::print_timestats_every = 5120 +CarpetLib::print_memstats_every = 5120 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 512 +NaNChecker::action_if_found = "just warn" +NaNChecker::check_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 + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift +" + + + +ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry90 SymBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmin = 0.00 +CoordBase::ymin = 0.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 + +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 = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +CartGrid3D::type = "coordbase" + +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_z = no + + + +ActiveThorns = "SphericalSurface" + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + + + +ActiveThorns = "CarpetMask" + +CarpetMask::excluded_surface [0] = 0 +CarpetMask::excluded_surface_factor[0] = 1.0 + + + +ActiveThorns = "CarpetRegrid2 CarpetTracker" + +CarpetTracker::surface[0] = 0 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating90 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 6 +CarpetRegrid2::radius_1[1] = 60.0 +CarpetRegrid2::radius_1[2] = 30.0 +CarpetRegrid2::radius_1[3] = 15.0 +CarpetRegrid2::radius_1[4] = 7.5 +CarpetRegrid2::radius_1[5] = 3.75 + + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::dtfac = 0.40 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" + +ADMMacros::spatial_order = 4 + + + +ActiveThorns = "Exact" + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.8 +Exact::Kerr_KerrSchild__parabolic = yes +Exact::Kerr_KerrSchild__power = 8 +Exact::Kerr_KerrSchild__epsilon = 1.0 + + + +ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad" + +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::timelevels = 3 + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::AlphaDriver = 1.0 +ML_BSSN::BetaDriver = 1.0 +ML_BSSN::LapseAdvectionCoeff = 0.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_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" + + + +ActiveThorns = "Dissipation" + +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 +" + + + +ActiveThorns = "ADMConstraints" + +ADMConstraints::constraints_persist = yes +ADMConstraints::constraints_timelevels = 3 + + + +ActiveThorns = "AHFinderDirect" + +AHFinderDirect::find_every = 1280 + +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +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::output_h_every = 0 + +AHFinderDirect::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no + + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 1280 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + ADMConstraints::ham +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 1280 +IOScalar::outScalar_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum + SphericalSurface::sf_radius +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 1280 +IOASCII::out0D_vars = " + Carpet::timing + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum + SphericalSurface::sf_active + SphericalSurface::sf_valid + SphericalSurface::sf_info + SphericalSurface::sf_radius + SphericalSurface::sf_origin + SphericalSurface::sf_coordinate_descriptors +" + +IOASCII::out1D_every = 1280 +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum + SphericalSurface::sf_radius +" + +IOASCII::out2D_every = 1280 +IOASCII::out2D_vars = " + SphericalSurface::sf_radius +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 5120 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + +#IOHDF5::checkpoint = yes +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = yes +#IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_on_terminate = yes +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + + +ActiveThorns = "Formaline" + + + +ActiveThorns = "TimerReport" + +TimerReport::out_every = 5120 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20 diff --git a/par/ks_bssn_0_8th.par b/par/ks_bssn_0_8th.par deleted file mode 100644 index fa8b6c6..0000000 --- a/par/ks_bssn_0_8th.par +++ /dev/null @@ -1,172 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.4 -CoordBase::dy = 0.4 -CoordBase::dz = 0.4 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 4 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 4 -CoordBase::boundary_size_y_lower = 4 -CoordBase::boundary_size_z_lower = 4 -CoordBase::boundary_size_x_upper = 4 -CoordBase::boundary_size_y_upper = 4 -CoordBase::boundary_size_z_upper = 4 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -ML_BSSN::log_confac_bound = "none" -ML_BSSN::metric_bound = "none" -ML_BSSN::Gamma_bound = "none" -ML_BSSN::trace_curv_bound = "none" -ML_BSSN::curv_bound = "none" -ML_BSSN::lapse_bound = "none" -ML_BSSN::dtlapse_bound = "none" -ML_BSSN::shift_bound = "none" -ML_BSSN::dtshift_bound = "none" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" diff --git a/par/ks_bssn_1.par b/par/ks_bssn_1.par deleted file mode 100644 index 4f993f4..0000000 --- a/par/ks_bssn_1.par +++ /dev/null @@ -1,173 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN - TmunuBase -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.2 -CoordBase::dy = 0.2 -CoordBase::dz = 0.2 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -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" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" diff --git a/par/ks_bssn_1_8th.par b/par/ks_bssn_1_8th.par deleted file mode 100644 index 64b0497..0000000 --- a/par/ks_bssn_1_8th.par +++ /dev/null @@ -1,172 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.2 -CoordBase::dy = 0.2 -CoordBase::dz = 0.2 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 4 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 4 -CoordBase::boundary_size_y_lower = 4 -CoordBase::boundary_size_z_lower = 4 -CoordBase::boundary_size_x_upper = 4 -CoordBase::boundary_size_y_upper = 4 -CoordBase::boundary_size_z_upper = 4 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -ML_BSSN::log_confac_bound = "none" -ML_BSSN::metric_bound = "none" -ML_BSSN::Gamma_bound = "none" -ML_BSSN::trace_curv_bound = "none" -ML_BSSN::curv_bound = "none" -ML_BSSN::lapse_bound = "none" -ML_BSSN::dtlapse_bound = "none" -ML_BSSN::shift_bound = "none" -ML_BSSN::dtshift_bound = "none" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" diff --git a/par/ks_bssn_2.par b/par/ks_bssn_2.par deleted file mode 100644 index 747dae3..0000000 --- a/par/ks_bssn_2.par +++ /dev/null @@ -1,173 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN - TmunuBase -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.1 -CoordBase::dy = 0.1 -CoordBase::dz = 0.1 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -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" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" diff --git a/par/ks_bssn_2_8th.par b/par/ks_bssn_2_8th.par deleted file mode 100644 index 4c757d8..0000000 --- a/par/ks_bssn_2_8th.par +++ /dev/null @@ -1,172 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.1 -CoordBase::dy = 0.1 -CoordBase::dz = 0.1 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 4 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 4 -CoordBase::boundary_size_y_lower = 4 -CoordBase::boundary_size_z_lower = 4 -CoordBase::boundary_size_x_upper = 4 -CoordBase::boundary_size_y_upper = 4 -CoordBase::boundary_size_z_upper = 4 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -ML_BSSN::log_confac_bound = "none" -ML_BSSN::metric_bound = "none" -ML_BSSN::Gamma_bound = "none" -ML_BSSN::trace_curv_bound = "none" -ML_BSSN::curv_bound = "none" -ML_BSSN::lapse_bound = "none" -ML_BSSN::dtlapse_bound = "none" -ML_BSSN::shift_bound = "none" -ML_BSSN::dtshift_bound = "none" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::log_confac - ML_BSSN::metric - ML_BSSN::trace_curv - ML_BSSN::curv - ML_BSSN::Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::log_confacrhs - ML_BSSN::metricrhs - ML_BSSN::Gammarhs - ML_BSSN::trace_curvrhs - ML_BSSN::curvrhs -" diff --git a/par/ks_bssn_3.par b/par/ks_bssn_3.par deleted file mode 100644 index 43a96bf..0000000 --- a/par/ks_bssn_3.par +++ /dev/null @@ -1,172 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - Exact - Formaline - GSL - - Slab - - GenericFD - - ML_BSSN -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.05 -CoordBase::dy = 0.05 -CoordBase::dz = 0.05 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -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" - -ADMMacros::spatial_order = 4 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" diff --git a/par/ks_bssn_evol.par b/par/ks_bssn_evol.par deleted file mode 100644 index 0c8852c..0000000 --- a/par/ks_bssn_evol.par +++ /dev/null @@ -1,327 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - Time - - ADMBase - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetInterp - CarpetLib - CarpetReduce - CarpetRegrid - LoopControl - - IsolatedHorizon - ReflectionSymmetry - - AEILocalInterp - AHFinderDirect - Dissipation - Exact - Formaline - GSL - SphericalSurface - - Slab - - GenericFD - - ML_BSSN - ML_BSSN_Helper - TmunuBase -" - -# output - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 32 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - ML_BSSN::ML_cons_detg - ML_BSSN::ML_cons_Gamma - ML_BSSN::ML_cons_traceA -" -# 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_dtlapse -# ML_BSSN::ML_shift -# ML_BSSN::ML_dtshift -# ML_BSSN::ML_log_confacrhs -# ML_BSSN::ML_metricrhs -# ML_BSSN::ML_Gammarhs -# ML_BSSN::ML_trace_curvrhs -# ML_BSSN::ML_curvrhs -# ML_BSSN::ML_curvrhs -# ML_BSSN::ML_lapserhs -# ML_BSSN::ML_dtlapserhs -# ML_BSSN::ML_shiftrhs -# ML_BSSN::ML_dtshiftrhs - -IOASCII::one_file_per_group = yes -IOASCII::out0D_every = 32 -IOASCII::out0D_vars = " - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_state - IsolatedHorizon::ih_grid_int - IsolatedHorizon::ih_grid_real - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_coordinates - IsolatedHorizon::ih_tetrad_l - IsolatedHorizon::ih_tetrad_n - IsolatedHorizon::ih_tetrad_m - IsolatedHorizon::ih_newman_penrose - IsolatedHorizon::ih_weyl_scalars - IsolatedHorizon::ih_twometric - IsolatedHorizon::ih_killing_vector - IsolatedHorizon::ih_killed_twometric - IsolatedHorizon::ih_scalars - IsolatedHorizon::ih_invariant_coordinates - IsolatedHorizon::ih_multipole_moments - IsolatedHorizon::ih_3determinant -" - -IOASCII::out1D_every = 32 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - ML_BSSN::ML_cons_detg - ML_BSSN::ML_cons_Gamma - ML_BSSN::ML_cons_traceA - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_newman_penrose - IsolatedHorizon::ih_weyl_scalars - IsolatedHorizon::ih_killing_vector - IsolatedHorizon::ih_killed_twometric - IsolatedHorizon::ih_3determinant -" -# 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_dtlapse -# ML_BSSN::ML_shift -# ML_BSSN::ML_dtshift -# ML_BSSN::ML_log_confacrhs -# ML_BSSN::ML_metricrhs -# ML_BSSN::ML_Gammarhs -# ML_BSSN::ML_trace_curvrhs -# ML_BSSN::ML_curvrhs -# ML_BSSN::ML_curvrhs -# ML_BSSN::ML_lapserhs -# ML_BSSN::ML_dtlapserhs -# ML_BSSN::ML_shiftrhs -# ML_BSSN::ML_dtshiftrhs - - -#--- driver - -Carpet::init_fill_timelevels = yes - -Cactus::cctk_itlast = 1024 -Carpet::use_buffer_zones = yes -Carpet::use_tapered_grids = yes - -CarpetLib::print_memstats_every = 1024 -CarpetLib::print_timestats_every = 1024 -Carpet::output_timers_every = 1024 -Carpet::print_timestats_every = 1024 - -carpet::verbose = yes -carpet::veryverbose = no - -#--- Carpet prolongation order settings - -Carpet::prolongation_order_space = 5 -Carpet::prolongation_order_time = 0 - - - -# MoL time integration - -MoL::ODE_Method = "RK4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -Time::dtfac = 0.25 - -#--- grid and symmetries - -driver::ghost_size = 3 -grid::type = "CoordBase" -CoordBase::domainsize = "minmax" -CoordBase::xmax = 96.0 -CoordBase::ymax = 96.0 -CoordBase::zmax = 96.0 -CoordBase::dx = 1.28 -CoordBase::dy = 1.28 -CoordBase::dz = 1.28 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes - -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 - - -carpetregrid::refined_regions = "manual-coordinate-list" - -carpetregrid::coordinates = " - [ [ ([0,0,0]:[32.0,32.0,32.0]:[0.64,0.64,0.64]) ], - [ ([0,0,0]:[16.0,16.0,16.0]:[0.32,0.32,0.32]) ], - [ ([0,0,0]:[8.0,8.0,8.0]:[0.16,0.16,0.16]) ], - [ ([0,0,0]:[4.0,4.0,4.0]:[0.08,0.08,0.08]) ], - [ ([0,0,0]:[2.0,2.0,2.0]:[0.04,0.04,0.04]) ], - [ ([0,0,0]:[1.0,1.0,1.0]:[0.02,0.02,0.02]) ]] -" - -carpetregrid::smart_outer_boundaries = "yes" -carpet::max_refinement_levels = 7 -carpetregrid::refinement_levels = 3 # 7 - -#--- initial data - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "exact" -ADMBase::initial_dtshift = "exact" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 -#Exact::Kerr_KerrSchild__epsilon = 1.0e-4 -Exact::Kerr_KerrSchild__parabolic = yes -Exact::Kerr_KerrSchild__power = 8 -Exact::Kerr_KerrSchild__epsilon = 1.0 - - -#--- ML_BSSN - -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" - -ML_BSSN::timelevels = 2 - -ML_BSSN::my_boundary_condition = "Minkowski" - -ML_BSSN::harmonicN = 1 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - -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" - -#--- Dissipation - -Dissipation::order = 5 -Dissipation::epsdis = 0.3 -Dissipation::vars = " - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv -" - -#--- Analysis - -nanchecker::check_every = 1 -nanchecker::check_vars = "admbase::gxx" -nanchecker::action_if_found = "terminate" - -#--- Horizons - -SphericalSurface::nsurfaces = 1 -SphericalSurface::maxntheta = 39 -SphericalSurface::maxnphi = 76 - -SphericalSurface::ntheta [0] = 39 -SphericalSurface::nphi [0] = 76 -SphericalSurface::nghoststheta[0] = 2 -SphericalSurface::nghostsphi [0] = 2 - -AHFinderDirect::find_every = 32 -AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no - -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::N_horizons = 1 - -AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 -AHFinderDirect::which_surface_to_store_info [1] = 0 -AHFinderDirect::reset_horizon_after_not_finding [1] = no - -IsolatedHorizon::verbose = yes -IsolatedHorizon::interpolator = "Lagrange polynomial interpolation" -IsolatedHorizon::interpolator_options = "order=4" -IsolatedHorizon::spatial_order = 4 - -IsolatedHorizon::num_horizons = 1 -IsolatedHorizon::surface_index[0] = 0 diff --git a/par/ks_ccatie_1.par b/par/ks_ccatie_1.par deleted file mode 100644 index 739fabb..0000000 --- a/par/ks_ccatie_1.par +++ /dev/null @@ -1,152 +0,0 @@ -ActiveThorns = " - Boundary - CartGrid3D - CoordBase - IOUtil - InitBase - MoL - SymBase - - ADMAnalysis - ADMBase - ADMConstraints - ADMCoupling - ADMMacros - CoordGauge - SpaceMask - StaticConformal - - NaNChecker - - Carpet - CarpetIOASCII - CarpetIOBasic - CarpetIOScalar - CarpetLib - CarpetReduce - LoopControl - - ReflectionSymmetry - - BSSN_MoL - Exact - Formaline - Fortran - GSL - - Slab -" - - - -Cactus::cctk_itlast = 0 - -grid::type = "CoordBase" - -CoordBase::domainsize = "minmax" - -CoordBase::xmax = 4.0 -CoordBase::ymax = 4.0 -CoordBase::zmax = 4.0 -CoordBase::dx = 0.2 -CoordBase::dy = 0.2 -CoordBase::dz = 0.2 - -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_x = no -ReflectionSymmetry::avoid_origin_y = no -ReflectionSymmetry::avoid_origin_z = no - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -Carpet::domain_from_coordbase = yes -driver::ghost_size = 2 -Carpet::use_buffer_zones = yes - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - - - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 - -ADM_BSSN::timelevels = 3 - -ADMBase::evolution_method = "ADM_BSSN" -ADM_BSSN::bound = "newrad" -ADM_BSSN::stencil_size = 2 - -ADMBase::lapse_evolution_method = "1+log" -ADM_BSSN::lapsesource = "modified" -ADM_BSSN::harmonic_f = 2.0 - -ADMBase::shift_evolution_method = "gamma2" -ADM_BSSN::ShiftGammaCoeff = 0.75 -ADM_BSSN::BetaDriver = 2 - -ADMMacros::spatial_order = 4 - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ADMConstraints::Hamiltonian - ADM_BSSN::ADM_BS_sK -" - -IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ADM_BSSN::ADM_BSSN_phi - ADM_BSSN::ADM_BSSN_metric - ADM_BSSN::ADM_BSSN_gamma - ADM_BSSN::ADM_BSSN_K - ADM_BSSN::ADM_BSSN_curv - ADM_BSSN::ADM_BSSN_dtlapse - ADM_BSSN::ADM_BSSN_B -" - -IOASCII::one_file_per_group = yes -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::lapse - ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ADM_BSSN::ADM_BSSN_phi - ADM_BSSN::ADM_BSSN_metric - ADM_BSSN::ADM_BSSN_gamma - ADM_BSSN::ADM_BSSN_K - ADM_BSSN::ADM_BSSN_curv - ADM_BSSN::ADM_BSSN_dtlapse - ADM_BSSN::ADM_BSSN_B -" diff --git a/par/m_bssn.par b/par/m_bssn.par deleted file mode 100644 index f32af64..0000000 --- a/par/m_bssn.par +++ /dev/null @@ -1,74 +0,0 @@ -ActiveThorns = " - Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase - ADMBase ADMCoupling ADMMacros StaticConformal - TmunuBase - NaNChecker - Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl - Formaline GSL - GenericFD - ML_BSSN ML_BSSN_Helper -" - - - -Cactus::cctk_itlast = 0 - -CoordBase::boundary_size_x_lower = 2 -CoordBase::boundary_size_y_lower = 2 -CoordBase::boundary_size_z_lower = 2 -CoordBase::boundary_size_x_upper = 2 -CoordBase::boundary_size_y_upper = 2 -CoordBase::boundary_size_z_upper = 2 - -driver::ghost_size = 2 - - - -ADMBase::initial_shift = "zero" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -ML_BSSN::my_initial_data = "Minkowski" - -ML_BSSN::harmonicN = 0 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 0.5 - - - -IO::out_dir = $parfile - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_vars = " - ML_BSSN::H - ML_BSSN::trKrhs -" - -IOScalar::outScalar_every = 1 -IOScalar::outScalar_vars = " - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" - -IOASCII::out1D_every = 1 -IOASCII::out1D_vars = " - ML_BSSN::Ham - ML_BSSN::mom - ML_BSSN::cons_detg - ML_BSSN::cons_Gamma - ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs -" diff --git a/par/old/bl-mclachlan-ab4.par b/par/old/bl-mclachlan-ab4.par new file mode 100644 index 0000000..56c7d42 --- /dev/null +++ b/par/old/bl-mclachlan-ab4.par @@ -0,0 +1,380 @@ +Cactus::cctk_run_title = "Brill-Lindquist" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 100.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "Fortran" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "Slab" + + + +ActiveThorns = "InitBase" + + + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" + +Carpet::verbose = no +Carpet::schedule_barriers = no +Carpet::veryverbose = no +Carpet::storage_verbose = no +Carpet::timers_verbose = no +CarpetLib::output_bboxes = no + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 10 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::init_fill_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 5120 +CarpetLib::print_timestats_every = 5120 +CarpetLib::print_memstats_every = 5120 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 1 # 512 +NaNChecker::verbose = "all" +NaNChecker::action_if_found = "just warn" +NaNChecker::check_vars = " + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift +" + + + +ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase" + +CoordBase::domainsize = "minmax" + +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 + +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_z_lower = 1 + +CartGrid3D::type = "coordbase" + +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_z = no + + + +ActiveThorns = "SphericalSurface" + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + + + +ActiveThorns = "CarpetMask" + +CarpetMask::excluded_surface [0] = 0 +CarpetMask::excluded_surface_factor[0] = 1.0 + + + +ActiveThorns = "CarpetRegrid2 CarpetTracker" + +CarpetTracker::surface[0] = 0 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 7 +CarpetRegrid2::radius_1[1] = 64.0 +CarpetRegrid2::radius_1[2] = 32.0 +CarpetRegrid2::radius_1[3] = 16.0 +CarpetRegrid2::radius_1[4] = 8.0 +CarpetRegrid2::radius_1[5] = 4.0 +CarpetRegrid2::radius_1[6] = 2.0 +CarpetRegrid2::radius_1[7] = 1.0 + + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "AB" +MoL::AB_Type = "4" +MoL::MoL_Intermediate_Steps = 1 + +Time::dtfac = 0.05 # 0.40 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" + +ADMMacros::spatial_order = 4 + + + +ActiveThorns = "Exact" + +ADMBase::metric_type = "physical" + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Schwarzschild/BL" +Exact::Schwarzschild_BL__mass = 1.0 + + + +ActiveThorns = "ML_BSSN ML_BSSN_Helper" + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_evolution_method = "ML_BSSN" + +ML_BSSN::timelevels = 3 +ML_BSSN::rhs_timelevels = 4 + +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::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.0 + +ML_BSSN::my_boundary_condition = "Minkowski" + +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" + + + +ActiveThorns = "Dissipation" + +Dissipation::order = 5 +Dissipation::vars = " + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv +" + + + +ActiveThorns = "ADMConstraints" + +ADMConstraints::constraints_persist = yes +ADMConstraints::constraints_timelevels = 3 + + + +ActiveThorns = "AHFinderDirect" + +AHFinderDirect::find_every = 1280 + +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +AHFinderDirect::move_origins = yes +AHFinderDirect::reshape_while_moving = yes +AHFinderDirect::predict_origin_movement = yes + +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::output_h_every = 0 + +AHFinderDirect::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no + + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 1280 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + ADMConstraints::ham +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 1280 +IOScalar::outScalar_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 1280 +IOASCII::out0D_vars = " + Carpet::timing + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + +IOASCII::out1D_every = 1280 +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + +IOASCII::out2D_every = 1280 +IOASCII::out2D_vars = " + SphericalSurface::sf_radius +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 5120 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::lapse + ADMConstraints::hamiltonian +" + +#IOHDF5::checkpoint = yes +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = yes +#IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_on_terminate = yes +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + + +ActiveThorns = "Formaline" + + + +ActiveThorns = "TimerReport" + +TimerReport::out_every = 5120 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20 diff --git a/par/old/bl_bssn_1.par b/par/old/bl_bssn_1.par new file mode 100644 index 0000000..862e994 --- /dev/null +++ b/par/old/bl_bssn_1.par @@ -0,0 +1,171 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.2 +CoordBase::dy = 0.2 +CoordBase::dz = 0.2 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Schwarzschild/BL" +Exact::Schwarzschild_BL__mass = 1.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +ML_BSSN::log_confac_bound = "none" +ML_BSSN::metric_bound = "none" +ML_BSSN::Gamma_bound = "none" +ML_BSSN::trace_curv_bound = "none" +ML_BSSN::curv_bound = "none" +ML_BSSN::lapse_bound = "none" +ML_BSSN::dtlapse_bound = "none" +ML_BSSN::shift_bound = "none" +ML_BSSN::dtshift_bound = "none" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" diff --git a/par/old/bl_bssn_2.par b/par/old/bl_bssn_2.par new file mode 100644 index 0000000..caefed4 --- /dev/null +++ b/par/old/bl_bssn_2.par @@ -0,0 +1,171 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.1 +CoordBase::dy = 0.1 +CoordBase::dz = 0.1 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Schwarzschild/BL" +Exact::Schwarzschild_BL__mass = 1.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +ML_BSSN::log_confac_bound = "none" +ML_BSSN::metric_bound = "none" +ML_BSSN::Gamma_bound = "none" +ML_BSSN::trace_curv_bound = "none" +ML_BSSN::curv_bound = "none" +ML_BSSN::lapse_bound = "none" +ML_BSSN::dtlapse_bound = "none" +ML_BSSN::shift_bound = "none" +ML_BSSN::dtshift_bound = "none" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" diff --git a/par/old/gw_bssn_1.par b/par/old/gw_bssn_1.par new file mode 100644 index 0000000..10d9255 --- /dev/null +++ b/par/old/gw_bssn_1.par @@ -0,0 +1,206 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + Time + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + Periodic + + Dissipation + Exact + Formaline + GSL + SphericalSurface + TmunuBase + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 100 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 1.0 +CoordBase::ymax = 0.0 +CoordBase::zmax = 0.0 +CoordBase::dx = 0.02 +CoordBase::dy = 0.02 +CoordBase::dz = 0.02 + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 1 +CoordBase::boundary_shiftout_z_upper = 1 + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + +Periodic::periodic = yes + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Minkowski/gauge wave" +Exact::Minkowski_gauge_wave__what_fn = "expsin" +Exact::Minkowski_gauge_wave__amplitude = 0.01 + +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" + +ADMMacros::spatial_order = 4 + + + +Time::dtfac = 0.5 + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Dissipation::epsdis = 0.25 +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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift +" + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + 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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs + ML_BSSN::ML_lapserhs + ML_BSSN::ML_dtlapserhs + ML_BSSN::ML_shiftrhs + ML_BSSN::ML_dtshiftrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + 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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs + ML_BSSN::ML_lapserhs + ML_BSSN::ML_dtlapserhs + ML_BSSN::ML_shiftrhs + ML_BSSN::ML_dtshiftrhs +" diff --git a/par/old/gw_bssn_2.par b/par/old/gw_bssn_2.par new file mode 100644 index 0000000..ff09c74 --- /dev/null +++ b/par/old/gw_bssn_2.par @@ -0,0 +1,206 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + Time + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + Periodic + + Dissipation + Exact + Formaline + GSL + SphericalSurface + TmunuBase + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 200 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 1.0 +CoordBase::ymax = 0.0 +CoordBase::zmax = 0.0 +CoordBase::dx = 0.01 +CoordBase::dy = 0.01 +CoordBase::dz = 0.01 + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 1 +CoordBase::boundary_shiftout_z_upper = 1 + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + +Periodic::periodic = yes + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Minkowski/gauge wave" +Exact::Minkowski_gauge_wave__what_fn = "expsin" +Exact::Minkowski_gauge_wave__amplitude = 0.01 + +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" + +ADMMacros::spatial_order = 4 + + + +Time::dtfac = 0.5 + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Dissipation::epsdis = 0.25 +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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift +" + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + 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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs + ML_BSSN::ML_lapserhs + ML_BSSN::ML_dtlapserhs + ML_BSSN::ML_shiftrhs + ML_BSSN::ML_dtshiftrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + 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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs + ML_BSSN::ML_lapserhs + ML_BSSN::ML_dtlapserhs + ML_BSSN::ML_shiftrhs + ML_BSSN::ML_dtshiftrhs +" diff --git a/par/old/head_on.par b/par/old/head_on.par new file mode 100644 index 0000000..e6c326d --- /dev/null +++ b/par/old/head_on.par @@ -0,0 +1,368 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + time + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + CarpetRegrid2 + CarpetInterp + LoopControl + + ReflectionSymmetry + + TwoPunctures + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN + ML_BSSN_Helper + SphericalSurface + TmunuBase + Dissipation + + carpetiohdf5 + + AHFinderDirect + CarpetTracker + + LocalInterp AEILocalInterp LocalReduce +" + +# output + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 64 +IOScalar::outScalar_vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 64 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + 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_dtlapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" + +### Checkpointing + +CarpetIOHDF5::checkpoint = yes +IO::checkpoint_ID = no +IO::recover = "autoprobe" +IO::checkpoint_every = 512 +IO::out_proc_every = 1 +IO::checkpoint_keep = 1 +IO::checkpoint_dir = $parfile +IO::recover_dir = $parfile +Carpet::regrid_during_recovery = no +CarpetIOHDF5::use_grid_structure_from_checkpoint = yes + +#--- driver + +#Cactus::cctk_itlast = 0 +Cactus::terminate = time +Cactus::cctk_final_time = 80 +Carpet::use_buffer_zones = yes + +CarpetLib::interleave_communications = yes +CarpetLib::vary_tags = yes +CarpetLib::combine_sends = yes +CarpetLib::reduce_mpi_waitall = yes +CarpetLib::print_memstats_every = 1024 +Carpet::regrid_in_level_mode = no +Carpet::output_timers_every = 1024 +CarpetLib::print_timestats_every = 1024 +Carpet::print_timestats_every = 1 +#Carpet::init_each_timelevel = yes + + +carpet::verbose = no +carpet::veryverbose = no + +#--- Carpet prolongation order settings + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + + +# MoL time integration + +MoL::ODE_Method = RK4 +Carpet::num_integrator_substeps = 4 +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 +time::dtfac = 0.25 + +#--- grid and symmetries + +driver::ghost_size = 3 +grid::type = "CoordBase" +CoordBase::domainsize = "minmax" +CoordBase::xmax = 65.536 +CoordBase::ymax = 65.536 +CoordBase::zmax = 65.536 +CoordBase::dx = 2.048 +CoordBase::dy = 2.048 +CoordBase::dz = 2.048 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes + +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 + + +carpet::max_refinement_levels = 7 +CarpetRegrid2::num_centres = 2 +CarpetRegrid2::num_levels_1 = 7 +CarpetRegrid2::position_x_1 = +3.0 +CarpetRegrid2::radius_1 [1] = 25.6 # 0.256 +CarpetRegrid2::radius_1 [2] = 12.8 # 0.128 +CarpetRegrid2::radius_1 [3] = 6.4 # 0.256 +CarpetRegrid2::radius_1 [4] = 3.2 # 0.128 +CarpetRegrid2::radius_1 [5] = 1.6 # 0.064 +CarpetRegrid2::radius_1 [6] = 0.8 # 0.032 + +CarpetRegrid2::num_levels_2 = 7 +CarpetRegrid2::position_x_2 = -3.0 +CarpetRegrid2::radius_2 [1] = 25.6 +CarpetRegrid2::radius_2 [2] = 12.8 +CarpetRegrid2::radius_2 [3] = 6.4 +CarpetRegrid2::radius_2 [4] = 3.2 +CarpetRegrid2::radius_2 [5] = 1.6 +CarpetRegrid2::radius_2 [6] = 0.8 + +### Refinement tracking + +CarpetTracker::surface [0] = 0 +CarpetTracker::surface [1] = 1 +CarpetRegrid2::regrid_every = 32 + + +#--- initial data + +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::par_b = 3.0 +TwoPunctures::TP_epsilon = 1e-6 +TwoPunctures::par_m_plus = 0.5 +TwoPunctures::par_m_minus = 0.5 +#TwoPunctures::grid_setup_method = evaluation +TwoPunctures::verbose = yes +TwoPunctures::do_residuum_debug_output = yes +TwoPunctures::do_initial_debug_output = yes +Carpet::init_3_timelevels = yes + + +#--- ML_BSSN + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_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::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.0 + +ML_BSSN::ML_log_confac_bound = "radiative" +ML_BSSN::ML_metric_bound = "radiative" +ML_BSSN::ML_Gamma_bound = "radiative" +ML_BSSN::ML_trace_curv_bound = "radiative" +ML_BSSN::ML_curv_bound = "radiative" +ML_BSSN::ML_lapse_bound = "radiative" +ML_BSSN::ML_dtlapse_bound = "radiative" +ML_BSSN::ML_shift_bound = "radiative" +ML_BSSN::ML_dtshift_bound = "radiative" + +ADMMacros::spatial_order = 4 + +#--- Dissipation + +Dissipation::vars = " + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv +" +Dissipation::order = 5 + +### Horizons + +AHFinderDirect::N_horizons = 3 +AHFinderDirect::find_every = 6 +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::verbose_level = "physics details" +AHFinderDirect::move_origins = yes + +AHFinderDirect::origin_x [1] = 3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center [1] = 3.0 +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::origin_x [2] =-3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center [2] =-3.0 +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::origin_x [3] = 0 +AHFinderDirect::find_after_individual [3] = 1024 +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 + + +# Horizon surfaces + +SphericalSurface::nsurfaces = 3 +SphericalSurface::maxntheta = 73 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 37 +SphericalSurface::nphi [0] = 76 +SphericalSurface::symmetric_z [0] = yes +SphericalSurface::nghoststheta [0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +SphericalSurface::ntheta [1] = 37 +SphericalSurface::nphi [1] = 76 +SphericalSurface::symmetric_z [1] = yes +SphericalSurface::nghoststheta [1] = 2 +SphericalSurface::nghostsphi [1] = 2 + +SphericalSurface::ntheta [2] = 37 +SphericalSurface::nphi [2] = 76 +SphericalSurface::symmetric_z [2] = yes +SphericalSurface::nghoststheta [2] = 2 +SphericalSurface::nghostsphi [2] = 2 + + +#--- Analysis + +admconstraints::constraints_timelevels = 3 +admconstraints::constraints_persist = yes +admconstraints::bound = "static" + +nanchecker::check_every = 1 +nanchecker::check_vars = " +ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv +" +nanchecker::action_if_found = "terminate" diff --git a/par/old/ks_bssn_0_8th.par b/par/old/ks_bssn_0_8th.par new file mode 100644 index 0000000..fa8b6c6 --- /dev/null +++ b/par/old/ks_bssn_0_8th.par @@ -0,0 +1,172 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.4 +CoordBase::dy = 0.4 +CoordBase::dz = 0.4 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 4 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 4 +CoordBase::boundary_size_y_lower = 4 +CoordBase::boundary_size_z_lower = 4 +CoordBase::boundary_size_x_upper = 4 +CoordBase::boundary_size_y_upper = 4 +CoordBase::boundary_size_z_upper = 4 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +ML_BSSN::log_confac_bound = "none" +ML_BSSN::metric_bound = "none" +ML_BSSN::Gamma_bound = "none" +ML_BSSN::trace_curv_bound = "none" +ML_BSSN::curv_bound = "none" +ML_BSSN::lapse_bound = "none" +ML_BSSN::dtlapse_bound = "none" +ML_BSSN::shift_bound = "none" +ML_BSSN::dtshift_bound = "none" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" diff --git a/par/old/ks_bssn_1.par b/par/old/ks_bssn_1.par new file mode 100644 index 0000000..4f993f4 --- /dev/null +++ b/par/old/ks_bssn_1.par @@ -0,0 +1,173 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN + TmunuBase +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.2 +CoordBase::dy = 0.2 +CoordBase::dz = 0.2 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +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" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" diff --git a/par/old/ks_bssn_1_8th.par b/par/old/ks_bssn_1_8th.par new file mode 100644 index 0000000..64b0497 --- /dev/null +++ b/par/old/ks_bssn_1_8th.par @@ -0,0 +1,172 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.2 +CoordBase::dy = 0.2 +CoordBase::dz = 0.2 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 4 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 4 +CoordBase::boundary_size_y_lower = 4 +CoordBase::boundary_size_z_lower = 4 +CoordBase::boundary_size_x_upper = 4 +CoordBase::boundary_size_y_upper = 4 +CoordBase::boundary_size_z_upper = 4 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +ML_BSSN::log_confac_bound = "none" +ML_BSSN::metric_bound = "none" +ML_BSSN::Gamma_bound = "none" +ML_BSSN::trace_curv_bound = "none" +ML_BSSN::curv_bound = "none" +ML_BSSN::lapse_bound = "none" +ML_BSSN::dtlapse_bound = "none" +ML_BSSN::shift_bound = "none" +ML_BSSN::dtshift_bound = "none" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" diff --git a/par/old/ks_bssn_2.par b/par/old/ks_bssn_2.par new file mode 100644 index 0000000..747dae3 --- /dev/null +++ b/par/old/ks_bssn_2.par @@ -0,0 +1,173 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN + TmunuBase +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.1 +CoordBase::dy = 0.1 +CoordBase::dz = 0.1 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +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" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" diff --git a/par/old/ks_bssn_2_8th.par b/par/old/ks_bssn_2_8th.par new file mode 100644 index 0000000..4c757d8 --- /dev/null +++ b/par/old/ks_bssn_2_8th.par @@ -0,0 +1,172 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.1 +CoordBase::dy = 0.1 +CoordBase::dz = 0.1 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 4 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 4 +CoordBase::boundary_size_y_lower = 4 +CoordBase::boundary_size_z_lower = 4 +CoordBase::boundary_size_x_upper = 4 +CoordBase::boundary_size_y_upper = 4 +CoordBase::boundary_size_z_upper = 4 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +ML_BSSN::log_confac_bound = "none" +ML_BSSN::metric_bound = "none" +ML_BSSN::Gamma_bound = "none" +ML_BSSN::trace_curv_bound = "none" +ML_BSSN::curv_bound = "none" +ML_BSSN::lapse_bound = "none" +ML_BSSN::dtlapse_bound = "none" +ML_BSSN::shift_bound = "none" +ML_BSSN::dtshift_bound = "none" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::log_confacrhs + ML_BSSN::metricrhs + ML_BSSN::Gammarhs + ML_BSSN::trace_curvrhs + ML_BSSN::curvrhs +" diff --git a/par/old/ks_bssn_3.par b/par/old/ks_bssn_3.par new file mode 100644 index 0000000..43a96bf --- /dev/null +++ b/par/old/ks_bssn_3.par @@ -0,0 +1,172 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + Exact + Formaline + GSL + + Slab + + GenericFD + + ML_BSSN +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.05 +CoordBase::dy = 0.05 +CoordBase::dz = 0.05 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +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" + +ADMMacros::spatial_order = 4 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv + ML_BSSN::ML_Gamma + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" diff --git a/par/old/ks_bssn_evol.par b/par/old/ks_bssn_evol.par new file mode 100644 index 0000000..0c8852c --- /dev/null +++ b/par/old/ks_bssn_evol.par @@ -0,0 +1,327 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + Time + + ADMBase + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetInterp + CarpetLib + CarpetReduce + CarpetRegrid + LoopControl + + IsolatedHorizon + ReflectionSymmetry + + AEILocalInterp + AHFinderDirect + Dissipation + Exact + Formaline + GSL + SphericalSurface + + Slab + + GenericFD + + ML_BSSN + ML_BSSN_Helper + TmunuBase +" + +# output + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 32 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + ML_BSSN::ML_cons_detg + ML_BSSN::ML_cons_Gamma + ML_BSSN::ML_cons_traceA +" +# 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_dtlapse +# ML_BSSN::ML_shift +# ML_BSSN::ML_dtshift +# ML_BSSN::ML_log_confacrhs +# ML_BSSN::ML_metricrhs +# ML_BSSN::ML_Gammarhs +# ML_BSSN::ML_trace_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_lapserhs +# ML_BSSN::ML_dtlapserhs +# ML_BSSN::ML_shiftrhs +# ML_BSSN::ML_dtshiftrhs + +IOASCII::one_file_per_group = yes +IOASCII::out0D_every = 32 +IOASCII::out0D_vars = " + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_state + IsolatedHorizon::ih_grid_int + IsolatedHorizon::ih_grid_real + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_coordinates + IsolatedHorizon::ih_tetrad_l + IsolatedHorizon::ih_tetrad_n + IsolatedHorizon::ih_tetrad_m + IsolatedHorizon::ih_newman_penrose + IsolatedHorizon::ih_weyl_scalars + IsolatedHorizon::ih_twometric + IsolatedHorizon::ih_killing_vector + IsolatedHorizon::ih_killed_twometric + IsolatedHorizon::ih_scalars + IsolatedHorizon::ih_invariant_coordinates + IsolatedHorizon::ih_multipole_moments + IsolatedHorizon::ih_3determinant +" + +IOASCII::out1D_every = 32 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + ML_BSSN::ML_cons_detg + ML_BSSN::ML_cons_Gamma + ML_BSSN::ML_cons_traceA + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_newman_penrose + IsolatedHorizon::ih_weyl_scalars + IsolatedHorizon::ih_killing_vector + IsolatedHorizon::ih_killed_twometric + IsolatedHorizon::ih_3determinant +" +# 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_dtlapse +# ML_BSSN::ML_shift +# ML_BSSN::ML_dtshift +# ML_BSSN::ML_log_confacrhs +# ML_BSSN::ML_metricrhs +# ML_BSSN::ML_Gammarhs +# ML_BSSN::ML_trace_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_lapserhs +# ML_BSSN::ML_dtlapserhs +# ML_BSSN::ML_shiftrhs +# ML_BSSN::ML_dtshiftrhs + + +#--- driver + +Carpet::init_fill_timelevels = yes + +Cactus::cctk_itlast = 1024 +Carpet::use_buffer_zones = yes +Carpet::use_tapered_grids = yes + +CarpetLib::print_memstats_every = 1024 +CarpetLib::print_timestats_every = 1024 +Carpet::output_timers_every = 1024 +Carpet::print_timestats_every = 1024 + +carpet::verbose = yes +carpet::veryverbose = no + +#--- Carpet prolongation order settings + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 0 + + + +# MoL time integration + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::dtfac = 0.25 + +#--- grid and symmetries + +driver::ghost_size = 3 +grid::type = "CoordBase" +CoordBase::domainsize = "minmax" +CoordBase::xmax = 96.0 +CoordBase::ymax = 96.0 +CoordBase::zmax = 96.0 +CoordBase::dx = 1.28 +CoordBase::dy = 1.28 +CoordBase::dz = 1.28 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes + +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 + + +carpetregrid::refined_regions = "manual-coordinate-list" + +carpetregrid::coordinates = " + [ [ ([0,0,0]:[32.0,32.0,32.0]:[0.64,0.64,0.64]) ], + [ ([0,0,0]:[16.0,16.0,16.0]:[0.32,0.32,0.32]) ], + [ ([0,0,0]:[8.0,8.0,8.0]:[0.16,0.16,0.16]) ], + [ ([0,0,0]:[4.0,4.0,4.0]:[0.08,0.08,0.08]) ], + [ ([0,0,0]:[2.0,2.0,2.0]:[0.04,0.04,0.04]) ], + [ ([0,0,0]:[1.0,1.0,1.0]:[0.02,0.02,0.02]) ]] +" + +carpetregrid::smart_outer_boundaries = "yes" +carpet::max_refinement_levels = 7 +carpetregrid::refinement_levels = 3 # 7 + +#--- initial data + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 +#Exact::Kerr_KerrSchild__epsilon = 1.0e-4 +Exact::Kerr_KerrSchild__parabolic = yes +Exact::Kerr_KerrSchild__power = 8 +Exact::Kerr_KerrSchild__epsilon = 1.0 + + +#--- ML_BSSN + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::shift_evolution_method = "ML_BSSN" + +ML_BSSN::timelevels = 2 + +ML_BSSN::my_boundary_condition = "Minkowski" + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + +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" + +#--- Dissipation + +Dissipation::order = 5 +Dissipation::epsdis = 0.3 +Dissipation::vars = " + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_trace_curv + ML_BSSN::ML_curv +" + +#--- Analysis + +nanchecker::check_every = 1 +nanchecker::check_vars = "admbase::gxx" +nanchecker::action_if_found = "terminate" + +#--- Horizons + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +AHFinderDirect::find_every = 32 +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +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::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no + +IsolatedHorizon::verbose = yes +IsolatedHorizon::interpolator = "Lagrange polynomial interpolation" +IsolatedHorizon::interpolator_options = "order=4" +IsolatedHorizon::spatial_order = 4 + +IsolatedHorizon::num_horizons = 1 +IsolatedHorizon::surface_index[0] = 0 diff --git a/par/old/ks_bssn_evol_mp.par b/par/old/ks_bssn_evol_mp.par new file mode 100644 index 0000000..93cd92d --- /dev/null +++ b/par/old/ks_bssn_evol_mp.par @@ -0,0 +1,309 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + Time + + Coordinates + Interpolate + + ADMBase + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetInterp + CarpetLib + CarpetReduce + CarpetRegrid2 + LoopControl + + AEILocalInterp + AHFinderDirect + Dissipation + Exact + Formaline + GSL + SphericalSurface + + Slab + + GenericFD + + ML_BSSN_MP_O2 + ML_BSSN_MP_O2_Helper + TmunuBase +" + +# output + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN_MP_O2::H + ML_BSSN_MP_O2::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 32 +IOScalar::outScalar_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out0D_every = 32 +IOASCII::out0D_vars = " +" + +IOASCII::out1D_every = 32 +IOASCII::out1D_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + +IOASCII::out2D_every = 32 +IOASCII::out2D_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + + +#--- driver + +Carpet::init_fill_timelevels = yes + +Cactus::cctk_itlast = 0 # 1024 +Carpet::use_buffer_zones = yes +Carpet::use_tapered_grids = yes + +CarpetLib::print_memstats_every = 1024 +CarpetLib::print_timestats_every = 1024 +Carpet::output_timers_every = 1024 +Carpet::print_timestats_every = 1024 + +carpet::verbose = yes +carpet::veryverbose = no + +#--- Carpet prolongation order settings + +Carpet::prolongation_order_space = 3 # 5 +Carpet::prolongation_order_time = 1 + + + +# MoL time integration + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::timestep_method = "given" +Time::timestep = 0.05 + +#--- grid + +driver::ghost_size = 2 # 3 +grid::type = "MultiPatch" + +Coordinates::coordinate_system = "Thornburg04" +Coordinates::h_cartesian = 0.2 +Coordinates::h_radial = 0.2 +Coordinates::sphere_inner_radius = 2.0 +Coordinates::sphere_outer_radius = 5.0 +Coordinates::n_angular = 11 + +Coordinates::patch_boundary_size = 2 +Coordinates::additional_overlap_size = 2 +Coordinates::outer_boundary_size = 2 + +Interpolate::interpolator_options = "order=4" + +Carpet::domain_from_multipatch = yes + +Carpet::max_refinement_levels = 9 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::num_levels_1 = 1 # 3 +CarpetRegrid2::radius_1[1] = 0.1 +CarpetRegrid2::radius_1[2] = 0.05 + +#--- initial data + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 +#Exact::Kerr_KerrSchild__epsilon = 1.0e-4 +Exact::Kerr_KerrSchild__parabolic = yes +Exact::Kerr_KerrSchild__power = 8 +Exact::Kerr_KerrSchild__epsilon = 1.0 + + +#--- ML_BSSN + +ADMBase::evolution_method = "ML_BSSN_MP_O2" +ADMBase::lapse_evolution_method = "ML_BSSN_MP_O2" +ADMBase::shift_evolution_method = "ML_BSSN_MP_O2" + +ML_BSSN_MP_O2::timelevels = 2 + +ML_BSSN_MP_O2::my_boundary_condition = "Minkowski" + +ML_BSSN_MP_O2::harmonicN = 1 # 1+log +ML_BSSN_MP_O2::harmonicF = 2.0 # 1+log +ML_BSSN_MP_O2::ShiftGammaCoeff = 0.75 +ML_BSSN_MP_O2::BetaDriver = 0.5 + +ML_BSSN_MP_O2::ML_log_confac_bound = "none" +ML_BSSN_MP_O2::ML_metric_bound = "none" +ML_BSSN_MP_O2::ML_Gamma_bound = "none" +ML_BSSN_MP_O2::ML_trace_curv_bound = "none" +ML_BSSN_MP_O2::ML_curv_bound = "none" +ML_BSSN_MP_O2::ML_lapse_bound = "none" +ML_BSSN_MP_O2::ML_dtlapse_bound = "none" +ML_BSSN_MP_O2::ML_shift_bound = "none" +ML_BSSN_MP_O2::ML_dtshift_bound = "none" + +#--- Dissipation + +Dissipation::order = 3 # 5 +Dissipation::epsdis = 0.3 +Dissipation::vars = " + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv +" + +#--- Analysis + +nanchecker::check_every = 1 +nanchecker::check_vars = "admbase::gxx" +nanchecker::action_if_found = "terminate" + +#--- Horizons + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +AHFinderDirect::find_every = 32 +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +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::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no diff --git a/par/old/ks_bssn_evol_mp_hi.par b/par/old/ks_bssn_evol_mp_hi.par new file mode 100644 index 0000000..d4ed7cc --- /dev/null +++ b/par/old/ks_bssn_evol_mp_hi.par @@ -0,0 +1,309 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + Time + + Coordinates + Interpolate + + ADMBase + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetInterp + CarpetLib + CarpetReduce + CarpetRegrid2 + LoopControl + + AEILocalInterp + AHFinderDirect + Dissipation + Exact + Formaline + GSL + SphericalSurface + + Slab + + GenericFD + + ML_BSSN_MP_O2 + ML_BSSN_MP_O2_Helper + TmunuBase +" + +# output + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN_MP_O2::H + ML_BSSN_MP_O2::trKrhs +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 64 +IOScalar::outScalar_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + +IOASCII::one_file_per_group = yes +IOASCII::out0D_every = 64 +IOASCII::out0D_vars = " +" + +IOASCII::out1D_every = 64 +IOASCII::out1D_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + +IOASCII::out2D_every = 64 +IOASCII::out2D_vars = " + grid::coordinates + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_BSSN_MP_O2::ML_Ham + ML_BSSN_MP_O2::ML_mom + ML_BSSN_MP_O2::ML_cons_detg + ML_BSSN_MP_O2::ML_cons_Gamma + ML_BSSN_MP_O2::ML_cons_traceA + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_dtlapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_dtshift + ML_BSSN_MP_O2::ML_log_confacrhs + ML_BSSN_MP_O2::ML_metricrhs + ML_BSSN_MP_O2::ML_Gammarhs + ML_BSSN_MP_O2::ML_trace_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_curvrhs + ML_BSSN_MP_O2::ML_lapserhs + ML_BSSN_MP_O2::ML_dtlapserhs + ML_BSSN_MP_O2::ML_shiftrhs + ML_BSSN_MP_O2::ML_dtshiftrhs +" + + +#--- driver + +Carpet::init_fill_timelevels = yes + +Cactus::cctk_itlast = 0 # 1024 +Carpet::use_buffer_zones = yes +Carpet::use_tapered_grids = yes + +CarpetLib::print_memstats_every = 1024 +CarpetLib::print_timestats_every = 1024 +Carpet::output_timers_every = 1024 +Carpet::print_timestats_every = 1024 + +carpet::verbose = yes +carpet::veryverbose = no + +#--- Carpet prolongation order settings + +Carpet::prolongation_order_space = 3 # 5 +Carpet::prolongation_order_time = 1 + + + +# MoL time integration + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Time::timestep_method = "given" +Time::timestep = 0.05 + +#--- grid + +driver::ghost_size = 2 # 3 +grid::type = "MultiPatch" + +Coordinates::coordinate_system = "Thornburg04" +Coordinates::h_cartesian = 0.1 +Coordinates::h_radial = 0.1 +Coordinates::sphere_inner_radius = 2.0 +Coordinates::sphere_outer_radius = 5.0 +Coordinates::n_angular = 21 + +Coordinates::patch_boundary_size = 2 +Coordinates::additional_overlap_size = 2 +Coordinates::outer_boundary_size = 2 + +Interpolate::interpolator_options = "order=4" + +Carpet::domain_from_multipatch = yes + +Carpet::max_refinement_levels = 9 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::num_levels_1 = 1 # 3 +CarpetRegrid2::radius_1[1] = 0.1 +CarpetRegrid2::radius_1[2] = 0.05 + +#--- initial data + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 +#Exact::Kerr_KerrSchild__epsilon = 1.0e-4 +Exact::Kerr_KerrSchild__parabolic = yes +Exact::Kerr_KerrSchild__power = 8 +Exact::Kerr_KerrSchild__epsilon = 1.0 + + +#--- ML_BSSN + +ADMBase::evolution_method = "ML_BSSN_MP_O2" +ADMBase::lapse_evolution_method = "ML_BSSN_MP_O2" +ADMBase::shift_evolution_method = "ML_BSSN_MP_O2" + +ML_BSSN_MP_O2::timelevels = 2 + +ML_BSSN_MP_O2::my_boundary_condition = "Minkowski" + +ML_BSSN_MP_O2::harmonicN = 1 # 1+log +ML_BSSN_MP_O2::harmonicF = 2.0 # 1+log +ML_BSSN_MP_O2::ShiftGammaCoeff = 0.75 +ML_BSSN_MP_O2::BetaDriver = 0.5 + +ML_BSSN_MP_O2::ML_log_confac_bound = "none" +ML_BSSN_MP_O2::ML_metric_bound = "none" +ML_BSSN_MP_O2::ML_Gamma_bound = "none" +ML_BSSN_MP_O2::ML_trace_curv_bound = "none" +ML_BSSN_MP_O2::ML_curv_bound = "none" +ML_BSSN_MP_O2::ML_lapse_bound = "none" +ML_BSSN_MP_O2::ML_dtlapse_bound = "none" +ML_BSSN_MP_O2::ML_shift_bound = "none" +ML_BSSN_MP_O2::ML_dtshift_bound = "none" + +#--- Dissipation + +Dissipation::order = 3 # 5 +Dissipation::epsdis = 0.3 +Dissipation::vars = " + ML_BSSN_MP_O2::ML_Gamma + ML_BSSN_MP_O2::ML_lapse + ML_BSSN_MP_O2::ML_shift + ML_BSSN_MP_O2::ML_log_confac + ML_BSSN_MP_O2::ML_metric + ML_BSSN_MP_O2::ML_trace_curv + ML_BSSN_MP_O2::ML_curv +" + +#--- Analysis + +nanchecker::check_every = 1 +nanchecker::check_vars = "admbase::gxx" +nanchecker::action_if_found = "terminate" + +#--- Horizons + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +AHFinderDirect::find_every = 64 +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +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::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no diff --git a/par/old/ks_ccatie_1.par b/par/old/ks_ccatie_1.par new file mode 100644 index 0000000..739fabb --- /dev/null +++ b/par/old/ks_ccatie_1.par @@ -0,0 +1,152 @@ +ActiveThorns = " + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + + NaNChecker + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + + ReflectionSymmetry + + BSSN_MoL + Exact + Formaline + Fortran + GSL + + Slab +" + + + +Cactus::cctk_itlast = 0 + +grid::type = "CoordBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmax = 4.0 +CoordBase::ymax = 4.0 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.2 +CoordBase::dy = 0.2 +CoordBase::dz = 0.2 + +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no + +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +Carpet::domain_from_coordbase = yes +driver::ghost_size = 2 +Carpet::use_buffer_zones = yes + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + + + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +Exact::exact_model = "Kerr/Kerr-Schild" +Exact::Kerr_KerrSchild__mass = 1.0 +Exact::Kerr_KerrSchild__spin = 0.0 + +ADM_BSSN::timelevels = 3 + +ADMBase::evolution_method = "ADM_BSSN" +ADM_BSSN::bound = "newrad" +ADM_BSSN::stencil_size = 2 + +ADMBase::lapse_evolution_method = "1+log" +ADM_BSSN::lapsesource = "modified" +ADM_BSSN::harmonic_f = 2.0 + +ADMBase::shift_evolution_method = "gamma2" +ADM_BSSN::ShiftGammaCoeff = 0.75 +ADM_BSSN::BetaDriver = 2 + +ADMMacros::spatial_order = 4 + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ADMConstraints::Hamiltonian + ADM_BSSN::ADM_BS_sK +" + +IOScalar::one_file_per_group = yes +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ADM_BSSN::ADM_BSSN_phi + ADM_BSSN::ADM_BSSN_metric + ADM_BSSN::ADM_BSSN_gamma + ADM_BSSN::ADM_BSSN_K + ADM_BSSN::ADM_BSSN_curv + ADM_BSSN::ADM_BSSN_dtlapse + ADM_BSSN::ADM_BSSN_B +" + +IOASCII::one_file_per_group = yes +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ADM_BSSN::ADM_BSSN_phi + ADM_BSSN::ADM_BSSN_metric + ADM_BSSN::ADM_BSSN_gamma + ADM_BSSN::ADM_BSSN_K + ADM_BSSN::ADM_BSSN_curv + ADM_BSSN::ADM_BSSN_dtlapse + ADM_BSSN::ADM_BSSN_B +" diff --git a/par/old/m_bssn.par b/par/old/m_bssn.par new file mode 100644 index 0000000..f32af64 --- /dev/null +++ b/par/old/m_bssn.par @@ -0,0 +1,74 @@ +ActiveThorns = " + Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase + ADMBase ADMCoupling ADMMacros StaticConformal + TmunuBase + NaNChecker + Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl + Formaline GSL + GenericFD + ML_BSSN ML_BSSN_Helper +" + + + +Cactus::cctk_itlast = 0 + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + +driver::ghost_size = 2 + + + +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +ML_BSSN::my_initial_data = "Minkowski" + +ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 0.5 + + + +IO::out_dir = $parfile + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_vars = " + ML_BSSN::H + ML_BSSN::trKrhs +" + +IOScalar::outScalar_every = 1 +IOScalar::outScalar_vars = " + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" + +IOASCII::out1D_every = 1 +IOASCII::out1D_vars = " + ML_BSSN::Ham + ML_BSSN::mom + ML_BSSN::cons_detg + ML_BSSN::cons_Gamma + ML_BSSN::cons_traceA + ML_BSSN::ML_log_confacrhs + ML_BSSN::ML_metricrhs + ML_BSSN::ML_Gammarhs + ML_BSSN::ML_trace_curvrhs + ML_BSSN::ML_curvrhs +" diff --git a/par/qc0-mclachlan.par b/par/qc0-mclachlan.par index 5ccdaed..ca2d880 100644 --- a/par/qc0-mclachlan.par +++ b/par/qc0-mclachlan.par @@ -4,8 +4,6 @@ Cactus::cctk_full_warnings = yes Cactus::highlight_warning_messages = no Cactus::terminate = "time" -#Cactus::cctk_final_time = 1.0 -#Cactus::cctk_final_time = 10.0 Cactus::cctk_final_time = 100.0 @@ -16,32 +14,19 @@ IO::out_dir = $parfile -ActiveThorns = "ManualTermination" - -ManualTermination::max_walltime = @WALLTIME_HOURS@ # hours -ManualTermination::on_remaining_walltime = 60 # minutes - -ManualTermination::termination_from_file = yes -ManualTermination::create_termination_file = yes -ManualTermination::termination_file = "../TERMINATE" - - - ActiveThorns = "AEILocalInterp" ActiveThorns = "BLAS LAPACK" ActiveThorns = "Fortran" -ActiveThorns = "LocalInterp" - ActiveThorns = "GenericFD" ActiveThorns = "HDF5" -ActiveThorns = "LoopControl" +ActiveThorns = "LocalInterp" -#ActiveThorns = "NaNCatcher" +ActiveThorns = "LoopControl" ActiveThorns = "Slab" @@ -64,8 +49,10 @@ ActiveThorns = "InitBase" ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" Carpet::verbose = no -Carpet::schedule_barriers = no Carpet::veryverbose = no +Carpet::schedule_barriers = no +Carpet::storage_verbose = no +Carpet::timers_verbose = no CarpetLib::output_bboxes = no Carpet::domain_from_coordbase = yes @@ -79,29 +66,39 @@ Carpet::prolongation_order_time = 2 Carpet::convergence_level = 0 -Carpet::init_3_timelevels = yes +Carpet::init_fill_timelevels = yes Carpet::poison_new_timelevels = yes CarpetLib::poison_new_memory = yes -Carpet::output_timers_every = 512 -CarpetLib::print_timestats_every = 512 -CarpetLib::print_memstats_every = 512 +Carpet::output_timers_every = 5120 +CarpetLib::print_timestats_every = 5120 +CarpetLib::print_memstats_every = 5120 ActiveThorns = "NaNChecker" -NaNChecker::check_every = 512 -NaNChecker::action_if_found = "just warn" +NaNChecker::check_every = 1 # 512 +#NaNChecker::verbose = "all" +#NaNChecker::action_if_found = "just warn" +NaNChecker::action_if_found = "terminate" NaNChecker::check_vars = " - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv + 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 + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift " @@ -177,7 +174,7 @@ SphericalSurface::nghostsphi [5] = 2 ActiveThorns = "CarpetMask" -CarpetMask::verbose = yes +CarpetMask::verbose = no CarpetMask::excluded_surface [0] = 0 CarpetMask::excluded_surface_factor[0] = 1.0 @@ -196,9 +193,10 @@ CarpetTracker::surface[0] = 0 CarpetTracker::surface[1] = 1 CarpetTracker::surface[2] = 2 -CarpetRegrid2::regrid_every = 128 - -CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::regrid_every = 128 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::verbose = yes CarpetRegrid2::num_centres = 3 @@ -240,7 +238,9 @@ MoL::ODE_Method = "RK4" MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1 -Time::dtfac = 0.40 +Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128, 256]" + +Time::dtfac = 0.25 @@ -268,7 +268,6 @@ TwoPunctures::par_P_minus[1] = -0.3331917498 #TODO# TwoPunctures::grid_setup_method = "evaluation" -#TwoPunctures::TP_epsilon = 1.0e-4 TwoPunctures::TP_epsilon = 1.0e-2 TwoPunctures::TP_Tiny = 1.0e-2 @@ -276,11 +275,13 @@ TwoPunctures::verbose = yes -ActiveThorns = "ML_BSSN ML_BSSN_Helper" +ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad" -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" +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 @@ -289,17 +290,11 @@ ML_BSSN::BetaDriver = 1.0 ML_BSSN::LapseAdvectionCoeff = 1.0 ML_BSSN::ShiftAdvectionCoeff = 1.0 -#ML_BSSN::ML_log_confac_bound = "radiative" -#ML_BSSN::ML_metric_bound = "radiative" -#ML_BSSN::ML_Gamma_bound = "radiative" -#ML_BSSN::ML_trace_curv_bound = "radiative" -#ML_BSSN::ML_curv_bound = "radiative" -#ML_BSSN::ML_lapse_bound = "radiative" -#ML_BSSN::ML_dtlapse_bound = "radiative" -#ML_BSSN::ML_shift_bound = "radiative" -#ML_BSSN::ML_dtshift_bound = "radiative" +ML_BSSN::MinimumLapse = 1.0e-8 -ML_BSSN::my_boundary_condition = "Minkowski" +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" +Boundary::radpower = 2 ML_BSSN::ML_log_confac_bound = "none" ML_BSSN::ML_metric_bound = "none" @@ -317,13 +312,15 @@ ActiveThorns = "Dissipation" Dissipation::order = 5 Dissipation::vars = " - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv + 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 " @@ -394,116 +391,13 @@ AHFinderDirect::find_after_individual_time [6] = 15.0 -ActiveThorns = "IsolatedHorizon" - -IsolatedHorizon::verbose = yes -IsolatedHorizon::interpolator = "Lagrange polynomial interpolation" -IsolatedHorizon::interpolator_options = "order=4" -IsolatedHorizon::spatial_order = 6 - -IsolatedHorizon::num_horizons = 6 - -IsolatedHorizon::surface_index [0] = 0 -IsolatedHorizon::companion_index[0] = 4 - -IsolatedHorizon::surface_index [1] = 1 - -IsolatedHorizon::surface_index [2] = 2 - -IsolatedHorizon::surface_index [3] = 3 -IsolatedHorizon::companion_index[3] = 5 - -IsolatedHorizon::surface_index [4] = 4 - -IsolatedHorizon::surface_index [5] = 5 - - - -ActiveThorns = "ProperDistance ProperTime" - -ProperDistance::number_geodesics = 512 -ProperDistance::direction = "x" -ProperDistance::plane = "xy" -ProperDistance::opening_angle = 180.0 -ProperDistance::step_size = 0.04 -ProperDistance::horizon_number = 1 -ProperDistance::interpolation_order = 3 -ProperDistance::integration_method = "rk4" -ProperDistance::eps = 0.00001 -ProperDistance::max_proper_distance = 15 -ProperDistance::calc_every = 128 -ProperDistance::use_second_horizon = yes -ProperDistance::second_horizon_number = 2 - - - -ActiveThorns = "PsiKadelia" - -PsiKadelia::psikadelia_persists = yes -PsiKadelia::calc_every = 256 # only on levels 0 and 1 -PsiKadelia::ricci_prolongation_type = "none" -PsiKadelia::weyl_timelevels = 3 -PsiKadelia::psif_vec = "standard-radial" - - - -ActiveThorns = "ModeDecomposition" - -ModeDecomposition::lmax = 8 -ModeDecomposition::ntheta = 35 -ModeDecomposition::nphi = 72 - -ModeDecomposition::num_vars = 5 -ModeDecomposition::variable_re[0] = "PsiKadelia::psi0re" -ModeDecomposition::variable_im[0] = "PsiKadelia::psi0im" -ModeDecomposition::spin_weight[0] = +2 -ModeDecomposition::variable_re[1] = "PsiKadelia::psi1re" -ModeDecomposition::variable_im[1] = "PsiKadelia::psi1im" -ModeDecomposition::spin_weight[1] = +1 -ModeDecomposition::variable_re[2] = "PsiKadelia::psi2re" -ModeDecomposition::variable_im[2] = "PsiKadelia::psi2im" -ModeDecomposition::spin_weight[2] = 0 -ModeDecomposition::variable_re[3] = "PsiKadelia::psi3re" -ModeDecomposition::variable_im[3] = "PsiKadelia::psi3im" -ModeDecomposition::spin_weight[3] = -1 -ModeDecomposition::variable_re[4] = "PsiKadelia::psi4re" -ModeDecomposition::variable_im[4] = "PsiKadelia::psi4im" -ModeDecomposition::spin_weight[4] = -2 - -ModeDecomposition::num_radii = 5 -ModeDecomposition::radius[0] = 30 -ModeDecomposition::radius[1] = 40 -ModeDecomposition::radius[2] = 50 -ModeDecomposition::radius[3] = 60 -ModeDecomposition::radius[4] = 70 - - - -ActiveThorns = "WaveExtract" - -WaveExtract::out_every = 128 -WaveExtract::maximum_detector_number = 5 -WaveExtract::switch_output_format = 100 -WaveExtract::rsch2_computation = "average Schwarzschild metric" -WaveExtract::l_mode = 8 -WaveExtract::m_mode = 8 -WaveExtract::detector_radius[0] = 30 -WaveExtract::detector_radius[1] = 40 -WaveExtract::detector_radius[2] = 50 -WaveExtract::detector_radius[3] = 60 -WaveExtract::detector_radius[4] = 70 - - - ActiveThorns = "CarpetIOBasic" IOBasic::outInfo_every = 128 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = " - Carpet::timing + Carpet::physical_time_per_hour ADMConstraints::ham - IsolatedHorizon::ih_spin[0] - IsolatedHorizon::ih_radius[0] " @@ -524,21 +418,6 @@ IOScalar::outScalar_vars = " ADMConstraints::hamiltonian ADMConstraints::momentum SphericalSurface::sf_radius - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_coordinates - IsolatedHorizon::ih_tetrad_l - IsolatedHorizon::ih_tetrad_n - IsolatedHorizon::ih_tetrad_m - IsolatedHorizon::ih_newman_penrose - IsolatedHorizon::ih_weyl_scalars - IsolatedHorizon::ih_twometric - IsolatedHorizon::ih_killing_vector - IsolatedHorizon::ih_killed_twometric - IsolatedHorizon::ih_invariant_coordinates - IsolatedHorizon::ih_multipole_moments - IsolatedHorizon::ih_3determinant - PsiKadelia::WeylComponents - ModeDecomposition::modes " @@ -568,28 +447,6 @@ IOASCII::out0D_vars = " SphericalSurface::sf_radius SphericalSurface::sf_origin SphericalSurface::sf_coordinate_descriptors - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_state - IsolatedHorizon::ih_grid_int - IsolatedHorizon::ih_grid_real - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_coordinates - IsolatedHorizon::ih_tetrad_l - IsolatedHorizon::ih_tetrad_n - IsolatedHorizon::ih_tetrad_m - IsolatedHorizon::ih_newman_penrose - IsolatedHorizon::ih_weyl_scalars - IsolatedHorizon::ih_twometric - IsolatedHorizon::ih_killing_vector - IsolatedHorizon::ih_killed_twometric - IsolatedHorizon::ih_scalars - IsolatedHorizon::ih_invariant_coordinates - IsolatedHorizon::ih_multipole_moments - IsolatedHorizon::ih_3determinant - PsiKadelia::WeylComponents - ModeDecomposition::modes - ProperDistance::pdistance - ProperTime::ptime " IOASCII::out1D_every = 128 @@ -604,29 +461,24 @@ IOASCII::out1D_vars = " ADMConstraints::hamiltonian ADMConstraints::momentum SphericalSurface::sf_radius - IsolatedHorizon::ih_shapes - IsolatedHorizon::ih_weyl_scalars - IsolatedHorizon::ih_killing_vector - IsolatedHorizon::ih_killed_twometric - IsolatedHorizon::ih_3determinant - PsiKadelia::WeylComponents - ModeDecomposition::modes " IOASCII::out2D_every = 128 IOASCII::out2D_vars = " SphericalSurface::sf_radius - PsiKadelia::WeylComponents " Activethorns = "CarpetIOHDF5" -IOHDF5::out_every = 512 -IOHDF5::one_file_per_group = yes -IOHDF5::compression_level = 1 -IOHDF5::out_vars = " +IOHDF5::out_every = 512 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " CarpetReduce::weight ADMBase::metric ADMBase::curv @@ -636,14 +488,13 @@ IOHDF5::out_vars = " ADMBase::dtshift ADMConstraints::hamiltonian ADMConstraints::momentum - PsiKadelia::WeylComponents " -IOHDF5::checkpoint = yes -IO::checkpoint_dir = $parfile -IO::checkpoint_ID = yes -IO::checkpoint_every = 512 -IO::checkpoint_on_terminate = yes +IOHDF5::checkpoint = yes +IO::checkpoint_dir = $parfile +IO::checkpoint_ID = yes +IO::checkpoint_every_walltime_hours = 6.0 +IO::checkpoint_on_terminate = yes IO::recover = "autoprobe" IO::recover_dir = $parfile @@ -652,14 +503,12 @@ IO::recover_dir = $parfile ActiveThorns = "Formaline" -#Formaline::publish_level = 2 -#Formaline::send_as_rdf = yes -#Formaline::rdf_hostname[0] = "buran.aei.mpg.de" -#Formaline::rdf_port [0] = 24997 - ActiveThorns = "TimerReport" -TimerReport::out_every = 512 -TimerReport::out_filename = "TimerReport" +TimerReport::out_every = 5120 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20 diff --git a/par/tov-mclachlan.par b/par/tov-mclachlan.par new file mode 100644 index 0000000..afb169b --- /dev/null +++ b/par/tov-mclachlan.par @@ -0,0 +1,404 @@ +Cactus::cctk_run_title = "Tolman-Oppenheimer-Volkoff" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 1500.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "Constants" + +ActiveThorns = "Fortran" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "Slab" + + + +ActiveThorns = "InitBase" + + + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" + +Carpet::verbose = no +Carpet::veryverbose = no +Carpet::schedule_barriers = no +Carpet::storage_verbose = no +Carpet::timers_verbose = no +CarpetLib::output_bboxes = no + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 10 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::init_fill_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 5120 +CarpetLib::print_timestats_every = 5120 +CarpetLib::print_memstats_every = 5120 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 1 # 512 +NaNChecker::verbose = "all" +NaNChecker::action_if_found = "just warn" +NaNChecker::check_vars = " + ML_BSSN_M::ML_log_confac + ML_BSSN_M::ML_metric + ML_BSSN_M::ML_trace_curv + ML_BSSN_M::ML_curv + ML_BSSN_M::ML_Gamma + ML_BSSN_M::ML_lapse + ML_BSSN_M::ML_shift + ML_BSSN_M::ML_dtlapse + ML_BSSN_M::ML_dtshift + Whisky::dens + Whisky::tau + Whisky::scon + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel +" + + + +ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry90 SymBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmin = 0.00 +CoordBase::ymin = 0.00 +CoordBase::zmin = 0.00 +CoordBase::xmax = +1800.00 +CoordBase::ymax = +1800.00 +CoordBase::zmax = +1800.00 +CoordBase::dx = 30.00 +CoordBase::dy = 30.00 +CoordBase::dz = 30.00 + +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 = 1 +CoordBase::boundary_shiftout_z_lower = 1 + +CartGrid3D::type = "coordbase" + +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_z = no + + + +ActiveThorns = "SphericalSurface" + + + +ActiveThorns = "CarpetRegrid2" + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating90 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 4 +CarpetRegrid2::radius_1[1] = 900.0 +CarpetRegrid2::radius_1[2] = 450.0 +CarpetRegrid2::radius_1[3] = 225.0 + + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "RK3" +MoL::MoL_Intermediate_Steps = 3 +MoL::MoL_Num_Scratch_Levels = 2 + +Time::dtfac = 0.25 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge HydroBase SpaceMask StaticConformal TmunuBase" + +ADMMacros::spatial_order = 4 + +HydroBase::timelevels = 3 +HydroBase::prolongation_type = "ENO" + +TmunuBase::stress_energy_storage = yes +TmunuBase::stress_energy_at_RHS = yes +TmunuBase::timelevels = 1 +TmunuBase::prolongation_type = "none" + +SpaceMask::use_mask = yes + + + +ActiveThorns = "EOS_Base EOS_Polytrope" + +EOS_Polytrope::eos_gamma = 1.3333333333 +EOS_Polytrope::eos_k = 0.4640517 + + + +ActiveThorns = "ML_BSSN_M ML_BSSN_M_Helper NewRad" + +ADMBase::evolution_method = "ML_BSSN_M" +ADMBase::lapse_evolution_method = "ML_BSSN_M" +ADMBase::shift_evolution_method = "ML_BSSN_M" +ADMBase::dtlapse_evolution_method = "ML_BSSN_M" +ADMBase::dtshift_evolution_method = "ML_BSSN_M" + +ML_BSSN_M::timelevels = 3 + +ML_BSSN_M::harmonicN = 1 # 1+log +ML_BSSN_M::harmonicF = 2.0 # 1+log +ML_BSSN_M::ShiftGammaCoeff = 0.75 +ML_BSSN_M::AlphaDriver = 0.0 +ML_BSSN_M::BetaDriver = 1.0 +ML_BSSN_M::LapseAdvectionCoeff = 1.0 +ML_BSSN_M::ShiftAdvectionCoeff = 1.0 + +ML_BSSN_M::MinimumLapse = 1.0e-8 + +ML_BSSN_M::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN_M::my_rhs_boundary_condition = "NewRad" +Boundary::radpower = 2 + +ML_BSSN_M::ML_log_confac_bound = "none" +ML_BSSN_M::ML_metric_bound = "none" +ML_BSSN_M::ML_Gamma_bound = "none" +ML_BSSN_M::ML_trace_curv_bound = "none" +ML_BSSN_M::ML_curv_bound = "none" +ML_BSSN_M::ML_lapse_bound = "none" +ML_BSSN_M::ML_dtlapse_bound = "none" +ML_BSSN_M::ML_shift_bound = "none" +ML_BSSN_M::ML_dtshift_bound = "none" + + + +ActiveThorns = "Whisky" + +HydroBase::evolution_method = "whisky" + +Whisky::riemann_solver = "Marquina" +Whisky::whisky_eos_type = "Polytype" +Whisky::whisky_eos_table = "2D_Polytrope" +Whisky::recon_method = "ppm" +Whisky::whisky_stencil = 3 +Whisky::bound = "none" +Whisky::rho_abs_min = 1.e-12 +Whisky::initial_atmosphere_factor = 0.999 + +Whisky::ppm_detect = no +Whisky::ppm_epsilon = 0.0 +Whisky::ppm_omega1 = 0.52 +Whisky::ppm_omega2 = 0.0 +Whisky::ppm_epsilon_shock = 0.1 +Whisky::ppm_eta1 = 5.0 +Whisky::ppm_k0 = 0.2 + + + +ActiveThorns = "Whisky_Init_Data Whisky_TOVSolverC" + +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +WhiskyTOVSolver::TOV_Rho_Central[0] = 1.61930347e-08 +WhiskyTOVSolver::TOV_Gamma[0] = 1.3333333333 +WhiskyTOVSolver::TOV_K[0] = 0.4640517 +WhiskyTOVSolver::TOV_Combine_Method = "maximum" +WhiskyTOVSolver::TOV_Num_Radial = 200000 +WhiskyTOVSolver::TOV_dr[0] = 0.05 + + + +ActiveThorns = "Dissipation" + +Dissipation::order = 5 +Dissipation::vars = " + ML_BSSN_M::ML_log_confac + ML_BSSN_M::ML_metric + ML_BSSN_M::ML_trace_curv + ML_BSSN_M::ML_curv + ML_BSSN_M::ML_Gamma + ML_BSSN_M::ML_lapse + ML_BSSN_M::ML_shift + ML_BSSN_M::ML_dtlapse + ML_BSSN_M::ML_dtshift +" + + + +ActiveThorns = "ADMConstraints" + +ADMConstraints::constraints_persist = yes +ADMConstraints::constraints_timelevels = 3 + + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 1280 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + ADMConstraints::ham + HydroBase::rho +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 1280 +IOScalar::outScalar_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 1280 +IOASCII::out0D_vars = " + Carpet::timing + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + +IOASCII::out1D_every = 1280 +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + HydroBase::rho + HydroBase::press + HydroBase::eps + HydroBase::vel + ADMConstraints::hamiltonian + ADMConstraints::momentum +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 5120 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::lapse + HydroBase::rho + ADMConstraints::hamiltonian +" + +#IOHDF5::checkpoint = yes +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = yes +#IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_on_terminate = yes +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + + +ActiveThorns = "Formaline" + + + +ActiveThorns = "TimerReport" + +TimerReport::out_every = 5120 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20 -- cgit v1.2.3