1 /* Ergo, version 3.8, a program for linear scaling electronic structure
2 * calculations.
3 * Copyright (C) 2019 Elias Rudberg, Emanuel H. Rubensson, Pawel Salek,
4 * and Anastasia Kruchinina.
5 *
6 * This program is free software: you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation, either version 3 of the License, or
9 * (at your option) any later version.
10 *
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
15 *
16 * You should have received a copy of the GNU General Public License
17 * along with this program. If not, see <http://www.gnu.org/licenses/>.
18 *
19 * Primary academic reference:
20 * Ergo: An open-source program for linear-scaling electronic structure
21 * calculations,
22 * Elias Rudberg, Emanuel H. Rubensson, Pawel Salek, and Anastasia
23 * Kruchinina,
24 * SoftwareX 7, 107 (2018),
25 * <http://dx.doi.org/10.1016/j.softx.2018.03.005>
26 *
27 * For further information about Ergo, see <http://www.ergoscf.org>.
28 */
29
30 /*-*-mode: C; c-indentation-style: "bsd"; c-basic-offset: 4; -*-*/
31 /** @file fun-pbec.c PBEc implementation.
32
33 Automatically generated code implementing PBEC functional and
34 its derivatives. It is generated by func-codegen.pl being a part of
35 a "Automatic code generation framework for analytical functional
36 derivative evaluation", Pawel Salek, 2005
37
38 This functional is connected by making following changes:
39 1. add "extern Functional pbecFunctional;" to 'functionals.h'
40 2. add "&pbecFunctional," to 'functionals.c'
41 3. add "fun-pbec.c" to 'Makefile.am', 'Makefile.in' or 'Makefile'.
42
43 This functional has been generated from following input:
44 ------ cut here -------
45 parameter fuer PW92C
46 c:1.709921;
47 T:[0.031091,0.015545,0.016887];
48 U:[0.21370,0.20548,0.11125];
49 V:[7.5957,14.1189,10.357];
50 W:[3.5876,6.1977,3.6231];
51 X:[1.6382,3.3662,0.88026];
52 Y:[0.49294,0.62517,0.49671];
53 P:[1,1,1];
54
55 funktionen fuer PW92C
56 r(a,b):=(3/(4*%PI*(a+b)))^(1/3);
57 zet(a,b):=(a-b)/(a+b);
58 omega(z):=((1+z)^(4/3)+(1-z)^(4/3)-2)/(2^(4/3)-2);
59 en(r,t,u,v,w,x,y,p):=-2*t*(1+u*r)
60 *log(1+1/(2*t*(v*sqrt(r)+w*r+x*r^(3/2)+y*r^(p+1))));
61
62 eps(a,b) := en(r(a,b),T[1],U[1],V[1],W[1],X[1],Y[1],P[1])
63 - en(r(a,b),T[3],U[3],V[3],W[3],X[3],Y[3],P[3])
64 *omega(zet(a,b))*(1-(zet(a,b))^4)/c
65 + (en(r(a,b),T[2],U[2],V[2],W[2],X[2],Y[2],P[2])
66 -en(r(a,b),T[1],U[1],V[1],W[1],X[1],Y[1],P[1]))
67 *omega(zet(a,b))*(zet(a,b))^4;
68
69 parameter fuer pbec
70 iota:0.0716;
71 nu:16/%PI*(3*%PI^2)^(1/3);
72 kappa:0.004235;
73 lamda:nu*kappa;
74 Z:-0.001667;
75 Xi:23.266;
76 Phi:0.007389;
77 Lamda:8.723;
78 Upsilon:0.472;
79
80 funktionen fuer pbec
81 u(a,b) := ((1+zet(a,b))^(2/3)+(1-zet(a,b)^(2/3)))/2;
82
83 u(a,b) := 0.5*((2*a)^(2/3)+(2*b)^(2/3))/(a+b)^(2/3);
84 theta(r):= 10^(-3)*(2.568+Xi*r+Phi*r^2)
85 /(1+Lamda*r+Upsilon*r^2+10*Phi*r^3);
86 phi(r) := theta(r)-Z;
87 sigma(ga,gb,gab) := ga*ga+gb*gb+2*gab;
88 d(a,b,ga,gb,gab) := sqrt(sigma(ga,gb,gab))
89 /(4*u(a,b)*(3/%PI)^(1/6)*(a+b)^(7/6));
90 A(a,b) := 2*iota/(lamda*(exp(-2*iota*eps(a,b)/(lamda^2*(u(a,b))^3))-1));
91 L(d,a,b) := (u(a,b))^3*lamda^2/(2*iota)
92 *log((1+2*(iota*(d^2+A(a,b)*d^4))
93 /(lamda*(1+A(a,b)*d^2+(A(a,b))^2*d^4))));
94
95 kernel
96 K(rhoa,grada,rhob,gradb,gradab):=
97 (rhoa+rhob)*(eps(rhoa,rhob)+
98 L(d(rhoa,rhob,grada,gradb,gradab),rhoa,rhob));
99
100 ------ cut here -------
101 */
102
103 #include <math.h>
104 #include <stddef.h>
105
106 #define __CVERSION__
107
108 #include "functionals.h"
109
110 /* INTERFACE PART */
pbec_isgga(void)111 static int pbec_isgga(void) { return 1; } /* FIXME: detect! */
112 static int pbec_read(const char *conf_line);
113 static real pbec_energy(const FunDensProp* dp);
114 static void pbec_first(FunFirstFuncDrv *ds, real factor,
115 const FunDensProp* dp);
116 static void pbec_second(FunSecondFuncDrv *ds, real factor,
117 const FunDensProp* dp);
118
119 Functional PbecFunctional = {
120 "PBEC", /* name */
121 pbec_isgga, /* gga-corrected */
122 pbec_read,
123 NULL,
124 pbec_energy,
125 pbec_first,
126 pbec_second /*,
127 pbec_third,
128 pbec_fourth */
129 };
130
131 /* IMPLEMENTATION PART */
132 static int
pbec_read(const char * conf_line)133 pbec_read(const char *conf_line)
134 {
135 fun_set_hf_weight(0);
136 return 1;
137 }
138
139 static real
pbec_energy(const FunDensProp * dp)140 pbec_energy(const FunDensProp *dp)
141 {
142 real res;
143 real rhoa = dp->rhoa, rhob = dp->rhob;
144 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
145
146 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
147 real t11, t12, t13, t14, t15, t16, t17, t18;
148 real t19, t20, t21, t22, t23, t24, t25, t26;
149 real t27, t28, t29, t30, t31, t32, t33, t34;
150 real t35, t36, t37, t38, t39, t40, t41, t42;
151 real t43;
152
153 t1 = rhob+rhoa;
154 t2 = POW(3.0,.3333333333333333);
155 t3 = 1/POW(4.0,.3333333333333333);
156 t4 = POW(3.141592653589793,.3333333333333333);
157 t5 = 1/t4;
158 t6 = 1/POW(t1,.3333333333333333);
159 t7 = 0.2137*t2*t3*t5*t6+1.0;
160 t8 = POW(3.0,.6666666666666666);
161 t9 = 1/POW(4.0,.6666666666666666);
162 t10 = POW(3.141592653589793,.6666666666666666);
163 t11 = 1/t10;
164 t12 = 1/POW(t1,.6666666666666666);
165 t13 = 1.732050807568877;
166 t14 = .5641895835477563;
167 t15 = 1/SQRT(t1);
168 t16 = POW(2.0,.3333333333333333);
169 t17 = 1/t16;
170 t18 = POW(3.0,.1666666666666667);
171 t19 = 1/POW(3.141592653589793,.1666666666666667);
172 t20 = 1/POW(t1,.1666666666666667);
173 t21 = LOG(16.0818243221511/(7.5957*t17*t18*t19*t20+3.5876*
174 t2*t3*t5*t6+0.8191*t13*t14*t15+0.49294*t8*t9*t11*t12)+1.0);
175 t22 = -
176 0.062182*t7*t21;
177 t23 = 1/(2.0*t16-2.0);
178 t24 = rhoa-1.0*rhob;
179 t25 = POW(t24,4.0);
180 t26 = 1/POW(t1,4.0);
181 t27 = 1/t1;
182 t28 = POW(t24*t27+1.0,1.333333333333333)+POW(1.0-1.0*
183 t24*t27,1.333333333333333)-2.0;
184 t29 = .0197517897025652*LOG(29.60857464321668/(10.357*
185 t17*t18*t19*t20+3.6231*t2*t3*t5*t6+0.44013*t13*t14*t15+0.49671*
186 t8*t9*t11*t12)+1.0)*t23*(1.0-1.0*t25*t26)*t28*(0.11125*t2*
187 t3*t5*t6+1.0);
188 t30 = t23*t25*t26*t28*(0.062182*t7*t21-0.03109*LOG(32.1646831778707/
189 (14.1189*t17*t18*t19*t20+6.1977*t2*t3*t5*t6+1.6831*t13*t14*
190 t15+0.62517*t8*t9*t11*t12)+1.0)*(0.20548*t2*t3*t5*t6+1.0));
191 t31 = POW(2.0,
192 .6666666666666666);
193 t32 = t31*POW(rhob,.6666666666666666)+t31*POW(rhoa,.6666666666666666);
194 t33 = POW(t32,
195 3.0);
196 t34 = POW(t1,2.0);
197 t35 = 1/t34;
198 t36 = 1/t2;
199 t37 = POW(gradb,2.0)+2.0*gradab+POW(grada,2.0);
200 t38 = 1/POW(t32,2.0);
201 t39 = POW(t37,2.0);
202 t40 = 1/POW(t32,4.0);
203 t41 = 1/t8;
204 t42 = 1/POW(2.718281828459045,249.5089969598932*t10*(t30+
205 t29+t22)*t34*t41/t33)-1.0;
206 t43 = 1/t42;
207
208 /* code */
209 res = t1*(.0040078715083799*t11*t33*t35*LOG(2.113341204250295*
210 t36*t4*(.1383178583708611*t35*t39*t40*t43+0.25*t36*t4*t37*
211 t38*t27)/(.5283353010625738*t41*t10*t37*t38*t27*t43+.0930460634496268*
212 POW(3.141592653589793,1.333333333333333)*t35*t36*t39*t40/POW(t42,
213 2.0)+1.0)+1.0)*t8+t30+t29+t22);
214
215 return res;
216 }
217
218 static void
pbec_first(FunFirstFuncDrv * ds,real factor,const FunDensProp * dp)219 pbec_first(FunFirstFuncDrv *ds, real factor, const FunDensProp *dp)
220 {
221 real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
222 real rhoa = dp->rhoa, rhob = dp->rhob;
223 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
224
225 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
226 real t11, t12, t13, t14, t15, t16, t17, t18;
227 real t19, t20, t21, t22, t23, t24, t25, t26;
228 real t27, t28, t29, t30, t31, t32, t33, t34;
229 real t35, t36, t37, t38, t39, t40, t41, t42;
230 real t43, t44, t45, t46, t47, t48, t49, t50;
231 real t51, t52, t53, t54, t55, t56, t57, t58;
232 real t59, t60, t61, t62, t63, t64, t65, t66;
233 real t67, t68, t69, t70, t71, t72, t73, t74;
234 real t75, t76, t77, t78, t79, t80, t81, t82;
235 real t83, t84, t85, t86, t87, t88, t89, t90;
236 real t91, t92, t93, t94, t95, t96, t97, t98;
237 real t99, t100, t101, t102, t103, t104, t105;
238 real t106, t107, t108, t109, t110, t111, t112;
239 real t113, t114, t115, t116;
240
241 t1 = POW(3.0,.3333333333333333);
242 t2 = 1/POW(4.0,.3333333333333333);
243 t3 = POW(3.141592653589793,.3333333333333333);
244 t4 = 1/t3;
245 t5 = rhob+rhoa;
246 t6 = 1/POW(t5,.3333333333333333);
247 t7 = 0.2137*t1*t2*t4*t6+1.0;
248 t8 = POW(3.0,.6666666666666666);
249 t9 = 1/POW(4.0,.6666666666666666);
250 t10 = POW(3.141592653589793,.6666666666666666);
251 t11 = 1/t10;
252 t12 = 1/POW(t5,.6666666666666666);
253 t13 = 1.732050807568877;
254 t14 = .5641895835477563;
255 t15 = SQRT(t5);
256 t16 = 1/t15;
257 t17 = POW(2.0,.3333333333333333);
258 t18 = 1/t17;
259 t19 = POW(3.0,.1666666666666667);
260 t20 = 1/POW(3.141592653589793,.1666666666666667);
261 t21 = 1/POW(t5,.1666666666666667);
262 t22 = 7.5957*t18*t19*t20*t21+3.5876*t1*t2*t4*t6+0.8191*
263 t13*t14*t16+0.49294*t8*t9*t11*t12;
264 t23 = 16.0818243221511/t22+1.0;
265 t24 = LOG(t23);
266 t25 = -0.062182*t7*t24;
267 t26 = 1/(2.0*t17-2.0);
268 t27 = rhoa-1.0*rhob;
269 t28 = POW(t27,4.0);
270 t29 = 1/POW(t5,4.0);
271 t30 = 1.0-1.0*t28*t29;
272 t31 = 0.11125*t1*t2*t4*t6+1.0;
273 t32 = 1/t5;
274 t33 = 1.0-1.0*t27*t32;
275 t34 = t27*t32+1.0;
276 t35 = POW(t34,1.333333333333333)+POW(t33,1.333333333333333)-
277 2.0;
278 t36 = 10.357*t18*t19*t20*t21+3.6231*t1*t2*t4*t6+0.44013*
279 t13*t14*t16+0.49671*t8*t9*t11*t12;
280 t37 = 29.60857464321668/t36+1.0;
281 t38 = LOG(t37);
282 t39 = .0197517897025652*t26*t30*t31*t35*t38;
283 t40 = 0.20548*t1*t2*t4*t6+1.0;
284 t41 = 14.1189*t18*t19*t20*t21+6.1977*t1*t2*t4*t6+1.6831*
285 t13*t14*t16+0.62517*t8*t9*t11*t12;
286 t42 = 32.1646831778707/t41+1.0;
287 t43 = LOG(t42);
288 t44 = 0.062182*t7*t24-0.03109*t40*t43;
289 t45 = t26*t28*t29*t35*t44;
290 t46 = POW(2.0,.6666666666666666);
291 t47 = t46*POW(rhob,.6666666666666666)+t46*POW(rhoa,.6666666666666666);
292 t48 = POW(t47,
293 3.0);
294 t49 = POW(t5,2.0);
295 t50 = 1/t49;
296 t51 = 1/t1;
297 t52 = POW(gradb,2.0)+2.0*gradab+POW(grada,2.0);
298 t53 = POW(t47,2.0);
299 t54 = 1/t53;
300 t55 = POW(t52,2.0);
301 t56 = 1/POW(t47,4.0);
302 t57 = 1/t8;
303 t58 = 1/t48;
304 t59 = t45+t39+t25;
305 t60 = 1/POW(2.718281828459045,249.5089969598932*t57*t10*
306 t58*t49*t59);
307 t61 = t60-1.0;
308 t62 = 1/t61;
309 t63 = .1383178583708611*t50*t55*t56*t62+0.25*t51*t3*t52*
310 t54*t32;
311 t64 = POW(3.141592653589793,1.333333333333333);
312 t65 = 1/POW(t61,2.0);
313 t66 = .0930460634496268*t51*t64*t55*t56*t50*t65+.5283353010625738*
314 t57*t10*t52*t54*t32*t62+1.0;
315 t67 = 1/t66;
316 t68 = 2.113341204250295*t51*t3*t63*t67+1.0;
317 t69 = LOG(t68);
318 t70 = .0040078715083799*t8*t11*t48*t50*t69;
319 t71 = 1/POW(t5,1.666666666666667);
320 t72 = 1/POW(t15,3.0);
321 t73 = 1/POW(t5,1.333333333333333);
322 t74 = 1/POW(t5,1.166666666666667);
323 t75 = -1.26595*t18*t19*t20*t74-1.195866666666667*t1*t2*
324 t4*t73-0.40955*t13*t14*t72-.3286266666666666*t8*t9*t11*t71;
325 t76 = 1/
326 POW(t22,2.0);
327 t77 = 1/t23;
328 t78 = .9999999999999999*t75*t7*t76*t77;
329 t79 = -.5848223397455204*t26*t30*t31*t35*(-1.726166666666666*
330 t18*t19*t20*t74-1.2077*t1*t2*t4*t73-0.220065*t13*t14*t72-0.33114*
331 t8*t9*t11*t71)/(POW(t36,2.0)*t37);
332 t80 = .0044294311333333*t1*t2*t4*t73*t24;
333 t81 = t27*t50;
334 t82 = -1.0*t32;
335 t83 = POW(t33,.3333333333333333);
336 t84 = -1.0*t27*t50;
337 t85 = POW(t34,.3333333333333333);
338 t86 = 1.333333333333333*(t32+t84)*t85+1.333333333333333*
339 (t82+t81)*t83;
340 t87 = .0197517897025652*t26*t30*t31*t86*t38;
341 t88 = -7.324622014701263e-4*t26*t1*t2*t4*t73*t30*t35*
342 t38;
343 t89 = 1/POW(t5,5.0);
344 t90 = 4.0*t28*t89;
345 t91 = POW(t27,3.0);
346 t92 = .0197517897025652*t26*t31*t35*t38*(t90-4.0*t29*
347 t91);
348 t93 = t26*t28*t29*t35*(-.9999999999999999*t75*t7*t76*
349 t77+t40*(-2.35315*t18*t19*t20*t74-2.0659*t1*t2*t4*t73-0.84155*
350 t13*t14*t72-0.41678*t8*t9*t11*t71)/(POW(t41,2.0)*t42)+.0021294577333333*
351 t1*t2*t4*t73*t43-.0044294311333333*t1*t2*t4*t73*t24);
352 t94 = t26*t28*t29*t86*t44;
353 t95 = -4.0*t26*t28*t35*t44*t89;
354 t96 = 4.0*t26*t29*t35*t44*t91;
355 t97 = 1/t68;
356 t98 = -0.25*t51*t3*t52*t54*t50;
357 t99 = 1/POW(rhoa,.3333333333333333);
358 t100 = 1/POW(t5,3.0);
359 t101 = -.2766357167417222*t100*t55*t56*t62;
360 t102 = 1/POW(t47,5.0);
361 t103 = -499.0179939197864*t57*t10*t58*t5*t59;
362 t104 = 499.0179939197864*t46*t57*t10*t99*t56*t49*t59+
363 t103-249.5089969598932*t57*t10*t58*t49*(t96+t95+t94+t93+t92+
364 t88+t87+t80+t79+t78);
365 t105 = 1/POW(t66,2.0);
366 t106 = -.1860921268992536*t51*t64*t55*t56*t100*t65;
367 t107 = -
368 .5283353010625738*t57*t10*t52*t54*t50*t62;
369 t108 = 1/POW(t61,3.0);
370 t109 = -.0080157430167598*t8*t11*t48*t100*t69;
371 t110 = 1.333333333333333*(t82+t84)*t85+1.333333333333333*
372 (t32+t81)*t83;
373 t111 = .0197517897025652*t26*t30*t31*t110*t38;
374 t112 = .0197517897025652*t26*t31*t35*t38*(4.0*t29*t91+
375 t90);
376 t113 = t26*t28*t29*t110*t44;
377 t114 = -4.0*t26*t29*t35*t44*t91;
378 t115 = 1/POW(rhob,.3333333333333333);
379 t116 = 499.0179939197864*t46*t57*t10*t56*t115*t49*t59+
380 t103-249.5089969598932*t57*t10*t58*t49*(t114+t95+t113+t93+
381 t112+t88+t111+t80+t79+t78);
382
383 /* code */
384 dfdra = t5*(.0040078715083799*t11*t48*t50*t8*t97*(2.113341204250295*
385 t3*t51*t67*(-.3688476223222962*t102*t46*t50*t55*t62*t99+t98-
386 .1383178583708611*t104*t50*t55*t56*t60*t65-.3333333333333333*
387 t46*t51*t3*t52*t99*t58*t32+t101)-2.113341204250295*t51*t3*
388 t63*t105*(-.5283353010625738*t57*t10*t52*t54*t32*t104*t65*
389 t60-.1860921268992536*t51*t64*t55*t56*t50*t104*t108*t60-.7044470680834316*
390 t46*t57*t10*t52*t99*t58*t32*t62+t107-.2481228358656715*t46*
391 t51*t64*t55*t99*t102*t50*t65+t106))+t96+t95+t94+t93+t92+t88+
392 t87+t80+t79+t78+.0080157430167598*t46*t8*t11*t99*t53*t50*t69+
393 t109)+t70+t45+t39+t25;
394 dfdrb = t5*(.0040078715083799*t11*t48*t50*t8*t97*(2.113341204250295*
395 t3*t51*t67*(t98-.1383178583708611*t116*t50*t55*t56*t60*t65-
396 .3688476223222962*t102*t115*t46*t50*t55*t62-.3333333333333333*
397 t46*t51*t3*t52*t58*t115*t32+t101)-2.113341204250295*t51*t3*
398 t63*t105*(-.5283353010625738*t57*t10*t52*t54*t32*t116*t65*
399 t60-.1860921268992536*t51*t64*t55*t56*t50*t116*t108*t60-.7044470680834316*
400 t46*t57*t10*t52*t58*t115*t32*t62+t107-.2481228358656715*t46*
401 t51*t64*t55*t102*t115*t50*t65+t106))+t95+t93+t88+t80+t79+t78+
402 .0080157430167598*t46*t8*t11*t53*t115*t50*t69+t114+t113+t112+
403 t111+t109)+t70+t45+t39+t25;
404 dfdga = .0040078715083799*t11*t32*t48*t8*t97*(2.113341204250295*
405 t3*t51*t67*(.5532714334834444*t50*t52*t56*t62*grada+0.5*t51*
406 t3*grada*t54*t32)-2.113341204250295*t51*t3*t63*(1.056670602125148*
407 t57*t10*grada*t54*t32*t62+.3721842537985072*t51*t64*grada*
408 t52*t56*t50*t65)*t105);
409 dfdgb = .0040078715083799*t11*t32*t48*t8*t97*(2.113341204250295*
410 t3*t51*t67*(.5532714334834444*t50*t52*t56*t62*gradb+0.5*t51*
411 t3*gradb*t54*t32)-2.113341204250295*t51*t3*t63*(1.056670602125148*
412 t57*t10*gradb*t54*t32*t62+.3721842537985072*t51*t64*gradb*
413 t52*t56*t50*t65)*t105);
414 dfdgab = .0040078715083799*t11*t32*t48*(2.113341204250295*
415 t3*t51*(.5532714334834444*t50*t52*t56*t62+0.5*t51*t3*t54*t32)*
416 t67-2.113341204250295*t51*t3*t63*(1.056670602125148*t57*t10*
417 t54*t32*t62+.3721842537985072*t51*t64*t52*t56*t50*t65)*t105)*
418 t8*t97;
419
420
421 ds->df1000 += factor*dfdra;
422 ds->df0100 += factor*dfdrb;
423 ds->df0010 += factor*dfdga;
424 ds->df0001 += factor*dfdgb;
425 ds->df00001+= factor*dfdgab;
426
427 }
428
429 static void
pbec_second(FunSecondFuncDrv * ds,real factor,const FunDensProp * dp)430 pbec_second(FunSecondFuncDrv *ds, real factor, const FunDensProp* dp)
431 {
432 real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
433 real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
434 d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
435 d2fdgbgb, d2fdgbgab, d2fdgabgab;
436 real rhoa = dp->rhoa, rhob = dp->rhob;
437 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
438
439 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
440 real t11, t12, t13, t14, t15, t16, t17, t18;
441 real t19, t20, t21, t22, t23, t24, t25, t26;
442 real t27, t28, t29, t30, t31, t32, t33, t34;
443 real t35, t36, t37, t38, t39, t40, t41, t42;
444 real t43, t44, t45, t46, t47, t48, t49, t50;
445 real t51, t52, t53, t54, t55, t56, t57, t58;
446 real t59, t60, t61, t62, t63, t64, t65, t66;
447 real t67, t68, t69, t70, t71, t72, t73, t74;
448 real t75, t76, t77, t78, t79, t80, t81, t82;
449 real t83, t84, t85, t86, t87, t88, t89, t90;
450 real t91, t92, t93, t94, t95, t96, t97, t98;
451 real t99, t100, t101, t102, t103, t104, t105;
452 real t106, t107, t108, t109, t110, t111, t112;
453 real t113, t114, t115, t116, t117, t118, t119;
454 real t120, t121, t122, t123, t124, t125, t126;
455 real t127, t128, t129, t130, t131, t132, t133;
456 real t134, t135, t136, t137, t138, t139, t140;
457 real t141, t142, t143, t144, t145, t146, t147;
458 real t148, t149, t150, t151, t152, t153, t154;
459 real t155, t156, t157, t158, t159, t160, t161;
460 real t162, t163, t164, t165, t166, t167, t168;
461 real t169, t170, t171, t172, t173, t174, t175;
462 real t176, t177, t178, t179, t180, t181, t182;
463 real t183, t184, t185, t186, t187, t188, t189;
464 real t190, t191, t192, t193, t194, t195, t196;
465 real t197, t198, t199, t200, t201, t202, t203;
466 real t204, t205, t206, t207, t208, t209, t210;
467 real t211, t212, t213, t214, t215, t216, t217;
468 real t218, t219, t220, t221, t222, t223, t224;
469 real t225, t226, t227, t228, t229, t230, t231;
470 real t232, t233, t234, t235, t236, t237, t238;
471 real t239, t240, t241, t242, t243, t244, t245;
472 real t246, t247, t248, t249, t250, t251, t252;
473 real t253, t254, t255, t256, t257, t258, t259;
474 real t260, t261, t262, t263, t264, t265, t266;
475 real t267, t268, t269, t270, t271, t272, t273;
476 real t274, t275, t276, t277, t278, t279, t280;
477 real t281, t282, t283, t284, t285, t286, t287;
478 real t288;
479
480 t1 = POW(3.0,.3333333333333333);
481 t2 = 1/POW(4.0,.3333333333333333);
482 t3 = POW(3.141592653589793,.3333333333333333);
483 t4 = 1/t3;
484 t5 = rhob+rhoa;
485 t6 = 1/POW(t5,.3333333333333333);
486 t7 = 0.2137*t1*t2*t4*t6+1.0;
487 t8 = POW(3.0,.6666666666666666);
488 t9 = 1/POW(4.0,.6666666666666666);
489 t10 = POW(3.141592653589793,.6666666666666666);
490 t11 = 1/t10;
491 t12 = 1/POW(t5,.6666666666666666);
492 t13 = 1.732050807568877;
493 t14 = .5641895835477563;
494 t15 = SQRT(t5);
495 t16 = 1/t15;
496 t17 = POW(2.0,.3333333333333333);
497 t18 = 1/t17;
498 t19 = POW(3.0,.1666666666666667);
499 t20 = 1/POW(3.141592653589793,.1666666666666667);
500 t21 = 1/POW(t5,.1666666666666667);
501 t22 = 7.5957*t18*t19*t20*t21+3.5876*t1*t2*t4*t6+0.8191*
502 t13*t14*t16+0.49294*t8*t9*t11*t12;
503 t23 = 16.0818243221511/t22+1.0;
504 t24 = LOG(t23);
505 t25 = -0.062182*t7*t24;
506 t26 = 1/(2.0*t17-2.0);
507 t27 = rhoa-1.0*rhob;
508 t28 = POW(t27,4.0);
509 t29 = 1/POW(t5,4.0);
510 t30 = 1.0-1.0*t28*t29;
511 t31 = 0.11125*t1*t2*t4*t6+1.0;
512 t32 = 1/t5;
513 t33 = 1.0-1.0*t27*t32;
514 t34 = t27*t32+1.0;
515 t35 = POW(t34,1.333333333333333)+POW(t33,1.333333333333333)-
516 2.0;
517 t36 = 10.357*t18*t19*t20*t21+3.6231*t1*t2*t4*t6+0.44013*
518 t13*t14*t16+0.49671*t8*t9*t11*t12;
519 t37 = 29.60857464321668/t36+1.0;
520 t38 = LOG(t37);
521 t39 = .0197517897025652*t26*t30*t31*t35*t38;
522 t40 = 0.20548*t1*t2*t4*t6+1.0;
523 t41 = 14.1189*t18*t19*t20*t21+6.1977*t1*t2*t4*t6+1.6831*
524 t13*t14*t16+0.62517*t8*t9*t11*t12;
525 t42 = 32.1646831778707/t41+1.0;
526 t43 = LOG(t42);
527 t44 = 0.062182*t7*t24-0.03109*t40*t43;
528 t45 = t26*t28*t29*t35*t44;
529 t46 = POW(2.0,.6666666666666666);
530 t47 = POW(rhoa,.6666666666666666);
531 t48 = POW(rhob,.6666666666666666);
532 t49 = t46*t48+t46*t47;
533 t50 = POW(t49,3.0);
534 t51 = POW(t5,2.0);
535 t52 = 1/t51;
536 t53 = 1/t1;
537 t54 = POW(grada,2.0);
538 t55 = POW(gradb,2.0);
539 t56 = 2.0*gradab+t55+t54;
540 t57 = POW(t49,2.0);
541 t58 = 1/t57;
542 t59 = POW(t56,2.0);
543 t60 = 1/POW(t49,4.0);
544 t61 = 1/t8;
545 t62 = 1/t50;
546 t63 = t45+t39+t25;
547 t64 = 1/POW(2.718281828459045,249.5089969598932*t61*t10*
548 t62*t51*t63);
549 t65 = t64-1.0;
550 t66 = 1/t65;
551 t67 = .1383178583708611*t52*t59*t60*t66+0.25*t53*t3*t56*
552 t58*t32;
553 t68 = POW(3.141592653589793,1.333333333333333);
554 t69 = 1/POW(t65,2.0);
555 t70 = .0930460634496268*t53*t68*t59*t60*t52*t69+.5283353010625738*
556 t61*t10*t56*t58*t32*t66+1.0;
557 t71 = 1/t70;
558 t72 = 2.113341204250295*t53*t3*t67*t71+1.0;
559 t73 = LOG(t72);
560 t74 = .0040078715083799*t8*t11*t50*t52*t73;
561 t75 = 1/POW(t5,1.666666666666667);
562 t76 = 1/POW(t15,3.0);
563 t77 = 1/POW(t5,1.333333333333333);
564 t78 = 1/POW(t5,1.166666666666667);
565 t79 = -1.26595*t18*t19*t20*t78-1.195866666666667*t1*t2*
566 t4*t77-0.40955*t13*t14*t76-.3286266666666666*t8*t9*t11*t75;
567 t80 = 1/
568 POW(t22,2.0);
569 t81 = 1/t23;
570 t82 = .9999999999999999*t79*t7*t80*t81;
571 t83 = -1.726166666666666*t18*t19*t20*t78-1.2077*t1*t2*
572 t4*t77-0.220065*t13*t14*t76-0.33114*t8*t9*t11*t75;
573 t84 = 1/POW(t36,2.0);
574 t85 = 1/t37;
575 t86 = -.5848223397455204*t26*t30*t83*t31*t84*t35*t85;
576 t87 = .0044294311333333*
577 t1*t2*t4*t77*t24;
578 t88 = t27*t52;
579 t89 = -1.0*t32;
580 t90 = t89+t88;
581 t91 = POW(t33,.3333333333333333);
582 t92 = -1.0*t27*t52;
583 t93 = t32+t92;
584 t94 = POW(t34,.3333333333333333);
585 t95 = 1.333333333333333*t93*t94+1.333333333333333*t90*
586 t91;
587 t96 = .0197517897025652*t26*t30*t31*t95*t38;
588 t97 = -7.324622014701263e-4*t26*t1*t2*t4*t77*t30*t35*
589 t38;
590 t98 = 1/POW(t5,5.0);
591 t99 = 4.0*t28*t98;
592 t100 = POW(t27,3.0);
593 t101 = t99-4.0*t100*t29;
594 t102 = .0197517897025652*t26*t101*t31*t35*t38;
595 t103 = -2.35315*t18*t19*t20*t78-2.0659*t1*t2*t4*t77-0.84155*
596 t13*t14*t76-0.41678*t8*t9*t11*t75;
597 t104 = 1/POW(t41,2.0);
598 t105 = 1/t42;
599 t106 = .0021294577333333*t1*t2*t4*t77*t43-.0044294311333333*
600 t1*t2*t4*t77*t24+t103*t40*t104*t105-.9999999999999999*t79*
601 t7*t80*t81;
602 t107 = t26*t28*t29*t35*t106;
603 t108 = t26*t28*t29*t95*t44;
604 t109 = -4.0*t26*t28*t35*t44*t98;
605 t110 = 4.0*t100*t26*t29*t35*t44;
606 t111 = 1/t72;
607 t112 = -0.25*t53*t3*t56*t58*t52;
608 t113 = 1/POW(rhoa,.3333333333333333);
609 t114 = 1/POW(t5,3.0);
610 t115 = -.2766357167417222*t114*t59*t60*t66;
611 t116 = 1/POW(t49,5.0);
612 t117 = t110+t109+t108+t107+t102+t97+t96+t87+t86+t82;
613 t118 = -
614 499.0179939197864*t61*t10*t62*t5*t63;
615 t119 = 499.0179939197864*t46*t61*t10*t113*t60*t51*t63+
616 t118-249.5089969598932*t61*t10*t62*t51*t117;
617 t120 = -.1383178583708611*t119*t52*t59*t60*t64*t69-.3688476223222962*
618 t113*t116*t46*t52*t59*t66-.3333333333333333*t46*t53*t3*t56*
619 t113*t62*t32+t115+t112;
620 t121 = 1/POW(t70,2.0);
621 t122 = -.1860921268992536*t53*t68*t59*t60*t114*t69;
622 t123 = -
623 .5283353010625738*t61*t10*t56*t58*t52*t66;
624 t124 = 1/POW(t65,3.0);
625 t125 = -.5283353010625738*t61*t10*t56*t58*t32*t119*t69*
626 t64-.1860921268992536*t53*t68*t59*t60*t52*t119*t124*t64-.7044470680834316*
627 t46*t61*t10*t56*t113*t62*t32*t66+t123-.2481228358656715*t46*
628 t53*t68*t59*t113*t116*t52*t69+t122;
629 t126 = 2.113341204250295*t53*t3*t71*t120-2.113341204250295*
630 t53*t3*t67*t121*t125;
631 t127 = .0040078715083799*t8*t11*t50*t52*t111*t126;
632 t128 = -.0080157430167598*t8*t11*t50*t114*t73;
633 t129 = .0080157430167598*t46*t8*t11*t113*t57*t52*t73;
634 t130 = t32+
635 t88;
636 t131 = t89+t92;
637 t132 = 1.333333333333333*t131*t94+1.333333333333333*t130*
638 t91;
639 t133 = .0197517897025652*t26*t30*t31*t132*t38;
640 t134 = t99+4.0*t100*t29;
641 t135 = .0197517897025652*t26*t134*t31*t35*t38;
642 t136 = t26*t28*t29*t132*t44;
643 t137 = -4.0*t100*t26*t29*t35*t44;
644 t138 = 1/POW(rhob,.3333333333333333);
645 t139 = t137+t109+t136+t107+t135+t97+t133+t87+t86+t82;
646 t140 = 499.0179939197864*
647 t46*t61*t10*t60*t138*t51*t63+t118-249.5089969598932*t61*t10*
648 t62*t51*t139;
649 t141 = -.1383178583708611*t140*t52*t59*t60*t64*t69-.3688476223222962*
650 t116*t138*t46*t52*t59*t66-.3333333333333333*t46*t53*t3*t56*
651 t62*t138*t32+t115+t112;
652 t142 = -.5283353010625738*t61*t10*t56*t58*t32*t140*t69*
653 t64-.1860921268992536*t53*t68*t59*t60*t52*t140*t124*t64-.7044470680834316*
654 t46*t61*t10*t56*t62*t138*t32*t66+t123-.2481228358656715*t46*
655 t53*t68*t59*t116*t138*t52*t69+t122;
656 t143 = 2.113341204250295*t53*t3*t71*t141-2.113341204250295*
657 t53*t3*t67*t121*t142;
658 t144 = .0040078715083799*t8*t11*t50*t52*t111*t143;
659 t145 = .0080157430167598*t46*t8*t11*t57*t138*t52*t73;
660 t146 = 1.056670602125148*
661 t61*t10*grada*t58*t32*t66+.3721842537985072*t53*t68*grada*
662 t56*t60*t52*t69;
663 t147 = .5532714334834444*t52*t56*t60*t66*grada+0.5*t53*
664 t3*grada*t58*t32;
665 t148 = 2.113341204250295*t53*t3*t147*t71-2.113341204250295*
666 t53*t3*t67*t146*t121;
667 t149 = 1.056670602125148*t61*t10*gradb*t58*t32*t66+.3721842537985072*
668 t53*t68*gradb*t56*t60*t52*t69;
669 t150 = .5532714334834444*t52*t56*t60*t66*gradb+0.5*t53*
670 t3*gradb*t58*t32;
671 t151 = 2.113341204250295*t53*t3*t150*t71-2.113341204250295*
672 t53*t3*t67*t149*t121;
673 t152 = .3721842537985072*t53*t68*t56*t60*t52*t69;
674 t153 = 1.056670602125148*t61*t10*t58*t32*t66;
675 t154 = t153+t152;
676 t155 = 0.5*t53*t3*t58*t32;
677 t156 = .5532714334834444*t52*t56*t60*t66;
678 t157 = t156+t155;
679 t158 = 2.113341204250295*t53*t3*t157*t71-2.113341204250295*
680 t53*t3*t67*t154*t121;
681 t159 = 2.0*t79*t7*t80*t81;
682 t160 = -1.169644679491041*t26*t30*t83*t31*t84*t35*t85;
683 t161 = .0088588622666667*
684 t1*t2*t4*t77*t24;
685 t162 = -.0014649244029403*t26*t1*t2*t4*t77*t30*t35*t38;
686 t163 = 2.0*
687 t106*t26*t28*t29*t35;
688 t164 = -8.0*t26*t28*t35*t44*t98;
689 t165 = -.0160314860335196*t8*t11*t50*t114*t73;
690 t166 = POW(t79,2.0);
691 t167 = 1/POW(t22,4.0);
692 t168 = 1/POW(t23,2.0);
693 t169 = 16.0818243221511*t166*t7*t167*t168;
694 t170 = 1/POW(t22,3.0);
695 t171 = -2.0*t166*t7*t170*t81;
696 t172 = -.1424666666666666*t1*t2*t4*t77*t79*t80*t81;
697 t173 = 1/
698 POW(t5,2.666666666666667);
699 t174 = 1/POW(t15,5.0);
700 t175 = 1/POW(t5,2.333333333333333);
701 t176 = 1/POW(t5,2.166666666666667);
702 t177 = 1.476941666666667*t18*t19*t20*t176+1.594488888888889*
703 t1*t2*t4*t175+0.614325*t13*t14*t174+0.547711111111111*t8*t9*
704 t11*t173;
705 t178 = .9999999999999999*t177*t7*t80*t81;
706 t179 = POW(t83,2.0);
707 t180 = -17.31575589937587*t179*t26*t30*t31*t35/(POW(t36,
708 4.0)*POW(t37,2.0));
709 t181 = -1.169644679491041*t26*t30*t83*t31*t84*t95*t85;
710 t182 = 1.169644679491041*
711 t179*t26*t30*t31*t35*t85/POW(t36,3.0);
712 t183 = .0433743235311261*t26*t1*t2*t4*t77*t30*t83*t84*
713 t35*t85;
714 t184 = -.5848223397455204*t26*t30*(2.013861111111111*
715 t18*t19*t20*t176+1.610266666666667*t1*t2*t4*t175+0.3300975*
716 t13*t14*t174+.5519000000000001*t8*t9*t11*t173)*t31*t84*t35*
717 t85;
718 t185 = -1.169644679491041*t26*t101*t83*t31*t84*t35*t85;
719 t186 = -
720 .0059059081777778*t1*t2*t4*t175*t24;
721 t187 = 1/POW(t33,.6666666666666666);
722 t188 = -2.0*t114*t27;
723 t189 = 2.0*t52;
724 t190 = 1/POW(t34,.6666666666666666);
725 t191 = 2.0*t114*t27;
726 t192 = -2.0*t52;
727 t193 = 1.333333333333333*(t192+t191)*t94+.4444444444444444*
728 t190*POW(t93,2.0)+1.333333333333333*(t189+t188)*t91+.4444444444444444*
729 t187*POW(t90,2.0);
730 t194 = .0197517897025652*t26*t30*t31*t193*t38;
731 t195 = -.0014649244029403*t26*t1*t2*t4*t77*t30*t95*t38;
732 t196 = .0395035794051304*
733 t26*t101*t31*t95*t38;
734 t197 = -.0014649244029403*t26*t1*t2*t4*t77*t101*t35*t38;
735 t198 = 9.766162686268349e-4*t26*t1*t2*t4*t175*t30*t35*
736 t38;
737 t199 = 1/POW(t5,6.0);
738 t200 = -20.0*t199*t28;
739 t201 = POW(t27,2.0);
740 t202 = -12.0*t201*t29;
741 t203 = .0197517897025652*t26*t31*t35*t38*(32.0*t100*t98+
742 t202+t200);
743 t204 = POW(t103,2.0);
744 t205 = t26*t28*t29*t35*(.1424666666666666*t1*t2*t4*t77*
745 t79*t80*t81-.9999999999999999*t177*t7*t80*t81+2.0*t166*t7*
746 t170*t81-.0028392769777778*t1*t2*t4*t175*t43+32.1646831778707*
747 t204*t40/(POW(t41,4.0)*POW(t42,2.0))-2.0*t105*t204*t40/POW(t41,
748 3.0)+.0059059081777778*t1*t2*t4*t175*t24-16.0818243221511*
749 t166*t7*t167*t168+(2.745341666666667*t18*t19*t20*t176+2.754533333333333*
750 t1*t2*t4*t175+1.262325*t13*t14*t174+.6946333333333333*t8*t9*
751 t11*t173)*t40*t104*t105-.1369866666666666*t1*t2*t4*t77*t103*
752 t104*t105);
753 t206 = 2.0*t106*t26*t28*t29*t95;
754 t207 = -8.0*t106*t26*t28*t35*t98;
755 t208 = 8.0*t100*t106*t26*t29*t35;
756 t209 = t26*t28*t29*t193*t44;
757 t210 = -8.0*t26*t28*t44*t95*t98;
758 t211 = 8.0*t100*t26*t29*t44*t95;
759 t212 = 20.0*t199*t26*t28*t35*t44;
760 t213 = -32.0*t100*t26*t35*t44*t98;
761 t214 = 12.0*t201*t26*t29*t35*t44;
762 t215 = 1/POW(t72,2.0);
763 t216 = 1/POW(t70,3.0);
764 t217 = 0.5*t53*t3*t56*t58*t114;
765 t218 = 1/t47;
766 t219 = 1/POW(rhoa,1.333333333333333);
767 t220 = POW(t119,2.0);
768 t221 = 1/POW(2.718281828459045,499.0179939197864*t61*
769 t10*t62*t51*t63);
770 t222 = .8299071502251667*t29*t59*t60*t66;
771 t223 = 1/POW(t49,6.0);
772 t224 = -499.0179939197864*t61*t10*t62*t63;
773 t225 = -166.3393313065955*t46*t61*t10*t219*t60*t51*t63-
774 2661.429300905527*t17*t61*t10*t218*t116*t51*t63+1996.071975679145*
775 t46*t61*t10*t113*t60*t5*t63+t224+998.0359878395727*t46*t61*
776 t10*t113*t60*t51*t117-998.0359878395727*t61*t10*t62*t5*t117-
777 249.5089969598932*t61*t10*t62*t51*(t214+t213+t212+t211+t210+
778 t209+t208+t207+t206+t205+t203+t198+t197+t196+t195+t194+t186+
779 t185+t184+t183+t182+t181+t180+t178+t172+t171+t169);
780 t226 = 1/POW(t65,4.0);
781 t227 = .5582763806977609*t53*t68*t59*t60*t29*t69;
782 t228 = 1.056670602125148*t61*t10*t56*t58*t114*t66;
783 t229 = .0240472290502793*t8*t11*t50*t29*t73;
784 t230 = -.5848223397455204*t26*t30*t83*t31*t84*t132*t85;
785 t231 = -
786 .5848223397455204*t26*t30*t83*t31*t84*t95*t85;
787 t232 = -.5848223397455204*t26*t101*t83*t31*t84*t35*t85;
788 t233 = -
789 .5848223397455204*t26*t134*t83*t31*t84*t35*t85;
790 t234 = 2.666666666666667*t114*t27*t94+.4444444444444444*
791 t131*t190*t93-2.666666666666667*t114*t27*t91+.4444444444444444*
792 t130*t187*t90;
793 t235 = .0197517897025652*t26*t30*t31*t234*t38;
794 t236 = -7.324622014701263e-4*t26*t1*t2*t4*t77*t30*t132*
795 t38;
796 t237 = .0197517897025652*t26*t101*t31*t132*t38;
797 t238 = -7.324622014701263e-4*t26*t1*t2*t4*t77*t30*t95*
798 t38;
799 t239 = .0197517897025652*t26*t134*t31*t95*t38;
800 t240 = -7.324622014701263e-4*t26*t1*t2*t4*t77*t101*t35*
801 t38;
802 t241 = -7.324622014701263e-4*t26*t1*t2*t4*t77*t134*t35*
803 t38;
804 t242 = .0197517897025652*t26*(12.0*t201*t29+t200)*t31*
805 t35*t38;
806 t243 = t26*t28*t29*t132*t106;
807 t244 = t26*t28*t29*t95*t106;
808 t245 = t26*t28*t29*t234*t44;
809 t246 = -4.0*t132*t26*t28*t44*t98;
810 t247 = 4.0*t100*t132*t26*t29*t44;
811 t248 = -4.0*t26*t28*t44*t95*t98;
812 t249 = -4.0*t100*t26*t29*t44*t95;
813 t250 = -12.0*t201*t26*t29*t35*t44;
814 t251 = -2661.429300905527*t17*t61*t10*t113*t116*t138*
815 t51*t63+998.0359878395727*t46*t61*t10*t60*t138*t5*t63+998.0359878395727*
816 t46*t61*t10*t113*t60*t5*t63+t224+499.0179939197864*t46*t61*
817 t10*t60*t138*t51*t117-499.0179939197864*t61*t10*t62*t5*t117+
818 499.0179939197864*t46*t61*t10*t113*t60*t51*t139-499.0179939197864*
819 t61*t10*t62*t5*t139-249.5089969598932*t61*t10*t62*t51*(t250+
820 t212+t249+t248+t247+t246+t245+t207+t244+t243+t205+t242+t198+
821 t241+t240+t239+t238+t237+t236+t235+t186+t233+t232+t184+t183+
822 t182+t231+t230+t180+t178+t172+t171+t169);
823 t252 = .0040078715083799*t8*t11*t50*t52*t148*t111;
824 t253 = -.0080157430167598*t8*t11*t50*t114*t148*t111;
825 t254 = -
826 0.5*t53*t3*grada*t58*t52;
827 t255 = -1.106542866966889*t114*t56*t60*t66*grada;
828 t256 = -.7443685075970145*t53*t68*grada*t56*t60*t114*
829 t69;
830 t257 = -1.056670602125148*t61*t10*grada*t58*t52*t66;
831 t258 = .0040078715083799*
832 t8*t11*t50*t52*t151*t111;
833 t259 = -.0080157430167598*t8*t11*t50*t114*t151*t111;
834 t260 = -
835 0.5*t53*t3*gradb*t58*t52;
836 t261 = -1.106542866966889*t114*t56*t60*t66*gradb;
837 t262 = -.7443685075970145*t53*t68*gradb*t56*t60*t114*
838 t69;
839 t263 = -1.056670602125148*t61*t10*gradb*t58*t52*t66;
840 t264 = .0040078715083799*
841 t8*t11*t50*t52*t158*t111;
842 t265 = -.0080157430167598*t8*t11*t50*t114*t158*t111;
843 t266 = -
844 0.5*t53*t3*t58*t52;
845 t267 = -1.106542866966889*t114*t56*t60*t66;
846 t268 = -.7443685075970145*t53*t68*t56*t60*t114*t69;
847 t269 = -
848 1.056670602125148*t61*t10*t58*t52*t66;
849 t270 = -1.169644679491041*t26*t30*t83*t31*t84*t132*t85;
850 t271 = -
851 1.169644679491041*t26*t134*t83*t31*t84*t35*t85;
852 t272 = 1.333333333333333*(t189+t191)*t94+1.333333333333333*
853 (t192+t188)*t91+.4444444444444444*POW(t131,2.0)*t190+.4444444444444444*
854 POW(t130,2.0)*t187;
855 t273 = .0197517897025652*t26*t30*t31*t272*t38;
856 t274 = -.0014649244029403*t26*t1*t2*t4*t77*t30*t132*t38;
857 t275 = .0395035794051304*
858 t26*t134*t31*t132*t38;
859 t276 = -.0014649244029403*t26*t1*t2*t4*t77*t134*t35*t38;
860 t277 = .0197517897025652*
861 t26*t31*t35*t38*(-32.0*t100*t98+t202+t200);
862 t278 = 2.0*t106*t132*t26*t28*t29;
863 t279 = -8.0*t100*t106*t26*t29*t35;
864 t280 = t26*t28*t29*t272*t44;
865 t281 = -8.0*t132*t26*t28*t44*t98;
866 t282 = -8.0*t100*t132*t26*t29*t44;
867 t283 = 32.0*t100*t26*t35*t44*t98;
868 t284 = 1/POW(rhob,1.333333333333333);
869 t285 = 1/t48;
870 t286 = POW(t140,2.0);
871 t287 = -2661.429300905527*t17*t61*t10*t116*t285*t51*t63-
872 166.3393313065955*t46*t61*t10*t60*t284*t51*t63+1996.071975679145*
873 t46*t61*t10*t60*t138*t5*t63+t224+998.0359878395727*t46*t61*
874 t10*t60*t138*t51*t139-998.0359878395727*t61*t10*t62*t5*t139-
875 249.5089969598932*t61*t10*t62*t51*(t214+t283+t212+t282+t281+
876 t280+t279+t207+t278+t205+t277+t198+t276+t275+t274+t273+t186+
877 t271+t184+t183+t182+t270+t180+t178+t172+t171+t169);
878 t288 = POW(3.141592653589793,1.666666666666667);
879
880 /* code */
881 dfdra = t5*(t129+t128+t127+t110+t109+t108+t107+t102+t97+
882 t96+t87+t86+t82)+t74+t45+t39+t25;
883 dfdrb = t5*(t145+t128+t144+t137+t109+t136+t107+t135+t97+
884 t133+t87+t86+t82)+t74+t45+t39+t25;
885 dfdga = .0040078715083799*t8*t11*t50*t32*t148*t111;
886 dfdgb = .0040078715083799*
887 t8*t11*t50*t32*t151*t111;
888 dfdgab = .0040078715083799*t8*t11*t50*t32*t158*t111;
889 d2fdrara = 2.0*
890 t26*t28*t29*t44*t95+t5*(.0040078715083799*t11*t111*t50*t52*
891 (2.113341204250295*t3*t53*(-.1383178583708611*t225*t52*t59*
892 t60*t64*t69-.1383178583708611*t220*t52*t59*t60*t64*t69+.5532714334834444*
893 t114*t119*t59*t60*t64*t69+.7376952446445925*t113*t116*t119*
894 t46*t52*t59*t64*t69+.1229492074407654*t116*t219*t46*t52*t59*
895 t66+2.458984148815308*t17*t218*t223*t52*t59*t66+1.475390489289185*
896 t113*t114*t116*t46*t59*t66+.2766357167417222*t124*t220*t221*
897 t52*t59*t60+.6666666666666666*t46*t53*t3*t56*t113*t62*t52+
898 .1111111111111111*t46*t53*t3*t56*t219*t62*t32+1.333333333333333*
899 t17*t53*t3*t56*t218*t60*t32+t222+t217)*t71+4.22668240850059*
900 POW(t125,2.0)*t216*t3*t53*t67-2.113341204250295*t53*t3*t67*
901 t121*(-.5283353010625738*t61*t10*t56*t58*t32*t220*t69*t64+
902 1.408894136166863*t46*t61*t10*t56*t113*t62*t32*t119*t69*t64+
903 1.056670602125148*t61*t10*t56*t58*t52*t119*t69*t64-.5283353010625738*
904 t61*t10*t56*t58*t32*t225*t69*t64-.1860921268992536*t53*t68*
905 t59*t60*t52*t220*t124*t64+.9924913434626859*t46*t53*t68*t59*
906 t113*t116*t52*t119*t124*t64+.7443685075970145*t53*t68*t59*
907 t60*t114*t119*t124*t64-.1860921268992536*t53*t68*t59*t60*t52*
908 t225*t124*t64+.2348156893611439*t46*t61*t10*t56*t219*t62*t32*
909 t66+2.817788272333726*t17*t61*t10*t56*t218*t60*t32*t66+1.408894136166863*
910 t46*t61*t10*t56*t113*t62*t52*t66+t228+.0827076119552238*t46*
911 t53*t68*t59*t219*t116*t52*t69+1.654152239104477*t17*t53*t68*
912 t59*t218*t223*t52*t69+.9924913434626859*t46*t53*t68*t59*t113*
913 t116*t114*t69+t227+1.056670602125148*t61*t10*t56*t58*t32*t220*
914 t221*t124+.5582763806977609*t53*t68*t59*t60*t52*t220*t221*
915 t226)-4.22668240850059*t53*t3*t121*t120*t125)*t8-.0040078715083799*
916 t11*POW(t126,2.0)*t215*t50*t52*t8-.0026719143389199*t46*t8*
917 t11*t219*t57*t52*t73+.0213753147113594*t17*t8*t11*t218*t49*
918 t52*t73-.0320629720670391*t46*t8*t11*t113*t57*t114*t73+t229+
919 t214+t213+t212+t211+t210+t209+t208+t207+t206+t205+t203+t198+
920 t197+t196+t195+t194+t186+t185+t184+t183+t182+t181+t180+t178+
921 t172+t171+t169+.0160314860335196*t46*t8*t11*t113*t57*t52*t111*
922 t126-.0160314860335196*t8*t11*t50*t114*t111*t126)+.0160314860335196*
923 t46*t8*t11*t113*t57*t52*t73+8.0*t100*t26*t29*t35*t44+.0395035794051304*
924 t26*t30*t31*t95*t38+.0395035794051304*t26*t101*t31*t35*t38+
925 t165+t164+t163+t162+t161+t160+t159+.0080157430167598*t8*t11*
926 t50*t52*t111*t126;
927 d2fdrarb = t96+t5*(.0040078715083799*t11*t111*t50*t52*
928 (2.113341204250295*t3*t53*(-.1383178583708611*t251*t52*t59*
929 t60*t64*t69-.1383178583708611*t119*t140*t52*t59*t60*t64*t69+
930 .2766357167417222*t114*t140*t59*t60*t64*t69+.2766357167417222*
931 t114*t119*t59*t60*t64*t69+.3688476223222962*t113*t116*t140*
932 t46*t52*t59*t64*t69+.3688476223222962*t116*t119*t138*t46*t52*
933 t59*t64*t69+2.458984148815308*t113*t138*t17*t223*t52*t59*t66+
934 .7376952446445925*t114*t116*t138*t46*t59*t66+.7376952446445925*
935 t113*t114*t116*t46*t59*t66+.2766357167417222*t119*t124*t140*
936 t221*t52*t59*t60+.3333333333333333*t46*t53*t3*t56*t113*t62*
937 t52+.3333333333333333*t46*t53*t3*t56*t62*t138*t52+1.333333333333333*
938 t17*t53*t3*t56*t113*t60*t138*t32+t222+t217)*t71-2.113341204250295*
939 t53*t3*t67*t121*(-.5283353010625738*t61*t10*t56*t58*t32*t119*
940 t140*t69*t64+.7044470680834316*t46*t61*t10*t56*t113*t62*t32*
941 t140*t69*t64+.5283353010625738*t61*t10*t56*t58*t52*t140*t69*
942 t64-.5283353010625738*t61*t10*t56*t58*t32*t251*t69*t64+.7044470680834316*
943 t46*t61*t10*t56*t62*t138*t32*t119*t69*t64+.5283353010625738*
944 t61*t10*t56*t58*t52*t119*t69*t64-.1860921268992536*t53*t68*
945 t59*t60*t52*t119*t140*t124*t64+.4962456717313429*t46*t53*t68*
946 t59*t113*t116*t52*t140*t124*t64+.3721842537985072*t53*t68*
947 t59*t60*t114*t140*t124*t64-.1860921268992536*t53*t68*t59*t60*
948 t52*t251*t124*t64+.4962456717313429*t46*t53*t68*t59*t116*t138*
949 t52*t119*t124*t64+.3721842537985072*t53*t68*t59*t60*t114*t119*
950 t124*t64+2.817788272333726*t17*t61*t10*t56*t113*t60*t138*t32*
951 t66+.7044470680834316*t46*t61*t10*t56*t62*t138*t52*t66+.7044470680834316*
952 t46*t61*t10*t56*t113*t62*t52*t66+t228+1.654152239104477*t17*
953 t53*t68*t59*t113*t223*t138*t52*t69+.4962456717313429*t46*t53*
954 t68*t59*t116*t138*t114*t69+.4962456717313429*t46*t53*t68*t59*
955 t113*t116*t114*t69+t227+1.056670602125148*t61*t10*t56*t58*
956 t32*t119*t140*t221*t124+.5582763806977609*t53*t68*t59*t60*
957 t52*t119*t140*t221*t226)+4.22668240850059*t53*t3*t67*t216*
958 t125*t142-2.113341204250295*t53*t3*t121*t120*t142-2.113341204250295*
959 t53*t3*t121*t125*t141)*t8+.0213753147113594*t17*t8*t11*t113*
960 t49*t138*t52*t73-.0160314860335196*t46*t8*t11*t113*t57*t114*
961 t73-.0160314860335196*t46*t8*t11*t57*t138*t114*t73+t250+t249+
962 t248+t247+t246+t245+t244+t243+t242+t241+t240+t239+t238+t237+
963 t236+t235+t233+t232+t231+t230+t229+t212+t207+t205+t198+t186+
964 t184+t183+t182+t180+t178+t172+t171+t169-.0040078715083799*
965 t8*t11*t50*t52*t215*t126*t143+.0080157430167598*t46*t8*t11*
966 t113*t57*t52*t111*t143-.0080157430167598*t8*t11*t50*t114*t111*
967 t143+.0080157430167598*t46*t8*t11*t57*t138*t52*t111*t126-.0080157430167598*
968 t8*t11*t50*t114*t111*t126)+t165+t164+t163+t162+t161+t160+t159+
969 t145+t144+t136+t135+t133+t129+t127+t108+t102;
970 d2fdraga = t5*(.0040078715083799*t11*t111*t50*t52*t8*
971 (2.113341204250295*t3*t53*t71*(-.5532714334834444*t119*t52*
972 t56*t60*t64*t69*grada-1.475390489289185*t113*t116*t46*t52*
973 t56*t66*grada-.6666666666666666*t46*t53*t3*grada*t113*t62*
974 t32+t255+t254)-2.113341204250295*t53*t3*t67*t121*(-1.056670602125148*
975 t61*t10*grada*t58*t32*t119*t69*t64-.7443685075970145*t53*t68*
976 grada*t56*t60*t52*t119*t124*t64-1.408894136166863*t46*t61*
977 t10*grada*t113*t62*t32*t66+t257-.9924913434626859*t46*t53*
978 t68*grada*t56*t113*t116*t52*t69+t256)+4.22668240850059*t53*
979 t3*t67*t146*t216*t125-2.113341204250295*t53*t3*t147*t121*t125-
980 2.113341204250295*t53*t3*t146*t121*t120)+t253-.0040078715083799*
981 t8*t11*t50*t52*t148*t215*t126+.0080157430167598*t46*t8*t11*
982 t113*t57*t52*t148*t111)+t252;
983 d2fdragb = t5*(.0040078715083799*t11*t111*t50*t52*t8*
984 (2.113341204250295*t3*t53*t71*(-.5532714334834444*t119*t52*
985 t56*t60*t64*t69*gradb-1.475390489289185*t113*t116*t46*t52*
986 t56*t66*gradb-.6666666666666666*t46*t53*t3*gradb*t113*t62*
987 t32+t261+t260)-2.113341204250295*t53*t3*t67*t121*(-1.056670602125148*
988 t61*t10*gradb*t58*t32*t119*t69*t64-.7443685075970145*t53*t68*
989 gradb*t56*t60*t52*t119*t124*t64-1.408894136166863*t46*t61*
990 t10*gradb*t113*t62*t32*t66+t263-.9924913434626859*t46*t53*
991 t68*gradb*t56*t113*t116*t52*t69+t262)+4.22668240850059*t53*
992 t3*t67*t149*t216*t125-2.113341204250295*t53*t3*t150*t121*t125-
993 2.113341204250295*t53*t3*t149*t121*t120)+t259-.0040078715083799*
994 t8*t11*t50*t52*t151*t215*t126+.0080157430167598*t46*t8*t11*
995 t113*t57*t52*t151*t111)+t258;
996 d2fdraab = t5*(.0040078715083799*t11*t111*t50*t52*(2.113341204250295*
997 t3*t53*(-.5532714334834444*t119*t52*t56*t60*t64*t69-1.475390489289185*
998 t113*t116*t46*t52*t56*t66-.6666666666666666*t46*t53*t3*t113*
999 t62*t32+t267+t266)*t71-2.113341204250295*t53*t3*t67*t121*(-
1000 1.056670602125148*t61*t10*t58*t32*t119*t69*t64-.7443685075970145*
1001 t53*t68*t56*t60*t52*t119*t124*t64-1.408894136166863*t46*t61*
1002 t10*t113*t62*t32*t66+t269-.9924913434626859*t46*t53*t68*t56*
1003 t113*t116*t52*t69+t268)+4.22668240850059*t53*t3*t67*t154*t216*
1004 t125-2.113341204250295*t53*t3*t157*t121*t125-2.113341204250295*
1005 t53*t3*t154*t121*t120)*t8+t265-.0040078715083799*t8*t11*t50*
1006 t52*t158*t215*t126+.0080157430167598*t46*t8*t11*t113*t57*t52*
1007 t158*t111)+t264;
1008 d2fdrbrb = t5*(.0040078715083799*t11*t111*t50*t52*(2.113341204250295*
1009 t3*t53*(-.1383178583708611*t287*t52*t59*t60*t64*t69-.1383178583708611*
1010 t286*t52*t59*t60*t64*t69+.5532714334834444*t114*t140*t59*t60*
1011 t64*t69+.7376952446445925*t116*t138*t140*t46*t52*t59*t64*t69+
1012 .1229492074407654*t116*t284*t46*t52*t59*t66+2.458984148815308*
1013 t17*t223*t285*t52*t59*t66+1.475390489289185*t114*t116*t138*
1014 t46*t59*t66+.2766357167417222*t124*t221*t286*t52*t59*t60+.6666666666666666*
1015 t46*t53*t3*t56*t62*t138*t52+1.333333333333333*t17*t53*t3*t56*
1016 t60*t285*t32+.1111111111111111*t46*t53*t3*t56*t62*t284*t32+
1017 t222+t217)*t71+4.22668240850059*POW(t142,2.0)*t216*t3*t53*
1018 t67-2.113341204250295*t53*t3*t67*t121*(-.5283353010625738*
1019 t61*t10*t56*t58*t32*t286*t69*t64+1.408894136166863*t46*t61*
1020 t10*t56*t62*t138*t32*t140*t69*t64+1.056670602125148*t61*t10*
1021 t56*t58*t52*t140*t69*t64-.5283353010625738*t61*t10*t56*t58*
1022 t32*t287*t69*t64-.1860921268992536*t53*t68*t59*t60*t52*t286*
1023 t124*t64+.9924913434626859*t46*t53*t68*t59*t116*t138*t52*t140*
1024 t124*t64+.7443685075970145*t53*t68*t59*t60*t114*t140*t124*
1025 t64-.1860921268992536*t53*t68*t59*t60*t52*t287*t124*t64+2.817788272333726*
1026 t17*t61*t10*t56*t60*t285*t32*t66+.2348156893611439*t46*t61*
1027 t10*t56*t62*t284*t32*t66+1.408894136166863*t46*t61*t10*t56*
1028 t62*t138*t52*t66+t228+1.654152239104477*t17*t53*t68*t59*t223*
1029 t285*t52*t69+.0827076119552238*t46*t53*t68*t59*t116*t284*t52*
1030 t69+.9924913434626859*t46*t53*t68*t59*t116*t138*t114*t69+t227+
1031 1.056670602125148*t61*t10*t56*t58*t32*t286*t221*t124+.5582763806977609*
1032 t53*t68*t59*t60*t52*t286*t221*t226)-4.22668240850059*t53*t3*
1033 t121*t141*t142)*t8-.0040078715083799*t11*POW(t143,2.0)*t215*
1034 t50*t52*t8+.0213753147113594*t17*t8*t11*t49*t285*t52*t73-.0026719143389199*
1035 t46*t8*t11*t57*t284*t52*t73-.0320629720670391*t46*t8*t11*t57*
1036 t138*t114*t73+t283+t282+t281+t280+t279+t278+t277+t276+t275+
1037 t274+t273+t271+t270+t229+t214+t212+t207+t205+t198+t186+t184+
1038 t183+t182+t180+t178+t172+t171+t169+.0160314860335196*t46*t8*
1039 t11*t57*t138*t52*t111*t143-.0160314860335196*t8*t11*t50*t114*
1040 t111*t143)+.0160314860335196*t46*t8*t11*t57*t138*t52*t73-8.0*
1041 t100*t26*t29*t35*t44+2.0*t132*t26*t28*t29*t44+.0395035794051304*
1042 t26*t134*t31*t35*t38+.0395035794051304*t26*t30*t31*t132*t38+
1043 t165+t164+t163+t162+t161+t160+t159+.0080157430167598*t8*t11*
1044 t50*t52*t111*t143;
1045 d2fdrbga = t5*(.0040078715083799*t11*t111*t50*t52*t8*
1046 (2.113341204250295*t3*t53*t71*(-.5532714334834444*t140*t52*
1047 t56*t60*t64*t69*grada-1.475390489289185*t116*t138*t46*t52*
1048 t56*t66*grada-.6666666666666666*t46*t53*t3*grada*t62*t138*
1049 t32+t255+t254)-2.113341204250295*t53*t3*t67*t121*(-1.056670602125148*
1050 t61*t10*grada*t58*t32*t140*t69*t64-.7443685075970145*t53*t68*
1051 grada*t56*t60*t52*t140*t124*t64-1.408894136166863*t46*t61*
1052 t10*grada*t62*t138*t32*t66+t257-.9924913434626859*t46*t53*
1053 t68*grada*t56*t116*t138*t52*t69+t256)+4.22668240850059*t53*
1054 t3*t67*t146*t216*t142-2.113341204250295*t53*t3*t147*t121*t142-
1055 2.113341204250295*t53*t3*t146*t121*t141)+t253-.0040078715083799*
1056 t8*t11*t50*t52*t148*t215*t143+.0080157430167598*t46*t8*t11*
1057 t57*t138*t52*t148*t111)+t252;
1058 d2fdrbgb = t5*(.0040078715083799*t11*t111*t50*t52*t8*
1059 (2.113341204250295*t3*t53*t71*(-.5532714334834444*t140*t52*
1060 t56*t60*t64*t69*gradb-1.475390489289185*t116*t138*t46*t52*
1061 t56*t66*gradb-.6666666666666666*t46*t53*t3*gradb*t62*t138*
1062 t32+t261+t260)-2.113341204250295*t53*t3*t67*t121*(-1.056670602125148*
1063 t61*t10*gradb*t58*t32*t140*t69*t64-.7443685075970145*t53*t68*
1064 gradb*t56*t60*t52*t140*t124*t64-1.408894136166863*t46*t61*
1065 t10*gradb*t62*t138*t32*t66+t263-.9924913434626859*t46*t53*
1066 t68*gradb*t56*t116*t138*t52*t69+t262)+4.22668240850059*t53*
1067 t3*t67*t149*t216*t142-2.113341204250295*t53*t3*t150*t121*t142-
1068 2.113341204250295*t53*t3*t149*t121*t141)+t259-.0040078715083799*
1069 t8*t11*t50*t52*t151*t215*t143+.0080157430167598*t46*t8*t11*
1070 t57*t138*t52*t151*t111)+t258;
1071 d2fdrbgab = t5*(.0040078715083799*t11*t111*t50*t52*(2.113341204250295*
1072 t3*t53*(-.5532714334834444*t140*t52*t56*t60*t64*t69-1.475390489289185*
1073 t116*t138*t46*t52*t56*t66-.6666666666666666*t46*t53*t3*t62*
1074 t138*t32+t267+t266)*t71-2.113341204250295*t53*t3*t67*t121*
1075 (-1.056670602125148*t61*t10*t58*t32*t140*t69*t64-.7443685075970145*
1076 t53*t68*t56*t60*t52*t140*t124*t64-1.408894136166863*t46*t61*
1077 t10*t62*t138*t32*t66+t269-.9924913434626859*t46*t53*t68*t56*
1078 t116*t138*t52*t69+t268)+4.22668240850059*t53*t3*t67*t154*t216*
1079 t142-2.113341204250295*t53*t3*t157*t121*t142-2.113341204250295*
1080 t53*t3*t154*t121*t141)*t8+t265-.0040078715083799*t8*t11*t50*
1081 t52*t158*t215*t143+.0080157430167598*t46*t8*t11*t57*t138*t52*
1082 t158*t111)+t264;
1083 d2fdgaga = .0040078715083799*t11*t111*t32*t50*(2.113341204250295*
1084 t3*t53*(1.106542866966889*t52*t54*t60*t66+t156+t155)*t71+4.22668240850059*
1085 POW(t146,2.0)*t216*t3*t53*t67-2.113341204250295*t53*t3*t67*
1086 (t153+t152+.7443685075970145*t53*t68*t54*t60*t52*t69)*t121-
1087 4.22668240850059*t53*t3*t147*t146*t121)*t8-.0040078715083799*
1088 t11*POW(t148,2.0)*t215*t32*t50*t8;
1089 d2fdgagb = .0040078715083799*t8*t11*t50*t32*t111*(.7443685075970145*
1090 t53*t68*grada*gradb*t60*t52*t71*t66-1.573104638251069*t61*
1091 t288*grada*gradb*t60*t52*t67*t121*t69-2.113341204250295*t53*
1092 t3*t147*t149*t121-2.113341204250295*t53*t3*t150*t146*t121+
1093 4.22668240850059*t53*t3*t67*t146*t149*t216)-.0040078715083799*
1094 t8*t11*t50*t32*t148*t151*t215;
1095 d2fdgagab = .0040078715083799*t8*t11*t50*t32*t111*(.7443685075970145*
1096 t53*t68*grada*t60*t52*t71*t66-1.573104638251069*t61*t288*grada*
1097 t60*t52*t67*t121*t69-2.113341204250295*t53*t3*t157*t146*t121-
1098 2.113341204250295*t53*t3*t147*t154*t121+4.22668240850059*t53*
1099 t3*t67*t154*t146*t216)-.0040078715083799*t8*t11*t50*t32*t158*
1100 t148*t215;
1101 d2fdgbgb = .0040078715083799*t11*t111*t32*t50*(2.113341204250295*
1102 t3*t53*(1.106542866966889*t52*t55*t60*t66+t156+t155)*t71+4.22668240850059*
1103 POW(t149,2.0)*t216*t3*t53*t67-2.113341204250295*t53*t3*t67*
1104 (t153+t152+.7443685075970145*t53*t68*t55*t60*t52*t69)*t121-
1105 4.22668240850059*t53*t3*t150*t149*t121)*t8-.0040078715083799*
1106 t11*POW(t151,2.0)*t215*t32*t50*t8;
1107 d2fdgbgab = .0040078715083799*t8*t11*t50*t32*t111*(.7443685075970145*
1108 t53*t68*gradb*t60*t52*t71*t66-1.573104638251069*t61*t288*gradb*
1109 t60*t52*t67*t121*t69-2.113341204250295*t53*t3*t157*t149*t121-
1110 2.113341204250295*t53*t3*t150*t154*t121+4.22668240850059*t53*
1111 t3*t67*t154*t149*t216)-.0040078715083799*t8*t11*t50*t32*t158*
1112 t151*t215;
1113 d2fdgabgab = .0040078715083799*t11*t111*t32*t50*(-1.573104638251069*
1114 t61*t288*t60*t52*t67*t121*t69+4.22668240850059*POW(t154,2.0)*
1115 t216*t3*t53*t67+.7443685075970145*t53*t68*t60*t52*t71*t66-
1116 4.22668240850059*t53*t3*t157*t154*t121)*t8-.0040078715083799*
1117 t11*POW(t158,2.0)*t215*t32*t50*t8;
1118
1119
1120 ds->df1000 += factor*dfdra;
1121 ds->df0100 += factor*dfdrb;
1122 ds->df0010 += factor*dfdga;
1123 ds->df0001 += factor*dfdgb;
1124 ds->df00001+= factor*dfdgab;
1125
1126 ds->df2000 += factor*d2fdrara;
1127 ds->df1100 += factor*d2fdrarb;
1128 ds->df1010 += factor*d2fdraga;
1129 ds->df1001 += factor*d2fdragb;
1130 ds->df10001+= factor*d2fdraab;
1131 ds->df0200 += factor*d2fdrbrb;
1132 ds->df0110 += factor*d2fdrbga;
1133 ds->df0101 += factor*d2fdrbgb;
1134 ds->df01001+= factor*d2fdrbgab;
1135 ds->df0020 += factor*d2fdgaga;
1136 ds->df0011 += factor*d2fdgagb;
1137 ds->df00101+= factor*d2fdgagab;
1138 ds->df0002 += factor*d2fdgbgb;
1139 ds->df00011+= factor*d2fdgbgab;
1140 ds->df00002+= factor*d2fdgabgab;
1141
1142 }
1143