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_c_wi.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, t2, t3, t4, t6, t7, t9, t12;
24   double t13, t15, t16, t17, t18, t22, t23, t25;
25   double t26, t27, t28, t31, t32, t33, t38, t42;
26   double t43;
27 
28 #ifndef XC_DONT_COMPILE_VXC
29   double t44, t46, t49, t50, t51, t53, t58, t59;
30   double t61, t62, t63, t71, t72, t73, t76, t77;
31   double t81, t83, t87, t88, t90, t91, t92, t93;
32   double t94, t95, t96, t97;
33 
34 #ifndef XC_DONT_COMPILE_FXC
35   double t102, t106, t110, t112, t117, t118, t119, t120;
36   double t121, t122, t127, t128, t133, t134, t135, t148;
37   double t151, t152, t153, t156, t157, t163, t164, t167;
38   double t168, t173, t174, t178, t179, t180, t183, t184;
39   double t187, t188, t190, t193, t194, t196, t197, t201;
40   double t204, t208, t209, t211, t212, t215, t216, t217;
41   double t218, t219, t220, t222, t225, t226, t227, t228;
42   double t229, t230, t233, t234, t235;
43 
44 #ifndef XC_DONT_COMPILE_KXC
45   double t240, t243, t249, t254, t259, t260, t265, t266;
46   double t267, t269, t270, t275, t276, t284, t285, t286;
47   double t287, t290, t291, t295, t312, t315, t316, t318;
48   double t319, t322, t325, t326, t328, t331, t332, t335;
49   double t336, t339, t341, t349, t355, t358, t360, t363;
50   double t365, t370, t371, t373, t374, t376, t379, t383;
51   double t389, t393, t399, t401, t405, t410, t411, t414;
52   double t415, t418, t419, t422, t425, t429, t430, t433;
53   double t434, t437, t438, t439, t441, t444, t447, t448;
54   double t451, t453, t457, t458, t459, t460, t461, t469;
55   double t473, t477, t481, t482, t484, t485, t488, t489;
56   double t492, t493, t498, t500, t501, t502, t506, t509;
57   double t510, t511, t512, t515, t516, t517, t521, t522;
58   double t523;
59 
60 #ifndef XC_DONT_COMPILE_LXC
61   double t562, t563, t580, t590, t595, t604, t608, t623;
62   double t624, t632, t725, t785, t872, t943, t1011, t1071;
63   double t1118;
64 #endif
65 
66 #endif
67 
68 #endif
69 
70 #endif
71 
72 #endif
73 
74 
75   gga_c_wi_params *params;
76 
77   assert(p->params != NULL);
78   params = (gga_c_wi_params * )(p->params);
79 
80   t1 = params->b * sigma[0];
81   t2 = rho[0] * rho[0];
82   t3 = POW_1_3(rho[0]);
83   t4 = t3 * t3;
84   t6 = 0.1e1 / t4 / t2;
85   t7 = params->k * sigma[0];
86   t9 = exp(-t7 * t6);
87   t12 = t1 * t6 * t9 + params->a;
88   t13 = M_CBRT3;
89   t15 = POW_1_3(0.1e1 / M_PI);
90   t16 = t13 * t15;
91   t17 = M_CBRT4;
92   t18 = t17 * t17;
93   t22 = t13 * t13;
94   t23 = M_CBRTPI;
95   t25 = sqrt(sigma[0]);
96   t26 = t25 * sigma[0];
97   t27 = t2 * t2;
98   t28 = 0.1e1 / t27;
99   t31 = 0.1e1 / t3 / rho[0];
100   t32 = t25 * t31;
101   t33 = sqrt(t32);
102   t38 = 0.1e1 + params->d * t17 * t22 * t23 * t33 * t26 * t28 / 0.3e1;
103   t42 = params->c + t16 * t18 / t3 * t38 / 0.4e1;
104   t43 = 0.1e1 / t42;
105   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
106     zk[0] = t12 * t43;
107 
108 #ifndef XC_DONT_COMPILE_VXC
109 
110   if(order < 1) return;
111 
112 
113   t44 = t2 * rho[0];
114   t46 = 0.1e1 / t4 / t44;
115   t49 = sigma[0] * sigma[0];
116   t50 = params->b * t49;
117   t51 = t27 * t2;
118   t53 = 0.1e1 / t3 / t51;
119   t58 = 0.8e1 / 0.3e1 * t50 * t53 * params->k * t9 - 0.8e1 / 0.3e1 * t1 * t46 * t9;
120   t59 = rho[0] * t58;
121   t61 = rho[0] * t12;
122   t62 = t42 * t42;
123   t63 = 0.1e1 / t62;
124   t71 = t33 * sigma[0] * t6;
125   t72 = t23 * t71;
126   t73 = t72 * t25;
127   t76 = -t16 * t18 * t31 * t38 / 0.12e2 - 0.14e2 / 0.3e1 * t15 * t6 * params->d * t73;
128   t77 = t63 * t76;
129   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
130     vrho[0] = t59 * t43 - t61 * t77 + (t12 * t43);
131 
132   t81 = t27 * rho[0];
133   t83 = 0.1e1 / t3 / t81;
134   t87 = -t1 * t83 * params->k * t9 + params->b * t6 * t9;
135   t88 = rho[0] * t87;
136   t90 = 0.1e1 / t4;
137   t91 = t90 * t12;
138   t92 = t63 * t15;
139   t93 = t91 * t92;
140   t94 = params->d * t23;
141   t95 = 0.1e1 / t25;
142   t96 = t71 * t95;
143   t97 = t94 * t96;
144   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
145     vsigma[0] = t88 * t43 - 0.7e1 / 0.4e1 * t93 * t97;
146 
147 #ifndef XC_DONT_COMPILE_FXC
148 
149   if(order < 2) return;
150 
151 
152   t102 = t12 * t63;
153   t106 = 0.1e1 / t4 / t27;
154   t110 = t27 * t44;
155   t112 = 0.1e1 / t3 / t110;
156   t117 = t49 * sigma[0];
157   t118 = params->b * t117;
158   t119 = t27 * t27;
159   t120 = t119 * t2;
160   t121 = 0.1e1 / t120;
161   t122 = params->k * params->k;
162   t127 = 0.88e2 / 0.9e1 * t1 * t106 * t9 - 0.24e2 * t50 * t112 * params->k * t9 + 0.64e2 / 0.9e1 * t118 * t121 * t122 * t9;
163   t128 = rho[0] * t127;
164   t133 = 0.1e1 / t62 / t42;
165   t134 = t76 * t76;
166   t135 = t133 * t134;
167   t148 = 0.1e1 / t81;
168   t151 = t33 * t32;
169   t152 = t23 * t151;
170   t153 = t152 * sigma[0];
171   t156 = t16 * t18 / t3 / t2 * t38 / 0.9e1 + 0.14e2 * t15 * t46 * params->d * t73 + 0.140e3 / 0.9e1 * t15 * t148 * params->d * t153;
172   t157 = t63 * t156;
173   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
174     v2rho2[0] = -0.2e1 * t102 * t76 + t128 * t43 + 0.2e1 * t61 * t135 - t61 * t157 + 0.2e1 * t58 * t43 - 0.2e1 * t59 * t77;
175 
176   t163 = params->b * t53;
177   t164 = t7 * t9;
178   t167 = t119 * rho[0];
179   t168 = 0.1e1 / t167;
180   t173 = -0.8e1 / 0.3e1 * params->b * t46 * t9 + 0.8e1 * t163 * t164 - 0.8e1 / 0.3e1 * t50 * t168 * t122 * t9;
181   t174 = rho[0] * t173;
182   t178 = 0.1e1 / t4 / rho[0];
183   t179 = t178 * t12;
184   t180 = t179 * t92;
185   t183 = t90 * t58;
186   t184 = t183 * t92;
187   t187 = t133 * t15;
188   t188 = t91 * t187;
189   t190 = t94 * t96 * t76;
190   t193 = 0.1e1 / t44;
191   t194 = t193 * t12;
192   t196 = t15 * params->d;
193   t197 = t196 * t152;
194   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
195     v2rhosigma[0] = t87 * t43 + t174 * t43 - t88 * t77 + 0.7e1 / 0.6e1 * t180 * t97 - 0.7e1 / 0.4e1 * t184 * t97 + 0.7e1 / 0.2e1 * t188 * t190 + 0.35e2 / 0.6e1 * t194 * t63 * t197;
196 
197   t201 = params->k * t9;
198   t204 = 0.1e1 / t119;
199   t208 = t1 * t204 * t122 * t9 - 0.2e1 * params->b * t83 * t201;
200   t209 = rho[0] * t208;
201   t211 = t90 * t87;
202   t212 = t211 * t92;
203   t215 = t168 * t12;
204   t216 = t215 * t133;
205   t217 = t15 * t15;
206   t218 = params->d * params->d;
207   t219 = t217 * t218;
208   t220 = t23 * t23;
209   t222 = t219 * t220 * t26;
210   t225 = 0.1e1 / t2;
211   t226 = t225 * t12;
212   t227 = t226 * t92;
213   t228 = 0.1e1 / sigma[0];
214   t229 = t151 * t228;
215   t230 = t94 * t229;
216   t233 = 0.1e1 / t26;
217   t234 = t71 * t233;
218   t235 = t94 * t234;
219   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
220     v2sigma2[0] = t209 * t43 - 0.7e1 / 0.2e1 * t212 * t97 + 0.49e2 / 0.8e1 * t216 * t222 - 0.35e2 / 0.16e2 * t227 * t230 + 0.7e1 / 0.8e1 * t93 * t235;
221 
222 #ifndef XC_DONT_COMPILE_KXC
223 
224   if(order < 3) return;
225 
226 
227   t240 = t58 * t63;
228   t243 = t12 * t133;
229   t249 = 0.1e1 / t4 / t81;
230   t254 = 0.1e1 / t3 / t119;
231   t259 = t119 * t44;
232   t260 = 0.1e1 / t259;
233   t265 = t49 * t49;
234   t266 = params->b * t265;
235   t267 = t119 * t81;
236   t269 = 0.1e1 / t4 / t267;
237   t270 = t122 * params->k;
238   t275 = -0.1232e4 / 0.27e2 * t1 * t249 * t9 + 0.5456e4 / 0.27e2 * t50 * t254 * params->k * t9 - 0.1216e4 / 0.9e1 * t118 * t260 * t122 * t9 + 0.512e3 / 0.27e2 * t266 * t269 * t270 * t9;
239   t276 = rho[0] * t275;
240   t284 = t62 * t62;
241   t285 = 0.1e1 / t284;
242   t286 = t134 * t76;
243   t287 = t285 * t286;
244   t290 = t133 * t76;
245   t291 = t290 * t156;
246   t295 = 0.1e1 / t3 / t44;
247   t312 = t23 * t33 * t26;
248   t315 = -0.7e1 / 0.27e2 * t16 * t18 * t295 * t38 - 0.1442e4 / 0.27e2 * t15 * t106 * params->d * t73 - 0.1120e4 / 0.9e1 * t15 / t51 * params->d * t153 - 0.280e3 / 0.9e1 * t15 * t112 * params->d * t312;
249   t316 = t63 * t315;
250   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
251     v3rho3[0] = -0.3e1 * t102 * t156 + 0.3e1 * t127 * t43 - 0.3e1 * t128 * t77 + 0.6e1 * t243 * t134 + 0.6e1 * t59 * t135 - 0.3e1 * t59 * t157 - 0.6e1 * t240 * t76 + t276 * t43 - 0.6e1 * t61 * t287 + 0.6e1 * t61 * t291 - t61 * t316;
252 
253   t318 = t6 * t12;
254   t319 = t318 * t92;
255   t322 = t179 * t187;
256   t325 = t285 * t15;
257   t326 = t91 * t325;
258   t328 = t94 * t96 * t134;
259   t331 = t90 * t127;
260   t332 = t331 * t92;
261   t335 = t178 * t58;
262   t336 = t335 * t92;
263   t339 = t194 * t187;
264   t341 = t94 * t151 * t76;
265   t349 = t87 * t63;
266   t355 = params->b * t112;
267   t358 = params->b * t121;
268   t360 = t122 * t49 * t9;
269   t363 = t119 * t27;
270   t365 = 0.1e1 / t4 / t363;
271   t370 = 0.88e2 / 0.9e1 * params->b * t106 * t9 - 0.520e3 / 0.9e1 * t355 * t164 + 0.136e3 / 0.3e1 * t358 * t360 - 0.64e2 / 0.9e1 * t118 * t365 * t270 * t9;
272   t371 = rho[0] * t370;
273   t373 = t83 * t12;
274   t374 = t373 * t92;
275   t376 = t94 * t33 * t25;
276   t379 = t28 * t12;
277   t383 = t193 * t58;
278   t389 = t183 * t187;
279   t393 = t94 * t96 * t156;
280   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
281     v3rho2sigma[0] = -0.35e2 / 0.18e2 * t319 * t97 - 0.14e2 / 0.3e1 * t322 * t190 - 0.21e2 / 0.2e1 * t326 * t328 - 0.7e1 / 0.4e1 * t332 * t97 + 0.7e1 / 0.3e1 * t336 * t97 - 0.70e2 / 0.3e1 * t339 * t341 + 0.2e1 * t88 * t135 - 0.2e1 * t174 * t77 - t88 * t157 - 0.2e1 * t349 * t76 + t371 * t43 - 0.35e2 / 0.3e1 * t374 * t376 - 0.385e3 / 0.18e2 * t379 * t63 * t197 + 0.35e2 / 0.3e1 * t383 * t63 * t197 + 0.2e1 * t173 * t43 + 0.7e1 * t389 * t190 + 0.7e1 / 0.2e1 * t188 * t393;
282 
283   t399 = params->b * t168;
284   t401 = t122 * sigma[0] * t9;
285   t405 = 0.1e1 / t4 / t259;
286   t410 = 0.32e2 / 0.3e1 * t163 * t201 - 0.40e2 / 0.3e1 * t399 * t401 + 0.8e1 / 0.3e1 * t50 * t405 * t270 * t9;
287   t411 = rho[0] * t410;
288   t414 = t178 * t87;
289   t415 = t414 * t92;
290   t418 = t90 * t173;
291   t419 = t418 * t92;
292   t422 = t211 * t187;
293   t425 = t193 * t87;
294   t429 = t121 * t12;
295   t430 = t429 * t133;
296   t433 = t168 * t58;
297   t434 = t433 * t133;
298   t437 = t285 * t217;
299   t438 = t215 * t437;
300   t439 = t218 * t220;
301   t441 = t439 * t26 * t76;
302   t444 = t194 * t92;
303   t447 = t225 * t58;
304   t448 = t447 * t92;
305   t451 = t226 * t187;
306   t453 = t94 * t229 * t76;
307   t457 = 0.1e1 / t3 / t27;
308   t458 = t457 * t12;
309   t459 = t458 * t92;
310   t460 = t33 * t95;
311   t461 = t94 * t460;
312   t469 = t94 * t234 * t76;
313   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
314     v3rhosigma2[0] = t208 * t43 + t411 * t43 - t209 * t77 + 0.7e1 / 0.3e1 * t415 * t97 - 0.7e1 / 0.2e1 * t419 * t97 + 0.7e1 * t422 * t190 + 0.35e2 / 0.3e1 * t425 * t63 * t197 - 0.441e3 / 0.8e1 * t430 * t222 + 0.49e2 / 0.8e1 * t434 * t222 - 0.147e3 / 0.8e1 * t438 * t441 + 0.35e2 / 0.24e2 * t444 * t230 - 0.35e2 / 0.16e2 * t448 * t230 + 0.35e2 / 0.8e1 * t451 * t453 + 0.35e2 / 0.8e1 * t459 * t461 - 0.7e1 / 0.12e2 * t180 * t235 + 0.7e1 / 0.8e1 * t184 * t235 - 0.7e1 / 0.4e1 * t188 * t469;
315 
316   t473 = t122 * t9;
317   t477 = 0.1e1 / t4 / t120;
318   t481 = -t1 * t477 * t270 * t9 + 0.3e1 * params->b * t204 * t473;
319   t482 = rho[0] * t481;
320   t484 = t90 * t208;
321   t485 = t484 * t92;
322   t488 = t168 * t87;
323   t489 = t488 * t133;
324   t492 = t225 * t87;
325   t493 = t492 * t92;
326   t498 = t477 * t12;
327   t500 = t218 * params->d;
328   t501 = t500 * sigma[0];
329   t502 = t501 * t71;
330   t506 = t219 * t220 * t25;
331   t509 = t295 * t12;
332   t510 = t509 * t92;
333   t511 = t33 * t233;
334   t512 = t94 * t511;
335   t515 = 0.1e1 / t49;
336   t516 = t151 * t515;
337   t517 = t94 * t516;
338   t521 = 0.1e1 / t25 / t49;
339   t522 = t71 * t521;
340   t523 = t94 * t522;
341   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
342     v3sigma3[0] = t482 * t43 - 0.21e2 / 0.4e1 * t485 * t97 + 0.147e3 / 0.8e1 * t489 * t222 - 0.105e3 / 0.16e2 * t493 * t230 + 0.21e2 / 0.8e1 * t212 * t235 - 0.1029e4 / 0.32e2 * t498 * t285 * t502 + 0.441e3 / 0.32e2 * t216 * t506 - 0.105e3 / 0.64e2 * t510 * t512 + 0.105e3 / 0.32e2 * t227 * t517 - 0.21e2 / 0.16e2 * t93 * t523;
343 
344 #ifndef XC_DONT_COMPILE_LXC
345 
346   if(order < 4) return;
347 
348 
349   t562 = 0.1e1 / t33;
350   t563 = t23 * t562;
351   t580 = 0.1e1 / t4 / t51;
352   t590 = 0.1e1 / t363;
353   t595 = t119 * t51;
354   t604 = t119 * t119;
355   t608 = t122 * t122;
356   t623 = 0.1e1 / t284 / t42;
357   t624 = t134 * t134;
358   t632 = t156 * t156;
359   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
360     v4rho4[0] = 0.24e2 * t58 * t133 * t134 - 0.24e2 * t12 * t285 * t286 + 0.24e2 * t243 * t76 * t156 - 0.4e1 * t276 * t77 - 0.6e1 * t128 * t157 - 0.4e1 * t59 * t316 - t61 * t63 * (0.70e2 / 0.81e2 * t16 * t18 * t457 * t38 + 0.6860e4 / 0.27e2 * t15 * t249 * params->d * t73 + 0.74900e5 / 0.81e2 * t15 / t110 * params->d * t153 + 0.12880e5 / 0.27e2 * t15 * t254 * params->d * t312 + 0.560e3 / 0.27e2 * t15 / t4 / t167 * params->d * t563 * t49) + 0.4e1 * t275 * t43 - 0.12e2 * t127 * t63 * t76 - 0.12e2 * t240 * t156 - 0.4e1 * t102 * t315 + rho[0] * (0.20944e5 / 0.81e2 * t1 * t580 * t9 - 0.48752e5 / 0.27e2 * t50 / t3 / t167 * params->k * t9 + 0.164032e6 / 0.81e2 * t118 * t590 * t122 * t9 - 0.50176e5 / 0.81e2 * t266 / t4 / t595 * t270 * t9 + 0.4096e4 / 0.81e2 * params->b * t265 * sigma[0] / t3 / t604 / rho[0] * t608 * t9) * t43 + 0.12e2 * t128 * t135 - 0.24e2 * t59 * t287 + 0.24e2 * t59 * t291 + 0.24e2 * t61 * t623 * t624 - 0.36e2 * t61 * t285 * t134 * t156 + 0.6e1 * t61 * t133 * t632 + 0.8e1 * t61 * t290 * t315;
361 
362   t725 = 0.6e1 * t88 * t291 - 0.3e1 * t349 * t156 + rho[0] * (-0.1232e4 / 0.27e2 * params->b * t249 * t9 + 0.4048e4 / 0.9e1 * params->b * t254 * t164 - 0.16400e5 / 0.27e2 * params->b * t260 * t360 + 0.5696e4 / 0.27e2 * params->b * t269 * t270 * t117 * t9 - 0.512e3 / 0.27e2 * t266 / t3 / t604 * t608 * t9) * t43 + 0.6e1 * t87 * t133 * t134 - 0.6e1 * t173 * t63 * t76 + 0.3e1 * t370 * t43 - 0.35e2 / 0.6e1 * t6 * t58 * t92 * t97 + 0.385e3 / 0.3e1 * t379 * t187 * t341 + 0.105e3 * t194 * t325 * t94 * t151 * t134 - 0.7e1 / 0.4e1 * t90 * t275 * t92 * t97 - 0.70e2 * t383 * t187 * t341 - 0.35e2 * t339 * t94 * t151 * t156 - 0.35e2 * t83 * t58 * t92 * t376 + 0.70e2 / 0.9e1 / t4 / t110 * t12 * t92 * t94 * t562 * sigma[0] + 0.140e3 / 0.27e2 * t46 * t12 * t92 * t97 + 0.7e1 / 0.2e1 * t178 * t127 * t92 * t97 + 0.105e3 * t53 * t12 * t92 * t376 + 0.2485e4 / 0.27e2 * t148 * t12 * t63 * t197;
363   t785 = 0.35e2 / 0.2e1 * t193 * t127 * t63 * t197 - 0.385e3 / 0.6e1 * t28 * t58 * t63 * t197 + 0.35e2 / 0.3e1 * t318 * t187 * t190 - 0.14e2 * t335 * t187 * t190 - 0.7e1 * t322 * t393 - 0.63e2 / 0.2e1 * t183 * t325 * t328 + 0.21e2 / 0.2e1 * t331 * t187 * t190 + 0.70e2 * t373 * t187 * t94 * t33 * t76 * t25 + 0.21e2 / 0.2e1 * t389 * t393 + 0.7e1 / 0.2e1 * t188 * t94 * t96 * t315 + 0.21e2 * t179 * t325 * t328 + 0.42e2 * t91 * t623 * t15 * t94 * t96 * t286 - 0.63e2 / 0.2e1 * t91 * t325 * params->d * t72 * t95 * t76 * t156 + 0.6e1 * t174 * t135 - 0.3e1 * t371 * t77 - 0.3e1 * t174 * t157 - t88 * t316 - 0.6e1 * t88 * t287;
364   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
365     v4rho3sigma[0] = t725 + t785;
366 
367   t872 = 0.147e3 / 0.2e1 * t215 * t623 * t217 * t439 * t26 * t134 - 0.175e3 / 0.72e2 * t379 * t92 * t230 + 0.35e2 / 0.36e2 * t319 * t235 + 0.14e2 / 0.3e1 * t178 * t173 * t92 * t97 - 0.7e1 / 0.2e1 * t90 * t370 * t92 * t97 - 0.140e3 / 0.3e1 * t425 * t187 * t341 - 0.70e2 / 0.3e1 * t83 * t87 * t92 * t376 + 0.1323e4 / 0.4e1 * t429 * t437 * t441 - 0.147e3 / 0.4e1 * t433 * t437 * t441 - 0.147e3 / 0.8e1 * t438 * t439 * t26 * t156 + 0.2205e4 / 0.4e1 * t260 * t12 * t133 * t222;
368   t943 = 0.35e2 / 0.8e1 * t451 * t94 * t229 * t156 - 0.35e2 / 0.2e1 * t458 * t187 * t94 * t460 * t76 + 0.7e1 / 0.3e1 * t322 * t469 - 0.7e1 / 0.2e1 * t389 * t469 - 0.7e1 / 0.4e1 * t188 * t94 * t234 * t156 - 0.21e2 * t211 * t325 * t328 - 0.105e3 / 0.8e1 * t226 * t325 * t94 * t229 * t134 + 0.21e2 / 0.4e1 * t326 * t94 * t234 * t134 - 0.2e1 * t411 * t77 - t209 * t157 + 0.2e1 * t209 * t135;
369   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
370     v4rho2sigma2[0] = -0.2e1 * t208 * t63 * t76 + rho[0] * (-0.608e3 / 0.9e1 * t355 * t201 + 0.1336e4 / 0.9e1 * t358 * t401 - 0.200e3 / 0.3e1 * params->b * t365 * t270 * t49 * t9 + 0.64e2 / 0.9e1 * t118 / t3 / t119 / t110 * t608 * t9) * t43 + 0.2e1 * t410 * t43 + 0.35e2 / 0.12e2 * t383 * t92 * t230 - 0.175e3 / 0.8e1 * t374 * t461 - 0.35e2 / 0.16e2 * t225 * t127 * t92 * t230 + 0.35e2 / 0.4e1 * t457 * t58 * t92 * t461 - 0.7e1 / 0.6e1 * t336 * t235 + 0.7e1 / 0.8e1 * t332 * t235 - 0.35e2 / 0.9e1 * t6 * t87 * t92 * t97 + t872 - 0.385e3 / 0.9e1 * t28 * t87 * t63 * t197 + 0.70e2 / 0.3e1 * t193 * t173 * t63 * t197 - 0.441e3 / 0.4e1 * t121 * t58 * t133 * t222 + 0.49e2 / 0.8e1 * t168 * t127 * t133 * t222 - 0.35e2 / 0.12e2 * t580 * t12 * t63 * t196 * t563 - 0.28e2 / 0.3e1 * t414 * t187 * t190 + 0.14e2 * t418 * t187 * t190 + 0.7e1 * t422 * t393 - 0.35e2 / 0.6e1 * t339 * t453 + 0.35e2 / 0.4e1 * t447 * t187 * t453 + t943;
371 
372   t1011 = t481 * t43 - 0.1029e4 / 0.32e2 * t477 * t58 * t285 * t502 + 0.1715e4 / 0.16e2 / t267 * t12 * t285 * t500 * t26 * t151 + 0.343e3 * t405 * t12 * t285 * t502 + rho[0] * (-0.24e2 * t399 * t473 + 0.56e2 / 0.3e1 * params->b * t405 * t270 * sigma[0] * t9 - 0.8e1 / 0.3e1 * t50 / t3 / t595 * t608 * t9) * t43 + 0.105e3 / 0.32e2 * t448 * t517 - 0.21e2 / 0.16e2 * t184 * t523 + 0.7e1 / 0.2e1 * t178 * t208 * t92 * t97 + 0.35e2 / 0.8e1 * t425 * t92 * t230 - 0.7e1 / 0.4e1 * t415 * t235 - 0.35e2 / 0.32e2 * t459 * t512 - 0.35e2 / 0.16e2 * t444 * t517 + 0.7e1 / 0.8e1 * t180 * t523 - 0.21e2 / 0.4e1 * t90 * t410 * t92 * t97 - 0.441e3 / 0.8e1 * t488 * t437 * t441 - 0.105e3 / 0.16e2 * t225 * t173 * t92 * t230 + 0.105e3 / 0.8e1 * t457 * t87 * t92 * t461;
373   t1071 = 0.21e2 / 0.8e1 * t419 * t235 - 0.1323e4 / 0.32e2 * t438 * t439 * t25 * t76 - 0.105e3 / 0.64e2 * t295 * t58 * t92 * t512 + 0.35e2 / 0.32e2 * t249 * t12 * t92 * t94 * t562 * t228 - 0.1323e4 / 0.8e1 * t121 * t87 * t133 * t222 - 0.3969e4 / 0.32e2 * t430 * t506 + 0.35e2 / 0.2e1 * t193 * t208 * t63 * t197 + 0.147e3 / 0.8e1 * t168 * t173 * t133 * t222 + 0.1029e4 / 0.8e1 * t498 * t623 * t501 * t71 * t76 + 0.441e3 / 0.32e2 * t434 * t506 + 0.21e2 / 0.2e1 * t484 * t187 * t190 + 0.105e3 / 0.8e1 * t492 * t187 * t453 - 0.21e2 / 0.4e1 * t422 * t469 + 0.105e3 / 0.32e2 * t509 * t187 * t94 * t511 * t76 - 0.105e3 / 0.16e2 * t451 * t94 * t516 * t76 + 0.21e2 / 0.8e1 * t188 * t94 * t522 * t76 - t482 * t77;
374   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
375     v4rhosigma3[0] = t1011 + t1071;
376 
377   t1118 = t218 * t218;
378   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
379     v4sigma4[0] = -0.7e1 * t90 * t481 * t92 * t97 - 0.105e3 / 0.8e1 * t225 * t208 * t92 * t230 + 0.441e3 / 0.8e1 * t489 * t506 - 0.105e3 / 0.16e2 * t295 * t87 * t92 * t512 + 0.735e3 / 0.128e3 * t216 * t219 * t220 * t95 - 0.105e3 / 0.256e3 * t106 * t12 * t92 * t94 * t562 * t515 - 0.1029e4 / 0.8e1 * t477 * t87 * t285 * t502 - 0.13377e5 / 0.128e3 * t498 * t285 * t500 * t71 - 0.5145e4 / 0.128e3 * t590 * t12 * t285 * t500 * t25 * t151 + 0.147e3 / 0.4e1 * t168 * t208 * t133 * t222 + 0.7203e4 / 0.32e2 / t604 / t44 * t12 * t623 * t1118 * t117 * t15 * t23 + 0.21e2 / 0.4e1 * t485 * t235 + 0.105e3 / 0.8e1 * t493 * t517 - 0.21e2 / 0.4e1 * t212 * t523 + 0.315e3 / 0.64e2 * t510 * t94 * t33 * t521 - 0.525e3 / 0.64e2 * t227 * t94 * t151 / t117 + 0.105e3 / 0.32e2 * t93 * t94 * t71 / t25 / t117 + rho[0] * (-0.4e1 * params->b * t477 * t270 * t9 + t1 / t3 / t267 * t608 * t9) * t43;
380 
381 #ifndef XC_DONT_COMPILE_MXC
382 
383   if(order < 5) return;
384 
385 
386 #endif
387 
388 #endif
389 
390 #endif
391 
392 #endif
393 
394 #endif
395 
396 
397 }
398 
399 
400 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 *))401 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
402 {
403 
404 #ifndef XC_DONT_COMPILE_EXC
405   double t1, t2, t3, t4, t6, t7, t9, t12;
406   double t13, t15, t16, t17, t18, t22, t23, t25;
407   double t26, t27, t28, t31, t32, t33, t38, t42;
408   double t43;
409 
410 #ifndef XC_DONT_COMPILE_VXC
411   double t44, t46, t49, t50, t51, t53, t58, t59;
412   double t61, t62, t63, t71, t72, t73, t76, t77;
413   double t81, t83, t87, t88, t90, t91, t92, t93;
414   double t94, t95, t96, t97;
415 
416 #ifndef XC_DONT_COMPILE_FXC
417   double t102, t106, t110, t112, t117, t118, t119, t120;
418   double t121, t122, t127, t128, t133, t134, t135, t148;
419   double t151, t152, t153, t156, t157, t163, t164, t167;
420   double t168, t173, t174, t178, t179, t180, t183, t184;
421   double t187, t188, t190, t193, t194, t196, t197, t201;
422   double t204, t208, t209, t211, t212, t215, t216, t217;
423   double t218, t219, t220, t222, t225, t226, t227, t228;
424   double t229, t230, t233, t234, t235;
425 
426 #ifndef XC_DONT_COMPILE_KXC
427   double t240, t243, t249, t254, t259, t260, t265, t266;
428   double t267, t269, t270, t275, t276, t284, t285, t286;
429   double t287, t290, t291, t295, t312, t315, t316, t318;
430   double t319, t322, t325, t326, t328, t333, t335, t343;
431   double t349, t352, t354, t357, t359, t364, t365, t367;
432   double t368, t371, t372, t375, t376, t378, t381, t385;
433   double t389, t393, t399, t401, t405, t410, t411, t414;
434   double t415, t418, t419, t422, t425, t429, t430, t433;
435   double t434, t437, t438, t439, t441, t444, t447, t448;
436   double t451, t453, t457, t458, t459, t460, t461, t469;
437   double t473, t477, t481, t482, t484, t485, t488, t489;
438   double t492, t493, t498, t500, t501, t502, t506, t509;
439   double t510, t511, t512, t515, t516, t517, t521, t522;
440   double t523;
441 
442 #ifndef XC_DONT_COMPILE_LXC
443   double t562, t563, t578, t588, t593, t602, t606, t623;
444   double t624, t632, t710, t785, t869, t943, t1016, t1071;
445   double t1118;
446 #endif
447 
448 #endif
449 
450 #endif
451 
452 #endif
453 
454 #endif
455 
456 
457   gga_c_wi_params *params;
458 
459   assert(p->params != NULL);
460   params = (gga_c_wi_params * )(p->params);
461 
462   t1 = params->b * sigma[0];
463   t2 = rho[0] * rho[0];
464   t3 = POW_1_3(rho[0]);
465   t4 = t3 * t3;
466   t6 = 0.1e1 / t4 / t2;
467   t7 = params->k * sigma[0];
468   t9 = exp(-t7 * t6);
469   t12 = t1 * t6 * t9 + params->a;
470   t13 = M_CBRT3;
471   t15 = POW_1_3(0.1e1 / M_PI);
472   t16 = t13 * t15;
473   t17 = M_CBRT4;
474   t18 = t17 * t17;
475   t22 = t13 * t13;
476   t23 = M_CBRTPI;
477   t25 = sqrt(sigma[0]);
478   t26 = t25 * sigma[0];
479   t27 = t2 * t2;
480   t28 = 0.1e1 / t27;
481   t31 = 0.1e1 / t3 / rho[0];
482   t32 = t25 * t31;
483   t33 = sqrt(t32);
484   t38 = 0.1e1 + params->d * t17 * t22 * t23 * t33 * t26 * t28 / 0.3e1;
485   t42 = params->c + t16 * t18 / t3 * t38 / 0.4e1;
486   t43 = 0.1e1 / t42;
487   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
488     zk[0] = t12 * t43;
489 
490 #ifndef XC_DONT_COMPILE_VXC
491 
492   if(order < 1) return;
493 
494 
495   t44 = t2 * rho[0];
496   t46 = 0.1e1 / t4 / t44;
497   t49 = sigma[0] * sigma[0];
498   t50 = params->b * t49;
499   t51 = t27 * t2;
500   t53 = 0.1e1 / t3 / t51;
501   t58 = 0.8e1 / 0.3e1 * t50 * t53 * params->k * t9 - 0.8e1 / 0.3e1 * t1 * t46 * t9;
502   t59 = rho[0] * t58;
503   t61 = rho[0] * t12;
504   t62 = t42 * t42;
505   t63 = 0.1e1 / t62;
506   t71 = t33 * sigma[0] * t6;
507   t72 = t23 * t71;
508   t73 = t72 * t25;
509   t76 = -t16 * t18 * t31 * t38 / 0.12e2 - 0.14e2 / 0.3e1 * t15 * t6 * params->d * t73;
510   t77 = t63 * t76;
511   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
512     vrho[0] = t59 * t43 - t61 * t77 + (t12 * t43);
513 
514   t81 = t27 * rho[0];
515   t83 = 0.1e1 / t3 / t81;
516   t87 = -t1 * t83 * params->k * t9 + params->b * t6 * t9;
517   t88 = rho[0] * t87;
518   t90 = 0.1e1 / t4;
519   t91 = t90 * t12;
520   t92 = t63 * t15;
521   t93 = t91 * t92;
522   t94 = params->d * t23;
523   t95 = 0.1e1 / t25;
524   t96 = t71 * t95;
525   t97 = t94 * t96;
526   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
527     vsigma[0] = t88 * t43 - 0.7e1 / 0.4e1 * t93 * t97;
528 
529 #ifndef XC_DONT_COMPILE_FXC
530 
531   if(order < 2) return;
532 
533 
534   t102 = t12 * t63;
535   t106 = 0.1e1 / t4 / t27;
536   t110 = t27 * t44;
537   t112 = 0.1e1 / t3 / t110;
538   t117 = t49 * sigma[0];
539   t118 = params->b * t117;
540   t119 = t27 * t27;
541   t120 = t119 * t2;
542   t121 = 0.1e1 / t120;
543   t122 = params->k * params->k;
544   t127 = 0.88e2 / 0.9e1 * t1 * t106 * t9 - 0.24e2 * t50 * t112 * params->k * t9 + 0.64e2 / 0.9e1 * t118 * t121 * t122 * t9;
545   t128 = rho[0] * t127;
546   t133 = 0.1e1 / t62 / t42;
547   t134 = t76 * t76;
548   t135 = t133 * t134;
549   t148 = 0.1e1 / t81;
550   t151 = t33 * t32;
551   t152 = t23 * t151;
552   t153 = t152 * sigma[0];
553   t156 = t16 * t18 / t3 / t2 * t38 / 0.9e1 + 0.14e2 * t15 * t46 * params->d * t73 + 0.140e3 / 0.9e1 * t15 * t148 * params->d * t153;
554   t157 = t63 * t156;
555   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
556     v2rho2[0] = -0.2e1 * t102 * t76 + t128 * t43 + 0.2e1 * t61 * t135 - t61 * t157 + 0.2e1 * t58 * t43 - 0.2e1 * t59 * t77;
557 
558   t163 = params->b * t53;
559   t164 = t7 * t9;
560   t167 = t119 * rho[0];
561   t168 = 0.1e1 / t167;
562   t173 = -0.8e1 / 0.3e1 * params->b * t46 * t9 + 0.8e1 * t163 * t164 - 0.8e1 / 0.3e1 * t50 * t168 * t122 * t9;
563   t174 = rho[0] * t173;
564   t178 = 0.1e1 / t4 / rho[0];
565   t179 = t178 * t12;
566   t180 = t179 * t92;
567   t183 = t90 * t58;
568   t184 = t183 * t92;
569   t187 = t133 * t15;
570   t188 = t91 * t187;
571   t190 = t94 * t96 * t76;
572   t193 = 0.1e1 / t44;
573   t194 = t193 * t12;
574   t196 = t15 * params->d;
575   t197 = t196 * t152;
576   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
577     v2rhosigma[0] = t87 * t43 + t174 * t43 - t88 * t77 + 0.7e1 / 0.6e1 * t180 * t97 - 0.7e1 / 0.4e1 * t184 * t97 + 0.7e1 / 0.2e1 * t188 * t190 + 0.35e2 / 0.6e1 * t194 * t63 * t197;
578 
579   t201 = params->k * t9;
580   t204 = 0.1e1 / t119;
581   t208 = t1 * t204 * t122 * t9 - 0.2e1 * params->b * t83 * t201;
582   t209 = rho[0] * t208;
583   t211 = t90 * t87;
584   t212 = t211 * t92;
585   t215 = t168 * t12;
586   t216 = t215 * t133;
587   t217 = t15 * t15;
588   t218 = params->d * params->d;
589   t219 = t217 * t218;
590   t220 = t23 * t23;
591   t222 = t219 * t220 * t26;
592   t225 = 0.1e1 / t2;
593   t226 = t225 * t12;
594   t227 = t226 * t92;
595   t228 = 0.1e1 / sigma[0];
596   t229 = t151 * t228;
597   t230 = t94 * t229;
598   t233 = 0.1e1 / t26;
599   t234 = t71 * t233;
600   t235 = t94 * t234;
601   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
602     v2sigma2[0] = t209 * t43 - 0.7e1 / 0.2e1 * t212 * t97 + 0.49e2 / 0.8e1 * t216 * t222 - 0.35e2 / 0.16e2 * t227 * t230 + 0.7e1 / 0.8e1 * t93 * t235;
603 
604 #ifndef XC_DONT_COMPILE_KXC
605 
606   if(order < 3) return;
607 
608 
609   t240 = t58 * t63;
610   t243 = t12 * t133;
611   t249 = 0.1e1 / t4 / t81;
612   t254 = 0.1e1 / t3 / t119;
613   t259 = t119 * t44;
614   t260 = 0.1e1 / t259;
615   t265 = t49 * t49;
616   t266 = params->b * t265;
617   t267 = t119 * t81;
618   t269 = 0.1e1 / t4 / t267;
619   t270 = t122 * params->k;
620   t275 = -0.1232e4 / 0.27e2 * t1 * t249 * t9 + 0.5456e4 / 0.27e2 * t50 * t254 * params->k * t9 - 0.1216e4 / 0.9e1 * t118 * t260 * t122 * t9 + 0.512e3 / 0.27e2 * t266 * t269 * t270 * t9;
621   t276 = rho[0] * t275;
622   t284 = t62 * t62;
623   t285 = 0.1e1 / t284;
624   t286 = t134 * t76;
625   t287 = t285 * t286;
626   t290 = t133 * t76;
627   t291 = t290 * t156;
628   t295 = 0.1e1 / t3 / t44;
629   t312 = t23 * t33 * t26;
630   t315 = -0.7e1 / 0.27e2 * t16 * t18 * t295 * t38 - 0.1442e4 / 0.27e2 * t15 * t106 * params->d * t73 - 0.1120e4 / 0.9e1 * t15 / t51 * params->d * t153 - 0.280e3 / 0.9e1 * t15 * t112 * params->d * t312;
631   t316 = t63 * t315;
632   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
633     v3rho3[0] = -0.3e1 * t102 * t156 + 0.3e1 * t127 * t43 - 0.3e1 * t128 * t77 + 0.6e1 * t243 * t134 + 0.6e1 * t59 * t135 - 0.3e1 * t59 * t157 - 0.6e1 * t240 * t76 + t276 * t43 - 0.6e1 * t61 * t287 + 0.6e1 * t61 * t291 - t61 * t316;
634 
635   t318 = t6 * t12;
636   t319 = t318 * t92;
637   t322 = t179 * t187;
638   t325 = t285 * t15;
639   t326 = t91 * t325;
640   t328 = t94 * t96 * t134;
641   t333 = t194 * t187;
642   t335 = t94 * t151 * t76;
643   t343 = t87 * t63;
644   t349 = params->b * t112;
645   t352 = params->b * t121;
646   t354 = t122 * t49 * t9;
647   t357 = t119 * t27;
648   t359 = 0.1e1 / t4 / t357;
649   t364 = 0.88e2 / 0.9e1 * params->b * t106 * t9 - 0.520e3 / 0.9e1 * t349 * t164 + 0.136e3 / 0.3e1 * t352 * t354 - 0.64e2 / 0.9e1 * t118 * t359 * t270 * t9;
650   t365 = rho[0] * t364;
651   t367 = t90 * t127;
652   t368 = t367 * t92;
653   t371 = t178 * t58;
654   t372 = t371 * t92;
655   t375 = t83 * t12;
656   t376 = t375 * t92;
657   t378 = t94 * t33 * t25;
658   t381 = t28 * t12;
659   t385 = t193 * t58;
660   t389 = t183 * t187;
661   t393 = t94 * t96 * t156;
662   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
663     v3rho2sigma[0] = -0.35e2 / 0.18e2 * t319 * t97 - 0.14e2 / 0.3e1 * t322 * t190 - 0.21e2 / 0.2e1 * t326 * t328 + 0.2e1 * t173 * t43 - 0.70e2 / 0.3e1 * t333 * t335 + 0.2e1 * t88 * t135 - 0.2e1 * t174 * t77 - t88 * t157 - 0.2e1 * t343 * t76 + t365 * t43 - 0.7e1 / 0.4e1 * t368 * t97 + 0.7e1 / 0.3e1 * t372 * t97 - 0.35e2 / 0.3e1 * t376 * t378 - 0.385e3 / 0.18e2 * t381 * t63 * t197 + 0.35e2 / 0.3e1 * t385 * t63 * t197 + 0.7e1 * t389 * t190 + 0.7e1 / 0.2e1 * t188 * t393;
664 
665   t399 = params->b * t168;
666   t401 = t122 * sigma[0] * t9;
667   t405 = 0.1e1 / t4 / t259;
668   t410 = 0.32e2 / 0.3e1 * t163 * t201 - 0.40e2 / 0.3e1 * t399 * t401 + 0.8e1 / 0.3e1 * t50 * t405 * t270 * t9;
669   t411 = rho[0] * t410;
670   t414 = t178 * t87;
671   t415 = t414 * t92;
672   t418 = t90 * t173;
673   t419 = t418 * t92;
674   t422 = t211 * t187;
675   t425 = t193 * t87;
676   t429 = t121 * t12;
677   t430 = t429 * t133;
678   t433 = t168 * t58;
679   t434 = t433 * t133;
680   t437 = t285 * t217;
681   t438 = t215 * t437;
682   t439 = t218 * t220;
683   t441 = t439 * t26 * t76;
684   t444 = t194 * t92;
685   t447 = t225 * t58;
686   t448 = t447 * t92;
687   t451 = t226 * t187;
688   t453 = t94 * t229 * t76;
689   t457 = 0.1e1 / t3 / t27;
690   t458 = t457 * t12;
691   t459 = t458 * t92;
692   t460 = t33 * t95;
693   t461 = t94 * t460;
694   t469 = t94 * t234 * t76;
695   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
696     v3rhosigma2[0] = t208 * t43 + t411 * t43 - t209 * t77 + 0.7e1 / 0.3e1 * t415 * t97 - 0.7e1 / 0.2e1 * t419 * t97 + 0.7e1 * t422 * t190 + 0.35e2 / 0.3e1 * t425 * t63 * t197 - 0.441e3 / 0.8e1 * t430 * t222 + 0.49e2 / 0.8e1 * t434 * t222 - 0.147e3 / 0.8e1 * t438 * t441 + 0.35e2 / 0.24e2 * t444 * t230 - 0.35e2 / 0.16e2 * t448 * t230 + 0.35e2 / 0.8e1 * t451 * t453 + 0.35e2 / 0.8e1 * t459 * t461 - 0.7e1 / 0.12e2 * t180 * t235 + 0.7e1 / 0.8e1 * t184 * t235 - 0.7e1 / 0.4e1 * t188 * t469;
697 
698   t473 = t122 * t9;
699   t477 = 0.1e1 / t4 / t120;
700   t481 = -t1 * t477 * t270 * t9 + 0.3e1 * params->b * t204 * t473;
701   t482 = rho[0] * t481;
702   t484 = t90 * t208;
703   t485 = t484 * t92;
704   t488 = t168 * t87;
705   t489 = t488 * t133;
706   t492 = t225 * t87;
707   t493 = t492 * t92;
708   t498 = t477 * t12;
709   t500 = t218 * params->d;
710   t501 = t500 * sigma[0];
711   t502 = t501 * t71;
712   t506 = t219 * t220 * t25;
713   t509 = t295 * t12;
714   t510 = t509 * t92;
715   t511 = t33 * t233;
716   t512 = t94 * t511;
717   t515 = 0.1e1 / t49;
718   t516 = t151 * t515;
719   t517 = t94 * t516;
720   t521 = 0.1e1 / t25 / t49;
721   t522 = t71 * t521;
722   t523 = t94 * t522;
723   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
724     v3sigma3[0] = t482 * t43 - 0.21e2 / 0.4e1 * t485 * t97 + 0.147e3 / 0.8e1 * t489 * t222 - 0.105e3 / 0.16e2 * t493 * t230 + 0.21e2 / 0.8e1 * t212 * t235 - 0.1029e4 / 0.32e2 * t498 * t285 * t502 + 0.441e3 / 0.32e2 * t216 * t506 - 0.105e3 / 0.64e2 * t510 * t512 + 0.105e3 / 0.32e2 * t227 * t517 - 0.21e2 / 0.16e2 * t93 * t523;
725 
726 #ifndef XC_DONT_COMPILE_LXC
727 
728   if(order < 4) return;
729 
730 
731   t562 = 0.1e1 / t33;
732   t563 = t23 * t562;
733   t578 = 0.1e1 / t4 / t51;
734   t588 = 0.1e1 / t357;
735   t593 = t119 * t51;
736   t602 = t119 * t119;
737   t606 = t122 * t122;
738   t623 = 0.1e1 / t284 / t42;
739   t624 = t134 * t134;
740   t632 = t156 * t156;
741   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
742     v4rho4[0] = 0.24e2 * t58 * t133 * t134 - 0.24e2 * t12 * t285 * t286 + 0.24e2 * t243 * t76 * t156 - 0.4e1 * t276 * t77 - 0.6e1 * t128 * t157 - 0.4e1 * t59 * t316 - t61 * t63 * (0.70e2 / 0.81e2 * t16 * t18 * t457 * t38 + 0.6860e4 / 0.27e2 * t15 * t249 * params->d * t73 + 0.74900e5 / 0.81e2 * t15 / t110 * params->d * t153 + 0.12880e5 / 0.27e2 * t15 * t254 * params->d * t312 + 0.560e3 / 0.27e2 * t15 / t4 / t167 * params->d * t563 * t49) - 0.12e2 * t127 * t63 * t76 - 0.12e2 * t240 * t156 - 0.4e1 * t102 * t315 + rho[0] * (0.20944e5 / 0.81e2 * t1 * t578 * t9 - 0.48752e5 / 0.27e2 * t50 / t3 / t167 * params->k * t9 + 0.164032e6 / 0.81e2 * t118 * t588 * t122 * t9 - 0.50176e5 / 0.81e2 * t266 / t4 / t593 * t270 * t9 + 0.4096e4 / 0.81e2 * params->b * t265 * sigma[0] / t3 / t602 / rho[0] * t606 * t9) * t43 + 0.4e1 * t275 * t43 + 0.12e2 * t128 * t135 - 0.24e2 * t59 * t287 + 0.24e2 * t59 * t291 + 0.24e2 * t61 * t623 * t624 - 0.36e2 * t61 * t285 * t134 * t156 + 0.6e1 * t61 * t133 * t632 + 0.8e1 * t61 * t290 * t315;
743 
744   t710 = 0.6e1 * t88 * t291 - 0.6e1 * t173 * t63 * t76 + rho[0] * (-0.1232e4 / 0.27e2 * params->b * t249 * t9 + 0.4048e4 / 0.9e1 * params->b * t254 * t164 - 0.16400e5 / 0.27e2 * params->b * t260 * t354 + 0.5696e4 / 0.27e2 * params->b * t269 * t270 * t117 * t9 - 0.512e3 / 0.27e2 * t266 / t3 / t602 * t606 * t9) * t43 + 0.3e1 * t364 * t43 + 0.6e1 * t174 * t135 - 0.3e1 * t365 * t77 - 0.3e1 * t174 * t157 - t88 * t316 - 0.6e1 * t88 * t287 + 0.2485e4 / 0.27e2 * t148 * t12 * t63 * t197 + 0.35e2 / 0.2e1 * t193 * t127 * t63 * t197 - 0.385e3 / 0.6e1 * t28 * t58 * t63 * t197 + 0.35e2 / 0.3e1 * t318 * t187 * t190 - 0.14e2 * t371 * t187 * t190 - 0.7e1 * t322 * t393 - 0.63e2 / 0.2e1 * t183 * t325 * t328 + 0.70e2 * t375 * t187 * t94 * t33 * t76 * t25 + 0.21e2 / 0.2e1 * t367 * t187 * t190;
745   t785 = 0.21e2 / 0.2e1 * t389 * t393 + 0.7e1 / 0.2e1 * t188 * t94 * t96 * t315 + 0.21e2 * t179 * t325 * t328 + 0.42e2 * t91 * t623 * t15 * t94 * t96 * t286 - 0.63e2 / 0.2e1 * t91 * t325 * params->d * t72 * t95 * t76 * t156 + 0.6e1 * t87 * t133 * t134 - 0.3e1 * t343 * t156 - 0.35e2 / 0.6e1 * t6 * t58 * t92 * t97 + 0.385e3 / 0.3e1 * t381 * t187 * t335 + 0.105e3 * t194 * t325 * t94 * t151 * t134 - 0.70e2 * t385 * t187 * t335 - 0.35e2 * t333 * t94 * t151 * t156 - 0.7e1 / 0.4e1 * t90 * t275 * t92 * t97 - 0.35e2 * t83 * t58 * t92 * t378 + 0.70e2 / 0.9e1 / t4 / t110 * t12 * t92 * t94 * t562 * sigma[0] + 0.140e3 / 0.27e2 * t46 * t12 * t92 * t97 + 0.7e1 / 0.2e1 * t178 * t127 * t92 * t97 + 0.105e3 * t53 * t12 * t92 * t378;
746   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
747     v4rho3sigma[0] = t710 + t785;
748 
749   t869 = -0.35e2 / 0.12e2 * t578 * t12 * t63 * t196 * t563 - 0.28e2 / 0.3e1 * t414 * t187 * t190 + 0.14e2 * t418 * t187 * t190 + 0.7e1 * t422 * t393 - 0.35e2 / 0.6e1 * t333 * t453 + 0.35e2 / 0.4e1 * t447 * t187 * t453 + 0.35e2 / 0.8e1 * t451 * t94 * t229 * t156 - 0.35e2 / 0.2e1 * t458 * t187 * t94 * t460 * t76 + 0.7e1 / 0.3e1 * t322 * t469 - 0.7e1 / 0.2e1 * t389 * t469 - 0.7e1 / 0.4e1 * t188 * t94 * t234 * t156;
750   t943 = -0.147e3 / 0.8e1 * t438 * t439 * t26 * t156 + 0.35e2 / 0.12e2 * t385 * t92 * t230 - 0.175e3 / 0.8e1 * t376 * t461 - 0.35e2 / 0.16e2 * t225 * t127 * t92 * t230 + 0.35e2 / 0.4e1 * t457 * t58 * t92 * t461 - 0.7e1 / 0.6e1 * t372 * t235 + 0.7e1 / 0.8e1 * t368 * t235 - 0.35e2 / 0.9e1 * t6 * t87 * t92 * t97 + 0.147e3 / 0.2e1 * t215 * t623 * t217 * t439 * t26 * t134 - 0.175e3 / 0.72e2 * t381 * t92 * t230 + 0.35e2 / 0.36e2 * t319 * t235;
751   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
752     v4rho2sigma2[0] = -0.2e1 * t208 * t63 * t76 + rho[0] * (-0.608e3 / 0.9e1 * t349 * t201 + 0.1336e4 / 0.9e1 * t352 * t401 - 0.200e3 / 0.3e1 * params->b * t359 * t270 * t49 * t9 + 0.64e2 / 0.9e1 * t118 / t3 / t119 / t110 * t606 * t9) * t43 - 0.2e1 * t411 * t77 - t209 * t157 + 0.2e1 * t209 * t135 + 0.2205e4 / 0.4e1 * t260 * t12 * t133 * t222 - 0.385e3 / 0.9e1 * t28 * t87 * t63 * t197 + 0.70e2 / 0.3e1 * t193 * t173 * t63 * t197 - 0.441e3 / 0.4e1 * t121 * t58 * t133 * t222 + 0.49e2 / 0.8e1 * t168 * t127 * t133 * t222 + t869 - 0.21e2 * t211 * t325 * t328 - 0.105e3 / 0.8e1 * t226 * t325 * t94 * t229 * t134 + 0.21e2 / 0.4e1 * t326 * t94 * t234 * t134 + 0.2e1 * t410 * t43 + 0.14e2 / 0.3e1 * t178 * t173 * t92 * t97 - 0.7e1 / 0.2e1 * t90 * t364 * t92 * t97 - 0.140e3 / 0.3e1 * t425 * t187 * t335 - 0.70e2 / 0.3e1 * t83 * t87 * t92 * t378 + 0.1323e4 / 0.4e1 * t429 * t437 * t441 - 0.147e3 / 0.4e1 * t433 * t437 * t441 + t943;
753 
754   t1016 = t481 * t43 - 0.1029e4 / 0.32e2 * t477 * t58 * t285 * t502 + 0.1715e4 / 0.16e2 / t267 * t12 * t285 * t500 * t26 * t151 + 0.343e3 * t405 * t12 * t285 * t502 + rho[0] * (-0.24e2 * t399 * t473 + 0.56e2 / 0.3e1 * params->b * t405 * t270 * sigma[0] * t9 - 0.8e1 / 0.3e1 * t50 / t3 / t593 * t606 * t9) * t43 - t482 * t77 - 0.1323e4 / 0.8e1 * t121 * t87 * t133 * t222 - 0.3969e4 / 0.32e2 * t430 * t506 + 0.35e2 / 0.2e1 * t193 * t208 * t63 * t197 + 0.147e3 / 0.8e1 * t168 * t173 * t133 * t222 + 0.1029e4 / 0.8e1 * t498 * t623 * t501 * t71 * t76 + 0.441e3 / 0.32e2 * t434 * t506 + 0.21e2 / 0.2e1 * t484 * t187 * t190 + 0.105e3 / 0.8e1 * t492 * t187 * t453 - 0.21e2 / 0.4e1 * t422 * t469 + 0.105e3 / 0.32e2 * t509 * t187 * t94 * t511 * t76 - 0.105e3 / 0.16e2 * t451 * t94 * t516 * t76;
755   t1071 = 0.21e2 / 0.8e1 * t188 * t94 * t522 * t76 - 0.21e2 / 0.4e1 * t90 * t410 * t92 * t97 - 0.441e3 / 0.8e1 * t488 * t437 * t441 - 0.105e3 / 0.16e2 * t225 * t173 * t92 * t230 + 0.105e3 / 0.8e1 * t457 * t87 * t92 * t461 + 0.21e2 / 0.8e1 * t419 * t235 - 0.1323e4 / 0.32e2 * t438 * t439 * t25 * t76 - 0.105e3 / 0.64e2 * t295 * t58 * t92 * t512 + 0.35e2 / 0.32e2 * t249 * t12 * t92 * t94 * t562 * t228 + 0.105e3 / 0.32e2 * t448 * t517 - 0.21e2 / 0.16e2 * t184 * t523 + 0.7e1 / 0.2e1 * t178 * t208 * t92 * t97 + 0.35e2 / 0.8e1 * t425 * t92 * t230 - 0.7e1 / 0.4e1 * t415 * t235 - 0.35e2 / 0.32e2 * t459 * t512 - 0.35e2 / 0.16e2 * t444 * t517 + 0.7e1 / 0.8e1 * t180 * t523;
756   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
757     v4rhosigma3[0] = t1016 + t1071;
758 
759   t1118 = t218 * t218;
760   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
761     v4sigma4[0] = -0.7e1 * t90 * t481 * t92 * t97 - 0.105e3 / 0.8e1 * t225 * t208 * t92 * t230 + 0.441e3 / 0.8e1 * t489 * t506 - 0.105e3 / 0.16e2 * t295 * t87 * t92 * t512 + 0.735e3 / 0.128e3 * t216 * t219 * t220 * t95 - 0.105e3 / 0.256e3 * t106 * t12 * t92 * t94 * t562 * t515 - 0.1029e4 / 0.8e1 * t477 * t87 * t285 * t502 - 0.13377e5 / 0.128e3 * t498 * t285 * t500 * t71 - 0.5145e4 / 0.128e3 * t588 * t12 * t285 * t500 * t25 * t151 + 0.147e3 / 0.4e1 * t168 * t208 * t133 * t222 + 0.7203e4 / 0.32e2 / t602 / t44 * t12 * t623 * t1118 * t117 * t15 * t23 + 0.21e2 / 0.4e1 * t485 * t235 + 0.105e3 / 0.8e1 * t493 * t517 - 0.21e2 / 0.4e1 * t212 * t523 + 0.315e3 / 0.64e2 * t510 * t94 * t33 * t521 - 0.525e3 / 0.64e2 * t227 * t94 * t151 / t117 + 0.105e3 / 0.32e2 * t93 * t94 * t71 / t25 / t117 + rho[0] * (-0.4e1 * params->b * t477 * t270 * t9 + t1 / t3 / t267 * t606 * t9) * t43;
762 
763 #ifndef XC_DONT_COMPILE_MXC
764 
765   if(order < 5) return;
766 
767 
768 #endif
769 
770 #endif
771 
772 #endif
773 
774 #endif
775 
776 #endif
777 
778 
779 }
780 
781 
782 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 *))783 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
784 {
785 
786 #ifndef XC_DONT_COMPILE_EXC
787   double t2, t3, t4, t5, t6, t7, t9, t10;
788   double t12, t15, t16, t18, t19, t20, t21, t25;
789   double t26, t28, t29, t30, t31, t34, t35, t36;
790   double t41, t45, t46;
791 
792 #ifndef XC_DONT_COMPILE_VXC
793   double t47, t49, t52, t53, t54, t56, t61, t62;
794   double t64, t65, t66, t74, t75, t76, t79, t80;
795   double t84, t86, t90, t91, t93, t94, t95, t96;
796   double t97, t98, t99, t100, t101, t103, t104;
797 
798 #ifndef XC_DONT_COMPILE_FXC
799   double t109, t113, t117, t119, t124, t125, t126, t127;
800   double t128, t129, t134, t135, t140, t141, t142, t155;
801   double t158, t159, t160, t163, t164, t168, t170, t171;
802   double t172, t174, t175, t178, t180, t181, t185, t186;
803   double t187, t188, t190, t191, t192, t194, t195, t197;
804   double t198, t200, t201, t203, t204, t205, t211, t212;
805   double t220, t221, t223, t226, t227, t228, t230, t231;
806   double t233, t234, t235, t236, t237, t238, t239, t241;
807   double t242, t244, t245, t246, t247, t248, t249, t250;
808   double t252, t253, t254, t255, t259, t260, t262, t263;
809   double t264, t271, t272;
810 
811 #ifndef XC_DONT_COMPILE_KXC
812   double t280, t283, t289, t294, t299, t300, t305, t306;
813   double t307, t309, t310, t315, t316, t324, t325, t326;
814   double t327, t330, t331, t335, t352, t355, t356, t358;
815   double t359, t361, t362, t369, t370, t371, t373, t374;
816   double t375, t377, t379, t380, t382, t383, t385, t386;
817   double t388, t390, t392, t394, t396, t397, t398, t400;
818   double t401, t403, t404, t407, t408, t412, t416, t418;
819   double t419, t421, t423, t424, t426, t428, t431, t433;
820   double t434, t454, t461, t462, t465, t467, t469, t470;
821   double t473, t476, t478, t479, t482, t483, t485, t486;
822   double t487, t489, t490, t492, t493, t496, t497, t498;
823   double t499, t501, t502, t503, t505, t506, t507, t509;
824   double t510, t512, t513, t515, t516, t517, t519, t521;
825   double t522, t525, t526, t527, t528, t529, t530, t532;
826   double t534, t537, t538, t544, t545, t548, t549, t550;
827   double t552, t553, t554, t556, t557, t559, t561, t574;
828   double t579, t580, t598, t599, t602, t605, t606, t607;
829   double t609, t610, t611, t613, t614, t615, t617, t618;
830   double t619, t621, t623, t625, t626, t627, t628, t631;
831   double t632, t634, t635, t636, t637, t638, t640, t641;
832   double t642, t643, t646, t647, t648, t649, t653, t654;
833   double t656, t657, t658, t660, t661, t662, t664, t665;
834   double t666, t668, t671, t672, t673, t681, t682, t684;
835   double t685, t686, t699, t700;
836 
837 #ifndef XC_DONT_COMPILE_LXC
838   double t747, t748, t765, t775, t780, t789, t793, t808;
839   double t809, t817, t825, t828, t831, t836, t842, t862;
840   double t865, t868, t874, t877, t879, t881, t885, t888;
841   double t894, t898, t900, t903, t907, t919, t923, t927;
842   double t930, t934, t938, t946, t950, t953, t957, t961;
843   double t968, t970, t1007, t1026, t1032, t1036, t1039, t1043;
844   double t1045, t1048, t1052, t1055, t1057, t1062, t1066, t1068;
845   double t1070, t1075, t1080, t1083, t1085, t1092, t1096, t1100;
846   double t1103, t1105, t1106, t1109, t1113, t1118, t1120, t1122;
847   double t1127, t1131, t1135, t1140, t1142, t1147, t1154, t1161;
848   double t1166, t1170, t1173, t1177, t1182, t1184, t1188, t1192;
849   double t1195, t1199, t1203, t1209, t1225, t1228, t1231, t1233;
850   double t1236, t1244, t1264, t1268, t1270, t1296, t1329, t1331;
851   double t1336, t1342, t1349, t1351, t1355, t1359, t1364, t1366;
852   double t1371, t1374, t1376, t1378, t1380, t1382, t1386, t1389;
853   double t1393, t1395, t1398, t1400, t1404, t1408, t1414, t1416;
854   double t1418, t1423, t1430, t1434, t1437, t1440, t1442, t1447;
855   double t1451, t1455, t1457, t1464, t1468, t1472, t1476, t1479;
856   double t1484, t1485, t1486, t1491, t1492, t1493, t1494, t1497;
857   double t1504, t1508, t1511, t1515, t1519, t1521, t1524, t1526;
858   double t1529, t1534, t1535, t1540, t1543, t1545, t1547, t1556;
859   double t1560, t1569, t1576, t1590, t1594, t1596, t1607, t1631;
860   double t1649, t1652, t1656, t1662, t1666, t1667, t1672, t1676;
861   double t1680, t1686, t1692, t1696, t1697, t1698, t1702, t1706;
862   double t1712, t1714, t1715, t1716, t1718, t1722, t1727, t1733;
863   double t1737, t1747, t1749, t1757, t1767, t1771, t1775, t1777;
864   double t1779, t1781, t1783, t1785, t1790, t1796, t1804, t1806;
865   double t1808, t1809, t1821, t1825, t1827, t1829, t1844, t1858;
866   double t1860;
867 #endif
868 
869 #endif
870 
871 #endif
872 
873 #endif
874 
875 #endif
876 
877 
878   gga_c_wi_params *params;
879 
880   assert(p->params != NULL);
881   params = (gga_c_wi_params * )(p->params);
882 
883   t2 = sigma[0] + 0.2e1 * sigma[1] + sigma[2];
884   t3 = params->b * t2;
885   t4 = rho[0] + rho[1];
886   t5 = t4 * t4;
887   t6 = POW_1_3(t4);
888   t7 = t6 * t6;
889   t9 = 0.1e1 / t7 / t5;
890   t10 = params->k * t2;
891   t12 = exp(-t10 * t9);
892   t15 = t3 * t9 * t12 + params->a;
893   t16 = M_CBRT3;
894   t18 = POW_1_3(0.1e1 / M_PI);
895   t19 = t16 * t18;
896   t20 = M_CBRT4;
897   t21 = t20 * t20;
898   t25 = t16 * t16;
899   t26 = M_CBRTPI;
900   t28 = sqrt(t2);
901   t29 = t28 * t2;
902   t30 = t5 * t5;
903   t31 = 0.1e1 / t30;
904   t34 = 0.1e1 / t6 / t4;
905   t35 = t28 * t34;
906   t36 = sqrt(t35);
907   t41 = 0.1e1 + params->d * t20 * t25 * t26 * t36 * t29 * t31 / 0.3e1;
908   t45 = params->c + t19 * t21 / t6 * t41 / 0.4e1;
909   t46 = 0.1e1 / t45;
910   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
911     zk[0] = t15 * t46;
912 
913 #ifndef XC_DONT_COMPILE_VXC
914 
915   if(order < 1) return;
916 
917 
918   t47 = t5 * t4;
919   t49 = 0.1e1 / t7 / t47;
920   t52 = t2 * t2;
921   t53 = params->b * t52;
922   t54 = t30 * t5;
923   t56 = 0.1e1 / t6 / t54;
924   t61 = 0.8e1 / 0.3e1 * t53 * t56 * params->k * t12 - 0.8e1 / 0.3e1 * t3 * t49 * t12;
925   t62 = t4 * t61;
926   t64 = t4 * t15;
927   t65 = t45 * t45;
928   t66 = 0.1e1 / t65;
929   t74 = t36 * t2 * t9;
930   t75 = t26 * t74;
931   t76 = t75 * t28;
932   t79 = -t19 * t21 * t34 * t41 / 0.12e2 - 0.14e2 / 0.3e1 * t18 * t9 * params->d * t76;
933   t80 = t66 * t79;
934   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
935     vrho[0] = t62 * t46 - t64 * t80 + (t15 * t46);
936 
937   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
938     vrho[1] = vrho[0];
939 
940   t84 = t30 * t4;
941   t86 = 0.1e1 / t6 / t84;
942   t90 = -t3 * t86 * params->k * t12 + params->b * t9 * t12;
943   t91 = t4 * t90;
944   t93 = 0.1e1 / t7;
945   t94 = t93 * t15;
946   t95 = t66 * t18;
947   t96 = t94 * t95;
948   t97 = params->d * t26;
949   t98 = 0.1e1 / t28;
950   t99 = t74 * t98;
951   t100 = t97 * t99;
952   t101 = t96 * t100;
953   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
954     vsigma[0] = t91 * t46 - 0.7e1 / 0.4e1 * t101;
955 
956   t103 = 0.2e1 * t90;
957   t104 = t4 * t103;
958   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
959     vsigma[1] = t104 * t46 - 0.7e1 / 0.2e1 * t101;
960 
961   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
962     vsigma[2] = vsigma[0];
963 
964 #ifndef XC_DONT_COMPILE_FXC
965 
966   if(order < 2) return;
967 
968 
969   t109 = t15 * t66;
970   t113 = 0.1e1 / t7 / t30;
971   t117 = t30 * t47;
972   t119 = 0.1e1 / t6 / t117;
973   t124 = t52 * t2;
974   t125 = params->b * t124;
975   t126 = t30 * t30;
976   t127 = t126 * t5;
977   t128 = 0.1e1 / t127;
978   t129 = params->k * params->k;
979   t134 = 0.88e2 / 0.9e1 * t3 * t113 * t12 - 0.24e2 * t53 * t119 * params->k * t12 + 0.64e2 / 0.9e1 * t125 * t128 * t129 * t12;
980   t135 = t4 * t134;
981   t140 = 0.1e1 / t65 / t45;
982   t141 = t79 * t79;
983   t142 = t140 * t141;
984   t155 = 0.1e1 / t84;
985   t158 = t36 * t35;
986   t159 = t26 * t158;
987   t160 = t159 * t2;
988   t163 = t19 * t21 / t6 / t5 * t41 / 0.9e1 + 0.14e2 * t18 * t49 * params->d * t76 + 0.140e3 / 0.9e1 * t18 * t155 * params->d * t160;
989   t164 = t66 * t163;
990   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
991     v2rho2[0] = -0.2e1 * t109 * t79 + t135 * t46 + 0.2e1 * t64 * t142 - t64 * t164 + 0.2e1 * t61 * t46 - 0.2e1 * t62 * t80;
992 
993   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
994     v2rho2[1] = v2rho2[0];
995 
996   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
997     v2rho2[2] = v2rho2[1];
998 
999   t168 = params->b * t49 * t12;
1000   t170 = params->b * t56;
1001   t171 = t10 * t12;
1002   t172 = t170 * t171;
1003   t174 = t126 * t4;
1004   t175 = 0.1e1 / t174;
1005   t178 = t53 * t175 * t129 * t12;
1006   t180 = -0.8e1 / 0.3e1 * t168 + 0.8e1 * t172 - 0.8e1 / 0.3e1 * t178;
1007   t181 = t4 * t180;
1008   t185 = 0.1e1 / t7 / t4;
1009   t186 = t185 * t15;
1010   t187 = t186 * t95;
1011   t188 = t187 * t100;
1012   t190 = t93 * t61;
1013   t191 = t190 * t95;
1014   t192 = t191 * t100;
1015   t194 = t140 * t18;
1016   t195 = t94 * t194;
1017   t197 = t97 * t99 * t79;
1018   t198 = t195 * t197;
1019   t200 = 0.1e1 / t47;
1020   t201 = t200 * t15;
1021   t203 = t18 * params->d;
1022   t204 = t203 * t159;
1023   t205 = t201 * t66 * t204;
1024   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1025     v2rhosigma[0] = t90 * t46 + t181 * t46 - t91 * t80 + 0.7e1 / 0.6e1 * t188 - 0.7e1 / 0.4e1 * t192 + 0.7e1 / 0.2e1 * t198 + 0.35e2 / 0.6e1 * t205;
1026 
1027   t211 = -0.16e2 / 0.3e1 * t168 + 0.16e2 * t172 - 0.16e2 / 0.3e1 * t178;
1028   t212 = t4 * t211;
1029   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1030     v2rhosigma[1] = t103 * t46 + t212 * t46 - t104 * t80 + 0.7e1 / 0.3e1 * t188 - 0.7e1 / 0.2e1 * t192 + 0.7e1 * t198 + 0.35e2 / 0.3e1 * t205;
1031 
1032   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1033     v2rhosigma[2] = v2rhosigma[0];
1034 
1035   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1036     v2rhosigma[3] = v2rhosigma[2];
1037 
1038   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1039     v2rhosigma[4] = v2rhosigma[1];
1040 
1041   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1042     v2rhosigma[5] = v2rhosigma[3];
1043 
1044   t220 = params->k * t12;
1045   t221 = params->b * t86 * t220;
1046   t223 = 0.1e1 / t126;
1047   t226 = t3 * t223 * t129 * t12;
1048   t227 = -0.2e1 * t221 + t226;
1049   t228 = t4 * t227;
1050   t230 = t93 * t90;
1051   t231 = t230 * t95;
1052   t233 = 0.7e1 / 0.2e1 * t231 * t100;
1053   t234 = t175 * t15;
1054   t235 = t234 * t140;
1055   t236 = t18 * t18;
1056   t237 = params->d * params->d;
1057   t238 = t236 * t237;
1058   t239 = t26 * t26;
1059   t241 = t238 * t239 * t29;
1060   t242 = t235 * t241;
1061   t244 = 0.1e1 / t5;
1062   t245 = t244 * t15;
1063   t246 = t245 * t95;
1064   t247 = 0.1e1 / t2;
1065   t248 = t158 * t247;
1066   t249 = t97 * t248;
1067   t250 = t246 * t249;
1068   t252 = 0.1e1 / t29;
1069   t253 = t74 * t252;
1070   t254 = t97 * t253;
1071   t255 = t96 * t254;
1072   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1073     v2sigma2[0] = t228 * t46 - t233 + 0.49e2 / 0.8e1 * t242 - 0.35e2 / 0.16e2 * t250 + 0.7e1 / 0.8e1 * t255;
1074 
1075   t259 = -0.4e1 * t221 + 0.2e1 * t226;
1076   t260 = t4 * t259;
1077   t262 = t93 * t103;
1078   t263 = t262 * t95;
1079   t264 = t263 * t100;
1080   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1081     v2sigma2[1] = t260 * t46 - 0.7e1 / 0.4e1 * t264 - t233 + 0.49e2 / 0.4e1 * t242 - 0.35e2 / 0.8e1 * t250 + 0.7e1 / 0.4e1 * t255;
1082 
1083   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1084     v2sigma2[2] = v2sigma2[0];
1085 
1086   t271 = -0.8e1 * t221 + 0.4e1 * t226;
1087   t272 = t4 * t271;
1088   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1089     v2sigma2[3] = t272 * t46 - 0.7e1 * t264 + 0.49e2 / 0.2e1 * t242 - 0.35e2 / 0.4e1 * t250 + 0.7e1 / 0.2e1 * t255;
1090 
1091   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1092     v2sigma2[4] = v2sigma2[1];
1093 
1094   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1095     v2sigma2[5] = v2sigma2[2];
1096 
1097 #ifndef XC_DONT_COMPILE_KXC
1098 
1099   if(order < 3) return;
1100 
1101 
1102   t280 = t61 * t66;
1103   t283 = t15 * t140;
1104   t289 = 0.1e1 / t7 / t84;
1105   t294 = 0.1e1 / t6 / t126;
1106   t299 = t126 * t47;
1107   t300 = 0.1e1 / t299;
1108   t305 = t52 * t52;
1109   t306 = params->b * t305;
1110   t307 = t126 * t84;
1111   t309 = 0.1e1 / t7 / t307;
1112   t310 = t129 * params->k;
1113   t315 = -0.1232e4 / 0.27e2 * t3 * t289 * t12 + 0.5456e4 / 0.27e2 * t53 * t294 * params->k * t12 - 0.1216e4 / 0.9e1 * t125 * t300 * t129 * t12 + 0.512e3 / 0.27e2 * t306 * t309 * t310 * t12;
1114   t316 = t4 * t315;
1115   t324 = t65 * t65;
1116   t325 = 0.1e1 / t324;
1117   t326 = t141 * t79;
1118   t327 = t325 * t326;
1119   t330 = t140 * t79;
1120   t331 = t330 * t163;
1121   t335 = 0.1e1 / t6 / t47;
1122   t352 = t26 * t36 * t29;
1123   t355 = -0.7e1 / 0.27e2 * t19 * t21 * t335 * t41 - 0.1442e4 / 0.27e2 * t18 * t113 * params->d * t76 - 0.1120e4 / 0.9e1 * t18 / t54 * params->d * t160 - 0.280e3 / 0.9e1 * t18 * t119 * params->d * t352;
1124   t356 = t66 * t355;
1125   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1126     v3rho3[0] = -0.3e1 * t109 * t163 + 0.3e1 * t134 * t46 - 0.3e1 * t135 * t80 + 0.6e1 * t283 * t141 + 0.6e1 * t62 * t142 - 0.3e1 * t62 * t164 - 0.6e1 * t280 * t79 + t316 * t46 - 0.6e1 * t64 * t327 + 0.6e1 * t64 * t331 - t64 * t356;
1127 
1128   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1129     v3rho3[1] = v3rho3[0];
1130 
1131   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1132     v3rho3[2] = v3rho3[1];
1133 
1134   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1135     v3rho3[3] = v3rho3[2];
1136 
1137   t358 = t325 * t18;
1138   t359 = t94 * t358;
1139   t361 = t97 * t99 * t141;
1140   t362 = t359 * t361;
1141   t369 = t185 * t61;
1142   t370 = t369 * t95;
1143   t371 = t370 * t100;
1144   t373 = t93 * t134;
1145   t374 = t373 * t95;
1146   t375 = t374 * t100;
1147   t377 = t201 * t194;
1148   t379 = t97 * t158 * t79;
1149   t380 = t377 * t379;
1150   t382 = t86 * t15;
1151   t383 = t382 * t95;
1152   t385 = t97 * t36 * t28;
1153   t386 = t383 * t385;
1154   t388 = t31 * t15;
1155   t390 = t388 * t66 * t204;
1156   t392 = t200 * t61;
1157   t394 = t392 * t66 * t204;
1158   t396 = t9 * t15;
1159   t397 = t396 * t95;
1160   t398 = t397 * t100;
1161   t400 = t186 * t194;
1162   t401 = t400 * t197;
1163   t403 = t190 * t194;
1164   t404 = t403 * t197;
1165   t407 = t97 * t99 * t163;
1166   t408 = t195 * t407;
1167   t412 = t90 * t66;
1168   t416 = params->b * t113 * t12;
1169   t418 = params->b * t119;
1170   t419 = t418 * t171;
1171   t421 = params->b * t128;
1172   t423 = t129 * t52 * t12;
1173   t424 = t421 * t423;
1174   t426 = t126 * t30;
1175   t428 = 0.1e1 / t7 / t426;
1176   t431 = t125 * t428 * t310 * t12;
1177   t433 = 0.88e2 / 0.9e1 * t416 - 0.520e3 / 0.9e1 * t419 + 0.136e3 / 0.3e1 * t424 - 0.64e2 / 0.9e1 * t431;
1178   t434 = t4 * t433;
1179   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1180     v3rho2sigma[0] = -0.21e2 / 0.2e1 * t362 + 0.2e1 * t91 * t142 - 0.2e1 * t181 * t80 - t91 * t164 + 0.7e1 / 0.3e1 * t371 - 0.7e1 / 0.4e1 * t375 - 0.70e2 / 0.3e1 * t380 - 0.35e2 / 0.3e1 * t386 - 0.385e3 / 0.18e2 * t390 + 0.35e2 / 0.3e1 * t394 - 0.35e2 / 0.18e2 * t398 - 0.14e2 / 0.3e1 * t401 + 0.7e1 * t404 + 0.7e1 / 0.2e1 * t408 + 0.2e1 * t180 * t46 - 0.2e1 * t412 * t79 + t434 * t46;
1181 
1182   t454 = t103 * t66;
1183   t461 = 0.176e3 / 0.9e1 * t416 - 0.1040e4 / 0.9e1 * t419 + 0.272e3 / 0.3e1 * t424 - 0.128e3 / 0.9e1 * t431;
1184   t462 = t4 * t461;
1185   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1186     v3rho2sigma[1] = -0.21e2 * t362 + 0.2e1 * t104 * t142 - 0.2e1 * t212 * t80 - t104 * t164 + 0.14e2 / 0.3e1 * t371 - 0.7e1 / 0.2e1 * t375 - 0.140e3 / 0.3e1 * t380 - 0.70e2 / 0.3e1 * t386 - 0.385e3 / 0.9e1 * t390 + 0.70e2 / 0.3e1 * t394 - 0.35e2 / 0.9e1 * t398 - 0.28e2 / 0.3e1 * t401 + 0.14e2 * t404 + 0.7e1 * t408 + 0.2e1 * t211 * t46 - 0.2e1 * t454 * t79 + t462 * t46;
1187 
1188   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1189     v3rho2sigma[2] = v3rho2sigma[0];
1190 
1191   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1192     v3rho2sigma[3] = v3rho2sigma[2];
1193 
1194   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1195     v3rho2sigma[4] = v3rho2sigma[1];
1196 
1197   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1198     v3rho2sigma[5] = v3rho2sigma[3];
1199 
1200   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1201     v3rho2sigma[6] = v3rho2sigma[5];
1202 
1203   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1204     v3rho2sigma[7] = v3rho2sigma[4];
1205 
1206   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1207     v3rho2sigma[8] = v3rho2sigma[6];
1208 
1209   t465 = t170 * t220;
1210   t467 = params->b * t175;
1211   t469 = t129 * t2 * t12;
1212   t470 = t467 * t469;
1213   t473 = 0.1e1 / t7 / t299;
1214   t476 = t53 * t473 * t310 * t12;
1215   t478 = 0.32e2 / 0.3e1 * t465 - 0.40e2 / 0.3e1 * t470 + 0.8e1 / 0.3e1 * t476;
1216   t479 = t4 * t478;
1217   t482 = t185 * t90;
1218   t483 = t482 * t95;
1219   t485 = 0.7e1 / 0.3e1 * t483 * t100;
1220   t486 = t93 * t180;
1221   t487 = t486 * t95;
1222   t489 = 0.7e1 / 0.2e1 * t487 * t100;
1223   t490 = t230 * t194;
1224   t492 = 0.7e1 * t490 * t197;
1225   t493 = t200 * t90;
1226   t496 = 0.35e2 / 0.3e1 * t493 * t66 * t204;
1227   t497 = t128 * t15;
1228   t498 = t497 * t140;
1229   t499 = t498 * t241;
1230   t501 = t175 * t61;
1231   t502 = t501 * t140;
1232   t503 = t502 * t241;
1233   t505 = t325 * t236;
1234   t506 = t234 * t505;
1235   t507 = t237 * t239;
1236   t509 = t507 * t29 * t79;
1237   t510 = t506 * t509;
1238   t512 = t201 * t95;
1239   t513 = t512 * t249;
1240   t515 = t244 * t61;
1241   t516 = t515 * t95;
1242   t517 = t516 * t249;
1243   t519 = t245 * t194;
1244   t521 = t97 * t248 * t79;
1245   t522 = t519 * t521;
1246   t525 = 0.1e1 / t6 / t30;
1247   t526 = t525 * t15;
1248   t527 = t526 * t95;
1249   t528 = t36 * t98;
1250   t529 = t97 * t528;
1251   t530 = t527 * t529;
1252   t532 = t187 * t254;
1253   t534 = t191 * t254;
1254   t537 = t97 * t253 * t79;
1255   t538 = t195 * t537;
1256   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1257     v3rhosigma2[0] = t227 * t46 + t479 * t46 - t228 * t80 + t485 - t489 + t492 + t496 - 0.441e3 / 0.8e1 * t499 + 0.49e2 / 0.8e1 * t503 - 0.147e3 / 0.8e1 * t510 + 0.35e2 / 0.24e2 * t513 - 0.35e2 / 0.16e2 * t517 + 0.35e2 / 0.8e1 * t522 + 0.35e2 / 0.8e1 * t530 - 0.7e1 / 0.12e2 * t532 + 0.7e1 / 0.8e1 * t534 - 0.7e1 / 0.4e1 * t538;
1258 
1259   t544 = 0.64e2 / 0.3e1 * t465 - 0.80e2 / 0.3e1 * t470 + 0.16e2 / 0.3e1 * t476;
1260   t545 = t4 * t544;
1261   t548 = t185 * t103;
1262   t549 = t548 * t95;
1263   t550 = t549 * t100;
1264   t552 = t93 * t211;
1265   t553 = t552 * t95;
1266   t554 = t553 * t100;
1267   t556 = t262 * t194;
1268   t557 = t556 * t197;
1269   t559 = t200 * t103;
1270   t561 = t559 * t66 * t204;
1271   t574 = t496 - 0.441e3 / 0.4e1 * t499 + 0.49e2 / 0.4e1 * t503 - 0.147e3 / 0.4e1 * t510 + 0.35e2 / 0.12e2 * t513 - 0.35e2 / 0.8e1 * t517 + 0.35e2 / 0.4e1 * t522 + 0.35e2 / 0.4e1 * t530 - 0.7e1 / 0.6e1 * t532 + 0.7e1 / 0.4e1 * t534 - 0.7e1 / 0.2e1 * t538;
1272   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1273     v3rhosigma2[1] = t259 * t46 + t545 * t46 - t260 * t80 + 0.7e1 / 0.6e1 * t550 - 0.7e1 / 0.4e1 * t554 + 0.7e1 / 0.2e1 * t557 + 0.35e2 / 0.6e1 * t561 + t485 - t489 + t492 + t574;
1274 
1275   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1276     v3rhosigma2[2] = v3rhosigma2[0];
1277 
1278   t579 = 0.128e3 / 0.3e1 * t465 - 0.160e3 / 0.3e1 * t470 + 0.32e2 / 0.3e1 * t476;
1279   t580 = t4 * t579;
1280   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1281     v3rhosigma2[3] = t271 * t46 + t580 * t46 - t272 * t80 + 0.14e2 / 0.3e1 * t550 - 0.7e1 * t554 + 0.14e2 * t557 + 0.70e2 / 0.3e1 * t561 - 0.441e3 / 0.2e1 * t499 + 0.49e2 / 0.2e1 * t503 - 0.147e3 / 0.2e1 * t510 + 0.35e2 / 0.6e1 * t513 - 0.35e2 / 0.4e1 * t517 + 0.35e2 / 0.2e1 * t522 + 0.35e2 / 0.2e1 * t530 - 0.7e1 / 0.3e1 * t532 + 0.7e1 / 0.2e1 * t534 - 0.7e1 * t538;
1282 
1283   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1284     v3rhosigma2[4] = v3rhosigma2[1];
1285 
1286   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1287     v3rhosigma2[5] = v3rhosigma2[2];
1288 
1289   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1290     v3rhosigma2[6] = v3rhosigma2[5];
1291 
1292   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1293     v3rhosigma2[7] = v3rhosigma2[4];
1294 
1295   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1296     v3rhosigma2[8] = v3rhosigma2[6];
1297 
1298   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1299     v3rhosigma2[9] = v3rhosigma2[3];
1300 
1301   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1302     v3rhosigma2[10] = v3rhosigma2[7];
1303 
1304   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1305     v3rhosigma2[11] = v3rhosigma2[8];
1306 
1307   t598 = t129 * t12;
1308   t599 = params->b * t223 * t598;
1309   t602 = 0.1e1 / t7 / t127;
1310   t605 = t3 * t602 * t310 * t12;
1311   t606 = 0.3e1 * t599 - t605;
1312   t607 = t4 * t606;
1313   t609 = t93 * t227;
1314   t610 = t609 * t95;
1315   t611 = t610 * t100;
1316   t613 = t175 * t90;
1317   t614 = t613 * t140;
1318   t615 = t614 * t241;
1319   t617 = t244 * t90;
1320   t618 = t617 * t95;
1321   t619 = t618 * t249;
1322   t621 = t231 * t254;
1323   t623 = t602 * t15;
1324   t625 = t237 * params->d;
1325   t626 = t625 * t2;
1326   t627 = t626 * t74;
1327   t628 = t623 * t325 * t627;
1328   t631 = t238 * t239 * t28;
1329   t632 = t235 * t631;
1330   t634 = t335 * t15;
1331   t635 = t634 * t95;
1332   t636 = t36 * t252;
1333   t637 = t97 * t636;
1334   t638 = t635 * t637;
1335   t640 = 0.1e1 / t52;
1336   t641 = t158 * t640;
1337   t642 = t97 * t641;
1338   t643 = t246 * t642;
1339   t646 = 0.1e1 / t28 / t52;
1340   t647 = t74 * t646;
1341   t648 = t97 * t647;
1342   t649 = t96 * t648;
1343   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1344     v3sigma3[0] = t607 * t46 - 0.21e2 / 0.4e1 * t611 + 0.147e3 / 0.8e1 * t615 - 0.105e3 / 0.16e2 * t619 + 0.21e2 / 0.8e1 * t621 - 0.1029e4 / 0.32e2 * t628 + 0.441e3 / 0.32e2 * t632 - 0.105e3 / 0.64e2 * t638 + 0.105e3 / 0.32e2 * t643 - 0.21e2 / 0.16e2 * t649;
1345 
1346   t653 = 0.6e1 * t599 - 0.2e1 * t605;
1347   t654 = t4 * t653;
1348   t656 = t93 * t259;
1349   t657 = t656 * t95;
1350   t658 = t657 * t100;
1351   t660 = t175 * t103;
1352   t661 = t660 * t140;
1353   t662 = t661 * t241;
1354   t664 = t244 * t103;
1355   t665 = t664 * t95;
1356   t666 = t665 * t249;
1357   t668 = t263 * t254;
1358   t671 = 0.49e2 / 0.2e1 * t615;
1359   t672 = 0.35e2 / 0.4e1 * t619;
1360   t673 = 0.7e1 / 0.2e1 * t621;
1361   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1362     v3sigma3[1] = t654 * t46 - 0.7e1 / 0.2e1 * t658 + 0.49e2 / 0.8e1 * t662 - 0.35e2 / 0.16e2 * t666 + 0.7e1 / 0.8e1 * t668 - 0.7e1 / 0.2e1 * t611 + t671 - t672 + t673 - 0.1029e4 / 0.16e2 * t628 + 0.441e3 / 0.16e2 * t632 - 0.105e3 / 0.32e2 * t638 + 0.105e3 / 0.16e2 * t643 - 0.21e2 / 0.8e1 * t649;
1363 
1364   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1365     v3sigma3[2] = v3sigma3[0];
1366 
1367   t681 = 0.12e2 * t599 - 0.4e1 * t605;
1368   t682 = t4 * t681;
1369   t684 = t93 * t271;
1370   t685 = t684 * t95;
1371   t686 = t685 * t100;
1372   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1373     v3sigma3[3] = t682 * t46 - 0.7e1 / 0.4e1 * t686 - 0.7e1 * t658 + 0.49e2 / 0.2e1 * t662 - 0.35e2 / 0.4e1 * t666 + 0.7e1 / 0.2e1 * t668 + t671 - 0.1029e4 / 0.8e1 * t628 + 0.441e3 / 0.8e1 * t632 - t672 - 0.105e3 / 0.16e2 * t638 + 0.105e3 / 0.8e1 * t643 + t673 - 0.21e2 / 0.4e1 * t649;
1374 
1375   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1376     v3sigma3[4] = v3sigma3[1];
1377 
1378   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1379     v3sigma3[5] = v3sigma3[2];
1380 
1381   t699 = 0.24e2 * t599 - 0.8e1 * t605;
1382   t700 = t4 * t699;
1383   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1384     v3sigma3[6] = t700 * t46 - 0.21e2 / 0.2e1 * t686 + 0.147e3 / 0.2e1 * t662 - 0.105e3 / 0.4e1 * t666 + 0.21e2 / 0.2e1 * t668 - 0.1029e4 / 0.4e1 * t628 + 0.441e3 / 0.4e1 * t632 - 0.105e3 / 0.8e1 * t638 + 0.105e3 / 0.4e1 * t643 - 0.21e2 / 0.2e1 * t649;
1385 
1386   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1387     v3sigma3[7] = v3sigma3[3];
1388 
1389   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1390     v3sigma3[8] = v3sigma3[4];
1391 
1392   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1393     v3sigma3[9] = v3sigma3[5];
1394 
1395 #ifndef XC_DONT_COMPILE_LXC
1396 
1397   if(order < 4) return;
1398 
1399 
1400   t747 = 0.1e1 / t36;
1401   t748 = t26 * t747;
1402   t765 = 0.1e1 / t7 / t54;
1403   t775 = 0.1e1 / t426;
1404   t780 = t126 * t54;
1405   t789 = t126 * t126;
1406   t793 = t129 * t129;
1407   t808 = 0.1e1 / t324 / t45;
1408   t809 = t141 * t141;
1409   t817 = t163 * t163;
1410   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1411     v4rho4[0] = 0.24e2 * t61 * t140 * t141 - 0.24e2 * t15 * t325 * t326 + 0.24e2 * t283 * t79 * t163 - 0.4e1 * t316 * t80 - 0.6e1 * t135 * t164 - 0.4e1 * t62 * t356 - t64 * t66 * (0.70e2 / 0.81e2 * t19 * t21 * t525 * t41 + 0.6860e4 / 0.27e2 * t18 * t289 * params->d * t76 + 0.74900e5 / 0.81e2 * t18 / t117 * params->d * t160 + 0.12880e5 / 0.27e2 * t18 * t294 * params->d * t352 + 0.560e3 / 0.27e2 * t18 / t7 / t174 * params->d * t748 * t52) + 0.4e1 * t315 * t46 - 0.12e2 * t134 * t66 * t79 - 0.12e2 * t280 * t163 - 0.4e1 * t109 * t355 + t4 * (0.20944e5 / 0.81e2 * t3 * t765 * t12 - 0.48752e5 / 0.27e2 * t53 / t6 / t174 * params->k * t12 + 0.164032e6 / 0.81e2 * t125 * t775 * t129 * t12 - 0.50176e5 / 0.81e2 * t306 / t7 / t780 * t310 * t12 + 0.4096e4 / 0.81e2 * params->b * t305 * t2 / t6 / t789 / t4 * t793 * t12) * t46 + 0.12e2 * t135 * t142 - 0.24e2 * t62 * t327 + 0.24e2 * t62 * t331 + 0.24e2 * t64 * t808 * t809 - 0.36e2 * t64 * t325 * t141 * t163 + 0.6e1 * t64 * t140 * t817 + 0.8e1 * t64 * t330 * t355;
1412 
1413   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1414     v4rho4[1] = v4rho4[0];
1415 
1416   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1417     v4rho4[2] = v4rho4[1];
1418 
1419   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1420     v4rho4[3] = v4rho4[2];
1421 
1422   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1423     v4rho4[4] = v4rho4[3];
1424 
1425   t825 = params->b * t289 * t12;
1426   t828 = params->b * t294 * t171;
1427   t831 = params->b * t300 * t423;
1428   t836 = params->b * t309 * t310 * t124 * t12;
1429   t842 = t306 / t6 / t789 * t793 * t12;
1430   t862 = t190 * t358 * t361;
1431   t865 = t369 * t194 * t197;
1432   t868 = t373 * t194 * t197;
1433   t874 = t382 * t194 * t97 * t36 * t79 * t28;
1434   t877 = t396 * t194 * t197;
1435   t879 = t400 * t407;
1436   t881 = t403 * t407;
1437   t885 = t195 * t97 * t99 * t355;
1438   t888 = t186 * t358 * t361;
1439   t894 = t94 * t808 * t18 * t97 * t99 * t326;
1440   t898 = t155 * t15 * t66 * t204;
1441   t900 = t4 * (-0.1232e4 / 0.27e2 * t825 + 0.4048e4 / 0.9e1 * t828 - 0.16400e5 / 0.27e2 * t831 + 0.5696e4 / 0.27e2 * t836 - 0.512e3 / 0.27e2 * t842) * t46 - 0.6e1 * t91 * t327 + 0.3e1 * t433 * t46 + 0.6e1 * t91 * t331 - 0.6e1 * t180 * t66 * t79 - 0.3e1 * t412 * t163 + 0.6e1 * t90 * t140 * t141 - 0.63e2 / 0.2e1 * t862 - 0.14e2 * t865 + 0.21e2 / 0.2e1 * t868 + 0.70e2 * t874 + 0.35e2 / 0.3e1 * t877 - 0.7e1 * t879 + 0.21e2 / 0.2e1 * t881 + 0.7e1 / 0.2e1 * t885 + 0.21e2 * t888 + 0.42e2 * t894 + 0.2485e4 / 0.27e2 * t898;
1442   t903 = t31 * t61 * t66 * t204;
1443   t907 = t200 * t134 * t66 * t204;
1444   t919 = t201 * t358 * t97 * t158 * t141;
1445   t923 = t185 * t134 * t95 * t100;
1446   t927 = t93 * t315 * t95 * t100;
1447   t930 = t392 * t194 * t379;
1448   t934 = t377 * t97 * t158 * t163;
1449   t938 = t86 * t61 * t95 * t385;
1450   t946 = 0.1e1 / t7 / t117 * t15 * t95 * t97 * t747 * t2;
1451   t950 = t9 * t61 * t95 * t100;
1452   t953 = t388 * t194 * t379;
1453   t957 = t56 * t15 * t95 * t385;
1454   t961 = t49 * t15 * t95 * t100;
1455   t968 = t94 * t358 * params->d * t75 * t98 * t79 * t163;
1456   t970 = -0.385e3 / 0.6e1 * t903 + 0.35e2 / 0.2e1 * t907 + 0.6e1 * t181 * t142 - 0.3e1 * t434 * t80 - 0.3e1 * t181 * t164 - t91 * t356 + 0.105e3 * t919 + 0.7e1 / 0.2e1 * t923 - 0.7e1 / 0.4e1 * t927 - 0.70e2 * t930 - 0.35e2 * t934 - 0.35e2 * t938 + 0.70e2 / 0.9e1 * t946 - 0.35e2 / 0.6e1 * t950 + 0.385e3 / 0.3e1 * t953 + 0.105e3 * t957 + 0.140e3 / 0.27e2 * t961 - 0.63e2 / 0.2e1 * t968;
1457   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1458     v4rho3sigma[0] = t900 + t970;
1459 
1460   t1007 = 0.3e1 * t461 * t46 + t4 * (-0.2464e4 / 0.27e2 * t825 + 0.8096e4 / 0.9e1 * t828 - 0.32800e5 / 0.27e2 * t831 + 0.11392e5 / 0.27e2 * t836 - 0.1024e4 / 0.27e2 * t842) * t46 - 0.3e1 * t454 * t163 + 0.6e1 * t103 * t140 * t141 - 0.6e1 * t211 * t66 * t79 - 0.6e1 * t104 * t327 + 0.6e1 * t212 * t142 - 0.3e1 * t462 * t80 - 0.3e1 * t212 * t164 - t104 * t356 + 0.6e1 * t104 * t331 - 0.63e2 * t862 - 0.28e2 * t865 + 0.21e2 * t868 + 0.140e3 * t874 + 0.70e2 / 0.3e1 * t877 - 0.14e2 * t879 + 0.21e2 * t881;
1461   t1026 = 0.7e1 * t885 + 0.42e2 * t888 + 0.84e2 * t894 + 0.4970e4 / 0.27e2 * t898 - 0.385e3 / 0.3e1 * t903 + 0.35e2 * t907 + 0.210e3 * t919 + 0.7e1 * t923 - 0.7e1 / 0.2e1 * t927 - 0.140e3 * t930 - 0.70e2 * t934 - 0.70e2 * t938 + 0.140e3 / 0.9e1 * t946 - 0.35e2 / 0.3e1 * t950 + 0.770e3 / 0.3e1 * t953 + 0.210e3 * t957 + 0.280e3 / 0.27e2 * t961 - 0.63e2 * t968;
1462   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1463     v4rho3sigma[1] = t1007 + t1026;
1464 
1465   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1466     v4rho3sigma[2] = v4rho3sigma[0];
1467 
1468   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1469     v4rho3sigma[3] = v4rho3sigma[2];
1470 
1471   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1472     v4rho3sigma[4] = v4rho3sigma[1];
1473 
1474   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1475     v4rho3sigma[5] = v4rho3sigma[3];
1476 
1477   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1478     v4rho3sigma[6] = v4rho3sigma[5];
1479 
1480   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1481     v4rho3sigma[7] = v4rho3sigma[4];
1482 
1483   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1484     v4rho3sigma[8] = v4rho3sigma[6];
1485 
1486   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1487     v4rho3sigma[9] = v4rho3sigma[8];
1488 
1489   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1490     v4rho3sigma[10] = v4rho3sigma[7];
1491 
1492   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1493     v4rho3sigma[11] = v4rho3sigma[9];
1494 
1495   t1032 = 0.14e2 / 0.3e1 * t185 * t180 * t95 * t100;
1496   t1036 = 0.7e1 / 0.2e1 * t93 * t433 * t95 * t100;
1497   t1039 = 0.140e3 / 0.3e1 * t493 * t194 * t379;
1498   t1043 = 0.70e2 / 0.3e1 * t86 * t90 * t95 * t385;
1499   t1045 = t497 * t505 * t509;
1500   t1048 = t501 * t505 * t509;
1501   t1052 = t506 * t507 * t29 * t163;
1502   t1055 = t392 * t95 * t249;
1503   t1057 = t383 * t529;
1504   t1062 = t244 * t134 * t95 * t249;
1505   t1066 = t525 * t61 * t95 * t529;
1506   t1068 = t370 * t254;
1507   t1070 = t374 * t254;
1508   t1075 = 0.35e2 / 0.9e1 * t9 * t90 * t95 * t100;
1509   t1080 = t234 * t808 * t236 * t507 * t29 * t141;
1510   t1083 = t388 * t95 * t249;
1511   t1085 = t397 * t254;
1512   t1092 = -0.35e2 / 0.16e2 * t1062 + 0.35e2 / 0.4e1 * t1066 - 0.7e1 / 0.6e1 * t1068 + 0.7e1 / 0.8e1 * t1070 - t1075 + 0.147e3 / 0.2e1 * t1080 - 0.175e3 / 0.72e2 * t1083 + 0.35e2 / 0.36e2 * t1085 + 0.2e1 * t228 * t142 - 0.2e1 * t479 * t80 - t228 * t164;
1513   t1096 = t359 * t97 * t253 * t141;
1514   t1100 = 0.28e2 / 0.3e1 * t482 * t194 * t197;
1515   t1103 = 0.14e2 * t486 * t194 * t197;
1516   t1105 = 0.7e1 * t490 * t407;
1517   t1106 = t377 * t521;
1518   t1109 = t515 * t194 * t521;
1519   t1113 = t519 * t97 * t248 * t163;
1520   t1118 = t526 * t194 * t97 * t528 * t79;
1521   t1120 = t400 * t537;
1522   t1122 = t403 * t537;
1523   t1127 = t195 * t97 * t253 * t163;
1524   t1131 = 0.21e2 * t230 * t358 * t361;
1525   t1135 = t245 * t358 * t97 * t248 * t141;
1526   t1140 = t418 * t220;
1527   t1142 = t421 * t469;
1528   t1147 = params->b * t428 * t310 * t52 * t12;
1529   t1154 = t125 / t6 / t126 / t117 * t793 * t12;
1530   t1161 = t300 * t15 * t140 * t241;
1531   t1166 = 0.385e3 / 0.9e1 * t31 * t90 * t66 * t204;
1532   t1170 = 0.70e2 / 0.3e1 * t200 * t180 * t66 * t204;
1533   t1173 = t128 * t61 * t140 * t241;
1534   t1177 = t175 * t134 * t140 * t241;
1535   t1182 = t765 * t15 * t66 * t203 * t748;
1536   t1184 = -0.7e1 / 0.4e1 * t1127 - t1131 - 0.105e3 / 0.8e1 * t1135 - 0.2e1 * t227 * t66 * t79 + t4 * (-0.608e3 / 0.9e1 * t1140 + 0.1336e4 / 0.9e1 * t1142 - 0.200e3 / 0.3e1 * t1147 + 0.64e2 / 0.9e1 * t1154) * t46 + 0.2205e4 / 0.4e1 * t1161 - t1166 + t1170 - 0.441e3 / 0.4e1 * t1173 + 0.49e2 / 0.8e1 * t1177 - 0.35e2 / 0.12e2 * t1182;
1537   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1538     v4rho2sigma2[0] = 0.2e1 * t478 * t46 + t1032 - t1036 - t1039 - t1043 + 0.1323e4 / 0.4e1 * t1045 - 0.147e3 / 0.4e1 * t1048 - 0.147e3 / 0.8e1 * t1052 + 0.35e2 / 0.12e2 * t1055 - 0.175e3 / 0.8e1 * t1057 + t1092 + 0.21e2 / 0.4e1 * t1096 - t1100 + t1103 + t1105 - 0.35e2 / 0.6e1 * t1106 + 0.35e2 / 0.4e1 * t1109 + 0.35e2 / 0.8e1 * t1113 - 0.35e2 / 0.2e1 * t1118 + 0.7e1 / 0.3e1 * t1120 - 0.7e1 / 0.2e1 * t1122 + t1184;
1539 
1540   t1188 = t185 * t211 * t95 * t100;
1541   t1192 = t93 * t461 * t95 * t100;
1542   t1195 = t559 * t194 * t379;
1543   t1199 = t86 * t103 * t95 * t385;
1544   t1203 = t9 * t103 * t95 * t100;
1545   t1209 = 0.7e1 / 0.3e1 * t1188 - 0.7e1 / 0.4e1 * t1192 - 0.70e2 / 0.3e1 * t1195 - 0.35e2 / 0.3e1 * t1199 - 0.35e2 / 0.18e2 * t1203 + t1032 - t1036 - t1039 - t1043 + 0.1323e4 / 0.2e1 * t1045 - 0.147e3 / 0.2e1 * t1048 - 0.147e3 / 0.4e1 * t1052 + 0.35e2 / 0.6e1 * t1055;
1546   t1225 = -0.175e3 / 0.4e1 * t1057 - 0.35e2 / 0.8e1 * t1062 + 0.35e2 / 0.2e1 * t1066 - 0.7e1 / 0.3e1 * t1068 + 0.7e1 / 0.4e1 * t1070 - t1075 + 0.147e3 * t1080 - 0.175e3 / 0.36e2 * t1083 + 0.35e2 / 0.18e2 * t1085 + 0.2e1 * t260 * t142 - 0.2e1 * t545 * t80 - t260 * t164 + 0.2e1 * t544 * t46;
1547   t1228 = t548 * t194 * t197;
1548   t1231 = t552 * t194 * t197;
1549   t1233 = t556 * t407;
1550   t1236 = t262 * t358 * t361;
1551   t1244 = -0.14e2 / 0.3e1 * t1228 + 0.7e1 * t1231 + 0.7e1 / 0.2e1 * t1233 - 0.21e2 / 0.2e1 * t1236 + 0.21e2 / 0.2e1 * t1096 - t1100 + t1103 + t1105 - 0.35e2 / 0.3e1 * t1106 + 0.35e2 / 0.2e1 * t1109 + 0.35e2 / 0.4e1 * t1113 - 0.35e2 * t1118 + 0.14e2 / 0.3e1 * t1120;
1552   t1264 = t31 * t103 * t66 * t204;
1553   t1268 = t200 * t211 * t66 * t204;
1554   t1270 = -0.7e1 * t1122 - 0.7e1 / 0.2e1 * t1127 - t1131 - 0.105e3 / 0.4e1 * t1135 - 0.2e1 * t259 * t66 * t79 + t4 * (-0.1216e4 / 0.9e1 * t1140 + 0.2672e4 / 0.9e1 * t1142 - 0.400e3 / 0.3e1 * t1147 + 0.128e3 / 0.9e1 * t1154) * t46 + 0.2205e4 / 0.2e1 * t1161 - t1166 + t1170 - 0.441e3 / 0.2e1 * t1173 + 0.49e2 / 0.4e1 * t1177 - 0.35e2 / 0.6e1 * t1182 - 0.385e3 / 0.18e2 * t1264 + 0.35e2 / 0.3e1 * t1268;
1555   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556     v4rho2sigma2[1] = t1209 + t1225 + t1244 + t1270;
1557 
1558   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559     v4rho2sigma2[2] = v4rho2sigma2[0];
1560 
1561   t1296 = -0.35e2 / 0.4e1 * t1062 + 0.35e2 * t1066 - 0.14e2 / 0.3e1 * t1068 + 0.7e1 / 0.2e1 * t1070 + 0.294e3 * t1080 - 0.175e3 / 0.18e2 * t1083 + 0.35e2 / 0.9e1 * t1085 + 0.2e1 * t272 * t142 - 0.2e1 * t580 * t80 - t272 * t164 - 0.56e2 / 0.3e1 * t1228;
1562   t1329 = -0.7e1 * t1127 - 0.105e3 / 0.2e1 * t1135 + 0.2e1 * t579 * t46 - 0.2e1 * t271 * t66 * t79 + t4 * (-0.2432e4 / 0.9e1 * t1140 + 0.5344e4 / 0.9e1 * t1142 - 0.800e3 / 0.3e1 * t1147 + 0.256e3 / 0.9e1 * t1154) * t46 + 0.2205e4 * t1161 - 0.441e3 * t1173 + 0.49e2 / 0.2e1 * t1177 - 0.35e2 / 0.3e1 * t1182 - 0.770e3 / 0.9e1 * t1264 + 0.140e3 / 0.3e1 * t1268;
1563   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1564     v4rho2sigma2[3] = 0.28e2 / 0.3e1 * t1188 - 0.7e1 * t1192 - 0.280e3 / 0.3e1 * t1195 - 0.140e3 / 0.3e1 * t1199 - 0.70e2 / 0.9e1 * t1203 + 0.1323e4 * t1045 - 0.147e3 * t1048 - 0.147e3 / 0.2e1 * t1052 + 0.35e2 / 0.3e1 * t1055 - 0.175e3 / 0.2e1 * t1057 + t1296 + 0.28e2 * t1231 + 0.14e2 * t1233 - 0.42e2 * t1236 + 0.21e2 * t1096 - 0.70e2 / 0.3e1 * t1106 + 0.35e2 * t1109 + 0.35e2 / 0.2e1 * t1113 - 0.70e2 * t1118 + 0.28e2 / 0.3e1 * t1120 - 0.14e2 * t1122 + t1329;
1565 
1566   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1567     v4rho2sigma2[4] = v4rho2sigma2[1];
1568 
1569   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1570     v4rho2sigma2[5] = v4rho2sigma2[2];
1571 
1572   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1573     v4rho2sigma2[6] = v4rho2sigma2[5];
1574 
1575   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1576     v4rho2sigma2[7] = v4rho2sigma2[4];
1577 
1578   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1579     v4rho2sigma2[8] = v4rho2sigma2[6];
1580 
1581   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1582     v4rho2sigma2[9] = v4rho2sigma2[3];
1583 
1584   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1585     v4rho2sigma2[10] = v4rho2sigma2[7];
1586 
1587   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1588     v4rho2sigma2[11] = v4rho2sigma2[8];
1589 
1590   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1591     v4rho2sigma2[12] = v4rho2sigma2[11];
1592 
1593   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594     v4rho2sigma2[13] = v4rho2sigma2[10];
1595 
1596   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597     v4rho2sigma2[14] = v4rho2sigma2[12];
1598 
1599   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600     v4rho2sigma2[15] = v4rho2sigma2[9];
1601 
1602   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603     v4rho2sigma2[16] = v4rho2sigma2[13];
1604 
1605   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606     v4rho2sigma2[17] = v4rho2sigma2[14];
1607 
1608   t1331 = t467 * t598;
1609   t1336 = params->b * t473 * t310 * t2 * t12;
1610   t1342 = t53 / t6 / t780 * t793 * t12;
1611   t1349 = t128 * t90 * t140 * t241;
1612   t1351 = t498 * t631;
1613   t1355 = t200 * t227 * t66 * t204;
1614   t1359 = t175 * t180 * t140 * t241;
1615   t1364 = t623 * t808 * t626 * t74 * t79;
1616   t1366 = t502 * t631;
1617   t1371 = t185 * t227 * t95 * t100;
1618   t1374 = t493 * t95 * t249;
1619   t1376 = t483 * t254;
1620   t1378 = t527 * t637;
1621   t1380 = t512 * t642;
1622   t1382 = t187 * t648;
1623   t1386 = t93 * t478 * t95 * t100;
1624   t1389 = t613 * t505 * t509;
1625   t1393 = t244 * t180 * t95 * t249;
1626   t1395 = t4 * (-0.24e2 * t1331 + 0.56e2 / 0.3e1 * t1336 - 0.8e1 / 0.3e1 * t1342) * t46 - 0.1323e4 / 0.8e1 * t1349 - 0.3969e4 / 0.32e2 * t1351 + 0.35e2 / 0.2e1 * t1355 + 0.147e3 / 0.8e1 * t1359 + 0.1029e4 / 0.8e1 * t1364 + 0.441e3 / 0.32e2 * t1366 + t606 * t46 + 0.7e1 / 0.2e1 * t1371 + 0.35e2 / 0.8e1 * t1374 - 0.7e1 / 0.4e1 * t1376 - 0.35e2 / 0.32e2 * t1378 - 0.35e2 / 0.16e2 * t1380 + 0.7e1 / 0.8e1 * t1382 - 0.21e2 / 0.4e1 * t1386 - 0.441e3 / 0.8e1 * t1389 - 0.105e3 / 0.16e2 * t1393;
1627   t1398 = t525 * t90 * t95 * t529;
1628   t1400 = t487 * t254;
1629   t1404 = t506 * t507 * t28 * t79;
1630   t1408 = t335 * t61 * t95 * t637;
1631   t1414 = t289 * t15 * t95 * t97 * t747 * t247;
1632   t1416 = t516 * t642;
1633   t1418 = t191 * t648;
1634   t1423 = t602 * t61 * t325 * t627;
1635   t1430 = 0.1e1 / t307 * t15 * t325 * t625 * t29 * t158;
1636   t1434 = t473 * t15 * t325 * t627;
1637   t1437 = t609 * t194 * t197;
1638   t1440 = t617 * t194 * t521;
1639   t1442 = t490 * t537;
1640   t1447 = t634 * t194 * t97 * t636 * t79;
1641   t1451 = t519 * t97 * t641 * t79;
1642   t1455 = t195 * t97 * t647 * t79;
1643   t1457 = 0.105e3 / 0.8e1 * t1398 + 0.21e2 / 0.8e1 * t1400 - 0.1323e4 / 0.32e2 * t1404 - 0.105e3 / 0.64e2 * t1408 + 0.35e2 / 0.32e2 * t1414 + 0.105e3 / 0.32e2 * t1416 - 0.21e2 / 0.16e2 * t1418 - t607 * t80 - 0.1029e4 / 0.32e2 * t1423 + 0.1715e4 / 0.16e2 * t1430 + 0.343e3 * t1434 + 0.21e2 / 0.2e1 * t1437 + 0.105e3 / 0.8e1 * t1440 - 0.21e2 / 0.4e1 * t1442 + 0.105e3 / 0.32e2 * t1447 - 0.105e3 / 0.16e2 * t1451 + 0.21e2 / 0.8e1 * t1455;
1644   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645     v4rhosigma3[0] = t1395 + t1457;
1646 
1647   t1464 = 0.441e3 / 0.2e1 * t1349;
1648   t1468 = t200 * t259 * t66 * t204;
1649   t1472 = t175 * t211 * t140 * t241;
1650   t1476 = t128 * t103 * t140 * t241;
1651   t1479 = 0.49e2 / 0.2e1 * t1359;
1652   t1484 = t4 * (-0.48e2 * t1331 + 0.112e3 / 0.3e1 * t1336 - 0.16e2 / 0.3e1 * t1342) * t46 - t1464 - 0.3969e4 / 0.16e2 * t1351 + 0.35e2 / 0.3e1 * t1468 + 0.49e2 / 0.8e1 * t1472 - 0.441e3 / 0.8e1 * t1476 + 0.35e2 / 0.3e1 * t1355 + t1479 + 0.1029e4 / 0.4e1 * t1364 + 0.441e3 / 0.16e2 * t1366 + t653 * t46 + 0.7e1 / 0.3e1 * t1371;
1653   t1485 = 0.35e2 / 0.6e1 * t1374;
1654   t1486 = 0.7e1 / 0.3e1 * t1376;
1655   t1491 = 0.147e3 / 0.2e1 * t1389;
1656   t1492 = 0.35e2 / 0.4e1 * t1393;
1657   t1493 = 0.35e2 / 0.2e1 * t1398;
1658   t1494 = 0.7e1 / 0.2e1 * t1400;
1659   t1497 = t1485 - t1486 - 0.35e2 / 0.16e2 * t1378 - 0.35e2 / 0.8e1 * t1380 + 0.7e1 / 0.4e1 * t1382 - 0.7e1 / 0.2e1 * t1386 - t1491 - t1492 + t1493 + t1494 - 0.1323e4 / 0.16e2 * t1404 - 0.105e3 / 0.32e2 * t1408;
1660   t1504 = t185 * t259 * t95 * t100;
1661   t1508 = t93 * t544 * t95 * t100;
1662   t1511 = t660 * t505 * t509;
1663   t1515 = t244 * t211 * t95 * t249;
1664   t1519 = t525 * t103 * t95 * t529;
1665   t1521 = t553 * t254;
1666   t1524 = t559 * t95 * t249;
1667   t1526 = t549 * t254;
1668   t1529 = 0.35e2 / 0.16e2 * t1414 + 0.105e3 / 0.16e2 * t1416 - 0.21e2 / 0.8e1 * t1418 + 0.7e1 / 0.3e1 * t1504 - 0.7e1 / 0.2e1 * t1508 - 0.147e3 / 0.8e1 * t1511 - 0.35e2 / 0.16e2 * t1515 + 0.35e2 / 0.8e1 * t1519 + 0.7e1 / 0.8e1 * t1521 + 0.35e2 / 0.24e2 * t1524 - 0.7e1 / 0.12e2 * t1526 - t654 * t80;
1669   t1534 = 0.35e2 / 0.2e1 * t1440;
1670   t1535 = 0.7e1 * t1442;
1671   t1540 = t656 * t194 * t197;
1672   t1543 = t664 * t194 * t521;
1673   t1545 = t556 * t537;
1674   t1547 = -0.1029e4 / 0.16e2 * t1423 + 0.1715e4 / 0.8e1 * t1430 + 0.686e3 * t1434 + 0.7e1 * t1437 + t1534 - t1535 + 0.105e3 / 0.16e2 * t1447 - 0.105e3 / 0.8e1 * t1451 + 0.21e2 / 0.4e1 * t1455 + 0.7e1 * t1540 + 0.35e2 / 0.8e1 * t1543 - 0.7e1 / 0.4e1 * t1545;
1675   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1676     v4rhosigma3[1] = t1484 + t1497 + t1529 + t1547;
1677 
1678   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1679     v4rhosigma3[2] = v4rhosigma3[0];
1680 
1681   t1556 = t684 * t194 * t197;
1682   t1560 = t200 * t271 * t66 * t204;
1683   t1569 = t4 * (-0.96e2 * t1331 + 0.224e3 / 0.3e1 * t1336 - 0.32e2 / 0.3e1 * t1342) * t46 + 0.7e1 / 0.2e1 * t1556 + 0.35e2 / 0.6e1 * t1560 - t1464 - 0.3969e4 / 0.8e1 * t1351 + 0.70e2 / 0.3e1 * t1468 + 0.49e2 / 0.2e1 * t1472 - 0.441e3 / 0.2e1 * t1476 + t1479 + 0.1029e4 / 0.2e1 * t1364 + 0.441e3 / 0.8e1 * t1366 + t681 * t46;
1684   t1576 = t1485 - t1486 - 0.35e2 / 0.8e1 * t1378 - 0.35e2 / 0.4e1 * t1380 + 0.7e1 / 0.2e1 * t1382 - t1491 - t1492 + t1493 + t1494 - 0.1323e4 / 0.8e1 * t1404 - 0.105e3 / 0.16e2 * t1408 + 0.35e2 / 0.8e1 * t1414;
1685   t1590 = t93 * t579 * t95 * t100;
1686   t1594 = t185 * t271 * t95 * t100;
1687   t1596 = 0.105e3 / 0.8e1 * t1416 - 0.21e2 / 0.4e1 * t1418 + 0.14e2 / 0.3e1 * t1504 - 0.7e1 * t1508 - 0.147e3 / 0.2e1 * t1511 - 0.35e2 / 0.4e1 * t1515 + 0.35e2 / 0.2e1 * t1519 + 0.7e1 / 0.2e1 * t1521 + 0.35e2 / 0.6e1 * t1524 - 0.7e1 / 0.3e1 * t1526 - 0.7e1 / 0.4e1 * t1590 + 0.7e1 / 0.6e1 * t1594;
1688   t1607 = -t682 * t80 - 0.1029e4 / 0.8e1 * t1423 + 0.1715e4 / 0.4e1 * t1430 + 0.1372e4 * t1434 + t1534 - t1535 + 0.105e3 / 0.8e1 * t1447 - 0.105e3 / 0.4e1 * t1451 + 0.21e2 / 0.2e1 * t1455 + 0.14e2 * t1540 + 0.35e2 / 0.2e1 * t1543 - 0.7e1 * t1545;
1689   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1690     v4rhosigma3[3] = t1569 + t1576 + t1596 + t1607;
1691 
1692   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1693     v4rhosigma3[4] = v4rhosigma3[1];
1694 
1695   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1696     v4rhosigma3[5] = v4rhosigma3[2];
1697 
1698   t1631 = 0.21e2 * t1556 + t699 * t46 + 0.35e2 * t1560 - 0.3969e4 / 0.4e1 * t1351 + 0.147e3 / 0.2e1 * t1472 - 0.1323e4 / 0.2e1 * t1476 + 0.1029e4 * t1364 + 0.441e3 / 0.4e1 * t1366 - t700 * t80 + t4 * (-0.192e3 * t1331 + 0.448e3 / 0.3e1 * t1336 - 0.64e2 / 0.3e1 * t1342) * t46 - 0.35e2 / 0.4e1 * t1378 - 0.35e2 / 0.2e1 * t1380 + 0.7e1 * t1382 - 0.1323e4 / 0.4e1 * t1404 - 0.105e3 / 0.8e1 * t1408 + 0.35e2 / 0.4e1 * t1414 + 0.105e3 / 0.4e1 * t1416;
1699   t1649 = -0.21e2 / 0.2e1 * t1418 - 0.441e3 / 0.2e1 * t1511 - 0.105e3 / 0.4e1 * t1515 + 0.105e3 / 0.2e1 * t1519 + 0.21e2 / 0.2e1 * t1521 + 0.35e2 / 0.2e1 * t1524 - 0.7e1 * t1526 - 0.21e2 / 0.2e1 * t1590 + 0.7e1 * t1594 - 0.1029e4 / 0.4e1 * t1423 + 0.1715e4 / 0.2e1 * t1430 + 0.2744e4 * t1434 + 0.105e3 / 0.4e1 * t1447 - 0.105e3 / 0.2e1 * t1451 + 0.21e2 * t1455 + 0.105e3 / 0.2e1 * t1543 - 0.21e2 * t1545;
1700   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1701     v4rhosigma3[6] = t1631 + t1649;
1702 
1703   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1704     v4rhosigma3[7] = v4rhosigma3[3];
1705 
1706   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1707     v4rhosigma3[8] = v4rhosigma3[4];
1708 
1709   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1710     v4rhosigma3[9] = v4rhosigma3[5];
1711 
1712   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1713     v4rhosigma3[10] = v4rhosigma3[9];
1714 
1715   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1716     v4rhosigma3[11] = v4rhosigma3[8];
1717 
1718   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1719     v4rhosigma3[12] = v4rhosigma3[10];
1720 
1721   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1722     v4rhosigma3[13] = v4rhosigma3[7];
1723 
1724   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1725     v4rhosigma3[14] = v4rhosigma3[11];
1726 
1727   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1728     v4rhosigma3[15] = v4rhosigma3[12];
1729 
1730   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1731     v4rhosigma3[16] = v4rhosigma3[6];
1732 
1733   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1734     v4rhosigma3[17] = v4rhosigma3[13];
1735 
1736   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1737     v4rhosigma3[18] = v4rhosigma3[14];
1738 
1739   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1740     v4rhosigma3[19] = v4rhosigma3[15];
1741 
1742   t1652 = t602 * t90 * t325 * t627;
1743   t1656 = t623 * t325 * t625 * t74;
1744   t1662 = t775 * t15 * t325 * t625 * t28 * t158;
1745   t1666 = t175 * t227 * t140 * t241;
1746   t1667 = 0.147e3 / 0.4e1 * t1666;
1747   t1672 = t237 * t237;
1748   t1676 = 0.1e1 / t789 / t47 * t15 * t808 * t1672 * t124 * t18 * t26;
1749   t1680 = params->b * t602 * t310 * t12;
1750   t1686 = t3 / t6 / t307 * t793 * t12;
1751   t1692 = t93 * t606 * t95 * t100;
1752   t1696 = t244 * t227 * t95 * t249;
1753   t1697 = 0.105e3 / 0.8e1 * t1696;
1754   t1698 = t614 * t631;
1755   t1702 = t335 * t90 * t95 * t637;
1756   t1706 = t235 * t238 * t239 * t98;
1757   t1712 = t113 * t15 * t95 * t97 * t747 * t640;
1758   t1714 = t610 * t254;
1759   t1715 = 0.21e2 / 0.4e1 * t1714;
1760   t1716 = t618 * t642;
1761   t1718 = t231 * t648;
1762   t1722 = t635 * t97 * t36 * t646;
1763   t1727 = t246 * t97 * t158 / t124;
1764   t1733 = t96 * t97 * t74 / t28 / t124;
1765   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1766     v4sigma4[0] = -0.1029e4 / 0.8e1 * t1652 - 0.13377e5 / 0.128e3 * t1656 - 0.5145e4 / 0.128e3 * t1662 + t1667 + 0.7203e4 / 0.32e2 * t1676 + t4 * (-0.4e1 * t1680 + t1686) * t46 - 0.7e1 * t1692 - t1697 + 0.441e3 / 0.8e1 * t1698 - 0.105e3 / 0.16e2 * t1702 + 0.735e3 / 0.128e3 * t1706 - 0.105e3 / 0.256e3 * t1712 + t1715 + 0.105e3 / 0.8e1 * t1716 - 0.21e2 / 0.4e1 * t1718 + 0.315e3 / 0.64e2 * t1722 - 0.525e3 / 0.64e2 * t1727 + 0.105e3 / 0.32e2 * t1733;
1767 
1768   t1737 = t602 * t103 * t325 * t627;
1769   t1747 = t175 * t259 * t140 * t241;
1770   t1749 = t661 * t631;
1771   t1757 = -0.1029e4 / 0.32e2 * t1737 - 0.3087e4 / 0.16e2 * t1652 - 0.5145e4 / 0.64e2 * t1662 - 0.13377e5 / 0.64e2 * t1656 + t1667 + 0.7203e4 / 0.16e2 * t1676 + 0.1323e4 / 0.16e2 * t1698 + 0.735e3 / 0.64e2 * t1706 + 0.147e3 / 0.8e1 * t1747 + 0.441e3 / 0.32e2 * t1749 + t4 * (-0.8e1 * t1680 + 0.2e1 * t1686) * t46 - 0.7e1 / 0.2e1 * t1692 - t1697;
1772   t1767 = t93 * t653 * t95 * t100;
1773   t1771 = t244 * t259 * t95 * t249;
1774   t1775 = t335 * t103 * t95 * t637;
1775   t1777 = t657 * t254;
1776   t1779 = t665 * t642;
1777   t1781 = t263 * t648;
1778   t1783 = -0.315e3 / 0.32e2 * t1702 - 0.105e3 / 0.128e3 * t1712 + t1715 + 0.315e3 / 0.16e2 * t1716 - 0.63e2 / 0.8e1 * t1718 + 0.315e3 / 0.32e2 * t1722 - 0.525e3 / 0.32e2 * t1727 + 0.105e3 / 0.16e2 * t1733 - 0.21e2 / 0.4e1 * t1767 - 0.105e3 / 0.16e2 * t1771 - 0.105e3 / 0.64e2 * t1775 + 0.21e2 / 0.8e1 * t1777 + 0.105e3 / 0.32e2 * t1779 - 0.21e2 / 0.16e2 * t1781;
1779   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1780     v4sigma4[1] = t1757 + t1783;
1781 
1782   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1783     v4sigma4[2] = v4sigma4[0];
1784 
1785   t1785 = 0.1029e4 / 0.4e1 * t1652;
1786   t1790 = 0.441e3 / 0.4e1 * t1698;
1787   t1796 = t175 * t271 * t140 * t241;
1788   t1804 = 0.105e3 / 0.8e1 * t1702;
1789   t1806 = -0.1029e4 / 0.8e1 * t1737 - t1785 - 0.5145e4 / 0.32e2 * t1662 - 0.13377e5 / 0.32e2 * t1656 + 0.49e2 / 0.2e1 * t1666 + 0.7203e4 / 0.8e1 * t1676 + t1790 + 0.735e3 / 0.32e2 * t1706 + 0.49e2 * t1747 + 0.441e3 / 0.8e1 * t1749 + 0.49e2 / 0.8e1 * t1796 + t4 * (-0.16e2 * t1680 + 0.4e1 * t1686) * t46 - 0.35e2 / 0.4e1 * t1696 - t1804 - 0.105e3 / 0.64e2 * t1712;
1790   t1808 = 0.105e3 / 0.4e1 * t1716;
1791   t1809 = 0.21e2 / 0.2e1 * t1718;
1792   t1821 = t93 * t681 * t95 * t100;
1793   t1825 = t244 * t271 * t95 * t249;
1794   t1827 = t685 * t254;
1795   t1829 = 0.7e1 / 0.2e1 * t1714 + t1808 - t1809 + 0.315e3 / 0.16e2 * t1722 - 0.525e3 / 0.16e2 * t1727 + 0.105e3 / 0.8e1 * t1733 - 0.7e1 * t1767 - 0.35e2 / 0.2e1 * t1771 - 0.105e3 / 0.16e2 * t1775 + 0.7e1 * t1777 + 0.105e3 / 0.8e1 * t1779 - 0.21e2 / 0.4e1 * t1781 - 0.7e1 / 0.2e1 * t1821 - 0.35e2 / 0.16e2 * t1825 + 0.7e1 / 0.8e1 * t1827;
1796   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1797     v4sigma4[3] = t1806 + t1829;
1798 
1799   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1800     v4sigma4[4] = v4sigma4[1];
1801 
1802   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1803     v4sigma4[5] = v4sigma4[2];
1804 
1805   t1844 = -0.3087e4 / 0.8e1 * t1737 - t1785 - 0.5145e4 / 0.16e2 * t1662 - 0.13377e5 / 0.16e2 * t1656 + 0.7203e4 / 0.4e1 * t1676 + t1790 + 0.735e3 / 0.16e2 * t1706 + 0.147e3 / 0.2e1 * t1747 + 0.1323e4 / 0.8e1 * t1749 + 0.147e3 / 0.4e1 * t1796 + t4 * (-0.32e2 * t1680 + 0.8e1 * t1686) * t46 - t1804 - 0.105e3 / 0.32e2 * t1712;
1806   t1858 = t93 * t699 * t95 * t100;
1807   t1860 = t1808 - t1809 + 0.315e3 / 0.8e1 * t1722 - 0.525e3 / 0.8e1 * t1727 + 0.105e3 / 0.4e1 * t1733 - 0.105e3 / 0.4e1 * t1771 - 0.315e3 / 0.16e2 * t1775 + 0.21e2 / 0.2e1 * t1777 + 0.315e3 / 0.8e1 * t1779 - 0.63e2 / 0.4e1 * t1781 - 0.21e2 / 0.2e1 * t1821 - 0.105e3 / 0.8e1 * t1825 + 0.21e2 / 0.4e1 * t1827 - 0.7e1 / 0.4e1 * t1858;
1808   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1809     v4sigma4[6] = t1844 + t1860;
1810 
1811   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1812     v4sigma4[7] = v4sigma4[3];
1813 
1814   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1815     v4sigma4[8] = v4sigma4[4];
1816 
1817   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1818     v4sigma4[9] = v4sigma4[5];
1819 
1820   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1821     v4sigma4[10] = -0.13377e5 / 0.8e1 * t1656 - 0.5145e4 / 0.8e1 * t1662 - 0.1029e4 * t1737 + 0.7203e4 / 0.2e1 * t1676 + 0.147e3 * t1796 + t4 * (-0.64e2 * t1680 + 0.16e2 * t1686) * t46 + 0.735e3 / 0.8e1 * t1706 - 0.105e3 / 0.16e2 * t1712 + 0.441e3 * t1749 - 0.105e3 / 0.2e1 * t1775 - 0.105e3 / 0.2e1 * t1825 - 0.14e2 * t1858 + 0.315e3 / 0.4e1 * t1722 - 0.525e3 / 0.4e1 * t1727 + 0.105e3 / 0.2e1 * t1733 + 0.105e3 * t1779 - 0.42e2 * t1781 + 0.21e2 * t1827;
1822 
1823   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1824     v4sigma4[11] = v4sigma4[6];
1825 
1826   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1827     v4sigma4[12] = v4sigma4[7];
1828 
1829   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1830     v4sigma4[13] = v4sigma4[8];
1831 
1832   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1833     v4sigma4[14] = v4sigma4[9];
1834 
1835 #ifndef XC_DONT_COMPILE_MXC
1836 
1837   if(order < 5) return;
1838 
1839 
1840 #endif
1841 
1842 #endif
1843 
1844 #endif
1845 
1846 #endif
1847 
1848 #endif
1849 
1850 
1851 }
1852 
1853