te[ 0 ] = n23 * n34 * n42 - n24 * n33 * n42 + n24 * n32 * n43 - n22 * n34 * n43 - n23 * n32 * n44 + n22 * n33 * n44;
        te[ 4 ] = n14 * n33 * n42 - n13 * n34 * n42 - n14 * n32 * n43 + n12 * n34 * n43 + n13 * n32 * n44 - n12 * n33 * n44;
        te[ 8 ] = n13 * n24 * n42 - n14 * n23 * n42 + n14 * n22 * n43 - n12 * n24 * n43 - n13 * n22 * n44 + n12 * n23 * n44;
        te[ 12 ] = n14 * n23 * n32 - n13 * n24 * n32 - n14 * n22 * n33 + n12 * n24 * n33 + n13 * n22 * n34 - n12 * n23 * n34;
        te[ 1 ] = n24 * n33 * n41 - n23 * n34 * n41 - n24 * n31 * n43 + n21 * n34 * n43 + n23 * n31 * n44 - n21 * n33 * n44;
        te[ 5 ] = n13 * n34 * n41 - n14 * n33 * n41 + n14 * n31 * n43 - n11 * n34 * n43 - n13 * n31 * n44 + n11 * n33 * n44;
        te[ 9 ] = n14 * n23 * n41 - n13 * n24 * n41 - n14 * n21 * n43 + n11 * n24 * n43 + n13 * n21 * n44 - n11 * n23 * n44;
        te[ 13 ] = n13 * n24 * n31 - n14 * n23 * n31 + n14 * n21 * n33 - n11 * n24 * n33 - n13 * n21 * n34 + n11 * n23 * n34;
        te[ 2 ] = n22 * n34 * n41 - n24 * n32 * n41 + n24 * n31 * n42 - n21 * n34 * n42 - n22 * n31 * n44 + n21 * n32 * n44;
        te[ 6 ] = n14 * n32 * n41 - n12 * n34 * n41 - n14 * n31 * n42 + n11 * n34 * n42 + n12 * n31 * n44 - n11 * n32 * n44;
        te[ 10 ] = n12 * n24 * n41 - n14 * n22 * n41 + n14 * n21 * n42 - n11 * n24 * n42 - n12 * n21 * n44 + n11 * n22 * n44;
        te[ 14 ] = n14 * n22 * n31 - n12 * n24 * n31 - n14 * n21 * n32 + n11 * n24 * n32 + n12 * n21 * n34 - n11 * n22 * n34;
        te[ 3 ] = n23 * n32 * n41 - n22 * n33 * n41 - n23 * n31 * n42 + n21 * n33 * n42 + n22 * n31 * n43 - n21 * n32 * n43;
        te[ 7 ] = n12 * n33 * n41 - n13 * n32 * n41 + n13 * n31 * n42 - n11 * n33 * n42 - n12 * n31 * n43 + n11 * n32 * n43;
        te[ 11 ] = n13 * n22 * n41 - n12 * n23 * n41 - n13 * n21 * n42 + n11 * n23 * n42 + n12 * n21 * n43 - n11 * n22 * n43;
        te[ 15 ] = n12 * n23 * n31 - n13 * n22 * n31 + n13 * n21 * n32 - n11 * n23 * n32 - n12 * n21 * n33 + n11 * n22 * n33;

