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
|
# Schedule definitions for thorn ADMConstraints
# $Header$
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_POSTINITIAL after (MoL_PostStep MoL_PostInitial)
{
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
schedule GROUP ADMConstraintsGroup at CCTK_POST_RECOVER_VARIABLES after (MoL_PostStep MoL_PostInitial)
{
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
schedule GROUP ADMConstraintsGroup at CCTK_EVOL after MoL_Evolution
{
} "Evaluate ADM constraints, and perform symmetry boundary conditions"
schedule GROUP ADMConstraintsBoundariesGroup at CCTK_POSTREGRID after MoL_PostStep
{
} "Set ADM constraints on the boundary"
schedule GROUP ADMConstraintsBoundariesGroup at CCTK_POSTRESTRICT after MoL_PostStep
{
} "Set ADM constraints on the boundary"
}
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 GROUP ADMConstraintsBoundariesGroup in ADMConstraintsGroup after ADMConstraints
{
} "Set ADM constraints on the boundary"
schedule ADMConstraints_Boundaries in ADMConstraintsBoundariesGroup
{
LANG: Fortran
OPTIONS: level
SYNC: hamiltonian, normalized_hamiltonian, momentum
} "Select boundary conditions for the ADM constraints"
schedule GROUP ApplyBCs as ADMConstraints_ApplyBCs in ADMConstraintsBoundariesGroup after ADMConstraints_Boundaries
{
} "Apply boundary conditions to the ADM constraints"
|