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
|
/*@@
@header KK_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Declarations for macro to calculate term K_{ik}K^{k}_{l}
Requires: Upper physical metric
Lower extrinsic curvature
Provides: Components of K_{ik}K^{k}_{l}
@enddesc
@@*/
#ifndef KK_DECLARE
#include "UPPERMET_declare.h"
#ifdef FCODE
/* Input variables */
#undef KK_HXX
#define KK_HXX kxx(i,j,k)
#undef KK_HXY
#define KK_HXY kxy(i,j,k)
#undef KK_HXZ
#define KK_HXZ kxz(i,j,k)
#undef KK_HYY
#define KK_HYY kyy(i,j,k)
#undef KK_HYZ
#define KK_HYZ kyz(i,j,k)
#undef KK_HZZ
#define KK_HZZ kzz(i,j,k)
/* Output variables */
#undef KK_KKXX
#define KK_KKXX KK11
#undef KK_KKXY
#define KK_KKXY KK12
#undef KK_KKXZ
#define KK_KKXZ KK13
#undef KK_KKYY
#define KK_KKYY KK22
#undef KK_KKYZ
#define KK_KKYZ KK23
#undef KK_KKZZ
#define KK_KKZZ KK33
/* Declare output variables */
CCTK_REAL KK_KKXX
CCTK_REAL KK_KKXY
CCTK_REAL KK_KKXZ
CCTK_REAL KK_KKYY
CCTK_REAL KK_KKYZ
CCTK_REAL KK_KKZZ
#endif
#ifdef CCODE
/* Input variables */
#undef KK_HXX
#define KK_HXX kxx[ijk]
#undef KK_HXY
#define KK_HXY kxy[ijk]
#undef KK_HXZ
#define KK_HXZ kxz[ijk]
#undef KK_HYY
#define KK_HYY kyy[ijk]
#undef KK_HYZ
#define KK_HYZ kyz[ijk]
#undef KK_HZZ
#define KK_HZZ kzz[ijk]
/* Output variables */
#undef KK_KKXX
#define KK_KKXX KK11
#undef KK_KKXY
#define KK_KKXY KK12
#undef KK_KKXZ
#define KK_KKXZ KK13
#undef KK_KKYY
#define KK_KKYY KK22
#undef KK_KKYZ
#define KK_KKYZ KK23
#undef KK_KKZZ
#define KK_KKZZ KK33
/* Declare output variables */
CCTK_REAL KK_KKXX;
CCTK_REAL KK_KKXY;
CCTK_REAL KK_KKXZ;
CCTK_REAL KK_KKYY;
CCTK_REAL KK_KKYZ;
CCTK_REAL KK_KKZZ;
#endif
#define KK_DECLARE
#endif
|
