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