1 /*
2 This file was generated automatically with ./scripts/maple2c.pl.
3 Do not edit this file directly as it can be overwritten!!
4
5 This Source Code Form is subject to the terms of the Mozilla Public
6 License, v. 2.0. If a copy of the MPL was not distributed with this
7 file, You can obtain one at http://mozilla.org/MPL/2.0/.
8
9 Maple version : Maple 2016 (X86 64 LINUX)
10 Maple source : ./maple/gga_exc/gga_c_wi.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t2, t3, t4, t6, t7, t9, t12;
24 double t13, t15, t16, t17, t18, t22, t23, t25;
25 double t26, t27, t28, t31, t32, t33, t38, t42;
26 double t43;
27
28 #ifndef XC_DONT_COMPILE_VXC
29 double t44, t46, t49, t50, t51, t53, t58, t59;
30 double t61, t62, t63, t71, t72, t73, t76, t77;
31 double t81, t83, t87, t88, t90, t91, t92, t93;
32 double t94, t95, t96, t97;
33
34 #ifndef XC_DONT_COMPILE_FXC
35 double t102, t106, t110, t112, t117, t118, t119, t120;
36 double t121, t122, t127, t128, t133, t134, t135, t148;
37 double t151, t152, t153, t156, t157, t163, t164, t167;
38 double t168, t173, t174, t178, t179, t180, t183, t184;
39 double t187, t188, t190, t193, t194, t196, t197, t201;
40 double t204, t208, t209, t211, t212, t215, t216, t217;
41 double t218, t219, t220, t222, t225, t226, t227, t228;
42 double t229, t230, t233, t234, t235;
43
44 #ifndef XC_DONT_COMPILE_KXC
45 double t240, t243, t249, t254, t259, t260, t265, t266;
46 double t267, t269, t270, t275, t276, t284, t285, t286;
47 double t287, t290, t291, t295, t312, t315, t316, t318;
48 double t319, t322, t325, t326, t328, t331, t332, t335;
49 double t336, t339, t341, t349, t355, t358, t360, t363;
50 double t365, t370, t371, t373, t374, t376, t379, t383;
51 double t389, t393, t399, t401, t405, t410, t411, t414;
52 double t415, t418, t419, t422, t425, t429, t430, t433;
53 double t434, t437, t438, t439, t441, t444, t447, t448;
54 double t451, t453, t457, t458, t459, t460, t461, t469;
55 double t473, t477, t481, t482, t484, t485, t488, t489;
56 double t492, t493, t498, t500, t501, t502, t506, t509;
57 double t510, t511, t512, t515, t516, t517, t521, t522;
58 double t523;
59
60 #ifndef XC_DONT_COMPILE_LXC
61 double t562, t563, t580, t590, t595, t604, t608, t623;
62 double t624, t632, t725, t785, t872, t943, t1011, t1071;
63 double t1118;
64 #endif
65
66 #endif
67
68 #endif
69
70 #endif
71
72 #endif
73
74
75 gga_c_wi_params *params;
76
77 assert(p->params != NULL);
78 params = (gga_c_wi_params * )(p->params);
79
80 t1 = params->b * sigma[0];
81 t2 = rho[0] * rho[0];
82 t3 = POW_1_3(rho[0]);
83 t4 = t3 * t3;
84 t6 = 0.1e1 / t4 / t2;
85 t7 = params->k * sigma[0];
86 t9 = exp(-t7 * t6);
87 t12 = t1 * t6 * t9 + params->a;
88 t13 = M_CBRT3;
89 t15 = POW_1_3(0.1e1 / M_PI);
90 t16 = t13 * t15;
91 t17 = M_CBRT4;
92 t18 = t17 * t17;
93 t22 = t13 * t13;
94 t23 = M_CBRTPI;
95 t25 = sqrt(sigma[0]);
96 t26 = t25 * sigma[0];
97 t27 = t2 * t2;
98 t28 = 0.1e1 / t27;
99 t31 = 0.1e1 / t3 / rho[0];
100 t32 = t25 * t31;
101 t33 = sqrt(t32);
102 t38 = 0.1e1 + params->d * t17 * t22 * t23 * t33 * t26 * t28 / 0.3e1;
103 t42 = params->c + t16 * t18 / t3 * t38 / 0.4e1;
104 t43 = 0.1e1 / t42;
105 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
106 zk[0] = t12 * t43;
107
108 #ifndef XC_DONT_COMPILE_VXC
109
110 if(order < 1) return;
111
112
113 t44 = t2 * rho[0];
114 t46 = 0.1e1 / t4 / t44;
115 t49 = sigma[0] * sigma[0];
116 t50 = params->b * t49;
117 t51 = t27 * t2;
118 t53 = 0.1e1 / t3 / t51;
119 t58 = 0.8e1 / 0.3e1 * t50 * t53 * params->k * t9 - 0.8e1 / 0.3e1 * t1 * t46 * t9;
120 t59 = rho[0] * t58;
121 t61 = rho[0] * t12;
122 t62 = t42 * t42;
123 t63 = 0.1e1 / t62;
124 t71 = t33 * sigma[0] * t6;
125 t72 = t23 * t71;
126 t73 = t72 * t25;
127 t76 = -t16 * t18 * t31 * t38 / 0.12e2 - 0.14e2 / 0.3e1 * t15 * t6 * params->d * t73;
128 t77 = t63 * t76;
129 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
130 vrho[0] = t59 * t43 - t61 * t77 + (t12 * t43);
131
132 t81 = t27 * rho[0];
133 t83 = 0.1e1 / t3 / t81;
134 t87 = -t1 * t83 * params->k * t9 + params->b * t6 * t9;
135 t88 = rho[0] * t87;
136 t90 = 0.1e1 / t4;
137 t91 = t90 * t12;
138 t92 = t63 * t15;
139 t93 = t91 * t92;
140 t94 = params->d * t23;
141 t95 = 0.1e1 / t25;
142 t96 = t71 * t95;
143 t97 = t94 * t96;
144 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
145 vsigma[0] = t88 * t43 - 0.7e1 / 0.4e1 * t93 * t97;
146
147 #ifndef XC_DONT_COMPILE_FXC
148
149 if(order < 2) return;
150
151
152 t102 = t12 * t63;
153 t106 = 0.1e1 / t4 / t27;
154 t110 = t27 * t44;
155 t112 = 0.1e1 / t3 / t110;
156 t117 = t49 * sigma[0];
157 t118 = params->b * t117;
158 t119 = t27 * t27;
159 t120 = t119 * t2;
160 t121 = 0.1e1 / t120;
161 t122 = params->k * params->k;
162 t127 = 0.88e2 / 0.9e1 * t1 * t106 * t9 - 0.24e2 * t50 * t112 * params->k * t9 + 0.64e2 / 0.9e1 * t118 * t121 * t122 * t9;
163 t128 = rho[0] * t127;
164 t133 = 0.1e1 / t62 / t42;
165 t134 = t76 * t76;
166 t135 = t133 * t134;
167 t148 = 0.1e1 / t81;
168 t151 = t33 * t32;
169 t152 = t23 * t151;
170 t153 = t152 * sigma[0];
171 t156 = t16 * t18 / t3 / t2 * t38 / 0.9e1 + 0.14e2 * t15 * t46 * params->d * t73 + 0.140e3 / 0.9e1 * t15 * t148 * params->d * t153;
172 t157 = t63 * t156;
173 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
174 v2rho2[0] = -0.2e1 * t102 * t76 + t128 * t43 + 0.2e1 * t61 * t135 - t61 * t157 + 0.2e1 * t58 * t43 - 0.2e1 * t59 * t77;
175
176 t163 = params->b * t53;
177 t164 = t7 * t9;
178 t167 = t119 * rho[0];
179 t168 = 0.1e1 / t167;
180 t173 = -0.8e1 / 0.3e1 * params->b * t46 * t9 + 0.8e1 * t163 * t164 - 0.8e1 / 0.3e1 * t50 * t168 * t122 * t9;
181 t174 = rho[0] * t173;
182 t178 = 0.1e1 / t4 / rho[0];
183 t179 = t178 * t12;
184 t180 = t179 * t92;
185 t183 = t90 * t58;
186 t184 = t183 * t92;
187 t187 = t133 * t15;
188 t188 = t91 * t187;
189 t190 = t94 * t96 * t76;
190 t193 = 0.1e1 / t44;
191 t194 = t193 * t12;
192 t196 = t15 * params->d;
193 t197 = t196 * t152;
194 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
195 v2rhosigma[0] = t87 * t43 + t174 * t43 - t88 * t77 + 0.7e1 / 0.6e1 * t180 * t97 - 0.7e1 / 0.4e1 * t184 * t97 + 0.7e1 / 0.2e1 * t188 * t190 + 0.35e2 / 0.6e1 * t194 * t63 * t197;
196
197 t201 = params->k * t9;
198 t204 = 0.1e1 / t119;
199 t208 = t1 * t204 * t122 * t9 - 0.2e1 * params->b * t83 * t201;
200 t209 = rho[0] * t208;
201 t211 = t90 * t87;
202 t212 = t211 * t92;
203 t215 = t168 * t12;
204 t216 = t215 * t133;
205 t217 = t15 * t15;
206 t218 = params->d * params->d;
207 t219 = t217 * t218;
208 t220 = t23 * t23;
209 t222 = t219 * t220 * t26;
210 t225 = 0.1e1 / t2;
211 t226 = t225 * t12;
212 t227 = t226 * t92;
213 t228 = 0.1e1 / sigma[0];
214 t229 = t151 * t228;
215 t230 = t94 * t229;
216 t233 = 0.1e1 / t26;
217 t234 = t71 * t233;
218 t235 = t94 * t234;
219 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
220 v2sigma2[0] = t209 * t43 - 0.7e1 / 0.2e1 * t212 * t97 + 0.49e2 / 0.8e1 * t216 * t222 - 0.35e2 / 0.16e2 * t227 * t230 + 0.7e1 / 0.8e1 * t93 * t235;
221
222 #ifndef XC_DONT_COMPILE_KXC
223
224 if(order < 3) return;
225
226
227 t240 = t58 * t63;
228 t243 = t12 * t133;
229 t249 = 0.1e1 / t4 / t81;
230 t254 = 0.1e1 / t3 / t119;
231 t259 = t119 * t44;
232 t260 = 0.1e1 / t259;
233 t265 = t49 * t49;
234 t266 = params->b * t265;
235 t267 = t119 * t81;
236 t269 = 0.1e1 / t4 / t267;
237 t270 = t122 * params->k;
238 t275 = -0.1232e4 / 0.27e2 * t1 * t249 * t9 + 0.5456e4 / 0.27e2 * t50 * t254 * params->k * t9 - 0.1216e4 / 0.9e1 * t118 * t260 * t122 * t9 + 0.512e3 / 0.27e2 * t266 * t269 * t270 * t9;
239 t276 = rho[0] * t275;
240 t284 = t62 * t62;
241 t285 = 0.1e1 / t284;
242 t286 = t134 * t76;
243 t287 = t285 * t286;
244 t290 = t133 * t76;
245 t291 = t290 * t156;
246 t295 = 0.1e1 / t3 / t44;
247 t312 = t23 * t33 * t26;
248 t315 = -0.7e1 / 0.27e2 * t16 * t18 * t295 * t38 - 0.1442e4 / 0.27e2 * t15 * t106 * params->d * t73 - 0.1120e4 / 0.9e1 * t15 / t51 * params->d * t153 - 0.280e3 / 0.9e1 * t15 * t112 * params->d * t312;
249 t316 = t63 * t315;
250 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
251 v3rho3[0] = -0.3e1 * t102 * t156 + 0.3e1 * t127 * t43 - 0.3e1 * t128 * t77 + 0.6e1 * t243 * t134 + 0.6e1 * t59 * t135 - 0.3e1 * t59 * t157 - 0.6e1 * t240 * t76 + t276 * t43 - 0.6e1 * t61 * t287 + 0.6e1 * t61 * t291 - t61 * t316;
252
253 t318 = t6 * t12;
254 t319 = t318 * t92;
255 t322 = t179 * t187;
256 t325 = t285 * t15;
257 t326 = t91 * t325;
258 t328 = t94 * t96 * t134;
259 t331 = t90 * t127;
260 t332 = t331 * t92;
261 t335 = t178 * t58;
262 t336 = t335 * t92;
263 t339 = t194 * t187;
264 t341 = t94 * t151 * t76;
265 t349 = t87 * t63;
266 t355 = params->b * t112;
267 t358 = params->b * t121;
268 t360 = t122 * t49 * t9;
269 t363 = t119 * t27;
270 t365 = 0.1e1 / t4 / t363;
271 t370 = 0.88e2 / 0.9e1 * params->b * t106 * t9 - 0.520e3 / 0.9e1 * t355 * t164 + 0.136e3 / 0.3e1 * t358 * t360 - 0.64e2 / 0.9e1 * t118 * t365 * t270 * t9;
272 t371 = rho[0] * t370;
273 t373 = t83 * t12;
274 t374 = t373 * t92;
275 t376 = t94 * t33 * t25;
276 t379 = t28 * t12;
277 t383 = t193 * t58;
278 t389 = t183 * t187;
279 t393 = t94 * t96 * t156;
280 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
281 v3rho2sigma[0] = -0.35e2 / 0.18e2 * t319 * t97 - 0.14e2 / 0.3e1 * t322 * t190 - 0.21e2 / 0.2e1 * t326 * t328 - 0.7e1 / 0.4e1 * t332 * t97 + 0.7e1 / 0.3e1 * t336 * t97 - 0.70e2 / 0.3e1 * t339 * t341 + 0.2e1 * t88 * t135 - 0.2e1 * t174 * t77 - t88 * t157 - 0.2e1 * t349 * t76 + t371 * t43 - 0.35e2 / 0.3e1 * t374 * t376 - 0.385e3 / 0.18e2 * t379 * t63 * t197 + 0.35e2 / 0.3e1 * t383 * t63 * t197 + 0.2e1 * t173 * t43 + 0.7e1 * t389 * t190 + 0.7e1 / 0.2e1 * t188 * t393;
282
283 t399 = params->b * t168;
284 t401 = t122 * sigma[0] * t9;
285 t405 = 0.1e1 / t4 / t259;
286 t410 = 0.32e2 / 0.3e1 * t163 * t201 - 0.40e2 / 0.3e1 * t399 * t401 + 0.8e1 / 0.3e1 * t50 * t405 * t270 * t9;
287 t411 = rho[0] * t410;
288 t414 = t178 * t87;
289 t415 = t414 * t92;
290 t418 = t90 * t173;
291 t419 = t418 * t92;
292 t422 = t211 * t187;
293 t425 = t193 * t87;
294 t429 = t121 * t12;
295 t430 = t429 * t133;
296 t433 = t168 * t58;
297 t434 = t433 * t133;
298 t437 = t285 * t217;
299 t438 = t215 * t437;
300 t439 = t218 * t220;
301 t441 = t439 * t26 * t76;
302 t444 = t194 * t92;
303 t447 = t225 * t58;
304 t448 = t447 * t92;
305 t451 = t226 * t187;
306 t453 = t94 * t229 * t76;
307 t457 = 0.1e1 / t3 / t27;
308 t458 = t457 * t12;
309 t459 = t458 * t92;
310 t460 = t33 * t95;
311 t461 = t94 * t460;
312 t469 = t94 * t234 * t76;
313 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
314 v3rhosigma2[0] = t208 * t43 + t411 * t43 - t209 * t77 + 0.7e1 / 0.3e1 * t415 * t97 - 0.7e1 / 0.2e1 * t419 * t97 + 0.7e1 * t422 * t190 + 0.35e2 / 0.3e1 * t425 * t63 * t197 - 0.441e3 / 0.8e1 * t430 * t222 + 0.49e2 / 0.8e1 * t434 * t222 - 0.147e3 / 0.8e1 * t438 * t441 + 0.35e2 / 0.24e2 * t444 * t230 - 0.35e2 / 0.16e2 * t448 * t230 + 0.35e2 / 0.8e1 * t451 * t453 + 0.35e2 / 0.8e1 * t459 * t461 - 0.7e1 / 0.12e2 * t180 * t235 + 0.7e1 / 0.8e1 * t184 * t235 - 0.7e1 / 0.4e1 * t188 * t469;
315
316 t473 = t122 * t9;
317 t477 = 0.1e1 / t4 / t120;
318 t481 = -t1 * t477 * t270 * t9 + 0.3e1 * params->b * t204 * t473;
319 t482 = rho[0] * t481;
320 t484 = t90 * t208;
321 t485 = t484 * t92;
322 t488 = t168 * t87;
323 t489 = t488 * t133;
324 t492 = t225 * t87;
325 t493 = t492 * t92;
326 t498 = t477 * t12;
327 t500 = t218 * params->d;
328 t501 = t500 * sigma[0];
329 t502 = t501 * t71;
330 t506 = t219 * t220 * t25;
331 t509 = t295 * t12;
332 t510 = t509 * t92;
333 t511 = t33 * t233;
334 t512 = t94 * t511;
335 t515 = 0.1e1 / t49;
336 t516 = t151 * t515;
337 t517 = t94 * t516;
338 t521 = 0.1e1 / t25 / t49;
339 t522 = t71 * t521;
340 t523 = t94 * t522;
341 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
342 v3sigma3[0] = t482 * t43 - 0.21e2 / 0.4e1 * t485 * t97 + 0.147e3 / 0.8e1 * t489 * t222 - 0.105e3 / 0.16e2 * t493 * t230 + 0.21e2 / 0.8e1 * t212 * t235 - 0.1029e4 / 0.32e2 * t498 * t285 * t502 + 0.441e3 / 0.32e2 * t216 * t506 - 0.105e3 / 0.64e2 * t510 * t512 + 0.105e3 / 0.32e2 * t227 * t517 - 0.21e2 / 0.16e2 * t93 * t523;
343
344 #ifndef XC_DONT_COMPILE_LXC
345
346 if(order < 4) return;
347
348
349 t562 = 0.1e1 / t33;
350 t563 = t23 * t562;
351 t580 = 0.1e1 / t4 / t51;
352 t590 = 0.1e1 / t363;
353 t595 = t119 * t51;
354 t604 = t119 * t119;
355 t608 = t122 * t122;
356 t623 = 0.1e1 / t284 / t42;
357 t624 = t134 * t134;
358 t632 = t156 * t156;
359 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
360 v4rho4[0] = 0.24e2 * t58 * t133 * t134 - 0.24e2 * t12 * t285 * t286 + 0.24e2 * t243 * t76 * t156 - 0.4e1 * t276 * t77 - 0.6e1 * t128 * t157 - 0.4e1 * t59 * t316 - t61 * t63 * (0.70e2 / 0.81e2 * t16 * t18 * t457 * t38 + 0.6860e4 / 0.27e2 * t15 * t249 * params->d * t73 + 0.74900e5 / 0.81e2 * t15 / t110 * params->d * t153 + 0.12880e5 / 0.27e2 * t15 * t254 * params->d * t312 + 0.560e3 / 0.27e2 * t15 / t4 / t167 * params->d * t563 * t49) + 0.4e1 * t275 * t43 - 0.12e2 * t127 * t63 * t76 - 0.12e2 * t240 * t156 - 0.4e1 * t102 * t315 + rho[0] * (0.20944e5 / 0.81e2 * t1 * t580 * t9 - 0.48752e5 / 0.27e2 * t50 / t3 / t167 * params->k * t9 + 0.164032e6 / 0.81e2 * t118 * t590 * t122 * t9 - 0.50176e5 / 0.81e2 * t266 / t4 / t595 * t270 * t9 + 0.4096e4 / 0.81e2 * params->b * t265 * sigma[0] / t3 / t604 / rho[0] * t608 * t9) * t43 + 0.12e2 * t128 * t135 - 0.24e2 * t59 * t287 + 0.24e2 * t59 * t291 + 0.24e2 * t61 * t623 * t624 - 0.36e2 * t61 * t285 * t134 * t156 + 0.6e1 * t61 * t133 * t632 + 0.8e1 * t61 * t290 * t315;
361
362 t725 = 0.6e1 * t88 * t291 - 0.3e1 * t349 * t156 + rho[0] * (-0.1232e4 / 0.27e2 * params->b * t249 * t9 + 0.4048e4 / 0.9e1 * params->b * t254 * t164 - 0.16400e5 / 0.27e2 * params->b * t260 * t360 + 0.5696e4 / 0.27e2 * params->b * t269 * t270 * t117 * t9 - 0.512e3 / 0.27e2 * t266 / t3 / t604 * t608 * t9) * t43 + 0.6e1 * t87 * t133 * t134 - 0.6e1 * t173 * t63 * t76 + 0.3e1 * t370 * t43 - 0.35e2 / 0.6e1 * t6 * t58 * t92 * t97 + 0.385e3 / 0.3e1 * t379 * t187 * t341 + 0.105e3 * t194 * t325 * t94 * t151 * t134 - 0.7e1 / 0.4e1 * t90 * t275 * t92 * t97 - 0.70e2 * t383 * t187 * t341 - 0.35e2 * t339 * t94 * t151 * t156 - 0.35e2 * t83 * t58 * t92 * t376 + 0.70e2 / 0.9e1 / t4 / t110 * t12 * t92 * t94 * t562 * sigma[0] + 0.140e3 / 0.27e2 * t46 * t12 * t92 * t97 + 0.7e1 / 0.2e1 * t178 * t127 * t92 * t97 + 0.105e3 * t53 * t12 * t92 * t376 + 0.2485e4 / 0.27e2 * t148 * t12 * t63 * t197;
363 t785 = 0.35e2 / 0.2e1 * t193 * t127 * t63 * t197 - 0.385e3 / 0.6e1 * t28 * t58 * t63 * t197 + 0.35e2 / 0.3e1 * t318 * t187 * t190 - 0.14e2 * t335 * t187 * t190 - 0.7e1 * t322 * t393 - 0.63e2 / 0.2e1 * t183 * t325 * t328 + 0.21e2 / 0.2e1 * t331 * t187 * t190 + 0.70e2 * t373 * t187 * t94 * t33 * t76 * t25 + 0.21e2 / 0.2e1 * t389 * t393 + 0.7e1 / 0.2e1 * t188 * t94 * t96 * t315 + 0.21e2 * t179 * t325 * t328 + 0.42e2 * t91 * t623 * t15 * t94 * t96 * t286 - 0.63e2 / 0.2e1 * t91 * t325 * params->d * t72 * t95 * t76 * t156 + 0.6e1 * t174 * t135 - 0.3e1 * t371 * t77 - 0.3e1 * t174 * t157 - t88 * t316 - 0.6e1 * t88 * t287;
364 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
365 v4rho3sigma[0] = t725 + t785;
366
367 t872 = 0.147e3 / 0.2e1 * t215 * t623 * t217 * t439 * t26 * t134 - 0.175e3 / 0.72e2 * t379 * t92 * t230 + 0.35e2 / 0.36e2 * t319 * t235 + 0.14e2 / 0.3e1 * t178 * t173 * t92 * t97 - 0.7e1 / 0.2e1 * t90 * t370 * t92 * t97 - 0.140e3 / 0.3e1 * t425 * t187 * t341 - 0.70e2 / 0.3e1 * t83 * t87 * t92 * t376 + 0.1323e4 / 0.4e1 * t429 * t437 * t441 - 0.147e3 / 0.4e1 * t433 * t437 * t441 - 0.147e3 / 0.8e1 * t438 * t439 * t26 * t156 + 0.2205e4 / 0.4e1 * t260 * t12 * t133 * t222;
368 t943 = 0.35e2 / 0.8e1 * t451 * t94 * t229 * t156 - 0.35e2 / 0.2e1 * t458 * t187 * t94 * t460 * t76 + 0.7e1 / 0.3e1 * t322 * t469 - 0.7e1 / 0.2e1 * t389 * t469 - 0.7e1 / 0.4e1 * t188 * t94 * t234 * t156 - 0.21e2 * t211 * t325 * t328 - 0.105e3 / 0.8e1 * t226 * t325 * t94 * t229 * t134 + 0.21e2 / 0.4e1 * t326 * t94 * t234 * t134 - 0.2e1 * t411 * t77 - t209 * t157 + 0.2e1 * t209 * t135;
369 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
370 v4rho2sigma2[0] = -0.2e1 * t208 * t63 * t76 + rho[0] * (-0.608e3 / 0.9e1 * t355 * t201 + 0.1336e4 / 0.9e1 * t358 * t401 - 0.200e3 / 0.3e1 * params->b * t365 * t270 * t49 * t9 + 0.64e2 / 0.9e1 * t118 / t3 / t119 / t110 * t608 * t9) * t43 + 0.2e1 * t410 * t43 + 0.35e2 / 0.12e2 * t383 * t92 * t230 - 0.175e3 / 0.8e1 * t374 * t461 - 0.35e2 / 0.16e2 * t225 * t127 * t92 * t230 + 0.35e2 / 0.4e1 * t457 * t58 * t92 * t461 - 0.7e1 / 0.6e1 * t336 * t235 + 0.7e1 / 0.8e1 * t332 * t235 - 0.35e2 / 0.9e1 * t6 * t87 * t92 * t97 + t872 - 0.385e3 / 0.9e1 * t28 * t87 * t63 * t197 + 0.70e2 / 0.3e1 * t193 * t173 * t63 * t197 - 0.441e3 / 0.4e1 * t121 * t58 * t133 * t222 + 0.49e2 / 0.8e1 * t168 * t127 * t133 * t222 - 0.35e2 / 0.12e2 * t580 * t12 * t63 * t196 * t563 - 0.28e2 / 0.3e1 * t414 * t187 * t190 + 0.14e2 * t418 * t187 * t190 + 0.7e1 * t422 * t393 - 0.35e2 / 0.6e1 * t339 * t453 + 0.35e2 / 0.4e1 * t447 * t187 * t453 + t943;
371
372 t1011 = t481 * t43 - 0.1029e4 / 0.32e2 * t477 * t58 * t285 * t502 + 0.1715e4 / 0.16e2 / t267 * t12 * t285 * t500 * t26 * t151 + 0.343e3 * t405 * t12 * t285 * t502 + rho[0] * (-0.24e2 * t399 * t473 + 0.56e2 / 0.3e1 * params->b * t405 * t270 * sigma[0] * t9 - 0.8e1 / 0.3e1 * t50 / t3 / t595 * t608 * t9) * t43 + 0.105e3 / 0.32e2 * t448 * t517 - 0.21e2 / 0.16e2 * t184 * t523 + 0.7e1 / 0.2e1 * t178 * t208 * t92 * t97 + 0.35e2 / 0.8e1 * t425 * t92 * t230 - 0.7e1 / 0.4e1 * t415 * t235 - 0.35e2 / 0.32e2 * t459 * t512 - 0.35e2 / 0.16e2 * t444 * t517 + 0.7e1 / 0.8e1 * t180 * t523 - 0.21e2 / 0.4e1 * t90 * t410 * t92 * t97 - 0.441e3 / 0.8e1 * t488 * t437 * t441 - 0.105e3 / 0.16e2 * t225 * t173 * t92 * t230 + 0.105e3 / 0.8e1 * t457 * t87 * t92 * t461;
373 t1071 = 0.21e2 / 0.8e1 * t419 * t235 - 0.1323e4 / 0.32e2 * t438 * t439 * t25 * t76 - 0.105e3 / 0.64e2 * t295 * t58 * t92 * t512 + 0.35e2 / 0.32e2 * t249 * t12 * t92 * t94 * t562 * t228 - 0.1323e4 / 0.8e1 * t121 * t87 * t133 * t222 - 0.3969e4 / 0.32e2 * t430 * t506 + 0.35e2 / 0.2e1 * t193 * t208 * t63 * t197 + 0.147e3 / 0.8e1 * t168 * t173 * t133 * t222 + 0.1029e4 / 0.8e1 * t498 * t623 * t501 * t71 * t76 + 0.441e3 / 0.32e2 * t434 * t506 + 0.21e2 / 0.2e1 * t484 * t187 * t190 + 0.105e3 / 0.8e1 * t492 * t187 * t453 - 0.21e2 / 0.4e1 * t422 * t469 + 0.105e3 / 0.32e2 * t509 * t187 * t94 * t511 * t76 - 0.105e3 / 0.16e2 * t451 * t94 * t516 * t76 + 0.21e2 / 0.8e1 * t188 * t94 * t522 * t76 - t482 * t77;
374 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
375 v4rhosigma3[0] = t1011 + t1071;
376
377 t1118 = t218 * t218;
378 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
379 v4sigma4[0] = -0.7e1 * t90 * t481 * t92 * t97 - 0.105e3 / 0.8e1 * t225 * t208 * t92 * t230 + 0.441e3 / 0.8e1 * t489 * t506 - 0.105e3 / 0.16e2 * t295 * t87 * t92 * t512 + 0.735e3 / 0.128e3 * t216 * t219 * t220 * t95 - 0.105e3 / 0.256e3 * t106 * t12 * t92 * t94 * t562 * t515 - 0.1029e4 / 0.8e1 * t477 * t87 * t285 * t502 - 0.13377e5 / 0.128e3 * t498 * t285 * t500 * t71 - 0.5145e4 / 0.128e3 * t590 * t12 * t285 * t500 * t25 * t151 + 0.147e3 / 0.4e1 * t168 * t208 * t133 * t222 + 0.7203e4 / 0.32e2 / t604 / t44 * t12 * t623 * t1118 * t117 * t15 * t23 + 0.21e2 / 0.4e1 * t485 * t235 + 0.105e3 / 0.8e1 * t493 * t517 - 0.21e2 / 0.4e1 * t212 * t523 + 0.315e3 / 0.64e2 * t510 * t94 * t33 * t521 - 0.525e3 / 0.64e2 * t227 * t94 * t151 / t117 + 0.105e3 / 0.32e2 * t93 * t94 * t71 / t25 / t117 + rho[0] * (-0.4e1 * params->b * t477 * t270 * t9 + t1 / t3 / t267 * t608 * t9) * t43;
380
381 #ifndef XC_DONT_COMPILE_MXC
382
383 if(order < 5) return;
384
385
386 #endif
387
388 #endif
389
390 #endif
391
392 #endif
393
394 #endif
395
396
397 }
398
399
400 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 *))401 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
402 {
403
404 #ifndef XC_DONT_COMPILE_EXC
405 double t1, t2, t3, t4, t6, t7, t9, t12;
406 double t13, t15, t16, t17, t18, t22, t23, t25;
407 double t26, t27, t28, t31, t32, t33, t38, t42;
408 double t43;
409
410 #ifndef XC_DONT_COMPILE_VXC
411 double t44, t46, t49, t50, t51, t53, t58, t59;
412 double t61, t62, t63, t71, t72, t73, t76, t77;
413 double t81, t83, t87, t88, t90, t91, t92, t93;
414 double t94, t95, t96, t97;
415
416 #ifndef XC_DONT_COMPILE_FXC
417 double t102, t106, t110, t112, t117, t118, t119, t120;
418 double t121, t122, t127, t128, t133, t134, t135, t148;
419 double t151, t152, t153, t156, t157, t163, t164, t167;
420 double t168, t173, t174, t178, t179, t180, t183, t184;
421 double t187, t188, t190, t193, t194, t196, t197, t201;
422 double t204, t208, t209, t211, t212, t215, t216, t217;
423 double t218, t219, t220, t222, t225, t226, t227, t228;
424 double t229, t230, t233, t234, t235;
425
426 #ifndef XC_DONT_COMPILE_KXC
427 double t240, t243, t249, t254, t259, t260, t265, t266;
428 double t267, t269, t270, t275, t276, t284, t285, t286;
429 double t287, t290, t291, t295, t312, t315, t316, t318;
430 double t319, t322, t325, t326, t328, t333, t335, t343;
431 double t349, t352, t354, t357, t359, t364, t365, t367;
432 double t368, t371, t372, t375, t376, t378, t381, t385;
433 double t389, t393, t399, t401, t405, t410, t411, t414;
434 double t415, t418, t419, t422, t425, t429, t430, t433;
435 double t434, t437, t438, t439, t441, t444, t447, t448;
436 double t451, t453, t457, t458, t459, t460, t461, t469;
437 double t473, t477, t481, t482, t484, t485, t488, t489;
438 double t492, t493, t498, t500, t501, t502, t506, t509;
439 double t510, t511, t512, t515, t516, t517, t521, t522;
440 double t523;
441
442 #ifndef XC_DONT_COMPILE_LXC
443 double t562, t563, t578, t588, t593, t602, t606, t623;
444 double t624, t632, t710, t785, t869, t943, t1016, t1071;
445 double t1118;
446 #endif
447
448 #endif
449
450 #endif
451
452 #endif
453
454 #endif
455
456
457 gga_c_wi_params *params;
458
459 assert(p->params != NULL);
460 params = (gga_c_wi_params * )(p->params);
461
462 t1 = params->b * sigma[0];
463 t2 = rho[0] * rho[0];
464 t3 = POW_1_3(rho[0]);
465 t4 = t3 * t3;
466 t6 = 0.1e1 / t4 / t2;
467 t7 = params->k * sigma[0];
468 t9 = exp(-t7 * t6);
469 t12 = t1 * t6 * t9 + params->a;
470 t13 = M_CBRT3;
471 t15 = POW_1_3(0.1e1 / M_PI);
472 t16 = t13 * t15;
473 t17 = M_CBRT4;
474 t18 = t17 * t17;
475 t22 = t13 * t13;
476 t23 = M_CBRTPI;
477 t25 = sqrt(sigma[0]);
478 t26 = t25 * sigma[0];
479 t27 = t2 * t2;
480 t28 = 0.1e1 / t27;
481 t31 = 0.1e1 / t3 / rho[0];
482 t32 = t25 * t31;
483 t33 = sqrt(t32);
484 t38 = 0.1e1 + params->d * t17 * t22 * t23 * t33 * t26 * t28 / 0.3e1;
485 t42 = params->c + t16 * t18 / t3 * t38 / 0.4e1;
486 t43 = 0.1e1 / t42;
487 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
488 zk[0] = t12 * t43;
489
490 #ifndef XC_DONT_COMPILE_VXC
491
492 if(order < 1) return;
493
494
495 t44 = t2 * rho[0];
496 t46 = 0.1e1 / t4 / t44;
497 t49 = sigma[0] * sigma[0];
498 t50 = params->b * t49;
499 t51 = t27 * t2;
500 t53 = 0.1e1 / t3 / t51;
501 t58 = 0.8e1 / 0.3e1 * t50 * t53 * params->k * t9 - 0.8e1 / 0.3e1 * t1 * t46 * t9;
502 t59 = rho[0] * t58;
503 t61 = rho[0] * t12;
504 t62 = t42 * t42;
505 t63 = 0.1e1 / t62;
506 t71 = t33 * sigma[0] * t6;
507 t72 = t23 * t71;
508 t73 = t72 * t25;
509 t76 = -t16 * t18 * t31 * t38 / 0.12e2 - 0.14e2 / 0.3e1 * t15 * t6 * params->d * t73;
510 t77 = t63 * t76;
511 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
512 vrho[0] = t59 * t43 - t61 * t77 + (t12 * t43);
513
514 t81 = t27 * rho[0];
515 t83 = 0.1e1 / t3 / t81;
516 t87 = -t1 * t83 * params->k * t9 + params->b * t6 * t9;
517 t88 = rho[0] * t87;
518 t90 = 0.1e1 / t4;
519 t91 = t90 * t12;
520 t92 = t63 * t15;
521 t93 = t91 * t92;
522 t94 = params->d * t23;
523 t95 = 0.1e1 / t25;
524 t96 = t71 * t95;
525 t97 = t94 * t96;
526 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
527 vsigma[0] = t88 * t43 - 0.7e1 / 0.4e1 * t93 * t97;
528
529 #ifndef XC_DONT_COMPILE_FXC
530
531 if(order < 2) return;
532
533
534 t102 = t12 * t63;
535 t106 = 0.1e1 / t4 / t27;
536 t110 = t27 * t44;
537 t112 = 0.1e1 / t3 / t110;
538 t117 = t49 * sigma[0];
539 t118 = params->b * t117;
540 t119 = t27 * t27;
541 t120 = t119 * t2;
542 t121 = 0.1e1 / t120;
543 t122 = params->k * params->k;
544 t127 = 0.88e2 / 0.9e1 * t1 * t106 * t9 - 0.24e2 * t50 * t112 * params->k * t9 + 0.64e2 / 0.9e1 * t118 * t121 * t122 * t9;
545 t128 = rho[0] * t127;
546 t133 = 0.1e1 / t62 / t42;
547 t134 = t76 * t76;
548 t135 = t133 * t134;
549 t148 = 0.1e1 / t81;
550 t151 = t33 * t32;
551 t152 = t23 * t151;
552 t153 = t152 * sigma[0];
553 t156 = t16 * t18 / t3 / t2 * t38 / 0.9e1 + 0.14e2 * t15 * t46 * params->d * t73 + 0.140e3 / 0.9e1 * t15 * t148 * params->d * t153;
554 t157 = t63 * t156;
555 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
556 v2rho2[0] = -0.2e1 * t102 * t76 + t128 * t43 + 0.2e1 * t61 * t135 - t61 * t157 + 0.2e1 * t58 * t43 - 0.2e1 * t59 * t77;
557
558 t163 = params->b * t53;
559 t164 = t7 * t9;
560 t167 = t119 * rho[0];
561 t168 = 0.1e1 / t167;
562 t173 = -0.8e1 / 0.3e1 * params->b * t46 * t9 + 0.8e1 * t163 * t164 - 0.8e1 / 0.3e1 * t50 * t168 * t122 * t9;
563 t174 = rho[0] * t173;
564 t178 = 0.1e1 / t4 / rho[0];
565 t179 = t178 * t12;
566 t180 = t179 * t92;
567 t183 = t90 * t58;
568 t184 = t183 * t92;
569 t187 = t133 * t15;
570 t188 = t91 * t187;
571 t190 = t94 * t96 * t76;
572 t193 = 0.1e1 / t44;
573 t194 = t193 * t12;
574 t196 = t15 * params->d;
575 t197 = t196 * t152;
576 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
577 v2rhosigma[0] = t87 * t43 + t174 * t43 - t88 * t77 + 0.7e1 / 0.6e1 * t180 * t97 - 0.7e1 / 0.4e1 * t184 * t97 + 0.7e1 / 0.2e1 * t188 * t190 + 0.35e2 / 0.6e1 * t194 * t63 * t197;
578
579 t201 = params->k * t9;
580 t204 = 0.1e1 / t119;
581 t208 = t1 * t204 * t122 * t9 - 0.2e1 * params->b * t83 * t201;
582 t209 = rho[0] * t208;
583 t211 = t90 * t87;
584 t212 = t211 * t92;
585 t215 = t168 * t12;
586 t216 = t215 * t133;
587 t217 = t15 * t15;
588 t218 = params->d * params->d;
589 t219 = t217 * t218;
590 t220 = t23 * t23;
591 t222 = t219 * t220 * t26;
592 t225 = 0.1e1 / t2;
593 t226 = t225 * t12;
594 t227 = t226 * t92;
595 t228 = 0.1e1 / sigma[0];
596 t229 = t151 * t228;
597 t230 = t94 * t229;
598 t233 = 0.1e1 / t26;
599 t234 = t71 * t233;
600 t235 = t94 * t234;
601 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
602 v2sigma2[0] = t209 * t43 - 0.7e1 / 0.2e1 * t212 * t97 + 0.49e2 / 0.8e1 * t216 * t222 - 0.35e2 / 0.16e2 * t227 * t230 + 0.7e1 / 0.8e1 * t93 * t235;
603
604 #ifndef XC_DONT_COMPILE_KXC
605
606 if(order < 3) return;
607
608
609 t240 = t58 * t63;
610 t243 = t12 * t133;
611 t249 = 0.1e1 / t4 / t81;
612 t254 = 0.1e1 / t3 / t119;
613 t259 = t119 * t44;
614 t260 = 0.1e1 / t259;
615 t265 = t49 * t49;
616 t266 = params->b * t265;
617 t267 = t119 * t81;
618 t269 = 0.1e1 / t4 / t267;
619 t270 = t122 * params->k;
620 t275 = -0.1232e4 / 0.27e2 * t1 * t249 * t9 + 0.5456e4 / 0.27e2 * t50 * t254 * params->k * t9 - 0.1216e4 / 0.9e1 * t118 * t260 * t122 * t9 + 0.512e3 / 0.27e2 * t266 * t269 * t270 * t9;
621 t276 = rho[0] * t275;
622 t284 = t62 * t62;
623 t285 = 0.1e1 / t284;
624 t286 = t134 * t76;
625 t287 = t285 * t286;
626 t290 = t133 * t76;
627 t291 = t290 * t156;
628 t295 = 0.1e1 / t3 / t44;
629 t312 = t23 * t33 * t26;
630 t315 = -0.7e1 / 0.27e2 * t16 * t18 * t295 * t38 - 0.1442e4 / 0.27e2 * t15 * t106 * params->d * t73 - 0.1120e4 / 0.9e1 * t15 / t51 * params->d * t153 - 0.280e3 / 0.9e1 * t15 * t112 * params->d * t312;
631 t316 = t63 * t315;
632 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
633 v3rho3[0] = -0.3e1 * t102 * t156 + 0.3e1 * t127 * t43 - 0.3e1 * t128 * t77 + 0.6e1 * t243 * t134 + 0.6e1 * t59 * t135 - 0.3e1 * t59 * t157 - 0.6e1 * t240 * t76 + t276 * t43 - 0.6e1 * t61 * t287 + 0.6e1 * t61 * t291 - t61 * t316;
634
635 t318 = t6 * t12;
636 t319 = t318 * t92;
637 t322 = t179 * t187;
638 t325 = t285 * t15;
639 t326 = t91 * t325;
640 t328 = t94 * t96 * t134;
641 t333 = t194 * t187;
642 t335 = t94 * t151 * t76;
643 t343 = t87 * t63;
644 t349 = params->b * t112;
645 t352 = params->b * t121;
646 t354 = t122 * t49 * t9;
647 t357 = t119 * t27;
648 t359 = 0.1e1 / t4 / t357;
649 t364 = 0.88e2 / 0.9e1 * params->b * t106 * t9 - 0.520e3 / 0.9e1 * t349 * t164 + 0.136e3 / 0.3e1 * t352 * t354 - 0.64e2 / 0.9e1 * t118 * t359 * t270 * t9;
650 t365 = rho[0] * t364;
651 t367 = t90 * t127;
652 t368 = t367 * t92;
653 t371 = t178 * t58;
654 t372 = t371 * t92;
655 t375 = t83 * t12;
656 t376 = t375 * t92;
657 t378 = t94 * t33 * t25;
658 t381 = t28 * t12;
659 t385 = t193 * t58;
660 t389 = t183 * t187;
661 t393 = t94 * t96 * t156;
662 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
663 v3rho2sigma[0] = -0.35e2 / 0.18e2 * t319 * t97 - 0.14e2 / 0.3e1 * t322 * t190 - 0.21e2 / 0.2e1 * t326 * t328 + 0.2e1 * t173 * t43 - 0.70e2 / 0.3e1 * t333 * t335 + 0.2e1 * t88 * t135 - 0.2e1 * t174 * t77 - t88 * t157 - 0.2e1 * t343 * t76 + t365 * t43 - 0.7e1 / 0.4e1 * t368 * t97 + 0.7e1 / 0.3e1 * t372 * t97 - 0.35e2 / 0.3e1 * t376 * t378 - 0.385e3 / 0.18e2 * t381 * t63 * t197 + 0.35e2 / 0.3e1 * t385 * t63 * t197 + 0.7e1 * t389 * t190 + 0.7e1 / 0.2e1 * t188 * t393;
664
665 t399 = params->b * t168;
666 t401 = t122 * sigma[0] * t9;
667 t405 = 0.1e1 / t4 / t259;
668 t410 = 0.32e2 / 0.3e1 * t163 * t201 - 0.40e2 / 0.3e1 * t399 * t401 + 0.8e1 / 0.3e1 * t50 * t405 * t270 * t9;
669 t411 = rho[0] * t410;
670 t414 = t178 * t87;
671 t415 = t414 * t92;
672 t418 = t90 * t173;
673 t419 = t418 * t92;
674 t422 = t211 * t187;
675 t425 = t193 * t87;
676 t429 = t121 * t12;
677 t430 = t429 * t133;
678 t433 = t168 * t58;
679 t434 = t433 * t133;
680 t437 = t285 * t217;
681 t438 = t215 * t437;
682 t439 = t218 * t220;
683 t441 = t439 * t26 * t76;
684 t444 = t194 * t92;
685 t447 = t225 * t58;
686 t448 = t447 * t92;
687 t451 = t226 * t187;
688 t453 = t94 * t229 * t76;
689 t457 = 0.1e1 / t3 / t27;
690 t458 = t457 * t12;
691 t459 = t458 * t92;
692 t460 = t33 * t95;
693 t461 = t94 * t460;
694 t469 = t94 * t234 * t76;
695 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
696 v3rhosigma2[0] = t208 * t43 + t411 * t43 - t209 * t77 + 0.7e1 / 0.3e1 * t415 * t97 - 0.7e1 / 0.2e1 * t419 * t97 + 0.7e1 * t422 * t190 + 0.35e2 / 0.3e1 * t425 * t63 * t197 - 0.441e3 / 0.8e1 * t430 * t222 + 0.49e2 / 0.8e1 * t434 * t222 - 0.147e3 / 0.8e1 * t438 * t441 + 0.35e2 / 0.24e2 * t444 * t230 - 0.35e2 / 0.16e2 * t448 * t230 + 0.35e2 / 0.8e1 * t451 * t453 + 0.35e2 / 0.8e1 * t459 * t461 - 0.7e1 / 0.12e2 * t180 * t235 + 0.7e1 / 0.8e1 * t184 * t235 - 0.7e1 / 0.4e1 * t188 * t469;
697
698 t473 = t122 * t9;
699 t477 = 0.1e1 / t4 / t120;
700 t481 = -t1 * t477 * t270 * t9 + 0.3e1 * params->b * t204 * t473;
701 t482 = rho[0] * t481;
702 t484 = t90 * t208;
703 t485 = t484 * t92;
704 t488 = t168 * t87;
705 t489 = t488 * t133;
706 t492 = t225 * t87;
707 t493 = t492 * t92;
708 t498 = t477 * t12;
709 t500 = t218 * params->d;
710 t501 = t500 * sigma[0];
711 t502 = t501 * t71;
712 t506 = t219 * t220 * t25;
713 t509 = t295 * t12;
714 t510 = t509 * t92;
715 t511 = t33 * t233;
716 t512 = t94 * t511;
717 t515 = 0.1e1 / t49;
718 t516 = t151 * t515;
719 t517 = t94 * t516;
720 t521 = 0.1e1 / t25 / t49;
721 t522 = t71 * t521;
722 t523 = t94 * t522;
723 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
724 v3sigma3[0] = t482 * t43 - 0.21e2 / 0.4e1 * t485 * t97 + 0.147e3 / 0.8e1 * t489 * t222 - 0.105e3 / 0.16e2 * t493 * t230 + 0.21e2 / 0.8e1 * t212 * t235 - 0.1029e4 / 0.32e2 * t498 * t285 * t502 + 0.441e3 / 0.32e2 * t216 * t506 - 0.105e3 / 0.64e2 * t510 * t512 + 0.105e3 / 0.32e2 * t227 * t517 - 0.21e2 / 0.16e2 * t93 * t523;
725
726 #ifndef XC_DONT_COMPILE_LXC
727
728 if(order < 4) return;
729
730
731 t562 = 0.1e1 / t33;
732 t563 = t23 * t562;
733 t578 = 0.1e1 / t4 / t51;
734 t588 = 0.1e1 / t357;
735 t593 = t119 * t51;
736 t602 = t119 * t119;
737 t606 = t122 * t122;
738 t623 = 0.1e1 / t284 / t42;
739 t624 = t134 * t134;
740 t632 = t156 * t156;
741 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
742 v4rho4[0] = 0.24e2 * t58 * t133 * t134 - 0.24e2 * t12 * t285 * t286 + 0.24e2 * t243 * t76 * t156 - 0.4e1 * t276 * t77 - 0.6e1 * t128 * t157 - 0.4e1 * t59 * t316 - t61 * t63 * (0.70e2 / 0.81e2 * t16 * t18 * t457 * t38 + 0.6860e4 / 0.27e2 * t15 * t249 * params->d * t73 + 0.74900e5 / 0.81e2 * t15 / t110 * params->d * t153 + 0.12880e5 / 0.27e2 * t15 * t254 * params->d * t312 + 0.560e3 / 0.27e2 * t15 / t4 / t167 * params->d * t563 * t49) - 0.12e2 * t127 * t63 * t76 - 0.12e2 * t240 * t156 - 0.4e1 * t102 * t315 + rho[0] * (0.20944e5 / 0.81e2 * t1 * t578 * t9 - 0.48752e5 / 0.27e2 * t50 / t3 / t167 * params->k * t9 + 0.164032e6 / 0.81e2 * t118 * t588 * t122 * t9 - 0.50176e5 / 0.81e2 * t266 / t4 / t593 * t270 * t9 + 0.4096e4 / 0.81e2 * params->b * t265 * sigma[0] / t3 / t602 / rho[0] * t606 * t9) * t43 + 0.4e1 * t275 * t43 + 0.12e2 * t128 * t135 - 0.24e2 * t59 * t287 + 0.24e2 * t59 * t291 + 0.24e2 * t61 * t623 * t624 - 0.36e2 * t61 * t285 * t134 * t156 + 0.6e1 * t61 * t133 * t632 + 0.8e1 * t61 * t290 * t315;
743
744 t710 = 0.6e1 * t88 * t291 - 0.6e1 * t173 * t63 * t76 + rho[0] * (-0.1232e4 / 0.27e2 * params->b * t249 * t9 + 0.4048e4 / 0.9e1 * params->b * t254 * t164 - 0.16400e5 / 0.27e2 * params->b * t260 * t354 + 0.5696e4 / 0.27e2 * params->b * t269 * t270 * t117 * t9 - 0.512e3 / 0.27e2 * t266 / t3 / t602 * t606 * t9) * t43 + 0.3e1 * t364 * t43 + 0.6e1 * t174 * t135 - 0.3e1 * t365 * t77 - 0.3e1 * t174 * t157 - t88 * t316 - 0.6e1 * t88 * t287 + 0.2485e4 / 0.27e2 * t148 * t12 * t63 * t197 + 0.35e2 / 0.2e1 * t193 * t127 * t63 * t197 - 0.385e3 / 0.6e1 * t28 * t58 * t63 * t197 + 0.35e2 / 0.3e1 * t318 * t187 * t190 - 0.14e2 * t371 * t187 * t190 - 0.7e1 * t322 * t393 - 0.63e2 / 0.2e1 * t183 * t325 * t328 + 0.70e2 * t375 * t187 * t94 * t33 * t76 * t25 + 0.21e2 / 0.2e1 * t367 * t187 * t190;
745 t785 = 0.21e2 / 0.2e1 * t389 * t393 + 0.7e1 / 0.2e1 * t188 * t94 * t96 * t315 + 0.21e2 * t179 * t325 * t328 + 0.42e2 * t91 * t623 * t15 * t94 * t96 * t286 - 0.63e2 / 0.2e1 * t91 * t325 * params->d * t72 * t95 * t76 * t156 + 0.6e1 * t87 * t133 * t134 - 0.3e1 * t343 * t156 - 0.35e2 / 0.6e1 * t6 * t58 * t92 * t97 + 0.385e3 / 0.3e1 * t381 * t187 * t335 + 0.105e3 * t194 * t325 * t94 * t151 * t134 - 0.70e2 * t385 * t187 * t335 - 0.35e2 * t333 * t94 * t151 * t156 - 0.7e1 / 0.4e1 * t90 * t275 * t92 * t97 - 0.35e2 * t83 * t58 * t92 * t378 + 0.70e2 / 0.9e1 / t4 / t110 * t12 * t92 * t94 * t562 * sigma[0] + 0.140e3 / 0.27e2 * t46 * t12 * t92 * t97 + 0.7e1 / 0.2e1 * t178 * t127 * t92 * t97 + 0.105e3 * t53 * t12 * t92 * t378;
746 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
747 v4rho3sigma[0] = t710 + t785;
748
749 t869 = -0.35e2 / 0.12e2 * t578 * t12 * t63 * t196 * t563 - 0.28e2 / 0.3e1 * t414 * t187 * t190 + 0.14e2 * t418 * t187 * t190 + 0.7e1 * t422 * t393 - 0.35e2 / 0.6e1 * t333 * t453 + 0.35e2 / 0.4e1 * t447 * t187 * t453 + 0.35e2 / 0.8e1 * t451 * t94 * t229 * t156 - 0.35e2 / 0.2e1 * t458 * t187 * t94 * t460 * t76 + 0.7e1 / 0.3e1 * t322 * t469 - 0.7e1 / 0.2e1 * t389 * t469 - 0.7e1 / 0.4e1 * t188 * t94 * t234 * t156;
750 t943 = -0.147e3 / 0.8e1 * t438 * t439 * t26 * t156 + 0.35e2 / 0.12e2 * t385 * t92 * t230 - 0.175e3 / 0.8e1 * t376 * t461 - 0.35e2 / 0.16e2 * t225 * t127 * t92 * t230 + 0.35e2 / 0.4e1 * t457 * t58 * t92 * t461 - 0.7e1 / 0.6e1 * t372 * t235 + 0.7e1 / 0.8e1 * t368 * t235 - 0.35e2 / 0.9e1 * t6 * t87 * t92 * t97 + 0.147e3 / 0.2e1 * t215 * t623 * t217 * t439 * t26 * t134 - 0.175e3 / 0.72e2 * t381 * t92 * t230 + 0.35e2 / 0.36e2 * t319 * t235;
751 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
752 v4rho2sigma2[0] = -0.2e1 * t208 * t63 * t76 + rho[0] * (-0.608e3 / 0.9e1 * t349 * t201 + 0.1336e4 / 0.9e1 * t352 * t401 - 0.200e3 / 0.3e1 * params->b * t359 * t270 * t49 * t9 + 0.64e2 / 0.9e1 * t118 / t3 / t119 / t110 * t606 * t9) * t43 - 0.2e1 * t411 * t77 - t209 * t157 + 0.2e1 * t209 * t135 + 0.2205e4 / 0.4e1 * t260 * t12 * t133 * t222 - 0.385e3 / 0.9e1 * t28 * t87 * t63 * t197 + 0.70e2 / 0.3e1 * t193 * t173 * t63 * t197 - 0.441e3 / 0.4e1 * t121 * t58 * t133 * t222 + 0.49e2 / 0.8e1 * t168 * t127 * t133 * t222 + t869 - 0.21e2 * t211 * t325 * t328 - 0.105e3 / 0.8e1 * t226 * t325 * t94 * t229 * t134 + 0.21e2 / 0.4e1 * t326 * t94 * t234 * t134 + 0.2e1 * t410 * t43 + 0.14e2 / 0.3e1 * t178 * t173 * t92 * t97 - 0.7e1 / 0.2e1 * t90 * t364 * t92 * t97 - 0.140e3 / 0.3e1 * t425 * t187 * t335 - 0.70e2 / 0.3e1 * t83 * t87 * t92 * t378 + 0.1323e4 / 0.4e1 * t429 * t437 * t441 - 0.147e3 / 0.4e1 * t433 * t437 * t441 + t943;
753
754 t1016 = t481 * t43 - 0.1029e4 / 0.32e2 * t477 * t58 * t285 * t502 + 0.1715e4 / 0.16e2 / t267 * t12 * t285 * t500 * t26 * t151 + 0.343e3 * t405 * t12 * t285 * t502 + rho[0] * (-0.24e2 * t399 * t473 + 0.56e2 / 0.3e1 * params->b * t405 * t270 * sigma[0] * t9 - 0.8e1 / 0.3e1 * t50 / t3 / t593 * t606 * t9) * t43 - t482 * t77 - 0.1323e4 / 0.8e1 * t121 * t87 * t133 * t222 - 0.3969e4 / 0.32e2 * t430 * t506 + 0.35e2 / 0.2e1 * t193 * t208 * t63 * t197 + 0.147e3 / 0.8e1 * t168 * t173 * t133 * t222 + 0.1029e4 / 0.8e1 * t498 * t623 * t501 * t71 * t76 + 0.441e3 / 0.32e2 * t434 * t506 + 0.21e2 / 0.2e1 * t484 * t187 * t190 + 0.105e3 / 0.8e1 * t492 * t187 * t453 - 0.21e2 / 0.4e1 * t422 * t469 + 0.105e3 / 0.32e2 * t509 * t187 * t94 * t511 * t76 - 0.105e3 / 0.16e2 * t451 * t94 * t516 * t76;
755 t1071 = 0.21e2 / 0.8e1 * t188 * t94 * t522 * t76 - 0.21e2 / 0.4e1 * t90 * t410 * t92 * t97 - 0.441e3 / 0.8e1 * t488 * t437 * t441 - 0.105e3 / 0.16e2 * t225 * t173 * t92 * t230 + 0.105e3 / 0.8e1 * t457 * t87 * t92 * t461 + 0.21e2 / 0.8e1 * t419 * t235 - 0.1323e4 / 0.32e2 * t438 * t439 * t25 * t76 - 0.105e3 / 0.64e2 * t295 * t58 * t92 * t512 + 0.35e2 / 0.32e2 * t249 * t12 * t92 * t94 * t562 * t228 + 0.105e3 / 0.32e2 * t448 * t517 - 0.21e2 / 0.16e2 * t184 * t523 + 0.7e1 / 0.2e1 * t178 * t208 * t92 * t97 + 0.35e2 / 0.8e1 * t425 * t92 * t230 - 0.7e1 / 0.4e1 * t415 * t235 - 0.35e2 / 0.32e2 * t459 * t512 - 0.35e2 / 0.16e2 * t444 * t517 + 0.7e1 / 0.8e1 * t180 * t523;
756 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
757 v4rhosigma3[0] = t1016 + t1071;
758
759 t1118 = t218 * t218;
760 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
761 v4sigma4[0] = -0.7e1 * t90 * t481 * t92 * t97 - 0.105e3 / 0.8e1 * t225 * t208 * t92 * t230 + 0.441e3 / 0.8e1 * t489 * t506 - 0.105e3 / 0.16e2 * t295 * t87 * t92 * t512 + 0.735e3 / 0.128e3 * t216 * t219 * t220 * t95 - 0.105e3 / 0.256e3 * t106 * t12 * t92 * t94 * t562 * t515 - 0.1029e4 / 0.8e1 * t477 * t87 * t285 * t502 - 0.13377e5 / 0.128e3 * t498 * t285 * t500 * t71 - 0.5145e4 / 0.128e3 * t588 * t12 * t285 * t500 * t25 * t151 + 0.147e3 / 0.4e1 * t168 * t208 * t133 * t222 + 0.7203e4 / 0.32e2 / t602 / t44 * t12 * t623 * t1118 * t117 * t15 * t23 + 0.21e2 / 0.4e1 * t485 * t235 + 0.105e3 / 0.8e1 * t493 * t517 - 0.21e2 / 0.4e1 * t212 * t523 + 0.315e3 / 0.64e2 * t510 * t94 * t33 * t521 - 0.525e3 / 0.64e2 * t227 * t94 * t151 / t117 + 0.105e3 / 0.32e2 * t93 * t94 * t71 / t25 / t117 + rho[0] * (-0.4e1 * params->b * t477 * t270 * t9 + t1 / t3 / t267 * t606 * t9) * t43;
762
763 #ifndef XC_DONT_COMPILE_MXC
764
765 if(order < 5) return;
766
767
768 #endif
769
770 #endif
771
772 #endif
773
774 #endif
775
776 #endif
777
778
779 }
780
781
782 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 *))783 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
784 {
785
786 #ifndef XC_DONT_COMPILE_EXC
787 double t2, t3, t4, t5, t6, t7, t9, t10;
788 double t12, t15, t16, t18, t19, t20, t21, t25;
789 double t26, t28, t29, t30, t31, t34, t35, t36;
790 double t41, t45, t46;
791
792 #ifndef XC_DONT_COMPILE_VXC
793 double t47, t49, t52, t53, t54, t56, t61, t62;
794 double t64, t65, t66, t74, t75, t76, t79, t80;
795 double t84, t86, t90, t91, t93, t94, t95, t96;
796 double t97, t98, t99, t100, t101, t103, t104;
797
798 #ifndef XC_DONT_COMPILE_FXC
799 double t109, t113, t117, t119, t124, t125, t126, t127;
800 double t128, t129, t134, t135, t140, t141, t142, t155;
801 double t158, t159, t160, t163, t164, t168, t170, t171;
802 double t172, t174, t175, t178, t180, t181, t185, t186;
803 double t187, t188, t190, t191, t192, t194, t195, t197;
804 double t198, t200, t201, t203, t204, t205, t211, t212;
805 double t220, t221, t223, t226, t227, t228, t230, t231;
806 double t233, t234, t235, t236, t237, t238, t239, t241;
807 double t242, t244, t245, t246, t247, t248, t249, t250;
808 double t252, t253, t254, t255, t259, t260, t262, t263;
809 double t264, t271, t272;
810
811 #ifndef XC_DONT_COMPILE_KXC
812 double t280, t283, t289, t294, t299, t300, t305, t306;
813 double t307, t309, t310, t315, t316, t324, t325, t326;
814 double t327, t330, t331, t335, t352, t355, t356, t358;
815 double t359, t361, t362, t369, t370, t371, t373, t374;
816 double t375, t377, t379, t380, t382, t383, t385, t386;
817 double t388, t390, t392, t394, t396, t397, t398, t400;
818 double t401, t403, t404, t407, t408, t412, t416, t418;
819 double t419, t421, t423, t424, t426, t428, t431, t433;
820 double t434, t454, t461, t462, t465, t467, t469, t470;
821 double t473, t476, t478, t479, t482, t483, t485, t486;
822 double t487, t489, t490, t492, t493, t496, t497, t498;
823 double t499, t501, t502, t503, t505, t506, t507, t509;
824 double t510, t512, t513, t515, t516, t517, t519, t521;
825 double t522, t525, t526, t527, t528, t529, t530, t532;
826 double t534, t537, t538, t544, t545, t548, t549, t550;
827 double t552, t553, t554, t556, t557, t559, t561, t574;
828 double t579, t580, t598, t599, t602, t605, t606, t607;
829 double t609, t610, t611, t613, t614, t615, t617, t618;
830 double t619, t621, t623, t625, t626, t627, t628, t631;
831 double t632, t634, t635, t636, t637, t638, t640, t641;
832 double t642, t643, t646, t647, t648, t649, t653, t654;
833 double t656, t657, t658, t660, t661, t662, t664, t665;
834 double t666, t668, t671, t672, t673, t681, t682, t684;
835 double t685, t686, t699, t700;
836
837 #ifndef XC_DONT_COMPILE_LXC
838 double t747, t748, t765, t775, t780, t789, t793, t808;
839 double t809, t817, t825, t828, t831, t836, t842, t862;
840 double t865, t868, t874, t877, t879, t881, t885, t888;
841 double t894, t898, t900, t903, t907, t919, t923, t927;
842 double t930, t934, t938, t946, t950, t953, t957, t961;
843 double t968, t970, t1007, t1026, t1032, t1036, t1039, t1043;
844 double t1045, t1048, t1052, t1055, t1057, t1062, t1066, t1068;
845 double t1070, t1075, t1080, t1083, t1085, t1092, t1096, t1100;
846 double t1103, t1105, t1106, t1109, t1113, t1118, t1120, t1122;
847 double t1127, t1131, t1135, t1140, t1142, t1147, t1154, t1161;
848 double t1166, t1170, t1173, t1177, t1182, t1184, t1188, t1192;
849 double t1195, t1199, t1203, t1209, t1225, t1228, t1231, t1233;
850 double t1236, t1244, t1264, t1268, t1270, t1296, t1329, t1331;
851 double t1336, t1342, t1349, t1351, t1355, t1359, t1364, t1366;
852 double t1371, t1374, t1376, t1378, t1380, t1382, t1386, t1389;
853 double t1393, t1395, t1398, t1400, t1404, t1408, t1414, t1416;
854 double t1418, t1423, t1430, t1434, t1437, t1440, t1442, t1447;
855 double t1451, t1455, t1457, t1464, t1468, t1472, t1476, t1479;
856 double t1484, t1485, t1486, t1491, t1492, t1493, t1494, t1497;
857 double t1504, t1508, t1511, t1515, t1519, t1521, t1524, t1526;
858 double t1529, t1534, t1535, t1540, t1543, t1545, t1547, t1556;
859 double t1560, t1569, t1576, t1590, t1594, t1596, t1607, t1631;
860 double t1649, t1652, t1656, t1662, t1666, t1667, t1672, t1676;
861 double t1680, t1686, t1692, t1696, t1697, t1698, t1702, t1706;
862 double t1712, t1714, t1715, t1716, t1718, t1722, t1727, t1733;
863 double t1737, t1747, t1749, t1757, t1767, t1771, t1775, t1777;
864 double t1779, t1781, t1783, t1785, t1790, t1796, t1804, t1806;
865 double t1808, t1809, t1821, t1825, t1827, t1829, t1844, t1858;
866 double t1860;
867 #endif
868
869 #endif
870
871 #endif
872
873 #endif
874
875 #endif
876
877
878 gga_c_wi_params *params;
879
880 assert(p->params != NULL);
881 params = (gga_c_wi_params * )(p->params);
882
883 t2 = sigma[0] + 0.2e1 * sigma[1] + sigma[2];
884 t3 = params->b * t2;
885 t4 = rho[0] + rho[1];
886 t5 = t4 * t4;
887 t6 = POW_1_3(t4);
888 t7 = t6 * t6;
889 t9 = 0.1e1 / t7 / t5;
890 t10 = params->k * t2;
891 t12 = exp(-t10 * t9);
892 t15 = t3 * t9 * t12 + params->a;
893 t16 = M_CBRT3;
894 t18 = POW_1_3(0.1e1 / M_PI);
895 t19 = t16 * t18;
896 t20 = M_CBRT4;
897 t21 = t20 * t20;
898 t25 = t16 * t16;
899 t26 = M_CBRTPI;
900 t28 = sqrt(t2);
901 t29 = t28 * t2;
902 t30 = t5 * t5;
903 t31 = 0.1e1 / t30;
904 t34 = 0.1e1 / t6 / t4;
905 t35 = t28 * t34;
906 t36 = sqrt(t35);
907 t41 = 0.1e1 + params->d * t20 * t25 * t26 * t36 * t29 * t31 / 0.3e1;
908 t45 = params->c + t19 * t21 / t6 * t41 / 0.4e1;
909 t46 = 0.1e1 / t45;
910 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
911 zk[0] = t15 * t46;
912
913 #ifndef XC_DONT_COMPILE_VXC
914
915 if(order < 1) return;
916
917
918 t47 = t5 * t4;
919 t49 = 0.1e1 / t7 / t47;
920 t52 = t2 * t2;
921 t53 = params->b * t52;
922 t54 = t30 * t5;
923 t56 = 0.1e1 / t6 / t54;
924 t61 = 0.8e1 / 0.3e1 * t53 * t56 * params->k * t12 - 0.8e1 / 0.3e1 * t3 * t49 * t12;
925 t62 = t4 * t61;
926 t64 = t4 * t15;
927 t65 = t45 * t45;
928 t66 = 0.1e1 / t65;
929 t74 = t36 * t2 * t9;
930 t75 = t26 * t74;
931 t76 = t75 * t28;
932 t79 = -t19 * t21 * t34 * t41 / 0.12e2 - 0.14e2 / 0.3e1 * t18 * t9 * params->d * t76;
933 t80 = t66 * t79;
934 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
935 vrho[0] = t62 * t46 - t64 * t80 + (t15 * t46);
936
937 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
938 vrho[1] = vrho[0];
939
940 t84 = t30 * t4;
941 t86 = 0.1e1 / t6 / t84;
942 t90 = -t3 * t86 * params->k * t12 + params->b * t9 * t12;
943 t91 = t4 * t90;
944 t93 = 0.1e1 / t7;
945 t94 = t93 * t15;
946 t95 = t66 * t18;
947 t96 = t94 * t95;
948 t97 = params->d * t26;
949 t98 = 0.1e1 / t28;
950 t99 = t74 * t98;
951 t100 = t97 * t99;
952 t101 = t96 * t100;
953 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
954 vsigma[0] = t91 * t46 - 0.7e1 / 0.4e1 * t101;
955
956 t103 = 0.2e1 * t90;
957 t104 = t4 * t103;
958 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
959 vsigma[1] = t104 * t46 - 0.7e1 / 0.2e1 * t101;
960
961 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
962 vsigma[2] = vsigma[0];
963
964 #ifndef XC_DONT_COMPILE_FXC
965
966 if(order < 2) return;
967
968
969 t109 = t15 * t66;
970 t113 = 0.1e1 / t7 / t30;
971 t117 = t30 * t47;
972 t119 = 0.1e1 / t6 / t117;
973 t124 = t52 * t2;
974 t125 = params->b * t124;
975 t126 = t30 * t30;
976 t127 = t126 * t5;
977 t128 = 0.1e1 / t127;
978 t129 = params->k * params->k;
979 t134 = 0.88e2 / 0.9e1 * t3 * t113 * t12 - 0.24e2 * t53 * t119 * params->k * t12 + 0.64e2 / 0.9e1 * t125 * t128 * t129 * t12;
980 t135 = t4 * t134;
981 t140 = 0.1e1 / t65 / t45;
982 t141 = t79 * t79;
983 t142 = t140 * t141;
984 t155 = 0.1e1 / t84;
985 t158 = t36 * t35;
986 t159 = t26 * t158;
987 t160 = t159 * t2;
988 t163 = t19 * t21 / t6 / t5 * t41 / 0.9e1 + 0.14e2 * t18 * t49 * params->d * t76 + 0.140e3 / 0.9e1 * t18 * t155 * params->d * t160;
989 t164 = t66 * t163;
990 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
991 v2rho2[0] = -0.2e1 * t109 * t79 + t135 * t46 + 0.2e1 * t64 * t142 - t64 * t164 + 0.2e1 * t61 * t46 - 0.2e1 * t62 * t80;
992
993 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
994 v2rho2[1] = v2rho2[0];
995
996 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
997 v2rho2[2] = v2rho2[1];
998
999 t168 = params->b * t49 * t12;
1000 t170 = params->b * t56;
1001 t171 = t10 * t12;
1002 t172 = t170 * t171;
1003 t174 = t126 * t4;
1004 t175 = 0.1e1 / t174;
1005 t178 = t53 * t175 * t129 * t12;
1006 t180 = -0.8e1 / 0.3e1 * t168 + 0.8e1 * t172 - 0.8e1 / 0.3e1 * t178;
1007 t181 = t4 * t180;
1008 t185 = 0.1e1 / t7 / t4;
1009 t186 = t185 * t15;
1010 t187 = t186 * t95;
1011 t188 = t187 * t100;
1012 t190 = t93 * t61;
1013 t191 = t190 * t95;
1014 t192 = t191 * t100;
1015 t194 = t140 * t18;
1016 t195 = t94 * t194;
1017 t197 = t97 * t99 * t79;
1018 t198 = t195 * t197;
1019 t200 = 0.1e1 / t47;
1020 t201 = t200 * t15;
1021 t203 = t18 * params->d;
1022 t204 = t203 * t159;
1023 t205 = t201 * t66 * t204;
1024 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1025 v2rhosigma[0] = t90 * t46 + t181 * t46 - t91 * t80 + 0.7e1 / 0.6e1 * t188 - 0.7e1 / 0.4e1 * t192 + 0.7e1 / 0.2e1 * t198 + 0.35e2 / 0.6e1 * t205;
1026
1027 t211 = -0.16e2 / 0.3e1 * t168 + 0.16e2 * t172 - 0.16e2 / 0.3e1 * t178;
1028 t212 = t4 * t211;
1029 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1030 v2rhosigma[1] = t103 * t46 + t212 * t46 - t104 * t80 + 0.7e1 / 0.3e1 * t188 - 0.7e1 / 0.2e1 * t192 + 0.7e1 * t198 + 0.35e2 / 0.3e1 * t205;
1031
1032 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1033 v2rhosigma[2] = v2rhosigma[0];
1034
1035 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1036 v2rhosigma[3] = v2rhosigma[2];
1037
1038 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1039 v2rhosigma[4] = v2rhosigma[1];
1040
1041 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1042 v2rhosigma[5] = v2rhosigma[3];
1043
1044 t220 = params->k * t12;
1045 t221 = params->b * t86 * t220;
1046 t223 = 0.1e1 / t126;
1047 t226 = t3 * t223 * t129 * t12;
1048 t227 = -0.2e1 * t221 + t226;
1049 t228 = t4 * t227;
1050 t230 = t93 * t90;
1051 t231 = t230 * t95;
1052 t233 = 0.7e1 / 0.2e1 * t231 * t100;
1053 t234 = t175 * t15;
1054 t235 = t234 * t140;
1055 t236 = t18 * t18;
1056 t237 = params->d * params->d;
1057 t238 = t236 * t237;
1058 t239 = t26 * t26;
1059 t241 = t238 * t239 * t29;
1060 t242 = t235 * t241;
1061 t244 = 0.1e1 / t5;
1062 t245 = t244 * t15;
1063 t246 = t245 * t95;
1064 t247 = 0.1e1 / t2;
1065 t248 = t158 * t247;
1066 t249 = t97 * t248;
1067 t250 = t246 * t249;
1068 t252 = 0.1e1 / t29;
1069 t253 = t74 * t252;
1070 t254 = t97 * t253;
1071 t255 = t96 * t254;
1072 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1073 v2sigma2[0] = t228 * t46 - t233 + 0.49e2 / 0.8e1 * t242 - 0.35e2 / 0.16e2 * t250 + 0.7e1 / 0.8e1 * t255;
1074
1075 t259 = -0.4e1 * t221 + 0.2e1 * t226;
1076 t260 = t4 * t259;
1077 t262 = t93 * t103;
1078 t263 = t262 * t95;
1079 t264 = t263 * t100;
1080 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1081 v2sigma2[1] = t260 * t46 - 0.7e1 / 0.4e1 * t264 - t233 + 0.49e2 / 0.4e1 * t242 - 0.35e2 / 0.8e1 * t250 + 0.7e1 / 0.4e1 * t255;
1082
1083 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1084 v2sigma2[2] = v2sigma2[0];
1085
1086 t271 = -0.8e1 * t221 + 0.4e1 * t226;
1087 t272 = t4 * t271;
1088 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1089 v2sigma2[3] = t272 * t46 - 0.7e1 * t264 + 0.49e2 / 0.2e1 * t242 - 0.35e2 / 0.4e1 * t250 + 0.7e1 / 0.2e1 * t255;
1090
1091 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1092 v2sigma2[4] = v2sigma2[1];
1093
1094 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1095 v2sigma2[5] = v2sigma2[2];
1096
1097 #ifndef XC_DONT_COMPILE_KXC
1098
1099 if(order < 3) return;
1100
1101
1102 t280 = t61 * t66;
1103 t283 = t15 * t140;
1104 t289 = 0.1e1 / t7 / t84;
1105 t294 = 0.1e1 / t6 / t126;
1106 t299 = t126 * t47;
1107 t300 = 0.1e1 / t299;
1108 t305 = t52 * t52;
1109 t306 = params->b * t305;
1110 t307 = t126 * t84;
1111 t309 = 0.1e1 / t7 / t307;
1112 t310 = t129 * params->k;
1113 t315 = -0.1232e4 / 0.27e2 * t3 * t289 * t12 + 0.5456e4 / 0.27e2 * t53 * t294 * params->k * t12 - 0.1216e4 / 0.9e1 * t125 * t300 * t129 * t12 + 0.512e3 / 0.27e2 * t306 * t309 * t310 * t12;
1114 t316 = t4 * t315;
1115 t324 = t65 * t65;
1116 t325 = 0.1e1 / t324;
1117 t326 = t141 * t79;
1118 t327 = t325 * t326;
1119 t330 = t140 * t79;
1120 t331 = t330 * t163;
1121 t335 = 0.1e1 / t6 / t47;
1122 t352 = t26 * t36 * t29;
1123 t355 = -0.7e1 / 0.27e2 * t19 * t21 * t335 * t41 - 0.1442e4 / 0.27e2 * t18 * t113 * params->d * t76 - 0.1120e4 / 0.9e1 * t18 / t54 * params->d * t160 - 0.280e3 / 0.9e1 * t18 * t119 * params->d * t352;
1124 t356 = t66 * t355;
1125 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1126 v3rho3[0] = -0.3e1 * t109 * t163 + 0.3e1 * t134 * t46 - 0.3e1 * t135 * t80 + 0.6e1 * t283 * t141 + 0.6e1 * t62 * t142 - 0.3e1 * t62 * t164 - 0.6e1 * t280 * t79 + t316 * t46 - 0.6e1 * t64 * t327 + 0.6e1 * t64 * t331 - t64 * t356;
1127
1128 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1129 v3rho3[1] = v3rho3[0];
1130
1131 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1132 v3rho3[2] = v3rho3[1];
1133
1134 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1135 v3rho3[3] = v3rho3[2];
1136
1137 t358 = t325 * t18;
1138 t359 = t94 * t358;
1139 t361 = t97 * t99 * t141;
1140 t362 = t359 * t361;
1141 t369 = t185 * t61;
1142 t370 = t369 * t95;
1143 t371 = t370 * t100;
1144 t373 = t93 * t134;
1145 t374 = t373 * t95;
1146 t375 = t374 * t100;
1147 t377 = t201 * t194;
1148 t379 = t97 * t158 * t79;
1149 t380 = t377 * t379;
1150 t382 = t86 * t15;
1151 t383 = t382 * t95;
1152 t385 = t97 * t36 * t28;
1153 t386 = t383 * t385;
1154 t388 = t31 * t15;
1155 t390 = t388 * t66 * t204;
1156 t392 = t200 * t61;
1157 t394 = t392 * t66 * t204;
1158 t396 = t9 * t15;
1159 t397 = t396 * t95;
1160 t398 = t397 * t100;
1161 t400 = t186 * t194;
1162 t401 = t400 * t197;
1163 t403 = t190 * t194;
1164 t404 = t403 * t197;
1165 t407 = t97 * t99 * t163;
1166 t408 = t195 * t407;
1167 t412 = t90 * t66;
1168 t416 = params->b * t113 * t12;
1169 t418 = params->b * t119;
1170 t419 = t418 * t171;
1171 t421 = params->b * t128;
1172 t423 = t129 * t52 * t12;
1173 t424 = t421 * t423;
1174 t426 = t126 * t30;
1175 t428 = 0.1e1 / t7 / t426;
1176 t431 = t125 * t428 * t310 * t12;
1177 t433 = 0.88e2 / 0.9e1 * t416 - 0.520e3 / 0.9e1 * t419 + 0.136e3 / 0.3e1 * t424 - 0.64e2 / 0.9e1 * t431;
1178 t434 = t4 * t433;
1179 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1180 v3rho2sigma[0] = -0.21e2 / 0.2e1 * t362 + 0.2e1 * t91 * t142 - 0.2e1 * t181 * t80 - t91 * t164 + 0.7e1 / 0.3e1 * t371 - 0.7e1 / 0.4e1 * t375 - 0.70e2 / 0.3e1 * t380 - 0.35e2 / 0.3e1 * t386 - 0.385e3 / 0.18e2 * t390 + 0.35e2 / 0.3e1 * t394 - 0.35e2 / 0.18e2 * t398 - 0.14e2 / 0.3e1 * t401 + 0.7e1 * t404 + 0.7e1 / 0.2e1 * t408 + 0.2e1 * t180 * t46 - 0.2e1 * t412 * t79 + t434 * t46;
1181
1182 t454 = t103 * t66;
1183 t461 = 0.176e3 / 0.9e1 * t416 - 0.1040e4 / 0.9e1 * t419 + 0.272e3 / 0.3e1 * t424 - 0.128e3 / 0.9e1 * t431;
1184 t462 = t4 * t461;
1185 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1186 v3rho2sigma[1] = -0.21e2 * t362 + 0.2e1 * t104 * t142 - 0.2e1 * t212 * t80 - t104 * t164 + 0.14e2 / 0.3e1 * t371 - 0.7e1 / 0.2e1 * t375 - 0.140e3 / 0.3e1 * t380 - 0.70e2 / 0.3e1 * t386 - 0.385e3 / 0.9e1 * t390 + 0.70e2 / 0.3e1 * t394 - 0.35e2 / 0.9e1 * t398 - 0.28e2 / 0.3e1 * t401 + 0.14e2 * t404 + 0.7e1 * t408 + 0.2e1 * t211 * t46 - 0.2e1 * t454 * t79 + t462 * t46;
1187
1188 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1189 v3rho2sigma[2] = v3rho2sigma[0];
1190
1191 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1192 v3rho2sigma[3] = v3rho2sigma[2];
1193
1194 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1195 v3rho2sigma[4] = v3rho2sigma[1];
1196
1197 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1198 v3rho2sigma[5] = v3rho2sigma[3];
1199
1200 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1201 v3rho2sigma[6] = v3rho2sigma[5];
1202
1203 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1204 v3rho2sigma[7] = v3rho2sigma[4];
1205
1206 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1207 v3rho2sigma[8] = v3rho2sigma[6];
1208
1209 t465 = t170 * t220;
1210 t467 = params->b * t175;
1211 t469 = t129 * t2 * t12;
1212 t470 = t467 * t469;
1213 t473 = 0.1e1 / t7 / t299;
1214 t476 = t53 * t473 * t310 * t12;
1215 t478 = 0.32e2 / 0.3e1 * t465 - 0.40e2 / 0.3e1 * t470 + 0.8e1 / 0.3e1 * t476;
1216 t479 = t4 * t478;
1217 t482 = t185 * t90;
1218 t483 = t482 * t95;
1219 t485 = 0.7e1 / 0.3e1 * t483 * t100;
1220 t486 = t93 * t180;
1221 t487 = t486 * t95;
1222 t489 = 0.7e1 / 0.2e1 * t487 * t100;
1223 t490 = t230 * t194;
1224 t492 = 0.7e1 * t490 * t197;
1225 t493 = t200 * t90;
1226 t496 = 0.35e2 / 0.3e1 * t493 * t66 * t204;
1227 t497 = t128 * t15;
1228 t498 = t497 * t140;
1229 t499 = t498 * t241;
1230 t501 = t175 * t61;
1231 t502 = t501 * t140;
1232 t503 = t502 * t241;
1233 t505 = t325 * t236;
1234 t506 = t234 * t505;
1235 t507 = t237 * t239;
1236 t509 = t507 * t29 * t79;
1237 t510 = t506 * t509;
1238 t512 = t201 * t95;
1239 t513 = t512 * t249;
1240 t515 = t244 * t61;
1241 t516 = t515 * t95;
1242 t517 = t516 * t249;
1243 t519 = t245 * t194;
1244 t521 = t97 * t248 * t79;
1245 t522 = t519 * t521;
1246 t525 = 0.1e1 / t6 / t30;
1247 t526 = t525 * t15;
1248 t527 = t526 * t95;
1249 t528 = t36 * t98;
1250 t529 = t97 * t528;
1251 t530 = t527 * t529;
1252 t532 = t187 * t254;
1253 t534 = t191 * t254;
1254 t537 = t97 * t253 * t79;
1255 t538 = t195 * t537;
1256 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1257 v3rhosigma2[0] = t227 * t46 + t479 * t46 - t228 * t80 + t485 - t489 + t492 + t496 - 0.441e3 / 0.8e1 * t499 + 0.49e2 / 0.8e1 * t503 - 0.147e3 / 0.8e1 * t510 + 0.35e2 / 0.24e2 * t513 - 0.35e2 / 0.16e2 * t517 + 0.35e2 / 0.8e1 * t522 + 0.35e2 / 0.8e1 * t530 - 0.7e1 / 0.12e2 * t532 + 0.7e1 / 0.8e1 * t534 - 0.7e1 / 0.4e1 * t538;
1258
1259 t544 = 0.64e2 / 0.3e1 * t465 - 0.80e2 / 0.3e1 * t470 + 0.16e2 / 0.3e1 * t476;
1260 t545 = t4 * t544;
1261 t548 = t185 * t103;
1262 t549 = t548 * t95;
1263 t550 = t549 * t100;
1264 t552 = t93 * t211;
1265 t553 = t552 * t95;
1266 t554 = t553 * t100;
1267 t556 = t262 * t194;
1268 t557 = t556 * t197;
1269 t559 = t200 * t103;
1270 t561 = t559 * t66 * t204;
1271 t574 = t496 - 0.441e3 / 0.4e1 * t499 + 0.49e2 / 0.4e1 * t503 - 0.147e3 / 0.4e1 * t510 + 0.35e2 / 0.12e2 * t513 - 0.35e2 / 0.8e1 * t517 + 0.35e2 / 0.4e1 * t522 + 0.35e2 / 0.4e1 * t530 - 0.7e1 / 0.6e1 * t532 + 0.7e1 / 0.4e1 * t534 - 0.7e1 / 0.2e1 * t538;
1272 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1273 v3rhosigma2[1] = t259 * t46 + t545 * t46 - t260 * t80 + 0.7e1 / 0.6e1 * t550 - 0.7e1 / 0.4e1 * t554 + 0.7e1 / 0.2e1 * t557 + 0.35e2 / 0.6e1 * t561 + t485 - t489 + t492 + t574;
1274
1275 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1276 v3rhosigma2[2] = v3rhosigma2[0];
1277
1278 t579 = 0.128e3 / 0.3e1 * t465 - 0.160e3 / 0.3e1 * t470 + 0.32e2 / 0.3e1 * t476;
1279 t580 = t4 * t579;
1280 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1281 v3rhosigma2[3] = t271 * t46 + t580 * t46 - t272 * t80 + 0.14e2 / 0.3e1 * t550 - 0.7e1 * t554 + 0.14e2 * t557 + 0.70e2 / 0.3e1 * t561 - 0.441e3 / 0.2e1 * t499 + 0.49e2 / 0.2e1 * t503 - 0.147e3 / 0.2e1 * t510 + 0.35e2 / 0.6e1 * t513 - 0.35e2 / 0.4e1 * t517 + 0.35e2 / 0.2e1 * t522 + 0.35e2 / 0.2e1 * t530 - 0.7e1 / 0.3e1 * t532 + 0.7e1 / 0.2e1 * t534 - 0.7e1 * t538;
1282
1283 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1284 v3rhosigma2[4] = v3rhosigma2[1];
1285
1286 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1287 v3rhosigma2[5] = v3rhosigma2[2];
1288
1289 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1290 v3rhosigma2[6] = v3rhosigma2[5];
1291
1292 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1293 v3rhosigma2[7] = v3rhosigma2[4];
1294
1295 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1296 v3rhosigma2[8] = v3rhosigma2[6];
1297
1298 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1299 v3rhosigma2[9] = v3rhosigma2[3];
1300
1301 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1302 v3rhosigma2[10] = v3rhosigma2[7];
1303
1304 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1305 v3rhosigma2[11] = v3rhosigma2[8];
1306
1307 t598 = t129 * t12;
1308 t599 = params->b * t223 * t598;
1309 t602 = 0.1e1 / t7 / t127;
1310 t605 = t3 * t602 * t310 * t12;
1311 t606 = 0.3e1 * t599 - t605;
1312 t607 = t4 * t606;
1313 t609 = t93 * t227;
1314 t610 = t609 * t95;
1315 t611 = t610 * t100;
1316 t613 = t175 * t90;
1317 t614 = t613 * t140;
1318 t615 = t614 * t241;
1319 t617 = t244 * t90;
1320 t618 = t617 * t95;
1321 t619 = t618 * t249;
1322 t621 = t231 * t254;
1323 t623 = t602 * t15;
1324 t625 = t237 * params->d;
1325 t626 = t625 * t2;
1326 t627 = t626 * t74;
1327 t628 = t623 * t325 * t627;
1328 t631 = t238 * t239 * t28;
1329 t632 = t235 * t631;
1330 t634 = t335 * t15;
1331 t635 = t634 * t95;
1332 t636 = t36 * t252;
1333 t637 = t97 * t636;
1334 t638 = t635 * t637;
1335 t640 = 0.1e1 / t52;
1336 t641 = t158 * t640;
1337 t642 = t97 * t641;
1338 t643 = t246 * t642;
1339 t646 = 0.1e1 / t28 / t52;
1340 t647 = t74 * t646;
1341 t648 = t97 * t647;
1342 t649 = t96 * t648;
1343 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1344 v3sigma3[0] = t607 * t46 - 0.21e2 / 0.4e1 * t611 + 0.147e3 / 0.8e1 * t615 - 0.105e3 / 0.16e2 * t619 + 0.21e2 / 0.8e1 * t621 - 0.1029e4 / 0.32e2 * t628 + 0.441e3 / 0.32e2 * t632 - 0.105e3 / 0.64e2 * t638 + 0.105e3 / 0.32e2 * t643 - 0.21e2 / 0.16e2 * t649;
1345
1346 t653 = 0.6e1 * t599 - 0.2e1 * t605;
1347 t654 = t4 * t653;
1348 t656 = t93 * t259;
1349 t657 = t656 * t95;
1350 t658 = t657 * t100;
1351 t660 = t175 * t103;
1352 t661 = t660 * t140;
1353 t662 = t661 * t241;
1354 t664 = t244 * t103;
1355 t665 = t664 * t95;
1356 t666 = t665 * t249;
1357 t668 = t263 * t254;
1358 t671 = 0.49e2 / 0.2e1 * t615;
1359 t672 = 0.35e2 / 0.4e1 * t619;
1360 t673 = 0.7e1 / 0.2e1 * t621;
1361 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1362 v3sigma3[1] = t654 * t46 - 0.7e1 / 0.2e1 * t658 + 0.49e2 / 0.8e1 * t662 - 0.35e2 / 0.16e2 * t666 + 0.7e1 / 0.8e1 * t668 - 0.7e1 / 0.2e1 * t611 + t671 - t672 + t673 - 0.1029e4 / 0.16e2 * t628 + 0.441e3 / 0.16e2 * t632 - 0.105e3 / 0.32e2 * t638 + 0.105e3 / 0.16e2 * t643 - 0.21e2 / 0.8e1 * t649;
1363
1364 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1365 v3sigma3[2] = v3sigma3[0];
1366
1367 t681 = 0.12e2 * t599 - 0.4e1 * t605;
1368 t682 = t4 * t681;
1369 t684 = t93 * t271;
1370 t685 = t684 * t95;
1371 t686 = t685 * t100;
1372 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1373 v3sigma3[3] = t682 * t46 - 0.7e1 / 0.4e1 * t686 - 0.7e1 * t658 + 0.49e2 / 0.2e1 * t662 - 0.35e2 / 0.4e1 * t666 + 0.7e1 / 0.2e1 * t668 + t671 - 0.1029e4 / 0.8e1 * t628 + 0.441e3 / 0.8e1 * t632 - t672 - 0.105e3 / 0.16e2 * t638 + 0.105e3 / 0.8e1 * t643 + t673 - 0.21e2 / 0.4e1 * t649;
1374
1375 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1376 v3sigma3[4] = v3sigma3[1];
1377
1378 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1379 v3sigma3[5] = v3sigma3[2];
1380
1381 t699 = 0.24e2 * t599 - 0.8e1 * t605;
1382 t700 = t4 * t699;
1383 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1384 v3sigma3[6] = t700 * t46 - 0.21e2 / 0.2e1 * t686 + 0.147e3 / 0.2e1 * t662 - 0.105e3 / 0.4e1 * t666 + 0.21e2 / 0.2e1 * t668 - 0.1029e4 / 0.4e1 * t628 + 0.441e3 / 0.4e1 * t632 - 0.105e3 / 0.8e1 * t638 + 0.105e3 / 0.4e1 * t643 - 0.21e2 / 0.2e1 * t649;
1385
1386 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1387 v3sigma3[7] = v3sigma3[3];
1388
1389 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1390 v3sigma3[8] = v3sigma3[4];
1391
1392 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1393 v3sigma3[9] = v3sigma3[5];
1394
1395 #ifndef XC_DONT_COMPILE_LXC
1396
1397 if(order < 4) return;
1398
1399
1400 t747 = 0.1e1 / t36;
1401 t748 = t26 * t747;
1402 t765 = 0.1e1 / t7 / t54;
1403 t775 = 0.1e1 / t426;
1404 t780 = t126 * t54;
1405 t789 = t126 * t126;
1406 t793 = t129 * t129;
1407 t808 = 0.1e1 / t324 / t45;
1408 t809 = t141 * t141;
1409 t817 = t163 * t163;
1410 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1411 v4rho4[0] = 0.24e2 * t61 * t140 * t141 - 0.24e2 * t15 * t325 * t326 + 0.24e2 * t283 * t79 * t163 - 0.4e1 * t316 * t80 - 0.6e1 * t135 * t164 - 0.4e1 * t62 * t356 - t64 * t66 * (0.70e2 / 0.81e2 * t19 * t21 * t525 * t41 + 0.6860e4 / 0.27e2 * t18 * t289 * params->d * t76 + 0.74900e5 / 0.81e2 * t18 / t117 * params->d * t160 + 0.12880e5 / 0.27e2 * t18 * t294 * params->d * t352 + 0.560e3 / 0.27e2 * t18 / t7 / t174 * params->d * t748 * t52) + 0.4e1 * t315 * t46 - 0.12e2 * t134 * t66 * t79 - 0.12e2 * t280 * t163 - 0.4e1 * t109 * t355 + t4 * (0.20944e5 / 0.81e2 * t3 * t765 * t12 - 0.48752e5 / 0.27e2 * t53 / t6 / t174 * params->k * t12 + 0.164032e6 / 0.81e2 * t125 * t775 * t129 * t12 - 0.50176e5 / 0.81e2 * t306 / t7 / t780 * t310 * t12 + 0.4096e4 / 0.81e2 * params->b * t305 * t2 / t6 / t789 / t4 * t793 * t12) * t46 + 0.12e2 * t135 * t142 - 0.24e2 * t62 * t327 + 0.24e2 * t62 * t331 + 0.24e2 * t64 * t808 * t809 - 0.36e2 * t64 * t325 * t141 * t163 + 0.6e1 * t64 * t140 * t817 + 0.8e1 * t64 * t330 * t355;
1412
1413 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1414 v4rho4[1] = v4rho4[0];
1415
1416 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1417 v4rho4[2] = v4rho4[1];
1418
1419 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1420 v4rho4[3] = v4rho4[2];
1421
1422 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1423 v4rho4[4] = v4rho4[3];
1424
1425 t825 = params->b * t289 * t12;
1426 t828 = params->b * t294 * t171;
1427 t831 = params->b * t300 * t423;
1428 t836 = params->b * t309 * t310 * t124 * t12;
1429 t842 = t306 / t6 / t789 * t793 * t12;
1430 t862 = t190 * t358 * t361;
1431 t865 = t369 * t194 * t197;
1432 t868 = t373 * t194 * t197;
1433 t874 = t382 * t194 * t97 * t36 * t79 * t28;
1434 t877 = t396 * t194 * t197;
1435 t879 = t400 * t407;
1436 t881 = t403 * t407;
1437 t885 = t195 * t97 * t99 * t355;
1438 t888 = t186 * t358 * t361;
1439 t894 = t94 * t808 * t18 * t97 * t99 * t326;
1440 t898 = t155 * t15 * t66 * t204;
1441 t900 = t4 * (-0.1232e4 / 0.27e2 * t825 + 0.4048e4 / 0.9e1 * t828 - 0.16400e5 / 0.27e2 * t831 + 0.5696e4 / 0.27e2 * t836 - 0.512e3 / 0.27e2 * t842) * t46 - 0.6e1 * t91 * t327 + 0.3e1 * t433 * t46 + 0.6e1 * t91 * t331 - 0.6e1 * t180 * t66 * t79 - 0.3e1 * t412 * t163 + 0.6e1 * t90 * t140 * t141 - 0.63e2 / 0.2e1 * t862 - 0.14e2 * t865 + 0.21e2 / 0.2e1 * t868 + 0.70e2 * t874 + 0.35e2 / 0.3e1 * t877 - 0.7e1 * t879 + 0.21e2 / 0.2e1 * t881 + 0.7e1 / 0.2e1 * t885 + 0.21e2 * t888 + 0.42e2 * t894 + 0.2485e4 / 0.27e2 * t898;
1442 t903 = t31 * t61 * t66 * t204;
1443 t907 = t200 * t134 * t66 * t204;
1444 t919 = t201 * t358 * t97 * t158 * t141;
1445 t923 = t185 * t134 * t95 * t100;
1446 t927 = t93 * t315 * t95 * t100;
1447 t930 = t392 * t194 * t379;
1448 t934 = t377 * t97 * t158 * t163;
1449 t938 = t86 * t61 * t95 * t385;
1450 t946 = 0.1e1 / t7 / t117 * t15 * t95 * t97 * t747 * t2;
1451 t950 = t9 * t61 * t95 * t100;
1452 t953 = t388 * t194 * t379;
1453 t957 = t56 * t15 * t95 * t385;
1454 t961 = t49 * t15 * t95 * t100;
1455 t968 = t94 * t358 * params->d * t75 * t98 * t79 * t163;
1456 t970 = -0.385e3 / 0.6e1 * t903 + 0.35e2 / 0.2e1 * t907 + 0.6e1 * t181 * t142 - 0.3e1 * t434 * t80 - 0.3e1 * t181 * t164 - t91 * t356 + 0.105e3 * t919 + 0.7e1 / 0.2e1 * t923 - 0.7e1 / 0.4e1 * t927 - 0.70e2 * t930 - 0.35e2 * t934 - 0.35e2 * t938 + 0.70e2 / 0.9e1 * t946 - 0.35e2 / 0.6e1 * t950 + 0.385e3 / 0.3e1 * t953 + 0.105e3 * t957 + 0.140e3 / 0.27e2 * t961 - 0.63e2 / 0.2e1 * t968;
1457 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1458 v4rho3sigma[0] = t900 + t970;
1459
1460 t1007 = 0.3e1 * t461 * t46 + t4 * (-0.2464e4 / 0.27e2 * t825 + 0.8096e4 / 0.9e1 * t828 - 0.32800e5 / 0.27e2 * t831 + 0.11392e5 / 0.27e2 * t836 - 0.1024e4 / 0.27e2 * t842) * t46 - 0.3e1 * t454 * t163 + 0.6e1 * t103 * t140 * t141 - 0.6e1 * t211 * t66 * t79 - 0.6e1 * t104 * t327 + 0.6e1 * t212 * t142 - 0.3e1 * t462 * t80 - 0.3e1 * t212 * t164 - t104 * t356 + 0.6e1 * t104 * t331 - 0.63e2 * t862 - 0.28e2 * t865 + 0.21e2 * t868 + 0.140e3 * t874 + 0.70e2 / 0.3e1 * t877 - 0.14e2 * t879 + 0.21e2 * t881;
1461 t1026 = 0.7e1 * t885 + 0.42e2 * t888 + 0.84e2 * t894 + 0.4970e4 / 0.27e2 * t898 - 0.385e3 / 0.3e1 * t903 + 0.35e2 * t907 + 0.210e3 * t919 + 0.7e1 * t923 - 0.7e1 / 0.2e1 * t927 - 0.140e3 * t930 - 0.70e2 * t934 - 0.70e2 * t938 + 0.140e3 / 0.9e1 * t946 - 0.35e2 / 0.3e1 * t950 + 0.770e3 / 0.3e1 * t953 + 0.210e3 * t957 + 0.280e3 / 0.27e2 * t961 - 0.63e2 * t968;
1462 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1463 v4rho3sigma[1] = t1007 + t1026;
1464
1465 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1466 v4rho3sigma[2] = v4rho3sigma[0];
1467
1468 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1469 v4rho3sigma[3] = v4rho3sigma[2];
1470
1471 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1472 v4rho3sigma[4] = v4rho3sigma[1];
1473
1474 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1475 v4rho3sigma[5] = v4rho3sigma[3];
1476
1477 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1478 v4rho3sigma[6] = v4rho3sigma[5];
1479
1480 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1481 v4rho3sigma[7] = v4rho3sigma[4];
1482
1483 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1484 v4rho3sigma[8] = v4rho3sigma[6];
1485
1486 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1487 v4rho3sigma[9] = v4rho3sigma[8];
1488
1489 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1490 v4rho3sigma[10] = v4rho3sigma[7];
1491
1492 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1493 v4rho3sigma[11] = v4rho3sigma[9];
1494
1495 t1032 = 0.14e2 / 0.3e1 * t185 * t180 * t95 * t100;
1496 t1036 = 0.7e1 / 0.2e1 * t93 * t433 * t95 * t100;
1497 t1039 = 0.140e3 / 0.3e1 * t493 * t194 * t379;
1498 t1043 = 0.70e2 / 0.3e1 * t86 * t90 * t95 * t385;
1499 t1045 = t497 * t505 * t509;
1500 t1048 = t501 * t505 * t509;
1501 t1052 = t506 * t507 * t29 * t163;
1502 t1055 = t392 * t95 * t249;
1503 t1057 = t383 * t529;
1504 t1062 = t244 * t134 * t95 * t249;
1505 t1066 = t525 * t61 * t95 * t529;
1506 t1068 = t370 * t254;
1507 t1070 = t374 * t254;
1508 t1075 = 0.35e2 / 0.9e1 * t9 * t90 * t95 * t100;
1509 t1080 = t234 * t808 * t236 * t507 * t29 * t141;
1510 t1083 = t388 * t95 * t249;
1511 t1085 = t397 * t254;
1512 t1092 = -0.35e2 / 0.16e2 * t1062 + 0.35e2 / 0.4e1 * t1066 - 0.7e1 / 0.6e1 * t1068 + 0.7e1 / 0.8e1 * t1070 - t1075 + 0.147e3 / 0.2e1 * t1080 - 0.175e3 / 0.72e2 * t1083 + 0.35e2 / 0.36e2 * t1085 + 0.2e1 * t228 * t142 - 0.2e1 * t479 * t80 - t228 * t164;
1513 t1096 = t359 * t97 * t253 * t141;
1514 t1100 = 0.28e2 / 0.3e1 * t482 * t194 * t197;
1515 t1103 = 0.14e2 * t486 * t194 * t197;
1516 t1105 = 0.7e1 * t490 * t407;
1517 t1106 = t377 * t521;
1518 t1109 = t515 * t194 * t521;
1519 t1113 = t519 * t97 * t248 * t163;
1520 t1118 = t526 * t194 * t97 * t528 * t79;
1521 t1120 = t400 * t537;
1522 t1122 = t403 * t537;
1523 t1127 = t195 * t97 * t253 * t163;
1524 t1131 = 0.21e2 * t230 * t358 * t361;
1525 t1135 = t245 * t358 * t97 * t248 * t141;
1526 t1140 = t418 * t220;
1527 t1142 = t421 * t469;
1528 t1147 = params->b * t428 * t310 * t52 * t12;
1529 t1154 = t125 / t6 / t126 / t117 * t793 * t12;
1530 t1161 = t300 * t15 * t140 * t241;
1531 t1166 = 0.385e3 / 0.9e1 * t31 * t90 * t66 * t204;
1532 t1170 = 0.70e2 / 0.3e1 * t200 * t180 * t66 * t204;
1533 t1173 = t128 * t61 * t140 * t241;
1534 t1177 = t175 * t134 * t140 * t241;
1535 t1182 = t765 * t15 * t66 * t203 * t748;
1536 t1184 = -0.7e1 / 0.4e1 * t1127 - t1131 - 0.105e3 / 0.8e1 * t1135 - 0.2e1 * t227 * t66 * t79 + t4 * (-0.608e3 / 0.9e1 * t1140 + 0.1336e4 / 0.9e1 * t1142 - 0.200e3 / 0.3e1 * t1147 + 0.64e2 / 0.9e1 * t1154) * t46 + 0.2205e4 / 0.4e1 * t1161 - t1166 + t1170 - 0.441e3 / 0.4e1 * t1173 + 0.49e2 / 0.8e1 * t1177 - 0.35e2 / 0.12e2 * t1182;
1537 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1538 v4rho2sigma2[0] = 0.2e1 * t478 * t46 + t1032 - t1036 - t1039 - t1043 + 0.1323e4 / 0.4e1 * t1045 - 0.147e3 / 0.4e1 * t1048 - 0.147e3 / 0.8e1 * t1052 + 0.35e2 / 0.12e2 * t1055 - 0.175e3 / 0.8e1 * t1057 + t1092 + 0.21e2 / 0.4e1 * t1096 - t1100 + t1103 + t1105 - 0.35e2 / 0.6e1 * t1106 + 0.35e2 / 0.4e1 * t1109 + 0.35e2 / 0.8e1 * t1113 - 0.35e2 / 0.2e1 * t1118 + 0.7e1 / 0.3e1 * t1120 - 0.7e1 / 0.2e1 * t1122 + t1184;
1539
1540 t1188 = t185 * t211 * t95 * t100;
1541 t1192 = t93 * t461 * t95 * t100;
1542 t1195 = t559 * t194 * t379;
1543 t1199 = t86 * t103 * t95 * t385;
1544 t1203 = t9 * t103 * t95 * t100;
1545 t1209 = 0.7e1 / 0.3e1 * t1188 - 0.7e1 / 0.4e1 * t1192 - 0.70e2 / 0.3e1 * t1195 - 0.35e2 / 0.3e1 * t1199 - 0.35e2 / 0.18e2 * t1203 + t1032 - t1036 - t1039 - t1043 + 0.1323e4 / 0.2e1 * t1045 - 0.147e3 / 0.2e1 * t1048 - 0.147e3 / 0.4e1 * t1052 + 0.35e2 / 0.6e1 * t1055;
1546 t1225 = -0.175e3 / 0.4e1 * t1057 - 0.35e2 / 0.8e1 * t1062 + 0.35e2 / 0.2e1 * t1066 - 0.7e1 / 0.3e1 * t1068 + 0.7e1 / 0.4e1 * t1070 - t1075 + 0.147e3 * t1080 - 0.175e3 / 0.36e2 * t1083 + 0.35e2 / 0.18e2 * t1085 + 0.2e1 * t260 * t142 - 0.2e1 * t545 * t80 - t260 * t164 + 0.2e1 * t544 * t46;
1547 t1228 = t548 * t194 * t197;
1548 t1231 = t552 * t194 * t197;
1549 t1233 = t556 * t407;
1550 t1236 = t262 * t358 * t361;
1551 t1244 = -0.14e2 / 0.3e1 * t1228 + 0.7e1 * t1231 + 0.7e1 / 0.2e1 * t1233 - 0.21e2 / 0.2e1 * t1236 + 0.21e2 / 0.2e1 * t1096 - t1100 + t1103 + t1105 - 0.35e2 / 0.3e1 * t1106 + 0.35e2 / 0.2e1 * t1109 + 0.35e2 / 0.4e1 * t1113 - 0.35e2 * t1118 + 0.14e2 / 0.3e1 * t1120;
1552 t1264 = t31 * t103 * t66 * t204;
1553 t1268 = t200 * t211 * t66 * t204;
1554 t1270 = -0.7e1 * t1122 - 0.7e1 / 0.2e1 * t1127 - t1131 - 0.105e3 / 0.4e1 * t1135 - 0.2e1 * t259 * t66 * t79 + t4 * (-0.1216e4 / 0.9e1 * t1140 + 0.2672e4 / 0.9e1 * t1142 - 0.400e3 / 0.3e1 * t1147 + 0.128e3 / 0.9e1 * t1154) * t46 + 0.2205e4 / 0.2e1 * t1161 - t1166 + t1170 - 0.441e3 / 0.2e1 * t1173 + 0.49e2 / 0.4e1 * t1177 - 0.35e2 / 0.6e1 * t1182 - 0.385e3 / 0.18e2 * t1264 + 0.35e2 / 0.3e1 * t1268;
1555 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556 v4rho2sigma2[1] = t1209 + t1225 + t1244 + t1270;
1557
1558 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559 v4rho2sigma2[2] = v4rho2sigma2[0];
1560
1561 t1296 = -0.35e2 / 0.4e1 * t1062 + 0.35e2 * t1066 - 0.14e2 / 0.3e1 * t1068 + 0.7e1 / 0.2e1 * t1070 + 0.294e3 * t1080 - 0.175e3 / 0.18e2 * t1083 + 0.35e2 / 0.9e1 * t1085 + 0.2e1 * t272 * t142 - 0.2e1 * t580 * t80 - t272 * t164 - 0.56e2 / 0.3e1 * t1228;
1562 t1329 = -0.7e1 * t1127 - 0.105e3 / 0.2e1 * t1135 + 0.2e1 * t579 * t46 - 0.2e1 * t271 * t66 * t79 + t4 * (-0.2432e4 / 0.9e1 * t1140 + 0.5344e4 / 0.9e1 * t1142 - 0.800e3 / 0.3e1 * t1147 + 0.256e3 / 0.9e1 * t1154) * t46 + 0.2205e4 * t1161 - 0.441e3 * t1173 + 0.49e2 / 0.2e1 * t1177 - 0.35e2 / 0.3e1 * t1182 - 0.770e3 / 0.9e1 * t1264 + 0.140e3 / 0.3e1 * t1268;
1563 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1564 v4rho2sigma2[3] = 0.28e2 / 0.3e1 * t1188 - 0.7e1 * t1192 - 0.280e3 / 0.3e1 * t1195 - 0.140e3 / 0.3e1 * t1199 - 0.70e2 / 0.9e1 * t1203 + 0.1323e4 * t1045 - 0.147e3 * t1048 - 0.147e3 / 0.2e1 * t1052 + 0.35e2 / 0.3e1 * t1055 - 0.175e3 / 0.2e1 * t1057 + t1296 + 0.28e2 * t1231 + 0.14e2 * t1233 - 0.42e2 * t1236 + 0.21e2 * t1096 - 0.70e2 / 0.3e1 * t1106 + 0.35e2 * t1109 + 0.35e2 / 0.2e1 * t1113 - 0.70e2 * t1118 + 0.28e2 / 0.3e1 * t1120 - 0.14e2 * t1122 + t1329;
1565
1566 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1567 v4rho2sigma2[4] = v4rho2sigma2[1];
1568
1569 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1570 v4rho2sigma2[5] = v4rho2sigma2[2];
1571
1572 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1573 v4rho2sigma2[6] = v4rho2sigma2[5];
1574
1575 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1576 v4rho2sigma2[7] = v4rho2sigma2[4];
1577
1578 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1579 v4rho2sigma2[8] = v4rho2sigma2[6];
1580
1581 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1582 v4rho2sigma2[9] = v4rho2sigma2[3];
1583
1584 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1585 v4rho2sigma2[10] = v4rho2sigma2[7];
1586
1587 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1588 v4rho2sigma2[11] = v4rho2sigma2[8];
1589
1590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1591 v4rho2sigma2[12] = v4rho2sigma2[11];
1592
1593 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594 v4rho2sigma2[13] = v4rho2sigma2[10];
1595
1596 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597 v4rho2sigma2[14] = v4rho2sigma2[12];
1598
1599 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600 v4rho2sigma2[15] = v4rho2sigma2[9];
1601
1602 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603 v4rho2sigma2[16] = v4rho2sigma2[13];
1604
1605 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606 v4rho2sigma2[17] = v4rho2sigma2[14];
1607
1608 t1331 = t467 * t598;
1609 t1336 = params->b * t473 * t310 * t2 * t12;
1610 t1342 = t53 / t6 / t780 * t793 * t12;
1611 t1349 = t128 * t90 * t140 * t241;
1612 t1351 = t498 * t631;
1613 t1355 = t200 * t227 * t66 * t204;
1614 t1359 = t175 * t180 * t140 * t241;
1615 t1364 = t623 * t808 * t626 * t74 * t79;
1616 t1366 = t502 * t631;
1617 t1371 = t185 * t227 * t95 * t100;
1618 t1374 = t493 * t95 * t249;
1619 t1376 = t483 * t254;
1620 t1378 = t527 * t637;
1621 t1380 = t512 * t642;
1622 t1382 = t187 * t648;
1623 t1386 = t93 * t478 * t95 * t100;
1624 t1389 = t613 * t505 * t509;
1625 t1393 = t244 * t180 * t95 * t249;
1626 t1395 = t4 * (-0.24e2 * t1331 + 0.56e2 / 0.3e1 * t1336 - 0.8e1 / 0.3e1 * t1342) * t46 - 0.1323e4 / 0.8e1 * t1349 - 0.3969e4 / 0.32e2 * t1351 + 0.35e2 / 0.2e1 * t1355 + 0.147e3 / 0.8e1 * t1359 + 0.1029e4 / 0.8e1 * t1364 + 0.441e3 / 0.32e2 * t1366 + t606 * t46 + 0.7e1 / 0.2e1 * t1371 + 0.35e2 / 0.8e1 * t1374 - 0.7e1 / 0.4e1 * t1376 - 0.35e2 / 0.32e2 * t1378 - 0.35e2 / 0.16e2 * t1380 + 0.7e1 / 0.8e1 * t1382 - 0.21e2 / 0.4e1 * t1386 - 0.441e3 / 0.8e1 * t1389 - 0.105e3 / 0.16e2 * t1393;
1627 t1398 = t525 * t90 * t95 * t529;
1628 t1400 = t487 * t254;
1629 t1404 = t506 * t507 * t28 * t79;
1630 t1408 = t335 * t61 * t95 * t637;
1631 t1414 = t289 * t15 * t95 * t97 * t747 * t247;
1632 t1416 = t516 * t642;
1633 t1418 = t191 * t648;
1634 t1423 = t602 * t61 * t325 * t627;
1635 t1430 = 0.1e1 / t307 * t15 * t325 * t625 * t29 * t158;
1636 t1434 = t473 * t15 * t325 * t627;
1637 t1437 = t609 * t194 * t197;
1638 t1440 = t617 * t194 * t521;
1639 t1442 = t490 * t537;
1640 t1447 = t634 * t194 * t97 * t636 * t79;
1641 t1451 = t519 * t97 * t641 * t79;
1642 t1455 = t195 * t97 * t647 * t79;
1643 t1457 = 0.105e3 / 0.8e1 * t1398 + 0.21e2 / 0.8e1 * t1400 - 0.1323e4 / 0.32e2 * t1404 - 0.105e3 / 0.64e2 * t1408 + 0.35e2 / 0.32e2 * t1414 + 0.105e3 / 0.32e2 * t1416 - 0.21e2 / 0.16e2 * t1418 - t607 * t80 - 0.1029e4 / 0.32e2 * t1423 + 0.1715e4 / 0.16e2 * t1430 + 0.343e3 * t1434 + 0.21e2 / 0.2e1 * t1437 + 0.105e3 / 0.8e1 * t1440 - 0.21e2 / 0.4e1 * t1442 + 0.105e3 / 0.32e2 * t1447 - 0.105e3 / 0.16e2 * t1451 + 0.21e2 / 0.8e1 * t1455;
1644 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645 v4rhosigma3[0] = t1395 + t1457;
1646
1647 t1464 = 0.441e3 / 0.2e1 * t1349;
1648 t1468 = t200 * t259 * t66 * t204;
1649 t1472 = t175 * t211 * t140 * t241;
1650 t1476 = t128 * t103 * t140 * t241;
1651 t1479 = 0.49e2 / 0.2e1 * t1359;
1652 t1484 = t4 * (-0.48e2 * t1331 + 0.112e3 / 0.3e1 * t1336 - 0.16e2 / 0.3e1 * t1342) * t46 - t1464 - 0.3969e4 / 0.16e2 * t1351 + 0.35e2 / 0.3e1 * t1468 + 0.49e2 / 0.8e1 * t1472 - 0.441e3 / 0.8e1 * t1476 + 0.35e2 / 0.3e1 * t1355 + t1479 + 0.1029e4 / 0.4e1 * t1364 + 0.441e3 / 0.16e2 * t1366 + t653 * t46 + 0.7e1 / 0.3e1 * t1371;
1653 t1485 = 0.35e2 / 0.6e1 * t1374;
1654 t1486 = 0.7e1 / 0.3e1 * t1376;
1655 t1491 = 0.147e3 / 0.2e1 * t1389;
1656 t1492 = 0.35e2 / 0.4e1 * t1393;
1657 t1493 = 0.35e2 / 0.2e1 * t1398;
1658 t1494 = 0.7e1 / 0.2e1 * t1400;
1659 t1497 = t1485 - t1486 - 0.35e2 / 0.16e2 * t1378 - 0.35e2 / 0.8e1 * t1380 + 0.7e1 / 0.4e1 * t1382 - 0.7e1 / 0.2e1 * t1386 - t1491 - t1492 + t1493 + t1494 - 0.1323e4 / 0.16e2 * t1404 - 0.105e3 / 0.32e2 * t1408;
1660 t1504 = t185 * t259 * t95 * t100;
1661 t1508 = t93 * t544 * t95 * t100;
1662 t1511 = t660 * t505 * t509;
1663 t1515 = t244 * t211 * t95 * t249;
1664 t1519 = t525 * t103 * t95 * t529;
1665 t1521 = t553 * t254;
1666 t1524 = t559 * t95 * t249;
1667 t1526 = t549 * t254;
1668 t1529 = 0.35e2 / 0.16e2 * t1414 + 0.105e3 / 0.16e2 * t1416 - 0.21e2 / 0.8e1 * t1418 + 0.7e1 / 0.3e1 * t1504 - 0.7e1 / 0.2e1 * t1508 - 0.147e3 / 0.8e1 * t1511 - 0.35e2 / 0.16e2 * t1515 + 0.35e2 / 0.8e1 * t1519 + 0.7e1 / 0.8e1 * t1521 + 0.35e2 / 0.24e2 * t1524 - 0.7e1 / 0.12e2 * t1526 - t654 * t80;
1669 t1534 = 0.35e2 / 0.2e1 * t1440;
1670 t1535 = 0.7e1 * t1442;
1671 t1540 = t656 * t194 * t197;
1672 t1543 = t664 * t194 * t521;
1673 t1545 = t556 * t537;
1674 t1547 = -0.1029e4 / 0.16e2 * t1423 + 0.1715e4 / 0.8e1 * t1430 + 0.686e3 * t1434 + 0.7e1 * t1437 + t1534 - t1535 + 0.105e3 / 0.16e2 * t1447 - 0.105e3 / 0.8e1 * t1451 + 0.21e2 / 0.4e1 * t1455 + 0.7e1 * t1540 + 0.35e2 / 0.8e1 * t1543 - 0.7e1 / 0.4e1 * t1545;
1675 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1676 v4rhosigma3[1] = t1484 + t1497 + t1529 + t1547;
1677
1678 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1679 v4rhosigma3[2] = v4rhosigma3[0];
1680
1681 t1556 = t684 * t194 * t197;
1682 t1560 = t200 * t271 * t66 * t204;
1683 t1569 = t4 * (-0.96e2 * t1331 + 0.224e3 / 0.3e1 * t1336 - 0.32e2 / 0.3e1 * t1342) * t46 + 0.7e1 / 0.2e1 * t1556 + 0.35e2 / 0.6e1 * t1560 - t1464 - 0.3969e4 / 0.8e1 * t1351 + 0.70e2 / 0.3e1 * t1468 + 0.49e2 / 0.2e1 * t1472 - 0.441e3 / 0.2e1 * t1476 + t1479 + 0.1029e4 / 0.2e1 * t1364 + 0.441e3 / 0.8e1 * t1366 + t681 * t46;
1684 t1576 = t1485 - t1486 - 0.35e2 / 0.8e1 * t1378 - 0.35e2 / 0.4e1 * t1380 + 0.7e1 / 0.2e1 * t1382 - t1491 - t1492 + t1493 + t1494 - 0.1323e4 / 0.8e1 * t1404 - 0.105e3 / 0.16e2 * t1408 + 0.35e2 / 0.8e1 * t1414;
1685 t1590 = t93 * t579 * t95 * t100;
1686 t1594 = t185 * t271 * t95 * t100;
1687 t1596 = 0.105e3 / 0.8e1 * t1416 - 0.21e2 / 0.4e1 * t1418 + 0.14e2 / 0.3e1 * t1504 - 0.7e1 * t1508 - 0.147e3 / 0.2e1 * t1511 - 0.35e2 / 0.4e1 * t1515 + 0.35e2 / 0.2e1 * t1519 + 0.7e1 / 0.2e1 * t1521 + 0.35e2 / 0.6e1 * t1524 - 0.7e1 / 0.3e1 * t1526 - 0.7e1 / 0.4e1 * t1590 + 0.7e1 / 0.6e1 * t1594;
1688 t1607 = -t682 * t80 - 0.1029e4 / 0.8e1 * t1423 + 0.1715e4 / 0.4e1 * t1430 + 0.1372e4 * t1434 + t1534 - t1535 + 0.105e3 / 0.8e1 * t1447 - 0.105e3 / 0.4e1 * t1451 + 0.21e2 / 0.2e1 * t1455 + 0.14e2 * t1540 + 0.35e2 / 0.2e1 * t1543 - 0.7e1 * t1545;
1689 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1690 v4rhosigma3[3] = t1569 + t1576 + t1596 + t1607;
1691
1692 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1693 v4rhosigma3[4] = v4rhosigma3[1];
1694
1695 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1696 v4rhosigma3[5] = v4rhosigma3[2];
1697
1698 t1631 = 0.21e2 * t1556 + t699 * t46 + 0.35e2 * t1560 - 0.3969e4 / 0.4e1 * t1351 + 0.147e3 / 0.2e1 * t1472 - 0.1323e4 / 0.2e1 * t1476 + 0.1029e4 * t1364 + 0.441e3 / 0.4e1 * t1366 - t700 * t80 + t4 * (-0.192e3 * t1331 + 0.448e3 / 0.3e1 * t1336 - 0.64e2 / 0.3e1 * t1342) * t46 - 0.35e2 / 0.4e1 * t1378 - 0.35e2 / 0.2e1 * t1380 + 0.7e1 * t1382 - 0.1323e4 / 0.4e1 * t1404 - 0.105e3 / 0.8e1 * t1408 + 0.35e2 / 0.4e1 * t1414 + 0.105e3 / 0.4e1 * t1416;
1699 t1649 = -0.21e2 / 0.2e1 * t1418 - 0.441e3 / 0.2e1 * t1511 - 0.105e3 / 0.4e1 * t1515 + 0.105e3 / 0.2e1 * t1519 + 0.21e2 / 0.2e1 * t1521 + 0.35e2 / 0.2e1 * t1524 - 0.7e1 * t1526 - 0.21e2 / 0.2e1 * t1590 + 0.7e1 * t1594 - 0.1029e4 / 0.4e1 * t1423 + 0.1715e4 / 0.2e1 * t1430 + 0.2744e4 * t1434 + 0.105e3 / 0.4e1 * t1447 - 0.105e3 / 0.2e1 * t1451 + 0.21e2 * t1455 + 0.105e3 / 0.2e1 * t1543 - 0.21e2 * t1545;
1700 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1701 v4rhosigma3[6] = t1631 + t1649;
1702
1703 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1704 v4rhosigma3[7] = v4rhosigma3[3];
1705
1706 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1707 v4rhosigma3[8] = v4rhosigma3[4];
1708
1709 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1710 v4rhosigma3[9] = v4rhosigma3[5];
1711
1712 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1713 v4rhosigma3[10] = v4rhosigma3[9];
1714
1715 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1716 v4rhosigma3[11] = v4rhosigma3[8];
1717
1718 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1719 v4rhosigma3[12] = v4rhosigma3[10];
1720
1721 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1722 v4rhosigma3[13] = v4rhosigma3[7];
1723
1724 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1725 v4rhosigma3[14] = v4rhosigma3[11];
1726
1727 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1728 v4rhosigma3[15] = v4rhosigma3[12];
1729
1730 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1731 v4rhosigma3[16] = v4rhosigma3[6];
1732
1733 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1734 v4rhosigma3[17] = v4rhosigma3[13];
1735
1736 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1737 v4rhosigma3[18] = v4rhosigma3[14];
1738
1739 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1740 v4rhosigma3[19] = v4rhosigma3[15];
1741
1742 t1652 = t602 * t90 * t325 * t627;
1743 t1656 = t623 * t325 * t625 * t74;
1744 t1662 = t775 * t15 * t325 * t625 * t28 * t158;
1745 t1666 = t175 * t227 * t140 * t241;
1746 t1667 = 0.147e3 / 0.4e1 * t1666;
1747 t1672 = t237 * t237;
1748 t1676 = 0.1e1 / t789 / t47 * t15 * t808 * t1672 * t124 * t18 * t26;
1749 t1680 = params->b * t602 * t310 * t12;
1750 t1686 = t3 / t6 / t307 * t793 * t12;
1751 t1692 = t93 * t606 * t95 * t100;
1752 t1696 = t244 * t227 * t95 * t249;
1753 t1697 = 0.105e3 / 0.8e1 * t1696;
1754 t1698 = t614 * t631;
1755 t1702 = t335 * t90 * t95 * t637;
1756 t1706 = t235 * t238 * t239 * t98;
1757 t1712 = t113 * t15 * t95 * t97 * t747 * t640;
1758 t1714 = t610 * t254;
1759 t1715 = 0.21e2 / 0.4e1 * t1714;
1760 t1716 = t618 * t642;
1761 t1718 = t231 * t648;
1762 t1722 = t635 * t97 * t36 * t646;
1763 t1727 = t246 * t97 * t158 / t124;
1764 t1733 = t96 * t97 * t74 / t28 / t124;
1765 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1766 v4sigma4[0] = -0.1029e4 / 0.8e1 * t1652 - 0.13377e5 / 0.128e3 * t1656 - 0.5145e4 / 0.128e3 * t1662 + t1667 + 0.7203e4 / 0.32e2 * t1676 + t4 * (-0.4e1 * t1680 + t1686) * t46 - 0.7e1 * t1692 - t1697 + 0.441e3 / 0.8e1 * t1698 - 0.105e3 / 0.16e2 * t1702 + 0.735e3 / 0.128e3 * t1706 - 0.105e3 / 0.256e3 * t1712 + t1715 + 0.105e3 / 0.8e1 * t1716 - 0.21e2 / 0.4e1 * t1718 + 0.315e3 / 0.64e2 * t1722 - 0.525e3 / 0.64e2 * t1727 + 0.105e3 / 0.32e2 * t1733;
1767
1768 t1737 = t602 * t103 * t325 * t627;
1769 t1747 = t175 * t259 * t140 * t241;
1770 t1749 = t661 * t631;
1771 t1757 = -0.1029e4 / 0.32e2 * t1737 - 0.3087e4 / 0.16e2 * t1652 - 0.5145e4 / 0.64e2 * t1662 - 0.13377e5 / 0.64e2 * t1656 + t1667 + 0.7203e4 / 0.16e2 * t1676 + 0.1323e4 / 0.16e2 * t1698 + 0.735e3 / 0.64e2 * t1706 + 0.147e3 / 0.8e1 * t1747 + 0.441e3 / 0.32e2 * t1749 + t4 * (-0.8e1 * t1680 + 0.2e1 * t1686) * t46 - 0.7e1 / 0.2e1 * t1692 - t1697;
1772 t1767 = t93 * t653 * t95 * t100;
1773 t1771 = t244 * t259 * t95 * t249;
1774 t1775 = t335 * t103 * t95 * t637;
1775 t1777 = t657 * t254;
1776 t1779 = t665 * t642;
1777 t1781 = t263 * t648;
1778 t1783 = -0.315e3 / 0.32e2 * t1702 - 0.105e3 / 0.128e3 * t1712 + t1715 + 0.315e3 / 0.16e2 * t1716 - 0.63e2 / 0.8e1 * t1718 + 0.315e3 / 0.32e2 * t1722 - 0.525e3 / 0.32e2 * t1727 + 0.105e3 / 0.16e2 * t1733 - 0.21e2 / 0.4e1 * t1767 - 0.105e3 / 0.16e2 * t1771 - 0.105e3 / 0.64e2 * t1775 + 0.21e2 / 0.8e1 * t1777 + 0.105e3 / 0.32e2 * t1779 - 0.21e2 / 0.16e2 * t1781;
1779 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1780 v4sigma4[1] = t1757 + t1783;
1781
1782 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1783 v4sigma4[2] = v4sigma4[0];
1784
1785 t1785 = 0.1029e4 / 0.4e1 * t1652;
1786 t1790 = 0.441e3 / 0.4e1 * t1698;
1787 t1796 = t175 * t271 * t140 * t241;
1788 t1804 = 0.105e3 / 0.8e1 * t1702;
1789 t1806 = -0.1029e4 / 0.8e1 * t1737 - t1785 - 0.5145e4 / 0.32e2 * t1662 - 0.13377e5 / 0.32e2 * t1656 + 0.49e2 / 0.2e1 * t1666 + 0.7203e4 / 0.8e1 * t1676 + t1790 + 0.735e3 / 0.32e2 * t1706 + 0.49e2 * t1747 + 0.441e3 / 0.8e1 * t1749 + 0.49e2 / 0.8e1 * t1796 + t4 * (-0.16e2 * t1680 + 0.4e1 * t1686) * t46 - 0.35e2 / 0.4e1 * t1696 - t1804 - 0.105e3 / 0.64e2 * t1712;
1790 t1808 = 0.105e3 / 0.4e1 * t1716;
1791 t1809 = 0.21e2 / 0.2e1 * t1718;
1792 t1821 = t93 * t681 * t95 * t100;
1793 t1825 = t244 * t271 * t95 * t249;
1794 t1827 = t685 * t254;
1795 t1829 = 0.7e1 / 0.2e1 * t1714 + t1808 - t1809 + 0.315e3 / 0.16e2 * t1722 - 0.525e3 / 0.16e2 * t1727 + 0.105e3 / 0.8e1 * t1733 - 0.7e1 * t1767 - 0.35e2 / 0.2e1 * t1771 - 0.105e3 / 0.16e2 * t1775 + 0.7e1 * t1777 + 0.105e3 / 0.8e1 * t1779 - 0.21e2 / 0.4e1 * t1781 - 0.7e1 / 0.2e1 * t1821 - 0.35e2 / 0.16e2 * t1825 + 0.7e1 / 0.8e1 * t1827;
1796 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1797 v4sigma4[3] = t1806 + t1829;
1798
1799 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1800 v4sigma4[4] = v4sigma4[1];
1801
1802 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1803 v4sigma4[5] = v4sigma4[2];
1804
1805 t1844 = -0.3087e4 / 0.8e1 * t1737 - t1785 - 0.5145e4 / 0.16e2 * t1662 - 0.13377e5 / 0.16e2 * t1656 + 0.7203e4 / 0.4e1 * t1676 + t1790 + 0.735e3 / 0.16e2 * t1706 + 0.147e3 / 0.2e1 * t1747 + 0.1323e4 / 0.8e1 * t1749 + 0.147e3 / 0.4e1 * t1796 + t4 * (-0.32e2 * t1680 + 0.8e1 * t1686) * t46 - t1804 - 0.105e3 / 0.32e2 * t1712;
1806 t1858 = t93 * t699 * t95 * t100;
1807 t1860 = t1808 - t1809 + 0.315e3 / 0.8e1 * t1722 - 0.525e3 / 0.8e1 * t1727 + 0.105e3 / 0.4e1 * t1733 - 0.105e3 / 0.4e1 * t1771 - 0.315e3 / 0.16e2 * t1775 + 0.21e2 / 0.2e1 * t1777 + 0.315e3 / 0.8e1 * t1779 - 0.63e2 / 0.4e1 * t1781 - 0.21e2 / 0.2e1 * t1821 - 0.105e3 / 0.8e1 * t1825 + 0.21e2 / 0.4e1 * t1827 - 0.7e1 / 0.4e1 * t1858;
1808 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1809 v4sigma4[6] = t1844 + t1860;
1810
1811 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1812 v4sigma4[7] = v4sigma4[3];
1813
1814 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1815 v4sigma4[8] = v4sigma4[4];
1816
1817 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1818 v4sigma4[9] = v4sigma4[5];
1819
1820 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1821 v4sigma4[10] = -0.13377e5 / 0.8e1 * t1656 - 0.5145e4 / 0.8e1 * t1662 - 0.1029e4 * t1737 + 0.7203e4 / 0.2e1 * t1676 + 0.147e3 * t1796 + t4 * (-0.64e2 * t1680 + 0.16e2 * t1686) * t46 + 0.735e3 / 0.8e1 * t1706 - 0.105e3 / 0.16e2 * t1712 + 0.441e3 * t1749 - 0.105e3 / 0.2e1 * t1775 - 0.105e3 / 0.2e1 * t1825 - 0.14e2 * t1858 + 0.315e3 / 0.4e1 * t1722 - 0.525e3 / 0.4e1 * t1727 + 0.105e3 / 0.2e1 * t1733 + 0.105e3 * t1779 - 0.42e2 * t1781 + 0.21e2 * t1827;
1822
1823 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1824 v4sigma4[11] = v4sigma4[6];
1825
1826 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1827 v4sigma4[12] = v4sigma4[7];
1828
1829 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1830 v4sigma4[13] = v4sigma4[8];
1831
1832 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1833 v4sigma4[14] = v4sigma4[9];
1834
1835 #ifndef XC_DONT_COMPILE_MXC
1836
1837 if(order < 5) return;
1838
1839
1840 #endif
1841
1842 #endif
1843
1844 #endif
1845
1846 #endif
1847
1848 #endif
1849
1850
1851 }
1852
1853