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_dk87.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t2, t3, t5, t6, t7, t8, t9;
24 double t10, t11, t12, t13, t14, t15, t16, t17;
25 double t18, t20, t22, t23, t24, t25, t26, t28;
26 double t29, t31, t34, t35, t36, t38, t39, t41;
27 double t42, t43, t47, t49;
28
29 #ifndef XC_DONT_COMPILE_VXC
30 double t53, t54, t55, t57, t59, t63, t64, t66;
31 double t67, t70, t74, t75, t76, t77, t79, t81;
32 double t82, t84, t88, t92, t96, t97, t98, t103;
33 double t105, t111;
34
35 #ifndef XC_DONT_COMPILE_FXC
36 double t114, t115, t123, t125, t129, t133, t135, t140;
37 double t142, t146, t148, t149, t150, t151, t152, t153;
38 double t156, t157, t158, t159, t160, t163, t165, t166;
39 double t170, t175, t182, t186, t188, t192, t193, t197;
40 double t199, t201, t203, t207, t208, t214, t217, t218;
41 double t219, t227, t236, t242;
42
43 #ifndef XC_DONT_COMPILE_KXC
44 double t246, t247, t258, t260, t264, t269, t277, t278;
45 double t282, t283, t288, t290, t294, t300, t301, t305;
46 double t306, t308, t311, t312, t313, t314, t315, t319;
47 double t324, t333, t345, t350, t351, t355, t356, t361;
48 double t366, t367, t368, t369, t373, t375, t382, t393;
49 double t401, t403, t407, t415, t416, t417, t422, t429;
50 double t432, t433, t434, t439, t440, t441, t449, t453;
51 double t454, t459, t460, t465, t472;
52
53 #ifndef XC_DONT_COMPILE_LXC
54 double t487, t494, t499, t501, t508, t516, t518, t522;
55 double t527, t538, t542, t553, t560, t571, t575, t588;
56 double t609, t630, t634, t640, t655, t661, t727, t789;
57 double t794, t846;
58 #endif
59
60 #endif
61
62 #endif
63
64 #endif
65
66 #endif
67
68
69 gga_x_dk87_params *params;
70
71 assert(p->params != NULL);
72 params = (gga_x_dk87_params * )(p->params);
73
74 t1 = M_CBRT3;
75 t2 = 0.1e1 / M_PI;
76 t3 = POW_1_3(t2);
77 t5 = M_CBRT4;
78 t6 = t5 * t5;
79 t7 = t1 * t3 * t6;
80 t8 = M_CBRT2;
81 t9 = t8 * t8;
82 t10 = POW_1_3(rho[0]);
83 t11 = t9 * t10;
84 t12 = M_CBRT6;
85 t13 = t12 * t12;
86 t14 = t2 * t13;
87 t15 = M_PI * M_PI;
88 t16 = POW_1_3(t15);
89 t17 = 0.1e1 / t16;
90 t18 = t1 * t1;
91 t20 = 0.1e1 / t3;
92 t22 = t14 * t17 * t18 * t20;
93 t23 = t5 * sigma[0];
94 t24 = t23 * t9;
95 t25 = rho[0] * rho[0];
96 t26 = t10 * t10;
97 t28 = 0.1e1 / t26 / t25;
98 t29 = sqrt(sigma[0]);
99 t31 = t10 * rho[0];
100 t34 = pow(t29 * t8 / t31, params->alpha);
101 t35 = params->a1 * t34;
102 t36 = 0.1e1 + t35;
103 t38 = params->b1 * sigma[0];
104 t39 = t9 * t28;
105 t41 = t38 * t39 + 0.1e1;
106 t42 = 0.1e1 / t41;
107 t43 = t28 * t36 * t42;
108 t47 = 0.1e1 + 0.7e1 / 0.11664e5 * t22 * t24 * t43;
109 t49 = t7 * t11 * t47;
110 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
111 zk[0] = -0.3e1 / 0.16e2 * t49;
112
113 #ifndef XC_DONT_COMPILE_VXC
114
115 if(order < 1) return;
116
117
118 t53 = t31 * t1 * t3;
119 t54 = t6 * t9;
120 t55 = t25 * rho[0];
121 t57 = 0.1e1 / t26 / t55;
122 t59 = t57 * t36 * t42;
123 t63 = t14 * t17;
124 t64 = t18 * t20;
125 t66 = t63 * t64 * t5;
126 t67 = sigma[0] * t9;
127 t70 = t35 * params->alpha * t42;
128 t74 = sigma[0] * sigma[0];
129 t75 = t74 * t8;
130 t76 = t25 * t25;
131 t77 = t76 * t25;
132 t79 = 0.1e1 / t10 / t77;
133 t81 = t41 * t41;
134 t82 = 0.1e1 / t81;
135 t84 = t36 * t82 * params->b1;
136 t88 = -0.7e1 / 0.4374e4 * t22 * t24 * t59 - 0.7e1 / 0.8748e4 * t66 * t67 * t57 * t70 + 0.7e1 / 0.2187e4 * t66 * t75 * t79 * t84;
137 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
138 vrho[0] = -t49 / 0.4e1 - 0.3e1 / 0.16e2 * t53 * t54 * t88;
139
140 t92 = t5 * t9;
141 t96 = t39 * params->a1;
142 t97 = t34 * params->alpha;
143 t98 = t97 * t42;
144 t103 = t76 * rho[0];
145 t105 = 0.1e1 / t10 / t103;
146 t111 = t54 * (0.7e1 / 0.11664e5 * t22 * t92 * t43 + 0.7e1 / 0.23328e5 * t66 * t96 * t98 - 0.7e1 / 0.5832e4 * t66 * sigma[0] * t8 * t105 * t84);
147 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
148 vsigma[0] = -0.3e1 / 0.16e2 * t53 * t111;
149
150 #ifndef XC_DONT_COMPILE_FXC
151
152 if(order < 2) return;
153
154
155 t114 = 0.1e1 / t26;
156 t115 = t9 * t114;
157 t123 = 0.1e1 / t26 / t76;
158 t125 = t123 * t36 * t42;
159 t129 = t67 * t123;
160 t133 = t76 * t55;
161 t135 = 0.1e1 / t10 / t133;
162 t140 = params->alpha * params->alpha;
163 t142 = t35 * t140 * t42;
164 t146 = t5 * t74;
165 t148 = t63 * t64 * t146;
166 t149 = t8 * t135;
167 t150 = t149 * params->a1;
168 t151 = t82 * params->b1;
169 t152 = t97 * t151;
170 t153 = t150 * t152;
171 t156 = t74 * sigma[0];
172 t157 = t5 * t156;
173 t158 = t76 * t76;
174 t159 = t158 * t25;
175 t160 = 0.1e1 / t159;
176 t163 = 0.1e1 / t81 / t41;
177 t165 = params->b1 * params->b1;
178 t166 = t36 * t163 * t165;
179 t170 = 0.77e2 / 0.13122e5 * t22 * t24 * t125 + 0.133e3 / 0.26244e5 * t66 * t129 * t70 - 0.7e1 / 0.243e3 * t66 * t75 * t135 * t84 + 0.7e1 / 0.6561e4 * t66 * t129 * t142 - 0.56e2 / 0.6561e4 * t148 * t153 + 0.224e3 / 0.6561e4 * t22 * t157 * t160 * t166;
180 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
181 v2rho2[0] = -t7 * t115 * t47 / 0.12e2 - t7 * t11 * t88 / 0.2e1 - 0.3e1 / 0.16e2 * t53 * t54 * t170;
182
183 t175 = t10 * t1 * t3;
184 t182 = t9 * t57 * params->a1;
185 t186 = t8 * t79;
186 t188 = t151 * sigma[0];
187 t192 = t34 * t140;
188 t193 = t192 * t42;
189 t197 = t5 * t8;
190 t199 = t63 * t64 * t197;
191 t201 = t79 * params->a1 * t34;
192 t203 = params->alpha * t82 * t38;
193 t207 = t158 * rho[0];
194 t208 = 0.1e1 / t207;
195 t214 = t54 * (-0.7e1 / 0.4374e4 * t22 * t92 * t59 - 0.7e1 / 0.4374e4 * t66 * t182 * t98 + 0.7e1 / 0.729e3 * t66 * t186 * t36 * t188 - 0.7e1 / 0.17496e5 * t66 * t182 * t193 + 0.7e1 / 0.2187e4 * t199 * t201 * t203 - 0.28e2 / 0.2187e4 * t22 * t146 * t208 * t166);
196 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
197 v2rhosigma[0] = -t175 * t111 / 0.4e1 - 0.3e1 / 0.16e2 * t53 * t214;
198
199 t217 = 0.1e1 / sigma[0];
200 t218 = t217 * t42;
201 t219 = t97 * t218;
202 t227 = t192 * t218;
203 t236 = 0.1e1 / t158;
204 t242 = t54 * (0.7e1 / 0.23328e5 * t66 * t96 * t219 - 0.7e1 / 0.2916e4 * t22 * t197 * t105 * t84 + 0.7e1 / 0.46656e5 * t66 * t96 * t227 - 0.7e1 / 0.5832e4 * t66 * t8 * t105 * params->a1 * t152 + 0.7e1 / 0.1458e4 * t22 * t23 * t236 * t166);
205 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
206 v2sigma2[0] = -0.3e1 / 0.16e2 * t53 * t242;
207
208 #ifndef XC_DONT_COMPILE_KXC
209
210 if(order < 3) return;
211
212
213 t246 = 0.1e1 / t26 / rho[0];
214 t247 = t9 * t246;
215 t258 = 0.1e1 / t26 / t103;
216 t260 = t258 * t36 * t42;
217 t264 = t67 * t258;
218 t269 = 0.1e1 / t10 / t158;
219 t277 = t8 * t269;
220 t278 = t277 * params->a1;
221 t282 = t158 * t55;
222 t283 = 0.1e1 / t282;
223 t288 = t140 * params->alpha;
224 t290 = t35 * t288 * t42;
225 t294 = t192 * t151;
226 t300 = t163 * t165;
227 t301 = t97 * t300;
228 t305 = t74 * t74;
229 t306 = t158 * t103;
230 t308 = 0.1e1 / t26 / t306;
231 t311 = t81 * t81;
232 t312 = 0.1e1 / t311;
233 t313 = t165 * params->b1;
234 t314 = t312 * t313;
235 t315 = t314 * t9;
236 t319 = -0.539e3 / 0.19683e5 * t22 * t24 * t260 - 0.413e3 / 0.13122e5 * t66 * t264 * t70 + 0.4774e4 / 0.19683e5 * t66 * t75 * t269 * t84 - 0.77e2 / 0.6561e4 * t66 * t264 * t142 + 0.280e3 / 0.2187e4 * t148 * t278 * t152 - 0.4256e4 / 0.6561e4 * t22 * t157 * t283 * t166 - 0.28e2 / 0.19683e5 * t66 * t264 * t290 + 0.112e3 / 0.6561e4 * t148 * t278 * t294 - 0.896e3 / 0.6561e4 * t66 * t156 * t283 * params->a1 * t301 + 0.1792e4 / 0.6561e4 * t66 * t305 * t308 * t36 * t315;
237 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
238 v3rho3[0] = t7 * t247 * t47 / 0.18e2 - t7 * t115 * t88 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t170 - 0.3e1 / 0.16e2 * t53 * t54 * t319;
239
240 t324 = t114 * t1 * t3;
241 t333 = t9 * t123 * params->a1;
242 t345 = t135 * params->a1 * t34;
243 t350 = t5 * t160 * t36;
244 t351 = t300 * t74;
245 t355 = t34 * t288;
246 t356 = t355 * t42;
247 t361 = t140 * t82 * t38;
248 t366 = t160 * params->a1 * t34;
249 t367 = params->alpha * t163;
250 t368 = t165 * t74;
251 t369 = t367 * t368;
252 t373 = t158 * t76;
253 t375 = 0.1e1 / t26 / t373;
254 t382 = t54 * (0.77e2 / 0.13122e5 * t22 * t92 * t125 + 0.35e2 / 0.4374e4 * t66 * t333 * t98 - 0.455e3 / 0.6561e4 * t66 * t149 * t36 * t188 + 0.7e1 / 0.1944e4 * t66 * t333 * t193 - 0.91e2 / 0.2187e4 * t199 * t345 * t203 + 0.476e3 / 0.2187e4 * t22 * t350 * t351 + 0.7e1 / 0.13122e5 * t66 * t333 * t356 - 0.14e2 / 0.2187e4 * t199 * t345 * t361 + 0.112e3 / 0.2187e4 * t66 * t366 * t369 - 0.224e3 / 0.2187e4 * t66 * t156 * t375 * t36 * t315);
255 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
256 v3rho2sigma[0] = -t324 * t111 / 0.12e2 - t175 * t214 / 0.2e1 - 0.3e1 / 0.16e2 * t53 * t382;
257
258 t393 = t186 * params->a1;
259 t401 = t5 * t208;
260 t403 = t300 * sigma[0];
261 t407 = t355 * t218;
262 t415 = t208 * params->a1 * t34;
263 t416 = t165 * sigma[0];
264 t417 = t367 * t416;
265 t422 = 0.1e1 / t26 / t282;
266 t429 = t54 * (-0.7e1 / 0.8748e4 * t66 * t182 * t219 - 0.7e1 / 0.8748e4 * t66 * t182 * t227 + 0.49e2 / 0.4374e4 * t66 * t393 * t152 + 0.28e2 / 0.2187e4 * t22 * t197 * t79 * t84 - 0.140e3 / 0.2187e4 * t22 * t401 * t36 * t403 - 0.7e1 / 0.34992e5 * t66 * t182 * t407 + 0.7e1 / 0.2916e4 * t66 * t393 * t294 - 0.14e2 / 0.729e3 * t66 * t415 * t417 + 0.28e2 / 0.729e3 * t66 * t74 * t422 * t36 * t315);
267 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
268 v3rhosigma2[0] = -t175 * t242 / 0.4e1 - 0.3e1 / 0.16e2 * t53 * t429;
269
270 t432 = 0.1e1 / t74;
271 t433 = t432 * t42;
272 t434 = t97 * t433;
273 t439 = t105 * params->a1 * t34;
274 t440 = params->alpha * t217;
275 t441 = t440 * t151;
276 t449 = t355 * t433;
277 t453 = t140 * t217;
278 t454 = t453 * t151;
279 t459 = t236 * params->a1 * t34;
280 t460 = t367 * t165;
281 t465 = 0.1e1 / t26 / t159;
282 t472 = t54 * (-0.7e1 / 0.23328e5 * t66 * t96 * t434 - 0.7e1 / 0.3888e4 * t199 * t439 * t441 + 0.7e1 / 0.486e3 * t22 * t5 * t236 * t166 + 0.7e1 / 0.93312e5 * t66 * t96 * t449 - 0.7e1 / 0.7776e4 * t199 * t439 * t454 + 0.7e1 / 0.972e3 * t66 * t459 * t460 - 0.7e1 / 0.486e3 * t66 * sigma[0] * t465 * t36 * t315);
283 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
284 v3sigma3[0] = -0.3e1 / 0.16e2 * t53 * t472;
285
286 #ifndef XC_DONT_COMPILE_LXC
287
288 if(order < 4) return;
289
290
291 t487 = 0.1e1 / t26 / t77;
292 t494 = 0.1e1 / t10 / t207;
293 t499 = t158 * t77;
294 t501 = 0.1e1 / t26 / t499;
295 t508 = t158 * t158;
296 t516 = t165 * t165;
297 t518 = 0.1e1 / t311 / t41 * t516 * t8;
298 t522 = 0.1e1 / t373;
299 t527 = t67 * t487;
300 t538 = t156 * t522 * params->a1;
301 t542 = t140 * t140;
302 t553 = t8 * t494 * params->a1;
303 t560 = t355 * t151;
304 t571 = params->alpha * t312 * t313 * t9;
305 t575 = 0.9163e4 / 0.59049e5 * t22 * t24 * t487 * t36 * t42 - 0.42658e5 / 0.19683e5 * t66 * t75 * t494 * t84 - 0.175616e6 / 0.19683e5 * t66 * t305 * t501 * t36 * t315 + 0.114688e6 / 0.19683e5 * t66 * t305 * sigma[0] / t10 / t508 / rho[0] * t36 * t518 + 0.574112e6 / 0.59049e5 * t22 * t157 * t522 * t166 + 0.25375e5 / 0.118098e6 * t66 * t527 * t70 + 0.2135e4 / 0.19683e5 * t66 * t527 * t142 + 0.1400e4 / 0.59049e5 * t66 * t527 * t290 + 0.73472e5 / 0.19683e5 * t66 * t538 * t301 + 0.112e3 / 0.59049e5 * t66 * t527 * t35 * t542 * t42 + 0.7168e4 / 0.19683e5 * t66 * t538 * t192 * t300 - 0.92008e5 / 0.59049e5 * t148 * t553 * t152 - 0.2464e4 / 0.6561e4 * t148 * t553 * t294 - 0.1792e4 / 0.59049e5 * t148 * t553 * t560 - 0.28672e5 / 0.19683e5 * t63 * t64 * t5 * t305 * t501 * params->a1 * t34 * t571;
306 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
307 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t39 * t47 + 0.2e1 / 0.9e1 * t7 * t247 * t88 - t7 * t115 * t170 / 0.2e1 - t7 * t11 * t319 - 0.3e1 / 0.16e2 * t53 * t54 * t575;
308
309 t588 = t9 * t258 * params->a1;
310 t609 = t34 * t542;
311 t630 = t283 * params->a1 * t34;
312 t634 = t140 * t163;
313 t640 = t269 * params->a1 * t34;
314 t655 = t35 * params->alpha;
315 t661 = -0.889e3 / 0.19683e5 * t66 * t588 * t98 + 0.3542e4 / 0.6561e4 * t66 * t277 * t36 * t188 - 0.721e3 / 0.26244e5 * t66 * t588 * t193 - 0.287e3 / 0.39366e5 * t66 * t588 * t356 + 0.19936e5 / 0.6561e4 * t66 * t308 * t36 * t312 * t313 * t156 * t9 - 0.14e2 / 0.19683e5 * t66 * t588 * t609 * t42 - 0.14336e5 / 0.6561e4 * t66 * t305 / t10 / t508 * t36 * t518 - 0.539e3 / 0.19683e5 * t22 * t92 * t260 - 0.57400e5 / 0.19683e5 * t22 * t5 * t283 * t36 * t351 - 0.8176e4 / 0.6561e4 * t66 * t630 * t369 - 0.896e3 / 0.6561e4 * t66 * t630 * t634 * t368 + 0.8666e4 / 0.19683e5 * t199 * t640 * t203 + 0.266e3 / 0.2187e4 * t199 * t640 * t361 + 0.224e3 / 0.19683e5 * t199 * t640 * t288 * t82 * t38 + 0.3584e4 / 0.6561e4 * t63 * t64 * t5 * t308 * t655 * t314 * t156 * t9;
316 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
317 v4rho3sigma[0] = t246 * t1 * t3 * t111 / 0.18e2 - t324 * t214 / 0.4e1 - 0.3e1 / 0.4e1 * t175 * t382 - 0.3e1 / 0.16e2 * t53 * t54 * t661;
318
319 t727 = 0.4676e4 / 0.6561e4 * t22 * t350 * t403 + 0.35e2 / 0.8748e4 * t66 * t333 * t227 - 0.1211e4 / 0.13122e5 * t66 * t153 + 0.7e1 / 0.3888e4 * t66 * t333 * t407 - 0.301e3 / 0.8748e4 * t66 * t150 * t294 + 0.826e3 / 0.2187e4 * t66 * t366 * t417 + 0.7e1 / 0.26244e5 * t66 * t333 * t609 * t218 - 0.28e2 / 0.6561e4 * t66 * t150 * t560 + 0.112e3 / 0.2187e4 * t66 * t366 * t634 * t416 + 0.77e2 / 0.26244e5 * t66 * t333 * t219 - 0.448e3 / 0.2187e4 * t63 * t64 * t5 * t375 * t655 * t314 * t74 * t9 - 0.532e3 / 0.6561e4 * t22 * t197 * t135 * t84 - 0.700e3 / 0.729e3 * t66 * t375 * t36 * t312 * t313 * t74 * t9 + 0.1792e4 / 0.2187e4 * t66 * t156 / t10 / t158 / t133 * t36 * t518;
320 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
321 v4rho2sigma2[0] = -t324 * t242 / 0.12e2 - t175 * t429 / 0.2e1 - 0.3e1 / 0.16e2 * t53 * t54 * t727;
322
323 t789 = 0.7e1 / 0.8748e4 * t66 * t182 * t434 + 0.7e1 / 0.17496e5 * t66 * t182 * t192 * t433 + 0.35e2 / 0.4374e4 * t199 * t201 * t441 + 0.7e1 / 0.972e3 * t199 * t201 * t454 - 0.70e2 / 0.729e3 * t66 * t415 * t460 - 0.28e2 / 0.243e3 * t22 * t401 * t166 + 0.196e3 / 0.729e3 * t66 * t422 * t36 * t312 * t313 * sigma[0] * t9 - 0.7e1 / 0.34992e5 * t66 * t182 * t449 - 0.7e1 / 0.69984e5 * t66 * t182 * t609 * t433 + 0.7e1 / 0.4374e4 * t199 * t201 * t288 * t217 * t151 - 0.14e2 / 0.729e3 * t66 * t415 * t634 * t165 + 0.56e2 / 0.729e3 * t63 * t64 * t5 * t422 * t655 * t314 * t67 - 0.224e3 / 0.729e3 * t66 * t74 / t10 / t499 * t36 * t518;
324 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
325 v4rhosigma3[0] = -t175 * t472 / 0.4e1 - 0.3e1 / 0.16e2 * t53 * t54 * t789;
326
327 t794 = 0.1e1 / t156 * t42;
328 t846 = -0.7e1 / 0.46656e5 * t66 * t96 * t192 * t794 + 0.7e1 / 0.11664e5 * t66 * t96 * t97 * t794 + 0.7e1 / 0.2916e4 * t199 * t439 * params->alpha * t432 * t151 + 0.7e1 / 0.486e3 * t66 * t459 * t440 * t300 - 0.14e2 / 0.243e3 * t22 * t5 * t465 * t36 * t315 + 0.7e1 / 0.186624e6 * t66 * t96 * t609 * t794 - 0.7e1 / 0.46656e5 * t66 * t96 * t355 * t794 - 0.7e1 / 0.11664e5 * t199 * t439 * t288 * t432 * t151 + 0.7e1 / 0.972e3 * t66 * t459 * t453 * t300 - 0.7e1 / 0.243e3 * t66 * t465 * params->a1 * t34 * t571 + 0.28e2 / 0.243e3 * t66 * sigma[0] / t10 / t306 * t36 * t518;
329 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
330 v4sigma4[0] = -0.3e1 / 0.16e2 * t53 * t54 * t846;
331
332 #ifndef XC_DONT_COMPILE_MXC
333
334 if(order < 5) return;
335
336
337 #endif
338
339 #endif
340
341 #endif
342
343 #endif
344
345 #endif
346
347
348 }
349
350
351 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 *))352 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
353 {
354
355 #ifndef XC_DONT_COMPILE_EXC
356 double t1, t2, t3, t4, t5, t6, t7, t8;
357 double t9, t10, t11, t12, t13, t14, t15, t16;
358 double t17, t19, t20, t21, t22, t24, t25, t26;
359 double t29, t31, t33, t35, t36, t37, t41, t43;
360
361 #ifndef XC_DONT_COMPILE_VXC
362 double t46, t47, t48, t50, t52, t59, t62, t63;
363 double t67, t68, t69, t70, t72, t74, t75, t77;
364 double t81, t85, t86, t95, t97, t103;
365
366 #ifndef XC_DONT_COMPILE_FXC
367 double t106, t107, t114, t116, t121, t125, t127, t132;
368 double t133, t134, t140, t141, t145, t146, t147, t148;
369 double t149, t152, t154, t155, t159, t163, t170, t174;
370 double t176, t184, t185, t186, t190, t191, t197, t201;
371 double t202, t203, t204, t212, t213, t218, t219, t223;
372 double t229;
373
374 #ifndef XC_DONT_COMPILE_KXC
375 double t233, t234, t245, t247, t252, t257, t266, t270;
376 double t271, t276, t277, t278, t282, t288, t289, t293;
377 double t294, t295, t297, t299, t300, t302, t303, t307;
378 double t311, t319, t323, t332, t337, t338, t345, t346;
379 double t351, t352, t353, t354, t358, t360, t366, t372;
380 double t385, t387, t391, t392, t396, t401, t402, t403;
381 double t408, t414, t417, t418, t419, t423, t431, t432;
382 double t436, t441, t442, t447, t453;
383
384 #ifndef XC_DONT_COMPILE_LXC
385 double t469, t476, t481, t486, t488, t495, t503, t504;
386 double t509, t519, t521, t526, t534, t548, t553, t565;
387 double t604, t612, t616, t621, t628, t634, t661, t700;
388 double t761, t765, t818;
389 #endif
390
391 #endif
392
393 #endif
394
395 #endif
396
397 #endif
398
399
400 gga_x_dk87_params *params;
401
402 assert(p->params != NULL);
403 params = (gga_x_dk87_params * )(p->params);
404
405 t1 = M_CBRT3;
406 t2 = 0.1e1 / M_PI;
407 t3 = POW_1_3(t2);
408 t4 = t1 * t3;
409 t5 = M_CBRT4;
410 t6 = t5 * t5;
411 t7 = POW_1_3(rho[0]);
412 t8 = t6 * t7;
413 t9 = M_CBRT6;
414 t10 = t9 * t9;
415 t11 = t2 * t10;
416 t12 = M_PI * M_PI;
417 t13 = POW_1_3(t12);
418 t14 = 0.1e1 / t13;
419 t15 = t1 * t1;
420 t16 = t14 * t15;
421 t17 = 0.1e1 / t3;
422 t19 = t11 * t16 * t17;
423 t20 = t5 * sigma[0];
424 t21 = rho[0] * rho[0];
425 t22 = t7 * t7;
426 t24 = 0.1e1 / t22 / t21;
427 t25 = sqrt(sigma[0]);
428 t26 = t7 * rho[0];
429 t29 = pow(t25 / t26, params->alpha);
430 t31 = params->a1 * t29 + 0.1e1;
431 t33 = params->b1 * sigma[0];
432 t35 = t33 * t24 + 0.1e1;
433 t36 = 0.1e1 / t35;
434 t37 = t24 * t31 * t36;
435 t41 = 0.1e1 + 0.7e1 / 0.11664e5 * t19 * t20 * t37;
436 t43 = t4 * t8 * t41;
437 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
438 zk[0] = -0.3e1 / 0.8e1 * t43;
439
440 #ifndef XC_DONT_COMPILE_VXC
441
442 if(order < 1) return;
443
444
445 t46 = t26 * t1;
446 t47 = t3 * t6;
447 t48 = t21 * rho[0];
448 t50 = 0.1e1 / t22 / t48;
449 t52 = t50 * t31 * t36;
450 t59 = t11 * t14 * t15 * t17 * t5;
451 t62 = t29 * params->alpha;
452 t63 = t62 * t36;
453 t67 = sigma[0] * sigma[0];
454 t68 = t5 * t67;
455 t69 = t21 * t21;
456 t70 = t69 * t21;
457 t72 = 0.1e1 / t7 / t70;
458 t74 = t35 * t35;
459 t75 = 0.1e1 / t74;
460 t77 = t31 * t75 * params->b1;
461 t81 = -0.7e1 / 0.4374e4 * t19 * t20 * t52 - 0.7e1 / 0.8748e4 * t59 * sigma[0] * t50 * params->a1 * t63 + 0.7e1 / 0.4374e4 * t19 * t68 * t72 * t77;
462 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
463 vrho[0] = -t43 / 0.2e1 - 0.3e1 / 0.8e1 * t46 * t47 * t81;
464
465 t85 = t11 * t16;
466 t86 = t17 * t5;
467 t95 = t69 * rho[0];
468 t97 = 0.1e1 / t7 / t95;
469 t103 = t47 * (0.7e1 / 0.11664e5 * t85 * t86 * t37 + 0.7e1 / 0.23328e5 * t19 * t5 * t24 * params->a1 * t63 - 0.7e1 / 0.11664e5 * t19 * t20 * t97 * t77);
470 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
471 vsigma[0] = -0.3e1 / 0.8e1 * t46 * t103;
472
473 #ifndef XC_DONT_COMPILE_FXC
474
475 if(order < 2) return;
476
477
478 t106 = 0.1e1 / t22;
479 t107 = t6 * t106;
480 t114 = 0.1e1 / t22 / t69;
481 t116 = t114 * t31 * t36;
482 t121 = sigma[0] * t114 * params->a1;
483 t125 = t69 * t48;
484 t127 = 0.1e1 / t7 / t125;
485 t132 = params->alpha * params->alpha;
486 t133 = t29 * t132;
487 t134 = t133 * t36;
488 t140 = t75 * params->b1;
489 t141 = t62 * t140;
490 t145 = t67 * sigma[0];
491 t146 = t5 * t145;
492 t147 = t69 * t69;
493 t148 = t147 * t21;
494 t149 = 0.1e1 / t148;
495 t152 = 0.1e1 / t74 / t35;
496 t154 = params->b1 * params->b1;
497 t155 = t31 * t152 * t154;
498 t159 = 0.77e2 / 0.13122e5 * t19 * t20 * t116 + 0.133e3 / 0.26244e5 * t59 * t121 * t63 - 0.7e1 / 0.486e3 * t19 * t68 * t127 * t77 + 0.7e1 / 0.6561e4 * t59 * t121 * t134 - 0.28e2 / 0.6561e4 * t59 * t67 * t127 * params->a1 * t141 + 0.56e2 / 0.6561e4 * t19 * t146 * t149 * t155;
499 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
500 v2rho2[0] = -t4 * t107 * t41 / 0.6e1 - t4 * t8 * t81 - 0.3e1 / 0.8e1 * t46 * t47 * t159;
501
502 t163 = t7 * t1;
503 t170 = t5 * t50 * params->a1;
504 t174 = t5 * t72;
505 t176 = t140 * sigma[0];
506 t184 = t72 * params->a1 * t29;
507 t185 = params->alpha * t75;
508 t186 = t185 * t33;
509 t190 = t147 * rho[0];
510 t191 = 0.1e1 / t190;
511 t197 = t47 * (-0.7e1 / 0.4374e4 * t85 * t86 * t52 - 0.7e1 / 0.4374e4 * t19 * t170 * t63 + 0.7e1 / 0.1458e4 * t19 * t174 * t31 * t176 - 0.7e1 / 0.17496e5 * t19 * t170 * t134 + 0.7e1 / 0.4374e4 * t59 * t184 * t186 - 0.7e1 / 0.2187e4 * t19 * t68 * t191 * t155);
512 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
513 v2rhosigma[0] = -t163 * t103 / 0.2e1 - 0.3e1 / 0.8e1 * t46 * t197;
514
515 t201 = t24 * params->a1 * t29;
516 t202 = 0.1e1 / sigma[0];
517 t203 = params->alpha * t202;
518 t204 = t203 * t36;
519 t212 = t132 * t202;
520 t213 = t212 * t36;
521 t218 = t97 * params->a1 * t29;
522 t219 = t185 * params->b1;
523 t223 = 0.1e1 / t147;
524 t229 = t47 * (0.7e1 / 0.23328e5 * t59 * t201 * t204 - 0.7e1 / 0.5832e4 * t19 * t5 * t97 * t77 + 0.7e1 / 0.46656e5 * t59 * t201 * t213 - 0.7e1 / 0.11664e5 * t59 * t218 * t219 + 0.7e1 / 0.5832e4 * t19 * t20 * t223 * t155);
525 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
526 v2sigma2[0] = -0.3e1 / 0.8e1 * t46 * t229;
527
528 #ifndef XC_DONT_COMPILE_KXC
529
530 if(order < 3) return;
531
532
533 t233 = 0.1e1 / t22 / rho[0];
534 t234 = t6 * t233;
535 t245 = 0.1e1 / t22 / t95;
536 t247 = t245 * t31 * t36;
537 t252 = sigma[0] * t245 * params->a1;
538 t257 = 0.1e1 / t7 / t147;
539 t266 = t67 * t257 * params->a1;
540 t270 = t147 * t48;
541 t271 = 0.1e1 / t270;
542 t276 = t132 * params->alpha;
543 t277 = t29 * t276;
544 t278 = t277 * t36;
545 t282 = t133 * t140;
546 t288 = t152 * t154;
547 t289 = t62 * t288;
548 t293 = t67 * t67;
549 t294 = t5 * t293;
550 t295 = t147 * t95;
551 t297 = 0.1e1 / t22 / t295;
552 t299 = t74 * t74;
553 t300 = 0.1e1 / t299;
554 t302 = t154 * params->b1;
555 t303 = t31 * t300 * t302;
556 t307 = -0.539e3 / 0.19683e5 * t19 * t20 * t247 - 0.413e3 / 0.13122e5 * t59 * t252 * t63 + 0.2387e4 / 0.19683e5 * t19 * t68 * t257 * t77 - 0.77e2 / 0.6561e4 * t59 * t252 * t134 + 0.140e3 / 0.2187e4 * t59 * t266 * t141 - 0.1064e4 / 0.6561e4 * t19 * t146 * t271 * t155 - 0.28e2 / 0.19683e5 * t59 * t252 * t278 + 0.56e2 / 0.6561e4 * t59 * t266 * t282 - 0.224e3 / 0.6561e4 * t59 * t145 * t271 * params->a1 * t289 + 0.448e3 / 0.6561e4 * t19 * t294 * t297 * t303;
557 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
558 v3rho3[0] = t4 * t234 * t41 / 0.9e1 - t4 * t107 * t81 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t159 - 0.3e1 / 0.8e1 * t46 * t47 * t307;
559
560 t311 = t106 * t1;
561 t319 = t5 * t114 * params->a1;
562 t323 = t5 * t127;
563 t332 = t127 * params->a1 * t29;
564 t337 = t5 * t149 * t31;
565 t338 = t288 * t67;
566 t345 = t132 * t75;
567 t346 = t345 * t33;
568 t351 = t149 * params->a1 * t29;
569 t352 = params->alpha * t152;
570 t353 = t154 * t67;
571 t354 = t352 * t353;
572 t358 = t147 * t69;
573 t360 = 0.1e1 / t22 / t358;
574 t366 = t47 * (0.77e2 / 0.13122e5 * t85 * t86 * t116 + 0.35e2 / 0.4374e4 * t19 * t319 * t63 - 0.455e3 / 0.13122e5 * t19 * t323 * t31 * t176 + 0.7e1 / 0.1944e4 * t19 * t319 * t134 - 0.91e2 / 0.4374e4 * t59 * t332 * t186 + 0.119e3 / 0.2187e4 * t19 * t337 * t338 + 0.7e1 / 0.13122e5 * t19 * t319 * t278 - 0.7e1 / 0.2187e4 * t59 * t332 * t346 + 0.28e2 / 0.2187e4 * t59 * t351 * t354 - 0.56e2 / 0.2187e4 * t19 * t146 * t360 * t303);
575 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
576 v3rho2sigma[0] = -t311 * t103 / 0.6e1 - t163 * t197 - 0.3e1 / 0.8e1 * t46 * t366;
577
578 t372 = t50 * params->a1 * t29;
579 t385 = t5 * t191;
580 t387 = t288 * sigma[0];
581 t391 = t276 * t202;
582 t392 = t391 * t36;
583 t396 = t345 * params->b1;
584 t401 = t191 * params->a1 * t29;
585 t402 = t154 * sigma[0];
586 t403 = t352 * t402;
587 t408 = 0.1e1 / t22 / t270;
588 t414 = t47 * (-0.7e1 / 0.8748e4 * t59 * t372 * t204 - 0.7e1 / 0.8748e4 * t59 * t372 * t213 + 0.49e2 / 0.8748e4 * t59 * t184 * t219 + 0.14e2 / 0.2187e4 * t19 * t174 * t77 - 0.35e2 / 0.2187e4 * t19 * t385 * t31 * t387 - 0.7e1 / 0.34992e5 * t59 * t372 * t392 + 0.7e1 / 0.5832e4 * t59 * t184 * t396 - 0.7e1 / 0.1458e4 * t59 * t401 * t403 + 0.7e1 / 0.729e3 * t19 * t68 * t408 * t303);
589 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
590 v3rhosigma2[0] = -t163 * t229 / 0.2e1 - 0.3e1 / 0.8e1 * t46 * t414;
591
592 t417 = 0.1e1 / t67;
593 t418 = params->alpha * t417;
594 t419 = t418 * t36;
595 t423 = t203 * t140;
596 t431 = t276 * t417;
597 t432 = t431 * t36;
598 t436 = t212 * t140;
599 t441 = t223 * params->a1 * t29;
600 t442 = t352 * t154;
601 t447 = 0.1e1 / t22 / t148;
602 t453 = t47 * (-0.7e1 / 0.23328e5 * t59 * t201 * t419 - 0.7e1 / 0.7776e4 * t59 * t218 * t423 + 0.7e1 / 0.1944e4 * t19 * t5 * t223 * t155 + 0.7e1 / 0.93312e5 * t59 * t201 * t432 - 0.7e1 / 0.15552e5 * t59 * t218 * t436 + 0.7e1 / 0.3888e4 * t59 * t441 * t442 - 0.7e1 / 0.1944e4 * t19 * t20 * t447 * t303);
603 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
604 v3sigma3[0] = -0.3e1 / 0.8e1 * t46 * t453;
605
606 #ifndef XC_DONT_COMPILE_LXC
607
608 if(order < 4) return;
609
610
611 t469 = 0.1e1 / t22 / t70;
612 t476 = 0.1e1 / t7 / t190;
613 t481 = 0.1e1 / t358;
614 t486 = t147 * t70;
615 t488 = 0.1e1 / t22 / t486;
616 t495 = t147 * t147;
617 t503 = t154 * t154;
618 t504 = t31 / t299 / t35 * t503;
619 t509 = sigma[0] * t469 * params->a1;
620 t519 = t132 * t132;
621 t521 = t29 * t519 * t36;
622 t526 = t67 * t476 * params->a1;
623 t534 = t145 * t481 * params->a1;
624 t548 = t300 * t302;
625 t553 = 0.9163e4 / 0.59049e5 * t19 * t20 * t469 * t31 * t36 - 0.21329e5 / 0.19683e5 * t19 * t68 * t476 * t77 + 0.143528e6 / 0.59049e5 * t19 * t146 * t481 * t155 - 0.43904e5 / 0.19683e5 * t19 * t294 * t488 * t303 + 0.14336e5 / 0.19683e5 * t19 * t5 * t293 * sigma[0] / t7 / t495 / rho[0] * t504 + 0.25375e5 / 0.118098e6 * t59 * t509 * t63 + 0.2135e4 / 0.19683e5 * t59 * t509 * t134 + 0.1400e4 / 0.59049e5 * t59 * t509 * t278 + 0.112e3 / 0.59049e5 * t59 * t509 * t521 - 0.46004e5 / 0.59049e5 * t59 * t526 * t141 - 0.1232e4 / 0.6561e4 * t59 * t526 * t282 + 0.18368e5 / 0.19683e5 * t59 * t534 * t289 - 0.896e3 / 0.59049e5 * t59 * t526 * t277 * t140 + 0.1792e4 / 0.19683e5 * t59 * t534 * t133 * t288 - 0.7168e4 / 0.19683e5 * t59 * t293 * t488 * params->a1 * t62 * t548;
626 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
627 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t24 * t41 + 0.4e1 / 0.9e1 * t4 * t234 * t81 - t4 * t107 * t159 - 0.2e1 * t4 * t8 * t307 - 0.3e1 / 0.8e1 * t46 * t47 * t553;
628
629 t565 = t5 * t245 * params->a1;
630 t604 = t257 * params->a1 * t29;
631 t612 = t271 * params->a1 * t29;
632 t616 = t276 * t75;
633 t621 = t132 * t152;
634 t628 = params->alpha * t300;
635 t634 = -0.889e3 / 0.19683e5 * t19 * t565 * t63 + 0.1771e4 / 0.6561e4 * t19 * t5 * t257 * t31 * t176 - 0.721e3 / 0.26244e5 * t19 * t565 * t134 - 0.14350e5 / 0.19683e5 * t19 * t5 * t271 * t31 * t338 - 0.287e3 / 0.39366e5 * t19 * t565 * t278 + 0.4984e4 / 0.6561e4 * t19 * t5 * t297 * t31 * t548 * t145 - 0.14e2 / 0.19683e5 * t19 * t565 * t521 - 0.1792e4 / 0.6561e4 * t19 * t294 / t7 / t495 * t504 - 0.539e3 / 0.19683e5 * t85 * t86 * t247 + 0.4333e4 / 0.19683e5 * t59 * t604 * t186 + 0.133e3 / 0.2187e4 * t59 * t604 * t346 - 0.2044e4 / 0.6561e4 * t59 * t612 * t354 + 0.112e3 / 0.19683e5 * t59 * t604 * t616 * t33 - 0.224e3 / 0.6561e4 * t59 * t612 * t621 * t353 + 0.896e3 / 0.6561e4 * t59 * t297 * params->a1 * t29 * t628 * t302 * t145;
636 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
637 v4rho3sigma[0] = t233 * t1 * t103 / 0.9e1 - t311 * t197 / 0.2e1 - 0.3e1 / 0.2e1 * t163 * t366 - 0.3e1 / 0.8e1 * t46 * t47 * t634;
638
639 t661 = t114 * params->a1 * t29;
640 t700 = -0.266e3 / 0.6561e4 * t19 * t323 * t77 + 0.1169e4 / 0.6561e4 * t19 * t337 * t387 - 0.175e3 / 0.729e3 * t19 * t5 * t360 * t31 * t548 * t67 + 0.224e3 / 0.2187e4 * t19 * t146 / t7 / t147 / t125 * t504 + 0.35e2 / 0.8748e4 * t59 * t661 * t213 - 0.1211e4 / 0.26244e5 * t59 * t332 * t219 + 0.7e1 / 0.3888e4 * t59 * t661 * t392 - 0.301e3 / 0.17496e5 * t59 * t332 * t396 + 0.7e1 / 0.26244e5 * t59 * t661 * t519 * t202 * t36 - 0.14e2 / 0.6561e4 * t59 * t332 * t616 * params->b1 + 0.77e2 / 0.26244e5 * t59 * t661 * t204 + 0.413e3 / 0.4374e4 * t59 * t351 * t403 + 0.28e2 / 0.2187e4 * t59 * t351 * t621 * t402 - 0.112e3 / 0.2187e4 * t59 * t360 * params->a1 * t29 * t628 * t302 * t67;
641 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
642 v4rho2sigma2[0] = -t311 * t229 / 0.6e1 - t163 * t414 - 0.3e1 / 0.8e1 * t46 * t47 * t700;
643
644 t761 = 0.7e1 / 0.8748e4 * t59 * t372 * t419 + 0.7e1 / 0.17496e5 * t59 * t372 * t132 * t417 * t36 + 0.35e2 / 0.8748e4 * t59 * t184 * t423 + 0.7e1 / 0.1944e4 * t59 * t184 * t436 - 0.35e2 / 0.1458e4 * t59 * t401 * t442 - 0.7e1 / 0.243e3 * t19 * t385 * t155 + 0.49e2 / 0.729e3 * t19 * t5 * t408 * t31 * t548 * sigma[0] - 0.7e1 / 0.34992e5 * t59 * t372 * t432 - 0.7e1 / 0.69984e5 * t59 * t372 * t519 * t417 * t36 + 0.7e1 / 0.8748e4 * t59 * t184 * t391 * t140 - 0.7e1 / 0.1458e4 * t59 * t401 * t621 * t154 + 0.14e2 / 0.729e3 * t59 * t408 * params->a1 * t29 * t628 * t302 * sigma[0] - 0.28e2 / 0.729e3 * t19 * t68 / t7 / t486 * t504;
645 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
646 v4rhosigma3[0] = -t163 * t453 / 0.2e1 - 0.3e1 / 0.8e1 * t46 * t47 * t761;
647
648 t765 = 0.1e1 / t145;
649 t818 = -0.7e1 / 0.46656e5 * t59 * t201 * t132 * t765 * t36 + 0.7e1 / 0.11664e5 * t59 * t201 * params->alpha * t765 * t36 + 0.7e1 / 0.5832e4 * t59 * t218 * t418 * t140 + 0.7e1 / 0.1944e4 * t59 * t441 * t203 * t288 - 0.7e1 / 0.486e3 * t19 * t5 * t447 * t303 + 0.7e1 / 0.186624e6 * t59 * t201 * t519 * t765 * t36 - 0.7e1 / 0.46656e5 * t59 * t201 * t276 * t765 * t36 - 0.7e1 / 0.23328e5 * t59 * t218 * t431 * t140 + 0.7e1 / 0.3888e4 * t59 * t441 * t212 * t288 - 0.7e1 / 0.972e3 * t59 * t447 * params->a1 * t29 * t628 * t302 + 0.7e1 / 0.486e3 * t19 * t20 / t7 / t295 * t504;
650 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
651 v4sigma4[0] = -0.3e1 / 0.8e1 * t46 * t47 * t818;
652
653 #ifndef XC_DONT_COMPILE_MXC
654
655 if(order < 5) return;
656
657
658 #endif
659
660 #endif
661
662 #endif
663
664 #endif
665
666 #endif
667
668
669 }
670
671
672 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 *))673 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
674 {
675
676 #ifndef XC_DONT_COMPILE_EXC
677 double t1, t2, t3, t4, t5, t6, t7, t8;
678 double t9, t10, t11, t13, t14, t15, t16, t17;
679 double t18, t19, t20, t21, t22, t23, t24, t25;
680 double t26, t28, t29, t30, t31, t32, t34, t35;
681 double t39, t41, t43, t45, t46, t47, t51, t53;
682 double t55, t56, t57, t58, t59, t60, t61, t62;
683 double t64, t65, t69, t71, t73, t75, t76, t77;
684 double t81, t83;
685
686 #ifndef XC_DONT_COMPILE_VXC
687 double t85, t86, t87, t88, t89, t90, t92, t93;
688 double t95, t96, t97, t98, t99, t101, t102, t103;
689 double t105, t107, t114, t117, t118, t122, t123, t124;
690 double t125, t127, t129, t130, t132, t136, t138, t139;
691 double t140, t141, t142, t144, t145, t146, t148, t149;
692 double t153, t154, t156, t157, t158, t159, t161, t162;
693 double t163, t165, t167, t173, t174, t178, t179, t180;
694 double t181, t183, t185, t186, t188, t192, t194, t195;
695 double t199, t200, t201, t202, t203, t212, t214, t219;
696 double t223, t232, t234, t239;
697
698 #ifndef XC_DONT_COMPILE_FXC
699 double t243, t245, t246, t247, t248, t249, t250, t252;
700 double t253, t254, t256, t257, t258, t260, t261, t262;
701 double t263, t264, t265, t267, t268, t270, t271, t273;
702 double t274, t276, t277, t279, t281, t286, t290, t292;
703 double t297, t298, t299, t305, t306, t310, t311, t312;
704 double t313, t314, t317, t319, t320, t324, t326, t327;
705 double t328, t329, t330, t331, t332, t334, t335, t336;
706 double t338, t339, t340, t341, t343, t344, t345, t347;
707 double t348, t349, t351, t352, t353, t354, t356, t359;
708 double t361, t363, t366, t367, t368, t370, t374, t375;
709 double t376, t377, t379, t382, t384, t385, t387, t390;
710 double t392, t395, t397, t400, t401, t403, t404, t405;
711 double t406, t407, t409, t410, t411, t412, t414, t415;
712 double t416, t417, t419, t420, t421, t423, t424, t425;
713 double t427, t429, t434, t438, t440, t445, t446, t452;
714 double t453, t457, t458, t459, t460, t461, t464, t466;
715 double t470, t472, t473, t474, t477, t479, t480, t488;
716 double t492, t494, t502, t503, t504, t508, t509, t514;
717 double t519, t520, t521, t526, t531, t538, t542, t544;
718 double t552, t553, t554, t558, t559, t564, t565, t569;
719 double t570, t571, t572, t580, t581, t586, t587, t591;
720 double t596, t597, t601, t602, t603, t604, t612, t613;
721 double t618, t619, t623, t628, t629;
722
723 #ifndef XC_DONT_COMPILE_KXC
724 double t636, t638, t642, t643, t644, t645, t647, t648;
725 double t649, t651, t653, t655, t656, t657, t659, t660;
726 double t662, t663, t664, t665, t667, t668, t669, t671;
727 double t673, t675, t676, t677, t679, t680, t681, t683;
728 double t684, t686, t689, t691, t693, t694, t695, t697;
729 double t698, t699, t700, t701, t703, t704, t706, t707;
730 double t708, t710, t711, t713, t715, t716, t718, t719;
731 double t720, t722, t723, t725, t726, t728, t731, t733;
732 double t738, t743, t752, t756, t757, t762, t763, t764;
733 double t768, t774, t775, t779, t780, t781, t783, t785;
734 double t786, t788, t789, t793, t795, t796, t797, t802;
735 double t803, t806, t807, t808, t810, t812, t814, t815;
736 double t816, t817, t819, t821, t822, t823, t825, t827;
737 double t829, t830, t831, t833, t835, t837, t839, t841;
738 double t843, t844, t846, t851, t854, t857, t859, t868;
739 double t871, t872, t874, t876, t878, t880, t882, t884;
740 double t887, t888, t890, t891, t892, t894, t896, t898;
741 double t900, t902, t904, t907, t912, t914, t916, t918;
742 double t920, t922, t924, t926, t928, t930, t932, t934;
743 double t936, t938, t943, t947, t953, t956, t959, t962;
744 double t965, t967, t969, t971, t973, t976, t977, t979;
745 double t981, t983, t984, t985, t987, t989, t991, t993;
746 double t996, t1004, t1005, t1010, t1011, t1013, t1014, t1015;
747 double t1017, t1018, t1019, t1020, t1022, t1023, t1025, t1027;
748 double t1028, t1029, t1032, t1033, t1035, t1036, t1037, t1039;
749 double t1040, t1041, t1043, t1044, t1045, t1047, t1048, t1049;
750 double t1050, t1052, t1053, t1055, t1057, t1062, t1067, t1076;
751 double t1080, t1081, t1086, t1087, t1091, t1097, t1098, t1102;
752 double t1103, t1104, t1106, t1108, t1109, t1111, t1115, t1117;
753 double t1118, t1119, t1123, t1127, t1129, t1130, t1134, t1137;
754 double t1145, t1149, t1158, t1163, t1164, t1171, t1172, t1177;
755 double t1178, t1179, t1180, t1184, t1186, t1191, t1196, t1200;
756 double t1201, t1205, t1211, t1212, t1213, t1215, t1216, t1217;
757 double t1218, t1222, t1225, t1226, t1228, t1232, t1234, t1235;
758 double t1239, t1242, t1245, t1250, t1254, t1259, t1263, t1266;
759 double t1273, t1277, t1286, t1291, t1292, t1299, t1300, t1305;
760 double t1306, t1307, t1308, t1312, t1314, t1319, t1320, t1323;
761 double t1326, t1330, t1343, t1345, t1349, t1350, t1354, t1359;
762 double t1360, t1361, t1366, t1371, t1372, t1375, t1377, t1381;
763 double t1386, t1390, t1403, t1405, t1409, t1410, t1414, t1419;
764 double t1420, t1421, t1426, t1431, t1432, t1435, t1436, t1437;
765 double t1441, t1449, t1450, t1454, t1459, t1460, t1465, t1470;
766 double t1471, t1474, t1475, t1476, t1480, t1488, t1489, t1493;
767 double t1498, t1499, t1504, t1509, t1510;
768
769 #ifndef XC_DONT_COMPILE_LXC
770 double t1524, t1528, t1531, t1534, t1536, t1538, t1540, t1542;
771 double t1548, t1555, t1566, t1576, t1579, t1582, t1585, t1592;
772 double t1597, t1602, t1604, t1611, t1619, t1620, t1625, t1633;
773 double t1647, t1653, t1663, t1665, t1669, t1676, t1680, t1686;
774 double t1688, t1690, t1693, t1696, t1702, t1705, t1708, t1722;
775 double t1724, t1726, t1732, t1739, t1741, t1742, t1743, t1748;
776 double t1750, t1752, t1754, t1756, t1762, t1764, t1766, t1767;
777 double t1768, t1773, t1776, t1786, t1789, t1792, t1794, t1799;
778 double t1806, t1809, t1810, t1812, t1813, t1815, t1823, t1826;
779 double t1829, t1835, t1836, t1863, t1864, t1865, t1867, t1879;
780 double t1880, t1886, t1888, t1890, t1894, t1901, t1907, t1909;
781 double t1912, t1918, t1925, t1931, t1938, t1940, t1946, t1949;
782 double t1962, t1967, t1969, t1972, t1978, t1982, t1989, t2000;
783 double t2003, t2006, t2012, t2015, t2022, t2028, t2030, t2036;
784 double t2039, t2042, t2050, t2053, t2063, t2105, t2108, t2113;
785 double t2126, t2133, t2140, t2150, t2153, t2160, t2169, t2184;
786 double t2191, t2207, t2217, t2232, t2241, t2253, t2263, t2266;
787 double t2269, t2272, t2275, t2278, t2281, t2304, t2319, t2374;
788 double t2380, t2382, t2395, t2398, t2406, t2409, t2416, t2421;
789 double t2434, t2442, t2447, t2452, t2462, t2487, t2500, t2532;
790 double t2539, t2544, t2549, t2551, t2558, t2566, t2571, t2579;
791 double t2593, t2599, t2610, t2614, t2618, t2621, t2626, t2629;
792 double t2632, t2637, t2676, t2677, t2682, t2689, t2701, t2706;
793 double t2709, t2713, t2731, t2734, t2741, t2744, t2747, t2751;
794 double t2764, t2765, t2767, t2769, t2771, t2772, t2786, t2790;
795 double t2797, t2810, t2813, t2814, t2816, t2829, t2832, t2845;
796 double t2857, t2874, t2882, t2891, t2900, t2932, t2971, t2979;
797 double t2983, t2988, t2995, t3001, t3024, t3025, t3027, t3059;
798 double t3098, t3103, t3104, t3116, t3131, t3141, t3188, t3227;
799 double t3231, t3292, t3296, t3365, t3369, t3422, t3426, t3479;
800 #endif
801
802 #endif
803
804 #endif
805
806 #endif
807
808 #endif
809
810
811 gga_x_dk87_params *params;
812
813 assert(p->params != NULL);
814 params = (gga_x_dk87_params * )(p->params);
815
816 t1 = M_CBRT3;
817 t2 = 0.1e1 / M_PI;
818 t3 = POW_1_3(t2);
819 t4 = t1 * t3;
820 t5 = M_CBRT4;
821 t6 = t5 * t5;
822 t7 = t4 * t6;
823 t8 = rho[0] - rho[1];
824 t9 = rho[0] + rho[1];
825 t10 = 0.1e1 / t9;
826 t11 = t8 * t10;
827 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
828 t14 = POW_1_3(t13);
829 t15 = t14 * t13;
830 t16 = POW_1_3(t9);
831 t17 = t15 * t16;
832 t18 = M_CBRT6;
833 t19 = t18 * t18;
834 t20 = t2 * t19;
835 t21 = M_PI * M_PI;
836 t22 = POW_1_3(t21);
837 t23 = 0.1e1 / t22;
838 t24 = t1 * t1;
839 t25 = t23 * t24;
840 t26 = 0.1e1 / t3;
841 t28 = t20 * t25 * t26;
842 t29 = t5 * sigma[0];
843 t30 = rho[0] * rho[0];
844 t31 = POW_1_3(rho[0]);
845 t32 = t31 * t31;
846 t34 = 0.1e1 / t32 / t30;
847 t35 = sqrt(sigma[0]);
848 t39 = pow(t35 / t31 / rho[0], params->alpha);
849 t41 = params->a1 * t39 + 0.1e1;
850 t43 = params->b1 * sigma[0];
851 t45 = t43 * t34 + 0.1e1;
852 t46 = 0.1e1 / t45;
853 t47 = t34 * t41 * t46;
854 t51 = 0.1e1 + 0.7e1 / 0.11664e5 * t28 * t29 * t47;
855 t53 = t7 * t17 * t51;
856 t55 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
857 t56 = POW_1_3(t55);
858 t57 = t56 * t55;
859 t58 = t57 * t16;
860 t59 = t5 * sigma[2];
861 t60 = rho[1] * rho[1];
862 t61 = POW_1_3(rho[1]);
863 t62 = t61 * t61;
864 t64 = 0.1e1 / t62 / t60;
865 t65 = sqrt(sigma[2]);
866 t69 = pow(t65 / t61 / rho[1], params->alpha);
867 t71 = params->a1 * t69 + 0.1e1;
868 t73 = params->b1 * sigma[2];
869 t75 = t73 * t64 + 0.1e1;
870 t76 = 0.1e1 / t75;
871 t77 = t64 * t71 * t76;
872 t81 = 0.1e1 + 0.7e1 / 0.11664e5 * t28 * t59 * t77;
873 t83 = t7 * t58 * t81;
874 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
875 zk[0] = -0.3e1 / 0.8e1 * t53 - 0.3e1 / 0.8e1 * t83;
876
877 #ifndef XC_DONT_COMPILE_VXC
878
879 if(order < 1) return;
880
881
882 t85 = 0.3e1 / 0.8e1 * t53;
883 t86 = 0.3e1 / 0.8e1 * t83;
884 t87 = t14 * t16;
885 t88 = t9 * t9;
886 t89 = 0.1e1 / t88;
887 t90 = t8 * t89;
888 t92 = t10 / 0.2e1 - t90 / 0.2e1;
889 t93 = t51 * t92;
890 t95 = t7 * t87 * t93;
891 t96 = t95 / 0.2e1;
892 t97 = t16 * t16;
893 t98 = 0.1e1 / t97;
894 t99 = t15 * t98;
895 t101 = t7 * t99 * t51;
896 t102 = t101 / 0.8e1;
897 t103 = t30 * rho[0];
898 t105 = 0.1e1 / t32 / t103;
899 t107 = t105 * t41 * t46;
900 t114 = t20 * t23 * t24 * t26 * t5;
901 t117 = t39 * params->alpha;
902 t118 = t117 * t46;
903 t122 = sigma[0] * sigma[0];
904 t123 = t5 * t122;
905 t124 = t30 * t30;
906 t125 = t124 * t30;
907 t127 = 0.1e1 / t31 / t125;
908 t129 = t45 * t45;
909 t130 = 0.1e1 / t129;
910 t132 = t41 * t130 * params->b1;
911 t136 = -0.7e1 / 0.4374e4 * t28 * t29 * t107 - 0.7e1 / 0.8748e4 * t114 * sigma[0] * t105 * params->a1 * t118 + 0.7e1 / 0.4374e4 * t28 * t123 * t127 * t132;
912 t138 = t7 * t17 * t136;
913 t139 = 0.3e1 / 0.8e1 * t138;
914 t140 = t56 * t16;
915 t141 = -t92;
916 t142 = t81 * t141;
917 t144 = t7 * t140 * t142;
918 t145 = t144 / 0.2e1;
919 t146 = t57 * t98;
920 t148 = t7 * t146 * t81;
921 t149 = t148 / 0.8e1;
922 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
923 vrho[0] = -t85 - t86 + t9 * (-t96 - t102 - t139 - t145 - t149);
924
925 t153 = -t10 / 0.2e1 - t90 / 0.2e1;
926 t154 = t51 * t153;
927 t156 = t7 * t87 * t154;
928 t157 = t156 / 0.2e1;
929 t158 = -t153;
930 t159 = t81 * t158;
931 t161 = t7 * t140 * t159;
932 t162 = t161 / 0.2e1;
933 t163 = t60 * rho[1];
934 t165 = 0.1e1 / t62 / t163;
935 t167 = t165 * t71 * t76;
936 t173 = t69 * params->alpha;
937 t174 = t173 * t76;
938 t178 = sigma[2] * sigma[2];
939 t179 = t5 * t178;
940 t180 = t60 * t60;
941 t181 = t180 * t60;
942 t183 = 0.1e1 / t61 / t181;
943 t185 = t75 * t75;
944 t186 = 0.1e1 / t185;
945 t188 = t71 * t186 * params->b1;
946 t192 = -0.7e1 / 0.4374e4 * t28 * t59 * t167 - 0.7e1 / 0.8748e4 * t114 * sigma[2] * t165 * params->a1 * t174 + 0.7e1 / 0.4374e4 * t28 * t179 * t183 * t188;
947 t194 = t7 * t58 * t192;
948 t195 = 0.3e1 / 0.8e1 * t194;
949 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
950 vrho[1] = -t85 - t86 + t9 * (-t157 - t102 - t162 - t149 - t195);
951
952 t199 = t16 * t9 * t1;
953 t200 = t199 * t3;
954 t201 = t6 * t15;
955 t202 = t20 * t25;
956 t203 = t26 * t5;
957 t212 = t124 * rho[0];
958 t214 = 0.1e1 / t31 / t212;
959 t219 = 0.7e1 / 0.11664e5 * t202 * t203 * t47 + 0.7e1 / 0.23328e5 * t28 * t5 * t34 * params->a1 * t118 - 0.7e1 / 0.11664e5 * t28 * t29 * t214 * t132;
960 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
961 vsigma[0] = -0.3e1 / 0.8e1 * t200 * t201 * t219;
962
963 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
964 vsigma[1] = 0.0e0;
965
966 t223 = t6 * t57;
967 t232 = t180 * rho[1];
968 t234 = 0.1e1 / t61 / t232;
969 t239 = 0.7e1 / 0.11664e5 * t202 * t203 * t77 + 0.7e1 / 0.23328e5 * t28 * t5 * t64 * params->a1 * t174 - 0.7e1 / 0.11664e5 * t28 * t59 * t234 * t188;
970 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
971 vsigma[2] = -0.3e1 / 0.8e1 * t200 * t223 * t239;
972
973 #ifndef XC_DONT_COMPILE_FXC
974
975 if(order < 2) return;
976
977
978 t243 = t101 / 0.4e1;
979 t245 = t148 / 0.4e1;
980 t246 = t14 * t14;
981 t247 = 0.1e1 / t246;
982 t248 = t247 * t16;
983 t249 = t92 * t92;
984 t250 = t51 * t249;
985 t252 = t7 * t248 * t250;
986 t253 = t252 / 0.6e1;
987 t254 = t14 * t98;
988 t256 = t7 * t254 * t93;
989 t257 = t256 / 0.3e1;
990 t258 = t136 * t92;
991 t260 = t7 * t87 * t258;
992 t261 = t88 * t9;
993 t262 = 0.1e1 / t261;
994 t263 = t8 * t262;
995 t264 = -t89 + t263;
996 t265 = t51 * t264;
997 t267 = t7 * t87 * t265;
998 t268 = t267 / 0.2e1;
999 t270 = 0.1e1 / t97 / t9;
1000 t271 = t15 * t270;
1001 t273 = t7 * t271 * t51;
1002 t274 = t273 / 0.12e2;
1003 t276 = t7 * t99 * t136;
1004 t277 = t276 / 0.4e1;
1005 t279 = 0.1e1 / t32 / t124;
1006 t281 = t279 * t41 * t46;
1007 t286 = sigma[0] * t279 * params->a1;
1008 t290 = t124 * t103;
1009 t292 = 0.1e1 / t31 / t290;
1010 t297 = params->alpha * params->alpha;
1011 t298 = t39 * t297;
1012 t299 = t298 * t46;
1013 t305 = t130 * params->b1;
1014 t306 = t117 * t305;
1015 t310 = t122 * sigma[0];
1016 t311 = t5 * t310;
1017 t312 = t124 * t124;
1018 t313 = t312 * t30;
1019 t314 = 0.1e1 / t313;
1020 t317 = 0.1e1 / t129 / t45;
1021 t319 = params->b1 * params->b1;
1022 t320 = t41 * t317 * t319;
1023 t324 = 0.77e2 / 0.13122e5 * t28 * t29 * t281 + 0.133e3 / 0.26244e5 * t114 * t286 * t118 - 0.7e1 / 0.486e3 * t28 * t123 * t292 * t132 + 0.7e1 / 0.6561e4 * t114 * t286 * t299 - 0.28e2 / 0.6561e4 * t114 * t122 * t292 * params->a1 * t306 + 0.56e2 / 0.6561e4 * t28 * t311 * t314 * t320;
1024 t326 = t7 * t17 * t324;
1025 t327 = 0.3e1 / 0.8e1 * t326;
1026 t328 = t56 * t56;
1027 t329 = 0.1e1 / t328;
1028 t330 = t329 * t16;
1029 t331 = t141 * t141;
1030 t332 = t81 * t331;
1031 t334 = t7 * t330 * t332;
1032 t335 = t334 / 0.6e1;
1033 t336 = t56 * t98;
1034 t338 = t7 * t336 * t142;
1035 t339 = t338 / 0.3e1;
1036 t340 = -t264;
1037 t341 = t81 * t340;
1038 t343 = t7 * t140 * t341;
1039 t344 = t343 / 0.2e1;
1040 t345 = t57 * t270;
1041 t347 = t7 * t345 * t81;
1042 t348 = t347 / 0.12e2;
1043 t349 = -t253 - t257 - t260 - t268 + t274 - t277 - t327 - t335 - t339 - t344 + t348;
1044 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1045 v2rho2[0] = -t95 - t243 - 0.3e1 / 0.4e1 * t138 - t144 - t245 + t9 * t349;
1046
1047 t351 = t6 * t247;
1048 t352 = t4 * t351;
1049 t353 = t16 * t51;
1050 t354 = t153 * t92;
1051 t356 = t352 * t353 * t354;
1052 t359 = t7 * t254 * t154;
1053 t361 = t136 * t153;
1054 t363 = t7 * t87 * t361;
1055 t366 = 0.1e1 / t97 / t88;
1056 t367 = t14 * t366;
1057 t368 = t51 * t8;
1058 t370 = t7 * t367 * t368;
1059 t374 = t6 * t329;
1060 t375 = t4 * t374;
1061 t376 = t16 * t81;
1062 t377 = t158 * t141;
1063 t379 = t375 * t376 * t377;
1064 t382 = t7 * t336 * t159;
1065 t384 = t56 * t366;
1066 t385 = t81 * t8;
1067 t387 = t7 * t384 * t385;
1068 t390 = t192 * t141;
1069 t392 = t7 * t140 * t390;
1070 t395 = t7 * t146 * t192;
1071 t397 = -t356 / 0.6e1 - t359 / 0.6e1 - t363 / 0.2e1 - t370 / 0.2e1 - t256 / 0.6e1 + t274 - t276 / 0.8e1 - t379 / 0.6e1 - t382 / 0.6e1 + t387 / 0.2e1 - t338 / 0.6e1 + t348 - t392 / 0.2e1 - t395 / 0.8e1;
1072 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1073 v2rho2[1] = t9 * t397 - t139 - t145 - t157 - t162 - t195 - t243 - t245 - t96;
1074
1075 t400 = t153 * t153;
1076 t401 = t51 * t400;
1077 t403 = t7 * t248 * t401;
1078 t404 = t403 / 0.6e1;
1079 t405 = t359 / 0.3e1;
1080 t406 = t89 + t263;
1081 t407 = t51 * t406;
1082 t409 = t7 * t87 * t407;
1083 t410 = t409 / 0.2e1;
1084 t411 = t158 * t158;
1085 t412 = t81 * t411;
1086 t414 = t7 * t330 * t412;
1087 t415 = t414 / 0.6e1;
1088 t416 = t382 / 0.3e1;
1089 t417 = t192 * t158;
1090 t419 = t7 * t140 * t417;
1091 t420 = -t406;
1092 t421 = t81 * t420;
1093 t423 = t7 * t140 * t421;
1094 t424 = t423 / 0.2e1;
1095 t425 = t395 / 0.4e1;
1096 t427 = 0.1e1 / t62 / t180;
1097 t429 = t427 * t71 * t76;
1098 t434 = sigma[2] * t427 * params->a1;
1099 t438 = t180 * t163;
1100 t440 = 0.1e1 / t61 / t438;
1101 t445 = t69 * t297;
1102 t446 = t445 * t76;
1103 t452 = t186 * params->b1;
1104 t453 = t173 * t452;
1105 t457 = t178 * sigma[2];
1106 t458 = t5 * t457;
1107 t459 = t180 * t180;
1108 t460 = t459 * t60;
1109 t461 = 0.1e1 / t460;
1110 t464 = 0.1e1 / t185 / t75;
1111 t466 = t71 * t464 * t319;
1112 t470 = 0.77e2 / 0.13122e5 * t28 * t59 * t429 + 0.133e3 / 0.26244e5 * t114 * t434 * t174 - 0.7e1 / 0.486e3 * t28 * t179 * t440 * t188 + 0.7e1 / 0.6561e4 * t114 * t434 * t446 - 0.28e2 / 0.6561e4 * t114 * t178 * t440 * params->a1 * t453 + 0.56e2 / 0.6561e4 * t28 * t458 * t461 * t466;
1113 t472 = t7 * t58 * t470;
1114 t473 = 0.3e1 / 0.8e1 * t472;
1115 t474 = -t404 - t405 - t410 + t274 - t415 - t416 - t419 - t424 + t348 - t425 - t473;
1116 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1117 v2rho2[2] = -t156 - t243 - t161 - t245 - 0.3e1 / 0.4e1 * t194 + t9 * t474;
1118
1119 t477 = t7 * t17 * t219;
1120 t479 = t6 * t14;
1121 t480 = t219 * t92;
1122 t488 = t5 * t105 * params->a1;
1123 t492 = t5 * t127;
1124 t494 = t305 * sigma[0];
1125 t502 = t127 * params->a1 * t39;
1126 t503 = params->alpha * t130;
1127 t504 = t503 * t43;
1128 t508 = t312 * rho[0];
1129 t509 = 0.1e1 / t508;
1130 t514 = -0.7e1 / 0.4374e4 * t202 * t203 * t107 - 0.7e1 / 0.4374e4 * t28 * t488 * t118 + 0.7e1 / 0.1458e4 * t28 * t492 * t41 * t494 - 0.7e1 / 0.17496e5 * t28 * t488 * t299 + 0.7e1 / 0.4374e4 * t114 * t502 * t504 - 0.7e1 / 0.2187e4 * t28 * t123 * t509 * t320;
1131 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1132 v2rhosigma[0] = -t477 / 0.2e1 - t200 * t479 * t480 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t201 * t514;
1133
1134 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1135 v2rhosigma[1] = 0.0e0;
1136
1137 t519 = t7 * t58 * t239;
1138 t520 = t6 * t56;
1139 t521 = t239 * t141;
1140 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1141 v2rhosigma[2] = -t200 * t520 * t521 / 0.2e1 - t519 / 0.2e1;
1142
1143 t526 = t479 * t219 * t153;
1144 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1145 v2rhosigma[3] = -t200 * t526 / 0.2e1 - t477 / 0.2e1;
1146
1147 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1148 v2rhosigma[4] = 0.0e0;
1149
1150 t531 = t520 * t239 * t158;
1151 t538 = t5 * t165 * params->a1;
1152 t542 = t5 * t183;
1153 t544 = t452 * sigma[2];
1154 t552 = t183 * params->a1 * t69;
1155 t553 = params->alpha * t186;
1156 t554 = t553 * t73;
1157 t558 = t459 * rho[1];
1158 t559 = 0.1e1 / t558;
1159 t564 = -0.7e1 / 0.4374e4 * t202 * t203 * t167 - 0.7e1 / 0.4374e4 * t28 * t538 * t174 + 0.7e1 / 0.1458e4 * t28 * t542 * t71 * t544 - 0.7e1 / 0.17496e5 * t28 * t538 * t446 + 0.7e1 / 0.4374e4 * t114 * t552 * t554 - 0.7e1 / 0.2187e4 * t28 * t179 * t559 * t466;
1160 t565 = t223 * t564;
1161 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1162 v2rhosigma[5] = -t519 / 0.2e1 - t200 * t531 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t565;
1163
1164 t569 = t34 * params->a1 * t39;
1165 t570 = 0.1e1 / sigma[0];
1166 t571 = params->alpha * t570;
1167 t572 = t571 * t46;
1168 t580 = t297 * t570;
1169 t581 = t580 * t46;
1170 t586 = t214 * params->a1 * t39;
1171 t587 = t503 * params->b1;
1172 t591 = 0.1e1 / t312;
1173 t596 = 0.7e1 / 0.23328e5 * t114 * t569 * t572 - 0.7e1 / 0.5832e4 * t28 * t5 * t214 * t132 + 0.7e1 / 0.46656e5 * t114 * t569 * t581 - 0.7e1 / 0.11664e5 * t114 * t586 * t587 + 0.7e1 / 0.5832e4 * t28 * t29 * t591 * t320;
1174 t597 = t201 * t596;
1175 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1176 v2sigma2[0] = -0.3e1 / 0.8e1 * t200 * t597;
1177
1178 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1179 v2sigma2[1] = 0.0e0;
1180
1181 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1182 v2sigma2[2] = 0.0e0;
1183
1184 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1185 v2sigma2[3] = 0.0e0;
1186
1187 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1188 v2sigma2[4] = 0.0e0;
1189
1190 t601 = t64 * params->a1 * t69;
1191 t602 = 0.1e1 / sigma[2];
1192 t603 = params->alpha * t602;
1193 t604 = t603 * t76;
1194 t612 = t297 * t602;
1195 t613 = t612 * t76;
1196 t618 = t234 * params->a1 * t69;
1197 t619 = t553 * params->b1;
1198 t623 = 0.1e1 / t459;
1199 t628 = 0.7e1 / 0.23328e5 * t114 * t601 * t604 - 0.7e1 / 0.5832e4 * t28 * t5 * t234 * t188 + 0.7e1 / 0.46656e5 * t114 * t601 * t613 - 0.7e1 / 0.11664e5 * t114 * t618 * t619 + 0.7e1 / 0.5832e4 * t28 * t59 * t623 * t466;
1200 t629 = t223 * t628;
1201 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1202 v2sigma2[5] = -0.3e1 / 0.8e1 * t200 * t629;
1203
1204 #ifndef XC_DONT_COMPILE_KXC
1205
1206 if(order < 3) return;
1207
1208
1209 t636 = t273 / 0.4e1;
1210 t638 = t347 / 0.4e1;
1211 t642 = 0.1e1 / t246 / t13;
1212 t643 = t642 * t16;
1213 t644 = t249 * t92;
1214 t645 = t51 * t644;
1215 t647 = t7 * t643 * t645;
1216 t648 = t647 / 0.9e1;
1217 t649 = t247 * t98;
1218 t651 = t7 * t649 * t250;
1219 t653 = t92 * t264;
1220 t655 = t352 * t353 * t653;
1221 t656 = t655 / 0.2e1;
1222 t657 = t14 * t270;
1223 t659 = t7 * t657 * t93;
1224 t660 = t659 / 0.3e1;
1225 t662 = 0.1e1 / t328 / t55;
1226 t663 = t662 * t16;
1227 t664 = t331 * t141;
1228 t665 = t81 * t664;
1229 t667 = t7 * t663 * t665;
1230 t668 = t667 / 0.9e1;
1231 t669 = t329 * t98;
1232 t671 = t7 * t669 * t332;
1233 t673 = t141 * t340;
1234 t675 = t375 * t376 * t673;
1235 t676 = t675 / 0.2e1;
1236 t677 = t56 * t270;
1237 t679 = t7 * t677 * t142;
1238 t680 = t679 / 0.3e1;
1239 t681 = t136 * t249;
1240 t683 = t7 * t248 * t681;
1241 t684 = t683 / 0.2e1;
1242 t686 = t7 * t254 * t258;
1243 t689 = t7 * t254 * t265;
1244 t691 = t324 * t92;
1245 t693 = t7 * t87 * t691;
1246 t694 = 0.3e1 / 0.2e1 * t693;
1247 t695 = t136 * t264;
1248 t697 = t7 * t87 * t695;
1249 t698 = 0.3e1 / 0.2e1 * t697;
1250 t699 = t88 * t88;
1251 t700 = 0.1e1 / t699;
1252 t701 = t8 * t700;
1253 t703 = 0.3e1 * t262 - 0.3e1 * t701;
1254 t704 = t51 * t703;
1255 t706 = t7 * t87 * t704;
1256 t707 = t706 / 0.2e1;
1257 t708 = t15 * t366;
1258 t710 = t7 * t708 * t51;
1259 t711 = 0.5e1 / 0.36e2 * t710;
1260 t713 = t7 * t336 * t341;
1261 t715 = -t703;
1262 t716 = t81 * t715;
1263 t718 = t7 * t140 * t716;
1264 t719 = t718 / 0.2e1;
1265 t720 = t57 * t366;
1266 t722 = t7 * t720 * t81;
1267 t723 = 0.5e1 / 0.36e2 * t722;
1268 t725 = t7 * t271 * t136;
1269 t726 = t725 / 0.4e1;
1270 t728 = t7 * t99 * t324;
1271 t731 = 0.1e1 / t32 / t212;
1272 t733 = t731 * t41 * t46;
1273 t738 = sigma[0] * t731 * params->a1;
1274 t743 = 0.1e1 / t31 / t312;
1275 t752 = t122 * t743 * params->a1;
1276 t756 = t312 * t103;
1277 t757 = 0.1e1 / t756;
1278 t762 = t297 * params->alpha;
1279 t763 = t39 * t762;
1280 t764 = t763 * t46;
1281 t768 = t298 * t305;
1282 t774 = t317 * t319;
1283 t775 = t117 * t774;
1284 t779 = t122 * t122;
1285 t780 = t5 * t779;
1286 t781 = t312 * t212;
1287 t783 = 0.1e1 / t32 / t781;
1288 t785 = t129 * t129;
1289 t786 = 0.1e1 / t785;
1290 t788 = t319 * params->b1;
1291 t789 = t41 * t786 * t788;
1292 t793 = -0.539e3 / 0.19683e5 * t28 * t29 * t733 - 0.413e3 / 0.13122e5 * t114 * t738 * t118 + 0.2387e4 / 0.19683e5 * t28 * t123 * t743 * t132 - 0.77e2 / 0.6561e4 * t114 * t738 * t299 + 0.140e3 / 0.2187e4 * t114 * t752 * t306 - 0.1064e4 / 0.6561e4 * t28 * t311 * t757 * t320 - 0.28e2 / 0.19683e5 * t114 * t738 * t764 + 0.56e2 / 0.6561e4 * t114 * t752 * t768 - 0.224e3 / 0.6561e4 * t114 * t310 * t757 * params->a1 * t775 + 0.448e3 / 0.6561e4 * t28 * t780 * t783 * t789;
1293 t795 = t7 * t17 * t793;
1294 t796 = 0.3e1 / 0.8e1 * t795;
1295 t797 = -t689 / 0.2e1 - t694 - t698 - t707 - t711 - t713 / 0.2e1 - t719 - t723 + t726 - 0.3e1 / 0.8e1 * t728 - t796;
1296 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1297 v3rho3[0] = -t256 - t252 / 0.2e1 - t334 / 0.2e1 - t338 - 0.3e1 * t260 - 0.3e1 / 0.2e1 * t267 + t636 - 0.3e1 / 0.2e1 * t343 + t638 - 0.3e1 / 0.4e1 * t276 - 0.9e1 / 0.8e1 * t326 + t9 * (t648 - t651 / 0.6e1 - t656 + t660 + t668 - t671 / 0.6e1 - t676 + t680 - t684 - t686 + t797);
1298
1299 t802 = t356 / 0.3e1;
1300 t803 = t379 / 0.3e1;
1301 t806 = t6 * t642;
1302 t807 = t4 * t806;
1303 t808 = t153 * t249;
1304 t810 = t807 * t353 * t808;
1305 t812 = t98 * t51;
1306 t814 = t352 * t812 * t354;
1307 t815 = t814 / 0.9e1;
1308 t816 = t366 * t51;
1309 t817 = t8 * t92;
1310 t819 = t352 * t816 * t817;
1311 t821 = t6 * t662;
1312 t822 = t4 * t821;
1313 t823 = t158 * t331;
1314 t825 = t822 * t376 * t823;
1315 t827 = t98 * t81;
1316 t829 = t375 * t827 * t377;
1317 t830 = t829 / 0.9e1;
1318 t831 = t16 * t136;
1319 t833 = t352 * t831 * t354;
1320 t835 = t153 * t264;
1321 t837 = t352 * t353 * t835;
1322 t839 = t158 * t340;
1323 t841 = t375 * t376 * t839;
1324 t843 = t366 * t81;
1325 t844 = t8 * t141;
1326 t846 = t375 * t843 * t844;
1327 t851 = t7 * t367 * t51;
1328 t854 = t7 * t384 * t81;
1329 t857 = t7 * t345 * t192;
1330 t859 = t810 / 0.9e1 - t815 - t819 / 0.3e1 + t825 / 0.9e1 - t830 - t833 / 0.3e1 - t837 / 0.6e1 - t841 / 0.6e1 + t846 / 0.3e1 - t711 - t723 + t725 / 0.6e1 - t728 / 0.8e1 - t851 / 0.2e1 + t854 / 0.2e1 + t857 / 0.12e2;
1331 t868 = t7 * t657 * t154;
1332 t871 = 0.1e1 / t97 / t261;
1333 t872 = t14 * t871;
1334 t874 = t7 * t872 * t368;
1335 t876 = t136 * t8;
1336 t878 = t7 * t367 * t876;
1337 t880 = t7 * t677 * t159;
1338 t882 = t192 * t331;
1339 t884 = t7 * t330 * t882;
1340 t887 = t7 * t336 * t390;
1341 t888 = t887 / 0.3e1;
1342 t890 = t7 * t254 * t361;
1343 t891 = t890 / 0.3e1;
1344 t892 = t324 * t153;
1345 t894 = t7 * t87 * t892;
1346 t896 = t56 * t871;
1347 t898 = t7 * t896 * t385;
1348 t900 = t192 * t340;
1349 t902 = t7 * t140 * t900;
1350 t904 = -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 + t868 / 0.9e1 + 0.7e1 / 0.6e1 * t874 - t878 + t880 / 0.9e1 - t884 / 0.6e1 - t888 - t891 - t894 / 0.2e1 - 0.7e1 / 0.6e1 * t898 - t902 / 0.2e1;
1351 t907 = -t268 + t636 - t344 + t638 - t363 + t387 - t392 - t425 - t276 / 0.2e1 - t327 + t9 * (t859 + t904);
1352 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1353 v3rho3[1] = -0.2e1 / 0.3e1 * t256 - t253 - t335 - 0.2e1 / 0.3e1 * t338 - t802 - t405 - t370 - t803 - t416 - t260 + t907;
1354
1355 t912 = t400 * t92;
1356 t914 = t807 * t353 * t912;
1357 t916 = t406 * t92;
1358 t918 = t352 * t353 * t916;
1359 t920 = t411 * t141;
1360 t922 = t822 * t376 * t920;
1361 t924 = t158 * t8;
1362 t926 = t375 * t843 * t924;
1363 t928 = t16 * t192;
1364 t930 = t375 * t928 * t377;
1365 t932 = t420 * t141;
1366 t934 = t375 * t376 * t932;
1367 t936 = t153 * t8;
1368 t938 = t352 * t816 * t936;
1369 t943 = t7 * t146 * t470;
1370 t947 = t914 / 0.9e1 - t918 / 0.6e1 + t922 / 0.9e1 + t926 / 0.3e1 - t930 / 0.3e1 - t934 / 0.6e1 - t938 / 0.3e1 - t815 - t830 - t711 - t723 + t725 / 0.12e2 + t857 / 0.6e1 - t943 / 0.8e1 + t659 / 0.9e1 + t679 / 0.9e1;
1371 t953 = t7 * t649 * t401;
1372 t956 = t7 * t254 * t407;
1373 t959 = t7 * t669 * t412;
1374 t962 = t7 * t336 * t417;
1375 t965 = t7 * t336 * t421;
1376 t967 = t136 * t400;
1377 t969 = t7 * t248 * t967;
1378 t971 = t136 * t406;
1379 t973 = t7 * t87 * t971;
1380 t976 = -t262 - 0.3e1 * t701;
1381 t977 = t51 * t976;
1382 t979 = t7 * t87 * t977;
1383 t981 = t192 * t8;
1384 t983 = t7 * t384 * t981;
1385 t984 = -t976;
1386 t985 = t81 * t984;
1387 t987 = t7 * t140 * t985;
1388 t989 = t470 * t141;
1389 t991 = t7 * t140 * t989;
1390 t993 = 0.2e1 / 0.9e1 * t868 - t874 / 0.3e1 + 0.2e1 / 0.9e1 * t880 - t888 - t891 + t898 / 0.3e1 - t953 / 0.18e2 - t956 / 0.6e1 - t959 / 0.18e2 - t962 / 0.3e1 - t965 / 0.6e1 - t969 / 0.6e1 - t973 / 0.2e1 - t979 / 0.2e1 + t983 - t987 / 0.2e1 - t991 / 0.2e1;
1391 t996 = t638 - t363 + t387 - t392 - t395 / 0.2e1 - t410 - t419 - t424 - t277 - t473 + t9 * (t947 + t993);
1392 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1393 v3rho3[2] = -t257 - t339 - t802 - 0.2e1 / 0.3e1 * t359 - t370 - t803 - 0.2e1 / 0.3e1 * t382 - t404 - t415 + t636 + t996;
1394
1395 t1004 = t868 / 0.3e1;
1396 t1005 = t880 / 0.3e1;
1397 t1010 = t400 * t153;
1398 t1011 = t51 * t1010;
1399 t1013 = t7 * t643 * t1011;
1400 t1014 = t1013 / 0.9e1;
1401 t1015 = t153 * t406;
1402 t1017 = t352 * t353 * t1015;
1403 t1018 = t1017 / 0.2e1;
1404 t1019 = t411 * t158;
1405 t1020 = t81 * t1019;
1406 t1022 = t7 * t663 * t1020;
1407 t1023 = t1022 / 0.9e1;
1408 t1025 = t158 * t420;
1409 t1027 = t375 * t376 * t1025;
1410 t1028 = t1027 / 0.2e1;
1411 t1029 = t857 / 0.4e1;
1412 t1032 = -0.3e1 * t262 - 0.3e1 * t701;
1413 t1033 = t51 * t1032;
1414 t1035 = t7 * t87 * t1033;
1415 t1036 = t1035 / 0.2e1;
1416 t1037 = t192 * t411;
1417 t1039 = t7 * t330 * t1037;
1418 t1040 = t1039 / 0.2e1;
1419 t1041 = t470 * t158;
1420 t1043 = t7 * t140 * t1041;
1421 t1044 = 0.3e1 / 0.2e1 * t1043;
1422 t1045 = t192 * t420;
1423 t1047 = t7 * t140 * t1045;
1424 t1048 = 0.3e1 / 0.2e1 * t1047;
1425 t1049 = -t1032;
1426 t1050 = t81 * t1049;
1427 t1052 = t7 * t140 * t1050;
1428 t1053 = t1052 / 0.2e1;
1429 t1055 = 0.1e1 / t62 / t232;
1430 t1057 = t1055 * t71 * t76;
1431 t1062 = sigma[2] * t1055 * params->a1;
1432 t1067 = 0.1e1 / t61 / t459;
1433 t1076 = t178 * t1067 * params->a1;
1434 t1080 = t459 * t163;
1435 t1081 = 0.1e1 / t1080;
1436 t1086 = t69 * t762;
1437 t1087 = t1086 * t76;
1438 t1091 = t445 * t452;
1439 t1097 = t464 * t319;
1440 t1098 = t173 * t1097;
1441 t1102 = t178 * t178;
1442 t1103 = t5 * t1102;
1443 t1104 = t459 * t232;
1444 t1106 = 0.1e1 / t62 / t1104;
1445 t1108 = t185 * t185;
1446 t1109 = 0.1e1 / t1108;
1447 t1111 = t71 * t1109 * t788;
1448 t1115 = -0.539e3 / 0.19683e5 * t28 * t59 * t1057 - 0.413e3 / 0.13122e5 * t114 * t1062 * t174 + 0.2387e4 / 0.19683e5 * t28 * t179 * t1067 * t188 - 0.77e2 / 0.6561e4 * t114 * t1062 * t446 + 0.140e3 / 0.2187e4 * t114 * t1076 * t453 - 0.1064e4 / 0.6561e4 * t28 * t458 * t1081 * t466 - 0.28e2 / 0.19683e5 * t114 * t1062 * t1087 + 0.56e2 / 0.6561e4 * t114 * t1076 * t1091 - 0.224e3 / 0.6561e4 * t114 * t457 * t1081 * params->a1 * t1098 + 0.448e3 / 0.6561e4 * t28 * t1103 * t1106 * t1111;
1449 t1117 = t7 * t58 * t1115;
1450 t1118 = 0.3e1 / 0.8e1 * t1117;
1451 t1119 = -t1028 - t711 - t723 + t1029 - 0.3e1 / 0.8e1 * t943 - t1036 - t1040 - t1044 - t1048 - t1053 - t1118;
1452 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1453 v3rho3[3] = -t359 - t382 - t403 / 0.2e1 - t414 / 0.2e1 + t636 + t638 - 0.3e1 / 0.4e1 * t395 - 0.3e1 / 0.2e1 * t409 - 0.3e1 * t419 - 0.3e1 / 0.2e1 * t423 - 0.9e1 / 0.8e1 * t472 + t9 * (t1004 + t1005 - t953 / 0.6e1 - t956 / 0.2e1 - t959 / 0.6e1 - t962 - t965 / 0.2e1 + t1014 - t1018 + t1023 + t1119);
1454
1455 t1123 = t7 * t87 * t480;
1456 t1127 = t7 * t99 * t219 / 0.6e1;
1457 t1129 = t7 * t17 * t514;
1458 t1130 = t219 * t249;
1459 t1134 = t514 * t92;
1460 t1137 = t219 * t264;
1461 t1145 = t5 * t279 * params->a1;
1462 t1149 = t5 * t292;
1463 t1158 = t292 * params->a1 * t39;
1464 t1163 = t5 * t314 * t41;
1465 t1164 = t774 * t122;
1466 t1171 = t297 * t130;
1467 t1172 = t1171 * t43;
1468 t1177 = t314 * params->a1 * t39;
1469 t1178 = params->alpha * t317;
1470 t1179 = t319 * t122;
1471 t1180 = t1178 * t1179;
1472 t1184 = t312 * t124;
1473 t1186 = 0.1e1 / t32 / t1184;
1474 t1191 = 0.77e2 / 0.13122e5 * t202 * t203 * t281 + 0.35e2 / 0.4374e4 * t28 * t1145 * t118 - 0.455e3 / 0.13122e5 * t28 * t1149 * t41 * t494 + 0.7e1 / 0.1944e4 * t28 * t1145 * t299 - 0.91e2 / 0.4374e4 * t114 * t1158 * t504 + 0.119e3 / 0.2187e4 * t28 * t1163 * t1164 + 0.7e1 / 0.13122e5 * t28 * t1145 * t764 - 0.7e1 / 0.2187e4 * t114 * t1158 * t1172 + 0.28e2 / 0.2187e4 * t114 * t1177 * t1180 - 0.56e2 / 0.2187e4 * t28 * t311 * t1186 * t789;
1475 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1476 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t1123 - t1127 - t1129 - t200 * t351 * t1130 / 0.6e1 - t200 * t479 * t1134 - t200 * t479 * t1137 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t201 * t1191;
1477
1478 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1479 v3rho2sigma[1] = 0.0e0;
1480
1481 t1196 = t7 * t140 * t521;
1482 t1200 = t7 * t146 * t239 / 0.6e1;
1483 t1201 = t239 * t331;
1484 t1205 = t239 * t340;
1485 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1486 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t1196 - t1200 - t200 * t374 * t1201 / 0.6e1 - t200 * t520 * t1205 / 0.2e1;
1487
1488 t1211 = t16 * t1;
1489 t1212 = t1211 * t3;
1490 t1213 = t1212 * t526;
1491 t1215 = t3 * t6;
1492 t1216 = t199 * t1215;
1493 t1217 = t247 * t219;
1494 t1218 = t1217 * t354;
1495 t1222 = t479 * t514 * t153;
1496 t1225 = t270 * t1;
1497 t1226 = t1225 * t3;
1498 t1228 = t479 * t219 * t8;
1499 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1500 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t1123 - t1127 - t1129 / 0.2e1 - 0.2e1 / 0.3e1 * t1213 - t1216 * t1218 / 0.6e1 - t200 * t1222 / 0.2e1 - t1226 * t1228 / 0.2e1;
1501
1502 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1503 v3rho2sigma[4] = 0.0e0;
1504
1505 t1232 = t1212 * t531;
1506 t1234 = t329 * t239;
1507 t1235 = t1234 * t377;
1508 t1239 = t520 * t239 * t8;
1509 t1242 = t1212 * t565;
1510 t1245 = t520 * t564 * t141;
1511 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1512 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t1196 - t1200 - 0.2e1 / 0.3e1 * t1232 - t1216 * t1235 / 0.6e1 + t1226 * t1239 / 0.2e1 - t1242 / 0.2e1 - t200 * t1245 / 0.2e1;
1513
1514 t1250 = t351 * t219 * t400;
1515 t1254 = t479 * t219 * t406;
1516 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1517 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t1213 - t1127 - t200 * t1250 / 0.6e1 - t200 * t1254 / 0.2e1;
1518
1519 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1520 v3rho2sigma[7] = 0.0e0;
1521
1522 t1259 = t374 * t239 * t411;
1523 t1263 = t520 * t564 * t158;
1524 t1266 = t520 * t239 * t420;
1525 t1273 = t5 * t427 * params->a1;
1526 t1277 = t5 * t440;
1527 t1286 = t440 * params->a1 * t69;
1528 t1291 = t5 * t461 * t71;
1529 t1292 = t1097 * t178;
1530 t1299 = t297 * t186;
1531 t1300 = t1299 * t73;
1532 t1305 = t461 * params->a1 * t69;
1533 t1306 = params->alpha * t464;
1534 t1307 = t319 * t178;
1535 t1308 = t1306 * t1307;
1536 t1312 = t459 * t180;
1537 t1314 = 0.1e1 / t62 / t1312;
1538 t1319 = 0.77e2 / 0.13122e5 * t202 * t203 * t429 + 0.35e2 / 0.4374e4 * t28 * t1273 * t174 - 0.455e3 / 0.13122e5 * t28 * t1277 * t71 * t544 + 0.7e1 / 0.1944e4 * t28 * t1273 * t446 - 0.91e2 / 0.4374e4 * t114 * t1286 * t554 + 0.119e3 / 0.2187e4 * t28 * t1291 * t1292 + 0.7e1 / 0.13122e5 * t28 * t1273 * t1087 - 0.7e1 / 0.2187e4 * t114 * t1286 * t1300 + 0.28e2 / 0.2187e4 * t114 * t1305 * t1308 - 0.56e2 / 0.2187e4 * t28 * t458 * t1314 * t1111;
1539 t1320 = t223 * t1319;
1540 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1541 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t1232 - t1200 - t1242 - t200 * t1259 / 0.6e1 - t200 * t1263 - t200 * t1266 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t1320;
1542
1543 t1323 = t1212 * t597;
1544 t1326 = t479 * t596 * t92;
1545 t1330 = t105 * params->a1 * t39;
1546 t1343 = t5 * t509;
1547 t1345 = t774 * sigma[0];
1548 t1349 = t762 * t570;
1549 t1350 = t1349 * t46;
1550 t1354 = t1171 * params->b1;
1551 t1359 = t509 * params->a1 * t39;
1552 t1360 = t319 * sigma[0];
1553 t1361 = t1178 * t1360;
1554 t1366 = 0.1e1 / t32 / t756;
1555 t1371 = -0.7e1 / 0.8748e4 * t114 * t1330 * t572 - 0.7e1 / 0.8748e4 * t114 * t1330 * t581 + 0.49e2 / 0.8748e4 * t114 * t502 * t587 + 0.14e2 / 0.2187e4 * t28 * t492 * t132 - 0.35e2 / 0.2187e4 * t28 * t1343 * t41 * t1345 - 0.7e1 / 0.34992e5 * t114 * t1330 * t1350 + 0.7e1 / 0.5832e4 * t114 * t502 * t1354 - 0.7e1 / 0.1458e4 * t114 * t1359 * t1361 + 0.7e1 / 0.729e3 * t28 * t123 * t1366 * t789;
1556 t1372 = t201 * t1371;
1557 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1558 v3rhosigma2[0] = -t1323 / 0.2e1 - t200 * t1326 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t1372;
1559
1560 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1561 v3rhosigma2[1] = 0.0e0;
1562
1563 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1564 v3rhosigma2[2] = 0.0e0;
1565
1566 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1567 v3rhosigma2[3] = 0.0e0;
1568
1569 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1570 v3rhosigma2[4] = 0.0e0;
1571
1572 t1375 = t1212 * t629;
1573 t1377 = t520 * t628 * t141;
1574 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1575 v3rhosigma2[5] = -t200 * t1377 / 0.2e1 - t1375 / 0.2e1;
1576
1577 t1381 = t479 * t596 * t153;
1578 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1579 v3rhosigma2[6] = -t200 * t1381 / 0.2e1 - t1323 / 0.2e1;
1580
1581 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1582 v3rhosigma2[7] = 0.0e0;
1583
1584 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1585 v3rhosigma2[8] = 0.0e0;
1586
1587 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1588 v3rhosigma2[9] = 0.0e0;
1589
1590 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1591 v3rhosigma2[10] = 0.0e0;
1592
1593 t1386 = t520 * t628 * t158;
1594 t1390 = t165 * params->a1 * t69;
1595 t1403 = t5 * t559;
1596 t1405 = t1097 * sigma[2];
1597 t1409 = t762 * t602;
1598 t1410 = t1409 * t76;
1599 t1414 = t1299 * params->b1;
1600 t1419 = t559 * params->a1 * t69;
1601 t1420 = t319 * sigma[2];
1602 t1421 = t1306 * t1420;
1603 t1426 = 0.1e1 / t62 / t1080;
1604 t1431 = -0.7e1 / 0.8748e4 * t114 * t1390 * t604 - 0.7e1 / 0.8748e4 * t114 * t1390 * t613 + 0.49e2 / 0.8748e4 * t114 * t552 * t619 + 0.14e2 / 0.2187e4 * t28 * t542 * t188 - 0.35e2 / 0.2187e4 * t28 * t1403 * t71 * t1405 - 0.7e1 / 0.34992e5 * t114 * t1390 * t1410 + 0.7e1 / 0.5832e4 * t114 * t552 * t1414 - 0.7e1 / 0.1458e4 * t114 * t1419 * t1421 + 0.7e1 / 0.729e3 * t28 * t179 * t1426 * t1111;
1605 t1432 = t223 * t1431;
1606 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1607 v3rhosigma2[11] = -t1375 / 0.2e1 - t200 * t1386 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t1432;
1608
1609 t1435 = 0.1e1 / t122;
1610 t1436 = params->alpha * t1435;
1611 t1437 = t1436 * t46;
1612 t1441 = t571 * t305;
1613 t1449 = t762 * t1435;
1614 t1450 = t1449 * t46;
1615 t1454 = t580 * t305;
1616 t1459 = t591 * params->a1 * t39;
1617 t1460 = t1178 * t319;
1618 t1465 = 0.1e1 / t32 / t313;
1619 t1470 = -0.7e1 / 0.23328e5 * t114 * t569 * t1437 - 0.7e1 / 0.7776e4 * t114 * t586 * t1441 + 0.7e1 / 0.1944e4 * t28 * t5 * t591 * t320 + 0.7e1 / 0.93312e5 * t114 * t569 * t1450 - 0.7e1 / 0.15552e5 * t114 * t586 * t1454 + 0.7e1 / 0.3888e4 * t114 * t1459 * t1460 - 0.7e1 / 0.1944e4 * t28 * t29 * t1465 * t789;
1620 t1471 = t201 * t1470;
1621 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1622 v3sigma3[0] = -0.3e1 / 0.8e1 * t200 * t1471;
1623
1624 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1625 v3sigma3[1] = 0.0e0;
1626
1627 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1628 v3sigma3[2] = 0.0e0;
1629
1630 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1631 v3sigma3[3] = 0.0e0;
1632
1633 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1634 v3sigma3[4] = 0.0e0;
1635
1636 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1637 v3sigma3[5] = 0.0e0;
1638
1639 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1640 v3sigma3[6] = 0.0e0;
1641
1642 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1643 v3sigma3[7] = 0.0e0;
1644
1645 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1646 v3sigma3[8] = 0.0e0;
1647
1648 t1474 = 0.1e1 / t178;
1649 t1475 = params->alpha * t1474;
1650 t1476 = t1475 * t76;
1651 t1480 = t603 * t452;
1652 t1488 = t762 * t1474;
1653 t1489 = t1488 * t76;
1654 t1493 = t612 * t452;
1655 t1498 = t623 * params->a1 * t69;
1656 t1499 = t1306 * t319;
1657 t1504 = 0.1e1 / t62 / t460;
1658 t1509 = -0.7e1 / 0.23328e5 * t114 * t601 * t1476 - 0.7e1 / 0.7776e4 * t114 * t618 * t1480 + 0.7e1 / 0.1944e4 * t28 * t5 * t623 * t466 + 0.7e1 / 0.93312e5 * t114 * t601 * t1489 - 0.7e1 / 0.15552e5 * t114 * t618 * t1493 + 0.7e1 / 0.3888e4 * t114 * t1498 * t1499 - 0.7e1 / 0.1944e4 * t28 * t59 * t1504 * t1111;
1659 t1510 = t223 * t1509;
1660 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1661 v3sigma3[9] = -0.3e1 / 0.8e1 * t200 * t1510;
1662
1663 #ifndef XC_DONT_COMPILE_LXC
1664
1665 if(order < 4) return;
1666
1667
1668 t1524 = 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;
1669 t1528 = 0.5e1 / 0.9e1 * t710;
1670 t1531 = 0.5e1 / 0.9e1 * t722;
1671 t1534 = t662 * t98;
1672 t1536 = t7 * t1534 * t665;
1673 t1538 = t329 * t270;
1674 t1540 = t7 * t1538 * t332;
1675 t1542 = t340 * t340;
1676 t1548 = t7 * t384 * t142;
1677 t1555 = t375 * t827 * t673;
1678 t1566 = t352 * t812 * t653;
1679 t1576 = t7 * t708 * t136;
1680 t1579 = t7 * t271 * t324;
1681 t1582 = t7 * t99 * t793;
1682 t1585 = 0.1e1 / t32 / t125;
1683 t1592 = 0.1e1 / t31 / t508;
1684 t1597 = 0.1e1 / t1184;
1685 t1602 = t312 * t125;
1686 t1604 = 0.1e1 / t32 / t1602;
1687 t1611 = t312 * t312;
1688 t1619 = t319 * t319;
1689 t1620 = t41 / t785 / t45 * t1619;
1690 t1625 = t122 * t1592 * params->a1;
1691 t1633 = t310 * t1597 * params->a1;
1692 t1647 = t786 * t788;
1693 t1653 = sigma[0] * t1585 * params->a1;
1694 t1663 = t297 * t297;
1695 t1665 = t39 * t1663 * t46;
1696 t1669 = 0.9163e4 / 0.59049e5 * t28 * t29 * t1585 * t41 * t46 - 0.21329e5 / 0.19683e5 * t28 * t123 * t1592 * t132 + 0.143528e6 / 0.59049e5 * t28 * t311 * t1597 * t320 - 0.43904e5 / 0.19683e5 * t28 * t780 * t1604 * t789 + 0.14336e5 / 0.19683e5 * t28 * t5 * t779 * sigma[0] / t31 / t1611 / rho[0] * t1620 - 0.46004e5 / 0.59049e5 * t114 * t1625 * t306 - 0.1232e4 / 0.6561e4 * t114 * t1625 * t768 + 0.18368e5 / 0.19683e5 * t114 * t1633 * t775 - 0.896e3 / 0.59049e5 * t114 * t1625 * t763 * t305 + 0.1792e4 / 0.19683e5 * t114 * t1633 * t298 * t774 - 0.7168e4 / 0.19683e5 * t114 * t779 * t1604 * params->a1 * t117 * t1647 + 0.25375e5 / 0.118098e6 * t114 * t1653 * t118 + 0.2135e4 / 0.19683e5 * t114 * t1653 * t299 + 0.1400e4 / 0.59049e5 * t114 * t1653 * t764 + 0.112e3 / 0.59049e5 * t114 * t1653 * t1665;
1697 t1676 = 0.10e2 / 0.27e2 * t7 * t15 * t871 * t51;
1698 t1680 = 0.10e2 / 0.27e2 * t7 * t57 * t871 * t81;
1699 t1686 = t7 * t649 * t681;
1700 t1688 = 0.4e1 / 0.27e2 * t1536 + 0.2e1 / 0.9e1 * t1540 - t7 * t330 * t81 * t1542 / 0.2e1 - 0.20e2 / 0.27e2 * t1548 + 0.2e1 / 0.3e1 * t822 * t376 * t331 * t340 - 0.2e1 / 0.3e1 * t1555 - 0.2e1 / 0.3e1 * t375 * t376 * t141 * t715 + 0.2e1 / 0.3e1 * t807 * t353 * t249 * t264 - 0.2e1 / 0.3e1 * t1566 - 0.2e1 * t352 * t831 * t653 - 0.2e1 / 0.3e1 * t352 * t353 * t92 * t703 - 0.5e1 / 0.9e1 * t1576 + t1579 / 0.2e1 - t1582 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * t1669 + t1676 + t1680 + 0.4e1 / 0.9e1 * t7 * t643 * t136 * t644 - 0.2e1 / 0.3e1 * t1686;
1701 t1690 = t7 * t657 * t258;
1702 t1693 = t7 * t657 * t265;
1703 t1696 = t7 * t677 * t341;
1704 t1702 = t7 * t254 * t691;
1705 t1705 = t7 * t254 * t695;
1706 t1708 = t7 * t254 * t704;
1707 t1722 = t699 * t9;
1708 t1724 = t8 / t1722;
1709 t1726 = -0.12e2 * t700 + 0.12e2 * t1724;
1710 t1732 = t7 * t336 * t716;
1711 t1739 = t13 * t13;
1712 t1741 = 0.1e1 / t246 / t1739;
1713 t1742 = t1741 * t16;
1714 t1743 = t249 * t249;
1715 t1748 = t642 * t98;
1716 t1750 = t7 * t1748 * t645;
1717 t1752 = t247 * t270;
1718 t1754 = t7 * t1752 * t250;
1719 t1756 = t264 * t264;
1720 t1762 = t7 * t367 * t93;
1721 t1764 = t55 * t55;
1722 t1766 = 0.1e1 / t328 / t1764;
1723 t1767 = t1766 * t16;
1724 t1768 = t331 * t331;
1725 t1773 = 0.4e1 / 0.3e1 * t1690 + 0.2e1 / 0.3e1 * t1693 + 0.2e1 / 0.3e1 * t1696 - t7 * t248 * t324 * t249 - 0.2e1 * t1702 - 0.2e1 * t1705 - 0.2e1 / 0.3e1 * t1708 - 0.2e1 * t7 * t87 * t793 * t92 - 0.3e1 * t7 * t87 * t324 * t264 - 0.2e1 * t7 * t87 * t136 * t703 - t7 * t87 * t51 * t1726 / 0.2e1 - 0.2e1 / 0.3e1 * t1732 + t7 * t140 * t81 * t1726 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1742 * t51 * t1743 + 0.4e1 / 0.27e2 * t1750 + 0.2e1 / 0.9e1 * t1754 - t7 * t248 * t51 * t1756 / 0.2e1 - 0.20e2 / 0.27e2 * t1762 - 0.5e1 / 0.27e2 * t7 * t1767 * t81 * t1768;
1726 t1776 = -0.6e1 * t693 - 0.6e1 * t697 - 0.2e1 * t706 - t1528 - 0.2e1 * t713 - 0.2e1 * t718 - t1531 + t725 - 0.3e1 / 0.2e1 * t728 - 0.3e1 / 0.2e1 * t795 + t9 * (t1688 + t1773);
1727 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1728 v4rho4[0] = t1524 + t1776;
1729
1730 t1786 = t7 * t367 * t154;
1731 t1789 = 0.1e1 / t97 / t699;
1732 t1792 = t7 * t14 * t1789 * t368;
1733 t1794 = t329 * t366;
1734 t1799 = t7 * t384 * t159;
1735 t1806 = t7 * t669 * t882;
1736 t1809 = t7 * t677 * t390;
1737 t1810 = t1809 / 0.3e1;
1738 t1812 = t7 * t657 * t361;
1739 t1813 = t1812 / 0.3e1;
1740 t1815 = t7 * t872 * t876;
1741 t1823 = t7 * t56 * t1789 * t385;
1742 t1826 = t7 * t336 * t900;
1743 t1829 = t7 * t254 * t892;
1744 t1835 = t807 * t353 * t835 * t92 / 0.3e1 + t822 * t376 * t839 * t141 / 0.3e1 - 0.5e1 / 0.27e2 * t1786 - 0.25e2 / 0.6e1 * t1792 + t7 * t1794 * t142 / 0.2e1 - 0.5e1 / 0.27e2 * t1799 + t7 * t663 * t192 * t664 / 0.9e1 - t1806 / 0.6e1 + t1810 + t1813 + 0.7e1 / 0.2e1 * t1815 - 0.3e1 / 0.2e1 * t7 * t367 * t324 * t8 + 0.25e2 / 0.6e1 * t1823 - t1826 / 0.2e1 - t1829 / 0.2e1 - t7 * t87 * t793 * t153 / 0.2e1;
1745 t1836 = t247 * t366;
1746 t1863 = t352 * t270 * t51 * t354;
1747 t1864 = t1863 / 0.9e1;
1748 t1865 = t871 * t51;
1749 t1867 = t352 * t1865 * t817;
1750 t1879 = t375 * t270 * t81 * t377;
1751 t1880 = t1879 / 0.9e1;
1752 t1886 = t352 * t98 * t136 * t354;
1753 t1888 = -t7 * t1836 * t93 / 0.2e1 - t7 * t140 * t192 * t715 / 0.2e1 + t1536 / 0.27e2 + t1540 / 0.9e1 - 0.5e1 / 0.9e1 * t1548 - t375 * t928 * t673 / 0.2e1 - t352 * t16 * t324 * t354 / 0.2e1 - t352 * t831 * t835 / 0.2e1 - t352 * t353 * t153 * t703 / 0.6e1 + t1864 + 0.7e1 / 0.6e1 * t1867 - t375 * t376 * t158 * t715 / 0.6e1 + t375 * t843 * t8 * t340 / 0.2e1 + t1880 + t807 * t831 * t808 / 0.3e1 - t1886 / 0.3e1;
1754 t1890 = t366 * t136;
1755 t1894 = t352 * t812 * t835;
1756 t1901 = t375 * t827 * t839;
1757 t1907 = t871 * t81;
1758 t1909 = t375 * t1907 * t844;
1759 t1912 = t4 * t6 * t1741;
1760 t1918 = t807 * t812 * t808;
1761 t1925 = t4 * t6 * t1766;
1762 t1931 = t822 * t827 * t823;
1763 t1938 = -t352 * t1890 * t817 - t1894 / 0.6e1 - t352 * t816 * t8 * t264 / 0.2e1 - t1901 / 0.6e1 - t822 * t843 * t8 * t331 / 0.3e1 - 0.7e1 / 0.6e1 * t1909 - 0.5e1 / 0.27e2 * t1912 * t353 * t153 * t644 + t1918 / 0.9e1 + t807 * t816 * t8 * t249 / 0.3e1 - 0.5e1 / 0.27e2 * t1925 * t376 * t158 * t664 + t1931 / 0.9e1 - t1555 / 0.6e1 - t1566 / 0.6e1 - 0.5e1 / 0.12e2 * t1576 + t1579 / 0.4e1 - t1582 / 0.8e1;
1764 t1940 = t7 * t872 * t51;
1765 t1946 = t7 * t896 * t81;
1766 t1949 = t7 * t720 * t192;
1767 t1962 = t1676 + t1680 + 0.5e1 / 0.2e1 * t1940 - 0.3e1 / 0.2e1 * t7 * t367 * t136 - 0.5e1 / 0.2e1 * t1946 - 0.5e1 / 0.36e2 * t1949 - t1686 / 0.6e1 + 0.2e1 / 0.3e1 * t1690 + t1693 / 0.3e1 + t1696 / 0.3e1 - t1702 / 0.2e1 - t1705 / 0.2e1 - t1708 / 0.6e1 - t1732 / 0.6e1 + t1750 / 0.27e2 + t1754 / 0.9e1 - 0.5e1 / 0.9e1 * t1762;
1768 t1967 = t814 / 0.3e1;
1769 t1969 = t829 / 0.3e1;
1770 t1972 = t9 * (t1835 + t1888 + t1938 + t1962) + t810 / 0.3e1 - t1967 - t819 + t825 / 0.3e1 - t1969 - t833 - t837 / 0.2e1 - t841 / 0.2e1 + t846 - t656;
1771 t1978 = -t676 - t1528 - t1531 + 0.3e1 / 0.4e1 * t725 - 0.3e1 / 0.4e1 * t728 - t796 - 0.3e1 / 0.2e1 * t851 + 0.3e1 / 0.2e1 * t854 + t1029 + t648 - t651 / 0.3e1;
1772 t1982 = t659 + t668 - t671 / 0.3e1 + t679 - t684 - 0.2e1 * t686 - t689 - t694 - t698 - t707 - t713;
1773 t1989 = -t719 + t1004 + 0.7e1 / 0.2e1 * t874 - 0.3e1 * t878 + t1005 - t884 / 0.2e1 - t887 - t890 - 0.3e1 / 0.2e1 * t894 - 0.7e1 / 0.2e1 * t898 - 0.3e1 / 0.2e1 * t902;
1774 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1775 v4rho4[1] = t1972 + t1978 + t1982 + t1989;
1776
1777 t2000 = t7 * t1752 * t401;
1778 t2003 = t7 * t657 * t407;
1779 t2006 = t7 * t1538 * t412;
1780 t2012 = t7 * t677 * t417;
1781 t2015 = t7 * t677 * t421;
1782 t2022 = t7 * t336 * t989;
1783 t2028 = 0.1e1 / t97 / t1722;
1784 t2030 = t8 * t8;
1785 t2036 = t7 * t649 * t967;
1786 t2039 = t7 * t254 * t971;
1787 t2042 = t7 * t254 * t977;
1788 t2050 = t7 * t896 * t981;
1789 t2053 = t7 * t336 * t985;
1790 t2063 = -0.4e1 / 0.9e1 * t822 * t843 * t377 * t8 + 0.4e1 / 0.9e1 * t807 * t816 * t354 * t8 + t2000 / 0.27e2 + t2003 / 0.9e1 + t2006 / 0.27e2 + t7 * t1794 * t159 / 0.3e1 + 0.2e1 / 0.9e1 * t2012 + t2015 / 0.9e1 - t7 * t330 * t470 * t331 / 0.6e1 - t2022 / 0.3e1 - t7 * t1836 * t154 / 0.3e1 - t7 * t247 * t2028 * t51 * t2030 / 0.3e1 - t2036 / 0.9e1 - t2039 / 0.3e1 - t2042 / 0.3e1 - t7 * t329 * t2028 * t81 * t2030 / 0.3e1 - 0.7e1 / 0.3e1 * t2050 - t2053 / 0.3e1 - t7 * t248 * t324 * t400 / 0.6e1 - t7 * t87 * t324 * t406 / 0.2e1;
1791 t2105 = -t7 * t87 * t136 * t976 - t7 * t140 * t470 * t340 / 0.2e1 - 0.10e2 / 0.27e2 * t1786 - 0.41e2 / 0.9e1 * t1792 - 0.10e2 / 0.27e2 * t1799 - t1806 / 0.9e1 + 0.4e1 / 0.9e1 * t1809 + 0.4e1 / 0.9e1 * t1812 - 0.2e1 / 0.3e1 * t1815 + 0.41e2 / 0.9e1 * t1823 - t1826 / 0.3e1 - t1829 / 0.3e1 + t1540 / 0.27e2 - 0.10e2 / 0.27e2 * t1548 + 0.2e1 / 0.9e1 * t807 * t831 * t912 + t807 * t353 * t400 * t264 / 0.9e1 - t352 * t831 * t916 / 0.3e1 - t352 * t353 * t976 * t92 / 0.3e1 - t352 * t353 * t406 * t264 / 0.6e1 + t822 * t376 * t411 * t340 / 0.9e1;
1792 t2108 = t375 * t1907 * t924;
1793 t2113 = t366 * t192;
1794 t2126 = t822 * t827 * t920;
1795 t2133 = t375 * t98 * t192 * t377;
1796 t2140 = t375 * t827 * t932;
1797 t2150 = t807 * t812 * t912;
1798 t2153 = t352 * t1865 * t936;
1799 t2160 = t352 * t812 * t916;
1800 t2169 = -0.7e1 / 0.9e1 * t2108 - t375 * t928 * t839 / 0.3e1 + 0.2e1 / 0.3e1 * t375 * t2113 * t844 - t375 * t376 * t984 * t141 / 0.3e1 - t375 * t376 * t420 * t340 / 0.6e1 + 0.2e1 / 0.27e2 * t2126 + 0.2e1 / 0.9e1 * t822 * t928 * t823 - 0.2e1 / 0.9e1 * t2133 + t822 * t376 * t420 * t331 / 0.9e1 - t2140 / 0.9e1 - 0.2e1 / 0.3e1 * t352 * t1890 * t936 - 0.5e1 / 0.27e2 * t1912 * t353 * t400 * t249 + 0.2e1 / 0.27e2 * t2150 + 0.7e1 / 0.9e1 * t2153 + t807 * t353 * t406 * t249 / 0.9e1 - t2160 / 0.9e1 - 0.5e1 / 0.27e2 * t1925 * t376 * t411 * t331 + 0.4e1 / 0.27e2 * t1863 - 0.2e1 / 0.9e1 * t1867 + 0.4e1 / 0.27e2 * t1879;
1801 t2184 = t7 * t345 * t470;
1802 t2191 = -0.2e1 / 0.9e1 * t1886 - t1894 / 0.9e1 - t1901 / 0.9e1 + 0.2e1 / 0.9e1 * t1909 + 0.2e1 / 0.27e2 * t1918 + 0.2e1 / 0.27e2 * t1931 - 0.5e1 / 0.18e2 * t1576 + t1579 / 0.12e2 + t1676 + t1680 - t1940 / 0.3e1 + t1946 / 0.3e1 - 0.5e1 / 0.18e2 * t1949 + t7 * t384 * t192 + t2184 / 0.12e2 + 0.2e1 / 0.9e1 * t1690 + t1693 / 0.9e1 + t1696 / 0.9e1 + t1754 / 0.27e2 - 0.10e2 / 0.27e2 * t1762;
1803 t2207 = t9 * (t2063 + t2105 + t2169 + t2191) + 0.2e1 / 0.9e1 * t914 - t918 / 0.3e1 + 0.2e1 / 0.9e1 * t922 + 0.2e1 / 0.3e1 * t926 - 0.2e1 / 0.3e1 * t930 - t934 / 0.3e1 - 0.2e1 / 0.3e1 * t938 + 0.2e1 / 0.9e1 * t810 - 0.4e1 / 0.9e1 * t814 - 0.2e1 / 0.3e1 * t819 + 0.2e1 / 0.9e1 * t825 - 0.4e1 / 0.9e1 * t829;
1804 t2217 = -0.2e1 / 0.3e1 * t833 - t837 / 0.3e1 - t841 / 0.3e1 + 0.2e1 / 0.3e1 * t846 - t1528 - t1531 + t725 / 0.2e1 - t728 / 0.4e1 - t851 + t854 + t857 / 0.2e1 - t943 / 0.4e1 - t651 / 0.9e1;
1805 t2232 = 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 * t868 + 0.5e1 / 0.3e1 * t874 - 0.2e1 * t878 + 0.2e1 / 0.3e1 * t880 - t884 / 0.3e1 - 0.4e1 / 0.3e1 * t887 - 0.4e1 / 0.3e1 * t890;
1806 t2241 = -t894 - 0.5e1 / 0.3e1 * t898 - t902 - t953 / 0.9e1 - t956 / 0.3e1 - t959 / 0.9e1 - 0.2e1 / 0.3e1 * t962 - t965 / 0.3e1 - t969 / 0.3e1 - t973 - t979 + 0.2e1 * t983 - t987 - t991;
1807 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1808 v4rho4[2] = t2207 + t2217 + t2232 + t2241;
1809
1810 t2253 = -0.6e1 * t700 - 0.12e2 * t1724;
1811 t2263 = t7 * t1748 * t1011;
1812 t2266 = t7 * t1534 * t1020;
1813 t2269 = t7 * t254 * t1033;
1814 t2272 = t7 * t669 * t1037;
1815 t2275 = t7 * t336 * t1041;
1816 t2278 = t7 * t336 * t1045;
1817 t2281 = t7 * t336 * t1050;
1818 t2304 = t807 * t353 * t1015 * t92 / 0.3e1 + t822 * t376 * t1025 * t141 / 0.3e1 - t7 * t140 * t81 * t2253 / 0.2e1 - t7 * t140 * t1115 * t141 / 0.2e1 + t2263 / 0.27e2 + t2266 / 0.27e2 - t2269 / 0.6e1 - t2272 / 0.6e1 - t2275 / 0.2e1 - t2278 / 0.2e1 - t2281 / 0.6e1 + t7 * t643 * t136 * t1010 / 0.9e1 - t7 * t87 * t136 * t1032 / 0.2e1 + t7 * t87 * t51 * t2253 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t384 * t470 * t8 - 0.3e1 / 0.2e1 * t7 * t140 * t192 * t984;
1819 t2319 = t2000 / 0.9e1 + t2003 / 0.3e1 + t2006 / 0.9e1 + 0.2e1 / 0.3e1 * t2012 + t2015 / 0.3e1 - t2022 / 0.2e1 - t2036 / 0.6e1 - t2039 / 0.2e1 - t2042 / 0.2e1 + t2050 - t2053 / 0.2e1 - 0.5e1 / 0.9e1 * t1786 + t1792 / 0.3e1 - 0.5e1 / 0.9e1 * t1799 + t1810 + t1813 - t1823 / 0.3e1;
1820 t2374 = t352 * t812 * t1015;
1821 t2380 = -0.5e1 / 0.27e2 * t1548 - 0.5e1 / 0.27e2 * t1912 * t353 * t1010 * t92 - t352 * t831 * t1015 / 0.2e1 - t352 * t353 * t153 * t976 / 0.2e1 - 0.5e1 / 0.27e2 * t1925 * t376 * t1019 * t141 - t822 * t843 * t411 * t8 / 0.3e1 - t375 * t376 * t158 * t984 / 0.2e1 + t375 * t843 * t8 * t420 / 0.2e1 - t352 * t353 * t1032 * t92 / 0.6e1 + t822 * t928 * t920 / 0.3e1 + t375 * t2113 * t924 - t375 * t16 * t470 * t377 / 0.2e1 - t375 * t928 * t932 / 0.2e1 - t375 * t376 * t1049 * t141 / 0.6e1 + t807 * t816 * t400 * t8 / 0.3e1 - t2374 / 0.6e1 - t352 * t816 * t8 * t406 / 0.2e1;
1822 t2382 = t375 * t827 * t1025;
1823 t2395 = t7 * t146 * t1115;
1824 t2398 = -t2382 / 0.6e1 + t2108 / 0.3e1 + t2126 / 0.9e1 - t2133 / 0.3e1 - t2140 / 0.6e1 + t2150 / 0.9e1 - t2153 / 0.3e1 - t2160 / 0.6e1 + t1864 + t1880 - 0.5e1 / 0.36e2 * t1576 + t1676 + t1680 - 0.5e1 / 0.12e2 * t1949 + t2184 / 0.4e1 - t2395 / 0.8e1 - 0.5e1 / 0.27e2 * t1762;
1825 t2406 = t9 * (t2304 + t2319 + t2380 + t2398) - t1018 - t1028 + t914 / 0.3e1 - t918 / 0.2e1 + t922 / 0.3e1 + t926 - t930 - t934 / 0.2e1 - t938 - t1967;
1826 t2409 = -t1969 - t1528 - t1531 + t726 + 0.3e1 / 0.4e1 * t857 - 0.3e1 / 0.4e1 * t943 - t1118 + t660 + t680 + t868 - t874;
1827 t2416 = t880 - t887 - t890 + t898 - t953 / 0.3e1 - t956 - t959 / 0.3e1 - 0.2e1 * t962 - t965 - t969 / 0.2e1 - 0.3e1 / 0.2e1 * t973;
1828 t2421 = -0.3e1 / 0.2e1 * t979 + 0.3e1 * t983 - 0.3e1 / 0.2e1 * t987 - 0.3e1 / 0.2e1 * t991 + t1014 + t1023 - t1036 - t1040 - t1044 - t1048 - t1053;
1829 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1830 v4rho4[3] = t2406 + t2409 + t2416 + t2421;
1831
1832 t2434 = 0.4e1 / 0.3e1 * t868 + 0.4e1 / 0.3e1 * t880 - 0.2e1 / 0.3e1 * t953 - 0.2e1 * t956 - 0.2e1 / 0.3e1 * t959 - 0.4e1 * t962 - 0.2e1 * t965 + 0.4e1 / 0.9e1 * t1013 - 0.2e1 * t1017 + 0.4e1 / 0.9e1 * t1022 - 0.2e1 * t1027;
1833 t2442 = t406 * t406;
1834 t2447 = t411 * t411;
1835 t2452 = t420 * t420;
1836 t2462 = 0.12e2 * t700 + 0.12e2 * t1724;
1837 t2487 = t400 * t400;
1838 t2500 = -t7 * t248 * t51 * t2442 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1767 * t81 * t2447 - t7 * t330 * t81 * t2452 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t663 * t192 * t1019 - t7 * t87 * t51 * t2462 / 0.2e1 - t7 * t330 * t470 * t411 - 0.2e1 * t7 * t140 * t1115 * t158 - 0.3e1 * t7 * t140 * t470 * t420 - 0.2e1 * t7 * t140 * t192 * t1049 + t7 * t140 * t81 * t2462 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1742 * t51 * t2487 + 0.4e1 / 0.27e2 * t2263 + 0.4e1 / 0.27e2 * t2266 - 0.2e1 / 0.3e1 * t2269 - 0.2e1 / 0.3e1 * t2272 - 0.2e1 * t2275 - 0.2e1 * t2278 - 0.2e1 / 0.3e1 * t2281 + 0.2e1 / 0.9e1 * t2000;
1839 t2532 = 0.1e1 / t62 / t181;
1840 t2539 = 0.1e1 / t61 / t558;
1841 t2544 = 0.1e1 / t1312;
1842 t2549 = t459 * t181;
1843 t2551 = 0.1e1 / t62 / t2549;
1844 t2558 = t459 * t459;
1845 t2566 = t71 / t1108 / t75 * t1619;
1846 t2571 = t178 * t2539 * params->a1;
1847 t2579 = t457 * t2544 * params->a1;
1848 t2593 = t1109 * t788;
1849 t2599 = sigma[2] * t2532 * params->a1;
1850 t2610 = t69 * t1663 * t76;
1851 t2614 = 0.9163e4 / 0.59049e5 * t28 * t59 * t2532 * t71 * t76 - 0.21329e5 / 0.19683e5 * t28 * t179 * t2539 * t188 + 0.143528e6 / 0.59049e5 * t28 * t458 * t2544 * t466 - 0.43904e5 / 0.19683e5 * t28 * t1103 * t2551 * t1111 + 0.14336e5 / 0.19683e5 * t28 * t5 * t1102 * sigma[2] / t61 / t2558 / rho[1] * t2566 - 0.46004e5 / 0.59049e5 * t114 * t2571 * t453 - 0.1232e4 / 0.6561e4 * t114 * t2571 * t1091 + 0.18368e5 / 0.19683e5 * t114 * t2579 * t1098 - 0.896e3 / 0.59049e5 * t114 * t2571 * t1086 * t452 + 0.1792e4 / 0.19683e5 * t114 * t2579 * t445 * t1097 - 0.7168e4 / 0.19683e5 * t114 * t1102 * t2551 * params->a1 * t173 * t2593 + 0.25375e5 / 0.118098e6 * t114 * t2599 * t174 + 0.2135e4 / 0.19683e5 * t114 * t2599 * t446 + 0.1400e4 / 0.59049e5 * t114 * t2599 * t1087 + 0.112e3 / 0.59049e5 * t114 * t2599 * t2610;
1852 t2618 = 0.2e1 / 0.3e1 * t2003 + 0.2e1 / 0.9e1 * t2006 + 0.4e1 / 0.3e1 * t2012 + 0.2e1 / 0.3e1 * t2015 - 0.20e2 / 0.27e2 * t1786 - 0.20e2 / 0.27e2 * t1799 + 0.2e1 / 0.3e1 * t807 * t353 * t400 * t406 - 0.2e1 / 0.3e1 * t352 * t353 * t153 * t1032 + 0.2e1 / 0.3e1 * t822 * t376 * t411 * t420 - 0.2e1 * t375 * t928 * t1025 - 0.2e1 / 0.3e1 * t375 * t376 * t158 * t1049 - 0.2e1 / 0.3e1 * t2374 - 0.2e1 / 0.3e1 * t2382 + t1676 + t1680 - 0.5e1 / 0.9e1 * t1949 + t2184 / 0.2e1 - t2395 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t58 * t2614;
1853 t2621 = -t1528 - t1531 + t857 - 0.3e1 / 0.2e1 * t943 - 0.2e1 * t1035 - 0.2e1 * t1039 - 0.6e1 * t1043 - 0.6e1 * t1047 - 0.2e1 * t1052 - 0.3e1 / 0.2e1 * t1117 + t9 * (t2500 + t2618);
1854 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1855 v4rho4[4] = t2434 + t2621;
1856
1857 t2626 = t7 * t99 * t514;
1858 t2629 = t7 * t17 * t1191;
1859 t2632 = t757 * params->a1 * t39;
1860 t2637 = t5 * t731 * params->a1;
1861 t2676 = t743 * params->a1 * t39;
1862 t2677 = t762 * t130;
1863 t2682 = t297 * t317;
1864 t2689 = params->alpha * t786;
1865 t2701 = -0.2044e4 / 0.6561e4 * t114 * t2632 * t1180 - 0.889e3 / 0.19683e5 * t28 * t2637 * t118 + 0.1771e4 / 0.6561e4 * t28 * t5 * t743 * t41 * t494 - 0.721e3 / 0.26244e5 * t28 * t2637 * t299 - 0.14350e5 / 0.19683e5 * t28 * t5 * t757 * t41 * t1164 - 0.287e3 / 0.39366e5 * t28 * t2637 * t764 + 0.4984e4 / 0.6561e4 * t28 * t5 * t783 * t41 * t1647 * t310 - 0.14e2 / 0.19683e5 * t28 * t2637 * t1665 - 0.1792e4 / 0.6561e4 * t28 * t780 / t31 / t1611 * t1620 - 0.539e3 / 0.19683e5 * t202 * t203 * t733 + 0.112e3 / 0.19683e5 * t114 * t2676 * t2677 * t43 - 0.224e3 / 0.6561e4 * t114 * t2632 * t2682 * t1179 + 0.896e3 / 0.6561e4 * t114 * t783 * params->a1 * t39 * t2689 * t788 * t310 + 0.4333e4 / 0.19683e5 * t114 * t2676 * t504 + 0.133e3 / 0.2187e4 * t114 * t2676 * t1172;
1866 t2706 = t7 * t87 * t1134;
1867 t2709 = t7 * t87 * t1137;
1868 t2713 = t7 * t271 * t219 / 0.9e1;
1869 t2731 = t7 * t248 * t1130;
1870 t2734 = t7 * t254 * t480;
1871 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1872 v4rho3sigma[0] = -t1216 * t1217 * t653 / 0.2e1 - t2626 / 0.2e1 - 0.3e1 / 0.2e1 * t2629 - 0.3e1 / 0.8e1 * t200 * t201 * t2701 - 0.4e1 * t2706 - 0.2e1 * t2709 + t2713 - t200 * t351 * t514 * t249 / 0.2e1 - 0.3e1 / 0.2e1 * t200 * t479 * t1191 * t92 - 0.3e1 / 0.2e1 * t200 * t479 * t514 * t264 - t200 * t479 * t219 * t703 / 0.2e1 - 0.2e1 / 0.3e1 * t2731 - 0.2e1 / 0.3e1 * t2734 + t200 * t806 * t219 * t644 / 0.9e1;
1873
1874 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1875 v4rho3sigma[1] = 0.0e0;
1876
1877 t2741 = t7 * t330 * t1201;
1878 t2744 = t7 * t336 * t521;
1879 t2747 = t7 * t140 * t1205;
1880 t2751 = t7 * t345 * t239 / 0.9e1;
1881 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1882 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2741 - 0.2e1 / 0.3e1 * t2744 - 0.2e1 * t2747 + t2751 + t200 * t821 * t239 * t664 / 0.9e1 - t1216 * t1234 * t673 / 0.2e1 - t200 * t520 * t239 * t715 / 0.2e1;
1883
1884 t2764 = t98 * t1 * t3;
1885 t2765 = t2764 * t526;
1886 t2767 = t1211 * t1215;
1887 t2769 = 0.4e1 / 0.9e1 * t2767 * t1218;
1888 t2771 = t366 * t1 * t3;
1889 t2772 = t2771 * t1228;
1890 t2786 = t642 * t219;
1891 t2790 = t1225 * t1215;
1892 t2797 = 0.4e1 / 0.3e1 * t1212 * t1222;
1893 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1894 v4rho3sigma[3] = -0.2e1 / 0.9e1 * t2765 - t2769 + t2772 / 0.6e1 - t1216 * t247 * t514 * t354 / 0.3e1 - t1216 * t1217 * t835 / 0.6e1 - t1226 * t479 * t514 * t8 - t2626 / 0.3e1 - t2629 / 0.2e1 + t1216 * t2786 * t808 / 0.9e1 - t2790 * t1217 * t817 / 0.3e1 - 0.4e1 / 0.3e1 * t2706 - 0.2e1 / 0.3e1 * t2709 + t2713 - t2797 - t200 * t479 * t1191 * t153 / 0.2e1 - t1226 * t479 * t219 / 0.2e1 - 0.2e1 / 0.9e1 * t2731 - 0.4e1 / 0.9e1 * t2734;
1895
1896 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1897 v4rho3sigma[4] = 0.0e0;
1898
1899 t2810 = t2764 * t531;
1900 t2813 = 0.4e1 / 0.9e1 * t2767 * t1235;
1901 t2814 = t2771 * t1239;
1902 t2816 = t662 * t239;
1903 t2829 = t2764 * t565;
1904 t2832 = 0.4e1 / 0.3e1 * t1212 * t1245;
1905 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1906 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2741 - 0.4e1 / 0.9e1 * t2744 - 0.2e1 / 0.3e1 * t2747 + t2751 - 0.2e1 / 0.9e1 * t2810 - t2813 - t2814 / 0.6e1 + t1216 * t2816 * t823 / 0.9e1 + t2790 * t1234 * t844 / 0.3e1 - t1216 * t1234 * t839 / 0.6e1 + t1226 * t520 * t239 / 0.2e1 - t2829 / 0.6e1 - t2832 - t200 * t374 * t564 * t331 / 0.6e1 - t200 * t520 * t564 * t340 / 0.2e1;
1907
1908 t2845 = t1212 * t1250;
1909 t2857 = t1212 * t1254;
1910 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1911 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2765 - t2769 - t2797 - 0.4e1 / 0.3e1 * t2772 - 0.2e1 / 0.9e1 * t2734 + t2713 - t2626 / 0.6e1 - 0.2e1 / 0.9e1 * t2845 + t1216 * t2786 * t912 / 0.9e1 - t200 * t351 * t514 * t400 / 0.6e1 - t2790 * t1217 * t936 / 0.3e1 - 0.2e1 / 0.3e1 * t2857 - t1216 * t1217 * t916 / 0.6e1 - t200 * t479 * t514 * t406 / 0.2e1 - t200 * t479 * t219 * t976 / 0.2e1;
1912
1913 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1914 v4rho3sigma[7] = 0.0e0;
1915
1916 t2874 = t1212 * t1259;
1917 t2882 = t1212 * t1263;
1918 t2891 = t1212 * t1266;
1919 t2900 = t1212 * t1320;
1920 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1921 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2810 - t2813 + 0.4e1 / 0.3e1 * t2814 - 0.2e1 / 0.9e1 * t2744 + t2751 - t2829 / 0.3e1 - t2832 - 0.2e1 / 0.9e1 * t2874 + t1216 * t2816 * t920 / 0.9e1 + t2790 * t1234 * t924 / 0.3e1 - 0.4e1 / 0.3e1 * t2882 - t1216 * t329 * t564 * t377 / 0.3e1 + t1226 * t520 * t564 * t8 - 0.2e1 / 0.3e1 * t2891 - t1216 * t1234 * t932 / 0.6e1 - t200 * t520 * t239 * t984 / 0.2e1 - t2900 / 0.2e1 - t200 * t520 * t1319 * t141 / 0.2e1;
1922
1923 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1924 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2765 - 0.2e1 / 0.3e1 * t2845 - 0.2e1 * t2857 + t2713 + t200 * t806 * t219 * t1010 / 0.9e1 - t1216 * t1217 * t1015 / 0.2e1 - t200 * t479 * t219 * t1032 / 0.2e1;
1925
1926 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1927 v4rho3sigma[10] = 0.0e0;
1928
1929 t2932 = t5 * t1055 * params->a1;
1930 t2971 = t1067 * params->a1 * t69;
1931 t2979 = t1081 * params->a1 * t69;
1932 t2983 = t762 * t186;
1933 t2988 = t297 * t464;
1934 t2995 = params->alpha * t1109;
1935 t3001 = -0.889e3 / 0.19683e5 * t28 * t2932 * t174 + 0.1771e4 / 0.6561e4 * t28 * t5 * t1067 * t71 * t544 - 0.721e3 / 0.26244e5 * t28 * t2932 * t446 - 0.14350e5 / 0.19683e5 * t28 * t5 * t1081 * t71 * t1292 - 0.287e3 / 0.39366e5 * t28 * t2932 * t1087 + 0.4984e4 / 0.6561e4 * t28 * t5 * t1106 * t71 * t2593 * t457 - 0.14e2 / 0.19683e5 * t28 * t2932 * t2610 - 0.1792e4 / 0.6561e4 * t28 * t1103 / t61 / t2558 * t2566 - 0.539e3 / 0.19683e5 * t202 * t203 * t1057 + 0.4333e4 / 0.19683e5 * t114 * t2971 * t554 + 0.133e3 / 0.2187e4 * t114 * t2971 * t1300 - 0.2044e4 / 0.6561e4 * t114 * t2979 * t1308 + 0.112e3 / 0.19683e5 * t114 * t2971 * t2983 * t73 - 0.224e3 / 0.6561e4 * t114 * t2979 * t2988 * t1307 + 0.896e3 / 0.6561e4 * t114 * t1106 * params->a1 * t69 * t2995 * t788 * t457;
1936 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1937 v4rho3sigma[11] = -0.2e1 / 0.3e1 * t2810 - 0.2e1 / 0.3e1 * t2874 - 0.4e1 * t2882 - 0.2e1 * t2891 + t200 * t821 * t239 * t1019 / 0.9e1 - t1216 * t1234 * t1025 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t223 * t3001 + t2751 - t2829 / 0.2e1 - 0.3e1 / 0.2e1 * t2900 - t200 * t374 * t564 * t411 / 0.2e1 - 0.3e1 / 0.2e1 * t200 * t520 * t1319 * t158 - 0.3e1 / 0.2e1 * t200 * t520 * t564 * t420 - t200 * t520 * t239 * t1049 / 0.2e1;
1938
1939 t3024 = t2764 * t597 / 0.6e1;
1940 t3025 = t1212 * t1326;
1941 t3027 = t1212 * t1372;
1942 t3059 = t279 * params->a1 * t39;
1943 t3098 = 0.224e3 / 0.2187e4 * t28 * t311 / t31 / t312 / t290 * t1620 - 0.266e3 / 0.6561e4 * t28 * t1149 * t132 + 0.1169e4 / 0.6561e4 * t28 * t1163 * t1345 - 0.175e3 / 0.729e3 * t28 * t5 * t1186 * t41 * t1647 * t122 + 0.77e2 / 0.26244e5 * t114 * t3059 * t572 + 0.413e3 / 0.4374e4 * t114 * t1177 * t1361 + 0.28e2 / 0.2187e4 * t114 * t1177 * t2682 * t1360 - 0.112e3 / 0.2187e4 * t114 * t1186 * params->a1 * t39 * t2689 * t788 * t122 + 0.35e2 / 0.8748e4 * t114 * t3059 * t581 - 0.1211e4 / 0.26244e5 * t114 * t1158 * t587 + 0.7e1 / 0.3888e4 * t114 * t3059 * t1350 - 0.301e3 / 0.17496e5 * t114 * t1158 * t1354 + 0.7e1 / 0.26244e5 * t114 * t3059 * t1663 * t570 * t46 - 0.14e2 / 0.6561e4 * t114 * t1158 * t2677 * params->b1;
1944 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1945 v4rho2sigma2[0] = -t3024 - 0.4e1 / 0.3e1 * t3025 - t3027 - t200 * t351 * t596 * t249 / 0.6e1 - t200 * t479 * t1371 * t92 - t200 * t479 * t596 * t264 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t201 * t3098;
1946
1947 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1948 v4rho2sigma2[1] = 0.0e0;
1949
1950 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1951 v4rho2sigma2[2] = 0.0e0;
1952
1953 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1954 v4rho2sigma2[3] = 0.0e0;
1955
1956 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1957 v4rho2sigma2[4] = 0.0e0;
1958
1959 t3103 = t2764 * t629 / 0.6e1;
1960 t3104 = t1212 * t1377;
1961 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1962 v4rho2sigma2[5] = -t3103 - 0.4e1 / 0.3e1 * t3104 - t200 * t374 * t628 * t331 / 0.6e1 - t200 * t520 * t628 * t340 / 0.2e1;
1963
1964 t3116 = t1212 * t1381;
1965 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1966 v4rho2sigma2[6] = -t3024 - 0.2e1 / 0.3e1 * t3025 - t3027 / 0.2e1 - 0.2e1 / 0.3e1 * t3116 - t1216 * t247 * t596 * t354 / 0.6e1 - t200 * t479 * t1371 * t153 / 0.2e1 - t1226 * t479 * t596 * t8 / 0.2e1;
1967
1968 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1969 v4rho2sigma2[7] = 0.0e0;
1970
1971 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1972 v4rho2sigma2[8] = 0.0e0;
1973
1974 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1975 v4rho2sigma2[9] = 0.0e0;
1976
1977 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1978 v4rho2sigma2[10] = 0.0e0;
1979
1980 t3131 = t1212 * t1386;
1981 t3141 = t1212 * t1432;
1982 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1983 v4rho2sigma2[11] = -t3103 - 0.2e1 / 0.3e1 * t3104 - 0.2e1 / 0.3e1 * t3131 - t1216 * t329 * t628 * t377 / 0.6e1 + t1226 * t520 * t628 * t8 / 0.2e1 - t3141 / 0.2e1 - t200 * t520 * t1431 * t141 / 0.2e1;
1984
1985 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1986 v4rho2sigma2[12] = -t3024 - 0.4e1 / 0.3e1 * t3116 - t200 * t351 * t596 * t400 / 0.6e1 - t200 * t479 * t596 * t406 / 0.2e1;
1987
1988 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1989 v4rho2sigma2[13] = 0.0e0;
1990
1991 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1992 v4rho2sigma2[14] = 0.0e0;
1993
1994 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1995 v4rho2sigma2[15] = 0.0e0;
1996
1997 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1998 v4rho2sigma2[16] = 0.0e0;
1999
2000 t3188 = t427 * params->a1 * t69;
2001 t3227 = -0.266e3 / 0.6561e4 * t28 * t1277 * t188 + 0.1169e4 / 0.6561e4 * t28 * t1291 * t1405 - 0.175e3 / 0.729e3 * t28 * t5 * t1314 * t71 * t2593 * t178 + 0.224e3 / 0.2187e4 * t28 * t458 / t61 / t459 / t438 * t2566 + 0.77e2 / 0.26244e5 * t114 * t3188 * t604 + 0.413e3 / 0.4374e4 * t114 * t1305 * t1421 + 0.28e2 / 0.2187e4 * t114 * t1305 * t2988 * t1420 - 0.112e3 / 0.2187e4 * t114 * t1314 * params->a1 * t69 * t2995 * t788 * t178 + 0.35e2 / 0.8748e4 * t114 * t3188 * t613 - 0.1211e4 / 0.26244e5 * t114 * t1286 * t619 + 0.7e1 / 0.3888e4 * t114 * t3188 * t1410 - 0.301e3 / 0.17496e5 * t114 * t1286 * t1414 + 0.7e1 / 0.26244e5 * t114 * t3188 * t1663 * t602 * t76 - 0.14e2 / 0.6561e4 * t114 * t1286 * t2983 * params->b1;
2002 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2003 v4rho2sigma2[17] = -t3103 - 0.4e1 / 0.3e1 * t3131 - t3141 - t200 * t374 * t628 * t411 / 0.6e1 - t200 * t520 * t1431 * t158 - t200 * t520 * t628 * t420 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t223 * t3227;
2004
2005 t3231 = t1212 * t1471;
2006 t3292 = 0.7e1 / 0.8748e4 * t114 * t1330 * t1437 + 0.7e1 / 0.17496e5 * t114 * t1330 * t297 * t1435 * t46 + 0.35e2 / 0.8748e4 * t114 * t502 * t1441 + 0.7e1 / 0.1944e4 * t114 * t502 * t1454 - 0.35e2 / 0.1458e4 * t114 * t1359 * t1460 - 0.7e1 / 0.243e3 * t28 * t1343 * t320 + 0.49e2 / 0.729e3 * t28 * t5 * t1366 * t41 * t1647 * sigma[0] - 0.7e1 / 0.34992e5 * t114 * t1330 * t1450 - 0.7e1 / 0.69984e5 * t114 * t1330 * t1663 * t1435 * t46 + 0.7e1 / 0.8748e4 * t114 * t502 * t1349 * t305 - 0.7e1 / 0.1458e4 * t114 * t1359 * t2682 * t319 + 0.14e2 / 0.729e3 * t114 * t1366 * params->a1 * t39 * t2689 * t788 * sigma[0] - 0.28e2 / 0.729e3 * t28 * t123 / t31 / t1602 * t1620;
2007 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2008 v4rhosigma3[0] = -t3231 / 0.2e1 - t200 * t479 * t1470 * t92 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t201 * t3292;
2009
2010 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2011 v4rhosigma3[1] = 0.0e0;
2012
2013 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2014 v4rhosigma3[2] = 0.0e0;
2015
2016 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2017 v4rhosigma3[3] = 0.0e0;
2018
2019 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2020 v4rhosigma3[4] = 0.0e0;
2021
2022 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2023 v4rhosigma3[5] = 0.0e0;
2024
2025 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2026 v4rhosigma3[6] = 0.0e0;
2027
2028 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2029 v4rhosigma3[7] = 0.0e0;
2030
2031 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2032 v4rhosigma3[8] = 0.0e0;
2033
2034 t3296 = t1212 * t1510;
2035 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2036 v4rhosigma3[9] = -t200 * t520 * t1509 * t141 / 0.2e1 - t3296 / 0.2e1;
2037
2038 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2039 v4rhosigma3[10] = -t200 * t479 * t1470 * t153 / 0.2e1 - t3231 / 0.2e1;
2040
2041 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2042 v4rhosigma3[11] = 0.0e0;
2043
2044 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2045 v4rhosigma3[12] = 0.0e0;
2046
2047 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2048 v4rhosigma3[13] = 0.0e0;
2049
2050 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2051 v4rhosigma3[14] = 0.0e0;
2052
2053 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2054 v4rhosigma3[15] = 0.0e0;
2055
2056 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2057 v4rhosigma3[16] = 0.0e0;
2058
2059 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2060 v4rhosigma3[17] = 0.0e0;
2061
2062 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2063 v4rhosigma3[18] = 0.0e0;
2064
2065 t3365 = 0.7e1 / 0.8748e4 * t114 * t1390 * t1476 + 0.7e1 / 0.17496e5 * t114 * t1390 * t297 * t1474 * t76 + 0.35e2 / 0.8748e4 * t114 * t552 * t1480 + 0.7e1 / 0.1944e4 * t114 * t552 * t1493 - 0.35e2 / 0.1458e4 * t114 * t1419 * t1499 - 0.7e1 / 0.243e3 * t28 * t1403 * t466 + 0.49e2 / 0.729e3 * t28 * t5 * t1426 * t71 * t2593 * sigma[2] - 0.7e1 / 0.34992e5 * t114 * t1390 * t1489 - 0.7e1 / 0.69984e5 * t114 * t1390 * t1663 * t1474 * t76 + 0.7e1 / 0.8748e4 * t114 * t552 * t1409 * t452 - 0.7e1 / 0.1458e4 * t114 * t1419 * t2988 * t319 + 0.14e2 / 0.729e3 * t114 * t1426 * params->a1 * t69 * t2995 * t788 * sigma[2] - 0.28e2 / 0.729e3 * t28 * t179 / t61 / t2549 * t2566;
2066 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2067 v4rhosigma3[19] = -t3296 / 0.2e1 - t200 * t520 * t1509 * t158 / 0.2e1 - 0.3e1 / 0.8e1 * t200 * t223 * t3365;
2068
2069 t3369 = 0.1e1 / t310;
2070 t3422 = -0.7e1 / 0.46656e5 * t114 * t569 * t297 * t3369 * t46 + 0.7e1 / 0.11664e5 * t114 * t569 * params->alpha * t3369 * t46 + 0.7e1 / 0.5832e4 * t114 * t586 * t1436 * t305 + 0.7e1 / 0.1944e4 * t114 * t1459 * t571 * t774 - 0.7e1 / 0.486e3 * t28 * t5 * t1465 * t789 + 0.7e1 / 0.186624e6 * t114 * t569 * t1663 * t3369 * t46 - 0.7e1 / 0.46656e5 * t114 * t569 * t762 * t3369 * t46 - 0.7e1 / 0.23328e5 * t114 * t586 * t1449 * t305 + 0.7e1 / 0.3888e4 * t114 * t1459 * t580 * t774 - 0.7e1 / 0.972e3 * t114 * t1465 * params->a1 * t39 * t2689 * t788 + 0.7e1 / 0.486e3 * t28 * t29 / t31 / t781 * t1620;
2071 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2072 v4sigma4[0] = -0.3e1 / 0.8e1 * t200 * t201 * t3422;
2073
2074 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2075 v4sigma4[1] = 0.0e0;
2076
2077 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2078 v4sigma4[2] = 0.0e0;
2079
2080 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2081 v4sigma4[3] = 0.0e0;
2082
2083 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2084 v4sigma4[4] = 0.0e0;
2085
2086 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2087 v4sigma4[5] = 0.0e0;
2088
2089 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2090 v4sigma4[6] = 0.0e0;
2091
2092 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2093 v4sigma4[7] = 0.0e0;
2094
2095 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2096 v4sigma4[8] = 0.0e0;
2097
2098 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2099 v4sigma4[9] = 0.0e0;
2100
2101 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2102 v4sigma4[10] = 0.0e0;
2103
2104 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2105 v4sigma4[11] = 0.0e0;
2106
2107 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2108 v4sigma4[12] = 0.0e0;
2109
2110 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2111 v4sigma4[13] = 0.0e0;
2112
2113 t3426 = 0.1e1 / t457;
2114 t3479 = -0.7e1 / 0.46656e5 * t114 * t601 * t297 * t3426 * t76 + 0.7e1 / 0.11664e5 * t114 * t601 * params->alpha * t3426 * t76 + 0.7e1 / 0.5832e4 * t114 * t618 * t1475 * t452 + 0.7e1 / 0.1944e4 * t114 * t1498 * t603 * t1097 - 0.7e1 / 0.486e3 * t28 * t5 * t1504 * t1111 + 0.7e1 / 0.186624e6 * t114 * t601 * t1663 * t3426 * t76 - 0.7e1 / 0.46656e5 * t114 * t601 * t762 * t3426 * t76 - 0.7e1 / 0.23328e5 * t114 * t618 * t1488 * t452 + 0.7e1 / 0.3888e4 * t114 * t1498 * t612 * t1097 - 0.7e1 / 0.972e3 * t114 * t1504 * params->a1 * t69 * t2995 * t788 + 0.7e1 / 0.486e3 * t28 * t59 / t61 / t1104 * t2566;
2115 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2116 v4sigma4[14] = -0.3e1 / 0.8e1 * t200 * t223 * t3479;
2117
2118 #ifndef XC_DONT_COMPILE_MXC
2119
2120 if(order < 5) return;
2121
2122
2123 #endif
2124
2125 #endif
2126
2127 #endif
2128
2129 #endif
2130
2131 #endif
2132
2133
2134 }
2135
2136