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_optx.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, t15, t16, t17, t18;
25   double t23, t26, t28, t29, t30, t31, t34, t36;
26 
27 #ifndef XC_DONT_COMPILE_VXC
28   double t41, t42, t43, t47, t51, t52, t53, t54;
29   double t55, t58, t62, t66, t69, t75;
30 
31 #ifndef XC_DONT_COMPILE_FXC
32   double t78, t79, t86, t87, t91, t94, t95, t100;
33   double t101, t102, t103, t104, t106, t107, t108, t110;
34   double t113, t118, t127, t128, t130, t132, t136, t145;
35   double t149, t153;
36 
37 #ifndef XC_DONT_COMPILE_KXC
38   double t157, t158, t171, t174, t179, t183, t187, t188;
39   double t189, t190, t194, t196, t199, t204, t217, t222;
40   double t226, t233, t239, t240, t244, t248, t257, t261;
41   double t265;
42 
43 #ifndef XC_DONT_COMPILE_LXC
44   double t296, t304, t310;
45 #endif
46 
47 #endif
48 
49 #endif
50 
51 #endif
52 
53 #endif
54 
55 
56   gga_x_optx_params *params;
57 
58   assert(p->params != NULL);
59   params = (gga_x_optx_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->gamma * params->gamma;
71   t13 = params->b * t12;
72   t14 = sigma[0] * sigma[0];
73   t15 = t13 * t14;
74   t16 = rho[0] * rho[0];
75   t17 = t16 * t16;
76   t18 = t17 * rho[0];
77   t23 = t10 * t10;
78   t26 = t9 / t23 / t16;
79   t28 = params->gamma * sigma[0] * t26 + 0.1e1;
80   t29 = t28 * t28;
81   t30 = 0.1e1 / t29;
82   t31 = t8 / t10 / t18 * t30;
83   t34 = 0.2e1 * t15 * t31 + params->a;
84   t36 = t7 * t11 * t34;
85   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
86     zk[0] = -0.3e1 / 0.16e2 * t36;
87 
88 #ifndef XC_DONT_COMPILE_VXC
89 
90   if(order < 1) return;
91 
92 
93   t41 = t10 * rho[0] * t1 * t3;
94   t42 = t6 * t9;
95   t43 = t17 * t16;
96   t47 = t8 / t10 / t43 * t30;
97   t51 = params->b * t12 * params->gamma;
98   t52 = t14 * sigma[0];
99   t53 = t17 * t17;
100   t54 = t53 * rho[0];
101   t55 = 0.1e1 / t54;
102   t58 = 0.1e1 / t29 / t28;
103   t62 = -0.32e2 / 0.3e1 * t15 * t47 + 0.64e2 / 0.3e1 * t51 * t52 * t55 * t58;
104   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
105     vrho[0] = -t36 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t42 * t62;
106 
107   t66 = t13 * sigma[0];
108   t69 = 0.1e1 / t53;
109   t75 = t42 * (-0.8e1 * t51 * t14 * t69 * t58 + 0.4e1 * t66 * t31);
110   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
111     vsigma[0] = -0.3e1 / 0.16e2 * t41 * t75;
112 
113 #ifndef XC_DONT_COMPILE_FXC
114 
115   if(order < 2) return;
116 
117 
118   t78 = 0.1e1 / t23;
119   t79 = t9 * t78;
120   t86 = t16 * rho[0];
121   t87 = t17 * t86;
122   t91 = t8 / t10 / t87 * t30;
123   t94 = t53 * t16;
124   t95 = 0.1e1 / t94;
125   t100 = t12 * t12;
126   t101 = params->b * t100;
127   t102 = t14 * t14;
128   t103 = t101 * t102;
129   t104 = t53 * t17;
130   t106 = 0.1e1 / t23 / t104;
131   t107 = t29 * t29;
132   t108 = 0.1e1 / t107;
133   t110 = t106 * t108 * t9;
134   t113 = 0.608e3 / 0.9e1 * t15 * t91 - 0.2752e4 / 0.9e1 * t51 * t52 * t95 * t58 + 0.512e3 / 0.3e1 * t103 * t110;
135   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
136     v2rho2[0] = -t7 * t79 * t34 / 0.12e2 - t7 * t11 * t62 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t42 * t113;
137 
138   t118 = t10 * t1 * t3;
139   t127 = t101 * t52;
140   t128 = t53 * t86;
141   t130 = 0.1e1 / t23 / t128;
142   t132 = t130 * t108 * t9;
143   t136 = t42 * (-0.64e2 / 0.3e1 * t66 * t47 + 0.320e3 / 0.3e1 * t51 * t14 * t55 * t58 - 0.64e2 * t127 * t132);
144   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
145     v2rhosigma[0] = -t118 * t75 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t136;
146 
147   t145 = t101 * t14;
148   t149 = 0.1e1 / t23 / t94 * t108 * t9;
149   t153 = t42 * (-0.32e2 * t51 * sigma[0] * t69 * t58 + 0.4e1 * t13 * t31 + 0.24e2 * t145 * t149);
150   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
151     v2sigma2[0] = -0.3e1 / 0.16e2 * t41 * t153;
152 
153 #ifndef XC_DONT_COMPILE_KXC
154 
155   if(order < 3) return;
156 
157 
158   t157 = 0.1e1 / t23 / rho[0];
159   t158 = t9 * t157;
160   t171 = t8 / t10 / t53 * t30;
161   t174 = 0.1e1 / t128;
162   t179 = t53 * t18;
163   t183 = 0.1e1 / t23 / t179 * t108 * t9;
164   t187 = params->b * t100 * params->gamma;
165   t188 = t102 * sigma[0];
166   t189 = t187 * t188;
167   t190 = t53 * t53;
168   t194 = 0.1e1 / t107 / t28;
169   t196 = 0.1e1 / t10 / t190 * t194 * t8;
170   t199 = -0.13376e5 / 0.27e2 * t15 * t171 + 0.102016e6 / 0.27e2 * t51 * t52 * t174 * t58 - 0.4608e4 * t103 * t183 + 0.32768e5 / 0.9e1 * t189 * t196;
171   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
172     v3rho3[0] = t7 * t158 * t34 / 0.18e2 - t7 * t79 * t62 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t113 - 0.3e1 / 0.16e2 * t41 * t42 * t199;
173 
174   t204 = t78 * t1 * t3;
175   t217 = t187 * t102;
176   t222 = 0.1e1 / t10 / t53 / t87 * t194 * t8;
177   t226 = t42 * (0.1216e4 / 0.9e1 * t66 * t91 - 0.10688e5 / 0.9e1 * t51 * t14 * t95 * t58 + 0.1600e4 * t127 * t110 - 0.4096e4 / 0.3e1 * t217 * t222);
178   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
179     v3rho2sigma[0] = -t204 * t75 / 0.12e2 - t118 * t136 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t226;
180 
181   t233 = t55 * t58;
182   t239 = t187 * t52;
183   t240 = t53 * t43;
184   t244 = 0.1e1 / t10 / t240 * t194 * t8;
185   t248 = t42 * (-0.64e2 / 0.3e1 * t13 * t47 + 0.896e3 / 0.3e1 * t51 * t233 * sigma[0] - 0.512e3 * t145 * t132 + 0.512e3 * t239 * t244);
186   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
187     v3rhosigma2[0] = -t118 * t153 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t248;
188 
189   t257 = t187 * t14;
190   t261 = 0.1e1 / t10 / t179 * t194 * t8;
191   t265 = t42 * (0.144e3 * t101 * sigma[0] * t149 - 0.48e2 * t51 * t69 * t58 - 0.192e3 * t257 * t261);
192   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
193     v3sigma3[0] = -0.3e1 / 0.16e2 * t41 * t265;
194 
195 #ifndef XC_DONT_COMPILE_LXC
196 
197   if(order < 4) return;
198 
199 
200   t296 = t190 * rho[0];
201   t304 = params->b * t100 * t12;
202   t310 = 0.1e1 / t107 / t29;
203   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
204     v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t26 * t34 + 0.2e1 / 0.9e1 * t7 * t158 * t62 - t7 * t79 * t113 / 0.2e1 - t7 * t11 * t199 - 0.3e1 / 0.16e2 * t41 * t42 * (0.334400e6 / 0.81e2 * t15 * t8 / t10 / t54 * t30 - 0.3794560e7 / 0.81e2 * t51 * t52 / t104 * t58 + 0.2516480e7 / 0.27e2 * t103 / t23 / t240 * t108 * t9 - 0.4259840e7 / 0.27e2 * t189 / t10 / t296 * t194 * t8 + 0.2621440e7 / 0.27e2 * t304 * t102 * t14 / t190 / t17 * t310);
205 
206   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
207     v4rho3sigma[0] = t157 * t1 * t3 * t75 / 0.18e2 - t204 * t136 / 0.4e1 - 0.3e1 / 0.4e1 * t118 * t226 - 0.3e1 / 0.16e2 * t41 * t42 * (-0.26752e5 / 0.27e2 * t66 * t171 + 0.359552e6 / 0.27e2 * t51 * t14 * t174 * t58 - 0.267904e6 / 0.9e1 * t127 * t183 + 0.495616e6 / 0.9e1 * t217 * t196 - 0.327680e6 / 0.9e1 * t304 * t188 / t190 / t86 * t310);
208 
209   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
210     v4rho2sigma2[0] = -t204 * t153 / 0.12e2 - t118 * t248 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t42 * (0.1216e4 / 0.9e1 * t13 * t91 - 0.26240e5 / 0.9e1 * t51 * t95 * t58 * sigma[0] + 0.25088e5 / 0.3e1 * t101 * t106 * t108 * t14 * t9 - 0.54784e5 / 0.3e1 * t239 * t222 + 0.40960e5 / 0.3e1 * t304 * t102 / t190 / t16 * t310);
211 
212   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
213     v4rhosigma3[0] = -t118 * t265 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t42 * (0.384e3 * t51 * t233 - 0.1920e4 * t101 * t130 * t108 * sigma[0] * t9 + 0.5632e4 * t257 * t244 - 0.5120e4 * t304 * t52 / t296 * t310);
214 
215   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
216     v4sigma4[0] = -0.3e1 / 0.16e2 * t41 * t42 * (0.288e3 * t101 * t149 - 0.1536e4 * t187 * sigma[0] * t261 + 0.1920e4 * t304 * t14 / t190 * t310);
217 
218 #ifndef XC_DONT_COMPILE_MXC
219 
220   if(order < 5) return;
221 
222 
223 #endif
224 
225 #endif
226 
227 #endif
228 
229 #endif
230 
231 #endif
232 
233 
234 }
235 
236 
237 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 *))238 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
239 {
240 
241 #ifndef XC_DONT_COMPILE_EXC
242   double t1, t3, t4, t5, t6, t7, t8, t9;
243   double t10, t11, t12, t13, t14, t16, t19, t21;
244   double t23, t24, t25, t28, t30;
245 
246 #ifndef XC_DONT_COMPILE_VXC
247   double t34, t35, t36, t38, t43, t44, t45, t46;
248   double t47, t50, t54, t61, t67;
249 
250 #ifndef XC_DONT_COMPILE_FXC
251   double t70, t71, t77, t78, t80, t85, t86, t91;
252   double t92, t93, t94, t96, t98, t99, t103, t107;
253   double t118, t120, t126, t137, t143;
254 
255 #ifndef XC_DONT_COMPILE_KXC
256   double t147, t148, t159, t164, t169, t171, t177, t178;
257   double t179, t181, t184, t188, t192, t210, t216, t224;
258   double t232, t234, t240, t251, t257;
259 
260 #ifndef XC_DONT_COMPILE_LXC
261   double t289, t297, t303;
262 #endif
263 
264 #endif
265 
266 #endif
267 
268 #endif
269 
270 #endif
271 
272 
273   gga_x_optx_params *params;
274 
275   assert(p->params != NULL);
276   params = (gga_x_optx_params * )(p->params);
277 
278   t1 = M_CBRT3;
279   t3 = POW_1_3(0.1e1 / M_PI);
280   t4 = t1 * t3;
281   t5 = M_CBRT4;
282   t6 = t5 * t5;
283   t7 = POW_1_3(rho[0]);
284   t8 = t6 * t7;
285   t9 = params->gamma * params->gamma;
286   t10 = params->b * t9;
287   t11 = sigma[0] * sigma[0];
288   t12 = rho[0] * rho[0];
289   t13 = t12 * t12;
290   t14 = t13 * rho[0];
291   t16 = 0.1e1 / t7 / t14;
292   t19 = t7 * t7;
293   t21 = 0.1e1 / t19 / t12;
294   t23 = params->gamma * sigma[0] * t21 + 0.1e1;
295   t24 = t23 * t23;
296   t25 = 0.1e1 / t24;
297   t28 = t10 * t11 * t16 * t25 + params->a;
298   t30 = t4 * t8 * t28;
299   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
300     zk[0] = -0.3e1 / 0.8e1 * t30;
301 
302 #ifndef XC_DONT_COMPILE_VXC
303 
304   if(order < 1) return;
305 
306 
307   t34 = t7 * rho[0] * t1;
308   t35 = t3 * t6;
309   t36 = t13 * t12;
310   t38 = 0.1e1 / t7 / t36;
311   t43 = params->b * t9 * params->gamma;
312   t44 = t11 * sigma[0];
313   t45 = t13 * t13;
314   t46 = t45 * rho[0];
315   t47 = 0.1e1 / t46;
316   t50 = 0.1e1 / t24 / t23;
317   t54 = -0.16e2 / 0.3e1 * t10 * t11 * t38 * t25 + 0.16e2 / 0.3e1 * t43 * t44 * t47 * t50;
318   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
319     vrho[0] = -t30 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t35 * t54;
320 
321   t61 = 0.1e1 / t45;
322   t67 = t35 * (0.2e1 * t10 * sigma[0] * t16 * t25 - 0.2e1 * t43 * t11 * t61 * t50);
323   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
324     vsigma[0] = -0.3e1 / 0.8e1 * t34 * t67;
325 
326 #ifndef XC_DONT_COMPILE_FXC
327 
328   if(order < 2) return;
329 
330 
331   t70 = 0.1e1 / t19;
332   t71 = t6 * t70;
333   t77 = t12 * rho[0];
334   t78 = t13 * t77;
335   t80 = 0.1e1 / t7 / t78;
336   t85 = t45 * t12;
337   t86 = 0.1e1 / t85;
338   t91 = t9 * t9;
339   t92 = params->b * t91;
340   t93 = t11 * t11;
341   t94 = t45 * t13;
342   t96 = 0.1e1 / t19 / t94;
343   t98 = t24 * t24;
344   t99 = 0.1e1 / t98;
345   t103 = 0.304e3 / 0.9e1 * t10 * t11 * t80 * t25 - 0.688e3 / 0.9e1 * t43 * t44 * t86 * t50 + 0.128e3 / 0.3e1 * t92 * t93 * t96 * t99;
346   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
347     v2rho2[0] = -t4 * t71 * t28 / 0.6e1 - t4 * t8 * t54 - 0.3e1 / 0.8e1 * t34 * t35 * t103;
348 
349   t107 = t7 * t1;
350   t118 = t45 * t77;
351   t120 = 0.1e1 / t19 / t118;
352   t126 = t35 * (-0.32e2 / 0.3e1 * t10 * sigma[0] * t38 * t25 + 0.80e2 / 0.3e1 * t43 * t11 * t47 * t50 - 0.16e2 * t92 * t44 * t120 * t99);
353   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
354     v2rhosigma[0] = -t107 * t67 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t126;
355 
356   t137 = 0.1e1 / t19 / t85;
357   t143 = t35 * (0.6e1 * t92 * t11 * t137 * t99 - 0.8e1 * t43 * sigma[0] * t61 * t50 + 0.2e1 * t10 * t16 * t25);
358   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
359     v2sigma2[0] = -0.3e1 / 0.8e1 * t34 * t143;
360 
361 #ifndef XC_DONT_COMPILE_KXC
362 
363   if(order < 3) return;
364 
365 
366   t147 = 0.1e1 / t19 / rho[0];
367   t148 = t6 * t147;
368   t159 = 0.1e1 / t7 / t45;
369   t164 = 0.1e1 / t118;
370   t169 = t45 * t14;
371   t171 = 0.1e1 / t19 / t169;
372   t177 = params->b * t91 * params->gamma;
373   t178 = t93 * sigma[0];
374   t179 = t45 * t45;
375   t181 = 0.1e1 / t7 / t179;
376   t184 = 0.1e1 / t98 / t23;
377   t188 = -0.6688e4 / 0.27e2 * t10 * t11 * t159 * t25 + 0.25504e5 / 0.27e2 * t43 * t44 * t164 * t50 - 0.1152e4 * t92 * t93 * t171 * t99 + 0.4096e4 / 0.9e1 * t177 * t178 * t181 * t184;
378   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
379     v3rho3[0] = t4 * t148 * t28 / 0.9e1 - t4 * t71 * t54 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t103 - 0.3e1 / 0.8e1 * t34 * t35 * t188;
380 
381   t192 = t70 * t1;
382   t210 = 0.1e1 / t7 / t45 / t78;
383   t216 = t35 * (0.608e3 / 0.9e1 * t10 * sigma[0] * t80 * t25 - 0.2672e4 / 0.9e1 * t43 * t11 * t86 * t50 + 0.400e3 * t92 * t44 * t96 * t99 - 0.512e3 / 0.3e1 * t177 * t93 * t210 * t184);
384   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
385     v3rho2sigma[0] = -t192 * t67 / 0.6e1 - t107 * t126 - 0.3e1 / 0.8e1 * t34 * t216;
386 
387   t224 = t47 * t50;
388   t232 = t45 * t36;
389   t234 = 0.1e1 / t7 / t232;
390   t240 = t35 * (-0.32e2 / 0.3e1 * t10 * t38 * t25 + 0.224e3 / 0.3e1 * t43 * t224 * sigma[0] - 0.128e3 * t92 * t11 * t120 * t99 + 0.64e2 * t177 * t44 * t234 * t184);
391   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
392     v3rhosigma2[0] = -t107 * t143 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t240;
393 
394   t251 = 0.1e1 / t7 / t169;
395   t257 = t35 * (-0.24e2 * t177 * t11 * t251 * t184 + 0.36e2 * t92 * sigma[0] * t137 * t99 - 0.12e2 * t43 * t61 * t50);
396   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
397     v3sigma3[0] = -0.3e1 / 0.8e1 * t34 * t257;
398 
399 #ifndef XC_DONT_COMPILE_LXC
400 
401   if(order < 4) return;
402 
403 
404   t289 = t179 * rho[0];
405   t297 = params->b * t91 * t9;
406   t303 = 0.1e1 / t98 / t24;
407   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
408     v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t21 * t28 + 0.4e1 / 0.9e1 * t4 * t148 * t54 - t4 * t71 * t103 - 0.2e1 * t4 * t8 * t188 - 0.3e1 / 0.8e1 * t34 * t35 * (0.167200e6 / 0.81e2 * t10 * t11 / t7 / t46 * t25 - 0.948640e6 / 0.81e2 * t43 * t44 / t94 * t50 + 0.629120e6 / 0.27e2 * t92 * t93 / t19 / t232 * t99 - 0.532480e6 / 0.27e2 * t177 * t178 / t7 / t289 * t184 + 0.163840e6 / 0.27e2 * t297 * t93 * t11 / t179 / t13 * t303);
409 
410   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
411     v4rho3sigma[0] = t147 * t1 * t67 / 0.9e1 - t192 * t126 / 0.2e1 - 0.3e1 / 0.2e1 * t107 * t216 - 0.3e1 / 0.8e1 * t34 * t35 * (-0.13376e5 / 0.27e2 * t10 * sigma[0] * t159 * t25 + 0.89888e5 / 0.27e2 * t43 * t11 * t164 * t50 - 0.66976e5 / 0.9e1 * t92 * t44 * t171 * t99 + 0.61952e5 / 0.9e1 * t177 * t93 * t181 * t184 - 0.20480e5 / 0.9e1 * t297 * t178 / t179 / t77 * t303);
412 
413   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
414     v4rho2sigma2[0] = -t192 * t143 / 0.6e1 - t107 * t240 - 0.3e1 / 0.8e1 * t34 * t35 * (0.608e3 / 0.9e1 * t10 * t80 * t25 - 0.6560e4 / 0.9e1 * t43 * t86 * t50 * sigma[0] + 0.6272e4 / 0.3e1 * t92 * t96 * t99 * t11 - 0.6848e4 / 0.3e1 * t177 * t44 * t210 * t184 + 0.2560e4 / 0.3e1 * t297 * t93 / t179 / t12 * t303);
415 
416   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
417     v4rhosigma3[0] = -t107 * t257 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t35 * (0.96e2 * t43 * t224 - 0.480e3 * t92 * t120 * t99 * sigma[0] + 0.704e3 * t177 * t11 * t234 * t184 - 0.320e3 * t297 * t44 / t289 * t303);
418 
419   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
420     v4sigma4[0] = -0.3e1 / 0.8e1 * t34 * t35 * (0.72e2 * t92 * t137 * t99 - 0.192e3 * t177 * sigma[0] * t251 * t184 + 0.120e3 * t297 * t11 / t179 * t303);
421 
422 #ifndef XC_DONT_COMPILE_MXC
423 
424   if(order < 5) return;
425 
426 
427 #endif
428 
429 #endif
430 
431 #endif
432 
433 #endif
434 
435 #endif
436 
437 
438 }
439 
440 
441 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 *))442 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
443 {
444 
445 #ifndef XC_DONT_COMPILE_EXC
446   double t1, t3, t4, t5, t6, t7, t8, t9;
447   double t10, t11, t13, t14, t15, t16, t17, t18;
448   double t19, t20, t21, t22, t23, t24, t26, t29;
449   double t33, t34, t35, t38, t40, t42, t43, t44;
450   double t45, t46, t47, t48, t49, t50, t52, t55;
451   double t59, t60, t61, t64, t66;
452 
453 #ifndef XC_DONT_COMPILE_VXC
454   double t68, t69, t70, t71, t72, t73, t75, t76;
455   double t78, t79, t80, t81, t82, t84, t85, t86;
456   double t88, t93, t94, t95, t96, t97, t100, t104;
457   double t106, t107, t108, t109, t110, t112, t113, t114;
458   double t116, t117, t121, t122, t124, t125, t126, t127;
459   double t129, t130, t131, t133, t137, t138, t139, t140;
460   double t143, t147, t149, t150, t154, t155, t156, t160;
461   double t165, t169, t173, t178;
462 
463 #ifndef XC_DONT_COMPILE_FXC
464   double t182, t184, t185, t186, t187, t188, t189, t191;
465   double t192, t193, t195, t196, t197, t199, t200, t201;
466   double t202, t203, t204, t206, t207, t209, t210, t212;
467   double t213, t215, t216, t217, t218, t220, t225, t226;
468   double t231, t232, t233, t234, t236, t238, t239, t243;
469   double t245, t246, t247, t248, t249, t250, t251, t253;
470   double t254, t255, t257, t258, t259, t260, t262, t263;
471   double t264, t266, t267, t268, t270, t271, t272, t273;
472   double t275, t278, t280, t282, t285, t286, t287, t289;
473   double t293, t294, t295, t296, t298, t301, t303, t304;
474   double t306, t309, t311, t314, t316, t319, t320, t322;
475   double t323, t324, t325, t326, t328, t329, t330, t331;
476   double t333, t334, t335, t336, t338, t339, t340, t342;
477   double t343, t344, t345, t346, t348, t353, t354, t359;
478   double t360, t362, t364, t365, t369, t371, t372, t373;
479   double t376, t378, t379, t391, t393, t398, t403, t404;
480   double t405, t410, t415, t426, t428, t433, t434, t445;
481   double t450, t451, t462, t467, t468;
482 
483 #ifndef XC_DONT_COMPILE_KXC
484   double t475, t477, t481, t482, t483, t484, t486, t487;
485   double t488, t490, t492, t494, t495, t496, t498, t499;
486   double t501, t502, t503, t504, t506, t507, t508, t510;
487   double t512, t514, t515, t516, t518, t519, t520, t522;
488   double t523, t525, t528, t530, t532, t533, t534, t536;
489   double t537, t538, t539, t540, t542, t543, t545, t546;
490   double t547, t549, t550, t552, t554, t555, t557, t558;
491   double t559, t561, t562, t564, t565, t567, t570, t575;
492   double t580, t582, t588, t589, t590, t592, t595, t599;
493   double t601, t602, t603, t608, t609, t612, t613, t614;
494   double t616, t618, t620, t621, t622, t623, t625, t627;
495   double t628, t629, t631, t633, t635, t636, t637, t639;
496   double t641, t643, t645, t647, t649, t650, t652, t661;
497   double t663, t666, t667, t669, t671, t673, t675, t677;
498   double t679, t682, t683, t685, t686, t687, t689, t691;
499   double t693, t695, t697, t702, t705, t708, t710, t713;
500   double t718, t720, t722, t724, t726, t728, t730, t732;
501   double t734, t736, t738, t740, t742, t744, t751, t754;
502   double t757, t760, t763, t766, t768, t770, t772, t774;
503   double t777, t778, t780, t782, t784, t785, t786, t788;
504   double t790, t792, t797, t799, t802, t810, t811, t816;
505   double t817, t819, t820, t821, t823, t824, t825, t826;
506   double t828, t829, t831, t833, t834, t835, t838, t839;
507   double t841, t842, t843, t845, t846, t847, t849, t850;
508   double t851, t853, t854, t855, t856, t858, t859, t861;
509   double t866, t871, t873, t878, t879, t881, t884, t888;
510   double t890, t891, t892, t896, t900, t902, t903, t907;
511   double t910, t928, t933, t938, t942, t943, t947, t953;
512   double t954, t955, t957, t958, t959, t960, t964, t967;
513   double t968, t970, t974, t976, t977, t981, t984, t987;
514   double t992, t996, t1001, t1005, t1008, t1025, t1030, t1031;
515   double t1034, t1037, t1043, t1051, t1053, t1058, t1059, t1062;
516   double t1064, t1068, t1073, t1079, t1087, t1089, t1094, t1095;
517   double t1106, t1111, t1112, t1123, t1128, t1129;
518 
519 #ifndef XC_DONT_COMPILE_LXC
520   double t1143, t1147, t1150, t1153, t1159, t1166, t1169, t1172;
521   double t1175, t1181, t1184, t1187, t1201, t1203, t1205, t1211;
522   double t1218, t1220, t1221, t1222, t1227, t1229, t1231, t1232;
523   double t1234, t1236, t1242, t1244, t1246, t1247, t1248, t1253;
524   double t1255, t1257, t1259, t1266, t1277, t1287, t1290, t1293;
525   double t1312, t1320, t1326, t1337, t1341, t1342, t1345, t1347;
526   double t1353, t1356, t1358, t1363, t1370, t1373, t1374, t1376;
527   double t1377, t1379, t1386, t1389, t1392, t1395, t1401, t1402;
528   double t1411, t1415, t1429, t1449, t1450, t1451, t1453, t1465;
529   double t1466, t1472, t1474, t1478, t1484, t1486, t1492, t1494;
530   double t1497, t1503, t1510, t1516, t1531, t1536, t1538, t1541;
531   double t1545, t1550, t1558, t1583, t1588, t1605, t1612, t1622;
532   double t1625, t1631, t1633, t1640, t1645, t1651, t1654, t1657;
533   double t1663, t1666, t1673, t1679, t1681, t1687, t1690, t1692;
534   double t1695, t1703, t1706, t1735, t1760, t1776, t1790, t1803;
535   double t1810, t1864, t1870, t1872, t1882, t1887, t1890, t1893;
536   double t1896, t1899, t1902, t1904, t1907, t1919, t1950, t1967;
537   double t1975, t1977, t1987, t1990, t2003, t2049, t2061, t2072;
538   double t2082, t2107, t2108, t2113, t2118, t2137, t2140, t2142;
539   double t2145, t2174, t2177, t2181, t2199, t2202, t2212, t2215;
540   double t2218, t2222, t2239, t2250, t2251, t2253, t2255, t2257;
541   double t2258, t2270, t2274, t2281, t2284, t2285, t2287, t2300;
542   double t2303, t2316, t2328, t2345, t2353, t2362, t2371, t2447;
543   double t2448, t2450, t2488, t2489, t2501, t2516, t2526, t2578;
544   double t2603;
545 #endif
546 
547 #endif
548 
549 #endif
550 
551 #endif
552 
553 #endif
554 
555 
556   gga_x_optx_params *params;
557 
558   assert(p->params != NULL);
559   params = (gga_x_optx_params * )(p->params);
560 
561   t1 = M_CBRT3;
562   t3 = POW_1_3(0.1e1 / M_PI);
563   t4 = t1 * t3;
564   t5 = M_CBRT4;
565   t6 = t5 * t5;
566   t7 = t4 * t6;
567   t8 = rho[0] - rho[1];
568   t9 = rho[0] + rho[1];
569   t10 = 0.1e1 / t9;
570   t11 = t8 * t10;
571   t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
572   t14 = POW_1_3(t13);
573   t15 = t14 * t13;
574   t16 = POW_1_3(t9);
575   t17 = t15 * t16;
576   t18 = params->gamma * params->gamma;
577   t19 = params->b * t18;
578   t20 = sigma[0] * sigma[0];
579   t21 = rho[0] * rho[0];
580   t22 = t21 * t21;
581   t23 = t22 * rho[0];
582   t24 = POW_1_3(rho[0]);
583   t26 = 0.1e1 / t24 / t23;
584   t29 = t24 * t24;
585   t33 = 0.1e1 + params->gamma * sigma[0] / t29 / t21;
586   t34 = t33 * t33;
587   t35 = 0.1e1 / t34;
588   t38 = t19 * t20 * t26 * t35 + params->a;
589   t40 = t7 * t17 * t38;
590   t42 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
591   t43 = POW_1_3(t42);
592   t44 = t43 * t42;
593   t45 = t44 * t16;
594   t46 = sigma[2] * sigma[2];
595   t47 = rho[1] * rho[1];
596   t48 = t47 * t47;
597   t49 = t48 * rho[1];
598   t50 = POW_1_3(rho[1]);
599   t52 = 0.1e1 / t50 / t49;
600   t55 = t50 * t50;
601   t59 = 0.1e1 + params->gamma * sigma[2] / t55 / t47;
602   t60 = t59 * t59;
603   t61 = 0.1e1 / t60;
604   t64 = t19 * t46 * t52 * t61 + params->a;
605   t66 = t7 * t45 * t64;
606   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
607     zk[0] = -0.3e1 / 0.8e1 * t40 - 0.3e1 / 0.8e1 * t66;
608 
609 #ifndef XC_DONT_COMPILE_VXC
610 
611   if(order < 1) return;
612 
613 
614   t68 = 0.3e1 / 0.8e1 * t40;
615   t69 = 0.3e1 / 0.8e1 * t66;
616   t70 = t14 * t16;
617   t71 = t9 * t9;
618   t72 = 0.1e1 / t71;
619   t73 = t8 * t72;
620   t75 = t10 / 0.2e1 - t73 / 0.2e1;
621   t76 = t38 * t75;
622   t78 = t7 * t70 * t76;
623   t79 = t78 / 0.2e1;
624   t80 = t16 * t16;
625   t81 = 0.1e1 / t80;
626   t82 = t15 * t81;
627   t84 = t7 * t82 * t38;
628   t85 = t84 / 0.8e1;
629   t86 = t22 * t21;
630   t88 = 0.1e1 / t24 / t86;
631   t93 = params->b * t18 * params->gamma;
632   t94 = t20 * sigma[0];
633   t95 = t22 * t22;
634   t96 = t95 * rho[0];
635   t97 = 0.1e1 / t96;
636   t100 = 0.1e1 / t34 / t33;
637   t104 = 0.16e2 / 0.3e1 * t93 * t94 * t97 * t100 - 0.16e2 / 0.3e1 * t19 * t20 * t88 * t35;
638   t106 = t7 * t17 * t104;
639   t107 = 0.3e1 / 0.8e1 * t106;
640   t108 = t43 * t16;
641   t109 = -t75;
642   t110 = t64 * t109;
643   t112 = t7 * t108 * t110;
644   t113 = t112 / 0.2e1;
645   t114 = t44 * t81;
646   t116 = t7 * t114 * t64;
647   t117 = t116 / 0.8e1;
648   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
649     vrho[0] = -t68 - t69 + t9 * (-t79 - t85 - t107 - t113 - t117);
650 
651   t121 = -t10 / 0.2e1 - t73 / 0.2e1;
652   t122 = t38 * t121;
653   t124 = t7 * t70 * t122;
654   t125 = t124 / 0.2e1;
655   t126 = -t121;
656   t127 = t64 * t126;
657   t129 = t7 * t108 * t127;
658   t130 = t129 / 0.2e1;
659   t131 = t48 * t47;
660   t133 = 0.1e1 / t50 / t131;
661   t137 = t46 * sigma[2];
662   t138 = t48 * t48;
663   t139 = t138 * rho[1];
664   t140 = 0.1e1 / t139;
665   t143 = 0.1e1 / t60 / t59;
666   t147 = -0.16e2 / 0.3e1 * t19 * t46 * t133 * t61 + 0.16e2 / 0.3e1 * t93 * t137 * t140 * t143;
667   t149 = t7 * t45 * t147;
668   t150 = 0.3e1 / 0.8e1 * t149;
669   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
670     vrho[1] = -t68 - t69 + t9 * (-t125 - t85 - t130 - t117 - t150);
671 
672   t154 = t16 * t9 * t1;
673   t155 = t154 * t3;
674   t156 = t6 * t15;
675   t160 = 0.1e1 / t95;
676   t165 = -0.2e1 * t93 * t20 * t160 * t100 + 0.2e1 * t19 * sigma[0] * t26 * t35;
677   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
678     vsigma[0] = -0.3e1 / 0.8e1 * t155 * t156 * t165;
679 
680   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
681     vsigma[1] = 0.0e0;
682 
683   t169 = t6 * t44;
684   t173 = 0.1e1 / t138;
685   t178 = -0.2e1 * t93 * t46 * t173 * t143 + 0.2e1 * t19 * sigma[2] * t52 * t61;
686   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
687     vsigma[2] = -0.3e1 / 0.8e1 * t155 * t169 * t178;
688 
689 #ifndef XC_DONT_COMPILE_FXC
690 
691   if(order < 2) return;
692 
693 
694   t182 = t84 / 0.4e1;
695   t184 = t116 / 0.4e1;
696   t185 = t14 * t14;
697   t186 = 0.1e1 / t185;
698   t187 = t186 * t16;
699   t188 = t75 * t75;
700   t189 = t38 * t188;
701   t191 = t7 * t187 * t189;
702   t192 = t191 / 0.6e1;
703   t193 = t14 * t81;
704   t195 = t7 * t193 * t76;
705   t196 = t195 / 0.3e1;
706   t197 = t104 * t75;
707   t199 = t7 * t70 * t197;
708   t200 = t71 * t9;
709   t201 = 0.1e1 / t200;
710   t202 = t8 * t201;
711   t203 = -t72 + t202;
712   t204 = t38 * t203;
713   t206 = t7 * t70 * t204;
714   t207 = t206 / 0.2e1;
715   t209 = 0.1e1 / t80 / t9;
716   t210 = t15 * t209;
717   t212 = t7 * t210 * t38;
718   t213 = t212 / 0.12e2;
719   t215 = t7 * t82 * t104;
720   t216 = t215 / 0.4e1;
721   t217 = t21 * rho[0];
722   t218 = t22 * t217;
723   t220 = 0.1e1 / t24 / t218;
724   t225 = t95 * t21;
725   t226 = 0.1e1 / t225;
726   t231 = t18 * t18;
727   t232 = params->b * t231;
728   t233 = t20 * t20;
729   t234 = t95 * t22;
730   t236 = 0.1e1 / t29 / t234;
731   t238 = t34 * t34;
732   t239 = 0.1e1 / t238;
733   t243 = 0.304e3 / 0.9e1 * t19 * t20 * t220 * t35 - 0.688e3 / 0.9e1 * t93 * t94 * t226 * t100 + 0.128e3 / 0.3e1 * t232 * t233 * t236 * t239;
734   t245 = t7 * t17 * t243;
735   t246 = 0.3e1 / 0.8e1 * t245;
736   t247 = t43 * t43;
737   t248 = 0.1e1 / t247;
738   t249 = t248 * t16;
739   t250 = t109 * t109;
740   t251 = t64 * t250;
741   t253 = t7 * t249 * t251;
742   t254 = t253 / 0.6e1;
743   t255 = t43 * t81;
744   t257 = t7 * t255 * t110;
745   t258 = t257 / 0.3e1;
746   t259 = -t203;
747   t260 = t64 * t259;
748   t262 = t7 * t108 * t260;
749   t263 = t262 / 0.2e1;
750   t264 = t44 * t209;
751   t266 = t7 * t264 * t64;
752   t267 = t266 / 0.12e2;
753   t268 = -t192 - t196 - t199 - t207 + t213 - t216 - t246 - t254 - t258 - t263 + t267;
754   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
755     v2rho2[0] = -t78 - t182 - 0.3e1 / 0.4e1 * t106 - t112 - t184 + t9 * t268;
756 
757   t270 = t6 * t186;
758   t271 = t4 * t270;
759   t272 = t16 * t38;
760   t273 = t121 * t75;
761   t275 = t271 * t272 * t273;
762   t278 = t7 * t193 * t122;
763   t280 = t104 * t121;
764   t282 = t7 * t70 * t280;
765   t285 = 0.1e1 / t80 / t71;
766   t286 = t14 * t285;
767   t287 = t38 * t8;
768   t289 = t7 * t286 * t287;
769   t293 = t6 * t248;
770   t294 = t4 * t293;
771   t295 = t16 * t64;
772   t296 = t126 * t109;
773   t298 = t294 * t295 * t296;
774   t301 = t7 * t255 * t127;
775   t303 = t43 * t285;
776   t304 = t64 * t8;
777   t306 = t7 * t303 * t304;
778   t309 = t147 * t109;
779   t311 = t7 * t108 * t309;
780   t314 = t7 * t114 * t147;
781   t316 = -t275 / 0.6e1 - t278 / 0.6e1 - t282 / 0.2e1 - t289 / 0.2e1 - t195 / 0.6e1 + t213 - t215 / 0.8e1 - t298 / 0.6e1 - t301 / 0.6e1 + t306 / 0.2e1 - t257 / 0.6e1 + t267 - t311 / 0.2e1 - t314 / 0.8e1;
782   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
783     v2rho2[1] = t9 * t316 - t107 - t113 - t125 - t130 - t150 - t182 - t184 - t79;
784 
785   t319 = t121 * t121;
786   t320 = t38 * t319;
787   t322 = t7 * t187 * t320;
788   t323 = t322 / 0.6e1;
789   t324 = t278 / 0.3e1;
790   t325 = t72 + t202;
791   t326 = t38 * t325;
792   t328 = t7 * t70 * t326;
793   t329 = t328 / 0.2e1;
794   t330 = t126 * t126;
795   t331 = t64 * t330;
796   t333 = t7 * t249 * t331;
797   t334 = t333 / 0.6e1;
798   t335 = t301 / 0.3e1;
799   t336 = t147 * t126;
800   t338 = t7 * t108 * t336;
801   t339 = -t325;
802   t340 = t64 * t339;
803   t342 = t7 * t108 * t340;
804   t343 = t342 / 0.2e1;
805   t344 = t314 / 0.4e1;
806   t345 = t47 * rho[1];
807   t346 = t48 * t345;
808   t348 = 0.1e1 / t50 / t346;
809   t353 = t138 * t47;
810   t354 = 0.1e1 / t353;
811   t359 = t46 * t46;
812   t360 = t138 * t48;
813   t362 = 0.1e1 / t55 / t360;
814   t364 = t60 * t60;
815   t365 = 0.1e1 / t364;
816   t369 = 0.304e3 / 0.9e1 * t19 * t46 * t348 * t61 - 0.688e3 / 0.9e1 * t93 * t137 * t354 * t143 + 0.128e3 / 0.3e1 * t232 * t359 * t362 * t365;
817   t371 = t7 * t45 * t369;
818   t372 = 0.3e1 / 0.8e1 * t371;
819   t373 = -t323 - t324 - t329 + t213 - t334 - t335 - t338 - t343 + t267 - t344 - t372;
820   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
821     v2rho2[2] = -t124 - t182 - t129 - t184 - 0.3e1 / 0.4e1 * t149 + t9 * t373;
822 
823   t376 = t7 * t17 * t165;
824   t378 = t6 * t14;
825   t379 = t165 * t75;
826   t391 = t95 * t217;
827   t393 = 0.1e1 / t29 / t391;
828   t398 = -0.32e2 / 0.3e1 * t19 * sigma[0] * t88 * t35 + 0.80e2 / 0.3e1 * t93 * t20 * t97 * t100 - 0.16e2 * t232 * t94 * t393 * t239;
829   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
830     v2rhosigma[0] = -t376 / 0.2e1 - t155 * t378 * t379 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * t398;
831 
832   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
833     v2rhosigma[1] = 0.0e0;
834 
835   t403 = t7 * t45 * t178;
836   t404 = t6 * t43;
837   t405 = t178 * t109;
838   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
839     v2rhosigma[2] = -t155 * t404 * t405 / 0.2e1 - t403 / 0.2e1;
840 
841   t410 = t378 * t165 * t121;
842   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
843     v2rhosigma[3] = -t155 * t410 / 0.2e1 - t376 / 0.2e1;
844 
845   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
846     v2rhosigma[4] = 0.0e0;
847 
848   t415 = t404 * t178 * t126;
849   t426 = t138 * t345;
850   t428 = 0.1e1 / t55 / t426;
851   t433 = -0.32e2 / 0.3e1 * t19 * sigma[2] * t133 * t61 + 0.80e2 / 0.3e1 * t93 * t46 * t140 * t143 - 0.16e2 * t232 * t137 * t428 * t365;
852   t434 = t169 * t433;
853   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
854     v2rhosigma[5] = -t403 / 0.2e1 - t155 * t415 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t434;
855 
856   t445 = 0.1e1 / t29 / t225;
857   t450 = -0.8e1 * t93 * sigma[0] * t160 * t100 + 0.6e1 * t232 * t20 * t445 * t239 + 0.2e1 * t19 * t26 * t35;
858   t451 = t156 * t450;
859   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
860     v2sigma2[0] = -0.3e1 / 0.8e1 * t155 * t451;
861 
862   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
863     v2sigma2[1] = 0.0e0;
864 
865   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
866     v2sigma2[2] = 0.0e0;
867 
868   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
869     v2sigma2[3] = 0.0e0;
870 
871   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
872     v2sigma2[4] = 0.0e0;
873 
874   t462 = 0.1e1 / t55 / t353;
875   t467 = -0.8e1 * t93 * sigma[2] * t173 * t143 + 0.6e1 * t232 * t46 * t462 * t365 + 0.2e1 * t19 * t52 * t61;
876   t468 = t169 * t467;
877   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
878     v2sigma2[5] = -0.3e1 / 0.8e1 * t155 * t468;
879 
880 #ifndef XC_DONT_COMPILE_KXC
881 
882   if(order < 3) return;
883 
884 
885   t475 = t212 / 0.4e1;
886   t477 = t266 / 0.4e1;
887   t481 = 0.1e1 / t185 / t13;
888   t482 = t481 * t16;
889   t483 = t188 * t75;
890   t484 = t38 * t483;
891   t486 = t7 * t482 * t484;
892   t487 = t486 / 0.9e1;
893   t488 = t186 * t81;
894   t490 = t7 * t488 * t189;
895   t492 = t75 * t203;
896   t494 = t271 * t272 * t492;
897   t495 = t494 / 0.2e1;
898   t496 = t14 * t209;
899   t498 = t7 * t496 * t76;
900   t499 = t498 / 0.3e1;
901   t501 = 0.1e1 / t247 / t42;
902   t502 = t501 * t16;
903   t503 = t250 * t109;
904   t504 = t64 * t503;
905   t506 = t7 * t502 * t504;
906   t507 = t506 / 0.9e1;
907   t508 = t248 * t81;
908   t510 = t7 * t508 * t251;
909   t512 = t109 * t259;
910   t514 = t294 * t295 * t512;
911   t515 = t514 / 0.2e1;
912   t516 = t43 * t209;
913   t518 = t7 * t516 * t110;
914   t519 = t518 / 0.3e1;
915   t520 = t104 * t188;
916   t522 = t7 * t187 * t520;
917   t523 = t522 / 0.2e1;
918   t525 = t7 * t193 * t197;
919   t528 = t7 * t193 * t204;
920   t530 = t243 * t75;
921   t532 = t7 * t70 * t530;
922   t533 = 0.3e1 / 0.2e1 * t532;
923   t534 = t104 * t203;
924   t536 = t7 * t70 * t534;
925   t537 = 0.3e1 / 0.2e1 * t536;
926   t538 = t71 * t71;
927   t539 = 0.1e1 / t538;
928   t540 = t8 * t539;
929   t542 = 0.3e1 * t201 - 0.3e1 * t540;
930   t543 = t38 * t542;
931   t545 = t7 * t70 * t543;
932   t546 = t545 / 0.2e1;
933   t547 = t15 * t285;
934   t549 = t7 * t547 * t38;
935   t550 = 0.5e1 / 0.36e2 * t549;
936   t552 = t7 * t255 * t260;
937   t554 = -t542;
938   t555 = t64 * t554;
939   t557 = t7 * t108 * t555;
940   t558 = t557 / 0.2e1;
941   t559 = t44 * t285;
942   t561 = t7 * t559 * t64;
943   t562 = 0.5e1 / 0.36e2 * t561;
944   t564 = t7 * t210 * t104;
945   t565 = t564 / 0.4e1;
946   t567 = t7 * t82 * t243;
947   t570 = 0.1e1 / t24 / t95;
948   t575 = 0.1e1 / t391;
949   t580 = t95 * t23;
950   t582 = 0.1e1 / t29 / t580;
951   t588 = params->b * t231 * params->gamma;
952   t589 = t233 * sigma[0];
953   t590 = t95 * t95;
954   t592 = 0.1e1 / t24 / t590;
955   t595 = 0.1e1 / t238 / t33;
956   t599 = -0.6688e4 / 0.27e2 * t19 * t20 * t570 * t35 + 0.25504e5 / 0.27e2 * t93 * t94 * t575 * t100 - 0.1152e4 * t232 * t233 * t582 * t239 + 0.4096e4 / 0.9e1 * t588 * t589 * t592 * t595;
957   t601 = t7 * t17 * t599;
958   t602 = 0.3e1 / 0.8e1 * t601;
959   t603 = -t528 / 0.2e1 - t533 - t537 - t546 - t550 - t552 / 0.2e1 - t558 - t562 + t565 - 0.3e1 / 0.8e1 * t567 - t602;
960   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
961     v3rho3[0] = -t195 - t191 / 0.2e1 - t253 / 0.2e1 - t257 - 0.3e1 * t199 - 0.3e1 / 0.2e1 * t206 + t475 - 0.3e1 / 0.2e1 * t262 + t477 - 0.3e1 / 0.4e1 * t215 - 0.9e1 / 0.8e1 * t245 + t9 * (t487 - t490 / 0.6e1 - t495 + t499 + t507 - t510 / 0.6e1 - t515 + t519 - t523 - t525 + t603);
962 
963   t608 = t275 / 0.3e1;
964   t609 = t298 / 0.3e1;
965   t612 = t6 * t481;
966   t613 = t4 * t612;
967   t614 = t121 * t188;
968   t616 = t613 * t272 * t614;
969   t618 = t81 * t38;
970   t620 = t271 * t618 * t273;
971   t621 = t620 / 0.9e1;
972   t622 = t285 * t38;
973   t623 = t8 * t75;
974   t625 = t271 * t622 * t623;
975   t627 = t6 * t501;
976   t628 = t4 * t627;
977   t629 = t126 * t250;
978   t631 = t628 * t295 * t629;
979   t633 = t81 * t64;
980   t635 = t294 * t633 * t296;
981   t636 = t635 / 0.9e1;
982   t637 = t16 * t104;
983   t639 = t271 * t637 * t273;
984   t641 = t121 * t203;
985   t643 = t271 * t272 * t641;
986   t645 = t126 * t259;
987   t647 = t294 * t295 * t645;
988   t649 = t285 * t64;
989   t650 = t8 * t109;
990   t652 = t294 * t649 * t650;
991   t661 = t616 / 0.9e1 - t621 - t625 / 0.3e1 + t631 / 0.9e1 - t636 - t639 / 0.3e1 - t643 / 0.6e1 - t647 / 0.6e1 + t652 / 0.3e1 - t490 / 0.18e2 + 0.2e1 / 0.9e1 * t498 - t510 / 0.18e2 + 0.2e1 / 0.9e1 * t518 - t525 / 0.3e1 - t528 / 0.6e1 - t552 / 0.6e1;
992   t663 = t7 * t496 * t122;
993   t666 = 0.1e1 / t80 / t200;
994   t667 = t14 * t666;
995   t669 = t7 * t667 * t287;
996   t671 = t104 * t8;
997   t673 = t7 * t286 * t671;
998   t675 = t7 * t516 * t127;
999   t677 = t147 * t250;
1000   t679 = t7 * t249 * t677;
1001   t682 = t7 * t255 * t309;
1002   t683 = t682 / 0.3e1;
1003   t685 = t7 * t193 * t280;
1004   t686 = t685 / 0.3e1;
1005   t687 = t243 * t121;
1006   t689 = t7 * t70 * t687;
1007   t691 = t43 * t666;
1008   t693 = t7 * t691 * t304;
1009   t695 = t147 * t259;
1010   t697 = t7 * t108 * t695;
1011   t702 = t7 * t286 * t38;
1012   t705 = t7 * t303 * t64;
1013   t708 = t7 * t264 * t147;
1014   t710 = t663 / 0.9e1 + 0.7e1 / 0.6e1 * t669 - t673 + t675 / 0.9e1 - t679 / 0.6e1 - t683 - t686 - t689 / 0.2e1 - 0.7e1 / 0.6e1 * t693 - t697 / 0.2e1 - t550 - t562 + t564 / 0.6e1 - t567 / 0.8e1 - t702 / 0.2e1 + t705 / 0.2e1 + t708 / 0.12e2;
1015   t713 = -t207 + t475 - t263 + t477 - t282 + t306 - t311 - t344 - t215 / 0.2e1 - t246 + t9 * (t661 + t710);
1016   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1017     v3rho3[1] = -0.2e1 / 0.3e1 * t195 - t192 - t254 - 0.2e1 / 0.3e1 * t257 - t608 - t324 - t289 - t609 - t335 - t199 + t713;
1018 
1019   t718 = t319 * t75;
1020   t720 = t613 * t272 * t718;
1021   t722 = t325 * t75;
1022   t724 = t271 * t272 * t722;
1023   t726 = t330 * t109;
1024   t728 = t628 * t295 * t726;
1025   t730 = t126 * t8;
1026   t732 = t294 * t649 * t730;
1027   t734 = t16 * t147;
1028   t736 = t294 * t734 * t296;
1029   t738 = t339 * t109;
1030   t740 = t294 * t295 * t738;
1031   t742 = t121 * t8;
1032   t744 = t271 * t622 * t742;
1033   t751 = t720 / 0.9e1 - t724 / 0.6e1 + t728 / 0.9e1 + t732 / 0.3e1 - t736 / 0.3e1 - t740 / 0.6e1 - t744 / 0.3e1 - t621 - t636 + t498 / 0.9e1 + t518 / 0.9e1 + 0.2e1 / 0.9e1 * t663 - t669 / 0.3e1 + 0.2e1 / 0.9e1 * t675 - t683 - t686;
1034   t754 = t7 * t488 * t320;
1035   t757 = t7 * t193 * t326;
1036   t760 = t7 * t508 * t331;
1037   t763 = t7 * t255 * t336;
1038   t766 = t7 * t255 * t340;
1039   t768 = t104 * t319;
1040   t770 = t7 * t187 * t768;
1041   t772 = t104 * t325;
1042   t774 = t7 * t70 * t772;
1043   t777 = -t201 - 0.3e1 * t540;
1044   t778 = t38 * t777;
1045   t780 = t7 * t70 * t778;
1046   t782 = t147 * t8;
1047   t784 = t7 * t303 * t782;
1048   t785 = -t777;
1049   t786 = t64 * t785;
1050   t788 = t7 * t108 * t786;
1051   t790 = t369 * t109;
1052   t792 = t7 * t108 * t790;
1053   t797 = t7 * t114 * t369;
1054   t799 = t693 / 0.3e1 - t754 / 0.18e2 - t757 / 0.6e1 - t760 / 0.18e2 - t763 / 0.3e1 - t766 / 0.6e1 - t770 / 0.6e1 - t774 / 0.2e1 - t780 / 0.2e1 + t784 - t788 / 0.2e1 - t792 / 0.2e1 - t550 - t562 + t564 / 0.12e2 + t708 / 0.6e1 - t797 / 0.8e1;
1055   t802 = t477 - t282 + t306 - t311 - t314 / 0.2e1 - t329 - t338 - t343 - t216 - t372 + t9 * (t751 + t799);
1056   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1057     v3rho3[2] = -t196 - t258 - t608 - 0.2e1 / 0.3e1 * t278 - t289 - t609 - 0.2e1 / 0.3e1 * t301 - t323 - t334 + t475 + t802;
1058 
1059   t810 = t663 / 0.3e1;
1060   t811 = t675 / 0.3e1;
1061   t816 = t319 * t121;
1062   t817 = t38 * t816;
1063   t819 = t7 * t482 * t817;
1064   t820 = t819 / 0.9e1;
1065   t821 = t121 * t325;
1066   t823 = t271 * t272 * t821;
1067   t824 = t823 / 0.2e1;
1068   t825 = t330 * t126;
1069   t826 = t64 * t825;
1070   t828 = t7 * t502 * t826;
1071   t829 = t828 / 0.9e1;
1072   t831 = t126 * t339;
1073   t833 = t294 * t295 * t831;
1074   t834 = t833 / 0.2e1;
1075   t835 = t708 / 0.4e1;
1076   t838 = -0.3e1 * t201 - 0.3e1 * t540;
1077   t839 = t38 * t838;
1078   t841 = t7 * t70 * t839;
1079   t842 = t841 / 0.2e1;
1080   t843 = t147 * t330;
1081   t845 = t7 * t249 * t843;
1082   t846 = t845 / 0.2e1;
1083   t847 = t369 * t126;
1084   t849 = t7 * t108 * t847;
1085   t850 = 0.3e1 / 0.2e1 * t849;
1086   t851 = t147 * t339;
1087   t853 = t7 * t108 * t851;
1088   t854 = 0.3e1 / 0.2e1 * t853;
1089   t855 = -t838;
1090   t856 = t64 * t855;
1091   t858 = t7 * t108 * t856;
1092   t859 = t858 / 0.2e1;
1093   t861 = 0.1e1 / t50 / t138;
1094   t866 = 0.1e1 / t426;
1095   t871 = t138 * t49;
1096   t873 = 0.1e1 / t55 / t871;
1097   t878 = t359 * sigma[2];
1098   t879 = t138 * t138;
1099   t881 = 0.1e1 / t50 / t879;
1100   t884 = 0.1e1 / t364 / t59;
1101   t888 = -0.6688e4 / 0.27e2 * t19 * t46 * t861 * t61 + 0.25504e5 / 0.27e2 * t93 * t137 * t866 * t143 - 0.1152e4 * t232 * t359 * t873 * t365 + 0.4096e4 / 0.9e1 * t588 * t878 * t881 * t884;
1102   t890 = t7 * t45 * t888;
1103   t891 = 0.3e1 / 0.8e1 * t890;
1104   t892 = -t834 - t550 - t562 + t835 - 0.3e1 / 0.8e1 * t797 - t842 - t846 - t850 - t854 - t859 - t891;
1105   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1106     v3rho3[3] = -t278 - t301 - t322 / 0.2e1 - t333 / 0.2e1 + t475 + t477 - 0.3e1 / 0.4e1 * t314 - 0.3e1 / 0.2e1 * t328 - 0.3e1 * t338 - 0.3e1 / 0.2e1 * t342 - 0.9e1 / 0.8e1 * t371 + t9 * (t810 + t811 - t754 / 0.6e1 - t757 / 0.2e1 - t760 / 0.6e1 - t763 - t766 / 0.2e1 + t820 - t824 + t829 + t892);
1107 
1108   t896 = t7 * t70 * t379;
1109   t900 = t7 * t82 * t165 / 0.6e1;
1110   t902 = t7 * t17 * t398;
1111   t903 = t165 * t188;
1112   t907 = t398 * t75;
1113   t910 = t165 * t203;
1114   t928 = 0.1e1 / t24 / t95 / t218;
1115   t933 = 0.608e3 / 0.9e1 * t19 * sigma[0] * t220 * t35 - 0.2672e4 / 0.9e1 * t93 * t20 * t226 * t100 + 0.400e3 * t232 * t94 * t236 * t239 - 0.512e3 / 0.3e1 * t588 * t233 * t928 * t595;
1116   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1117     v3rho2sigma[0] = -0.4e1 / 0.3e1 * t896 - t900 - t902 - t155 * t270 * t903 / 0.6e1 - t155 * t378 * t907 - t155 * t378 * t910 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * t933;
1118 
1119   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1120     v3rho2sigma[1] = 0.0e0;
1121 
1122   t938 = t7 * t108 * t405;
1123   t942 = t7 * t114 * t178 / 0.6e1;
1124   t943 = t178 * t250;
1125   t947 = t178 * t259;
1126   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1127     v3rho2sigma[2] = -0.4e1 / 0.3e1 * t938 - t942 - t155 * t293 * t943 / 0.6e1 - t155 * t404 * t947 / 0.2e1;
1128 
1129   t953 = t16 * t1;
1130   t954 = t953 * t3;
1131   t955 = t954 * t410;
1132   t957 = t3 * t6;
1133   t958 = t154 * t957;
1134   t959 = t186 * t165;
1135   t960 = t959 * t273;
1136   t964 = t378 * t398 * t121;
1137   t967 = t209 * t1;
1138   t968 = t967 * t3;
1139   t970 = t378 * t165 * t8;
1140   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1141     v3rho2sigma[3] = -0.2e1 / 0.3e1 * t896 - t900 - t902 / 0.2e1 - 0.2e1 / 0.3e1 * t955 - t958 * t960 / 0.6e1 - t155 * t964 / 0.2e1 - t968 * t970 / 0.2e1;
1142 
1143   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1144     v3rho2sigma[4] = 0.0e0;
1145 
1146   t974 = t954 * t415;
1147   t976 = t248 * t178;
1148   t977 = t976 * t296;
1149   t981 = t404 * t178 * t8;
1150   t984 = t954 * t434;
1151   t987 = t404 * t433 * t109;
1152   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1153     v3rho2sigma[5] = -0.2e1 / 0.3e1 * t938 - t942 - 0.2e1 / 0.3e1 * t974 - t958 * t977 / 0.6e1 + t968 * t981 / 0.2e1 - t984 / 0.2e1 - t155 * t987 / 0.2e1;
1154 
1155   t992 = t270 * t165 * t319;
1156   t996 = t378 * t165 * t325;
1157   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1158     v3rho2sigma[6] = -0.4e1 / 0.3e1 * t955 - t900 - t155 * t992 / 0.6e1 - t155 * t996 / 0.2e1;
1159 
1160   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1161     v3rho2sigma[7] = 0.0e0;
1162 
1163   t1001 = t293 * t178 * t330;
1164   t1005 = t404 * t433 * t126;
1165   t1008 = t404 * t178 * t339;
1166   t1025 = 0.1e1 / t50 / t138 / t346;
1167   t1030 = 0.608e3 / 0.9e1 * t19 * sigma[2] * t348 * t61 - 0.2672e4 / 0.9e1 * t93 * t46 * t354 * t143 + 0.400e3 * t232 * t137 * t362 * t365 - 0.512e3 / 0.3e1 * t588 * t359 * t1025 * t884;
1168   t1031 = t169 * t1030;
1169   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1170     v3rho2sigma[8] = -0.4e1 / 0.3e1 * t974 - t942 - t984 - t155 * t1001 / 0.6e1 - t155 * t1005 - t155 * t1008 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1031;
1171 
1172   t1034 = t954 * t451;
1173   t1037 = t378 * t450 * t75;
1174   t1043 = t97 * t100;
1175   t1051 = t95 * t86;
1176   t1053 = 0.1e1 / t24 / t1051;
1177   t1058 = -0.32e2 / 0.3e1 * t19 * t88 * t35 + 0.224e3 / 0.3e1 * t93 * t1043 * sigma[0] - 0.128e3 * t232 * t20 * t393 * t239 + 0.64e2 * t588 * t94 * t1053 * t595;
1178   t1059 = t156 * t1058;
1179   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1180     v3rhosigma2[0] = -t1034 / 0.2e1 - t155 * t1037 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1059;
1181 
1182   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1183     v3rhosigma2[1] = 0.0e0;
1184 
1185   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1186     v3rhosigma2[2] = 0.0e0;
1187 
1188   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1189     v3rhosigma2[3] = 0.0e0;
1190 
1191   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1192     v3rhosigma2[4] = 0.0e0;
1193 
1194   t1062 = t954 * t468;
1195   t1064 = t404 * t467 * t109;
1196   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1197     v3rhosigma2[5] = -t155 * t1064 / 0.2e1 - t1062 / 0.2e1;
1198 
1199   t1068 = t378 * t450 * t121;
1200   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1201     v3rhosigma2[6] = -t155 * t1068 / 0.2e1 - t1034 / 0.2e1;
1202 
1203   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1204     v3rhosigma2[7] = 0.0e0;
1205 
1206   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1207     v3rhosigma2[8] = 0.0e0;
1208 
1209   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1210     v3rhosigma2[9] = 0.0e0;
1211 
1212   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1213     v3rhosigma2[10] = 0.0e0;
1214 
1215   t1073 = t404 * t467 * t126;
1216   t1079 = t140 * t143;
1217   t1087 = t138 * t131;
1218   t1089 = 0.1e1 / t50 / t1087;
1219   t1094 = -0.32e2 / 0.3e1 * t19 * t133 * t61 + 0.224e3 / 0.3e1 * t93 * t1079 * sigma[2] - 0.128e3 * t232 * t46 * t428 * t365 + 0.64e2 * t588 * t137 * t1089 * t884;
1220   t1095 = t169 * t1094;
1221   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1222     v3rhosigma2[11] = -t1062 / 0.2e1 - t155 * t1073 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1095;
1223 
1224   t1106 = 0.1e1 / t24 / t580;
1225   t1111 = -0.24e2 * t588 * t20 * t1106 * t595 + 0.36e2 * t232 * sigma[0] * t445 * t239 - 0.12e2 * t93 * t160 * t100;
1226   t1112 = t156 * t1111;
1227   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1228     v3sigma3[0] = -0.3e1 / 0.8e1 * t155 * t1112;
1229 
1230   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1231     v3sigma3[1] = 0.0e0;
1232 
1233   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1234     v3sigma3[2] = 0.0e0;
1235 
1236   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1237     v3sigma3[3] = 0.0e0;
1238 
1239   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1240     v3sigma3[4] = 0.0e0;
1241 
1242   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1243     v3sigma3[5] = 0.0e0;
1244 
1245   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1246     v3sigma3[6] = 0.0e0;
1247 
1248   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1249     v3sigma3[7] = 0.0e0;
1250 
1251   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1252     v3sigma3[8] = 0.0e0;
1253 
1254   t1123 = 0.1e1 / t50 / t871;
1255   t1128 = -0.24e2 * t588 * t46 * t1123 * t884 + 0.36e2 * t232 * sigma[2] * t462 * t365 - 0.12e2 * t93 * t173 * t143;
1256   t1129 = t169 * t1128;
1257   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1258     v3sigma3[9] = -0.3e1 / 0.8e1 * t155 * t1129;
1259 
1260 #ifndef XC_DONT_COMPILE_LXC
1261 
1262   if(order < 4) return;
1263 
1264 
1265   t1143 = 0.4e1 / 0.9e1 * t486 - 0.2e1 / 0.3e1 * t490 - 0.2e1 * t494 + 0.4e1 / 0.3e1 * t498 + 0.4e1 / 0.9e1 * t506 - 0.2e1 / 0.3e1 * t510 - 0.2e1 * t514 + 0.4e1 / 0.3e1 * t518 - 0.2e1 * t522 - 0.4e1 * t525 - 0.2e1 * t528;
1266   t1147 = 0.5e1 / 0.9e1 * t549;
1267   t1150 = 0.5e1 / 0.9e1 * t561;
1268   t1153 = t259 * t259;
1269   t1159 = t7 * t303 * t110;
1270   t1166 = t7 * t488 * t520;
1271   t1169 = t7 * t496 * t197;
1272   t1172 = t7 * t496 * t204;
1273   t1175 = t7 * t516 * t260;
1274   t1181 = t7 * t193 * t530;
1275   t1184 = t7 * t193 * t534;
1276   t1187 = t7 * t193 * t543;
1277   t1201 = t538 * t9;
1278   t1203 = t8 / t1201;
1279   t1205 = -0.12e2 * t539 + 0.12e2 * t1203;
1280   t1211 = t7 * t255 * t555;
1281   t1218 = t13 * t13;
1282   t1220 = 0.1e1 / t185 / t1218;
1283   t1221 = t1220 * t16;
1284   t1222 = t188 * t188;
1285   t1227 = t481 * t81;
1286   t1229 = t7 * t1227 * t484;
1287   t1231 = -t7 * t249 * t64 * t1153 / 0.2e1 - 0.20e2 / 0.27e2 * t1159 + 0.4e1 / 0.9e1 * t7 * t482 * t104 * t483 - 0.2e1 / 0.3e1 * t1166 + 0.4e1 / 0.3e1 * t1169 + 0.2e1 / 0.3e1 * t1172 + 0.2e1 / 0.3e1 * t1175 - t7 * t187 * t243 * t188 - 0.2e1 * t1181 - 0.2e1 * t1184 - 0.2e1 / 0.3e1 * t1187 - 0.2e1 * t7 * t70 * t599 * t75 - 0.3e1 * t7 * t70 * t243 * t203 - 0.2e1 * t7 * t70 * t104 * t542 - t7 * t70 * t38 * t1205 / 0.2e1 - 0.2e1 / 0.3e1 * t1211 + t7 * t108 * t64 * t1205 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1221 * t38 * t1222 + 0.4e1 / 0.27e2 * t1229;
1288   t1232 = t186 * t209;
1289   t1234 = t7 * t1232 * t189;
1290   t1236 = t203 * t203;
1291   t1242 = t7 * t286 * t76;
1292   t1244 = t42 * t42;
1293   t1246 = 0.1e1 / t247 / t1244;
1294   t1247 = t1246 * t16;
1295   t1248 = t250 * t250;
1296   t1253 = t501 * t81;
1297   t1255 = t7 * t1253 * t504;
1298   t1257 = t248 * t209;
1299   t1259 = t7 * t1257 * t251;
1300   t1266 = t294 * t633 * t512;
1301   t1277 = t271 * t618 * t492;
1302   t1287 = t7 * t547 * t104;
1303   t1290 = t7 * t210 * t243;
1304   t1293 = t7 * t82 * t599;
1305   t1312 = t590 * rho[0];
1306   t1320 = params->b * t231 * t18;
1307   t1326 = 0.1e1 / t238 / t34;
1308   t1337 = 0.10e2 / 0.27e2 * t7 * t15 * t666 * t38;
1309   t1341 = 0.10e2 / 0.27e2 * t7 * t44 * t666 * t64;
1310   t1342 = 0.2e1 / 0.9e1 * t1234 - t7 * t187 * t38 * t1236 / 0.2e1 - 0.20e2 / 0.27e2 * t1242 - 0.5e1 / 0.27e2 * t7 * t1247 * t64 * t1248 + 0.4e1 / 0.27e2 * t1255 + 0.2e1 / 0.9e1 * t1259 + 0.2e1 / 0.3e1 * t628 * t295 * t250 * t259 - 0.2e1 / 0.3e1 * t1266 - 0.2e1 / 0.3e1 * t294 * t295 * t109 * t554 + 0.2e1 / 0.3e1 * t613 * t272 * t188 * t203 - 0.2e1 / 0.3e1 * t1277 - 0.2e1 * t271 * t637 * t492 - 0.2e1 / 0.3e1 * t271 * t272 * t75 * t542 - 0.5e1 / 0.9e1 * t1287 + t1290 / 0.2e1 - t1293 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * (0.167200e6 / 0.81e2 * t19 * t20 / t24 / t96 * t35 - 0.948640e6 / 0.81e2 * t93 * t94 / t234 * t100 + 0.629120e6 / 0.27e2 * t232 * t233 / t29 / t1051 * t239 - 0.532480e6 / 0.27e2 * t588 * t589 / t24 / t1312 * t595 + 0.163840e6 / 0.27e2 * t1320 * t233 * t20 / t590 / t22 * t1326) + t1337 + t1341;
1311   t1345 = -0.6e1 * t532 - 0.6e1 * t536 - 0.2e1 * t545 - t1147 - 0.2e1 * t552 - 0.2e1 * t557 - t1150 + t564 - 0.3e1 / 0.2e1 * t567 - 0.3e1 / 0.2e1 * t601 + t9 * (t1231 + t1342);
1312   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1313     v4rho4[0] = t1143 + t1345;
1314 
1315   t1347 = t7 * t667 * t38;
1316   t1353 = t7 * t691 * t64;
1317   t1356 = t7 * t559 * t147;
1318   t1358 = t248 * t285;
1319   t1363 = t7 * t303 * t127;
1320   t1370 = t7 * t508 * t677;
1321   t1373 = t7 * t516 * t309;
1322   t1374 = t1373 / 0.3e1;
1323   t1376 = t7 * t496 * t280;
1324   t1377 = t1376 / 0.3e1;
1325   t1379 = t7 * t667 * t671;
1326   t1386 = 0.1e1 / t80 / t538;
1327   t1389 = t7 * t43 * t1386 * t304;
1328   t1392 = t7 * t255 * t695;
1329   t1395 = t7 * t193 * t687;
1330   t1401 = 0.5e1 / 0.2e1 * t1347 - 0.3e1 / 0.2e1 * t7 * t286 * t104 - 0.5e1 / 0.2e1 * t1353 - 0.5e1 / 0.36e2 * t1356 + t7 * t1358 * t110 / 0.2e1 - 0.5e1 / 0.27e2 * t1363 + t7 * t502 * t147 * t503 / 0.9e1 - t1370 / 0.6e1 + t1374 + t1377 + 0.7e1 / 0.2e1 * t1379 - 0.3e1 / 0.2e1 * t7 * t286 * t243 * t8 + 0.25e2 / 0.6e1 * t1389 - t1392 / 0.2e1 - t1395 / 0.2e1 - t7 * t70 * t599 * t121 / 0.2e1;
1331   t1402 = t186 * t285;
1332   t1411 = t7 * t286 * t122;
1333   t1415 = t7 * t14 * t1386 * t287;
1334   t1429 = -t7 * t1402 * t76 / 0.2e1 - t7 * t108 * t147 * t554 / 0.2e1 - 0.5e1 / 0.27e2 * t1411 - 0.25e2 / 0.6e1 * t1415 - 0.5e1 / 0.9e1 * t1159 - t1166 / 0.6e1 + 0.2e1 / 0.3e1 * t1169 + t1172 / 0.3e1 + t1175 / 0.3e1 - t1181 / 0.2e1 - t1184 / 0.2e1 - t1187 / 0.6e1 - t1211 / 0.6e1 + t1229 / 0.27e2 + t1234 / 0.9e1 - 0.5e1 / 0.9e1 * t1242;
1335   t1449 = t271 * t209 * t38 * t273;
1336   t1450 = t1449 / 0.9e1;
1337   t1451 = t666 * t38;
1338   t1453 = t271 * t1451 * t623;
1339   t1465 = t294 * t209 * t64 * t296;
1340   t1466 = t1465 / 0.9e1;
1341   t1472 = t271 * t81 * t104 * t273;
1342   t1474 = t285 * t104;
1343   t1478 = t271 * t618 * t641;
1344   t1484 = t1255 / 0.27e2 + t1259 / 0.9e1 - t294 * t734 * t512 / 0.2e1 - t271 * t16 * t243 * t273 / 0.2e1 - t271 * t637 * t641 / 0.2e1 - t271 * t272 * t121 * t542 / 0.6e1 + t1450 + 0.7e1 / 0.6e1 * t1453 - t294 * t295 * t126 * t554 / 0.6e1 + t294 * t649 * t8 * t259 / 0.2e1 + t1466 + t613 * t637 * t614 / 0.3e1 - t1472 / 0.3e1 - t271 * t1474 * t623 - t1478 / 0.6e1 - t271 * t622 * t8 * t203 / 0.2e1;
1345   t1486 = t294 * t633 * t645;
1346   t1492 = t666 * t64;
1347   t1494 = t294 * t1492 * t650;
1348   t1497 = t4 * t6 * t1220;
1349   t1503 = t613 * t618 * t614;
1350   t1510 = t4 * t6 * t1246;
1351   t1516 = t628 * t633 * t629;
1352   t1531 = -t1486 / 0.6e1 - t628 * t649 * t8 * t250 / 0.3e1 - 0.7e1 / 0.6e1 * t1494 - 0.5e1 / 0.27e2 * t1497 * t272 * t121 * t483 + t1503 / 0.9e1 + t613 * t622 * t8 * t188 / 0.3e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t126 * t503 + t1516 / 0.9e1 - t1266 / 0.6e1 - t1277 / 0.6e1 + t613 * t272 * t641 * t75 / 0.3e1 + t628 * t295 * t645 * t109 / 0.3e1 - 0.5e1 / 0.12e2 * t1287 + t1290 / 0.4e1 - t1293 / 0.8e1 + t1337 + t1341;
1353   t1536 = t620 / 0.3e1;
1354   t1538 = t635 / 0.3e1;
1355   t1541 = t9 * (t1401 + t1429 + t1484 + t1531) + t616 / 0.3e1 - t1536 - t625 + t631 / 0.3e1 - t1538 - t639 - t643 / 0.2e1 - t647 / 0.2e1 + t652 - t495;
1356   t1545 = -t515 + t487 - t490 / 0.3e1 + t498 + t507 - t510 / 0.3e1 + t518 - t523 - 0.2e1 * t525 - t528 - t533;
1357   t1550 = -t537 - t546 - t552 - t558 + t810 + 0.7e1 / 0.2e1 * t669 - 0.3e1 * t673 + t811 - t679 / 0.2e1 - t682 - t685;
1358   t1558 = -0.3e1 / 0.2e1 * t689 - 0.7e1 / 0.2e1 * t693 - 0.3e1 / 0.2e1 * t697 - t1147 - t1150 + 0.3e1 / 0.4e1 * t564 - 0.3e1 / 0.4e1 * t567 - t602 - 0.3e1 / 0.2e1 * t702 + 0.3e1 / 0.2e1 * t705 + t835;
1359   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1360     v4rho4[1] = t1541 + t1545 + t1550 + t1558;
1361 
1362   t1583 = t294 * t1492 * t730;
1363   t1588 = t285 * t147;
1364   t1605 = t294 * t81 * t147 * t296;
1365   t1612 = t294 * t633 * t738;
1366   t1622 = t613 * t618 * t718;
1367   t1625 = t271 * t1451 * t742;
1368   t1631 = 0.2e1 / 0.9e1 * t613 * t637 * t718 + t613 * t272 * t319 * t203 / 0.9e1 - t271 * t637 * t722 / 0.3e1 - t271 * t272 * t777 * t75 / 0.3e1 - t271 * t272 * t325 * t203 / 0.6e1 + t628 * t295 * t330 * t259 / 0.9e1 - 0.7e1 / 0.9e1 * t1583 - t294 * t734 * t645 / 0.3e1 + 0.2e1 / 0.3e1 * t294 * t1588 * t650 - t294 * t295 * t785 * t109 / 0.3e1 - t294 * t295 * t339 * t259 / 0.6e1 + 0.2e1 / 0.9e1 * t628 * t734 * t629 - 0.2e1 / 0.9e1 * t1605 + t628 * t295 * t339 * t250 / 0.9e1 - t1612 / 0.9e1 - 0.2e1 / 0.3e1 * t271 * t1474 * t742 - 0.5e1 / 0.27e2 * t1497 * t272 * t319 * t188 + 0.2e1 / 0.27e2 * t1622 + 0.7e1 / 0.9e1 * t1625 + t613 * t272 * t325 * t188 / 0.9e1;
1369   t1633 = t271 * t618 * t722;
1370   t1640 = t628 * t633 * t726;
1371   t1645 = t7 * t264 * t369;
1372   t1651 = t7 * t1232 * t320;
1373   t1654 = t7 * t496 * t326;
1374   t1657 = t7 * t1257 * t331;
1375   t1663 = t7 * t516 * t336;
1376   t1666 = t7 * t516 * t340;
1377   t1673 = t7 * t255 * t790;
1378   t1679 = 0.1e1 / t80 / t1201;
1379   t1681 = t8 * t8;
1380   t1687 = t7 * t488 * t768;
1381   t1690 = t7 * t193 * t772;
1382   t1692 = -t1633 / 0.9e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t330 * t250 + 0.2e1 / 0.27e2 * t1640 + t7 * t303 * t147 + t1645 / 0.12e2 - t1347 / 0.3e1 + t1353 / 0.3e1 - 0.5e1 / 0.18e2 * t1356 + t1651 / 0.27e2 + t1654 / 0.9e1 + t1657 / 0.27e2 + t7 * t1358 * t127 / 0.3e1 + 0.2e1 / 0.9e1 * t1663 + t1666 / 0.9e1 - t7 * t249 * t369 * t250 / 0.6e1 - t1673 / 0.3e1 - t7 * t1402 * t122 / 0.3e1 - t7 * t186 * t1679 * t38 * t1681 / 0.3e1 - t1687 / 0.9e1 - t1690 / 0.3e1;
1383   t1695 = t7 * t193 * t778;
1384   t1703 = t7 * t691 * t782;
1385   t1706 = t7 * t255 * t786;
1386   t1735 = -t1695 / 0.3e1 - t7 * t248 * t1679 * t64 * t1681 / 0.3e1 - 0.7e1 / 0.3e1 * t1703 - t1706 / 0.3e1 - t7 * t187 * t243 * t319 / 0.6e1 - t7 * t70 * t243 * t325 / 0.2e1 - t7 * t70 * t104 * t777 - t7 * t108 * t369 * t259 / 0.2e1 - 0.10e2 / 0.27e2 * t1363 - t1370 / 0.9e1 + 0.4e1 / 0.9e1 * t1373 + 0.4e1 / 0.9e1 * t1376 - 0.2e1 / 0.3e1 * t1379 + 0.41e2 / 0.9e1 * t1389 - t1392 / 0.3e1 - t1395 / 0.3e1 - 0.10e2 / 0.27e2 * t1411 - 0.41e2 / 0.9e1 * t1415 - 0.10e2 / 0.27e2 * t1159 + 0.2e1 / 0.9e1 * t1169;
1387   t1760 = t1172 / 0.9e1 + t1175 / 0.9e1 + t1234 / 0.27e2 - 0.10e2 / 0.27e2 * t1242 + t1259 / 0.27e2 + 0.4e1 / 0.27e2 * t1449 - 0.2e1 / 0.9e1 * t1453 + 0.4e1 / 0.27e2 * t1465 - 0.2e1 / 0.9e1 * t1472 - t1478 / 0.9e1 - t1486 / 0.9e1 + 0.2e1 / 0.9e1 * t1494 + 0.2e1 / 0.27e2 * t1503 + 0.2e1 / 0.27e2 * t1516 - 0.4e1 / 0.9e1 * t628 * t649 * t296 * t8 + 0.4e1 / 0.9e1 * t613 * t622 * t273 * t8 - 0.5e1 / 0.18e2 * t1287 + t1290 / 0.12e2 + t1337 + t1341;
1388   t1776 = t9 * (t1631 + t1692 + t1735 + t1760) + 0.2e1 / 0.9e1 * t720 - t724 / 0.3e1 + 0.2e1 / 0.9e1 * t728 + 0.2e1 / 0.3e1 * t732 - 0.2e1 / 0.3e1 * t736 - t740 / 0.3e1 - 0.2e1 / 0.3e1 * t744 + 0.2e1 / 0.9e1 * t616 - 0.4e1 / 0.9e1 * t620 - 0.2e1 / 0.3e1 * t625 + 0.2e1 / 0.9e1 * t631 - 0.4e1 / 0.9e1 * t635;
1389   t1790 = -0.2e1 / 0.3e1 * t639 - t643 / 0.3e1 - t647 / 0.3e1 + 0.2e1 / 0.3e1 * t652 - t490 / 0.9e1 + 0.2e1 / 0.3e1 * t498 - t510 / 0.9e1 + 0.2e1 / 0.3e1 * t518 - 0.2e1 / 0.3e1 * t525 - t528 / 0.3e1 - t552 / 0.3e1 + 0.2e1 / 0.3e1 * t663 + 0.5e1 / 0.3e1 * t669;
1390   t1803 = -0.2e1 * t673 + 0.2e1 / 0.3e1 * t675 - t679 / 0.3e1 - 0.4e1 / 0.3e1 * t682 - 0.4e1 / 0.3e1 * t685 - t689 - 0.5e1 / 0.3e1 * t693 - t697 - t754 / 0.9e1 - t757 / 0.3e1 - t760 / 0.9e1 - 0.2e1 / 0.3e1 * t763 - t766 / 0.3e1;
1391   t1810 = -t770 / 0.3e1 - t774 - t780 + 0.2e1 * t784 - t788 - t792 - t1147 - t1150 + t564 / 0.2e1 - t567 / 0.4e1 - t702 + t705 + t708 / 0.2e1 - t797 / 0.4e1;
1392   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1393     v4rho4[2] = t1776 + t1790 + t1803 + t1810;
1394 
1395   t1864 = t271 * t618 * t821;
1396   t1870 = -0.5e1 / 0.27e2 * t1497 * t272 * t816 * t75 - t271 * t637 * t821 / 0.2e1 - t271 * t272 * t121 * t777 / 0.2e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t825 * t109 - t628 * t649 * t330 * t8 / 0.3e1 - t294 * t295 * t126 * t785 / 0.2e1 + t294 * t649 * t8 * t339 / 0.2e1 - t271 * t272 * t838 * t75 / 0.6e1 + t628 * t734 * t726 / 0.3e1 + t294 * t1588 * t730 - t294 * t16 * t369 * t296 / 0.2e1 - t294 * t734 * t738 / 0.2e1 - t294 * t295 * t855 * t109 / 0.6e1 + t613 * t622 * t319 * t8 / 0.3e1 - t1864 / 0.6e1 - t271 * t622 * t8 * t325 / 0.2e1;
1397   t1872 = t294 * t633 * t831;
1398   t1882 = t7 * t114 * t888;
1399   t1887 = t7 * t1227 * t817;
1400   t1890 = t7 * t1253 * t826;
1401   t1893 = t7 * t193 * t839;
1402   t1896 = t7 * t508 * t843;
1403   t1899 = t7 * t255 * t847;
1404   t1902 = t7 * t255 * t851;
1405   t1904 = -t1872 / 0.6e1 + t1583 / 0.3e1 - t1605 / 0.3e1 - t1612 / 0.6e1 + t1622 / 0.9e1 - t1625 / 0.3e1 - t1633 / 0.6e1 + t1640 / 0.9e1 - t1882 / 0.8e1 + t1645 / 0.4e1 - 0.5e1 / 0.12e2 * t1356 + t1887 / 0.27e2 + t1890 / 0.27e2 - t1893 / 0.6e1 - t1896 / 0.6e1 - t1899 / 0.2e1 - t1902 / 0.2e1;
1406   t1907 = t7 * t255 * t856;
1407   t1919 = 0.6e1 * t539 + 0.12e2 * t1203;
1408   t1950 = -t1907 / 0.6e1 + t7 * t482 * t104 * t816 / 0.9e1 - t7 * t70 * t104 * t838 / 0.2e1 - t7 * t70 * t38 * t1919 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t303 * t369 * t8 - 0.3e1 / 0.2e1 * t7 * t108 * t147 * t785 + t7 * t108 * t64 * t1919 / 0.2e1 - t7 * t108 * t888 * t109 / 0.2e1 + t1651 / 0.9e1 + t1654 / 0.3e1 + t1657 / 0.9e1 + 0.2e1 / 0.3e1 * t1663 + t1666 / 0.3e1 - t1673 / 0.2e1 - t1687 / 0.6e1 - t1690 / 0.2e1 - t1695 / 0.2e1;
1409   t1967 = t1703 - t1706 / 0.2e1 - 0.5e1 / 0.9e1 * t1363 + t1374 + t1377 - t1389 / 0.3e1 - 0.5e1 / 0.9e1 * t1411 + t1415 / 0.3e1 - 0.5e1 / 0.27e2 * t1159 - 0.5e1 / 0.27e2 * t1242 + t1450 + t1466 + t613 * t272 * t821 * t75 / 0.3e1 + t628 * t295 * t831 * t109 / 0.3e1 - 0.5e1 / 0.36e2 * t1287 + t1337 + t1341;
1410   t1975 = t9 * (t1870 + t1904 + t1950 + t1967) - t824 - t834 + t720 / 0.3e1 - t724 / 0.2e1 + t728 / 0.3e1 + t732 - t736 - t740 / 0.2e1 - t744 - t1536;
1411   t1977 = -t1538 + t499 + t519 + t663 - t669 + t675 - t682 - t685 + t693 - t754 / 0.3e1 - t757;
1412   t1987 = -t760 / 0.3e1 - 0.2e1 * t763 - t766 - t770 / 0.2e1 - 0.3e1 / 0.2e1 * t774 - 0.3e1 / 0.2e1 * t780 + 0.3e1 * t784 - 0.3e1 / 0.2e1 * t788 - 0.3e1 / 0.2e1 * t792 + t820 + t829;
1413   t1990 = -t842 - t846 - t850 - t854 - t859 - t1147 - t1150 + t565 + 0.3e1 / 0.4e1 * t708 - 0.3e1 / 0.4e1 * t797 - t891;
1414   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1415     v4rho4[3] = t1975 + t1977 + t1987 + t1990;
1416 
1417   t2003 = 0.4e1 / 0.3e1 * t663 + 0.4e1 / 0.3e1 * t675 - 0.2e1 / 0.3e1 * t754 - 0.2e1 * t757 - 0.2e1 / 0.3e1 * t760 - 0.4e1 * t763 - 0.2e1 * t766 + 0.4e1 / 0.9e1 * t819 - 0.2e1 * t823 + 0.4e1 / 0.9e1 * t828 - 0.2e1 * t833;
1418   t2049 = t879 * rho[1];
1419   t2061 = 0.1e1 / t364 / t60;
1420   t2072 = t339 * t339;
1421   t2082 = 0.12e2 * t539 + 0.12e2 * t1203;
1422   t2107 = 0.2e1 / 0.3e1 * t613 * t272 * t319 * t325 - 0.2e1 / 0.3e1 * t271 * t272 * t121 * t838 + 0.2e1 / 0.3e1 * t628 * t295 * t330 * t339 - 0.2e1 * t294 * t734 * t831 - 0.2e1 / 0.3e1 * t294 * t295 * t126 * t855 - 0.2e1 / 0.3e1 * t1864 - 0.2e1 / 0.3e1 * t1872 - 0.3e1 / 0.8e1 * t7 * t45 * (0.167200e6 / 0.81e2 * t19 * t46 / t50 / t139 * t61 - 0.948640e6 / 0.81e2 * t93 * t137 / t360 * t143 + 0.629120e6 / 0.27e2 * t232 * t359 / t55 / t1087 * t365 - 0.532480e6 / 0.27e2 * t588 * t878 / t50 / t2049 * t884 + 0.163840e6 / 0.27e2 * t1320 * t359 * t46 / t879 / t48 * t2061) - t1882 / 0.2e1 + t1645 / 0.2e1 - 0.5e1 / 0.9e1 * t1356 - t7 * t249 * t64 * t2072 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t502 * t147 * t825 - t7 * t70 * t38 * t2082 / 0.2e1 - t7 * t249 * t369 * t330 - 0.2e1 * t7 * t108 * t888 * t126 - 0.3e1 * t7 * t108 * t369 * t339 - 0.2e1 * t7 * t108 * t147 * t855 + t7 * t108 * t64 * t2082 / 0.2e1;
1423   t2108 = t319 * t319;
1424   t2113 = t325 * t325;
1425   t2118 = t330 * t330;
1426   t2137 = -0.5e1 / 0.27e2 * t7 * t1221 * t38 * t2108 - t7 * t187 * t38 * t2113 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1247 * t64 * t2118 + 0.4e1 / 0.27e2 * t1887 + 0.4e1 / 0.27e2 * t1890 - 0.2e1 / 0.3e1 * t1893 - 0.2e1 / 0.3e1 * t1896 - 0.2e1 * t1899 - 0.2e1 * t1902 - 0.2e1 / 0.3e1 * t1907 + 0.2e1 / 0.9e1 * t1651 + 0.2e1 / 0.3e1 * t1654 + 0.2e1 / 0.9e1 * t1657 + 0.4e1 / 0.3e1 * t1663 + 0.2e1 / 0.3e1 * t1666 - 0.20e2 / 0.27e2 * t1363 - 0.20e2 / 0.27e2 * t1411 + t1337 + t1341;
1427   t2140 = -t1147 - t1150 + t708 - 0.3e1 / 0.2e1 * t797 - 0.2e1 * t841 - 0.2e1 * t845 - 0.6e1 * t849 - 0.6e1 * t853 - 0.2e1 * t858 - 0.3e1 / 0.2e1 * t890 + t9 * (t2107 + t2137);
1428   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1429     v4rho4[4] = t2003 + t2140;
1430 
1431   t2142 = t7 * t82 * t398;
1432   t2145 = t7 * t17 * t933;
1433   t2174 = t7 * t70 * t907;
1434   t2177 = t7 * t70 * t910;
1435   t2181 = t7 * t210 * t165 / 0.9e1;
1436   t2199 = t7 * t187 * t903;
1437   t2202 = t7 * t193 * t379;
1438   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1439     v4rho3sigma[0] = -t2142 / 0.2e1 - 0.3e1 / 0.2e1 * t2145 - 0.3e1 / 0.8e1 * t155 * t156 * (-0.13376e5 / 0.27e2 * t19 * sigma[0] * t570 * t35 + 0.89888e5 / 0.27e2 * t93 * t20 * t575 * t100 - 0.66976e5 / 0.9e1 * t232 * t94 * t582 * t239 + 0.61952e5 / 0.9e1 * t588 * t233 * t592 * t595 - 0.20480e5 / 0.9e1 * t1320 * t589 / t590 / t217 * t1326) - 0.4e1 * t2174 - 0.2e1 * t2177 + t2181 - t155 * t270 * t398 * t188 / 0.2e1 - 0.3e1 / 0.2e1 * t155 * t378 * t933 * t75 - 0.3e1 / 0.2e1 * t155 * t378 * t398 * t203 - t155 * t378 * t165 * t542 / 0.2e1 - 0.2e1 / 0.3e1 * t2199 - 0.2e1 / 0.3e1 * t2202 + t155 * t612 * t165 * t483 / 0.9e1 - t958 * t959 * t492 / 0.2e1;
1440 
1441   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1442     v4rho3sigma[1] = 0.0e0;
1443 
1444   t2212 = t7 * t249 * t943;
1445   t2215 = t7 * t255 * t405;
1446   t2218 = t7 * t108 * t947;
1447   t2222 = t7 * t264 * t178 / 0.9e1;
1448   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1449     v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2212 - 0.2e1 / 0.3e1 * t2215 - 0.2e1 * t2218 + t2222 + t155 * t627 * t178 * t503 / 0.9e1 - t958 * t976 * t512 / 0.2e1 - t155 * t404 * t178 * t554 / 0.2e1;
1450 
1451   t2239 = 0.4e1 / 0.3e1 * t954 * t964;
1452   t2250 = t81 * t1 * t3;
1453   t2251 = t2250 * t410;
1454   t2253 = t953 * t957;
1455   t2255 = 0.4e1 / 0.9e1 * t2253 * t960;
1456   t2257 = t285 * t1 * t3;
1457   t2258 = t2257 * t970;
1458   t2270 = t481 * t165;
1459   t2274 = t967 * t957;
1460   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1461     v4rho3sigma[3] = -t2142 / 0.3e1 - t2145 / 0.2e1 - 0.4e1 / 0.3e1 * t2174 - 0.2e1 / 0.3e1 * t2177 + t2181 - t2239 - t155 * t378 * t933 * t121 / 0.2e1 - t968 * t378 * t165 / 0.2e1 - 0.2e1 / 0.9e1 * t2199 - 0.4e1 / 0.9e1 * t2202 - 0.2e1 / 0.9e1 * t2251 - t2255 + t2258 / 0.6e1 - t958 * t186 * t398 * t273 / 0.3e1 - t958 * t959 * t641 / 0.6e1 - t968 * t378 * t398 * t8 + t958 * t2270 * t614 / 0.9e1 - t2274 * t959 * t623 / 0.3e1;
1462 
1463   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1464     v4rho3sigma[4] = 0.0e0;
1465 
1466   t2281 = t2250 * t415;
1467   t2284 = 0.4e1 / 0.9e1 * t2253 * t977;
1468   t2285 = t2257 * t981;
1469   t2287 = t501 * t178;
1470   t2300 = t2250 * t434;
1471   t2303 = 0.4e1 / 0.3e1 * t954 * t987;
1472   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1473     v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2212 - 0.4e1 / 0.9e1 * t2215 - 0.2e1 / 0.3e1 * t2218 + t2222 - 0.2e1 / 0.9e1 * t2281 - t2284 - t2285 / 0.6e1 + t958 * t2287 * t629 / 0.9e1 + t2274 * t976 * t650 / 0.3e1 - t958 * t976 * t645 / 0.6e1 + t968 * t404 * t178 / 0.2e1 - t2300 / 0.6e1 - t2303 - t155 * t293 * t433 * t250 / 0.6e1 - t155 * t404 * t433 * t259 / 0.2e1;
1474 
1475   t2316 = t954 * t992;
1476   t2328 = t954 * t996;
1477   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1478     v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2251 - t2255 - t2239 - 0.4e1 / 0.3e1 * t2258 - 0.2e1 / 0.9e1 * t2202 + t2181 - t2142 / 0.6e1 - 0.2e1 / 0.9e1 * t2316 + t958 * t2270 * t718 / 0.9e1 - t155 * t270 * t398 * t319 / 0.6e1 - t2274 * t959 * t742 / 0.3e1 - 0.2e1 / 0.3e1 * t2328 - t958 * t959 * t722 / 0.6e1 - t155 * t378 * t398 * t325 / 0.2e1 - t155 * t378 * t165 * t777 / 0.2e1;
1479 
1480   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1481     v4rho3sigma[7] = 0.0e0;
1482 
1483   t2345 = t954 * t1001;
1484   t2353 = t954 * t1005;
1485   t2362 = t954 * t1008;
1486   t2371 = t954 * t1031;
1487   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1488     v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2281 - t2284 + 0.4e1 / 0.3e1 * t2285 - 0.2e1 / 0.9e1 * t2215 + t2222 - t2300 / 0.3e1 - t2303 - 0.2e1 / 0.9e1 * t2345 + t958 * t2287 * t726 / 0.9e1 + t2274 * t976 * t730 / 0.3e1 - 0.4e1 / 0.3e1 * t2353 - t958 * t248 * t433 * t296 / 0.3e1 + t968 * t404 * t433 * t8 - 0.2e1 / 0.3e1 * t2362 - t958 * t976 * t738 / 0.6e1 - t155 * t404 * t178 * t785 / 0.2e1 - t2371 / 0.2e1 - t155 * t404 * t1030 * t109 / 0.2e1;
1489 
1490   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1491     v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2251 - 0.2e1 / 0.3e1 * t2316 - 0.2e1 * t2328 + t2181 + t155 * t612 * t165 * t816 / 0.9e1 - t958 * t959 * t821 / 0.2e1 - t155 * t378 * t165 * t838 / 0.2e1;
1492 
1493   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1494     v4rho3sigma[10] = 0.0e0;
1495 
1496   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1497     v4rho3sigma[11] = -0.3e1 / 0.8e1 * t155 * t169 * (-0.13376e5 / 0.27e2 * t19 * sigma[2] * t861 * t61 + 0.89888e5 / 0.27e2 * t93 * t46 * t866 * t143 - 0.66976e5 / 0.9e1 * t232 * t137 * t873 * t365 + 0.61952e5 / 0.9e1 * t588 * t359 * t881 * t884 - 0.20480e5 / 0.9e1 * t1320 * t878 / t879 / t345 * t2061) + t2222 - t2300 / 0.2e1 - 0.3e1 / 0.2e1 * t2371 - t155 * t293 * t433 * t330 / 0.2e1 - 0.3e1 / 0.2e1 * t155 * t404 * t1030 * t126 - 0.3e1 / 0.2e1 * t155 * t404 * t433 * t339 - t155 * t404 * t178 * t855 / 0.2e1 - 0.2e1 / 0.3e1 * t2281 - 0.2e1 / 0.3e1 * t2345 - 0.4e1 * t2353 - 0.2e1 * t2362 + t155 * t627 * t178 * t825 / 0.9e1 - t958 * t976 * t831 / 0.2e1;
1498 
1499   t2447 = t2250 * t451 / 0.6e1;
1500   t2448 = t954 * t1037;
1501   t2450 = t954 * t1059;
1502   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1503     v4rho2sigma2[0] = -t2447 - 0.4e1 / 0.3e1 * t2448 - t2450 - t155 * t270 * t450 * t188 / 0.6e1 - t155 * t378 * t1058 * t75 - t155 * t378 * t450 * t203 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * (0.608e3 / 0.9e1 * t19 * t220 * t35 - 0.6560e4 / 0.9e1 * t93 * t226 * t100 * sigma[0] + 0.6272e4 / 0.3e1 * t232 * t236 * t239 * t20 - 0.6848e4 / 0.3e1 * t588 * t94 * t928 * t595 + 0.2560e4 / 0.3e1 * t1320 * t233 / t590 / t21 * t1326);
1504 
1505   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1506     v4rho2sigma2[1] = 0.0e0;
1507 
1508   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1509     v4rho2sigma2[2] = 0.0e0;
1510 
1511   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1512     v4rho2sigma2[3] = 0.0e0;
1513 
1514   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1515     v4rho2sigma2[4] = 0.0e0;
1516 
1517   t2488 = t2250 * t468 / 0.6e1;
1518   t2489 = t954 * t1064;
1519   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1520     v4rho2sigma2[5] = -t2488 - 0.4e1 / 0.3e1 * t2489 - t155 * t293 * t467 * t250 / 0.6e1 - t155 * t404 * t467 * t259 / 0.2e1;
1521 
1522   t2501 = t954 * t1068;
1523   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1524     v4rho2sigma2[6] = -t2447 - 0.2e1 / 0.3e1 * t2448 - t2450 / 0.2e1 - 0.2e1 / 0.3e1 * t2501 - t958 * t186 * t450 * t273 / 0.6e1 - t155 * t378 * t1058 * t121 / 0.2e1 - t968 * t378 * t450 * t8 / 0.2e1;
1525 
1526   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1527     v4rho2sigma2[7] = 0.0e0;
1528 
1529   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1530     v4rho2sigma2[8] = 0.0e0;
1531 
1532   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1533     v4rho2sigma2[9] = 0.0e0;
1534 
1535   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1536     v4rho2sigma2[10] = 0.0e0;
1537 
1538   t2516 = t954 * t1073;
1539   t2526 = t954 * t1095;
1540   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1541     v4rho2sigma2[11] = -t2488 - 0.2e1 / 0.3e1 * t2489 - 0.2e1 / 0.3e1 * t2516 - t958 * t248 * t467 * t296 / 0.6e1 + t968 * t404 * t467 * t8 / 0.2e1 - t2526 / 0.2e1 - t155 * t404 * t1094 * t109 / 0.2e1;
1542 
1543   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1544     v4rho2sigma2[12] = -t2447 - 0.4e1 / 0.3e1 * t2501 - t155 * t270 * t450 * t319 / 0.6e1 - t155 * t378 * t450 * t325 / 0.2e1;
1545 
1546   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1547     v4rho2sigma2[13] = 0.0e0;
1548 
1549   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1550     v4rho2sigma2[14] = 0.0e0;
1551 
1552   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1553     v4rho2sigma2[15] = 0.0e0;
1554 
1555   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556     v4rho2sigma2[16] = 0.0e0;
1557 
1558   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559     v4rho2sigma2[17] = -t2488 - 0.4e1 / 0.3e1 * t2516 - t2526 - t155 * t293 * t467 * t330 / 0.6e1 - t155 * t404 * t1094 * t126 - t155 * t404 * t467 * t339 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t169 * (0.608e3 / 0.9e1 * t19 * t348 * t61 - 0.6560e4 / 0.9e1 * t93 * t354 * t143 * sigma[2] + 0.6272e4 / 0.3e1 * t232 * t362 * t365 * t46 - 0.6848e4 / 0.3e1 * t588 * t137 * t1025 * t884 + 0.2560e4 / 0.3e1 * t1320 * t359 / t879 / t47 * t2061);
1560 
1561   t2578 = t954 * t1112;
1562   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1563     v4rhosigma3[0] = -t2578 / 0.2e1 - t155 * t378 * t1111 * t75 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * (0.96e2 * t93 * t1043 - 0.480e3 * t232 * t393 * t239 * sigma[0] + 0.704e3 * t588 * t20 * t1053 * t595 - 0.320e3 * t1320 * t94 / t1312 * t1326);
1564 
1565   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1566     v4rhosigma3[1] = 0.0e0;
1567 
1568   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1569     v4rhosigma3[2] = 0.0e0;
1570 
1571   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1572     v4rhosigma3[3] = 0.0e0;
1573 
1574   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1575     v4rhosigma3[4] = 0.0e0;
1576 
1577   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1578     v4rhosigma3[5] = 0.0e0;
1579 
1580   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1581     v4rhosigma3[6] = 0.0e0;
1582 
1583   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1584     v4rhosigma3[7] = 0.0e0;
1585 
1586   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1587     v4rhosigma3[8] = 0.0e0;
1588 
1589   t2603 = t954 * t1129;
1590   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1591     v4rhosigma3[9] = -t155 * t404 * t1128 * t109 / 0.2e1 - t2603 / 0.2e1;
1592 
1593   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594     v4rhosigma3[10] = -t155 * t378 * t1111 * t121 / 0.2e1 - t2578 / 0.2e1;
1595 
1596   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597     v4rhosigma3[11] = 0.0e0;
1598 
1599   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600     v4rhosigma3[12] = 0.0e0;
1601 
1602   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603     v4rhosigma3[13] = 0.0e0;
1604 
1605   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606     v4rhosigma3[14] = 0.0e0;
1607 
1608   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1609     v4rhosigma3[15] = 0.0e0;
1610 
1611   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1612     v4rhosigma3[16] = 0.0e0;
1613 
1614   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1615     v4rhosigma3[17] = 0.0e0;
1616 
1617   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1618     v4rhosigma3[18] = 0.0e0;
1619 
1620   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1621     v4rhosigma3[19] = -t2603 / 0.2e1 - t155 * t404 * t1128 * t126 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t169 * (0.96e2 * t93 * t1079 - 0.480e3 * t232 * t428 * t365 * sigma[2] + 0.704e3 * t588 * t46 * t1089 * t884 - 0.320e3 * t1320 * t137 / t2049 * t2061);
1622 
1623   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1624     v4sigma4[0] = -0.3e1 / 0.8e1 * t155 * t156 * (0.72e2 * t232 * t445 * t239 - 0.192e3 * t588 * sigma[0] * t1106 * t595 + 0.120e3 * t1320 * t20 / t590 * t1326);
1625 
1626   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1627     v4sigma4[1] = 0.0e0;
1628 
1629   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1630     v4sigma4[2] = 0.0e0;
1631 
1632   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1633     v4sigma4[3] = 0.0e0;
1634 
1635   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1636     v4sigma4[4] = 0.0e0;
1637 
1638   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1639     v4sigma4[5] = 0.0e0;
1640 
1641   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1642     v4sigma4[6] = 0.0e0;
1643 
1644   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645     v4sigma4[7] = 0.0e0;
1646 
1647   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1648     v4sigma4[8] = 0.0e0;
1649 
1650   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1651     v4sigma4[9] = 0.0e0;
1652 
1653   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1654     v4sigma4[10] = 0.0e0;
1655 
1656   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1657     v4sigma4[11] = 0.0e0;
1658 
1659   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1660     v4sigma4[12] = 0.0e0;
1661 
1662   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1663     v4sigma4[13] = 0.0e0;
1664 
1665   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1666     v4sigma4[14] = -0.3e1 / 0.8e1 * t155 * t169 * (0.72e2 * t232 * t462 * t365 - 0.192e3 * t588 * sigma[2] * t1123 * t884 + 0.120e3 * t1320 * t46 / t879 * t2061);
1667 
1668 #ifndef XC_DONT_COMPILE_MXC
1669 
1670   if(order < 5) return;
1671 
1672 
1673 #endif
1674 
1675 #endif
1676 
1677 #endif
1678 
1679 #endif
1680 
1681 #endif
1682 
1683 
1684 }
1685 
1686