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
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
|
/*@@
@header DXZDG_declare.h
@date Jun 98
@author Gabrielle Allen
@desc
Declaration for macro to calculate the (first and) second derivatives
of the physical metric with respect to x,z
The macro uses @seefile DXDG_guts.h , @seefile DZDG_guts.h and
@seefile DXDG_declare.h , @seefile DZDG_declare.h
@enddesc
@@*/
#ifndef DXZDG_DECLARE
#define DXZDG_DECLARE
#include "DXDG_declare.h"
#include "DZDG_declare.h"
#ifdef FCODE
/* Input variables */
#undef DXZDG_GXX_IPKP
#define DXZDG_GXX_IPKP gxx(i+1,j,k+1)
#undef DXZDG_GXX_IPKM
#define DXZDG_GXX_IPKM gxx(i+1,j,k-1)
#undef DXZDG_GXX_IMKP
#define DXZDG_GXX_IMKP gxx(i-1,j,k+1)
#undef DXZDG_GXX_IMKM
#define DXZDG_GXX_IMKM gxx(i-1,j,k-1)
#undef DXZDG_GXY_IPKP
#define DXZDG_GXY_IPKP gxy(i+1,j,k+1)
#undef DXZDG_GXY_IPKM
#define DXZDG_GXY_IPKM gxy(i+1,j,k-1)
#undef DXZDG_GXY_IMKP
#define DXZDG_GXY_IMKP gxy(i-1,j,k+1)
#undef DXZDG_GXY_IMKM
#define DXZDG_GXY_IMKM gxy(i-1,j,k-1)
#undef DXZDG_GXZ_IPKP
#define DXZDG_GXZ_IPKP gxz(i+1,j,k+1)
#undef DXZDG_GXZ_IPKM
#define DXZDG_GXZ_IPKM gxz(i+1,j,k-1)
#undef DXZDG_GXZ_IMKP
#define DXZDG_GXZ_IMKP gxz(i-1,j,k+1)
#undef DXZDG_GXZ_IMKM
#define DXZDG_GXZ_IMKM gxz(i-1,j,k-1)
#undef DXZDG_GYY_IPKP
#define DXZDG_GYY_IPKP gyy(i+1,j,k+1)
#undef DXZDG_GYY_IPKM
#define DXZDG_GYY_IPKM gyy(i+1,j,k-1)
#undef DXZDG_GYY_IMKP
#define DXZDG_GYY_IMKP gyy(i-1,j,k+1)
#undef DXZDG_GYY_IMKM
#define DXZDG_GYY_IMKM gyy(i-1,j,k-1)
#undef DXZDG_GYZ_IPKP
#define DXZDG_GYZ_IPKP gyz(i+1,j,k+1)
#undef DXZDG_GYZ_IPKM
#define DXZDG_GYZ_IPKM gyz(i+1,j,k-1)
#undef DXZDG_GYZ_IMKP
#define DXZDG_GYZ_IMKP gyz(i-1,j,k+1)
#undef DXZDG_GYZ_IMKM
#define DXZDG_GYZ_IMKM gyz(i-1,j,k-1)
#undef DXZDG_GZZ_IPKP
#define DXZDG_GZZ_IPKP gzz(i+1,j,k+1)
#undef DXZDG_GZZ_IPKM
#define DXZDG_GZZ_IPKM gzz(i+1,j,k-1)
#undef DXZDG_GZZ_IMKP
#define DXZDG_GZZ_IMKP gzz(i-1,j,k+1)
#undef DXZDG_GZZ_IMKM
#define DXZDG_GZZ_IMKM gzz(i-1,j,k-1)
/* Output variables */
#undef DXZDG_DXZDGXX
#define DXZDG_DXZDGXX deldelg1311
#undef DXZDG_DXZDGXY
#define DXZDG_DXZDGXY deldelg1312
#undef DXZDG_DXZDGXZ
#define DXZDG_DXZDGXZ deldelg1313
#undef DXZDG_DXZDGYY
#define DXZDG_DXZDGYY deldelg1322
#undef DXZDG_DXZDGYZ
#define DXZDG_DXZDGYZ deldelg1323
#undef DXZDG_DXZDGZZ
#define DXZDG_DXZDGZZ deldelg1333
/* Internal variables */
#undef DXZDG_FAC
#define DXZDG_FAC dxzdg_fac
#undef DXZDG_OO4DXDZ
#define DXZDG_OO4DXDZ dxzdg_oo4dxdz
#undef DXZDG_DXZDPSI_O_PSI
#define DXZDG_DXZDPSI_O_PSI psixz(i,j,k)
/* Declare internal variables */
CCTK_REAL DXZDG_FAC
CCTK_REAL DXZDG_OO4DXDZ
/* Declare output variables */
CCTK_REAL DXZDG_DXZDGXX
CCTK_REAL DXZDG_DXZDGXY
CCTK_REAL DXZDG_DXZDGXZ
CCTK_REAL DXZDG_DXZDGYY
CCTK_REAL DXZDG_DXZDGYZ
CCTK_REAL DXZDG_DXZDGZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DXZDG_GXX_IPKP
#define DXZDG_GXX_IPKP gxx[ di + dk + ijk]
#undef DXZDG_GXX_IPKM
#define DXZDG_GXX_IPKM gxx[ di - dk + ijk]
#undef DXZDG_GXX_IMKP
#define DXZDG_GXX_IMKP gxx[-di + dk + ijk]
#undef DXZDG_GXX_IMKM
#define DXZDG_GXX_IMKM gxx[-di - dk + ijk]
#undef DXZDG_GXY_IPKP
#define DXZDG_GXY_IPKP gxy[ di + dk + ijk]
#undef DXZDG_GXY_IPKM
#define DXZDG_GXY_IPKM gxy[ di - dk + ijk]
#undef DXZDG_GXY_IMKP
#define DXZDG_GXY_IMKP gxy[-di + dk + ijk]
#undef DXZDG_GXY_IMKM
#define DXZDG_GXY_IMKM gxy[-di - dk + ijk]
#undef DXZDG_GXZ_IPKP
#define DXZDG_GXZ_IPKP gxz[ di + dk + ijk]
#undef DXZDG_GXZ_IPKM
#define DXZDG_GXZ_IPKM gxz[ di - dk + ijk]
#undef DXZDG_GXZ_IMKP
#define DXZDG_GXZ_IMKP gxz[-di + dk + ijk]
#undef DXZDG_GXZ_IMKM
#define DXZDG_GXZ_IMKM gxz[-di - dk + ijk]
#undef DXZDG_GYY_IPKP
#define DXZDG_GYY_IPKP gyy[ di + dk + ijk]
#undef DXZDG_GYY_IPKM
#define DXZDG_GYY_IPKM gyy[ di - dk + ijk]
#undef DXZDG_GYY_IMKP
#define DXZDG_GYY_IMKP gyy[-di + dk + ijk]
#undef DXZDG_GYY_IMKM
#define DXZDG_GYY_IMKM gyy[-di - dk + ijk]
#undef DXZDG_GYZ_IPKP
#define DXZDG_GYZ_IPKP gyz[ di + dk + ijk]
#undef DXZDG_GYZ_IPKM
#define DXZDG_GYZ_IPKM gyz[ di - dk + ijk]
#undef DXZDG_GYZ_IMKP
#define DXZDG_GYZ_IMKP gyz[-di + dk + ijk]
#undef DXZDG_GYZ_IMKM
#define DXZDG_GYZ_IMKM gyz[-di - dk + ijk]
#undef DXZDG_GZZ_IPKP
#define DXZDG_GZZ_IPKP gzz[ di + dk + ijk]
#undef DXZDG_GZZ_IPKM
#define DXZDG_GZZ_IPKM gzz[ di - dk + ijk]
#undef DXZDG_GZZ_IMKP
#define DXZDG_GZZ_IMKP gzz[-di + dk + ijk]
#undef DXZDG_GZZ_IMKM
#define DXZDG_GZZ_IMKM gzz[-di - dk + ijk]
/* Output variables */
#undef DXZDG_DXZDGXX
#define DXZDG_DXZDGXX deldelg1311
#undef DXZDG_DXZDGXY
#define DXZDG_DXZDGXY deldelg1312
#undef DXZDG_DXZDGXZ
#define DXZDG_DXZDGXZ deldelg1313
#undef DXZDG_DXZDGYY
#define DXZDG_DXZDGYY deldelg1322
#undef DXZDG_DXZDGYZ
#define DXZDG_DXZDGYZ deldelg1323
#undef DXZDG_DXZDGZZ
#define DXZDG_DXZDGZZ deldelg1333
/* Internal variables */
#undef DXZDG_FAC
#define DXZDG_FAC dxzdg_fac
#undef DXZDG_OO4DXDZ
#define DXZDG_OO4DXDZ dxzdg_oo4dxdz
#undef DXZDG_DXZDPSI_O_PSI
#define DXZDG_DXZDPSI_O_PSI psixz[ijk]
/* Declare internal variables */
CCTK_REAL DXZDG_FAC;
CCTK_REAL DXZDG_OO4DXDZ = 1/(4*cctkGH->cctk_delta_space[0]*cctkGH->cctk_delta_space[2]);
/* Declare output variables */
CCTK_REAL DXZDG_DXZDGXX;
CCTK_REAL DXZDG_DXZDGXY;
CCTK_REAL DXZDG_DXZDGXZ;
CCTK_REAL DXZDG_DXZDGYY;
CCTK_REAL DXZDG_DXZDGYZ;
CCTK_REAL DXZDG_DXZDGZZ;
#endif
#endif
|
