1 /*
2   This file was generated automatically with ./scripts/maple2c.pl.
3   Do not edit this file directly as it can be overwritten!!
4 
5   This Source Code Form is subject to the terms of the Mozilla Public
6   License, v. 2.0. If a copy of the MPL was not distributed with this
7   file, You can obtain one at http://mozilla.org/MPL/2.0/.
8 
9   Maple version     : Maple 2016 (X86 64 LINUX)
10   Maple source      : ./maple/gga_exc/gga_x_ol2.mpl
11   Type of functional: gga_exc
12 */
13 
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16 
17 
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21 
22 #ifndef XC_DONT_COMPILE_EXC
23   double t1, t3, t5, t6, t7, t8, t9, t10;
24   double t11, t12, t13, t14, t16, t17, t20, t21;
25   double t22, t23, t28, t29, t30, t32, t34;
26 
27 #ifndef XC_DONT_COMPILE_VXC
28   double t38, t39, t40, t42, t43, t49, t52, t53;
29   double t54, t55, t58, t62, t65, t66, t69, t74;
30 
31 #ifndef XC_DONT_COMPILE_FXC
32   double t77, t78, t85, t87, t88, t94, t97, t100;
33   double t101, t102, t103, t105, t106, t109, t114, t124;
34   double t125, t127, t131, t134, t135, t138, t139, t140;
35   double t142, t148;
36 
37 #ifndef XC_DONT_COMPILE_KXC
38   double t152, t153, t164, t165, t171, t177, t178, t182;
39   double t183, t184, t186, t187, t188, t193, t198, t215;
40   double t217, t221, t230, t235, t236, t240, t243, t245;
41   double t248, t249, t253, t254, t256, t260;
42 
43 #ifndef XC_DONT_COMPILE_LXC
44   double t275, t276, t290, t302, t387;
45 #endif
46 
47 #endif
48 
49 #endif
50 
51 #endif
52 
53 #endif
54 
55 
56   gga_x_ol2_params *params;
57 
58   assert(p->params != NULL);
59   params = (gga_x_ol2_params * )(p->params);
60 
61   t1 = M_CBRT3;
62   t3 = POW_1_3(0.1e1 / M_PI);
63   t5 = M_CBRT4;
64   t6 = t5 * t5;
65   t7 = t1 * t3 * t6;
66   t8 = M_CBRT2;
67   t9 = t8 * t8;
68   t10 = POW_1_3(rho[0]);
69   t11 = t9 * t10;
70   t12 = params->bb * sigma[0];
71   t13 = rho[0] * rho[0];
72   t14 = t10 * t10;
73   t16 = 0.1e1 / t14 / t13;
74   t17 = t9 * t16;
75   t20 = sqrt(sigma[0]);
76   t21 = params->cc * t20;
77   t22 = t10 * rho[0];
78   t23 = 0.1e1 / t22;
79   t28 = 0.4e1 * t20 * t8 * t23 + t8;
80   t29 = 0.1e1 / t28;
81   t30 = t8 * t23 * t29;
82   t32 = params->aa + 0.13888888888888888889e-1 * t12 * t17 + t21 * t30;
83   t34 = t7 * t11 * t32;
84   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
85     zk[0] = -0.3e1 / 0.16e2 * t34;
86 
87 #ifndef XC_DONT_COMPILE_VXC
88 
89   if(order < 1) return;
90 
91 
92   t38 = t22 * t1 * t3;
93   t39 = t6 * t9;
94   t40 = t13 * rho[0];
95   t42 = 0.1e1 / t14 / t40;
96   t43 = t9 * t42;
97   t49 = t8 / t10 / t13 * t29;
98   t52 = params->cc * sigma[0];
99   t53 = t28 * t28;
100   t54 = 0.1e1 / t53;
101   t55 = t43 * t54;
102   t58 = -0.37037037037037037037e-1 * t12 * t43 - 0.4e1 / 0.3e1 * t21 * t49 + 0.16e2 / 0.3e1 * t52 * t55;
103   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
104     vrho[0] = -t34 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t39 * t58;
105 
106   t62 = params->bb * t9;
107   t65 = 0.1e1 / t20;
108   t66 = params->cc * t65;
109   t69 = params->cc * t9;
110   t74 = t39 * (0.13888888888888888889e-1 * t62 * t16 + t66 * t30 / 0.2e1 - 0.2e1 * t69 * t16 * t54);
111   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
112     vsigma[0] = -0.3e1 / 0.16e2 * t38 * t74;
113 
114 #ifndef XC_DONT_COMPILE_FXC
115 
116   if(order < 2) return;
117 
118 
119   t77 = 0.1e1 / t14;
120   t78 = t9 * t77;
121   t85 = t13 * t13;
122   t87 = 0.1e1 / t14 / t85;
123   t88 = t9 * t87;
124   t94 = t8 / t10 / t40 * t29;
125   t97 = t88 * t54;
126   t100 = t20 * sigma[0];
127   t101 = params->cc * t100;
128   t102 = t85 * t13;
129   t103 = 0.1e1 / t102;
130   t105 = 0.1e1 / t53 / t28;
131   t106 = t103 * t105;
132   t109 = 0.13580246913580246914e0 * t12 * t88 + 0.28e2 / 0.9e1 * t21 * t94 - 0.80e2 / 0.3e1 * t52 * t97 + 0.1024e4 / 0.9e1 * t101 * t106;
133   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
134     v2rho2[0] = -t7 * t78 * t32 / 0.12e2 - t7 * t11 * t58 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t39 * t109;
135 
136   t114 = t10 * t1 * t3;
137   t124 = t85 * rho[0];
138   t125 = 0.1e1 / t124;
139   t127 = t105 * t20;
140   t131 = t39 * (-0.37037037037037037037e-1 * t62 * t42 - 0.2e1 / 0.3e1 * t66 * t49 + 0.8e1 * t69 * t42 * t54 - 0.128e3 / 0.3e1 * params->cc * t125 * t127);
141   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
142     v2rhosigma[0] = -t114 * t74 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t131;
143 
144   t134 = 0.1e1 / t100;
145   t135 = params->cc * t134;
146   t138 = 0.1e1 / sigma[0];
147   t139 = params->cc * t138;
148   t140 = t17 * t54;
149   t142 = 0.1e1 / t85;
150   t148 = t39 * (-t135 * t30 / 0.4e1 - t139 * t140 + 0.16e2 * params->cc * t142 * t105 * t65);
151   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
152     v2sigma2[0] = -0.3e1 / 0.16e2 * t38 * t148;
153 
154 #ifndef XC_DONT_COMPILE_KXC
155 
156   if(order < 3) return;
157 
158 
159   t152 = 0.1e1 / t14 / rho[0];
160   t153 = t9 * t152;
161   t164 = 0.1e1 / t14 / t124;
162   t165 = t9 * t164;
163   t171 = t8 / t10 / t85 * t29;
164   t177 = t85 * t40;
165   t178 = 0.1e1 / t177;
166   t182 = sigma[0] * sigma[0];
167   t183 = params->cc * t182;
168   t184 = t85 * t85;
169   t186 = 0.1e1 / t10 / t184;
170   t187 = t53 * t53;
171   t188 = 0.1e1 / t187;
172   t193 = -0.63374485596707818932e0 * t12 * t165 - 0.280e3 / 0.27e2 * t21 * t171 + 0.3808e4 / 0.27e2 * t52 * t165 * t54 - 0.11264e5 / 0.9e1 * t101 * t178 * t105 + 0.16384e5 / 0.9e1 * t183 * t186 * t188 * t8;
173   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
174     v3rho3[0] = t7 * t153 * t32 / 0.18e2 - t7 * t78 * t58 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t109 - 0.3e1 / 0.16e2 * t38 * t39 * t193;
175 
176   t198 = t77 * t1 * t3;
177   t215 = params->cc / t10 / t177;
178   t217 = t188 * sigma[0] * t8;
179   t221 = t39 * (0.13580246913580246914e0 * t62 * t87 + 0.14e2 / 0.9e1 * t66 * t94 - 0.296e3 / 0.9e1 * t69 * t87 * t54 + 0.384e3 * params->cc * t103 * t127 - 0.2048e4 / 0.3e1 * t215 * t217);
180   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
181     v3rho2sigma[0] = -t198 * t74 / 0.12e2 - t114 * t131 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t221;
182 
183   t230 = t125 * t105;
184   t235 = params->cc / t10 / t102;
185   t236 = t188 * t8;
186   t240 = t39 * (t135 * t49 / 0.3e1 + 0.4e1 / 0.3e1 * t139 * t55 - 0.256e3 / 0.3e1 * t66 * t230 + 0.256e3 * t235 * t236);
187   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
188     v3rhosigma2[0] = -t114 * t148 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t240;
189 
190   t243 = t20 * t182;
191   t245 = params->cc / t243;
192   t248 = 0.1e1 / t182;
193   t249 = params->cc * t248;
194   t253 = 0.1e1 / t10 / t124;
195   t254 = params->cc * t253;
196   t256 = t188 * t138 * t8;
197   t260 = t39 * (0.3e1 / 0.8e1 * t245 * t30 + 0.3e1 / 0.2e1 * t249 * t140 - 0.96e2 * t254 * t256);
198   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
199     v3sigma3[0] = -0.3e1 / 0.16e2 * t38 * t260;
200 
201 #ifndef XC_DONT_COMPILE_LXC
202 
203   if(order < 4) return;
204 
205 
206   t275 = 0.1e1 / t14 / t102;
207   t276 = t9 * t275;
208   t290 = t184 * rho[0];
209   t302 = 0.1e1 / t187 / t28;
210   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
211     v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t17 * t32 + 0.2e1 / 0.9e1 * t7 * t153 * t58 - t7 * t78 * t109 / 0.2e1 - t7 * t11 * t193 - 0.3e1 / 0.16e2 * t38 * t39 * (0.35912208504801097395e1 * t12 * t276 + 0.3640e4 / 0.81e2 * t21 * t8 * t253 * t29 - 0.23072e5 / 0.27e2 * t52 * t276 * t54 + 0.953344e6 / 0.81e2 * t101 / t184 * t105 - 0.950272e6 / 0.27e2 * t183 / t10 / t290 * t188 * t8 + 0.1048576e7 / 0.27e2 * params->cc * t243 / t14 / t184 / t13 * t302 * t9);
212 
213   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
214     v4rho3sigma[0] = t152 * t1 * t3 * t74 / 0.18e2 - t198 * t131 / 0.4e1 - 0.3e1 / 0.4e1 * t114 * t221 - 0.3e1 / 0.16e2 * t38 * t39 * (-0.63374485596707818932e0 * t62 * t164 - 0.140e3 / 0.27e2 * t66 * t171 + 0.1456e4 / 0.9e1 * t69 * t164 * t54 - 0.81152e5 / 0.27e2 * params->cc * t178 * t127 + 0.100352e6 / 0.9e1 * params->cc * t186 * t217 - 0.131072e6 / 0.9e1 * params->cc / t14 / t290 * t302 * t100 * t9);
215 
216   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
217     v4rho2sigma2[0] = -t198 * t148 / 0.12e2 - t114 * t240 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t39 * (-0.7e1 / 0.9e1 * t135 * t94 - 0.28e2 / 0.9e1 * t139 * t97 + 0.4096e4 / 0.9e1 * t66 * t106 - 0.8960e4 / 0.3e1 * t215 * t236 + 0.16384e5 / 0.3e1 * params->cc / t14 / t184 * t302 * t9 * t20);
218 
219   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
220     v4rhosigma3[0] = -t114 * t260 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t39 * (-t245 * t49 / 0.2e1 - 0.2e1 * t249 * t55 + 0.32e2 * t135 * t230 + 0.512e3 * t235 * t256 - 0.2048e4 * params->cc / t14 / t177 * t302 * t65 * t9);
221 
222   t387 = t182 * sigma[0];
223   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
224     v4sigma4[0] = -0.3e1 / 0.16e2 * t38 * t39 * (-0.15e2 / 0.16e2 * params->cc / t20 / t387 * t30 - 0.15e2 / 0.4e1 * params->cc / t387 * t140 - 0.12e2 * t245 * t142 * t105 + 0.768e3 * params->cc * t275 * t302 * t134 * t9 + 0.96e2 * t254 * t188 * t248 * t8);
225 
226 #ifndef XC_DONT_COMPILE_MXC
227 
228   if(order < 5) return;
229 
230 
231 #endif
232 
233 #endif
234 
235 #endif
236 
237 #endif
238 
239 #endif
240 
241 
242 }
243 
244 
245 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))246 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
247 {
248 
249 #ifndef XC_DONT_COMPILE_EXC
250   double t1, t3, t4, t5, t6, t7, t8, t9;
251   double t10, t11, t13, t16, t17, t18, t19, t20;
252   double t23, t24, t25, t27, t29;
253 
254 #ifndef XC_DONT_COMPILE_VXC
255   double t32, t33, t34, t36, t41, t44, t45, t46;
256   double t47, t50, t56, t57, t64;
257 
258 #ifndef XC_DONT_COMPILE_FXC
259   double t67, t68, t74, t76, t81, t84, t87, t88;
260   double t89, t90, t92, t93, t96, t100, t110, t111;
261   double t113, t117, t120, t121, t124, t125, t126, t128;
262   double t134;
263 
264 #ifndef XC_DONT_COMPILE_KXC
265   double t138, t139, t150, t155, t161, t162, t166, t167;
266   double t168, t170, t171, t172, t176, t180, t196, t197;
267   double t201, t210, t215, t219, t222, t224, t227, t228;
268   double t232, t233, t234, t238;
269 
270 #ifndef XC_DONT_COMPILE_LXC
271   double t254, t267, t278, t357;
272 #endif
273 
274 #endif
275 
276 #endif
277 
278 #endif
279 
280 #endif
281 
282 
283   gga_x_ol2_params *params;
284 
285   assert(p->params != NULL);
286   params = (gga_x_ol2_params * )(p->params);
287 
288   t1 = M_CBRT3;
289   t3 = POW_1_3(0.1e1 / M_PI);
290   t4 = t1 * t3;
291   t5 = M_CBRT4;
292   t6 = t5 * t5;
293   t7 = POW_1_3(rho[0]);
294   t8 = t6 * t7;
295   t9 = params->bb * sigma[0];
296   t10 = rho[0] * rho[0];
297   t11 = t7 * t7;
298   t13 = 0.1e1 / t11 / t10;
299   t16 = sqrt(sigma[0]);
300   t17 = params->cc * t16;
301   t18 = t7 * rho[0];
302   t19 = 0.1e1 / t18;
303   t20 = M_CBRT2;
304   t23 = 0.4e1 * t16 * t19 + t20;
305   t24 = 0.1e1 / t23;
306   t25 = t19 * t24;
307   t27 = params->aa + 0.13888888888888888889e-1 * t9 * t13 + t17 * t25;
308   t29 = t4 * t8 * t27;
309   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
310     zk[0] = -0.3e1 / 0.8e1 * t29;
311 
312 #ifndef XC_DONT_COMPILE_VXC
313 
314   if(order < 1) return;
315 
316 
317   t32 = t18 * t1;
318   t33 = t3 * t6;
319   t34 = t10 * rho[0];
320   t36 = 0.1e1 / t11 / t34;
321   t41 = 0.1e1 / t7 / t10 * t24;
322   t44 = params->cc * sigma[0];
323   t45 = t23 * t23;
324   t46 = 0.1e1 / t45;
325   t47 = t36 * t46;
326   t50 = -0.37037037037037037037e-1 * t9 * t36 - 0.4e1 / 0.3e1 * t17 * t41 + 0.16e2 / 0.3e1 * t44 * t47;
327   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
328     vrho[0] = -t29 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t33 * t50;
329 
330   t56 = 0.1e1 / t16;
331   t57 = params->cc * t56;
332   t64 = t33 * (0.13888888888888888889e-1 * params->bb * t13 + t57 * t25 / 0.2e1 - 0.2e1 * params->cc * t13 * t46);
333   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
334     vsigma[0] = -0.3e1 / 0.8e1 * t32 * t64;
335 
336 #ifndef XC_DONT_COMPILE_FXC
337 
338   if(order < 2) return;
339 
340 
341   t67 = 0.1e1 / t11;
342   t68 = t6 * t67;
343   t74 = t10 * t10;
344   t76 = 0.1e1 / t11 / t74;
345   t81 = 0.1e1 / t7 / t34 * t24;
346   t84 = t76 * t46;
347   t87 = t16 * sigma[0];
348   t88 = params->cc * t87;
349   t89 = t74 * t10;
350   t90 = 0.1e1 / t89;
351   t92 = 0.1e1 / t45 / t23;
352   t93 = t90 * t92;
353   t96 = 0.13580246913580246914e0 * t9 * t76 + 0.28e2 / 0.9e1 * t17 * t81 - 0.80e2 / 0.3e1 * t44 * t84 + 0.512e3 / 0.9e1 * t88 * t93;
354   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
355     v2rho2[0] = -t4 * t68 * t27 / 0.6e1 - t4 * t8 * t50 - 0.3e1 / 0.8e1 * t32 * t33 * t96;
356 
357   t100 = t7 * t1;
358   t110 = t74 * rho[0];
359   t111 = 0.1e1 / t110;
360   t113 = t92 * t16;
361   t117 = t33 * (-0.37037037037037037037e-1 * params->bb * t36 - 0.2e1 / 0.3e1 * t57 * t41 + 0.8e1 * params->cc * t36 * t46 - 0.64e2 / 0.3e1 * params->cc * t111 * t113);
362   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
363     v2rhosigma[0] = -t100 * t64 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t117;
364 
365   t120 = 0.1e1 / t87;
366   t121 = params->cc * t120;
367   t124 = 0.1e1 / sigma[0];
368   t125 = params->cc * t124;
369   t126 = t13 * t46;
370   t128 = 0.1e1 / t74;
371   t134 = t33 * (-t121 * t25 / 0.4e1 - t125 * t126 + 0.8e1 * params->cc * t128 * t92 * t56);
372   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
373     v2sigma2[0] = -0.3e1 / 0.8e1 * t32 * t134;
374 
375 #ifndef XC_DONT_COMPILE_KXC
376 
377   if(order < 3) return;
378 
379 
380   t138 = 0.1e1 / t11 / rho[0];
381   t139 = t6 * t138;
382   t150 = 0.1e1 / t11 / t110;
383   t155 = 0.1e1 / t7 / t74 * t24;
384   t161 = t74 * t34;
385   t162 = 0.1e1 / t161;
386   t166 = sigma[0] * sigma[0];
387   t167 = params->cc * t166;
388   t168 = t74 * t74;
389   t170 = 0.1e1 / t7 / t168;
390   t171 = t45 * t45;
391   t172 = 0.1e1 / t171;
392   t176 = -0.63374485596707818932e0 * t9 * t150 - 0.280e3 / 0.27e2 * t17 * t155 + 0.3808e4 / 0.27e2 * t44 * t150 * t46 - 0.5632e4 / 0.9e1 * t88 * t162 * t92 + 0.8192e4 / 0.9e1 * t167 * t170 * t172;
393   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
394     v3rho3[0] = t4 * t139 * t27 / 0.9e1 - t4 * t68 * t50 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t96 - 0.3e1 / 0.8e1 * t32 * t33 * t176;
395 
396   t180 = t67 * t1;
397   t196 = params->cc / t7 / t161;
398   t197 = t172 * sigma[0];
399   t201 = t33 * (0.13580246913580246914e0 * params->bb * t76 + 0.14e2 / 0.9e1 * t57 * t81 - 0.296e3 / 0.9e1 * params->cc * t76 * t46 + 0.192e3 * params->cc * t90 * t113 - 0.1024e4 / 0.3e1 * t196 * t197);
400   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
401     v3rho2sigma[0] = -t180 * t64 / 0.6e1 - t100 * t117 - 0.3e1 / 0.8e1 * t32 * t201;
402 
403   t210 = t111 * t92;
404   t215 = params->cc / t7 / t89;
405   t219 = t33 * (t121 * t41 / 0.3e1 + 0.4e1 / 0.3e1 * t125 * t47 - 0.128e3 / 0.3e1 * t57 * t210 + 0.128e3 * t215 * t172);
406   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
407     v3rhosigma2[0] = -t100 * t134 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t219;
408 
409   t222 = t16 * t166;
410   t224 = params->cc / t222;
411   t227 = 0.1e1 / t166;
412   t228 = params->cc * t227;
413   t232 = 0.1e1 / t7 / t110;
414   t233 = params->cc * t232;
415   t234 = t172 * t124;
416   t238 = t33 * (0.3e1 / 0.8e1 * t224 * t25 + 0.3e1 / 0.2e1 * t228 * t126 - 0.48e2 * t233 * t234);
417   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
418     v3sigma3[0] = -0.3e1 / 0.8e1 * t32 * t238;
419 
420 #ifndef XC_DONT_COMPILE_LXC
421 
422   if(order < 4) return;
423 
424 
425   t254 = 0.1e1 / t11 / t89;
426   t267 = t168 * rho[0];
427   t278 = 0.1e1 / t171 / t23;
428   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
429     v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t13 * t27 + 0.4e1 / 0.9e1 * t4 * t139 * t50 - t4 * t68 * t96 - 0.2e1 * t4 * t8 * t176 - 0.3e1 / 0.8e1 * t32 * t33 * (0.35912208504801097395e1 * t9 * t254 + 0.3640e4 / 0.81e2 * t17 * t232 * t24 - 0.23072e5 / 0.27e2 * t44 * t254 * t46 + 0.476672e6 / 0.81e2 * t88 / t168 * t92 - 0.475136e6 / 0.27e2 * t167 / t7 / t267 * t172 + 0.524288e6 / 0.27e2 * params->cc * t222 / t11 / t168 / t10 * t278);
430 
431   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
432     v4rho3sigma[0] = t138 * t1 * t64 / 0.9e1 - t180 * t117 / 0.2e1 - 0.3e1 / 0.2e1 * t100 * t201 - 0.3e1 / 0.8e1 * t32 * t33 * (-0.63374485596707818932e0 * params->bb * t150 - 0.140e3 / 0.27e2 * t57 * t155 + 0.1456e4 / 0.9e1 * params->cc * t150 * t46 - 0.40576e5 / 0.27e2 * params->cc * t162 * t113 + 0.50176e5 / 0.9e1 * params->cc * t170 * t197 - 0.65536e5 / 0.9e1 * params->cc / t11 / t267 * t278 * t87);
433 
434   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
435     v4rho2sigma2[0] = -t180 * t134 / 0.6e1 - t100 * t219 - 0.3e1 / 0.8e1 * t32 * t33 * (-0.7e1 / 0.9e1 * t121 * t81 - 0.28e2 / 0.9e1 * t125 * t84 + 0.2048e4 / 0.9e1 * t57 * t93 - 0.4480e4 / 0.3e1 * t196 * t172 + 0.8192e4 / 0.3e1 * params->cc / t11 / t168 * t278 * t16);
436 
437   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
438     v4rhosigma3[0] = -t100 * t238 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t33 * (-t224 * t41 / 0.2e1 - 0.2e1 * t228 * t47 + 0.16e2 * t121 * t210 + 0.256e3 * t215 * t234 - 0.1024e4 * params->cc / t11 / t161 * t278 * t56);
439 
440   t357 = t166 * sigma[0];
441   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
442     v4sigma4[0] = -0.3e1 / 0.8e1 * t32 * t33 * (-0.15e2 / 0.16e2 * params->cc / t16 / t357 * t25 - 0.15e2 / 0.4e1 * params->cc / t357 * t126 - 0.6e1 * t224 * t128 * t92 + 0.384e3 * params->cc * t254 * t278 * t120 + 0.48e2 * t233 * t172 * t227);
443 
444 #ifndef XC_DONT_COMPILE_MXC
445 
446   if(order < 5) return;
447 
448 
449 #endif
450 
451 #endif
452 
453 #endif
454 
455 #endif
456 
457 #endif
458 
459 
460 }
461 
462 
463 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))464 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
465 {
466 
467 #ifndef XC_DONT_COMPILE_EXC
468   double t1, t3, t4, t5, t6, t7, t8, t9;
469   double t10, t11, t13, t14, t15, t16, t17, t18;
470   double t19, t20, t21, t23, t26, t27, t29, t30;
471   double t33, t34, t35, t37, t39, t41, t42, t43;
472   double t44, t45, t46, t47, t48, t50, t53, t54;
473   double t56, t59, t60, t61, t63, t65;
474 
475 #ifndef XC_DONT_COMPILE_VXC
476   double t67, t68, t69, t70, t71, t72, t74, t75;
477   double t77, t78, t79, t80, t81, t83, t84, t85;
478   double t87, t92, t95, t96, t97, t98, t101, t103;
479   double t104, t105, t106, t107, t109, t110, t111, t113;
480   double t114, t118, t119, t121, t122, t123, t124, t126;
481   double t127, t128, t130, t135, t138, t139, t140, t141;
482   double t144, t146, t147, t151, t152, t153, t156, t157;
483   double t163, t167, t170, t171, t177;
484 
485 #ifndef XC_DONT_COMPILE_FXC
486   double t181, t183, t184, t185, t186, t187, t188, t190;
487   double t191, t192, t194, t195, t196, t198, t199, t200;
488   double t201, t202, t203, t205, t206, t208, t209, t211;
489   double t212, t214, t215, t216, t218, t223, t226, t229;
490   double t230, t231, t232, t234, t235, t238, t240, t241;
491   double t242, t243, t244, t245, t246, t248, t249, t250;
492   double t252, t253, t254, t255, t257, t258, t259, t261;
493   double t262, t263, t265, t266, t267, t268, t270, t273;
494   double t275, t277, t280, t281, t282, t284, t288, t289;
495   double t290, t291, t293, t296, t298, t299, t301, t304;
496   double t306, t309, t311, t314, t315, t317, t318, t319;
497   double t320, t321, t323, t324, t325, t326, t328, t329;
498   double t330, t331, t333, t334, t335, t337, t338, t339;
499   double t340, t342, t347, t350, t353, t354, t355, t356;
500   double t358, t359, t362, t364, t365, t366, t369, t371;
501   double t372, t383, t384, t386, t389, t394, t395, t396;
502   double t401, t406, t416, t417, t419, t422, t423, t426;
503   double t427, t430, t431, t432, t434, t439, t440, t443;
504   double t444, t447, t448, t449, t451, t456, t457;
505 
506 #ifndef XC_DONT_COMPILE_KXC
507   double t464, t466, t470, t471, t472, t473, t475, t476;
508   double t477, t479, t481, t483, t484, t485, t487, t488;
509   double t490, t491, t492, t493, t495, t496, t497, t499;
510   double t501, t503, t504, t505, t507, t508, t509, t511;
511   double t512, t514, t517, t519, t521, t522, t523, t525;
512   double t526, t527, t528, t529, t531, t532, t534, t535;
513   double t536, t538, t539, t541, t543, t544, t546, t547;
514   double t548, t550, t551, t553, t554, t556, t559, t564;
515   double t570, t571, t575, t576, t577, t579, t580, t581;
516   double t585, t587, t588, t589, t594, t595, t598, t599;
517   double t600, t602, t604, t606, t607, t608, t609, t611;
518   double t613, t614, t615, t617, t619, t621, t622, t623;
519   double t625, t627, t629, t631, t633, t635, t636, t638;
520   double t643, t646, t649, t651, t653, t656, t657, t659;
521   double t661, t663, t665, t667, t669, t672, t673, t675;
522   double t676, t677, t679, t681, t683, t685, t687, t696;
523   double t699, t704, t706, t708, t710, t712, t714, t716;
524   double t718, t720, t722, t724, t726, t728, t730, t735;
525   double t738, t741, t743, t745, t748, t750, t752, t754;
526   double t756, t759, t760, t762, t764, t766, t767, t768;
527   double t770, t772, t774, t777, t785, t788, t796, t797;
528   double t802, t807, t813, t814, t818, t819, t820, t822;
529   double t823, t824, t828, t830, t831, t832, t836, t837;
530   double t839, t840, t841, t843, t844, t845, t847, t848;
531   double t849, t851, t852, t853, t854, t856, t857, t859;
532   double t860, t862, t863, t864, t866, t867, t868, t869;
533   double t871, t872, t873, t875, t876, t877, t881, t885;
534   double t887, t888, t892, t895, t911, t912, t915, t920;
535   double t924, t925, t929, t935, t936, t937, t939, t940;
536   double t941, t942, t946, t949, t950, t952, t956, t958;
537   double t959, t963, t966, t969, t974, t978, t983, t987;
538   double t990, t1005, t1006, t1009, t1010, t1013, t1016, t1023;
539   double t1028, t1031, t1032, t1035, t1037, t1041, t1046, t1053;
540   double t1058, t1061, t1062, t1065, t1067, t1070, t1071, t1075;
541   double t1076, t1077, t1080, t1081, t1084, t1086, t1089, t1090;
542   double t1094, t1095, t1096, t1099, t1100;
543 
544 #ifndef XC_DONT_COMPILE_LXC
545   double t1114, t1118, t1121, t1129, t1143, t1152, t1156, t1158;
546   double t1161, t1164, t1167, t1180, t1191, t1199, t1205, t1207;
547   double t1209, t1210, t1211, t1216, t1218, t1220, t1222, t1224;
548   double t1229, t1231, t1238, t1241, t1244, t1247, t1253, t1256;
549   double t1259, t1273, t1275, t1277, t1283, t1290, t1292, t1293;
550   double t1294, t1299, t1301, t1303, t1305, t1307, t1310, t1319;
551   double t1328, t1331, t1334, t1336, t1341, t1348, t1351, t1352;
552   double t1366, t1367, t1368, t1370, t1372, t1383, t1384, t1393;
553   double t1395, t1399, t1406, t1412, t1414, t1417, t1423, t1430;
554   double t1435, t1438, t1446, t1452, t1455, t1458, t1459, t1461;
555   double t1469, t1471, t1473, t1476, t1496, t1501, t1503, t1506;
556   double t1512, t1520, t1523, t1537, t1540, t1547, t1553, t1556;
557   double t1559, t1561, t1567, t1570, t1573, t1596, t1602, t1605;
558   double t1608, t1615, t1618, t1625, t1632, t1639, t1655, t1666;
559   double t1671, t1694, t1704, t1706, t1725, t1741, t1751, t1762;
560   double t1775, t1790, t1793, t1796, t1799, t1802, t1805, t1808;
561   double t1816, t1838, t1858, t1899, t1906, t1917, t1925, t1932;
562   double t1939, t1943, t1951, t1955, t1986, t1991, t1996, t2001;
563   double t2007, t2035, t2050, t2063, t2074, t2085, t2089, t2101;
564   double t2103, t2106, t2124, t2127, t2131, t2141, t2144, t2170;
565   double t2173, t2176, t2180, t2193, t2194, t2207, t2208, t2210;
566   double t2212, t2223, t2226, t2230, t2239, t2242, t2243, t2245;
567   double t2258, t2261, t2274, t2286, t2303, t2311, t2320, t2329;
568   double t2402, t2403, t2405, t2436, t2437, t2449, t2464, t2474;
569   double t2519, t2543, t2575, t2599;
570 #endif
571 
572 #endif
573 
574 #endif
575 
576 #endif
577 
578 #endif
579 
580 
581   gga_x_ol2_params *params;
582 
583   assert(p->params != NULL);
584   params = (gga_x_ol2_params * )(p->params);
585 
586   t1 = M_CBRT3;
587   t3 = POW_1_3(0.1e1 / M_PI);
588   t4 = t1 * t3;
589   t5 = M_CBRT4;
590   t6 = t5 * t5;
591   t7 = t4 * t6;
592   t8 = rho[0] - rho[1];
593   t9 = rho[0] + rho[1];
594   t10 = 0.1e1 / t9;
595   t11 = t8 * t10;
596   t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
597   t14 = POW_1_3(t13);
598   t15 = t14 * t13;
599   t16 = POW_1_3(t9);
600   t17 = t15 * t16;
601   t18 = params->bb * sigma[0];
602   t19 = rho[0] * rho[0];
603   t20 = POW_1_3(rho[0]);
604   t21 = t20 * t20;
605   t23 = 0.1e1 / t21 / t19;
606   t26 = sqrt(sigma[0]);
607   t27 = params->cc * t26;
608   t29 = 0.1e1 / t20 / rho[0];
609   t30 = M_CBRT2;
610   t33 = 0.4e1 * t26 * t29 + t30;
611   t34 = 0.1e1 / t33;
612   t35 = t29 * t34;
613   t37 = params->aa + 0.13888888888888888889e-1 * t18 * t23 + t27 * t35;
614   t39 = t7 * t17 * t37;
615   t41 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
616   t42 = POW_1_3(t41);
617   t43 = t42 * t41;
618   t44 = t43 * t16;
619   t45 = params->bb * sigma[2];
620   t46 = rho[1] * rho[1];
621   t47 = POW_1_3(rho[1]);
622   t48 = t47 * t47;
623   t50 = 0.1e1 / t48 / t46;
624   t53 = sqrt(sigma[2]);
625   t54 = params->cc * t53;
626   t56 = 0.1e1 / t47 / rho[1];
627   t59 = 0.4e1 * t53 * t56 + t30;
628   t60 = 0.1e1 / t59;
629   t61 = t56 * t60;
630   t63 = params->aa + 0.13888888888888888889e-1 * t45 * t50 + t54 * t61;
631   t65 = t7 * t44 * t63;
632   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
633     zk[0] = -0.3e1 / 0.8e1 * t39 - 0.3e1 / 0.8e1 * t65;
634 
635 #ifndef XC_DONT_COMPILE_VXC
636 
637   if(order < 1) return;
638 
639 
640   t67 = 0.3e1 / 0.8e1 * t39;
641   t68 = 0.3e1 / 0.8e1 * t65;
642   t69 = t14 * t16;
643   t70 = t9 * t9;
644   t71 = 0.1e1 / t70;
645   t72 = t8 * t71;
646   t74 = t10 / 0.2e1 - t72 / 0.2e1;
647   t75 = t37 * t74;
648   t77 = t7 * t69 * t75;
649   t78 = t77 / 0.2e1;
650   t79 = t16 * t16;
651   t80 = 0.1e1 / t79;
652   t81 = t15 * t80;
653   t83 = t7 * t81 * t37;
654   t84 = t83 / 0.8e1;
655   t85 = t19 * rho[0];
656   t87 = 0.1e1 / t21 / t85;
657   t92 = 0.1e1 / t20 / t19 * t34;
658   t95 = params->cc * sigma[0];
659   t96 = t33 * t33;
660   t97 = 0.1e1 / t96;
661   t98 = t87 * t97;
662   t101 = -0.37037037037037037037e-1 * t18 * t87 - 0.4e1 / 0.3e1 * t27 * t92 + 0.16e2 / 0.3e1 * t95 * t98;
663   t103 = t7 * t17 * t101;
664   t104 = 0.3e1 / 0.8e1 * t103;
665   t105 = t42 * t16;
666   t106 = -t74;
667   t107 = t63 * t106;
668   t109 = t7 * t105 * t107;
669   t110 = t109 / 0.2e1;
670   t111 = t43 * t80;
671   t113 = t7 * t111 * t63;
672   t114 = t113 / 0.8e1;
673   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
674     vrho[0] = -t67 - t68 + t9 * (-t78 - t84 - t104 - t110 - t114);
675 
676   t118 = -t10 / 0.2e1 - t72 / 0.2e1;
677   t119 = t37 * t118;
678   t121 = t7 * t69 * t119;
679   t122 = t121 / 0.2e1;
680   t123 = -t118;
681   t124 = t63 * t123;
682   t126 = t7 * t105 * t124;
683   t127 = t126 / 0.2e1;
684   t128 = t46 * rho[1];
685   t130 = 0.1e1 / t48 / t128;
686   t135 = 0.1e1 / t47 / t46 * t60;
687   t138 = params->cc * sigma[2];
688   t139 = t59 * t59;
689   t140 = 0.1e1 / t139;
690   t141 = t130 * t140;
691   t144 = -0.37037037037037037037e-1 * t45 * t130 - 0.4e1 / 0.3e1 * t54 * t135 + 0.16e2 / 0.3e1 * t138 * t141;
692   t146 = t7 * t44 * t144;
693   t147 = 0.3e1 / 0.8e1 * t146;
694   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
695     vrho[1] = -t67 - t68 + t9 * (-t122 - t84 - t127 - t114 - t147);
696 
697   t151 = t16 * t9 * t1;
698   t152 = t151 * t3;
699   t153 = t6 * t15;
700   t156 = 0.1e1 / t26;
701   t157 = params->cc * t156;
702   t163 = 0.13888888888888888889e-1 * params->bb * t23 + t157 * t35 / 0.2e1 - 0.2e1 * params->cc * t23 * t97;
703   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
704     vsigma[0] = -0.3e1 / 0.8e1 * t152 * t153 * t163;
705 
706   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
707     vsigma[1] = 0.0e0;
708 
709   t167 = t6 * t43;
710   t170 = 0.1e1 / t53;
711   t171 = params->cc * t170;
712   t177 = 0.13888888888888888889e-1 * params->bb * t50 + t171 * t61 / 0.2e1 - 0.2e1 * params->cc * t50 * t140;
713   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
714     vsigma[2] = -0.3e1 / 0.8e1 * t152 * t167 * t177;
715 
716 #ifndef XC_DONT_COMPILE_FXC
717 
718   if(order < 2) return;
719 
720 
721   t181 = t83 / 0.4e1;
722   t183 = t113 / 0.4e1;
723   t184 = t14 * t14;
724   t185 = 0.1e1 / t184;
725   t186 = t185 * t16;
726   t187 = t74 * t74;
727   t188 = t37 * t187;
728   t190 = t7 * t186 * t188;
729   t191 = t190 / 0.6e1;
730   t192 = t14 * t80;
731   t194 = t7 * t192 * t75;
732   t195 = t194 / 0.3e1;
733   t196 = t101 * t74;
734   t198 = t7 * t69 * t196;
735   t199 = t70 * t9;
736   t200 = 0.1e1 / t199;
737   t201 = t8 * t200;
738   t202 = -t71 + t201;
739   t203 = t37 * t202;
740   t205 = t7 * t69 * t203;
741   t206 = t205 / 0.2e1;
742   t208 = 0.1e1 / t79 / t9;
743   t209 = t15 * t208;
744   t211 = t7 * t209 * t37;
745   t212 = t211 / 0.12e2;
746   t214 = t7 * t81 * t101;
747   t215 = t214 / 0.4e1;
748   t216 = t19 * t19;
749   t218 = 0.1e1 / t21 / t216;
750   t223 = 0.1e1 / t20 / t85 * t34;
751   t226 = t218 * t97;
752   t229 = t26 * sigma[0];
753   t230 = params->cc * t229;
754   t231 = t216 * t19;
755   t232 = 0.1e1 / t231;
756   t234 = 0.1e1 / t96 / t33;
757   t235 = t232 * t234;
758   t238 = 0.13580246913580246914e0 * t18 * t218 + 0.28e2 / 0.9e1 * t27 * t223 - 0.80e2 / 0.3e1 * t95 * t226 + 0.512e3 / 0.9e1 * t230 * t235;
759   t240 = t7 * t17 * t238;
760   t241 = 0.3e1 / 0.8e1 * t240;
761   t242 = t42 * t42;
762   t243 = 0.1e1 / t242;
763   t244 = t243 * t16;
764   t245 = t106 * t106;
765   t246 = t63 * t245;
766   t248 = t7 * t244 * t246;
767   t249 = t248 / 0.6e1;
768   t250 = t42 * t80;
769   t252 = t7 * t250 * t107;
770   t253 = t252 / 0.3e1;
771   t254 = -t202;
772   t255 = t63 * t254;
773   t257 = t7 * t105 * t255;
774   t258 = t257 / 0.2e1;
775   t259 = t43 * t208;
776   t261 = t7 * t259 * t63;
777   t262 = t261 / 0.12e2;
778   t263 = -t191 - t195 - t198 - t206 + t212 - t215 - t241 - t249 - t253 - t258 + t262;
779   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
780     v2rho2[0] = -t77 - t181 - 0.3e1 / 0.4e1 * t103 - t109 - t183 + t9 * t263;
781 
782   t265 = t6 * t185;
783   t266 = t4 * t265;
784   t267 = t16 * t37;
785   t268 = t118 * t74;
786   t270 = t266 * t267 * t268;
787   t273 = t7 * t192 * t119;
788   t275 = t101 * t118;
789   t277 = t7 * t69 * t275;
790   t280 = 0.1e1 / t79 / t70;
791   t281 = t14 * t280;
792   t282 = t37 * t8;
793   t284 = t7 * t281 * t282;
794   t288 = t6 * t243;
795   t289 = t4 * t288;
796   t290 = t16 * t63;
797   t291 = t123 * t106;
798   t293 = t289 * t290 * t291;
799   t296 = t7 * t250 * t124;
800   t298 = t42 * t280;
801   t299 = t63 * t8;
802   t301 = t7 * t298 * t299;
803   t304 = t144 * t106;
804   t306 = t7 * t105 * t304;
805   t309 = t7 * t111 * t144;
806   t311 = -t270 / 0.6e1 - t273 / 0.6e1 - t277 / 0.2e1 - t284 / 0.2e1 - t194 / 0.6e1 + t212 - t214 / 0.8e1 - t293 / 0.6e1 - t296 / 0.6e1 + t301 / 0.2e1 - t252 / 0.6e1 + t262 - t306 / 0.2e1 - t309 / 0.8e1;
807   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
808     v2rho2[1] = t9 * t311 - t104 - t110 - t122 - t127 - t147 - t181 - t183 - t78;
809 
810   t314 = t118 * t118;
811   t315 = t37 * t314;
812   t317 = t7 * t186 * t315;
813   t318 = t317 / 0.6e1;
814   t319 = t273 / 0.3e1;
815   t320 = t71 + t201;
816   t321 = t37 * t320;
817   t323 = t7 * t69 * t321;
818   t324 = t323 / 0.2e1;
819   t325 = t123 * t123;
820   t326 = t63 * t325;
821   t328 = t7 * t244 * t326;
822   t329 = t328 / 0.6e1;
823   t330 = t296 / 0.3e1;
824   t331 = t144 * t123;
825   t333 = t7 * t105 * t331;
826   t334 = -t320;
827   t335 = t63 * t334;
828   t337 = t7 * t105 * t335;
829   t338 = t337 / 0.2e1;
830   t339 = t309 / 0.4e1;
831   t340 = t46 * t46;
832   t342 = 0.1e1 / t48 / t340;
833   t347 = 0.1e1 / t47 / t128 * t60;
834   t350 = t342 * t140;
835   t353 = t53 * sigma[2];
836   t354 = params->cc * t353;
837   t355 = t340 * t46;
838   t356 = 0.1e1 / t355;
839   t358 = 0.1e1 / t139 / t59;
840   t359 = t356 * t358;
841   t362 = 0.13580246913580246914e0 * t45 * t342 + 0.28e2 / 0.9e1 * t54 * t347 - 0.80e2 / 0.3e1 * t138 * t350 + 0.512e3 / 0.9e1 * t354 * t359;
842   t364 = t7 * t44 * t362;
843   t365 = 0.3e1 / 0.8e1 * t364;
844   t366 = -t318 - t319 - t324 + t212 - t329 - t330 - t333 - t338 + t262 - t339 - t365;
845   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
846     v2rho2[2] = -t121 - t181 - t126 - t183 - 0.3e1 / 0.4e1 * t146 + t9 * t366;
847 
848   t369 = t7 * t17 * t163;
849   t371 = t6 * t14;
850   t372 = t163 * t74;
851   t383 = t216 * rho[0];
852   t384 = 0.1e1 / t383;
853   t386 = t234 * t26;
854   t389 = -0.37037037037037037037e-1 * params->bb * t87 - 0.2e1 / 0.3e1 * t157 * t92 + 0.8e1 * params->cc * t87 * t97 - 0.64e2 / 0.3e1 * params->cc * t384 * t386;
855   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
856     v2rhosigma[0] = -t369 / 0.2e1 - t152 * t371 * t372 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * t389;
857 
858   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
859     v2rhosigma[1] = 0.0e0;
860 
861   t394 = t7 * t44 * t177;
862   t395 = t6 * t42;
863   t396 = t177 * t106;
864   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
865     v2rhosigma[2] = -t152 * t395 * t396 / 0.2e1 - t394 / 0.2e1;
866 
867   t401 = t371 * t163 * t118;
868   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
869     v2rhosigma[3] = -t152 * t401 / 0.2e1 - t369 / 0.2e1;
870 
871   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
872     v2rhosigma[4] = 0.0e0;
873 
874   t406 = t395 * t177 * t123;
875   t416 = t340 * rho[1];
876   t417 = 0.1e1 / t416;
877   t419 = t358 * t53;
878   t422 = -0.37037037037037037037e-1 * params->bb * t130 - 0.2e1 / 0.3e1 * t171 * t135 + 0.8e1 * params->cc * t130 * t140 - 0.64e2 / 0.3e1 * params->cc * t417 * t419;
879   t423 = t167 * t422;
880   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
881     v2rhosigma[5] = -t394 / 0.2e1 - t152 * t406 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t423;
882 
883   t426 = 0.1e1 / t229;
884   t427 = params->cc * t426;
885   t430 = 0.1e1 / sigma[0];
886   t431 = params->cc * t430;
887   t432 = t23 * t97;
888   t434 = 0.1e1 / t216;
889   t439 = -t427 * t35 / 0.4e1 - t431 * t432 + 0.8e1 * params->cc * t434 * t234 * t156;
890   t440 = t153 * t439;
891   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
892     v2sigma2[0] = -0.3e1 / 0.8e1 * t152 * t440;
893 
894   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
895     v2sigma2[1] = 0.0e0;
896 
897   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
898     v2sigma2[2] = 0.0e0;
899 
900   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
901     v2sigma2[3] = 0.0e0;
902 
903   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
904     v2sigma2[4] = 0.0e0;
905 
906   t443 = 0.1e1 / t353;
907   t444 = params->cc * t443;
908   t447 = 0.1e1 / sigma[2];
909   t448 = params->cc * t447;
910   t449 = t50 * t140;
911   t451 = 0.1e1 / t340;
912   t456 = -t444 * t61 / 0.4e1 - t448 * t449 + 0.8e1 * params->cc * t451 * t358 * t170;
913   t457 = t167 * t456;
914   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
915     v2sigma2[5] = -0.3e1 / 0.8e1 * t152 * t457;
916 
917 #ifndef XC_DONT_COMPILE_KXC
918 
919   if(order < 3) return;
920 
921 
922   t464 = t211 / 0.4e1;
923   t466 = t261 / 0.4e1;
924   t470 = 0.1e1 / t184 / t13;
925   t471 = t470 * t16;
926   t472 = t187 * t74;
927   t473 = t37 * t472;
928   t475 = t7 * t471 * t473;
929   t476 = t475 / 0.9e1;
930   t477 = t185 * t80;
931   t479 = t7 * t477 * t188;
932   t481 = t74 * t202;
933   t483 = t266 * t267 * t481;
934   t484 = t483 / 0.2e1;
935   t485 = t14 * t208;
936   t487 = t7 * t485 * t75;
937   t488 = t487 / 0.3e1;
938   t490 = 0.1e1 / t242 / t41;
939   t491 = t490 * t16;
940   t492 = t245 * t106;
941   t493 = t63 * t492;
942   t495 = t7 * t491 * t493;
943   t496 = t495 / 0.9e1;
944   t497 = t243 * t80;
945   t499 = t7 * t497 * t246;
946   t501 = t106 * t254;
947   t503 = t289 * t290 * t501;
948   t504 = t503 / 0.2e1;
949   t505 = t42 * t208;
950   t507 = t7 * t505 * t107;
951   t508 = t507 / 0.3e1;
952   t509 = t101 * t187;
953   t511 = t7 * t186 * t509;
954   t512 = t511 / 0.2e1;
955   t514 = t7 * t192 * t196;
956   t517 = t7 * t192 * t203;
957   t519 = t238 * t74;
958   t521 = t7 * t69 * t519;
959   t522 = 0.3e1 / 0.2e1 * t521;
960   t523 = t101 * t202;
961   t525 = t7 * t69 * t523;
962   t526 = 0.3e1 / 0.2e1 * t525;
963   t527 = t70 * t70;
964   t528 = 0.1e1 / t527;
965   t529 = t8 * t528;
966   t531 = 0.3e1 * t200 - 0.3e1 * t529;
967   t532 = t37 * t531;
968   t534 = t7 * t69 * t532;
969   t535 = t534 / 0.2e1;
970   t536 = t15 * t280;
971   t538 = t7 * t536 * t37;
972   t539 = 0.5e1 / 0.36e2 * t538;
973   t541 = t7 * t250 * t255;
974   t543 = -t531;
975   t544 = t63 * t543;
976   t546 = t7 * t105 * t544;
977   t547 = t546 / 0.2e1;
978   t548 = t43 * t280;
979   t550 = t7 * t548 * t63;
980   t551 = 0.5e1 / 0.36e2 * t550;
981   t553 = t7 * t209 * t101;
982   t554 = t553 / 0.4e1;
983   t556 = t7 * t81 * t238;
984   t559 = 0.1e1 / t21 / t383;
985   t564 = 0.1e1 / t20 / t216 * t34;
986   t570 = t216 * t85;
987   t571 = 0.1e1 / t570;
988   t575 = sigma[0] * sigma[0];
989   t576 = params->cc * t575;
990   t577 = t216 * t216;
991   t579 = 0.1e1 / t20 / t577;
992   t580 = t96 * t96;
993   t581 = 0.1e1 / t580;
994   t585 = -0.63374485596707818932e0 * t18 * t559 - 0.280e3 / 0.27e2 * t27 * t564 + 0.3808e4 / 0.27e2 * t95 * t559 * t97 - 0.5632e4 / 0.9e1 * t230 * t571 * t234 + 0.8192e4 / 0.9e1 * t576 * t579 * t581;
995   t587 = t7 * t17 * t585;
996   t588 = 0.3e1 / 0.8e1 * t587;
997   t589 = -t517 / 0.2e1 - t522 - t526 - t535 - t539 - t541 / 0.2e1 - t547 - t551 + t554 - 0.3e1 / 0.8e1 * t556 - t588;
998   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
999     v3rho3[0] = -t190 / 0.2e1 - t194 - t248 / 0.2e1 - t252 - 0.3e1 * t198 - 0.3e1 / 0.2e1 * t205 + t464 - 0.3e1 / 0.2e1 * t257 + t466 - 0.3e1 / 0.4e1 * t214 - 0.9e1 / 0.8e1 * t240 + t9 * (t476 - t479 / 0.6e1 - t484 + t488 + t496 - t499 / 0.6e1 - t504 + t508 - t512 - t514 + t589);
1000 
1001   t594 = t270 / 0.3e1;
1002   t595 = t293 / 0.3e1;
1003   t598 = t6 * t470;
1004   t599 = t4 * t598;
1005   t600 = t118 * t187;
1006   t602 = t599 * t267 * t600;
1007   t604 = t80 * t37;
1008   t606 = t266 * t604 * t268;
1009   t607 = t606 / 0.9e1;
1010   t608 = t280 * t37;
1011   t609 = t8 * t74;
1012   t611 = t266 * t608 * t609;
1013   t613 = t6 * t490;
1014   t614 = t4 * t613;
1015   t615 = t123 * t245;
1016   t617 = t614 * t290 * t615;
1017   t619 = t80 * t63;
1018   t621 = t289 * t619 * t291;
1019   t622 = t621 / 0.9e1;
1020   t623 = t16 * t101;
1021   t625 = t266 * t623 * t268;
1022   t627 = t118 * t202;
1023   t629 = t266 * t267 * t627;
1024   t631 = t123 * t254;
1025   t633 = t289 * t290 * t631;
1026   t635 = t280 * t63;
1027   t636 = t8 * t106;
1028   t638 = t289 * t635 * t636;
1029   t643 = t7 * t281 * t37;
1030   t646 = t7 * t298 * t63;
1031   t649 = t7 * t259 * t144;
1032   t651 = t602 / 0.9e1 - t607 - t611 / 0.3e1 + t617 / 0.9e1 - t622 - t625 / 0.3e1 - t629 / 0.6e1 - t633 / 0.6e1 + t638 / 0.3e1 + t553 / 0.6e1 - t556 / 0.8e1 - t643 / 0.2e1 + t646 / 0.2e1 + t649 / 0.12e2 - t539 - t551;
1033   t653 = t7 * t485 * t119;
1034   t656 = 0.1e1 / t79 / t199;
1035   t657 = t14 * t656;
1036   t659 = t7 * t657 * t282;
1037   t661 = t101 * t8;
1038   t663 = t7 * t281 * t661;
1039   t665 = t7 * t505 * t124;
1040   t667 = t144 * t245;
1041   t669 = t7 * t244 * t667;
1042   t672 = t7 * t250 * t304;
1043   t673 = t672 / 0.3e1;
1044   t675 = t7 * t192 * t275;
1045   t676 = t675 / 0.3e1;
1046   t677 = t238 * t118;
1047   t679 = t7 * t69 * t677;
1048   t681 = t42 * t656;
1049   t683 = t7 * t681 * t299;
1050   t685 = t144 * t254;
1051   t687 = t7 * t105 * t685;
1052   t696 = t653 / 0.9e1 + 0.7e1 / 0.6e1 * t659 - t663 + t665 / 0.9e1 - t669 / 0.6e1 - t673 - t676 - t679 / 0.2e1 - 0.7e1 / 0.6e1 * t683 - t687 / 0.2e1 - t479 / 0.18e2 + 0.2e1 / 0.9e1 * t487 - t499 / 0.18e2 + 0.2e1 / 0.9e1 * t507 - t514 / 0.3e1 - t517 / 0.6e1 - t541 / 0.6e1;
1053   t699 = -t206 + t464 - t258 + t466 - t277 + t301 - t306 - t339 - t214 / 0.2e1 - t241 + t9 * (t651 + t696);
1054   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1055     v3rho3[1] = -t191 - 0.2e1 / 0.3e1 * t194 - t249 - 0.2e1 / 0.3e1 * t252 - t594 - t319 - t284 - t595 - t330 - t198 + t699;
1056 
1057   t704 = t320 * t74;
1058   t706 = t266 * t267 * t704;
1059   t708 = t325 * t106;
1060   t710 = t614 * t290 * t708;
1061   t712 = t118 * t8;
1062   t714 = t266 * t608 * t712;
1063   t716 = t123 * t8;
1064   t718 = t289 * t635 * t716;
1065   t720 = t16 * t144;
1066   t722 = t289 * t720 * t291;
1067   t724 = t334 * t106;
1068   t726 = t289 * t290 * t724;
1069   t728 = t314 * t74;
1070   t730 = t599 * t267 * t728;
1071   t735 = t7 * t111 * t362;
1072   t738 = t7 * t250 * t335;
1073   t741 = t7 * t477 * t315;
1074   t743 = -t706 / 0.6e1 + t710 / 0.9e1 - t714 / 0.3e1 + t718 / 0.3e1 - t722 / 0.3e1 - t726 / 0.6e1 - t607 - t622 + t730 / 0.9e1 + t553 / 0.12e2 + t649 / 0.6e1 - t735 / 0.8e1 - t539 - t551 - t738 / 0.6e1 - t741 / 0.18e2;
1075   t745 = t7 * t192 * t321;
1076   t748 = t7 * t497 * t326;
1077   t750 = t101 * t314;
1078   t752 = t7 * t186 * t750;
1079   t754 = t101 * t320;
1080   t756 = t7 * t69 * t754;
1081   t759 = -t200 - 0.3e1 * t529;
1082   t760 = t37 * t759;
1083   t762 = t7 * t69 * t760;
1084   t764 = t144 * t8;
1085   t766 = t7 * t298 * t764;
1086   t767 = -t759;
1087   t768 = t63 * t767;
1088   t770 = t7 * t105 * t768;
1089   t772 = t362 * t106;
1090   t774 = t7 * t105 * t772;
1091   t777 = t7 * t250 * t331;
1092   t785 = -t745 / 0.6e1 - t748 / 0.18e2 - t752 / 0.6e1 - t756 / 0.2e1 - t762 / 0.2e1 + t766 - t770 / 0.2e1 - t774 / 0.2e1 - t777 / 0.3e1 + 0.2e1 / 0.9e1 * t653 - t659 / 0.3e1 + 0.2e1 / 0.9e1 * t665 - t673 - t676 + t683 / 0.3e1 + t487 / 0.9e1 + t507 / 0.9e1;
1093   t788 = t466 - t277 + t301 - t306 - t309 / 0.2e1 - t324 - t333 - t338 - t215 - t365 + t9 * (t743 + t785);
1094   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1095     v3rho3[2] = -t195 - t253 - t594 - 0.2e1 / 0.3e1 * t273 - t284 - t595 - 0.2e1 / 0.3e1 * t296 - t318 - t329 + t464 + t788;
1096 
1097   t796 = t653 / 0.3e1;
1098   t797 = t665 / 0.3e1;
1099   t802 = 0.1e1 / t48 / t416;
1100   t807 = 0.1e1 / t47 / t340 * t60;
1101   t813 = t340 * t128;
1102   t814 = 0.1e1 / t813;
1103   t818 = sigma[2] * sigma[2];
1104   t819 = params->cc * t818;
1105   t820 = t340 * t340;
1106   t822 = 0.1e1 / t47 / t820;
1107   t823 = t139 * t139;
1108   t824 = 0.1e1 / t823;
1109   t828 = -0.63374485596707818932e0 * t45 * t802 - 0.280e3 / 0.27e2 * t54 * t807 + 0.3808e4 / 0.27e2 * t138 * t802 * t140 - 0.5632e4 / 0.9e1 * t354 * t814 * t358 + 0.8192e4 / 0.9e1 * t819 * t822 * t824;
1110   t830 = t7 * t44 * t828;
1111   t831 = 0.3e1 / 0.8e1 * t830;
1112   t832 = t649 / 0.4e1;
1113   t836 = -0.3e1 * t200 - 0.3e1 * t529;
1114   t837 = t37 * t836;
1115   t839 = t7 * t69 * t837;
1116   t840 = t839 / 0.2e1;
1117   t841 = t144 * t325;
1118   t843 = t7 * t244 * t841;
1119   t844 = t843 / 0.2e1;
1120   t845 = t362 * t123;
1121   t847 = t7 * t105 * t845;
1122   t848 = 0.3e1 / 0.2e1 * t847;
1123   t849 = t144 * t334;
1124   t851 = t7 * t105 * t849;
1125   t852 = 0.3e1 / 0.2e1 * t851;
1126   t853 = -t836;
1127   t854 = t63 * t853;
1128   t856 = t7 * t105 * t854;
1129   t857 = t856 / 0.2e1;
1130   t859 = t314 * t118;
1131   t860 = t37 * t859;
1132   t862 = t7 * t471 * t860;
1133   t863 = t862 / 0.9e1;
1134   t864 = t118 * t320;
1135   t866 = t266 * t267 * t864;
1136   t867 = t866 / 0.2e1;
1137   t868 = t325 * t123;
1138   t869 = t63 * t868;
1139   t871 = t7 * t491 * t869;
1140   t872 = t871 / 0.9e1;
1141   t873 = t123 * t334;
1142   t875 = t289 * t290 * t873;
1143   t876 = t875 / 0.2e1;
1144   t877 = -t840 - t844 - t848 - t852 - t857 - t777 - t738 / 0.2e1 + t863 - t867 + t872 - t876;
1145   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1146     v3rho3[3] = -t273 - t296 - t317 / 0.2e1 - t328 / 0.2e1 + t464 + t466 - 0.3e1 / 0.4e1 * t309 - 0.3e1 / 0.2e1 * t323 - 0.3e1 * t333 - 0.3e1 / 0.2e1 * t337 - 0.9e1 / 0.8e1 * t364 + t9 * (t796 + t797 - t741 / 0.6e1 - t745 / 0.2e1 - t748 / 0.6e1 - t539 - t551 - t831 + t832 - 0.3e1 / 0.8e1 * t735 + t877);
1147 
1148   t881 = t7 * t69 * t372;
1149   t885 = t7 * t81 * t163 / 0.6e1;
1150   t887 = t7 * t17 * t389;
1151   t888 = t163 * t187;
1152   t892 = t389 * t74;
1153   t895 = t163 * t202;
1154   t911 = params->cc / t20 / t570;
1155   t912 = t581 * sigma[0];
1156   t915 = 0.13580246913580246914e0 * params->bb * t218 + 0.14e2 / 0.9e1 * t157 * t223 - 0.296e3 / 0.9e1 * params->cc * t218 * t97 + 0.192e3 * params->cc * t232 * t386 - 0.1024e4 / 0.3e1 * t911 * t912;
1157   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1158     v3rho2sigma[0] = -0.4e1 / 0.3e1 * t881 - t885 - t887 - t152 * t265 * t888 / 0.6e1 - t152 * t371 * t892 - t152 * t371 * t895 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * t915;
1159 
1160   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1161     v3rho2sigma[1] = 0.0e0;
1162 
1163   t920 = t7 * t105 * t396;
1164   t924 = t7 * t111 * t177 / 0.6e1;
1165   t925 = t177 * t245;
1166   t929 = t177 * t254;
1167   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1168     v3rho2sigma[2] = -0.4e1 / 0.3e1 * t920 - t924 - t152 * t288 * t925 / 0.6e1 - t152 * t395 * t929 / 0.2e1;
1169 
1170   t935 = t16 * t1;
1171   t936 = t935 * t3;
1172   t937 = t936 * t401;
1173   t939 = t3 * t6;
1174   t940 = t151 * t939;
1175   t941 = t185 * t163;
1176   t942 = t941 * t268;
1177   t946 = t371 * t389 * t118;
1178   t949 = t208 * t1;
1179   t950 = t949 * t3;
1180   t952 = t371 * t163 * t8;
1181   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1182     v3rho2sigma[3] = -0.2e1 / 0.3e1 * t881 - t885 - t887 / 0.2e1 - 0.2e1 / 0.3e1 * t937 - t940 * t942 / 0.6e1 - t152 * t946 / 0.2e1 - t950 * t952 / 0.2e1;
1183 
1184   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1185     v3rho2sigma[4] = 0.0e0;
1186 
1187   t956 = t936 * t406;
1188   t958 = t243 * t177;
1189   t959 = t958 * t291;
1190   t963 = t395 * t177 * t8;
1191   t966 = t936 * t423;
1192   t969 = t395 * t422 * t106;
1193   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1194     v3rho2sigma[5] = -0.2e1 / 0.3e1 * t920 - t924 - 0.2e1 / 0.3e1 * t956 - t940 * t959 / 0.6e1 + t950 * t963 / 0.2e1 - t966 / 0.2e1 - t152 * t969 / 0.2e1;
1195 
1196   t974 = t265 * t163 * t314;
1197   t978 = t371 * t163 * t320;
1198   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1199     v3rho2sigma[6] = -0.4e1 / 0.3e1 * t937 - t885 - t152 * t974 / 0.6e1 - t152 * t978 / 0.2e1;
1200 
1201   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1202     v3rho2sigma[7] = 0.0e0;
1203 
1204   t983 = t288 * t177 * t325;
1205   t987 = t395 * t422 * t123;
1206   t990 = t395 * t177 * t334;
1207   t1005 = params->cc / t47 / t813;
1208   t1006 = t824 * sigma[2];
1209   t1009 = 0.13580246913580246914e0 * params->bb * t342 + 0.14e2 / 0.9e1 * t171 * t347 - 0.296e3 / 0.9e1 * params->cc * t342 * t140 + 0.192e3 * params->cc * t356 * t419 - 0.1024e4 / 0.3e1 * t1005 * t1006;
1210   t1010 = t167 * t1009;
1211   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1212     v3rho2sigma[8] = -0.4e1 / 0.3e1 * t956 - t924 - t966 - t152 * t983 / 0.6e1 - t152 * t987 - t152 * t990 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1010;
1213 
1214   t1013 = t936 * t440;
1215   t1016 = t371 * t439 * t74;
1216   t1023 = t384 * t234;
1217   t1028 = params->cc / t20 / t231;
1218   t1031 = t427 * t92 / 0.3e1 + 0.4e1 / 0.3e1 * t431 * t98 - 0.128e3 / 0.3e1 * t157 * t1023 + 0.128e3 * t1028 * t581;
1219   t1032 = t153 * t1031;
1220   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1221     v3rhosigma2[0] = -t1013 / 0.2e1 - t152 * t1016 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1032;
1222 
1223   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1224     v3rhosigma2[1] = 0.0e0;
1225 
1226   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1227     v3rhosigma2[2] = 0.0e0;
1228 
1229   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1230     v3rhosigma2[3] = 0.0e0;
1231 
1232   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1233     v3rhosigma2[4] = 0.0e0;
1234 
1235   t1035 = t936 * t457;
1236   t1037 = t395 * t456 * t106;
1237   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1238     v3rhosigma2[5] = -t152 * t1037 / 0.2e1 - t1035 / 0.2e1;
1239 
1240   t1041 = t371 * t439 * t118;
1241   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1242     v3rhosigma2[6] = -t152 * t1041 / 0.2e1 - t1013 / 0.2e1;
1243 
1244   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1245     v3rhosigma2[7] = 0.0e0;
1246 
1247   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1248     v3rhosigma2[8] = 0.0e0;
1249 
1250   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1251     v3rhosigma2[9] = 0.0e0;
1252 
1253   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1254     v3rhosigma2[10] = 0.0e0;
1255 
1256   t1046 = t395 * t456 * t123;
1257   t1053 = t417 * t358;
1258   t1058 = params->cc / t47 / t355;
1259   t1061 = t444 * t135 / 0.3e1 + 0.4e1 / 0.3e1 * t448 * t141 - 0.128e3 / 0.3e1 * t171 * t1053 + 0.128e3 * t1058 * t824;
1260   t1062 = t167 * t1061;
1261   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1262     v3rhosigma2[11] = -t1035 / 0.2e1 - t152 * t1046 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1062;
1263 
1264   t1065 = t26 * t575;
1265   t1067 = params->cc / t1065;
1266   t1070 = 0.1e1 / t575;
1267   t1071 = params->cc * t1070;
1268   t1075 = 0.1e1 / t20 / t383;
1269   t1076 = params->cc * t1075;
1270   t1077 = t581 * t430;
1271   t1080 = 0.3e1 / 0.8e1 * t1067 * t35 + 0.3e1 / 0.2e1 * t1071 * t432 - 0.48e2 * t1076 * t1077;
1272   t1081 = t153 * t1080;
1273   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1274     v3sigma3[0] = -0.3e1 / 0.8e1 * t152 * t1081;
1275 
1276   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1277     v3sigma3[1] = 0.0e0;
1278 
1279   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1280     v3sigma3[2] = 0.0e0;
1281 
1282   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1283     v3sigma3[3] = 0.0e0;
1284 
1285   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1286     v3sigma3[4] = 0.0e0;
1287 
1288   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1289     v3sigma3[5] = 0.0e0;
1290 
1291   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1292     v3sigma3[6] = 0.0e0;
1293 
1294   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1295     v3sigma3[7] = 0.0e0;
1296 
1297   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1298     v3sigma3[8] = 0.0e0;
1299 
1300   t1084 = t53 * t818;
1301   t1086 = params->cc / t1084;
1302   t1089 = 0.1e1 / t818;
1303   t1090 = params->cc * t1089;
1304   t1094 = 0.1e1 / t47 / t416;
1305   t1095 = params->cc * t1094;
1306   t1096 = t824 * t447;
1307   t1099 = 0.3e1 / 0.8e1 * t1086 * t61 + 0.3e1 / 0.2e1 * t1090 * t449 - 0.48e2 * t1095 * t1096;
1308   t1100 = t167 * t1099;
1309   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1310     v3sigma3[9] = -0.3e1 / 0.8e1 * t152 * t1100;
1311 
1312 #ifndef XC_DONT_COMPILE_LXC
1313 
1314   if(order < 4) return;
1315 
1316 
1317   t1114 = 0.4e1 / 0.9e1 * t475 - 0.2e1 / 0.3e1 * t479 - 0.2e1 * t483 + 0.4e1 / 0.3e1 * t487 + 0.4e1 / 0.9e1 * t495 - 0.2e1 / 0.3e1 * t499 - 0.2e1 * t503 + 0.4e1 / 0.3e1 * t507 - 0.2e1 * t511 - 0.4e1 * t514 - 0.2e1 * t517;
1318   t1118 = 0.5e1 / 0.9e1 * t538;
1319   t1121 = 0.5e1 / 0.9e1 * t550;
1320   t1129 = t266 * t604 * t481;
1321   t1143 = t289 * t619 * t501;
1322   t1152 = 0.10e2 / 0.27e2 * t7 * t15 * t656 * t37;
1323   t1156 = 0.10e2 / 0.27e2 * t7 * t43 * t656 * t63;
1324   t1158 = t7 * t536 * t101;
1325   t1161 = t7 * t209 * t238;
1326   t1164 = t7 * t81 * t585;
1327   t1167 = 0.1e1 / t21 / t231;
1328   t1180 = t577 * rho[0];
1329   t1191 = 0.1e1 / t580 / t33;
1330   t1199 = t202 * t202;
1331   t1205 = t7 * t281 * t75;
1332   t1207 = t41 * t41;
1333   t1209 = 0.1e1 / t242 / t1207;
1334   t1210 = t1209 * t16;
1335   t1211 = t245 * t245;
1336   t1216 = t490 * t80;
1337   t1218 = t7 * t1216 * t493;
1338   t1220 = t243 * t208;
1339   t1222 = t7 * t1220 * t246;
1340   t1224 = t254 * t254;
1341   t1229 = 0.2e1 / 0.3e1 * t599 * t267 * t187 * t202 - 0.2e1 / 0.3e1 * t1129 - 0.2e1 * t266 * t623 * t481 - 0.2e1 / 0.3e1 * t266 * t267 * t74 * t531 + 0.2e1 / 0.3e1 * t614 * t290 * t245 * t254 - 0.2e1 / 0.3e1 * t1143 - 0.2e1 / 0.3e1 * t289 * t290 * t106 * t543 + t1152 + t1156 - 0.5e1 / 0.9e1 * t1158 + t1161 / 0.2e1 - t1164 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * (0.35912208504801097395e1 * t18 * t1167 + 0.3640e4 / 0.81e2 * t27 * t1075 * t34 - 0.23072e5 / 0.27e2 * t95 * t1167 * t97 + 0.476672e6 / 0.81e2 * t230 / t577 * t234 - 0.475136e6 / 0.27e2 * t576 / t20 / t1180 * t581 + 0.524288e6 / 0.27e2 * params->cc * t1065 / t21 / t577 / t19 * t1191) - t7 * t186 * t37 * t1199 / 0.2e1 - 0.20e2 / 0.27e2 * t1205 - 0.5e1 / 0.27e2 * t7 * t1210 * t63 * t1211 + 0.4e1 / 0.27e2 * t1218 + 0.2e1 / 0.9e1 * t1222 - t7 * t244 * t63 * t1224 / 0.2e1;
1342   t1231 = t7 * t298 * t107;
1343   t1238 = t7 * t477 * t509;
1344   t1241 = t7 * t485 * t196;
1345   t1244 = t7 * t485 * t203;
1346   t1247 = t7 * t505 * t255;
1347   t1253 = t7 * t192 * t519;
1348   t1256 = t7 * t192 * t523;
1349   t1259 = t7 * t192 * t532;
1350   t1273 = t527 * t9;
1351   t1275 = t8 / t1273;
1352   t1277 = -0.12e2 * t528 + 0.12e2 * t1275;
1353   t1283 = t7 * t250 * t544;
1354   t1290 = t13 * t13;
1355   t1292 = 0.1e1 / t184 / t1290;
1356   t1293 = t1292 * t16;
1357   t1294 = t187 * t187;
1358   t1299 = t470 * t80;
1359   t1301 = t7 * t1299 * t473;
1360   t1303 = t185 * t208;
1361   t1305 = t7 * t1303 * t188;
1362   t1307 = -0.20e2 / 0.27e2 * t1231 + 0.4e1 / 0.9e1 * t7 * t471 * t101 * t472 - 0.2e1 / 0.3e1 * t1238 + 0.4e1 / 0.3e1 * t1241 + 0.2e1 / 0.3e1 * t1244 + 0.2e1 / 0.3e1 * t1247 - t7 * t186 * t238 * t187 - 0.2e1 * t1253 - 0.2e1 * t1256 - 0.2e1 / 0.3e1 * t1259 - 0.2e1 * t7 * t69 * t585 * t74 - 0.3e1 * t7 * t69 * t238 * t202 - 0.2e1 * t7 * t69 * t101 * t531 - t7 * t69 * t37 * t1277 / 0.2e1 - 0.2e1 / 0.3e1 * t1283 + t7 * t105 * t63 * t1277 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1293 * t37 * t1294 + 0.4e1 / 0.27e2 * t1301 + 0.2e1 / 0.9e1 * t1305;
1363   t1310 = -0.6e1 * t521 - 0.6e1 * t525 - 0.2e1 * t534 - t1118 - 0.2e1 * t541 - 0.2e1 * t546 - t1121 + t553 - 0.3e1 / 0.2e1 * t556 - 0.3e1 / 0.2e1 * t587 + t9 * (t1229 + t1307);
1364   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1365     v4rho4[0] = t1114 + t1310;
1366 
1367   t1319 = t185 * t280;
1368   t1328 = t7 * t281 * t119;
1369   t1331 = 0.1e1 / t79 / t527;
1370   t1334 = t7 * t14 * t1331 * t282;
1371   t1336 = t243 * t280;
1372   t1341 = t7 * t298 * t124;
1373   t1348 = t7 * t497 * t667;
1374   t1351 = t7 * t505 * t304;
1375   t1352 = t1351 / 0.3e1;
1376   t1366 = t266 * t208 * t37 * t268;
1377   t1367 = t1366 / 0.9e1;
1378   t1368 = t656 * t37;
1379   t1370 = t266 * t1368 * t609;
1380   t1372 = t614 * t290 * t631 * t106 / 0.3e1 - t7 * t69 * t585 * t118 / 0.2e1 - t7 * t1319 * t75 / 0.2e1 - t7 * t105 * t144 * t543 / 0.2e1 - 0.5e1 / 0.27e2 * t1328 - 0.25e2 / 0.6e1 * t1334 + t7 * t1336 * t107 / 0.2e1 - 0.5e1 / 0.27e2 * t1341 + t7 * t491 * t144 * t492 / 0.9e1 - t1348 / 0.6e1 + t1352 - t266 * t16 * t238 * t268 / 0.2e1 - t266 * t623 * t627 / 0.2e1 - t266 * t267 * t118 * t531 / 0.6e1 + t1367 + 0.7e1 / 0.6e1 * t1370;
1381   t1383 = t289 * t208 * t63 * t291;
1382   t1384 = t1383 / 0.9e1;
1383   t1393 = t266 * t80 * t101 * t268;
1384   t1395 = t280 * t101;
1385   t1399 = t266 * t604 * t627;
1386   t1406 = t289 * t619 * t631;
1387   t1412 = t656 * t63;
1388   t1414 = t289 * t1412 * t636;
1389   t1417 = t4 * t6 * t1292;
1390   t1423 = t599 * t604 * t600;
1391   t1430 = t4 * t6 * t1209;
1392   t1435 = -t289 * t290 * t123 * t543 / 0.6e1 + t289 * t635 * t8 * t254 / 0.2e1 + t1384 - t289 * t720 * t501 / 0.2e1 + t599 * t623 * t600 / 0.3e1 - t1393 / 0.3e1 - t266 * t1395 * t609 - t1399 / 0.6e1 - t266 * t608 * t8 * t202 / 0.2e1 - t1406 / 0.6e1 - t614 * t635 * t8 * t245 / 0.3e1 - 0.7e1 / 0.6e1 * t1414 - 0.5e1 / 0.27e2 * t1417 * t267 * t118 * t472 + t1423 / 0.9e1 + t599 * t608 * t8 * t187 / 0.3e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t123 * t492;
1393   t1438 = t614 * t619 * t615;
1394   t1446 = t7 * t657 * t37;
1395   t1452 = t7 * t681 * t63;
1396   t1455 = t7 * t548 * t144;
1397   t1458 = t7 * t485 * t275;
1398   t1459 = t1458 / 0.3e1;
1399   t1461 = t7 * t657 * t661;
1400   t1469 = t7 * t42 * t1331 * t299;
1401   t1471 = t1438 / 0.9e1 - t1129 / 0.6e1 - t1143 / 0.6e1 + t1152 + t1156 - 0.5e1 / 0.12e2 * t1158 + t1161 / 0.4e1 - t1164 / 0.8e1 + 0.5e1 / 0.2e1 * t1446 - 0.3e1 / 0.2e1 * t7 * t281 * t101 - 0.5e1 / 0.2e1 * t1452 - 0.5e1 / 0.36e2 * t1455 + t1459 + 0.7e1 / 0.2e1 * t1461 - 0.3e1 / 0.2e1 * t7 * t281 * t238 * t8 + 0.25e2 / 0.6e1 * t1469;
1402   t1473 = t7 * t250 * t685;
1403   t1476 = t7 * t192 * t677;
1404   t1496 = -t1473 / 0.2e1 - t1476 / 0.2e1 - 0.5e1 / 0.9e1 * t1205 + t1218 / 0.27e2 + t1222 / 0.9e1 - 0.5e1 / 0.9e1 * t1231 - t1238 / 0.6e1 + 0.2e1 / 0.3e1 * t1241 + t1244 / 0.3e1 + t1247 / 0.3e1 - t1253 / 0.2e1 - t1256 / 0.2e1 - t1259 / 0.6e1 - t1283 / 0.6e1 + t1301 / 0.27e2 + t1305 / 0.9e1 + t599 * t267 * t627 * t74 / 0.3e1;
1405   t1501 = t606 / 0.3e1;
1406   t1503 = t621 / 0.3e1;
1407   t1506 = t9 * (t1372 + t1435 + t1471 + t1496) + t602 / 0.3e1 - t1501 - t611 + t617 / 0.3e1 - t1503 - t625 - t629 / 0.2e1 - t633 / 0.2e1 + t638 - t484;
1408   t1512 = -t504 + 0.3e1 / 0.4e1 * t553 - 0.3e1 / 0.4e1 * t556 - t588 - 0.3e1 / 0.2e1 * t643 + 0.3e1 / 0.2e1 * t646 + t832 - t1118 - t1121 + t796 + 0.7e1 / 0.2e1 * t659;
1409   t1520 = -0.3e1 * t663 + t797 - t669 / 0.2e1 - t672 - t675 - 0.3e1 / 0.2e1 * t679 - 0.7e1 / 0.2e1 * t683 - 0.3e1 / 0.2e1 * t687 + t476 - t479 / 0.3e1 + t487;
1410   t1523 = t496 - t499 / 0.3e1 + t507 - t512 - 0.2e1 * t514 - t517 - t522 - t526 - t535 - t541 - t547;
1411   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1412     v4rho4[1] = t1506 + t1512 + t1520 + t1523;
1413 
1414   t1537 = t7 * t681 * t764;
1415   t1540 = t7 * t250 * t768;
1416   t1547 = t7 * t250 * t772;
1417   t1553 = t7 * t505 * t331;
1418   t1556 = t7 * t505 * t335;
1419   t1559 = 0.1e1 / t79 / t1273;
1420   t1561 = t8 * t8;
1421   t1567 = t7 * t477 * t750;
1422   t1570 = t7 * t192 * t754;
1423   t1573 = t7 * t192 * t760;
1424   t1596 = -0.4e1 / 0.9e1 * t614 * t635 * t291 * t8 + 0.4e1 / 0.9e1 * t599 * t608 * t268 * t8 - t7 * t1319 * t119 / 0.3e1 - 0.7e1 / 0.3e1 * t1537 - t1540 / 0.3e1 - t7 * t244 * t362 * t245 / 0.6e1 - t1547 / 0.3e1 + t7 * t1336 * t124 / 0.3e1 + 0.2e1 / 0.9e1 * t1553 + t1556 / 0.9e1 - t7 * t185 * t1559 * t37 * t1561 / 0.3e1 - t1567 / 0.9e1 - t1570 / 0.3e1 - t1573 / 0.3e1 - t7 * t186 * t238 * t314 / 0.6e1 - t7 * t69 * t238 * t320 / 0.2e1 - t7 * t69 * t101 * t759 - t7 * t105 * t362 * t254 / 0.2e1 - t7 * t243 * t1559 * t63 * t1561 / 0.3e1 - 0.10e2 / 0.27e2 * t1328;
1425   t1602 = t7 * t1303 * t315;
1426   t1605 = t7 * t485 * t321;
1427   t1608 = t7 * t1220 * t326;
1428   t1615 = t599 * t604 * t728;
1429   t1618 = t266 * t1368 * t712;
1430   t1625 = t266 * t604 * t704;
1431   t1632 = t614 * t619 * t708;
1432   t1639 = t289 * t619 * t724;
1433   t1655 = -0.41e2 / 0.9e1 * t1334 - 0.10e2 / 0.27e2 * t1341 - t1348 / 0.9e1 + 0.4e1 / 0.9e1 * t1351 + t1602 / 0.27e2 + t1605 / 0.9e1 + t1608 / 0.27e2 - 0.5e1 / 0.27e2 * t1417 * t267 * t314 * t187 + 0.2e1 / 0.27e2 * t1615 + 0.7e1 / 0.9e1 * t1618 + t599 * t267 * t320 * t187 / 0.9e1 - t1625 / 0.9e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t325 * t245 + 0.2e1 / 0.27e2 * t1632 + t614 * t290 * t334 * t245 / 0.9e1 - t1639 / 0.9e1 + 0.2e1 / 0.9e1 * t599 * t623 * t728 + t599 * t267 * t314 * t202 / 0.9e1 - t266 * t623 * t704 / 0.3e1 - t266 * t267 * t759 * t74 / 0.3e1;
1434   t1666 = t289 * t1412 * t716;
1435   t1671 = t280 * t144;
1436   t1694 = t289 * t80 * t144 * t291;
1437   t1704 = -t266 * t267 * t320 * t202 / 0.6e1 + t614 * t290 * t325 * t254 / 0.9e1 - 0.7e1 / 0.9e1 * t1666 - 0.2e1 / 0.3e1 * t266 * t1395 * t712 + 0.2e1 / 0.3e1 * t289 * t1671 * t636 - t289 * t290 * t767 * t106 / 0.3e1 - t289 * t720 * t631 / 0.3e1 - t289 * t290 * t334 * t254 / 0.6e1 + 0.4e1 / 0.27e2 * t1366 - 0.2e1 / 0.9e1 * t1370 + 0.4e1 / 0.27e2 * t1383 + 0.2e1 / 0.9e1 * t614 * t720 * t615 - 0.2e1 / 0.9e1 * t1694 - 0.2e1 / 0.9e1 * t1393 - t1399 / 0.9e1 - t1406 / 0.9e1 + 0.2e1 / 0.9e1 * t1414 + 0.2e1 / 0.27e2 * t1423 + 0.2e1 / 0.27e2 * t1438 + t7 * t298 * t144;
1438   t1706 = t7 * t259 * t362;
1439   t1725 = t1706 / 0.12e2 + t1152 + t1156 - 0.5e1 / 0.18e2 * t1158 + t1161 / 0.12e2 - t1446 / 0.3e1 + t1452 / 0.3e1 - 0.5e1 / 0.18e2 * t1455 + 0.4e1 / 0.9e1 * t1458 - 0.2e1 / 0.3e1 * t1461 + 0.41e2 / 0.9e1 * t1469 - t1473 / 0.3e1 - t1476 / 0.3e1 - 0.10e2 / 0.27e2 * t1205 + t1222 / 0.27e2 - 0.10e2 / 0.27e2 * t1231 + 0.2e1 / 0.9e1 * t1241 + t1244 / 0.9e1 + t1247 / 0.9e1 + t1305 / 0.27e2;
1440   t1741 = t9 * (t1596 + t1655 + t1704 + t1725) - t706 / 0.3e1 + 0.2e1 / 0.9e1 * t710 - 0.2e1 / 0.3e1 * t714 + 0.2e1 / 0.3e1 * t718 - 0.2e1 / 0.3e1 * t722 - t726 / 0.3e1 + 0.2e1 / 0.9e1 * t602 - 0.4e1 / 0.9e1 * t606 - 0.2e1 / 0.3e1 * t611 + 0.2e1 / 0.9e1 * t617 - 0.4e1 / 0.9e1 * t621 + 0.2e1 / 0.9e1 * t730;
1441   t1751 = -0.2e1 / 0.3e1 * t625 - t629 / 0.3e1 - t633 / 0.3e1 + 0.2e1 / 0.3e1 * t638 + t553 / 0.2e1 - t556 / 0.4e1 - t643 + t646 + t649 / 0.2e1 - t735 / 0.4e1 - t1118 - t1121 - t738 / 0.3e1;
1442   t1762 = -t741 / 0.9e1 - t745 / 0.3e1 - t748 / 0.9e1 - t752 / 0.3e1 - t756 - t762 + 0.2e1 * t766 - t770 - t774 - 0.2e1 / 0.3e1 * t777 + 0.2e1 / 0.3e1 * t653 + 0.5e1 / 0.3e1 * t659 - 0.2e1 * t663;
1443   t1775 = 0.2e1 / 0.3e1 * t665 - t669 / 0.3e1 - 0.4e1 / 0.3e1 * t672 - 0.4e1 / 0.3e1 * t675 - t679 - 0.5e1 / 0.3e1 * t683 - t687 - t479 / 0.9e1 + 0.2e1 / 0.3e1 * t487 - t499 / 0.9e1 + 0.2e1 / 0.3e1 * t507 - 0.2e1 / 0.3e1 * t514 - t517 / 0.3e1 - t541 / 0.3e1;
1444   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1445     v4rho4[2] = t1741 + t1751 + t1762 + t1775;
1446 
1447   t1790 = t7 * t192 * t837;
1448   t1793 = t7 * t497 * t841;
1449   t1796 = t7 * t250 * t845;
1450   t1799 = t7 * t250 * t849;
1451   t1802 = t7 * t250 * t854;
1452   t1805 = t7 * t1299 * t860;
1453   t1808 = t7 * t1216 * t869;
1454   t1816 = 0.6e1 * t528 + 0.12e2 * t1275;
1455   t1838 = t599 * t267 * t864 * t74 / 0.3e1 + t614 * t290 * t873 * t106 / 0.3e1 - t7 * t105 * t828 * t106 / 0.2e1 - t1790 / 0.6e1 - t1793 / 0.6e1 - t1796 / 0.2e1 - t1799 / 0.2e1 - t1802 / 0.6e1 + t1805 / 0.27e2 + t1808 / 0.27e2 - t7 * t69 * t101 * t836 / 0.2e1 - t7 * t69 * t37 * t1816 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t298 * t362 * t8 - 0.3e1 / 0.2e1 * t7 * t105 * t144 * t767 + t7 * t105 * t63 * t1816 / 0.2e1 + t7 * t471 * t101 * t859 / 0.9e1;
1456   t1858 = t1537 - t1540 / 0.2e1 - t1547 / 0.2e1 + 0.2e1 / 0.3e1 * t1553 + t1556 / 0.3e1 - t1567 / 0.6e1 - t1570 / 0.2e1 - t1573 / 0.2e1 - 0.5e1 / 0.9e1 * t1328 + t1334 / 0.3e1 - 0.5e1 / 0.9e1 * t1341 + t1352 + t1602 / 0.9e1 + t1605 / 0.3e1 + t1608 / 0.9e1 + t289 * t1671 * t716 - t289 * t16 * t362 * t291 / 0.2e1;
1457   t1899 = t289 * t619 * t873;
1458   t1906 = t266 * t604 * t864;
1459   t1917 = -t289 * t720 * t724 / 0.2e1 - t289 * t290 * t853 * t106 / 0.6e1 - 0.5e1 / 0.27e2 * t1417 * t267 * t859 * t74 - t266 * t623 * t864 / 0.2e1 - t266 * t267 * t118 * t759 / 0.2e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t868 * t106 - t614 * t635 * t325 * t8 / 0.3e1 - t289 * t290 * t123 * t767 / 0.2e1 + t289 * t635 * t8 * t334 / 0.2e1 - t266 * t608 * t8 * t320 / 0.2e1 - t1899 / 0.6e1 + t599 * t608 * t314 * t8 / 0.3e1 - t1906 / 0.6e1 - t266 * t267 * t836 * t74 / 0.6e1 + t614 * t720 * t708 / 0.3e1 + t1615 / 0.9e1 - t1618 / 0.3e1;
1460   t1925 = t7 * t111 * t828;
1461   t1932 = -t1625 / 0.6e1 + t1632 / 0.9e1 - t1639 / 0.6e1 + t1666 / 0.3e1 + t1367 + t1384 - t1694 / 0.3e1 + t1706 / 0.4e1 - t1925 / 0.8e1 + t1152 + t1156 - 0.5e1 / 0.36e2 * t1158 - 0.5e1 / 0.12e2 * t1455 + t1459 - t1469 / 0.3e1 - 0.5e1 / 0.27e2 * t1205 - 0.5e1 / 0.27e2 * t1231;
1462   t1939 = t9 * (t1838 + t1858 + t1917 + t1932) - t867 - t876 - t706 / 0.2e1 + t710 / 0.3e1 - t714 + t718 - t722 - t726 / 0.2e1 - t1501 - t1503;
1463   t1943 = t730 / 0.3e1 + t554 + 0.3e1 / 0.4e1 * t649 - 0.3e1 / 0.4e1 * t735 - t831 - t1118 - t1121 - t738 - t840 - t844 - t848;
1464   t1951 = -t852 - t857 + t863 + t872 - t741 / 0.3e1 - t745 - t748 / 0.3e1 - t752 / 0.2e1 - 0.3e1 / 0.2e1 * t756 - 0.3e1 / 0.2e1 * t762 + 0.3e1 * t766;
1465   t1955 = -0.3e1 / 0.2e1 * t770 - 0.3e1 / 0.2e1 * t774 - 0.2e1 * t777 + t653 - t659 + t665 - t672 - t675 + t683 + t488 + t508;
1466   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1467     v4rho4[3] = t1939 + t1943 + t1951 + t1955;
1468 
1469   t1986 = t320 * t320;
1470   t1991 = t325 * t325;
1471   t1996 = t334 * t334;
1472   t2001 = t314 * t314;
1473   t2007 = 0.12e2 * t528 + 0.12e2 * t1275;
1474   t2035 = 0.2e1 / 0.3e1 * t614 * t290 * t325 * t334 - 0.2e1 / 0.3e1 * t289 * t290 * t853 * t123 + 0.2e1 / 0.3e1 * t599 * t267 * t314 * t320 - 0.2e1 / 0.3e1 * t266 * t267 * t836 * t118 - 0.2e1 * t289 * t720 * t873 + 0.4e1 / 0.9e1 * t7 * t491 * t144 * t868 - t7 * t186 * t37 * t1986 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1210 * t63 * t1991 - t7 * t244 * t63 * t1996 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1293 * t37 * t2001 - t7 * t69 * t37 * t2007 / 0.2e1 - t7 * t244 * t362 * t325 - 0.3e1 * t7 * t105 * t362 * t334 - 0.2e1 * t7 * t105 * t144 * t853 + t7 * t105 * t63 * t2007 / 0.2e1 - 0.2e1 * t7 * t105 * t828 * t123 - 0.2e1 / 0.3e1 * t1790 - 0.2e1 / 0.3e1 * t1793 - 0.2e1 * t1796;
1475   t2050 = 0.1e1 / t48 / t355;
1476   t2063 = t820 * rho[1];
1477   t2074 = 0.1e1 / t823 / t59;
1478   t2085 = -0.2e1 * t1799 - 0.2e1 / 0.3e1 * t1802 + 0.4e1 / 0.27e2 * t1805 + 0.4e1 / 0.27e2 * t1808 + 0.4e1 / 0.3e1 * t1553 + 0.2e1 / 0.3e1 * t1556 - 0.20e2 / 0.27e2 * t1328 - 0.20e2 / 0.27e2 * t1341 + 0.2e1 / 0.9e1 * t1602 + 0.2e1 / 0.3e1 * t1605 + 0.2e1 / 0.9e1 * t1608 - 0.2e1 / 0.3e1 * t1899 - 0.2e1 / 0.3e1 * t1906 - 0.3e1 / 0.8e1 * t7 * t44 * (0.35912208504801097395e1 * t45 * t2050 + 0.3640e4 / 0.81e2 * t54 * t1094 * t60 - 0.23072e5 / 0.27e2 * t138 * t2050 * t140 + 0.476672e6 / 0.81e2 * t354 / t820 * t358 - 0.475136e6 / 0.27e2 * t819 / t47 / t2063 * t824 + 0.524288e6 / 0.27e2 * params->cc * t1084 / t48 / t820 / t46 * t2074) + t1706 / 0.2e1 - t1925 / 0.2e1 + t1152 + t1156 - 0.5e1 / 0.9e1 * t1455;
1479   t2089 = 0.4e1 / 0.3e1 * t653 + 0.4e1 / 0.3e1 * t665 - 0.2e1 / 0.3e1 * t741 - 0.2e1 * t745 - 0.2e1 / 0.3e1 * t748 - t1118 - t1121 - 0.3e1 / 0.2e1 * t830 + t9 * (t2035 + t2085) + t649 - 0.3e1 / 0.2e1 * t735;
1480   t2101 = -0.2e1 * t839 - 0.2e1 * t843 - 0.6e1 * t847 - 0.6e1 * t851 - 0.2e1 * t856 - 0.4e1 * t777 - 0.2e1 * t738 + 0.4e1 / 0.9e1 * t862 - 0.2e1 * t866 + 0.4e1 / 0.9e1 * t871 - 0.2e1 * t875;
1481   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1482     v4rho4[4] = t2089 + t2101;
1483 
1484   t2103 = t7 * t186 * t888;
1485   t2106 = t7 * t192 * t372;
1486   t2124 = t7 * t69 * t892;
1487   t2127 = t7 * t69 * t895;
1488   t2131 = t7 * t209 * t163 / 0.9e1;
1489   t2141 = t7 * t81 * t389;
1490   t2144 = t7 * t17 * t915;
1491   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1492     v4rho3sigma[0] = -0.2e1 / 0.3e1 * t2103 - 0.2e1 / 0.3e1 * t2106 + t152 * t598 * t163 * t472 / 0.9e1 - t940 * t941 * t481 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t371 * t389 * t202 - t152 * t371 * t163 * t531 / 0.2e1 - 0.4e1 * t2124 - 0.2e1 * t2127 + t2131 - t152 * t265 * t389 * t187 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t371 * t915 * t74 - t2141 / 0.2e1 - 0.3e1 / 0.2e1 * t2144 - 0.3e1 / 0.8e1 * t152 * t153 * (-0.63374485596707818932e0 * params->bb * t559 - 0.140e3 / 0.27e2 * t157 * t564 + 0.1456e4 / 0.9e1 * params->cc * t559 * t97 - 0.40576e5 / 0.27e2 * params->cc * t571 * t386 + 0.50176e5 / 0.9e1 * params->cc * t579 * t912 - 0.65536e5 / 0.9e1 * params->cc / t21 / t1180 * t1191 * t229);
1493 
1494   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1495     v4rho3sigma[1] = 0.0e0;
1496 
1497   t2170 = t7 * t244 * t925;
1498   t2173 = t7 * t250 * t396;
1499   t2176 = t7 * t105 * t929;
1500   t2180 = t7 * t259 * t177 / 0.9e1;
1501   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1502     v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2170 - 0.2e1 / 0.3e1 * t2173 - 0.2e1 * t2176 + t2180 + t152 * t613 * t177 * t492 / 0.9e1 - t940 * t958 * t501 / 0.2e1 - t152 * t395 * t177 * t543 / 0.2e1;
1503 
1504   t2193 = t280 * t1 * t3;
1505   t2194 = t2193 * t952;
1506   t2207 = t80 * t1 * t3;
1507   t2208 = t2207 * t401;
1508   t2210 = t935 * t939;
1509   t2212 = 0.4e1 / 0.9e1 * t2210 * t942;
1510   t2223 = 0.4e1 / 0.3e1 * t936 * t946;
1511   t2226 = t470 * t163;
1512   t2230 = t949 * t939;
1513   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1514     v4rho3sigma[3] = t2194 / 0.6e1 - t940 * t185 * t389 * t268 / 0.3e1 - t940 * t941 * t627 / 0.6e1 - t950 * t371 * t389 * t8 - 0.2e1 / 0.9e1 * t2208 - t2212 - 0.2e1 / 0.9e1 * t2103 - 0.4e1 / 0.9e1 * t2106 - t152 * t371 * t915 * t118 / 0.2e1 - t950 * t371 * t163 / 0.2e1 - t2223 - 0.4e1 / 0.3e1 * t2124 - 0.2e1 / 0.3e1 * t2127 + t2131 + t940 * t2226 * t600 / 0.9e1 - t2230 * t941 * t609 / 0.3e1 - t2141 / 0.3e1 - t2144 / 0.2e1;
1515 
1516   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1517     v4rho3sigma[4] = 0.0e0;
1518 
1519   t2239 = t2207 * t406;
1520   t2242 = 0.4e1 / 0.9e1 * t2210 * t959;
1521   t2243 = t2193 * t963;
1522   t2245 = t490 * t177;
1523   t2258 = t2207 * t423;
1524   t2261 = 0.4e1 / 0.3e1 * t936 * t969;
1525   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1526     v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2170 - 0.4e1 / 0.9e1 * t2173 - 0.2e1 / 0.3e1 * t2176 + t2180 - 0.2e1 / 0.9e1 * t2239 - t2242 - t2243 / 0.6e1 + t940 * t2245 * t615 / 0.9e1 + t2230 * t958 * t636 / 0.3e1 - t940 * t958 * t631 / 0.6e1 + t950 * t395 * t177 / 0.2e1 - t2258 / 0.6e1 - t2261 - t152 * t288 * t422 * t245 / 0.6e1 - t152 * t395 * t422 * t254 / 0.2e1;
1527 
1528   t2274 = t936 * t974;
1529   t2286 = t936 * t978;
1530   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1531     v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2208 - t2212 - t2223 - 0.4e1 / 0.3e1 * t2194 - 0.2e1 / 0.9e1 * t2106 + t2131 - t2141 / 0.6e1 - 0.2e1 / 0.9e1 * t2274 + t940 * t2226 * t728 / 0.9e1 - t152 * t265 * t389 * t314 / 0.6e1 - t2230 * t941 * t712 / 0.3e1 - 0.2e1 / 0.3e1 * t2286 - t940 * t941 * t704 / 0.6e1 - t152 * t371 * t389 * t320 / 0.2e1 - t152 * t371 * t163 * t759 / 0.2e1;
1532 
1533   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1534     v4rho3sigma[7] = 0.0e0;
1535 
1536   t2303 = t936 * t983;
1537   t2311 = t936 * t987;
1538   t2320 = t936 * t990;
1539   t2329 = t936 * t1010;
1540   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1541     v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2239 - t2242 + 0.4e1 / 0.3e1 * t2243 - 0.2e1 / 0.9e1 * t2173 + t2180 - t2258 / 0.3e1 - t2261 - 0.2e1 / 0.9e1 * t2303 + t940 * t2245 * t708 / 0.9e1 + t2230 * t958 * t716 / 0.3e1 - 0.4e1 / 0.3e1 * t2311 - t940 * t243 * t422 * t291 / 0.3e1 + t950 * t395 * t422 * t8 - 0.2e1 / 0.3e1 * t2320 - t940 * t958 * t724 / 0.6e1 - t152 * t395 * t177 * t767 / 0.2e1 - t2329 / 0.2e1 - t152 * t395 * t1009 * t106 / 0.2e1;
1542 
1543   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1544     v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2208 - 0.2e1 / 0.3e1 * t2274 - 0.2e1 * t2286 + t2131 + t152 * t598 * t163 * t859 / 0.9e1 - t940 * t941 * t864 / 0.2e1 - t152 * t371 * t163 * t836 / 0.2e1;
1545 
1546   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1547     v4rho3sigma[10] = 0.0e0;
1548 
1549   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1550     v4rho3sigma[11] = -t152 * t288 * t422 * t325 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t395 * t1009 * t123 - 0.3e1 / 0.2e1 * t152 * t395 * t422 * t334 - t152 * t395 * t177 * t853 / 0.2e1 - t940 * t958 * t873 / 0.2e1 + t152 * t613 * t177 * t868 / 0.9e1 - 0.2e1 / 0.3e1 * t2303 - 0.4e1 * t2311 - 0.2e1 * t2320 - 0.2e1 / 0.3e1 * t2239 - 0.3e1 / 0.2e1 * t2329 - t2258 / 0.2e1 + t2180 - 0.3e1 / 0.8e1 * t152 * t167 * (-0.63374485596707818932e0 * params->bb * t802 - 0.140e3 / 0.27e2 * t171 * t807 + 0.1456e4 / 0.9e1 * params->cc * t802 * t140 - 0.40576e5 / 0.27e2 * params->cc * t814 * t419 + 0.50176e5 / 0.9e1 * params->cc * t822 * t1006 - 0.65536e5 / 0.9e1 * params->cc / t48 / t2063 * t2074 * t353);
1551 
1552   t2402 = t2207 * t440 / 0.6e1;
1553   t2403 = t936 * t1016;
1554   t2405 = t936 * t1032;
1555   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556     v4rho2sigma2[0] = -t2402 - 0.4e1 / 0.3e1 * t2403 - t2405 - t152 * t265 * t439 * t187 / 0.6e1 - t152 * t371 * t1031 * t74 - t152 * t371 * t439 * t202 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * (-0.7e1 / 0.9e1 * t427 * t223 - 0.28e2 / 0.9e1 * t431 * t226 + 0.2048e4 / 0.9e1 * t157 * t235 - 0.4480e4 / 0.3e1 * t911 * t581 + 0.8192e4 / 0.3e1 * params->cc / t21 / t577 * t1191 * t26);
1557 
1558   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559     v4rho2sigma2[1] = 0.0e0;
1560 
1561   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1562     v4rho2sigma2[2] = 0.0e0;
1563 
1564   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1565     v4rho2sigma2[3] = 0.0e0;
1566 
1567   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1568     v4rho2sigma2[4] = 0.0e0;
1569 
1570   t2436 = t2207 * t457 / 0.6e1;
1571   t2437 = t936 * t1037;
1572   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1573     v4rho2sigma2[5] = -t2436 - 0.4e1 / 0.3e1 * t2437 - t152 * t288 * t456 * t245 / 0.6e1 - t152 * t395 * t456 * t254 / 0.2e1;
1574 
1575   t2449 = t936 * t1041;
1576   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1577     v4rho2sigma2[6] = -t2402 - 0.2e1 / 0.3e1 * t2403 - t2405 / 0.2e1 - 0.2e1 / 0.3e1 * t2449 - t940 * t185 * t439 * t268 / 0.6e1 - t152 * t371 * t1031 * t118 / 0.2e1 - t950 * t371 * t439 * t8 / 0.2e1;
1578 
1579   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1580     v4rho2sigma2[7] = 0.0e0;
1581 
1582   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1583     v4rho2sigma2[8] = 0.0e0;
1584 
1585   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1586     v4rho2sigma2[9] = 0.0e0;
1587 
1588   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1589     v4rho2sigma2[10] = 0.0e0;
1590 
1591   t2464 = t936 * t1046;
1592   t2474 = t936 * t1062;
1593   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594     v4rho2sigma2[11] = -t2436 - 0.2e1 / 0.3e1 * t2437 - 0.2e1 / 0.3e1 * t2464 - t940 * t243 * t456 * t291 / 0.6e1 + t950 * t395 * t456 * t8 / 0.2e1 - t2474 / 0.2e1 - t152 * t395 * t1061 * t106 / 0.2e1;
1595 
1596   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597     v4rho2sigma2[12] = -t2402 - 0.4e1 / 0.3e1 * t2449 - t152 * t265 * t439 * t314 / 0.6e1 - t152 * t371 * t439 * t320 / 0.2e1;
1598 
1599   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600     v4rho2sigma2[13] = 0.0e0;
1601 
1602   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603     v4rho2sigma2[14] = 0.0e0;
1604 
1605   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606     v4rho2sigma2[15] = 0.0e0;
1607 
1608   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1609     v4rho2sigma2[16] = 0.0e0;
1610 
1611   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1612     v4rho2sigma2[17] = -t2436 - 0.4e1 / 0.3e1 * t2464 - t2474 - t152 * t288 * t456 * t325 / 0.6e1 - t152 * t395 * t1061 * t123 - t152 * t395 * t456 * t334 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t167 * (-0.7e1 / 0.9e1 * t444 * t347 - 0.28e2 / 0.9e1 * t448 * t350 + 0.2048e4 / 0.9e1 * t171 * t359 - 0.4480e4 / 0.3e1 * t1005 * t824 + 0.8192e4 / 0.3e1 * params->cc / t48 / t820 * t2074 * t53);
1613 
1614   t2519 = t936 * t1081;
1615   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1616     v4rhosigma3[0] = -t2519 / 0.2e1 - t152 * t371 * t1080 * t74 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * (-t1067 * t92 / 0.2e1 - 0.2e1 * t1071 * t98 + 0.16e2 * t427 * t1023 + 0.256e3 * t1028 * t1077 - 0.1024e4 * params->cc / t21 / t570 * t1191 * t156);
1617 
1618   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1619     v4rhosigma3[1] = 0.0e0;
1620 
1621   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1622     v4rhosigma3[2] = 0.0e0;
1623 
1624   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1625     v4rhosigma3[3] = 0.0e0;
1626 
1627   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1628     v4rhosigma3[4] = 0.0e0;
1629 
1630   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1631     v4rhosigma3[5] = 0.0e0;
1632 
1633   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1634     v4rhosigma3[6] = 0.0e0;
1635 
1636   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1637     v4rhosigma3[7] = 0.0e0;
1638 
1639   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1640     v4rhosigma3[8] = 0.0e0;
1641 
1642   t2543 = t936 * t1100;
1643   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1644     v4rhosigma3[9] = -t152 * t395 * t1099 * t106 / 0.2e1 - t2543 / 0.2e1;
1645 
1646   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1647     v4rhosigma3[10] = -t152 * t371 * t1080 * t118 / 0.2e1 - t2519 / 0.2e1;
1648 
1649   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1650     v4rhosigma3[11] = 0.0e0;
1651 
1652   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1653     v4rhosigma3[12] = 0.0e0;
1654 
1655   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1656     v4rhosigma3[13] = 0.0e0;
1657 
1658   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1659     v4rhosigma3[14] = 0.0e0;
1660 
1661   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1662     v4rhosigma3[15] = 0.0e0;
1663 
1664   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1665     v4rhosigma3[16] = 0.0e0;
1666 
1667   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1668     v4rhosigma3[17] = 0.0e0;
1669 
1670   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1671     v4rhosigma3[18] = 0.0e0;
1672 
1673   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1674     v4rhosigma3[19] = -t2543 / 0.2e1 - t152 * t395 * t1099 * t123 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t167 * (-t1086 * t135 / 0.2e1 - 0.2e1 * t1090 * t141 + 0.16e2 * t444 * t1053 + 0.256e3 * t1058 * t1096 - 0.1024e4 * params->cc / t48 / t813 * t2074 * t170);
1675 
1676   t2575 = t575 * sigma[0];
1677   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1678     v4sigma4[0] = -0.3e1 / 0.8e1 * t152 * t153 * (-0.15e2 / 0.16e2 * params->cc / t26 / t2575 * t35 - 0.15e2 / 0.4e1 * params->cc / t2575 * t432 - 0.6e1 * t1067 * t434 * t234 + 0.384e3 * params->cc * t1167 * t1191 * t426 + 0.48e2 * t1076 * t581 * t1070);
1679 
1680   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1681     v4sigma4[1] = 0.0e0;
1682 
1683   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1684     v4sigma4[2] = 0.0e0;
1685 
1686   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1687     v4sigma4[3] = 0.0e0;
1688 
1689   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1690     v4sigma4[4] = 0.0e0;
1691 
1692   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1693     v4sigma4[5] = 0.0e0;
1694 
1695   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1696     v4sigma4[6] = 0.0e0;
1697 
1698   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1699     v4sigma4[7] = 0.0e0;
1700 
1701   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1702     v4sigma4[8] = 0.0e0;
1703 
1704   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1705     v4sigma4[9] = 0.0e0;
1706 
1707   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1708     v4sigma4[10] = 0.0e0;
1709 
1710   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1711     v4sigma4[11] = 0.0e0;
1712 
1713   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1714     v4sigma4[12] = 0.0e0;
1715 
1716   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1717     v4sigma4[13] = 0.0e0;
1718 
1719   t2599 = t818 * sigma[2];
1720   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1721     v4sigma4[14] = -0.3e1 / 0.8e1 * t152 * t167 * (-0.15e2 / 0.16e2 * params->cc / t53 / t2599 * t61 - 0.15e2 / 0.4e1 * params->cc / t2599 * t449 - 0.6e1 * t1086 * t451 * t358 + 0.384e3 * params->cc * t2050 * t2074 * t443 + 0.48e2 * t1095 * t824 * t1089);
1722 
1723 #ifndef XC_DONT_COMPILE_MXC
1724 
1725   if(order < 5) return;
1726 
1727 
1728 #endif
1729 
1730 #endif
1731 
1732 #endif
1733 
1734 #endif
1735 
1736 #endif
1737 
1738 
1739 }
1740 
1741