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