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 DYDCG_declare.h
@date Jun 98
@author Gabrielle Allen
@desc
Declarations for macro @seefile DYDCG_guts.h to compute first
derivatives of the conformal metric with respect to y
@enddesc
@@*/
#ifndef DYDCG_DECLARE
#ifdef FCODE
/* Input variables */
#undef DYDCG_GXX_JP
#define DYDCG_GXX_JP gxx(i,j+1,k)
#undef DYDCG_GXX_JM
#define DYDCG_GXX_JM gxx(i,j-1,k)
#undef DYDCG_GXY_JP
#define DYDCG_GXY_JP gxy(i,j+1,k)
#undef DYDCG_GXY_JM
#define DYDCG_GXY_JM gxy(i,j-1,k)
#undef DYDCG_GXZ_JP
#define DYDCG_GXZ_JP gxz(i,j+1,k)
#undef DYDCG_GXZ_JM
#define DYDCG_GXZ_JM gxz(i,j-1,k)
#undef DYDCG_GYY_JP
#define DYDCG_GYY_JP gyy(i,j+1,k)
#undef DYDCG_GYY_JM
#define DYDCG_GYY_JM gyy(i,j-1,k)
#undef DYDCG_GYZ_JP
#define DYDCG_GYZ_JP gyz(i,j+1,k)
#undef DYDCG_GYZ_JM
#define DYDCG_GYZ_JM gyz(i,j-1,k)
#undef DYDCG_GZZ_JP
#define DYDCG_GZZ_JP gzz(i,j+1,k)
#undef DYDCG_GZZ_JM
#define DYDCG_GZZ_JM gzz(i,j-1,k)
/* Output variables */
#undef DYDCG_DYDCGXX
#define DYDCG_DYDCGXX delgb211
#undef DYDCG_DYDCGXY
#define DYDCG_DYDCGXY delgb212
#undef DYDCG_DYDCGXZ
#define DYDCG_DYDCGXZ delgb213
#undef DYDCG_DYDCGYY
#define DYDCG_DYDCGYY delgb222
#undef DYDCG_DYDCGYZ
#define DYDCG_DYDCGYZ delgb223
#undef DYDCG_DYDCGZZ
#define DYDCG_DYDCGZZ delgb233
/* Internal variables */
#undef DYDCG_DY
#define DYDCG_DY dy
#undef DYDCG_OO2DY
#define DYDCG_OO2DY cdydg_oo2dy
/* Declare internal variables */
REAL DYDCG_OO2DY
/* Declare output variables */
REAL DYDCG_DYDCGXX
REAL DYDCG_DYDCGXY
REAL DYDCG_DYDCGXZ
REAL DYDCG_DYDCGYY
REAL DYDCG_DYDCGYZ
REAL DYDCG_DYDCGZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DYDCG_GXX_JP
#define DYDCG_GXX_JP gxx[ dj+ijk]
#undef DYDCG_GXX_JM
#define DYDCG_GXX_JM gxx[-dj+ijk]
#undef DYDCG_GXY_JP
#define DYDCG_GXY_JP gxy[ dj+ijk]
#undef DYDCG_GXY_JM
#define DYDCG_GXY_JM gxy[-dj+ijk]
#undef DYDCG_GXZ_JP
#define DYDCG_GXZ_JP gxz[ dj+ijk]
#undef DYDCG_GXZ_JM
#define DYDCG_GXZ_JM gxz[-dj+ijk]
#undef DYDCG_GYY_JP
#define DYDCG_GYY_JP gyy[ dj+ijk]
#undef DYDCG_GYY_JM
#define DYDCG_GYY_JM gyy[-dj+ijk]
#undef DYDCG_GYZ_JP
#define DYDCG_GYZ_JP gyz[ dj+ijk]
#undef DYDCG_GYZ_JM
#define DYDCG_GYZ_JM gyz[-dj+ijk]
#undef DYDCG_GZZ_JP
#define DYDCG_GZZ_JP gzz[ dj+ijk]
#undef DYDCG_GZZ_JM
#define DYDCG_GZZ_JM gzz[-dj+ijk]
/* Output variables */
#undef DYDCG_DYDCGXX
#define DYDCG_DYDCGXX delgb211
#undef DYDCG_DYDCGXY
#define DYDCG_DYDCGXY delgb212
#undef DYDCG_DYDCGXZ
#define DYDCG_DYDCGXZ delgb213
#undef DYDCG_DYDCGYY
#define DYDCG_DYDCGYY delgb222
#undef DYDCG_DYDCGYZ
#define DYDCG_DYDCGYZ delgb223
#undef DYDCG_DYDCGZZ
#define DYDCG_DYDCGZZ delgb233
/* Internal variables */
#undef DYDCG_OO2DY
#define DYDCG_OO2DY cdydg_oo2dy
/* Declare internal variables */
double DYDCG_OO2DY;
/* Declare output variables */
double DYDCG_DYDCGXX;
double DYDCG_DYDCGXY;
double DYDCG_DYDCGXZ;
double DYDCG_DYDCGYY;
double DYDCG_DYDCGYZ;
double DYDCG_DYDCGZZ;
#endif
#define DYDCG_DECLARE
#endif
|