1 /*
2 This file was generated automatically with ./scripts/maple2c.pl.
3 Do not edit this file directly as it can be overwritten!!
4
5 This Source Code Form is subject to the terms of the Mozilla Public
6 License, v. 2.0. If a copy of the MPL was not distributed with this
7 file, You can obtain one at http://mozilla.org/MPL/2.0/.
8
9 Maple version : Maple 2016 (X86 64 LINUX)
10 Maple source : ./maple/gga_exc/gga_x_cap.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t3, t5, t6, t7, t8, t9, t10;
24 double t11, t12, t13, t14, t15, t16, t17, t18;
25 double t20, t21, t22, t24, t29, t30, t32, t33;
26 double t34, t35, t38, t40;
27
28 #ifndef XC_DONT_COMPILE_VXC
29 double t44, t45, t46, t50, t53, t54, t55, t56;
30 double t57, t58, t59, t61, t63, t64, t65, t69;
31 double t71, t72, t73, t74, t75, t78, t82, t84;
32 double t87, t89, t90, t91, t94, t95, t97, t101;
33
34 #ifndef XC_DONT_COMPILE_FXC
35 double t104, t105, t115, t118, t120, t122, t126, t130;
36 double t131, t132, t133, t134, t135, t136, t137, t138;
37 double t141, t145, t147, t149, t151, t156, t159, t164;
38 double t173, t174, t175, t176, t177, t180, t182, t186;
39 double t187, t189, t195, t198, t200, t203, t204, t205;
40 double t209, t212, t213, t218, t222, t230;
41
42 #ifndef XC_DONT_COMPILE_KXC
43 double t234, t235, t248, t252, t254, t258, t261, t262;
44 double t263, t271, t276, t277, t279, t280, t283, t285;
45 double t286, t290, t291, t294, t298, t299, t300, t301;
46 double t302, t304, t308, t312, t315, t320, t331, t337;
47 double t343, t344, t345, t347, t353, t359, t362, t364;
48 double t365, t367, t371, t372, t374, t378, t379, t381;
49 double t384, t385, t396, t397, t398, t401, t405, t411;
50 double t412, t413, t416, t423, t425, t426, t430, t431;
51 double t435, t436, t439, t440, t443, t444, t446, t449;
52 double t450, t451, t455, t459, t460, t461, t463, t467;
53 double t469, t473, t475, t478, t480, t482, t486, t488;
54 double t492, t494, t498;
55
56 #ifndef XC_DONT_COMPILE_LXC
57 double t516, t519, t526, t528, t529, t531, t532, t536;
58 double t537, t545, t546, t551, t557, t559, t571, t579;
59 double t581, t582, t605, t610, t627, t632, t633, t641;
60 double t673, t686, t689, t696, t703, t710, t716, t721;
61 double t743, t744, t755, t765, t794, t812, t819, t835;
62 double t836, t864, t892, t904, t909, t911, t920, t925;
63 double t926, t950, t1012;
64 #endif
65
66 #endif
67
68 #endif
69
70 #endif
71
72 #endif
73
74
75 gga_x_cap_params *params;
76
77 assert(p->params != NULL);
78 params = (gga_x_cap_params * )(p->params);
79
80 t1 = M_CBRT3;
81 t3 = POW_1_3(0.1e1 / M_PI);
82 t5 = M_CBRT4;
83 t6 = t5 * t5;
84 t7 = t1 * t3 * t6;
85 t8 = M_CBRT2;
86 t9 = t8 * t8;
87 t10 = POW_1_3(rho[0]);
88 t11 = t9 * t10;
89 t12 = M_CBRT6;
90 t13 = t12 * t12;
91 t14 = params->alphaoAx * t13;
92 t15 = M_PI * M_PI;
93 t16 = POW_1_3(t15);
94 t17 = 0.1e1 / t16;
95 t18 = sqrt(sigma[0]);
96 t20 = t14 * t17 * t18;
97 t21 = t10 * rho[0];
98 t22 = 0.1e1 / t21;
99 t24 = t13 * t17;
100 t29 = 0.1e1 + t24 * t18 * t8 * t22 / 0.12e2;
101 t30 = log(t29);
102 t32 = params->c * t30 + 0.1e1;
103 t33 = 0.1e1 / t32;
104 t34 = t30 * t33;
105 t35 = t8 * t22 * t34;
106 t38 = 0.1e1 - t20 * t35 / 0.12e2;
107 t40 = t7 * t11 * t38;
108 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
109 zk[0] = -0.3e1 / 0.16e2 * t40;
110
111 #ifndef XC_DONT_COMPILE_VXC
112
113 if(order < 1) return;
114
115
116 t44 = t21 * t1 * t3;
117 t45 = t6 * t9;
118 t46 = rho[0] * rho[0];
119 t50 = t8 / t10 / t46 * t34;
120 t53 = params->alphaoAx * t12;
121 t54 = t16 * t16;
122 t55 = 0.1e1 / t54;
123 t56 = t55 * sigma[0];
124 t57 = t53 * t56;
125 t58 = t46 * rho[0];
126 t59 = t10 * t10;
127 t61 = 0.1e1 / t59 / t58;
128 t63 = 0.1e1 / t29;
129 t64 = t63 * t33;
130 t65 = t9 * t61 * t64;
131 t69 = t53 * t56 * t9;
132 t71 = t32 * t32;
133 t72 = 0.1e1 / t71;
134 t73 = t72 * params->c;
135 t74 = t73 * t63;
136 t75 = t61 * t30 * t74;
137 t78 = t20 * t50 / 0.9e1 + t57 * t65 / 0.18e2 - t69 * t75 / 0.18e2;
138 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
139 vrho[0] = -t40 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t45 * t78;
140
141 t82 = 0.1e1 / t18;
142 t84 = t14 * t17 * t82;
143 t87 = t53 * t55;
144 t89 = 0.1e1 / t59 / t46;
145 t90 = t9 * t89;
146 t91 = t90 * t64;
147 t94 = t55 * t9;
148 t95 = t53 * t94;
149 t97 = t89 * t30 * t74;
150 t101 = t45 * (-t84 * t35 / 0.24e2 - t87 * t91 / 0.48e2 + t95 * t97 / 0.48e2);
151 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
152 vsigma[0] = -0.3e1 / 0.16e2 * t44 * t101;
153
154 #ifndef XC_DONT_COMPILE_FXC
155
156 if(order < 2) return;
157
158
159 t104 = 0.1e1 / t59;
160 t105 = t9 * t104;
161 t115 = t8 / t10 / t58 * t34;
162 t118 = t46 * t46;
163 t120 = 0.1e1 / t59 / t118;
164 t122 = t9 * t120 * t64;
165 t126 = t120 * t30 * t74;
166 t130 = params->alphaoAx / t15;
167 t131 = t18 * sigma[0];
168 t132 = t130 * t131;
169 t133 = t118 * t46;
170 t134 = 0.1e1 / t133;
171 t135 = t29 * t29;
172 t136 = 0.1e1 / t135;
173 t137 = t134 * t136;
174 t138 = t137 * t33;
175 t141 = t137 * t73;
176 t145 = t130 * t131 * t134;
177 t147 = 0.1e1 / t71 / t32;
178 t149 = params->c * params->c;
179 t151 = t30 * t147 * t149 * t136;
180 t156 = t30 * t72 * params->c * t136;
181 t159 = -0.7e1 / 0.27e2 * t20 * t115 - 0.5e1 / 0.18e2 * t57 * t122 + 0.5e1 / 0.18e2 * t69 * t126 + 0.2e1 / 0.27e2 * t132 * t138 + 0.4e1 / 0.27e2 * t132 * t141 - 0.4e1 / 0.27e2 * t145 * t151 - 0.2e1 / 0.27e2 * t145 * t156;
182 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
183 v2rho2[0] = -t7 * t105 * t38 / 0.12e2 - t7 * t11 * t78 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t45 * t159;
184
185 t164 = t10 * t1 * t3;
186 t173 = t118 * rho[0];
187 t174 = 0.1e1 / t173;
188 t175 = t130 * t174;
189 t176 = t136 * t33;
190 t177 = t176 * t18;
191 t180 = t136 * t72;
192 t182 = t180 * params->c * t18;
193 t186 = t130 * t174 * t30;
194 t187 = t147 * t149;
195 t189 = t187 * t136 * t18;
196 t195 = t45 * (t84 * t50 / 0.18e2 + t87 * t65 / 0.12e2 - t95 * t75 / 0.12e2 - t175 * t177 / 0.36e2 - t175 * t182 / 0.18e2 + t186 * t189 / 0.18e2 + t186 * t182 / 0.36e2);
197 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
198 v2rhosigma[0] = -t164 * t101 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t195;
199
200 t198 = 0.1e1 / t131;
201 t200 = t14 * t17 * t198;
202 t203 = 0.1e1 / sigma[0];
203 t204 = t55 * t203;
204 t205 = t53 * t204;
205 t209 = t53 * t204 * t9;
206 t212 = 0.1e1 / t118;
207 t213 = t130 * t212;
208 t218 = t180 * params->c * t82;
209 t222 = t130 * t212 * t30;
210 t230 = t45 * (t200 * t35 / 0.48e2 - t205 * t91 / 0.96e2 + t209 * t97 / 0.96e2 + t213 * t176 * t82 / 0.96e2 + t213 * t218 / 0.48e2 - t222 * t187 * t136 * t82 / 0.48e2 - t222 * t218 / 0.96e2);
211 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
212 v2sigma2[0] = -0.3e1 / 0.16e2 * t44 * t230;
213
214 #ifndef XC_DONT_COMPILE_KXC
215
216 if(order < 3) return;
217
218
219 t234 = 0.1e1 / t59 / rho[0];
220 t235 = t9 * t234;
221 t248 = t8 / t10 / t118 * t34;
222 t252 = 0.1e1 / t59 / t173;
223 t254 = t9 * t252 * t64;
224 t258 = t252 * t30 * t74;
225 t261 = t118 * t58;
226 t262 = 0.1e1 / t261;
227 t263 = t262 * t136;
228 t271 = t130 * t131 * t262;
229 t276 = sigma[0] * sigma[0];
230 t277 = t118 * t118;
231 t279 = 0.1e1 / t10 / t277;
232 t280 = t276 * t279;
233 t283 = 0.1e1 / t135 / t29;
234 t285 = t24 * t8;
235 t286 = t283 * t33 * t285;
236 t290 = t130 * t280 * t283;
237 t291 = t73 * t285;
238 t294 = t187 * t285;
239 t298 = t130 * t280 * t30;
240 t299 = t71 * t71;
241 t300 = 0.1e1 / t299;
242 t301 = t149 * params->c;
243 t302 = t300 * t301;
244 t304 = t302 * t283 * t285;
245 t308 = t187 * t283 * t285;
246 t312 = t73 * t283 * t285;
247 t315 = 0.70e2 / 0.81e2 * t20 * t248 + 0.119e3 / 0.81e2 * t57 * t254 - 0.119e3 / 0.81e2 * t69 * t258 - 0.22e2 / 0.27e2 * t132 * t263 * t33 - 0.44e2 / 0.27e2 * t132 * t263 * t73 + 0.44e2 / 0.27e2 * t271 * t151 + 0.22e2 / 0.27e2 * t271 * t156 + 0.4e1 / 0.243e3 * t130 * t280 * t286 + 0.4e1 / 0.81e2 * t290 * t291 + 0.4e1 / 0.81e2 * t290 * t294 - 0.4e1 / 0.81e2 * t298 * t304 - 0.4e1 / 0.81e2 * t298 * t308 - 0.4e1 / 0.243e3 * t298 * t312;
248 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
249 v3rho3[0] = t7 * t235 * t38 / 0.18e2 - t7 * t105 * t78 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t159 - 0.3e1 / 0.16e2 * t44 * t45 * t315;
250
251 t320 = t104 * t1 * t3;
252 t331 = t130 * t134;
253 t337 = t130 * t134 * t30;
254 t343 = 0.1e1 / t10 / t261;
255 t344 = t343 * t283;
256 t345 = t130 * t344;
257 t347 = t33 * sigma[0] * t285;
258 t353 = sigma[0] * params->c * t285;
259 t359 = t149 * sigma[0] * t285;
260 t362 = t343 * t30;
261 t364 = t130 * t362 * t300;
262 t365 = t301 * t283;
263 t367 = t365 * sigma[0] * t285;
264 t371 = t130 * t362 * t147;
265 t372 = t149 * t283;
266 t374 = t372 * sigma[0] * t285;
267 t378 = t130 * t362 * t72;
268 t379 = params->c * t283;
269 t381 = t379 * sigma[0] * t285;
270 t384 = -0.7e1 / 0.54e2 * t84 * t115 - 0.37e2 / 0.108e3 * t87 * t122 + 0.37e2 / 0.108e3 * t95 * t126 + t331 * t177 / 0.4e1 + t331 * t182 / 0.2e1 - t337 * t189 / 0.2e1 - t337 * t182 / 0.4e1 - t345 * t347 / 0.162e3 - t130 * t344 * t72 * t353 / 0.54e2 - t130 * t344 * t147 * t359 / 0.54e2 + t364 * t367 / 0.54e2 + t371 * t374 / 0.54e2 + t378 * t381 / 0.162e3;
271 t385 = t45 * t384;
272 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
273 v3rho2sigma[0] = -t320 * t101 / 0.12e2 - t164 * t195 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t385;
274
275 t396 = t130 * t82;
276 t397 = t174 * t136;
277 t398 = t397 * t33;
278 t401 = t397 * t73;
279 t405 = t130 * t82 * t174;
280 t411 = 0.1e1 / t10 / t133;
281 t412 = t411 * t283;
282 t413 = t130 * t412;
283 t416 = t33 * t13 * t17 * t8;
284 t423 = t411 * t30;
285 t425 = t130 * t423 * t300;
286 t426 = t365 * t285;
287 t430 = t130 * t423 * t147;
288 t431 = t372 * t285;
289 t435 = t130 * t423 * t72;
290 t436 = t379 * t285;
291 t439 = -t200 * t50 / 0.36e2 + t205 * t65 / 0.72e2 - t209 * t75 / 0.72e2 - t396 * t398 / 0.18e2 - t396 * t401 / 0.9e1 + t405 * t151 / 0.9e1 + t405 * t156 / 0.18e2 + t413 * t416 / 0.432e3 + t413 * t291 / 0.144e3 + t413 * t294 / 0.144e3 - t425 * t426 / 0.144e3 - t430 * t431 / 0.144e3 - t435 * t436 / 0.432e3;
292 t440 = t45 * t439;
293 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
294 v3rhosigma2[0] = -t164 * t230 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t440;
295
296 t443 = t18 * t276;
297 t444 = 0.1e1 / t443;
298 t446 = t14 * t17 * t444;
299 t449 = 0.1e1 / t276;
300 t450 = t55 * t449;
301 t451 = t53 * t450;
302 t455 = t53 * t450 * t9;
303 t459 = 0.1e1 / t10 / t173;
304 t460 = t459 * t283;
305 t461 = t130 * t460;
306 t463 = t33 * t203 * t285;
307 t467 = t130 * t460 * t72;
308 t469 = t203 * params->c * t285;
309 t473 = t130 * t460 * t147;
310 t475 = t149 * t203 * t285;
311 t478 = t459 * t30;
312 t480 = t130 * t478 * t300;
313 t482 = t365 * t203 * t285;
314 t486 = t130 * t478 * t147;
315 t488 = t372 * t203 * t285;
316 t492 = t130 * t478 * t72;
317 t494 = t379 * t203 * t285;
318 t498 = t45 * (-t446 * t35 / 0.32e2 + t451 * t91 / 0.64e2 - t455 * t97 / 0.64e2 - t461 * t463 / 0.1152e4 - t467 * t469 / 0.384e3 - t473 * t475 / 0.384e3 + t480 * t482 / 0.384e3 + t486 * t488 / 0.384e3 + t492 * t494 / 0.1152e4);
319 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
320 v3sigma3[0] = -0.3e1 / 0.16e2 * t44 * t498;
321
322 #ifndef XC_DONT_COMPILE_LXC
323
324 if(order < 4) return;
325
326
327 t516 = t277 * rho[0];
328 t519 = t276 / t10 / t516;
329 t526 = t443 / t59 / t277 / t46;
330 t528 = t135 * t135;
331 t529 = 0.1e1 / t528;
332 t531 = t12 * t55;
333 t532 = t531 * t9;
334 t536 = 0.1e1 / t277;
335 t537 = t536 * t136;
336 t545 = 0.1e1 / t59 / t133;
337 t546 = t545 * t30;
338 t551 = t130 * t519 * t283;
339 t557 = t130 * t526 * t529;
340 t559 = t94 * params->c;
341 t571 = t130 * t519 * t30;
342 t579 = t130 * t526 * t30;
343 t581 = 0.1e1 / t299 / t32;
344 t582 = t149 * t149;
345 t605 = t130 * t131 * t536;
346 t610 = 0.32e2 / 0.243e3 * t557 * t302 * t532 + 0.232e3 / 0.243e3 * t571 * t304 + 0.232e3 / 0.243e3 * t571 * t308 + 0.232e3 / 0.729e3 * t571 * t312 - 0.32e2 / 0.243e3 * t579 * t581 * t582 * t529 * t532 - 0.16e2 / 0.81e2 * t579 * t302 * t529 * t532 - 0.88e2 / 0.729e3 * t579 * t187 * t529 * t532 - 0.8e1 / 0.243e3 * t579 * t73 * t529 * t532 - 0.721e3 / 0.81e2 * t57 * t9 * t545 * t64 - 0.3724e4 / 0.243e3 * t605 * t151 - 0.1862e4 / 0.243e3 * t605 * t156;
347 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
348 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t90 * t38 + 0.2e1 / 0.9e1 * t7 * t235 * t78 - t7 * t105 * t159 / 0.2e1 - t7 * t11 * t315 - 0.3e1 / 0.16e2 * t44 * t45 * (-0.910e3 / 0.243e3 * t20 * t8 * t459 * t34 - 0.232e3 / 0.729e3 * t130 * t519 * t286 + 0.8e1 / 0.243e3 * t130 * t526 * t529 * t33 * t532 + 0.3724e4 / 0.243e3 * t132 * t537 * t73 + 0.1862e4 / 0.243e3 * t132 * t537 * t33 + 0.721e3 / 0.81e2 * t69 * t546 * t74 - 0.232e3 / 0.243e3 * t551 * t291 - 0.232e3 / 0.243e3 * t551 * t294 + 0.88e2 / 0.729e3 * t557 * t72 * t12 * t559 + 0.16e2 / 0.81e2 * t557 * t187 * t532 + t610);
349
350 t627 = t279 * t283;
351 t632 = 0.1e1 / t59 / t516;
352 t633 = t632 * t529;
353 t641 = t130 * t262;
354 t673 = t279 * t30;
355 t686 = t632 * t30;
356 t689 = t582 * t529;
357 t696 = t301 * t529;
358 t703 = t149 * t529;
359 t710 = params->c * t529;
360 t716 = t130 * t262 * t30;
361 t721 = -0.2e1 / 0.27e2 * t130 * t633 * t147 * t131 * t149 * t532 - 0.4e1 / 0.81e2 * t130 * t633 * t300 * t301 * t131 * t532 - 0.49e2 / 0.162e3 * t130 * t673 * t300 * t367 - 0.49e2 / 0.162e3 * t130 * t673 * t147 * t374 - 0.49e2 / 0.486e3 * t130 * t673 * t72 * t381 + 0.4e1 / 0.81e2 * t130 * t686 * t581 * t689 * t131 * t532 + 0.2e1 / 0.27e2 * t130 * t686 * t300 * t696 * t131 * t532 + 0.11e2 / 0.243e3 * t130 * t686 * t147 * t703 * t131 * t532 + t130 * t686 * t72 * t710 * t131 * t532 / 0.81e2 + 0.317e3 / 0.81e2 * t716 * t189 + 0.317e3 / 0.162e3 * t716 * t182;
362 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
363 v4rho3sigma[0] = t234 * t1 * t3 * t101 / 0.18e2 - t320 * t195 / 0.4e1 - 0.3e1 / 0.4e1 * t164 * t385 - 0.3e1 / 0.16e2 * t44 * t45 * (0.35e2 / 0.81e2 * t84 * t248 - 0.91e2 / 0.54e2 * t95 * t258 + 0.49e2 / 0.486e3 * t130 * t627 * t347 - t130 * t633 * t33 * t131 * t532 / 0.81e2 + 0.91e2 / 0.54e2 * t87 * t254 - 0.317e3 / 0.81e2 * t641 * t182 - 0.317e3 / 0.162e3 * t641 * t177 + 0.49e2 / 0.162e3 * t130 * t627 * t72 * t353 + 0.49e2 / 0.162e3 * t130 * t627 * t147 * t359 - 0.11e2 / 0.243e3 * t130 * t633 * t72 * t131 * t12 * t559 + t721);
364
365 t743 = 0.1e1 / t59 / t277;
366 t744 = t743 * t529;
367 t755 = t94 * t18;
368 t765 = t743 * t30;
369 t794 = t130 * t82 * t134;
370 t812 = -t130 * t765 * t300 * t696 * t12 * t755 / 0.36e2 - 0.11e2 / 0.648e3 * t130 * t765 * t147 * t703 * t12 * t755 - t130 * t765 * t72 * t710 * t12 * t755 / 0.216e3 - 0.7e1 / 0.216e3 * t205 * t122 - 0.16e2 / 0.27e2 * t794 * t151 - 0.8e1 / 0.27e2 * t794 * t156 - 0.35e2 / 0.1296e4 * t345 * t416 + 0.7e1 / 0.108e3 * t200 * t115 - 0.35e2 / 0.432e3 * t345 * t291 - 0.35e2 / 0.432e3 * t345 * t294 + t130 * t744 * t33 * t12 * t755 / 0.216e3;
371 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
372 v4rho2sigma2[0] = -t320 * t230 / 0.12e2 - t164 * t440 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t45 * (0.16e2 / 0.27e2 * t396 * t141 + 0.8e1 / 0.27e2 * t396 * t138 + 0.7e1 / 0.216e3 * t209 * t126 + 0.35e2 / 0.432e3 * t364 * t426 + 0.35e2 / 0.432e3 * t371 * t431 + 0.35e2 / 0.1296e4 * t378 * t436 + 0.11e2 / 0.648e3 * t130 * t744 * t72 * t531 * t9 * params->c * t18 + t130 * t744 * t147 * t149 * t12 * t755 / 0.36e2 + t130 * t744 * t300 * t301 * t12 * t755 / 0.54e2 - t130 * t765 * t581 * t689 * t12 * t755 / 0.54e2 + t812);
373
374 t819 = t130 * t198;
375 t835 = 0.1e1 / t59 / t261;
376 t836 = t835 * t529;
377 t864 = t835 * t30;
378 t892 = t130 * t198 * t174;
379 t904 = t130 * t412 * t72 * t469 / 0.72e2 + t130 * t412 * t147 * t475 / 0.72e2 + t130 * t864 * t581 * t689 * t82 * t532 / 0.144e3 + t130 * t864 * t300 * t696 * t82 * t532 / 0.96e2 + 0.11e2 / 0.1728e4 * t130 * t864 * t147 * t703 * t82 * t532 + t130 * t864 * t72 * t710 * t82 * t532 / 0.576e3 - t451 * t65 / 0.48e2 - t892 * t151 / 0.24e2 - t892 * t156 / 0.48e2 + t446 * t50 / 0.24e2 - t130 * t836 * t33 * t82 * t532 / 0.576e3;
380 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
381 v4rhosigma3[0] = -t164 * t498 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t45 * (t819 * t401 / 0.24e2 - t425 * t482 / 0.72e2 - t430 * t488 / 0.72e2 - t435 * t494 / 0.216e3 + t819 * t398 / 0.48e2 + t455 * t75 / 0.48e2 + t413 * t463 / 0.216e3 - 0.11e2 / 0.1728e4 * t130 * t836 * t72 * t82 * t12 * t559 - t130 * t836 * t147 * t82 * t149 * t532 / 0.96e2 - t130 * t836 * t300 * t301 * t82 * t532 / 0.144e3 + t904);
382
383 t909 = t276 * sigma[0];
384 t911 = t55 / t909;
385 t920 = t130 * t444 * t212;
386 t925 = t130 * t444;
387 t926 = t212 * t136;
388 t950 = t545 * t529;
389 t1012 = 0.11e2 / 0.4608e4 * t130 * t950 * t72 * t198 * t12 * t559 + t130 * t950 * t147 * t198 * t149 * t532 / 0.256e3 + t130 * t950 * t300 * t301 * t198 * t532 / 0.384e3 + t467 * t449 * params->c * t285 / 0.384e3 + t473 * t149 * t449 * t285 / 0.384e3 - t130 * t546 * t581 * t689 * t198 * t532 / 0.384e3 - t130 * t546 * t300 * t696 * t198 * t532 / 0.256e3 - 0.11e2 / 0.4608e4 * t130 * t546 * t147 * t703 * t198 * t532 - t130 * t546 * t72 * t710 * t198 * t532 / 0.1536e4 + 0.5e1 / 0.64e2 * t14 * t17 / t18 / t909 * t35 + t130 * t950 * t33 * t198 * t532 / 0.1536e4;
390 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
391 v4sigma4[0] = -0.3e1 / 0.16e2 * t44 * t45 * (0.5e1 / 0.128e3 * t53 * t911 * t9 * t97 - 0.5e1 / 0.128e3 * t53 * t911 * t91 + t920 * t151 / 0.64e2 + t920 * t156 / 0.128e3 - t925 * t926 * t73 / 0.64e2 - t480 * t365 * t449 * t285 / 0.384e3 - t486 * t372 * t449 * t285 / 0.384e3 - t492 * t379 * t449 * t285 / 0.1152e4 - t925 * t926 * t33 / 0.128e3 + t461 * t33 * t449 * t285 / 0.1152e4 + t1012);
392
393 #ifndef XC_DONT_COMPILE_MXC
394
395 if(order < 5) return;
396
397
398 #endif
399
400 #endif
401
402 #endif
403
404 #endif
405
406 #endif
407
408
409 }
410
411
412 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 *))413 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
414 {
415
416 #ifndef XC_DONT_COMPILE_EXC
417 double t1, t3, t4, t5, t6, t7, t8, t9;
418 double t10, t12, t13, t14, t15, t16, t17, t18;
419 double t19, t20, t23, t24, t26, t27, t28, t32;
420 double t34;
421
422 #ifndef XC_DONT_COMPILE_VXC
423 double t37, t38, t39, t41, t46, t47, t48, t49;
424 double t50, t51, t53, t55, t56, t61, t63, t64;
425 double t65, t66, t67, t70, t74, t80, t87, t89;
426 double t93;
427
428 #ifndef XC_DONT_COMPILE_FXC
429 double t96, t97, t104, t109, t111, t117, t121, t122;
430 double t123, t124, t125, t126, t127, t128, t129, t132;
431 double t136, t138, t140, t142, t146, t149, t153, t168;
432 double t169, t170, t171, t172, t175, t177, t181, t182;
433 double t184, t190, t193, t198, t204, t206, t209, t210;
434 double t215, t219, t227;
435
436 #ifndef XC_DONT_COMPILE_KXC
437 double t231, t232, t243, t249, t258, t259, t260, t268;
438 double t273, t274, t276, t277, t278, t280, t282, t287;
439 double t288, t294, t298, t299, t300, t301, t302, t305;
440 double t312, t316, t332, t338, t344, t346, t348, t352;
441 double t356, t357, t360, t363, t364, t369, t370, t375;
442 double t376, t379, t380, t395, t396, t397, t400, t404;
443 double t410, t415, t416, t419, t422, t423, t430, t431;
444 double t434, t435, t440, t446, t450, t452, t454, t458;
445 double t462, t463, t466, t468, t469, t473, t474, t478;
446 double t479, t483;
447
448 #ifndef XC_DONT_COMPILE_LXC
449 double t502, t504, t509, t512, t513, t519, t520, t521;
450 double t522, t524, t529, t534, t538, t548, t549, t554;
451 double t555, t561, t568, t576, t578, t579, t592, t613;
452 double t619, t623, t625, t630, t647, t657, t670, t673;
453 double t679, t685, t691, t695, t708, t728, t737, t739;
454 double t746, t755, t781, t793, t795, t805, t809, t829;
455 double t855, t870, t875, t876, t886, t896, t900, t901;
456 double t906, t951, t967;
457 #endif
458
459 #endif
460
461 #endif
462
463 #endif
464
465 #endif
466
467
468 gga_x_cap_params *params;
469
470 assert(p->params != NULL);
471 params = (gga_x_cap_params * )(p->params);
472
473 t1 = M_CBRT3;
474 t3 = POW_1_3(0.1e1 / M_PI);
475 t4 = t1 * t3;
476 t5 = M_CBRT4;
477 t6 = t5 * t5;
478 t7 = POW_1_3(rho[0]);
479 t8 = t6 * t7;
480 t9 = M_CBRT6;
481 t10 = t9 * t9;
482 t12 = M_PI * M_PI;
483 t13 = POW_1_3(t12);
484 t14 = 0.1e1 / t13;
485 t15 = params->alphaoAx * t10 * t14;
486 t16 = sqrt(sigma[0]);
487 t17 = t7 * rho[0];
488 t18 = 0.1e1 / t17;
489 t19 = t16 * t18;
490 t20 = t10 * t14;
491 t23 = 0.1e1 + t20 * t19 / 0.12e2;
492 t24 = log(t23);
493 t26 = params->c * t24 + 0.1e1;
494 t27 = 0.1e1 / t26;
495 t28 = t24 * t27;
496 t32 = 0.1e1 - t15 * t19 * t28 / 0.12e2;
497 t34 = t4 * t8 * t32;
498 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
499 zk[0] = -0.3e1 / 0.8e1 * t34;
500
501 #ifndef XC_DONT_COMPILE_VXC
502
503 if(order < 1) return;
504
505
506 t37 = t17 * t1;
507 t38 = t3 * t6;
508 t39 = rho[0] * rho[0];
509 t41 = 0.1e1 / t7 / t39;
510 t46 = params->alphaoAx * t9;
511 t47 = t13 * t13;
512 t48 = 0.1e1 / t47;
513 t49 = t46 * t48;
514 t50 = t39 * rho[0];
515 t51 = t7 * t7;
516 t53 = 0.1e1 / t51 / t50;
517 t55 = 0.1e1 / t23;
518 t56 = t55 * t27;
519 t61 = t46 * t48 * sigma[0];
520 t63 = t26 * t26;
521 t64 = 0.1e1 / t63;
522 t65 = t64 * params->c;
523 t66 = t65 * t55;
524 t67 = t53 * t24 * t66;
525 t70 = t15 * t16 * t41 * t28 / 0.9e1 + t49 * sigma[0] * t53 * t56 / 0.18e2 - t61 * t67 / 0.18e2;
526 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
527 vrho[0] = -t34 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * t70;
528
529 t74 = 0.1e1 / t16;
530 t80 = 0.1e1 / t51 / t39;
531 t87 = t24 * t64;
532 t89 = t87 * params->c * t55;
533 t93 = t38 * (-t15 * t74 * t18 * t28 / 0.24e2 - t49 * t80 * t55 * t27 / 0.48e2 + t46 * t48 * t80 * t89 / 0.48e2);
534 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
535 vsigma[0] = -0.3e1 / 0.8e1 * t37 * t93;
536
537 #ifndef XC_DONT_COMPILE_FXC
538
539 if(order < 2) return;
540
541
542 t96 = 0.1e1 / t51;
543 t97 = t6 * t96;
544 t104 = 0.1e1 / t7 / t50;
545 t109 = t39 * t39;
546 t111 = 0.1e1 / t51 / t109;
547 t117 = t111 * t24 * t66;
548 t121 = params->alphaoAx / t12;
549 t122 = t16 * sigma[0];
550 t123 = t121 * t122;
551 t124 = t109 * t39;
552 t125 = 0.1e1 / t124;
553 t126 = t23 * t23;
554 t127 = 0.1e1 / t126;
555 t128 = t125 * t127;
556 t129 = t128 * t27;
557 t132 = t128 * t65;
558 t136 = t121 * t122 * t125;
559 t138 = 0.1e1 / t63 / t26;
560 t140 = params->c * params->c;
561 t142 = t24 * t138 * t140 * t127;
562 t146 = t87 * params->c * t127;
563 t149 = -0.7e1 / 0.27e2 * t15 * t16 * t104 * t28 - 0.5e1 / 0.18e2 * t49 * sigma[0] * t111 * t56 + 0.5e1 / 0.18e2 * t61 * t117 + t123 * t129 / 0.27e2 + 0.2e1 / 0.27e2 * t123 * t132 - 0.2e1 / 0.27e2 * t136 * t142 - t136 * t146 / 0.27e2;
564 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
565 v2rho2[0] = -t4 * t97 * t32 / 0.6e1 - t4 * t8 * t70 - 0.3e1 / 0.8e1 * t37 * t38 * t149;
566
567 t153 = t7 * t1;
568 t168 = t109 * rho[0];
569 t169 = 0.1e1 / t168;
570 t170 = t121 * t169;
571 t171 = t127 * t27;
572 t172 = t171 * t16;
573 t175 = t127 * t64;
574 t177 = t175 * params->c * t16;
575 t181 = t121 * t169 * t24;
576 t182 = t138 * t140;
577 t184 = t182 * t127 * t16;
578 t190 = t38 * (t15 * t74 * t41 * t28 / 0.18e2 + t49 * t53 * t55 * t27 / 0.12e2 - t46 * t48 * t53 * t89 / 0.12e2 - t170 * t172 / 0.72e2 - t170 * t177 / 0.36e2 + t181 * t184 / 0.36e2 + t181 * t177 / 0.72e2);
579 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
580 v2rhosigma[0] = -t153 * t93 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t190;
581
582 t193 = 0.1e1 / t122;
583 t198 = 0.1e1 / sigma[0];
584 t204 = t46 * t48 * t198;
585 t206 = t80 * t24 * t66;
586 t209 = 0.1e1 / t109;
587 t210 = t121 * t209;
588 t215 = t175 * params->c * t74;
589 t219 = t121 * t209 * t24;
590 t227 = t38 * (t15 * t193 * t18 * t28 / 0.48e2 - t49 * t198 * t80 * t56 / 0.96e2 + t204 * t206 / 0.96e2 + t210 * t171 * t74 / 0.192e3 + t210 * t215 / 0.96e2 - t219 * t182 * t127 * t74 / 0.96e2 - t219 * t215 / 0.192e3);
591 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
592 v2sigma2[0] = -0.3e1 / 0.8e1 * t37 * t227;
593
594 #ifndef XC_DONT_COMPILE_KXC
595
596 if(order < 3) return;
597
598
599 t231 = 0.1e1 / t51 / rho[0];
600 t232 = t6 * t231;
601 t243 = 0.1e1 / t7 / t109;
602 t249 = 0.1e1 / t51 / t168;
603 t258 = t109 * t50;
604 t259 = 0.1e1 / t258;
605 t260 = t259 * t127;
606 t268 = t121 * t122 * t259;
607 t273 = sigma[0] * sigma[0];
608 t274 = t109 * t109;
609 t276 = 0.1e1 / t7 / t274;
610 t277 = t273 * t276;
611 t278 = t121 * t277;
612 t280 = 0.1e1 / t126 / t23;
613 t282 = t280 * t27 * t20;
614 t287 = params->c * t10 * t14;
615 t288 = t280 * t64 * t287;
616 t294 = t280 * t138 * t140 * t10 * t14;
617 t298 = t121 * t277 * t24;
618 t299 = t63 * t63;
619 t300 = 0.1e1 / t299;
620 t301 = t140 * params->c;
621 t302 = t300 * t301;
622 t305 = t302 * t280 * t10 * t14;
623 t312 = 0.70e2 / 0.81e2 * t15 * t16 * t243 * t28 + 0.119e3 / 0.81e2 * t49 * sigma[0] * t249 * t56 - 0.119e3 / 0.81e2 * t61 * t249 * t24 * t66 - 0.11e2 / 0.27e2 * t123 * t260 * t27 - 0.22e2 / 0.27e2 * t123 * t260 * t65 + 0.22e2 / 0.27e2 * t268 * t142 + 0.11e2 / 0.27e2 * t268 * t146 + 0.2e1 / 0.243e3 * t278 * t282 + 0.2e1 / 0.81e2 * t278 * t288 + 0.2e1 / 0.81e2 * t278 * t294 - 0.2e1 / 0.81e2 * t298 * t305 - 0.2e1 / 0.81e2 * t298 * t294 - 0.2e1 / 0.243e3 * t298 * t288;
624 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
625 v3rho3[0] = t4 * t232 * t32 / 0.9e1 - t4 * t97 * t70 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t149 - 0.3e1 / 0.8e1 * t37 * t38 * t312;
626
627 t316 = t96 * t1;
628 t332 = t121 * t125;
629 t338 = t121 * t125 * t24;
630 t344 = 0.1e1 / t7 / t258;
631 t346 = t121 * t344 * t280;
632 t348 = t27 * sigma[0] * t20;
633 t352 = t64 * sigma[0] * t287;
634 t356 = sigma[0] * t10 * t14;
635 t357 = t182 * t356;
636 t360 = t344 * t24;
637 t363 = t301 * t280;
638 t364 = t363 * t356;
639 t369 = t140 * t280;
640 t370 = t369 * t356;
641 t375 = params->c * t280;
642 t376 = t375 * t356;
643 t379 = -0.7e1 / 0.54e2 * t15 * t74 * t104 * t28 - 0.37e2 / 0.108e3 * t49 * t111 * t55 * t27 + 0.37e2 / 0.108e3 * t46 * t48 * t111 * t89 + t332 * t172 / 0.8e1 + t332 * t177 / 0.4e1 - t338 * t184 / 0.4e1 - t338 * t177 / 0.8e1 - t346 * t348 / 0.324e3 - t346 * t352 / 0.108e3 - t346 * t357 / 0.108e3 + t121 * t360 * t300 * t364 / 0.108e3 + t121 * t360 * t138 * t370 / 0.108e3 + t121 * t360 * t64 * t376 / 0.324e3;
644 t380 = t38 * t379;
645 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
646 v3rho2sigma[0] = -t316 * t93 / 0.6e1 - t153 * t190 - 0.3e1 / 0.8e1 * t37 * t380;
647
648 t395 = t121 * t74;
649 t396 = t169 * t127;
650 t397 = t396 * t27;
651 t400 = t396 * t65;
652 t404 = t121 * t74 * t169;
653 t410 = 0.1e1 / t7 / t124;
654 t415 = t121 * t410 * t280;
655 t416 = t65 * t20;
656 t419 = t182 * t20;
657 t422 = t410 * t24;
658 t423 = t121 * t422;
659 t430 = -t15 * t193 * t41 * t28 / 0.36e2 + t49 * t198 * t53 * t56 / 0.72e2 - t204 * t67 / 0.72e2 - t395 * t397 / 0.36e2 - t395 * t400 / 0.18e2 + t404 * t142 / 0.18e2 + t404 * t146 / 0.36e2 + t121 * t410 * t282 / 0.864e3 + t415 * t416 / 0.288e3 + t415 * t419 / 0.288e3 - t423 * t305 / 0.288e3 - t423 * t294 / 0.288e3 - t423 * t288 / 0.864e3;
660 t431 = t38 * t430;
661 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
662 v3rhosigma2[0] = -t153 * t227 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t431;
663
664 t434 = t16 * t273;
665 t435 = 0.1e1 / t434;
666 t440 = 0.1e1 / t273;
667 t446 = t46 * t48 * t440;
668 t450 = 0.1e1 / t7 / t168;
669 t452 = t121 * t450 * t280;
670 t454 = t27 * t198 * t20;
671 t458 = t64 * t198 * t287;
672 t462 = t198 * t10 * t14;
673 t463 = t182 * t462;
674 t466 = t450 * t24;
675 t468 = t121 * t466 * t300;
676 t469 = t363 * t462;
677 t473 = t121 * t466 * t138;
678 t474 = t369 * t462;
679 t478 = t121 * t466 * t64;
680 t479 = t375 * t462;
681 t483 = t38 * (-t15 * t435 * t18 * t28 / 0.32e2 + t49 * t440 * t80 * t56 / 0.64e2 - t446 * t206 / 0.64e2 - t452 * t454 / 0.2304e4 - t452 * t458 / 0.768e3 - t452 * t463 / 0.768e3 + t468 * t469 / 0.768e3 + t473 * t474 / 0.768e3 + t478 * t479 / 0.2304e4);
682 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
683 v3sigma3[0] = -0.3e1 / 0.8e1 * t37 * t483;
684
685 #ifndef XC_DONT_COMPILE_LXC
686
687 if(order < 4) return;
688
689
690 t502 = 0.1e1 / t274;
691 t504 = t121 * t122 * t502;
692 t509 = t274 * rho[0];
693 t512 = t273 / t7 / t509;
694 t513 = t121 * t512;
695 t519 = t434 / t51 / t274 / t39;
696 t520 = t121 * t519;
697 t521 = t126 * t126;
698 t522 = 0.1e1 / t521;
699 t524 = t9 * t48;
700 t529 = 0.1e1 / t51 / t124;
701 t534 = t502 * t127;
702 t538 = t529 * t24;
703 t548 = t524 * params->c;
704 t549 = t522 * t64 * t548;
705 t554 = t140 * t9 * t48;
706 t555 = t522 * t138 * t554;
707 t561 = t522 * t300 * t301 * t9 * t48;
708 t568 = t121 * t512 * t24;
709 t576 = t121 * t519 * t24;
710 t578 = 0.1e1 / t299 / t26;
711 t579 = t140 * t140;
712 t592 = 0.44e2 / 0.729e3 * t520 * t549 + 0.8e1 / 0.81e2 * t520 * t555 + 0.16e2 / 0.243e3 * t520 * t561 + 0.931e3 / 0.243e3 * t123 * t534 * t27 + 0.116e3 / 0.243e3 * t568 * t305 + 0.116e3 / 0.243e3 * t568 * t294 + 0.116e3 / 0.729e3 * t568 * t288 - 0.16e2 / 0.243e3 * t576 * t578 * t579 * t522 * t9 * t48 - 0.8e1 / 0.81e2 * t576 * t561 - 0.44e2 / 0.729e3 * t576 * t555 - 0.4e1 / 0.243e3 * t576 * t549;
713 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
714 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t80 * t32 + 0.4e1 / 0.9e1 * t4 * t232 * t70 - t4 * t97 * t149 - 0.2e1 * t4 * t8 * t312 - 0.3e1 / 0.8e1 * t37 * t38 * (-0.910e3 / 0.243e3 * t15 * t16 * t450 * t28 - 0.1862e4 / 0.243e3 * t504 * t142 - 0.931e3 / 0.243e3 * t504 * t146 - 0.116e3 / 0.729e3 * t513 * t282 + 0.4e1 / 0.243e3 * t520 * t522 * t27 * t524 - 0.721e3 / 0.81e2 * t49 * sigma[0] * t529 * t56 + 0.1862e4 / 0.243e3 * t123 * t534 * t65 + 0.721e3 / 0.81e2 * t61 * t538 * t66 - 0.116e3 / 0.243e3 * t513 * t288 - 0.116e3 / 0.243e3 * t513 * t294 + t592);
715
716 t613 = t121 * t259 * t24;
717 t619 = t121 * t276 * t280;
718 t623 = 0.1e1 / t51 / t509;
719 t625 = t121 * t623 * t522;
720 t630 = t121 * t259;
721 t647 = t122 * t9 * t48;
722 t657 = t276 * t24;
723 t670 = t623 * t24;
724 t673 = t579 * t522;
725 t679 = t301 * t522;
726 t685 = t140 * t522;
727 t691 = params->c * t522;
728 t695 = -t625 * t138 * t122 * t554 / 0.27e2 - 0.2e1 / 0.81e2 * t625 * t302 * t647 + 0.91e2 / 0.54e2 * t49 * t249 * t55 * t27 - 0.317e3 / 0.324e3 * t630 * t172 - 0.49e2 / 0.324e3 * t121 * t657 * t300 * t364 - 0.49e2 / 0.324e3 * t121 * t657 * t138 * t370 - 0.49e2 / 0.972e3 * t121 * t657 * t64 * t376 + 0.2e1 / 0.81e2 * t121 * t670 * t578 * t673 * t647 + t121 * t670 * t300 * t679 * t647 / 0.27e2 + 0.11e2 / 0.486e3 * t121 * t670 * t138 * t685 * t647 + t121 * t670 * t64 * t691 * t647 / 0.162e3;
729 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
730 v4rho3sigma[0] = t231 * t1 * t93 / 0.9e1 - t316 * t190 / 0.2e1 - 0.3e1 / 0.2e1 * t153 * t380 - 0.3e1 / 0.8e1 * t37 * t38 * (0.35e2 / 0.81e2 * t15 * t74 * t243 * t28 - 0.91e2 / 0.54e2 * t46 * t48 * t249 * t89 + 0.317e3 / 0.162e3 * t613 * t184 + 0.317e3 / 0.324e3 * t613 * t177 + 0.49e2 / 0.972e3 * t619 * t348 - t625 * t27 * t122 * t524 / 0.162e3 - 0.317e3 / 0.162e3 * t630 * t177 + 0.49e2 / 0.324e3 * t619 * t352 + 0.49e2 / 0.324e3 * t619 * t357 - 0.11e2 / 0.486e3 * t625 * t64 * t122 * t548 + t695);
731
732 t708 = t121 * t74 * t125;
733 t728 = t121 * t360;
734 t737 = 0.1e1 / t51 / t274;
735 t739 = t121 * t737 * t522;
736 t746 = t524 * t16;
737 t755 = t737 * t24;
738 t781 = 0.35e2 / 0.864e3 * t728 * t294 + 0.35e2 / 0.2592e4 * t728 * t288 + 0.11e2 / 0.1296e4 * t739 * t64 * t9 * t48 * params->c * t16 + t739 * t182 * t746 / 0.72e2 + t739 * t302 * t746 / 0.108e3 + 0.4e1 / 0.27e2 * t395 * t129 - t121 * t755 * t578 * t673 * t746 / 0.108e3 - t121 * t755 * t300 * t679 * t746 / 0.72e2 - 0.11e2 / 0.1296e4 * t121 * t755 * t138 * t685 * t746 - t121 * t755 * t64 * t691 * t746 / 0.432e3 + t739 * t27 * t9 * t48 * t16 / 0.432e3;
739 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
740 v4rho2sigma2[0] = -t316 * t227 / 0.6e1 - t153 * t431 - 0.3e1 / 0.8e1 * t37 * t38 * (0.7e1 / 0.108e3 * t15 * t193 * t104 * t28 - 0.8e1 / 0.27e2 * t708 * t142 - 0.4e1 / 0.27e2 * t708 * t146 - 0.35e2 / 0.864e3 * t346 * t416 - 0.35e2 / 0.864e3 * t346 * t419 - 0.7e1 / 0.216e3 * t49 * t198 * t111 * t56 + 0.8e1 / 0.27e2 * t395 * t132 - 0.35e2 / 0.2592e4 * t121 * t344 * t282 + 0.7e1 / 0.216e3 * t204 * t117 + 0.35e2 / 0.864e3 * t728 * t305 + t781);
741
742 t793 = 0.1e1 / t51 / t258;
743 t795 = t121 * t793 * t522;
744 t805 = t74 * t9 * t48;
745 t809 = t121 * t193;
746 t829 = t793 * t24;
747 t855 = t121 * t193 * t169;
748 t870 = t415 * t463 / 0.144e3 + t121 * t829 * t578 * t673 * t805 / 0.288e3 + t121 * t829 * t300 * t679 * t805 / 0.192e3 + 0.11e2 / 0.3456e4 * t121 * t829 * t138 * t685 * t805 + t121 * t829 * t64 * t691 * t805 / 0.1152e4 + t15 * t435 * t41 * t28 / 0.24e2 - t855 * t142 / 0.48e2 - t855 * t146 / 0.96e2 - t795 * t27 * t74 * t524 / 0.1152e4 - t49 * t440 * t53 * t56 / 0.48e2 + t809 * t400 / 0.48e2;
749 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
750 v4rhosigma3[0] = -t153 * t483 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * (t446 * t67 / 0.48e2 + t415 * t454 / 0.432e3 - 0.11e2 / 0.3456e4 * t795 * t64 * t74 * t548 - t795 * t138 * t74 * t554 / 0.192e3 - t795 * t302 * t805 / 0.288e3 + t809 * t397 / 0.96e2 - t121 * t422 * t300 * t469 / 0.144e3 - t121 * t422 * t138 * t474 / 0.144e3 - t121 * t422 * t64 * t479 / 0.432e3 + t415 * t458 / 0.144e3 + t870);
751
752 t875 = t273 * sigma[0];
753 t876 = 0.1e1 / t875;
754 t886 = t121 * t529 * t522;
755 t896 = t193 * t9 * t48;
756 t900 = t121 * t435;
757 t901 = t209 * t127;
758 t906 = t440 * t10 * t14;
759 t951 = t121 * t435 * t209;
760 t967 = t452 * t182 * t906 / 0.768e3 - t121 * t538 * t578 * t673 * t896 / 0.768e3 - t121 * t538 * t300 * t679 * t896 / 0.512e3 - 0.11e2 / 0.9216e4 * t121 * t538 * t138 * t685 * t896 - t121 * t538 * t64 * t691 * t896 / 0.3072e4 + 0.5e1 / 0.64e2 * t15 / t16 / t875 * t18 * t28 + t951 * t142 / 0.128e3 + t951 * t146 / 0.256e3 + t886 * t27 * t193 * t524 / 0.3072e4 - 0.5e1 / 0.128e3 * t49 * t876 * t80 * t56 - t900 * t901 * t65 / 0.128e3;
761 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
762 v4sigma4[0] = -0.3e1 / 0.8e1 * t37 * t38 * (0.5e1 / 0.128e3 * t46 * t48 * t876 * t206 + t452 * t27 * t440 * t20 / 0.2304e4 + 0.11e2 / 0.9216e4 * t886 * t64 * t193 * t548 + t886 * t138 * t193 * t554 / 0.512e3 + t886 * t302 * t896 / 0.768e3 - t900 * t901 * t27 / 0.256e3 - t468 * t363 * t906 / 0.768e3 - t473 * t369 * t906 / 0.768e3 - t478 * t375 * t906 / 0.2304e4 + t452 * t64 * t440 * t287 / 0.768e3 + t967);
763
764 #ifndef XC_DONT_COMPILE_MXC
765
766 if(order < 5) return;
767
768
769 #endif
770
771 #endif
772
773 #endif
774
775 #endif
776
777 #endif
778
779
780 }
781
782
783 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 *))784 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
785 {
786
787 #ifndef XC_DONT_COMPILE_EXC
788 double t1, t3, t4, t5, t6, t7, t8, t9;
789 double t10, t11, t13, t14, t15, t16, t17, t18;
790 double t19, t21, t22, t23, t24, t25, t26, t28;
791 double t29, t30, t33, t34, t36, t37, t38, t42;
792 double t44, t46, t47, t48, t49, t50, t51, t53;
793 double t54, t57, t58, t60, t61, t62, t66, t68;
794
795 #ifndef XC_DONT_COMPILE_VXC
796 double t70, t71, t72, t73, t74, t75, t77, t78;
797 double t80, t81, t82, t83, t84, t86, t87, t88;
798 double t90, t95, t96, t97, t98, t99, t100, t102;
799 double t104, t105, t110, t112, t113, t114, t115, t116;
800 double t119, t121, t122, t123, t124, t125, t127, t128;
801 double t129, t131, t132, t136, t137, t139, t140, t141;
802 double t142, t144, t145, t146, t148, t153, t154, t156;
803 double t158, t159, t164, t166, t167, t168, t169, t170;
804 double t173, t175, t176, t180, t181, t182, t183, t189;
805 double t196, t198, t201, t205, t206, t212, t219, t221;
806 double t224;
807
808 #ifndef XC_DONT_COMPILE_FXC
809 double t228, t230, t231, t232, t233, t234, t235, t237;
810 double t238, t239, t241, t242, t243, t245, t246, t247;
811 double t248, t249, t250, t252, t253, t255, t256, t258;
812 double t259, t261, t262, t264, t269, t271, t277, t281;
813 double t282, t283, t284, t285, t286, t287, t288, t289;
814 double t292, t296, t298, t300, t302, t306, t309, t311;
815 double t312, t313, t314, t315, t316, t317, t319, t320;
816 double t321, t323, t324, t325, t326, t328, t329, t330;
817 double t332, t333, t334, t336, t337, t338, t339, t341;
818 double t344, t346, t348, t351, t352, t353, t355, t359;
819 double t360, t361, t362, t364, t367, t369, t370, t372;
820 double t375, t377, t380, t382, t385, t386, t388, t389;
821 double t390, t391, t392, t394, t395, t396, t397, t399;
822 double t400, t401, t402, t404, t405, t406, t408, t409;
823 double t410, t412, t417, t419, t425, t428, t429, t430;
824 double t431, t432, t433, t434, t435, t438, t442, t444;
825 double t447, t451, t454, t456, t457, t458, t461, t463;
826 double t464, t480, t481, t482, t483, t484, t487, t489;
827 double t493, t494, t496, t501, t506, t507, t508, t513;
828 double t518, t533, t534, t535, t536, t537, t540, t542;
829 double t546, t547, t549, t554, t555, t558, t563, t569;
830 double t571, t574, t575, t580, t584, t591, t592, t595;
831 double t600, t606, t608, t611, t612, t617, t621, t628;
832 double t629;
833
834 #ifndef XC_DONT_COMPILE_KXC
835 double t636, t638, t642, t643, t644, t645, t647, t648;
836 double t649, t651, t653, t655, t656, t657, t659, t660;
837 double t662, t663, t664, t665, t667, t668, t669, t671;
838 double t673, t675, t676, t677, t679, t680, t681, t683;
839 double t684, t686, t689, t691, t693, t694, t695, t697;
840 double t698, t699, t700, t701, t703, t704, t706, t707;
841 double t708, t710, t711, t713, t715, t716, t718, t719;
842 double t720, t722, t723, t725, t726, t728, t731, t737;
843 double t746, t747, t748, t756, t761, t762, t764, t765;
844 double t766, t768, t770, t775, t776, t781, t782, t786;
845 double t787, t788, t789, t790, t793, t800, t802, t803;
846 double t804, t809, t810, t813, t814, t815, t817, t819;
847 double t821, t822, t823, t824, t826, t828, t829, t830;
848 double t832, t834, t836, t837, t838, t840, t842, t843;
849 double t845, t847, t849, t851, t853, t858, t861, t864;
850 double t866, t875, t878, t879, t881, t883, t885, t887;
851 double t889, t891, t894, t895, t897, t898, t899, t901;
852 double t903, t905, t907, t909, t911, t914, t919, t921;
853 double t923, t925, t927, t929, t931, t933, t935, t937;
854 double t939, t941, t943, t945, t948, t954, t960, t963;
855 double t965, t967, t969, t971, t974, t975, t977, t979;
856 double t981, t982, t983, t985, t987, t989, t992, t995;
857 double t998, t1000, t1003, t1011, t1012, t1015, t1018, t1019;
858 double t1021, t1022, t1023, t1025, t1026, t1028, t1030, t1031;
859 double t1032, t1034, t1035, t1036, t1037, t1039, t1040, t1042;
860 double t1048, t1057, t1058, t1059, t1067, t1072, t1073, t1075;
861 double t1076, t1077, t1079, t1081, t1085, t1089, t1093, t1094;
862 double t1095, t1096, t1099, t1106, t1108, t1109, t1112, t1113;
863 double t1115, t1116, t1117, t1119, t1120, t1121, t1122, t1124;
864 double t1125, t1126, t1128, t1129, t1130, t1134, t1138, t1140;
865 double t1141, t1145, t1148, t1164, t1170, t1176, t1178, t1180;
866 double t1184, t1188, t1189, t1192, t1195, t1196, t1201, t1202;
867 double t1207, t1208, t1211, t1216, t1220, t1221, t1225, t1231;
868 double t1232, t1233, t1235, t1236, t1237, t1238, t1242, t1245;
869 double t1246, t1248, t1252, t1254, t1255, t1259, t1262, t1265;
870 double t1270, t1274, t1279, t1283, t1286, t1301, t1307, t1313;
871 double t1315, t1317, t1321, t1325, t1326, t1329, t1332, t1333;
872 double t1338, t1339, t1344, t1345, t1348, t1349, t1352, t1355;
873 double t1368, t1369, t1370, t1373, t1377, t1383, t1388, t1389;
874 double t1392, t1395, t1396, t1403, t1404, t1407, t1409, t1413;
875 double t1418, t1431, t1432, t1433, t1436, t1440, t1446, t1451;
876 double t1452, t1455, t1458, t1459, t1466, t1467, t1470, t1471;
877 double t1476, t1482, t1486, t1488, t1490, t1494, t1498, t1499;
878 double t1502, t1504, t1505, t1509, t1510, t1514, t1515, t1518;
879 double t1519, t1522, t1523, t1528, t1534, t1538, t1540, t1542;
880 double t1546, t1550, t1551, t1554, t1556, t1557, t1561, t1562;
881 double t1566, t1567, t1570, t1571;
882
883 #ifndef XC_DONT_COMPILE_LXC
884 double t1585, t1589, t1592, t1596, t1598, t1600, t1602, t1608;
885 double t1610, t1612, t1613, t1614, t1619, t1621, t1623, t1625;
886 double t1627, t1633, t1636, t1642, t1645, t1648, t1667, t1669;
887 double t1670, t1671, t1676, t1678, t1680, t1681, t1683, t1685;
888 double t1695, t1698, t1701, t1708, t1723, t1729, t1732, t1735;
889 double t1737, t1738, t1742, t1745, t1747, t1757, t1759, t1761;
890 double t1762, t1764, t1765, t1767, t1771, t1774, t1777, t1781;
891 double t1790, t1794, t1799, t1811, t1819, t1823, t1831, t1835;
892 double t1836, t1839, t1856, t1857, t1858, t1860, t1872, t1873;
893 double t1883, t1884, t1886, t1889, t1892, t1894, t1899, t1901;
894 double t1907, t1909, t1918, t1919, t1921, t1929, t1932, t1935;
895 double t1946, t1958, t1965, t1971, t1973, t1976, t1982, t1984;
896 double t1990, t1996, t2003, t2005, t2014, t2020, t2023, t2025;
897 double t2030, t2032, t2035, t2041, t2045, t2052, t2059, t2065;
898 double t2071, t2074, t2077, t2080, t2083, t2086, t2089, t2091;
899 double t2097, t2100, t2123, t2125, t2137, t2146, t2171, t2181;
900 double t2188, t2195, t2197, t2200, t2207, t2232, t2252, t2254;
901 double t2270, t2280, t2295, t2304, t2358, t2364, t2366, t2369;
902 double t2372, t2375, t2378, t2381, t2384, t2387, t2391, t2423;
903 double t2440, t2459, t2461, t2469, t2472, t2480, t2484, t2494;
904 double t2518, t2523, t2528, t2533, t2538, t2577, t2594, t2595;
905 double t2599, t2602, t2604, t2614, t2616, t2618, t2620, t2621;
906 double t2623, t2627, t2630, t2633, t2637, t2638, t2642, t2648;
907 double t2660, t2678, t2683, t2686, t2688, t2691, t2697, t2700;
908 double t2714, t2715, t2718, t2720, t2726, t2732, t2738, t2743;
909 double t2750, t2752, t2758, t2776, t2787, t2793, t2796, t2800;
910 double t2818, t2821, t2831, t2834, t2837, t2841, t2855, t2859;
911 double t2866, t2877, t2878, t2880, t2882, t2884, t2885, t2900;
912 double t2903, t2904, t2906, t2919, t2922, t2935, t2947, t2964;
913 double t2972, t2981, t2990, t3011, t3012, t3015, t3017, t3023;
914 double t3028, t3034, t3050, t3053, t3069, t3075, t3090, t3101;
915 double t3136, t3137, t3139, t3152, t3162, t3164, t3172, t3180;
916 double t3186, t3195, t3229, t3235, t3236, t3248, t3263, t3273;
917 double t3305, t3315, t3317, t3325, t3338, t3347, t3377, t3382;
918 double t3393, t3399, t3401, t3419, t3423, t3431, t3470, t3475;
919 double t3494, t3500, t3502, t3511, t3528, t3538, t3571, t3576;
920 double t3584, t3590, t3595, t3600, t3601, t3623, t3635, t3668;
921 double t3676, t3685, t3704, t3712, t3718, t3724, t3729, t3730;
922 double t3765;
923 #endif
924
925 #endif
926
927 #endif
928
929 #endif
930
931 #endif
932
933
934 gga_x_cap_params *params;
935
936 assert(p->params != NULL);
937 params = (gga_x_cap_params * )(p->params);
938
939 t1 = M_CBRT3;
940 t3 = POW_1_3(0.1e1 / M_PI);
941 t4 = t1 * t3;
942 t5 = M_CBRT4;
943 t6 = t5 * t5;
944 t7 = t4 * t6;
945 t8 = rho[0] - rho[1];
946 t9 = rho[0] + rho[1];
947 t10 = 0.1e1 / t9;
948 t11 = t8 * t10;
949 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
950 t14 = POW_1_3(t13);
951 t15 = t14 * t13;
952 t16 = POW_1_3(t9);
953 t17 = t15 * t16;
954 t18 = M_CBRT6;
955 t19 = t18 * t18;
956 t21 = M_PI * M_PI;
957 t22 = POW_1_3(t21);
958 t23 = 0.1e1 / t22;
959 t24 = params->alphaoAx * t19 * t23;
960 t25 = sqrt(sigma[0]);
961 t26 = POW_1_3(rho[0]);
962 t28 = 0.1e1 / t26 / rho[0];
963 t29 = t25 * t28;
964 t30 = t19 * t23;
965 t33 = 0.1e1 + t30 * t29 / 0.12e2;
966 t34 = log(t33);
967 t36 = params->c * t34 + 0.1e1;
968 t37 = 0.1e1 / t36;
969 t38 = t34 * t37;
970 t42 = 0.1e1 - t24 * t29 * t38 / 0.12e2;
971 t44 = t7 * t17 * t42;
972 t46 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
973 t47 = POW_1_3(t46);
974 t48 = t47 * t46;
975 t49 = t48 * t16;
976 t50 = sqrt(sigma[2]);
977 t51 = POW_1_3(rho[1]);
978 t53 = 0.1e1 / t51 / rho[1];
979 t54 = t50 * t53;
980 t57 = 0.1e1 + t30 * t54 / 0.12e2;
981 t58 = log(t57);
982 t60 = params->c * t58 + 0.1e1;
983 t61 = 0.1e1 / t60;
984 t62 = t58 * t61;
985 t66 = 0.1e1 - t24 * t54 * t62 / 0.12e2;
986 t68 = t7 * t49 * t66;
987 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
988 zk[0] = -0.3e1 / 0.8e1 * t44 - 0.3e1 / 0.8e1 * t68;
989
990 #ifndef XC_DONT_COMPILE_VXC
991
992 if(order < 1) return;
993
994
995 t70 = 0.3e1 / 0.8e1 * t44;
996 t71 = 0.3e1 / 0.8e1 * t68;
997 t72 = t14 * t16;
998 t73 = t9 * t9;
999 t74 = 0.1e1 / t73;
1000 t75 = t8 * t74;
1001 t77 = t10 / 0.2e1 - t75 / 0.2e1;
1002 t78 = t42 * t77;
1003 t80 = t7 * t72 * t78;
1004 t81 = t80 / 0.2e1;
1005 t82 = t16 * t16;
1006 t83 = 0.1e1 / t82;
1007 t84 = t15 * t83;
1008 t86 = t7 * t84 * t42;
1009 t87 = t86 / 0.8e1;
1010 t88 = rho[0] * rho[0];
1011 t90 = 0.1e1 / t26 / t88;
1012 t95 = params->alphaoAx * t18;
1013 t96 = t22 * t22;
1014 t97 = 0.1e1 / t96;
1015 t98 = t95 * t97;
1016 t99 = t88 * rho[0];
1017 t100 = t26 * t26;
1018 t102 = 0.1e1 / t100 / t99;
1019 t104 = 0.1e1 / t33;
1020 t105 = t104 * t37;
1021 t110 = t95 * t97 * sigma[0];
1022 t112 = t36 * t36;
1023 t113 = 0.1e1 / t112;
1024 t114 = t113 * params->c;
1025 t115 = t114 * t104;
1026 t116 = t102 * t34 * t115;
1027 t119 = t24 * t25 * t90 * t38 / 0.9e1 + t98 * sigma[0] * t102 * t105 / 0.18e2 - t110 * t116 / 0.18e2;
1028 t121 = t7 * t17 * t119;
1029 t122 = 0.3e1 / 0.8e1 * t121;
1030 t123 = t47 * t16;
1031 t124 = -t77;
1032 t125 = t66 * t124;
1033 t127 = t7 * t123 * t125;
1034 t128 = t127 / 0.2e1;
1035 t129 = t48 * t83;
1036 t131 = t7 * t129 * t66;
1037 t132 = t131 / 0.8e1;
1038 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1039 vrho[0] = -t70 - t71 + t9 * (-t81 - t87 - t122 - t128 - t132);
1040
1041 t136 = -t10 / 0.2e1 - t75 / 0.2e1;
1042 t137 = t42 * t136;
1043 t139 = t7 * t72 * t137;
1044 t140 = t139 / 0.2e1;
1045 t141 = -t136;
1046 t142 = t66 * t141;
1047 t144 = t7 * t123 * t142;
1048 t145 = t144 / 0.2e1;
1049 t146 = rho[1] * rho[1];
1050 t148 = 0.1e1 / t51 / t146;
1051 t153 = t146 * rho[1];
1052 t154 = t51 * t51;
1053 t156 = 0.1e1 / t154 / t153;
1054 t158 = 0.1e1 / t57;
1055 t159 = t158 * t61;
1056 t164 = t95 * t97 * sigma[2];
1057 t166 = t60 * t60;
1058 t167 = 0.1e1 / t166;
1059 t168 = t167 * params->c;
1060 t169 = t168 * t158;
1061 t170 = t156 * t58 * t169;
1062 t173 = t24 * t50 * t148 * t62 / 0.9e1 + t98 * sigma[2] * t156 * t159 / 0.18e2 - t164 * t170 / 0.18e2;
1063 t175 = t7 * t49 * t173;
1064 t176 = 0.3e1 / 0.8e1 * t175;
1065 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1066 vrho[1] = -t70 - t71 + t9 * (-t140 - t87 - t145 - t132 - t176);
1067
1068 t180 = t16 * t9 * t1;
1069 t181 = t180 * t3;
1070 t182 = t6 * t15;
1071 t183 = 0.1e1 / t25;
1072 t189 = 0.1e1 / t100 / t88;
1073 t196 = t34 * t113;
1074 t198 = t196 * params->c * t104;
1075 t201 = -t24 * t183 * t28 * t38 / 0.24e2 - t98 * t189 * t104 * t37 / 0.48e2 + t95 * t97 * t189 * t198 / 0.48e2;
1076 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1077 vsigma[0] = -0.3e1 / 0.8e1 * t181 * t182 * t201;
1078
1079 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1080 vsigma[1] = 0.0e0;
1081
1082 t205 = t6 * t48;
1083 t206 = 0.1e1 / t50;
1084 t212 = 0.1e1 / t154 / t146;
1085 t219 = t58 * t167;
1086 t221 = t219 * params->c * t158;
1087 t224 = -t24 * t206 * t53 * t62 / 0.24e2 - t98 * t212 * t158 * t61 / 0.48e2 + t95 * t97 * t212 * t221 / 0.48e2;
1088 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
1089 vsigma[2] = -0.3e1 / 0.8e1 * t181 * t205 * t224;
1090
1091 #ifndef XC_DONT_COMPILE_FXC
1092
1093 if(order < 2) return;
1094
1095
1096 t228 = t86 / 0.4e1;
1097 t230 = t131 / 0.4e1;
1098 t231 = t14 * t14;
1099 t232 = 0.1e1 / t231;
1100 t233 = t232 * t16;
1101 t234 = t77 * t77;
1102 t235 = t42 * t234;
1103 t237 = t7 * t233 * t235;
1104 t238 = t237 / 0.6e1;
1105 t239 = t14 * t83;
1106 t241 = t7 * t239 * t78;
1107 t242 = t241 / 0.3e1;
1108 t243 = t119 * t77;
1109 t245 = t7 * t72 * t243;
1110 t246 = t73 * t9;
1111 t247 = 0.1e1 / t246;
1112 t248 = t8 * t247;
1113 t249 = -t74 + t248;
1114 t250 = t42 * t249;
1115 t252 = t7 * t72 * t250;
1116 t253 = t252 / 0.2e1;
1117 t255 = 0.1e1 / t82 / t9;
1118 t256 = t15 * t255;
1119 t258 = t7 * t256 * t42;
1120 t259 = t258 / 0.12e2;
1121 t261 = t7 * t84 * t119;
1122 t262 = t261 / 0.4e1;
1123 t264 = 0.1e1 / t26 / t99;
1124 t269 = t88 * t88;
1125 t271 = 0.1e1 / t100 / t269;
1126 t277 = t271 * t34 * t115;
1127 t281 = params->alphaoAx / t21;
1128 t282 = t25 * sigma[0];
1129 t283 = t281 * t282;
1130 t284 = t269 * t88;
1131 t285 = 0.1e1 / t284;
1132 t286 = t33 * t33;
1133 t287 = 0.1e1 / t286;
1134 t288 = t285 * t287;
1135 t289 = t288 * t37;
1136 t292 = t288 * t114;
1137 t296 = t281 * t282 * t285;
1138 t298 = 0.1e1 / t112 / t36;
1139 t300 = params->c * params->c;
1140 t302 = t34 * t298 * t300 * t287;
1141 t306 = t196 * params->c * t287;
1142 t309 = -0.7e1 / 0.27e2 * t24 * t25 * t264 * t38 - 0.5e1 / 0.18e2 * t98 * sigma[0] * t271 * t105 + 0.5e1 / 0.18e2 * t110 * t277 + t283 * t289 / 0.27e2 + 0.2e1 / 0.27e2 * t283 * t292 - 0.2e1 / 0.27e2 * t296 * t302 - t296 * t306 / 0.27e2;
1143 t311 = t7 * t17 * t309;
1144 t312 = 0.3e1 / 0.8e1 * t311;
1145 t313 = t47 * t47;
1146 t314 = 0.1e1 / t313;
1147 t315 = t314 * t16;
1148 t316 = t124 * t124;
1149 t317 = t66 * t316;
1150 t319 = t7 * t315 * t317;
1151 t320 = t319 / 0.6e1;
1152 t321 = t47 * t83;
1153 t323 = t7 * t321 * t125;
1154 t324 = t323 / 0.3e1;
1155 t325 = -t249;
1156 t326 = t66 * t325;
1157 t328 = t7 * t123 * t326;
1158 t329 = t328 / 0.2e1;
1159 t330 = t48 * t255;
1160 t332 = t7 * t330 * t66;
1161 t333 = t332 / 0.12e2;
1162 t334 = -t238 - t242 - t245 - t253 + t259 - t262 - t312 - t320 - t324 - t329 + t333;
1163 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1164 v2rho2[0] = -t80 - t228 - 0.3e1 / 0.4e1 * t121 - t127 - t230 + t9 * t334;
1165
1166 t336 = t6 * t232;
1167 t337 = t4 * t336;
1168 t338 = t16 * t42;
1169 t339 = t136 * t77;
1170 t341 = t337 * t338 * t339;
1171 t344 = t7 * t239 * t137;
1172 t346 = t119 * t136;
1173 t348 = t7 * t72 * t346;
1174 t351 = 0.1e1 / t82 / t73;
1175 t352 = t14 * t351;
1176 t353 = t42 * t8;
1177 t355 = t7 * t352 * t353;
1178 t359 = t6 * t314;
1179 t360 = t4 * t359;
1180 t361 = t16 * t66;
1181 t362 = t141 * t124;
1182 t364 = t360 * t361 * t362;
1183 t367 = t7 * t321 * t142;
1184 t369 = t47 * t351;
1185 t370 = t66 * t8;
1186 t372 = t7 * t369 * t370;
1187 t375 = t173 * t124;
1188 t377 = t7 * t123 * t375;
1189 t380 = t7 * t129 * t173;
1190 t382 = -t341 / 0.6e1 - t344 / 0.6e1 - t348 / 0.2e1 - t355 / 0.2e1 - t241 / 0.6e1 + t259 - t261 / 0.8e1 - t364 / 0.6e1 - t367 / 0.6e1 + t372 / 0.2e1 - t323 / 0.6e1 + t333 - t377 / 0.2e1 - t380 / 0.8e1;
1191 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1192 v2rho2[1] = t9 * t382 - t122 - t128 - t140 - t145 - t176 - t228 - t230 - t81;
1193
1194 t385 = t136 * t136;
1195 t386 = t42 * t385;
1196 t388 = t7 * t233 * t386;
1197 t389 = t388 / 0.6e1;
1198 t390 = t344 / 0.3e1;
1199 t391 = t74 + t248;
1200 t392 = t42 * t391;
1201 t394 = t7 * t72 * t392;
1202 t395 = t394 / 0.2e1;
1203 t396 = t141 * t141;
1204 t397 = t66 * t396;
1205 t399 = t7 * t315 * t397;
1206 t400 = t399 / 0.6e1;
1207 t401 = t367 / 0.3e1;
1208 t402 = t173 * t141;
1209 t404 = t7 * t123 * t402;
1210 t405 = -t391;
1211 t406 = t66 * t405;
1212 t408 = t7 * t123 * t406;
1213 t409 = t408 / 0.2e1;
1214 t410 = t380 / 0.4e1;
1215 t412 = 0.1e1 / t51 / t153;
1216 t417 = t146 * t146;
1217 t419 = 0.1e1 / t154 / t417;
1218 t425 = t419 * t58 * t169;
1219 t428 = t50 * sigma[2];
1220 t429 = t281 * t428;
1221 t430 = t417 * t146;
1222 t431 = 0.1e1 / t430;
1223 t432 = t57 * t57;
1224 t433 = 0.1e1 / t432;
1225 t434 = t431 * t433;
1226 t435 = t434 * t61;
1227 t438 = t434 * t168;
1228 t442 = t281 * t428 * t431;
1229 t444 = 0.1e1 / t166 / t60;
1230 t447 = t58 * t444 * t300 * t433;
1231 t451 = t219 * params->c * t433;
1232 t454 = -0.7e1 / 0.27e2 * t24 * t50 * t412 * t62 - 0.5e1 / 0.18e2 * t98 * sigma[2] * t419 * t159 + 0.5e1 / 0.18e2 * t164 * t425 + t429 * t435 / 0.27e2 + 0.2e1 / 0.27e2 * t429 * t438 - 0.2e1 / 0.27e2 * t442 * t447 - t442 * t451 / 0.27e2;
1233 t456 = t7 * t49 * t454;
1234 t457 = 0.3e1 / 0.8e1 * t456;
1235 t458 = -t389 - t390 - t395 + t259 - t400 - t401 - t404 - t409 + t333 - t410 - t457;
1236 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1237 v2rho2[2] = -t139 - t228 - t144 - t230 - 0.3e1 / 0.4e1 * t175 + t9 * t458;
1238
1239 t461 = t7 * t17 * t201;
1240 t463 = t6 * t14;
1241 t464 = t201 * t77;
1242 t480 = t269 * rho[0];
1243 t481 = 0.1e1 / t480;
1244 t482 = t281 * t481;
1245 t483 = t287 * t37;
1246 t484 = t483 * t25;
1247 t487 = t287 * t113;
1248 t489 = t487 * params->c * t25;
1249 t493 = t281 * t481 * t34;
1250 t494 = t298 * t300;
1251 t496 = t494 * t287 * t25;
1252 t501 = t24 * t183 * t90 * t38 / 0.18e2 + t98 * t102 * t104 * t37 / 0.12e2 - t95 * t97 * t102 * t198 / 0.12e2 - t482 * t484 / 0.72e2 - t482 * t489 / 0.36e2 + t493 * t496 / 0.36e2 + t493 * t489 / 0.72e2;
1253 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1254 v2rhosigma[0] = -t461 / 0.2e1 - t181 * t463 * t464 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t182 * t501;
1255
1256 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1257 v2rhosigma[1] = 0.0e0;
1258
1259 t506 = t7 * t49 * t224;
1260 t507 = t6 * t47;
1261 t508 = t224 * t124;
1262 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1263 v2rhosigma[2] = -t181 * t507 * t508 / 0.2e1 - t506 / 0.2e1;
1264
1265 t513 = t463 * t201 * t136;
1266 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1267 v2rhosigma[3] = -t181 * t513 / 0.2e1 - t461 / 0.2e1;
1268
1269 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1270 v2rhosigma[4] = 0.0e0;
1271
1272 t518 = t507 * t224 * t141;
1273 t533 = t417 * rho[1];
1274 t534 = 0.1e1 / t533;
1275 t535 = t281 * t534;
1276 t536 = t433 * t61;
1277 t537 = t536 * t50;
1278 t540 = t433 * t167;
1279 t542 = t540 * params->c * t50;
1280 t546 = t281 * t534 * t58;
1281 t547 = t444 * t300;
1282 t549 = t547 * t433 * t50;
1283 t554 = t24 * t206 * t148 * t62 / 0.18e2 + t98 * t156 * t158 * t61 / 0.12e2 - t95 * t97 * t156 * t221 / 0.12e2 - t535 * t537 / 0.72e2 - t535 * t542 / 0.36e2 + t546 * t549 / 0.36e2 + t546 * t542 / 0.72e2;
1284 t555 = t205 * t554;
1285 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1286 v2rhosigma[5] = -t506 / 0.2e1 - t181 * t518 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t555;
1287
1288 t558 = 0.1e1 / t282;
1289 t563 = 0.1e1 / sigma[0];
1290 t569 = t95 * t97 * t563;
1291 t571 = t189 * t34 * t115;
1292 t574 = 0.1e1 / t269;
1293 t575 = t281 * t574;
1294 t580 = t487 * params->c * t183;
1295 t584 = t281 * t574 * t34;
1296 t591 = t24 * t558 * t28 * t38 / 0.48e2 - t98 * t563 * t189 * t105 / 0.96e2 + t569 * t571 / 0.96e2 + t575 * t483 * t183 / 0.192e3 + t575 * t580 / 0.96e2 - t584 * t494 * t287 * t183 / 0.96e2 - t584 * t580 / 0.192e3;
1297 t592 = t182 * t591;
1298 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1299 v2sigma2[0] = -0.3e1 / 0.8e1 * t181 * t592;
1300
1301 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1302 v2sigma2[1] = 0.0e0;
1303
1304 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1305 v2sigma2[2] = 0.0e0;
1306
1307 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1308 v2sigma2[3] = 0.0e0;
1309
1310 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1311 v2sigma2[4] = 0.0e0;
1312
1313 t595 = 0.1e1 / t428;
1314 t600 = 0.1e1 / sigma[2];
1315 t606 = t95 * t97 * t600;
1316 t608 = t212 * t58 * t169;
1317 t611 = 0.1e1 / t417;
1318 t612 = t281 * t611;
1319 t617 = t540 * params->c * t206;
1320 t621 = t281 * t611 * t58;
1321 t628 = t24 * t595 * t53 * t62 / 0.48e2 - t98 * t600 * t212 * t159 / 0.96e2 + t606 * t608 / 0.96e2 + t612 * t536 * t206 / 0.192e3 + t612 * t617 / 0.96e2 - t621 * t547 * t433 * t206 / 0.96e2 - t621 * t617 / 0.192e3;
1322 t629 = t205 * t628;
1323 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1324 v2sigma2[5] = -0.3e1 / 0.8e1 * t181 * t629;
1325
1326 #ifndef XC_DONT_COMPILE_KXC
1327
1328 if(order < 3) return;
1329
1330
1331 t636 = t258 / 0.4e1;
1332 t638 = t332 / 0.4e1;
1333 t642 = 0.1e1 / t231 / t13;
1334 t643 = t642 * t16;
1335 t644 = t234 * t77;
1336 t645 = t42 * t644;
1337 t647 = t7 * t643 * t645;
1338 t648 = t647 / 0.9e1;
1339 t649 = t232 * t83;
1340 t651 = t7 * t649 * t235;
1341 t653 = t77 * t249;
1342 t655 = t337 * t338 * t653;
1343 t656 = t655 / 0.2e1;
1344 t657 = t14 * t255;
1345 t659 = t7 * t657 * t78;
1346 t660 = t659 / 0.3e1;
1347 t662 = 0.1e1 / t313 / t46;
1348 t663 = t662 * t16;
1349 t664 = t316 * t124;
1350 t665 = t66 * t664;
1351 t667 = t7 * t663 * t665;
1352 t668 = t667 / 0.9e1;
1353 t669 = t314 * t83;
1354 t671 = t7 * t669 * t317;
1355 t673 = t124 * t325;
1356 t675 = t360 * t361 * t673;
1357 t676 = t675 / 0.2e1;
1358 t677 = t47 * t255;
1359 t679 = t7 * t677 * t125;
1360 t680 = t679 / 0.3e1;
1361 t681 = t119 * t234;
1362 t683 = t7 * t233 * t681;
1363 t684 = t683 / 0.2e1;
1364 t686 = t7 * t239 * t243;
1365 t689 = t7 * t239 * t250;
1366 t691 = t309 * t77;
1367 t693 = t7 * t72 * t691;
1368 t694 = 0.3e1 / 0.2e1 * t693;
1369 t695 = t119 * t249;
1370 t697 = t7 * t72 * t695;
1371 t698 = 0.3e1 / 0.2e1 * t697;
1372 t699 = t73 * t73;
1373 t700 = 0.1e1 / t699;
1374 t701 = t8 * t700;
1375 t703 = 0.3e1 * t247 - 0.3e1 * t701;
1376 t704 = t42 * t703;
1377 t706 = t7 * t72 * t704;
1378 t707 = t706 / 0.2e1;
1379 t708 = t15 * t351;
1380 t710 = t7 * t708 * t42;
1381 t711 = 0.5e1 / 0.36e2 * t710;
1382 t713 = t7 * t321 * t326;
1383 t715 = -t703;
1384 t716 = t66 * t715;
1385 t718 = t7 * t123 * t716;
1386 t719 = t718 / 0.2e1;
1387 t720 = t48 * t351;
1388 t722 = t7 * t720 * t66;
1389 t723 = 0.5e1 / 0.36e2 * t722;
1390 t725 = t7 * t256 * t119;
1391 t726 = t725 / 0.4e1;
1392 t728 = t7 * t84 * t309;
1393 t731 = 0.1e1 / t26 / t269;
1394 t737 = 0.1e1 / t100 / t480;
1395 t746 = t269 * t99;
1396 t747 = 0.1e1 / t746;
1397 t748 = t747 * t287;
1398 t756 = t281 * t282 * t747;
1399 t761 = sigma[0] * sigma[0];
1400 t762 = t269 * t269;
1401 t764 = 0.1e1 / t26 / t762;
1402 t765 = t761 * t764;
1403 t766 = t281 * t765;
1404 t768 = 0.1e1 / t286 / t33;
1405 t770 = t768 * t37 * t30;
1406 t775 = params->c * t19 * t23;
1407 t776 = t768 * t113 * t775;
1408 t781 = t300 * t19 * t23;
1409 t782 = t768 * t298 * t781;
1410 t786 = t281 * t765 * t34;
1411 t787 = t112 * t112;
1412 t788 = 0.1e1 / t787;
1413 t789 = t300 * params->c;
1414 t790 = t788 * t789;
1415 t793 = t790 * t768 * t19 * t23;
1416 t800 = 0.70e2 / 0.81e2 * t24 * t25 * t731 * t38 + 0.119e3 / 0.81e2 * t98 * sigma[0] * t737 * t105 - 0.119e3 / 0.81e2 * t110 * t737 * t34 * t115 - 0.11e2 / 0.27e2 * t283 * t748 * t37 - 0.22e2 / 0.27e2 * t283 * t748 * t114 + 0.22e2 / 0.27e2 * t756 * t302 + 0.11e2 / 0.27e2 * t756 * t306 + 0.2e1 / 0.243e3 * t766 * t770 + 0.2e1 / 0.81e2 * t766 * t776 + 0.2e1 / 0.81e2 * t766 * t782 - 0.2e1 / 0.81e2 * t786 * t793 - 0.2e1 / 0.81e2 * t786 * t782 - 0.2e1 / 0.243e3 * t786 * t776;
1417 t802 = t7 * t17 * t800;
1418 t803 = 0.3e1 / 0.8e1 * t802;
1419 t804 = -t689 / 0.2e1 - t694 - t698 - t707 - t711 - t713 / 0.2e1 - t719 - t723 + t726 - 0.3e1 / 0.8e1 * t728 - t803;
1420 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1421 v3rho3[0] = -t237 / 0.2e1 - t241 - t319 / 0.2e1 - t323 - 0.3e1 * t245 - 0.3e1 / 0.2e1 * t252 + t636 - 0.3e1 / 0.2e1 * t328 + t638 - 0.3e1 / 0.4e1 * t261 - 0.9e1 / 0.8e1 * t311 + t9 * (t648 - t651 / 0.6e1 - t656 + t660 + t668 - t671 / 0.6e1 - t676 + t680 - t684 - t686 + t804);
1422
1423 t809 = t341 / 0.3e1;
1424 t810 = t364 / 0.3e1;
1425 t813 = t6 * t642;
1426 t814 = t4 * t813;
1427 t815 = t136 * t234;
1428 t817 = t814 * t338 * t815;
1429 t819 = t83 * t42;
1430 t821 = t337 * t819 * t339;
1431 t822 = t821 / 0.9e1;
1432 t823 = t351 * t42;
1433 t824 = t8 * t77;
1434 t826 = t337 * t823 * t824;
1435 t828 = t6 * t662;
1436 t829 = t4 * t828;
1437 t830 = t141 * t316;
1438 t832 = t829 * t361 * t830;
1439 t834 = t83 * t66;
1440 t836 = t360 * t834 * t362;
1441 t837 = t836 / 0.9e1;
1442 t838 = t141 * t325;
1443 t840 = t360 * t361 * t838;
1444 t842 = t351 * t66;
1445 t843 = t8 * t124;
1446 t845 = t360 * t842 * t843;
1447 t847 = t16 * t119;
1448 t849 = t337 * t847 * t339;
1449 t851 = t136 * t249;
1450 t853 = t337 * t338 * t851;
1451 t858 = t7 * t352 * t42;
1452 t861 = t7 * t369 * t66;
1453 t864 = t7 * t330 * t173;
1454 t866 = t817 / 0.9e1 - t822 - t826 / 0.3e1 + t832 / 0.9e1 - t837 - t840 / 0.6e1 + t845 / 0.3e1 - t849 / 0.3e1 - t853 / 0.6e1 - t711 - t723 + t725 / 0.6e1 - t728 / 0.8e1 - t858 / 0.2e1 + t861 / 0.2e1 + t864 / 0.12e2;
1455 t875 = t7 * t657 * t137;
1456 t878 = 0.1e1 / t82 / t246;
1457 t879 = t14 * t878;
1458 t881 = t7 * t879 * t353;
1459 t883 = t119 * t8;
1460 t885 = t7 * t352 * t883;
1461 t887 = t7 * t677 * t142;
1462 t889 = t173 * t316;
1463 t891 = t7 * t315 * t889;
1464 t894 = t7 * t321 * t375;
1465 t895 = t894 / 0.3e1;
1466 t897 = t7 * t239 * t346;
1467 t898 = t897 / 0.3e1;
1468 t899 = t309 * t136;
1469 t901 = t7 * t72 * t899;
1470 t903 = t47 * t878;
1471 t905 = t7 * t903 * t370;
1472 t907 = t173 * t325;
1473 t909 = t7 * t123 * t907;
1474 t911 = -t651 / 0.18e2 + 0.2e1 / 0.9e1 * t659 - t671 / 0.18e2 + 0.2e1 / 0.9e1 * t679 - t686 / 0.3e1 - t689 / 0.6e1 - t713 / 0.6e1 + t875 / 0.9e1 + 0.7e1 / 0.6e1 * t881 - t885 + t887 / 0.9e1 - t891 / 0.6e1 - t895 - t898 - t901 / 0.2e1 - 0.7e1 / 0.6e1 * t905 - t909 / 0.2e1;
1475 t914 = -t253 + t636 - t329 + t638 - t348 + t372 - t377 - t410 - t261 / 0.2e1 - t312 + t9 * (t866 + t911);
1476 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1477 v3rho3[1] = -t238 - 0.2e1 / 0.3e1 * t241 - t320 - 0.2e1 / 0.3e1 * t323 - t809 - t390 - t355 - t810 - t401 - t245 + t914;
1478
1479 t919 = t136 * t8;
1480 t921 = t337 * t823 * t919;
1481 t923 = t141 * t8;
1482 t925 = t360 * t842 * t923;
1483 t927 = t16 * t173;
1484 t929 = t360 * t927 * t362;
1485 t931 = t405 * t124;
1486 t933 = t360 * t361 * t931;
1487 t935 = t396 * t124;
1488 t937 = t829 * t361 * t935;
1489 t939 = t385 * t77;
1490 t941 = t814 * t338 * t939;
1491 t943 = t391 * t77;
1492 t945 = t337 * t338 * t943;
1493 t948 = t7 * t129 * t454;
1494 t954 = -t921 / 0.3e1 + t925 / 0.3e1 - t929 / 0.3e1 - t933 / 0.6e1 + t937 / 0.9e1 + t941 / 0.9e1 - t945 / 0.6e1 - t822 - t837 - t948 / 0.8e1 - t711 - t723 + t725 / 0.12e2 + t864 / 0.6e1 + t659 / 0.9e1 + t679 / 0.9e1;
1495 t960 = t7 * t649 * t386;
1496 t963 = t7 * t239 * t392;
1497 t965 = t119 * t385;
1498 t967 = t7 * t233 * t965;
1499 t969 = t119 * t391;
1500 t971 = t7 * t72 * t969;
1501 t974 = -t247 - 0.3e1 * t701;
1502 t975 = t42 * t974;
1503 t977 = t7 * t72 * t975;
1504 t979 = t173 * t8;
1505 t981 = t7 * t369 * t979;
1506 t982 = -t974;
1507 t983 = t66 * t982;
1508 t985 = t7 * t123 * t983;
1509 t987 = t454 * t124;
1510 t989 = t7 * t123 * t987;
1511 t992 = t7 * t669 * t397;
1512 t995 = t7 * t321 * t402;
1513 t998 = t7 * t321 * t406;
1514 t1000 = 0.2e1 / 0.9e1 * t875 - t881 / 0.3e1 + 0.2e1 / 0.9e1 * t887 - t895 - t898 + t905 / 0.3e1 - t960 / 0.18e2 - t963 / 0.6e1 - t967 / 0.6e1 - t971 / 0.2e1 - t977 / 0.2e1 + t981 - t985 / 0.2e1 - t989 / 0.2e1 - t992 / 0.18e2 - t995 / 0.3e1 - t998 / 0.6e1;
1515 t1003 = t638 - t348 + t372 - t377 - t380 / 0.2e1 - t395 - t404 - t409 - t262 - t457 + t9 * (t954 + t1000);
1516 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1517 v3rho3[2] = -t242 - t324 - t809 - 0.2e1 / 0.3e1 * t344 - t355 - t810 - 0.2e1 / 0.3e1 * t367 - t389 - t400 + t636 + t1003;
1518
1519 t1011 = t875 / 0.3e1;
1520 t1012 = t887 / 0.3e1;
1521 t1015 = t864 / 0.4e1;
1522 t1018 = -0.3e1 * t247 - 0.3e1 * t701;
1523 t1019 = t42 * t1018;
1524 t1021 = t7 * t72 * t1019;
1525 t1022 = t1021 / 0.2e1;
1526 t1023 = t173 * t396;
1527 t1025 = t7 * t315 * t1023;
1528 t1026 = t1025 / 0.2e1;
1529 t1028 = t454 * t141;
1530 t1030 = t7 * t123 * t1028;
1531 t1031 = 0.3e1 / 0.2e1 * t1030;
1532 t1032 = t173 * t405;
1533 t1034 = t7 * t123 * t1032;
1534 t1035 = 0.3e1 / 0.2e1 * t1034;
1535 t1036 = -t1018;
1536 t1037 = t66 * t1036;
1537 t1039 = t7 * t123 * t1037;
1538 t1040 = t1039 / 0.2e1;
1539 t1042 = 0.1e1 / t51 / t417;
1540 t1048 = 0.1e1 / t154 / t533;
1541 t1057 = t417 * t153;
1542 t1058 = 0.1e1 / t1057;
1543 t1059 = t1058 * t433;
1544 t1067 = t281 * t428 * t1058;
1545 t1072 = sigma[2] * sigma[2];
1546 t1073 = t417 * t417;
1547 t1075 = 0.1e1 / t51 / t1073;
1548 t1076 = t1072 * t1075;
1549 t1077 = t281 * t1076;
1550 t1079 = 0.1e1 / t432 / t57;
1551 t1081 = t1079 * t61 * t30;
1552 t1085 = t1079 * t167 * t775;
1553 t1089 = t1079 * t444 * t781;
1554 t1093 = t281 * t1076 * t58;
1555 t1094 = t166 * t166;
1556 t1095 = 0.1e1 / t1094;
1557 t1096 = t1095 * t789;
1558 t1099 = t1096 * t1079 * t19 * t23;
1559 t1106 = 0.70e2 / 0.81e2 * t24 * t50 * t1042 * t62 + 0.119e3 / 0.81e2 * t98 * sigma[2] * t1048 * t159 - 0.119e3 / 0.81e2 * t164 * t1048 * t58 * t169 - 0.11e2 / 0.27e2 * t429 * t1059 * t61 - 0.22e2 / 0.27e2 * t429 * t1059 * t168 + 0.22e2 / 0.27e2 * t1067 * t447 + 0.11e2 / 0.27e2 * t1067 * t451 + 0.2e1 / 0.243e3 * t1077 * t1081 + 0.2e1 / 0.81e2 * t1077 * t1085 + 0.2e1 / 0.81e2 * t1077 * t1089 - 0.2e1 / 0.81e2 * t1093 * t1099 - 0.2e1 / 0.81e2 * t1093 * t1089 - 0.2e1 / 0.243e3 * t1093 * t1085;
1560 t1108 = t7 * t49 * t1106;
1561 t1109 = 0.3e1 / 0.8e1 * t1108;
1562 t1112 = t385 * t136;
1563 t1113 = t42 * t1112;
1564 t1115 = t7 * t643 * t1113;
1565 t1116 = t1115 / 0.9e1;
1566 t1117 = t136 * t391;
1567 t1119 = t337 * t338 * t1117;
1568 t1120 = t1119 / 0.2e1;
1569 t1121 = t396 * t141;
1570 t1122 = t66 * t1121;
1571 t1124 = t7 * t663 * t1122;
1572 t1125 = t1124 / 0.9e1;
1573 t1126 = t141 * t405;
1574 t1128 = t360 * t361 * t1126;
1575 t1129 = t1128 / 0.2e1;
1576 t1130 = -t1031 - t1035 - t1040 - t1109 - t992 / 0.6e1 - t995 - t998 / 0.2e1 + t1116 - t1120 + t1125 - t1129;
1577 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1578 v3rho3[3] = -t344 - t367 - t388 / 0.2e1 - t399 / 0.2e1 + t636 + t638 - 0.3e1 / 0.4e1 * t380 - 0.3e1 / 0.2e1 * t394 - 0.3e1 * t404 - 0.3e1 / 0.2e1 * t408 - 0.9e1 / 0.8e1 * t456 + t9 * (t1011 + t1012 - t960 / 0.6e1 - t963 / 0.2e1 - t711 - t723 + t1015 - 0.3e1 / 0.8e1 * t948 - t1022 - t1026 + t1130);
1579
1580 t1134 = t7 * t72 * t464;
1581 t1138 = t7 * t84 * t201 / 0.6e1;
1582 t1140 = t7 * t17 * t501;
1583 t1141 = t201 * t234;
1584 t1145 = t501 * t77;
1585 t1148 = t201 * t249;
1586 t1164 = t281 * t285;
1587 t1170 = t281 * t285 * t34;
1588 t1176 = 0.1e1 / t26 / t746;
1589 t1178 = t281 * t1176 * t768;
1590 t1180 = t37 * sigma[0] * t30;
1591 t1184 = t113 * sigma[0] * t775;
1592 t1188 = sigma[0] * t19 * t23;
1593 t1189 = t494 * t1188;
1594 t1192 = t1176 * t34;
1595 t1195 = t789 * t768;
1596 t1196 = t1195 * t1188;
1597 t1201 = t300 * t768;
1598 t1202 = t1201 * t1188;
1599 t1207 = params->c * t768;
1600 t1208 = t1207 * t1188;
1601 t1211 = -0.7e1 / 0.54e2 * t24 * t183 * t264 * t38 - 0.37e2 / 0.108e3 * t98 * t271 * t104 * t37 + 0.37e2 / 0.108e3 * t95 * t97 * t271 * t198 + t1164 * t484 / 0.8e1 + t1164 * t489 / 0.4e1 - t1170 * t496 / 0.4e1 - t1170 * t489 / 0.8e1 - t1178 * t1180 / 0.324e3 - t1178 * t1184 / 0.108e3 - t1178 * t1189 / 0.108e3 + t281 * t1192 * t788 * t1196 / 0.108e3 + t281 * t1192 * t298 * t1202 / 0.108e3 + t281 * t1192 * t113 * t1208 / 0.324e3;
1602 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1603 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t1134 - t1138 - t1140 - t181 * t336 * t1141 / 0.6e1 - t181 * t463 * t1145 - t181 * t463 * t1148 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t182 * t1211;
1604
1605 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1606 v3rho2sigma[1] = 0.0e0;
1607
1608 t1216 = t7 * t123 * t508;
1609 t1220 = t7 * t129 * t224 / 0.6e1;
1610 t1221 = t224 * t316;
1611 t1225 = t224 * t325;
1612 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1613 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t1216 - t1220 - t181 * t359 * t1221 / 0.6e1 - t181 * t507 * t1225 / 0.2e1;
1614
1615 t1231 = t16 * t1;
1616 t1232 = t1231 * t3;
1617 t1233 = t1232 * t513;
1618 t1235 = t3 * t6;
1619 t1236 = t180 * t1235;
1620 t1237 = t232 * t201;
1621 t1238 = t1237 * t339;
1622 t1242 = t463 * t501 * t136;
1623 t1245 = t255 * t1;
1624 t1246 = t1245 * t3;
1625 t1248 = t463 * t201 * t8;
1626 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1627 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t1134 - t1138 - t1140 / 0.2e1 - 0.2e1 / 0.3e1 * t1233 - t1236 * t1238 / 0.6e1 - t181 * t1242 / 0.2e1 - t1246 * t1248 / 0.2e1;
1628
1629 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1630 v3rho2sigma[4] = 0.0e0;
1631
1632 t1252 = t1232 * t518;
1633 t1254 = t314 * t224;
1634 t1255 = t1254 * t362;
1635 t1259 = t507 * t224 * t8;
1636 t1262 = t1232 * t555;
1637 t1265 = t507 * t554 * t124;
1638 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1639 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t1216 - t1220 - 0.2e1 / 0.3e1 * t1252 - t1236 * t1255 / 0.6e1 + t1246 * t1259 / 0.2e1 - t1262 / 0.2e1 - t181 * t1265 / 0.2e1;
1640
1641 t1270 = t336 * t201 * t385;
1642 t1274 = t463 * t201 * t391;
1643 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1644 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t1233 - t1138 - t181 * t1270 / 0.6e1 - t181 * t1274 / 0.2e1;
1645
1646 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1647 v3rho2sigma[7] = 0.0e0;
1648
1649 t1279 = t359 * t224 * t396;
1650 t1283 = t507 * t554 * t141;
1651 t1286 = t507 * t224 * t405;
1652 t1301 = t281 * t431;
1653 t1307 = t281 * t431 * t58;
1654 t1313 = 0.1e1 / t51 / t1057;
1655 t1315 = t281 * t1313 * t1079;
1656 t1317 = t61 * sigma[2] * t30;
1657 t1321 = t167 * sigma[2] * t775;
1658 t1325 = sigma[2] * t19 * t23;
1659 t1326 = t547 * t1325;
1660 t1329 = t1313 * t58;
1661 t1332 = t789 * t1079;
1662 t1333 = t1332 * t1325;
1663 t1338 = t300 * t1079;
1664 t1339 = t1338 * t1325;
1665 t1344 = params->c * t1079;
1666 t1345 = t1344 * t1325;
1667 t1348 = -0.7e1 / 0.54e2 * t24 * t206 * t412 * t62 - 0.37e2 / 0.108e3 * t98 * t419 * t158 * t61 + 0.37e2 / 0.108e3 * t95 * t97 * t419 * t221 + t1301 * t537 / 0.8e1 + t1301 * t542 / 0.4e1 - t1307 * t549 / 0.4e1 - t1307 * t542 / 0.8e1 - t1315 * t1317 / 0.324e3 - t1315 * t1321 / 0.108e3 - t1315 * t1326 / 0.108e3 + t281 * t1329 * t1095 * t1333 / 0.108e3 + t281 * t1329 * t444 * t1339 / 0.108e3 + t281 * t1329 * t167 * t1345 / 0.324e3;
1668 t1349 = t205 * t1348;
1669 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1670 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t1252 - t1220 - t1262 - t181 * t1279 / 0.6e1 - t181 * t1283 - t181 * t1286 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t1349;
1671
1672 t1352 = t1232 * t592;
1673 t1355 = t463 * t591 * t77;
1674 t1368 = t281 * t183;
1675 t1369 = t481 * t287;
1676 t1370 = t1369 * t37;
1677 t1373 = t1369 * t114;
1678 t1377 = t281 * t183 * t481;
1679 t1383 = 0.1e1 / t26 / t284;
1680 t1388 = t281 * t1383 * t768;
1681 t1389 = t114 * t30;
1682 t1392 = t494 * t30;
1683 t1395 = t1383 * t34;
1684 t1396 = t281 * t1395;
1685 t1403 = -t24 * t558 * t90 * t38 / 0.36e2 + t98 * t563 * t102 * t105 / 0.72e2 - t569 * t116 / 0.72e2 - t1368 * t1370 / 0.36e2 - t1368 * t1373 / 0.18e2 + t1377 * t302 / 0.18e2 + t1377 * t306 / 0.36e2 + t281 * t1383 * t770 / 0.864e3 + t1388 * t1389 / 0.288e3 + t1388 * t1392 / 0.288e3 - t1396 * t793 / 0.288e3 - t1396 * t782 / 0.288e3 - t1396 * t776 / 0.864e3;
1686 t1404 = t182 * t1403;
1687 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1688 v3rhosigma2[0] = -t1352 / 0.2e1 - t181 * t1355 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t1404;
1689
1690 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1691 v3rhosigma2[1] = 0.0e0;
1692
1693 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1694 v3rhosigma2[2] = 0.0e0;
1695
1696 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1697 v3rhosigma2[3] = 0.0e0;
1698
1699 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1700 v3rhosigma2[4] = 0.0e0;
1701
1702 t1407 = t1232 * t629;
1703 t1409 = t507 * t628 * t124;
1704 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1705 v3rhosigma2[5] = -t181 * t1409 / 0.2e1 - t1407 / 0.2e1;
1706
1707 t1413 = t463 * t591 * t136;
1708 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1709 v3rhosigma2[6] = -t181 * t1413 / 0.2e1 - t1352 / 0.2e1;
1710
1711 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1712 v3rhosigma2[7] = 0.0e0;
1713
1714 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1715 v3rhosigma2[8] = 0.0e0;
1716
1717 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1718 v3rhosigma2[9] = 0.0e0;
1719
1720 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1721 v3rhosigma2[10] = 0.0e0;
1722
1723 t1418 = t507 * t628 * t141;
1724 t1431 = t281 * t206;
1725 t1432 = t534 * t433;
1726 t1433 = t1432 * t61;
1727 t1436 = t1432 * t168;
1728 t1440 = t281 * t206 * t534;
1729 t1446 = 0.1e1 / t51 / t430;
1730 t1451 = t281 * t1446 * t1079;
1731 t1452 = t168 * t30;
1732 t1455 = t547 * t30;
1733 t1458 = t1446 * t58;
1734 t1459 = t281 * t1458;
1735 t1466 = -t24 * t595 * t148 * t62 / 0.36e2 + t98 * t600 * t156 * t159 / 0.72e2 - t606 * t170 / 0.72e2 - t1431 * t1433 / 0.36e2 - t1431 * t1436 / 0.18e2 + t1440 * t447 / 0.18e2 + t1440 * t451 / 0.36e2 + t281 * t1446 * t1081 / 0.864e3 + t1451 * t1452 / 0.288e3 + t1451 * t1455 / 0.288e3 - t1459 * t1099 / 0.288e3 - t1459 * t1089 / 0.288e3 - t1459 * t1085 / 0.864e3;
1736 t1467 = t205 * t1466;
1737 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1738 v3rhosigma2[11] = -t1407 / 0.2e1 - t181 * t1418 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t1467;
1739
1740 t1470 = t25 * t761;
1741 t1471 = 0.1e1 / t1470;
1742 t1476 = 0.1e1 / t761;
1743 t1482 = t95 * t97 * t1476;
1744 t1486 = 0.1e1 / t26 / t480;
1745 t1488 = t281 * t1486 * t768;
1746 t1490 = t37 * t563 * t30;
1747 t1494 = t113 * t563 * t775;
1748 t1498 = t563 * t19 * t23;
1749 t1499 = t494 * t1498;
1750 t1502 = t1486 * t34;
1751 t1504 = t281 * t1502 * t788;
1752 t1505 = t1195 * t1498;
1753 t1509 = t281 * t1502 * t298;
1754 t1510 = t1201 * t1498;
1755 t1514 = t281 * t1502 * t113;
1756 t1515 = t1207 * t1498;
1757 t1518 = -t24 * t1471 * t28 * t38 / 0.32e2 + t98 * t1476 * t189 * t105 / 0.64e2 - t1482 * t571 / 0.64e2 - t1488 * t1490 / 0.2304e4 - t1488 * t1494 / 0.768e3 - t1488 * t1499 / 0.768e3 + t1504 * t1505 / 0.768e3 + t1509 * t1510 / 0.768e3 + t1514 * t1515 / 0.2304e4;
1758 t1519 = t182 * t1518;
1759 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1760 v3sigma3[0] = -0.3e1 / 0.8e1 * t181 * t1519;
1761
1762 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1763 v3sigma3[1] = 0.0e0;
1764
1765 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1766 v3sigma3[2] = 0.0e0;
1767
1768 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1769 v3sigma3[3] = 0.0e0;
1770
1771 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1772 v3sigma3[4] = 0.0e0;
1773
1774 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1775 v3sigma3[5] = 0.0e0;
1776
1777 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1778 v3sigma3[6] = 0.0e0;
1779
1780 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1781 v3sigma3[7] = 0.0e0;
1782
1783 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1784 v3sigma3[8] = 0.0e0;
1785
1786 t1522 = t50 * t1072;
1787 t1523 = 0.1e1 / t1522;
1788 t1528 = 0.1e1 / t1072;
1789 t1534 = t95 * t97 * t1528;
1790 t1538 = 0.1e1 / t51 / t533;
1791 t1540 = t281 * t1538 * t1079;
1792 t1542 = t61 * t600 * t30;
1793 t1546 = t167 * t600 * t775;
1794 t1550 = t600 * t19 * t23;
1795 t1551 = t547 * t1550;
1796 t1554 = t1538 * t58;
1797 t1556 = t281 * t1554 * t1095;
1798 t1557 = t1332 * t1550;
1799 t1561 = t281 * t1554 * t444;
1800 t1562 = t1338 * t1550;
1801 t1566 = t281 * t1554 * t167;
1802 t1567 = t1344 * t1550;
1803 t1570 = -t24 * t1523 * t53 * t62 / 0.32e2 + t98 * t1528 * t212 * t159 / 0.64e2 - t1534 * t608 / 0.64e2 - t1540 * t1542 / 0.2304e4 - t1540 * t1546 / 0.768e3 - t1540 * t1551 / 0.768e3 + t1556 * t1557 / 0.768e3 + t1561 * t1562 / 0.768e3 + t1566 * t1567 / 0.2304e4;
1804 t1571 = t205 * t1570;
1805 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1806 v3sigma3[9] = -0.3e1 / 0.8e1 * t181 * t1571;
1807
1808 #ifndef XC_DONT_COMPILE_LXC
1809
1810 if(order < 4) return;
1811
1812
1813 t1585 = 0.4e1 / 0.9e1 * t647 - 0.2e1 / 0.3e1 * t651 - 0.2e1 * t655 + 0.4e1 / 0.3e1 * t659 + 0.4e1 / 0.9e1 * t667 - 0.2e1 / 0.3e1 * t671 - 0.2e1 * t675 + 0.4e1 / 0.3e1 * t679 - 0.2e1 * t683 - 0.4e1 * t686 - 0.2e1 * t689;
1814 t1589 = 0.5e1 / 0.9e1 * t710;
1815 t1592 = 0.5e1 / 0.9e1 * t722;
1816 t1596 = t7 * t321 * t716;
1817 t1598 = t699 * t9;
1818 t1600 = t8 / t1598;
1819 t1602 = 0.12e2 * t700 - 0.12e2 * t1600;
1820 t1608 = t7 * t352 * t78;
1821 t1610 = t46 * t46;
1822 t1612 = 0.1e1 / t313 / t1610;
1823 t1613 = t1612 * t16;
1824 t1614 = t316 * t316;
1825 t1619 = t662 * t83;
1826 t1621 = t7 * t1619 * t665;
1827 t1623 = t314 * t255;
1828 t1625 = t7 * t1623 * t317;
1829 t1627 = t325 * t325;
1830 t1633 = t7 * t369 * t125;
1831 t1636 = t7 * t677 * t326;
1832 t1642 = t7 * t239 * t691;
1833 t1645 = t7 * t239 * t695;
1834 t1648 = t7 * t239 * t704;
1835 t1667 = t13 * t13;
1836 t1669 = 0.1e1 / t231 / t1667;
1837 t1670 = t1669 * t16;
1838 t1671 = t234 * t234;
1839 t1676 = t642 * t83;
1840 t1678 = t7 * t1676 * t645;
1841 t1680 = -0.2e1 / 0.3e1 * t1596 - t7 * t123 * t66 * t1602 / 0.2e1 - 0.20e2 / 0.27e2 * t1608 - 0.5e1 / 0.27e2 * t7 * t1613 * t66 * t1614 + 0.4e1 / 0.27e2 * t1621 + 0.2e1 / 0.9e1 * t1625 - t7 * t315 * t66 * t1627 / 0.2e1 - 0.20e2 / 0.27e2 * t1633 + 0.2e1 / 0.3e1 * t1636 - t7 * t233 * t309 * t234 - 0.2e1 * t1642 - 0.2e1 * t1645 - 0.2e1 / 0.3e1 * t1648 - 0.2e1 * t7 * t72 * t800 * t77 - 0.3e1 * t7 * t72 * t309 * t249 - 0.2e1 * t7 * t72 * t119 * t703 + t7 * t72 * t42 * t1602 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1670 * t42 * t1671 + 0.4e1 / 0.27e2 * t1678;
1842 t1681 = t232 * t255;
1843 t1683 = t7 * t1681 * t235;
1844 t1685 = t249 * t249;
1845 t1695 = t7 * t649 * t681;
1846 t1698 = t7 * t657 * t243;
1847 t1701 = t7 * t657 * t250;
1848 t1708 = t360 * t834 * t673;
1849 t1723 = t337 * t819 * t653;
1850 t1729 = t7 * t708 * t119;
1851 t1732 = t7 * t256 * t309;
1852 t1735 = t7 * t84 * t800;
1853 t1737 = 0.1e1 / t762;
1854 t1738 = t1737 * t287;
1855 t1742 = t762 * rho[0];
1856 t1745 = t761 / t26 / t1742;
1857 t1747 = t281 * t1745 * t34;
1858 t1757 = t1470 / t100 / t762 / t88;
1859 t1759 = t281 * t1757 * t34;
1860 t1761 = 0.1e1 / t787 / t36;
1861 t1762 = t300 * t300;
1862 t1764 = t286 * t286;
1863 t1765 = 0.1e1 / t1764;
1864 t1767 = t1765 * t18 * t97;
1865 t1771 = t790 * t1767;
1866 t1774 = t494 * t1767;
1867 t1777 = t114 * t1767;
1868 t1781 = 0.1e1 / t100 / t284;
1869 t1790 = t1781 * t34;
1870 t1794 = t281 * t1745;
1871 t1799 = t281 * t1757;
1872 t1811 = t281 * t282 * t1737;
1873 t1819 = t18 * t97;
1874 t1823 = 0.721e3 / 0.81e2 * t110 * t1790 * t115 - 0.116e3 / 0.243e3 * t1794 * t776 - 0.116e3 / 0.243e3 * t1794 * t782 + 0.44e2 / 0.729e3 * t1799 * t1777 + 0.8e1 / 0.81e2 * t1799 * t1774 + 0.16e2 / 0.243e3 * t1799 * t1771 - 0.910e3 / 0.243e3 * t24 * t25 * t1486 * t38 - 0.1862e4 / 0.243e3 * t1811 * t302 - 0.931e3 / 0.243e3 * t1811 * t306 - 0.116e3 / 0.729e3 * t1794 * t770 + 0.4e1 / 0.243e3 * t1799 * t1765 * t37 * t1819;
1875 t1831 = 0.10e2 / 0.27e2 * t7 * t15 * t878 * t42;
1876 t1835 = 0.10e2 / 0.27e2 * t7 * t48 * t878 * t66;
1877 t1836 = 0.2e1 / 0.9e1 * t1683 - t7 * t233 * t42 * t1685 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t643 * t119 * t644 - 0.2e1 / 0.3e1 * t1695 + 0.4e1 / 0.3e1 * t1698 + 0.2e1 / 0.3e1 * t1701 + 0.2e1 / 0.3e1 * t829 * t361 * t316 * t325 - 0.2e1 / 0.3e1 * t1708 - 0.2e1 / 0.3e1 * t360 * t361 * t124 * t715 - 0.2e1 / 0.3e1 * t337 * t338 * t77 * t703 + 0.2e1 / 0.3e1 * t814 * t338 * t234 * t249 - 0.2e1 / 0.3e1 * t1723 - 0.2e1 * t337 * t847 * t653 - 0.5e1 / 0.9e1 * t1729 + t1732 / 0.2e1 - t1735 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * (0.931e3 / 0.243e3 * t283 * t1738 * t37 + 0.116e3 / 0.243e3 * t1747 * t793 + 0.116e3 / 0.243e3 * t1747 * t782 + 0.116e3 / 0.729e3 * t1747 * t776 - 0.16e2 / 0.243e3 * t1759 * t1761 * t1762 * t1767 - 0.8e1 / 0.81e2 * t1759 * t1771 - 0.44e2 / 0.729e3 * t1759 * t1774 - 0.4e1 / 0.243e3 * t1759 * t1777 - 0.721e3 / 0.81e2 * t98 * sigma[0] * t1781 * t105 + 0.1862e4 / 0.243e3 * t283 * t1738 * t114 + t1823) + t1831 + t1835;
1878 t1839 = -0.6e1 * t693 - 0.6e1 * t697 - 0.2e1 * t706 - t1589 - 0.2e1 * t713 - 0.2e1 * t718 - t1592 + t725 - 0.3e1 / 0.2e1 * t728 - 0.3e1 / 0.2e1 * t802 + t9 * (t1680 + t1836);
1879 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1880 v4rho4[0] = t1585 + t1839;
1881
1882 t1856 = t337 * t255 * t42 * t339;
1883 t1857 = t1856 / 0.9e1;
1884 t1858 = t878 * t42;
1885 t1860 = t337 * t1858 * t824;
1886 t1872 = t360 * t255 * t66 * t362;
1887 t1873 = t1872 / 0.9e1;
1888 t1883 = t7 * t677 * t375;
1889 t1884 = t1883 / 0.3e1;
1890 t1886 = t7 * t352 * t137;
1891 t1889 = 0.1e1 / t82 / t699;
1892 t1892 = t7 * t14 * t1889 * t353;
1893 t1894 = t314 * t351;
1894 t1899 = t7 * t369 * t142;
1895 t1901 = -t360 * t927 * t673 / 0.2e1 - t337 * t16 * t309 * t339 / 0.2e1 - t337 * t847 * t851 / 0.2e1 - t337 * t338 * t136 * t703 / 0.6e1 + t1857 + 0.7e1 / 0.6e1 * t1860 - t360 * t361 * t141 * t715 / 0.6e1 + t360 * t842 * t8 * t325 / 0.2e1 + t1873 + t814 * t338 * t851 * t77 / 0.3e1 + t829 * t361 * t838 * t124 / 0.3e1 + t1884 - 0.5e1 / 0.27e2 * t1886 - 0.25e2 / 0.6e1 * t1892 + t7 * t1894 * t125 / 0.2e1 - 0.5e1 / 0.27e2 * t1899;
1896 t1907 = t7 * t669 * t889;
1897 t1909 = t232 * t351;
1898 t1918 = t7 * t657 * t346;
1899 t1919 = t1918 / 0.3e1;
1900 t1921 = t7 * t879 * t883;
1901 t1929 = t7 * t47 * t1889 * t370;
1902 t1932 = t7 * t321 * t907;
1903 t1935 = t7 * t239 * t899;
1904 t1946 = t7 * t663 * t173 * t664 / 0.9e1 - t1907 / 0.6e1 - t7 * t1909 * t78 / 0.2e1 - t7 * t123 * t173 * t715 / 0.2e1 + t1919 + 0.7e1 / 0.2e1 * t1921 - 0.3e1 / 0.2e1 * t7 * t352 * t309 * t8 + 0.25e2 / 0.6e1 * t1929 - t1932 / 0.2e1 - t1935 / 0.2e1 - t7 * t72 * t800 * t136 / 0.2e1 - t1596 / 0.6e1 - 0.5e1 / 0.9e1 * t1608 + t1621 / 0.27e2 + t1625 / 0.9e1 - 0.5e1 / 0.9e1 * t1633;
1905 t1958 = t337 * t819 * t851;
1906 t1965 = t360 * t834 * t838;
1907 t1971 = t878 * t66;
1908 t1973 = t360 * t1971 * t843;
1909 t1976 = t4 * t6 * t1669;
1910 t1982 = t814 * t819 * t815;
1911 t1984 = t1636 / 0.3e1 - t1642 / 0.2e1 - t1645 / 0.2e1 - t1648 / 0.6e1 + t1678 / 0.27e2 + t1683 / 0.9e1 - t1695 / 0.6e1 + 0.2e1 / 0.3e1 * t1698 + t1701 / 0.3e1 - t1958 / 0.6e1 - t337 * t823 * t8 * t249 / 0.2e1 - t1965 / 0.6e1 - t829 * t842 * t8 * t316 / 0.3e1 - 0.7e1 / 0.6e1 * t1973 - 0.5e1 / 0.27e2 * t1976 * t338 * t136 * t644 + t1982 / 0.9e1;
1912 t1990 = t4 * t6 * t1612;
1913 t1996 = t829 * t834 * t830;
1914 t2003 = t337 * t83 * t119 * t339;
1915 t2005 = t351 * t119;
1916 t2014 = t7 * t879 * t42;
1917 t2020 = t7 * t903 * t66;
1918 t2023 = t7 * t720 * t173;
1919 t2025 = t814 * t823 * t8 * t234 / 0.3e1 - 0.5e1 / 0.27e2 * t1990 * t361 * t141 * t664 + t1996 / 0.9e1 + t814 * t847 * t815 / 0.3e1 - t2003 / 0.3e1 - t337 * t2005 * t824 - t1708 / 0.6e1 - t1723 / 0.6e1 - 0.5e1 / 0.12e2 * t1729 + t1732 / 0.4e1 - t1735 / 0.8e1 + t1831 + t1835 + 0.5e1 / 0.2e1 * t2014 - 0.3e1 / 0.2e1 * t7 * t352 * t119 - 0.5e1 / 0.2e1 * t2020 - 0.5e1 / 0.36e2 * t2023;
1920 t2030 = t821 / 0.3e1;
1921 t2032 = t836 / 0.3e1;
1922 t2035 = t9 * (t1901 + t1946 + t1984 + t2025) + t817 / 0.3e1 - t2030 - t826 + t832 / 0.3e1 - t2032 - t840 / 0.2e1 + t845 - t849 - t853 / 0.2e1 - t656;
1923 t2041 = -t676 - t1589 - t1592 + 0.3e1 / 0.4e1 * t725 - 0.3e1 / 0.4e1 * t728 - t803 - 0.3e1 / 0.2e1 * t858 + 0.3e1 / 0.2e1 * t861 + t1015 + t648 - t651 / 0.3e1;
1924 t2045 = t659 + t668 - t671 / 0.3e1 + t679 - t684 - 0.2e1 * t686 - t689 - t694 - t698 - t707 - t713;
1925 t2052 = -t719 + t1011 + 0.7e1 / 0.2e1 * t881 - 0.3e1 * t885 + t1012 - t891 / 0.2e1 - t894 - t897 - 0.3e1 / 0.2e1 * t901 - 0.7e1 / 0.2e1 * t905 - 0.3e1 / 0.2e1 * t909;
1926 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1927 v4rho4[1] = t2035 + t2041 + t2045 + t2052;
1928
1929 t2059 = t7 * t321 * t987;
1930 t2065 = t7 * t1623 * t397;
1931 t2071 = t7 * t677 * t402;
1932 t2074 = t7 * t677 * t406;
1933 t2077 = t7 * t1681 * t386;
1934 t2080 = t7 * t657 * t392;
1935 t2083 = t7 * t903 * t979;
1936 t2086 = t7 * t321 * t983;
1937 t2089 = 0.1e1 / t82 / t1598;
1938 t2091 = t8 * t8;
1939 t2097 = t7 * t239 * t969;
1940 t2100 = t7 * t239 * t975;
1941 t2123 = t7 * t649 * t965;
1942 t2125 = -t7 * t315 * t454 * t316 / 0.6e1 - t2059 / 0.3e1 - t7 * t1909 * t137 / 0.3e1 + t2065 / 0.27e2 + t7 * t1894 * t142 / 0.3e1 + 0.2e1 / 0.9e1 * t2071 + t2074 / 0.9e1 + t2077 / 0.27e2 + t2080 / 0.9e1 - 0.7e1 / 0.3e1 * t2083 - t2086 / 0.3e1 - t7 * t232 * t2089 * t42 * t2091 / 0.3e1 - t2097 / 0.3e1 - t2100 / 0.3e1 - t7 * t233 * t309 * t385 / 0.6e1 - t7 * t72 * t309 * t391 / 0.2e1 - t7 * t72 * t119 * t974 - t7 * t123 * t454 * t325 / 0.2e1 - t7 * t314 * t2089 * t66 * t2091 / 0.3e1 - t2123 / 0.9e1;
1943 t2137 = t351 * t173;
1944 t2146 = t360 * t1971 * t923;
1945 t2171 = t829 * t834 * t935;
1946 t2181 = t360 * t83 * t173 * t362;
1947 t2188 = t360 * t834 * t931;
1948 t2195 = t814 * t819 * t939;
1949 t2197 = -t360 * t927 * t838 / 0.3e1 - t360 * t361 * t405 * t325 / 0.6e1 + t829 * t361 * t396 * t325 / 0.9e1 + 0.2e1 / 0.3e1 * t360 * t2137 * t843 - t360 * t361 * t982 * t124 / 0.3e1 - 0.7e1 / 0.9e1 * t2146 + 0.2e1 / 0.9e1 * t814 * t847 * t939 + t814 * t338 * t385 * t249 / 0.9e1 - t337 * t847 * t943 / 0.3e1 - t337 * t338 * t974 * t77 / 0.3e1 - t337 * t338 * t391 * t249 / 0.6e1 - 0.5e1 / 0.27e2 * t1990 * t361 * t396 * t316 + 0.2e1 / 0.27e2 * t2171 - 0.2e1 / 0.3e1 * t337 * t2005 * t919 + 0.2e1 / 0.9e1 * t829 * t927 * t830 - 0.2e1 / 0.9e1 * t2181 + t829 * t361 * t405 * t316 / 0.9e1 - t2188 / 0.9e1 - 0.5e1 / 0.27e2 * t1976 * t338 * t385 * t234 + 0.2e1 / 0.27e2 * t2195;
1950 t2200 = t337 * t1858 * t919;
1951 t2207 = t337 * t819 * t943;
1952 t2232 = 0.7e1 / 0.9e1 * t2200 + t814 * t338 * t391 * t234 / 0.9e1 - t2207 / 0.9e1 + 0.4e1 / 0.27e2 * t1856 - 0.2e1 / 0.9e1 * t1860 + 0.4e1 / 0.27e2 * t1872 - 0.4e1 / 0.9e1 * t829 * t842 * t362 * t8 + 0.4e1 / 0.9e1 * t814 * t823 * t339 * t8 + 0.4e1 / 0.9e1 * t1883 - 0.10e2 / 0.27e2 * t1886 - 0.41e2 / 0.9e1 * t1892 - 0.10e2 / 0.27e2 * t1899 - t1907 / 0.9e1 + 0.4e1 / 0.9e1 * t1918 - 0.2e1 / 0.3e1 * t1921 + 0.41e2 / 0.9e1 * t1929 - t1932 / 0.3e1 - t1935 / 0.3e1 - 0.10e2 / 0.27e2 * t1608 + t1625 / 0.27e2;
1953 t2252 = t7 * t330 * t454;
1954 t2254 = -0.10e2 / 0.27e2 * t1633 + t1636 / 0.9e1 + t1683 / 0.27e2 + 0.2e1 / 0.9e1 * t1698 + t1701 / 0.9e1 - t1958 / 0.9e1 - t1965 / 0.9e1 + 0.2e1 / 0.9e1 * t1973 + 0.2e1 / 0.27e2 * t1982 + 0.2e1 / 0.27e2 * t1996 - 0.2e1 / 0.9e1 * t2003 - 0.5e1 / 0.18e2 * t1729 + t1732 / 0.12e2 + t1831 + t1835 - t2014 / 0.3e1 + t2020 / 0.3e1 - 0.5e1 / 0.18e2 * t2023 + t7 * t369 * t173 + t2252 / 0.12e2;
1955 t2270 = t9 * (t2125 + t2197 + t2232 + t2254) - 0.2e1 / 0.3e1 * t921 + 0.2e1 / 0.3e1 * t925 - 0.2e1 / 0.3e1 * t929 - t933 / 0.3e1 + 0.2e1 / 0.9e1 * t937 + 0.2e1 / 0.9e1 * t941 - t945 / 0.3e1 + 0.2e1 / 0.9e1 * t817 - 0.4e1 / 0.9e1 * t821 - 0.2e1 / 0.3e1 * t826 + 0.2e1 / 0.9e1 * t832 - 0.4e1 / 0.9e1 * t836;
1956 t2280 = -t840 / 0.3e1 + 0.2e1 / 0.3e1 * t845 - 0.2e1 / 0.3e1 * t849 - t853 / 0.3e1 - t948 / 0.4e1 - t1589 - t1592 + t725 / 0.2e1 - t728 / 0.4e1 - t858 + t861 + t864 / 0.2e1 - t651 / 0.9e1;
1957 t2295 = 0.2e1 / 0.3e1 * t659 - t671 / 0.9e1 + 0.2e1 / 0.3e1 * t679 - 0.2e1 / 0.3e1 * t686 - t689 / 0.3e1 - t713 / 0.3e1 + 0.2e1 / 0.3e1 * t875 + 0.5e1 / 0.3e1 * t881 - 0.2e1 * t885 + 0.2e1 / 0.3e1 * t887 - t891 / 0.3e1 - 0.4e1 / 0.3e1 * t894 - 0.4e1 / 0.3e1 * t897;
1958 t2304 = -t901 - 0.5e1 / 0.3e1 * t905 - t909 - t960 / 0.9e1 - t963 / 0.3e1 - t967 / 0.3e1 - t971 - t977 + 0.2e1 * t981 - t985 - t989 - t992 / 0.9e1 - 0.2e1 / 0.3e1 * t995 - t998 / 0.3e1;
1959 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1960 v4rho4[2] = t2270 + t2280 + t2295 + t2304;
1961
1962 t2358 = t337 * t819 * t1117;
1963 t2364 = -0.5e1 / 0.27e2 * t1976 * t338 * t1112 * t77 - t337 * t847 * t1117 / 0.2e1 - t337 * t338 * t136 * t974 / 0.2e1 - 0.5e1 / 0.27e2 * t1990 * t361 * t1121 * t124 - t829 * t842 * t396 * t8 / 0.3e1 - t360 * t361 * t141 * t982 / 0.2e1 + t360 * t842 * t8 * t405 / 0.2e1 + t829 * t927 * t935 / 0.3e1 + t360 * t2137 * t923 - t360 * t16 * t454 * t362 / 0.2e1 - t360 * t927 * t931 / 0.2e1 - t360 * t361 * t1036 * t124 / 0.6e1 - t337 * t338 * t1018 * t77 / 0.6e1 + t814 * t823 * t385 * t8 / 0.3e1 - t2358 / 0.6e1 - t337 * t823 * t8 * t391 / 0.2e1;
1964 t2366 = t360 * t834 * t1126;
1965 t2369 = t7 * t1619 * t1122;
1966 t2372 = t7 * t669 * t1023;
1967 t2375 = t7 * t321 * t1028;
1968 t2378 = t7 * t321 * t1032;
1969 t2381 = t7 * t321 * t1037;
1970 t2384 = t7 * t1676 * t1113;
1971 t2387 = t7 * t239 * t1019;
1972 t2391 = -0.6e1 * t700 - 0.12e2 * t1600;
1973 t2423 = -t2366 / 0.6e1 + t2369 / 0.27e2 - t2372 / 0.6e1 - t2375 / 0.2e1 - t2378 / 0.2e1 - t2381 / 0.6e1 + t2384 / 0.27e2 - t2387 / 0.6e1 - t7 * t123 * t66 * t2391 / 0.2e1 - t7 * t123 * t1106 * t124 / 0.2e1 + t7 * t643 * t119 * t1112 / 0.9e1 - t7 * t72 * t119 * t1018 / 0.2e1 + t7 * t72 * t42 * t2391 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t369 * t454 * t8 - 0.3e1 / 0.2e1 * t7 * t123 * t173 * t982 - t2059 / 0.2e1 + t2065 / 0.9e1;
1974 t2440 = 0.2e1 / 0.3e1 * t2071 + t2074 / 0.3e1 + t2077 / 0.9e1 + t2080 / 0.3e1 + t2083 - t2086 / 0.2e1 - t2097 / 0.2e1 - t2100 / 0.2e1 - t2123 / 0.6e1 + t2146 / 0.3e1 + t2171 / 0.9e1 - t2181 / 0.3e1 - t2188 / 0.6e1 + t2195 / 0.9e1 - t2200 / 0.3e1 - t2207 / 0.6e1 + t1857;
1975 t2459 = t7 * t129 * t1106;
1976 t2461 = t1873 + t829 * t361 * t1126 * t124 / 0.3e1 + t814 * t338 * t1117 * t77 / 0.3e1 + t1884 - 0.5e1 / 0.9e1 * t1886 + t1892 / 0.3e1 - 0.5e1 / 0.9e1 * t1899 + t1919 - t1929 / 0.3e1 - 0.5e1 / 0.27e2 * t1608 - 0.5e1 / 0.27e2 * t1633 - 0.5e1 / 0.36e2 * t1729 + t1831 + t1835 - 0.5e1 / 0.12e2 * t2023 + t2252 / 0.4e1 - t2459 / 0.8e1;
1977 t2469 = t9 * (t2364 + t2423 + t2440 + t2461) - t1120 - t1129 - t921 + t925 - t929 - t933 / 0.2e1 + t937 / 0.3e1 + t941 / 0.3e1 - t945 / 0.2e1 - t2030;
1978 t2472 = -t2032 - 0.3e1 / 0.4e1 * t948 - t1109 - t1589 - t1592 + t726 + 0.3e1 / 0.4e1 * t864 + t660 + t680 + t875 - t881;
1979 t2480 = t887 - t894 - t897 + t905 - t960 / 0.3e1 - t963 - t967 / 0.2e1 - 0.3e1 / 0.2e1 * t971 - 0.3e1 / 0.2e1 * t977 + 0.3e1 * t981 - 0.3e1 / 0.2e1 * t985;
1980 t2484 = -0.3e1 / 0.2e1 * t989 - t992 / 0.3e1 - 0.2e1 * t995 - t998 - t1022 - t1026 - t1031 - t1035 - t1040 + t1116 + t1125;
1981 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1982 v4rho4[3] = t2469 + t2472 + t2480 + t2484;
1983
1984 t2494 = 0.4e1 / 0.3e1 * t875 + 0.4e1 / 0.3e1 * t887 - 0.2e1 / 0.3e1 * t960 - 0.2e1 * t963 - t1589 - t1592 + t864 - 0.3e1 / 0.2e1 * t948 - 0.2e1 * t1021 - 0.2e1 * t1025 - 0.6e1 * t1030;
1985 t2518 = -0.12e2 * t700 - 0.12e2 * t1600;
1986 t2523 = t385 * t385;
1987 t2528 = t391 * t391;
1988 t2533 = t396 * t396;
1989 t2538 = t405 * t405;
1990 t2577 = -0.2e1 * t7 * t123 * t1106 * t141 - 0.3e1 * t7 * t123 * t454 * t405 - 0.2e1 * t7 * t123 * t173 * t1036 - t7 * t123 * t66 * t2518 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1670 * t42 * t2523 - t7 * t233 * t42 * t2528 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1613 * t66 * t2533 - t7 * t315 * t66 * t2538 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t663 * t173 * t1121 + t7 * t72 * t42 * t2518 / 0.2e1 - t7 * t315 * t454 * t396 + 0.2e1 / 0.3e1 * t829 * t361 * t396 * t405 + 0.2e1 / 0.3e1 * t814 * t338 * t385 * t391 - 0.2e1 / 0.3e1 * t337 * t338 * t1018 * t136 - 0.2e1 * t360 * t927 * t1126 - 0.2e1 / 0.3e1 * t360 * t361 * t1036 * t141 - 0.2e1 / 0.3e1 * t2358 - 0.2e1 / 0.3e1 * t2366 + 0.4e1 / 0.27e2 * t2369;
1991 t2594 = 0.1e1 / t1073;
1992 t2595 = t2594 * t433;
1993 t2599 = t1073 * rho[1];
1994 t2602 = t1072 / t51 / t2599;
1995 t2604 = t281 * t2602 * t58;
1996 t2614 = t1522 / t154 / t1073 / t146;
1997 t2616 = t281 * t2614 * t58;
1998 t2618 = 0.1e1 / t1094 / t60;
1999 t2620 = t432 * t432;
2000 t2621 = 0.1e1 / t2620;
2001 t2623 = t2621 * t18 * t97;
2002 t2627 = t1096 * t2623;
2003 t2630 = t547 * t2623;
2004 t2633 = t168 * t2623;
2005 t2637 = 0.1e1 / t154 / t430;
2006 t2638 = t2637 * t58;
2007 t2642 = t281 * t2602;
2008 t2648 = t281 * t2614;
2009 t2660 = t281 * t428 * t2594;
2010 t2678 = -0.116e3 / 0.243e3 * t2642 * t1089 + 0.44e2 / 0.729e3 * t2648 * t2633 + 0.8e1 / 0.81e2 * t2648 * t2630 + 0.16e2 / 0.243e3 * t2648 * t2627 - 0.910e3 / 0.243e3 * t24 * t50 * t1538 * t62 - 0.1862e4 / 0.243e3 * t2660 * t447 - 0.931e3 / 0.243e3 * t2660 * t451 - 0.116e3 / 0.729e3 * t2642 * t1081 + 0.4e1 / 0.243e3 * t2648 * t2621 * t61 * t1819 - 0.721e3 / 0.81e2 * t98 * sigma[2] * t2637 * t159 + 0.1862e4 / 0.243e3 * t429 * t2595 * t168;
2011 t2683 = -0.2e1 / 0.3e1 * t2372 - 0.2e1 * t2375 - 0.2e1 * t2378 - 0.2e1 / 0.3e1 * t2381 + 0.4e1 / 0.27e2 * t2384 - 0.2e1 / 0.3e1 * t2387 + 0.2e1 / 0.9e1 * t2065 + 0.4e1 / 0.3e1 * t2071 + 0.2e1 / 0.3e1 * t2074 + 0.2e1 / 0.9e1 * t2077 + 0.2e1 / 0.3e1 * t2080 - 0.20e2 / 0.27e2 * t1886 - 0.20e2 / 0.27e2 * t1899 + t1831 + t1835 - 0.5e1 / 0.9e1 * t2023 + t2252 / 0.2e1 - t2459 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t49 * (0.931e3 / 0.243e3 * t429 * t2595 * t61 + 0.116e3 / 0.243e3 * t2604 * t1099 + 0.116e3 / 0.243e3 * t2604 * t1089 + 0.116e3 / 0.729e3 * t2604 * t1085 - 0.16e2 / 0.243e3 * t2616 * t2618 * t1762 * t2623 - 0.8e1 / 0.81e2 * t2616 * t2627 - 0.44e2 / 0.729e3 * t2616 * t2630 - 0.4e1 / 0.243e3 * t2616 * t2633 + 0.721e3 / 0.81e2 * t164 * t2638 * t169 - 0.116e3 / 0.243e3 * t2642 * t1085 + t2678);
2012 t2686 = -0.6e1 * t1034 - 0.2e1 * t1039 - 0.3e1 / 0.2e1 * t1108 - 0.2e1 / 0.3e1 * t992 - 0.4e1 * t995 - 0.2e1 * t998 + 0.4e1 / 0.9e1 * t1115 - 0.2e1 * t1119 + 0.4e1 / 0.9e1 * t1124 - 0.2e1 * t1128 + t9 * (t2577 + t2683);
2013 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2014 v4rho4[4] = t2494 + t2686;
2015
2016 t2688 = t7 * t84 * t501;
2017 t2691 = t7 * t17 * t1211;
2018 t2697 = t281 * t747;
2019 t2700 = t764 * t34;
2020 t2714 = 0.1e1 / t100 / t1742;
2021 t2715 = t2714 * t34;
2022 t2718 = t1762 * t1765;
2023 t2720 = t282 * t18 * t97;
2024 t2726 = t789 * t1765;
2025 t2732 = t300 * t1765;
2026 t2738 = params->c * t1765;
2027 t2743 = t281 * t764 * t768;
2028 t2750 = t281 * t2714 * t1765;
2029 t2752 = t1819 * params->c;
2030 t2758 = t300 * t18 * t97;
2031 t2776 = t281 * t747 * t34;
2032 t2787 = 0.49e2 / 0.324e3 * t2743 * t1189 - 0.11e2 / 0.486e3 * t2750 * t113 * t282 * t2752 - t2750 * t298 * t282 * t2758 / 0.27e2 - 0.2e1 / 0.81e2 * t2750 * t790 * t2720 - 0.317e3 / 0.162e3 * t2697 * t489 + 0.35e2 / 0.81e2 * t24 * t183 * t731 * t38 - 0.91e2 / 0.54e2 * t95 * t97 * t737 * t198 + 0.317e3 / 0.162e3 * t2776 * t496 + 0.317e3 / 0.324e3 * t2776 * t489 + 0.49e2 / 0.972e3 * t2743 * t1180 - t2750 * t37 * t282 * t1819 / 0.162e3;
2033 t2793 = t7 * t72 * t1145;
2034 t2796 = t7 * t72 * t1148;
2035 t2800 = t7 * t256 * t201 / 0.9e1;
2036 t2818 = t7 * t233 * t1141;
2037 t2821 = t7 * t239 * t464;
2038 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2039 v4rho3sigma[0] = -t2688 / 0.2e1 - 0.3e1 / 0.2e1 * t2691 - 0.3e1 / 0.8e1 * t181 * t182 * (0.91e2 / 0.54e2 * t98 * t737 * t104 * t37 - 0.317e3 / 0.324e3 * t2697 * t484 - 0.49e2 / 0.324e3 * t281 * t2700 * t788 * t1196 - 0.49e2 / 0.324e3 * t281 * t2700 * t298 * t1202 - 0.49e2 / 0.972e3 * t281 * t2700 * t113 * t1208 + 0.2e1 / 0.81e2 * t281 * t2715 * t1761 * t2718 * t2720 + t281 * t2715 * t788 * t2726 * t2720 / 0.27e2 + 0.11e2 / 0.486e3 * t281 * t2715 * t298 * t2732 * t2720 + t281 * t2715 * t113 * t2738 * t2720 / 0.162e3 + 0.49e2 / 0.324e3 * t2743 * t1184 + t2787) - 0.4e1 * t2793 - 0.2e1 * t2796 + t2800 - t181 * t336 * t501 * t234 / 0.2e1 - 0.3e1 / 0.2e1 * t181 * t463 * t1211 * t77 - 0.3e1 / 0.2e1 * t181 * t463 * t501 * t249 - t181 * t463 * t201 * t703 / 0.2e1 - 0.2e1 / 0.3e1 * t2818 - 0.2e1 / 0.3e1 * t2821 + t181 * t813 * t201 * t644 / 0.9e1 - t1236 * t1237 * t653 / 0.2e1;
2040
2041 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2042 v4rho3sigma[1] = 0.0e0;
2043
2044 t2831 = t7 * t315 * t1221;
2045 t2834 = t7 * t321 * t508;
2046 t2837 = t7 * t123 * t1225;
2047 t2841 = t7 * t330 * t224 / 0.9e1;
2048 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2049 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2831 - 0.2e1 / 0.3e1 * t2834 - 0.2e1 * t2837 + t2841 + t181 * t828 * t224 * t664 / 0.9e1 - t1236 * t1254 * t673 / 0.2e1 - t181 * t507 * t224 * t715 / 0.2e1;
2050
2051 t2855 = t642 * t201;
2052 t2859 = t1245 * t1235;
2053 t2866 = 0.4e1 / 0.3e1 * t1232 * t1242;
2054 t2877 = t83 * t1 * t3;
2055 t2878 = t2877 * t513;
2056 t2880 = t1231 * t1235;
2057 t2882 = 0.4e1 / 0.9e1 * t2880 * t1238;
2058 t2884 = t351 * t1 * t3;
2059 t2885 = t2884 * t1248;
2060 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2061 v4rho3sigma[3] = -t2688 / 0.3e1 - t2691 / 0.2e1 + t1236 * t2855 * t815 / 0.9e1 - t2859 * t1237 * t824 / 0.3e1 - 0.4e1 / 0.3e1 * t2793 - 0.2e1 / 0.3e1 * t2796 + t2800 - t2866 - t181 * t463 * t1211 * t136 / 0.2e1 - t1246 * t463 * t201 / 0.2e1 - 0.2e1 / 0.9e1 * t2818 - 0.4e1 / 0.9e1 * t2821 - 0.2e1 / 0.9e1 * t2878 - t2882 + t2885 / 0.6e1 - t1236 * t232 * t501 * t339 / 0.3e1 - t1236 * t1237 * t851 / 0.6e1 - t1246 * t463 * t501 * t8;
2062
2063 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2064 v4rho3sigma[4] = 0.0e0;
2065
2066 t2900 = t2877 * t518;
2067 t2903 = 0.4e1 / 0.9e1 * t2880 * t1255;
2068 t2904 = t2884 * t1259;
2069 t2906 = t662 * t224;
2070 t2919 = t2877 * t555;
2071 t2922 = 0.4e1 / 0.3e1 * t1232 * t1265;
2072 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2073 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2831 - 0.4e1 / 0.9e1 * t2834 - 0.2e1 / 0.3e1 * t2837 + t2841 - 0.2e1 / 0.9e1 * t2900 - t2903 - t2904 / 0.6e1 + t1236 * t2906 * t830 / 0.9e1 + t2859 * t1254 * t843 / 0.3e1 - t1236 * t1254 * t838 / 0.6e1 + t1246 * t507 * t224 / 0.2e1 - t2919 / 0.6e1 - t2922 - t181 * t359 * t554 * t316 / 0.6e1 - t181 * t507 * t554 * t325 / 0.2e1;
2074
2075 t2935 = t1232 * t1270;
2076 t2947 = t1232 * t1274;
2077 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2078 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2878 - t2882 - t2866 - 0.4e1 / 0.3e1 * t2885 - 0.2e1 / 0.9e1 * t2821 + t2800 - t2688 / 0.6e1 - 0.2e1 / 0.9e1 * t2935 + t1236 * t2855 * t939 / 0.9e1 - t181 * t336 * t501 * t385 / 0.6e1 - t2859 * t1237 * t919 / 0.3e1 - 0.2e1 / 0.3e1 * t2947 - t1236 * t1237 * t943 / 0.6e1 - t181 * t463 * t501 * t391 / 0.2e1 - t181 * t463 * t201 * t974 / 0.2e1;
2079
2080 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2081 v4rho3sigma[7] = 0.0e0;
2082
2083 t2964 = t1232 * t1279;
2084 t2972 = t1232 * t1283;
2085 t2981 = t1232 * t1286;
2086 t2990 = t1232 * t1349;
2087 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2088 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2900 - t2903 + 0.4e1 / 0.3e1 * t2904 - 0.2e1 / 0.9e1 * t2834 + t2841 - t2919 / 0.3e1 - t2922 - 0.2e1 / 0.9e1 * t2964 + t1236 * t2906 * t935 / 0.9e1 + t2859 * t1254 * t923 / 0.3e1 - 0.4e1 / 0.3e1 * t2972 - t1236 * t314 * t554 * t362 / 0.3e1 + t1246 * t507 * t554 * t8 - 0.2e1 / 0.3e1 * t2981 - t1236 * t1254 * t931 / 0.6e1 - t181 * t507 * t224 * t982 / 0.2e1 - t2990 / 0.2e1 - t181 * t507 * t1348 * t124 / 0.2e1;
2089
2090 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2091 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2878 - 0.2e1 / 0.3e1 * t2935 - 0.2e1 * t2947 + t2800 + t181 * t813 * t201 * t1112 / 0.9e1 - t1236 * t1237 * t1117 / 0.2e1 - t181 * t463 * t201 * t1018 / 0.2e1;
2092
2093 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2094 v4rho3sigma[10] = 0.0e0;
2095
2096 t3011 = 0.1e1 / t154 / t2599;
2097 t3012 = t3011 * t58;
2098 t3015 = t300 * t2621;
2099 t3017 = t428 * t18 * t97;
2100 t3023 = params->c * t2621;
2101 t3028 = t281 * t1075 * t1079;
2102 t3034 = t281 * t3011 * t2621;
2103 t3050 = t281 * t1058;
2104 t3053 = t1075 * t58;
2105 t3069 = t1762 * t2621;
2106 t3075 = t789 * t2621;
2107 t3090 = t281 * t1058 * t58;
2108 t3101 = -0.49e2 / 0.324e3 * t281 * t3053 * t444 * t1339 - 0.49e2 / 0.972e3 * t281 * t3053 * t167 * t1345 + 0.2e1 / 0.81e2 * t281 * t3012 * t2618 * t3069 * t3017 + t281 * t3012 * t1095 * t3075 * t3017 / 0.27e2 - 0.317e3 / 0.162e3 * t3050 * t542 + 0.35e2 / 0.81e2 * t24 * t206 * t1042 * t62 - 0.91e2 / 0.54e2 * t95 * t97 * t1048 * t221 + 0.317e3 / 0.162e3 * t3090 * t549 + 0.317e3 / 0.324e3 * t3090 * t542 + 0.49e2 / 0.972e3 * t3028 * t1317 - t3034 * t61 * t428 * t1819 / 0.162e3;
2109 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2110 v4rho3sigma[11] = -0.3e1 / 0.8e1 * t181 * t205 * (0.11e2 / 0.486e3 * t281 * t3012 * t444 * t3015 * t3017 + t281 * t3012 * t167 * t3023 * t3017 / 0.162e3 + 0.49e2 / 0.324e3 * t3028 * t1321 + 0.49e2 / 0.324e3 * t3028 * t1326 - 0.11e2 / 0.486e3 * t3034 * t167 * t428 * t2752 - t3034 * t444 * t428 * t2758 / 0.27e2 - 0.2e1 / 0.81e2 * t3034 * t1096 * t3017 + 0.91e2 / 0.54e2 * t98 * t1048 * t158 * t61 - 0.317e3 / 0.324e3 * t3050 * t537 - 0.49e2 / 0.324e3 * t281 * t3053 * t1095 * t1333 + t3101) + t2841 - t2919 / 0.2e1 - 0.3e1 / 0.2e1 * t2990 - t181 * t359 * t554 * t396 / 0.2e1 - 0.3e1 / 0.2e1 * t181 * t507 * t1348 * t141 - 0.3e1 / 0.2e1 * t181 * t507 * t554 * t405 - t181 * t507 * t224 * t1036 / 0.2e1 - 0.2e1 / 0.3e1 * t2900 - 0.2e1 / 0.3e1 * t2964 - 0.4e1 * t2972 - 0.2e1 * t2981 + t181 * t828 * t224 * t1121 / 0.9e1 - t1236 * t1254 * t1126 / 0.2e1;
2111
2112 t3136 = t2877 * t592 / 0.6e1;
2113 t3137 = t1232 * t1355;
2114 t3139 = t1232 * t1404;
2115 t3152 = t281 * t183 * t285;
2116 t3162 = 0.1e1 / t100 / t762;
2117 t3164 = t281 * t3162 * t1765;
2118 t3172 = t281 * t1192;
2119 t3180 = t97 * params->c;
2120 t3186 = t1819 * t25;
2121 t3195 = t3162 * t34;
2122 t3229 = t3164 * t494 * t3186 / 0.72e2 + t3164 * t790 * t3186 / 0.108e3 + 0.4e1 / 0.27e2 * t1368 * t289 - t281 * t3195 * t1761 * t2718 * t3186 / 0.108e3 - t281 * t3195 * t788 * t2726 * t3186 / 0.72e2 - 0.11e2 / 0.1296e4 * t281 * t3195 * t298 * t2732 * t3186 - t281 * t3195 * t113 * t2738 * t3186 / 0.432e3 - 0.7e1 / 0.216e3 * t98 * t563 * t271 * t105 + 0.8e1 / 0.27e2 * t1368 * t292 - 0.35e2 / 0.2592e4 * t281 * t1176 * t770 + 0.7e1 / 0.108e3 * t24 * t558 * t264 * t38;
2123 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2124 v4rho2sigma2[0] = -t3136 - 0.4e1 / 0.3e1 * t3137 - t3139 - t181 * t336 * t591 * t234 / 0.6e1 - t181 * t463 * t1403 * t77 - t181 * t463 * t591 * t249 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t182 * (-0.8e1 / 0.27e2 * t3152 * t302 - 0.4e1 / 0.27e2 * t3152 * t306 - 0.35e2 / 0.864e3 * t1178 * t1389 - 0.35e2 / 0.864e3 * t1178 * t1392 + t3164 * t37 * t18 * t97 * t25 / 0.432e3 + 0.7e1 / 0.216e3 * t569 * t277 + 0.35e2 / 0.864e3 * t3172 * t793 + 0.35e2 / 0.864e3 * t3172 * t782 + 0.35e2 / 0.2592e4 * t3172 * t776 + 0.11e2 / 0.1296e4 * t3164 * t113 * t18 * t3180 * t25 + t3229);
2125
2126 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2127 v4rho2sigma2[1] = 0.0e0;
2128
2129 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2130 v4rho2sigma2[2] = 0.0e0;
2131
2132 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2133 v4rho2sigma2[3] = 0.0e0;
2134
2135 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2136 v4rho2sigma2[4] = 0.0e0;
2137
2138 t3235 = t2877 * t629 / 0.6e1;
2139 t3236 = t1232 * t1409;
2140 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2141 v4rho2sigma2[5] = -t3235 - 0.4e1 / 0.3e1 * t3236 - t181 * t359 * t628 * t316 / 0.6e1 - t181 * t507 * t628 * t325 / 0.2e1;
2142
2143 t3248 = t1232 * t1413;
2144 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2145 v4rho2sigma2[6] = -t3136 - 0.2e1 / 0.3e1 * t3137 - t3139 / 0.2e1 - 0.2e1 / 0.3e1 * t3248 - t1236 * t232 * t591 * t339 / 0.6e1 - t181 * t463 * t1403 * t136 / 0.2e1 - t1246 * t463 * t591 * t8 / 0.2e1;
2146
2147 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2148 v4rho2sigma2[7] = 0.0e0;
2149
2150 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2151 v4rho2sigma2[8] = 0.0e0;
2152
2153 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2154 v4rho2sigma2[9] = 0.0e0;
2155
2156 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2157 v4rho2sigma2[10] = 0.0e0;
2158
2159 t3263 = t1232 * t1418;
2160 t3273 = t1232 * t1467;
2161 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2162 v4rho2sigma2[11] = -t3235 - 0.2e1 / 0.3e1 * t3236 - 0.2e1 / 0.3e1 * t3263 - t1236 * t314 * t628 * t362 / 0.6e1 + t1246 * t507 * t628 * t8 / 0.2e1 - t3273 / 0.2e1 - t181 * t507 * t1466 * t124 / 0.2e1;
2163
2164 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2165 v4rho2sigma2[12] = -t3136 - 0.4e1 / 0.3e1 * t3248 - t181 * t336 * t591 * t385 / 0.6e1 - t181 * t463 * t591 * t391 / 0.2e1;
2166
2167 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2168 v4rho2sigma2[13] = 0.0e0;
2169
2170 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2171 v4rho2sigma2[14] = 0.0e0;
2172
2173 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2174 v4rho2sigma2[15] = 0.0e0;
2175
2176 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2177 v4rho2sigma2[16] = 0.0e0;
2178
2179 t3305 = t281 * t206 * t431;
2180 t3315 = 0.1e1 / t154 / t1073;
2181 t3317 = t281 * t3315 * t2621;
2182 t3325 = t281 * t1329;
2183 t3338 = t1819 * t50;
2184 t3347 = t3315 * t58;
2185 t3377 = 0.11e2 / 0.1296e4 * t3317 * t167 * t18 * t3180 * t50 + t3317 * t547 * t3338 / 0.72e2 + t3317 * t1096 * t3338 / 0.108e3 + 0.4e1 / 0.27e2 * t1431 * t435 - t281 * t3347 * t2618 * t3069 * t3338 / 0.108e3 - t281 * t3347 * t1095 * t3075 * t3338 / 0.72e2 - 0.11e2 / 0.1296e4 * t281 * t3347 * t444 * t3015 * t3338 - t281 * t3347 * t167 * t3023 * t3338 / 0.432e3 - 0.7e1 / 0.216e3 * t98 * t600 * t419 * t159 + 0.8e1 / 0.27e2 * t1431 * t438 - 0.35e2 / 0.2592e4 * t281 * t1313 * t1081;
2186 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2187 v4rho2sigma2[17] = -t3235 - 0.4e1 / 0.3e1 * t3263 - t3273 - t181 * t359 * t628 * t396 / 0.6e1 - t181 * t507 * t1466 * t141 - t181 * t507 * t628 * t405 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t205 * (0.7e1 / 0.108e3 * t24 * t595 * t412 * t62 - 0.8e1 / 0.27e2 * t3305 * t447 - 0.4e1 / 0.27e2 * t3305 * t451 - 0.35e2 / 0.864e3 * t1315 * t1452 - 0.35e2 / 0.864e3 * t1315 * t1455 + t3317 * t61 * t18 * t97 * t50 / 0.432e3 + 0.7e1 / 0.216e3 * t606 * t425 + 0.35e2 / 0.864e3 * t3325 * t1099 + 0.35e2 / 0.864e3 * t3325 * t1089 + 0.35e2 / 0.2592e4 * t3325 * t1085 + t3377);
2188
2189 t3382 = t1232 * t1519;
2190 t3393 = t281 * t558 * t481;
2191 t3399 = 0.1e1 / t100 / t746;
2192 t3401 = t281 * t3399 * t1765;
2193 t3419 = t183 * t18 * t97;
2194 t3423 = t281 * t558;
2195 t3431 = t3399 * t34;
2196 t3470 = t1388 * t1494 / 0.144e3 + t1388 * t1499 / 0.144e3 + t281 * t3431 * t1761 * t2718 * t3419 / 0.288e3 + t281 * t3431 * t788 * t2726 * t3419 / 0.192e3 + 0.11e2 / 0.3456e4 * t281 * t3431 * t298 * t2732 * t3419 + t281 * t3431 * t113 * t2738 * t3419 / 0.1152e4 - t281 * t1395 * t788 * t1505 / 0.144e3 - t281 * t1395 * t298 * t1510 / 0.144e3 - t281 * t1395 * t113 * t1515 / 0.432e3 - t98 * t1476 * t102 * t105 / 0.48e2 + t3423 * t1373 / 0.48e2;
2197 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2198 v4rhosigma3[0] = -t3382 / 0.2e1 - t181 * t463 * t1518 * t77 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t182 * (t24 * t1471 * t90 * t38 / 0.24e2 - t3393 * t302 / 0.48e2 - t3393 * t306 / 0.96e2 - t3401 * t37 * t183 * t1819 / 0.1152e4 + t1482 * t116 / 0.48e2 + t1388 * t1490 / 0.432e3 - 0.11e2 / 0.3456e4 * t3401 * t113 * t183 * t2752 - t3401 * t298 * t183 * t2758 / 0.192e3 - t3401 * t790 * t3419 / 0.288e3 + t3423 * t1370 / 0.96e2 + t3470);
2199
2200 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2201 v4rhosigma3[1] = 0.0e0;
2202
2203 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2204 v4rhosigma3[2] = 0.0e0;
2205
2206 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2207 v4rhosigma3[3] = 0.0e0;
2208
2209 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2210 v4rhosigma3[4] = 0.0e0;
2211
2212 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2213 v4rhosigma3[5] = 0.0e0;
2214
2215 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2216 v4rhosigma3[6] = 0.0e0;
2217
2218 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2219 v4rhosigma3[7] = 0.0e0;
2220
2221 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2222 v4rhosigma3[8] = 0.0e0;
2223
2224 t3475 = t1232 * t1571;
2225 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2226 v4rhosigma3[9] = -t181 * t507 * t1570 * t124 / 0.2e1 - t3475 / 0.2e1;
2227
2228 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2229 v4rhosigma3[10] = -t181 * t463 * t1518 * t136 / 0.2e1 - t3382 / 0.2e1;
2230
2231 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2232 v4rhosigma3[11] = 0.0e0;
2233
2234 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2235 v4rhosigma3[12] = 0.0e0;
2236
2237 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2238 v4rhosigma3[13] = 0.0e0;
2239
2240 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2241 v4rhosigma3[14] = 0.0e0;
2242
2243 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2244 v4rhosigma3[15] = 0.0e0;
2245
2246 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2247 v4rhosigma3[16] = 0.0e0;
2248
2249 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2250 v4rhosigma3[17] = 0.0e0;
2251
2252 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2253 v4rhosigma3[18] = 0.0e0;
2254
2255 t3494 = t281 * t595 * t534;
2256 t3500 = 0.1e1 / t154 / t1057;
2257 t3502 = t281 * t3500 * t2621;
2258 t3511 = t281 * t595;
2259 t3528 = t206 * t18 * t97;
2260 t3538 = t3500 * t58;
2261 t3571 = -t3502 * t1096 * t3528 / 0.288e3 + t3511 * t1433 / 0.96e2 + t1451 * t1546 / 0.144e3 + t1451 * t1551 / 0.144e3 + t281 * t3538 * t2618 * t3069 * t3528 / 0.288e3 + t281 * t3538 * t1095 * t3075 * t3528 / 0.192e3 + 0.11e2 / 0.3456e4 * t281 * t3538 * t444 * t3015 * t3528 + t281 * t3538 * t167 * t3023 * t3528 / 0.1152e4 - t281 * t1458 * t1095 * t1557 / 0.144e3 - t281 * t1458 * t444 * t1562 / 0.144e3 - t281 * t1458 * t167 * t1567 / 0.432e3;
2262 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2263 v4rhosigma3[19] = -t3475 / 0.2e1 - t181 * t507 * t1570 * t141 / 0.2e1 - 0.3e1 / 0.8e1 * t181 * t205 * (t24 * t1523 * t148 * t62 / 0.24e2 - t3494 * t447 / 0.48e2 - t3494 * t451 / 0.96e2 - t3502 * t61 * t206 * t1819 / 0.1152e4 - t98 * t1528 * t156 * t159 / 0.48e2 + t3511 * t1436 / 0.48e2 + t1534 * t170 / 0.48e2 + t1451 * t1542 / 0.432e3 - 0.11e2 / 0.3456e4 * t3502 * t167 * t206 * t2752 - t3502 * t444 * t206 * t2758 / 0.192e3 + t3571);
2264
2265 t3576 = t761 * sigma[0];
2266 t3584 = t281 * t1471 * t574;
2267 t3590 = t281 * t1781 * t1765;
2268 t3595 = 0.1e1 / t3576;
2269 t3600 = t281 * t1471;
2270 t3601 = t574 * t287;
2271 t3623 = t558 * t18 * t97;
2272 t3635 = t1476 * t19 * t23;
2273 t3668 = t3590 * t790 * t3623 / 0.768e3 - t3600 * t3601 * t37 / 0.256e3 + t1488 * t113 * t1476 * t775 / 0.768e3 + t1488 * t494 * t3635 / 0.768e3 - t281 * t1790 * t1761 * t2718 * t3623 / 0.768e3 - t281 * t1790 * t788 * t2726 * t3623 / 0.512e3 - 0.11e2 / 0.9216e4 * t281 * t1790 * t298 * t2732 * t3623 - t281 * t1790 * t113 * t2738 * t3623 / 0.3072e4 - t1504 * t1195 * t3635 / 0.768e3 - t1509 * t1201 * t3635 / 0.768e3 - t1514 * t1207 * t3635 / 0.2304e4;
2274 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2275 v4sigma4[0] = -0.3e1 / 0.8e1 * t181 * t182 * (0.5e1 / 0.64e2 * t24 / t25 / t3576 * t28 * t38 + t3584 * t302 / 0.128e3 + t3584 * t306 / 0.256e3 + t3590 * t37 * t558 * t1819 / 0.3072e4 - 0.5e1 / 0.128e3 * t98 * t3595 * t189 * t105 - t3600 * t3601 * t114 / 0.128e3 + 0.5e1 / 0.128e3 * t95 * t97 * t3595 * t571 + t1488 * t37 * t1476 * t30 / 0.2304e4 + 0.11e2 / 0.9216e4 * t3590 * t113 * t558 * t2752 + t3590 * t298 * t558 * t2758 / 0.512e3 + t3668);
2276
2277 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2278 v4sigma4[1] = 0.0e0;
2279
2280 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2281 v4sigma4[2] = 0.0e0;
2282
2283 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2284 v4sigma4[3] = 0.0e0;
2285
2286 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2287 v4sigma4[4] = 0.0e0;
2288
2289 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2290 v4sigma4[5] = 0.0e0;
2291
2292 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2293 v4sigma4[6] = 0.0e0;
2294
2295 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2296 v4sigma4[7] = 0.0e0;
2297
2298 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2299 v4sigma4[8] = 0.0e0;
2300
2301 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2302 v4sigma4[9] = 0.0e0;
2303
2304 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2305 v4sigma4[10] = 0.0e0;
2306
2307 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2308 v4sigma4[11] = 0.0e0;
2309
2310 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2311 v4sigma4[12] = 0.0e0;
2312
2313 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2314 v4sigma4[13] = 0.0e0;
2315
2316 t3676 = t595 * t18 * t97;
2317 t3685 = t1528 * t19 * t23;
2318 t3704 = t1072 * sigma[2];
2319 t3712 = t281 * t1523 * t611;
2320 t3718 = t281 * t2637 * t2621;
2321 t3724 = 0.1e1 / t3704;
2322 t3729 = t281 * t1523;
2323 t3730 = t611 * t433;
2324 t3765 = -0.5e1 / 0.128e3 * t98 * t3724 * t212 * t159 - t3729 * t3730 * t168 / 0.128e3 - t3729 * t3730 * t61 / 0.256e3 - t1556 * t1332 * t3685 / 0.768e3 - t1561 * t1338 * t3685 / 0.768e3 - t1566 * t1344 * t3685 / 0.2304e4 + 0.5e1 / 0.128e3 * t95 * t97 * t3724 * t608 + t1540 * t61 * t1528 * t30 / 0.2304e4 + 0.11e2 / 0.9216e4 * t3718 * t167 * t595 * t2752 + t3718 * t444 * t595 * t2758 / 0.512e3 + t3718 * t1096 * t3676 / 0.768e3;
2325 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2326 v4sigma4[14] = -0.3e1 / 0.8e1 * t181 * t205 * (-t281 * t2638 * t167 * t3023 * t3676 / 0.3072e4 + t1540 * t167 * t1528 * t775 / 0.768e3 + t1540 * t547 * t3685 / 0.768e3 - t281 * t2638 * t2618 * t3069 * t3676 / 0.768e3 - t281 * t2638 * t1095 * t3075 * t3676 / 0.512e3 - 0.11e2 / 0.9216e4 * t281 * t2638 * t444 * t3015 * t3676 + 0.5e1 / 0.64e2 * t24 / t50 / t3704 * t53 * t62 + t3712 * t447 / 0.128e3 + t3712 * t451 / 0.256e3 + t3718 * t61 * t595 * t1819 / 0.3072e4 + t3765);
2327
2328 #ifndef XC_DONT_COMPILE_MXC
2329
2330 if(order < 5) return;
2331
2332
2333 #endif
2334
2335 #endif
2336
2337 #endif
2338
2339 #endif
2340
2341 #endif
2342
2343
2344 }
2345
2346