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
|
/*@@
@header LIEG_declare.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the Lie derivative of the lower
physical metric
@enddesc
@@*/
#ifndef LIEG_DECLARE
#define LIEG_DECLARE
#include "DB_declare.h"
#include "DG_declare.h"
#ifdef FCODE
/* Input variables */
#undef LIEG_BX
#define LIEG_BX betax(i,j,k)
#undef LIEG_BY
#define LIEG_BY betay(i,j,k)
#undef LIEG_BZ
#define LIEG_BZ betaz(i,j,k)
#undef LIEG_GXX
#define LIEG_GXX gxx(i,j,k)
#undef LIEG_GXY
#define LIEG_GXY gxy(i,j,k)
#undef LIEG_GXZ
#define LIEG_GXZ gxz(i,j,k)
#undef LIEG_GYY
#define LIEG_GYY gyy(i,j,k)
#undef LIEG_GYZ
#define LIEG_GYZ gyz(i,j,k)
#undef LIEG_GZZ
#define LIEG_GZZ gzz(i,j,k)
/* Output variables */
#undef LIEG_LGXX
#define LIEG_LGXX lieg_lgxx
#undef LIEG_LGXY
#define LIEG_LGXY lieg_lgxy
#undef LIEG_LGXZ
#define LIEG_LGXZ lieg_lgxz
#undef LIEG_LGYY
#define LIEG_LGYY lieg_lgyy
#undef LIEG_LGYZ
#define LIEG_LGYZ lieg_lgyz
#undef LIEG_LGZZ
#define LIEG_LGZZ lieg_lgzz
/* Declare output variables */
CCTK_REAL LIEG_LGXX
CCTK_REAL LIEG_LGXY
CCTK_REAL LIEG_LGXZ
CCTK_REAL LIEG_LGYY
CCTK_REAL LIEG_LGYZ
CCTK_REAL LIEG_LGZZ
#ifndef ADMMACROS_ADVECTION
integer admmacros_advectionx
integer admmacros_advectiony
integer admmacros_advectionz
#define ADMMACROS_ADVECTION
#endif
#endif
#ifdef CCODE
/* Input variables */
#undef LIEG_BX
#define LIEG_BX betax[ijk]
#undef LIEG_BY
#define LIEG_BY betay[ijk]
#undef LIEG_BZ
#define LIEG_BZ betaz[ijk]
#undef LIEG_GXX
#define LIEG_GXX gxx[ijk]
#undef LIEG_GXY
#define LIEG_GXY gxy[ijk]
#undef LIEG_GXZ
#define LIEG_GXZ gxz[ijk]
#undef LIEG_GYY
#define LIEG_GYY gyy[ijk]
#undef LIEG_GYZ
#define LIEG_GYZ gyz[ijk]
#undef LIEG_GZZ
#define LIEG_GZZ gzz[ijk]
/* Output variables */
#undef LIEG_LGXX
#define LIEG_LGXX lieg_lgxx
#undef LIEG_LGXY
#define LIEG_LGXY lieg_lgxy
#undef LIEG_LGXZ
#define LIEG_LGXZ lieg_lgxz
#undef LIEG_LGYY
#define LIEG_LGYY lieg_lgyy
#undef LIEG_LGYZ
#define LIEG_LGYZ lieg_lgyz
#undef LIEG_LGZZ
#define LIEG_LGZZ lieg_lgzz
/* Declare output variables */
CCTK_REAL LIEG_LGXX;
CCTK_REAL LIEG_LGXY;
CCTK_REAL LIEG_LGXZ;
CCTK_REAL LIEG_LGYY;
CCTK_REAL LIEG_LGYZ;
CCTK_REAL LIEG_LGZZ;
#endif
#endif
|
