blob: c822b1ece9f23a3fac9edf3a9ad4ebf606e8be48 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
ActiveThorns = "CoordBase SymBase Boundary NanChecker CartGrid3d Time MoL CarpetIOBasic CarpetIOScalar IOUtil Carpet CarpetLib CarpetReduce CarpetInterp CarpetSlab CarpetIOASCII ADMBase StaticConformal SpaceMask Slab Periodic Exact GenericFD CoordGauge ADMCoupling LoopControl ML_BSSN ML_BSSN_Helper GenericFD TMuNuBase SphericalSurface ADMMacros TimerReport"
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"
# We use the BBH gauge to ensure it is regression-tested. Also, there
# is no harmonic shift in McLachlan
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::MinimumLapse = 1.0e-8
ML_BSSN::conformalMethod = 1 # 1 for W, 0 for phi
#ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition = "static"
ML_BSSN::apply_dissipation = "always"
ML_BSSN::EpsDiss = 0.1
ML_BSSN::fdOrder = 8
Boundary::radpower = 2
#############################################################
# Grid
#############################################################
CoordBase::domainsize = minmax
CoordBase::boundary_size_x_lower = 5
CoordBase::boundary_size_y_lower = 5
CoordBase::boundary_size_z_lower = 5
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_y_lower = 1
CoordBase::boundary_shiftout_z_lower = 1
CoordBase::boundary_size_x_upper = 5
CoordBase::boundary_size_y_upper = 5
CoordBase::boundary_size_z_upper = 5
CoordBase::boundary_shiftout_x_upper = 0
CoordBase::boundary_shiftout_y_upper = 0
CoordBase::boundary_shiftout_z_upper = 0
CartGrid3D::type = "coordbase"
CartGrid3D::domain = "full"
CartGrid3D::avoid_origin = "no"
Periodic::periodic = "yes"
CoordBase::xmin = 0
CoordBase::ymin = 0
CoordBase::zmin = 0
CoordBase::xmax = 1
CoordBase::ymax = 1
CoordBase::zmax = 1
CoordBase::dx = 0.1
CoordBase::dy = 0.1
CoordBase::dz = 0.1
#############################################################
# Carpet
#############################################################
Carpet::ghost_size = 5
Carpet::domain_from_coordbase = "yes"
Carpet::init_3_timelevels = "no"
Carpet::poison_new_timelevels = yes
#############################################################
# Timers
#############################################################
# Do not commit these as testsuite output as they will change on
# different machines. They are useful for performance monitoring in
# automated build and test systems.
TimerReport::output_all_timers_readable = yes
TimerReport::n_top_timers = 40
TimerReport::output_schedule_timers = no
#############################################################
# Time integration
#############################################################
Cactus::terminate = "iteration"
Cactus::cctk_itlast = 8
# Use 0.5 for RK4 and 0.25 for ICN
Time::dtfac = 0.5
MethodOfLines::ode_method = "RK4"
MethodOfLines::MoL_Intermediate_Steps = 4
#MethodOfLines::ode_method = "generic"
#MethodOfLines::generic_type = "RK"
#MethodOfLines::MoL_Intermediate_Steps = 1
MethodOfLines::MoL_Num_Scratch_Levels = 1
MethodOfLines::MoL_NaN_Check = "no"
MethodOfLines::initial_data_is_crap = "no"
#############################################################
# 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 = "Minkowski/shifted gauge wave"
Exact::Minkowski_gauge_wave__what_fn = "sin"
Exact::Minkowski_gauge_wave__amplitude = 0.1
Exact::Minkowski_gauge_wave__omega = 1
Exact::Minkowski_gauge_wave__lambda = 0.577350269189626
Exact::rotation_euler_psi = -1.9216757376671543544
Exact::rotation_euler_theta = 0.66214523564555227398
Exact::rotation_euler_phi = 1.2199169159226388841
Exact::exact_order = 4
Exact::exact_eps = 1e-2 # This is large to ensure regression tests are insensitive to roundoff
#############################################################
# Output
#############################################################
IO::out_dir = $parfile
IO::out_fileinfo = "none"
IO::parfile_write = "no"
CarpetIOBasic::outInfo_every = 1
CarpetIOBasic::outInfo_vars = "ADMBase::alp"
IOASCII::out1D_every = 8
IOASCII::out_precision = 19
IOASCII::out1D_x = "yes"
IOASCII::out1D_y = "no"
IOASCII::out1D_z = "no"
IOASCII::out1D_d = "no"
IOASCII::out1D_vars = "ADMBase::gxx ADMBase::kxx ML_BSSN::ML_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA ML_BSSN::ML_curv ML_BSSN::ML_dtlapse ML_BSSN::ML_dtshift ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_log_confac ML_BSSN::ML_metric ML_BSSN::ML_shift ML_BSSN::ML_trace_curv"
|
