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
|
# Schedule definitions for thorn ADMAnalysis
# $Header$
SCHEDULE ADMAnalysis_ParamCheck AT CCTK_PARAMCHECK
{
LANG: C
OPTIONS: global
} "Check that the metric_type is recognised"
SCHEDULE ADMAnalysis_RegisterSymmetry AT CCTK_WRAGH
{
LANG: C
OPTIONS: global
} "Register symmetry of Ricci tensor and scalar"
SCHEDULE ADMAnalysis_EvaltrK AT CCTK_ANALYSIS
{
STORAGE: trace_of_K,detofg
LANG: C
TRIGGERS: trace_of_K,detofg
SYNC: trace_of_K,detofg
} "Compute the trace of the extrinsic curvature and the determinant of the metric"
SCHEDULE ADMAnalysis_MetricCartToSphere AT CCTK_ANALYSIS
{
STORAGE: spherical_metric
LANG: C
TRIGGERS: spherical_metric
SYNC: spherical_metric
} "Calculate the spherical metric in r,theta(q), phi(p)"
SCHEDULE ADMAnalysis_CurvCartToSphere AT CCTK_ANALYSIS
{
STORAGE: spherical_curv
LANG: C
TRIGGERS: spherical_curv
SYNC:spherical_curv
} "Calculate the spherical ex. curvature in r, theta(q), phi(p)"
if (ricci_persist)
{
if (ricci_timelevels == 1)
{
STORAGE: ricci_tensor[1], ricci_scalar[1]
}
else if (ricci_timelevels == 2)
{
STORAGE: ricci_tensor[2], ricci_scalar[2]
}
else if (ricci_timelevels == 3)
{
STORAGE: ricci_tensor[3], ricci_scalar[3]
}
SCHEDULE GROUP RicciGroup at CCTK_POSTINITIAL after (MoL_PostStep MoL_PostInitial)
{
STORAGE: detofg
} "Calculate Ricci tensor, with boundary conditions"
SCHEDULE GROUP RicciGroup at CCTK_POST_RECOVER_VARIABLES after (MoL_PostStep MoL_PostInitial)
{
STORAGE: detofg
} "Calculate Ricci tensor, with boundary conditions"
SCHEDULE GROUP RicciGroup at CCTK_EVOL after MoL_Evolution
{
STORAGE: detofg
} "Calculate Ricci tensor, with boundary conditions"
SCHEDULE GROUP RicciBoundariesGroup at CCTK_POSTREGRID
{
STORAGE: detofg
} "Set Ricci tensor on the boundary"
SCHEDULE GROUP RicciBoundariesGroup at CCTK_POSTRESTRICT
{
} "Set Ricci tensor on the boundary"
}
else
{
SCHEDULE GROUP RicciGroup at CCTK_ANALYSIS
{
STORAGE: ricci_tensor[1], ricci_scalar[1], detofg
TRIGGERS: ricci_tensor, ricci_scalar
} "Calculate Ricci tensor, with boundary conditions"
}
SCHEDULE ADMAnalysis_Ricci in RicciGroup
{
LANG: C
} "Calculate Ricci tensor, with boundary conditions"
schedule GROUP RicciBoundariesGroup in RicciGroup after ADMAnalysis_Ricci
{
} "Set Ricci tensor on the boundary"
SCHEDULE ADMAnalysis_Ricci_Boundaries in RicciBoundariesGroup
{
LANG: C
OPTIONS: level
SYNC: ricci_tensor, ricci_scalar
} "Select boundary conditions for the Ricci tensor"
SCHEDULE GROUP ApplyBCs as ADMAnalysis_ApplyBCs in RicciBoundariesGroup after ADMAnalysis_Ricci
{
} "Apply boundary conditions to the Ricci tensor"
|
