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 hxx(i+1,j,k)
#undef DXDK_KXX_IM
#define DXDK_KXX_IM hxx(i-1,j,k)
#undef DXDK_KXY_IP
#define DXDK_KXY_IP hxy(i+1,j,k)
#undef DXDK_KXY_IM
#define DXDK_KXY_IM hxy(i-1,j,k)
#undef DXDK_KXZ_IP
#define DXDK_KXZ_IP hxz(i+1,j,k)
#undef DXDK_KXZ_IM
#define DXDK_KXZ_IM hxz(i-1,j,k)
#undef DXDK_KYY_IP
#define DXDK_KYY_IP hyy(i+1,j,k)
#undef DXDK_KYY_IM
#define DXDK_KYY_IM hyy(i-1,j,k)
#undef DXDK_KYZ_IP
#define DXDK_KYZ_IP hyz(i+1,j,k)
#undef DXDK_KYZ_IM
#define DXDK_KYZ_IM hyz(i-1,j,k)
#undef DXDK_KZZ_IP
#define DXDK_KZZ_IP hzz(i+1,j,k)
#undef DXDK_KZZ_IM
#define DXDK_KZZ_IM hzz(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 */
REAL DXDK_OO2DX
/* Declare output variables */
REAL DXDK_DXDKXX
REAL DXDK_DXDKXY
REAL DXDK_DXDKXZ
REAL DXDK_DXDKYY
REAL DXDK_DXDKYZ
REAL DXDK_DXDKZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DXDK_KXX_IP
#define DXDK_KXX_IP hxx[ di+ijk]
#undef DXDK_KXX_IM
#define DXDK_KXX_IM hxx[-di+ijk]
#undef DXDK_KXY_IP
#define DXDK_KXY_IP hxy[ di+ijk]
#undef DXDK_KXY_IM
#define DXDK_KXY_IM hxy[-di+ijk]
#undef DXDK_KXZ_IP
#define DXDK_KXZ_IP hxz[ di+ijk]
#undef DXDK_KXZ_IM
#define DXDK_KXZ_IM hxz[-di+ijk]
#undef DXDK_KYY_IP
#define DXDK_KYY_IP hyy[ di+ijk]
#undef DXDK_KYY_IM
#define DXDK_KYY_IM hyy[-di+ijk]
#undef DXDK_KYZ_IP
#define DXDK_KYZ_IP hyz[ di+ijk]
#undef DXDK_KYZ_IM
#define DXDK_KYZ_IM hyz[-di+ijk]
#undef DXDK_KZZ_IP
#define DXDK_KZZ_IP hzz[ di+ijk]
#undef DXDK_KZZ_IM
#define DXDK_KZZ_IM hzz[-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
|