1 /*
2 This file was generated automatically with ./scripts/maple2c.pl.
3 Do not edit this file directly as it can be overwritten!!
4
5 This Source Code Form is subject to the terms of the Mozilla Public
6 License, v. 2.0. If a copy of the MPL was not distributed with this
7 file, You can obtain one at http://mozilla.org/MPL/2.0/.
8
9 Maple version : Maple 2016 (X86 64 LINUX)
10 Maple source : ./maple/gga_exc/gga_c_p86.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t2, t3, t4, t5, t6, t7, t9;
24 double t10, t11, t12, t13, t16, t19, t26, t27;
25 double t29, t30, t32, t33, t34, t35, t38, t40;
26 double t43, t45, t46, t48, t49, t50, t51, t52;
27 double t54, t57, t58, t59;
28
29 #ifndef XC_DONT_COMPILE_VXC
30 double t60, t61, t63, t64, t66, t70, t71, t73;
31 double t82, t83, t85, t86, t87, t89, t90, t91;
32 double t96, t98, t100, t101, t102, t105, t107, t109;
33 double t110, t114, t117, t118, t119, t120, t121, t122;
34 double t125, t126, t127, t129, t130, t132;
35
36 #ifndef XC_DONT_COMPILE_FXC
37 double t139, t140, t145, t146, t148, t155, t156, t158;
38 double t168, t169, t171, t172, t173, t175, t177, t180;
39 double t181, t182, t183, t186, t191, t193, t195, t199;
40 double t200, t201, t206, t208, t210, t215, t218, t219;
41 double t220, t221, t222, t223, t224, t225, t226, t227;
42 double t229, t230, t233, t234, t236, t237, t238, t240;
43 double t241, t242, t244, t248, t249, t251, t253, t254;
44 double t256;
45
46 #ifndef XC_DONT_COMPILE_KXC
47 double t266, t268, t269, t271, t272, t273, t274, t275;
48 double t277, t278, t279, t281, t283, t285, t287, t289;
49 double t290, t306, t307, t310, t312, t314, t316, t319;
50 double t324, t325, t326, t327, t330, t335, t337, t339;
51 double t344, t347, t348, t349, t350, t351, t352, t354;
52 double t355, t357, t358, t360, t362, t364, t365, t366;
53 double t367, t371, t377, t387, t397, t398, t405, t406;
54 double t408, t409, t411, t413, t414, t415, t416, t417;
55 double t419, t421, t422, t423, t425, t427, t429, t431;
56 double t433, t434, t436, t439, t440, t441, t443, t444;
57 double t446, t451, t452, t454, t455, t456, t457, t459;
58 double t460, t461, t462, t464, t465, t466, t468;
59
60 #ifndef XC_DONT_COMPILE_LXC
61 double t470, t474, t481, t504, t510, t524, t533, t534;
62 double t546, t549, t553, t569, t587, t609, t615, t626;
63 double t634, t670, t739, t789, t799, t822;
64 #endif
65
66 #endif
67
68 #endif
69
70 #endif
71
72 #endif
73
74
75
76 t1 = M_CBRT3;
77 t2 = 0.1e1 / M_PI;
78 t3 = POW_1_3(t2);
79 t4 = t1 * t3;
80 t5 = M_CBRT4;
81 t6 = t5 * t5;
82 t7 = POW_1_3(rho[0]);
83 t9 = t6 / t7;
84 t10 = t4 * t9;
85 t11 = t10 / 0.4e1;
86 t12 = 0.1e1 <= t11;
87 t13 = sqrt(t10);
88 t16 = 0.1e1 + 0.52645000000000000000e0 * t13 + 0.83350000000000000000e-1 * t10;
89 t19 = log(t11);
90 t26 = my_piecewise3(t12, -0.1423e0 / t16, 0.311e-1 * t19 - 0.48e-1 + 0.50000000000000000000e-3 * t4 * t9 * t19 - 0.29000000000000000000e-2 * t10);
91 t27 = rho[0] * rho[0];
92 t29 = 0.1e1 / t7 / t27;
93 t30 = sigma[0] * t29;
94 t32 = t1 * t1;
95 t33 = t3 * t3;
96 t34 = t32 * t33;
97 t35 = t7 * t7;
98 t38 = t34 * t5 / t35;
99 t40 = 0.2568e-2 + 0.58165000000000000000e-2 * t10 + 0.18472500000000000000e-5 * t38;
100 t43 = 0.1e1 / rho[0];
101 t45 = 0.1e1 + 0.21807500000000000000e1 * t10 + 0.11800000000000000000e0 * t38 + 0.17639938117590219540e-1 * t43;
102 t46 = 0.1e1 / t45;
103 t48 = 0.1667e-2 + t40 * t46;
104 t49 = 0.1e1 / t48;
105 t50 = sqrt(sigma[0]);
106 t51 = t49 * t50;
107 t52 = pow(rho[0], 0.1e1 / 0.6e1);
108 t54 = 0.1e1 / t52 / rho[0];
109 t57 = exp(-0.81290825e-3 * t51 * t54);
110 t58 = t57 * t48;
111 t59 = t30 * t58;
112 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
113 zk[0] = t26 + t59;
114
115 #ifndef XC_DONT_COMPILE_VXC
116
117 if(order < 1) return;
118
119
120 t60 = t16 * t16;
121 t61 = 0.1e1 / t60;
122 t63 = 0.1e1 / t13 * t1;
123 t64 = t3 * t6;
124 t66 = 0.1e1 / t7 / rho[0];
125 t70 = t6 * t66;
126 t71 = t4 * t70;
127 t73 = -0.87741666666666666667e-1 * t63 * t64 * t66 - 0.27783333333333333333e-1 * t71;
128 t82 = my_piecewise3(t12, 0.1423e0 * t61 * t73, -0.10366666666666666667e-1 * t43 - 0.16666666666666666667e-3 * t4 * t70 * t19 + 0.80000000000000000000e-3 * t71);
129 t83 = t27 * rho[0];
130 t85 = 0.1e1 / t7 / t83;
131 t86 = sigma[0] * t85;
132 t87 = t86 * t58;
133 t89 = t48 * t48;
134 t90 = 0.1e1 / t89;
135 t91 = t90 * t50;
136 t96 = t34 * t5 / t35 / rho[0];
137 t98 = -0.19388333333333333333e-2 * t71 - 0.12315000000000000000e-5 * t96;
138 t100 = t45 * t45;
139 t101 = 0.1e1 / t100;
140 t102 = t40 * t101;
141 t105 = 0.1e1 / t27;
142 t107 = -0.72691666666666666667e0 * t71 - 0.78666666666666666667e-1 * t96 - 0.17639938117590219540e-1 * t105;
143 t109 = -t102 * t107 + t98 * t46;
144 t110 = t54 * t109;
145 t114 = 0.1e1 / t52 / t27;
146 t117 = 0.81290825e-3 * t91 * t110 + 0.94839295833333333334e-3 * t51 * t114;
147 t118 = t117 * t57;
148 t119 = t118 * t48;
149 t120 = t30 * t119;
150 t121 = t57 * t109;
151 t122 = t30 * t121;
152 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
153 vrho[0] = t26 + t59 + rho[0] * (t82 - 0.7e1 / 0.3e1 * t87 + t120 + t122);
154
155 t125 = t29 * t57;
156 t126 = t125 * t48;
157 t127 = sqrt(rho[0]);
158 t129 = 0.1e1 / t127 / t83;
159 t130 = t50 * t129;
160 t132 = 0.40645412500000000000e-3 * t130 * t57;
161 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
162 vsigma[0] = rho[0] * (t126 - t132);
163
164 #ifndef XC_DONT_COMPILE_FXC
165
166 if(order < 2) return;
167
168
169 t139 = 0.1e1 / t60 / t16;
170 t140 = t73 * t73;
171 t145 = 0.1e1 / t13 / t10 * t32;
172 t146 = t33 * t5;
173 t148 = 0.1e1 / t35 / t27;
174 t155 = t6 * t29;
175 t156 = t4 * t155;
176 t158 = -0.58494444444444444445e-1 * t145 * t146 * t148 + 0.11698888888888888889e0 * t63 * t64 * t29 + 0.37044444444444444444e-1 * t156;
177 t168 = my_piecewise3(t12, -0.2846e0 * t139 * t140 + 0.1423e0 * t61 * t158, 0.10366666666666666667e-1 * t105 + 0.22222222222222222223e-3 * t4 * t155 * t19 - 0.10111111111111111111e-2 * t156);
178 t169 = t27 * t27;
179 t171 = 0.1e1 / t7 / t169;
180 t172 = sigma[0] * t171;
181 t173 = t172 * t58;
182 t175 = t86 * t119;
183 t177 = t86 * t121;
184 t180 = 0.1e1 / t89 / t48;
185 t181 = t180 * t50;
186 t182 = t109 * t109;
187 t183 = t54 * t182;
188 t186 = t114 * t109;
189 t191 = t34 * t5 * t148;
190 t193 = 0.25851111111111111111e-2 * t156 + 0.20525000000000000000e-5 * t191;
191 t195 = t98 * t101;
192 t199 = 0.1e1 / t100 / t45;
193 t200 = t40 * t199;
194 t201 = t107 * t107;
195 t206 = 0.1e1 / t83;
196 t208 = 0.96922222222222222222e0 * t156 + 0.13111111111111111111e0 * t191 + 0.35279876235180439080e-1 * t206;
197 t210 = -t102 * t208 - 0.2e1 * t195 * t107 + t193 * t46 + 0.2e1 * t200 * t201;
198 t215 = 0.1e1 / t52 / t83;
199 t218 = -0.162581650e-2 * t181 * t183 - 0.18967859166666666667e-2 * t91 * t186 + 0.81290825e-3 * t91 * t54 * t210 - 0.20548514097222222222e-2 * t51 * t215;
200 t219 = t218 * t57;
201 t220 = t219 * t48;
202 t221 = t30 * t220;
203 t222 = t117 * t117;
204 t223 = t222 * t57;
205 t224 = t223 * t48;
206 t225 = t30 * t224;
207 t226 = t118 * t109;
208 t227 = t30 * t226;
209 t229 = t57 * t210;
210 t230 = t30 * t229;
211 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
212 v2rho2[0] = 0.2e1 * t82 - 0.14e2 / 0.3e1 * t87 + 0.2e1 * t120 + 0.2e1 * t122 + rho[0] * (t168 + 0.70e2 / 0.9e1 * t173 - 0.14e2 / 0.3e1 * t175 - 0.14e2 / 0.3e1 * t177 + t221 + t225 + 0.2e1 * t227 + t230);
213
214 t233 = t85 * t57;
215 t234 = t233 * t48;
216 t236 = t29 * t117;
217 t237 = t236 * t58;
218 t238 = t125 * t109;
219 t240 = 0.1e1 / t127 / t169;
220 t241 = t50 * t240;
221 t242 = t241 * t57;
222 t244 = t130 * t118;
223 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
224 v2rhosigma[0] = t126 - t132 + rho[0] * (-0.7e1 / 0.3e1 * t234 + t237 + t238 + 0.14225894375000000000e-2 * t242 - 0.40645412500000000000e-3 * t244);
225
226 t248 = 0.1e1 / t50;
227 t249 = t129 * t248;
228 t251 = 0.60968118750000000000e-3 * t249 * t57;
229 t253 = 0.1e1 / t35 / t169;
230 t254 = t253 * t49;
231 t256 = 0.16520495572951562500e-6 * t254 * t57;
232 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
233 v2sigma2[0] = rho[0] * (-t251 + t256);
234
235 #ifndef XC_DONT_COMPILE_KXC
236
237 if(order < 3) return;
238
239
240 t266 = t86 * t224;
241 t268 = t30 * t218;
242 t269 = t268 * t119;
243 t271 = t222 * t117;
244 t272 = t271 * t57;
245 t273 = t272 * t48;
246 t274 = t30 * t273;
247 t275 = t169 * rho[0];
248 t277 = 0.1e1 / t7 / t275;
249 t278 = sigma[0] * t277;
250 t279 = t278 * t58;
251 t281 = t172 * t119;
252 t283 = t86 * t220;
253 t285 = t86 * t226;
254 t287 = t89 * t89;
255 t289 = 0.1e1 / t287 * t50;
256 t290 = t182 * t109;
257 t306 = t6 * t85;
258 t307 = t4 * t306;
259 t310 = 0.1e1 / t35 / t83;
260 t312 = t34 * t5 * t310;
261 t314 = -0.60319259259259259259e-2 * t307 - 0.54733333333333333333e-5 * t312;
262 t316 = t193 * t101;
263 t319 = t98 * t199;
264 t324 = t100 * t100;
265 t325 = 0.1e1 / t324;
266 t326 = t40 * t325;
267 t327 = t201 * t107;
268 t330 = t107 * t208;
269 t335 = 0.1e1 / t169;
270 t337 = -0.22615185185185185185e1 * t307 - 0.34962962962962962963e0 * t312 - 0.10583962870554131724e0 * t335;
271 t339 = -t102 * t337 - 0.3e1 * t316 * t107 - 0.3e1 * t195 * t208 + 0.6e1 * t200 * t330 + 0.6e1 * t319 * t201 + t314 * t46 - 0.6e1 * t326 * t327;
272 t344 = 0.1e1 / t52 / t169;
273 t347 = 0.487744950e-2 * t289 * t54 * t290 + 0.56903577500000000001e-2 * t181 * t114 * t182 - 0.487744950e-2 * t181 * t110 * t210 + 0.61645542291666666667e-2 * t91 * t215 * t109 - 0.28451788750000000000e-2 * t91 * t114 * t210 + 0.81290825e-3 * t91 * t54 * t339 + 0.65070294641203703703e-2 * t51 * t344;
274 t348 = t347 * t57;
275 t349 = t348 * t48;
276 t350 = t30 * t349;
277 t351 = t219 * t109;
278 t352 = t30 * t351;
279 t354 = t223 * t109;
280 t355 = t30 * t354;
281 t357 = t118 * t210;
282 t358 = t30 * t357;
283 t360 = t172 * t121;
284 t362 = t86 * t229;
285 t364 = t57 * t339;
286 t365 = t30 * t364;
287 t366 = t60 * t60;
288 t367 = 0.1e1 / t366;
289 t371 = t139 * t73;
290 t377 = 0.1e1 / t13 / t38 * t2 / 0.4e1;
291 t387 = -0.35096666666666666667e0 * t377 * t335 + 0.23397777777777777778e0 * t145 * t146 * t310 - 0.27297407407407407408e0 * t63 * t64 * t85 - 0.86437037037037037036e-1 * t307;
292 t397 = my_piecewise3(t12, 0.8538e0 * t367 * t140 * t73 - 0.8538e0 * t371 * t158 + 0.1423e0 * t61 * t387, -0.20733333333333333334e-1 * t206 - 0.51851851851851851854e-3 * t4 * t306 * t19 + 0.22851851851851851851e-2 * t307);
293 t398 = -0.7e1 * t266 + 0.3e1 * t269 + t274 - 0.910e3 / 0.27e2 * t279 + 0.70e2 / 0.3e1 * t281 - 0.7e1 * t283 - 0.14e2 * t285 + t350 + 0.3e1 * t352 + 0.3e1 * t355 + 0.3e1 * t358 + 0.70e2 / 0.3e1 * t360 - 0.7e1 * t362 + t365 + t397;
294 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
295 v3rho3[0] = 0.3e1 * t168 + 0.70e2 / 0.3e1 * t173 - 0.14e2 * t175 - 0.14e2 * t177 + 0.3e1 * t221 + 0.3e1 * t225 + 0.6e1 * t227 + 0.3e1 * t230 + rho[0] * t398;
296
297 t405 = t171 * t57;
298 t406 = t405 * t48;
299 t408 = t85 * t117;
300 t409 = t408 * t58;
301 t411 = t233 * t109;
302 t413 = t29 * t218;
303 t414 = t413 * t58;
304 t415 = t29 * t222;
305 t416 = t415 * t58;
306 t417 = t236 * t121;
307 t419 = t125 * t210;
308 t421 = 0.1e1 / t127 / t275;
309 t422 = t50 * t421;
310 t423 = t422 * t57;
311 t425 = t241 * t118;
312 t427 = t130 * t219;
313 t429 = t130 * t223;
314 t431 = 0.70e2 / 0.9e1 * t406 - 0.14e2 / 0.3e1 * t409 - 0.14e2 / 0.3e1 * t411 + t414 + t416 + 0.2e1 * t417 + t419 - 0.64016524687500000000e-2 * t423 + 0.28451788750000000000e-2 * t425 - 0.40645412500000000000e-3 * t427 - 0.40645412500000000000e-3 * t429;
315 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
316 v3rho2sigma[0] = -0.14e2 / 0.3e1 * t234 + 0.2e1 * t237 + 0.2e1 * t238 + 0.28451788750000000000e-2 * t242 - 0.81290825000000000000e-3 * t244 + rho[0] * t431;
317
318 t433 = t240 * t248;
319 t434 = t433 * t57;
320 t436 = t249 * t118;
321 t439 = 0.1e1 / t35 / t275;
322 t440 = t439 * t49;
323 t441 = t440 * t57;
324 t443 = t253 * t90;
325 t444 = t443 * t121;
326 t446 = t254 * t118;
327 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
328 v3rhosigma2[0] = -t251 + t256 + rho[0] * (0.21338841562500000000e-2 * t434 - 0.60968118750000000000e-3 * t436 - 0.77095646007107291667e-6 * t441 - 0.16520495572951562500e-6 * t444 + 0.16520495572951562500e-6 * t446);
329
330 t451 = 0.1e1 / t50 / sigma[0];
331 t452 = t129 * t451;
332 t454 = 0.30484059375000000000e-3 * t452 * t57;
333 t455 = 0.1e1 / sigma[0];
334 t456 = t253 * t455;
335 t457 = t49 * t57;
336 t459 = 0.24780743359427343750e-6 * t456 * t457;
337 t460 = t52 * t52;
338 t461 = t460 * t460;
339 t462 = t461 * t52;
340 t464 = 0.1e1 / t462 / t275;
341 t465 = t464 * t90;
342 t466 = t248 * t57;
343 t468 = 0.67148235726704010033e-10 * t465 * t466;
344 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
345 v3sigma3[0] = rho[0] * (t454 + t459 - t468);
346
347 #ifndef XC_DONT_COMPILE_LXC
348
349 if(order < 4) return;
350
351
352 t470 = t222 * t222;
353 t474 = t218 * t218;
354 t481 = t169 * t27;
355 t504 = t140 * t140;
356 t510 = t158 * t158;
357 t524 = 0.1e1 / t275;
358 t533 = t6 * t171;
359 t534 = t4 * t533;
360 t546 = my_piecewise3(t12, -0.34152e1 / t366 / t16 * t504 + 0.51228e1 * t367 * t140 * t158 - 0.8538e0 * t139 * t510 - 0.11384e1 * t371 * t387 + 0.1423e0 * t61 * (-0.60931712962962962962e-2 / t13 / t43 * t277 * t4 * t6 + 0.28077333333333333334e1 * t377 * t524 - 0.10399012345679012346e1 * t145 * t146 * t253 + 0.90991358024691358027e0 * t63 * t64 * t171 + 0.28812345679012345679e0 * t534), 0.62200000000000000002e-1 * t335 + 0.17283950617283950618e-2 * t4 * t533 * t19 - 0.74444444444444444441e-2 * t534);
361 t549 = t30 * t470 * t57 * t48 + 0.3e1 * t30 * t474 * t57 * t48 - 0.28e2 / 0.3e1 * t86 * t273 + 0.14560e5 / 0.81e2 * sigma[0] / t7 / t481 * t58 - 0.28e2 * t86 * t354 + 0.140e3 / 0.3e1 * t172 * t224 - 0.28e2 / 0.3e1 * t86 * t349 + 0.280e3 / 0.3e1 * t172 * t226 + 0.140e3 / 0.3e1 * t172 * t220 + 0.4e1 * t30 * t272 * t109 - 0.3640e4 / 0.27e2 * t278 * t119 + t546 - 0.28e2 * t86 * t357;
362 t553 = t182 * t182;
363 t569 = t210 * t210;
364 t587 = t34 * t5 * t253;
365 t609 = t201 * t201;
366 t615 = t208 * t208;
367 t626 = (0.20106419753086419753e-1 * t534 + 0.20068888888888888889e-4 * t587) * t46 - 0.4e1 * t314 * t101 * t107 + 0.12e2 * t193 * t199 * t201 - 0.6e1 * t316 * t208 - 0.24e2 * t98 * t325 * t327 + 0.24e2 * t319 * t330 - 0.4e1 * t195 * t337 + 0.24e2 * t40 / t324 / t45 * t609 - 0.36e2 * t326 * t201 * t208 + 0.6e1 * t200 * t615 + 0.8e1 * t200 * t107 * t337 - t102 * (0.75383950617283950617e1 * t534 + 0.12819753086419753086e1 * t587 + 0.42335851482216526896e0 * t524);
368 t634 = -0.1950979800e-1 / t287 / t48 * t50 * t54 * t553 - 0.22761431000000000000e-1 * t289 * t114 * t290 + 0.2926469700e-1 * t289 * t183 * t210 - 0.24658216916666666667e-1 * t181 * t215 * t182 + 0.22761431000000000000e-1 * t181 * t186 * t210 - 0.487744950e-2 * t181 * t54 * t569 - 0.650326600e-2 * t181 * t110 * t339 - 0.26028117856481481481e-1 * t91 * t344 * t109 + 0.12329108458333333333e-1 * t91 * t215 * t210 - 0.37935718333333333333e-2 * t91 * t114 * t339 + 0.81290825e-3 * t91 * t54 * t626 - 0.27112622767168209876e-1 * t51 / t52 / t275;
369 t670 = t30 * t634 * t57 * t48 - 0.28e2 * t86 * t351 + 0.6e1 * t30 * t223 * t210 + 0.6e1 * t30 * t219 * t210 + 0.4e1 * t30 * t348 * t109 + 0.4e1 * t30 * t118 * t339 - 0.3640e4 / 0.27e2 * t278 * t121 - 0.28e2 / 0.3e1 * t86 * t364 + 0.140e3 / 0.3e1 * t172 * t229 + t30 * t57 * t626 - 0.28e2 * t86 * t117 * t220 + 0.4e1 * t30 * t347 * t119 + 0.12e2 * t268 * t226 + 0.6e1 * t268 * t224;
370 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
371 v4rho4[0] = rho[0] * (t549 + t670) - 0.28e2 * t266 + 0.12e2 * t269 + 0.4e1 * t274 - 0.3640e4 / 0.27e2 * t279 + 0.280e3 / 0.3e1 * t281 - 0.28e2 * t283 - 0.56e2 * t285 + 0.4e1 * t350 + 0.12e2 * t352 + 0.12e2 * t355 + 0.12e2 * t358 + 0.280e3 / 0.3e1 * t360 - 0.28e2 * t362 + 0.4e1 * t365 + 0.4e1 * t397;
372
373 t739 = t29 * t347 * t58 + 0.70e2 / 0.3e1 * t171 * t117 * t58 - 0.7e1 * t85 * t218 * t58 - 0.14e2 * t408 * t121 - 0.910e3 / 0.27e2 * t277 * t57 * t48 - 0.12193623750000000000e-2 * t130 * t218 * t117 * t57 - 0.40645412500000000000e-3 * t130 * t272 + 0.3e1 * t413 * t119 + t29 * t271 * t58 + 0.42677683125000000000e-2 * t241 * t223 - 0.7e1 * t85 * t222 * t58;
374 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
375 v4rho3sigma[0] = rho[0] * (t125 * t339 - 0.7e1 * t233 * t210 + 0.70e2 / 0.3e1 * t405 * t109 + 0.42677683125000000000e-2 * t241 * t219 - 0.40645412500000000000e-3 * t130 * t348 + 0.35209088578125000000e-1 * t50 / t127 / t481 * t57 - 0.19204957406250000000e-1 * t422 * t118 + 0.3e1 * t236 * t229 + 0.3e1 * t413 * t121 + 0.3e1 * t415 * t121 + t739) - 0.14e2 * t411 + 0.3e1 * t419 + 0.70e2 / 0.3e1 * t406 - 0.14e2 * t409 + 0.3e1 * t414 + 0.6e1 * t417 - 0.19204957406250000000e-1 * t423 + 0.85355366250000000000e-2 * t425 - 0.12193623750000000000e-2 * t427 + 0.3e1 * t416 - 0.12193623750000000000e-2 * t429;
376
377 t789 = -0.96024787031250000000e-2 * t421 * t248 * t57 + 0.42677683125000000000e-2 * t433 * t118 - 0.60968118750000000000e-3 * t249 * t219 - 0.60968118750000000000e-3 * t249 * t223 + 0.43687532737360798611e-5 / t35 / t481 * t49 * t57 + 0.15419129201421458333e-5 * t439 * t90 * t121 - 0.15419129201421458333e-5 * t440 * t118 + 0.33040991145903125000e-6 * t253 * t180 * t57 * t182 - 0.33040991145903125000e-6 * t443 * t226 - 0.16520495572951562500e-6 * t443 * t229 + 0.16520495572951562500e-6 * t254 * t219 + 0.16520495572951562500e-6 * t254 * t223;
378 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
379 v4rho2sigma2[0] = 0.42677683125000000000e-2 * t434 - 0.12193623750000000000e-2 * t436 - 0.15419129201421458333e-5 * t441 - 0.33040991145903125000e-6 * t444 + 0.33040991145903125000e-6 * t446 + rho[0] * t789;
380
381 t799 = t90 * t57;
382 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
383 v4rhosigma3[0] = t454 + t459 - t468 + rho[0] * (-0.10669420781250000000e-2 * t240 * t451 * t57 + 0.30484059375000000000e-3 * t452 * t118 - 0.11564346901066093750e-5 * t439 * t455 * t457 - 0.24780743359427343750e-6 * t456 * t799 * t109 + 0.24780743359427343750e-6 * t456 * t49 * t117 * t57 + 0.39169804173910672519e-9 / t462 / t481 * t90 * t466 + 0.13429647145340802007e-9 * t464 * t180 * t466 * t109 - 0.67148235726704010033e-10 * t465 * t248 * t117 * t57);
384
385 t822 = sigma[0] * sigma[0];
386 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
387 v4sigma4[0] = rho[0] * (-0.45726089062500000000e-3 * t129 / t50 / t822 * t57 - 0.37171115039141015625e-6 * t253 / t822 * t457 - 0.67148235726704010034e-10 * t464 * t451 * t799 + 0.27292677397591217532e-13 / t169 / t83 * t180 * t455 * t57);
388
389 #ifndef XC_DONT_COMPILE_MXC
390
391 if(order < 5) return;
392
393
394 #endif
395
396 #endif
397
398 #endif
399
400 #endif
401
402 #endif
403
404
405 }
406
407
408 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 *))409 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
410 {
411
412 #ifndef XC_DONT_COMPILE_EXC
413 double t1, t2, t3, t4, t5, t6, t7, t9;
414 double t10, t11, t12, t13, t16, t19, t26, t27;
415 double t29, t30, t32, t33, t34, t35, t38, t40;
416 double t43, t45, t46, t48, t49, t50, t51, t52;
417 double t54, t57, t59, t60, t61, t63;
418
419 #ifndef XC_DONT_COMPILE_VXC
420 double t64, t65, t67, t68, t70, t74, t75, t77;
421 double t86, t87, t89, t90, t91, t93, t94, t95;
422 double t100, t102, t104, t105, t106, t109, t111, t113;
423 double t114, t118, t121, t122, t123, t126, t127, t131;
424 double t132, t134, t135, t137, t138, t139, t141;
425
426 #ifndef XC_DONT_COMPILE_FXC
427 double t146, t147, t152, t153, t155, t162, t163, t165;
428 double t175, t176, t178, t179, t180, t182, t183, t185;
429 double t188, t189, t190, t191, t194, t199, t201, t203;
430 double t207, t208, t209, t214, t216, t218, t223, t226;
431 double t227, t228, t230, t231, t232, t234, t236, t237;
432 double t241, t242, t244, t245, t247, t248, t251, t252;
433 double t253, t255, t256, t257, t261, t262, t264, t266;
434 double t267, t269;
435
436 #ifndef XC_DONT_COMPILE_KXC
437 double t279, t281, t282, t283, t285, t286, t288, t289;
438 double t291, t293, t295, t296, t312, t313, t316, t318;
439 double t320, t322, t325, t330, t331, t332, t333, t336;
440 double t341, t343, t345, t350, t353, t354, t355, t357;
441 double t359, t361, t363, t365, t368, t369, t371, t372;
442 double t374, t375, t377, t378, t379, t381, t382, t386;
443 double t392, t402, t412, t413, t418, t419, t421, t422;
444 double t424, t426, t427, t429, t430, t432, t433, t434;
445 double t437, t438, t439, t441, t444, t445, t447, t448;
446 double t449, t451, t453, t454, t456, t459, t460, t461;
447 double t463, t464, t466, t471, t472, t474, t475, t476;
448 double t478, t480, t481, t482, t483, t485, t486, t488;
449 double t490;
450
451 #ifndef XC_DONT_COMPILE_LXC
452 double t492, t493, t496, t518, t524, t532, t536, t549;
453 double t555, t588, t589, t612, t620, t636, t659, t675;
454 double t679, t693, t770, t821, t851;
455 #endif
456
457 #endif
458
459 #endif
460
461 #endif
462
463 #endif
464
465
466
467 t1 = M_CBRT3;
468 t2 = 0.1e1 / M_PI;
469 t3 = POW_1_3(t2);
470 t4 = t1 * t3;
471 t5 = M_CBRT4;
472 t6 = t5 * t5;
473 t7 = POW_1_3(rho[0]);
474 t9 = t6 / t7;
475 t10 = t4 * t9;
476 t11 = t10 / 0.4e1;
477 t12 = 0.1e1 <= t11;
478 t13 = sqrt(t10);
479 t16 = 0.1e1 + 0.69905000000000000000e0 * t13 + 0.65275000000000000000e-1 * t10;
480 t19 = log(t11);
481 t26 = my_piecewise3(t12, -0.843e-1 / t16, 0.1555e-1 * t19 - 0.269e-1 + 0.17500000000000000000e-3 * t4 * t9 * t19 - 0.12000000000000000000e-2 * t10);
482 t27 = rho[0] * rho[0];
483 t29 = 0.1e1 / t7 / t27;
484 t30 = sigma[0] * t29;
485 t32 = t1 * t1;
486 t33 = t3 * t3;
487 t34 = t32 * t33;
488 t35 = t7 * t7;
489 t38 = t34 * t5 / t35;
490 t40 = 0.2568e-2 + 0.58165000000000000000e-2 * t10 + 0.18472500000000000000e-5 * t38;
491 t43 = 0.1e1 / rho[0];
492 t45 = 0.1e1 + 0.21807500000000000000e1 * t10 + 0.11800000000000000000e0 * t38 + 0.17639938117590219540e-1 * t43;
493 t46 = 0.1e1 / t45;
494 t48 = 0.1667e-2 + t40 * t46;
495 t49 = 0.1e1 / t48;
496 t50 = sqrt(sigma[0]);
497 t51 = t49 * t50;
498 t52 = pow(rho[0], 0.1e1 / 0.6e1);
499 t54 = 0.1e1 / t52 / rho[0];
500 t57 = exp(-0.81290825e-3 * t51 * t54);
501 t59 = M_CBRT2;
502 t60 = t59 * t59;
503 t61 = t57 * t48 * t60;
504 t63 = t30 * t61 / 0.2e1;
505 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
506 zk[0] = t26 + t63;
507
508 #ifndef XC_DONT_COMPILE_VXC
509
510 if(order < 1) return;
511
512
513 t64 = t16 * t16;
514 t65 = 0.1e1 / t64;
515 t67 = 0.1e1 / t13 * t1;
516 t68 = t3 * t6;
517 t70 = 0.1e1 / t7 / rho[0];
518 t74 = t6 * t70;
519 t75 = t4 * t74;
520 t77 = -0.11650833333333333333e0 * t67 * t68 * t70 - 0.21758333333333333333e-1 * t75;
521 t86 = my_piecewise3(t12, 0.843e-1 * t65 * t77, -0.51833333333333333333e-2 * t43 - 0.58333333333333333333e-4 * t4 * t74 * t19 + 0.34166666666666666667e-3 * t75);
522 t87 = t27 * rho[0];
523 t89 = 0.1e1 / t7 / t87;
524 t90 = sigma[0] * t89;
525 t91 = t90 * t61;
526 t93 = t48 * t48;
527 t94 = 0.1e1 / t93;
528 t95 = t94 * t50;
529 t100 = t34 * t5 / t35 / rho[0];
530 t102 = -0.19388333333333333333e-2 * t75 - 0.12315000000000000000e-5 * t100;
531 t104 = t45 * t45;
532 t105 = 0.1e1 / t104;
533 t106 = t40 * t105;
534 t109 = 0.1e1 / t27;
535 t111 = -0.72691666666666666667e0 * t75 - 0.78666666666666666667e-1 * t100 - 0.17639938117590219540e-1 * t109;
536 t113 = t102 * t46 - t106 * t111;
537 t114 = t54 * t113;
538 t118 = 0.1e1 / t52 / t27;
539 t121 = 0.81290825e-3 * t95 * t114 + 0.94839295833333333334e-3 * t51 * t118;
540 t122 = t30 * t121;
541 t123 = t122 * t61;
542 t126 = t57 * t113 * t60;
543 t127 = t30 * t126;
544 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
545 vrho[0] = t26 + t63 + rho[0] * (t86 - 0.7e1 / 0.6e1 * t91 + t123 / 0.2e1 + t127 / 0.2e1);
546
547 t131 = t29 * t57;
548 t132 = t48 * t60;
549 t134 = t131 * t132 / 0.2e1;
550 t135 = sqrt(rho[0]);
551 t137 = 0.1e1 / t135 / t87;
552 t138 = t50 * t137;
553 t139 = t57 * t60;
554 t141 = 0.20322706250000000000e-3 * t138 * t139;
555 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
556 vsigma[0] = rho[0] * (t134 - t141);
557
558 #ifndef XC_DONT_COMPILE_FXC
559
560 if(order < 2) return;
561
562
563 t146 = 0.1e1 / t64 / t16;
564 t147 = t77 * t77;
565 t152 = 0.1e1 / t13 / t10 * t32;
566 t153 = t33 * t5;
567 t155 = 0.1e1 / t35 / t27;
568 t162 = t6 * t29;
569 t163 = t4 * t162;
570 t165 = -0.77672222222222222220e-1 * t152 * t153 * t155 + 0.15534444444444444444e0 * t67 * t68 * t29 + 0.29011111111111111111e-1 * t163;
571 t175 = my_piecewise3(t12, -0.1686e0 * t146 * t147 + 0.843e-1 * t65 * t165, 0.51833333333333333333e-2 * t109 + 0.77777777777777777777e-4 * t4 * t162 * t19 - 0.43611111111111111112e-3 * t163);
572 t176 = t27 * t27;
573 t178 = 0.1e1 / t7 / t176;
574 t179 = sigma[0] * t178;
575 t180 = t179 * t61;
576 t182 = t90 * t121;
577 t183 = t182 * t61;
578 t185 = t90 * t126;
579 t188 = 0.1e1 / t93 / t48;
580 t189 = t188 * t50;
581 t190 = t113 * t113;
582 t191 = t54 * t190;
583 t194 = t118 * t113;
584 t199 = t34 * t5 * t155;
585 t201 = 0.25851111111111111111e-2 * t163 + 0.20525000000000000000e-5 * t199;
586 t203 = t102 * t105;
587 t207 = 0.1e1 / t104 / t45;
588 t208 = t40 * t207;
589 t209 = t111 * t111;
590 t214 = 0.1e1 / t87;
591 t216 = 0.96922222222222222222e0 * t163 + 0.13111111111111111111e0 * t199 + 0.35279876235180439080e-1 * t214;
592 t218 = -t106 * t216 - 0.2e1 * t203 * t111 + t201 * t46 + 0.2e1 * t208 * t209;
593 t223 = 0.1e1 / t52 / t87;
594 t226 = -0.162581650e-2 * t189 * t191 - 0.18967859166666666667e-2 * t95 * t194 + 0.81290825e-3 * t95 * t54 * t218 - 0.20548514097222222222e-2 * t51 * t223;
595 t227 = t30 * t226;
596 t228 = t227 * t61;
597 t230 = t121 * t121;
598 t231 = t30 * t230;
599 t232 = t231 * t61;
600 t234 = t122 * t126;
601 t236 = t57 * t218 * t60;
602 t237 = t30 * t236;
603 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
604 v2rho2[0] = 0.2e1 * t86 - 0.7e1 / 0.3e1 * t91 + t123 + t127 + rho[0] * (t175 + 0.35e2 / 0.9e1 * t180 - 0.7e1 / 0.3e1 * t183 - 0.7e1 / 0.3e1 * t185 + t228 / 0.2e1 + t232 / 0.2e1 + t234 + t237 / 0.2e1);
605
606 t241 = t89 * t57;
607 t242 = t241 * t132;
608 t244 = t29 * t121;
609 t245 = t244 * t61;
610 t247 = t113 * t60;
611 t248 = t131 * t247;
612 t251 = 0.1e1 / t135 / t176;
613 t252 = t50 * t251;
614 t253 = t252 * t139;
615 t255 = t121 * t57;
616 t256 = t255 * t60;
617 t257 = t138 * t256;
618 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
619 v2rhosigma[0] = t134 - t141 + rho[0] * (-0.7e1 / 0.6e1 * t242 + t245 / 0.2e1 + t248 / 0.2e1 + 0.71129471875000000000e-3 * t253 - 0.20322706250000000000e-3 * t257);
620
621 t261 = 0.1e1 / t50;
622 t262 = t137 * t261;
623 t264 = 0.30484059375000000000e-3 * t262 * t139;
624 t266 = 0.1e1 / t35 / t176;
625 t267 = t266 * t49;
626 t269 = 0.82602477864757812500e-7 * t267 * t139;
627 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
628 v2sigma2[0] = rho[0] * (-t264 + t269);
629
630 #ifndef XC_DONT_COMPILE_KXC
631
632 if(order < 3) return;
633
634
635 t279 = t176 * rho[0];
636 t281 = 0.1e1 / t7 / t279;
637 t282 = sigma[0] * t281;
638 t283 = t282 * t61;
639 t285 = t179 * t121;
640 t286 = t285 * t61;
641 t288 = t90 * t226;
642 t289 = t288 * t61;
643 t291 = t182 * t126;
644 t293 = t93 * t93;
645 t295 = 0.1e1 / t293 * t50;
646 t296 = t190 * t113;
647 t312 = t6 * t89;
648 t313 = t4 * t312;
649 t316 = 0.1e1 / t35 / t87;
650 t318 = t34 * t5 * t316;
651 t320 = -0.60319259259259259259e-2 * t313 - 0.54733333333333333333e-5 * t318;
652 t322 = t201 * t105;
653 t325 = t102 * t207;
654 t330 = t104 * t104;
655 t331 = 0.1e1 / t330;
656 t332 = t40 * t331;
657 t333 = t209 * t111;
658 t336 = t111 * t216;
659 t341 = 0.1e1 / t176;
660 t343 = -0.22615185185185185185e1 * t313 - 0.34962962962962962963e0 * t318 - 0.10583962870554131724e0 * t341;
661 t345 = -t106 * t343 - 0.3e1 * t322 * t111 - 0.3e1 * t203 * t216 + 0.6e1 * t208 * t336 + 0.6e1 * t325 * t209 + t320 * t46 - 0.6e1 * t332 * t333;
662 t350 = 0.1e1 / t52 / t176;
663 t353 = 0.487744950e-2 * t295 * t54 * t296 + 0.56903577500000000001e-2 * t189 * t118 * t190 - 0.487744950e-2 * t189 * t114 * t218 + 0.61645542291666666667e-2 * t95 * t223 * t113 - 0.28451788750000000000e-2 * t95 * t118 * t218 + 0.81290825e-3 * t95 * t54 * t345 + 0.65070294641203703703e-2 * t51 * t350;
664 t354 = t30 * t353;
665 t355 = t354 * t61;
666 t357 = t227 * t126;
667 t359 = t231 * t126;
668 t361 = t122 * t236;
669 t363 = t179 * t126;
670 t365 = t90 * t236;
671 t368 = t57 * t345 * t60;
672 t369 = t30 * t368;
673 t371 = t90 * t230;
674 t372 = t371 * t61;
675 t374 = t255 * t132;
676 t375 = t227 * t374;
677 t377 = t230 * t121;
678 t378 = t30 * t377;
679 t379 = t378 * t61;
680 t381 = t64 * t64;
681 t382 = 0.1e1 / t381;
682 t386 = t146 * t77;
683 t392 = 0.1e1 / t13 / t38 * t2 / 0.4e1;
684 t402 = -0.46603333333333333332e0 * t392 * t341 + 0.31068888888888888888e0 * t152 * t153 * t316 - 0.36247037037037037036e0 * t67 * t68 * t89 - 0.67692592592592592592e-1 * t313;
685 t412 = my_piecewise3(t12, 0.5058e0 * t382 * t147 * t77 - 0.5058e0 * t386 * t165 + 0.843e-1 * t65 * t402, -0.10366666666666666667e-1 * t214 - 0.18148148148148148148e-3 * t4 * t312 * t19 + 0.99166666666666666667e-3 * t313);
686 t413 = -0.455e3 / 0.27e2 * t283 + 0.35e2 / 0.3e1 * t286 - 0.7e1 / 0.2e1 * t289 - 0.7e1 * t291 + t355 / 0.2e1 + 0.3e1 / 0.2e1 * t357 + 0.3e1 / 0.2e1 * t359 + 0.3e1 / 0.2e1 * t361 + 0.35e2 / 0.3e1 * t363 - 0.7e1 / 0.2e1 * t365 + t369 / 0.2e1 - 0.7e1 / 0.2e1 * t372 + 0.3e1 / 0.2e1 * t375 + t379 / 0.2e1 + t412;
687 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
688 v3rho3[0] = 0.3e1 * t175 + 0.35e2 / 0.3e1 * t180 - 0.7e1 * t183 - 0.7e1 * t185 + 0.3e1 / 0.2e1 * t228 + 0.3e1 / 0.2e1 * t232 + 0.3e1 * t234 + 0.3e1 / 0.2e1 * t237 + rho[0] * t413;
689
690 t418 = t178 * t57;
691 t419 = t418 * t132;
692 t421 = t89 * t121;
693 t422 = t421 * t61;
694 t424 = t241 * t247;
695 t426 = t29 * t226;
696 t427 = t426 * t61;
697 t429 = t29 * t230;
698 t430 = t429 * t61;
699 t432 = t244 * t126;
700 t433 = t218 * t60;
701 t434 = t131 * t433;
702 t437 = 0.1e1 / t135 / t279;
703 t438 = t50 * t437;
704 t439 = t438 * t139;
705 t441 = t252 * t256;
706 t444 = t226 * t57 * t60;
707 t445 = t138 * t444;
708 t447 = t230 * t57;
709 t448 = t447 * t60;
710 t449 = t138 * t448;
711 t451 = 0.35e2 / 0.9e1 * t419 - 0.7e1 / 0.3e1 * t422 - 0.7e1 / 0.3e1 * t424 + t427 / 0.2e1 + t430 / 0.2e1 + t432 + t434 / 0.2e1 - 0.32008262343750000000e-2 * t439 + 0.14225894375000000000e-2 * t441 - 0.20322706250000000000e-3 * t445 - 0.20322706250000000000e-3 * t449;
712 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
713 v3rho2sigma[0] = -0.7e1 / 0.3e1 * t242 + t245 + t248 + 0.14225894375000000000e-2 * t253 - 0.40645412500000000000e-3 * t257 + rho[0] * t451;
714
715 t453 = t251 * t261;
716 t454 = t453 * t139;
717 t456 = t262 * t256;
718 t459 = 0.1e1 / t35 / t279;
719 t460 = t459 * t49;
720 t461 = t460 * t139;
721 t463 = t266 * t94;
722 t464 = t463 * t126;
723 t466 = t267 * t256;
724 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
725 v3rhosigma2[0] = -t264 + t269 + rho[0] * (0.10669420781250000000e-2 * t454 - 0.30484059375000000000e-3 * t456 - 0.38547823003553645833e-6 * t461 - 0.82602477864757812500e-7 * t464 + 0.82602477864757812500e-7 * t466);
726
727 t471 = 0.1e1 / t50 / sigma[0];
728 t472 = t137 * t471;
729 t474 = 0.15242029687500000000e-3 * t472 * t139;
730 t475 = 0.1e1 / sigma[0];
731 t476 = t266 * t475;
732 t478 = t49 * t57 * t60;
733 t480 = 0.12390371679713671875e-6 * t476 * t478;
734 t481 = t52 * t52;
735 t482 = t481 * t481;
736 t483 = t482 * t52;
737 t485 = 0.1e1 / t483 / t279;
738 t486 = t485 * t94;
739 t488 = t261 * t57 * t60;
740 t490 = 0.33574117863352005017e-10 * t486 * t488;
741 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
742 v3sigma3[0] = rho[0] * (t474 + t480 - t490);
743
744 #ifndef XC_DONT_COMPILE_LXC
745
746 if(order < 4) return;
747
748
749 t492 = t6 * t178;
750 t493 = t4 * t492;
751 t496 = t34 * t5 * t266;
752 t518 = t209 * t209;
753 t524 = t216 * t216;
754 t532 = 0.1e1 / t279;
755 t536 = (0.20106419753086419753e-1 * t493 + 0.20068888888888888889e-4 * t496) * t46 - 0.4e1 * t320 * t105 * t111 + 0.12e2 * t201 * t207 * t209 - 0.6e1 * t322 * t216 - 0.24e2 * t102 * t331 * t333 + 0.24e2 * t325 * t336 - 0.4e1 * t203 * t343 + 0.24e2 * t40 / t330 / t45 * t518 - 0.36e2 * t332 * t209 * t216 + 0.6e1 * t208 * t524 + 0.8e1 * t208 * t111 * t343 - t106 * (0.75383950617283950617e1 * t493 + 0.12819753086419753086e1 * t496 + 0.42335851482216526896e0 * t532);
756 t549 = t147 * t147;
757 t555 = t165 * t165;
758 t588 = my_piecewise3(t12, -0.20232e1 / t381 / t16 * t549 + 0.30348e1 * t382 * t147 * t165 - 0.5058e0 * t146 * t555 - 0.6744e0 * t386 * t402 + 0.843e-1 * t65 * (-0.80908564814814814812e-2 / t13 / t43 * t281 * t4 * t6 + 0.37282666666666666666e1 * t392 * t532 - 0.13808395061728395062e1 * t152 * t153 * t266 + 0.12082345679012345679e1 * t67 * t68 * t178 + 0.22564197530864197531e0 * t493), 0.31100000000000000001e-1 * t341 + 0.60493827160493827160e-3 * t4 * t492 * t19 - 0.32450617283950617284e-2 * t493);
759 t589 = t176 * t27;
760 t612 = t30 * t57 * t536 * t60 / 0.2e1 - 0.1820e4 / 0.27e2 * t282 * t126 - 0.14e2 / 0.3e1 * t90 * t368 + 0.70e2 / 0.3e1 * t179 * t236 + t588 + 0.7280e4 / 0.81e2 * sigma[0] / t7 / t589 * t61 + 0.3e1 * t227 * t236 + 0.3e1 * t231 * t236 + 0.2e1 * t122 * t368 - 0.1820e4 / 0.27e2 * t282 * t121 * t61 + 0.70e2 / 0.3e1 * t179 * t226 * t61 + 0.140e3 / 0.3e1 * t285 * t126 - 0.14e2 / 0.3e1 * t90 * t353 * t61;
761 t620 = t190 * t190;
762 t636 = t218 * t218;
763 t659 = -0.1950979800e-1 / t293 / t48 * t50 * t54 * t620 - 0.22761431000000000000e-1 * t295 * t118 * t296 + 0.2926469700e-1 * t295 * t191 * t218 - 0.24658216916666666667e-1 * t189 * t223 * t190 + 0.22761431000000000000e-1 * t189 * t194 * t218 - 0.487744950e-2 * t189 * t54 * t636 - 0.650326600e-2 * t189 * t114 * t345 - 0.26028117856481481481e-1 * t95 * t350 * t113 + 0.12329108458333333333e-1 * t95 * t223 * t218 - 0.37935718333333333333e-2 * t95 * t118 * t345 + 0.81290825e-3 * t95 * t54 * t536 - 0.27112622767168209876e-1 * t51 / t52 / t279;
764 t675 = t226 * t226;
765 t679 = t230 * t230;
766 t693 = -0.14e2 * t288 * t126 - 0.14e2 * t182 * t236 + t30 * t659 * t61 / 0.2e1 + 0.2e1 * t354 * t126 + 0.70e2 / 0.3e1 * t179 * t230 * t61 - 0.14e2 * t371 * t126 + 0.2e1 * t378 * t126 - 0.14e2 / 0.3e1 * t90 * t377 * t61 + 0.3e1 / 0.2e1 * t30 * t675 * t61 + t30 * t679 * t61 / 0.2e1 + 0.3e1 * t227 * t447 * t132 - 0.14e2 * t288 * t374 + 0.2e1 * t354 * t374 + 0.6e1 * t227 * t255 * t247;
767 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
768 v4rho4[0] = rho[0] * (t612 + t693) - 0.1820e4 / 0.27e2 * t283 + 0.140e3 / 0.3e1 * t286 - 0.14e2 * t289 - 0.28e2 * t291 + 0.2e1 * t355 + 0.6e1 * t357 + 0.6e1 * t359 + 0.6e1 * t361 + 0.140e3 / 0.3e1 * t363 - 0.14e2 * t365 + 0.2e1 * t369 - 0.14e2 * t372 + 0.6e1 * t375 + 0.2e1 * t379 + 0.4e1 * t412;
769
770 t770 = -0.96024787031250000000e-2 * t438 * t256 + 0.21338841562500000000e-2 * t252 * t444 - 0.20322706250000000000e-3 * t138 * t353 * t57 * t60 - 0.7e1 * t421 * t126 + t29 * t353 * t61 / 0.2e1 + 0.3e1 / 0.2e1 * t426 * t126 + 0.3e1 / 0.2e1 * t429 * t126 + 0.3e1 / 0.2e1 * t244 * t236 - 0.455e3 / 0.27e2 * t281 * t57 * t132 + 0.35e2 / 0.3e1 * t178 * t121 * t61 - 0.7e1 / 0.2e1 * t89 * t226 * t61;
771 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
772 v4rho3sigma[0] = rho[0] * (t131 * t345 * t60 / 0.2e1 - 0.7e1 / 0.2e1 * t241 * t433 + 0.35e2 / 0.3e1 * t418 * t247 + 0.3e1 / 0.2e1 * t426 * t121 * t61 + t29 * t377 * t61 / 0.2e1 + 0.21338841562500000000e-2 * t252 * t448 - 0.60968118750000000000e-3 * t138 * t226 * t256 - 0.20322706250000000000e-3 * t138 * t377 * t57 * t60 - 0.7e1 / 0.2e1 * t89 * t230 * t61 + 0.17604544289062500000e-1 * t50 / t135 / t589 * t139 + t770) + 0.3e1 / 0.2e1 * t430 - 0.60968118750000000000e-3 * t449 - 0.7e1 * t424 + 0.3e1 / 0.2e1 * t434 + 0.35e2 / 0.3e1 * t419 - 0.7e1 * t422 + 0.3e1 / 0.2e1 * t427 + 0.3e1 * t432 - 0.96024787031250000000e-2 * t439 + 0.42677683125000000000e-2 * t441 - 0.60968118750000000000e-3 * t445;
773
774 t821 = -0.48012393515625000000e-2 * t437 * t261 * t139 + 0.21338841562500000000e-2 * t453 * t256 - 0.30484059375000000000e-3 * t262 * t444 - 0.30484059375000000000e-3 * t262 * t448 + 0.21843766368680399305e-5 / t35 / t589 * t49 * t139 + 0.77095646007107291666e-6 * t459 * t94 * t126 - 0.77095646007107291666e-6 * t460 * t256 + 0.16520495572951562500e-6 * t266 * t188 * t139 * t190 - 0.16520495572951562500e-6 * t463 * t121 * t126 - 0.82602477864757812500e-7 * t463 * t236 + 0.82602477864757812500e-7 * t267 * t444 + 0.82602477864757812500e-7 * t267 * t448;
775 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
776 v4rho2sigma2[0] = 0.21338841562500000000e-2 * t454 - 0.60968118750000000000e-3 * t456 - 0.77095646007107291666e-6 * t461 - 0.16520495572951562500e-6 * t464 + 0.16520495572951562500e-6 * t466 + rho[0] * t821;
777
778 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
779 v4rhosigma3[0] = t474 + t480 - t490 + rho[0] * (-0.53347103906250000000e-3 * t251 * t471 * t139 + 0.15242029687500000000e-3 * t472 * t256 - 0.57821734505330468750e-6 * t459 * t475 * t478 - 0.12390371679713671875e-6 * t476 * t94 * t126 + 0.12390371679713671875e-6 * t476 * t49 * t256 + 0.19584902086955336260e-9 / t483 / t589 * t94 * t488 + 0.67148235726704010034e-10 * t485 * t188 * t261 * t126 - 0.33574117863352005017e-10 * t486 * t261 * t256);
780
781 t851 = sigma[0] * sigma[0];
782 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
783 v4sigma4[0] = rho[0] * (-0.22863044531250000000e-3 * t137 / t50 / t851 * t139 - 0.18585557519570507812e-6 * t266 / t851 * t478 - 0.33574117863352005017e-10 * t485 * t471 * t94 * t57 * t60 + 0.13646338698795608766e-13 / t176 / t87 * t188 * t475 * t57 * t60);
784
785 #ifndef XC_DONT_COMPILE_MXC
786
787 if(order < 5) return;
788
789
790 #endif
791
792 #endif
793
794 #endif
795
796 #endif
797
798 #endif
799
800
801 }
802
803
804 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 *))805 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
806 {
807
808 #ifndef XC_DONT_COMPILE_EXC
809 double t1, t2, t3, t4, t5, t6, t7, t8;
810 double t10, t11, t12, t13, t14, t17, t20, t23;
811 double t27, t30, t37, t38, t39, t40, t41, t42;
812 double t43, t44, t45, t46, t47, t48, t50, t53;
813 double t54, t56, t57, t59, t60, t62, t63, t64;
814 double t65, t68, t70, t74, t75, t77, t78, t79;
815 double t80, t81, t83, t86, t87, t88, t89, t90;
816 double t91, t92, t93, t94, t95, t96, t97, t98;
817
818 #ifndef XC_DONT_COMPILE_VXC
819 double t99, t100, t102, t103, t105, t107, t109, t110;
820 double t112, t117, t121, t122, t123, t126, t133, t134;
821 double t136, t137, t138, t139, t141, t144, t146, t147;
822 double t149, t150, t151, t152, t153, t154, t155, t156;
823 double t161, t163, t165, t166, t167, t171, t173, t174;
824 double t178, t181, t182, t183, t184, t185, t186, t188;
825 double t189, t191, t192, t196, t197, t198, t199, t200;
826 double t203, t205, t208, t210, t214, t215, t216, t217;
827 double t218, t221, t222, t223, t225, t226, t228, t229;
828 double t230, t232, t233;
829
830 #ifndef XC_DONT_COMPILE_FXC
831 double t235, t236, t238, t239, t240, t242, t243, t248;
832 double t249, t251, t253, t256, t258, t259, t261, t267;
833 double t271, t273, t274, t280, t288, t289, t291, t293;
834 double t294, t295, t296, t299, t300, t302, t305, t306;
835 double t309, t312, t314, t315, t317, t318, t319, t320;
836 double t321, t322, t323, t324, t325, t326, t327, t328;
837 double t330, t331, t332, t333, t336, t341, t343, t345;
838 double t349, t350, t351, t357, t359, t364, t367, t368;
839 double t369, t370, t371, t372, t373, t374, t375, t376;
840 double t377, t378, t379, t381, t382, t383, t384, t385;
841 double t386, t388, t389, t390, t391, t392, t393, t394;
842 double t395, t400, t405, t406, t407, t408, t409, t410;
843 double t413, t414, t417, t420, t423, t426, t428, t432;
844 double t434, t436, t437, t439, t440, t441, t442, t443;
845 double t445, t448, t451, t454, t457, t458, t459, t460;
846 double t462, t465, t466, t470, t473, t476, t479, t481;
847 double t482, t483, t484, t485, t486, t487, t496, t497;
848 double t498, t499, t500, t501, t503, t504, t505, t506;
849 double t507, t508, t509, t510, t511, t512, t513, t514;
850 double t516, t517, t518, t519, t520, t521, t522, t523;
851 double t524, t525, t528, t529, t530, t531, t532, t533;
852 double t536, t537, t538, t539, t540, t543, t546, t547;
853 double t548, t549, t551, t552, t553, t554, t556, t557;
854 double t559, t560;
855
856 #ifndef XC_DONT_COMPILE_KXC
857 double t562, t564, t565, t566, t567, t568, t569, t570;
858 double t573, t575, t577, t579, t582, t587, t588, t589;
859 double t590, t593, t598, t600, t602, t604, t605, t606;
860 double t607, t608, t609, t610, t611, t613, t614, t617;
861 double t620, t622, t625, t626, t629, t632, t635, t636;
862 double t637, t638, t639, t640, t641, t642, t643, t645;
863 double t646, t647, t648, t649, t650, t654, t655, t656;
864 double t657, t658, t659, t660, t661, t662, t663, t664;
865 double t665, t666, t667, t668, t669, t670, t671, t672;
866 double t673, t675, t676, t678, t679, t699, t702, t703;
867 double t704, t705, t706, t707, t708, t709, t710, t711;
868 double t712, t713, t714, t715, t716, t717, t718, t722;
869 double t728, t729, t732, t735, t738, t744, t748, t749;
870 double t750, t751, t752, t754, t755, t756, t757, t758;
871 double t759, t760, t761, t762, t763, t766, t767, t768;
872 double t769, t771, t772, t773, t774, t775, t779, t786;
873 double t794, t795, t797, t799, t800, t802, t804, t807;
874 double t812, t815, t820, t822, t823, t827, t829, t832;
875 double t833, t834, t836, t837, t838, t841, t842, t843;
876 double t844, t845, t846, t848, t849, t850, t851, t853;
877 double t854, t855, t858, t859, t868, t871, t872, t881;
878 double t882, t883, t884, t886, t889, t895, t897, t899;
879 double t901, t902, t904, t905, t906, t907, t909, t911;
880 double t912, t913, t914, t915, t917, t918, t920, t921;
881 double t922, t924, t931, t933, t934, t935, t938, t947;
882 double t950, t959, t961, t962, t967, t969, t970, t973;
883 double t974, t975, t977, t978, t980, t985, t990, t993;
884 double t998, t1001, t1004, t1005, t1006, t1007, t1009, t1011;
885 double t1014, t1018, t1022, t1023, t1026, t1027, t1029, t1032;
886 double t1033, t1035, t1039, t1045, t1046, t1051, t1056, t1061;
887 double t1066, t1068, t1069, t1076, t1078, t1079, t1085, t1088;
888 double t1093, t1096, t1097, t1098, t1099, t1100, t1101, t1102;
889 double t1103, t1104, t1105, t1110, t1111, t1112, t1113, t1114;
890 double t1115, t1116, t1117, t1120, t1121, t1122, t1123, t1124;
891 double t1138, t1140, t1141, t1151, t1152, t1153, t1154, t1155;
892 double t1156, t1157, t1158, t1159, t1160, t1161, t1162, t1163;
893 double t1164, t1165, t1166, t1167, t1168, t1169, t1170, t1172;
894 double t1173, t1174, t1175, t1176, t1177, t1178, t1180, t1181;
895 double t1182, t1183, t1184, t1185, t1186, t1187, t1188, t1189;
896 double t1190, t1191, t1193, t1194, t1195, t1196, t1197, t1198;
897 double t1199, t1200, t1201, t1202, t1203, t1204, t1205, t1206;
898 double t1207, t1209, t1210, t1211, t1213, t1214, t1216, t1217;
899 double t1218, t1219, t1220, t1221, t1222, t1223, t1224, t1225;
900 double t1226, t1227, t1228, t1229, t1230, t1231, t1232, t1233;
901 double t1234, t1235, t1237, t1239, t1241, t1242, t1244, t1246;
902 double t1251, t1252, t1253, t1255, t1257, t1261, t1263, t1266;
903 double t1267, t1268, t1269, t1270, t1271, t1272, t1273, t1276;
904 double t1277, t1278, t1279, t1280, t1281, t1282, t1283, t1284;
905 double t1285, t1286, t1290, t1291, t1292, t1293, t1294, t1295;
906 double t1296, t1297, t1298, t1300, t1301, t1302, t1303, t1304;
907 double t1305, t1306, t1307, t1308, t1310, t1311, t1312, t1313;
908 double t1314, t1315, t1316, t1317, t1318, t1319, t1320, t1321;
909 double t1322, t1323, t1326, t1327, t1328, t1329, t1330, t1331;
910 double t1332, t1335, t1336, t1337, t1338, t1339, t1340, t1341;
911 double t1344, t1345, t1346, t1347, t1350, t1351, t1354, t1355;
912 double t1359, t1360, t1361, t1362, t1363, t1364, t1365, t1366;
913 double t1367, t1368, t1369, t1370, t1372, t1374, t1375, t1376;
914 double t1378, t1379, t1380, t1382, t1383, t1384, t1386, t1387;
915 double t1388;
916
917 #ifndef XC_DONT_COMPILE_LXC
918 double t1390, t1391, t1392, t1393, t1396, t1398, t1402, t1403;
919 double t1404, t1407, t1408, t1409, t1410, t1411, t1413, t1415;
920 double t1422, t1427, t1431, t1435, t1441, t1453, t1455, t1456;
921 double t1459, t1462, t1464, t1465, t1473, t1477, t1479, t1480;
922 double t1483, t1485, t1488, t1490, t1492, t1493, t1499, t1504;
923 double t1506, t1509, t1511, t1512, t1518, t1523, t1531, t1533;
924 double t1537, t1540, t1541, t1543, t1546, t1548, t1550, t1551;
925 double t1552, t1556, t1559, t1561, t1565, t1581, t1599, t1621;
926 double t1627, t1638, t1646, t1648, t1649, t1653, t1655, t1656;
927 double t1659, t1660, t1662, t1664, t1666, t1671, t1672, t1674;
928 double t1676, t1678, t1680, t1681, t1682, t1687, t1692, t1695;
929 double t1696, t1698, t1699, t1701, t1703, t1709, t1711, t1712;
930 double t1714, t1716, t1721, t1722, t1723, t1725, t1727, t1728;
931 double t1730, t1735, t1736, t1739, t1740, t1742, t1744, t1746;
932 double t1748, t1750, t1751, t1753, t1755, t1760, t1762, t1763;
933 double t1765, t1767, t1774, t1777, t1790, t1806, t1809, t1813;
934 double t1819, t1837, t1840, t1841, t1847, t1850, t1852, t1854;
935 double t1856, t1857, t1861, t1862, t1863, t1864, t1867, t1873;
936 double t1877, t1878, t1879, t1883, t1891, t1895, t1897, t1898;
937 double t1899, t1902, t1905, t1906, t1908, t1910, t1913, t1917;
938 double t1921, t1925, t1931, t1935, t1939, t1943, t1949, t1960;
939 double t1965, t1967, t1969, t1970, t1971, t1980, t1988, t1989;
940 double t1994, t1996, t1997, t1998, t2004, t2006, t2008, t2018;
941 double t2020, t2022, t2024, t2026, t2028, t2031, t2032, t2034;
942 double t2035, t2037, t2038, t2039, t2041, t2043, t2044, t2046;
943 double t2047, t2048, t2050, t2052, t2053, t2058, t2060, t2064;
944 double t2078, t2087, t2090, t2091, t2094, t2102, t2104, t2105;
945 double t2106, t2107, t2108, t2112, t2127, t2142, t2156, t2160;
946 double t2161, t2163, t2166, t2168, t2169, t2170, t2175, t2179;
947 double t2185, t2190, t2192, t2211, t2231, t2235, t2237, t2240;
948 double t2243, t2288, t2293, t2295, t2297, t2302, t2304, t2306;
949 double t2310, t2320, t2327, t2345, t2359, t2372, t2374, t2380;
950 double t2382, t2384, t2386, t2388, t2398, t2400, t2402, t2404;
951 double t2406, t2408, t2413, t2418, t2424, t2425, t2435, t2443;
952 double t2447, t2449, t2462, t2469, t2485, t2499, t2507, t2510;
953 double t2521, t2534, t2537, t2544, t2547, t2557, t2565, t2568;
954 double t2576, t2579, t2585, t2587, t2589, t2594, t2605, t2609;
955 double t2618, t2630, t2643, t2652, t2662, t2668, t2673, t2679;
956 double t2682, t2715, t2721, t2726, t2737, t2738, t2743, t2749;
957 double t2752, t2757, t2762, t2795, t2800, t2816, t2820, t2835;
958 double t2855, t2869, t2873, t2874, t2875, t2876, t2879, t2881;
959 double t2883, t2885, t2886, t2887, t2888, t2890, t2892, t2894;
960 double t2896, t2898, t2900, t2901, t2903, t2905, t2907, t2910;
961 double t2911, t2912, t2913, t2914, t2915, t2916, t2917, t2919;
962 double t2920, t2921, t2922, t2925, t2926, t2927, t2928, t2930;
963 double t2932, t2933, t2935, t2939, t2940, t2943, t2944, t2945;
964 double t2946, t2949, t2950, t2952, t2954, t2955, t2957, t2960;
965 double t2961, t2963, t2964, t2965, t2966, t2968, t2970, t2971;
966 double t2974, t2975, t2979, t2980, t2981, t2984, t2985, t2986;
967 double t2987, t2989, t2990, t2992, t2993, t2995, t2996, t2997;
968 double t3002, t3003, t3005, t3006, t3009, t3012, t3013, t3014;
969 double t3015, t3017, t3020, t3022, t3023, t3024, t3029, t3035;
970 double t3040, t3041, t3042, t3044, t3046, t3048, t3052, t3053;
971 double t3054, t3055, t3060, t3061, t3062, t3065, t3066, t3067;
972 double t3069, t3071, t3072, t3073, t3074, t3079, t3080, t3082;
973 double t3083, t3086, t3089, t3090, t3091, t3092, t3094, t3097;
974 double t3099, t3100, t3101, t3103, t3105, t3107, t3109, t3111;
975 double t3113, t3115, t3117, t3118, t3119, t3121, t3123, t3124;
976 double t3125, t3128, t3129, t3130, t3131, t3132, t3133, t3134;
977 double t3136, t3138, t3139, t3141, t3142, t3144, t3145, t3146;
978 double t3149, t3151, t3152, t3154, t3155, t3157, t3159, t3160;
979 double t3162, t3164, t3166, t3167, t3168, t3170, t3171, t3172;
980 double t3174, t3176, t3181, t3182, t3186, t3187, t3188, t3192;
981 double t3193, t3196, t3197, t3200, t3202, t3204, t3205, t3207;
982 double t3208, t3209, t3210, t3212, t3213, t3214, t3225, t3228;
983 double t3229, t3230, t3231, t3232, t3233, t3237, t3244, t3248;
984 double t3249, t3250, t3251, t3254, t3256, t3258, t3259, t3261;
985 double t3263, t3265, t3267, t3269, t3272, t3280, t3282, t3284;
986 double t3286, t3288, t3290, t3291, t3293, t3295, t3297, t3300;
987 double t3304, t3308, t3312, t3325, t3332, t3344, t3345, t3349;
988 double t3360, t3362, t3363, t3365, t3367, t3368, t3370, t3375;
989 double t3384, t3385, t3387, t3389, t3403, t3414, t3422, t3431;
990 double t3434, t3448, t3449, t3450, t3452, t3453, t3454, t3459;
991 double t3460, t3461, t3462, t3464, t3465, t3466, t3467, t3469;
992 double t3470, t3471, t3473, t3477, t3478, t3479, t3480, t3483;
993 double t3484, t3485, t3486, t3488, t3489, t3490, t3495, t3496;
994 double t3499, t3500, t3502, t3503, t3504, t3505, t3506, t3508;
995 double t3509, t3510, t3513, t3514, t3515, t3516, t3517, t3518;
996 double t3520, t3521, t3522, t3523, t3524, t3528, t3530, t3531;
997 double t3532, t3534, t3535, t3536, t3537, t3538, t3539, t3542;
998 double t3543, t3545, t3546, t3547, t3550, t3552, t3554, t3556;
999 double t3560, t3562, t3563, t3564, t3566, t3567, t3568, t3570;
1000 double t3571, t3574, t3576, t3578, t3580, t3583, t3585, t3587;
1001 double t3589, t3591, t3593, t3598, t3601, t3602, t3605, t3608;
1002 double t3609, t3611, t3612, t3615, t3616, t3619, t3622, t3623;
1003 double t3625, t3626, t3630, t3632, t3634, t3636, t3638, t3646;
1004 double t3661, t3665, t3666, t3668, t3669, t3670, t3671, t3675;
1005 double t3676, t3678, t3679, t3680, t3681, t3682, t3688, t3689;
1006 double t3692, t3693, t3694, t3695, t3696, t3698, t3700, t3701;
1007 double t3703, t3704, t3705, t3707, t3708, t3709, t3711, t3713;
1008 double t3714, t3716, t3717, t3718, t3720, t3721, t3722, t3724;
1009 double t3726, t3727, t3729, t3730, t3731, t3733, t3734, t3735;
1010 double t3737, t3739, t3741, t3743, t3745, t3750, t3755, t3760;
1011 double t3762, t3766, t3771, t3775, t3781, t3784, t3785, t3789;
1012 double t3790, t3792, t3794, t3795;
1013 #endif
1014
1015 #endif
1016
1017 #endif
1018
1019 #endif
1020
1021 #endif
1022
1023
1024
1025 t1 = M_CBRT3;
1026 t2 = 0.1e1 / M_PI;
1027 t3 = POW_1_3(t2);
1028 t4 = t1 * t3;
1029 t5 = M_CBRT4;
1030 t6 = t5 * t5;
1031 t7 = rho[0] + rho[1];
1032 t8 = POW_1_3(t7);
1033 t10 = t6 / t8;
1034 t11 = t4 * t10;
1035 t12 = t11 / 0.4e1;
1036 t13 = 0.1e1 <= t12;
1037 t14 = sqrt(t11);
1038 t17 = 0.1e1 + 0.52645000000000000000e0 * t14 + 0.83350000000000000000e-1 * t11;
1039 t20 = log(t12);
1040 t23 = t4 * t10 * t20;
1041 t27 = my_piecewise3(t13, -0.1423e0 / t17, 0.311e-1 * t20 - 0.48e-1 + 0.50000000000000000000e-3 * t23 - 0.29000000000000000000e-2 * t11);
1042 t30 = 0.1e1 + 0.69905000000000000000e0 * t14 + 0.65275000000000000000e-1 * t11;
1043 t37 = my_piecewise3(t13, -0.843e-1 / t30, 0.1555e-1 * t20 - 0.269e-1 + 0.17500000000000000000e-3 * t23 - 0.12000000000000000000e-2 * t11);
1044 t38 = t37 - t27;
1045 t39 = rho[0] - rho[1];
1046 t40 = 0.1e1 / t7;
1047 t41 = t39 * t40;
1048 t42 = 0.1e1 + t41;
1049 t43 = POW_1_3(t42);
1050 t44 = t43 * t42;
1051 t45 = 0.1e1 - t41;
1052 t46 = POW_1_3(t45);
1053 t47 = t46 * t45;
1054 t48 = t44 + t47 - 0.2e1;
1055 t50 = M_CBRT2;
1056 t53 = 0.1e1 / (0.2e1 * t50 - 0.2e1);
1057 t54 = t38 * t48 * t53;
1058 t56 = sigma[0] + 0.2e1 * sigma[1] + sigma[2];
1059 t57 = t7 * t7;
1060 t59 = 0.1e1 / t8 / t57;
1061 t60 = t56 * t59;
1062 t62 = t1 * t1;
1063 t63 = t3 * t3;
1064 t64 = t62 * t63;
1065 t65 = t8 * t8;
1066 t68 = t64 * t5 / t65;
1067 t70 = 0.2568e-2 + 0.58165000000000000000e-2 * t11 + 0.18472500000000000000e-5 * t68;
1068 t74 = 0.1e1 + 0.21807500000000000000e1 * t11 + 0.11800000000000000000e0 * t68 + 0.17639938117590219540e-1 * t40;
1069 t75 = 0.1e1 / t74;
1070 t77 = 0.1667e-2 + t70 * t75;
1071 t78 = 0.1e1 / t77;
1072 t79 = sqrt(t56);
1073 t80 = t78 * t79;
1074 t81 = pow(t7, 0.1e1 / 0.6e1);
1075 t83 = 0.1e1 / t81 / t7;
1076 t86 = exp(-0.81290825e-3 * t80 * t83);
1077 t87 = t60 * t86;
1078 t88 = t43 * t43;
1079 t89 = t88 * t42;
1080 t90 = t46 * t46;
1081 t91 = t90 * t45;
1082 t92 = t89 + t91;
1083 t93 = sqrt(t92);
1084 t94 = 0.1e1 / t93;
1085 t95 = t77 * t94;
1086 t96 = M_SQRT2;
1087 t97 = t95 * t96;
1088 t98 = t87 * t97;
1089 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
1090 zk[0] = t27 + t54 + t98;
1091
1092 #ifndef XC_DONT_COMPILE_VXC
1093
1094 if(order < 1) return;
1095
1096
1097 t99 = t17 * t17;
1098 t100 = 0.1e1 / t99;
1099 t102 = 0.1e1 / t14 * t1;
1100 t103 = t3 * t6;
1101 t105 = 0.1e1 / t8 / t7;
1102 t107 = t102 * t103 * t105;
1103 t109 = t6 * t105;
1104 t110 = t4 * t109;
1105 t112 = -0.87741666666666666667e-1 * t107 - 0.27783333333333333333e-1 * t110;
1106 t117 = t4 * t109 * t20;
1107 t121 = my_piecewise3(t13, 0.1423e0 * t100 * t112, -0.10366666666666666667e-1 * t40 - 0.16666666666666666667e-3 * t117 + 0.80000000000000000000e-3 * t110);
1108 t122 = t30 * t30;
1109 t123 = 0.1e1 / t122;
1110 t126 = -0.11650833333333333333e0 * t107 - 0.21758333333333333333e-1 * t110;
1111 t133 = my_piecewise3(t13, 0.843e-1 * t123 * t126, -0.51833333333333333333e-2 * t40 - 0.58333333333333333333e-4 * t117 + 0.34166666666666666667e-3 * t110);
1112 t134 = t133 - t121;
1113 t136 = t134 * t48 * t53;
1114 t137 = 0.1e1 / t57;
1115 t138 = t39 * t137;
1116 t139 = t40 - t138;
1117 t141 = -t139;
1118 t144 = 0.4e1 / 0.3e1 * t43 * t139 + 0.4e1 / 0.3e1 * t46 * t141;
1119 t146 = t38 * t144 * t53;
1120 t147 = t57 * t7;
1121 t149 = 0.1e1 / t8 / t147;
1122 t150 = t56 * t149;
1123 t151 = t150 * t86;
1124 t152 = t151 * t97;
1125 t153 = 0.7e1 / 0.3e1 * t152;
1126 t154 = t77 * t77;
1127 t155 = 0.1e1 / t154;
1128 t156 = t155 * t79;
1129 t161 = t64 * t5 / t65 / t7;
1130 t163 = -0.19388333333333333333e-2 * t110 - 0.12315000000000000000e-5 * t161;
1131 t165 = t74 * t74;
1132 t166 = 0.1e1 / t165;
1133 t167 = t70 * t166;
1134 t171 = -0.72691666666666666667e0 * t110 - 0.78666666666666666667e-1 * t161 - 0.17639938117590219540e-1 * t137;
1135 t173 = t163 * t75 - t167 * t171;
1136 t174 = t83 * t173;
1137 t178 = 0.1e1 / t81 / t57;
1138 t181 = 0.81290825e-3 * t156 * t174 + 0.94839295833333333334e-3 * t80 * t178;
1139 t182 = t60 * t181;
1140 t183 = t86 * t77;
1141 t184 = t94 * t96;
1142 t185 = t183 * t184;
1143 t186 = t182 * t185;
1144 t188 = t173 * t94 * t96;
1145 t189 = t87 * t188;
1146 t191 = 0.1e1 / t93 / t92;
1147 t192 = t77 * t191;
1148 t196 = 0.5e1 / 0.3e1 * t88 * t139 + 0.5e1 / 0.3e1 * t90 * t141;
1149 t197 = t96 * t196;
1150 t198 = t192 * t197;
1151 t199 = t87 * t198;
1152 t200 = t199 / 0.2e1;
1153 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1154 vrho[0] = t27 + t54 + t98 + t7 * (t121 + t136 + t146 - t153 + t186 + t189 - t200);
1155
1156 t203 = -t40 - t138;
1157 t205 = -t203;
1158 t208 = 0.4e1 / 0.3e1 * t43 * t203 + 0.4e1 / 0.3e1 * t46 * t205;
1159 t210 = t38 * t208 * t53;
1160 t214 = 0.5e1 / 0.3e1 * t88 * t203 + 0.5e1 / 0.3e1 * t90 * t205;
1161 t215 = t96 * t214;
1162 t216 = t192 * t215;
1163 t217 = t87 * t216;
1164 t218 = t217 / 0.2e1;
1165 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1166 vrho[1] = t27 + t54 + t98 + t7 * (t121 + t136 + t210 - t153 + t186 + t189 - t218);
1167
1168 t221 = t59 * t86;
1169 t222 = t221 * t97;
1170 t223 = sqrt(t7);
1171 t225 = 0.1e1 / t223 / t147;
1172 t226 = t79 * t225;
1173 t228 = t86 * t94 * t96;
1174 t229 = t226 * t228;
1175 t230 = 0.40645412500000000000e-3 * t229;
1176 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1177 vsigma[0] = t7 * (t222 - t230);
1178
1179 t232 = 0.2e1 * t222;
1180 t233 = 0.81290825e-3 * t229;
1181 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1182 vsigma[1] = t7 * (t232 - t233);
1183
1184 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1185 vsigma[2] = vsigma[0];
1186
1187 #ifndef XC_DONT_COMPILE_FXC
1188
1189 if(order < 2) return;
1190
1191
1192 t235 = 0.2e1 * t121;
1193 t236 = 0.2e1 * t136;
1194 t238 = 0.14e2 / 0.3e1 * t152;
1195 t239 = 0.2e1 * t186;
1196 t240 = 0.2e1 * t189;
1197 t242 = 0.1e1 / t99 / t17;
1198 t243 = t112 * t112;
1199 t248 = 0.1e1 / t14 / t11 * t62;
1200 t249 = t63 * t5;
1201 t251 = 0.1e1 / t65 / t57;
1202 t253 = t248 * t249 * t251;
1203 t256 = t102 * t103 * t59;
1204 t258 = t6 * t59;
1205 t259 = t4 * t258;
1206 t261 = -0.58494444444444444445e-1 * t253 + 0.11698888888888888889e0 * t256 + 0.37044444444444444444e-1 * t259;
1207 t267 = t4 * t258 * t20;
1208 t271 = my_piecewise3(t13, -0.2846e0 * t242 * t243 + 0.1423e0 * t100 * t261, 0.10366666666666666667e-1 * t137 + 0.22222222222222222223e-3 * t267 - 0.10111111111111111111e-2 * t259);
1209 t273 = 0.1e1 / t122 / t30;
1210 t274 = t126 * t126;
1211 t280 = -0.77672222222222222220e-1 * t253 + 0.15534444444444444444e0 * t256 + 0.29011111111111111111e-1 * t259;
1212 t288 = my_piecewise3(t13, -0.1686e0 * t273 * t274 + 0.843e-1 * t123 * t280, 0.51833333333333333333e-2 * t137 + 0.77777777777777777777e-4 * t267 - 0.43611111111111111112e-3 * t259);
1213 t289 = t288 - t271;
1214 t291 = t289 * t48 * t53;
1215 t293 = t134 * t144 * t53;
1216 t294 = 0.2e1 * t293;
1217 t295 = 0.1e1 / t88;
1218 t296 = t139 * t139;
1219 t299 = 0.1e1 / t147;
1220 t300 = t39 * t299;
1221 t302 = -0.2e1 * t137 + 0.2e1 * t300;
1222 t305 = 0.1e1 / t90;
1223 t306 = t141 * t141;
1224 t309 = -t302;
1225 t312 = 0.4e1 / 0.9e1 * t295 * t296 + 0.4e1 / 0.3e1 * t43 * t302 + 0.4e1 / 0.9e1 * t305 * t306 + 0.4e1 / 0.3e1 * t46 * t309;
1226 t314 = t38 * t312 * t53;
1227 t315 = t57 * t57;
1228 t317 = 0.1e1 / t8 / t315;
1229 t318 = t56 * t317;
1230 t319 = t318 * t86;
1231 t320 = t319 * t97;
1232 t321 = 0.70e2 / 0.9e1 * t320;
1233 t322 = t150 * t181;
1234 t323 = t322 * t185;
1235 t324 = 0.14e2 / 0.3e1 * t323;
1236 t325 = t151 * t188;
1237 t326 = 0.14e2 / 0.3e1 * t325;
1238 t327 = t151 * t198;
1239 t328 = 0.7e1 / 0.3e1 * t327;
1240 t330 = 0.1e1 / t154 / t77;
1241 t331 = t330 * t79;
1242 t332 = t173 * t173;
1243 t333 = t83 * t332;
1244 t336 = t178 * t173;
1245 t341 = t64 * t5 * t251;
1246 t343 = 0.25851111111111111111e-2 * t259 + 0.20525000000000000000e-5 * t341;
1247 t345 = t163 * t166;
1248 t349 = 0.1e1 / t165 / t74;
1249 t350 = t70 * t349;
1250 t351 = t171 * t171;
1251 t357 = 0.96922222222222222222e0 * t259 + 0.13111111111111111111e0 * t341 + 0.35279876235180439080e-1 * t299;
1252 t359 = -t167 * t357 - 0.2e1 * t345 * t171 + t343 * t75 + 0.2e1 * t350 * t351;
1253 t364 = 0.1e1 / t81 / t147;
1254 t367 = -0.162581650e-2 * t331 * t333 - 0.18967859166666666667e-2 * t156 * t336 + 0.81290825e-3 * t156 * t83 * t359 - 0.20548514097222222222e-2 * t80 * t364;
1255 t368 = t60 * t367;
1256 t369 = t368 * t185;
1257 t370 = t181 * t181;
1258 t371 = t60 * t370;
1259 t372 = t371 * t185;
1260 t373 = t86 * t173;
1261 t374 = t373 * t184;
1262 t375 = t182 * t374;
1263 t376 = 0.2e1 * t375;
1264 t377 = t181 * t86;
1265 t378 = t60 * t377;
1266 t379 = t378 * t198;
1267 t381 = t359 * t94 * t96;
1268 t382 = t87 * t381;
1269 t383 = t173 * t191;
1270 t384 = t383 * t197;
1271 t385 = t87 * t384;
1272 t386 = t92 * t92;
1273 t388 = 0.1e1 / t93 / t386;
1274 t389 = t77 * t388;
1275 t390 = t196 * t196;
1276 t391 = t96 * t390;
1277 t392 = t389 * t391;
1278 t393 = t87 * t392;
1279 t394 = 0.3e1 / 0.4e1 * t393;
1280 t395 = 0.1e1 / t43;
1281 t400 = 0.1e1 / t46;
1282 t405 = 0.10e2 / 0.9e1 * t395 * t296 + 0.5e1 / 0.3e1 * t88 * t302 + 0.10e2 / 0.9e1 * t400 * t306 + 0.5e1 / 0.3e1 * t90 * t309;
1283 t406 = t96 * t405;
1284 t407 = t192 * t406;
1285 t408 = t87 * t407;
1286 t409 = t408 / 0.2e1;
1287 t410 = t271 + t291 + t294 + t314 + t321 - t324 - t326 + t328 + t369 + t372 + t376 - t379 + t382 - t385 + t394 - t409;
1288 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1289 v2rho2[0] = t7 * t410 + 0.2e1 * t146 - t199 + t235 + t236 - t238 + t239 + t240;
1290
1291 t413 = t134 * t208 * t53;
1292 t414 = t295 * t203;
1293 t417 = t43 * t39;
1294 t420 = t305 * t205;
1295 t423 = t46 * t39;
1296 t426 = 0.4e1 / 0.9e1 * t414 * t139 + 0.8e1 / 0.3e1 * t417 * t299 + 0.4e1 / 0.9e1 * t420 * t141 - 0.8e1 / 0.3e1 * t423 * t299;
1297 t428 = t38 * t426 * t53;
1298 t432 = t151 * t216;
1299 t434 = t378 * t216;
1300 t436 = t383 * t215;
1301 t437 = t87 * t436;
1302 t439 = t60 * t183;
1303 t440 = t388 * t96;
1304 t441 = t214 * t196;
1305 t442 = t440 * t441;
1306 t443 = t439 * t442;
1307 t445 = t395 * t203;
1308 t448 = t88 * t39;
1309 t451 = t400 * t205;
1310 t454 = t90 * t39;
1311 t457 = 0.10e2 / 0.9e1 * t445 * t139 + 0.10e2 / 0.3e1 * t448 * t299 + 0.10e2 / 0.9e1 * t451 * t141 - 0.10e2 / 0.3e1 * t454 * t299;
1312 t458 = t96 * t457;
1313 t459 = t192 * t458;
1314 t460 = t87 * t459;
1315 t462 = t271 + t291 + t293 + t413 + t428 + t321 - t324 - t326 + 0.7e1 / 0.6e1 * t327 + t369 + t372 + t376 - t379 / 0.2e1 + t382 - t385 / 0.2e1 + 0.7e1 / 0.6e1 * t432 - t434 / 0.2e1 - t437 / 0.2e1 + 0.3e1 / 0.4e1 * t443 - t460 / 0.2e1;
1316 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1317 v2rho2[1] = t7 * t462 + t146 - t200 + t210 - t218 + t235 + t236 - t238 + t239 + t240;
1318
1319 t465 = 0.2e1 * t413;
1320 t466 = t203 * t203;
1321 t470 = 0.2e1 * t137 + 0.2e1 * t300;
1322 t473 = t205 * t205;
1323 t476 = -t470;
1324 t479 = 0.4e1 / 0.9e1 * t295 * t466 + 0.4e1 / 0.3e1 * t43 * t470 + 0.4e1 / 0.9e1 * t305 * t473 + 0.4e1 / 0.3e1 * t46 * t476;
1325 t481 = t38 * t479 * t53;
1326 t482 = 0.7e1 / 0.3e1 * t432;
1327 t483 = t214 * t214;
1328 t484 = t96 * t483;
1329 t485 = t389 * t484;
1330 t486 = t87 * t485;
1331 t487 = 0.3e1 / 0.4e1 * t486;
1332 t496 = 0.10e2 / 0.9e1 * t395 * t466 + 0.5e1 / 0.3e1 * t88 * t470 + 0.10e2 / 0.9e1 * t400 * t473 + 0.5e1 / 0.3e1 * t90 * t476;
1333 t497 = t96 * t496;
1334 t498 = t192 * t497;
1335 t499 = t87 * t498;
1336 t500 = t499 / 0.2e1;
1337 t501 = t271 + t291 + t465 + t481 + t321 - t324 - t326 + t482 + t369 + t372 + t376 - t434 + t382 - t437 + t487 - t500;
1338 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1339 v2rho2[2] = t7 * t501 + 0.2e1 * t210 - t217 + t235 + t236 - t238 + t239 + t240;
1340
1341 t503 = t149 * t86;
1342 t504 = t503 * t97;
1343 t505 = 0.7e1 / 0.3e1 * t504;
1344 t506 = t59 * t181;
1345 t507 = t506 * t86;
1346 t508 = t507 * t97;
1347 t509 = t221 * t188;
1348 t510 = t221 * t77;
1349 t511 = t191 * t96;
1350 t512 = t511 * t196;
1351 t513 = t510 * t512;
1352 t514 = t513 / 0.2e1;
1353 t516 = 0.1e1 / t223 / t315;
1354 t517 = t79 * t516;
1355 t518 = t517 * t228;
1356 t519 = 0.14225894375000000000e-2 * t518;
1357 t520 = t226 * t181;
1358 t521 = t520 * t228;
1359 t522 = 0.40645412500000000000e-3 * t521;
1360 t523 = t226 * t86;
1361 t524 = t523 * t512;
1362 t525 = 0.20322706250000000000e-3 * t524;
1363 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1364 v2rhosigma[0] = t222 - t230 + t7 * (-t505 + t508 + t509 - t514 + t519 - t522 + t525);
1365
1366 t528 = 0.14e2 / 0.3e1 * t504;
1367 t529 = 0.2e1 * t508;
1368 t530 = 0.2e1 * t509;
1369 t531 = 0.28451788750000000000e-2 * t518;
1370 t532 = 0.81290825e-3 * t521;
1371 t533 = 0.40645412500000000000e-3 * t524;
1372 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1373 v2rhosigma[1] = t232 - t233 + t7 * (-t528 + t529 + t530 - t513 + t531 - t532 + t533);
1374
1375 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1376 v2rhosigma[2] = v2rhosigma[0];
1377
1378 t536 = t511 * t214;
1379 t537 = t510 * t536;
1380 t538 = t537 / 0.2e1;
1381 t539 = t523 * t536;
1382 t540 = 0.20322706250000000000e-3 * t539;
1383 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1384 v2rhosigma[3] = t222 - t230 + t7 * (-t505 + t508 + t509 - t538 + t519 - t522 + t540);
1385
1386 t543 = 0.40645412500000000000e-3 * t539;
1387 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1388 v2rhosigma[4] = t232 - t233 + t7 * (-t528 + t529 + t530 - t537 + t531 - t532 + t543);
1389
1390 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1391 v2rhosigma[5] = v2rhosigma[3];
1392
1393 t546 = 0.1e1 / t79;
1394 t547 = t225 * t546;
1395 t548 = t547 * t228;
1396 t549 = 0.60968118750000000000e-3 * t548;
1397 t551 = 0.1e1 / t65 / t315;
1398 t552 = t551 * t78;
1399 t553 = t552 * t228;
1400 t554 = 0.16520495572951562500e-6 * t553;
1401 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1402 v2sigma2[0] = t7 * (-t549 + t554);
1403
1404 t556 = 0.12193623750000000000e-2 * t548;
1405 t557 = 0.33040991145903125000e-6 * t553;
1406 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1407 v2sigma2[1] = t7 * (-t556 + t557);
1408
1409 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1410 v2sigma2[2] = v2sigma2[0];
1411
1412 t559 = 0.243872475e-2 * t548;
1413 t560 = 0.6608198229180625e-6 * t553;
1414 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1415 v2sigma2[3] = t7 * (-t559 + t560);
1416
1417 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1418 v2sigma2[4] = v2sigma2[1];
1419
1420 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1421 v2sigma2[5] = v2sigma2[2];
1422
1423 #ifndef XC_DONT_COMPILE_KXC
1424
1425 if(order < 3) return;
1426
1427
1428 t562 = t315 * t7;
1429 t564 = 0.1e1 / t8 / t562;
1430 t565 = t56 * t564;
1431 t566 = t565 * t86;
1432 t567 = t566 * t97;
1433 t568 = 0.910e3 / 0.27e2 * t567;
1434 t569 = t6 * t149;
1435 t570 = t4 * t569;
1436 t573 = 0.1e1 / t65 / t147;
1437 t575 = t64 * t5 * t573;
1438 t577 = -0.60319259259259259259e-2 * t570 - 0.54733333333333333333e-5 * t575;
1439 t579 = t343 * t166;
1440 t582 = t163 * t349;
1441 t587 = t165 * t165;
1442 t588 = 0.1e1 / t587;
1443 t589 = t70 * t588;
1444 t590 = t351 * t171;
1445 t593 = t171 * t357;
1446 t598 = 0.1e1 / t315;
1447 t600 = -0.22615185185185185185e1 * t570 - 0.34962962962962962963e0 * t575 - 0.10583962870554131724e0 * t598;
1448 t602 = -t167 * t600 - 0.3e1 * t579 * t171 - 0.3e1 * t345 * t357 + 0.6e1 * t350 * t593 + 0.6e1 * t582 * t351 + t577 * t75 - 0.6e1 * t589 * t590;
1449 t604 = t602 * t94 * t96;
1450 t605 = t87 * t604;
1451 t606 = t151 * t381;
1452 t607 = 0.7e1 * t606;
1453 t608 = t319 * t188;
1454 t609 = 0.70e2 / 0.3e1 * t608;
1455 t610 = t383 * t406;
1456 t611 = t87 * t610;
1457 t613 = 0.1e1 / t44;
1458 t614 = t296 * t139;
1459 t617 = t395 * t139;
1460 t620 = t39 * t598;
1461 t622 = 0.6e1 * t299 - 0.6e1 * t620;
1462 t625 = 0.1e1 / t47;
1463 t626 = t306 * t141;
1464 t629 = t400 * t141;
1465 t632 = -t622;
1466 t635 = -0.10e2 / 0.27e2 * t613 * t614 + 0.10e2 / 0.3e1 * t617 * t302 + 0.5e1 / 0.3e1 * t88 * t622 - 0.10e2 / 0.27e2 * t625 * t626 + 0.10e2 / 0.3e1 * t629 * t309 + 0.5e1 / 0.3e1 * t90 * t632;
1467 t636 = t96 * t635;
1468 t637 = t192 * t636;
1469 t638 = t87 * t637;
1470 t639 = t638 / 0.2e1;
1471 t640 = t150 * t370;
1472 t641 = t640 * t185;
1473 t642 = 0.7e1 * t641;
1474 t643 = t151 * t392;
1475 t645 = t370 * t181;
1476 t646 = t60 * t645;
1477 t647 = t646 * t185;
1478 t648 = t173 * t388;
1479 t649 = t648 * t391;
1480 t650 = t87 * t649;
1481 t654 = 0.1e1 / t93 / t386 / t92;
1482 t655 = t77 * t654;
1483 t656 = t390 * t196;
1484 t657 = t96 * t656;
1485 t658 = t655 * t657;
1486 t659 = t87 * t658;
1487 t660 = 0.15e2 / 0.8e1 * t659;
1488 t661 = t319 * t198;
1489 t662 = 0.35e2 / 0.3e1 * t661;
1490 t663 = t318 * t181;
1491 t664 = t663 * t185;
1492 t665 = 0.70e2 / 0.3e1 * t664;
1493 t666 = t150 * t367;
1494 t667 = t666 * t185;
1495 t668 = 0.7e1 * t667;
1496 t669 = t322 * t374;
1497 t670 = 0.14e2 * t669;
1498 t671 = t151 * t384;
1499 t672 = 0.7e1 * t671;
1500 t673 = t151 * t407;
1501 t675 = -t568 + t605 - t607 + t609 - 0.3e1 / 0.2e1 * t611 - t639 - t642 - 0.21e2 / 0.4e1 * t643 + t647 + 0.9e1 / 0.4e1 * t650 - t660 - t662 + t665 - t668 - t670 + t672 + 0.7e1 / 0.2e1 * t673;
1502 t676 = t154 * t154;
1503 t678 = 0.1e1 / t676 * t79;
1504 t679 = t332 * t173;
1505 t699 = 0.1e1 / t81 / t315;
1506 t702 = 0.487744950e-2 * t678 * t83 * t679 + 0.56903577500000000001e-2 * t331 * t178 * t332 - 0.487744950e-2 * t331 * t174 * t359 + 0.61645542291666666667e-2 * t156 * t364 * t173 - 0.28451788750000000000e-2 * t156 * t178 * t359 + 0.81290825e-3 * t156 * t83 * t602 + 0.65070294641203703703e-2 * t80 * t699;
1507 t703 = t60 * t702;
1508 t704 = t703 * t185;
1509 t705 = t368 * t374;
1510 t706 = 0.3e1 * t705;
1511 t707 = t371 * t374;
1512 t708 = 0.3e1 * t707;
1513 t709 = t86 * t359;
1514 t710 = t709 * t184;
1515 t711 = t182 * t710;
1516 t712 = 0.3e1 * t711;
1517 t713 = t359 * t191;
1518 t714 = t713 * t197;
1519 t715 = t87 * t714;
1520 t716 = 0.3e1 / 0.2e1 * t715;
1521 t717 = t99 * t99;
1522 t718 = 0.1e1 / t717;
1523 t722 = t242 * t112;
1524 t728 = 0.1e1 / t14 / t68 * t2 / 0.4e1;
1525 t729 = t728 * t598;
1526 t732 = t248 * t249 * t573;
1527 t735 = t102 * t103 * t149;
1528 t738 = -0.35096666666666666667e0 * t729 + 0.23397777777777777778e0 * t732 - 0.27297407407407407408e0 * t735 - 0.86437037037037037036e-1 * t570;
1529 t744 = t4 * t569 * t20;
1530 t748 = my_piecewise3(t13, 0.8538e0 * t718 * t243 * t112 - 0.8538e0 * t722 * t261 + 0.1423e0 * t100 * t738, -0.20733333333333333334e-1 * t299 - 0.51851851851851851854e-3 * t744 + 0.22851851851851851851e-2 * t570);
1531 t749 = t367 * t86;
1532 t750 = t60 * t749;
1533 t751 = t750 * t198;
1534 t752 = 0.3e1 / 0.2e1 * t751;
1535 t754 = t60 * t367 * t181;
1536 t755 = t754 * t185;
1537 t756 = 0.3e1 * t755;
1538 t757 = t370 * t86;
1539 t758 = t60 * t757;
1540 t759 = t758 * t198;
1541 t760 = 0.3e1 / 0.2e1 * t759;
1542 t761 = t378 * t384;
1543 t762 = 0.3e1 * t761;
1544 t763 = t378 * t407;
1545 t766 = t440 * t196 * t405;
1546 t767 = t439 * t766;
1547 t768 = 0.9e1 / 0.4e1 * t767;
1548 t769 = t378 * t392;
1549 t771 = t150 * t377;
1550 t772 = t771 * t198;
1551 t773 = 0.7e1 * t772;
1552 t774 = t122 * t122;
1553 t775 = 0.1e1 / t774;
1554 t779 = t273 * t126;
1555 t786 = -0.46603333333333333332e0 * t729 + 0.31068888888888888888e0 * t732 - 0.36247037037037037036e0 * t735 - 0.67692592592592592592e-1 * t570;
1556 t794 = my_piecewise3(t13, 0.5058e0 * t775 * t274 * t126 - 0.5058e0 * t779 * t280 + 0.843e-1 * t123 * t786, -0.10366666666666666667e-1 * t299 - 0.18148148148148148148e-3 * t744 + 0.99166666666666666667e-3 * t570);
1557 t795 = t794 - t748;
1558 t797 = t795 * t48 * t53;
1559 t799 = t289 * t144 * t53;
1560 t800 = 0.3e1 * t799;
1561 t802 = t134 * t312 * t53;
1562 t804 = 0.1e1 / t89;
1563 t807 = t295 * t139;
1564 t812 = 0.1e1 / t91;
1565 t815 = t305 * t141;
1566 t820 = -0.8e1 / 0.27e2 * t804 * t614 + 0.4e1 / 0.3e1 * t807 * t302 + 0.4e1 / 0.3e1 * t43 * t622 - 0.8e1 / 0.27e2 * t812 * t626 + 0.4e1 / 0.3e1 * t815 * t309 + 0.4e1 / 0.3e1 * t46 * t632;
1567 t822 = t38 * t820 * t53;
1568 t823 = t704 + t706 + t708 + t712 - t716 + t748 - t752 + t756 - t760 - t762 - 0.3e1 / 0.2e1 * t763 + t768 + 0.9e1 / 0.4e1 * t769 + t773 + t797 + t800 + 0.3e1 * t802 + t822;
1569 t827 = 0.3e1 * t372;
1570 t829 = 0.3e1 * t291;
1571 t832 = 0.14e2 * t325;
1572 t833 = 0.3e1 * t382;
1573 t834 = 0.70e2 / 0.3e1 * t320;
1574 t836 = 0.14e2 * t323;
1575 t837 = 0.3e1 * t369;
1576 t838 = 0.6e1 * t375;
1577 t841 = 0.3e1 * t271;
1578 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1579 v3rho3[0] = t7 * (t675 + t823) - 0.3e1 * t379 + t827 + 0.9e1 / 0.4e1 * t393 + t829 + 0.6e1 * t293 + 0.3e1 * t314 - t832 + t833 + t834 + 0.7e1 * t327 - t836 + t837 + t838 - 0.3e1 * t385 - 0.3e1 / 0.2e1 * t408 + t841;
1580
1581 t842 = t215 * t196;
1582 t843 = t389 * t842;
1583 t844 = t378 * t843;
1584 t845 = 0.3e1 / 0.2e1 * t844;
1585 t846 = t151 * t436;
1586 t848 = t151 * t459;
1587 t849 = 0.7e1 / 0.3e1 * t848;
1588 t850 = t713 * t215;
1589 t851 = t87 * t850;
1590 t853 = t383 * t458;
1591 t854 = t87 * t853;
1592 t855 = t613 * t203;
1593 t858 = t395 * t39;
1594 t859 = t299 * t139;
1595 t868 = t625 * t205;
1596 t871 = t400 * t39;
1597 t872 = t299 * t141;
1598 t881 = -0.10e2 / 0.27e2 * t855 * t296 + 0.40e2 / 0.9e1 * t858 * t859 + 0.10e2 / 0.9e1 * t445 * t302 + 0.10e2 / 0.3e1 * t88 * t299 - 0.10e2 * t448 * t598 - 0.10e2 / 0.27e2 * t868 * t306 - 0.40e2 / 0.9e1 * t871 * t872 + 0.10e2 / 0.9e1 * t451 * t309 - 0.10e2 / 0.3e1 * t90 * t299 + 0.10e2 * t454 * t598;
1599 t882 = t96 * t881;
1600 t883 = t192 * t882;
1601 t884 = t87 * t883;
1602 t886 = t319 * t216;
1603 t889 = t845 - t568 + t605 - t607 + t609 + 0.7e1 / 0.3e1 * t846 + t849 - t851 / 0.2e1 - t854 - t884 / 0.2e1 - 0.35e2 / 0.9e1 * t886 - t611 / 0.2e1;
1604 t895 = -t642 - 0.7e1 / 0.4e1 * t643 + t647 + 0.3e1 / 0.4e1 * t650 - 0.70e2 / 0.9e1 * t661 + t665 - t668 - t670 + 0.14e2 / 0.3e1 * t671 + 0.7e1 / 0.6e1 * t673 + t704 + t706 + t708;
1605 t897 = t758 * t216;
1606 t899 = t654 * t96;
1607 t901 = t899 * t214 * t390;
1608 t902 = t439 * t901;
1609 t904 = t150 * t183;
1610 t905 = t904 * t442;
1611 t906 = 0.7e1 / 0.2e1 * t905;
1612 t907 = t771 * t216;
1613 t909 = t750 * t216;
1614 t911 = t378 * t436;
1615 t912 = t378 * t459;
1616 t913 = t60 * t373;
1617 t914 = t913 * t442;
1618 t915 = 0.3e1 / 0.2e1 * t914;
1619 t917 = t440 * t457 * t196;
1620 t918 = t439 * t917;
1621 t920 = t214 * t405;
1622 t921 = t440 * t920;
1623 t922 = t439 * t921;
1624 t924 = t712 - t715 + t748 - t897 / 0.2e1 - 0.15e2 / 0.8e1 * t902 - t906 + 0.7e1 / 0.3e1 * t907 - t909 / 0.2e1 - t911 - t912 + t915 + 0.3e1 / 0.2e1 * t918 + 0.3e1 / 0.4e1 * t922;
1625 t931 = t289 * t208 * t53;
1626 t933 = t134 * t426 * t53;
1627 t934 = 0.2e1 * t933;
1628 t935 = t804 * t203;
1629 t938 = t295 * t39;
1630 t947 = t812 * t205;
1631 t950 = t305 * t39;
1632 t959 = -0.8e1 / 0.27e2 * t935 * t296 + 0.16e2 / 0.9e1 * t938 * t859 + 0.4e1 / 0.9e1 * t414 * t302 + 0.8e1 / 0.3e1 * t43 * t299 - 0.8e1 * t417 * t598 - 0.8e1 / 0.27e2 * t947 * t306 - 0.16e2 / 0.9e1 * t950 * t872 + 0.4e1 / 0.9e1 * t420 * t309 - 0.8e1 / 0.3e1 * t46 * t299 + 0.8e1 * t423 * t598;
1633 t961 = t38 * t959 * t53;
1634 t962 = -t751 + t756 - t759 - 0.2e1 * t761 - t763 / 0.2e1 + 0.3e1 / 0.4e1 * t769 + 0.14e2 / 0.3e1 * t772 + t797 + 0.2e1 * t799 + t802 + t931 + t934 + t961;
1635 t967 = 0.3e1 / 0.2e1 * t443;
1636 t969 = 0.2e1 * t428;
1637 t970 = t7 * (t889 + t895 + t924 + t962) - 0.2e1 * t379 + t827 + t394 + t482 + t967 - t434 + t829 + 0.4e1 * t293 + t314 + t465 + t969;
1638 t973 = -t832 + t833 + t834 + 0.14e2 / 0.3e1 * t327 - t836 + t837 + t838 - 0.2e1 * t385 - t409 - t437 - t460 + t841;
1639 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1640 v3rho3[1] = t970 + t973;
1641
1642 t974 = t648 * t484;
1643 t975 = t87 * t974;
1644 t977 = t383 * t497;
1645 t978 = t87 * t977;
1646 t980 = t613 * t466;
1647 t985 = t395 * t470;
1648 t990 = -0.2e1 * t299 - 0.6e1 * t620;
1649 t993 = t625 * t473;
1650 t998 = t400 * t476;
1651 t1001 = -t990;
1652 t1004 = -0.10e2 / 0.27e2 * t980 * t139 + 0.40e2 / 0.9e1 * t445 * t300 + 0.10e2 / 0.9e1 * t985 * t139 + 0.5e1 / 0.3e1 * t88 * t990 - 0.10e2 / 0.27e2 * t993 * t141 - 0.40e2 / 0.9e1 * t451 * t300 + 0.10e2 / 0.9e1 * t998 * t141 + 0.5e1 / 0.3e1 * t90 * t1001;
1653 t1005 = t96 * t1004;
1654 t1006 = t192 * t1005;
1655 t1007 = t87 * t1006;
1656 t1009 = t151 * t485;
1657 t1011 = t151 * t498;
1658 t1014 = t845 - t568 + t605 - t607 + t609 + 0.3e1 / 0.4e1 * t975 - t978 / 0.2e1 - t1007 / 0.2e1 - 0.7e1 / 0.4e1 * t1009 + 0.7e1 / 0.6e1 * t1011 + 0.14e2 / 0.3e1 * t846 + t849;
1659 t1018 = -t851 - t854 - 0.70e2 / 0.9e1 * t886 - t642 + t647 - 0.35e2 / 0.9e1 * t661 + t665 - t668 - t670 + 0.7e1 / 0.3e1 * t671 + t704 + t706 + t708;
1660 t1022 = t899 * t483 * t196;
1661 t1023 = t439 * t1022;
1662 t1026 = t440 * t214 * t457;
1663 t1027 = t439 * t1026;
1664 t1029 = t378 * t485;
1665 t1032 = t440 * t496 * t196;
1666 t1033 = t439 * t1032;
1667 t1035 = t378 * t498;
1668 t1039 = t712 - t715 / 0.2e1 + t748 - 0.15e2 / 0.8e1 * t1023 + 0.3e1 / 0.2e1 * t1027 + 0.3e1 / 0.4e1 * t1029 + 0.3e1 / 0.4e1 * t1033 - t1035 / 0.2e1 - t897 - t906 + 0.14e2 / 0.3e1 * t907 - t909 - 0.2e1 * t911;
1669 t1045 = t134 * t479 * t53;
1670 t1046 = t804 * t466;
1671 t1051 = t295 * t470;
1672 t1056 = t812 * t473;
1673 t1061 = t305 * t476;
1674 t1066 = -0.8e1 / 0.27e2 * t1046 * t139 + 0.16e2 / 0.9e1 * t414 * t300 + 0.4e1 / 0.9e1 * t1051 * t139 + 0.4e1 / 0.3e1 * t43 * t990 - 0.8e1 / 0.27e2 * t1056 * t141 - 0.16e2 / 0.9e1 * t420 * t300 + 0.4e1 / 0.9e1 * t1061 * t141 + 0.4e1 / 0.3e1 * t46 * t1001;
1675 t1068 = t38 * t1066 * t53;
1676 t1069 = -t912 + t915 - t751 / 0.2e1 + t756 - t759 / 0.2e1 - t761 + 0.7e1 / 0.3e1 * t772 + t797 + t799 + 0.2e1 * t931 + t934 + t1045 + t1068;
1677 t1076 = t7 * (t1014 + t1018 + t1039 + t1069) + t481 - t379 + t827 + 0.14e2 / 0.3e1 * t432 + t967 - 0.2e1 * t434 + t487 + t829 + t294 + 0.4e1 * t413 + t969;
1678 t1078 = -t832 + t833 + t834 + t328 - t836 + t837 + t838 - t385 - 0.2e1 * t437 - t460 - t500 + t841;
1679 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1680 v3rho3[2] = t1076 + t1078;
1681
1682 t1079 = t466 * t203;
1683 t1085 = -0.6e1 * t299 - 0.6e1 * t620;
1684 t1088 = t473 * t205;
1685 t1093 = -t1085;
1686 t1096 = -0.10e2 / 0.27e2 * t613 * t1079 + 0.10e2 / 0.3e1 * t445 * t470 + 0.5e1 / 0.3e1 * t88 * t1085 - 0.10e2 / 0.27e2 * t625 * t1088 + 0.10e2 / 0.3e1 * t451 * t476 + 0.5e1 / 0.3e1 * t90 * t1093;
1687 t1097 = t96 * t1096;
1688 t1098 = t192 * t1097;
1689 t1099 = t87 * t1098;
1690 t1100 = t1099 / 0.2e1;
1691 t1101 = t483 * t214;
1692 t1102 = t96 * t1101;
1693 t1103 = t655 * t1102;
1694 t1104 = t87 * t1103;
1695 t1105 = 0.15e2 / 0.8e1 * t1104;
1696 t1110 = 0.7e1 * t846;
1697 t1111 = 0.3e1 / 0.2e1 * t851;
1698 t1112 = 0.35e2 / 0.3e1 * t886;
1699 t1113 = -t568 + t605 - t607 + t609 - t1100 - t1105 + 0.9e1 / 0.4e1 * t975 - 0.3e1 / 0.2e1 * t978 - 0.21e2 / 0.4e1 * t1009 + 0.7e1 / 0.2e1 * t1011 + t1110 - t1111 - t1112 - t642 + t647 + t665 - t668;
1700 t1114 = t214 * t496;
1701 t1115 = t440 * t1114;
1702 t1116 = t439 * t1115;
1703 t1117 = 0.9e1 / 0.4e1 * t1116;
1704 t1120 = 0.3e1 / 0.2e1 * t897;
1705 t1121 = 0.7e1 * t907;
1706 t1122 = 0.3e1 / 0.2e1 * t909;
1707 t1123 = 0.3e1 * t911;
1708 t1124 = 0.3e1 * t931;
1709 t1138 = -0.8e1 / 0.27e2 * t804 * t1079 + 0.4e1 / 0.3e1 * t414 * t470 + 0.4e1 / 0.3e1 * t43 * t1085 - 0.8e1 / 0.27e2 * t812 * t1088 + 0.4e1 / 0.3e1 * t420 * t476 + 0.4e1 / 0.3e1 * t46 * t1093;
1710 t1140 = t38 * t1138 * t53;
1711 t1141 = -t670 + t704 + t706 + t708 + t712 + t748 + t1117 + 0.9e1 / 0.4e1 * t1029 - 0.3e1 / 0.2e1 * t1035 - t1120 + t1121 - t1122 - t1123 + t756 + t797 + t1124 + 0.3e1 * t1045 + t1140;
1712 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1713 v3rho3[3] = t7 * (t1113 + t1141) + 0.3e1 * t481 + t827 + 0.7e1 * t432 - 0.3e1 * t434 + 0.9e1 / 0.4e1 * t486 + t829 + 0.6e1 * t413 - t832 + t833 + t834 - t836 + t837 + t838 - 0.3e1 * t437 - 0.3e1 / 0.2e1 * t499 + t841;
1714
1715 t1151 = 0.81290825000000000000e-3 * t521;
1716 t1152 = t511 * t405;
1717 t1153 = t523 * t1152;
1718 t1154 = 0.20322706250000000000e-3 * t1153;
1719 t1155 = t221 * t381;
1720 t1156 = t503 * t188;
1721 t1157 = 0.14e2 / 0.3e1 * t1156;
1722 t1158 = t517 * t86;
1723 t1159 = t1158 * t512;
1724 t1160 = 0.14225894375000000000e-2 * t1159;
1725 t1161 = t517 * t181;
1726 t1162 = t1161 * t228;
1727 t1163 = 0.28451788750000000000e-2 * t1162;
1728 t1164 = t226 * t367;
1729 t1165 = t1164 * t228;
1730 t1166 = 0.40645412500000000000e-3 * t1165;
1731 t1167 = t221 * t173;
1732 t1168 = t1167 * t512;
1733 t1169 = t510 * t1152;
1734 t1170 = t1169 / 0.2e1;
1735 t1172 = 0.1e1 / t223 / t562;
1736 t1173 = t79 * t1172;
1737 t1174 = t1173 * t228;
1738 t1175 = 0.64016524687500000000e-2 * t1174;
1739 t1176 = t503 * t77;
1740 t1177 = t1176 * t512;
1741 t1178 = 0.7e1 / 0.3e1 * t1177;
1742 t1180 = t149 * t181 * t86;
1743 t1181 = t1180 * t97;
1744 t1182 = 0.14e2 / 0.3e1 * t1181;
1745 t1183 = t59 * t367;
1746 t1184 = t1183 * t86;
1747 t1185 = t1184 * t97;
1748 t1186 = t507 * t188;
1749 t1187 = 0.2e1 * t1186;
1750 t1188 = t317 * t86;
1751 t1189 = t1188 * t97;
1752 t1190 = 0.70e2 / 0.9e1 * t1189;
1753 t1191 = t507 * t198;
1754 t1193 = t59 * t370 * t86;
1755 t1194 = t1193 * t97;
1756 t1195 = t440 * t390;
1757 t1196 = t510 * t1195;
1758 t1197 = 0.3e1 / 0.4e1 * t1196;
1759 t1198 = t86 * t191;
1760 t1199 = t1198 * t197;
1761 t1200 = t520 * t1199;
1762 t1201 = 0.40645412500000000000e-3 * t1200;
1763 t1202 = t226 * t370;
1764 t1203 = t1202 * t228;
1765 t1204 = 0.40645412500000000000e-3 * t1203;
1766 t1205 = t523 * t1195;
1767 t1206 = 0.30484059375000000000e-3 * t1205;
1768 t1207 = t1154 + t1155 - t1157 - t1160 + t1163 - t1166 - t1168 - t1170 - t1175 + t1178 - t1182 + t1185 + t1187 + t1190 - t1191 + t1194 + t1197 + t1201 - t1204 - t1206;
1769 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1770 v3rho2sigma[0] = t7 * t1207 - t1151 - t513 - t528 + t529 + t530 + t531 + t533;
1771
1772 t1209 = 0.28e2 / 0.3e1 * t504;
1773 t1210 = 0.4e1 * t508;
1774 t1211 = 0.4e1 * t509;
1775 t1213 = 0.56903577500000000000e-2 * t518;
1776 t1214 = 0.162581650e-2 * t521;
1777 t1216 = 0.40645412500000000000e-3 * t1153;
1778 t1217 = 0.2e1 * t1155;
1779 t1218 = 0.28e2 / 0.3e1 * t1156;
1780 t1219 = 0.28451788750000000000e-2 * t1159;
1781 t1220 = 0.56903577500000000000e-2 * t1162;
1782 t1221 = 0.81290825e-3 * t1165;
1783 t1222 = 0.2e1 * t1168;
1784 t1223 = 0.12803304937500000000e-1 * t1174;
1785 t1224 = 0.14e2 / 0.3e1 * t1177;
1786 t1225 = 0.28e2 / 0.3e1 * t1181;
1787 t1226 = 0.2e1 * t1185;
1788 t1227 = 0.4e1 * t1186;
1789 t1228 = 0.140e3 / 0.9e1 * t1189;
1790 t1229 = 0.2e1 * t1191;
1791 t1230 = 0.2e1 * t1194;
1792 t1231 = 0.3e1 / 0.2e1 * t1196;
1793 t1232 = 0.81290825000000000000e-3 * t1200;
1794 t1233 = 0.81290825e-3 * t1203;
1795 t1234 = 0.60968118750000000000e-3 * t1205;
1796 t1235 = t1216 + t1217 - t1218 - t1219 + t1220 - t1221 - t1222 - t1169 - t1223 + t1224 - t1225 + t1226 + t1227 + t1228 - t1229 + t1230 + t1231 + t1232 - t1233 - t1234;
1797 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1798 v3rho2sigma[1] = -t1209 + t1210 + t1211 - 0.2e1 * t513 + t1213 - t1214 + 0.81290825000000000000e-3 * t524 + t7 * t1235;
1799
1800 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1801 v3rho2sigma[2] = v3rho2sigma[0];
1802
1803 t1237 = t1176 * t536;
1804 t1239 = t1158 * t536;
1805 t1241 = t511 * t457;
1806 t1242 = t523 * t1241;
1807 t1244 = t1167 * t536;
1808 t1246 = t510 * t1241;
1809 t1251 = -t1175 + 0.7e1 / 0.6e1 * t1237 - 0.71129471875000000000e-3 * t1239 + 0.20322706250000000000e-3 * t1242 - t1244 / 0.2e1 - t1246 / 0.2e1 - 0.71129471875000000000e-3 * t1159 + t1163 - t1166 - t1168 / 0.2e1 + 0.7e1 / 0.6e1 * t1177 - t1182;
1810 t1252 = t1198 * t215;
1811 t1253 = t520 * t1252;
1812 t1255 = t510 * t442;
1813 t1257 = t507 * t216;
1814 t1261 = t523 * t442;
1815 t1263 = t1185 + t1187 + t1194 - t1204 + t1155 - t1157 + 0.20322706250000000000e-3 * t1253 + 0.3e1 / 0.4e1 * t1255 - t1257 / 0.2e1 - t1191 / 0.2e1 + 0.20322706250000000000e-3 * t1200 + t1190 - 0.30484059375000000000e-3 * t1261;
1816 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1817 v3rho2sigma[3] = -t528 + t529 + t530 - t514 + t531 - t1151 + t525 - t538 + t540 + t7 * (t1251 + t1263);
1818
1819 t1266 = 0.7e1 / 0.3e1 * t1237;
1820 t1267 = 0.14225894375000000000e-2 * t1239;
1821 t1268 = 0.40645412500000000000e-3 * t1242;
1822 t1269 = -t1223 + t1266 - t1267 + t1268 - t1244 - t1246 - t1160 + t1220 - t1221 - t1168 + t1178 - t1225;
1823 t1270 = 0.40645412500000000000e-3 * t1253;
1824 t1271 = 0.3e1 / 0.2e1 * t1255;
1825 t1272 = 0.60968118750000000000e-3 * t1261;
1826 t1273 = t1226 + t1227 + t1230 - t1233 + t1217 - t1218 + t1270 + t1271 - t1257 - t1191 + t1201 + t1228 - t1272;
1827 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1828 v3rho2sigma[4] = -t1209 + t1210 + t1211 - t513 + t1213 - t1214 + t533 - t537 + t543 + t7 * (t1269 + t1273);
1829
1830 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1831 v3rho2sigma[5] = v3rho2sigma[3];
1832
1833 t1276 = t511 * t496;
1834 t1277 = t523 * t1276;
1835 t1278 = 0.20322706250000000000e-3 * t1277;
1836 t1279 = t510 * t1276;
1837 t1280 = t1279 / 0.2e1;
1838 t1281 = t440 * t483;
1839 t1282 = t510 * t1281;
1840 t1283 = 0.3e1 / 0.4e1 * t1282;
1841 t1284 = t523 * t1281;
1842 t1285 = 0.30484059375000000000e-3 * t1284;
1843 t1286 = t1278 - t1280 - t1244 + t1155 - t1157 + t1163 - t1166 - t1175 - t1182 + t1185 + t1187 + t1190 + t1283 - t1285 + t1270 + t1266 - t1257 - t1267 + t1194 - t1204;
1844 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1845 v3rho2sigma[6] = t7 * t1286 - t1151 - t528 + t529 + t530 + t531 - t537 + t543;
1846
1847 t1290 = 0.40645412500000000000e-3 * t1277;
1848 t1291 = 0.2e1 * t1244;
1849 t1292 = 0.3e1 / 0.2e1 * t1282;
1850 t1293 = 0.60968118750000000000e-3 * t1284;
1851 t1294 = 0.81290825000000000000e-3 * t1253;
1852 t1295 = 0.14e2 / 0.3e1 * t1237;
1853 t1296 = 0.2e1 * t1257;
1854 t1297 = 0.28451788750000000000e-2 * t1239;
1855 t1298 = t1290 - t1279 - t1291 + t1217 - t1218 + t1220 - t1221 - t1223 - t1225 + t1226 + t1227 + t1228 + t1292 - t1293 + t1294 + t1295 - t1296 - t1297 + t1230 - t1233;
1856 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1857 v3rho2sigma[7] = -t1209 + t1210 + t1211 - 0.2e1 * t537 + t1213 - t1214 + 0.81290825000000000000e-3 * t539 + t7 * t1298;
1858
1859 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1860 v3rho2sigma[8] = v3rho2sigma[6];
1861
1862 t1300 = t516 * t546;
1863 t1301 = t1300 * t228;
1864 t1302 = 0.21338841562500000000e-2 * t1301;
1865 t1303 = t547 * t181;
1866 t1304 = t1303 * t228;
1867 t1305 = 0.60968118750000000000e-3 * t1304;
1868 t1306 = t547 * t86;
1869 t1307 = t1306 * t512;
1870 t1308 = 0.30484059375000000000e-3 * t1307;
1871 t1310 = 0.1e1 / t65 / t562;
1872 t1311 = t1310 * t78;
1873 t1312 = t1311 * t228;
1874 t1313 = 0.77095646007107291667e-6 * t1312;
1875 t1314 = t551 * t155;
1876 t1315 = t1314 * t86;
1877 t1316 = t1315 * t188;
1878 t1317 = 0.16520495572951562500e-6 * t1316;
1879 t1318 = t552 * t181;
1880 t1319 = t1318 * t228;
1881 t1320 = 0.16520495572951562500e-6 * t1319;
1882 t1321 = t552 * t86;
1883 t1322 = t1321 * t512;
1884 t1323 = 0.82602477864757812500e-7 * t1322;
1885 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1886 v3rhosigma2[0] = -t549 + t554 + t7 * (t1302 - t1305 + t1308 - t1313 - t1317 + t1320 - t1323);
1887
1888 t1326 = 0.42677683125000000000e-2 * t1301;
1889 t1327 = 0.12193623750000000000e-2 * t1304;
1890 t1328 = 0.60968118750000000000e-3 * t1307;
1891 t1329 = 0.15419129201421458333e-5 * t1312;
1892 t1330 = 0.33040991145903125000e-6 * t1316;
1893 t1331 = 0.33040991145903125000e-6 * t1319;
1894 t1332 = 0.16520495572951562500e-6 * t1322;
1895 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1896 v3rhosigma2[1] = -t556 + t557 + t7 * (t1326 - t1327 + t1328 - t1329 - t1330 + t1331 - t1332);
1897
1898 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1899 v3rhosigma2[2] = v3rhosigma2[0];
1900
1901 t1335 = 0.85355366250000000000e-2 * t1301;
1902 t1336 = 0.243872475e-2 * t1304;
1903 t1337 = 0.12193623750000000000e-2 * t1307;
1904 t1338 = 0.30838258402842916667e-5 * t1312;
1905 t1339 = 0.6608198229180625e-6 * t1316;
1906 t1340 = 0.6608198229180625e-6 * t1319;
1907 t1341 = 0.33040991145903125000e-6 * t1322;
1908 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1909 v3rhosigma2[3] = -t559 + t560 + t7 * (t1335 - t1336 + t1337 - t1338 - t1339 + t1340 - t1341);
1910
1911 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1912 v3rhosigma2[4] = v3rhosigma2[1];
1913
1914 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1915 v3rhosigma2[5] = v3rhosigma2[2];
1916
1917 t1344 = t1306 * t536;
1918 t1345 = 0.30484059375000000000e-3 * t1344;
1919 t1346 = t1321 * t536;
1920 t1347 = 0.82602477864757812500e-7 * t1346;
1921 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1922 v3rhosigma2[6] = -t549 + t554 + t7 * (t1302 - t1305 + t1345 - t1313 - t1317 + t1320 - t1347);
1923
1924 t1350 = 0.60968118750000000000e-3 * t1344;
1925 t1351 = 0.16520495572951562500e-6 * t1346;
1926 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1927 v3rhosigma2[7] = -t556 + t557 + t7 * (t1326 - t1327 + t1350 - t1329 - t1330 + t1331 - t1351);
1928
1929 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1930 v3rhosigma2[8] = v3rhosigma2[6];
1931
1932 t1354 = 0.12193623750000000000e-2 * t1344;
1933 t1355 = 0.33040991145903125000e-6 * t1346;
1934 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1935 v3rhosigma2[9] = -t559 + t560 + t7 * (t1335 - t1336 + t1354 - t1338 - t1339 + t1340 - t1355);
1936
1937 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1938 v3rhosigma2[10] = v3rhosigma2[7];
1939
1940 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1941 v3rhosigma2[11] = v3rhosigma2[8];
1942
1943 t1359 = 0.1e1 / t79 / t56;
1944 t1360 = t225 * t1359;
1945 t1361 = t1360 * t228;
1946 t1362 = 0.30484059375000000000e-3 * t1361;
1947 t1363 = 0.1e1 / t56;
1948 t1364 = t551 * t1363;
1949 t1365 = t1364 * t78;
1950 t1366 = t1365 * t228;
1951 t1367 = 0.24780743359427343750e-6 * t1366;
1952 t1368 = t81 * t81;
1953 t1369 = t1368 * t1368;
1954 t1370 = t1369 * t81;
1955 t1372 = 0.1e1 / t1370 / t562;
1956 t1374 = t1372 * t155 * t546;
1957 t1375 = t1374 * t228;
1958 t1376 = 0.67148235726704010033e-10 * t1375;
1959 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1960 v3sigma3[0] = t7 * (t1362 + t1367 - t1376);
1961
1962 t1378 = 0.60968118750000000000e-3 * t1361;
1963 t1379 = 0.49561486718854687500e-6 * t1366;
1964 t1380 = 0.13429647145340802007e-9 * t1375;
1965 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1966 v3sigma3[1] = t7 * (t1378 + t1379 - t1380);
1967
1968 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1969 v3sigma3[2] = v3sigma3[0];
1970
1971 t1382 = 0.12193623750000000000e-2 * t1361;
1972 t1383 = 0.99122973437709375000e-6 * t1366;
1973 t1384 = 0.26859294290681604013e-9 * t1375;
1974 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1975 v3sigma3[3] = t7 * (t1382 + t1383 - t1384);
1976
1977 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1978 v3sigma3[4] = v3sigma3[1];
1979
1980 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1981 v3sigma3[5] = v3sigma3[2];
1982
1983 t1386 = 0.243872475e-2 * t1361;
1984 t1387 = 0.19824594687541875e-5 * t1366;
1985 t1388 = 0.53718588581363208027e-9 * t1375;
1986 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1987 v3sigma3[6] = t7 * (t1386 + t1387 - t1388);
1988
1989 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1990 v3sigma3[7] = v3sigma3[3];
1991
1992 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1993 v3sigma3[8] = v3sigma3[4];
1994
1995 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1996 v3sigma3[9] = v3sigma3[5];
1997
1998 #ifndef XC_DONT_COMPILE_LXC
1999
2000 if(order < 4) return;
2001
2002
2003 t1390 = 0.3640e4 / 0.27e2 * t567;
2004 t1391 = 0.4e1 * t605;
2005 t1392 = 0.28e2 * t606;
2006 t1393 = 0.280e3 / 0.3e1 * t608;
2007 t1396 = 0.28e2 * t641;
2008 t1398 = 0.4e1 * t647;
2009 t1402 = 0.280e3 / 0.3e1 * t664;
2010 t1403 = 0.28e2 * t667;
2011 t1404 = 0.56e2 * t669;
2012 t1407 = 0.4e1 * t704;
2013 t1408 = -t1390 + t1391 - t1392 + t1393 - 0.6e1 * t611 - 0.2e1 * t638 - t1396 - 0.21e2 * t643 + t1398 + 0.9e1 * t650 - 0.15e2 / 0.2e1 * t659 - 0.140e3 / 0.3e1 * t661 + t1402 - t1403 - t1404 + 0.28e2 * t671 + 0.14e2 * t673 + t1407;
2014 t1409 = 0.12e2 * t705;
2015 t1410 = 0.12e2 * t707;
2016 t1411 = 0.12e2 * t711;
2017 t1413 = 0.4e1 * t748;
2018 t1415 = 0.12e2 * t755;
2019 t1422 = 0.4e1 * t797;
2020 t1427 = t754 * t183 * t512;
2021 t1431 = t378 * t389 * t406 * t196;
2022 t1435 = t243 * t243;
2023 t1441 = t261 * t261;
2024 t1453 = 0.1e1 / t14 / t40 * t564 * t4 * t6 / 0.48e2;
2025 t1455 = 0.1e1 / t562;
2026 t1456 = t728 * t1455;
2027 t1459 = t248 * t249 * t551;
2028 t1462 = t102 * t103 * t317;
2029 t1464 = t6 * t317;
2030 t1465 = t4 * t1464;
2031 t1473 = t4 * t1464 * t20;
2032 t1477 = my_piecewise3(t13, -0.34152e1 / t717 / t17 * t1435 + 0.51228e1 * t718 * t243 * t261 - 0.8538e0 * t242 * t1441 - 0.11384e1 * t722 * t738 + 0.1423e0 * t100 * (-0.29247222222222222222e0 * t1453 + 0.28077333333333333334e1 * t1456 - 0.10399012345679012346e1 * t1459 + 0.90991358024691358027e0 * t1462 + 0.28812345679012345679e0 * t1465), 0.62200000000000000002e-1 * t598 + 0.17283950617283950618e-2 * t1473 - 0.74444444444444444441e-2 * t1465);
2033 t1479 = 0.28e2 * t322 * t710;
2034 t1480 = t319 * t384;
2035 t1483 = 0.280e3 / 0.3e1 * t663 * t374;
2036 t1485 = t87 * t713 * t406;
2037 t1488 = t87 * t383 * t636;
2038 t1490 = t42 * t42;
2039 t1492 = 0.1e1 / t43 / t1490;
2040 t1493 = t296 * t296;
2041 t1499 = t302 * t302;
2042 t1504 = t39 * t1455;
2043 t1506 = -0.24e2 * t598 + 0.24e2 * t1504;
2044 t1509 = t45 * t45;
2045 t1511 = 0.1e1 / t46 / t1509;
2046 t1512 = t306 * t306;
2047 t1518 = t309 * t309;
2048 t1523 = -t1506;
2049 t1531 = t151 * t610;
2050 t1533 = t151 * t637;
2051 t1537 = 0.140e3 / 0.3e1 * t318 * t370 * t185;
2052 t1540 = 0.28e2 / 0.3e1 * t150 * t645 * t185;
2053 t1541 = t319 * t392;
2054 t1543 = t151 * t658;
2055 t1546 = 0.4e1 * t703 * t374;
2056 t1548 = 0.6e1 * t368 * t710;
2057 t1550 = 0.6e1 * t371 * t710;
2058 t1551 = -0.6e1 * t1427 + 0.9e1 * t1431 + t1477 - t1479 - 0.140e3 / 0.3e1 * t1480 + t1483 - 0.3e1 * t1485 - 0.2e1 * t1488 - t87 * t192 * t96 * (0.40e2 / 0.81e2 * t1492 * t1493 - 0.20e2 / 0.9e1 * t613 * t296 * t302 + 0.10e2 / 0.3e1 * t395 * t1499 + 0.40e2 / 0.9e1 * t617 * t622 + 0.5e1 / 0.3e1 * t88 * t1506 + 0.40e2 / 0.81e2 * t1511 * t1512 - 0.20e2 / 0.9e1 * t625 * t306 * t309 + 0.10e2 / 0.3e1 * t400 * t1518 + 0.40e2 / 0.9e1 * t629 * t632 + 0.5e1 / 0.3e1 * t90 * t1523) / 0.2e1 + 0.14e2 * t1531 + 0.14e2 / 0.3e1 * t1533 + t1537 - t1540 + 0.35e2 * t1541 + 0.35e2 / 0.2e1 * t1543 + t1546 + t1548 + t1550;
2059 t1552 = t319 * t407;
2060 t1556 = 0.140e3 / 0.3e1 * t318 * t367 * t185;
2061 t1559 = 0.28e2 / 0.3e1 * t150 * t702 * t185;
2062 t1561 = 0.28e2 * t666 * t374;
2063 t1565 = t332 * t332;
2064 t1581 = t359 * t359;
2065 t1599 = t64 * t5 * t551;
2066 t1621 = t351 * t351;
2067 t1627 = t357 * t357;
2068 t1638 = (0.20106419753086419753e-1 * t1465 + 0.20068888888888888889e-4 * t1599) * t75 - 0.4e1 * t577 * t166 * t171 + 0.12e2 * t343 * t349 * t351 - 0.6e1 * t579 * t357 - 0.24e2 * t163 * t588 * t590 + 0.24e2 * t582 * t593 - 0.4e1 * t345 * t600 + 0.24e2 * t70 / t587 / t74 * t1621 - 0.36e2 * t589 * t351 * t357 + 0.6e1 * t350 * t1627 + 0.8e1 * t350 * t171 * t600 - t167 * (0.75383950617283950617e1 * t1465 + 0.12819753086419753086e1 * t1599 + 0.42335851482216526896e0 * t1455);
2069 t1646 = -0.1950979800e-1 / t676 / t77 * t79 * t83 * t1565 - 0.22761431000000000000e-1 * t678 * t178 * t679 + 0.2926469700e-1 * t678 * t333 * t359 - 0.24658216916666666667e-1 * t331 * t364 * t332 + 0.22761431000000000000e-1 * t331 * t336 * t359 - 0.487744950e-2 * t331 * t83 * t1581 - 0.650326600e-2 * t331 * t174 * t602 - 0.26028117856481481481e-1 * t156 * t699 * t173 + 0.12329108458333333333e-1 * t156 * t364 * t359 - 0.37935718333333333333e-2 * t156 * t178 * t602 + 0.81290825e-3 * t156 * t83 * t1638 - 0.27112622767168209876e-1 * t80 / t81 / t562;
2070 t1648 = t60 * t1646 * t185;
2071 t1649 = t566 * t198;
2072 t1653 = 0.3640e4 / 0.27e2 * t565 * t181 * t185;
2073 t1655 = 0.28e2 * t640 * t374;
2074 t1656 = t151 * t649;
2075 t1659 = 0.4e1 * t646 * t374;
2076 t1660 = t359 * t388;
2077 t1662 = t87 * t1660 * t391;
2078 t1664 = t602 * t191;
2079 t1666 = t87 * t1664 * t197;
2080 t1671 = 0.4e1 * t182 * t86 * t602 * t184;
2081 t1672 = t151 * t714;
2082 t1674 = t173 * t654;
2083 t1676 = t87 * t1674 * t657;
2084 t1678 = t386 * t386;
2085 t1680 = 0.1e1 / t93 / t1678;
2086 t1681 = t77 * t1680;
2087 t1682 = t390 * t390;
2088 t1687 = t405 * t405;
2089 t1692 = t367 * t367;
2090 t1695 = 0.3e1 * t60 * t1692 * t185;
2091 t1696 = t370 * t370;
2092 t1698 = t60 * t1696 * t185;
2093 t1699 = -0.70e2 / 0.3e1 * t1552 + t1556 - t1559 - t1561 + t1648 + 0.1820e4 / 0.27e2 * t1649 - t1653 - t1655 - 0.21e2 * t1656 + t1659 + 0.9e1 / 0.2e1 * t1662 - 0.2e1 * t1666 + t1671 + 0.14e2 * t1672 - 0.15e2 / 0.2e1 * t1676 + 0.105e3 / 0.16e2 * t87 * t1681 * t96 * t1682 + 0.9e1 / 0.4e1 * t87 * t389 * t96 * t1687 + t1695 + t1698;
2094 t1701 = t913 * t766;
2095 t1703 = t378 * t610;
2096 t1709 = t378 * t637;
2097 t1711 = t150 * t749;
2098 t1712 = t1711 * t198;
2099 t1714 = t771 * t384;
2100 t1716 = t771 * t407;
2101 t1721 = 0.28e2 * t771 * t95 * t96 * t367;
2102 t1722 = t318 * t377;
2103 t1723 = t1722 * t198;
2104 t1725 = t771 * t392;
2105 t1727 = t150 * t757;
2106 t1728 = t1727 * t198;
2107 t1730 = t750 * t392;
2108 t1735 = 0.6e1 * t60 * t367 * t370 * t185;
2109 t1736 = t758 * t392;
2110 t1739 = t60 * t645 * t86;
2111 t1740 = t1739 * t198;
2112 t1742 = t378 * t649;
2113 t1744 = t378 * t658;
2114 t1746 = t904 * t766;
2115 t1748 = 0.9e1 * t1701 - 0.6e1 * t1703 + 0.3e1 * t439 * t440 * t635 * t196 - 0.2e1 * t1709 + 0.14e2 * t1712 + 0.28e2 * t1714 + 0.14e2 * t1716 - t1721 - 0.140e3 / 0.3e1 * t1723 - 0.21e2 * t1725 + 0.14e2 * t1728 + 0.9e1 / 0.2e1 * t1730 + t1735 + 0.9e1 / 0.2e1 * t1736 - 0.2e1 * t1740 + 0.9e1 * t1742 - 0.15e2 / 0.2e1 * t1744 - 0.21e2 * t1746;
2116 t1750 = t60 * t702 * t86;
2117 t1751 = t1750 * t198;
2118 t1753 = t750 * t384;
2119 t1755 = t750 * t407;
2120 t1760 = 0.4e1 * t60 * t702 * t181 * t185;
2121 t1762 = 0.12e2 * t754 * t374;
2122 t1763 = t758 * t384;
2123 t1765 = t758 * t407;
2124 t1767 = t378 * t714;
2125 t1774 = t134 * t820 * t53;
2126 t1777 = 0.1e1 / t88 / t1490;
2127 t1790 = 0.1e1 / t90 / t1509;
2128 t1806 = t795 * t144 * t53;
2129 t1809 = t289 * t312 * t53;
2130 t1813 = t274 * t274;
2131 t1819 = t280 * t280;
2132 t1837 = my_piecewise3(t13, -0.20232e1 / t774 / t30 * t1813 + 0.30348e1 * t775 * t274 * t280 - 0.5058e0 * t273 * t1819 - 0.6744e0 * t779 * t786 + 0.843e-1 * t123 * (-0.38836111111111111110e0 * t1453 + 0.37282666666666666666e1 * t1456 - 0.13808395061728395062e1 * t1459 + 0.12082345679012345679e1 * t1462 + 0.22564197530864197531e0 * t1465), 0.31100000000000000001e-1 * t598 + 0.60493827160493827160e-3 * t1473 - 0.32450617283950617284e-2 * t1465);
2133 t1840 = (t1837 - t1477) * t48 * t53;
2134 t1841 = t315 * t57;
2135 t1847 = 0.14560e5 / 0.81e2 * t56 / t8 / t1841 * t86 * t97;
2136 t1850 = t87 * t1638 * t94 * t96;
2137 t1852 = 0.28e2 / 0.3e1 * t151 * t604;
2138 t1854 = 0.140e3 / 0.3e1 * t319 * t381;
2139 t1856 = 0.3640e4 / 0.27e2 * t566 * t188;
2140 t1857 = -0.2e1 * t1751 - 0.6e1 * t1753 - 0.3e1 * t1755 + t1760 + t1762 - 0.6e1 * t1763 - 0.3e1 * t1765 - 0.6e1 * t1767 - 0.45e2 / 0.4e1 * t439 * t899 * t390 * t405 + 0.4e1 * t1774 + t38 * (0.40e2 / 0.81e2 * t1777 * t1493 - 0.16e2 / 0.9e1 * t804 * t296 * t302 + 0.4e1 / 0.3e1 * t295 * t1499 + 0.16e2 / 0.9e1 * t807 * t622 + 0.4e1 / 0.3e1 * t43 * t1506 + 0.40e2 / 0.81e2 * t1790 * t1512 - 0.16e2 / 0.9e1 * t812 * t306 * t309 + 0.4e1 / 0.3e1 * t305 * t1518 + 0.16e2 / 0.9e1 * t815 * t632 + 0.4e1 / 0.3e1 * t46 * t1523) * t53 + 0.4e1 * t1806 + 0.6e1 * t1809 + t1840 + t1847 + t1850 - t1852 + t1854 - t1856;
2141 t1861 = t1409 + t1410 + t1411 - 0.6e1 * t715 + t1413 - 0.6e1 * t751 + t1415 - 0.6e1 * t759 - 0.12e2 * t761 - 0.6e1 * t763 + 0.9e1 * t767 + 0.9e1 * t769 + 0.28e2 * t772 + t1422 + 0.12e2 * t799 + 0.12e2 * t802 + 0.4e1 * t822 + t7 * (t1551 + t1699 + t1748 + t1857);
2142 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2143 v4rho4[0] = t1408 + t1861;
2144
2145 t1862 = 0.9e1 / 0.2e1 * t844;
2146 t1863 = 0.7e1 * t848;
2147 t1864 = 0.3e1 * t854;
2148 t1867 = t1862 - t1390 + t1391 - t1392 + t1393 + t1110 + t1863 - t1111 - t1864 - 0.3e1 / 0.2e1 * t884 - t1112 - 0.3e1 * t611 - t639 - t1396;
2149 t1873 = -0.21e2 / 0.2e1 * t643 + t1398 + 0.9e1 / 0.2e1 * t650 - t660 - 0.35e2 * t661 + t1402 - t1403 - t1404 + 0.21e2 * t671 + 0.7e1 * t673 + t1407 + t1409 + t1410 + t1411;
2150 t1877 = 0.21e2 / 0.2e1 * t905;
2151 t1878 = 0.3e1 * t912;
2152 t1879 = 0.9e1 / 0.2e1 * t914;
2153 t1883 = -0.9e1 / 0.2e1 * t715 + t1413 - t1120 - 0.45e2 / 0.8e1 * t902 - t1877 + t1121 - t1122 - t1123 - t1878 + t1879 + 0.9e1 / 0.2e1 * t918 + 0.9e1 / 0.4e1 * t922 - 0.9e1 / 0.2e1 * t751 + t1415;
2154 t1891 = 0.6e1 * t933;
2155 t1895 = t151 * t850;
2156 t1897 = t151 * t853;
2157 t1898 = 0.7e1 * t1897;
2158 t1899 = t151 * t883;
2159 t1902 = t87 * t1664 * t215;
2160 t1905 = t87 * t713 * t458;
2161 t1906 = 0.3e1 / 0.2e1 * t1905;
2162 t1908 = t87 * t383 * t882;
2163 t1910 = t598 * t141;
2164 t1913 = t598 * t139;
2165 t1917 = 0.40e2 * t454 * t1455;
2166 t1921 = t139 * t302;
2167 t1925 = t299 * t296;
2168 t1931 = t299 * t302;
2169 t1935 = 0.40e2 * t448 * t1455;
2170 t1939 = t141 * t309;
2171 t1943 = t299 * t306;
2172 t1949 = t299 * t309;
2173 t1960 = 0.20e2 * t871 * t1910 - 0.20e2 * t858 * t1913 - t1917 + 0.40e2 / 0.81e2 * t1492 * t203 * t614 - 0.10e2 / 0.9e1 * t855 * t1921 - 0.20e2 / 0.9e1 * t613 * t39 * t1925 + 0.20e2 / 0.3e1 * t395 * t299 * t139 + 0.20e2 / 0.3e1 * t858 * t1931 + t1935 + 0.40e2 / 0.81e2 * t1511 * t205 * t626 - 0.10e2 / 0.9e1 * t868 * t1939 + 0.20e2 / 0.9e1 * t625 * t39 * t1943 - 0.20e2 / 0.3e1 * t400 * t299 * t141 - 0.20e2 / 0.3e1 * t871 * t1949 + 0.10e2 / 0.9e1 * t451 * t632 + 0.20e2 * t90 * t598 + 0.10e2 / 0.9e1 * t445 * t622 - 0.20e2 * t88 * t598;
2174 t1965 = t566 * t216;
2175 t1967 = t319 * t436;
2176 t1969 = t319 * t459;
2177 t1970 = 0.35e2 / 0.3e1 * t1969;
2178 t1971 = -0.9e1 / 0.2e1 * t1427 + 0.9e1 / 0.4e1 * t1431 + t1477 + 0.7e1 / 0.2e1 * t1895 + t1898 + 0.7e1 / 0.2e1 * t1899 - t1902 / 0.2e1 - t1906 - 0.3e1 / 0.2e1 * t1908 - t87 * t192 * t96 * t1960 / 0.2e1 + 0.455e3 / 0.27e2 * t1965 - 0.35e2 / 0.3e1 * t1967 - t1970 - t1479;
2179 t1980 = -0.35e2 * t1480 + t1483 - 0.3e1 / 0.2e1 * t1485 - t1488 / 0.2e1 + 0.7e1 * t1531 + 0.7e1 / 0.6e1 * t1533 + t1537 - t1540 + 0.35e2 / 0.2e1 * t1541 + 0.35e2 / 0.8e1 * t1543 + t1546 + t1548 + t1550 - 0.35e2 / 0.3e1 * t1552 + t1556;
2180 t1988 = -t1559 - t1561 + t1648 + 0.455e3 / 0.9e1 * t1649 - t1653 - t1655 - 0.21e2 / 0.2e1 * t1656 + t1659 + 0.9e1 / 0.4e1 * t1662 - 0.3e1 / 0.2e1 * t1666 + t1671 + 0.21e2 / 0.2e1 * t1672 - 0.15e2 / 0.8e1 * t1676 + t1695 + t1698;
2181 t1989 = t1680 * t96;
2182 t1994 = t758 * t436;
2183 t1996 = t758 * t459;
2184 t1997 = 0.3e1 / 0.2e1 * t1996;
2185 t1998 = t913 * t901;
2186 t2004 = t904 * t901;
2187 t2006 = t1727 * t216;
2188 t2008 = t913 * t921;
2189 t2018 = t1722 * t216;
2190 t2020 = t904 * t921;
2191 t2022 = t1711 * t216;
2192 t2024 = t1750 * t216;
2193 t2026 = t750 * t436;
2194 t2028 = 0.105e3 / 0.16e2 * t439 * t1989 * t214 * t656 - 0.3e1 / 0.2e1 * t1994 - t1997 - 0.45e2 / 0.8e1 * t1998 - 0.45e2 / 0.8e1 * t439 * t899 * t457 * t390 + 0.105e3 / 0.8e1 * t2004 + 0.7e1 / 0.2e1 * t2006 + 0.9e1 / 0.4e1 * t2008 + 0.9e1 / 0.4e1 * t439 * t440 * t457 * t405 + 0.3e1 / 0.4e1 * t439 * t440 * t214 * t635 - 0.35e2 / 0.3e1 * t2018 - 0.21e2 / 0.4e1 * t2020 + 0.7e1 / 0.2e1 * t2022 - t2024 / 0.2e1 - 0.3e1 / 0.2e1 * t2026;
2195 t2031 = t750 * t459;
2196 t2032 = 0.3e1 / 0.2e1 * t2031;
2197 t2034 = t318 * t183 * t442;
2198 t2035 = 0.35e2 / 0.2e1 * t2034;
2199 t2037 = t150 * t373 * t442;
2200 t2038 = 0.21e2 / 0.2e1 * t2037;
2201 t2039 = t771 * t436;
2202 t2041 = t904 * t917;
2203 t2043 = t771 * t459;
2204 t2044 = 0.7e1 * t2043;
2205 t2046 = t60 * t709 * t442;
2206 t2047 = 0.9e1 / 0.4e1 * t2046;
2207 t2048 = t378 * t850;
2208 t2050 = t913 * t917;
2209 t2052 = t378 * t853;
2210 t2053 = 0.3e1 * t2052;
2211 t2058 = t378 * t883;
2212 t2060 = t1739 * t216;
2213 t2064 = -t2032 + t2035 - t2038 + 0.7e1 * t2039 - 0.21e2 / 0.2e1 * t2041 + t2044 + t2047 - 0.3e1 / 0.2e1 * t2048 + 0.9e1 / 0.2e1 * t2050 - t2053 + 0.9e1 / 0.4e1 * t439 * t440 * t881 * t196 - 0.3e1 / 0.2e1 * t2058 - t2060 / 0.2e1 + 0.9e1 / 0.4e1 * t1701 - 0.3e1 * t1703;
2214 t2078 = -t1709 / 0.2e1 + 0.21e2 / 0.2e1 * t1712 + 0.21e2 * t1714 + 0.7e1 * t1716 - t1721 - 0.35e2 * t1723 - 0.21e2 / 0.2e1 * t1725 + 0.21e2 / 0.2e1 * t1728 + 0.9e1 / 0.4e1 * t1730 + t1735 + 0.9e1 / 0.4e1 * t1736 - 0.3e1 / 0.2e1 * t1740 + 0.9e1 / 0.2e1 * t1742 - 0.15e2 / 0.8e1 * t1744 - 0.21e2 / 0.4e1 * t1746;
2215 t2087 = t754 * t183 * t536;
2216 t2090 = t378 * t648 * t842;
2217 t2091 = 0.9e1 / 0.2e1 * t2090;
2218 t2094 = t378 * t389 * t458 * t196;
2219 t2102 = t378 * t389 * t215 * t405;
2220 t2104 = t771 * t843;
2221 t2105 = 0.21e2 / 0.2e1 * t2104;
2222 t2106 = t750 * t843;
2223 t2107 = 0.9e1 / 0.4e1 * t2106;
2224 t2108 = -0.3e1 / 0.2e1 * t1751 - 0.9e1 / 0.2e1 * t1753 - 0.3e1 / 0.2e1 * t1755 + t1760 + t1762 - 0.9e1 / 0.2e1 * t1763 - 0.3e1 / 0.2e1 * t1765 - 0.9e1 / 0.2e1 * t1767 - 0.3e1 / 0.2e1 * t2087 + t2091 + 0.9e1 / 0.2e1 * t2094 - 0.45e2 / 0.8e1 * t439 * t899 * t920 * t196 + 0.9e1 / 0.4e1 * t2102 - t2105 + t2107;
2225 t2112 = t134 * t959 * t53;
2226 t2127 = 0.32e2 * t417 * t1455;
2227 t2142 = 0.32e2 * t423 * t1455;
2228 t2156 = 0.4e1 / 0.9e1 * t420 * t632 + 0.16e2 * t46 * t598 + 0.4e1 / 0.9e1 * t414 * t622 - 0.16e2 * t43 * t598 - 0.8e1 * t938 * t1913 + 0.8e1 * t950 * t1910 + t2127 + 0.40e2 / 0.81e2 * t1790 * t205 * t626 - 0.8e1 / 0.9e1 * t947 * t1939 + 0.16e2 / 0.9e1 * t812 * t39 * t1943 - 0.8e1 / 0.3e1 * t305 * t299 * t141 - 0.8e1 / 0.3e1 * t950 * t1949 - t2142 + 0.40e2 / 0.81e2 * t1777 * t203 * t614 - 0.8e1 / 0.9e1 * t935 * t1921 - 0.16e2 / 0.9e1 * t804 * t39 * t1925 + 0.8e1 / 0.3e1 * t295 * t299 * t139 + 0.8e1 / 0.3e1 * t938 * t1931;
2229 t2160 = t289 * t426 * t53;
2230 t2161 = 0.3e1 * t2160;
2231 t2163 = t795 * t208 * t53;
2232 t2166 = t378 * t655 * t215 * t390;
2233 t2168 = t758 * t843;
2234 t2169 = 0.9e1 / 0.4e1 * t2168;
2235 t2170 = t1774 + 0.3e1 * t1806 + 0.3e1 * t1809 + t1840 + 0.3e1 * t2112 + t38 * t2156 * t53 + t2161 + t2163 + t1847 + t1850 - t1852 + t1854 - t1856 - 0.45e2 / 0.8e1 * t2166 + t2169;
2236 t2175 = -0.9e1 / 0.2e1 * t759 - 0.9e1 * t761 - 0.3e1 * t763 + t768 + 0.9e1 / 0.2e1 * t769 + 0.21e2 * t772 + t1422 + 0.9e1 * t799 + 0.6e1 * t802 + t822 + t1124 + t1891 + 0.3e1 * t961 + t7 * (t1971 + t1980 + t1988 + t2028 + t2064 + t2078 + t2108 + t2170);
2237 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2238 v4rho4[1] = t1867 + t1873 + t1883 + t2175;
2239
2240 t2179 = t134 * t1066 * t53;
2241 t2185 = t139 * t39 * t299;
2242 t2190 = t39 * t39;
2243 t2192 = 0.1e1 / t1841;
2244 t2211 = t141 * t39 * t299;
2245 t2231 = 0.40e2 / 0.81e2 * t1777 * t466 * t296 - 0.64e2 / 0.27e2 * t935 * t2185 - 0.8e1 / 0.27e2 * t1046 * t302 + 0.32e2 / 0.9e1 * t295 * t2190 * t2192 + 0.16e2 / 0.9e1 * t414 * t299 - 0.16e2 / 0.3e1 * t414 * t620 - 0.8e1 / 0.27e2 * t804 * t470 * t296 + 0.8e1 / 0.9e1 * t295 * t990 * t139 + 0.4e1 / 0.9e1 * t1051 * t302 + t2127 + 0.40e2 / 0.81e2 * t1790 * t473 * t306 + 0.64e2 / 0.27e2 * t947 * t2211 - 0.8e1 / 0.27e2 * t1056 * t309 + 0.32e2 / 0.9e1 * t305 * t2190 * t2192 - 0.16e2 / 0.9e1 * t420 * t299 + 0.16e2 / 0.3e1 * t420 * t620 - 0.8e1 / 0.27e2 * t812 * t476 * t306 + 0.8e1 / 0.9e1 * t305 * t1001 * t141 + 0.4e1 / 0.9e1 * t1061 * t309 - t2142;
2246 t2235 = t289 * t479 * t53;
2247 t2237 = t87 * t1660 * t484;
2248 t2240 = t87 * t713 * t497;
2249 t2243 = t87 * t383 * t1005;
2250 t2288 = 0.40e2 / 0.81e2 * t1492 * t466 * t296 - 0.80e2 / 0.27e2 * t855 * t2185 - 0.10e2 / 0.27e2 * t980 * t302 + 0.80e2 / 0.9e1 * t395 * t2190 * t2192 + 0.40e2 / 0.9e1 * t445 * t299 - 0.40e2 / 0.3e1 * t445 * t620 - 0.10e2 / 0.27e2 * t613 * t470 * t296 + 0.20e2 / 0.9e1 * t395 * t990 * t139 + 0.10e2 / 0.9e1 * t985 * t302 + t1935 + 0.40e2 / 0.81e2 * t1511 * t473 * t306 + 0.80e2 / 0.27e2 * t868 * t2211 - 0.10e2 / 0.27e2 * t993 * t309 + 0.80e2 / 0.9e1 * t400 * t2190 * t2192 - 0.40e2 / 0.9e1 * t451 * t299 + 0.40e2 / 0.3e1 * t451 * t620 - 0.10e2 / 0.27e2 * t625 * t476 * t306 + 0.20e2 / 0.9e1 * t400 * t1001 * t141 + 0.10e2 / 0.9e1 * t998 * t309 - t1917;
2251 t2293 = t319 * t498;
2252 t2295 = t319 * t485;
2253 t2297 = t457 * t457;
2254 t2302 = t151 * t974;
2255 t2304 = t151 * t977;
2256 t2306 = t151 * t1006;
2257 t2310 = -0.3e1 * t1427 + 0.2e1 * t2179 + t38 * t2231 * t53 + t2235 + t1477 + 0.3e1 / 0.4e1 * t2237 - t2240 / 0.2e1 - t2243 - t87 * t192 * t96 * t2288 / 0.2e1 - 0.35e2 / 0.9e1 * t2293 + 0.35e2 / 0.6e1 * t2295 + 0.3e1 / 0.2e1 * t87 * t389 * t96 * t2297 - 0.7e1 / 0.2e1 * t2302 + 0.7e1 / 0.3e1 * t2304 + 0.7e1 / 0.3e1 * t2306 + 0.7e1 * t1895 + 0.28e2 / 0.3e1 * t1897;
2258 t2320 = 0.7e1 / 0.3e1 * t1899 - t1902 - 0.2e1 * t1905 - t1908 + 0.910e3 / 0.27e2 * t1965 - 0.70e2 / 0.3e1 * t1967 - 0.140e3 / 0.9e1 * t1969 - t1479 - 0.70e2 / 0.3e1 * t1480 + t1483 - t1485 / 0.2e1 + 0.7e1 / 0.3e1 * t1531 + t1537 - t1540 + 0.35e2 / 0.6e1 * t1541 + t1546 + t1548;
2259 t2327 = t1550 - 0.35e2 / 0.9e1 * t1552 + t1556 - t1559 - t1561 + t1648 + 0.910e3 / 0.27e2 * t1649 - t1653 - t1655 - 0.7e1 / 0.2e1 * t1656 + t1659 + 0.3e1 / 0.4e1 * t1662 - t1666 + t1671 + 0.7e1 * t1672 + t1695 + t1698;
2260 t2345 = -0.3e1 * t1994 - 0.2e1 * t1996 - 0.15e2 / 0.4e1 * t1998 + 0.35e2 / 0.4e1 * t2004 + 0.7e1 * t2006 + 0.3e1 / 0.2e1 * t2008 - 0.70e2 / 0.3e1 * t2018 - 0.7e1 / 0.2e1 * t2020 + 0.7e1 * t2022 - t2024 - 0.3e1 * t2026 - 0.2e1 * t2031 + 0.70e2 / 0.3e1 * t2034 - 0.14e2 * t2037 + 0.14e2 * t2039 - 0.7e1 * t2041 + 0.28e2 / 0.3e1 * t2043 + 0.3e1 * t2046;
2261 t2359 = -0.3e1 * t2048 + 0.3e1 * t2050 - 0.4e1 * t2052 - t2058 - t2060 - t1703 + 0.7e1 * t1712 + 0.14e2 * t1714 + 0.7e1 / 0.3e1 * t1716 - t1721 - 0.70e2 / 0.3e1 * t1723 - 0.7e1 / 0.2e1 * t1725 + 0.7e1 * t1728 + 0.3e1 / 0.4e1 * t1730 + t1735 + 0.3e1 / 0.4e1 * t1736 - t1740;
2262 t2372 = t904 * t1032;
2263 t2374 = t771 * t498;
2264 t2380 = 0.3e1 / 0.2e1 * t1742 - t1751 - 0.3e1 * t1753 - t1755 / 0.2e1 + t1760 + t1762 - 0.3e1 * t1763 - t1765 / 0.2e1 - 0.3e1 * t1767 - 0.3e1 * t2087 + 0.6e1 * t2090 + 0.3e1 * t2094 + 0.3e1 / 0.2e1 * t2102 - 0.14e2 * t2104 + 0.3e1 * t2106 - 0.7e1 / 0.2e1 * t2372 + 0.7e1 / 0.3e1 * t2374 + 0.3e1 / 0.4e1 * t439 * t440 * t496 * t405;
2265 t2382 = t750 * t498;
2266 t2384 = t904 * t1022;
2267 t2386 = t904 * t1026;
2268 t2388 = t771 * t485;
2269 t2398 = t750 * t485;
2270 t2400 = t913 * t1022;
2271 t2402 = t913 * t1026;
2272 t2404 = t378 * t974;
2273 t2406 = t913 * t1032;
2274 t2408 = t378 * t977;
2275 t2413 = t378 * t1006;
2276 t2418 = t758 * t498;
2277 t2424 = -t2382 / 0.2e1 + 0.35e2 / 0.4e1 * t2384 - 0.7e1 * t2386 - 0.7e1 / 0.2e1 * t2388 - 0.15e2 / 0.8e1 * t439 * t899 * t483 * t405 + 0.3e1 / 0.2e1 * t439 * t440 * t214 * t881 + 0.3e1 / 0.4e1 * t2398 - 0.15e2 / 0.4e1 * t2400 + 0.3e1 * t2402 + 0.3e1 / 0.2e1 * t2404 + 0.3e1 / 0.2e1 * t2406 - t2408 + 0.3e1 / 0.2e1 * t439 * t440 * t1004 * t196 - t2413 - 0.15e2 / 0.8e1 * t439 * t899 * t496 * t390 - t2418 / 0.2e1 + 0.105e3 / 0.16e2 * t439 * t1989 * t483 * t390;
2278 t2425 = t758 * t485;
2279 t2435 = t378 * t389 * t497 * t196;
2280 t2443 = t378 * t655 * t484 * t196;
2281 t2447 = t378 * t389 * t215 * t457;
2282 t2449 = 0.3e1 / 0.4e1 * t2425 + 0.2e1 * t1806 + t1809 + t1840 + 0.2e1 * t2112 + 0.4e1 * t2160 + 0.2e1 * t2163 + t1847 + t1850 - t1852 + t1854 - t1856 - 0.15e2 / 0.4e1 * t2166 + 0.3e1 * t2168 + 0.3e1 / 0.2e1 * t2435 - 0.15e2 / 0.2e1 * t439 * t899 * t441 * t457 - 0.15e2 / 0.4e1 * t2443 + 0.3e1 * t2447;
2283 t2462 = t7 * (t2310 + t2320 + t2327 + t2345 + t2359 + t2380 + t2424 + t2449) + 0.6e1 * t844 - t1390 + t1391 - t1392 + t1393 + 0.3e1 / 0.2e1 * t975 - t978 - t1007 - 0.7e1 / 0.2e1 * t1009 + 0.7e1 / 0.3e1 * t1011 + 0.14e2 * t846 + 0.28e2 / 0.3e1 * t848 - 0.3e1 * t851 - 0.4e1 * t854 - t884;
2284 t2469 = -0.70e2 / 0.3e1 * t886 - t611 - t1396 - 0.7e1 / 0.2e1 * t643 + t1398 + 0.3e1 / 0.2e1 * t650 - 0.70e2 / 0.3e1 * t661 + t1402 - t1403 - t1404 + 0.14e2 * t671 + 0.7e1 / 0.3e1 * t673 + t1407 + t1409 + t1410 + t1411;
2285 t2485 = -0.3e1 * t715 + t1413 - 0.15e2 / 0.4e1 * t1023 + 0.3e1 * t1027 + 0.3e1 / 0.2e1 * t1029 + 0.3e1 / 0.2e1 * t1033 - t1035 - 0.3e1 * t897 - 0.15e2 / 0.4e1 * t902 - 0.14e2 * t905 + 0.14e2 * t907 - 0.3e1 * t909 - 0.6e1 * t911 - 0.4e1 * t912 + 0.6e1 * t914 + 0.3e1 * t918;
2286 t2499 = 0.3e1 / 0.2e1 * t922 - 0.3e1 * t751 + t1415 - 0.3e1 * t759 - 0.6e1 * t761 - t763 + 0.3e1 / 0.2e1 * t769 + 0.14e2 * t772 + t1422 + 0.6e1 * t799 + 0.2e1 * t802 + 0.6e1 * t931 + 0.8e1 * t933 + 0.2e1 * t961 + 0.2e1 * t1045 + 0.2e1 * t1068;
2287 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2288 v4rho4[2] = t2462 + t2469 + t2485 + t2499;
2289
2290 t2507 = t1862 - t1390 + t1391 - t1392 + t1393 - t1100 - t1105 + 0.9e1 / 0.2e1 * t975 - 0.3e1 * t978 - 0.3e1 / 0.2e1 * t1007 - 0.21e2 / 0.2e1 * t1009 + 0.7e1 * t1011 + 0.21e2 * t846 + t1863;
2291 t2510 = -0.9e1 / 0.2e1 * t851 - t1864 - 0.35e2 * t886 - t1396 + t1398 - t662 + t1402 - t1403 - t1404 + t672 + t1407 + t1409 + t1410 + t1411;
2292 t2521 = -t716 + t1413 + t1117 - 0.45e2 / 0.8e1 * t1023 + 0.9e1 / 0.2e1 * t1027 + 0.9e1 / 0.2e1 * t1029 + 0.9e1 / 0.4e1 * t1033 - 0.3e1 * t1035 - 0.9e1 / 0.2e1 * t897 - t1877 + 0.21e2 * t907 - 0.9e1 / 0.2e1 * t909 - 0.9e1 * t911 - t1878;
2293 t2534 = t378 * t389 * t215 * t496;
2294 t2537 = t87 * t1674 * t1102;
2295 t2544 = t470 * t139;
2296 t2547 = t299 * t470;
2297 t2557 = 0.12e2 * t598 + 0.24e2 * t1504;
2298 t2565 = t476 * t141;
2299 t2568 = t299 * t476;
2300 t2576 = -t2557;
2301 t2579 = 0.40e2 / 0.81e2 * t1492 * t1079 * t139 - 0.20e2 / 0.9e1 * t980 * t300 - 0.10e2 / 0.9e1 * t855 * t2544 + 0.20e2 / 0.3e1 * t858 * t2547 + 0.10e2 / 0.3e1 * t445 * t990 + 0.10e2 / 0.9e1 * t395 * t1085 * t139 + 0.5e1 / 0.3e1 * t88 * t2557 + 0.40e2 / 0.81e2 * t1511 * t1088 * t141 + 0.20e2 / 0.9e1 * t993 * t300 - 0.10e2 / 0.9e1 * t868 * t2565 - 0.20e2 / 0.3e1 * t871 * t2568 + 0.10e2 / 0.3e1 * t451 * t1001 + 0.10e2 / 0.9e1 * t400 * t1093 * t141 + 0.5e1 / 0.3e1 * t90 * t2576;
2302 t2585 = t87 * t383 * t1097;
2303 t2587 = t151 * t1103;
2304 t2589 = t151 * t1098;
2305 t2594 = -0.3e1 / 0.2e1 * t1427 + 0.3e1 * t2179 + 0.3e1 * t2235 + t1477 - 0.45e2 / 0.8e1 * t439 * t899 * t1114 * t196 + 0.9e1 / 0.4e1 * t2534 - 0.15e2 / 0.8e1 * t2537 - t87 * t192 * t96 * t2579 / 0.2e1 - t2585 / 0.2e1 + 0.35e2 / 0.8e1 * t2587 + 0.7e1 / 0.6e1 * t2589 + 0.9e1 / 0.4e1 * t2237 - 0.3e1 / 0.2e1 * t2240 - 0.3e1 / 0.2e1 * t2243;
2306 t2605 = -0.35e2 / 0.3e1 * t2293 + 0.35e2 / 0.2e1 * t2295 - 0.21e2 / 0.2e1 * t2302 + 0.7e1 * t2304 + 0.7e1 / 0.2e1 * t2306 + 0.21e2 / 0.2e1 * t1895 + t1898 - 0.3e1 / 0.2e1 * t1902 - t1906 + 0.455e3 / 0.9e1 * t1965 - 0.35e2 * t1967 - t1970 - t1479 - 0.35e2 / 0.3e1 * t1480 + t1483;
2307 t2609 = t1537 - t1540 + t1546 + t1548 + t1550 + t1556 - t1559 - t1561 + t1648 + 0.455e3 / 0.27e2 * t1649 - t1653 - t1655 + t1659 - t1666 / 0.2e1 + t1671;
2308 t2618 = 0.7e1 / 0.2e1 * t1672 + t1695 + t1698 - 0.9e1 / 0.2e1 * t1994 - t1997 + 0.21e2 / 0.2e1 * t2006 - 0.35e2 * t2018 + 0.21e2 / 0.2e1 * t2022 - 0.3e1 / 0.2e1 * t2024 - 0.9e1 / 0.2e1 * t2026 - t2032 + t2035 - t2038 + 0.21e2 * t2039 + t2044;
2309 t2630 = t2047 - 0.9e1 / 0.2e1 * t2048 - t2053 - 0.3e1 / 0.2e1 * t2060 + 0.7e1 / 0.2e1 * t1712 + 0.7e1 * t1714 - t1721 - 0.35e2 / 0.3e1 * t1723 + 0.7e1 / 0.2e1 * t1728 + t1735 - t1740 / 0.2e1 - t1751 / 0.2e1 - 0.3e1 / 0.2e1 * t1753 + t1760 + t1762;
2310 t2643 = -0.3e1 / 0.2e1 * t1763 - 0.3e1 / 0.2e1 * t1767 - 0.9e1 / 0.2e1 * t2087 + t2091 - t2105 + t2107 - 0.21e2 / 0.4e1 * t2372 + 0.7e1 * t2374 - 0.3e1 / 0.2e1 * t2382 + 0.105e3 / 0.8e1 * t2384 - 0.21e2 / 0.2e1 * t2386 - 0.21e2 / 0.2e1 * t2388 + 0.9e1 / 0.4e1 * t2398 - 0.45e2 / 0.8e1 * t2400 + 0.9e1 / 0.2e1 * t2402;
2311 t2652 = t913 * t1115;
2312 t2662 = t378 * t1098;
2313 t2668 = 0.9e1 / 0.2e1 * t2404 + 0.9e1 / 0.4e1 * t2406 - 0.3e1 * t2408 - 0.3e1 / 0.2e1 * t2413 - 0.3e1 / 0.2e1 * t2418 + 0.9e1 / 0.4e1 * t2425 + t1806 + t1840 + t2161 + 0.3e1 * t2163 + 0.9e1 / 0.4e1 * t2652 + 0.9e1 / 0.4e1 * t439 * t440 * t457 * t496 + 0.9e1 / 0.4e1 * t439 * t440 * t214 * t1004 - t2662 / 0.2e1 + 0.105e3 / 0.16e2 * t439 * t1989 * t1101 * t196;
2314 t2673 = t378 * t1103;
2315 t2679 = t904 * t1115;
2316 t2682 = t134 * t1138 * t53;
2317 t2715 = 0.40e2 / 0.81e2 * t1777 * t1079 * t139 - 0.16e2 / 0.9e1 * t1046 * t300 - 0.8e1 / 0.9e1 * t935 * t2544 + 0.8e1 / 0.3e1 * t938 * t2547 + 0.4e1 / 0.3e1 * t414 * t990 + 0.4e1 / 0.9e1 * t295 * t1085 * t139 + 0.4e1 / 0.3e1 * t43 * t2557 + 0.40e2 / 0.81e2 * t1790 * t1088 * t141 + 0.16e2 / 0.9e1 * t1056 * t300 - 0.8e1 / 0.9e1 * t947 * t2565 - 0.8e1 / 0.3e1 * t950 * t2568 + 0.4e1 / 0.3e1 * t420 * t1001 + 0.4e1 / 0.9e1 * t305 * t1093 * t141 + 0.4e1 / 0.3e1 * t46 * t2576;
2318 t2721 = -0.45e2 / 0.8e1 * t439 * t899 * t483 * t457 - 0.15e2 / 0.8e1 * t2673 + 0.3e1 / 0.4e1 * t439 * t440 * t1096 * t196 - 0.21e2 / 0.4e1 * t2679 + t2682 + t38 * t2715 * t53 + t1847 + t1850 - t1852 + t1854 - t1856 + t2169 + 0.9e1 / 0.4e1 * t2435 - 0.45e2 / 0.8e1 * t2443 + 0.9e1 / 0.2e1 * t2447;
2319 t2726 = t1879 - t752 + t1415 - t760 - t762 + t773 + t1422 + t800 + 0.9e1 * t931 + t1891 + 0.6e1 * t1045 + 0.3e1 * t1068 + t1140 + t7 * (t2594 + t2605 + t2609 + t2618 + t2630 + t2643 + t2668 + t2721);
2320 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2321 v4rho4[3] = t2507 + t2510 + t2521 + t2726;
2322
2323 t2737 = -t1390 + t1391 - t1392 + t1393 - 0.2e1 * t1099 - 0.15e2 / 0.2e1 * t1104 + 0.9e1 * t975 - 0.6e1 * t978 - 0.21e2 * t1009 + 0.14e2 * t1011 + 0.28e2 * t846 - 0.6e1 * t851 - 0.140e3 / 0.3e1 * t886 - t1396 + t1398 + t1402 - t1403 - t1404;
2324 t2738 = t466 * t466;
2325 t2743 = t470 * t470;
2326 t2749 = 0.24e2 * t598 + 0.24e2 * t1504;
2327 t2752 = t473 * t473;
2328 t2757 = t476 * t476;
2329 t2762 = -t2749;
2330 t2795 = t483 * t483;
2331 t2800 = t496 * t496;
2332 t2816 = t38 * (0.40e2 / 0.81e2 * t1777 * t2738 - 0.16e2 / 0.9e1 * t1046 * t470 + 0.4e1 / 0.3e1 * t295 * t2743 + 0.16e2 / 0.9e1 * t414 * t1085 + 0.4e1 / 0.3e1 * t43 * t2749 + 0.40e2 / 0.81e2 * t1790 * t2752 - 0.16e2 / 0.9e1 * t1056 * t476 + 0.4e1 / 0.3e1 * t305 * t2757 + 0.16e2 / 0.9e1 * t420 * t1093 + 0.4e1 / 0.3e1 * t46 * t2762) * t53 + 0.6e1 * t2235 + t1477 + 0.9e1 * t2534 - t87 * t192 * t96 * (0.40e2 / 0.81e2 * t1492 * t2738 - 0.20e2 / 0.9e1 * t980 * t470 + 0.10e2 / 0.3e1 * t395 * t2743 + 0.40e2 / 0.9e1 * t445 * t1085 + 0.5e1 / 0.3e1 * t88 * t2749 + 0.40e2 / 0.81e2 * t1511 * t2752 - 0.20e2 / 0.9e1 * t993 * t476 + 0.10e2 / 0.3e1 * t400 * t2757 + 0.40e2 / 0.9e1 * t451 * t1093 + 0.5e1 / 0.3e1 * t90 * t2762) / 0.2e1 + 0.105e3 / 0.16e2 * t87 * t1681 * t96 * t2795 + 0.9e1 / 0.4e1 * t87 * t389 * t96 * t2800 - 0.15e2 / 0.2e1 * t2537 - 0.2e1 * t2585 + 0.35e2 / 0.2e1 * t2587 + 0.14e2 / 0.3e1 * t2589 + 0.9e1 / 0.2e1 * t2237 - 0.3e1 * t2240 - 0.70e2 / 0.3e1 * t2293 + 0.35e2 * t2295 - 0.21e2 * t2302 + 0.14e2 * t2304 + 0.14e2 * t1895;
2333 t2820 = -0.2e1 * t1902 + 0.1820e4 / 0.27e2 * t1965 - 0.140e3 / 0.3e1 * t1967 - t1479 + t1483 + t1537 - t1540 + t1546 + t1548 + t1550 + t1556 - t1559 - t1561 + t1648 - t1653 - t1655 + t1659 + t1671 + t1695;
2334 t2835 = t1698 - 0.6e1 * t1994 + 0.14e2 * t2006 - 0.140e3 / 0.3e1 * t2018 + 0.14e2 * t2022 - 0.2e1 * t2024 - 0.6e1 * t2026 + 0.28e2 * t2039 - 0.6e1 * t2048 - 0.2e1 * t2060 - t1721 + t1735 + t1760 + t1762 - 0.6e1 * t2087 + 0.14e2 * t2374 - 0.3e1 * t2382 - 0.21e2 * t2388;
2335 t2855 = 0.9e1 / 0.2e1 * t2398 + 0.9e1 * t2404 - 0.6e1 * t2408 - 0.3e1 * t2418 + 0.9e1 / 0.2e1 * t2425 + t1840 + 0.4e1 * t2163 + 0.9e1 * t2652 - 0.2e1 * t2662 - 0.15e2 / 0.2e1 * t2673 - 0.21e2 * t2679 - 0.45e2 / 0.4e1 * t439 * t899 * t483 * t496 + 0.3e1 * t439 * t440 * t1096 * t214 + 0.4e1 * t2682 + t1847 + t1850 - t1852 + t1854 - t1856;
2336 t2869 = t1407 + t1409 + t1410 + t1411 + t1413 + t7 * (t2816 + t2820 + t2835 + t2855) + 0.9e1 * t1116 + 0.9e1 * t1029 - 0.6e1 * t1035 - 0.6e1 * t897 + 0.28e2 * t907 - 0.6e1 * t909 - 0.12e2 * t911 + t1415 + t1422 + 0.12e2 * t931 + 0.12e2 * t1045 + 0.4e1 * t1140;
2337 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2338 v4rho4[4] = t2737 + t2869;
2339
2340 t2873 = t79 / t223 / t1841 * t228;
2341 t2874 = 0.35209088578125000000e-1 * t2873;
2342 t2875 = t221 * t604;
2343 t2876 = t523 * t766;
2344 t2879 = t520 * t1198 * t406;
2345 t2881 = t1161 * t1199;
2346 t2883 = t1164 * t1199;
2347 t2885 = t377 * t184;
2348 t2886 = t1164 * t2885;
2349 t2887 = 0.12193623750000000000e-2 * t2886;
2350 t2888 = t507 * t384;
2351 t2890 = t510 * t766;
2352 t2892 = t507 * t392;
2353 t2894 = t1193 * t198;
2354 t2896 = t86 * t388;
2355 t2898 = t520 * t2896 * t391;
2356 t2900 = t2874 + t2875 - 0.91452178125000000000e-3 * t2876 + 0.60968118750000000000e-3 * t2879 - 0.42677683125000000000e-2 * t2881 + 0.60968118750000000000e-3 * t2883 - t2887 - 0.3e1 * t2888 + 0.9e1 / 0.4e1 * t2890 + 0.9e1 / 0.4e1 * t2892 - 0.3e1 / 0.2e1 * t2894 - 0.91452178125000000000e-3 * t2898;
2357 t2901 = t1202 * t1199;
2358 t2903 = t507 * t407;
2359 t2905 = t1180 * t198;
2360 t2907 = t1184 * t198;
2361 t2910 = t1183 * t181 * t185;
2362 t2911 = 0.3e1 * t2910;
2363 t2912 = t1188 * t188;
2364 t2913 = 0.70e2 / 0.3e1 * t2912;
2365 t2914 = t503 * t381;
2366 t2915 = 0.7e1 * t2914;
2367 t2916 = t221 * t359;
2368 t2917 = t2916 * t512;
2369 t2919 = t507 * t381;
2370 t2920 = 0.3e1 * t2919;
2371 t2921 = t511 * t635;
2372 t2922 = t523 * t2921;
2373 t2925 = t226 * t645 * t228;
2374 t2926 = 0.40645412500000000000e-3 * t2925;
2375 t2927 = t899 * t656;
2376 t2928 = t523 * t2927;
2377 t2930 = 0.60968118750000000000e-3 * t2901 - 0.3e1 / 0.2e1 * t2903 + 0.7e1 * t2905 - 0.3e1 / 0.2e1 * t2907 + t2911 + t2913 - t2915 - 0.3e1 / 0.2e1 * t2917 + t2920 + 0.20322706250000000000e-3 * t2922 - t2926 + 0.76210148437500000000e-3 * t2928;
2378 t2932 = t1188 * t77;
2379 t2933 = t2932 * t512;
2380 t2935 = t1176 * t1195;
2381 t2939 = t317 * t181 * t86 * t97;
2382 t2940 = 0.70e2 / 0.3e1 * t2939;
2383 t2943 = t149 * t370 * t86 * t97;
2384 t2944 = 0.7e1 * t2943;
2385 t2945 = t1193 * t188;
2386 t2946 = 0.3e1 * t2945;
2387 t2949 = t59 * t645 * t86 * t97;
2388 t2950 = t510 * t2927;
2389 t2952 = t1158 * t1152;
2390 t2954 = t1173 * t86;
2391 t2955 = t2954 * t512;
2392 t2957 = t1158 * t1195;
2393 t2960 = t1173 * t181 * t228;
2394 t2961 = 0.19204957406250000000e-1 * t2960;
2395 t2963 = t517 * t370 * t228;
2396 t2964 = 0.42677683125000000000e-2 * t2963;
2397 t2965 = -0.35e2 / 0.3e1 * t2933 - 0.21e2 / 0.4e1 * t2935 + t2940 - t2944 + t2946 + t2949 - 0.15e2 / 0.8e1 * t2950 - 0.21338841562500000000e-2 * t2952 + 0.96024787031250000000e-2 * t2955 + 0.32008262343750000000e-2 * t2957 - t2961 + t2964;
2398 t2966 = t1167 * t1195;
2399 t2968 = t1176 * t1152;
2400 t2970 = t1184 * t188;
2401 t2971 = 0.3e1 * t2970;
2402 t2974 = t59 * t702 * t86 * t97;
2403 t2975 = t510 * t2921;
2404 t2979 = t149 * t367 * t86 * t97;
2405 t2980 = 0.7e1 * t2979;
2406 t2981 = t1167 * t1152;
2407 t2984 = t226 * t702 * t228;
2408 t2985 = 0.40645412500000000000e-3 * t2984;
2409 t2986 = t503 * t173;
2410 t2987 = t2986 * t512;
2411 t2989 = t1180 * t188;
2412 t2990 = 0.14e2 * t2989;
2413 t2992 = t517 * t367 * t228;
2414 t2993 = 0.42677683125000000000e-2 * t2992;
2415 t2995 = t564 * t86 * t97;
2416 t2996 = 0.910e3 / 0.27e2 * t2995;
2417 t2997 = 0.9e1 / 0.4e1 * t2966 + 0.7e1 / 0.2e1 * t2968 + t2971 + t2974 - t2975 / 0.2e1 - t2980 - 0.3e1 / 0.2e1 * t2981 - t2985 + 0.7e1 * t2987 - t2990 + t2993 - t2996;
2418 t3002 = 0.3e1 * t1155;
2419 t3003 = 0.14e2 * t1156;
2420 t3005 = 0.85355366250000000000e-2 * t1162;
2421 t3006 = 0.12193623750000000000e-2 * t1165;
2422 t3009 = 0.19204957406250000000e-1 * t1174;
2423 t3012 = 0.14e2 * t1181;
2424 t3013 = 0.3e1 * t1185;
2425 t3014 = 0.6e1 * t1186;
2426 t3015 = 0.70e2 / 0.3e1 * t1189;
2427 t3017 = 0.3e1 * t1194;
2428 t3020 = 0.12193623750000000000e-2 * t1203;
2429 t3022 = 0.7e1 * t1177 - t3012 + t3013 + t3014 + t3015 - 0.3e1 * t1191 + t3017 + 0.9e1 / 0.4e1 * t1196 + 0.12193623750000000000e-2 * t1200 - t3020 - 0.91452178125000000000e-3 * t1205;
2430 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2431 v4rho3sigma[0] = t7 * (t2900 + t2930 + t2965 + t2997) + 0.60968118750000000000e-3 * t1153 + t3002 - t3003 - 0.42677683125000000000e-2 * t1159 + t3005 - t3006 - 0.3e1 * t1168 - 0.3e1 / 0.2e1 * t1169 - t3009 + t3022;
2432
2433 t3023 = 0.70418177156250000000e-1 * t2873;
2434 t3024 = 0.2e1 * t2875;
2435 t3029 = 0.243872475e-2 * t2886;
2436 t3035 = t3023 + t3024 - 0.18290435625000000000e-2 * t2876 + 0.12193623750000000000e-2 * t2879 - 0.85355366250000000000e-2 * t2881 + 0.12193623750000000000e-2 * t2883 - t3029 - 0.6e1 * t2888 + 0.9e1 / 0.2e1 * t2890 + 0.9e1 / 0.2e1 * t2892 - 0.3e1 * t2894 - 0.18290435625000000000e-2 * t2898;
2437 t3040 = 0.6e1 * t2910;
2438 t3041 = 0.140e3 / 0.3e1 * t2912;
2439 t3042 = 0.14e2 * t2914;
2440 t3044 = 0.6e1 * t2919;
2441 t3046 = 0.81290825e-3 * t2925;
2442 t3048 = 0.12193623750000000000e-2 * t2901 - 0.3e1 * t2903 + 0.14e2 * t2905 - 0.3e1 * t2907 + t3040 + t3041 - t3042 - 0.3e1 * t2917 + t3044 + 0.40645412500000000000e-3 * t2922 - t3046 + 0.15242029687500000000e-2 * t2928;
2443 t3052 = 0.140e3 / 0.3e1 * t2939;
2444 t3053 = 0.14e2 * t2943;
2445 t3054 = 0.6e1 * t2945;
2446 t3055 = 0.2e1 * t2949;
2447 t3060 = 0.38409914812500000000e-1 * t2960;
2448 t3061 = 0.85355366250000000000e-2 * t2963;
2449 t3062 = -0.70e2 / 0.3e1 * t2933 - 0.21e2 / 0.2e1 * t2935 + t3052 - t3053 + t3054 + t3055 - 0.15e2 / 0.4e1 * t2950 - 0.42677683125000000000e-2 * t2952 + 0.19204957406250000000e-1 * t2955 + 0.64016524687500000000e-2 * t2957 - t3060 + t3061;
2450 t3065 = 0.6e1 * t2970;
2451 t3066 = 0.2e1 * t2974;
2452 t3067 = 0.14e2 * t2979;
2453 t3069 = 0.81290825e-3 * t2984;
2454 t3071 = 0.28e2 * t2989;
2455 t3072 = 0.85355366250000000000e-2 * t2992;
2456 t3073 = 0.1820e4 / 0.27e2 * t2995;
2457 t3074 = 0.9e1 / 0.2e1 * t2966 + 0.7e1 * t2968 + t3065 + t3066 - t2975 - t3067 - 0.3e1 * t2981 - t3069 + 0.14e2 * t2987 - t3071 + t3072 - t3073;
2458 t3079 = 0.6e1 * t1155;
2459 t3080 = 0.28e2 * t1156;
2460 t3082 = 0.17071073250000000000e-1 * t1162;
2461 t3083 = 0.243872475e-2 * t1165;
2462 t3086 = 0.38409914812500000000e-1 * t1174;
2463 t3089 = 0.28e2 * t1181;
2464 t3090 = 0.6e1 * t1185;
2465 t3091 = 0.12e2 * t1186;
2466 t3092 = 0.140e3 / 0.3e1 * t1189;
2467 t3094 = 0.6e1 * t1194;
2468 t3097 = 0.243872475e-2 * t1203;
2469 t3099 = 0.14e2 * t1177 - t3089 + t3090 + t3091 + t3092 - 0.6e1 * t1191 + t3094 + 0.9e1 / 0.2e1 * t1196 + 0.24387247500000000000e-2 * t1200 - t3097 - 0.18290435625000000000e-2 * t1205;
2470 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2471 v4rho3sigma[1] = t7 * (t3035 + t3048 + t3062 + t3074) + 0.12193623750000000000e-2 * t1153 + t3079 - t3080 - 0.85355366250000000000e-2 * t1159 + t3082 - t3083 - 0.6e1 * t1168 - 0.3e1 * t1169 - t3086 + t3099;
2472
2473 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2474 v4rho3sigma[2] = v4rho3sigma[0];
2475
2476 t3100 = -t3009 + t1266 - t1267 + t1268 - t1244 - t1246 + t1154 - t1219 + t3005 - t3006 - t1222 - t1170 + t1224 - t3012 + t3013;
2477 t3101 = t1161 * t1252;
2478 t3103 = t1202 * t1252;
2479 t3105 = t510 * t901;
2480 t3107 = t1193 * t216;
2481 t3109 = t523 * t901;
2482 t3111 = t1180 * t216;
2483 t3113 = t510 * t921;
2484 t3115 = t1184 * t216;
2485 t3117 = t1158 * t442;
2486 t3118 = 0.21338841562500000000e-2 * t3117;
2487 t3119 = t523 * t921;
2488 t3121 = t1164 * t1252;
2489 t3123 = t1176 * t442;
2490 t3124 = 0.7e1 / 0.2e1 * t3123;
2491 t3125 = t523 * t917;
2492 t3128 = t520 * t1198 * t458;
2493 t3129 = 0.40645412500000000000e-3 * t3128;
2494 t3130 = t1167 * t442;
2495 t3131 = 0.3e1 / 0.2e1 * t3130;
2496 t3132 = t507 * t436;
2497 t3133 = -0.14225894375000000000e-2 * t3101 + 0.20322706250000000000e-3 * t3103 - 0.15e2 / 0.8e1 * t3105 - t3107 / 0.2e1 + 0.76210148437500000000e-3 * t3109 + t2874 + 0.7e1 / 0.3e1 * t3111 + 0.3e1 / 0.4e1 * t3113 - t3115 / 0.2e1 + t3118 - 0.30484059375000000000e-3 * t3119 + 0.20322706250000000000e-3 * t3121 - t3124 - 0.60968118750000000000e-3 * t3125 + t3129 + t3131 - t3132;
2498 t3134 = t510 * t917;
2499 t3136 = t507 * t459;
2500 t3138 = t226 * t377 * t442;
2501 t3139 = 0.60968118750000000000e-3 * t3138;
2502 t3141 = t506 * t183 * t442;
2503 t3142 = 0.3e1 / 0.2e1 * t3141;
2504 t3144 = 0.28451788750000000000e-2 * t2881;
2505 t3145 = 0.40645412500000000000e-3 * t2883;
2506 t3146 = 0.2e1 * t2888;
2507 t3149 = 0.40645412500000000000e-3 * t2901;
2508 t3151 = 0.14e2 / 0.3e1 * t2905;
2509 t3152 = 0.3e1 / 0.2e1 * t3134 - t3136 + t2875 - t3139 + t3142 + 0.20322706250000000000e-3 * t2879 - t3144 + t3145 - t2887 - t3146 + 0.3e1 / 0.4e1 * t2892 - t2894 - 0.30484059375000000000e-3 * t2898 + t3149 - t2903 / 0.2e1 + t3151 - t2907 + t2911;
2510 t3154 = t511 * t881;
2511 t3155 = t523 * t3154;
2512 t3157 = t2916 * t536;
2513 t3159 = t1167 * t1241;
2514 t3160 = t510 * t3154;
2515 t3162 = t2954 * t536;
2516 t3164 = t2932 * t536;
2517 t3166 = t1158 * t1241;
2518 t3167 = 0.14225894375000000000e-2 * t3166;
2519 t3168 = t2986 * t536;
2520 t3170 = t1176 * t1241;
2521 t3171 = 0.7e1 / 0.3e1 * t3170;
2522 t3172 = 0.70e2 / 0.9e1 * t2933;
2523 t3174 = t2913 - t2915 + 0.20322706250000000000e-3 * t3155 - t3157 / 0.2e1 - t3159 - t3160 / 0.2e1 + 0.32008262343750000000e-2 * t3162 - 0.35e2 / 0.9e1 * t3164 - t3167 + 0.7e1 / 0.3e1 * t3168 + t3171 - t2917 + t2920 - t2926 - t3172 - 0.7e1 / 0.4e1 * t2935 + t2940 - t2944;
2524 t3176 = 0.64016524687500000000e-2 * t2955;
2525 t3181 = 0.14e2 / 0.3e1 * t2987;
2526 t3182 = t2946 + t2949 - 0.71129471875000000000e-3 * t2952 + t3176 + 0.10669420781250000000e-2 * t2957 - t2961 + t2964 + 0.3e1 / 0.4e1 * t2966 + 0.7e1 / 0.6e1 * t2968 + t2971 + t2974 - t2980 - t2981 / 0.2e1 - t2985 + t3181 - t2990 + t2993 - t2996;
2527 t3186 = t3014 + t3017 + t1197 - t3020 - t1206 + t3002 - t3003 + t1270 + t1271 - t1257 - t1229 + t1232 + t7 * (t3133 + t3152 + t3174 + t3182) + t3015 - t1272;
2528 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2529 v4rho3sigma[3] = t3100 + t3186;
2530
2531 t3187 = 0.81290825000000000000e-3 * t1242;
2532 t3188 = 0.2e1 * t1246;
2533 t3192 = -t3086 + t1295 - t1297 + t3187 - t1291 - t3188 + t1216 - 0.56903577500000000000e-2 * t1159 + t3082 - t3083 - 0.4e1 * t1168 - t1169 + 0.28e2 / 0.3e1 * t1177 - t3089 + t3090;
2534 t3193 = 0.3e1 * t1255;
2535 t3196 = 0.28451788750000000000e-2 * t3101;
2536 t3197 = 0.40645412500000000000e-3 * t3103;
2537 t3200 = 0.14e2 / 0.3e1 * t3111;
2538 t3202 = 0.42677683125000000000e-2 * t3117;
2539 t3204 = 0.40645412500000000000e-3 * t3121;
2540 t3205 = 0.7e1 * t3123;
2541 t3207 = 0.81290825000000000000e-3 * t3128;
2542 t3208 = 0.3e1 * t3130;
2543 t3209 = 0.2e1 * t3132;
2544 t3210 = -t3196 + t3197 - 0.15e2 / 0.4e1 * t3105 - t3107 + 0.15242029687500000000e-2 * t3109 + t3023 + t3200 + 0.3e1 / 0.2e1 * t3113 - t3115 + t3202 - 0.60968118750000000000e-3 * t3119 + t3204 - t3205 - 0.12193623750000000000e-2 * t3125 + t3207 + t3208 - t3209;
2545 t3212 = 0.2e1 * t3136;
2546 t3213 = 0.12193623750000000000e-2 * t3138;
2547 t3214 = 0.3e1 * t3141;
2548 t3225 = 0.3e1 * t3134 - t3212 + t3024 - t3213 + t3214 + 0.40645412500000000000e-3 * t2879 - 0.56903577500000000000e-2 * t2881 + 0.81290825000000000000e-3 * t2883 - t3029 - 0.4e1 * t2888 + 0.3e1 / 0.2e1 * t2892 - 0.2e1 * t2894 - 0.60968118750000000000e-3 * t2898 + 0.81290825000000000000e-3 * t2901 - t2903 + 0.28e2 / 0.3e1 * t2905 - 0.2e1 * t2907 + t3040;
2549 t3228 = 0.2e1 * t3159;
2550 t3229 = 0.64016524687500000000e-2 * t3162;
2551 t3230 = 0.70e2 / 0.9e1 * t3164;
2552 t3231 = 0.28451788750000000000e-2 * t3166;
2553 t3232 = 0.14e2 / 0.3e1 * t3168;
2554 t3233 = 0.14e2 / 0.3e1 * t3170;
2555 t3237 = t3041 - t3042 + 0.40645412500000000000e-3 * t3155 - t3157 - t3228 - t3160 + t3229 - t3230 - t3231 + t3232 + t3233 - 0.2e1 * t2917 + t3044 - t3046 - 0.140e3 / 0.9e1 * t2933 - 0.7e1 / 0.2e1 * t2935 + t3052 - t3053;
2556 t3244 = t3054 + t3055 - 0.14225894375000000000e-2 * t2952 + 0.12803304937500000000e-1 * t2955 + 0.21338841562500000000e-2 * t2957 - t3060 + t3061 + 0.3e1 / 0.2e1 * t2966 + 0.7e1 / 0.3e1 * t2968 + t3065 + t3066 - t3067 - t2981 - t3069 + 0.28e2 / 0.3e1 * t2987 - t3071 + t3072 - t3073;
2557 t3248 = 0.12193623750000000000e-2 * t1261;
2558 t3249 = t3091 + t3094 + t1231 - t3097 - t1234 + t3079 - t3080 + t1294 + t3193 - t1296 - 0.4e1 * t1191 + 0.16258165000000000000e-2 * t1200 + t7 * (t3210 + t3225 + t3237 + t3244) + t3092 - t3248;
2559 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2560 v4rho3sigma[4] = t3192 + t3249;
2561
2562 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2563 v4rho3sigma[5] = v4rho3sigma[3];
2564
2565 t3250 = -t3009 - t1280 + t1283 - t1285 + t1278 + t1295 - t1297 + t1268 - t1291 - t1246 - t1160 + t3005 - t3006 - t1168 + t1178;
2566 t3251 = t523 * t1032;
2567 t3254 = t520 * t1198 * t497;
2568 t3256 = t510 * t1032;
2569 t3258 = -t3196 + t3197 - t3107 + t2874 + t3200 - t3115 + t3118 + t3204 - t3124 + t3129 + t3131 - t3209 - t3136 + t2875 - 0.30484059375000000000e-3 * t3251 + 0.20322706250000000000e-3 * t3254 + 0.3e1 / 0.4e1 * t3256;
2570 t3259 = t507 * t498;
2571 t3261 = t510 * t1022;
2572 t3263 = t510 * t1026;
2573 t3265 = t507 * t485;
2574 t3267 = t523 * t1022;
2575 t3269 = t523 * t1026;
2576 t3272 = t520 * t2896 * t484;
2577 t3280 = -t3259 / 0.2e1 - 0.15e2 / 0.8e1 * t3261 + 0.3e1 / 0.2e1 * t3263 + 0.3e1 / 0.4e1 * t3265 + 0.76210148437500000000e-3 * t3267 - 0.60968118750000000000e-3 * t3269 - 0.30484059375000000000e-3 * t3272 - t3139 + t3142 - 0.14225894375000000000e-2 * t2881 + 0.20322706250000000000e-3 * t2883 - t2887 - t2888 - t2894 / 0.2e1 + 0.20322706250000000000e-3 * t2901 + 0.7e1 / 0.3e1 * t2905 - t2907 / 0.2e1 + t2911;
2578 t3282 = t1176 * t1276;
2579 t3284 = t1176 * t1281;
2580 t3286 = t1158 * t1281;
2581 t3288 = t1167 * t1281;
2582 t3290 = t511 * t1004;
2583 t3291 = t523 * t3290;
2584 t3293 = t1167 * t1276;
2585 t3295 = t510 * t3290;
2586 t3297 = t1158 * t1276;
2587 t3300 = t2913 - t2915 + 0.7e1 / 0.6e1 * t3282 - 0.7e1 / 0.4e1 * t3284 + 0.10669420781250000000e-2 * t3286 + 0.3e1 / 0.4e1 * t3288 + 0.20322706250000000000e-3 * t3291 - t3293 / 0.2e1 - t3295 / 0.2e1 - 0.71129471875000000000e-3 * t3297 - t3157 - t3159 + t3229 - t3230 - t3167 + t3232 + t3171 - t2917 / 0.2e1;
2588 t3304 = t2920 - t2926 - 0.35e2 / 0.9e1 * t2933 + t2940 - t2944 + t2946 + t2949 + 0.32008262343750000000e-2 * t2955 - t2961 + t2964 + t2971 + t2974 - t2980 - t2985 + 0.7e1 / 0.3e1 * t2987 - t2990 + t2993 - t2996;
2589 t3308 = -t3012 + t3013 + t3014 + t3017 - t3020 + t3002 - t3003 + t1294 + t1271 - t1296 - t1191 + t1201 + t7 * (t3258 + t3280 + t3300 + t3304) + t3015 - t1272;
2590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2591 v4rho3sigma[6] = t3250 + t3308;
2592
2593 t3312 = -t3086 - t1279 + t1292 - t1293 + t1290 + 0.28e2 / 0.3e1 * t1237 - 0.56903577500000000000e-2 * t1239 + t3187 - 0.4e1 * t1244 - t3188 - t1219 + t3082 - t3083 - t1222 + t1224;
2594 t3325 = -0.56903577500000000000e-2 * t3101 + 0.81290825000000000000e-3 * t3103 - 0.2e1 * t3107 + t3023 + 0.28e2 / 0.3e1 * t3111 - 0.2e1 * t3115 + t3202 + 0.81290825000000000000e-3 * t3121 - t3205 + t3207 + t3208 - 0.4e1 * t3132 - t3212 + t3024 - 0.60968118750000000000e-3 * t3251 + 0.40645412500000000000e-3 * t3254 + 0.3e1 / 0.2e1 * t3256;
2595 t3332 = -t3259 - 0.15e2 / 0.4e1 * t3261 + 0.3e1 * t3263 + 0.3e1 / 0.2e1 * t3265 + 0.15242029687500000000e-2 * t3267 - 0.12193623750000000000e-2 * t3269 - 0.60968118750000000000e-3 * t3272 - t3213 + t3214 - t3144 + t3145 - t3029 - t3146 - t2894 + t3149 + t3151 - t2907 + t3040;
2596 t3344 = t3041 - t3042 + 0.7e1 / 0.3e1 * t3282 - 0.7e1 / 0.2e1 * t3284 + 0.21338841562500000000e-2 * t3286 + 0.3e1 / 0.2e1 * t3288 + 0.40645412500000000000e-3 * t3291 - t3293 - t3295 - 0.14225894375000000000e-2 * t3297 - 0.2e1 * t3157 - t3228 + 0.12803304937500000000e-1 * t3162 - 0.140e3 / 0.9e1 * t3164 - t3231 + 0.28e2 / 0.3e1 * t3168 + t3233 - t2917;
2597 t3345 = t3044 - t3046 - t3172 + t3052 - t3053 + t3054 + t3055 + t3176 - t3060 + t3061 + t3065 + t3066 - t3067 - t3069 + t3181 - t3071 + t3072 - t3073;
2598 t3349 = -t3089 + t3090 + t3091 + t3094 - t3097 + t3079 - t3080 + 0.16258165000000000000e-2 * t1253 + t3193 - 0.4e1 * t1257 - t1229 + t1232 + t7 * (t3325 + t3332 + t3344 + t3345) + t3092 - t3248;
2599 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2600 v4rho3sigma[7] = t3312 + t3349;
2601
2602 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2603 v4rho3sigma[8] = v4rho3sigma[6];
2604
2605 t3360 = -0.42677683125000000000e-2 * t3101 + 0.60968118750000000000e-3 * t3103 - 0.3e1 / 0.2e1 * t3107 + t2874 + 0.7e1 * t3111 - 0.3e1 / 0.2e1 * t3115 + 0.60968118750000000000e-3 * t3121 - 0.3e1 * t3132 + t2875 + 0.60968118750000000000e-3 * t3254 - 0.3e1 / 0.2e1 * t3259 + 0.9e1 / 0.4e1 * t3265;
2606 t3362 = t511 * t1096;
2607 t3363 = t523 * t3362;
2608 t3365 = t510 * t3362;
2609 t3367 = t899 * t1101;
2610 t3368 = t510 * t3367;
2611 t3370 = t523 * t3367;
2612 t3375 = -0.91452178125000000000e-3 * t3272 - t2887 + t2911 + t2913 - t2915 + 0.20322706250000000000e-3 * t3363 - t3365 / 0.2e1 - 0.15e2 / 0.8e1 * t3368 + 0.76210148437500000000e-3 * t3370 + 0.7e1 / 0.2e1 * t3282 - 0.21e2 / 0.4e1 * t3284 + 0.32008262343750000000e-2 * t3286;
2613 t3384 = 0.9e1 / 0.4e1 * t3288 - 0.3e1 / 0.2e1 * t3293 - 0.21338841562500000000e-2 * t3297 - 0.3e1 / 0.2e1 * t3157 + 0.96024787031250000000e-2 * t3162 - 0.35e2 / 0.3e1 * t3164 + 0.7e1 * t3168 + t2920 - t2926 + t2940 - t2944 + t2946;
2614 t3385 = t523 * t1115;
2615 t3387 = t510 * t1115;
2616 t3389 = t2949 - t2961 + t2964 + t2971 + t2974 - t2980 - t2985 - t2990 + t2993 - 0.91452178125000000000e-3 * t3385 + 0.9e1 / 0.4e1 * t3387 - t2996;
2617 t3403 = t3013 + t3014 + t3015 + 0.9e1 / 0.4e1 * t1282 - 0.91452178125000000000e-3 * t1284 + 0.12193623750000000000e-2 * t1253 + 0.7e1 * t1237 - 0.3e1 * t1257 - 0.42677683125000000000e-2 * t1239 + t3017 - t3020;
2618 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2619 v4rho3sigma[9] = t7 * (t3360 + t3375 + t3384 + t3389) + 0.60968118750000000000e-3 * t1277 - 0.3e1 / 0.2e1 * t1279 - 0.3e1 * t1244 + t3002 - t3003 + t3005 - t3006 - t3009 - t3012 + t3403;
2620
2621 t3414 = -0.85355366250000000000e-2 * t3101 + 0.12193623750000000000e-2 * t3103 - 0.3e1 * t3107 + t3023 + 0.14e2 * t3111 - 0.3e1 * t3115 + 0.12193623750000000000e-2 * t3121 - 0.6e1 * t3132 + t3024 + 0.12193623750000000000e-2 * t3254 - 0.3e1 * t3259 + 0.9e1 / 0.2e1 * t3265;
2622 t3422 = -0.18290435625000000000e-2 * t3272 - t3029 + t3040 + t3041 - t3042 + 0.40645412500000000000e-3 * t3363 - t3365 - 0.15e2 / 0.4e1 * t3368 + 0.15242029687500000000e-2 * t3370 + 0.7e1 * t3282 - 0.21e2 / 0.2e1 * t3284 + 0.64016524687500000000e-2 * t3286;
2623 t3431 = 0.9e1 / 0.2e1 * t3288 - 0.3e1 * t3293 - 0.42677683125000000000e-2 * t3297 - 0.3e1 * t3157 + 0.19204957406250000000e-1 * t3162 - 0.70e2 / 0.3e1 * t3164 + 0.14e2 * t3168 + t3044 - t3046 + t3052 - t3053 + t3054;
2624 t3434 = t3055 - t3060 + t3061 + t3065 + t3066 - t3067 - t3069 - t3071 + t3072 - 0.18290435625000000000e-2 * t3385 + 0.9e1 / 0.2e1 * t3387 - t3073;
2625 t3448 = t3090 + t3091 + t3092 + 0.9e1 / 0.2e1 * t1282 - 0.18290435625000000000e-2 * t1284 + 0.24387247500000000000e-2 * t1253 + 0.14e2 * t1237 - 0.6e1 * t1257 - 0.85355366250000000000e-2 * t1239 + t3094 - t3097;
2626 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2627 v4rho3sigma[10] = t7 * (t3414 + t3422 + t3431 + t3434) + 0.12193623750000000000e-2 * t1277 - 0.3e1 * t1279 - 0.6e1 * t1244 + t3079 - t3080 + t3082 - t3083 - t3086 - t3089 + t3448;
2628
2629 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2630 v4rho3sigma[11] = v4rho3sigma[9];
2631
2632 t3449 = t1303 * t1199;
2633 t3450 = 0.60968118750000000000e-3 * t3449;
2634 t3452 = t547 * t370 * t228;
2635 t3453 = 0.60968118750000000000e-3 * t3452;
2636 t3454 = t1306 * t1195;
2637 t3459 = t551 * t330 * t86 * t184 * t332;
2638 t3460 = 0.33040991145903125000e-6 * t3459;
2639 t3461 = t1315 * t384;
2640 t3462 = 0.16520495572951562500e-6 * t3461;
2641 t3464 = t1314 * t181 * t374;
2642 t3465 = 0.33040991145903125000e-6 * t3464;
2643 t3466 = t1318 * t1199;
2644 t3467 = 0.16520495572951562500e-6 * t3466;
2645 t3469 = t552 * t370 * t228;
2646 t3470 = 0.16520495572951562500e-6 * t3469;
2647 t3471 = t1321 * t1195;
2648 t3473 = t1321 * t1152;
2649 t3477 = t552 * t367 * t228;
2650 t3478 = 0.16520495572951562500e-6 * t3477;
2651 t3479 = t1311 * t86;
2652 t3480 = t3479 * t512;
2653 t3483 = t1311 * t181 * t228;
2654 t3484 = 0.15419129201421458333e-5 * t3483;
2655 t3485 = t1315 * t381;
2656 t3486 = 0.16520495572951562500e-6 * t3485;
2657 t3488 = t547 * t367 * t228;
2658 t3489 = 0.60968118750000000000e-3 * t3488;
2659 t3490 = t1306 * t1152;
2660 t3495 = 0.1e1 / t65 / t1841 * t78 * t228;
2661 t3496 = 0.43687532737360798611e-5 * t3495;
2662 t3499 = t1310 * t155 * t86 * t188;
2663 t3500 = 0.15419129201421458333e-5 * t3499;
2664 t3502 = t1172 * t546 * t228;
2665 t3503 = 0.96024787031250000000e-2 * t3502;
2666 t3504 = t1300 * t86;
2667 t3505 = t3504 * t512;
2668 t3506 = 0.21338841562500000000e-2 * t3505;
2669 t3508 = t1300 * t181 * t228;
2670 t3509 = 0.42677683125000000000e-2 * t3508;
2671 t3510 = t3478 + 0.77095646007107291667e-6 * t3480 - t3484 - t3486 - t3489 + 0.30484059375000000000e-3 * t3490 + t3496 + t3500 - t3503 - t3506 + t3509;
2672 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2673 v4rho2sigma2[0] = t1326 - t1327 + t1328 - t1329 - t1330 + t1331 - t1332 + t7 * (t3450 - t3453 - 0.45726089062500000000e-3 * t3454 + t3460 + t3462 - t3465 - t3467 + t3470 + 0.12390371679713671875e-6 * t3471 - 0.82602477864757812500e-7 * t3473 + t3510);
2674
2675 t3513 = 0.24387247500000000000e-2 * t1304;
2676 t3514 = 0.30838258402842916666e-5 * t1312;
2677 t3515 = 0.66081982291806250000e-6 * t1316;
2678 t3516 = 0.66081982291806250000e-6 * t1319;
2679 t3517 = 0.12193623750000000000e-2 * t3449;
2680 t3518 = 0.12193623750000000000e-2 * t3452;
2681 t3520 = 0.66081982291806250000e-6 * t3459;
2682 t3521 = 0.33040991145903125000e-6 * t3461;
2683 t3522 = 0.66081982291806250000e-6 * t3464;
2684 t3523 = 0.33040991145903125000e-6 * t3466;
2685 t3524 = 0.33040991145903125000e-6 * t3469;
2686 t3528 = 0.33040991145903125000e-6 * t3477;
2687 t3530 = 0.30838258402842916666e-5 * t3483;
2688 t3531 = 0.33040991145903125000e-6 * t3485;
2689 t3532 = 0.12193623750000000000e-2 * t3488;
2690 t3534 = 0.87375065474721597220e-5 * t3495;
2691 t3535 = 0.30838258402842916666e-5 * t3499;
2692 t3536 = 0.19204957406250000000e-1 * t3502;
2693 t3537 = 0.42677683125000000000e-2 * t3505;
2694 t3538 = 0.85355366250000000000e-2 * t3508;
2695 t3539 = t3528 + 0.15419129201421458333e-5 * t3480 - t3530 - t3531 - t3532 + 0.60968118750000000000e-3 * t3490 + t3534 + t3535 - t3536 - t3537 + t3538;
2696 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2697 v4rho2sigma2[1] = t1335 - t3513 + t1337 - t3514 - t3515 + t3516 - t1341 + t7 * (t3517 - t3518 - 0.91452178125000000000e-3 * t3454 + t3520 + t3521 - t3522 - t3523 + t3524 + 0.24780743359427343750e-6 * t3471 - 0.16520495572951562500e-6 * t3473 + t3539);
2698
2699 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2700 v4rho2sigma2[2] = v4rho2sigma2[0];
2701
2702 t3542 = 0.17071073250000000000e-1 * t1301;
2703 t3543 = 0.487744950e-2 * t1304;
2704 t3545 = 0.61676516805685833334e-5 * t1312;
2705 t3546 = 0.13216396458361250e-5 * t1316;
2706 t3547 = 0.13216396458361250e-5 * t1319;
2707 t3550 = 0.243872475e-2 * t3452;
2708 t3552 = 0.13216396458361250e-5 * t3459;
2709 t3554 = 0.13216396458361250e-5 * t3464;
2710 t3556 = 0.6608198229180625e-6 * t3469;
2711 t3560 = 0.6608198229180625e-6 * t3477;
2712 t3562 = 0.61676516805685833334e-5 * t3483;
2713 t3563 = 0.6608198229180625e-6 * t3485;
2714 t3564 = 0.243872475e-2 * t3488;
2715 t3566 = 0.17475013094944319445e-4 * t3495;
2716 t3567 = 0.61676516805685833334e-5 * t3499;
2717 t3568 = 0.38409914812500000000e-1 * t3502;
2718 t3570 = 0.17071073250000000000e-1 * t3508;
2719 t3571 = t3560 + 0.30838258402842916667e-5 * t3480 - t3562 - t3563 - t3564 + 0.12193623750000000000e-2 * t3490 + t3566 + t3567 - t3568 - 0.85355366250000000000e-2 * t3505 + t3570;
2720 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2721 v4rho2sigma2[3] = t3542 - t3543 + 0.24387247500000000000e-2 * t1307 - t3545 - t3546 + t3547 - 0.66081982291806250000e-6 * t1322 + t7 * (0.24387247500000000000e-2 * t3449 - t3550 - 0.18290435625000000000e-2 * t3454 + t3552 + 0.66081982291806250000e-6 * t3461 - t3554 - 0.66081982291806250000e-6 * t3466 + t3556 + 0.49561486718854687500e-6 * t3471 - 0.33040991145903125000e-6 * t3473 + t3571);
2722
2723 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2724 v4rho2sigma2[4] = v4rho2sigma2[1];
2725
2726 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2727 v4rho2sigma2[5] = v4rho2sigma2[2];
2728
2729 t3574 = t3504 * t536;
2730 t3576 = t3479 * t536;
2731 t3578 = t1321 * t1241;
2732 t3580 = t1306 * t1241;
2733 t3583 = t3496 - t3503 - 0.10669420781250000000e-2 * t3574 + 0.38547823003553645833e-6 * t3576 - 0.82602477864757812500e-7 * t3578 + 0.30484059375000000000e-3 * t3580 - t3489 + t3500 - 0.10669420781250000000e-2 * t3505 + t3509 - t3453 + t3460 + t3470;
2734 t3585 = t1306 * t442;
2735 t3587 = t1303 * t1252;
2736 t3589 = t1315 * t436;
2737 t3591 = t1321 * t442;
2738 t3593 = t1318 * t1252;
2739 t3598 = t3478 + 0.38547823003553645834e-6 * t3480 - t3484 - t3486 - 0.45726089062500000000e-3 * t3585 + 0.30484059375000000000e-3 * t3587 + 0.82602477864757812500e-7 * t3589 + 0.12390371679713671875e-6 * t3591 - 0.82602477864757812500e-7 * t3593 + 0.30484059375000000000e-3 * t3449 + 0.82602477864757812500e-7 * t3461 - t3465 - 0.82602477864757812500e-7 * t3466;
2740 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2741 v4rho2sigma2[6] = t1326 - t1327 + t1308 - t1329 - t1330 + t1331 - t1323 + t1345 - t1347 + t7 * (t3583 + t3598);
2742
2743 t3601 = 0.21338841562500000000e-2 * t3574;
2744 t3602 = 0.77095646007107291667e-6 * t3576;
2745 t3605 = t3534 - t3536 - t3601 + t3602 - 0.16520495572951562500e-6 * t3578 + 0.60968118750000000000e-3 * t3580 - t3532 + t3535 - t3506 + t3538 - t3518 + t3520 + t3524;
2746 t3608 = 0.60968118750000000000e-3 * t3587;
2747 t3609 = 0.16520495572951562500e-6 * t3589;
2748 t3611 = 0.16520495572951562500e-6 * t3593;
2749 t3612 = t3528 + 0.77095646007107291665e-6 * t3480 - t3530 - t3531 - 0.91452178125000000000e-3 * t3585 + t3608 + t3609 + 0.24780743359427343750e-6 * t3591 - t3611 + t3450 + t3462 - t3522 - t3467;
2750 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2751 v4rho2sigma2[7] = t1335 - t3513 + t1328 - t3514 - t3515 + t3516 - t1332 + t1350 - t1351 + t7 * (t3605 + t3612);
2752
2753 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2754 v4rho2sigma2[8] = v4rho2sigma2[6];
2755
2756 t3615 = 0.42677683125000000000e-2 * t3574;
2757 t3616 = 0.15419129201421458333e-5 * t3576;
2758 t3619 = t3566 - t3568 - t3615 + t3616 - 0.33040991145903125000e-6 * t3578 + 0.12193623750000000000e-2 * t3580 - t3564 + t3567 - t3537 + t3570 - t3550 + t3552 + t3556;
2759 t3622 = 0.12193623750000000000e-2 * t3587;
2760 t3623 = 0.33040991145903125000e-6 * t3589;
2761 t3625 = 0.33040991145903125000e-6 * t3593;
2762 t3626 = t3560 + 0.15419129201421458334e-5 * t3480 - t3562 - t3563 - 0.18290435625000000000e-2 * t3585 + t3622 + t3623 + 0.49561486718854687500e-6 * t3591 - t3625 + t3517 + t3521 - t3554 - t3523;
2763 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2764 v4rho2sigma2[9] = t3542 - t3543 + t1337 - t3545 - t3546 + t3547 - t1341 + t1354 - t1355 + t7 * (t3619 + t3626);
2765
2766 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2767 v4rho2sigma2[10] = v4rho2sigma2[7];
2768
2769 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2770 v4rho2sigma2[11] = v4rho2sigma2[8];
2771
2772 t3630 = t1321 * t1276;
2773 t3632 = t1306 * t1276;
2774 t3634 = t1306 * t1281;
2775 t3636 = t1321 * t1281;
2776 t3638 = -t3484 - t3486 - 0.82602477864757812500e-7 * t3630 + 0.30484059375000000000e-3 * t3632 - 0.45726089062500000000e-3 * t3634 + 0.12390371679713671875e-6 * t3636 - t3489 + t3496 + t3500 - t3503 + t3509;
2777 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2778 v4rho2sigma2[12] = t1326 - t1327 + t1350 - t1329 - t1330 + t1331 - t1351 + t7 * (-t3601 + t3608 + t3602 + t3609 - t3611 - t3453 + t3460 - t3465 + t3470 + t3478 + t3638);
2779
2780 t3646 = -t3530 - t3531 - 0.16520495572951562500e-6 * t3630 + 0.60968118750000000000e-3 * t3632 - 0.91452178125000000000e-3 * t3634 + 0.24780743359427343750e-6 * t3636 - t3532 + t3534 + t3535 - t3536 + t3538;
2781 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2782 v4rho2sigma2[13] = t1335 - t3513 + t1354 - t3514 - t3515 + t3516 - t1355 + t7 * (-t3615 + t3622 + t3616 + t3623 - t3625 - t3518 + t3520 - t3522 + t3524 + t3528 + t3646);
2783
2784 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2785 v4rho2sigma2[14] = v4rho2sigma2[12];
2786
2787 t3661 = -t3562 - t3563 - 0.33040991145903125000e-6 * t3630 + 0.12193623750000000000e-2 * t3632 - 0.18290435625000000000e-2 * t3634 + 0.49561486718854687500e-6 * t3636 - t3564 + t3566 + t3567 - t3568 + t3570;
2788 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2789 v4rho2sigma2[15] = t3542 - t3543 + 0.24387247500000000000e-2 * t1344 - t3545 - t3546 + t3547 - 0.66081982291806250000e-6 * t1346 + t7 * (-0.85355366250000000000e-2 * t3574 + 0.24387247500000000000e-2 * t3587 + 0.30838258402842916667e-5 * t3576 + 0.66081982291806250000e-6 * t3589 - 0.66081982291806250000e-6 * t3593 - t3550 + t3552 - t3554 + t3556 + t3560 + t3661);
2790
2791 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2792 v4rho2sigma2[16] = v4rho2sigma2[13];
2793
2794 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2795 v4rho2sigma2[17] = v4rho2sigma2[14];
2796
2797 t3665 = t516 * t1359 * t228;
2798 t3666 = 0.10669420781250000000e-2 * t3665;
2799 t3668 = t1360 * t181 * t228;
2800 t3669 = 0.30484059375000000000e-3 * t3668;
2801 t3670 = t1360 * t86;
2802 t3671 = t3670 * t512;
2803 t3675 = t1310 * t1363 * t78 * t228;
2804 t3676 = 0.11564346901066093750e-5 * t3675;
2805 t3678 = t1364 * t155 * t374;
2806 t3679 = 0.24780743359427343750e-6 * t3678;
2807 t3680 = t1365 * t2885;
2808 t3681 = 0.24780743359427343750e-6 * t3680;
2809 t3682 = t1365 * t1199;
2810 t3688 = 0.1e1 / t1370 / t1841 * t155 * t546 * t228;
2811 t3689 = 0.39169804173910672519e-9 * t3688;
2812 t3692 = t1372 * t330 * t546 * t374;
2813 t3693 = 0.13429647145340802007e-9 * t3692;
2814 t3694 = t1374 * t2885;
2815 t3695 = 0.67148235726704010033e-10 * t3694;
2816 t3696 = t1374 * t1199;
2817 t3698 = -t3666 + t3669 - 0.15242029687500000000e-3 * t3671 - t3676 - t3679 + t3681 - 0.12390371679713671875e-6 * t3682 + t3689 + t3693 - t3695 + 0.33574117863352005016e-10 * t3696;
2818 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2819 v4rhosigma3[0] = t7 * t3698 + t1362 + t1367 - t1376;
2820
2821 t3700 = 0.21338841562500000000e-2 * t3665;
2822 t3701 = 0.60968118750000000000e-3 * t3668;
2823 t3703 = 0.23128693802132187500e-5 * t3675;
2824 t3704 = 0.49561486718854687500e-6 * t3678;
2825 t3705 = 0.49561486718854687500e-6 * t3680;
2826 t3707 = 0.78339608347821345041e-9 * t3688;
2827 t3708 = 0.26859294290681604014e-9 * t3692;
2828 t3709 = 0.13429647145340802007e-9 * t3694;
2829 t3711 = -t3700 + t3701 - 0.30484059375000000000e-3 * t3671 - t3703 - t3704 + t3705 - 0.24780743359427343750e-6 * t3682 + t3707 + t3708 - t3709 + 0.67148235726704010035e-10 * t3696;
2830 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2831 v4rhosigma3[1] = t7 * t3711 + t1378 + t1379 - t1380;
2832
2833 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2834 v4rhosigma3[2] = v4rhosigma3[0];
2835
2836 t3713 = 0.42677683125000000000e-2 * t3665;
2837 t3714 = 0.12193623750000000000e-2 * t3668;
2838 t3716 = 0.46257387604264375000e-5 * t3675;
2839 t3717 = 0.99122973437709375000e-6 * t3678;
2840 t3718 = 0.99122973437709375000e-6 * t3680;
2841 t3720 = 0.15667921669564269008e-8 * t3688;
2842 t3721 = 0.53718588581363208026e-9 * t3692;
2843 t3722 = 0.26859294290681604013e-9 * t3694;
2844 t3724 = -t3713 + t3714 - 0.60968118750000000000e-3 * t3671 - t3716 - t3717 + t3718 - 0.49561486718854687500e-6 * t3682 + t3720 + t3721 - t3722 + 0.13429647145340802006e-9 * t3696;
2845 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2846 v4rhosigma3[3] = t7 * t3724 + t1382 + t1383 - t1384;
2847
2848 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2849 v4rhosigma3[4] = v4rhosigma3[1];
2850
2851 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2852 v4rhosigma3[5] = v4rhosigma3[2];
2853
2854 t3726 = 0.85355366250000000000e-2 * t3665;
2855 t3727 = 0.243872475e-2 * t3668;
2856 t3729 = 0.92514775208528750000e-5 * t3675;
2857 t3730 = 0.19824594687541875e-5 * t3678;
2858 t3731 = 0.19824594687541875e-5 * t3680;
2859 t3733 = 0.31335843339128538016e-8 * t3688;
2860 t3734 = 0.10743717716272641605e-8 * t3692;
2861 t3735 = 0.53718588581363208027e-9 * t3694;
2862 t3737 = -t3726 + t3727 - 0.12193623750000000000e-2 * t3671 - t3729 - t3730 + t3731 - 0.99122973437709375000e-6 * t3682 + t3733 + t3734 - t3735 + 0.26859294290681604014e-9 * t3696;
2863 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2864 v4rhosigma3[6] = t7 * t3737 + t1386 + t1387 - t1388;
2865
2866 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2867 v4rhosigma3[7] = v4rhosigma3[3];
2868
2869 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2870 v4rhosigma3[8] = v4rhosigma3[4];
2871
2872 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2873 v4rhosigma3[9] = v4rhosigma3[5];
2874
2875 t3739 = t3670 * t536;
2876 t3741 = t1365 * t1252;
2877 t3743 = t1374 * t1252;
2878 t3745 = -t3666 + t3669 - 0.15242029687500000000e-3 * t3739 - t3676 - t3679 + t3681 - 0.12390371679713671875e-6 * t3741 + t3689 + t3693 - t3695 + 0.33574117863352005016e-10 * t3743;
2879 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2880 v4rhosigma3[10] = t7 * t3745 + t1362 + t1367 - t1376;
2881
2882 t3750 = -t3700 + t3701 - 0.30484059375000000000e-3 * t3739 - t3703 - t3704 + t3705 - 0.24780743359427343750e-6 * t3741 + t3707 + t3708 - t3709 + 0.67148235726704010035e-10 * t3743;
2883 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2884 v4rhosigma3[11] = t7 * t3750 + t1378 + t1379 - t1380;
2885
2886 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2887 v4rhosigma3[12] = v4rhosigma3[10];
2888
2889 t3755 = -t3713 + t3714 - 0.60968118750000000000e-3 * t3739 - t3716 - t3717 + t3718 - 0.49561486718854687500e-6 * t3741 + t3720 + t3721 - t3722 + 0.13429647145340802006e-9 * t3743;
2890 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2891 v4rhosigma3[13] = t7 * t3755 + t1382 + t1383 - t1384;
2892
2893 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2894 v4rhosigma3[14] = v4rhosigma3[11];
2895
2896 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2897 v4rhosigma3[15] = v4rhosigma3[12];
2898
2899 t3760 = -t3726 + t3727 - 0.12193623750000000000e-2 * t3739 - t3729 - t3730 + t3731 - 0.99122973437709375000e-6 * t3741 + t3733 + t3734 - t3735 + 0.26859294290681604014e-9 * t3743;
2900 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2901 v4rhosigma3[16] = t7 * t3760 + t1386 + t1387 - t1388;
2902
2903 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2904 v4rhosigma3[17] = v4rhosigma3[13];
2905
2906 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2907 v4rhosigma3[18] = v4rhosigma3[14];
2908
2909 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2910 v4rhosigma3[19] = v4rhosigma3[15];
2911
2912 t3762 = t56 * t56;
2913 t3766 = t225 / t79 / t3762 * t228;
2914 t3771 = t551 / t3762 * t78 * t228;
2915 t3775 = t1372 * t1359 * t155 * t228;
2916 t3781 = 0.1e1 / t315 / t147 * t330 * t1363 * t228;
2917 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2918 v4sigma4[0] = t7 * (-0.45726089062500000000e-3 * t3766 - 0.37171115039141015625e-6 * t3771 - 0.67148235726704010034e-10 * t3775 + 0.27292677397591217532e-13 * t3781);
2919
2920 t3784 = 0.91452178125000000000e-3 * t3766;
2921 t3785 = 0.74342230078282031250e-6 * t3771;
2922 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2923 v4sigma4[1] = t7 * (-t3784 - t3785 - 0.13429647145340802006e-9 * t3775 + 0.54585354795182435066e-13 * t3781);
2924
2925 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2926 v4sigma4[2] = v4sigma4[0];
2927
2928 t3789 = 0.18290435625000000000e-2 * t3766;
2929 t3790 = 0.14868446015656406250e-5 * t3771;
2930 t3792 = 0.10917070959036487013e-12 * t3781;
2931 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2932 v4sigma4[3] = t7 * (-t3789 - t3790 - 0.26859294290681604014e-9 * t3775 + t3792);
2933
2934 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2935 v4sigma4[4] = v4sigma4[1];
2936
2937 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2938 v4sigma4[5] = v4sigma4[2];
2939
2940 t3794 = 0.36580871250000000000e-2 * t3766;
2941 t3795 = 0.29736892031312812500e-5 * t3771;
2942 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2943 v4sigma4[6] = t7 * (-t3794 - t3795 - 0.53718588581363208026e-9 * t3775 + 0.21834141918072974026e-12 * t3781);
2944
2945 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2946 v4sigma4[7] = v4sigma4[3];
2947
2948 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2949 v4sigma4[8] = v4sigma4[4];
2950
2951 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2952 v4sigma4[9] = v4sigma4[5];
2953
2954 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2955 v4sigma4[10] = t7 * (-0.731617425e-2 * t3766 - 0.59473784062625625e-5 * t3771 - 0.10743717716272641605e-8 * t3775 + 0.43668283836145948052e-12 * t3781);
2956
2957 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2958 v4sigma4[11] = t7 * (-t3794 - t3795 - 0.53718588581363208027e-9 * t3775 + 0.21834141918072974025e-12 * t3781);
2959
2960 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2961 v4sigma4[12] = t7 * (-t3789 - t3790 - 0.26859294290681604013e-9 * t3775 + t3792);
2962
2963 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2964 v4sigma4[13] = t7 * (-t3784 - t3785 - 0.13429647145340802007e-9 * t3775 + 0.54585354795182435064e-13 * t3781);
2965
2966 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2967 v4sigma4[14] = v4sigma4[9];
2968
2969 #ifndef XC_DONT_COMPILE_MXC
2970
2971 if(order < 5) return;
2972
2973
2974 #endif
2975
2976 #endif
2977
2978 #endif
2979
2980 #endif
2981
2982 #endif
2983
2984
2985 }
2986
2987