1. /*! Calculate the Christoffel symbols of the second kind at position 'pos'.
2. *
3. *  \param pos : pointer to position.
4. */
5. bool MetricHartleThorneGB::calculateChristoffels( const double* pos ) {
6.     double M = mMass;
7.     double a = mAngmom;
8.     double eta = mEta;
10.     double r = pos[1];
11.     double theta = pos[2];
12.     calcFuncDiff(pos);
14.     double t1 = sigma;      // Sigma(r,theta);
15.     double t2 = t1*t1;
16.     double t3 = t2*t2;
17.     double t4 = 1/t3;
18.     double t5 = delta;      // Delta(r);
19.     double t6 = t5*t5;
20.     double t7 = t4*t6;
21.     double t8 = a*a;
22.     double t9 = eta*t8;
23.     double t10 = hrr;       // hrr(r,theta);
24.     double t11 = t1*t10;
25.     double t14 = 1/(t5+t9*t11);
26.     double t15 = M*t2;
27.     double t17 = M*r;
28.     double t18 = dsigmadr;      // diff(Sigma(r,theta),r);
29.     double t19 = t18*t1;
30.     double t22 = t2*t1;
31.     double t23 = dhttdr;        // diff(htt(r,theta),r);
32.     double t24 = t22*t23;
33.     double t26 = t9*t2;
34.     double t27 = t23*M;
35.     double t28 = t27*r;
36.     double t31 = t9*t1;
37.     double t32 = htt;           // htt(r,theta);
38.     double t33 = t32*t18;
39.     double t34 = t33*t17;
40.     double t37 = t2*t32;
41.     double t41 = M*M;
42.     double t42 = t23*t41;
43.     double t43 = r*r;
44.     double t44 = t42*t43;
45.     double t47 = t32*t41;
46.     double t48 = t47*r;
47.     double t51 = t8*t8;
48.     double t52 = eta*t51;
49.     double t53 = t52*t1;
50.     double t54 = t41*r;
51.     double t55 = sin(theta);
52.     double t56 = t55*t55;
53.     double t57 = t56*t56;
54.     double t58 = hphph;           // hpp(r,theta);
55.     double t59 = t57*t58;
56.     double t60 = t54*t59;
57.     double t63 = t41*t43;
58.     double t64 = dhppdr;        // diff(hpp(r,theta),r);
59.     double t65 = t57*t64;
60.     double t70 = t47*t43;
61.     double t74 = t63*t59;
62.     double t77 = 2.0*t15-2.0*t17*t19+t9*t24-4.0*t26*t28+4.0*t31*t34-4.0*t9*t37*M+4.0*t31*t44+8.0*t31*t48+8.0*t53*t60+4.0*t53*t63*t65-8.0*t9*t18*t70-8.0*t52*t18*t74;
63.     double t81 = hthth;           // haa(r,theta);
64.     double t82 = t1*t81;
65.     double t85 = 1/(1.0+t9*t82);
66.     double t86 = t4*t85;
67.     double t87 = dsigmadth;     // diff(Sigma(r,theta),theta);
68.     double t88 = t87*t1;
69.     double t89 = t17*t88;
70.     double t91 = dhttdth;       // diff(htt(r,theta),theta);
71.     double t92 = t22*t91;
72.     double t94 = t91*M;
73.     double t95 = t94*r;
74.     double t98 = t32*t87;
75.     double t99 = t98*t17;
76.     double t102 = t91*t41;
77.     double t103 = t102*t43;
78.     double t106 = t1*t41;
79.     double t107 = t52*t106;
80.     double t108 = t56*t55;
81.     double t109 = t43*t108;
82.     double t110 = cos(theta);
83.     double t111 = t58*t110;
84.     double t112 = t109*t111;
85.     double t115 = dhppdth;      // diff(hpp(r,theta),theta);
86.     double t116 = t57*t115;
87.     double t129 = 1/t2;
88.     double t130 = t8*t22;
89.     double t133 = t43*t22;
90.     double t134 = t133*M;
91.     double t136 = t56*t58;
92.     double t137 = t52*t136;
93.     double t138 = t43*r;
94.     double t139 = t138*t2;
95.     double t140 = M*t18;
96.     double t141 = t139*t140;
97.     double t144 = eta*eta;
98.     double t145 = t51*t8;
99.     double t146 = t144*t145;
100.     double t147 = t146*t136;
101.     double t148 = t33*M;
102.     double t152 = t138*t22;
103.     double t153 = t152*t27;
104.     double t158 = 16.0*t147*t139*t47;
105.     double t159 = t8*t2;
106.     double t160 = t17*t18;
107.     double t163 = t146*t56;
108.     double t164 = t58*t43;
109.     double t165 = t3*t23;
110.     double t170 = t8*t32;
111.     double t172 = t133*eta*t170*M;
112.     double t174 = t51*t22;
113.     double t175 = t174*eta;
114.     double t179 = t51*t51;
115.     double t180 = t144*t179;
116.     double t181 = t180*t56;
117.     double t182 = t58*t22;
118.     double t183 = t32*M;
119.     double t187 = t51*t2;
120.     double t188 = t187*eta;
121.     double t191 = eta*t32;
122.     double t195 = eta*t145;
123.     double t196 = t195*t136;
124.     double t197 = r*t18;
125.     double t201 = t58*t3;
126.     double t204 = t57*t56;
127.     double t205 = t58*t58;
128.     double t206 = t204*t205;
129.     double t207 = t180*t206;
130.     double t208 = t43*t43;
131.     double t209 = t208*t1;
132.     double t210 = t18*t41;
133.     double t213 = 8.0*t207*t209*t210;
134.     double t214 = t180*t136;
135.     double t218 = t146*t59;
136.     double t219 = t208*t2;
137.     double t220 = t219*t42;
138.     double t221 = t218*t220;
139.     double t223 = t47*t18;
140.     double t224 = t209*t223;
141.     double t226 = 8.0*t218*t224;
142.     double t228 = t51*t3;
143.     double t229 = eta*t23;
144.     double t231 = -4.0*t181*t182*t183+8.0*t188*t48-4.0*t174*t191*M-2.0*t196*t15*t197+t181*t201*t23-t213+4.0*t214*t37*t160-12.0*t221+t226-2.0*t141+t228*t229;
145.     double t233 = t208*r;
146.     double t234 = t233*t145;
147.     double t235 = t41*M;
148.     double t236 = t235*t57;
149.     double t237 = t234*t236;
150.     double t238 = t144*t32;
151.     double t239 = t1*t64;
152.     double t240 = t238*t239;
153.     double t242 = 8.0*t237*t240;
154.     double t245 = t235*t32;
155.     double t248 = 8.0*t218*t209*t245;
156.     double t249 = t144*t51;
157.     double t250 = t249*t136;
158.     double t251 = t208*t43;
159.     double t252 = t251*t1;
160.     double t256 = t152*eta;
161.     double t257 = t8*t23;
162.     double t258 = t257*M;
163.     double t261 = t219*eta;
164.     double t265 = t139*eta;
165.     double t268 = 8.0*t265*t170*t41;
166.     double t273 = t146*t204;
167.     double t274 = t205*t233;
168.     double t278 = 4.0*t273*t274*t2*t41;
169.     double t280 = t233*t2;
170.     double t282 = t250*t280*t47;
171.     double t284 = t146*t206;
172.     double t287 = 4.0*t284*t252*t210;
173.     double t288 = t204*t58;
174.     double t289 = t180*t288;
175.     double t290 = t235*t138;
176.     double t291 = t290*t33;
177.     double t293 = 8.0*t289*t291;
178.     double t294 = t9*t136;
179.     double t298 = t249*t56;
180.     double t299 = t58*t208;
181.     double t302 = t209*eta;
182.     double t307 = t33*t63;
183.     double t310 = t233*t22;
184.     double t314 = t145*t1;
185.     double t319 = t56*t41;
186.     double t320 = t319*t51;
187.     double t321 = t43*t2;
188.     double t322 = t321*t229;
189.     double t325 = t56*M;
190.     double t326 = t325*t51;
191.     double t327 = r*t22;
192.     double t331 = 8.0*t282-t287-t293-2.0*t294*t280*t140+t298*t299*t165-8.0*t302*t8*t18*t47-8.0*t53*t307-4.0*t250*t310*t27-8.0*t314*eta*t18*t74-4.0*t320*t322+2.0*t326*t327*t229;
193.     double t334 = t43*t145;
194.     double t335 = t334*t319;
195.     double t336 = t32*t32;
196.     double t337 = t144*t336;
197.     double t340 = 4.0*t335*t337*t19;
198.     double t341 = t52*t56;
199.     double t342 = t22*M;
200.     double t343 = t164*t342;
201.     double t344 = t341*t343;
202.     double t346 = t138*t179;
203.     double t347 = t235*t204;
204.     double t348 = t346*t347;
205.     double t350 = 8.0*t348*t240;
206.     double t351 = t208*t145;
207.     double t352 = t41*t57;
208.     double t353 = t351*t352;
209.     double t354 = t2*t64;
210.     double t355 = t238*t354;
211.     double t356 = t353*t355;
212.     double t358 = t195*t56;
213.     double t359 = t182*M;
214.     double t362 = t43*t179;
215.     double t363 = t362*t352;
216.     double t364 = t363*t355;
217.     double t366 = t346*t236;
218.     double t368 = 8.0*t366*t240;
219.     double t369 = t179*t8;
220.     double t370 = t144*t369;
221.     double t371 = t57*t57;
222.     double t372 = t370*t371;
223.     double t373 = t205*t235;
224.     double t374 = t138*t18;
225.     double t377 = 8.0*t372*t373*t374;
226.     double t380 = t56*t144;
227.     double t384 = 4.0*r*t145*t41*t380*t336*t2;
228.     double t385 = t43*t3;
229.     double t386 = t9*t23;
230.     double t389 = t208*t22;
231.     double t391 = t250*t389*t183;
232.     double t393 = t251*t2;
233.     double t397 = t205*t1;
234.     double t398 = t235*t43;
235.     double t401 = 8.0*t372*t397*t398;
236.     double t402 = t180*t59;
237.     double t403 = t1*t235;
238.     double t405 = t403*t138*t23;
239.     double t407 = 8.0*t402*t405;
240.     double t408 = t19*t70;
241.     double t411 = t370*t204;
242.     double t412 = t205*t2;
243.     double t415 = 4.0*t411*t412*t54;
244.     double t416 = t37*t54;
245.     double t421 = t170*t140;
246.     double t425 = 8.0*t289*t405;
247.     double t427 = 8.0*t402*t291;
248.     double t428 = -4.0*t391+4.0*t250*t393*t42+t401+t407-8.0*t214*t408+t415+8.0*t214*t416-16.0*t147*t224+4.0*t265*t421-t425-t427;
249.     double t431 = 8.0*t402*t408;
250.     double t433 = t2*t23*t63;
251.     double t437 = 8.0*t402*t416;
252.     double t438 = t209*t52;
253.     double t442 = t1*t32;
254.     double t443 = t398*t442;
255.     double t445 = 8.0*t289*t443;
256.     double t446 = t402*t433;
257.     double t451 = 8.0*t402*t443;
258.     double t452 = t370*t206;
259.     double t455 = 4.0*t452*t19*t63;
260.     double t456 = t233*t1;
261.     double t460 = 8.0*t218*t456*t235*t23;
262.     double t461 = t233*t235;
263.     double t464 = 8.0*t218*t461*t33;
264.     double t465 = t431+4.0*t214*t433-t437-8.0*t438*t210*t59+t445-12.0*t446+4.0*t289*t433+t451-t455+t460-t464;
265.     double t466 = t180*t371;
266.     double t467 = t205*t208;
267.     double t470 = 8.0*t466*t467*t403;
268.     double t474 = 8.0*t466*t274*t235*t18;
269.     double t475 = t24*t17;
270.     double t481 = t147*t133*t183;
271.     double t488 = t145*t2;
272.     double t489 = t488*eta;
273.     double t492 = t9*t56;
274.     double t494 = t492*t299*t342;
275.     double t496 = t336*t1;
276.     double t499 = 8.0*t163*t398*t496;
277.     double t500 = t336*t18;
278.     double t503 = 8.0*t163*t290*t500;
279.     double t504 = -t470-t474-4.0*t214*t475+4.0*t402*t475-8.0*t481+4.0*t218*t153+4.0*t250*t280*t148+8.0*t489*t60+2.0*t494+t499-t503;
280.     double t509 = M*t138;
281.     double t510 = t509*t1;
282.     double t514 = t52*t57;
283.     double t515 = t58*t138;
284.     double t516 = t1*M;
285.     double t520 = t146*t32;
286.     double t521 = t41*t208;
287.     double t525 = t180*t32;
288.     double t529 = M*t233;
289.     double t538 = t249*t32;
290.     double t546 = t51*t57;
291.     double t557 = t52*t32;
292.     double t558 = t1*t56;
293.     double t562 = t146*t442;
294.     double t573 = 2.0*t510-t159-4.0*t341*t70-4.0*t514*t515*t516-t321+8.0*t520*t521*t136+4.0*t525*t63*t136+2.0*t529*t1*t294+4.0*t9*t47*t208+4.0*t52*t70+4.0*t538*t41*t251*t136+4.0*t510*t137+4.0*t521*eta*t546*t58+4.0*t63*eta*t145*t57*t58+t9*t37*t43+2.0*t557*t558*t17-8.0*t562*t509*t136-4.0*t249*t442*t529*t136-8.0*t520*t521*t59;
295.     double t574 = t180*t442;
296.     double t583 = t9*t32;
297.     double t586 = t516*r;
298.     double t589 = t136*t2;
299.     double t591 = t195*t57;
300.     double t592 = t58*t1;
301.     double t596 = t195*t204;
302.     double t598 = t58*t41*t43;
303.     double t604 = t59*t17;
304.     double t609 = t164*t2;
305.     double t616 = t2*t56;
306.     double t617 = t616*t164;
307.     double t625 = t52*t37;
308.     double t629 = -4.0*t574*t17*t136-8.0*t525*t74+4.0*t525*t288*t63-4.0*t583*t510-4.0*t557*t586-t195*t589-4.0*t591*t592*t17-4.0*t596*t598+4.0*t562*t59*t509+4.0*t574*t604-t492*t299*t2-2.0*t341*t609-2.0*t325*t8*r*t1+2.0*t520*t617+t538*t616*t299+t525*t589+2.0*t586*t196+t625+2.0*t17*t1*t8;
309.     double t631 = 1/(t573+t629);
310.     double t633 = t129*(4.0*t261*t257*t41+2.0*t130*M-8.0*t250*t252*t223-t242-4.0*t172-4.0*t175*t28+8.0*t147*t139*t148+2.0*t163*t164*t165-4.0*t256*t258-t248-t278+t331+t158-2.0*t159*t160+2.0*t358*t359-8.0*t147*t153+t350+4.0*t356+t340+4.0*t344+8.0*t147*t220-t384+t385*t386-t368+t428-4.0*t137*t141+t465+t504+t268+4.0*t188*t44+4.0*t188*t34+2.0*t134+4.0*t364-t377+t231)*t631/2.0;
311.     double t634 = a*M;
312.     double t635 = t515*t24;
313.     double t637 = M*t43;
314.     double t639 = t56*t64;
315.     double t640 = t52*t639;
316.     double t643 = t456*t42;
317.     double t646 = t249*t47;
318.     double t650 = M*t208;
319.     double t652 = t9*t639;
320.     double t658 = t467*t106;
321.     double t665 = t274*t210;
322.     double t668 = t180*t204;
323.     double t670 = t41*t138;
324.     double t671 = t205*t18*t670;
325.     double t674 = t515*t15;
326.     double t677 = t146*t47;
327.     double t678 = t138*t1;
328.     double t682 = t136*t22;
329.     double t684 = t1*t23;
330.     double t685 = t670*t684;
331.     double t688 = t298*t635-2.0*t637*t2*t640+4.0*t250*t643-4.0*t646*t456*t639-2.0*t650*t2*t652-4.0*t250*t209*t47-4.0*t273*t658-4.0*t250*t233*t18*t47-4.0*t273*t665-4.0*t668*t671-4.0*t492*t674-4.0*t677*t678*t639+t52*t682-4.0*t218*t685;
332.     double t689 = t321*t27;
333.     double t692 = t43*t18;
334.     double t697 = t146*t37;
335.     double t698 = r*t56;
336.     double t699 = t58*t18;
337.     double t700 = t698*t699;
338.     double t702 = t43*t57;
339.     double t703 = t64*M;
340.     double t704 = t702*t703;
341.     double t707 = t249*t336;
342.     double t708 = t1*t43;
343.     double t709 = t708*t140;
344.     double t715 = t327*t639;
345.     double t720 = t146*t183;
346.     double t721 = t321*t639;
347.     double t724 = t249*t183;
348.     double t725 = t138*t56;
349.     double t726 = t58*t2;
350.     double t732 = t219*t27;
351.     double t737 = t33*t670;
352.     double t740 = -4.0*t147*t689+4.0*t562*t692*t59*M+t697*t700-2.0*t697*t704+4.0*t707*t709-4.0*t677*t374*t59-t520*t715+4.0*t677*t708*t59+4.0*t720*t721+8.0*t724*t725*t726+2.0*t218*t689-4.0*t250*t732-4.0*t697*t604-4.0*t147*t737;
353.     double t742 = t183*r;
354.     double t745 = t32*t22;
355.     double t748 = t219*t639;
356.     double t751 = t63*t1;
357.     double t760 = t9*t138;
358.     double t761 = t22*t56;
359.     double t762 = t761*t64;
360.     double t769 = t761*t164;
361.     double t772 = t52*r;
362.     double t774 = t152*t639;
363.     double t776 = t249*t37;
364.     double t779 = t52*t43;
365.     double t780 = t2*t57;
366.     double t784 = 4.0*t26*t742-2.0*t9*t745+4.0*t724*t748+4.0*t707*t751-4.0*t707*t670*t18+4.0*t677*t678*t65+t760*t762+4.0*t188*t604+3.0*t492*t164*t22-3.0*t538*t769+t772*t762-t538*t774+t776*t725*t699+2.0*t779*t780*t703;
367.     double t785 = t197*t2;
368.     double t788 = t397*t63;
369.     double t791 = t15*r;
370.     double t794 = t616*t699;
371.     double t796 = t9*r;
372.     double t802 = t182*r*t23;
373.     double t806 = t442*t63;
374.     double t809 = t1*eta;
375.     double t810 = t51*t43;
376.     double t813 = t58*M;
377.     double t814 = t18*t57*t813;
378.     double t817 = t692*t183;
379.     double t820 = t336*t22;
380.     double t822 = -t707*t785-t520*t682+4.0*t668*t788-4.0*t707*t791-t772*t794+2.0*t796*t37*t18-t760*t794+t163*t802+4.0*t147*t685+4.0*t147*t806-4.0*t809*t810*t814-4.0*t31*t817+t22+t249*t820-t785;
381.     double t827 = t634*t129*(t688+t740+t784+t822)*t631;
382.     double t828 = t41*t108;
383.     double t830 = t2*t110;
384.     double t831 = t830*t58;
385.     double t832 = t238*t831;
386.     double t835 = t57*t108;
387.     double t836 = t835*t205;
388.     double t837 = t370*t836;
389.     double t838 = t138*t110;
390.     double t846 = t1*t91;
391.     double t847 = t290*t846;
392.     double t850 = t2*t91;
393.     double t851 = t63*t850;
394.     double t854 = t371*t55;
395.     double t855 = t370*t854;
396.     double t856 = t41*t41;
397.     double t858 = t208*t110;
398.     double t864 = t98*t63;
399.     double t871 = t9*t91;
400.     double t873 = M*t87;
401.     double t874 = t139*t873;
402.     double t876 = t290*t98;
403.     double t879 = t98*t41;
404.     double t888 = t98*M;
405.     double t889 = t139*t888;
406.     double t892 = t152*t94;
407.     double t896 = t43*t8*t41;
408.     double t904 = 24.0*t351*t828*t832-16.0*t837*t403*t838-8.0*t372*t373*t138*t87-8.0*t289*t847+4.0*t289*t851-32.0*t855*t205*t856*t858+8.0*t402*t847-8.0*t53*t864-8.0*t314*eta*t87*t74+t385*t871-2.0*t874-8.0*t402*t876-8.0*t250*t252*t879+t181*t201*t91+4.0*t250*t393*t102+8.0*t147*t889-8.0*t147*t892-8.0*t896*t55*t110*t2+4.0*t250*t280*t888;
409.     double t905 = t219*t102;
410.     double t912 = t17*t87;
411.     double t921 = t87*t41;
412.     double t929 = t209*t879;
413.     double t937 = t1*t115;
414.     double t938 = t238*t937;
415.     double t941 = t2*t115;
416.     double t942 = t238*t941;
417.     double t945 = t41*t55;
418.     double t954 = t351*t856;
419.     double t956 = t336*t110;
420.     double t963 = -12.0*t218*t905+8.0*t218*t456*t235*t91-2.0*t159*t912+4.0*t188*t99+4.0*t188*t103-8.0*t289*t876-8.0*t207*t209*t921+4.0*t218*t892+8.0*t147*t905-16.0*t147*t929-4.0*t284*t252*t921-4.0*t137*t874+8.0*t348*t938+4.0*t353*t942-8.0*t334*t945*t337*t830-8.0*t366*t938+4.0*t363*t942-32.0*t954*t55*t144*t956-8.0*t438*t921*t59;
421.     double t972 = eta*t91;
422.     double t976 = t57*t55;
423.     double t979 = t208*t32;
424.     double t984 = r*t87;
425.     double t988 = t92*t17;
426.     double t1009 = t3*t91;
427.     double t1013 = t170*t873;
428.     double t1020 = t8*t91;
429.     double t1024 = t43*t32;
430.     double t1026 = t106*t1024*t87;
431.     double t1029 = -8.0*t237*t938-4.0*t250*t310*t94+8.0*t218*t929+2.0*t326*t327*t972-64.0*t180*t976*t856*t979*t111+t228*t972-2.0*t196*t15*t984-4.0*t214*t988-4.0*t320*t321*t972-4.0*t452*t88*t63+4.0*t402*t988-12.0*t402*t851-8.0*t219*t52*t828*t111-8.0*t302*t8*t87*t47+2.0*t163*t164*t1009+4.0*t265*t1013-8.0*t466*t274*t235*t87+4.0*t261*t1020*t41+8.0*t402*t1026;
432.     double t1030 = t976*t235;
433.     double t1031 = t145*t138;
434.     double t1033 = t809*t111;
435.     double t1036 = t336*t87;
436.     double t1037 = t290*t1036;
437.     double t1045 = t235*t108;
438.     double t1047 = t110*t1;
439.     double t1048 = t1047*t58;
440.     double t1049 = t238*t1048;
441.     double t1052 = t810*t945;
442.     double t1053 = eta*t110;
443.     double t1054 = t1053*t37;
444.     double t1057 = t138*t51;
445.     double t1058 = t235*t55;
446.     double t1060 = t1053*t442;
447.     double t1063 = t180*t836;
448.     double t1071 = t346*t1030;
449.     double t1077 = t1020*M;
450.     double t1094 = t37*t912;
451.     double t1103 = -32.0*t1030*t1031*t1033-8.0*t163*t1037-2.0*t294*t280*t873+t298*t299*t1009-48.0*t346*t1045*t1049+16.0*t1052*t1054-32.0*t1057*t1058*t1060-16.0*t1063*t456*t235*t110+4.0*t335*t337*t88+64.0*t1071*t1049+24.0*t362*t828*t832-4.0*t256*t1077-8.0*t214*t1026-8.0*t488*eta*t41*t112+4.0*t214*t851-48.0*t234*t1045*t1049-8.0*t218*t461*t98+4.0*t214*t1094+32.0*t1031*t1058*t337*t1047-4.0*t175*t95;
452.     double t1108 = t129*(t904+t963+t1029+t1103)*t631/2.0;
453.     double t1110 = t55*t87;
454.     double t1111 = t1110*t2;
455.     double t1112 = r*t110;
456.     double t1116 = t63*t1047;
457.     double t1122 = t110*t22;
458.     double t1126 = t110*t57;
459.     double t1127 = t813*r;
460.     double t1131 = t976*t115;
461.     double t1132 = t1131*t17;
462.     double t1139 = t976*t87;
463.     double t1143 = t2*t108;
464.     double t1144 = t58*t87;
465.     double t1147 = t22*t108;
466.     double t1148 = t115*t43;
467.     double t1151 = t108*t58;
468.     double t1152 = t43*t87;
469.     double t1158 = t52*t110;
470.     double t1161 = t249*t1151;
471.     double t1162 = t139*t94;
472.     double t1165 = -t1111+16.0*t720*t1112*t589+24.0*t707*t1116+16.0*t724*t838*t589-4.0*t9*t1122*t32-12.0*t697*t1126*t1127-2.0*t697*t1132+4.0*t249*t496*t1110*t17+4.0*t562*t1139*t1127+t520*t1143*t1144-t538*t1147*t1148+t776*t1151*t1152-4.0*t677*t109*t937+4.0*t1158*t682-4.0*t1161*t1162;
473.     double t1166 = t1147*t115;
474.     double t1168 = t9*t110;
475.     double t1172 = t138*t57;
476.     double t1176 = t146*t1151;
477.     double t1177 = t846*t63;
478.     double t1180 = t108*t115;
479.     double t1181 = t52*t1180;
480.     double t1184 = t110*t32;
481.     double t1189 = t9*t1180;
482.     double t1192 = t2*eta;
483.     double t1195 = t51*t110*t136;
484.     double t1198 = t209*t102;
485.     double t1201 = t43*t110;
486.     double t1202 = t558*t58;
487.     double t1213 = t8*t110;
488.     double t1214 = t1213*t136;
489.     double t1217 = t850*t17;
490.     double t1220 = t1*t57;
491.     double t1225 = t146*t108;
492.     double t1226 = t182*t91;
493.     double t1228 = -t520*t1166+4.0*t1168*t769-32.0*t146*t245*t1172*t111+4.0*t1176*t1177-2.0*t791*t1181-16.0*t751*t9*t1184-2.0*t509*t2*t1189-8.0*t17*t1192*t1195+4.0*t1161*t1198-16.0*t677*t1201*t1202-16.0*t646*t858*t1202-16.0*t707*t290*t110-8.0*t509*t1192*t1214-4.0*t1176*t1217+40.0*t677*t1201*t1220*t58+t1225*t1226;
494.     double t1237 = t976*t58;
495.     double t1238 = t146*t1237;
496.     double t1242 = t205*t87;
497.     double t1243 = t1242*t63;
498.     double t1246 = t37*t17;
499.     double t1251 = t373*t838;
500.     double t1258 = t51*t976;
501.     double t1259 = t809*t1258;
502.     double t1260 = t1144*t17;
503.     double t1263 = t138*t108;
504.     double t1272 = t43*t976;
505.     double t1276 = r*t108;
506.     double t1283 = t17*t830;
507.     double t1285 = -16.0*t107*t702*t111-4.0*t1161*t208*t87*t47+2.0*t1238*t1217-4.0*t180*t835*t1243+16.0*t1168*t1246-4.0*t1176*t864-16.0*t466*t1251-4.0*t809*t8*t55*t99-4.0*t1259*t1260+4.0*t724*t1263*t941-4.0*t1238*t1177-4.0*t707*t63*t1110-4.0*t677*t1272*t1144+4.0*t720*t1276*t941+4.0*t677*t1272*t937-4.0*t1283;
508.     double t1287 = t467*t921;
509.     double t1290 = t208*t108;
510.     double t1298 = t110*t56;
511.     double t1305 = t9*t55;
512.     double t1309 = t9*t108;
513.     double t1310 = t22*t115;
514.     double t1311 = t1310*t43;
515.     double t1315 = t52*t108;
516.     double t1318 = t52*t976;
517.     double t1319 = t941*t17;
518.     double t1331 = t164*t92;
519.     double t1333 = -4.0*t146*t835*t1287-4.0*t646*t1290*t937-4.0*t520*t1122*t136-4.0*t249*t745*t1298*t164+2.0*t249*t820*t110+2.0*t1305*t37*t87+t1309*t1311-t1309*t726*t1152-t1315*t726*t87+2.0*t1318*t1319+2.0*t1122+12.0*t52*t830*t604-12.0*t707*t1283-t707*t1111+t52*t1166+t249*t108*t1331;
520.     double t1338 = t631/t55;
521.     double t1340 = t17*a*t129*(t1165+t1228+t1285+t1333)*t1338;
522.     double t1342 = t197*t1;
523.     double t1351 = r*t2;
524.     double t1356 = t9*t43;
525.     double t1360 = t616*t64;
526.     double t1362 = t558*t699;
527.     double t1373 = -t2+t1342+t9*t37-4.0*t583*t586-3.0*t492*t609-t52*t589-4.0*t53*t604+t9*t1351*t23-t796*t442*t18-2.0*t1356*t684*M-t760*t1360+t760*t1362-t772*t1360+t772*t1362-2.0*t779*t1220*t703+4.0*t1356*t148+4.0*t779*t814;
528.     double t1376 = t7*t14*t634*t56*t1373;
529.     double t1381 = t830*t32;
530.     double t1408 = t592*t87;
531.     double t1417 = -2.0*t830+t1110*t1+2.0*t9*t1381-4.0*t1168*t442*t17-4.0*t1168*t617-4.0*t1158*t589-12.0*t52*t1047*t604+t9*t55*t2*t91-t1305*t442*t87-2.0*t1305*t846*t17-t1309*t941*t43+t1309*t592*t1152-t52*t1143*t115+t1315*t1408-2.0*t1318*t937*t17+4.0*t1305*t99+4.0*t1318*t1260;
532.     double t1420 = t86*M*a*r*t55*t1417;
533.     double t1421 = 1/t1;
534.     double t1425 = deltadr;         // diff(Delta(r),r);
535.     double t1436 = dhrrdr;          // diff(hrr(r,theta),r);
536.     double t1449 = dhrrdth;         // diff(hrr(r,theta),theta);
537.     double t1452 = t87*t5+2.0*t9*t11*t87+t9*t2*t1449;
538.     double t1457 = t1421*t14;
539.     double t1459 = t1457*t1452/2.0;
540.     double t1463 = dhaadr;        // diff(haa(r,theta),r);
541.     double t1466 = t18+2.0*t9*t82*t18+t9*t2*t1463;
542.     double t1469 = t1466*t85*t1421/2.0;
543.     double t1470 = t634*t56;
544.     double t1472 = t146*t352;
545.     double t1473 = t233*t32;
546.     double t1492 = t370*t57;
547.     double t1513 = t1492*t412*t197+4.0*t214*t708*t148-t214*t1351*t33+t181*t802-4.0*t411*t412*t17+4.0*t466*t665+t181*t182*t32-4.0*t214*t1246+4.0*t218*t732+t174*t191-t130;
548.     double t1515 = t180*t352;
549.     double t1522 = t138*t32;
550.     double t1530 = t8*t56;
551.     double t1559 = r*t32;
552.     double t1567 = t43*t23;
553.     double t1575 = -t370*t57*t205*t22-t133*t583+4.0*t265*t546*t813+t152*t386-t265*t170*t18-t188*t1559*t18+t174*eta*r*t23-4.0*t489*t604-2.0*t188*t1567*M+4.0*t438*t814+4.0*t302*t421;
554.     double t1580 = t145*t22;
555.     double t1585 = t180*t183;
556.     double t1590 = t204*t64;
557.     double t1603 = t761*t64*t233;
558.     double t1609 = t180*t37;
559.     double t1612 = t251*t56;
560.     double t1613 = t1612*t354;
561.     double t1626 = t51*t56;
562.     double t1641 = -4.0*t697*t65*t650-4.0*t1609*t704+2.0*t724*t1613-r*t144*t179*t745*t639-2.0*t261*t258+2.0*t280*eta*t1530*t699+4.0*t265*t1626*t699+4.0*t466*t658+8.0*t147*t209*t148+8.0*t207*t209*t140-2.0*t147*t139*t33;
563.     double t1651 = t180*t57;
564.     double t1652 = t205*t43;
565.     double t1661 = t2*t18;
566.     double t1668 = t249*t57;
567.     double t1685 = t146*t57;
568.     double t1695 = t205*t251;
569.     double t1706 = t298*t58*t233*t24-t250*t280*t33-2.0*t250*t393*t27+3.0*t1685*t274*t1661-3.0*t298*t299*t745+4.0*t250*t280*t183+t1668*t1695*t22+t1685*t467*t22+4.0*t163*t670*t500+t159*t197+2.0*t489*t700;
570.     double t1712 = t1470*t129*(-8.0*t218*t139*t183+4.0*t250*t252*t148-t1651*t1652*t22+t133-2.0*t163*t164*t745+t1668*t205*t208*t138*t1661+4.0*t1472*t1473*t699+4.0*t1515*t1522*t239+4.0*t1515*t1522*t699-4.0*t574*t1590*t670+3.0*t1651*t205*t138*t1661-4.0*t1515*t1024*t592+4.0*t1472*t979*t592+4.0*t1472*t1473*t239+4.0*t284*t252*t140+t1706+4.0*t273*t274*t15-4.0*t163*t63*t496-4.0*t147*t732+4.0*t402*t689-4.0*t188*t742+4.0*t720*t748+4.0*t372*t671-4.0*t372*t788+2.0*t326*t322+2.0*t389*eta*t1530*t58+4.0*t314*eta*t43*t814-2.0*t1580*eta*t56*t58+4.0*t452*t709-4.0*t402*t685+4.0*t289*t685-2.0*t214*t689+4.0*t289*t737-4.0*t289*t806-t538*t1603+t1641+2.0*t163*t635+4.0*t53*t817-4.0*t218*t643+t1575-2.0*t520*t774+t1513+2.0*t1585*t721+t139*t18)*t631;
571.     double t1718 = t146*t745;
572.     double t1723 = t3*t56;
573.     double t1724 = t1723*t64;
574.     double t1727 = t509*t22;
575.     double t1731 = -4.0*t596*t354*t63-2.0*r*t3+4.0*t1718*t65*t509+4.0*t134+t158+t525*t1724-8.0*t172+4.0*t1727*t640+t213+4.0*t221+t226;
576.     double t1736 = -12.0*t514*t343+t242+t248+t268+t278+16.0*t282+t287+t293-t340+8.0*t344-t350;
577.     double t1748 = t180*t745;
578.     double t1752 = -12.0*t356-12.0*t364+t368+t377+t384-16.0*t391+4.0*t525*t22*t57*t813+t538*t1723*t64*t208-t401+4.0*t1748*t65*t17-t407;
579.     double t1778 = t208*t51;
580.     double t1779 = t1778*t41;
581.     double t1780 = t57*eta;
582.     double t1781 = t1780*t354;
583.     double t1784 = t334*t41;
584.     double t1791 = -8.0*t720*t774-32.0*t677*t1172*t726-4.0*t1609*t59*t160+4.0*t538*t1723*t515+2.0*t520*t1723*t64*t43+2.0*t557*t761*M+8.0*t1609*t288*t54-8.0*t574*t204*t18*t598+4.0*t1779*t1781+4.0*t1784*t1781-4.0*t1585*t715+8.0*t677*t748;
585.     double t1797 = t325*t8;
586.     double t1800 = t180*t47;
587.     double t1805 = t699*t63;
588.     double t1808 = t17*t22;
589.     double t1814 = t529*t22;
590.     double t1821 = -t415+8.0*t809*t1626*t307+2.0*t1797*t785+4.0*t1800*t721+t425+t427+8.0*t809*t145*t204*t1805+2.0*t1808*t195*t639+4.0*t646*t1613+2.0*t1814*t652+16.0*t1057*t41*t1780*t726;
591.     double t1840 = t3*t64;
592.     double t1843 = t22*t64;
593.     double t1853 = 4.0*t1609*t1590*t63-4.0*t724*t1603+2.0*t9*t32*t3*r+t431+12.0*t1718*t59*t637-4.0*t697*t59*t374*M-8.0*t341*t416-t492*t1840*t208-4.0*t591*t1843*t17-8.0*t596*t726*t54-4.0*t514*t1843*t509;
594.     double t1856 = t515*t140;
595.     double t1879 = t699*t17;
596.     double t1882 = 4.0*t52*t780*t1856-t195*t1724-4.0*t492*t201*t138-2.0*t341*t1840*t43-4.0*t591*t359-4.0*t341*t201*r+4.0*t520*t1723*t58*r-2.0*t625*t325*t197+4.0*t195*t780*t1879-t437-t445;
597.     double t1888 = 4.0*t446-t451+t455-t460+t464+t470+t474-2.0*t325*t130-16.0*t481+8.0*t494-t499+t503;
598.     double t1894 = t129*(t1731+t1736+t1752+t1791+t1821+t1853+t1882+t1888)*t631/2.0;
599.     double t1901 = dhaadth;         // diff(haa(r,theta),theta);
600.     double t1917 = t43*t336;
601.     double t1942 = t146*t828;
602.     double t1946 = M*t110;
603.     double t1947 = t1946*t32;
604.     double t1948 = t139*t1947;
605.     double t1955 = t180*t1237;
606.     double t1968 = -t321*t87-t159*t87-t181*t1226+t525*t761*t115+4.0*t147*t1162+8.0*t146*t945*t1917*t1047+4.0*t697*t116*t509-2.0*t1585*t698*t941-4.0*t720*t725*t941-4.0*t163*t63*t1036-4.0*t289*t1177-4.0*t289*t864-4.0*t372*t1243-4.0*t678*eta*t1013-32.0*t1942*t979*t1048+16.0*t1176*t1948-32.0*t180*t828*t1024*t1048+24.0*t1955*t63*t1047*t32+4.0*t574*t204*t115*t63-2.0*t724*t233*t56*t941;
607.     double t1971 = t139*t52;
608.     double t1977 = t321*eta;
609.     double t1980 = t976*t41;
610.     double t2003 = t55*t58;
611.     double t2005 = t830*t742;
612.     double t2008 = t976*t205;
613.     double t2015 = t55*M;
614.     double t2017 = t1351*t110;
615.     double t2022 = t314*t1780;
616.     double t2034 = -4.0*t466*t1287+2.0*t180*t55*t182*t1184-4.0*t180*t2003*t2005-4.0*t370*t2008*t1283+2.0*t520*t761*t1148-4.0*t2015*t8*t2017-4.0*t53*t99-4.0*t2022*t1260-2.0*t489*t136*t87-8.0*t488*t1053*t1151*t17+t187*t191*t87;
617.     double t2046 = t51*r;
618.     double t2069 = t146*t2003;
619.     double t2089 = t2*t87;
620.     double t2102 = -4.0*t1977*t1626*t1144+2.0*t188*t95-16.0*t855*t1251-4.0*t1472*t979*t937+t250*t219*t98+2.0*t250*t280*t94-t1668*t1695*t2089-3.0*t1685*t467*t2089-t298*t299*t92-2.0*t163*t1331+2.0*t147*t321*t98;
621.     double t2117 = t115*t208;
622.     double t2133 = t249*t2003;
623.     double t2162 = 2.0*t2133*t389*t1184-4.0*t214*t88*t742-4.0*t452*t89+16.0*t180*t1151*t2005+t181*t726*t98+2.0*t214*t1217-2.0*t326*t1351*t972-4.0*t2133*t280*t1947-4.0*t146*t2008*t280*t1946+2.0*t265*t1077-t174*t972;
624.     double t2168 = t1470*r*t129*(-8.0*t1971*t110*t108*t813+t1977*t170*t87-16.0*t1980*t334*t1033-4.0*t250*t456*t888-3.0*t1651*t1652*t2089-8.0*t180*t2008*t139*t1946-32.0*t180*t1030*t1522*t111-16.0*t146*t55*t290*t956-16.0*t837*t1116-16.0*t1063*t209*t41*t110+t538*t761*t2117-8.0*t207*t678*t873-8.0*t147*t678*t888+4.0*t2069*t133*t1184-4.0*t284*t456*t873-t1492*t412*t87+4.0*t1609*t116*t17+4.0*t2015*t2046*t1054-4.0*t1515*t1024*t1144-4.0*t1515*t1024*t937-4.0*t1472*t979*t1144-2.0*t261*t1530*t1144+4.0*t218*t1198+t2034+t1968-8.0*t2069*t1948-16.0*t1052*t1060-4.0*t678*t52*t57*t87*t813+t2102-t133*t871+t2162-4.0*t402*t1217-4.0*t218*t1162+4.0*t402*t1177)*t631;
625.     double t2169 = t3*t108;
626.     double t2180 = t108*M;
627.     double t2181 = t2180*t8;
628.     double t2190 = t835*t41;
629.     double t2199 = t385*eta;
630.     double t2203 = t976*eta;
631.     double t2204 = t2203*t941;
632.     double t2212 = 2.0*t520*t2169*t1148+4.0*t181*t201*t1184-2.0*t385*t110+t538*t2169*t2117+2.0*t2181*t1351*t87-4.0*t1580*eta*t1132+8.0*t1225*t1037-4.0*t2190*t145*t321*eta*t115-8.0*t175*t1184*t17+2.0*t2199*t1213*t32+4.0*t1779*t2204+2.0*t1808*t195*t1180+4.0*t1727*t1181;
633.     double t2226 = t208*t3;
634.     double t2227 = t2226*eta;
635.     double t2238 = t138*t369;
636.     double t2241 = t854*t144*t1242;
637.     double t2258 = 32.0*t954*t380*t956+32.0*t208*t369*t856*t371*t56*t144*t205*t110+4.0*t509*t1122-4.0*t2227*t1214-8.0*t2199*t1195-4.0*t256*t1258*t115*M+2.0*t1814*t1189+8.0*t2238*t235*t2241-2.0*t8*t3*t110-8.0*t256*t1213*t183+8.0*t583*t521*t830+8.0*t557*t63*t830+4.0*t1784*t2204;
638.     double t2260 = t233*t179;
639.     double t2274 = t208*t179;
640.     double t2284 = t1053*t726;
641.     double t2299 = t342*t1112*t32;
642.     double t2304 = 8.0*t2260*t235*t2241-2.0*t2199*t51*t108*t115+8.0*t896*t1298*t2-t2227*t8*t108*t115+64.0*t2274*t856*t204*t238*t111-4.0*t1942*t1917*t88+32.0*t1778*t352*t2284+4.0*t1238*t905-8.0*t720*t1263*t1310+4.0*t1800*t109*t941-4.0*t1585*t1276*t1310+24.0*t402*t2299+4.0*t1748*t1132;
643.     double t2314 = t180*t835*t58;
644.     double t2322 = t152*t1947;
645.     double t2346 = 4.0*t1609*t835*t115*t63-4.0*t724*t233*t108*t1310-8.0*t2314*t1026-72.0*t1515*t1024*t831+8.0*t2314*t847+24.0*t218*t2322-16.0*t214*t2299-4.0*t1238*t889+8.0*t146*t319*t1917*t830+32.0*t677*t208*t56*t831+8.0*t2190*t334*t809*t1144+4.0*t1718*t1131*t509+t525*t2169*t115;
646.     double t2366 = t234*t1030;
647.     double t2367 = t238*t1144;
648.     double t2378 = t22*eta;
649.     double t2391 = -72.0*t1472*t979*t831-4.0*t1955*t1094+32.0*t289*t63*t1381+16.0*t646*t1612*t831+32.0*t1031*t347*t1033+8.0*t677*t1290*t941+8.0*t2366*t2367-24.0*t152*t52*t1126*t813-32.0*t319*t810*t1054+8.0*t325*t2046*t2378*t1184+4.0*t488*t2203*t1260-24.0*t1580*t1053*t604+4.0*t1971*t1139*t813;
650.     double t2392 = t56*t235;
651.     double t2398 = t145*t3;
652.     double t2409 = t1213*t22;
653.     double t2432 = t235*t371;
654.     double t2434 = t144*t205;
655.     double t2435 = t2434*t1047;
656.     double t2442 = -32.0*t146*t2392*t138*t336*t1047-t2398*eta*t108*t115-8.0*t1797*t327*t110-4.0*t2398*eta*t1298*t58+4.0*t17*t2409+2.0*t228*t1053*t32+8.0*t2366*t938+8.0*t17*t2378*t145*t56*t111+16.0*t509*t2378*t1195+8.0*t529*t2378*t1214+16.0*t1800*t43*t56*t831+16.0*t2238*t2432*t2435+4.0*t646*t251*t108*t941;
657.     double t2448 = t144*t58*t846;
658.     double t2463 = t2434*t88;
659.     double t2485 = 32.0*t1057*t2392*t1060-8.0*t2366*t2448+16.0*t2260*t2432*t2435+8.0
660. *t1071*t2367+8.0*t1071*t938+32.0*t334*t352*t2284+4.0*t43*t369*t2190*t2463-8.0*t1071*t2448+8.0*t2274*t2190*t2463+4.0*t251*t145*t2190*t2463+4.0*t1955*t851+8.0*t828*t810*t809*t98-2.0*t2180*t2046*t1192*t98;
661.     double t2506 = t180*t835*t235;
662.     double t2513 = t180*t1980;
663.     double t2520 = t146*t1980;
664.     double t2531 = -32.0*t204*t41*t334*t2284+48.0*t366*t1049+48.0*t237*t1049-32.0*t147*t2322-16.0*t250*t310*t1947+8.0*t147*t385*t1184+4.0*t250*t2226*t1184+8.0*t2506*t1522*t1144-8.0*t2506*t1522*t937+8.0*t2513*t1024*t1408-12.0*t2513*t1024*t941+8.0*t2520*t979*t1408-12.0*t2520*t979*t941-64.0*t180*t347*t1522*t1048;
665.     double t2537 = t129*(t2212+t2258+t2304+t2346+t2391+t2442+t2485+t2531)*t1338/2.0;
666.     double t2589 = 2.0*t327+2.0*t325*t159-2.0*t1797*t1342+8.0*t341*t106*t1559+4.0*t341*t106*t1567+t492*t1843*t208+4.0*t492*t182*t138+2.0*t341*t1843*t43+4.0*t341*t182*r+4.0*t514*t354*t509+12.0*t514*t726*t637-4.0*t52*t1220*t1856+t195*t762+4.0*t591*t354*t17-4.0*t2022*t1879+4.0*t591*t726*M+4.0*t596*t239*t63+8.0*t596*t592*t54-8.0*t341*t307-8.0*t596*t1805;
667.     double t2625 = 2.0*t1201*t22+2.0*t2409+8.0*t1797*t2017-2.0*t2181*t984*t1+16.0*t52*t1298*t806+4.0*t492*t1122*t299+8.0*t341*t1122*t164+24.0*t52*t1126*t674+4.0*t358*t1122*t58+24.0*t195*t1126*t726*t17+32.0*t195*t204*t110*t592*t63;
668.     double t2646 = t195*t835;
669.     double t2655 = 4.0*t1315*t1177+t1309*t1310*t208+2.0*t1315*t1311+4.0*t1318*t941*t509-4.0*t1259*t515*t873+t195*t1166+4.0*t195*t976*t1319-4.0*t195*t976*t1*t1260+
670. 4.0*t2646*t937*t63-8.0*t1315*t864-8.0*t2646*t1144*t63;
671.     christoffel[0][0][0] = 0.0;
672.     christoffel[0][0][1] = -t7*t14*t77/2.0;
673.     christoffel[0][0][2] = -t86*(-2.0*t89+t9*t92-4.0*t26*t95+4.0*t31*t99+4.0*t31*t103+16.0*t107*t112+4.0*t53*t63*t116-8.0*t9*t87*t70-8.0*t52*t87*t74)/2.0;
674.     christoffel[0][0][3] = 0.0;
675.     christoffel[0][1][0] = t633;
676.     christoffel[0][1][1] = 0.0;
677.     christoffel[0][1][2] = 0.0;
678.     christoffel[0][1][3] = t827;
679.     christoffel[0][2][0] = t1108;
680.     christoffel[0][2][1] = 0.0;
681.     christoffel[0][2][2] = 0.0;
682.     christoffel[0][2][3] = t1340;
683.     christoffel[0][3][0] = 0.0;
684.     christoffel[0][3][1] = -t1376;
685.     christoffel[0][3][2] = -t1420;
686.     christoffel[0][3][3] = 0.0;
687.     christoffel[1][0][0] = t633;
688.     christoffel[1][0][1] = 0.0;
689.     christoffel[1][0][2] = 0.0;
690.     christoffel[1][0][3] = t827;
691.     christoffel[1][1][0] = 0.0;
692.     christoffel[1][1][1] = t1421/t5*t14*(t18*t6-t1*t1425*t5+2.0*t31*t10*t18*t5-2.0*t9*t2*t10*t1425+t9*t2*t1436*t5)/2.0;
693.     christoffel[1][1][2] = -t1421*t85*t1452/t6/2.0;
694.     christoffel[1][1][3] = 0.0;
695.     christoffel[1][2][0] = 0.0;
696.     christoffel[1][2][1] = t1459;
697.     christoffel[1][2][2] = t1469;
698.     christoffel[1][2][3] = 0.0;
699.     christoffel[1][3][0] = t1712;
700.     christoffel[1][3][1] = 0.0;
701.     christoffel[1][3][2] = 0.0;
702.     christoffel[1][3][3] = t1894;
703.     christoffel[2][0][0] = t1108;
704.     christoffel[2][0][1] = 0.0;
705.     christoffel[2][0][2] = 0.0;
706.     christoffel[2][0][3] = t1340;
707.     christoffel[2][1][0] = 0.0;
708.     christoffel[2][1][1] = t1459;
709.     christoffel[2][1][2] = t1469;
710.     christoffel[2][1][3] = 0.0;
711.     christoffel[2][2][0] = 0.0;
712.     christoffel[2][2][1] = -t6*t1466*t1457/2.0;
713.     christoffel[2][2][2] = (t87+2.0*t9*t82*t87+t9*t2*t1901)*t1421*t85/2.0;
714.     christoffel[2][2][3] = 0.0;
715.     christoffel[2][3][0] = -t2168;
716.     christoffel[2][3][1] = 0.0;
717.     christoffel[2][3][2] = 0.0;
718.     christoffel[2][3][3] = t2537;
719.     christoffel[3][0][0] = 0.0;
720.     christoffel[3][0][1] = -t1376;
721.     christoffel[3][0][2] = -t1420;
722.     christoffel[3][0][3] = 0.0;
723.     christoffel[3][1][0] = t1712;
724.     christoffel[3][1][1] = 0.0;
725.     christoffel[3][1][2] = 0.0;
726.     christoffel[3][1][3] = t1894;
727.     christoffel[3][2][0] = -t2168;
728.     christoffel[3][2][1] = 0.0;
729.     christoffel[3][2][2] = 0.0;
730.     christoffel[3][2][3] = t2537;
731.     christoffel[3][3][0] = 0.0;
732.     christoffel[3][3][1] = -t7*t14*t56*t2589/2.0;
733.     christoffel[3][3][2] = -t86*t55*(t2625+t2655)/2.0;
734.     christoffel[3][3][3] = 0.0;
735.     return true;
736. }