1 /*
2
3
4 !
5 ! Dalton, a molecular electronic structure program
6 ! Copyright (C) by the authors of Dalton.
7 !
8 ! This program is free software; you can redistribute it and/or
9 ! modify it under the terms of the GNU Lesser General Public
10 ! License version 2.1 as published by the Free Software Foundation.
11 !
12 ! This program is distributed in the hope that it will be useful,
13 ! but WITHOUT ANY WARRANTY; without even the implied warranty of
14 ! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 ! Lesser General Public License for more details.
16 !
17 ! If a copy of the GNU LGPL v2.1 was not distributed with this
18 ! code, you can obtain one at https://www.gnu.org/licenses/old-licenses/lgpl-2.1.en.html.
19 !
20
21 !
22
23 */
24 /*-*-mode: C; c-indentation-style: "bsd"; c-basic-offset: 4; -*-*/
25 /* fun-pw91nc.c:
26
27 Automatically generated code implementing PW91NC functional and
28 its derivatives. It is generated by func-codegen.pl being a part of
29 a "Automatic code generation framework for analytical functional
30 derivative evaluation", Pawel Salek, 2005
31
32 This functional is connected by making following changes:
33 1. add "extern Functional pw91ncFunctional;" to 'functionals.h'
34 2. add "&pw91ncFunctional," to 'functionals.c'
35 3. add "fun-pw91nc.c" to 'Makefile.am', 'Makefile.in' or 'Makefile'.
36
37 This functional has been generated from following input:
38 ------ cut here -------
39 rho: rhoa + rhob;
40 grad: sqrt(grada*grada + gradb*gradb + 2*gradab);
41 zeta: (rhoa-rhob)/(rhoa+rhob);
42
43 c:1.709921;
44 T:[0.031091,0.015545,0.016887];
45 U:[0.21370,0.20548,0.11125];
46 V:[7.5957,14.1189,10.357];
47 W:[3.5876,6.1977,3.6231];
48 X:[1.6382,3.3662,0.88026];
49 Y:[0.49294,0.62517,0.49671];
50 P:[1,1,1];
51
52 r(a,b):=(3/(4*%PI*(a+b)))^(1/3);
53 zet(a,b):=(a-b)/(a+b);
54 omega(z):=((1+z)^(4/3)+(1-z)^(4/3)-2)/(2^(4/3)-2);
55 en(r,t,u,v,w,x,y,p):=-2*t*(1+u*r)*log(1+1/(2*t*(v*sqrt(r)+w*r+x*r^(3/2)+y*r^(p+1))));
56
57 eps(a,b):=en(r(a,b),T[1],U[1],V[1],W[1],X[1],Y[1],P[1]);
58
59 iota:0.09;
60 nu:16/%PI*(3*%PI^2)^(1/3);
61 kappa:0.004235;
62 lamda:nu*kappa;
63 Z:-0.001667;
64 Xi:23.266;
65 Phi:0.007389;
66 Lamda:8.723;
67 Upsilon:0.472;
68
69 theta(r):=10^(-3)*(2.568+Xi*r+Phi*r^2)/(1+Lamda*r+Upsilon*r^2+10*Phi*r^3);
70 phi(r):=theta(r)-Z;
71 sigma(ga,gb,gab):=ga*ga+gb*gb+2*gab;
72 d(a,b,ga,gb,gab):=sqrt(sigma(ga,gb,gab))/4/(3/%PI)^(1/6)/(a+b)^(7/6);
73 A(a,b):=2*iota/lamda*(1/(exp(-2*iota*eps(a,b)/lamda^2)-1));
74 L(d,a,b):=lamda^2/(2*iota)*log((1+2*(iota*(d^2+A(a,b)*d^4))/(lamda*(1+A(a,b)*d^2+(A(a,b))^2*d^4))));
75 J(d,a,b):=nu*(phi(r(a,b))-kappa-3*Z/7)*d^2*exp(-100*4/(3*%PI^5*(a+b))^(1/3)*d^2);
76 H(d,a,b):=L(d,a,b)+J(d,a,b);
77
78 K(rhoa,grada,rhob,gradb,gradab):=(rhoa+rhob)*(H(d(rhoa,rhob,grada,gradb,gradab),rhoa,rhob));
79
80
81 ------ cut here -------
82 */
83
84
85 /* strictly conform to XOPEN ANSI C standard */
86 #if !defined(SYS_DEC)
87 /* XOPEN compliance is missing on old Tru64 4.0E Alphas and pow() prototype
88 * is not specified. */
89 #define _XOPEN_SOURCE 500
90 #define _XOPEN_SOURCE_EXTENDED 1
91 #endif
92 #include <math.h>
93 #include <stddef.h>
94 #include "general.h"
95
96 #define __CVERSION__
97
98 #include "functionals.h"
99
100 /* INTERFACE PART */
pw91nc_isgga(void)101 static integer pw91nc_isgga(void) { return 1; } /* FIXME: detect! */
102 static integer pw91nc_read(const char *conf_line);
103 static real pw91nc_energy(const FunDensProp* dp);
104 static void pw91nc_first(FunFirstFuncDrv *ds, real factor,
105 const FunDensProp* dp);
106 static void pw91nc_second(FunSecondFuncDrv *ds, real factor,
107 const FunDensProp* dp);
108 static void pw91nc_third(FunThirdFuncDrv *ds, real factor,
109 const FunDensProp* dp);
110
111 Functional PW91ncFunctional = {
112 "PW91nc", /* name */
113 pw91nc_isgga, /* gga-corrected */
114 1,
115 pw91nc_read,
116 NULL,
117 pw91nc_energy,
118 pw91nc_first,
119 pw91nc_second,
120 pw91nc_third
121 };
122
123 /* IMPLEMENTATION PART */
124 static integer
pw91nc_read(const char * conf_line)125 pw91nc_read(const char *conf_line)
126 {
127 fun_set_hf_weight(0);
128 return 1;
129 }
130
131 static real
pw91nc_energy(const FunDensProp * dp)132 pw91nc_energy(const FunDensProp *dp)
133 {
134 real res;
135 real rhoa = dp->rhoa, rhob = dp->rhob;
136 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
137
138 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
139 real t11, t12, t13, t14, t15, t16, t17, t18;
140 real t19, t20;
141
142 t1 = rhob+rhoa;
143 t2 = pow(gradb,2.0)+2.0*gradab+pow(grada,2.0);
144 t3 = 1/pow(t1,2.333333333333333);
145 t4 = pow(3.0,0.666666666666667);
146 t5 = 1/pow(4.0,0.666666666666667);
147 t6 = pow(3.141592653589793,0.666666666666667);
148 t7 = 1/t6;
149 t8 = 1/pow(t1,0.666666666666667);
150 t9 = pow(3.0,0.333333333333333);
151 t10 = 1/pow(4.0,0.333333333333333);
152 t11 = pow(3.141592653589793,0.333333333333333);
153 t12 = 1/t11;
154 t13 = 1/pow(t1,0.333333333333333);
155 t14 = 1/t4;
156 t15 = pow(3.141592653589793,1.333333333333333);
157 t16 = 1/t9;
158 t17 = pow(t2,2.0);
159 t18 = 1/pow(t1,4.666666666666667);
160 t19 = pow(2.718281828459045,2.43775691411733*(0.2137*
161 t9*t10*t12*t13+1.0)*t14*t6*log(16.0818243221511/(0.49294*t4*
162 t5*t7*t8+3.5876*t9*t10*t12*t13+7.5957*pow(3.0,0.166666666666667)/
163 (pow(2.0,0.333333333333333)*pow(3.141592653589793,0.166666666666667)*
164 pow(t1,0.166666666666667))+0.800428634999363/sqrt(t1))+1.0))-
165 1.0;
166 t20 = 1/t19;
167
168 /* code */
169 res = t1*(t2*t3*(0.001*(23.266*t9*t10*t12*t13+0.007389*
170 t4*t5*t7*t8+2.568)/(0.472*t4*t5*t7*t8+8.723*t9*t10*t12*t13+
171 0.01763993811759/t1+1.0)-0.001853571428571)/pow(2.718281828459045,
172 25.0*t14*t2/(pow(t1,2.666666666666667)*t15))+0.025507875555556*
173 log(2.656434474616293*t11*t16*(0.0625*t11*t16*t2*t3+0.010866451862236*
174 t17*t18*t20)/(0.166027154663518*t14*t6*t2*t3*t20+0.009188338695221*
175 t15*t16*t17*t18/pow(t19,2.0)+1.0)+1.0)*t4*t7);
176
177 return res;
178 }
179
180 static void
pw91nc_first(FunFirstFuncDrv * ds,real factor,const FunDensProp * dp)181 pw91nc_first(FunFirstFuncDrv *ds, real factor, const FunDensProp *dp)
182 {
183 real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
184 real rhoa = dp->rhoa, rhob = dp->rhob;
185 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
186
187 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
188 real t11, t12, t13, t14, t15, t16, t17, t18;
189 real t19, t20, t21, t22, t23, t24, t25, t26;
190 real t27, t28, t29, t30, t31, t32, t33, t34;
191 real t35, t36, t37, t38, t39, t40, t41, t42;
192 real t43, t44, t45, t46, t47, t48, t49, t50;
193 real t51, t52;
194
195 t1 = pow(gradb,2.0)+2.0*gradab+pow(grada,2.0);
196 t2 = rhob+rhoa;
197 t3 = 1/pow(t2,2.333333333333333);
198 t4 = 0.318309886183791;
199 t5 = pow(3.0,0.666666666666667);
200 t6 = 1/pow(4.0,0.666666666666667);
201 t7 = pow(3.141592653589793,0.666666666666667);
202 t8 = 1/t7;
203 t9 = 1/pow(t2,0.666666666666667);
204 t10 = pow(3.0,0.333333333333333);
205 t11 = 1/pow(4.0,0.333333333333333);
206 t12 = pow(3.141592653589793,0.333333333333333);
207 t13 = 1/t12;
208 t14 = 1/pow(t2,0.333333333333333);
209 t15 = 0.472*t5*t6*t8*t9+0.0554175*t4/t2+8.723*t10*t11*
210 t13*t14+1.0;
211 t16 = 1/t15;
212 t17 = 23.266*t10*t11*t13*t14+0.007389*t5*t6*t8*t9+2.568;
213 t18 = 0.001*
214 t16*t17-0.001853571428571;
215 t19 = 1/t5;
216 t20 = pow(3.141592653589793,1.333333333333333);
217 t21 = 1/t20;
218 t22 = 1/pow(2.718281828459045,25.0*t1*t19*t21/pow(t2,
219 2.666666666666667));
220 t23 = 1/pow(t2,1.666666666666667);
221 t24 = 1/pow(t2,1.333333333333333);
222 t25 = pow(t1,2.0);
223 t26 = 1/pow(t2,3.333333333333333);
224 t27 = 1/t10;
225 t28 = 1/pow(t2,4.666666666666667);
226 t29 = 0.2137*t10*t11*t13*t14+1.0;
227 t30 = 1.732050807568877;
228 t31 = 0.564189583547756;
229 t32 = sqrt(t2);
230 t33 = 1/pow(2.0,0.333333333333333);
231 t34 = pow(3.0,0.166666666666667);
232 t35 = 1/pow(3.141592653589793,0.166666666666667);
233 t36 = 0.49294*t5*t6*t8*t9+7.5957*t33*t34*t35/pow(t2,0.166666666666667)+
234 0.8191*t30*t31/t32+3.5876*t10*t11*t13*t14;
235 t37 = 16.0818243221511/t36+1.0;
236 t38 = log(t37);
237 t39 = pow(2.718281828459045,2.43775691411733*t19*t7*t29*
238 t38);
239 t40 = t39-1.0;
240 t41 = 1/t40;
241 t42 = 0.010866451862236*t25*t28*t41+0.0625*t1*t12*t27*
242 t3;
243 t43 = 1/pow(t40,2.0);
244 t44 = 0.009188338695221*t27*t20*t25*t28*t43+0.166027154663518*
245 t19*t7*t1*t3*t41+1.0;
246 t45 = 1/t44;
247 t46 = 2.656434474616293*t27*t12*t42*t45+1.0;
248 t47 = 1/t46;
249 t48 = 1/pow(t2,5.666666666666667);
250 t49 = -39.2035784329441*t19*t29*(-1.26595*t33*t34*t35/
251 pow(t2,1.166666666666667)-0.40955*t30*t31/pow(t32,3.0)-1.195866666666667*
252 t10*t11*t13*t24-0.328626666666667*t5*t6*t8*t23)*t7/(pow(t36,
253 2.0)*t37)-0.173649550848958*t27*t11*t12*t24*t38;
254 t50 = 1/pow(t44,2.0);
255 t51 = t2*(0.025507875555556*t47*t5*(2.656434474616293*
256 t12*t27*t45*(-0.010866451862236*t25*t28*t39*t43*t49-0.050710108690434*
257 t25*t41*t48-0.145833333333333*t1*t12*t26*t27)-2.656434474616293*
258 t12*t27*t42*(-0.018376677390443*t20*t25*t27*t28*t39*t49/pow(t40,
259 3.0)-0.042878913911033*t27*t20*t25*t48*t43-0.387396694214876*
260 t19*t7*t1*t26*t41-0.166027154663518*t19*t7*t1*t3*t49*t43*t39)*
261 t50)*t8+t1*t22*t3*(0.001*t16*(-7.755333333333333*t10*t11*t13*
262 t24-0.004926*t5*t6*t8*t23)-0.001*t17*(-0.0554175*t4/pow(t2,
263 2.0)-2.907666666666667*t10*t11*t13*t24-0.314666666666667*t5*
264 t6*t8*t23)/pow(t15,2.0))-2.333333333333333*t1*t18*t22*t26+
265 66.66666666666667*t18*t19*t21*t22*t25/pow(t2,6.0))+0.025507875555556*
266 t5*t8*log(t46)+t1*t3*t18*t22;
267 t52 = 1/pow(t2,5.0);
268
269 /* code */
270 dfdra = t51;
271 dfdrb = t51;
272 dfdga = t2*(0.025507875555556*t47*t5*t8*(2.656434474616293*
273 t12*t27*t45*(0.043465807448944*t1*t28*t41*grada+0.125*t12*
274 t27*t3*grada)-2.656434474616293*t27*t12*t42*(0.332054309327037*
275 t19*t7*grada*t3*t41+0.036753354780885*t27*t20*grada*t1*t28*
276 t43)*t50)-50.0*t1*t18*t19*t21*t22*t52*grada+2.0*t18*t22*t3*
277 grada);
278 dfdgb = t2*(0.025507875555556*t47*t5*t8*(2.656434474616293*
279 t12*t27*t45*(0.043465807448944*t1*t28*t41*gradb+0.125*t12*
280 t27*t3*gradb)-2.656434474616293*t27*t12*t42*(0.332054309327037*
281 t19*t7*gradb*t3*t41+0.036753354780885*t27*t20*gradb*t1*t28*
282 t43)*t50)-50.0*t1*t18*t19*t21*t22*t52*gradb+2.0*t18*t22*t3*
283 gradb);
284 dfdgab = t2*(0.025507875555556*t47*t5*(2.656434474616293*
285 t12*t27*(0.043465807448944*t1*t28*t41+0.125*t12*t27*t3)*t45-
286 2.656434474616293*t27*t12*t42*(0.332054309327037*t19*t7*t3*
287 t41+0.036753354780885*t27*t20*t1*t28*t43)*t50)*t8-50.0*t1*
288 t18*t19*t21*t22*t52+2.0*t18*t22*t3);
289
290
291 ds->df1000 += factor*dfdra;
292 ds->df0100 += factor*dfdrb;
293 ds->df0010 += factor*dfdga;
294 ds->df0001 += factor*dfdgb;
295 ds->df00001+= factor*dfdgab;
296
297 }
298
299 static void
pw91nc_second(FunSecondFuncDrv * ds,real factor,const FunDensProp * dp)300 pw91nc_second(FunSecondFuncDrv *ds, real factor, const FunDensProp* dp)
301 {
302 real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
303 real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
304 d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
305 d2fdgbgb, d2fdgbgab, d2fdgabgab;
306 real rhoa = dp->rhoa, rhob = dp->rhob;
307 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
308
309 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
310 real t11, t12, t13, t14, t15, t16, t17, t18;
311 real t19, t20, t21, t22, t23, t24, t25, t26;
312 real t27, t28, t29, t30, t31, t32, t33, t34;
313 real t35, t36, t37, t38, t39, t40, t41, t42;
314 real t43, t44, t45, t46, t47, t48, t49, t50;
315 real t51, t52, t53, t54, t55, t56, t57, t58;
316 real t59, t60, t61, t62, t63, t64, t65, t66;
317 real t67, t68, t69, t70, t71, t72, t73, t74;
318 real t75, t76, t77, t78, t79, t80, t81, t82;
319 real t83, t84, t85, t86, t87, t88, t89, t90;
320 real t91, t92, t93, t94, t95, t96, t97, t98;
321 real t99, t100, t101, t102, t103, t104, t105;
322
323 t1 = pow(grada,2.0);
324 t2 = pow(gradb,2.0);
325 t3 = 2.0*gradab+t2+t1;
326 t4 = rhob+rhoa;
327 t5 = 1/pow(t4,2.333333333333333);
328 t6 = 0.318309886183791;
329 t7 = pow(3.0,0.666666666666667);
330 t8 = 1/pow(4.0,0.666666666666667);
331 t9 = pow(3.141592653589793,0.666666666666667);
332 t10 = 1/t9;
333 t11 = 1/pow(t4,0.666666666666667);
334 t12 = pow(3.0,0.333333333333333);
335 t13 = 1/pow(4.0,0.333333333333333);
336 t14 = pow(3.141592653589793,0.333333333333333);
337 t15 = 1/t14;
338 t16 = 1/pow(t4,0.333333333333333);
339 t17 = 0.0554175*t6/t4+8.723*t12*t13*t15*t16+0.472*t7*
340 t8*t10*t11+1.0;
341 t18 = 1/t17;
342 t19 = 23.266*t12*t13*t15*t16+0.007389*t7*t8*t10*t11+2.568;
343 t20 = 0.001*
344 t18*t19-0.001853571428571;
345 t21 = 1/t7;
346 t22 = pow(3.141592653589793,1.333333333333333);
347 t23 = 1/t22;
348 t24 = 1/pow(t4,2.666666666666667);
349 t25 = 1/pow(2.718281828459045,25.0*t21*t23*t24*t3);
350 t26 = 1/
351 pow(t4,1.666666666666667);
352 t27 = 1/pow(t4,1.333333333333333);
353 t28 = -7.755333333333333*t12*t13*t15*t27-0.004926*t7*
354 t8*t10*t26;
355 t29 = -0.0554175*t6/pow(t4,2.0)-2.907666666666667*t12*
356 t13*t15*t27-0.314666666666667*t7*t8*t10*t26;
357 t30 = 1/pow(t17,2.0);
358 t31 = 0.001*t18*t28-0.001*t19*t29*t30;
359 t32 = pow(t3,2.0);
360 t33 = 1/pow(t4,6.0);
361 t34 = 1/pow(t4,3.333333333333333);
362 t35 = 1/t12;
363 t36 = 1/pow(t4,4.666666666666667);
364 t37 = 0.2137*t12*t13*t15*t16+1.0;
365 t38 = 1.732050807568877;
366 t39 = 0.564189583547756;
367 t40 = sqrt(t4);
368 t41 = 1/pow(2.0,0.333333333333333);
369 t42 = pow(3.0,0.166666666666667);
370 t43 = 1/pow(3.141592653589793,0.166666666666667);
371 t44 = 7.5957*t41*t42*t43/pow(t4,0.166666666666667)+0.8191*
372 t38*t39/t40+3.5876*t12*t13*t15*t16+0.49294*t7*t8*t10*t11;
373 t45 = 16.0818243221511/
374 t44+1.0;
375 t46 = log(t45);
376 t47 = pow(2.718281828459045,2.43775691411733*t21*t9*t37*
377 t46);
378 t48 = t47-1.0;
379 t49 = 1/t48;
380 t50 = 0.0625*t14*t3*t35*t5+0.010866451862236*t32*t36*
381 t49;
382 t51 = 1/pow(t48,2.0);
383 t52 = 0.009188338695221*t35*t22*t32*t36*t51+0.166027154663518*
384 t21*t9*t3*t5*t49+1.0;
385 t53 = 1/t52;
386 t54 = 2.656434474616293*t35*t14*t50*t53+1.0;
387 t55 = 1/t54;
388 t56 = 1/pow(t4,5.666666666666667);
389 t57 = -1.26595*t41*t42*t43/pow(t4,1.166666666666667)-
390 0.40955*t38*t39/pow(t40,3.0)-1.195866666666667*t12*t13*t15*
391 t27-0.328626666666667*t7*t8*t10*t26;
392 t58 = 1/pow(t44,2.0);
393 t59 = 1/t45;
394 t60 = -0.173649550848958*t35*t13*t14*t27*t46-39.2035784329441*
395 t21*t9*t57*t37*t58*t59;
396 t61 = -0.010866451862236*t32*t36*t47*t51*t60-0.050710108690434*
397 t32*t49*t56-0.145833333333333*t14*t3*t34*t35;
398 t62 = 1/pow(t52,2.0);
399 t63 = 1/pow(t48,3.0);
400 t64 = -0.166027154663518*t21*t9*t3*t5*t60*t51*t47-0.018376677390443*
401 t35*t22*t32*t36*t60*t63*t47-0.387396694214876*t21*t9*t3*t34*
402 t49-0.042878913911033*t35*t22*t32*t56*t51;
403 t65 = 2.656434474616293*t35*t14*t53*t61-2.656434474616293*
404 t35*t14*t50*t62*t64;
405 t66 = t4*(0.025507875555556*t7*t10*t55*t65-2.333333333333333*
406 t20*t25*t3*t34+66.66666666666667*t20*t21*t23*t25*t32*t33+t3*
407 t5*t31*t25)+0.025507875555556*t7*t10*log(t54)+t3*t5*t20*t25;
408 t67 = 1/
409 pow(t4,5.0);
410 t68 = -50.0*t20*t21*t23*t25*t3*t67*grada;
411 t69 = 2.0*t20*t25*t5*grada;
412 t70 = 0.332054309327037*t21*t9*grada*t5*t49+0.036753354780885*
413 t35*t22*grada*t3*t36*t51;
414 t71 = 0.125*t14*t35*t5*grada+0.043465807448944*t3*t36*
415 t49*grada;
416 t72 = 2.656434474616293*t35*t14*t71*t53-2.656434474616293*
417 t35*t14*t50*t70*t62;
418 t73 = 0.025507875555556*t7*t10*t72*t55;
419 t74 = -50.0*t20*t21*t23*t25*t3*t67*gradb;
420 t75 = 2.0*t20*t25*t5*gradb;
421 t76 = 0.332054309327037*t21*t9*gradb*t5*t49+0.036753354780885*
422 t35*t22*gradb*t3*t36*t51;
423 t77 = 0.125*t14*t35*t5*gradb+0.043465807448944*t3*t36*
424 t49*gradb;
425 t78 = 2.656434474616293*t35*t14*t77*t53-2.656434474616293*
426 t35*t14*t50*t76*t62;
427 t79 = 0.025507875555556*t7*t10*t78*t55;
428 t80 = -50.0*t20*t21*t23*t25*t3*t67;
429 t81 = 2.0*t20*t25*t5;
430 t82 = 0.036753354780885*t35*t22*t3*t36*t51;
431 t83 = 0.332054309327037*t21*t9*t5*t49;
432 t84 = t83+t82;
433 t85 = 0.125*t14*t35*t5;
434 t86 = 0.043465807448944*t3*t36*t49;
435 t87 = t86+t85;
436 t88 = 2.656434474616293*t35*t14*t87*t53-2.656434474616293*
437 t35*t14*t50*t84*t62;
438 t89 = 0.025507875555556*t7*t10*t88*t55;
439 t90 = 1/pow(3.141592653589793,2.666666666666667);
440 t91 = 1/pow(t4,4.333333333333333);
441 t92 = 1/pow(t54,2.0);
442 t93 = 1/pow(t52,3.0);
443 t94 = 1/pow(t4,6.666666666666667);
444 t95 = pow(t57,2.0);
445 t96 = 78.40715686588821*t21*t37*t59*t9*t95/pow(t44,3.0)-
446 630.4650611582789*t21*t37*t9*t95/(pow(t44,4.0)*pow(t45,2.0))-
447 39.2035784329441*t21*t37*(1.594488888888889*t12*t13*t15*t5+
448 1.476941666666667*t41*t42*t43/pow(t4,2.166666666666667)+0.614325*
449 t38*t39/pow(t40,5.0)+0.547711111111111*t7*t8*t10*t24)*t58*
450 t59*t9+5.58520314074677*t35*t13*t14*t27*t57*t58*t59+0.231532734465277*
451 t35*t13*t14*t5*t46;
452 t97 = pow(t60,2.0);
453 t98 = pow(2.718281828459045,4.875513828234661*t21*t9*
454 t37*t46);
455 t99 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
456 t14*t35*t53*(0.021732903724472*t32*t36*t63*t97*t98-0.010866451862236*
457 t32*t36*t47*t51*t97-0.010866451862236*t32*t36*t47*t51*t96+
458 0.287357282579128*t32*t49*t94+0.486111111111111*t14*t3*t35*
459 t91+0.101420217380869*t32*t47*t51*t56*t60)-2.656434474616293*
460 t14*t35*t50*t62*(0.055130032171328*t22*t32*t35*t36*t97*t98/
461 pow(t48,4.0)+0.332054309327037*t21*t9*t3*t5*t97*t63*t98+0.242980512162518*
462 t35*t22*t32*t94*t51+1.291322314049587*t21*t9*t3*t91*t49-0.018376677390443*
463 t35*t22*t32*t36*t97*t63*t47-0.018376677390443*t35*t22*t32*
464 t36*t96*t63*t47+0.17151565564413*t35*t22*t32*t56*t60*t63*t47-
465 0.166027154663518*t21*t9*t3*t5*t97*t51*t47-0.166027154663518*
466 t21*t9*t3*t5*t96*t51*t47+0.774793388429752*t21*t9*t3*t34*t60*
467 t51*t47)+5.312868949232586*t14*t35*t50*pow(t64,2.0)*t93-5.312868949232586*
468 t35*t14*t62*t61*t64)-0.025507875555556*t10*pow(t65,2.0)*t7*
469 t92+7.777777777777778*t20*t25*t3*t91+1481.481481481482*t20*
470 t25*pow(t3,3.0)*t35*t90/pow(t4,9.666666666666666)+t25*t3*t5*
471 (-0.001*t19*t30*(0.110835*t6/pow(t4,3.0)+3.876888888888889*
472 t12*t13*t15*t5+0.524444444444444*t7*t8*t10*t24)+0.001*t18*
473 (10.34044444444444*t12*t13*t15*t5+0.00821*t7*t8*t10*t24)-0.002*
474 t28*t29*t30+0.002*t19*pow(t29,2.0)/pow(t17,3.0))-555.5555555555555*
475 t20*t21*t23*t25*t32/pow(t4,7.0)-4.666666666666667*t25*t3*t31*
476 t34+133.3333333333333*t21*t23*t25*t31*t32*t33)+0.051015751111111*
477 t7*t10*t55*t65+2.0*t25*t3*t31*t5-4.666666666666667*t20*t25*
478 t3*t34+133.3333333333333*t20*t21*t23*t25*t32*t33;
479 t100 = 1/pow(t4,8.666666666666666);
480 t101 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
481 t14*t35*t53*(-0.043465807448944*t3*t36*t47*t51*t60*grada-0.202840434761738*
482 t3*t49*t56*grada-0.291666666666667*t14*t34*t35*grada)+5.312868949232586*
483 t35*t14*t50*t70*t93*t64-2.656434474616293*t35*t14*t71*t62*
484 t64-2.656434474616293*t35*t14*t70*t62*t61-2.656434474616293*
485 t35*t14*t50*t62*(-0.332054309327037*t21*t9*grada*t5*t60*t51*
486 t47-0.07350670956177*t35*t22*grada*t3*t36*t60*t63*t47-0.774793388429752*
487 t21*t9*grada*t34*t49-0.17151565564413*t35*t22*grada*t3*t56*
488 t51))-1111.111111111111*t100*t20*t25*t32*t35*t90*grada-50.0*
489 t21*t23*t25*t3*t31*t67*grada+2.0*t25*t31*t5*grada-4.666666666666667*
490 t20*t25*t34*grada+383.3333333333333*t20*t21*t23*t25*t3*t33*
491 grada-0.025507875555556*t7*t10*t72*t92*t65)+t73+t69+t68;
492 t102 = t4*
493 (0.025507875555556*t10*t55*t7*(2.656434474616293*t14*t35*t53*
494 (-0.043465807448944*t3*t36*t47*t51*t60*gradb-0.202840434761738*
495 t3*t49*t56*gradb-0.291666666666667*t14*t34*t35*gradb)+5.312868949232586*
496 t35*t14*t50*t76*t93*t64-2.656434474616293*t35*t14*t77*t62*
497 t64-2.656434474616293*t35*t14*t76*t62*t61-2.656434474616293*
498 t35*t14*t50*t62*(-0.332054309327037*t21*t9*gradb*t5*t60*t51*
499 t47-0.07350670956177*t35*t22*gradb*t3*t36*t60*t63*t47-0.774793388429752*
500 t21*t9*gradb*t34*t49-0.17151565564413*t35*t22*gradb*t3*t56*
501 t51))-1111.111111111111*t100*t20*t25*t32*t35*t90*gradb-50.0*
502 t21*t23*t25*t3*t31*t67*gradb+2.0*t25*t31*t5*gradb-4.666666666666667*
503 t20*t25*t34*gradb+383.3333333333333*t20*t21*t23*t25*t3*t33*
504 gradb-0.025507875555556*t7*t10*t78*t92*t65)+t79+t75+t74;
505 t103 = t4*
506 (-1111.111111111111*t100*t20*t25*t32*t35*t90+0.025507875555556*
507 t10*t55*(5.312868949232586*t35*t14*t50*t84*t93*t64-2.656434474616293*
508 t35*t14*t87*t62*t64-2.656434474616293*t35*t14*t84*t62*t61+
509 2.656434474616293*t14*t35*t53*(-0.043465807448944*t3*t36*t47*
510 t51*t60-0.202840434761738*t3*t49*t56-0.291666666666667*t14*
511 t34*t35)-2.656434474616293*t35*t14*t50*t62*(-0.332054309327037*
512 t21*t9*t5*t60*t51*t47-0.07350670956177*t35*t22*t3*t36*t60*
513 t63*t47-0.774793388429752*t21*t9*t34*t49-0.17151565564413*
514 t35*t22*t3*t56*t51))*t7-50.0*t21*t23*t25*t3*t31*t67-0.025507875555556*
515 t7*t10*t88*t92*t65+2.0*t25*t31*t5-4.666666666666667*t20*t25*
516 t34+383.3333333333333*t20*t21*t23*t25*t3*t33)+t89+t81+t80;
517 t104 = 1/
518 pow(t4,7.666666666666667);
519 t105 = pow(3.141592653589793,1.666666666666667);
520
521 /* code */
522 dfdra = t66;
523 dfdrb = t66;
524 dfdga = t4*(t73+t69+t68);
525 dfdgb = t4*(t79+t75+t74);
526 dfdgab = t4*(t89+t81+t80);
527 d2fdrara = t99;
528 d2fdrarb = t99;
529 d2fdraga = t101;
530 d2fdragb = t102;
531 d2fdraab = t103;
532 d2fdrbrb = t99;
533 d2fdrbga = t101;
534 d2fdrbgb = t102;
535 d2fdrbgab = t103;
536 d2fdgaga = t4*(0.025507875555556*t10*t55*t7*(5.312868949232586*
537 t14*t35*t50*pow(t70,2.0)*t93+2.656434474616293*t14*t35*t53*
538 (t86+t85+0.086931614897888*t1*t36*t49)-2.656434474616293*t35*
539 t14*t50*(t83+t82+0.07350670956177*t35*t22*t1*t36*t51)*t62-
540 5.312868949232586*t35*t14*t71*t70*t62)-0.025507875555556*t10*
541 t7*pow(t72,2.0)*t92+833.3333333333334*t1*t104*t20*t25*t3*t35*
542 t90+t81+t80-200.0*t1*t20*t21*t23*t25*t67);
543 d2fdgagb = t4*(833.3333333333334*t104*t20*t25*t3*t35*
544 t90*grada*gradb-200.0*t20*t21*t23*t25*t67*grada*gradb-0.025507875555556*
545 t7*t10*t72*t78*t92+0.025507875555556*t7*t10*t55*(0.07350670956177*
546 t35*t22*grada*gradb*t36*t53*t49-0.195265757395493*t21*t105*
547 grada*gradb*t36*t50*t62*t51-2.656434474616293*t35*t14*t71*
548 t76*t62-2.656434474616293*t35*t14*t77*t70*t62+5.312868949232586*
549 t35*t14*t50*t70*t76*t93));
550 d2fdgagab = t4*(833.3333333333334*t104*t20*t25*t3*t35*
551 t90*grada-200.0*t20*t21*t23*t25*t67*grada-0.025507875555556*
552 t7*t10*t88*t72*t92+0.025507875555556*t7*t10*t55*(0.07350670956177*
553 t35*t22*grada*t36*t53*t49-0.195265757395493*t21*t105*grada*
554 t36*t50*t62*t51-2.656434474616293*t35*t14*t87*t70*t62-2.656434474616293*
555 t35*t14*t71*t84*t62+5.312868949232586*t35*t14*t50*t84*t70*
556 t93));
557 d2fdgbgb = t4*(0.025507875555556*t10*t55*t7*(5.312868949232586*
558 t14*t35*t50*pow(t76,2.0)*t93+2.656434474616293*t14*t35*t53*
559 (t86+t85+0.086931614897888*t2*t36*t49)-2.656434474616293*t35*
560 t14*t50*(t83+t82+0.07350670956177*t35*t22*t2*t36*t51)*t62-
561 5.312868949232586*t35*t14*t77*t76*t62)-0.025507875555556*t10*
562 t7*pow(t78,2.0)*t92+833.3333333333334*t104*t2*t20*t25*t3*t35*
563 t90+t81+t80-200.0*t2*t20*t21*t23*t25*t67);
564 d2fdgbgab = t4*(833.3333333333334*t104*t20*t25*t3*t35*
565 t90*gradb-200.0*t20*t21*t23*t25*t67*gradb-0.025507875555556*
566 t7*t10*t88*t78*t92+0.025507875555556*t7*t10*t55*(0.07350670956177*
567 t35*t22*gradb*t36*t53*t49-0.195265757395493*t21*t105*gradb*
568 t36*t50*t62*t51-2.656434474616293*t35*t14*t87*t76*t62-2.656434474616293*
569 t35*t14*t77*t84*t62+5.312868949232586*t35*t14*t50*t84*t76*
570 t93));
571 d2fdgabgab = t4*(0.025507875555556*t10*t55*t7*(5.312868949232586*
572 t14*t35*t50*pow(t84,2.0)*t93-5.312868949232586*t35*t14*t87*
573 t84*t62-0.195265757395493*t21*t105*t36*t50*t62*t51+0.07350670956177*
574 t35*t22*t36*t53*t49)-0.025507875555556*t10*t7*pow(t88,2.0)*
575 t92+833.3333333333334*t104*t20*t25*t3*t35*t90-200.0*t20*t21*
576 t23*t25*t67);
577
578
579 ds->df1000 += factor*dfdra;
580 ds->df0100 += factor*dfdrb;
581 ds->df0010 += factor*dfdga;
582 ds->df0001 += factor*dfdgb;
583 ds->df00001+= factor*dfdgab;
584
585 ds->df2000 += factor*d2fdrara;
586 ds->df1100 += factor*d2fdrarb;
587 ds->df1010 += factor*d2fdraga;
588 ds->df1001 += factor*d2fdragb;
589 ds->df10001+= factor*d2fdraab;
590 ds->df0200 += factor*d2fdrbrb;
591 ds->df0110 += factor*d2fdrbga;
592 ds->df0101 += factor*d2fdrbgb;
593 ds->df01001+= factor*d2fdrbgab;
594 ds->df0020 += factor*d2fdgaga;
595 ds->df0011 += factor*d2fdgagb;
596 ds->df00101+= factor*d2fdgagab;
597 ds->df0002 += factor*d2fdgbgb;
598 ds->df00011+= factor*d2fdgbgab;
599 ds->df00002+= factor*d2fdgabgab;
600
601 }
602
603 static void
pw91nc_third(FunThirdFuncDrv * ds,real factor,const FunDensProp * dp)604 pw91nc_third(FunThirdFuncDrv *ds, real factor, const FunDensProp* dp)
605 {
606 real dfdra, dfdrb, dfdga, dfdgb, dfdgab;
607 real d2fdrara, d2fdrarb, d2fdraga, d2fdragb, d2fdraab, d2fdrbrb,
608 d2fdrbga, d2fdrbgb, d2fdrbgab, d2fdgaga, d2fdgagb, d2fdgagab,
609 d2fdgbgb, d2fdgbgab, d2fdgabgab;
610 real d3fdrarara, d3fdrararb, d3fdraraga, d3fdraragb, d3fdraraab,
611 d3fdrarbrb, d3fdrarbga, d3fdrarbgb, d3fdrarbab, d3fdragaga,
612 d3fdragagb, d3fdragaab, d3fdragbgb, d3fdragbab, d3fdraabab,
613 d3fdrbrbrb, d3fdrbrbga, d3fdrbrbgb, d3fdrbrbab, d3fdrbgaga,
614 d3fdrbgagb, d3fdrbgaab, d3fdrbgbgb, d3fdrbgbab, d3fdrbabab,
615 d3fdgagaga, d3fdgagagb, d3fdgagaab, d3fdgagbgb, d3fdgagbab,
616 d3fdgaabab, d3fdgbgbgb, d3fdgbgbab, d3fdgbabab, d3fdababab;
617 real rhoa = dp->rhoa, rhob = dp->rhob;
618 real grada = dp->grada, gradb = dp->gradb, gradab = dp->gradab;
619
620 real t1, t2, t3, t4, t5, t6, t7, t8, t9, t10;
621 real t11, t12, t13, t14, t15, t16, t17, t18;
622 real t19, t20, t21, t22, t23, t24, t25, t26;
623 real t27, t28, t29, t30, t31, t32, t33, t34;
624 real t35, t36, t37, t38, t39, t40, t41, t42;
625 real t43, t44, t45, t46, t47, t48, t49, t50;
626 real t51, t52, t53, t54, t55, t56, t57, t58;
627 real t59, t60, t61, t62, t63, t64, t65, t66;
628 real t67, t68, t69, t70, t71, t72, t73, t74;
629 real t75, t76, t77, t78, t79, t80, t81, t82;
630 real t83, t84, t85, t86, t87, t88, t89, t90;
631 real t91, t92, t93, t94, t95, t96, t97, t98;
632 real t99, t100, t101, t102, t103, t104, t105;
633 real t106, t107, t108, t109, t110, t111, t112;
634 real t113, t114, t115, t116, t117, t118, t119;
635 real t120, t121, t122, t123, t124, t125, t126;
636 real t127, t128, t129, t130, t131, t132, t133;
637 real t134, t135, t136, t137, t138, t139, t140;
638 real t141, t142, t143, t144, t145, t146, t147;
639 real t148, t149, t150, t151, t152, t153, t154;
640 real t155, t156, t157, t158, t159, t160, t161;
641 real t162, t163, t164, t165, t166, t167, t168;
642 real t169, t170, t171, t172, t173, t174, t175;
643 real t176, t177, t178, t179, t180, t181, t182;
644 real t183, t184, t185, t186, t187, t188, t189;
645 real t190, t191, t192, t193, t194, t195, t196;
646 real t197, t198, t199, t200, t201, t202, t203;
647 real t204, t205, t206, t207, t208, t209, t210;
648 real t211, t212, t213, t214, t215;
649
650 t1 = pow(grada,2.0);
651 t2 = pow(gradb,2.0);
652 t3 = 2.0*gradab+t2+t1;
653 t4 = rhob+rhoa;
654 t5 = 1/pow(t4,2.333333333333333);
655 t6 = 0.318309886183791;
656 t7 = pow(3.0,0.666666666666667);
657 t8 = 1/pow(4.0,0.666666666666667);
658 t9 = pow(3.141592653589793,0.666666666666667);
659 t10 = 1/t9;
660 t11 = 1/pow(t4,0.666666666666667);
661 t12 = pow(3.0,0.333333333333333);
662 t13 = 1/pow(4.0,0.333333333333333);
663 t14 = pow(3.141592653589793,0.333333333333333);
664 t15 = 1/t14;
665 t16 = 1/pow(t4,0.333333333333333);
666 t17 = 0.0554175*t6/t4+8.723*t12*t13*t15*t16+0.472*t7*
667 t8*t10*t11+1.0;
668 t18 = 1/t17;
669 t19 = 23.266*t12*t13*t15*t16+0.007389*t7*t8*t10*t11+2.568;
670 t20 = 0.001*
671 t18*t19-0.001853571428571;
672 t21 = 1/t7;
673 t22 = pow(3.141592653589793,1.333333333333333);
674 t23 = 1/t22;
675 t24 = 1/pow(t4,2.666666666666667);
676 t25 = 1/pow(2.718281828459045,25.0*t21*t23*t24*t3);
677 t26 = 1/
678 pow(t4,1.666666666666667);
679 t27 = 1/pow(t4,1.333333333333333);
680 t28 = -7.755333333333333*t12*t13*t15*t27-0.004926*t7*
681 t8*t10*t26;
682 t29 = -0.0554175*t6/pow(t4,2.0)-2.907666666666667*t12*
683 t13*t15*t27-0.314666666666667*t7*t8*t10*t26;
684 t30 = 1/pow(t17,2.0);
685 t31 = 0.001*t18*t28-0.001*t19*t29*t30;
686 t32 = pow(t3,2.0);
687 t33 = 1/pow(t4,6.0);
688 t34 = 1/pow(t4,3.333333333333333);
689 t35 = 1/t12;
690 t36 = 1/pow(t4,4.666666666666667);
691 t37 = 0.2137*t12*t13*t15*t16+1.0;
692 t38 = 1.732050807568877;
693 t39 = 0.564189583547756;
694 t40 = sqrt(t4);
695 t41 = 1/pow(2.0,0.333333333333333);
696 t42 = pow(3.0,0.166666666666667);
697 t43 = 1/pow(3.141592653589793,0.166666666666667);
698 t44 = 7.5957*t41*t42*t43/pow(t4,0.166666666666667)+0.8191*
699 t38*t39/t40+3.5876*t12*t13*t15*t16+0.49294*t7*t8*t10*t11;
700 t45 = 16.0818243221511/
701 t44+1.0;
702 t46 = log(t45);
703 t47 = pow(2.718281828459045,2.43775691411733*t21*t9*t37*
704 t46);
705 t48 = t47-1.0;
706 t49 = 1/t48;
707 t50 = 0.0625*t14*t3*t35*t5+0.010866451862236*t32*t36*
708 t49;
709 t51 = 1/pow(t48,2.0);
710 t52 = 0.009188338695221*t35*t22*t32*t36*t51+0.166027154663518*
711 t21*t9*t3*t5*t49+1.0;
712 t53 = 1/t52;
713 t54 = 2.656434474616293*t35*t14*t50*t53+1.0;
714 t55 = 1/t54;
715 t56 = 1/pow(t4,5.666666666666667);
716 t57 = -1.26595*t41*t42*t43/pow(t4,1.166666666666667)-
717 0.40955*t38*t39/pow(t40,3.0)-1.195866666666667*t12*t13*t15*
718 t27-0.328626666666667*t7*t8*t10*t26;
719 t58 = 1/pow(t44,2.0);
720 t59 = 1/t45;
721 t60 = -0.173649550848958*t35*t13*t14*t27*t46-39.2035784329441*
722 t21*t9*t57*t37*t58*t59;
723 t61 = -0.010866451862236*t32*t36*t47*t51*t60-0.050710108690434*
724 t32*t49*t56-0.145833333333333*t14*t3*t34*t35;
725 t62 = 1/pow(t52,2.0);
726 t63 = 1/pow(t48,3.0);
727 t64 = -0.166027154663518*t21*t9*t3*t5*t60*t51*t47-0.018376677390443*
728 t35*t22*t32*t36*t60*t63*t47-0.387396694214876*t21*t9*t3*t34*
729 t49-0.042878913911033*t35*t22*t32*t56*t51;
730 t65 = 2.656434474616293*t35*t14*t53*t61-2.656434474616293*
731 t35*t14*t50*t62*t64;
732 t66 = t4*(0.025507875555556*t7*t10*t55*t65-2.333333333333333*
733 t20*t25*t3*t34+66.66666666666667*t20*t21*t23*t25*t32*t33+t3*
734 t5*t31*t25)+0.025507875555556*t7*t10*log(t54)+t3*t5*t20*t25;
735 t67 = 1/
736 pow(t4,5.0);
737 t68 = -50.0*t20*t21*t23*t25*t3*t67*grada;
738 t69 = 2.0*t20*t25*t5*grada;
739 t70 = 0.332054309327037*t21*t9*grada*t5*t49+0.036753354780885*
740 t35*t22*grada*t3*t36*t51;
741 t71 = 0.125*t14*t35*t5*grada+0.043465807448944*t3*t36*
742 t49*grada;
743 t72 = 2.656434474616293*t35*t14*t71*t53-2.656434474616293*
744 t35*t14*t50*t70*t62;
745 t73 = 0.025507875555556*t7*t10*t72*t55;
746 t74 = -50.0*t20*t21*t23*t25*t3*t67*gradb;
747 t75 = 2.0*t20*t25*t5*gradb;
748 t76 = 0.332054309327037*t21*t9*gradb*t5*t49+0.036753354780885*
749 t35*t22*gradb*t3*t36*t51;
750 t77 = 0.125*t14*t35*t5*gradb+0.043465807448944*t3*t36*
751 t49*gradb;
752 t78 = 2.656434474616293*t35*t14*t77*t53-2.656434474616293*
753 t35*t14*t50*t76*t62;
754 t79 = 0.025507875555556*t7*t10*t78*t55;
755 t80 = -50.0*t20*t21*t23*t25*t3*t67;
756 t81 = 2.0*t20*t25*t5;
757 t82 = 0.036753354780885*t35*t22*t3*t36*t51;
758 t83 = 0.332054309327037*t21*t9*t5*t49;
759 t84 = t83+t82;
760 t85 = 0.125*t14*t35*t5;
761 t86 = 0.043465807448944*t3*t36*t49;
762 t87 = t86+t85;
763 t88 = 2.656434474616293*t35*t14*t87*t53-2.656434474616293*
764 t35*t14*t50*t84*t62;
765 t89 = 0.025507875555556*t7*t10*t88*t55;
766 t90 = 10.34044444444444*t12*t13*t15*t5+0.00821*t7*t8*
767 t10*t24;
768 t91 = pow(t29,2.0);
769 t92 = 1/pow(t17,3.0);
770 t93 = 0.110835*t6/pow(t4,3.0)+3.876888888888889*t12*t13*
771 t15*t5+0.524444444444444*t7*t8*t10*t24;
772 t94 = -0.001*t19*t30*t93+0.002*t19*t91*t92+0.001*t18*
773 t90-0.002*t28*t29*t30;
774 t95 = 1/pow(3.141592653589793,2.666666666666667);
775 t96 = pow(t3,3.0);
776 t97 = 1/pow(t4,9.666666666666666);
777 t98 = 1/pow(t4,7.0);
778 t99 = 1/pow(t4,4.333333333333333);
779 t100 = 1/pow(t54,2.0);
780 t101 = pow(t65,2.0);
781 t102 = 1/pow(t52,3.0);
782 t103 = pow(t64,2.0);
783 t104 = 1/pow(t4,6.666666666666667);
784 t105 = pow(t57,2.0);
785 t106 = 1/pow(t44,4.0);
786 t107 = 1/pow(t45,2.0);
787 t108 = 1/pow(t44,3.0);
788 t109 = 1.594488888888889*t12*t13*t15*t5+1.476941666666667*
789 t41*t42*t43/pow(t4,2.166666666666667)+0.614325*t38*t39/pow(t40,
790 5.0)+0.547711111111111*t7*t8*t10*t24;
791 t110 = 0.231532734465277*t35*t13*t14*t5*t46-39.2035784329441*
792 t21*t9*t109*t37*t58*t59+5.58520314074677*t35*t13*t14*t27*t57*
793 t58*t59+78.40715686588821*t21*t9*t105*t37*t108*t59-630.4650611582789*
794 t21*t9*t105*t37*t106*t107;
795 t111 = pow(t60,2.0);
796 t112 = pow(2.718281828459045,4.875513828234661*t21*t9*
797 t37*t46);
798 t113 = 0.486111111111111*t14*t3*t35*t99+0.021732903724472*
799 t111*t112*t32*t36*t63+0.101420217380869*t32*t47*t51*t56*t60-
800 0.010866451862236*t111*t32*t36*t47*t51-0.010866451862236*t110*
801 t32*t36*t47*t51+0.287357282579128*t104*t32*t49;
802 t114 = 1/pow(t48,4.0);
803 t115 = 0.332054309327037*t21*t9*t3*t5*t111*t63*t112+0.055130032171328*
804 t35*t22*t32*t36*t111*t114*t112-0.166027154663518*t21*t9*t3*
805 t5*t111*t51*t47+0.774793388429752*t21*t9*t3*t34*t60*t51*t47-
806 0.166027154663518*t21*t9*t3*t5*t110*t51*t47-0.018376677390443*
807 t35*t22*t32*t36*t111*t63*t47+0.17151565564413*t35*t22*t32*
808 t56*t60*t63*t47-0.018376677390443*t35*t22*t32*t36*t110*t63*
809 t47+1.291322314049587*t21*t9*t3*t99*t49+0.242980512162518*
810 t35*t22*t32*t104*t51;
811 t116 = -2.656434474616293*t35*t14*t50*t62*t115+2.656434474616293*
812 t35*t14*t53*t113+5.312868949232586*t35*t14*t50*t102*t103-5.312868949232586*
813 t35*t14*t62*t61*t64;
814 t117 = t4*(7.777777777777778*t20*t25*t3*t99-555.5555555555555*
815 t20*t21*t23*t25*t32*t98+1481.481481481482*t20*t25*t35*t95*
816 t96*t97-4.666666666666667*t25*t3*t31*t34+133.3333333333333*
817 t21*t23*t25*t31*t32*t33+t3*t5*t94*t25+0.025507875555556*t7*
818 t10*t55*t116-0.025507875555556*t7*t10*t100*t101)+0.051015751111111*
819 t7*t10*t55*t65+2.0*t25*t3*t31*t5-4.666666666666667*t20*t25*
820 t3*t34+133.3333333333333*t20*t21*t23*t25*t32*t33;
821 t118 = 1/pow(t4,8.666666666666666);
822 t119 = -0.043465807448944*t3*t36*t47*t51*t60*grada-0.202840434761738*
823 t3*t49*t56*grada-0.291666666666667*t14*t34*t35*grada;
824 t120 = -0.332054309327037*t21*t9*grada*t5*t60*t51*t47-
825 0.07350670956177*t35*t22*grada*t3*t36*t60*t63*t47-0.774793388429752*
826 t21*t9*grada*t34*t49-0.17151565564413*t35*t22*grada*t3*t56*
827 t51;
828 t121 = -2.656434474616293*t35*t14*t71*t62*t64+5.312868949232586*
829 t35*t14*t50*t70*t102*t64-2.656434474616293*t35*t14*t50*t62*
830 t120-2.656434474616293*t35*t14*t70*t62*t61+2.656434474616293*
831 t35*t14*t53*t119;
832 t122 = t4*(-1111.111111111111*t118*t20*t25*t32*t35*t95*
833 grada-50.0*t21*t23*t25*t3*t31*t67*grada+2.0*t25*t31*t5*grada-
834 4.666666666666667*t20*t25*t34*grada+383.3333333333333*t20*
835 t21*t23*t25*t3*t33*grada-0.025507875555556*t7*t10*t72*t100*
836 t65+0.025507875555556*t7*t10*t55*t121)+t73+t69+t68;
837 t123 = -0.043465807448944*t3*t36*t47*t51*t60*gradb-0.202840434761738*
838 t3*t49*t56*gradb-0.291666666666667*t14*t34*t35*gradb;
839 t124 = -0.332054309327037*t21*t9*gradb*t5*t60*t51*t47-
840 0.07350670956177*t35*t22*gradb*t3*t36*t60*t63*t47-0.774793388429752*
841 t21*t9*gradb*t34*t49-0.17151565564413*t35*t22*gradb*t3*t56*
842 t51;
843 t125 = -2.656434474616293*t35*t14*t77*t62*t64+5.312868949232586*
844 t35*t14*t50*t76*t102*t64-2.656434474616293*t35*t14*t50*t62*
845 t124-2.656434474616293*t35*t14*t76*t62*t61+2.656434474616293*
846 t35*t14*t53*t123;
847 t126 = t4*(-1111.111111111111*t118*t20*t25*t32*t35*t95*
848 gradb-50.0*t21*t23*t25*t3*t31*t67*gradb+2.0*t25*t31*t5*gradb-
849 4.666666666666667*t20*t25*t34*gradb+383.3333333333333*t20*
850 t21*t23*t25*t3*t33*gradb-0.025507875555556*t7*t10*t78*t100*
851 t65+0.025507875555556*t7*t10*t55*t125)+t79+t75+t74;
852 t127 = -50.0*t21*t23*t25*t3*t31*t67;
853 t128 = 2.0*t25*t31*t5;
854 t129 = -1111.111111111111*t118*t20*t25*t32*t35*t95;
855 t130 = 383.3333333333333*
856 t20*t21*t23*t25*t3*t33;
857 t131 = -4.666666666666667*t20*t25*t34;
858 t132 = -0.291666666666667*t14*t34*t35;
859 t133 = -0.202840434761738*t3*t49*t56;
860 t134 = -0.043465807448944*t3*t36*t47*t51*t60;
861 t135 = t134+t133+t132;
862 t136 = -0.17151565564413*t35*t22*t3*t56*t51;
863 t137 = -0.774793388429752*t21*t9*t34*t49;
864 t138 = -0.07350670956177*t35*t22*t3*t36*t60*t63*t47;
865 t139 = -
866 0.332054309327037*t21*t9*t5*t60*t51*t47;
867 t140 = t139+t138+t137+t136;
868 t141 = -2.656434474616293*t35*t14*t87*t62*t64+5.312868949232586*
869 t35*t14*t50*t84*t102*t64-2.656434474616293*t35*t14*t50*t62*
870 t140-2.656434474616293*t35*t14*t84*t62*t61+2.656434474616293*
871 t35*t14*t53*t135;
872 t142 = t4*(-0.025507875555556*t7*t10*t88*t100*t65+0.025507875555556*
873 t7*t10*t55*t141+t131+t130+t129+t128+t127)+t89+t81+t80;
874 t143 = 1/pow(t4,7.666666666666667);
875 t144 = 833.3333333333334*t1*t143*t20*t25*t3*t35*t95;
876 t145 = -
877 200.0*t1*t20*t21*t23*t25*t67;
878 t146 = pow(t72,2.0);
879 t147 = -0.025507875555556*t7*t10*t146*t100;
880 t148 = pow(t70,2.0);
881 t149 = t83+t82+0.07350670956177*t35*t22*t1*t36*t51;
882 t150 = t86+
883 t85+0.086931614897888*t1*t36*t49;
884 t151 = 2.656434474616293*t35*t14*t150*t53-5.312868949232586*
885 t35*t14*t71*t70*t62-2.656434474616293*t35*t14*t50*t149*t62+
886 5.312868949232586*t35*t14*t50*t148*t102;
887 t152 = 0.025507875555556*t7*t10*t151*t55;
888 t153 = 833.3333333333334*t143*t20*t25*t3*t35*t95*grada*
889 gradb;
890 t154 = -200.0*t20*t21*t23*t25*t67*grada*gradb;
891 t155 = -0.025507875555556*t7*t10*t72*t78*t100;
892 t156 = pow(3.141592653589793,1.666666666666667);
893 t157 = 0.07350670956177*t35*t22*grada*gradb*t36*t53*t49-
894 0.195265757395493*t21*t156*grada*gradb*t36*t50*t62*t51-2.656434474616293*
895 t35*t14*t71*t76*t62-2.656434474616293*t35*t14*t77*t70*t62+
896 5.312868949232586*t35*t14*t50*t70*t76*t102;
897 t158 = 0.025507875555556*t7*t10*t55*t157;
898 t159 = 833.3333333333334*t143*t20*t25*t3*t35*t95*grada;
899 t160 = -
900 200.0*t20*t21*t23*t25*t67*grada;
901 t161 = -0.025507875555556*t7*t10*t88*t72*t100;
902 t162 = -0.195265757395493*t21*t156*grada*t36*t50*t62*
903 t51;
904 t163 = 0.07350670956177*t35*t22*grada*t36*t53*t49;
905 t164 = t163+t162-2.656434474616293*t35*t14*t87*t70*t62-
906 2.656434474616293*t35*t14*t71*t84*t62+5.312868949232586*t35*
907 t14*t50*t84*t70*t102;
908 t165 = 0.025507875555556*t7*t10*t55*t164;
909 t166 = 833.3333333333334*t143*t2*t20*t25*t3*t35*t95;
910 t167 = -
911 200.0*t2*t20*t21*t23*t25*t67;
912 t168 = pow(t78,2.0);
913 t169 = -0.025507875555556*t7*t10*t168*t100;
914 t170 = pow(t76,2.0);
915 t171 = t83+t82+0.07350670956177*t35*t22*t2*t36*t51;
916 t172 = t86+
917 t85+0.086931614897888*t2*t36*t49;
918 t173 = 2.656434474616293*t35*t14*t172*t53-5.312868949232586*
919 t35*t14*t77*t76*t62-2.656434474616293*t35*t14*t50*t171*t62+
920 5.312868949232586*t35*t14*t50*t170*t102;
921 t174 = 0.025507875555556*t7*t10*t173*t55;
922 t175 = 833.3333333333334*t143*t20*t25*t3*t35*t95*gradb;
923 t176 = -
924 200.0*t20*t21*t23*t25*t67*gradb;
925 t177 = -0.025507875555556*t7*t10*t88*t78*t100;
926 t178 = -0.195265757395493*t21*t156*gradb*t36*t50*t62*
927 t51;
928 t179 = 0.07350670956177*t35*t22*gradb*t36*t53*t49;
929 t180 = t179+t178-2.656434474616293*t35*t14*t87*t76*t62-
930 2.656434474616293*t35*t14*t77*t84*t62+5.312868949232586*t35*
931 t14*t50*t84*t76*t102;
932 t181 = 0.025507875555556*t7*t10*t55*t180;
933 t182 = 833.3333333333334*t143*t20*t25*t3*t35*t95;
934 t183 = -200.0*t20*t21*t23*t25*t67;
935 t184 = pow(t88,2.0);
936 t185 = -0.025507875555556*t7*t10*t184*t100;
937 t186 = pow(t84,2.0);
938 t187 = -0.195265757395493*t21*t156*t36*t50*t62*t51;
939 t188 = 0.07350670956177*
940 t35*t22*t36*t53*t49;
941 t189 = t188+t187-5.312868949232586*t35*t14*t87*t84*t62+
942 5.312868949232586*t35*t14*t50*t186*t102;
943 t190 = 0.025507875555556*t7*t10*t55*t189;
944 t191 = 1/pow(t4,3.666666666666667);
945 t192 = 1/pow(3.141592653589793,4.0);
946 t193 = 1/pow(t4,5.333333333333333);
947 t194 = 1/pow(t54,3.0);
948 t195 = 1/pow(t52,4.0);
949 t196 = pow(t57,3.0);
950 t197 = -39.2035784329441*t21*t37*(-3.200040277777777*
951 t41*t42*t43/pow(t4,3.166666666666667)-1.5358125*t38*t39/pow(t40,
952 7.0)-3.720474074074074*t12*t13*t15*t34-1.460562962962963*t7*
953 t8*t10*t191)*t58*t59*t9-20278.05670960338*t196*t21*t37*t9/
954 (pow(t44,6.0)*pow(t45,3.0))+3782.790366949674*t107*t196*t21*
955 t37*t9/pow(t44,5.0)-11.17040628149354*t35*t13*t14*t5*t57*t58*
956 t59+8.377804711120154*t35*t13*t14*t27*t109*t58*t59+235.2214705976646*
957 t21*t9*t109*t57*t37*t108*t59-16.75560942224031*t35*t13*t14*
958 t27*t105*t108*t59-235.2214705976646*t21*t9*t196*t37*t106*t59-
959 0.540243047085647*t35*t13*t14*t34*t46-1891.395183474837*t21*
960 t9*t109*t57*t37*t106*t107+134.7303835695242*t35*t13*t14*t27*
961 t105*t106*t107;
962 t198 = pow(t60,3.0);
963 t199 = -0.347299101697916*t35*t13*t14*t27*t46-78.40715686588821*
964 t21*t9*t57*t37*t58*t59;
965 t200 = pow(2.718281828459045,7.313270742351992*t21*t9*
966 t37*t46);
967 t201 = t4*(23.33333333333333*t25*t3*t31*t99-1666.666666666667*
968 t21*t23*t25*t31*t32*t98+4444.444444444444*t25*t31*t35*t95*
969 t96*t97-26666.66666666667*t20*t25*t35*t95*t96/pow(t4,10.66666666666667)-
970 7.0*t25*t3*t34*t94+200.0*t21*t23*t25*t32*t33*t94+t25*t3*t5*
971 (0.006*t19*t29*t92*t93-0.003*t28*t30*t93+0.006*t28*t91*t92-
972 0.003*t29*t30*t90-0.001*t19*t30*(-0.332505*t6/pow(t4,4.0)-
973 9.046074074074074*t12*t13*t15*t34-1.398518518518518*t7*t8*
974 t10*t191)+0.001*t18*(-24.1277037037037*t12*t13*t15*t34-0.021893333333333*
975 t7*t8*t10*t191)-0.006*t19*pow(t29,3.0)/pow(t17,4.0))+0.051015751111111*
976 t10*t194*pow(t65,3.0)*t7+0.025507875555556*t10*t55*(-15.93860684769776*
977 t14*t195*t35*t50*pow(t64,3.0)+2.656434474616293*t14*t35*t53*
978 (0.065198711173416*t110*t112*t32*t36*t60*t63-0.304260652142607*
979 t111*t112*t32*t56*t63+0.021732903724472*t111*t112*t199*t32*
980 t36*t63+0.021732903724472*t112*t198*t32*t36*t63-0.032599355586708*
981 t110*t32*t36*t47*t51*t60-0.862071847737385*t104*t32*t47*t51*
982 t60+0.152130326071303*t111*t32*t47*t51*t56+0.152130326071303*
983 t110*t32*t47*t51*t56-0.010866451862236*t198*t32*t36*t47*t51-
984 0.010866451862236*t197*t32*t36*t47*t51-1.91571521719419*t143*
985 t32*t49-0.065198711173416*t114*t198*t200*t32*t36-2.106481481481481*
986 t14*t193*t3*t35)-2.656434474616293*t14*t35*t50*(-1.619870081083455*
987 t35*t22*t32*t143*t51-5.59573002754821*t21*t9*t3*t193*t49-0.220520128685311*
988 t198*t200*t22*t32*t35*t36/pow(t48,5.0)-0.055130032171328*t35*
989 t22*t32*t36*t110*t60*t63*t47-1.457883072975109*t35*t22*t32*
990 t104*t60*t63*t47-0.018376677390443*t35*t22*t32*t36*t198*t63*
991 t47-0.018376677390443*t35*t22*t32*t36*t197*t63*t47+0.257273483466196*
992 t35*t22*t32*t56*t111*t63*t47+0.257273483466196*t35*t22*t32*
993 t56*t110*t63*t47-3.87396694214876*t21*t9*t3*t99*t60*t51*t47-
994 0.498081463990555*t21*t9*t3*t5*t110*t60*t51*t47-0.166027154663518*
995 t21*t9*t3*t5*t198*t51*t47-0.166027154663518*t21*t9*t3*t5*t197*
996 t51*t47+1.162190082644628*t21*t9*t3*t34*t111*t51*t47+1.162190082644628*
997 t21*t9*t3*t34*t110*t51*t47-0.99616292798111*t21*t9*t3*t5*t198*
998 t114*t200+0.99616292798111*t21*t9*t3*t5*t110*t60*t63*t112+
999 0.332054309327037*t21*t9*t3*t5*t198*t63*t112-2.324380165289256*
1000 t21*t9*t3*t34*t111*t63*t112+0.332054309327037*t21*t9*t3*t5*
1001 t199*t111*t63*t112+0.165390096513983*t35*t22*t32*t36*t110*
1002 t60*t114*t112+0.055130032171328*t35*t22*t32*t36*t198*t114*
1003 t112-0.771820450398587*t35*t22*t32*t56*t111*t114*t112+0.055130032171328*
1004 t35*t22*t32*t36*t199*t111*t114*t112)*t62+15.93860684769776*
1005 t35*t14*t50*t102*t64*t115-7.969303423848879*t35*t14*t62*t61*
1006 t115-7.969303423848879*t35*t14*t62*t64*t113+15.93860684769776*
1007 t35*t14*t102*t61*t103)*t7+32921.81069958848*t192*t20*t25*pow(t3,
1008 4.0)/pow(t4,13.33333333333333)+4407.407407407408*t20*t21*t23*
1009 t25*t32/pow(t4,8.0)-33.7037037037037*t193*t20*t25*t3-0.076523626666667*
1010 t7*t10*t100*t65*t116)+23.33333333333333*t20*t25*t3*t99-1666.666666666667*
1011 t20*t21*t23*t25*t32*t98+4444.444444444444*t20*t25*t35*t95*
1012 t96*t97+3.0*t25*t3*t5*t94-14.0*t25*t3*t31*t34+400.0*t21*t23*
1013 t25*t31*t32*t33+0.076523626666667*t7*t10*t55*t116-0.076523626666667*
1014 t7*t10*t100*t101;
1015 t202 = 1/pow(t4,12.33333333333333);
1016 t203 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1017 t14*t35*t53*(0.972222222222222*t14*t35*t99*grada+0.086931614897888*
1018 t111*t112*t3*t36*t63*grada+0.405680869523475*t3*t47*t51*t56*
1019 t60*grada-0.043465807448944*t111*t3*t36*t47*t51*grada-0.043465807448944*
1020 t110*t3*t36*t47*t51*grada+1.149429130316514*t104*t3*t49*grada)+
1021 10.62573789846517*t35*t14*t70*t102*t61*t64+10.62573789846517*
1022 t35*t14*t50*t102*t120*t64-5.312868949232586*t35*t14*t62*t119*
1023 t64-5.312868949232586*t35*t14*t62*t61*t120-2.656434474616293*
1024 t35*t14*t71*t62*t115+5.312868949232586*t35*t14*t50*t70*t102*
1025 t115-2.656434474616293*t35*t14*t70*t62*t113-2.656434474616293*
1026 t35*t14*t50*t62*(0.664108618654073*t21*t9*grada*t5*t111*t63*
1027 t112+0.220520128685311*t35*t22*grada*t3*t36*t111*t114*t112-
1028 0.332054309327037*t21*t9*grada*t5*t111*t51*t47+1.549586776859504*
1029 t21*t9*grada*t34*t60*t51*t47-0.332054309327037*t21*t9*grada*
1030 t5*t110*t51*t47-0.07350670956177*t35*t22*grada*t3*t36*t111*
1031 t63*t47+0.686062622576522*t35*t22*grada*t3*t56*t60*t63*t47-
1032 0.07350670956177*t35*t22*grada*t3*t36*t110*t63*t47+2.582644628099174*
1033 t21*t9*grada*t99*t49+0.971922048650073*t35*t22*grada*t3*t104*
1034 t51)-15.93860684769776*t35*t14*t50*t70*t195*t103+5.312868949232586*
1035 t35*t14*t71*t102*t103)+15.55555555555556*t20*t25*t99*grada-
1036 2611.111111111111*t20*t21*t23*t25*t3*t98*grada+18148.14814814815*
1037 t20*t25*t32*t35*t95*t97*grada-24691.35802469136*t192*t20*t202*
1038 t25*t96*grada-2222.222222222222*t118*t25*t31*t32*t35*t95*grada-
1039 50.0*t21*t23*t25*t3*t67*t94*grada+2.0*t25*t5*t94*grada-9.333333333333334*
1040 t25*t31*t34*grada+766.6666666666666*t21*t23*t25*t3*t31*t33*
1041 grada-0.051015751111111*t7*t10*t100*t121*t65-0.025507875555556*
1042 t7*t10*t72*t100*t116+0.051015751111111*t7*t10*t72*t194*t101)-
1043 2222.222222222222*t118*t20*t25*t32*t35*t95*grada-100.0*t21*
1044 t23*t25*t3*t31*t67*grada+4.0*t25*t31*t5*grada-9.333333333333334*
1045 t20*t25*t34*grada+766.6666666666666*t20*t21*t23*t25*t3*t33*
1046 grada-0.051015751111111*t7*t10*t72*t100*t65+0.051015751111111*
1047 t7*t10*t55*t121;
1048 t204 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1049 t14*t35*t53*(0.972222222222222*t14*t35*t99*gradb+0.086931614897888*
1050 t111*t112*t3*t36*t63*gradb+0.405680869523475*t3*t47*t51*t56*
1051 t60*gradb-0.043465807448944*t111*t3*t36*t47*t51*gradb-0.043465807448944*
1052 t110*t3*t36*t47*t51*gradb+1.149429130316514*t104*t3*t49*gradb)+
1053 10.62573789846517*t35*t14*t76*t102*t61*t64+10.62573789846517*
1054 t35*t14*t50*t102*t124*t64-5.312868949232586*t35*t14*t62*t123*
1055 t64-5.312868949232586*t35*t14*t62*t61*t124-2.656434474616293*
1056 t35*t14*t77*t62*t115+5.312868949232586*t35*t14*t50*t76*t102*
1057 t115-2.656434474616293*t35*t14*t76*t62*t113-2.656434474616293*
1058 t35*t14*t50*t62*(0.664108618654073*t21*t9*gradb*t5*t111*t63*
1059 t112+0.220520128685311*t35*t22*gradb*t3*t36*t111*t114*t112-
1060 0.332054309327037*t21*t9*gradb*t5*t111*t51*t47+1.549586776859504*
1061 t21*t9*gradb*t34*t60*t51*t47-0.332054309327037*t21*t9*gradb*
1062 t5*t110*t51*t47-0.07350670956177*t35*t22*gradb*t3*t36*t111*
1063 t63*t47+0.686062622576522*t35*t22*gradb*t3*t56*t60*t63*t47-
1064 0.07350670956177*t35*t22*gradb*t3*t36*t110*t63*t47+2.582644628099174*
1065 t21*t9*gradb*t99*t49+0.971922048650073*t35*t22*gradb*t3*t104*
1066 t51)-15.93860684769776*t35*t14*t50*t76*t195*t103+5.312868949232586*
1067 t35*t14*t77*t102*t103)+15.55555555555556*t20*t25*t99*gradb-
1068 2611.111111111111*t20*t21*t23*t25*t3*t98*gradb+18148.14814814815*
1069 t20*t25*t32*t35*t95*t97*gradb-24691.35802469136*t192*t20*t202*
1070 t25*t96*gradb-2222.222222222222*t118*t25*t31*t32*t35*t95*gradb-
1071 50.0*t21*t23*t25*t3*t67*t94*gradb+2.0*t25*t5*t94*gradb-9.333333333333334*
1072 t25*t31*t34*gradb+766.6666666666666*t21*t23*t25*t3*t31*t33*
1073 gradb-0.051015751111111*t7*t10*t100*t125*t65-0.025507875555556*
1074 t7*t10*t78*t100*t116+0.051015751111111*t7*t10*t78*t194*t101)-
1075 2222.222222222222*t118*t20*t25*t32*t35*t95*gradb-100.0*t21*
1076 t23*t25*t3*t31*t67*gradb+4.0*t25*t31*t5*gradb-9.333333333333334*
1077 t20*t25*t34*gradb+766.6666666666666*t20*t21*t23*t25*t3*t33*
1078 gradb-0.051015751111111*t7*t10*t78*t100*t65+0.051015751111111*
1079 t7*t10*t55*t125;
1080 t205 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1081 t14*t35*t53*(0.972222222222222*t14*t35*t99+0.086931614897888*
1082 t111*t112*t3*t36*t63+0.405680869523475*t3*t47*t51*t56*t60-
1083 0.043465807448944*t111*t3*t36*t47*t51-0.043465807448944*t110*
1084 t3*t36*t47*t51+1.149429130316514*t104*t3*t49)+10.62573789846517*
1085 t35*t14*t84*t102*t61*t64+10.62573789846517*t35*t14*t50*t102*
1086 t140*t64-5.312868949232586*t35*t14*t62*t135*t64-5.312868949232586*
1087 t35*t14*t62*t61*t140-2.656434474616293*t35*t14*t87*t62*t115+
1088 5.312868949232586*t35*t14*t50*t84*t102*t115-2.656434474616293*
1089 t35*t14*t84*t62*t113-2.656434474616293*t35*t14*t50*t62*(0.664108618654073*
1090 t21*t9*t5*t111*t63*t112+0.220520128685311*t35*t22*t3*t36*t111*
1091 t114*t112-0.332054309327037*t21*t9*t5*t111*t51*t47+1.549586776859504*
1092 t21*t9*t34*t60*t51*t47-0.332054309327037*t21*t9*t5*t110*t51*
1093 t47-0.07350670956177*t35*t22*t3*t36*t111*t63*t47+0.686062622576522*
1094 t35*t22*t3*t56*t60*t63*t47-0.07350670956177*t35*t22*t3*t36*
1095 t110*t63*t47+2.582644628099174*t21*t9*t99*t49+0.971922048650073*
1096 t35*t22*t3*t104*t51)-15.93860684769776*t35*t14*t50*t84*t195*
1097 t103+5.312868949232586*t35*t14*t87*t102*t103)+15.55555555555556*
1098 t20*t25*t99-2611.111111111111*t20*t21*t23*t25*t3*t98+18148.14814814815*
1099 t20*t25*t32*t35*t95*t97-24691.35802469136*t192*t20*t202*t25*
1100 t96-2222.222222222222*t118*t25*t31*t32*t35*t95-50.0*t21*t23*
1101 t25*t3*t67*t94+2.0*t25*t5*t94-0.051015751111111*t7*t10*t100*
1102 t141*t65-9.333333333333334*t25*t31*t34+766.6666666666666*t21*
1103 t23*t25*t3*t31*t33-0.025507875555556*t7*t10*t88*t100*t116+
1104 0.051015751111111*t7*t10*t88*t194*t101)-2222.222222222222*
1105 t118*t20*t25*t32*t35*t95-100.0*t21*t23*t25*t3*t31*t67-0.051015751111111*
1106 t7*t10*t88*t100*t65+4.0*t25*t31*t5-9.333333333333334*t20*t25*
1107 t34+766.6666666666666*t20*t21*t23*t25*t3*t33+0.051015751111111*
1108 t7*t10*t55*t141;
1109 t206 = 1/pow(t4,11.33333333333333);
1110 t207 = t4*(833.3333333333334*t1*t143*t25*t3*t31*t35*t95-
1111 10833.33333333333*t1*t118*t20*t25*t3*t35*t95+0.025507875555556*
1112 t10*t55*(-2.656434474616293*t35*t14*t150*t62*t64-15.93860684769776*
1113 t35*t14*t50*t148*t195*t64+10.62573789846517*t35*t14*t71*t70*
1114 t102*t64+5.312868949232586*t35*t14*t50*t149*t102*t64-2.656434474616293*
1115 t35*t14*t149*t62*t61+5.312868949232586*t35*t14*t148*t102*t61+
1116 2.656434474616293*t14*t35*t53*(-0.086931614897888*t1*t36*t47*
1117 t51*t60-0.405680869523475*t1*t49*t56+t134+t133+t132)-2.656434474616293*
1118 t35*t14*t50*t62*(t139+t138-0.14701341912354*t35*t22*t1*t36*
1119 t60*t63*t47+t137+t136-0.343031311288261*t35*t22*t1*t56*t51)-
1120 5.312868949232586*t35*t14*t71*t62*t120+10.62573789846517*t35*
1121 t14*t50*t70*t102*t120-5.312868949232586*t35*t14*t70*t62*t119)*
1122 t7-200.0*t1*t21*t23*t25*t31*t67+0.051015751111111*t7*t10*t146*
1123 t194*t65-0.025507875555556*t7*t10*t151*t100*t65+1000.0*t1*
1124 t20*t21*t23*t25*t33+18518.51851851852*t1*t192*t20*t206*t25*
1125 t32+t131+t130+t129+t128+t127-0.051015751111111*t7*t10*t72*
1126 t100*t121)+t81+t80+t152+t147+t145+t144;
1127 t208 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1128 t14*t35*t53*(-0.086931614897888*t36*t47*t51*t60*grada*gradb-
1129 0.405680869523475*t49*t56*grada*gradb)+0.390531514790987*t21*
1130 t156*grada*gradb*t36*t50*t102*t51*t64-0.07350670956177*t35*
1131 t22*grada*gradb*t36*t62*t49*t64-15.93860684769776*t35*t14*
1132 t50*t70*t76*t195*t64+5.312868949232586*t35*t14*t71*t76*t102*
1133 t64+5.312868949232586*t35*t14*t77*t70*t102*t64-0.195265757395493*
1134 t21*t156*grada*gradb*t36*t62*t51*t61+5.312868949232586*t35*
1135 t14*t70*t76*t102*t61-2.656434474616293*t35*t14*t50*t62*(-0.14701341912354*
1136 t35*t22*grada*gradb*t36*t60*t63*t47-0.343031311288261*t35*
1137 t22*grada*gradb*t56*t51)-2.656434474616293*t35*t14*t71*t62*
1138 t124+5.312868949232586*t35*t14*t50*t70*t102*t124-2.656434474616293*
1139 t35*t14*t70*t62*t123-2.656434474616293*t35*t14*t77*t62*t120+
1140 5.312868949232586*t35*t14*t50*t76*t102*t120-2.656434474616293*
1141 t35*t14*t76*t62*t119)+833.3333333333334*t143*t25*t3*t31*t35*
1142 t95*grada*gradb-10833.33333333333*t118*t20*t25*t3*t35*t95*
1143 grada*gradb-200.0*t21*t23*t25*t31*t67*grada*gradb+1000.0*t20*
1144 t21*t23*t25*t33*grada*gradb+18518.51851851852*t192*t20*t206*
1145 t25*t32*grada*gradb+0.051015751111111*t7*t10*t72*t78*t194*
1146 t65-0.025507875555556*t7*t10*t100*t157*t65-0.025507875555556*
1147 t7*t10*t72*t100*t125-0.025507875555556*t7*t10*t78*t100*t121)+
1148 t158+t155+t154+t153;
1149 t209 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1150 t14*t35*t53*(-0.086931614897888*t36*t47*t51*t60*grada-0.405680869523475*
1151 t49*t56*grada)+0.390531514790987*t21*t156*grada*t36*t50*t102*
1152 t51*t64-0.07350670956177*t35*t22*grada*t36*t62*t49*t64-15.93860684769776*
1153 t35*t14*t50*t84*t70*t195*t64+5.312868949232586*t35*t14*t71*
1154 t84*t102*t64+5.312868949232586*t35*t14*t87*t70*t102*t64-0.195265757395493*
1155 t21*t156*grada*t36*t62*t51*t61+5.312868949232586*t35*t14*t84*
1156 t70*t102*t61-2.656434474616293*t35*t14*t50*t62*(-0.14701341912354*
1157 t35*t22*grada*t36*t60*t63*t47-0.343031311288261*t35*t22*grada*
1158 t56*t51)-2.656434474616293*t35*t14*t71*t62*t140+5.312868949232586*
1159 t35*t14*t50*t70*t102*t140-2.656434474616293*t35*t14*t70*t62*
1160 t135-2.656434474616293*t35*t14*t87*t62*t120+5.312868949232586*
1161 t35*t14*t50*t84*t102*t120-2.656434474616293*t35*t14*t84*t62*
1162 t119)+833.3333333333334*t143*t25*t3*t31*t35*t95*grada-10833.33333333333*
1163 t118*t20*t25*t3*t35*t95*grada-200.0*t21*t23*t25*t31*t67*grada+
1164 1000.0*t20*t21*t23*t25*t33*grada+18518.51851851852*t192*t20*
1165 t206*t25*t32*grada+0.051015751111111*t7*t10*t88*t72*t194*t65-
1166 0.025507875555556*t7*t10*t100*t164*t65-0.025507875555556*t7*
1167 t10*t72*t100*t141-0.025507875555556*t7*t10*t88*t100*t121)+
1168 t165+t161+t160+t159;
1169 t210 = t4*(833.3333333333334*t143*t2*t25*t3*t31*t35*t95-
1170 10833.33333333333*t118*t2*t20*t25*t3*t35*t95+0.025507875555556*
1171 t10*t55*(-2.656434474616293*t35*t14*t172*t62*t64-15.93860684769776*
1172 t35*t14*t50*t170*t195*t64+10.62573789846517*t35*t14*t77*t76*
1173 t102*t64+5.312868949232586*t35*t14*t50*t171*t102*t64-2.656434474616293*
1174 t35*t14*t171*t62*t61+5.312868949232586*t35*t14*t170*t102*t61+
1175 2.656434474616293*t14*t35*t53*(-0.086931614897888*t2*t36*t47*
1176 t51*t60-0.405680869523475*t2*t49*t56+t134+t133+t132)-2.656434474616293*
1177 t35*t14*t50*t62*(t139+t138-0.14701341912354*t35*t22*t2*t36*
1178 t60*t63*t47+t137+t136-0.343031311288261*t35*t22*t2*t56*t51)-
1179 5.312868949232586*t35*t14*t77*t62*t124+10.62573789846517*t35*
1180 t14*t50*t76*t102*t124-5.312868949232586*t35*t14*t76*t62*t123)*
1181 t7-200.0*t2*t21*t23*t25*t31*t67+0.051015751111111*t7*t10*t168*
1182 t194*t65-0.025507875555556*t7*t10*t173*t100*t65+1000.0*t2*
1183 t20*t21*t23*t25*t33+18518.51851851852*t192*t2*t20*t206*t25*
1184 t32+t131+t130+t129+t128+t127-0.051015751111111*t7*t10*t78*
1185 t100*t125)+t81+t80+t174+t169+t167+t166;
1186 t211 = t4*(0.025507875555556*t10*t55*t7*(2.656434474616293*
1187 t14*t35*t53*(-0.086931614897888*t36*t47*t51*t60*gradb-0.405680869523475*
1188 t49*t56*gradb)+0.390531514790987*t21*t156*gradb*t36*t50*t102*
1189 t51*t64-0.07350670956177*t35*t22*gradb*t36*t62*t49*t64-15.93860684769776*
1190 t35*t14*t50*t84*t76*t195*t64+5.312868949232586*t35*t14*t77*
1191 t84*t102*t64+5.312868949232586*t35*t14*t87*t76*t102*t64-0.195265757395493*
1192 t21*t156*gradb*t36*t62*t51*t61+5.312868949232586*t35*t14*t84*
1193 t76*t102*t61-2.656434474616293*t35*t14*t50*t62*(-0.14701341912354*
1194 t35*t22*gradb*t36*t60*t63*t47-0.343031311288261*t35*t22*gradb*
1195 t56*t51)-2.656434474616293*t35*t14*t77*t62*t140+5.312868949232586*
1196 t35*t14*t50*t76*t102*t140-2.656434474616293*t35*t14*t76*t62*
1197 t135-2.656434474616293*t35*t14*t87*t62*t124+5.312868949232586*
1198 t35*t14*t50*t84*t102*t124-2.656434474616293*t35*t14*t84*t62*
1199 t123)+833.3333333333334*t143*t25*t3*t31*t35*t95*gradb-10833.33333333333*
1200 t118*t20*t25*t3*t35*t95*gradb-200.0*t21*t23*t25*t31*t67*gradb+
1201 1000.0*t20*t21*t23*t25*t33*gradb+18518.51851851852*t192*t20*
1202 t206*t25*t32*gradb+0.051015751111111*t7*t10*t88*t78*t194*t65-
1203 0.025507875555556*t7*t10*t100*t180*t65-0.025507875555556*t7*
1204 t10*t78*t100*t141-0.025507875555556*t7*t10*t88*t100*t125)+
1205 t181+t177+t176+t175;
1206 t212 = t4*(833.3333333333334*t143*t25*t3*t31*t35*t95-
1207 10833.33333333333*t118*t20*t25*t3*t35*t95+0.025507875555556*
1208 t10*t55*(0.390531514790987*t21*t156*t36*t50*t102*t51*t64-0.07350670956177*
1209 t35*t22*t36*t62*t49*t64-15.93860684769776*t35*t14*t50*t186*
1210 t195*t64+10.62573789846517*t35*t14*t87*t84*t102*t64-0.195265757395493*
1211 t21*t156*t36*t62*t51*t61+5.312868949232586*t35*t14*t186*t102*
1212 t61+2.656434474616293*t14*t35*t53*(-0.086931614897888*t36*
1213 t47*t51*t60-0.405680869523475*t49*t56)-2.656434474616293*t35*
1214 t14*t50*t62*(-0.14701341912354*t35*t22*t36*t60*t63*t47-0.343031311288261*
1215 t35*t22*t56*t51)-5.312868949232586*t35*t14*t87*t62*t140+10.62573789846517*
1216 t35*t14*t50*t84*t102*t140-5.312868949232586*t35*t14*t84*t62*
1217 t135)*t7-200.0*t21*t23*t25*t31*t67+0.051015751111111*t7*t10*
1218 t184*t194*t65-0.025507875555556*t7*t10*t100*t189*t65+1000.0*
1219 t20*t21*t23*t25*t33+18518.51851851852*t192*t20*t206*t25*t32-
1220 0.051015751111111*t7*t10*t88*t100*t141)+t190+t185+t183+t182;
1221 t213 = pow(grada,
1222 3.0);
1223 t214 = 1/pow(t4,10.33333333333333);
1224 t215 = pow(gradb,3.0);
1225
1226 /* code */
1227 dfdra = t66;
1228 dfdrb = t66;
1229 dfdga = t4*(t73+t69+t68);
1230 dfdgb = t4*(t79+t75+t74);
1231 dfdgab = t4*(t89+t81+t80);
1232 d2fdrara = t117;
1233 d2fdrarb = t117;
1234 d2fdraga = t122;
1235 d2fdragb = t126;
1236 d2fdraab = t142;
1237 d2fdrbrb = t117;
1238 d2fdrbga = t122;
1239 d2fdrbgb = t126;
1240 d2fdrbgab = t142;
1241 d2fdgaga = t4*(t152+t147+t81+t80+t145+t144);
1242 d2fdgagb = t4*(t158+t155+t154+t153);
1243 d2fdgagab = t4*(t165+t161+t160+t159);
1244 d2fdgbgb = t4*(t174+t169+t81+t80+t167+t166);
1245 d2fdgbgab = t4*(t181+t177+t176+t175);
1246 d2fdgabgab = t4*(t190+t185+t183+t182);
1247 d3fdrarara = t201;
1248 d3fdrararb = t201;
1249 d3fdraraga = t203;
1250 d3fdraragb = t204;
1251 d3fdraraab = t205;
1252 d3fdrarbrb = t201;
1253 d3fdrarbga = t203;
1254 d3fdrarbgb = t204;
1255 d3fdrarbab = t205;
1256 d3fdragaga = t207;
1257 d3fdragagb = t208;
1258 d3fdragaab = t209;
1259 d3fdragbgb = t210;
1260 d3fdragbab = t211;
1261 d3fdraabab = t212;
1262 d3fdrbrbrb = t201;
1263 d3fdrbrbga = t203;
1264 d3fdrbrbgb = t204;
1265 d3fdrbrbab = t205;
1266 d3fdrbgaga = t207;
1267 d3fdrbgagb = t208;
1268 d3fdrbgaab = t209;
1269 d3fdrbgbgb = t210;
1270 d3fdrbgbab = t211;
1271 d3fdrbabab = t212;
1272 d3fdgagaga = t4*(2500.0*t143*t20*t25*t3*t35*t95*grada-
1273 600.0*t20*t21*t23*t25*t67*grada+5000.0*t143*t20*t213*t25*t35*
1274 t95+0.051015751111111*t10*t194*t7*pow(t72,3.0)+0.025507875555556*
1275 t10*t55*t7*(-15.93860684769776*t14*t195*t35*t50*pow(t70,3.0)-
1276 7.969303423848879*t35*t14*t150*t70*t62-7.969303423848879*t35*
1277 t14*t71*t149*t62-0.58579727218648*t21*t156*grada*t36*t50*t62*
1278 t51+0.220520128685311*t35*t22*grada*t36*t53*t49+15.93860684769776*
1279 t35*t14*t50*t149*t70*t102+15.93860684769776*t35*t14*t71*t148*
1280 t102)-13888.88888888889*t192*t20*t213*t214*t25*t3-0.076523626666667*
1281 t7*t10*t151*t72*t100);
1282 d3fdgagagb = t4*(5000.0*t1*t143*t20*t25*t35*t95*gradb-
1283 13888.88888888889*t1*t192*t20*t214*t25*t3*gradb+0.051015751111111*
1284 t7*t10*t146*t78*t194+0.025507875555556*t7*t10*t55*(t179-0.14701341912354*
1285 t35*t22*grada*gradb*t36*t70*t62*t49+t178-0.390531514790987*
1286 t21*t156*grada*gradb*t36*t71*t62*t51+0.781063029581974*t21*
1287 t156*grada*gradb*t36*t50*t70*t102*t51-2.656434474616293*t35*
1288 t14*t150*t76*t62-2.656434474616293*t35*t14*t77*t149*t62+10.62573789846517*
1289 t35*t14*t71*t70*t76*t102+5.312868949232586*t35*t14*t50*t149*
1290 t76*t102+5.312868949232586*t35*t14*t77*t148*t102-15.93860684769776*
1291 t35*t14*t50*t148*t76*t195)+t176+t175-0.051015751111111*t7*
1292 t10*t72*t100*t157-0.025507875555556*t7*t10*t151*t78*t100);
1293 d3fdgagaab = t4*
1294 (5000.0*t1*t143*t20*t25*t35*t95-13888.88888888889*t1*t192*
1295 t20*t214*t25*t3+0.051015751111111*t7*t10*t88*t146*t194+0.025507875555556*
1296 t7*t10*t55*(t188-0.14701341912354*t35*t22*grada*t36*t70*t62*
1297 t49+t187-0.390531514790987*t21*t156*grada*t36*t71*t62*t51+
1298 0.781063029581974*t21*t156*grada*t36*t50*t70*t102*t51-2.656434474616293*
1299 t35*t14*t87*t149*t62-2.656434474616293*t35*t14*t150*t84*t62+
1300 5.312868949232586*t35*t14*t87*t148*t102+10.62573789846517*
1301 t35*t14*t71*t84*t70*t102+5.312868949232586*t35*t14*t50*t84*
1302 t149*t102-15.93860684769776*t35*t14*t50*t84*t148*t195)+t183+
1303 t182-0.051015751111111*t7*t10*t72*t100*t164-0.025507875555556*
1304 t7*t10*t88*t151*t100);
1305 d3fdgagbgb = t4*(5000.0*t143*t2*t20*t25*t35*t95*grada-
1306 13888.88888888889*t192*t2*t20*t214*t25*t3*grada+0.051015751111111*
1307 t7*t10*t72*t168*t194+0.025507875555556*t7*t10*t55*(t163-0.14701341912354*
1308 t35*t22*grada*gradb*t36*t76*t62*t49+t162-0.390531514790987*
1309 t21*t156*grada*gradb*t36*t77*t62*t51+0.781063029581974*t21*
1310 t156*grada*gradb*t36*t50*t76*t102*t51-2.656434474616293*t35*
1311 t14*t172*t70*t62-2.656434474616293*t35*t14*t71*t171*t62+5.312868949232586*
1312 t35*t14*t71*t170*t102+10.62573789846517*t35*t14*t77*t70*t76*
1313 t102+5.312868949232586*t35*t14*t50*t171*t70*t102-15.93860684769776*
1314 t35*t14*t50*t70*t170*t195)+t160+t159-0.051015751111111*t7*
1315 t10*t78*t100*t157-0.025507875555556*t7*t10*t173*t72*t100);
1316 d3fdgagbab = t4*
1317 (5000.0*t143*t20*t25*t35*t95*grada*gradb-13888.88888888889*
1318 t192*t20*t214*t25*t3*grada*gradb+0.025507875555556*t7*t10*
1319 t55*(-0.07350670956177*t35*t22*grada*t36*t76*t62*t49-0.07350670956177*
1320 t35*t22*gradb*t36*t70*t62*t49-0.07350670956177*t35*t22*grada*
1321 gradb*t36*t84*t62*t49-0.195265757395493*t21*t156*grada*t36*
1322 t77*t62*t51-0.195265757395493*t21*t156*gradb*t36*t71*t62*t51-
1323 0.195265757395493*t21*t156*grada*gradb*t36*t87*t62*t51+0.390531514790987*
1324 t21*t156*grada*t36*t50*t76*t102*t51+0.390531514790987*t21*
1325 t156*gradb*t36*t50*t70*t102*t51+0.390531514790987*t21*t156*
1326 grada*gradb*t36*t50*t84*t102*t51+5.312868949232586*t35*t14*
1327 t87*t70*t76*t102+5.312868949232586*t35*t14*t71*t84*t76*t102+
1328 5.312868949232586*t35*t14*t77*t84*t70*t102-15.93860684769776*
1329 t35*t14*t50*t84*t70*t76*t195)+0.051015751111111*t7*t10*t88*
1330 t72*t78*t194-0.025507875555556*t7*t10*t72*t100*t180-0.025507875555556*
1331 t7*t10*t78*t100*t164-0.025507875555556*t7*t10*t88*t100*t157);
1332 d3fdgaabab = t4*
1333 (5000.0*t143*t20*t25*t35*t95*grada-13888.88888888889*t192*
1334 t20*t214*t25*t3*grada+0.025507875555556*t7*t10*t55*(-0.07350670956177*
1335 t35*t22*t36*t70*t62*t49-0.14701341912354*t35*t22*grada*t36*
1336 t84*t62*t49-0.195265757395493*t21*t156*t36*t71*t62*t51-0.390531514790987*
1337 t21*t156*grada*t36*t87*t62*t51+0.390531514790987*t21*t156*
1338 t36*t50*t70*t102*t51+0.781063029581974*t21*t156*grada*t36*
1339 t50*t84*t102*t51+10.62573789846517*t35*t14*t87*t84*t70*t102+
1340 5.312868949232586*t35*t14*t71*t186*t102-15.93860684769776*
1341 t35*t14*t50*t186*t70*t195)+0.051015751111111*t7*t10*t184*t72*
1342 t194-0.025507875555556*t7*t10*t72*t100*t189-0.051015751111111*
1343 t7*t10*t88*t100*t164);
1344 d3fdgbgbgb = t4*(2500.0*t143*t20*t25*t3*t35*t95*gradb-
1345 600.0*t20*t21*t23*t25*t67*gradb+5000.0*t143*t20*t215*t25*t35*
1346 t95+0.051015751111111*t10*t194*t7*pow(t78,3.0)+0.025507875555556*
1347 t10*t55*t7*(-15.93860684769776*t14*t195*t35*t50*pow(t76,3.0)-
1348 7.969303423848879*t35*t14*t172*t76*t62-7.969303423848879*t35*
1349 t14*t77*t171*t62-0.58579727218648*t21*t156*gradb*t36*t50*t62*
1350 t51+0.220520128685311*t35*t22*gradb*t36*t53*t49+15.93860684769776*
1351 t35*t14*t50*t171*t76*t102+15.93860684769776*t35*t14*t77*t170*
1352 t102)-13888.88888888889*t192*t20*t214*t215*t25*t3-0.076523626666667*
1353 t7*t10*t173*t78*t100);
1354 d3fdgbgbab = t4*(5000.0*t143*t2*t20*t25*t35*t95-13888.88888888889*
1355 t192*t2*t20*t214*t25*t3+0.051015751111111*t7*t10*t88*t168*
1356 t194+0.025507875555556*t7*t10*t55*(t188-0.14701341912354*t35*
1357 t22*gradb*t36*t76*t62*t49+t187-0.390531514790987*t21*t156*
1358 gradb*t36*t77*t62*t51+0.781063029581974*t21*t156*gradb*t36*
1359 t50*t76*t102*t51-2.656434474616293*t35*t14*t87*t171*t62-2.656434474616293*
1360 t35*t14*t172*t84*t62+5.312868949232586*t35*t14*t87*t170*t102+
1361 10.62573789846517*t35*t14*t77*t84*t76*t102+5.312868949232586*
1362 t35*t14*t50*t84*t171*t102-15.93860684769776*t35*t14*t50*t84*
1363 t170*t195)+t183+t182-0.051015751111111*t7*t10*t78*t100*t180-
1364 0.025507875555556*t7*t10*t88*t173*t100);
1365 d3fdgbabab = t4*(5000.0*t143*t20*t25*t35*t95*gradb-13888.88888888889*
1366 t192*t20*t214*t25*t3*gradb+0.025507875555556*t7*t10*t55*(-
1367 0.07350670956177*t35*t22*t36*t76*t62*t49-0.14701341912354*
1368 t35*t22*gradb*t36*t84*t62*t49-0.195265757395493*t21*t156*t36*
1369 t77*t62*t51-0.390531514790987*t21*t156*gradb*t36*t87*t62*t51+
1370 0.390531514790987*t21*t156*t36*t50*t76*t102*t51+0.781063029581974*
1371 t21*t156*gradb*t36*t50*t84*t102*t51+10.62573789846517*t35*
1372 t14*t87*t84*t76*t102+5.312868949232586*t35*t14*t77*t186*t102-
1373 15.93860684769776*t35*t14*t50*t186*t76*t195)+0.051015751111111*
1374 t7*t10*t184*t78*t194-0.025507875555556*t7*t10*t78*t100*t189-
1375 0.051015751111111*t7*t10*t88*t100*t180);
1376 d3fdababab = t4*(5000.0*t143*t20*t25*t35*t95+0.051015751111111*
1377 t10*t194*t7*pow(t88,3.0)+0.025507875555556*t10*t55*t7*(-15.93860684769776*
1378 t14*t195*t35*t50*pow(t84,3.0)-0.58579727218648*t21*t156*t36*
1379 t87*t62*t51+1.171594544372961*t21*t156*t36*t50*t84*t102*t51-
1380 0.220520128685311*t35*t22*t36*t84*t62*t49+15.93860684769776*
1381 t35*t14*t87*t186*t102)-13888.88888888889*t192*t20*t214*t25*
1382 t3-0.076523626666667*t7*t10*t88*t100*t189);
1383
1384
1385 ds->df1000 += factor*dfdra;
1386 ds->df0100 += factor*dfdrb;
1387 ds->df0010 += factor*dfdga;
1388 ds->df0001 += factor*dfdgb;
1389 ds->df00001+= factor*dfdgab;
1390
1391 ds->df2000 += factor*d2fdrara;
1392 ds->df1100 += factor*d2fdrarb;
1393 ds->df1010 += factor*d2fdraga;
1394 ds->df1001 += factor*d2fdragb;
1395 ds->df10001+= factor*d2fdraab;
1396 ds->df0200 += factor*d2fdrbrb;
1397 ds->df0110 += factor*d2fdrbga;
1398 ds->df0101 += factor*d2fdrbgb;
1399 ds->df01001+= factor*d2fdrbgab;
1400 ds->df0020 += factor*d2fdgaga;
1401 ds->df0011 += factor*d2fdgagb;
1402 ds->df00101+= factor*d2fdgagab;
1403 ds->df0002 += factor*d2fdgbgb;
1404 ds->df00011+= factor*d2fdgbgab;
1405 ds->df00002+= factor*d2fdgabgab;
1406
1407 ds->df3000 += factor*d3fdrarara;
1408 ds->df2100 += factor*d3fdrararb;
1409 ds->df2010 += factor*d3fdraraga;
1410 ds->df2001 += factor*d3fdraragb;
1411 ds->df20001 += factor*d3fdraraab;
1412 ds->df1200 += factor*d3fdrarbrb;
1413 ds->df1110 += factor*d3fdrarbga;
1414 ds->df1101 += factor*d3fdrarbgb;
1415 ds->df11001 += factor*d3fdrarbab;
1416 ds->df1020 += factor*d3fdragaga;
1417 ds->df1011 += factor*d3fdragagb;
1418 ds->df10101 += factor*d3fdragaab;
1419 ds->df1002 += factor*d3fdragbgb;
1420 ds->df10011 += factor*d3fdragbab;
1421 ds->df10002 += factor*d3fdraabab;
1422 ds->df0300 += factor*d3fdrbrbrb;
1423 ds->df0210 += factor*d3fdrbrbga;
1424 ds->df0201 += factor*d3fdrbrbgb;
1425 ds->df02001 += factor*d3fdrbrbab;
1426 ds->df0120 += factor*d3fdrbgaga;
1427 ds->df0111 += factor*d3fdrbgagb;
1428 ds->df01101 += factor*d3fdrbgaab;
1429 ds->df0102 += factor*d3fdrbgbgb;
1430 ds->df01011 += factor*d3fdrbgbab;
1431 ds->df01002 += factor*d3fdrbabab;
1432 ds->df0030 += factor*d3fdgagaga;
1433 ds->df0021 += factor*d3fdgagagb;
1434 ds->df00201 += factor*d3fdgagaab;
1435 ds->df0012 += factor*d3fdgagbgb;
1436 ds->df00111 += factor*d3fdgagbab;
1437 ds->df00102 += factor*d3fdgaabab;
1438 ds->df0003 += factor*d3fdgbgbgb;
1439 ds->df00021 += factor*d3fdgbgbab;
1440 ds->df00012 += factor*d3fdgbabab;
1441 ds->df00003 += factor*d3fdababab;
1442
1443 }
1444