blob: 14e4369e40f7707462d2af294ed88fd495e32b15 (
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
|
Comparison between BSSN_MoL (ADM_BSSN) and McLachlan (ML_BSSN) parameters
2009-11-14: Written by Erik Schnetter, not yet trustworthy
XXXX-XX-XX: Cross-checked by NNN
Formulation:
ADMBase::evolution_method = "ADM_BSSN"
ADMBase::evolution_method = "ML_BSSN"
ADM_BSSN::RfromH = yes
[always for ML_BSSN]
ADM_BSSN::bound = "radiative"
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition = "NewRad"
# These are necessary
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"
Boundary::radpower
Boundary::radpower
Lapse:
ADMBase::lapse_evolution_method = "1+log"
ADMBase::lapse_evolution_method = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ML_BSSN::harmonicN = 1 # 1+log
ADM_BSSN::lapsesource = "straight"
ML_BSSN::AlphaDriver = 0
ADM_BSSN::lapsesource = "modified"
ML_BSSN::AlphaDriver = 1
ADM_BSSN::harmonic_f
ML_BSSN::harmonicF
ADM_BSSN::force_lapse_positive = yes
[always for ML_BSSN, but MinimumLapse can be negative]
ADM_BSSN::force_lapse_minimum
ML_BSSN::MinimumLapse
ADM_BSSN::lapse_advection_coeff
ML_BSSN::LapseAdvectionCoeff
[should be 0 or 1]
Shift:
ADMBase::shift_evolution_method = "gamma0"
ADMBase::shift_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"
ML_BSSN::ShiftGammaCoeff = 0
ADMBase::shift_evolution_method = "gamma2"
ADMBase::shift_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"
ML_BSSN::ShiftGammaCoeff = 1
ADM_BSSN::ShiftGammaCoeff
ML_BSSN::ShiftGammaCoeff
ADM_BSSN::ShiftAlpPower = 0
[always for ML_BSSN]
[There is a parameter ML_BSSN::ShiftAlphaPower, but it is unused]
ADM_BSSN::GammaDriver
[has no effect in ADM_BSSN]
ADM_BSSN::BetaDriver
ML_BSSN::BetaDriver
ADM_BSSN::gamma_driver_advection_coeff
ML_BSSN::ShiftAdvectionCoeff
[should be 0 or 1]
ADM_BSSN::ApplyShiftBoundary = no
[not possible with ML_BSSN; shift boundary is always applied]
Discretisation:
ADM_BSSN::timelevels
ML_BSSN::timelevels
ADM_BSSN::advection = "upwind4"
[always for ML_BSSN when 4th order is used]
|
