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
|
/*@@
@header DXDK_declare.h
@date Jul 98
@author Gabrielle Allen
@desc
Declarations for macro @seefile DXDK_guts.h to compute first
derivatives of the extrinsic curvature with respect to x
@enddesc
@@*/
#ifndef DXDK_DECLARE
#define DXDK_DECLARE
#ifdef FCODE
/* Input variables */
#undef DXDK_KXX_IP
#define DXDK_KXX_IP kxx(i+1,j,k)
#undef DXDK_KXX_IM
#define DXDK_KXX_IM kxx(i-1,j,k)
#undef DXDK_KXY_IP
#define DXDK_KXY_IP kxy(i+1,j,k)
#undef DXDK_KXY_IM
#define DXDK_KXY_IM kxy(i-1,j,k)
#undef DXDK_KXZ_IP
#define DXDK_KXZ_IP kxz(i+1,j,k)
#undef DXDK_KXZ_IM
#define DXDK_KXZ_IM kxz(i-1,j,k)
#undef DXDK_KYY_IP
#define DXDK_KYY_IP kyy(i+1,j,k)
#undef DXDK_KYY_IM
#define DXDK_KYY_IM kyy(i-1,j,k)
#undef DXDK_KYZ_IP
#define DXDK_KYZ_IP kyz(i+1,j,k)
#undef DXDK_KYZ_IM
#define DXDK_KYZ_IM kyz(i-1,j,k)
#undef DXDK_KZZ_IP
#define DXDK_KZZ_IP kzz(i+1,j,k)
#undef DXDK_KZZ_IM
#define DXDK_KZZ_IM kzz(i-1,j,k)
/* Output variables */
#undef DXDK_DXDKXX
#define DXDK_DXDKXX dxdk_dxdkxx
#undef DXDK_DXDKXY
#define DXDK_DXDKXY dxdk_dxdkxy
#undef DXDK_DXDKXZ
#define DXDK_DXDKXZ dxdk_dxdkxz
#undef DXDK_DXDKYY
#define DXDK_DXDKYY dxdk_dxdkyy
#undef DXDK_DXDKYZ
#define DXDK_DXDKYZ dxdk_dxdkyz
#undef DXDK_DXDKZZ
#define DXDK_DXDKZZ dxdk_dxdkzz
/* Internal variables */
#undef DXDK_DX
#define DXDK_DX dx
#undef DXDK_OO2DX
#define DXDK_OO2DX dxdk_oo2dx
/* Declare internal variables */
CCTK_REAL DXDK_OO2DX
/* Declare output variables */
CCTK_REAL DXDK_DXDKXX
CCTK_REAL DXDK_DXDKXY
CCTK_REAL DXDK_DXDKXZ
CCTK_REAL DXDK_DXDKYY
CCTK_REAL DXDK_DXDKYZ
CCTK_REAL DXDK_DXDKZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DXDK_KXX_IP
#define DXDK_KXX_IP kxx[ di+ijk]
#undef DXDK_KXX_IM
#define DXDK_KXX_IM kxx[-di+ijk]
#undef DXDK_KXY_IP
#define DXDK_KXY_IP kxy[ di+ijk]
#undef DXDK_KXY_IM
#define DXDK_KXY_IM kxy[-di+ijk]
#undef DXDK_KXZ_IP
#define DXDK_KXZ_IP kxz[ di+ijk]
#undef DXDK_KXZ_IM
#define DXDK_KXZ_IM kxz[-di+ijk]
#undef DXDK_KYY_IP
#define DXDK_KYY_IP kyy[ di+ijk]
#undef DXDK_KYY_IM
#define DXDK_KYY_IM kyy[-di+ijk]
#undef DXDK_KYZ_IP
#define DXDK_KYZ_IP kyz[ di+ijk]
#undef DXDK_KYZ_IM
#define DXDK_KYZ_IM kyz[-di+ijk]
#undef DXDK_KZZ_IP
#define DXDK_KZZ_IP kzz[ di+ijk]
#undef DXDK_KZZ_IM
#define DXDK_KZZ_IM kzz[-di+ijk]
/* Output variables */
#undef DXDK_DXDKXX
#define DXDK_DXDKXX dxdk_dxdkxx
#undef DXDK_DXDKXY
#define DXDK_DXDKXY dxdk_dxdkxy
#undef DXDK_DXDKXZ
#define DXDK_DXDKXZ dxdk_dxdkxz
#undef DXDK_DXDKYY
#define DXDK_DXDKYY dxdk_dxdkyy
#undef DXDK_DXDKYZ
#define DXDK_DXDKYZ dxdk_dxdkyz
#undef DXDK_DXDKZZ
#define DXDK_DXDKZZ dxdk_dxdkzz
/* Internal variables */
#undef DXDK_DX
#define DXDK_DX dx0
#undef DXDK_OO2DX
#define DXDK_OO2DX dxdk_oo2dx
/* Declare internal variables */
double DXDK_DX;
double DXDK_OO2DX;
/* Declare output variables */
double DXDK_DXDKXX;
double DXDK_DXDKXY;
double DXDK_DXDKXZ;
double DXDK_DXDKYY;
double DXDK_DXDKYZ;
double DXDK_DXDKZZ;
#endif
#endif
|