1 /*
2   This file was generated automatically with ./scripts/maple2c.pl.
3   Do not edit this file directly as it can be overwritten!!
4 
5   This Source Code Form is subject to the terms of the Mozilla Public
6   License, v. 2.0. If a copy of the MPL was not distributed with this
7   file, You can obtain one at http://mozilla.org/MPL/2.0/.
8 
9   Maple version     : Maple 2016 (X86 64 LINUX)
10   Maple source      : ./maple/gga_exc/gga_x_rge2.mpl
11   Type of functional: gga_exc
12 */
13 
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16 
17 
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21 
22 #ifndef XC_DONT_COMPILE_EXC
23   double t1, t3, t4, t5, t6, t7, t8, t9;
24   double t10, t11, t12, t13, t14, t15, t17, t18;
25   double t19, t20, t22, t26, t28, t29, t30, t31;
26   double t32, t33, t35, t39, t42, t44;
27 
28 #ifndef XC_DONT_COMPILE_VXC
29   double t48, t49, t50, t51, t52, t53, t55, t59;
30   double t61, t65, t66, t70, t73, t77, t79;
31 
32 #ifndef XC_DONT_COMPILE_FXC
33   double t82, t83, t91, t92, t93, t98, t104, t108;
34   double t109, t113, t114, t117, t118, t119, t120, t121;
35   double t130, t132, t135, t137, t141, t144;
36 
37 #ifndef XC_DONT_COMPILE_KXC
38   double t148, t149, t162, t163, t164, t165, t169, t174;
39   double t178, t180, t184, t185, t189, t190, t193, t198;
40   double t200, t204, t208, t217, t219, t225, t229, t233;
41   double t237, t238, t240, t243, t245, t249;
42 
43 #ifndef XC_DONT_COMPILE_LXC
44   double t276, t277, t286, t302, t313, t382, t396, t429;
45   double t441;
46 #endif
47 
48 #endif
49 
50 #endif
51 
52 #endif
53 
54 #endif
55 
56 
57 
58   t1 = M_CBRT3;
59   t3 = POW_1_3(0.1e1 / M_PI);
60   t4 = t1 * t3;
61   t5 = M_CBRT4;
62   t6 = t5 * t5;
63   t7 = t4 * t6;
64   t8 = M_CBRT2;
65   t9 = t8 * t8;
66   t10 = POW_1_3(rho[0]);
67   t11 = t9 * t10;
68   t12 = M_CBRT6;
69   t13 = M_PI * M_PI;
70   t14 = POW_1_3(t13);
71   t15 = t14 * t14;
72   t17 = t12 / t15;
73   t18 = sigma[0] * t9;
74   t19 = rho[0] * rho[0];
75   t20 = t10 * t10;
76   t22 = 0.1e1 / t20 / t19;
77   t26 = t12 * t12;
78   t28 = 0.1e1 / t14 / t13;
79   t29 = t26 * t28;
80   t30 = sigma[0] * sigma[0];
81   t31 = t30 * t8;
82   t32 = t19 * t19;
83   t33 = t32 * rho[0];
84   t35 = 0.1e1 / t10 / t33;
85   t39 = 0.8040e0 + 0.5e1 / 0.972e3 * t17 * t18 * t22 + 0.65823568907145082055e-4 * t29 * t31 * t35;
86   t42 = 0.18040e1 - 0.64641600e0 / t39;
87   t44 = t7 * t11 * t42;
88   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
89     zk[0] = -0.3e1 / 0.16e2 * t44;
90 
91 #ifndef XC_DONT_COMPILE_VXC
92 
93   if(order < 1) return;
94 
95 
96   t48 = t10 * rho[0] * t1;
97   t49 = t48 * t3;
98   t50 = t6 * t9;
99   t51 = t39 * t39;
100   t52 = 0.1e1 / t51;
101   t53 = t19 * rho[0];
102   t55 = 0.1e1 / t20 / t53;
103   t59 = t32 * t19;
104   t61 = 0.1e1 / t10 / t59;
105   t65 = -0.10e2 / 0.729e3 * t17 * t18 * t55 - 0.35105903417144043763e-3 * t29 * t31 * t61;
106   t66 = t52 * t65;
107   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
108     vrho[0] = -t44 / 0.4e1 - 0.12120300000000000000e0 * t49 * t50 * t66;
109 
110   t70 = t9 * t22;
111   t73 = sigma[0] * t8;
112   t77 = 0.5e1 / 0.972e3 * t17 * t70 + 0.13164713781429016411e-3 * t29 * t73 * t35;
113   t79 = t50 * t52 * t77;
114   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
115     vsigma[0] = -0.12120300000000000000e0 * t49 * t79;
116 
117 #ifndef XC_DONT_COMPILE_FXC
118 
119   if(order < 2) return;
120 
121 
122   t82 = 0.1e1 / t20;
123   t83 = t9 * t82;
124   t91 = 0.1e1 / t51 / t39;
125   t92 = t65 * t65;
126   t93 = t91 * t92;
127   t98 = 0.1e1 / t20 / t32;
128   t104 = 0.1e1 / t10 / t32 / t53;
129   t108 = 0.110e3 / 0.2187e4 * t17 * t18 * t98 + 0.22233738830857894383e-2 * t29 * t31 * t104;
130   t109 = t52 * t108;
131   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
132     v2rho2[0] = -t7 * t83 * t42 / 0.12e2 - 0.32320800000000000000e0 * t7 * t11 * t66 + 0.24240600000000000000e0 * t49 * t50 * t93 - 0.12120300000000000000e0 * t49 * t50 * t109;
133 
134   t113 = t10 * t1;
135   t114 = t113 * t3;
136   t117 = t3 * t6;
137   t118 = t48 * t117;
138   t119 = t9 * t91;
139   t120 = t77 * t65;
140   t121 = t119 * t120;
141   t130 = -0.10e2 / 0.729e3 * t17 * t9 * t55 - 0.70211806834288087525e-3 * t29 * t73 * t61;
142   t132 = t50 * t52 * t130;
143   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
144     v2rhosigma[0] = -0.16160400000000000000e0 * t114 * t79 + 0.24240600000000000000e0 * t118 * t121 - 0.12120300000000000000e0 * t49 * t132;
145 
146   t135 = t77 * t77;
147   t137 = t50 * t91 * t135;
148   t141 = 0.1e1 / t32 * t1;
149   t144 = t6 * t52 * t29;
150   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
151     v2sigma2[0] = 0.24240600000000000000e0 * t49 * t137 - 0.31912056089010821521e-4 * t141 * t3 * t144;
152 
153 #ifndef XC_DONT_COMPILE_KXC
154 
155   if(order < 3) return;
156 
157 
158   t148 = 0.1e1 / t20 / rho[0];
159   t149 = t9 * t148;
160   t162 = t51 * t51;
161   t163 = 0.1e1 / t162;
162   t164 = t92 * t65;
163   t165 = t163 * t164;
164   t169 = t65 * t108;
165   t174 = 0.1e1 / t20 / t33;
166   t178 = t32 * t32;
167   t180 = 0.1e1 / t10 / t178;
168   t184 = -0.1540e4 / 0.6561e4 * t17 * t18 * t174 - 0.16304741809295789214e-1 * t29 * t31 * t180;
169   t185 = t52 * t184;
170   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
171     v3rho3[0] = t7 * t149 * t42 / 0.18e2 - 0.16160400000000000000e0 * t7 * t83 * t66 + 0.96962400000000000000e0 * t7 * t11 * t93 - 0.48481200000000000000e0 * t7 * t11 * t109 - 0.72721800000000000000e0 * t49 * t50 * t165 + 0.72721800000000000000e0 * t118 * t119 * t169 - 0.12120300000000000000e0 * t49 * t50 * t185;
172 
173   t189 = t82 * t1;
174   t190 = t189 * t3;
175   t193 = t113 * t117;
176   t198 = t9 * t163;
177   t200 = t198 * t77 * t92;
178   t204 = t119 * t130 * t65;
179   t208 = t119 * t77 * t108;
180   t217 = 0.110e3 / 0.2187e4 * t17 * t9 * t98 + 0.44467477661715788766e-2 * t29 * t73 * t104;
181   t219 = t50 * t52 * t217;
182   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
183     v3rho2sigma[0] = -0.53868000000000000000e-1 * t190 * t79 + 0.64641600000000000000e0 * t193 * t121 - 0.32320800000000000000e0 * t114 * t132 - 0.72721800000000000000e0 * t118 * t200 + 0.48481200000000000000e0 * t118 * t204 + 0.24240600000000000000e0 * t118 * t208 - 0.12120300000000000000e0 * t49 * t219;
184 
185   t225 = t198 * t135 * t65;
186   t229 = t119 * t77 * t130;
187   t233 = 0.1e1 / t33 * t1;
188   t237 = t141 * t117;
189   t238 = t91 * t26;
190   t240 = t238 * t28 * t65;
191   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
192     v3rhosigma2[0] = 0.32320800000000000000e0 * t114 * t137 - 0.72721800000000000000e0 * t118 * t225 + 0.48481200000000000000e0 * t118 * t229 + 0.12764822435604328608e-3 * t233 * t3 * t144 + 0.63824112178021643042e-4 * t237 * t240;
193 
194   t243 = t135 * t77;
195   t245 = t50 * t163 * t243;
196   t249 = t91 * t77 * t29;
197   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
198     v3sigma3[0] = -0.72721800000000000000e0 * t49 * t245 + 0.19147233653406492912e-3 * t237 * t249;
199 
200 #ifndef XC_DONT_COMPILE_LXC
201 
202   if(order < 4) return;
203 
204 
205   t276 = 0.1e1 / t162 / t39;
206   t277 = t92 * t92;
207   t286 = t108 * t108;
208   t302 = 0.1e1 / t10 / t178 / rho[0];
209   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
210     v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t70 * t42 + 0.14364800000000000000e0 * t7 * t149 * t66 + 0.64641600000000000000e0 * t7 * t83 * t93 - 0.32320800000000000000e0 * t7 * t83 * t109 - 0.38784960000000000000e1 * t7 * t11 * t165 + 0.38784960000000000000e1 * t4 * t50 * t10 * t91 * t169 - 0.64641600000000000000e0 * t7 * t11 * t185 + 0.29088720000000000000e1 * t49 * t50 * t276 * t277 - 0.43633080000000000000e1 * t118 * t198 * t92 * t108 + 0.72721800000000000000e0 * t49 * t50 * t91 * t286 + 0.96962400000000000000e0 * t118 * t119 * t65 * t184 - 0.12120300000000000000e0 * t49 * t50 * t52 * (0.26180e5 / 0.19683e5 * t17 * t18 / t20 / t59 + 0.13587284841079824345e0 * t29 * t31 * t302);
211 
212   t313 = t9 * t276;
213   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
214     v4rho3sigma[0] = -0.29088720000000000000e1 * t193 * t200 + 0.29088720000000000000e1 * t118 * t313 * t77 * t164 - 0.21816540000000000000e1 * t118 * t198 * t120 * t108 + 0.35912000000000000000e-1 * t148 * t1 * t3 * t79 + 0.32320800000000000000e0 * t189 * t117 * t121 + 0.19392480000000000000e1 * t193 * t204 + 0.96962400000000000000e0 * t193 * t208 - 0.21816540000000000000e1 * t118 * t198 * t130 * t92 + 0.72721800000000000000e0 * t118 * t119 * t217 * t65 + 0.72721800000000000000e0 * t118 * t119 * t130 * t108 + 0.24240600000000000000e0 * t118 * t119 * t77 * t184 - 0.16160400000000000000e0 * t190 * t132 - 0.48481200000000000000e0 * t114 * t219 - 0.12120300000000000000e0 * t49 * t50 * t52 * (-0.1540e4 / 0.6561e4 * t17 * t9 * t174 - 0.32609483618591578428e-1 * t29 * t73 * t180);
215 
216   t382 = t130 * t130;
217   t396 = t233 * t117;
218   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
219     v4rho2sigma2[0] = 0.10773600000000000000e0 * t190 * t137 - 0.19392480000000000000e1 * t193 * t225 + 0.12928320000000000000e1 * t193 * t229 + 0.29088720000000000000e1 * t118 * t313 * t135 * t92 - 0.29088720000000000000e1 * t118 * t198 * t120 * t130 - 0.72721800000000000000e0 * t118 * t198 * t135 * t108 + 0.48481200000000000000e0 * t49 * t50 * t91 * t382 + 0.48481200000000000000e0 * t118 * t119 * t77 * t217 - 0.63824112178021643040e-3 / t59 * t1 * t3 * t144 - 0.51059289742417314433e-3 * t396 * t240 - 0.19147233653406492913e-3 * t237 * t163 * t26 * t28 * t92 + 0.63824112178021643042e-4 * t237 * t238 * t28 * t108;
220 
221   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
222     v4rhosigma3[0] = -0.96962400000000000000e0 * t114 * t245 + 0.29088720000000000000e1 * t118 * t313 * t243 * t65 - 0.21816540000000000000e1 * t118 * t198 * t135 * t130 - 0.76588934613625971648e-3 * t396 * t249 - 0.57441700960219478736e-3 * t237 * t163 * t77 * t29 * t65 + 0.19147233653406492912e-3 * t237 * t91 * t130 * t29;
223 
224   t429 = t135 * t135;
225   t441 = t13 * t13;
226   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
227     v4sigma4[0] = 0.29088720000000000000e1 * t49 * t50 * t276 * t429 - 0.11488340192043895747e-2 * t237 * t163 * t135 * t29 + 0.15124071045194514738e-6 * t302 * t1 * t117 * t91 * t12 / t15 / t441 * t8;
228 
229 #ifndef XC_DONT_COMPILE_MXC
230 
231   if(order < 5) return;
232 
233 
234 #endif
235 
236 #endif
237 
238 #endif
239 
240 #endif
241 
242 #endif
243 
244 
245 }
246 
247 
248 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 *))249 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
250 {
251 
252 #ifndef XC_DONT_COMPILE_EXC
253   double t1, t3, t4, t5, t6, t7, t9, t10;
254   double t11, t12, t14, t15, t16, t18, t22, t24;
255   double t25, t26, t27, t28, t30, t34, t37, t39;
256 
257 #ifndef XC_DONT_COMPILE_VXC
258   double t43, t44, t45, t46, t47, t48, t50, t54;
259   double t56, t60, t69, t70;
260 
261 #ifndef XC_DONT_COMPILE_FXC
262   double t73, t78, t79, t84, t85, t86, t91, t97;
263   double t101, t106, t110, t118, t119, t122, t123, t127;
264   double t129;
265 
266 #ifndef XC_DONT_COMPILE_KXC
267   double t133, t138, t142, t149, t150, t151, t152, t156;
268   double t161, t165, t167, t171, t176, t184, t187, t188;
269   double t192, t200, t201, t207, t211, t215, t219, t220;
270   double t221, t223, t226, t227, t231;
271 
272 #ifndef XC_DONT_COMPILE_LXC
273   double t261, t262, t270, t285, t299, t300, t360, t373;
274   double t405, t415;
275 #endif
276 
277 #endif
278 
279 #endif
280 
281 #endif
282 
283 #endif
284 
285 
286 
287   t1 = M_CBRT3;
288   t3 = POW_1_3(0.1e1 / M_PI);
289   t4 = t1 * t3;
290   t5 = M_CBRT4;
291   t6 = t5 * t5;
292   t7 = POW_1_3(rho[0]);
293   t9 = M_CBRT6;
294   t10 = M_PI * M_PI;
295   t11 = POW_1_3(t10);
296   t12 = t11 * t11;
297   t14 = t9 / t12;
298   t15 = rho[0] * rho[0];
299   t16 = t7 * t7;
300   t18 = 0.1e1 / t16 / t15;
301   t22 = t9 * t9;
302   t24 = 0.1e1 / t11 / t10;
303   t25 = t22 * t24;
304   t26 = sigma[0] * sigma[0];
305   t27 = t15 * t15;
306   t28 = t27 * rho[0];
307   t30 = 0.1e1 / t7 / t28;
308   t34 = 0.8040e0 + 0.5e1 / 0.972e3 * t14 * sigma[0] * t18 + 0.32911784453572541027e-4 * t25 * t26 * t30;
309   t37 = 0.18040e1 - 0.64641600e0 / t34;
310   t39 = t4 * t6 * t7 * t37;
311   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
312     zk[0] = -0.3e1 / 0.8e1 * t39;
313 
314 #ifndef XC_DONT_COMPILE_VXC
315 
316   if(order < 1) return;
317 
318 
319   t43 = t7 * rho[0] * t1;
320   t44 = t43 * t3;
321   t45 = t34 * t34;
322   t46 = 0.1e1 / t45;
323   t47 = t6 * t46;
324   t48 = t15 * rho[0];
325   t50 = 0.1e1 / t16 / t48;
326   t54 = t27 * t15;
327   t56 = 0.1e1 / t7 / t54;
328   t60 = -0.10e2 / 0.729e3 * t14 * sigma[0] * t50 - 0.17552951708572021881e-3 * t25 * t26 * t56;
329   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
330     vrho[0] = -t39 / 0.2e1 - 0.24240600000000000000e0 * t44 * t47 * t60;
331 
332   t69 = 0.5e1 / 0.972e3 * t14 * t18 + 0.65823568907145082054e-4 * t25 * sigma[0] * t30;
333   t70 = t47 * t69;
334   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
335     vsigma[0] = -0.24240600000000000000e0 * t44 * t70;
336 
337 #ifndef XC_DONT_COMPILE_FXC
338 
339   if(order < 2) return;
340 
341 
342   t73 = 0.1e1 / t16;
343   t78 = t4 * t6;
344   t79 = t7 * t46;
345   t84 = 0.1e1 / t45 / t34;
346   t85 = t6 * t84;
347   t86 = t60 * t60;
348   t91 = 0.1e1 / t16 / t27;
349   t97 = 0.1e1 / t7 / t27 / t48;
350   t101 = 0.110e3 / 0.2187e4 * t14 * sigma[0] * t91 + 0.11116869415428947191e-2 * t25 * t26 * t97;
351   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
352     v2rho2[0] = -t4 * t6 * t73 * t37 / 0.6e1 - 0.64641600000000000000e0 * t78 * t79 * t60 + 0.48481200000000000000e0 * t44 * t85 * t86 - 0.24240600000000000000e0 * t44 * t47 * t101;
353 
354   t106 = t7 * t1 * t3;
355   t110 = t85 * t69 * t60;
356   t118 = -0.10e2 / 0.729e3 * t14 * t50 - 0.35105903417144043762e-3 * t25 * sigma[0] * t56;
357   t119 = t47 * t118;
358   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
359     v2rhosigma[0] = -0.32320800000000000000e0 * t106 * t70 + 0.48481200000000000000e0 * t44 * t110 - 0.24240600000000000000e0 * t44 * t119;
360 
361   t122 = t69 * t69;
362   t123 = t85 * t122;
363   t127 = 0.1e1 / t27 * t1;
364   t129 = t47 * t25;
365   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
366     v2sigma2[0] = 0.48481200000000000000e0 * t44 * t123 - 0.15956028044505410760e-4 * t127 * t3 * t129;
367 
368 #ifndef XC_DONT_COMPILE_KXC
369 
370   if(order < 3) return;
371 
372 
373   t133 = 0.1e1 / t16 / rho[0];
374   t138 = t73 * t46;
375   t142 = t7 * t84;
376   t149 = t45 * t45;
377   t150 = 0.1e1 / t149;
378   t151 = t6 * t150;
379   t152 = t86 * t60;
380   t156 = t60 * t101;
381   t161 = 0.1e1 / t16 / t28;
382   t165 = t27 * t27;
383   t167 = 0.1e1 / t7 / t165;
384   t171 = -0.1540e4 / 0.6561e4 * t14 * sigma[0] * t161 - 0.81523709046478946067e-2 * t25 * t26 * t167;
385   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
386     v3rho3[0] = t4 * t6 * t133 * t37 / 0.9e1 - 0.32320800000000000000e0 * t78 * t138 * t60 + 0.19392480000000000000e1 * t78 * t142 * t86 - 0.96962400000000000000e0 * t78 * t79 * t101 - 0.14544360000000000000e1 * t44 * t151 * t152 + 0.14544360000000000000e1 * t44 * t85 * t156 - 0.24240600000000000000e0 * t44 * t47 * t171;
387 
388   t176 = t73 * t1 * t3;
389   t184 = t151 * t69 * t86;
390   t187 = t118 * t60;
391   t188 = t85 * t187;
392   t192 = t85 * t69 * t101;
393   t200 = 0.110e3 / 0.2187e4 * t14 * t91 + 0.22233738830857894383e-2 * t25 * sigma[0] * t97;
394   t201 = t47 * t200;
395   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
396     v3rho2sigma[0] = -0.10773600000000000000e0 * t176 * t70 + 0.12928320000000000000e1 * t106 * t110 - 0.64641600000000000000e0 * t106 * t119 - 0.14544360000000000000e1 * t44 * t184 + 0.96962400000000000000e0 * t44 * t188 + 0.48481200000000000000e0 * t44 * t192 - 0.24240600000000000000e0 * t44 * t201;
397 
398   t207 = t151 * t122 * t60;
399   t211 = t85 * t69 * t118;
400   t215 = 0.1e1 / t28 * t1;
401   t219 = t3 * t6;
402   t220 = t127 * t219;
403   t221 = t84 * t22;
404   t223 = t221 * t24 * t60;
405   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
406     v3rhosigma2[0] = 0.64641600000000000000e0 * t106 * t123 - 0.14544360000000000000e1 * t44 * t207 + 0.96962400000000000000e0 * t44 * t211 + 0.63824112178021643040e-4 * t215 * t3 * t129 + 0.31912056089010821520e-4 * t220 * t223;
407 
408   t226 = t122 * t69;
409   t227 = t151 * t226;
410   t231 = t84 * t69 * t25;
411   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
412     v3sigma3[0] = -0.14544360000000000000e1 * t44 * t227 + 0.95736168267032464562e-4 * t220 * t231;
413 
414 #ifndef XC_DONT_COMPILE_LXC
415 
416   if(order < 4) return;
417 
418 
419   t261 = t6 / t149 / t34;
420   t262 = t86 * t86;
421   t270 = t101 * t101;
422   t285 = 0.1e1 / t7 / t165 / rho[0];
423   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
424     v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t18 * t37 + 0.28729600000000000000e0 * t78 * t133 * t46 * t60 + 0.12928320000000000000e1 * t78 * t73 * t84 * t86 - 0.64641600000000000000e0 * t78 * t138 * t101 - 0.77569920000000000000e1 * t78 * t7 * t150 * t152 + 0.77569920000000000000e1 * t78 * t142 * t156 - 0.12928320000000000000e1 * t78 * t79 * t171 + 0.58177440000000000000e1 * t44 * t261 * t262 - 0.87266160000000000000e1 * t44 * t151 * t86 * t101 + 0.14544360000000000000e1 * t44 * t85 * t270 + 0.19392480000000000000e1 * t44 * t85 * t60 * t171 - 0.24240600000000000000e0 * t44 * t47 * (0.26180e5 / 0.19683e5 * t14 * sigma[0] / t16 / t54 + 0.67936424205399121722e-1 * t25 * t26 * t285);
425 
426   t299 = t43 * t219;
427   t300 = t150 * t69;
428   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
429     v4rho3sigma[0] = -0.58177440000000000000e1 * t106 * t184 + 0.58177440000000000000e1 * t44 * t261 * t69 * t152 - 0.43633080000000000000e1 * t299 * t300 * t156 + 0.71824000000000000000e-1 * t133 * t1 * t3 * t70 + 0.64641600000000000000e0 * t176 * t110 + 0.38784960000000000000e1 * t106 * t188 + 0.19392480000000000000e1 * t106 * t192 - 0.43633080000000000000e1 * t44 * t151 * t118 * t86 + 0.14544360000000000000e1 * t44 * t85 * t200 * t60 + 0.14544360000000000000e1 * t44 * t85 * t118 * t101 + 0.48481200000000000000e0 * t44 * t85 * t69 * t171 - 0.32320800000000000000e0 * t176 * t119 - 0.96962400000000000000e0 * t106 * t201 - 0.24240600000000000000e0 * t44 * t47 * (-0.1540e4 / 0.6561e4 * t14 * t161 - 0.16304741809295789214e-1 * t25 * sigma[0] * t167);
430 
431   t360 = t118 * t118;
432   t373 = t215 * t219;
433   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
434     v4rho2sigma2[0] = 0.21547200000000000000e0 * t176 * t123 - 0.38784960000000000000e1 * t106 * t207 + 0.25856640000000000000e1 * t106 * t211 + 0.58177440000000000000e1 * t44 * t261 * t122 * t86 - 0.58177440000000000000e1 * t299 * t300 * t187 - 0.14544360000000000000e1 * t44 * t151 * t122 * t101 + 0.96962400000000000000e0 * t44 * t85 * t360 + 0.96962400000000000000e0 * t44 * t85 * t69 * t200 - 0.31912056089010821520e-3 / t54 * t1 * t3 * t129 - 0.25529644871208657216e-3 * t373 * t223 - 0.95736168267032464560e-4 * t220 * t150 * t22 * t24 * t86 + 0.31912056089010821520e-4 * t220 * t221 * t24 * t101;
435 
436   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
437     v4rhosigma3[0] = -0.19392480000000000000e1 * t106 * t227 + 0.58177440000000000000e1 * t44 * t261 * t226 * t60 - 0.43633080000000000000e1 * t44 * t151 * t122 * t118 - 0.38294467306812985825e-3 * t373 * t231 - 0.28720850480109739369e-3 * t220 * t300 * t25 * t60 + 0.95736168267032464562e-4 * t220 * t84 * t118 * t25;
438 
439   t405 = t122 * t122;
440   t415 = t10 * t10;
441   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
442     v4sigma4[0] = 0.58177440000000000000e1 * t44 * t261 * t405 - 0.57441700960219478737e-3 * t220 * t150 * t122 * t25 + 0.37810177612986286844e-7 * t285 * t1 * t3 * t85 * t9 / t12 / t415;
443 
444 #ifndef XC_DONT_COMPILE_MXC
445 
446   if(order < 5) return;
447 
448 
449 #endif
450 
451 #endif
452 
453 #endif
454 
455 #endif
456 
457 #endif
458 
459 
460 }
461 
462 
463 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))464 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
465 {
466 
467 #ifndef XC_DONT_COMPILE_EXC
468   double t1, t3, t4, t5, t6, t7, t8, t9;
469   double t10, t11, t13, t14, t15, t16, t17, t18;
470   double t19, t20, t21, t23, t24, t25, t26, t28;
471   double t32, t34, t35, t36, t37, t38, t40, t44;
472   double t47, t49, t51, t52, t53, t54, t55, t56;
473   double t57, t59, t63, t64, t65, t67, t71, t74;
474   double t76;
475 
476 #ifndef XC_DONT_COMPILE_VXC
477   double t78, t79, t80, t81, t82, t83, t85, t86;
478   double t88, t89, t90, t91, t92, t94, t95, t96;
479   double t97, t98, t100, t104, t106, t110, t111, t113;
480   double t114, t115, t116, t117, t119, t120, t121, t123;
481   double t124, t128, t129, t131, t132, t133, t134, t136;
482   double t137, t138, t139, t140, t142, t146, t148, t152;
483   double t153, t155, t156, t160, t161, t162, t168, t169;
484   double t173, t179, t180;
485 
486 #ifndef XC_DONT_COMPILE_FXC
487   double t184, t186, t187, t188, t189, t190, t191, t193;
488   double t194, t195, t197, t198, t199, t200, t201, t202;
489   double t204, t205, t206, t207, t208, t209, t210, t212;
490   double t213, t215, t216, t218, t219, t221, t222, t224;
491   double t225, t226, t228, t229, t231, t237, t241, t242;
492   double t244, t245, t246, t247, t248, t249, t250, t252;
493   double t253, t254, t256, t257, t258, t259, t261, t262;
494   double t263, t265, t266, t267, t270, t271, t272, t274;
495   double t277, t279, t281, t284, t285, t286, t288, t293;
496   double t294, t295, t297, t300, t302, t303, t305, t308;
497   double t309, t310, t311, t313, t316, t318, t321, t322;
498   double t324, t325, t326, t327, t328, t330, t331, t332;
499   double t333, t335, t336, t337, t338, t340, t341, t342;
500   double t343, t345, t346, t347, t349, t350, t351, t353;
501   double t354, t356, t362, t366, t367, t369, t370, t371;
502   double t375, t376, t377, t378, t379, t383, t384, t393;
503   double t394, t400, t401, t402, t406, t407, t410, t411;
504   double t414, t415, t416, t424, t426, t429, t431, t434;
505   double t435, t436, t439, t441, t444, t445, t446;
506 
507 #ifndef XC_DONT_COMPILE_KXC
508   double t454, t458, t459, t461, t462, t463, t465, t466;
509   double t467, t469, t471, t473, t474, t475, t477, t478;
510   double t479, t480, t481, t483, t484, t485, t487, t488;
511   double t489, t491, t492, t494, t495, t496, t497, t499;
512   double t500, t501, t503, t505, t507, t508, t510, t512;
513   double t513, t514, t515, t516, t518, t519, t521, t522;
514   double t523, t524, t526, t527, t528, t530, t532, t534;
515   double t535, t537, t539, t540, t541, t543, t544, t546;
516   double t547, t549, t550, t552, t555, t559, t561, t565;
517   double t566, t568, t569, t571, t573, t574, t576, t577;
518   double t578, t580, t581, t582, t584, t585, t586, t591;
519   double t593, t594, t595, t598, t599, t601, t603, t605;
520   double t606, t607, t609, t611, t613, t614, t616, t617;
521   double t618, t620, t622, t624, t626, t627, t629, t631;
522   double t633, t635, t637, t640, t641, t643, t645, t647;
523   double t648, t649, t650, t652, t654, t656, t658, t659;
524   double t661, t665, t666, t668, t670, t672, t674, t676;
525   double t679, t682, t693, t695, t698, t701, t703, t706;
526   double t711, t712, t714, t717, t719, t721, t723, t725;
527   double t727, t729, t731, t733, t735, t737, t739, t741;
528   double t743, t744, t746, t748, t750, t752, t753, t755;
529   double t757, t759, t762, t764, t766, t769, t772, t775;
530   double t778, t781, t782, t784, t786, t787, t789, t799;
531   double t801, t803, t806, t815, t817, t818, t819, t820;
532   double t822, t823, t824, t826, t827, t828, t830, t831;
533   double t832, t834, t835, t836, t838, t839, t840, t842;
534   double t843, t845, t846, t848, t849, t850, t851, t853;
535   double t854, t855, t856, t857, t858, t860, t861, t863;
536   double t864, t866, t867, t868, t869, t870, t872, t873;
537   double t875, t879, t881, t885, t886, t888, t889, t896;
538   double t897, t898, t899, t903, t907, t909, t912, t914;
539   double t915, t919, t920, t924, t928, t932, t933, t937;
540   double t941, t950, t951, t956, t960, t961, t962, t966;
541   double t973, t974, t975, t978, t982, t985, t986, t989;
542   double t990, t991, t992, t996, t999, t1002, t1003, t1006;
543   double t1008, t1010, t1013, t1014, t1017, t1021, t1022, t1025;
544   double t1026, t1032, t1033, t1037, t1040, t1041, t1044, t1045;
545   double t1048, t1049, t1050, t1053, t1054, t1057, t1058, t1066;
546   double t1068, t1072, t1073, t1074, t1078, t1082, t1086, t1087;
547   double t1088, t1089, t1090, t1091, t1092, t1095, t1096, t1097;
548   double t1098, t1099, t1102, t1103, t1107, t1109, t1113, t1114;
549   double t1115, t1116, t1117, t1118, t1119, t1122, t1123, t1126;
550   double t1127, t1130, t1131, t1134, t1135, t1138, t1139, t1142;
551   double t1143, t1146, t1147, t1148, t1149, t1150, t1153, t1154;
552   double t1157, t1159, t1162, t1163, t1166, t1168, t1171, t1172;
553 
554 #ifndef XC_DONT_COMPILE_LXC
555   double t1183, t1185, t1187, t1205, t1216, t1222, t1232, t1235;
556   double t1245, t1247, t1250, t1253, t1259, t1261, t1263, t1268;
557   double t1270, t1271, t1272, t1277, t1279, t1281, t1283, t1285;
558   double t1287, t1293, t1296, t1299, t1302, t1303, t1308, t1313;
559   double t1315, t1316, t1317, t1322, t1324, t1326, t1328, t1330;
560   double t1336, t1338, t1340, t1343, t1346, t1354, t1357, t1360;
561   double t1379, t1392, t1396, t1397, t1413, t1425, t1426, t1427;
562   double t1429, t1432, t1435, t1436, t1438, t1457, t1465, t1466;
563   double t1472, t1474, t1476, t1477, t1479, t1492, t1494, t1496;
564   double t1502, t1505, t1511, t1519, t1520, t1521, t1523, t1529;
565   double t1531, t1549, t1552, t1555, t1558, t1562, t1565, t1568;
566   double t1570, t1577, t1591, t1618, t1637, t1640, t1642, t1647;
567   double t1649, t1656, t1657, t1661, t1668, t1674, t1694, t1701;
568   double t1716, t1718, t1732, t1738, t1740, t1747, t1750, t1757;
569   double t1759, t1770, t1785, t1788, t1791, t1795, t1807, t1814;
570   double t1816, t1822, t1833, t1836, t1839, t1844, t1849, t1852;
571   double t1855, t1863, t1866, t1881, t1887, t1915, t1933, t1948;
572   double t1962, t1974, t1977, t1980, t1983, t1986, t1989, t1992;
573   double t1996, t2016, t2019, t2022, t2025, t2028, t2035, t2042;
574   double t2044, t2049, t2104, t2121, t2166, t2173, t2182, t2188;
575   double t2193, t2203, t2209, t2210, t2215, t2220, t2225, t2230;
576   double t2236, t2263, t2301, t2333, t2345, t2361, t2373, t2374;
577   double t2379, t2405, t2409, t2421, t2424, t2427, t2430, t2432;
578   double t2437, t2471, t2474, t2477, t2480, t2483, t2485, t2491;
579   double t2493, t2495, t2498, t2501, t2505, t2506, t2536, t2538;
580   double t2558, t2559, t2564, t2565, t2579, t2584, t2585, t2586;
581   double t2590, t2594, t2598, t2599, t2602, t2607, t2609, t2612;
582   double t2621, t2624, t2636, t2650, t2668, t2698, t2707, t2722;
583   double t2724, t2725, t2727, t2729, t2731, t2733, t2735, t2754;
584   double t2758, t2774, t2835, t2865, t2866, t2867, t2869, t2870;
585   double t2884, t2914, t2915, t2917, t2919, t2926, t2927, t2946;
586   double t2947, t2952, t2958, t2960, t2961, t2972, t2973, t2974;
587   double t2987, t3008, t3011, t3030, t3040, t3046, t3053, t3056;
588   double t3066, t3067, t3069, t3075, t3077, t3140, t3168, t3177;
589   double t3179, t3193, t3194, t3211, t3217, t3218, t3256, t3265;
590   double t3268, t3276;
591 #endif
592 
593 #endif
594 
595 #endif
596 
597 #endif
598 
599 #endif
600 
601 
602 
603   t1 = M_CBRT3;
604   t3 = POW_1_3(0.1e1 / M_PI);
605   t4 = t1 * t3;
606   t5 = M_CBRT4;
607   t6 = t5 * t5;
608   t7 = t4 * t6;
609   t8 = rho[0] - rho[1];
610   t9 = rho[0] + rho[1];
611   t10 = 0.1e1 / t9;
612   t11 = t8 * t10;
613   t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
614   t14 = POW_1_3(t13);
615   t15 = t14 * t13;
616   t16 = POW_1_3(t9);
617   t17 = t15 * t16;
618   t18 = M_CBRT6;
619   t19 = M_PI * M_PI;
620   t20 = POW_1_3(t19);
621   t21 = t20 * t20;
622   t23 = t18 / t21;
623   t24 = rho[0] * rho[0];
624   t25 = POW_1_3(rho[0]);
625   t26 = t25 * t25;
626   t28 = 0.1e1 / t26 / t24;
627   t32 = t18 * t18;
628   t34 = 0.1e1 / t20 / t19;
629   t35 = t32 * t34;
630   t36 = sigma[0] * sigma[0];
631   t37 = t24 * t24;
632   t38 = t37 * rho[0];
633   t40 = 0.1e1 / t25 / t38;
634   t44 = 0.8040e0 + 0.5e1 / 0.972e3 * t23 * sigma[0] * t28 + 0.32911784453572541027e-4 * t35 * t36 * t40;
635   t47 = 0.18040e1 - 0.64641600e0 / t44;
636   t49 = t7 * t17 * t47;
637   t51 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
638   t52 = POW_1_3(t51);
639   t53 = t52 * t51;
640   t54 = t53 * t16;
641   t55 = rho[1] * rho[1];
642   t56 = POW_1_3(rho[1]);
643   t57 = t56 * t56;
644   t59 = 0.1e1 / t57 / t55;
645   t63 = sigma[2] * sigma[2];
646   t64 = t55 * t55;
647   t65 = t64 * rho[1];
648   t67 = 0.1e1 / t56 / t65;
649   t71 = 0.8040e0 + 0.5e1 / 0.972e3 * t23 * sigma[2] * t59 + 0.32911784453572541027e-4 * t35 * t63 * t67;
650   t74 = 0.18040e1 - 0.64641600e0 / t71;
651   t76 = t7 * t54 * t74;
652   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
653     zk[0] = -0.3e1 / 0.8e1 * t49 - 0.3e1 / 0.8e1 * t76;
654 
655 #ifndef XC_DONT_COMPILE_VXC
656 
657   if(order < 1) return;
658 
659 
660   t78 = 0.3e1 / 0.8e1 * t49;
661   t79 = 0.3e1 / 0.8e1 * t76;
662   t80 = t14 * t16;
663   t81 = t9 * t9;
664   t82 = 0.1e1 / t81;
665   t83 = t8 * t82;
666   t85 = t10 / 0.2e1 - t83 / 0.2e1;
667   t86 = t47 * t85;
668   t88 = t7 * t80 * t86;
669   t89 = t88 / 0.2e1;
670   t90 = t16 * t16;
671   t91 = 0.1e1 / t90;
672   t92 = t15 * t91;
673   t94 = t7 * t92 * t47;
674   t95 = t94 / 0.8e1;
675   t96 = t44 * t44;
676   t97 = 0.1e1 / t96;
677   t98 = t24 * rho[0];
678   t100 = 0.1e1 / t26 / t98;
679   t104 = t37 * t24;
680   t106 = 0.1e1 / t25 / t104;
681   t110 = -0.10e2 / 0.729e3 * t23 * sigma[0] * t100 - 0.17552951708572021881e-3 * t35 * t36 * t106;
682   t111 = t97 * t110;
683   t113 = t7 * t17 * t111;
684   t114 = 0.24240600000000000000e0 * t113;
685   t115 = t52 * t16;
686   t116 = -t85;
687   t117 = t74 * t116;
688   t119 = t7 * t115 * t117;
689   t120 = t119 / 0.2e1;
690   t121 = t53 * t91;
691   t123 = t7 * t121 * t74;
692   t124 = t123 / 0.8e1;
693   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
694     vrho[0] = -t78 - t79 + t9 * (-t89 - t95 - t114 - t120 - t124);
695 
696   t128 = -t10 / 0.2e1 - t83 / 0.2e1;
697   t129 = t47 * t128;
698   t131 = t7 * t80 * t129;
699   t132 = t131 / 0.2e1;
700   t133 = -t128;
701   t134 = t74 * t133;
702   t136 = t7 * t115 * t134;
703   t137 = t136 / 0.2e1;
704   t138 = t71 * t71;
705   t139 = 0.1e1 / t138;
706   t140 = t55 * rho[1];
707   t142 = 0.1e1 / t57 / t140;
708   t146 = t64 * t55;
709   t148 = 0.1e1 / t56 / t146;
710   t152 = -0.10e2 / 0.729e3 * t23 * sigma[2] * t142 - 0.17552951708572021881e-3 * t35 * t63 * t148;
711   t153 = t139 * t152;
712   t155 = t7 * t54 * t153;
713   t156 = 0.24240600000000000000e0 * t155;
714   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
715     vrho[1] = -t78 - t79 + t9 * (-t132 - t95 - t137 - t124 - t156);
716 
717   t160 = t16 * t9 * t1;
718   t161 = t160 * t3;
719   t162 = t6 * t15;
720   t168 = 0.5e1 / 0.972e3 * t23 * t28 + 0.65823568907145082054e-4 * t35 * sigma[0] * t40;
721   t169 = t97 * t168;
722   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
723     vsigma[0] = -0.24240600000000000000e0 * t161 * t162 * t169;
724 
725   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
726     vsigma[1] = 0.0e0;
727 
728   t173 = t6 * t53;
729   t179 = 0.5e1 / 0.972e3 * t23 * t59 + 0.65823568907145082054e-4 * t35 * sigma[2] * t67;
730   t180 = t139 * t179;
731   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
732     vsigma[2] = -0.24240600000000000000e0 * t161 * t173 * t180;
733 
734 #ifndef XC_DONT_COMPILE_FXC
735 
736   if(order < 2) return;
737 
738 
739   t184 = t94 / 0.4e1;
740   t186 = t123 / 0.4e1;
741   t187 = t14 * t14;
742   t188 = 0.1e1 / t187;
743   t189 = t188 * t16;
744   t190 = t85 * t85;
745   t191 = t47 * t190;
746   t193 = t7 * t189 * t191;
747   t194 = t193 / 0.6e1;
748   t195 = t14 * t91;
749   t197 = t7 * t195 * t86;
750   t198 = t197 / 0.3e1;
751   t199 = t6 * t14;
752   t200 = t4 * t199;
753   t201 = t16 * t97;
754   t202 = t110 * t85;
755   t204 = t200 * t201 * t202;
756   t205 = 0.64641600000000000000e0 * t204;
757   t206 = t81 * t9;
758   t207 = 0.1e1 / t206;
759   t208 = t8 * t207;
760   t209 = -t82 + t208;
761   t210 = t47 * t209;
762   t212 = t7 * t80 * t210;
763   t213 = t212 / 0.2e1;
764   t215 = 0.1e1 / t90 / t9;
765   t216 = t15 * t215;
766   t218 = t7 * t216 * t47;
767   t219 = t218 / 0.12e2;
768   t221 = t7 * t92 * t111;
769   t222 = 0.16160400000000000000e0 * t221;
770   t224 = 0.1e1 / t96 / t44;
771   t225 = t110 * t110;
772   t226 = t224 * t225;
773   t228 = t7 * t17 * t226;
774   t229 = 0.48481200000000000000e0 * t228;
775   t231 = 0.1e1 / t26 / t37;
776   t237 = 0.1e1 / t25 / t37 / t98;
777   t241 = 0.110e3 / 0.2187e4 * t23 * sigma[0] * t231 + 0.11116869415428947191e-2 * t35 * t36 * t237;
778   t242 = t97 * t241;
779   t244 = t7 * t17 * t242;
780   t245 = 0.24240600000000000000e0 * t244;
781   t246 = t52 * t52;
782   t247 = 0.1e1 / t246;
783   t248 = t247 * t16;
784   t249 = t116 * t116;
785   t250 = t74 * t249;
786   t252 = t7 * t248 * t250;
787   t253 = t252 / 0.6e1;
788   t254 = t52 * t91;
789   t256 = t7 * t254 * t117;
790   t257 = t256 / 0.3e1;
791   t258 = -t209;
792   t259 = t74 * t258;
793   t261 = t7 * t115 * t259;
794   t262 = t261 / 0.2e1;
795   t263 = t53 * t215;
796   t265 = t7 * t263 * t74;
797   t266 = t265 / 0.12e2;
798   t267 = -t194 - t198 - t205 - t213 + t219 - t222 + t229 - t245 - t253 - t257 - t262 + t266;
799   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
800     v2rho2[0] = -t88 - t184 - 0.48481200000000000000e0 * t113 - t119 - t186 + t9 * t267;
801 
802   t270 = t4 * t6 * t188;
803   t271 = t16 * t47;
804   t272 = t128 * t85;
805   t274 = t270 * t271 * t272;
806   t277 = t7 * t195 * t129;
807   t279 = t110 * t128;
808   t281 = t200 * t201 * t279;
809   t284 = 0.1e1 / t90 / t81;
810   t285 = t14 * t284;
811   t286 = t47 * t8;
812   t288 = t7 * t285 * t286;
813   t293 = t4 * t6 * t247;
814   t294 = t16 * t74;
815   t295 = t133 * t116;
816   t297 = t293 * t294 * t295;
817   t300 = t7 * t254 * t134;
818   t302 = t52 * t284;
819   t303 = t74 * t8;
820   t305 = t7 * t302 * t303;
821   t308 = t6 * t52;
822   t309 = t4 * t308;
823   t310 = t16 * t139;
824   t311 = t152 * t116;
825   t313 = t309 * t310 * t311;
826   t316 = t7 * t121 * t153;
827   t318 = -t274 / 0.6e1 - t277 / 0.6e1 - 0.32320800000000000000e0 * t281 - t288 / 0.2e1 - t197 / 0.6e1 + t219 - 0.80802000000000000000e-1 * t221 - t297 / 0.6e1 - t300 / 0.6e1 + t305 / 0.2e1 - t256 / 0.6e1 + t266 - 0.32320800000000000000e0 * t313 - 0.80802000000000000000e-1 * t316;
828   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
829     v2rho2[1] = t9 * t318 - t114 - t120 - t132 - t137 - t156 - t184 - t186 - t89;
830 
831   t321 = t128 * t128;
832   t322 = t47 * t321;
833   t324 = t7 * t189 * t322;
834   t325 = t324 / 0.6e1;
835   t326 = t277 / 0.3e1;
836   t327 = t82 + t208;
837   t328 = t47 * t327;
838   t330 = t7 * t80 * t328;
839   t331 = t330 / 0.2e1;
840   t332 = t133 * t133;
841   t333 = t74 * t332;
842   t335 = t7 * t248 * t333;
843   t336 = t335 / 0.6e1;
844   t337 = t300 / 0.3e1;
845   t338 = t152 * t133;
846   t340 = t309 * t310 * t338;
847   t341 = 0.64641600000000000000e0 * t340;
848   t342 = -t327;
849   t343 = t74 * t342;
850   t345 = t7 * t115 * t343;
851   t346 = t345 / 0.2e1;
852   t347 = 0.16160400000000000000e0 * t316;
853   t349 = 0.1e1 / t138 / t71;
854   t350 = t152 * t152;
855   t351 = t349 * t350;
856   t353 = t7 * t54 * t351;
857   t354 = 0.48481200000000000000e0 * t353;
858   t356 = 0.1e1 / t57 / t64;
859   t362 = 0.1e1 / t56 / t64 / t140;
860   t366 = 0.110e3 / 0.2187e4 * t23 * sigma[2] * t356 + 0.11116869415428947191e-2 * t35 * t63 * t362;
861   t367 = t139 * t366;
862   t369 = t7 * t54 * t367;
863   t370 = 0.24240600000000000000e0 * t369;
864   t371 = -t325 - t326 - t331 + t219 - t336 - t337 - t341 - t346 + t266 - t347 + t354 - t370;
865   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
866     v2rho2[2] = -t131 - t184 - t136 - t186 - 0.48481200000000000000e0 * t155 + t9 * t371;
867 
868   t375 = 0.32320800000000000000e0 * t7 * t17 * t169;
869   t376 = t3 * t6;
870   t377 = t160 * t376;
871   t378 = t14 * t97;
872   t379 = t168 * t85;
873   t383 = t15 * t224;
874   t384 = t168 * t110;
875   t393 = -0.10e2 / 0.729e3 * t23 * t100 - 0.35105903417144043762e-3 * t35 * sigma[0] * t106;
876   t394 = t97 * t393;
877   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
878     v2rhosigma[0] = -t375 - 0.32320800000000000000e0 * t377 * t378 * t379 + 0.48481200000000000000e0 * t377 * t383 * t384 - 0.24240600000000000000e0 * t161 * t162 * t394;
879 
880   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
881     v2rhosigma[1] = 0.0e0;
882 
883   t400 = 0.32320800000000000000e0 * t7 * t54 * t180;
884   t401 = t52 * t139;
885   t402 = t179 * t116;
886   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
887     v2rhosigma[2] = -t400 - 0.32320800000000000000e0 * t377 * t401 * t402;
888 
889   t406 = t168 * t128;
890   t407 = t378 * t406;
891   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
892     v2rhosigma[3] = -t375 - 0.32320800000000000000e0 * t377 * t407;
893 
894   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
895     v2rhosigma[4] = 0.0e0;
896 
897   t410 = t179 * t133;
898   t411 = t401 * t410;
899   t414 = t53 * t349;
900   t415 = t179 * t152;
901   t416 = t414 * t415;
902   t424 = -0.10e2 / 0.729e3 * t23 * t142 - 0.35105903417144043762e-3 * t35 * sigma[2] * t148;
903   t426 = t173 * t139 * t424;
904   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
905     v2rhosigma[5] = -t400 - 0.32320800000000000000e0 * t377 * t411 + 0.48481200000000000000e0 * t377 * t416 - 0.24240600000000000000e0 * t161 * t426;
906 
907   t429 = t168 * t168;
908   t431 = t162 * t224 * t429;
909   t434 = t15 * t97;
910   t435 = t35 * t40;
911   t436 = t434 * t435;
912   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
913     v2sigma2[0] = 0.48481200000000000000e0 * t161 * t431 - 0.15956028044505410760e-4 * t377 * t436;
914 
915   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
916     v2sigma2[1] = 0.0e0;
917 
918   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
919     v2sigma2[2] = 0.0e0;
920 
921   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
922     v2sigma2[3] = 0.0e0;
923 
924   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
925     v2sigma2[4] = 0.0e0;
926 
927   t439 = t179 * t179;
928   t441 = t173 * t349 * t439;
929   t444 = t53 * t139;
930   t445 = t35 * t67;
931   t446 = t444 * t445;
932   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
933     v2sigma2[5] = 0.48481200000000000000e0 * t161 * t441 - 0.15956028044505410760e-4 * t377 * t446;
934 
935 #ifndef XC_DONT_COMPILE_KXC
936 
937   if(order < 3) return;
938 
939 
940   t454 = t218 / 0.4e1;
941   t458 = t265 / 0.4e1;
942   t459 = t85 * t209;
943   t461 = t270 * t271 * t459;
944   t462 = t461 / 0.2e1;
945   t463 = t110 * t190;
946   t465 = t270 * t201 * t463;
947   t466 = 0.32320800000000000000e0 * t465;
948   t467 = t91 * t97;
949   t469 = t200 * t467 * t202;
950   t471 = t241 * t85;
951   t473 = t200 * t201 * t471;
952   t474 = 0.96962400000000000000e0 * t473;
953   t475 = t110 * t209;
954   t477 = t200 * t201 * t475;
955   t478 = 0.96962400000000000000e0 * t477;
956   t479 = t4 * t162;
957   t480 = t16 * t224;
958   t481 = t110 * t241;
959   t483 = t479 * t480 * t481;
960   t484 = 0.14544360000000000000e1 * t483;
961   t485 = t116 * t258;
962   t487 = t293 * t294 * t485;
963   t488 = t487 / 0.2e1;
964   t489 = t225 * t85;
965   t491 = t200 * t480 * t489;
966   t492 = 0.19392480000000000000e1 * t491;
967   t494 = 0.1e1 / t187 / t13;
968   t495 = t494 * t16;
969   t496 = t190 * t85;
970   t497 = t47 * t496;
971   t499 = t7 * t495 * t497;
972   t500 = t499 / 0.9e1;
973   t501 = t188 * t91;
974   t503 = t7 * t501 * t191;
975   t505 = t14 * t215;
976   t507 = t7 * t505 * t86;
977   t508 = t507 / 0.3e1;
978   t510 = t7 * t92 * t226;
979   t512 = -t462 - t466 - 0.64641600000000000000e0 * t469 - t474 - t478 + t484 - t488 + t492 + t500 - t503 / 0.6e1 + t508 + 0.48481200000000000000e0 * t510;
980   t513 = t96 * t96;
981   t514 = 0.1e1 / t513;
982   t515 = t225 * t110;
983   t516 = t514 * t515;
984   t518 = t7 * t17 * t516;
985   t519 = 0.14544360000000000000e1 * t518;
986   t521 = 0.1e1 / t246 / t51;
987   t522 = t521 * t16;
988   t523 = t249 * t116;
989   t524 = t74 * t523;
990   t526 = t7 * t522 * t524;
991   t527 = t526 / 0.9e1;
992   t528 = t247 * t91;
993   t530 = t7 * t528 * t250;
994   t532 = t52 * t215;
995   t534 = t7 * t532 * t117;
996   t535 = t534 / 0.3e1;
997   t537 = t7 * t195 * t210;
998   t539 = t81 * t81;
999   t540 = 0.1e1 / t539;
1000   t541 = t8 * t540;
1001   t543 = 0.3e1 * t207 - 0.3e1 * t541;
1002   t544 = t47 * t543;
1003   t546 = t7 * t80 * t544;
1004   t547 = t546 / 0.2e1;
1005   t549 = t7 * t216 * t111;
1006   t550 = 0.16160400000000000000e0 * t549;
1007   t552 = t7 * t92 * t242;
1008   t555 = 0.1e1 / t26 / t38;
1009   t559 = t37 * t37;
1010   t561 = 0.1e1 / t25 / t559;
1011   t565 = -0.1540e4 / 0.6561e4 * t23 * sigma[0] * t555 - 0.81523709046478946067e-2 * t35 * t36 * t561;
1012   t566 = t97 * t565;
1013   t568 = t7 * t17 * t566;
1014   t569 = 0.24240600000000000000e0 * t568;
1015   t571 = t7 * t254 * t259;
1016   t573 = -t543;
1017   t574 = t74 * t573;
1018   t576 = t7 * t115 * t574;
1019   t577 = t576 / 0.2e1;
1020   t578 = t15 * t284;
1021   t580 = t7 * t578 * t47;
1022   t581 = 0.5e1 / 0.36e2 * t580;
1023   t582 = t53 * t284;
1024   t584 = t7 * t582 * t74;
1025   t585 = 0.5e1 / 0.36e2 * t584;
1026   t586 = -t519 + t527 - t530 / 0.6e1 + t535 - t537 / 0.2e1 - t547 + t550 - 0.24240600000000000000e0 * t552 - t569 - t571 / 0.2e1 - t577 - t581 - t585;
1027   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1028     v3rho3[0] = -t193 / 0.2e1 - t197 - 0.19392480000000000000e1 * t204 + 0.14544360000000000000e1 * t228 - t252 / 0.2e1 - t256 - 0.3e1 / 0.2e1 * t212 + t454 - 0.48481200000000000000e0 * t221 - 0.72721800000000000000e0 * t244 - 0.3e1 / 0.2e1 * t261 + t458 + t9 * (t512 + t586);
1029 
1030   t591 = t274 / 0.3e1;
1031   t593 = t297 / 0.3e1;
1032   t594 = -t194 - 0.2e1 / 0.3e1 * t197 - t205 + t229 - t253 - 0.2e1 / 0.3e1 * t256 - t591 - t326 - 0.64641600000000000000e0 * t281 - t288 - t593;
1033   t595 = 0.64641600000000000000e0 * t313;
1034   t598 = t4 * t6 * t521;
1035   t599 = t133 * t249;
1036   t601 = t598 * t294 * t599;
1037   t603 = t91 * t74;
1038   t605 = t293 * t603 * t295;
1039   t606 = t605 / 0.9e1;
1040   t607 = t152 * t249;
1041   t609 = t293 * t310 * t607;
1042   t611 = t91 * t139;
1043   t613 = t309 * t611 * t311;
1044   t614 = 0.21547200000000000000e0 * t613;
1045   t616 = t200 * t467 * t279;
1046   t617 = 0.21547200000000000000e0 * t616;
1047   t618 = t241 * t128;
1048   t620 = t200 * t201 * t618;
1049   t622 = t133 * t258;
1050   t624 = t293 * t294 * t622;
1051   t626 = t284 * t74;
1052   t627 = t8 * t116;
1053   t629 = t293 * t626 * t627;
1054   t631 = t152 * t258;
1055   t633 = t309 * t310 * t631;
1056   t635 = t128 * t209;
1057   t637 = t270 * t271 * t635;
1058   t640 = t4 * t6 * t494;
1059   t641 = t128 * t190;
1060   t643 = t640 * t271 * t641;
1061   t645 = t91 * t47;
1062   t647 = t270 * t645 * t272;
1063   t648 = t647 / 0.9e1;
1064   t649 = t284 * t47;
1065   t650 = t8 * t85;
1066   t652 = t270 * t649 * t650;
1067   t654 = t225 * t128;
1068   t656 = t200 * t480 * t654;
1069   t658 = t284 * t97;
1070   t659 = t110 * t8;
1071   t661 = t200 * t658 * t659;
1072   t665 = 0.1e1 / t90 / t206;
1073   t666 = t52 * t665;
1074   t668 = t7 * t666 * t303;
1075   t670 = t601 / 0.9e1 - t606 - 0.10773600000000000000e0 * t609 - t614 - t617 - 0.32320800000000000000e0 * t620 - t624 / 0.6e1 + t629 / 0.3e1 - 0.32320800000000000000e0 * t633 - t637 / 0.6e1 + t643 / 0.9e1 - t648 - t652 / 0.3e1 + 0.64641600000000000000e0 * t656 - 0.64641600000000000000e0 * t661 - 0.21547200000000000000e0 * t469 - 0.7e1 / 0.6e1 * t668;
1076   t672 = t7 * t505 * t129;
1077   t674 = t14 * t665;
1078   t676 = t7 * t674 * t286;
1079   t679 = t7 * t532 * t134;
1080   t682 = t7 * t263 * t153;
1081   t693 = t279 * t85;
1082   t695 = t270 * t201 * t693;
1083   t698 = t7 * t285 * t47;
1084   t701 = t7 * t302 * t74;
1085   t703 = t672 / 0.9e1 + 0.7e1 / 0.6e1 * t676 + t679 / 0.9e1 + 0.53868000000000000000e-1 * t682 - t503 / 0.18e2 + 0.2e1 / 0.9e1 * t507 + 0.16160400000000000000e0 * t510 - t530 / 0.18e2 + 0.2e1 / 0.9e1 * t534 - t537 / 0.6e1 + 0.10773600000000000000e0 * t549 - 0.80802000000000000000e-1 * t552 - t571 / 0.6e1 - 0.21547200000000000000e0 * t695 - t581 - t585 - t698 / 0.2e1 + t701 / 0.2e1;
1086   t706 = -t337 - t595 - t347 - t213 + t454 - 0.32320800000000000000e0 * t221 - t245 - t262 + t458 + t305 + t9 * (t670 + t703);
1087   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1088     v3rho3[1] = t594 + t706;
1089 
1090   t711 = -t198 - t257 - t591 - 0.2e1 / 0.3e1 * t277 - 0.64641600e0 * t281 - t288 - t593 - 0.2e1 / 0.3e1 * t300 - t595 - 0.32320800000000000000e0 * t316 - t325;
1091   t712 = t128 * t8;
1092   t714 = t270 * t649 * t712;
1093   t717 = t309 * t611 * t338;
1094   t719 = t321 * t85;
1095   t721 = t640 * t271 * t719;
1096   t723 = t110 * t321;
1097   t725 = t270 * t201 * t723;
1098   t727 = t327 * t85;
1099   t729 = t270 * t271 * t727;
1100   t731 = t110 * t327;
1101   t733 = t200 * t201 * t731;
1102   t735 = t332 * t116;
1103   t737 = t598 * t294 * t735;
1104   t739 = t133 * t8;
1105   t741 = t293 * t626 * t739;
1106   t743 = t284 * t139;
1107   t744 = t152 * t8;
1108   t746 = t309 * t743 * t744;
1109   t748 = t342 * t116;
1110   t750 = t293 * t294 * t748;
1111   t752 = t16 * t349;
1112   t753 = t350 * t116;
1113   t755 = t309 * t752 * t753;
1114   t757 = t366 * t116;
1115   t759 = t309 * t310 * t757;
1116   t762 = t7 * t501 * t322;
1117   t764 = -t714 / 0.3e1 - 0.21547200000000000000e0 * t717 + t721 / 0.9e1 - 0.10773600000000000000e0 * t725 - t729 / 0.6e1 - 0.32320800000000000000e0 * t733 + t737 / 0.9e1 + t741 / 0.3e1 + 0.64641600000000000000e0 * t746 - t750 / 0.6e1 + 0.64641600000000000000e0 * t755 - 0.32320800000000000000e0 * t759 - t606 - t614 - t617 - t648 - t762 / 0.18e2;
1118   t766 = t7 * t195 * t328;
1119   t769 = t7 * t528 * t333;
1120   t772 = t7 * t254 * t343;
1121   t775 = t7 * t121 * t351;
1122   t778 = t7 * t121 * t367;
1123   t781 = -t207 - 0.3e1 * t541;
1124   t782 = t47 * t781;
1125   t784 = t7 * t80 * t782;
1126   t786 = -t781;
1127   t787 = t74 * t786;
1128   t789 = t7 * t115 * t787;
1129   t799 = t338 * t116;
1130   t801 = t293 * t310 * t799;
1131   t803 = -t766 / 0.6e1 - t769 / 0.18e2 - t772 / 0.6e1 + 0.16160400000000000000e0 * t775 - 0.80802000000000000000e-1 * t778 - t784 / 0.2e1 - t789 / 0.2e1 + t668 / 0.3e1 + 0.2e1 / 0.9e1 * t672 - t676 / 0.3e1 + 0.2e1 / 0.9e1 * t679 + 0.10773600000000000000e0 * t682 + t507 / 0.9e1 + t534 / 0.9e1 + 0.53868000000000000000e-1 * t549 - 0.21547200000000000000e0 * t801 - t581 - t585;
1132   t806 = -t336 - t341 + t354 + t454 - t222 + t458 + t305 - t331 - t346 - t370 + t9 * (t764 + t803);
1133   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1134     v3rho3[2] = t711 + t806;
1135 
1136   t815 = t152 * t342;
1137   t817 = t309 * t310 * t815;
1138   t818 = 0.96962400000000000000e0 * t817;
1139   t819 = t4 * t173;
1140   t820 = t152 * t366;
1141   t822 = t819 * t752 * t820;
1142   t823 = 0.14544360000000000000e1 * t822;
1143   t824 = t350 * t133;
1144   t826 = t309 * t752 * t824;
1145   t827 = 0.19392480000000000000e1 * t826;
1146   t828 = t128 * t327;
1147   t830 = t270 * t271 * t828;
1148   t831 = t830 / 0.2e1;
1149   t832 = t133 * t342;
1150   t834 = t293 * t294 * t832;
1151   t835 = t834 / 0.2e1;
1152   t836 = t152 * t332;
1153   t838 = t293 * t310 * t836;
1154   t839 = 0.32320800000000000000e0 * t838;
1155   t840 = t366 * t133;
1156   t842 = t309 * t310 * t840;
1157   t843 = 0.96962400000000000000e0 * t842;
1158   t845 = t321 * t128;
1159   t846 = t47 * t845;
1160   t848 = t7 * t495 * t846;
1161   t849 = t848 / 0.9e1;
1162   t850 = t332 * t133;
1163   t851 = t74 * t850;
1164   t853 = t7 * t522 * t851;
1165   t854 = t853 / 0.9e1;
1166   t855 = t138 * t138;
1167   t856 = 0.1e1 / t855;
1168   t857 = t350 * t152;
1169   t858 = t856 * t857;
1170   t860 = t7 * t54 * t858;
1171   t861 = 0.14544360000000000000e1 * t860;
1172   t863 = -0.3e1 * t207 - 0.3e1 * t541;
1173   t864 = t47 * t863;
1174   t866 = t7 * t80 * t864;
1175   t867 = t866 / 0.2e1;
1176   t868 = -t818 + t823 + t827 - t831 - t835 - t839 - t843 - 0.64641600000000000000e0 * t717 + t849 + t854 - t861 - t867;
1177   t869 = -t863;
1178   t870 = t74 * t869;
1179   t872 = t7 * t115 * t870;
1180   t873 = t872 / 0.2e1;
1181   t875 = 0.1e1 / t57 / t65;
1182   t879 = t64 * t64;
1183   t881 = 0.1e1 / t56 / t879;
1184   t885 = -0.1540e4 / 0.6561e4 * t23 * sigma[2] * t875 - 0.81523709046478946067e-2 * t35 * t63 * t881;
1185   t886 = t139 * t885;
1186   t888 = t7 * t54 * t886;
1187   t889 = 0.24240600000000000000e0 * t888;
1188   t896 = t672 / 0.3e1;
1189   t897 = t679 / 0.3e1;
1190   t898 = 0.16160400000000000000e0 * t682;
1191   t899 = -t873 - t889 - t762 / 0.6e1 - t766 / 0.2e1 - t769 / 0.6e1 - t772 / 0.2e1 + 0.48481200000000000000e0 * t775 - 0.24240600000000000000e0 * t778 + t896 + t897 + t898 - t581 - t585;
1192   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1193     v3rho3[3] = -t277 - t300 - 0.48481200000000000000e0 * t316 - t324 / 0.2e1 - t335 / 0.2e1 - 0.19392480000000000000e1 * t340 + 0.14544360000000000000e1 * t353 + t454 + t458 - 0.3e1 / 0.2e1 * t330 - 0.3e1 / 0.2e1 * t345 - 0.72721800000000000000e0 * t369 + t9 * (t868 + t899);
1194 
1195   t903 = t200 * t201 * t379;
1196   t907 = 0.10773600000000000000e0 * t7 * t92 * t169;
1197   t909 = t479 * t480 * t384;
1198   t912 = t7 * t17 * t394;
1199   t914 = t188 * t97;
1200   t915 = t168 * t190;
1201   t919 = t14 * t224;
1202   t920 = t379 * t110;
1203   t924 = t393 * t85;
1204   t928 = t168 * t209;
1205   t932 = t15 * t514;
1206   t933 = t168 * t225;
1207   t937 = t393 * t110;
1208   t941 = t168 * t241;
1209   t950 = 0.110e3 / 0.2187e4 * t23 * t231 + 0.22233738830857894383e-2 * t35 * sigma[0] * t237;
1210   t951 = t97 * t950;
1211   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1212     v3rho2sigma[0] = -0.86188800000000000000e0 * t903 - t907 + 0.12928320000000000000e1 * t909 - 0.64641600000000000000e0 * t912 - 0.10773600000000000000e0 * t377 * t914 * t915 + 0.12928320000000000000e1 * t377 * t919 * t920 - 0.64641600000000000000e0 * t377 * t378 * t924 - 0.32320800000000000000e0 * t377 * t378 * t928 - 0.14544360000000000000e1 * t377 * t932 * t933 + 0.96962400000000000000e0 * t377 * t383 * t937 + 0.48481200000000000000e0 * t377 * t383 * t941 - 0.24240600000000000000e0 * t161 * t162 * t951;
1213 
1214   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1215     v3rho2sigma[1] = 0.0e0;
1216 
1217   t956 = t309 * t310 * t402;
1218   t960 = 0.10773600000000000000e0 * t7 * t121 * t180;
1219   t961 = t247 * t139;
1220   t962 = t179 * t249;
1221   t966 = t179 * t258;
1222   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1223     v3rho2sigma[2] = -0.86188800000000000000e0 * t956 - t960 - 0.10773600000000000000e0 * t377 * t961 * t962 - 0.32320800000000000000e0 * t377 * t401 * t966;
1224 
1225   t973 = t16 * t1;
1226   t974 = t973 * t376;
1227   t975 = t974 * t407;
1228   t978 = t914 * t406 * t85;
1229   t982 = t919 * t406 * t110;
1230   t985 = t393 * t128;
1231   t986 = t378 * t985;
1232   t989 = t215 * t1;
1233   t990 = t989 * t376;
1234   t991 = t168 * t8;
1235   t992 = t378 * t991;
1236   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1237     v3rho2sigma[3] = -0.43094400000000000000e0 * t903 - t907 + 0.64641600000000000000e0 * t909 - 0.32320800000000000000e0 * t912 - 0.43094400000000000000e0 * t975 - 0.10773600000000000000e0 * t377 * t978 + 0.64641600000000000000e0 * t377 * t982 - 0.32320800000000000000e0 * t377 * t986 - 0.32320800000000000000e0 * t990 * t992;
1238 
1239   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1240     v3rho2sigma[4] = 0.0e0;
1241 
1242   t996 = t974 * t411;
1243   t999 = t961 * t410 * t116;
1244   t1002 = t179 * t8;
1245   t1003 = t401 * t1002;
1246   t1006 = t974 * t416;
1247   t1008 = t52 * t349;
1248   t1010 = t1008 * t415 * t116;
1249   t1013 = t973 * t3;
1250   t1014 = t1013 * t426;
1251   t1017 = t401 * t424 * t116;
1252   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1253     v3rho2sigma[5] = -0.43094400000000000000e0 * t956 - t960 - 0.43094400000000000000e0 * t996 - 0.10773600000000000000e0 * t377 * t999 + 0.32320800000000000000e0 * t990 * t1003 + 0.64641600000000000000e0 * t1006 + 0.64641600000000000000e0 * t377 * t1010 - 0.32320800000000000000e0 * t1014 - 0.32320800000000000000e0 * t377 * t1017;
1254 
1255   t1021 = t168 * t321;
1256   t1022 = t914 * t1021;
1257   t1025 = t168 * t327;
1258   t1026 = t378 * t1025;
1259   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1260     v3rho2sigma[6] = -0.86188800000000000000e0 * t975 - t907 - 0.10773600000000000000e0 * t377 * t1022 - 0.32320800000000000000e0 * t377 * t1026;
1261 
1262   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1263     v3rho2sigma[7] = 0.0e0;
1264 
1265   t1032 = t179 * t332;
1266   t1033 = t961 * t1032;
1267   t1037 = t1008 * t410 * t152;
1268   t1040 = t424 * t133;
1269   t1041 = t401 * t1040;
1270   t1044 = t179 * t342;
1271   t1045 = t401 * t1044;
1272   t1048 = t53 * t856;
1273   t1049 = t179 * t350;
1274   t1050 = t1048 * t1049;
1275   t1053 = t424 * t152;
1276   t1054 = t414 * t1053;
1277   t1057 = t179 * t366;
1278   t1058 = t414 * t1057;
1279   t1066 = 0.110e3 / 0.2187e4 * t23 * t356 + 0.22233738830857894383e-2 * t35 * sigma[2] * t362;
1280   t1068 = t173 * t139 * t1066;
1281   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1282     v3rho2sigma[8] = -0.86188800000000000000e0 * t996 - t960 + 0.12928320000000000000e1 * t1006 - 0.64641600000000000000e0 * t1014 - 0.10773600000000000000e0 * t377 * t1033 + 0.12928320000000000000e1 * t377 * t1037 - 0.64641600000000000000e0 * t377 * t1041 - 0.32320800000000000000e0 * t377 * t1045 - 0.14544360000000000000e1 * t377 * t1050 + 0.96962400000000000000e0 * t377 * t1054 + 0.48481200000000000000e0 * t377 * t1058 - 0.24240600000000000000e0 * t161 * t1068;
1283 
1284   t1072 = 0.64641600000000000000e0 * t1013 * t431;
1285   t1073 = t429 * t85;
1286   t1074 = t919 * t1073;
1287   t1078 = t932 * t429 * t110;
1288   t1082 = t383 * t168 * t393;
1289   t1086 = 0.21274704059340547680e-4 * t974 * t436;
1290   t1087 = t376 * t14;
1291   t1088 = t160 * t1087;
1292   t1089 = t97 * t32;
1293   t1090 = t34 * t40;
1294   t1091 = t1090 * t85;
1295   t1092 = t1089 * t1091;
1296   t1095 = t376 * t15;
1297   t1096 = t160 * t1095;
1298   t1097 = t224 * t32;
1299   t1098 = t1090 * t110;
1300   t1099 = t1097 * t1098;
1301   t1102 = t35 * t106;
1302   t1103 = t434 * t1102;
1303   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1304     v3rhosigma2[0] = t1072 + 0.64641600000000000000e0 * t377 * t1074 - 0.14544360000000000000e1 * t377 * t1078 + 0.96962400000000000000e0 * t377 * t1082 - t1086 - 0.21274704059340547680e-4 * t1088 * t1092 + 0.31912056089010821520e-4 * t1096 * t1099 + 0.85098816237362190720e-4 * t377 * t1103;
1305 
1306   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1307     v3rhosigma2[1] = 0.0e0;
1308 
1309   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1310     v3rhosigma2[2] = 0.0e0;
1311 
1312   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1313     v3rhosigma2[3] = 0.0e0;
1314 
1315   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1316     v3rhosigma2[4] = 0.0e0;
1317 
1318   t1107 = 0.64641600000000000000e0 * t1013 * t441;
1319   t1109 = t1008 * t439 * t116;
1320   t1113 = 0.21274704059340547680e-4 * t974 * t446;
1321   t1114 = t376 * t52;
1322   t1115 = t160 * t1114;
1323   t1116 = t139 * t32;
1324   t1117 = t34 * t67;
1325   t1118 = t1117 * t116;
1326   t1119 = t1116 * t1118;
1327   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1328     v3rhosigma2[5] = t1107 + 0.64641600000000000000e0 * t377 * t1109 - t1113 - 0.21274704059340547680e-4 * t1115 * t1119;
1329 
1330   t1122 = t429 * t128;
1331   t1123 = t919 * t1122;
1332   t1126 = t1090 * t128;
1333   t1127 = t1089 * t1126;
1334   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1335     v3rhosigma2[6] = t1072 + 0.64641600000000000000e0 * t377 * t1123 - t1086 - 0.21274704059340547680e-4 * t1088 * t1127;
1336 
1337   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1338     v3rhosigma2[7] = 0.0e0;
1339 
1340   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1341     v3rhosigma2[8] = 0.0e0;
1342 
1343   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1344     v3rhosigma2[9] = 0.0e0;
1345 
1346   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1347     v3rhosigma2[10] = 0.0e0;
1348 
1349   t1130 = t439 * t133;
1350   t1131 = t1008 * t1130;
1351   t1134 = t439 * t152;
1352   t1135 = t1048 * t1134;
1353   t1138 = t179 * t424;
1354   t1139 = t414 * t1138;
1355   t1142 = t1117 * t133;
1356   t1143 = t1116 * t1142;
1357   t1146 = t376 * t53;
1358   t1147 = t160 * t1146;
1359   t1148 = t349 * t32;
1360   t1149 = t1117 * t152;
1361   t1150 = t1148 * t1149;
1362   t1153 = t35 * t148;
1363   t1154 = t444 * t1153;
1364   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1365     v3rhosigma2[11] = t1107 + 0.64641600000000000000e0 * t377 * t1131 - 0.14544360000000000000e1 * t377 * t1135 + 0.96962400000000000000e0 * t377 * t1139 - t1113 - 0.21274704059340547680e-4 * t1115 * t1143 + 0.31912056089010821520e-4 * t1147 * t1150 + 0.85098816237362190720e-4 * t377 * t1154;
1366 
1367   t1157 = t429 * t168;
1368   t1159 = t162 * t514 * t1157;
1369   t1162 = t224 * t168;
1370   t1163 = t1162 * t435;
1371   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1372     v3sigma3[0] = -0.14544360000000000000e1 * t161 * t1159 + 0.95736168267032464562e-4 * t1096 * t1163;
1373 
1374   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1375     v3sigma3[1] = 0.0e0;
1376 
1377   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1378     v3sigma3[2] = 0.0e0;
1379 
1380   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1381     v3sigma3[3] = 0.0e0;
1382 
1383   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1384     v3sigma3[4] = 0.0e0;
1385 
1386   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1387     v3sigma3[5] = 0.0e0;
1388 
1389   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1390     v3sigma3[6] = 0.0e0;
1391 
1392   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1393     v3sigma3[7] = 0.0e0;
1394 
1395   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1396     v3sigma3[8] = 0.0e0;
1397 
1398   t1166 = t439 * t179;
1399   t1168 = t173 * t856 * t1166;
1400   t1171 = t349 * t179;
1401   t1172 = t1171 * t445;
1402   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1403     v3sigma3[9] = -0.14544360000000000000e1 * t161 * t1168 + 0.95736168267032464562e-4 * t1147 * t1172;
1404 
1405 #ifndef XC_DONT_COMPILE_LXC
1406 
1407   if(order < 4) return;
1408 
1409 
1410   t1183 = t91 * t224;
1411   t1185 = t479 * t1183 * t481;
1412   t1187 = t16 * t514;
1413   t1205 = t293 * t603 * t485;
1414   t1216 = t200 * t1183 * t489;
1415   t1222 = -0.19392480000000000000e1 * t200 * t201 * t241 * t209 - 0.12928320000000000000e1 * t200 * t201 * t110 * t543 + 0.19392480000000000000e1 * t1185 - 0.77569920000000000000e1 * t200 * t1187 * t515 * t85 - 0.87266160000000000000e1 * t479 * t1187 * t225 * t241 + 0.19392480000000000000e1 * t479 * t480 * t110 * t565 + 0.2e1 / 0.3e1 * t598 * t294 * t249 * t258 - 0.2e1 / 0.3e1 * t1205 - 0.2e1 / 0.3e1 * t293 * t294 * t116 * t573 + 0.12928320000000000000e1 * t270 * t480 * t225 * t190 + 0.25856640000000000000e1 * t1216 + 0.38784960000000000000e1 * t200 * t480 * t225 * t209;
1416   t1232 = t270 * t645 * t459;
1417   t1235 = t270 * t467 * t463;
1418   t1245 = t215 * t97;
1419   t1247 = t200 * t1245 * t202;
1420   t1250 = t200 * t467 * t471;
1421   t1253 = t200 * t467 * t475;
1422   t1259 = t539 * t9;
1423   t1261 = t8 / t1259;
1424   t1263 = 0.12e2 * t540 - 0.12e2 * t1261;
1425   t1268 = t13 * t13;
1426   t1270 = 0.1e1 / t187 / t1268;
1427   t1271 = t1270 * t16;
1428   t1272 = t190 * t190;
1429   t1277 = t494 * t91;
1430   t1279 = t7 * t1277 * t497;
1431   t1281 = 0.2e1 / 0.3e1 * t640 * t271 * t190 * t209 + 0.28729600000000000000e0 * t640 * t201 * t110 * t496 - 0.2e1 / 0.3e1 * t1232 - 0.43094400000000000000e0 * t1235 - 0.2e1 / 0.3e1 * t270 * t271 * t85 * t543 - 0.64641600000000000000e0 * t270 * t201 * t241 * t190 + 0.86188800000000000000e0 * t1247 - 0.12928320000000000000e1 * t1250 - 0.12928320000000000000e1 * t1253 - 0.12928320000000000000e1 * t200 * t201 * t565 * t85 - t7 * t115 * t74 * t1263 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1271 * t47 * t1272 + 0.4e1 / 0.27e2 * t1279;
1432   t1283 = t188 * t215;
1433   t1285 = t7 * t1283 * t191;
1434   t1287 = t209 * t209;
1435   t1293 = t7 * t285 * t86;
1436   t1296 = t7 * t216 * t226;
1437   t1299 = t7 * t92 * t516;
1438   t1302 = 0.1e1 / t513 / t44;
1439   t1303 = t225 * t225;
1440   t1308 = t241 * t241;
1441   t1313 = t51 * t51;
1442   t1315 = 0.1e1 / t246 / t1313;
1443   t1316 = t1315 * t16;
1444   t1317 = t249 * t249;
1445   t1322 = t521 * t91;
1446   t1324 = t7 * t1322 * t524;
1447   t1326 = t247 * t215;
1448   t1328 = t7 * t1326 * t250;
1449   t1330 = t258 * t258;
1450   t1336 = t7 * t302 * t117;
1451   t1338 = 0.2e1 / 0.9e1 * t1285 - t7 * t189 * t47 * t1287 / 0.2e1 - 0.20e2 / 0.27e2 * t1293 - 0.64641600000000000000e0 * t1296 - 0.19392480000000000000e1 * t1299 + 0.58177440000000000000e1 * t7 * t17 * t1302 * t1303 + 0.14544360000000000000e1 * t7 * t17 * t224 * t1308 - 0.5e1 / 0.27e2 * t7 * t1316 * t74 * t1317 + 0.4e1 / 0.27e2 * t1324 + 0.2e1 / 0.9e1 * t1328 - t7 * t248 * t74 * t1330 / 0.2e1 - 0.20e2 / 0.27e2 * t1336;
1452   t1340 = t7 * t505 * t210;
1453   t1343 = t7 * t532 * t259;
1454   t1346 = t7 * t195 * t544;
1455   t1354 = t7 * t578 * t111;
1456   t1357 = t7 * t216 * t242;
1457   t1360 = t7 * t92 * t566;
1458   t1379 = t7 * t254 * t574;
1459   t1392 = 0.10e2 / 0.27e2 * t7 * t15 * t665 * t47;
1460   t1396 = 0.10e2 / 0.27e2 * t7 * t53 * t665 * t74;
1461   t1397 = 0.2e1 / 0.3e1 * t1340 + 0.2e1 / 0.3e1 * t1343 - 0.2e1 / 0.3e1 * t1346 + t7 * t80 * t47 * t1263 / 0.2e1 - 0.35912000000000000000e0 * t1354 + 0.32320800000000000000e0 * t1357 - 0.32320800000000000000e0 * t1360 - 0.24240600000000000000e0 * t7 * t17 * t97 * (0.26180e5 / 0.19683e5 * t23 * sigma[0] / t26 / t104 + 0.67936424205399121722e-1 * t35 * t36 / t25 / t559 / rho[0]) - 0.2e1 / 0.3e1 * t1379 - 0.12928320000000000000e1 * t270 * t201 * t202 * t209 + 0.77569920000000000000e1 * t200 * t480 * t471 * t110 + t1392 + t1396;
1462   t1413 = t9 * (t1222 + t1281 + t1338 + t1397) - 0.2e1 * t461 - 0.12928320000000000000e1 * t465 - 0.25856640000000000000e1 * t469 - 0.38784960000000000000e1 * t473 - 0.38784960000000000000e1 * t477 + 0.58177440000000000000e1 * t483 - 0.2e1 * t487 + 0.77569920000000000000e1 * t491 + 0.4e1 / 0.9e1 * t499 - 0.2e1 / 0.3e1 * t503 + 0.4e1 / 0.3e1 * t507 + 0.19392480000000000000e1 * t510;
1463   t1425 = 0.5e1 / 0.9e1 * t580;
1464   t1426 = 0.5e1 / 0.9e1 * t584;
1465   t1427 = -0.58177440000000000000e1 * t518 + 0.4e1 / 0.9e1 * t526 - 0.2e1 / 0.3e1 * t530 + 0.4e1 / 0.3e1 * t534 - 0.2e1 * t537 - 0.2e1 * t546 + 0.64641600000000000000e0 * t549 - 0.96962400000000000000e0 * t552 - 0.96962400000000000000e0 * t568 - 0.2e1 * t571 - 0.2e1 * t576 - t1425 - t1426;
1466   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1467     v4rho4[0] = t1413 + t1427;
1468 
1469   t1429 = t293 * t603 * t622;
1470   t1432 = t309 * t611 * t631;
1471   t1435 = t200 * t1245 * t279;
1472   t1436 = 0.21547200000000000000e0 * t1435;
1473   t1438 = t200 * t467 * t618;
1474   t1457 = t270 * t645 * t635;
1475   t1465 = t293 * t215 * t74 * t295;
1476   t1466 = t1465 / 0.9e1;
1477   t1472 = t293 * t611 * t607;
1478   t1474 = t215 * t139;
1479   t1476 = t309 * t1474 * t311;
1480   t1477 = 0.21547200000000000000e0 * t1476;
1481   t1479 = t200 * t1183 * t654;
1482   t1492 = t200 * t665 * t97 * t659;
1483   t1494 = -t1429 / 0.6e1 - 0.32320800000000000000e0 * t1432 + t1436 - 0.32320800000000000000e0 * t1438 - 0.32320800000000000000e0 * t200 * t201 * t565 * t128 - t293 * t294 * t133 * t573 / 0.6e1 + t293 * t626 * t8 * t258 / 0.2e1 - 0.32320800000000000000e0 * t309 * t310 * t152 * t573 - t1457 / 0.6e1 - t270 * t649 * t8 * t209 / 0.2e1 + t1466 + 0.71824000000000000000e-1 * t598 * t310 * t152 * t523 - 0.10773600000000000000e0 * t1472 + t1477 + 0.64641600000000000000e0 * t1479 - 0.19392480000000000000e1 * t200 * t1187 * t515 * t128 + 0.19392480000000000000e1 * t200 * t284 * t224 * t225 * t8 + 0.22624560000000000000e1 * t1492;
1484   t1496 = t4 * t6 * t1315;
1485   t1502 = t598 * t603 * t599;
1486   t1505 = t4 * t6 * t1270;
1487   t1511 = t640 * t645 * t641;
1488   t1519 = t270 * t215 * t47 * t272;
1489   t1520 = t1519 / 0.9e1;
1490   t1521 = t665 * t47;
1491   t1523 = t270 * t1521 * t650;
1492   t1529 = t665 * t74;
1493   t1531 = t293 * t1529 * t627;
1494   t1549 = -0.5e1 / 0.27e2 * t1496 * t294 * t133 * t523 + t1502 / 0.9e1 - 0.5e1 / 0.27e2 * t1505 * t271 * t128 * t496 + t1511 / 0.9e1 + t640 * t649 * t8 * t190 / 0.3e1 + t1520 + 0.7e1 / 0.6e1 * t1523 - t270 * t271 * t128 * t543 / 0.6e1 - 0.7e1 / 0.6e1 * t1531 - 0.96962400000000000000e0 * t200 * t658 * t241 * t8 - t598 * t626 * t8 * t249 / 0.3e1 + 0.48481200000000000000e0 * t1185 - t1205 / 0.6e1 + 0.64641600000000000000e0 * t1216 - t1232 / 0.6e1 - 0.10773600000000000000e0 * t1235 + 0.43094400000000000000e0 * t1247 - 0.32320800000000000000e0 * t1250 - 0.32320800000000000000e0 * t1253;
1495   t1552 = 0.1e1 / t90 / t539;
1496   t1555 = t7 * t52 * t1552 * t303;
1497   t1558 = t7 * t285 * t129;
1498   t1562 = t7 * t14 * t1552 * t286;
1499   t1565 = t7 * t302 * t134;
1500   t1568 = t7 * t582 * t153;
1501   t1570 = t188 * t284;
1502   t1577 = t247 * t284;
1503   t1591 = 0.25e2 / 0.6e1 * t1555 - 0.5e1 / 0.27e2 * t1558 - 0.25e2 / 0.6e1 * t1562 - 0.5e1 / 0.27e2 * t1565 - 0.89780000000000000000e-1 * t1568 - t7 * t1570 * t86 / 0.2e1 - 0.96962400000000000000e0 * t7 * t285 * t111 + t7 * t1577 * t117 / 0.2e1 + t1279 / 0.27e2 + t1285 / 0.9e1 - 0.5e1 / 0.9e1 * t1293 - 0.32320800000000000000e0 * t1296 - 0.48481200000000000000e0 * t1299 + t1324 / 0.27e2 + t1328 / 0.9e1 - 0.5e1 / 0.9e1 * t1336 + t1340 / 0.3e1 + t1343 / 0.3e1;
1504   t1618 = t270 * t467 * t693;
1505   t1637 = t7 * t674 * t47;
1506   t1640 = t7 * t666 * t74;
1507   t1642 = -t1346 / 0.6e1 - 0.26934000000000000000e0 * t1354 + 0.16160400000000000000e0 * t1357 - 0.80802000000000000000e-1 * t1360 - t1379 / 0.6e1 - 0.32320800000000000000e0 * t270 * t201 * t279 * t209 + 0.64641600000000000000e0 * t270 * t480 * t654 * t85 - 0.64641600000000000000e0 * t270 * t658 * t659 * t85 + t598 * t294 * t295 * t258 / 0.3e1 + 0.21547200000000000000e0 * t640 * t201 * t279 * t190 - 0.21547200000000000000e0 * t1618 - 0.32320800000000000000e0 * t270 * t201 * t618 * t85 + t640 * t271 * t635 * t85 / 0.3e1 - 0.32320800000000000000e0 * t293 * t310 * t311 * t258 + 0.19392480000000000000e1 * t200 * t480 * t279 * t241 + 0.5e1 / 0.2e1 * t1637 - 0.5e1 / 0.2e1 * t1640 + t1392 + t1396;
1508   t1647 = t605 / 0.3e1;
1509   t1649 = 0.64641600000000000000e0 * t613;
1510   t1656 = t9 * (t1494 + t1549 + t1591 + t1642) + t601 / 0.3e1 - t1647 - 0.32320800000000000000e0 * t609 - t1649 - 0.64641600000000000000e0 * t616 - 0.96962400000000000000e0 * t620 - t624 / 0.2e1 + t629 - 0.96962400000000000000e0 * t633 - t637 / 0.2e1 + t643 / 0.3e1;
1511   t1657 = t647 / 0.3e1;
1512   t1661 = -t1657 - t652 + 0.19392480000000000000e1 * t656 - 0.19392480000000000000e1 * t661 - t462 - t466 - 0.12928320000000000000e1 * t469 - t474 - t478 + t484 - t488 + t492;
1513   t1668 = -0.7e1 / 0.2e1 * t668 + t896 + 0.7e1 / 0.2e1 * t676 + t897 + t898 + t500 - t503 / 0.3e1 + t507 + 0.96962400000000000000e0 * t510 - t519 + t527 - t530 / 0.3e1;
1514   t1674 = t534 - t537 - t547 + 0.48481200000000000000e0 * t549 - 0.48481200000000000000e0 * t552 - t569 - t571 - t577 - 0.64641600000000000000e0 * t695 - t1425 - t1426 - 0.3e1 / 0.2e1 * t698 + 0.3e1 / 0.2e1 * t701;
1515   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1516     v4rho4[1] = t1656 + t1661 + t1668 + t1674;
1517 
1518   t1694 = t309 * t665 * t139 * t744;
1519   t1701 = t270 * t467 * t723;
1520   t1716 = t200 * t467 * t731;
1521   t1718 = -t293 * t294 * t786 * t116 / 0.3e1 - t293 * t294 * t342 * t258 / 0.6e1 + 0.64641600000000000000e0 * t309 * t752 * t350 * t258 - 0.32320800000000000000e0 * t309 * t310 * t366 * t258 - 0.15083040000000000000e1 * t1694 + t640 * t271 * t321 * t209 / 0.9e1 - 0.71824000000000000000e-1 * t1701 - 0.10773600000000000000e0 * t270 * t201 * t241 * t321 - t270 * t271 * t781 * t85 / 0.3e1 - t270 * t271 * t327 * t209 / 0.6e1 - 0.21547200000000000000e0 * t1716;
1522   t1732 = t293 * t1529 * t739;
1523   t1738 = t91 * t349;
1524   t1740 = t309 * t1738 * t753;
1525   t1747 = t309 * t611 * t757;
1526   t1750 = t598 * t603 * t735;
1527   t1757 = t293 * t603 * t748;
1528   t1759 = -0.32320800000000000000e0 * t200 * t201 * t241 * t327 - 0.64641600000000000000e0 * t200 * t201 * t110 * t781 + t598 * t294 * t332 * t258 / 0.9e1 - 0.7e1 / 0.9e1 * t1732 + 0.21547200000000000000e0 * t293 * t752 * t350 * t249 + 0.43094400000000000000e0 * t1740 - 0.10773600000000000000e0 * t293 * t310 * t366 * t249 - 0.21547200000000000000e0 * t1747 + 0.2e1 / 0.27e2 * t1750 + t598 * t294 * t342 * t249 / 0.9e1 - t1757 / 0.9e1;
1529   t1770 = t270 * t645 * t727;
1530   t1785 = t640 * t645 * t719;
1531   t1788 = t270 * t1521 * t712;
1532   t1791 = t309 * t1474 * t338;
1533   t1795 = 0.21547200000000000000e0 * t270 * t480 * t225 * t321 + t640 * t271 * t327 * t190 / 0.9e1 - t1770 / 0.9e1 + 0.64641600000000000000e0 * t200 * t480 * t225 * t327 - 0.5e1 / 0.27e2 * t1496 * t294 * t332 * t249 - 0.5e1 / 0.27e2 * t1505 * t271 * t321 * t190 + 0.2e1 / 0.27e2 * t1785 + 0.7e1 / 0.9e1 * t1788 + 0.14364800000000000000e0 * t1791 - t1429 / 0.9e1 - 0.21547200000000000000e0 * t1432;
1534   t1807 = 0.28729600000000000000e0 * t1435 - 0.21547200000000000000e0 * t1438 - t1457 / 0.9e1 + 0.4e1 / 0.27e2 * t1465 - 0.71824000000000000000e-1 * t1472 + 0.28729600000000000000e0 * t1476 + 0.43094400000000000000e0 * t1479 - 0.43094400000000000000e0 * t1492 + 0.2e1 / 0.27e2 * t1502 + 0.2e1 / 0.27e2 * t1511 + 0.4e1 / 0.27e2 * t1519;
1535   t1814 = 0.1e1 / t90 / t1259;
1536   t1816 = t8 * t8;
1537   t1822 = t7 * t254 * t787;
1538   t1833 = t7 * t1283 * t322;
1539   t1836 = t7 * t505 * t328;
1540   t1839 = t7 * t1326 * t333;
1541   t1844 = -0.2e1 / 0.9e1 * t1523 + 0.2e1 / 0.9e1 * t1531 + 0.14364800000000000000e0 * t1247 - t7 * t247 * t1814 * t74 * t1816 / 0.3e1 - t1822 / 0.3e1 - t7 * t188 * t1814 * t47 * t1816 / 0.3e1 - t7 * t1570 * t129 / 0.3e1 + t1833 / 0.27e2 + t1836 / 0.9e1 + t1839 / 0.27e2 + t7 * t1577 * t134 / 0.3e1;
1542   t1849 = t7 * t532 * t343;
1543   t1852 = t7 * t263 * t351;
1544   t1855 = t7 * t263 * t367;
1545   t1863 = t7 * t195 * t782;
1546   t1866 = 0.64641600000000000000e0 * t7 * t302 * t153 + t1849 / 0.9e1 - 0.10773600000000000000e0 * t1852 + 0.53868000000000000000e-1 * t1855 + 0.41e2 / 0.9e1 * t1555 - 0.10e2 / 0.27e2 * t1558 - 0.41e2 / 0.9e1 * t1562 - 0.10e2 / 0.27e2 * t1565 - 0.17956000000000000000e0 * t1568 - t1863 / 0.3e1 + t1285 / 0.27e2;
1547   t1881 = t293 * t611 * t799;
1548   t1887 = -0.10e2 / 0.27e2 * t1293 - 0.10773600000000000000e0 * t1296 + t1328 / 0.27e2 - 0.10e2 / 0.27e2 * t1336 + t1340 / 0.9e1 + t1343 / 0.9e1 - 0.17956000000000000000e0 * t1354 + 0.53868000000000000000e-1 * t1357 + 0.14364800000000000000e0 * t598 * t310 * t338 * t249 - 0.14364800000000000000e0 * t1881 - 0.4e1 / 0.9e1 * t598 * t626 * t295 * t8;
1549   t1915 = -0.21547200000000000000e0 * t270 * t201 * t731 * t85 + 0.14364800000000000000e0 * t640 * t201 * t723 * t85 + 0.4e1 / 0.9e1 * t640 * t649 * t712 * t85 - 0.43094400000000000000e0 * t270 * t658 * t279 * t8 - 0.21547200000000000000e0 * t293 * t310 * t338 * t258 + 0.43094400000000000000e0 * t293 * t743 * t744 * t116 - 0.14364800000000000000e0 * t1618 - t1637 / 0.3e1 + t1640 / 0.3e1 + t1392 + t1396;
1550   t1933 = t9 * (t1718 + t1759 + t1795 + t1807 + t1844 + t1866 + t1887 + t1915) - 0.2e1 / 0.3e1 * t714 - 0.43094400000000000000e0 * t717 + 0.2e1 / 0.9e1 * t721 - 0.21547200000000000000e0 * t725 - t729 / 0.3e1 - 0.64641600000000000000e0 * t733 + 0.2e1 / 0.9e1 * t737 + 0.2e1 / 0.3e1 * t741 + 0.12928320000000000000e1 * t746 - t750 / 0.3e1 + 0.12928320000000000000e1 * t755 - 0.64641600000000000000e0 * t759 + 0.2e1 / 0.9e1 * t601;
1551   t1948 = -0.4e1 / 0.9e1 * t605 - 0.21547200000000000000e0 * t609 - 0.86188800000000000000e0 * t613 - 0.86188800000000000000e0 * t616 - 0.64641600e0 * t620 - t624 / 0.3e1 + 0.2e1 / 0.3e1 * t629 - 0.64641600000000000000e0 * t633 - t637 / 0.3e1 + 0.2e1 / 0.9e1 * t643 - 0.4e1 / 0.9e1 * t647 - 0.2e1 / 0.3e1 * t652 + 0.129283200e1 * t656 - 0.129283200e1 * t661;
1552   t1962 = -0.43094400000000000000e0 * t469 - t762 / 0.9e1 - t766 / 0.3e1 - t769 / 0.9e1 - t772 / 0.3e1 + 0.32320800000000000000e0 * t775 - 0.16160400000000000000e0 * t778 - t784 - t789 - 0.5e1 / 0.3e1 * t668 + 0.2e1 / 0.3e1 * t672 + 0.5e1 / 0.3e1 * t676 + 0.2e1 / 0.3e1 * t679 + 0.32320800000000000000e0 * t682;
1553   t1974 = -t503 / 0.9e1 + 0.2e1 / 0.3e1 * t507 + 0.32320800000000000000e0 * t510 - t530 / 0.9e1 + 0.2e1 / 0.3e1 * t534 - t537 / 0.3e1 + 0.32320800000000000000e0 * t549 - 0.16160400000000000000e0 * t552 - t571 / 0.3e1 - 0.43094400000000000000e0 * t801 - 0.43094400000000000000e0 * t695 - t1425 - t1426 - t698 + t701;
1554   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1555     v4rho4[2] = t1933 + t1948 + t1962 + t1974;
1556 
1557   t1977 = t7 * t121 * t858;
1558   t1980 = t7 * t195 * t864;
1559   t1983 = t7 * t254 * t870;
1560   t1986 = t7 * t121 * t886;
1561   t1989 = t7 * t1277 * t846;
1562   t1992 = t7 * t1322 * t851;
1563   t1996 = 0.6e1 * t540 + 0.12e2 * t1261;
1564   t2016 = t309 * t1738 * t824;
1565   t2019 = t819 * t1738 * t820;
1566   t2022 = t293 * t611 * t836;
1567   t2025 = t309 * t611 * t840;
1568   t2028 = t309 * t611 * t815;
1569   t2035 = t270 * t645 * t828;
1570   t2042 = t293 * t603 * t832;
1571   t2044 = -0.48481200000000000000e0 * t1977 - t1980 / 0.6e1 - t1983 / 0.6e1 - 0.80802000000000000000e-1 * t1986 + t1989 / 0.27e2 + t1992 / 0.27e2 - t7 * t80 * t47 * t1996 / 0.2e1 + t7 * t115 * t74 * t1996 / 0.2e1 - 0.19392480000000000000e1 * t309 * t284 * t349 * t350 * t8 + 0.96962400000000000000e0 * t309 * t743 * t366 * t8 + 0.64641600000000000000e0 * t2016 + 0.48481200000000000000e0 * t2019 - 0.10773600000000000000e0 * t2022 - 0.32320800000000000000e0 * t2025 - 0.32320800000000000000e0 * t2028 + t640 * t649 * t321 * t8 / 0.3e1 - t2035 / 0.6e1 - t270 * t649 * t8 * t327 / 0.2e1 - t2042 / 0.6e1;
1572   t2049 = t16 * t856;
1573   t2104 = -0.96962400000000000000e0 * t309 * t310 * t152 * t786 - 0.19392480000000000000e1 * t309 * t2049 * t857 * t116 - t270 * t271 * t863 * t85 / 0.6e1 - 0.32320800000000000000e0 * t200 * t201 * t110 * t863 - t293 * t294 * t869 * t116 / 0.6e1 - 0.32320800000000000000e0 * t309 * t310 * t885 * t116 - 0.5e1 / 0.27e2 * t1505 * t271 * t845 * t85 + 0.71824000000000000000e-1 * t640 * t201 * t110 * t845 - t270 * t271 * t128 * t781 / 0.2e1 - 0.5e1 / 0.27e2 * t1496 * t294 * t850 * t116 - t598 * t626 * t332 * t8 / 0.3e1 - t293 * t294 * t133 * t786 / 0.2e1 + t293 * t626 * t8 * t342 / 0.2e1 + 0.64641600000000000000e0 * t1694 - 0.10773600000000000000e0 * t1701 - 0.32320800000000000000e0 * t1716 + t1732 / 0.3e1 + 0.64641600000000000000e0 * t1740 - 0.32320800000000000000e0 * t1747;
1574   t2121 = t1750 / 0.9e1 - t1757 / 0.6e1 - t1770 / 0.6e1 + t1785 / 0.9e1 - t1788 / 0.3e1 + 0.43094400000000000000e0 * t1791 + t1436 + t1466 + t1477 + t1520 - t1822 / 0.2e1 + t1833 / 0.9e1 + t1836 / 0.3e1 + t1839 / 0.9e1 + t1849 / 0.3e1 - 0.32320800000000000000e0 * t1852 + 0.16160400000000000000e0 * t1855 - t1555 / 0.3e1 - 0.5e1 / 0.9e1 * t1558;
1575   t2166 = t1562 / 0.3e1 - 0.5e1 / 0.9e1 * t1565 - 0.26934000000000000000e0 * t1568 - t1863 / 0.2e1 - 0.5e1 / 0.27e2 * t1293 - 0.5e1 / 0.27e2 * t1336 - 0.89780000000000000000e-1 * t1354 + 0.64641600000000000000e0 * t293 * t752 * t824 * t116 + 0.19392480000000000000e1 * t309 * t752 * t820 * t116 + 0.21547200000000000000e0 * t598 * t310 * t836 * t116 + 0.64641600000000000000e0 * t293 * t743 * t338 * t8 - 0.32320800000000000000e0 * t293 * t310 * t840 * t116 - 0.32320800000000000000e0 * t293 * t310 * t815 * t116 + t640 * t271 * t828 * t85 / 0.3e1 - 0.32320800000000000000e0 * t270 * t201 * t279 * t327 + t598 * t294 * t832 * t116 / 0.3e1 - 0.21547200000000000000e0 * t1881 + t1392 + t1396;
1576   t2173 = t9 * (t2044 + t2104 + t2121 + t2166) - t818 + t823 + t827 - t831 - t835 - t839 - t843 - t714 - 0.12928320000000000000e1 * t717 + t721 / 0.3e1 - 0.32320800000000000000e0 * t725;
1577   t2182 = -t729 / 0.2e1 - 0.96962400000000000000e0 * t733 + t737 / 0.3e1 + t741 + 0.19392480000000000000e1 * t746 - t750 / 0.2e1 + 0.19392480000000000000e1 * t755 - 0.96962400000000000000e0 * t759 - t1647 - t1649 - 0.64641600e0 * t616 - t1657;
1578   t2188 = t849 + t854 - t861 - t867 - t873 - t889 - t762 / 0.3e1 - t766 - t769 / 0.3e1 - t772 + 0.96962400000000000000e0 * t775 - 0.48481200000000000000e0 * t778;
1579   t2193 = -0.3e1 / 0.2e1 * t784 - 0.3e1 / 0.2e1 * t789 + t668 + t672 - t676 + t679 + 0.48481200000000000000e0 * t682 + t508 + t535 + t550 - 0.64641600000000000000e0 * t801 - t1425 - t1426;
1580   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1581     v4rho4[3] = t2173 + t2182 + t2188 + t2193;
1582 
1583   t2203 = t342 * t342;
1584   t2209 = 0.1e1 / t855 / t71;
1585   t2210 = t350 * t350;
1586   t2215 = t366 * t366;
1587   t2220 = t321 * t321;
1588   t2225 = t327 * t327;
1589   t2230 = t332 * t332;
1590   t2236 = 0.12e2 * t540 + 0.12e2 * t1261;
1591   t2263 = -0.12928320000000000000e1 * t293 * t310 * t338 * t342 + 0.77569920000000000000e1 * t309 * t752 * t338 * t366 - t7 * t248 * t74 * t2203 / 0.2e1 + 0.58177440000000000000e1 * t7 * t54 * t2209 * t2210 + 0.14544360000000000000e1 * t7 * t54 * t349 * t2215 - 0.5e1 / 0.27e2 * t7 * t1271 * t47 * t2220 - t7 * t189 * t47 * t2225 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1316 * t74 * t2230 - t7 * t80 * t47 * t2236 / 0.2e1 + t7 * t115 * t74 * t2236 / 0.2e1 - 0.24240600000000000000e0 * t7 * t54 * t139 * (0.26180e5 / 0.19683e5 * t23 * sigma[2] / t57 / t146 + 0.67936424205399121722e-1 * t35 * t63 / t56 / t879 / rho[1]) - 0.19392480000000000000e1 * t1977;
1592   t2301 = -0.2e1 / 0.3e1 * t1980 - 0.2e1 / 0.3e1 * t1983 - 0.32320800000000000000e0 * t1986 + 0.4e1 / 0.27e2 * t1989 + 0.4e1 / 0.27e2 * t1992 - 0.2e1 / 0.3e1 * t293 * t294 * t133 * t869 - 0.64641600000000000000e0 * t293 * t310 * t366 * t332 - 0.12928320000000000000e1 * t309 * t310 * t885 * t133 - 0.19392480000000000000e1 * t309 * t310 * t366 * t342 - 0.12928320000000000000e1 * t309 * t310 * t152 * t869 - 0.87266160000000000000e1 * t819 * t2049 * t350 * t366 + 0.19392480000000000000e1 * t819 * t752 * t152 * t885 + 0.2e1 / 0.3e1 * t640 * t271 * t321 * t327;
1593   t2333 = -0.2e1 / 0.3e1 * t270 * t271 * t128 * t863 + 0.38784960000000000000e1 * t309 * t752 * t350 * t342 + 0.12928320000000000000e1 * t293 * t752 * t350 * t332 - 0.77569920000000000000e1 * t309 * t2049 * t857 * t133 + 0.2e1 / 0.3e1 * t598 * t294 * t332 * t342 + 0.28729600000000000000e0 * t598 * t310 * t152 * t850 + 0.25856640000000000000e1 * t2016 + 0.19392480000000000000e1 * t2019 - 0.43094400000000000000e0 * t2022 - 0.12928320000000000000e1 * t2025 - 0.12928320000000000000e1 * t2028 - 0.2e1 / 0.3e1 * t2035;
1594   t2345 = -0.2e1 / 0.3e1 * t2042 + 0.86188800000000000000e0 * t1791 + 0.2e1 / 0.9e1 * t1833 + 0.2e1 / 0.3e1 * t1836 + 0.2e1 / 0.9e1 * t1839 + 0.2e1 / 0.3e1 * t1849 - 0.64641600000000000000e0 * t1852 + 0.32320800000000000000e0 * t1855 - 0.20e2 / 0.27e2 * t1558 - 0.20e2 / 0.27e2 * t1565 - 0.35912000000000000000e0 * t1568 + t1392 + t1396;
1595   t2361 = t9 * (t2263 + t2301 + t2333 + t2345) - 0.38784960000000000000e1 * t817 + 0.58177440000000000000e1 * t822 + 0.77569920000000000000e1 * t826 - 0.2e1 * t830 - 0.2e1 * t834 - 0.12928320000000000000e1 * t838 - 0.38784960000000000000e1 * t842 - 0.25856640000000000000e1 * t717 + 0.4e1 / 0.9e1 * t848 + 0.4e1 / 0.9e1 * t853 - 0.58177440000000000000e1 * t860 - 0.2e1 * t866;
1596   t2373 = -0.2e1 * t872 - 0.96962400000000000000e0 * t888 - 0.2e1 / 0.3e1 * t762 - 0.2e1 * t766 - 0.2e1 / 0.3e1 * t769 - 0.2e1 * t772 + 0.19392480000000000000e1 * t775 - 0.96962400000000000000e0 * t778 + 0.4e1 / 0.3e1 * t672 + 0.4e1 / 0.3e1 * t679 + 0.64641600000000000000e0 * t682 - t1425 - t1426;
1597   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1598     v4rho4[4] = t2361 + t2373;
1599 
1600   t2374 = t14 * t514;
1601   t2379 = t188 * t224;
1602   t2405 = t200 * t480 * t920;
1603   t2409 = 0.71824000000000000000e-1 * t7 * t216 * t169;
1604   t2421 = t7 * t92 * t394;
1605   t2424 = t7 * t17 * t951;
1606   t2427 = t200 * t467 * t379;
1607   t2430 = t479 * t1187 * t933;
1608   t2432 = t494 * t97;
1609   t2437 = -0.58177440000000000000e1 * t377 * t2374 * t379 * t225 + 0.64641600000000000000e0 * t377 * t2379 * t915 * t110 - 0.32320800000000000000e0 * t377 * t914 * t379 * t209 + 0.38784960000000000000e1 * t377 * t919 * t924 * t110 + 0.19392480000000000000e1 * t377 * t919 * t928 * t110 + 0.19392480000000000000e1 * t377 * t919 * t379 * t241 - 0.43633080000000000000e1 * t377 * t932 * t384 * t241 + 0.51713280000000000000e1 * t2405 + t2409 - 0.24240600000000000000e0 * t161 * t162 * t97 * (-0.1540e4 / 0.6561e4 * t23 * t555 - 0.16304741809295789214e-1 * t35 * sigma[0] * t561) - 0.32320800000000000000e0 * t2421 - 0.96962400000000000000e0 * t2424 - 0.43094400000000000000e0 * t2427 - 0.58177440000000000000e1 * t2430 + 0.71824000000000000000e-1 * t377 * t2432 * t168 * t496;
1610   t2471 = t200 * t201 * t924;
1611   t2474 = t200 * t201 * t928;
1612   t2477 = t479 * t1183 * t384;
1613   t2480 = t479 * t480 * t937;
1614   t2483 = t479 * t480 * t941;
1615   t2485 = t15 * t1302;
1616   t2491 = t270 * t201 * t915;
1617   t2493 = -0.32320800000000000000e0 * t377 * t914 * t393 * t190 - 0.96962400000000000000e0 * t377 * t378 * t950 * t85 - 0.96962400000000000000e0 * t377 * t378 * t393 * t209 - 0.32320800000000000000e0 * t377 * t378 * t168 * t543 - 0.43633080000000000000e1 * t377 * t932 * t393 * t225 + 0.14544360000000000000e1 * t377 * t383 * t950 * t110 + 0.14544360000000000000e1 * t377 * t383 * t393 * t241 + 0.48481200000000000000e0 * t377 * t383 * t168 * t565 - 0.25856640000000000000e1 * t2471 - 0.12928320000000000000e1 * t2474 + 0.64641600000000000000e0 * t2477 + 0.38784960000000000000e1 * t2480 + 0.19392480000000000000e1 * t2483 + 0.58177440000000000000e1 * t377 * t2485 * t168 * t515 - 0.43094400000000000000e0 * t2491;
1618   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1619     v4rho3sigma[0] = t2437 + t2493;
1620 
1621   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1622     v4rho3sigma[1] = 0.0e0;
1623 
1624   t2495 = t293 * t310 * t962;
1625   t2498 = t309 * t611 * t402;
1626   t2501 = t309 * t310 * t966;
1627   t2505 = 0.71824000000000000000e-1 * t7 * t263 * t180;
1628   t2506 = t521 * t139;
1629   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1630     v4rho3sigma[2] = -0.43094400000000000000e0 * t2495 - 0.43094400000000000000e0 * t2498 - 0.12928320000000000000e1 * t2501 + t2505 + 0.71824000000000000000e-1 * t377 * t2506 * t179 * t523 - 0.32320800000000000000e0 * t377 * t961 * t402 * t258 - 0.32320800000000000000e0 * t377 * t401 * t179 * t573;
1631 
1632   t2536 = 0.28729600000000000000e0 * t974 * t978;
1633   t2538 = 0.17237760000000000000e1 * t974 * t982;
1634   t2558 = 0.71824000000000000000e-1 * t377 * t2432 * t406 * t190 - 0.21547200000000000000e0 * t990 * t914 * t991 * t85 - 0.19392480000000000000e1 * t377 * t2374 * t406 * t225 + 0.12928320000000000000e1 * t990 * t919 * t991 * t110 - t2536 + t2538 - 0.21547200000000000000e0 * t377 * t914 * t985 * t85 - 0.10773600000000000000e0 * t377 * t914 * t406 * t209 + 0.12928320000000000000e1 * t377 * t919 * t985 * t110 + 0.64641600000000000000e0 * t377 * t919 * t406 * t241 + 0.17237760000000000000e1 * t2405 + t2409 - 0.21547200000000000000e0 * t2421 - 0.32320800000000000000e0 * t2424;
1635   t2559 = t989 * t3;
1636   t2564 = t284 * t1 * t376;
1637   t2565 = t2564 * t992;
1638   t2579 = 0.86188800000000000000e0 * t974 * t986;
1639   t2584 = t91 * t1;
1640   t2585 = t2584 * t376;
1641   t2586 = t2585 * t407;
1642   t2590 = t160 * t376 * t188;
1643   t2594 = -0.32320800000000000000e0 * t2559 * t199 * t169 + 0.10773600000000000000e0 * t2565 - 0.64641600000000000000e0 * t990 * t378 * t393 * t8 - 0.28729600000000000000e0 * t2427 - 0.19392480000000000000e1 * t2430 - 0.86188800000000000000e0 * t2471 - 0.43094400000000000000e0 * t2474 + 0.43094400000000000000e0 * t2477 + 0.12928320000000000000e1 * t2480 + 0.64641600000000000000e0 * t2483 - t2579 - 0.32320800000000000000e0 * t377 * t378 * t950 * t128 - 0.14364800000000000000e0 * t2586 - 0.14364800000000000000e0 * t2491 + 0.43094400000000000000e0 * t2590 * t1162 * t693;
1644   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645     v4rho3sigma[3] = t2558 + t2594;
1646 
1647   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1648     v4rho3sigma[4] = 0.0e0;
1649 
1650   t2598 = t2584 * t3;
1651   t2599 = t2598 * t426;
1652   t2602 = 0.86188800000000000000e0 * t974 * t1017;
1653   t2607 = t2564 * t1003;
1654   t2609 = t2585 * t416;
1655   t2612 = 0.17237760000000000000e1 * t974 * t1010;
1656   t2621 = t2585 * t411;
1657   t2624 = 0.28729600000000000000e0 * t974 * t999;
1658   t2636 = t247 * t349;
1659   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1660     v4rho3sigma[5] = 0.32320800000000000000e0 * t2559 * t308 * t180 - 0.10773600000000000000e0 * t2599 - t2602 - 0.32320800000000000000e0 * t377 * t401 * t424 * t258 - 0.10773600000000000000e0 * t2607 + 0.21547200000000000000e0 * t2609 + t2612 + 0.64641600000000000000e0 * t377 * t1008 * t415 * t258 - 0.10773600000000000000e0 * t377 * t961 * t424 * t249 - 0.14364800000000000000e0 * t2621 - t2624 - 0.10773600000000000000e0 * t377 * t961 * t410 * t258 + 0.21547200000000000000e0 * t990 * t961 * t1002 * t116 - 0.14364800000000000000e0 * t2495 - 0.28729600000000000000e0 * t2498 - 0.43094400000000000000e0 * t2501 + t2505 + 0.21547200000000000000e0 * t377 * t2636 * t415 * t249 + 0.71824000000000000000e-1 * t377 * t2506 * t410 * t249;
1661 
1662   t2650 = t974 * t1022;
1663   t2668 = t974 * t1026;
1664   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1665     v4rho3sigma[6] = -0.28729600000000000000e0 * t2586 - t2536 + t2538 - t2579 - 0.86188800000000000000e0 * t2565 - 0.14364800000000000000e0 * t2427 + t2409 + 0.21547200000000000000e0 * t2477 - 0.10773600000000000000e0 * t2421 - 0.14364800000000000000e0 * t2650 + 0.71824000000000000000e-1 * t377 * t2432 * t1021 * t85 + 0.21547200000000000000e0 * t377 * t2379 * t1021 * t110 - 0.10773600000000000000e0 * t377 * t914 * t393 * t321 - 0.21547200000000000000e0 * t990 * t914 * t406 * t8 - 0.43094400000000000000e0 * t2668 - 0.10773600000000000000e0 * t377 * t914 * t1025 * t85 + 0.64641600000000000000e0 * t377 * t919 * t1025 * t110 - 0.32320800000000000000e0 * t377 * t378 * t393 * t327 - 0.32320800000000000000e0 * t377 * t378 * t168 * t781;
1666 
1667   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1668     v4rho3sigma[7] = 0.0e0;
1669 
1670   t2698 = t52 * t856;
1671   t2707 = t974 * t1037;
1672   t2722 = t1013 * t1068;
1673   t2724 = -0.12928320000000000000e1 * t990 * t1008 * t1002 * t152 - 0.21547200000000000000e0 * t377 * t961 * t1040 * t116 - 0.10773600000000000000e0 * t377 * t961 * t1044 * t116 - 0.19392480000000000000e1 * t377 * t2698 * t1049 * t116 + 0.71824000000000000000e-1 * t377 * t2506 * t1032 * t116 + 0.17237760000000000000e1 * t2707 + 0.12928320000000000000e1 * t377 * t1008 * t1053 * t116 + 0.64641600000000000000e0 * t377 * t1008 * t1057 * t116 + 0.21547200000000000000e0 * t990 * t961 * t410 * t8 + t2612 - t2624 + t2505 - 0.21547200000000000000e0 * t2599 - 0.32320800000000000000e0 * t2722;
1674   t2725 = t974 * t1054;
1675   t2727 = t974 * t1058;
1676   t2729 = t974 * t1041;
1677   t2731 = t974 * t1045;
1678   t2733 = t974 * t1050;
1679   t2735 = t974 * t1033;
1680   t2754 = t160 * t376 * t247;
1681   t2758 = 0.12928320000000000000e1 * t2725 + 0.64641600000000000000e0 * t2727 - 0.86188800000000000000e0 * t2729 - 0.43094400000000000000e0 * t2731 - 0.19392480000000000000e1 * t2733 - 0.14364800000000000000e0 * t2735 - 0.32320800000000000000e0 * t377 * t401 * t1066 * t116 + 0.64641600000000000000e0 * t990 * t401 * t424 * t8 - 0.32320800000000000000e0 * t377 * t401 * t179 * t786 - t2602 + 0.86188800000000000000e0 * t2607 + 0.43094400000000000000e0 * t2609 - 0.28729600000000000000e0 * t2621 - 0.14364800000000000000e0 * t2498 + 0.43094400000000000000e0 * t2754 * t1171 * t799;
1682   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1683     v4rho3sigma[8] = t2724 + t2758;
1684 
1685   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1686     v4rho3sigma[9] = -0.43094400000000000000e0 * t2586 - 0.43094400000000000000e0 * t2650 - 0.12928320000000000000e1 * t2668 + t2409 + 0.71824000000000000000e-1 * t377 * t2432 * t168 * t845 - 0.32320800000000000000e0 * t377 * t914 * t406 * t327 - 0.32320800000000000000e0 * t377 * t378 * t168 * t863;
1687 
1688   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1689     v4rho3sigma[10] = 0.0e0;
1690 
1691   t2774 = t53 * t2209;
1692   t2835 = 0.58177440000000000000e1 * t377 * t2774 * t179 * t857 + 0.71824000000000000000e-1 * t377 * t2506 * t179 * t850 - 0.32320800000000000000e0 * t377 * t401 * t179 * t869 - 0.43633080000000000000e1 * t377 * t1048 * t424 * t350 + 0.14544360000000000000e1 * t377 * t414 * t1066 * t152 + 0.14544360000000000000e1 * t377 * t414 * t424 * t366 + 0.48481200000000000000e0 * t377 * t414 * t179 * t885 - 0.32320800000000000000e0 * t377 * t961 * t424 * t332 - 0.96962400000000000000e0 * t377 * t401 * t1066 * t133 - 0.96962400000000000000e0 * t377 * t401 * t424 * t342 - 0.43633080000000000000e1 * t377 * t1048 * t415 * t366 + 0.64641600000000000000e0 * t377 * t2636 * t1032 * t152 - 0.32320800000000000000e0 * t377 * t961 * t410 * t342 + 0.38784960000000000000e1 * t377 * t1008 * t1040 * t152 + 0.19392480000000000000e1 * t377 * t1008 * t1044 * t152;
1693   t2865 = 0.19392480000000000000e1 * t377 * t1008 * t410 * t366 - 0.58177440000000000000e1 * t377 * t2698 * t410 * t350 + 0.51713280000000000000e1 * t2707 + t2505 - 0.32320800000000000000e0 * t2599 - 0.96962400000000000000e0 * t2722 - 0.24240600000000000000e0 * t161 * t173 * t139 * (-0.1540e4 / 0.6561e4 * t23 * t875 - 0.16304741809295789214e-1 * t35 * sigma[2] * t881) + 0.38784960000000000000e1 * t2725 + 0.19392480000000000000e1 * t2727 - 0.25856640000000000000e1 * t2729 - 0.12928320000000000000e1 * t2731 - 0.58177440000000000000e1 * t2733 - 0.43094400000000000000e0 * t2735 + 0.64641600000000000000e0 * t2609 - 0.43094400000000000000e0 * t2621;
1694   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1695     v4rho3sigma[11] = t2835 + t2865;
1696 
1697   t2866 = t973 * t1087;
1698   t2867 = t2866 * t1092;
1699   t2869 = t973 * t1095;
1700   t2870 = t2869 * t1099;
1701   t2884 = t34 * t106;
1702   t2914 = -0.56732544158241460480e-4 * t2867 + 0.85098816237362190720e-4 * t2870 - 0.21274704059340547680e-4 * t1088 * t1089 * t1090 * t209 + 0.31912056089010821520e-4 * t1096 * t1097 * t1090 * t241 - 0.70915680197801825600e-5 * t2590 * t1089 * t1090 * t190 + 0.22693017663296584192e-3 * t1088 * t1089 * t2884 * t85 - 0.95736168267032464560e-4 * t1096 * t514 * t32 * t1090 * t225 - 0.34039526494944876288e-3 * t1096 * t1097 * t2884 * t110 - 0.14544360000000000000e1 * t377 * t932 * t429 * t241 + 0.96962400000000000000e0 * t377 * t383 * t168 * t950 + 0.21547200000000000000e0 * t377 * t2379 * t429 * t190 + 0.58177440000000000000e1 * t377 * t2485 * t429 * t225;
1703   t2915 = t974 * t1074;
1704   t2917 = t974 * t1078;
1705   t2919 = t974 * t1082;
1706   t2926 = 0.70915680197801825600e-5 * t2585 * t436;
1707   t2927 = t974 * t1103;
1708   t2946 = 0.21547200000000000000e0 * t2598 * t431;
1709   t2947 = t393 * t393;
1710   t2952 = t1097 * t34;
1711   t2958 = 0.17237760000000000000e1 * t2915 - 0.38784960000000000000e1 * t2917 + 0.25856640000000000000e1 * t2919 + 0.64641600000000000000e0 * t377 * t919 * t429 * t209 - t2926 + 0.22693017663296584192e-3 * t2927 - 0.53895916950329387456e-3 * t377 * t434 * t35 * t237 - 0.58177440000000000000e1 * t377 * t932 * t384 * t393 - 0.38784960000000000000e1 * t377 * t2374 * t1073 * t110 + 0.25856640000000000000e1 * t377 * t919 * t379 * t393 + t2946 + 0.96962400000000000000e0 * t161 * t162 * t224 * t2947 + 0.85098816237362190720e-4 * t1088 * t2952 * t40 * t85 * t110;
1712   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1713     v4rho2sigma2[0] = t2914 + t2958;
1714 
1715   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1716     v4rho2sigma2[1] = 0.0e0;
1717 
1718   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1719     v4rho2sigma2[2] = 0.0e0;
1720 
1721   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1722     v4rho2sigma2[3] = 0.0e0;
1723 
1724   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1725     v4rho2sigma2[4] = 0.0e0;
1726 
1727   t2960 = 0.21547200000000000000e0 * t2598 * t441;
1728   t2961 = t974 * t1109;
1729   t2972 = 0.70915680197801825600e-5 * t2585 * t446;
1730   t2973 = t973 * t1114;
1731   t2974 = t2973 * t1119;
1732   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1733     v4rho2sigma2[5] = t2960 + 0.17237760000000000000e1 * t2961 + 0.21547200000000000000e0 * t377 * t2636 * t439 * t249 + 0.64641600000000000000e0 * t377 * t1008 * t439 * t258 - t2972 - 0.56732544158241460480e-4 * t2974 - 0.70915680197801825600e-5 * t2754 * t1116 * t1117 * t249 - 0.21274704059340547680e-4 * t1115 * t1116 * t1117 * t258;
1734 
1735   t2987 = t974 * t1123;
1736   t3008 = t2866 * t1127;
1737   t3011 = t40 * t128;
1738   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1739     v4rho2sigma2[6] = t2946 + 0.86188800000000000000e0 * t2915 - 0.19392480000000000000e1 * t2917 + 0.12928320000000000000e1 * t2919 + 0.86188800000000000000e0 * t2987 + 0.21547200000000000000e0 * t377 * t2379 * t1122 * t85 - 0.19392480000000000000e1 * t377 * t2374 * t1122 * t110 + 0.12928320000000000000e1 * t377 * t919 * t406 * t393 + 0.64641600000000000000e0 * t990 * t919 * t429 * t8 - t2926 - 0.28366272079120730240e-4 * t2867 + 0.42549408118681095360e-4 * t2870 + 0.11346508831648292096e-3 * t2927 - 0.28366272079120730240e-4 * t3008 - 0.70915680197801825600e-5 * t2590 * t1089 * t34 * t3011 * t85 + 0.42549408118681095360e-4 * t1088 * t2952 * t3011 * t110 + 0.11346508831648292096e-3 * t1088 * t1089 * t2884 * t128 - 0.21274704059340547680e-4 * t989 * t1087 * t1089 * t1090 * t8;
1740 
1741   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1742     v4rho2sigma2[7] = 0.0e0;
1743 
1744   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1745     v4rho2sigma2[8] = 0.0e0;
1746 
1747   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1748     v4rho2sigma2[9] = 0.0e0;
1749 
1750   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1751     v4rho2sigma2[10] = 0.0e0;
1752 
1753   t3030 = t974 * t1131;
1754   t3040 = t974 * t1135;
1755   t3046 = t974 * t1139;
1756   t3053 = t2973 * t1143;
1757   t3056 = t67 * t133;
1758   t3066 = t973 * t1146;
1759   t3067 = t3066 * t1150;
1760   t3069 = t1148 * t34;
1761   t3075 = t974 * t1154;
1762   t3077 = t34 * t148;
1763   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1764     v4rho2sigma2[11] = t2960 + 0.86188800000000000000e0 * t2961 + 0.86188800000000000000e0 * t3030 + 0.21547200000000000000e0 * t377 * t2636 * t1130 * t116 - 0.64641600000000000000e0 * t990 * t1008 * t439 * t8 - 0.19392480000000000000e1 * t3040 - 0.19392480000000000000e1 * t377 * t2698 * t1134 * t116 + 0.12928320000000000000e1 * t3046 + 0.12928320000000000000e1 * t377 * t1008 * t1138 * t116 - t2972 - 0.28366272079120730240e-4 * t2974 - 0.28366272079120730240e-4 * t3053 - 0.70915680197801825600e-5 * t2754 * t1116 * t34 * t3056 * t116 + 0.21274704059340547680e-4 * t989 * t1114 * t1116 * t1117 * t8 + 0.42549408118681095360e-4 * t3067 + 0.42549408118681095360e-4 * t1115 * t3069 * t67 * t152 * t116 + 0.11346508831648292096e-3 * t3075 + 0.11346508831648292096e-3 * t1115 * t1116 * t3077 * t116;
1765 
1766   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1767     v4rho2sigma2[12] = t2946 + 0.17237760000000000000e1 * t2987 + 0.21547200000000000000e0 * t377 * t2379 * t429 * t321 + 0.64641600000000000000e0 * t377 * t919 * t429 * t327 - t2926 - 0.56732544158241460480e-4 * t3008 - 0.70915680197801825600e-5 * t2590 * t1089 * t1090 * t321 - 0.21274704059340547680e-4 * t1088 * t1089 * t1090 * t327;
1768 
1769   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1770     v4rho2sigma2[13] = 0.0e0;
1771 
1772   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1773     v4rho2sigma2[14] = 0.0e0;
1774 
1775   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1776     v4rho2sigma2[15] = 0.0e0;
1777 
1778   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1779     v4rho2sigma2[16] = 0.0e0;
1780 
1781   t3140 = -0.70915680197801825600e-5 * t2754 * t1116 * t1117 * t332 + 0.22693017663296584192e-3 * t1115 * t1116 * t3077 * t133 - 0.95736168267032464560e-4 * t1147 * t856 * t32 * t1117 * t350 - 0.34039526494944876288e-3 * t1147 * t1148 * t3077 * t152 - 0.56732544158241460480e-4 * t3053 + 0.85098816237362190720e-4 * t3067 - 0.21274704059340547680e-4 * t1115 * t1116 * t1117 * t342 + 0.31912056089010821520e-4 * t1147 * t1148 * t1117 * t366 + 0.96962400000000000000e0 * t377 * t414 * t179 * t1066 + 0.21547200000000000000e0 * t377 * t2636 * t439 * t332 + 0.58177440000000000000e1 * t377 * t2774 * t439 * t350 + 0.17237760000000000000e1 * t3030;
1782   t3168 = t424 * t424;
1783   t3177 = -0.38784960000000000000e1 * t3040 + 0.25856640000000000000e1 * t3046 + 0.64641600000000000000e0 * t377 * t1008 * t439 * t342 - 0.14544360000000000000e1 * t377 * t1048 * t439 * t366 - 0.58177440000000000000e1 * t377 * t1048 * t415 * t424 + 0.22693017663296584192e-3 * t3075 - 0.53895916950329387456e-3 * t377 * t444 * t35 * t362 - 0.38784960000000000000e1 * t377 * t2698 * t1130 * t152 + 0.25856640000000000000e1 * t377 * t1008 * t410 * t424 - t2972 + t2960 + 0.96962400000000000000e0 * t161 * t173 * t349 * t3168 + 0.85098816237362190720e-4 * t1115 * t3069 * t3056 * t152;
1784   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1785     v4rho2sigma2[17] = t3140 + t3177;
1786 
1787   t3179 = 0.19392480000000000000e1 * t1013 * t1159;
1788   t3193 = 0.12764822435604328608e-3 * t2869 * t1163;
1789   t3194 = t1162 * t32;
1790   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1791     v4rhosigma3[0] = -t3179 - 0.19392480000000000000e1 * t377 * t2374 * t1157 * t85 + 0.58177440000000000000e1 * t377 * t2485 * t1157 * t110 - 0.43633080000000000000e1 * t377 * t932 * t429 * t393 + t3193 + 0.12764822435604328608e-3 * t1088 * t3194 * t1091 - 0.28720850480109739369e-3 * t1096 * t514 * t168 * t32 * t1098 + 0.95736168267032464562e-4 * t1096 * t224 * t393 * t435 - 0.51059289742417314433e-3 * t1096 * t1162 * t1102;
1792 
1793   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1794     v4rhosigma3[1] = 0.0e0;
1795 
1796   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1797     v4rhosigma3[2] = 0.0e0;
1798 
1799   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1800     v4rhosigma3[3] = 0.0e0;
1801 
1802   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1803     v4rhosigma3[4] = 0.0e0;
1804 
1805   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1806     v4rhosigma3[5] = 0.0e0;
1807 
1808   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1809     v4rhosigma3[6] = 0.0e0;
1810 
1811   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1812     v4rhosigma3[7] = 0.0e0;
1813 
1814   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1815     v4rhosigma3[8] = 0.0e0;
1816 
1817   t3211 = 0.19392480000000000000e1 * t1013 * t1168;
1818   t3217 = 0.12764822435604328608e-3 * t3066 * t1172;
1819   t3218 = t1171 * t32;
1820   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1821     v4rhosigma3[9] = -t3211 - 0.19392480000000000000e1 * t377 * t2698 * t1166 * t116 + t3217 + 0.12764822435604328608e-3 * t1115 * t3218 * t1118;
1822 
1823   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1824     v4rhosigma3[10] = -t3179 - 0.19392480000000000000e1 * t377 * t2374 * t1157 * t128 + t3193 + 0.12764822435604328608e-3 * t1088 * t3194 * t1126;
1825 
1826   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1827     v4rhosigma3[11] = 0.0e0;
1828 
1829   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1830     v4rhosigma3[12] = 0.0e0;
1831 
1832   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1833     v4rhosigma3[13] = 0.0e0;
1834 
1835   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1836     v4rhosigma3[14] = 0.0e0;
1837 
1838   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1839     v4rhosigma3[15] = 0.0e0;
1840 
1841   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1842     v4rhosigma3[16] = 0.0e0;
1843 
1844   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1845     v4rhosigma3[17] = 0.0e0;
1846 
1847   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1848     v4rhosigma3[18] = 0.0e0;
1849 
1850   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1851     v4rhosigma3[19] = -t3211 - 0.19392480000000000000e1 * t377 * t2698 * t1166 * t133 + 0.58177440000000000000e1 * t377 * t2774 * t1166 * t152 - 0.43633080000000000000e1 * t377 * t1048 * t439 * t424 + t3217 + 0.12764822435604328608e-3 * t1115 * t3218 * t1142 - 0.28720850480109739369e-3 * t1147 * t856 * t179 * t32 * t1149 + 0.95736168267032464562e-4 * t1147 * t349 * t424 * t445 - 0.51059289742417314433e-3 * t1147 * t1171 * t1153;
1852 
1853   t3256 = t429 * t429;
1854   t3265 = t19 * t19;
1855   t3268 = t18 / t21 / t3265;
1856   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1857     v4sigma4[0] = 0.58177440000000000000e1 * t161 * t162 * t1302 * t3256 - 0.57441700960219478737e-3 * t1096 * t514 * t429 * t435 + 0.37810177612986286844e-7 * t377 * t383 * t3268 / t26 / t559 / t24;
1858 
1859   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1860     v4sigma4[1] = 0.0e0;
1861 
1862   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1863     v4sigma4[2] = 0.0e0;
1864 
1865   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1866     v4sigma4[3] = 0.0e0;
1867 
1868   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1869     v4sigma4[4] = 0.0e0;
1870 
1871   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1872     v4sigma4[5] = 0.0e0;
1873 
1874   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1875     v4sigma4[6] = 0.0e0;
1876 
1877   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1878     v4sigma4[7] = 0.0e0;
1879 
1880   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1881     v4sigma4[8] = 0.0e0;
1882 
1883   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1884     v4sigma4[9] = 0.0e0;
1885 
1886   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1887     v4sigma4[10] = 0.0e0;
1888 
1889   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1890     v4sigma4[11] = 0.0e0;
1891 
1892   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1893     v4sigma4[12] = 0.0e0;
1894 
1895   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1896     v4sigma4[13] = 0.0e0;
1897 
1898   t3276 = t439 * t439;
1899   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1900     v4sigma4[14] = 0.58177440000000000000e1 * t161 * t173 * t2209 * t3276 - 0.57441700960219478737e-3 * t1147 * t856 * t439 * t445 + 0.37810177612986286844e-7 * t377 * t414 * t3268 / t57 / t879 / t55;
1901 
1902 #ifndef XC_DONT_COMPILE_MXC
1903 
1904   if(order < 5) return;
1905 
1906 
1907 #endif
1908 
1909 #endif
1910 
1911 #endif
1912 
1913 #endif
1914 
1915 #endif
1916 
1917 
1918 }
1919 
1920