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
|
// cg.hh -- interface between machine-generated and other code
// $Id$
//
// This header file defines the "virtual machine" used by machine-generated
// code.
//
// prerequisites
// gfn.hh
//
// The machine-generated code uses the following variables:
// patch& p // current patch
// int irho,isigma // current generic integer coords within patch
//
//
// partial derivatives
//
#define PARTIAL_RHO(ghosted_gridfn_name) \
p.partial_rho(gfn__ ## ghosted_gridfn_name, irho,isigma)
#define PARTIAL_SIGMA(ghosted_gridfn_name) \
p.partial_sigma(gfn__ ## ghosted_gridfn_name, irho,isigma)
#define PARTIAL_RHO_RHO(ghosted_gridfn_name) \
p.partial_rho_rho(gfn__ ## ghosted_gridfn_name, irho,isigma)
#define PARTIAL_RHO_SIGMA(ghosted_gridfn_name) \
p.partial_rho_sigma(gfn__ ## ghosted_gridfn_name, irho,isigma)
#define PARTIAL_SIGMA_SIGMA(ghosted_gridfn_name) \
p.partial_sigma_sigma(gfn__ ## ghosted_gridfn_name, irho,isigma)
//
// ghosted-grid gridfns
//
#define h p.ghosted_gridfn(gfn__h, irho,isigma)
//
// nominal-grid gridfns
//
#define xx p.gridfn(gfn__xx, irho,isigma)
#define yy p.gridfn(gfn__yy, irho,isigma)
#define zz p.gridfn(gfn__zz, irho,isigma)
#define X_ud_11 p.gridfn(gfn__X_ud_11, irho,isigma)
#define X_ud_12 p.gridfn(gfn__X_ud_12, irho,isigma)
#define X_ud_13 p.gridfn(gfn__X_ud_13, irho,isigma)
#define X_ud_21 p.gridfn(gfn__X_ud_21, irho,isigma)
#define X_ud_22 p.gridfn(gfn__X_ud_22, irho,isigma)
#define X_ud_23 p.gridfn(gfn__X_ud_23, irho,isigma)
#define X_udd_111 p.gridfn(gfn__X_udd_111, irho,isigma)
#define X_udd_112 p.gridfn(gfn__X_udd_112, irho,isigma)
#define X_udd_113 p.gridfn(gfn__X_udd_113, irho,isigma)
#define X_udd_122 p.gridfn(gfn__X_udd_122, irho,isigma)
#define X_udd_123 p.gridfn(gfn__X_udd_123, irho,isigma)
#define X_udd_133 p.gridfn(gfn__X_udd_133, irho,isigma)
#define X_udd_211 p.gridfn(gfn__X_udd_211, irho,isigma)
#define X_udd_212 p.gridfn(gfn__X_udd_212, irho,isigma)
#define X_udd_213 p.gridfn(gfn__X_udd_213, irho,isigma)
#define X_udd_222 p.gridfn(gfn__X_udd_222, irho,isigma)
#define X_udd_223 p.gridfn(gfn__X_udd_223, irho,isigma)
#define X_udd_233 p.gridfn(gfn__X_udd_233, irho,isigma)
#define g_dd_11 p.gridfn(gfn__g_dd_11, irho,isigma)
#define g_dd_12 p.gridfn(gfn__g_dd_12, irho,isigma)
#define g_dd_13 p.gridfn(gfn__g_dd_13, irho,isigma)
#define g_dd_22 p.gridfn(gfn__g_dd_22, irho,isigma)
#define g_dd_23 p.gridfn(gfn__g_dd_23, irho,isigma)
#define g_dd_33 p.gridfn(gfn__g_dd_33, irho,isigma)
#define K_dd_11 p.gridfn(gfn__K_dd_11, irho,isigma)
#define K_dd_12 p.gridfn(gfn__K_dd_12, irho,isigma)
#define K_dd_13 p.gridfn(gfn__K_dd_13, irho,isigma)
#define K_dd_22 p.gridfn(gfn__K_dd_22, irho,isigma)
#define K_dd_23 p.gridfn(gfn__K_dd_23, irho,isigma)
#define K_dd_33 p.gridfn(gfn__K_dd_33, irho,isigma)
#define g_uu_11 p.gridfn(gfn__g_uu_11, irho,isigma)
#define g_uu_12 p.gridfn(gfn__g_uu_12, irho,isigma)
#define g_uu_13 p.gridfn(gfn__g_uu_13, irho,isigma)
#define g_uu_22 p.gridfn(gfn__g_uu_22, irho,isigma)
#define g_uu_23 p.gridfn(gfn__g_uu_23, irho,isigma)
#define g_uu_33 p.gridfn(gfn__g_uu_33, irho,isigma)
#define K_uu_11 p.gridfn(gfn__K_uu_11, irho,isigma)
#define K_uu_12 p.gridfn(gfn__K_uu_12, irho,isigma)
#define K_uu_13 p.gridfn(gfn__K_uu_13, irho,isigma)
#define K_uu_22 p.gridfn(gfn__K_uu_22, irho,isigma)
#define K_uu_23 p.gridfn(gfn__K_uu_23, irho,isigma)
#define K_uu_33 p.gridfn(gfn__K_uu_33, irho,isigma)
#define K p.gridfn(gfn__K, irho,isigma)
#define partial_d_g_dd_111 p.gridfn(gfn__partial_d_g_dd_111, irho,isigma)
#define partial_d_g_dd_112 p.gridfn(gfn__partial_d_g_dd_112, irho,isigma)
#define partial_d_g_dd_113 p.gridfn(gfn__partial_d_g_dd_113, irho,isigma)
#define partial_d_g_dd_122 p.gridfn(gfn__partial_d_g_dd_122, irho,isigma)
#define partial_d_g_dd_123 p.gridfn(gfn__partial_d_g_dd_123, irho,isigma)
#define partial_d_g_dd_133 p.gridfn(gfn__partial_d_g_dd_133, irho,isigma)
#define partial_d_g_dd_211 p.gridfn(gfn__partial_d_g_dd_211, irho,isigma)
#define partial_d_g_dd_212 p.gridfn(gfn__partial_d_g_dd_212, irho,isigma)
#define partial_d_g_dd_213 p.gridfn(gfn__partial_d_g_dd_213, irho,isigma)
#define partial_d_g_dd_222 p.gridfn(gfn__partial_d_g_dd_222, irho,isigma)
#define partial_d_g_dd_223 p.gridfn(gfn__partial_d_g_dd_223, irho,isigma)
#define partial_d_g_dd_233 p.gridfn(gfn__partial_d_g_dd_233, irho,isigma)
#define partial_d_g_dd_311 p.gridfn(gfn__partial_d_g_dd_311, irho,isigma)
#define partial_d_g_dd_312 p.gridfn(gfn__partial_d_g_dd_312, irho,isigma)
#define partial_d_g_dd_313 p.gridfn(gfn__partial_d_g_dd_313, irho,isigma)
#define partial_d_g_dd_322 p.gridfn(gfn__partial_d_g_dd_322, irho,isigma)
#define partial_d_g_dd_323 p.gridfn(gfn__partial_d_g_dd_323, irho,isigma)
#define partial_d_g_dd_333 p.gridfn(gfn__partial_d_g_dd_333, irho,isigma)
#define partial_d_ln_sqrt_g_1 p.gridfn(gfn__partial_d_ln_sqrt_g_1, irho,isigma)
#define partial_d_ln_sqrt_g_2 p.gridfn(gfn__partial_d_ln_sqrt_g_2, irho,isigma)
#define partial_d_ln_sqrt_g_3 p.gridfn(gfn__partial_d_ln_sqrt_g_3, irho,isigma)
#define partial_d_g_uu_111 p.gridfn(gfn__partial_d_g_uu_111, irho,isigma)
#define partial_d_g_uu_112 p.gridfn(gfn__partial_d_g_uu_112, irho,isigma)
#define partial_d_g_uu_113 p.gridfn(gfn__partial_d_g_uu_113, irho,isigma)
#define partial_d_g_uu_122 p.gridfn(gfn__partial_d_g_uu_122, irho,isigma)
#define partial_d_g_uu_123 p.gridfn(gfn__partial_d_g_uu_123, irho,isigma)
#define partial_d_g_uu_133 p.gridfn(gfn__partial_d_g_uu_133, irho,isigma)
#define partial_d_g_uu_211 p.gridfn(gfn__partial_d_g_uu_211, irho,isigma)
#define partial_d_g_uu_212 p.gridfn(gfn__partial_d_g_uu_212, irho,isigma)
#define partial_d_g_uu_213 p.gridfn(gfn__partial_d_g_uu_213, irho,isigma)
#define partial_d_g_uu_222 p.gridfn(gfn__partial_d_g_uu_222, irho,isigma)
#define partial_d_g_uu_223 p.gridfn(gfn__partial_d_g_uu_223, irho,isigma)
#define partial_d_g_uu_233 p.gridfn(gfn__partial_d_g_uu_233, irho,isigma)
#define partial_d_g_uu_311 p.gridfn(gfn__partial_d_g_uu_311, irho,isigma)
#define partial_d_g_uu_312 p.gridfn(gfn__partial_d_g_uu_312, irho,isigma)
#define partial_d_g_uu_313 p.gridfn(gfn__partial_d_g_uu_313, irho,isigma)
#define partial_d_g_uu_322 p.gridfn(gfn__partial_d_g_uu_322, irho,isigma)
#define partial_d_g_uu_323 p.gridfn(gfn__partial_d_g_uu_323, irho,isigma)
#define partial_d_g_uu_333 p.gridfn(gfn__partial_d_g_uu_333, irho,isigma)
#define H p.gridfn(gfn__H, irho,isigma)
|
