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 hxx(i,j+1,k)
#undef DYDK_KXX_JM
#define DYDK_KXX_JM hxx(i,j-1,k)
#undef DYDK_KXY_JP
#define DYDK_KXY_JP hxy(i,j+1,k)
#undef DYDK_KXY_JM
#define DYDK_KXY_JM hxy(i,j-1,k)
#undef DYDK_KXZ_JP
#define DYDK_KXZ_JP hxz(i,j+1,k)
#undef DYDK_KXZ_JM
#define DYDK_KXZ_JM hxz(i,j-1,k)
#undef DYDK_KYY_JP
#define DYDK_KYY_JP hyy(i,j+1,k)
#undef DYDK_KYY_JM
#define DYDK_KYY_JM hyy(i,j-1,k)
#undef DYDK_KYZ_JP
#define DYDK_KYZ_JP hyz(i,j+1,k)
#undef DYDK_KYZ_JM
#define DYDK_KYZ_JM hyz(i,j-1,k)
#undef DYDK_KZZ_JP
#define DYDK_KZZ_JP hzz(i,j+1,k)
#undef DYDK_KZZ_JM
#define DYDK_KZZ_JM hzz(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 */
REAL DYDK_OO2DY
/* Declare output variables */
REAL DYDK_DYDKXX
REAL DYDK_DYDKXY
REAL DYDK_DYDKXZ
REAL DYDK_DYDKYY
REAL DYDK_DYDKYZ
REAL DYDK_DYDKZZ
#endif
#ifdef CCODE
/* Input variables */
#undef DYDK_KXX_JP
#define DYDK_KXX_JP hxx[ dj+ijk]
#undef DYDK_KXX_JM
#define DYDK_KXX_JM hxx[-dj+ijk]
#undef DYDK_KXY_JP
#define DYDK_KXY_JP hxy[ dj+ijk]
#undef DYDK_KXY_JM
#define DYDK_KXY_JM hxy[-dj+ijk]
#undef DYDK_KXZ_JP
#define DYDK_KXZ_JP hxz[ dj+ijk]
#undef DYDK_KXZ_JM
#define DYDK_KXZ_JM hxz[-dj+ijk]
#undef DYDK_KYY_JP
#define DYDK_KYY_JP hyy[ dj+ijk]
#undef DYDK_KYY_JM
#define DYDK_KYY_JM hyy[-dj+ijk]
#undef DYDK_KYZ_JP
#define DYDK_KYZ_JP hyz[ dj+ijk]
#undef DYDK_KYZ_JM
#define DYDK_KYZ_JM hyz[-dj+ijk]
#undef DYDK_KZZ_JP
#define DYDK_KZZ_JP hzz[ dj+ijk]
#undef DYDK_KZZ_JM
#define DYDK_KZZ_JM hzz[-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
|