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
|
/*@@
@header DYDK_declare.h
@date Jul 98
@author Gabrielle Allen
@desc
Declarations for macro @seefile DYDG_guts.h to compute first
derivatives of the extrinsic curvature with respect to y
@enddesc
@@*/
#ifndef DYDK_DECLARE
#ifdef FCODE
/* Input variables */
#undef DYDK_KXX_JP
#define DYDK_KXX_JP kxx(i,j+1,k)
#undef DYDK_KXX_JM
#define DYDK_KXX_JM kxx(i,j-1,k)
#undef DYDK_KXY_JP
#define DYDK_KXY_JP kxy(i,j+1,k)
#undef DYDK_KXY_JM
#define DYDK_KXY_JM kxy(i,j-1,k)
#undef DYDK_KXZ_JP
#define DYDK_KXZ_JP kxz(i,j+1,k)
#undef DYDK_KXZ_JM
#define DYDK_KXZ_JM kxz(i,j-1,k)
#undef DYDK_KYY_JP
#define DYDK_KYY_JP kyy(i,j+1,k)
#undef DYDK_KYY_JM
#define DYDK_KYY_JM kyy(i,j-1,k)
#undef DYDK_KYZ_JP
#define DYDK_KYZ_JP kyz(i,j+1,k)
#undef DYDK_KYZ_JM
#define DYDK_KYZ_JM kyz(i,j-1,k)
#undef DYDK_KZZ_JP
#define DYDK_KZZ_JP kzz(i,j+1,k)
#undef DYDK_KZZ_JM
#define DYDK_KZZ_JM kzz(i,j-1,k)
/* Output variables */
#undef DYDK_DYDKXX
#define DYDK_DYDKXX dydk_dydkxx
#undef DYDK_DYDKXY
#define DYDK_DYDKXY dydk_dydkxy
#undef DYDK_DYDKXZ
#define DYDK_DYDKXZ dydk_dydkxz
#undef DYDK_DYDKYY
#define DYDK_DYDKYY dydk_dydkyy
#undef DYDK_DYDKYZ
#define DYDK_DYDKYZ dydk_dydkyz
#undef DYDK_DYDKZZ
#define DYDK_DYDKZZ dydk_dydkzz
/* Internal variables */
#undef DYDK_DY
#define DYDK_DY dy
#undef DYDK_OO2DY
#define DYDK_OO2DY dydk_oo2dy
/* Declare internal variables */
CCTK_REAL DYDK_OO2DY
/* Declare output variables */
CCTK_REAL DYDK_DYDKXX
CCTK_REAL DYDK_DYDKXY
CCTK_REAL DYDK_DYDKXZ
CCTK_REAL DYDK_DYDKYY
CCTK_REAL DYDK_DYDKYZ
CCTK_REAL DYDK_DYDKZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DYDK_KXX_JP
#define DYDK_KXX_JP kxx[ dj+ijk]
#undef DYDK_KXX_JM
#define DYDK_KXX_JM kxx[-dj+ijk]
#undef DYDK_KXY_JP
#define DYDK_KXY_JP kxy[ dj+ijk]
#undef DYDK_KXY_JM
#define DYDK_KXY_JM kxy[-dj+ijk]
#undef DYDK_KXZ_JP
#define DYDK_KXZ_JP kxz[ dj+ijk]
#undef DYDK_KXZ_JM
#define DYDK_KXZ_JM kxz[-dj+ijk]
#undef DYDK_KYY_JP
#define DYDK_KYY_JP kyy[ dj+ijk]
#undef DYDK_KYY_JM
#define DYDK_KYY_JM kyy[-dj+ijk]
#undef DYDK_KYZ_JP
#define DYDK_KYZ_JP kyz[ dj+ijk]
#undef DYDK_KYZ_JM
#define DYDK_KYZ_JM kyz[-dj+ijk]
#undef DYDK_KZZ_JP
#define DYDK_KZZ_JP kzz[ dj+ijk]
#undef DYDK_KZZ_JM
#define DYDK_KZZ_JM kzz[-dj+ijk]
/* Output variables */
#undef DYDK_DYDKXX
#define DYDK_DYDKXX dydk_dydkxx
#undef DYDK_DYDKXY
#define DYDK_DYDKXY dydk_dydkxy
#undef DYDK_DYDKXZ
#define DYDK_DYDKXZ dydk_dydkxz
#undef DYDK_DYDKYY
#define DYDK_DYDKYY dydk_dydkyy
#undef DYDK_DYDKYZ
#define DYDK_DYDKYZ dydk_dydkyz
#undef DYDK_DYDKZZ
#define DYDK_DYDKZZ dydk_dydkzz
/* Internal variables */
#undef DYDK_OO2DY
#define DYDK_OO2DY dydk_oo2dy
/* Declare internal variables */
double DYDK_OO2DY;
/* Declare output variables */
double DYDK_DYDKXX;
double DYDK_DYDKXY;
double DYDK_DYDKXZ;
double DYDK_DYDKYY;
double DYDK_DYDKYZ;
double DYDK_DYDKZZ;
#endif
#define DYDK_DECLARE
#endif
|