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
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
|
/*
* inputs = {r, partial_d_ln_sqrt_g, partial_d_g_uu, X_ud, X_udd, g_uu, K_uu, h}
* outputs = {HA, HB, HC, HD}
* cost = 134*assignments+3*divisions+5*functions+401*multiplications+173*additions
*/
fp t1, t2, t4, t5, t7, t8, t10, t11, t12, t14;
fp t16, t18, t19, t20, t23, t24, t25, t26, t29, t30;
fp t31, t33, t35, t37, t38, t40, t41, t42, t43, t45;
fp t46, t47, t48, t50, t51, t52, t53, t55, t58, t59;
fp t61, t64, t67, t68, t70, t72, t73, t75, t77, t83;
fp t84, t87, t90, t91, t93, t95, t112, t130, t134, t137;
fp t139, t143, t147, t150, t152, t159, t164, t166, t168, t176;
fp t181, t183, t199, t203, t206, t209, t213, t216, t217, t219;
fp t228, t229, t232, t235, t236, t237, t240, t243, t247, t251;
fp t254, t259, t260, t262, t263, t269, t272, t273, t276, t279;
fp t282, t283, t290, t293, t296, t299, t302, t305, t312, t315;
fp t318, t324, t325, t329, t330, t338, t358, t381, t392, t395;
fp t396, t397, t408, t411, t431, t440, t444, t447, t450, t465;
t1 = g_uu_33;
t2 = 1/r;
t4 = X_ud_13;
t5 = PARTIAL_RHO(h);
t7 = X_ud_23;
t8 = PARTIAL_SIGMA(h);
t10 = zz*t2-t4*t5-t7*t8;
t11 = t1*t10;
t12 = partial_d_g_uu_311;
t14 = X_ud_11;
t16 = X_ud_21;
t18 = xx*t2-t14*t5-t16*t8;
t19 = t18*t18;
t20 = t12*t19;
t23 = g_uu_23;
t24 = t23*t10;
t25 = partial_d_g_uu_211;
t26 = t25*t19;
t29 = g_uu_13;
t30 = t29*t19;
t31 = partial_d_g_uu_312;
t33 = X_ud_12;
t35 = X_ud_22;
t37 = yy*t2-t33*t5-t35*t8;
t38 = t31*t37;
t40 = g_uu_12;
t41 = t40*t19;
t42 = partial_d_g_uu_212;
t43 = t42*t37;
t45 = g_uu_11;
t46 = t45*t19;
t47 = partial_d_g_uu_112;
t48 = t47*t37;
t50 = g_uu_22;
t51 = t37*t37;
t52 = t50*t51;
t53 = t42*t18;
t55 = t40*t51;
t58 = partial_d_g_uu_113;
t59 = t58*t10;
t61 = t23*t51;
t64 = t29*t10;
t67 = partial_d_g_uu_213;
t68 = t67*t10;
t70 = t45*t45;
t72 = yy*yy;
t73 = zz*zz;
t75 = r*r;
t77 = 1/t75/r;
t83 = t14*t14;
t84 = PARTIAL_RHO_RHO(h);
t87 = PARTIAL_RHO_SIGMA(h);
t90 = t16*t16;
t91 = PARTIAL_SIGMA_SIGMA(h);
t93 = (t72+t73)*t77-X_udd_111*t5-X_udd_211*t8-t83*t84-2.0*t16*t14*t87-t90
*t91;
t95 = RATIONAL(-1.0,2.0)*t11*t20+RATIONAL(-1.0,2.0)*t24*t26-t30*t38-t41*
t43-t46*t48-t52*t53-t55*t47*t18-t46*t59-t61*t31*t18-t64*t48*t18-t41*t68-t70*t19
*t93;
t112 = -xx*yy*t77-X_udd_112*t5-X_udd_212*t8-t14*t33*t84-t16*t33*t87-t14*
t35*t87-t16*t35*t91;
t130 = -xx*zz*t77-X_udd_113*t5-X_udd_213*t8-t14*t4*t84-t16*t4*t87-t14*t7*
t87-t16*t7*t91;
t134 = t40*t37;
t137 = t51*t37;
t139 = partial_d_g_uu_322;
t143 = partial_d_g_uu_122;
t147 = partial_d_g_uu_222;
t150 = t40*t40;
t152 = xx*xx;
t159 = t33*t33;
t164 = t35*t35;
t166 = (t152+t73)*t77-X_udd_122*t5-X_udd_222*t8-t159*t84-2.0*t35*t33*t87-
t164*t91;
t168 = t29*t29;
t176 = t4*t4;
t181 = t7*t7;
t183 = (t152+t72)*t77-X_udd_133*t5-X_udd_233*t8-t176*t84-2.0*t7*t4*t87-
t181*t91;
t199 = -yy*zz*t77-X_udd_123*t5-X_udd_223*t8-t33*t4*t84-t35*t4*t87-t33*t7*
t87-t35*t7*t91;
t203 = t40*t112;
t206 = t50*t37;
t209 = -t24*t43*t18-2.0*t41*t45*t112-2.0*t30*t45*t130-t134*t59*t18+
RATIONAL(-1.0,2.0)*t23*t137*t139+RATIONAL(-1.0,2.0)*t40*t137*t143+RATIONAL(-1.0
,2.0)*t50*t137*t147-t150*t19*t166-t168*t19*t183-2.0*t30*t40*t199-2.0*t52*t203-
t206*t68*t18;
t213 = t50*t50;
t216 = t10*t10;
t217 = t216*t10;
t219 = partial_d_g_uu_333;
t228 = partial_d_g_uu_123;
t229 = t228*t10;
t232 = partial_d_g_uu_233;
t235 = t23*t37;
t236 = partial_d_g_uu_313;
t237 = t236*t10;
t240 = t23*t23;
t243 = t29*t18;
t247 = t243*t183;
t251 = partial_d_g_uu_133;
t254 = -t150*t51*t93-t213*t51*t166+RATIONAL(-1.0,2.0)*t1*t217*t219-2.0*
t61*t40*t130-2.0*t61*t50*t199-t55*t229+RATIONAL(-1.0,2.0)*t23*t217*t232-t235*
t237*t18-t240*t51*t183-2.0*t134*t243*t130-2.0*t235*t247+RATIONAL(-1.0,2.0)*t29*
t217*t251;
t259 = partial_d_g_uu_323;
t260 = t259*t10;
t262 = partial_d_g_uu_223;
t263 = t262*t10;
t269 = t243*t199;
t272 = t40*t18;
t273 = t272*t166;
t276 = t272*t199;
t279 = t219*t216;
t282 = t45*t18;
t283 = t251*t216;
t290 = t232*t216;
t293 = t23*t216;
t296 = -2.0*t150*t37*t18*t112-t61*t260-t52*t263-2.0*t168*t10*t18*t130-2.0
*t206*t269-2.0*t206*t273-2.0*t235*t276+RATIONAL(-1.0,2.0)*t243*t279+RATIONAL(
-1.0,2.0)*t282*t283+RATIONAL(-1.0,2.0)*t235*t279+RATIONAL(-1.0,2.0)*t134*t283+
RATIONAL(-1.0,2.0)*t206*t290-t293*t262*t37;
t299 = t29*t216;
t302 = t147*t51;
t305 = t139*t51;
t312 = t282*t93;
t315 = t282*t130;
t318 = t282*t112;
t324 = t1*t216;
t325 = t236*t18;
t329 = -t299*t228*t37+RATIONAL(-1.0,2.0)*t272*t302+RATIONAL(-1.0,2.0)*t11
*t305+RATIONAL(-1.0,2.0)*t24*t302+RATIONAL(-1.0,2.0)*t272*t290-2.0*t134*t312
-2.0*t235*t315-2.0*t206*t318-t30*t237-2.0*t24*t269-t324*t325-2.0*t11*t247;
t330 = t259*t37;
t338 = t143*t51;
t358 = -t324*t330-2.0*t11*t315-t272*t263*t37-2.0*t11*t276+RATIONAL(-1.0,
2.0)*t282*t338-t243*t260*t37-2.0*t64*t312-t282*t229*t37-t299*t58*t18-2.0*t24*
t273-2.0*t24*t318+RATIONAL(-1.0,2.0)*t64*t338-2.0*t64*t235*t130;
t381 = t29*t130;
t392 = t23*t199;
t395 = -2.0*t240*t10*t37*t199-2.0*t11*t235*t183-2.0*t11*t134*t130-2.0*t64
*t134*t93-2.0*t24*t206*t166-t11*t38*t18-2.0*t293*t29*t112-2.0*t324*t381-2.0*t11
*t206*t199-2.0*t64*t206*t112-t168*t216*t93-2.0*t324*t392;
t396 = partial_d_g_uu_111;
t397 = t396*t19;
t408 = t1*t1;
t411 = t19*t18;
t431 = RATIONAL(-1.0,2.0)*t134*t397+RATIONAL(-1.0,2.0)*t206*t26+RATIONAL(
-1.0,2.0)*t235*t20+RATIONAL(-1.0,2.0)*t64*t397-t240*t216*t166-t408*t216*t183+
RATIONAL(-1.0,2.0)*t45*t411*t396+RATIONAL(-1.0,2.0)*t40*t411*t25+RATIONAL(-1.0,
2.0)*t29*t411*t12-2.0*t24*t134*t112+RATIONAL(-1.0,2.0)*t243*t305-t293*t67*t18
-2.0*t64*t272*t112;
HA = t95+t209+t254+t296+t329+t358+t395+t431;
t440 = t396*t18+t53+t325+t48+t147*t37+t330+t59+t263+t219*t10+t45*t93+2.0*
t203+2.0*t381;
t444 = partial_d_ln_sqrt_g_1;
t447 = partial_d_ln_sqrt_g_2;
t450 = partial_d_ln_sqrt_g_3;
t465 = t50*t166+2.0*t392+t1*t183+t45*t444*t18+t40*t447*t18+t29*t450*t18+
t40*t444*t37+t50*t447*t37+t23*t450*t37+t29*t444*t10+t23*t447*t10+t1*t450*t10;
HB = t440+t465;
HC = K_uu_11*t19+2.0*K_uu_12*t37*t18+2.0*K_uu_13*t10*t18+K_uu_22*t51+2.0*
K_uu_23*t10*t37+K_uu_33*t216;
HD = t46+2.0*t134*t18+2.0*t64*t18+t52+2.0*t24*t37+t324;
|
