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
|
# Make sure there is storage for the conformal state flag
# so we can use it irrespective of metric_type
storage: conformal_state
schedule ADMConstraints_ParamCheck at CCTK_PARAMCHECK
{
LANG: C
OPTIONS: global
} "Check that we can deal with this metric_type and have enough conformal derivatives"
schedule ADMConstraints_ConformalCheck at CCTK_POSTINITIAL
{
LANG: C
} "Set conformal_state to 0 for physical, or check that it is at least 3"
schedule ADMConstraint_InitSymBound at CCTK_WRAGH
{
LANG: Fortran
OPTIONS: global
} "Register GF symmetries for ADM Constraints"
if (constraints_persist)
{
if (constraints_timelevels == 1)
{
STORAGE: hamiltonian[1], normalized_hamiltonian[1], momentum[1]
}
else if (constraints_timelevels == 2)
{
STORAGE: hamiltonian[2], normalized_hamiltonian[2], momentum[2]
}
else if (constraints_timelevels == 3)
{
STORAGE: hamiltonian[3], normalized_hamiltonian[3], momentum[3]
}
schedule GROUP ADMConstraintsGroup at CCTK_EVOL after MoL_Evolution
{
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
schedule GROUP ADMConstraintsGroup at CCTK_POSTREGRID
{
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
}
else
{
schedule GROUP ADMConstraintsGroup at CCTK_ANALYSIS
{
STORAGE: hamiltonian[1], normalized_hamiltonian[1], momentum[1]
TRIGGERS: hamiltonian, normalized_hamiltonian, momentum
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
}
schedule ADMConstraints in ADMConstraintsGroup
{
LANG: Fortran
} "Evaluate ADM constraints"
schedule ADMConstraints_Boundaries in ADMConstraintsGroup after ADMConstraints
{
LANG: Fortran
OPTIONS: level
SYNC: hamiltonian, normalized_hamiltonian, momentum
} "Select boundary conditions for the ADM constraints"
schedule GROUP ApplyBCs as ADMConstraints_ApplyBCs in ADMConstraintsGroup after ADMConstraints_Boundaries
{
} "Apply boundary conditions to the ADM constraints"
|
