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_optx.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t3, t5, t6, t7, t8, t9, t10;
24 double t11, t12, t13, t14, t15, t16, t17, t18;
25 double t23, t26, t28, t29, t30, t31, t34, t36;
26
27 #ifndef XC_DONT_COMPILE_VXC
28 double t41, t42, t43, t47, t51, t52, t53, t54;
29 double t55, t58, t62, t66, t69, t75;
30
31 #ifndef XC_DONT_COMPILE_FXC
32 double t78, t79, t86, t87, t91, t94, t95, t100;
33 double t101, t102, t103, t104, t106, t107, t108, t110;
34 double t113, t118, t127, t128, t130, t132, t136, t145;
35 double t149, t153;
36
37 #ifndef XC_DONT_COMPILE_KXC
38 double t157, t158, t171, t174, t179, t183, t187, t188;
39 double t189, t190, t194, t196, t199, t204, t217, t222;
40 double t226, t233, t239, t240, t244, t248, t257, t261;
41 double t265;
42
43 #ifndef XC_DONT_COMPILE_LXC
44 double t296, t304, t310;
45 #endif
46
47 #endif
48
49 #endif
50
51 #endif
52
53 #endif
54
55
56 gga_x_optx_params *params;
57
58 assert(p->params != NULL);
59 params = (gga_x_optx_params * )(p->params);
60
61 t1 = M_CBRT3;
62 t3 = POW_1_3(0.1e1 / M_PI);
63 t5 = M_CBRT4;
64 t6 = t5 * t5;
65 t7 = t1 * t3 * t6;
66 t8 = M_CBRT2;
67 t9 = t8 * t8;
68 t10 = POW_1_3(rho[0]);
69 t11 = t9 * t10;
70 t12 = params->gamma * params->gamma;
71 t13 = params->b * t12;
72 t14 = sigma[0] * sigma[0];
73 t15 = t13 * t14;
74 t16 = rho[0] * rho[0];
75 t17 = t16 * t16;
76 t18 = t17 * rho[0];
77 t23 = t10 * t10;
78 t26 = t9 / t23 / t16;
79 t28 = params->gamma * sigma[0] * t26 + 0.1e1;
80 t29 = t28 * t28;
81 t30 = 0.1e1 / t29;
82 t31 = t8 / t10 / t18 * t30;
83 t34 = 0.2e1 * t15 * t31 + params->a;
84 t36 = t7 * t11 * t34;
85 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
86 zk[0] = -0.3e1 / 0.16e2 * t36;
87
88 #ifndef XC_DONT_COMPILE_VXC
89
90 if(order < 1) return;
91
92
93 t41 = t10 * rho[0] * t1 * t3;
94 t42 = t6 * t9;
95 t43 = t17 * t16;
96 t47 = t8 / t10 / t43 * t30;
97 t51 = params->b * t12 * params->gamma;
98 t52 = t14 * sigma[0];
99 t53 = t17 * t17;
100 t54 = t53 * rho[0];
101 t55 = 0.1e1 / t54;
102 t58 = 0.1e1 / t29 / t28;
103 t62 = -0.32e2 / 0.3e1 * t15 * t47 + 0.64e2 / 0.3e1 * t51 * t52 * t55 * t58;
104 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
105 vrho[0] = -t36 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t42 * t62;
106
107 t66 = t13 * sigma[0];
108 t69 = 0.1e1 / t53;
109 t75 = t42 * (-0.8e1 * t51 * t14 * t69 * t58 + 0.4e1 * t66 * t31);
110 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
111 vsigma[0] = -0.3e1 / 0.16e2 * t41 * t75;
112
113 #ifndef XC_DONT_COMPILE_FXC
114
115 if(order < 2) return;
116
117
118 t78 = 0.1e1 / t23;
119 t79 = t9 * t78;
120 t86 = t16 * rho[0];
121 t87 = t17 * t86;
122 t91 = t8 / t10 / t87 * t30;
123 t94 = t53 * t16;
124 t95 = 0.1e1 / t94;
125 t100 = t12 * t12;
126 t101 = params->b * t100;
127 t102 = t14 * t14;
128 t103 = t101 * t102;
129 t104 = t53 * t17;
130 t106 = 0.1e1 / t23 / t104;
131 t107 = t29 * t29;
132 t108 = 0.1e1 / t107;
133 t110 = t106 * t108 * t9;
134 t113 = 0.608e3 / 0.9e1 * t15 * t91 - 0.2752e4 / 0.9e1 * t51 * t52 * t95 * t58 + 0.512e3 / 0.3e1 * t103 * t110;
135 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
136 v2rho2[0] = -t7 * t79 * t34 / 0.12e2 - t7 * t11 * t62 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t42 * t113;
137
138 t118 = t10 * t1 * t3;
139 t127 = t101 * t52;
140 t128 = t53 * t86;
141 t130 = 0.1e1 / t23 / t128;
142 t132 = t130 * t108 * t9;
143 t136 = t42 * (-0.64e2 / 0.3e1 * t66 * t47 + 0.320e3 / 0.3e1 * t51 * t14 * t55 * t58 - 0.64e2 * t127 * t132);
144 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
145 v2rhosigma[0] = -t118 * t75 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t136;
146
147 t145 = t101 * t14;
148 t149 = 0.1e1 / t23 / t94 * t108 * t9;
149 t153 = t42 * (-0.32e2 * t51 * sigma[0] * t69 * t58 + 0.4e1 * t13 * t31 + 0.24e2 * t145 * t149);
150 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
151 v2sigma2[0] = -0.3e1 / 0.16e2 * t41 * t153;
152
153 #ifndef XC_DONT_COMPILE_KXC
154
155 if(order < 3) return;
156
157
158 t157 = 0.1e1 / t23 / rho[0];
159 t158 = t9 * t157;
160 t171 = t8 / t10 / t53 * t30;
161 t174 = 0.1e1 / t128;
162 t179 = t53 * t18;
163 t183 = 0.1e1 / t23 / t179 * t108 * t9;
164 t187 = params->b * t100 * params->gamma;
165 t188 = t102 * sigma[0];
166 t189 = t187 * t188;
167 t190 = t53 * t53;
168 t194 = 0.1e1 / t107 / t28;
169 t196 = 0.1e1 / t10 / t190 * t194 * t8;
170 t199 = -0.13376e5 / 0.27e2 * t15 * t171 + 0.102016e6 / 0.27e2 * t51 * t52 * t174 * t58 - 0.4608e4 * t103 * t183 + 0.32768e5 / 0.9e1 * t189 * t196;
171 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
172 v3rho3[0] = t7 * t158 * t34 / 0.18e2 - t7 * t79 * t62 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t113 - 0.3e1 / 0.16e2 * t41 * t42 * t199;
173
174 t204 = t78 * t1 * t3;
175 t217 = t187 * t102;
176 t222 = 0.1e1 / t10 / t53 / t87 * t194 * t8;
177 t226 = t42 * (0.1216e4 / 0.9e1 * t66 * t91 - 0.10688e5 / 0.9e1 * t51 * t14 * t95 * t58 + 0.1600e4 * t127 * t110 - 0.4096e4 / 0.3e1 * t217 * t222);
178 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
179 v3rho2sigma[0] = -t204 * t75 / 0.12e2 - t118 * t136 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t226;
180
181 t233 = t55 * t58;
182 t239 = t187 * t52;
183 t240 = t53 * t43;
184 t244 = 0.1e1 / t10 / t240 * t194 * t8;
185 t248 = t42 * (-0.64e2 / 0.3e1 * t13 * t47 + 0.896e3 / 0.3e1 * t51 * t233 * sigma[0] - 0.512e3 * t145 * t132 + 0.512e3 * t239 * t244);
186 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
187 v3rhosigma2[0] = -t118 * t153 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t248;
188
189 t257 = t187 * t14;
190 t261 = 0.1e1 / t10 / t179 * t194 * t8;
191 t265 = t42 * (0.144e3 * t101 * sigma[0] * t149 - 0.48e2 * t51 * t69 * t58 - 0.192e3 * t257 * t261);
192 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
193 v3sigma3[0] = -0.3e1 / 0.16e2 * t41 * t265;
194
195 #ifndef XC_DONT_COMPILE_LXC
196
197 if(order < 4) return;
198
199
200 t296 = t190 * rho[0];
201 t304 = params->b * t100 * t12;
202 t310 = 0.1e1 / t107 / t29;
203 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
204 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t26 * t34 + 0.2e1 / 0.9e1 * t7 * t158 * t62 - t7 * t79 * t113 / 0.2e1 - t7 * t11 * t199 - 0.3e1 / 0.16e2 * t41 * t42 * (0.334400e6 / 0.81e2 * t15 * t8 / t10 / t54 * t30 - 0.3794560e7 / 0.81e2 * t51 * t52 / t104 * t58 + 0.2516480e7 / 0.27e2 * t103 / t23 / t240 * t108 * t9 - 0.4259840e7 / 0.27e2 * t189 / t10 / t296 * t194 * t8 + 0.2621440e7 / 0.27e2 * t304 * t102 * t14 / t190 / t17 * t310);
205
206 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
207 v4rho3sigma[0] = t157 * t1 * t3 * t75 / 0.18e2 - t204 * t136 / 0.4e1 - 0.3e1 / 0.4e1 * t118 * t226 - 0.3e1 / 0.16e2 * t41 * t42 * (-0.26752e5 / 0.27e2 * t66 * t171 + 0.359552e6 / 0.27e2 * t51 * t14 * t174 * t58 - 0.267904e6 / 0.9e1 * t127 * t183 + 0.495616e6 / 0.9e1 * t217 * t196 - 0.327680e6 / 0.9e1 * t304 * t188 / t190 / t86 * t310);
208
209 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
210 v4rho2sigma2[0] = -t204 * t153 / 0.12e2 - t118 * t248 / 0.2e1 - 0.3e1 / 0.16e2 * t41 * t42 * (0.1216e4 / 0.9e1 * t13 * t91 - 0.26240e5 / 0.9e1 * t51 * t95 * t58 * sigma[0] + 0.25088e5 / 0.3e1 * t101 * t106 * t108 * t14 * t9 - 0.54784e5 / 0.3e1 * t239 * t222 + 0.40960e5 / 0.3e1 * t304 * t102 / t190 / t16 * t310);
211
212 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
213 v4rhosigma3[0] = -t118 * t265 / 0.4e1 - 0.3e1 / 0.16e2 * t41 * t42 * (0.384e3 * t51 * t233 - 0.1920e4 * t101 * t130 * t108 * sigma[0] * t9 + 0.5632e4 * t257 * t244 - 0.5120e4 * t304 * t52 / t296 * t310);
214
215 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
216 v4sigma4[0] = -0.3e1 / 0.16e2 * t41 * t42 * (0.288e3 * t101 * t149 - 0.1536e4 * t187 * sigma[0] * t261 + 0.1920e4 * t304 * t14 / t190 * t310);
217
218 #ifndef XC_DONT_COMPILE_MXC
219
220 if(order < 5) return;
221
222
223 #endif
224
225 #endif
226
227 #endif
228
229 #endif
230
231 #endif
232
233
234 }
235
236
237 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 *))238 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
239 {
240
241 #ifndef XC_DONT_COMPILE_EXC
242 double t1, t3, t4, t5, t6, t7, t8, t9;
243 double t10, t11, t12, t13, t14, t16, t19, t21;
244 double t23, t24, t25, t28, t30;
245
246 #ifndef XC_DONT_COMPILE_VXC
247 double t34, t35, t36, t38, t43, t44, t45, t46;
248 double t47, t50, t54, t61, t67;
249
250 #ifndef XC_DONT_COMPILE_FXC
251 double t70, t71, t77, t78, t80, t85, t86, t91;
252 double t92, t93, t94, t96, t98, t99, t103, t107;
253 double t118, t120, t126, t137, t143;
254
255 #ifndef XC_DONT_COMPILE_KXC
256 double t147, t148, t159, t164, t169, t171, t177, t178;
257 double t179, t181, t184, t188, t192, t210, t216, t224;
258 double t232, t234, t240, t251, t257;
259
260 #ifndef XC_DONT_COMPILE_LXC
261 double t289, t297, t303;
262 #endif
263
264 #endif
265
266 #endif
267
268 #endif
269
270 #endif
271
272
273 gga_x_optx_params *params;
274
275 assert(p->params != NULL);
276 params = (gga_x_optx_params * )(p->params);
277
278 t1 = M_CBRT3;
279 t3 = POW_1_3(0.1e1 / M_PI);
280 t4 = t1 * t3;
281 t5 = M_CBRT4;
282 t6 = t5 * t5;
283 t7 = POW_1_3(rho[0]);
284 t8 = t6 * t7;
285 t9 = params->gamma * params->gamma;
286 t10 = params->b * t9;
287 t11 = sigma[0] * sigma[0];
288 t12 = rho[0] * rho[0];
289 t13 = t12 * t12;
290 t14 = t13 * rho[0];
291 t16 = 0.1e1 / t7 / t14;
292 t19 = t7 * t7;
293 t21 = 0.1e1 / t19 / t12;
294 t23 = params->gamma * sigma[0] * t21 + 0.1e1;
295 t24 = t23 * t23;
296 t25 = 0.1e1 / t24;
297 t28 = t10 * t11 * t16 * t25 + params->a;
298 t30 = t4 * t8 * t28;
299 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
300 zk[0] = -0.3e1 / 0.8e1 * t30;
301
302 #ifndef XC_DONT_COMPILE_VXC
303
304 if(order < 1) return;
305
306
307 t34 = t7 * rho[0] * t1;
308 t35 = t3 * t6;
309 t36 = t13 * t12;
310 t38 = 0.1e1 / t7 / t36;
311 t43 = params->b * t9 * params->gamma;
312 t44 = t11 * sigma[0];
313 t45 = t13 * t13;
314 t46 = t45 * rho[0];
315 t47 = 0.1e1 / t46;
316 t50 = 0.1e1 / t24 / t23;
317 t54 = -0.16e2 / 0.3e1 * t10 * t11 * t38 * t25 + 0.16e2 / 0.3e1 * t43 * t44 * t47 * t50;
318 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
319 vrho[0] = -t30 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t35 * t54;
320
321 t61 = 0.1e1 / t45;
322 t67 = t35 * (0.2e1 * t10 * sigma[0] * t16 * t25 - 0.2e1 * t43 * t11 * t61 * t50);
323 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
324 vsigma[0] = -0.3e1 / 0.8e1 * t34 * t67;
325
326 #ifndef XC_DONT_COMPILE_FXC
327
328 if(order < 2) return;
329
330
331 t70 = 0.1e1 / t19;
332 t71 = t6 * t70;
333 t77 = t12 * rho[0];
334 t78 = t13 * t77;
335 t80 = 0.1e1 / t7 / t78;
336 t85 = t45 * t12;
337 t86 = 0.1e1 / t85;
338 t91 = t9 * t9;
339 t92 = params->b * t91;
340 t93 = t11 * t11;
341 t94 = t45 * t13;
342 t96 = 0.1e1 / t19 / t94;
343 t98 = t24 * t24;
344 t99 = 0.1e1 / t98;
345 t103 = 0.304e3 / 0.9e1 * t10 * t11 * t80 * t25 - 0.688e3 / 0.9e1 * t43 * t44 * t86 * t50 + 0.128e3 / 0.3e1 * t92 * t93 * t96 * t99;
346 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
347 v2rho2[0] = -t4 * t71 * t28 / 0.6e1 - t4 * t8 * t54 - 0.3e1 / 0.8e1 * t34 * t35 * t103;
348
349 t107 = t7 * t1;
350 t118 = t45 * t77;
351 t120 = 0.1e1 / t19 / t118;
352 t126 = t35 * (-0.32e2 / 0.3e1 * t10 * sigma[0] * t38 * t25 + 0.80e2 / 0.3e1 * t43 * t11 * t47 * t50 - 0.16e2 * t92 * t44 * t120 * t99);
353 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
354 v2rhosigma[0] = -t107 * t67 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t126;
355
356 t137 = 0.1e1 / t19 / t85;
357 t143 = t35 * (0.6e1 * t92 * t11 * t137 * t99 - 0.8e1 * t43 * sigma[0] * t61 * t50 + 0.2e1 * t10 * t16 * t25);
358 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
359 v2sigma2[0] = -0.3e1 / 0.8e1 * t34 * t143;
360
361 #ifndef XC_DONT_COMPILE_KXC
362
363 if(order < 3) return;
364
365
366 t147 = 0.1e1 / t19 / rho[0];
367 t148 = t6 * t147;
368 t159 = 0.1e1 / t7 / t45;
369 t164 = 0.1e1 / t118;
370 t169 = t45 * t14;
371 t171 = 0.1e1 / t19 / t169;
372 t177 = params->b * t91 * params->gamma;
373 t178 = t93 * sigma[0];
374 t179 = t45 * t45;
375 t181 = 0.1e1 / t7 / t179;
376 t184 = 0.1e1 / t98 / t23;
377 t188 = -0.6688e4 / 0.27e2 * t10 * t11 * t159 * t25 + 0.25504e5 / 0.27e2 * t43 * t44 * t164 * t50 - 0.1152e4 * t92 * t93 * t171 * t99 + 0.4096e4 / 0.9e1 * t177 * t178 * t181 * t184;
378 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
379 v3rho3[0] = t4 * t148 * t28 / 0.9e1 - t4 * t71 * t54 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t103 - 0.3e1 / 0.8e1 * t34 * t35 * t188;
380
381 t192 = t70 * t1;
382 t210 = 0.1e1 / t7 / t45 / t78;
383 t216 = t35 * (0.608e3 / 0.9e1 * t10 * sigma[0] * t80 * t25 - 0.2672e4 / 0.9e1 * t43 * t11 * t86 * t50 + 0.400e3 * t92 * t44 * t96 * t99 - 0.512e3 / 0.3e1 * t177 * t93 * t210 * t184);
384 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
385 v3rho2sigma[0] = -t192 * t67 / 0.6e1 - t107 * t126 - 0.3e1 / 0.8e1 * t34 * t216;
386
387 t224 = t47 * t50;
388 t232 = t45 * t36;
389 t234 = 0.1e1 / t7 / t232;
390 t240 = t35 * (-0.32e2 / 0.3e1 * t10 * t38 * t25 + 0.224e3 / 0.3e1 * t43 * t224 * sigma[0] - 0.128e3 * t92 * t11 * t120 * t99 + 0.64e2 * t177 * t44 * t234 * t184);
391 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
392 v3rhosigma2[0] = -t107 * t143 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t240;
393
394 t251 = 0.1e1 / t7 / t169;
395 t257 = t35 * (-0.24e2 * t177 * t11 * t251 * t184 + 0.36e2 * t92 * sigma[0] * t137 * t99 - 0.12e2 * t43 * t61 * t50);
396 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
397 v3sigma3[0] = -0.3e1 / 0.8e1 * t34 * t257;
398
399 #ifndef XC_DONT_COMPILE_LXC
400
401 if(order < 4) return;
402
403
404 t289 = t179 * rho[0];
405 t297 = params->b * t91 * t9;
406 t303 = 0.1e1 / t98 / t24;
407 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
408 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t21 * t28 + 0.4e1 / 0.9e1 * t4 * t148 * t54 - t4 * t71 * t103 - 0.2e1 * t4 * t8 * t188 - 0.3e1 / 0.8e1 * t34 * t35 * (0.167200e6 / 0.81e2 * t10 * t11 / t7 / t46 * t25 - 0.948640e6 / 0.81e2 * t43 * t44 / t94 * t50 + 0.629120e6 / 0.27e2 * t92 * t93 / t19 / t232 * t99 - 0.532480e6 / 0.27e2 * t177 * t178 / t7 / t289 * t184 + 0.163840e6 / 0.27e2 * t297 * t93 * t11 / t179 / t13 * t303);
409
410 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
411 v4rho3sigma[0] = t147 * t1 * t67 / 0.9e1 - t192 * t126 / 0.2e1 - 0.3e1 / 0.2e1 * t107 * t216 - 0.3e1 / 0.8e1 * t34 * t35 * (-0.13376e5 / 0.27e2 * t10 * sigma[0] * t159 * t25 + 0.89888e5 / 0.27e2 * t43 * t11 * t164 * t50 - 0.66976e5 / 0.9e1 * t92 * t44 * t171 * t99 + 0.61952e5 / 0.9e1 * t177 * t93 * t181 * t184 - 0.20480e5 / 0.9e1 * t297 * t178 / t179 / t77 * t303);
412
413 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
414 v4rho2sigma2[0] = -t192 * t143 / 0.6e1 - t107 * t240 - 0.3e1 / 0.8e1 * t34 * t35 * (0.608e3 / 0.9e1 * t10 * t80 * t25 - 0.6560e4 / 0.9e1 * t43 * t86 * t50 * sigma[0] + 0.6272e4 / 0.3e1 * t92 * t96 * t99 * t11 - 0.6848e4 / 0.3e1 * t177 * t44 * t210 * t184 + 0.2560e4 / 0.3e1 * t297 * t93 / t179 / t12 * t303);
415
416 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
417 v4rhosigma3[0] = -t107 * t257 / 0.2e1 - 0.3e1 / 0.8e1 * t34 * t35 * (0.96e2 * t43 * t224 - 0.480e3 * t92 * t120 * t99 * sigma[0] + 0.704e3 * t177 * t11 * t234 * t184 - 0.320e3 * t297 * t44 / t289 * t303);
418
419 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
420 v4sigma4[0] = -0.3e1 / 0.8e1 * t34 * t35 * (0.72e2 * t92 * t137 * t99 - 0.192e3 * t177 * sigma[0] * t251 * t184 + 0.120e3 * t297 * t11 / t179 * t303);
421
422 #ifndef XC_DONT_COMPILE_MXC
423
424 if(order < 5) return;
425
426
427 #endif
428
429 #endif
430
431 #endif
432
433 #endif
434
435 #endif
436
437
438 }
439
440
441 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 *))442 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
443 {
444
445 #ifndef XC_DONT_COMPILE_EXC
446 double t1, t3, t4, t5, t6, t7, t8, t9;
447 double t10, t11, t13, t14, t15, t16, t17, t18;
448 double t19, t20, t21, t22, t23, t24, t26, t29;
449 double t33, t34, t35, t38, t40, t42, t43, t44;
450 double t45, t46, t47, t48, t49, t50, t52, t55;
451 double t59, t60, t61, t64, t66;
452
453 #ifndef XC_DONT_COMPILE_VXC
454 double t68, t69, t70, t71, t72, t73, t75, t76;
455 double t78, t79, t80, t81, t82, t84, t85, t86;
456 double t88, t93, t94, t95, t96, t97, t100, t104;
457 double t106, t107, t108, t109, t110, t112, t113, t114;
458 double t116, t117, t121, t122, t124, t125, t126, t127;
459 double t129, t130, t131, t133, t137, t138, t139, t140;
460 double t143, t147, t149, t150, t154, t155, t156, t160;
461 double t165, t169, t173, t178;
462
463 #ifndef XC_DONT_COMPILE_FXC
464 double t182, t184, t185, t186, t187, t188, t189, t191;
465 double t192, t193, t195, t196, t197, t199, t200, t201;
466 double t202, t203, t204, t206, t207, t209, t210, t212;
467 double t213, t215, t216, t217, t218, t220, t225, t226;
468 double t231, t232, t233, t234, t236, t238, t239, t243;
469 double t245, t246, t247, t248, t249, t250, t251, t253;
470 double t254, t255, t257, t258, t259, t260, t262, t263;
471 double t264, t266, t267, t268, t270, t271, t272, t273;
472 double t275, t278, t280, t282, t285, t286, t287, t289;
473 double t293, t294, t295, t296, t298, t301, t303, t304;
474 double t306, t309, t311, t314, t316, t319, t320, t322;
475 double t323, t324, t325, t326, t328, t329, t330, t331;
476 double t333, t334, t335, t336, t338, t339, t340, t342;
477 double t343, t344, t345, t346, t348, t353, t354, t359;
478 double t360, t362, t364, t365, t369, t371, t372, t373;
479 double t376, t378, t379, t391, t393, t398, t403, t404;
480 double t405, t410, t415, t426, t428, t433, t434, t445;
481 double t450, t451, t462, t467, t468;
482
483 #ifndef XC_DONT_COMPILE_KXC
484 double t475, t477, t481, t482, t483, t484, t486, t487;
485 double t488, t490, t492, t494, t495, t496, t498, t499;
486 double t501, t502, t503, t504, t506, t507, t508, t510;
487 double t512, t514, t515, t516, t518, t519, t520, t522;
488 double t523, t525, t528, t530, t532, t533, t534, t536;
489 double t537, t538, t539, t540, t542, t543, t545, t546;
490 double t547, t549, t550, t552, t554, t555, t557, t558;
491 double t559, t561, t562, t564, t565, t567, t570, t575;
492 double t580, t582, t588, t589, t590, t592, t595, t599;
493 double t601, t602, t603, t608, t609, t612, t613, t614;
494 double t616, t618, t620, t621, t622, t623, t625, t627;
495 double t628, t629, t631, t633, t635, t636, t637, t639;
496 double t641, t643, t645, t647, t649, t650, t652, t661;
497 double t663, t666, t667, t669, t671, t673, t675, t677;
498 double t679, t682, t683, t685, t686, t687, t689, t691;
499 double t693, t695, t697, t702, t705, t708, t710, t713;
500 double t718, t720, t722, t724, t726, t728, t730, t732;
501 double t734, t736, t738, t740, t742, t744, t751, t754;
502 double t757, t760, t763, t766, t768, t770, t772, t774;
503 double t777, t778, t780, t782, t784, t785, t786, t788;
504 double t790, t792, t797, t799, t802, t810, t811, t816;
505 double t817, t819, t820, t821, t823, t824, t825, t826;
506 double t828, t829, t831, t833, t834, t835, t838, t839;
507 double t841, t842, t843, t845, t846, t847, t849, t850;
508 double t851, t853, t854, t855, t856, t858, t859, t861;
509 double t866, t871, t873, t878, t879, t881, t884, t888;
510 double t890, t891, t892, t896, t900, t902, t903, t907;
511 double t910, t928, t933, t938, t942, t943, t947, t953;
512 double t954, t955, t957, t958, t959, t960, t964, t967;
513 double t968, t970, t974, t976, t977, t981, t984, t987;
514 double t992, t996, t1001, t1005, t1008, t1025, t1030, t1031;
515 double t1034, t1037, t1043, t1051, t1053, t1058, t1059, t1062;
516 double t1064, t1068, t1073, t1079, t1087, t1089, t1094, t1095;
517 double t1106, t1111, t1112, t1123, t1128, t1129;
518
519 #ifndef XC_DONT_COMPILE_LXC
520 double t1143, t1147, t1150, t1153, t1159, t1166, t1169, t1172;
521 double t1175, t1181, t1184, t1187, t1201, t1203, t1205, t1211;
522 double t1218, t1220, t1221, t1222, t1227, t1229, t1231, t1232;
523 double t1234, t1236, t1242, t1244, t1246, t1247, t1248, t1253;
524 double t1255, t1257, t1259, t1266, t1277, t1287, t1290, t1293;
525 double t1312, t1320, t1326, t1337, t1341, t1342, t1345, t1347;
526 double t1353, t1356, t1358, t1363, t1370, t1373, t1374, t1376;
527 double t1377, t1379, t1386, t1389, t1392, t1395, t1401, t1402;
528 double t1411, t1415, t1429, t1449, t1450, t1451, t1453, t1465;
529 double t1466, t1472, t1474, t1478, t1484, t1486, t1492, t1494;
530 double t1497, t1503, t1510, t1516, t1531, t1536, t1538, t1541;
531 double t1545, t1550, t1558, t1583, t1588, t1605, t1612, t1622;
532 double t1625, t1631, t1633, t1640, t1645, t1651, t1654, t1657;
533 double t1663, t1666, t1673, t1679, t1681, t1687, t1690, t1692;
534 double t1695, t1703, t1706, t1735, t1760, t1776, t1790, t1803;
535 double t1810, t1864, t1870, t1872, t1882, t1887, t1890, t1893;
536 double t1896, t1899, t1902, t1904, t1907, t1919, t1950, t1967;
537 double t1975, t1977, t1987, t1990, t2003, t2049, t2061, t2072;
538 double t2082, t2107, t2108, t2113, t2118, t2137, t2140, t2142;
539 double t2145, t2174, t2177, t2181, t2199, t2202, t2212, t2215;
540 double t2218, t2222, t2239, t2250, t2251, t2253, t2255, t2257;
541 double t2258, t2270, t2274, t2281, t2284, t2285, t2287, t2300;
542 double t2303, t2316, t2328, t2345, t2353, t2362, t2371, t2447;
543 double t2448, t2450, t2488, t2489, t2501, t2516, t2526, t2578;
544 double t2603;
545 #endif
546
547 #endif
548
549 #endif
550
551 #endif
552
553 #endif
554
555
556 gga_x_optx_params *params;
557
558 assert(p->params != NULL);
559 params = (gga_x_optx_params * )(p->params);
560
561 t1 = M_CBRT3;
562 t3 = POW_1_3(0.1e1 / M_PI);
563 t4 = t1 * t3;
564 t5 = M_CBRT4;
565 t6 = t5 * t5;
566 t7 = t4 * t6;
567 t8 = rho[0] - rho[1];
568 t9 = rho[0] + rho[1];
569 t10 = 0.1e1 / t9;
570 t11 = t8 * t10;
571 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
572 t14 = POW_1_3(t13);
573 t15 = t14 * t13;
574 t16 = POW_1_3(t9);
575 t17 = t15 * t16;
576 t18 = params->gamma * params->gamma;
577 t19 = params->b * t18;
578 t20 = sigma[0] * sigma[0];
579 t21 = rho[0] * rho[0];
580 t22 = t21 * t21;
581 t23 = t22 * rho[0];
582 t24 = POW_1_3(rho[0]);
583 t26 = 0.1e1 / t24 / t23;
584 t29 = t24 * t24;
585 t33 = 0.1e1 + params->gamma * sigma[0] / t29 / t21;
586 t34 = t33 * t33;
587 t35 = 0.1e1 / t34;
588 t38 = t19 * t20 * t26 * t35 + params->a;
589 t40 = t7 * t17 * t38;
590 t42 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
591 t43 = POW_1_3(t42);
592 t44 = t43 * t42;
593 t45 = t44 * t16;
594 t46 = sigma[2] * sigma[2];
595 t47 = rho[1] * rho[1];
596 t48 = t47 * t47;
597 t49 = t48 * rho[1];
598 t50 = POW_1_3(rho[1]);
599 t52 = 0.1e1 / t50 / t49;
600 t55 = t50 * t50;
601 t59 = 0.1e1 + params->gamma * sigma[2] / t55 / t47;
602 t60 = t59 * t59;
603 t61 = 0.1e1 / t60;
604 t64 = t19 * t46 * t52 * t61 + params->a;
605 t66 = t7 * t45 * t64;
606 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
607 zk[0] = -0.3e1 / 0.8e1 * t40 - 0.3e1 / 0.8e1 * t66;
608
609 #ifndef XC_DONT_COMPILE_VXC
610
611 if(order < 1) return;
612
613
614 t68 = 0.3e1 / 0.8e1 * t40;
615 t69 = 0.3e1 / 0.8e1 * t66;
616 t70 = t14 * t16;
617 t71 = t9 * t9;
618 t72 = 0.1e1 / t71;
619 t73 = t8 * t72;
620 t75 = t10 / 0.2e1 - t73 / 0.2e1;
621 t76 = t38 * t75;
622 t78 = t7 * t70 * t76;
623 t79 = t78 / 0.2e1;
624 t80 = t16 * t16;
625 t81 = 0.1e1 / t80;
626 t82 = t15 * t81;
627 t84 = t7 * t82 * t38;
628 t85 = t84 / 0.8e1;
629 t86 = t22 * t21;
630 t88 = 0.1e1 / t24 / t86;
631 t93 = params->b * t18 * params->gamma;
632 t94 = t20 * sigma[0];
633 t95 = t22 * t22;
634 t96 = t95 * rho[0];
635 t97 = 0.1e1 / t96;
636 t100 = 0.1e1 / t34 / t33;
637 t104 = 0.16e2 / 0.3e1 * t93 * t94 * t97 * t100 - 0.16e2 / 0.3e1 * t19 * t20 * t88 * t35;
638 t106 = t7 * t17 * t104;
639 t107 = 0.3e1 / 0.8e1 * t106;
640 t108 = t43 * t16;
641 t109 = -t75;
642 t110 = t64 * t109;
643 t112 = t7 * t108 * t110;
644 t113 = t112 / 0.2e1;
645 t114 = t44 * t81;
646 t116 = t7 * t114 * t64;
647 t117 = t116 / 0.8e1;
648 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
649 vrho[0] = -t68 - t69 + t9 * (-t79 - t85 - t107 - t113 - t117);
650
651 t121 = -t10 / 0.2e1 - t73 / 0.2e1;
652 t122 = t38 * t121;
653 t124 = t7 * t70 * t122;
654 t125 = t124 / 0.2e1;
655 t126 = -t121;
656 t127 = t64 * t126;
657 t129 = t7 * t108 * t127;
658 t130 = t129 / 0.2e1;
659 t131 = t48 * t47;
660 t133 = 0.1e1 / t50 / t131;
661 t137 = t46 * sigma[2];
662 t138 = t48 * t48;
663 t139 = t138 * rho[1];
664 t140 = 0.1e1 / t139;
665 t143 = 0.1e1 / t60 / t59;
666 t147 = -0.16e2 / 0.3e1 * t19 * t46 * t133 * t61 + 0.16e2 / 0.3e1 * t93 * t137 * t140 * t143;
667 t149 = t7 * t45 * t147;
668 t150 = 0.3e1 / 0.8e1 * t149;
669 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
670 vrho[1] = -t68 - t69 + t9 * (-t125 - t85 - t130 - t117 - t150);
671
672 t154 = t16 * t9 * t1;
673 t155 = t154 * t3;
674 t156 = t6 * t15;
675 t160 = 0.1e1 / t95;
676 t165 = -0.2e1 * t93 * t20 * t160 * t100 + 0.2e1 * t19 * sigma[0] * t26 * t35;
677 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
678 vsigma[0] = -0.3e1 / 0.8e1 * t155 * t156 * t165;
679
680 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
681 vsigma[1] = 0.0e0;
682
683 t169 = t6 * t44;
684 t173 = 0.1e1 / t138;
685 t178 = -0.2e1 * t93 * t46 * t173 * t143 + 0.2e1 * t19 * sigma[2] * t52 * t61;
686 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
687 vsigma[2] = -0.3e1 / 0.8e1 * t155 * t169 * t178;
688
689 #ifndef XC_DONT_COMPILE_FXC
690
691 if(order < 2) return;
692
693
694 t182 = t84 / 0.4e1;
695 t184 = t116 / 0.4e1;
696 t185 = t14 * t14;
697 t186 = 0.1e1 / t185;
698 t187 = t186 * t16;
699 t188 = t75 * t75;
700 t189 = t38 * t188;
701 t191 = t7 * t187 * t189;
702 t192 = t191 / 0.6e1;
703 t193 = t14 * t81;
704 t195 = t7 * t193 * t76;
705 t196 = t195 / 0.3e1;
706 t197 = t104 * t75;
707 t199 = t7 * t70 * t197;
708 t200 = t71 * t9;
709 t201 = 0.1e1 / t200;
710 t202 = t8 * t201;
711 t203 = -t72 + t202;
712 t204 = t38 * t203;
713 t206 = t7 * t70 * t204;
714 t207 = t206 / 0.2e1;
715 t209 = 0.1e1 / t80 / t9;
716 t210 = t15 * t209;
717 t212 = t7 * t210 * t38;
718 t213 = t212 / 0.12e2;
719 t215 = t7 * t82 * t104;
720 t216 = t215 / 0.4e1;
721 t217 = t21 * rho[0];
722 t218 = t22 * t217;
723 t220 = 0.1e1 / t24 / t218;
724 t225 = t95 * t21;
725 t226 = 0.1e1 / t225;
726 t231 = t18 * t18;
727 t232 = params->b * t231;
728 t233 = t20 * t20;
729 t234 = t95 * t22;
730 t236 = 0.1e1 / t29 / t234;
731 t238 = t34 * t34;
732 t239 = 0.1e1 / t238;
733 t243 = 0.304e3 / 0.9e1 * t19 * t20 * t220 * t35 - 0.688e3 / 0.9e1 * t93 * t94 * t226 * t100 + 0.128e3 / 0.3e1 * t232 * t233 * t236 * t239;
734 t245 = t7 * t17 * t243;
735 t246 = 0.3e1 / 0.8e1 * t245;
736 t247 = t43 * t43;
737 t248 = 0.1e1 / t247;
738 t249 = t248 * t16;
739 t250 = t109 * t109;
740 t251 = t64 * t250;
741 t253 = t7 * t249 * t251;
742 t254 = t253 / 0.6e1;
743 t255 = t43 * t81;
744 t257 = t7 * t255 * t110;
745 t258 = t257 / 0.3e1;
746 t259 = -t203;
747 t260 = t64 * t259;
748 t262 = t7 * t108 * t260;
749 t263 = t262 / 0.2e1;
750 t264 = t44 * t209;
751 t266 = t7 * t264 * t64;
752 t267 = t266 / 0.12e2;
753 t268 = -t192 - t196 - t199 - t207 + t213 - t216 - t246 - t254 - t258 - t263 + t267;
754 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
755 v2rho2[0] = -t78 - t182 - 0.3e1 / 0.4e1 * t106 - t112 - t184 + t9 * t268;
756
757 t270 = t6 * t186;
758 t271 = t4 * t270;
759 t272 = t16 * t38;
760 t273 = t121 * t75;
761 t275 = t271 * t272 * t273;
762 t278 = t7 * t193 * t122;
763 t280 = t104 * t121;
764 t282 = t7 * t70 * t280;
765 t285 = 0.1e1 / t80 / t71;
766 t286 = t14 * t285;
767 t287 = t38 * t8;
768 t289 = t7 * t286 * t287;
769 t293 = t6 * t248;
770 t294 = t4 * t293;
771 t295 = t16 * t64;
772 t296 = t126 * t109;
773 t298 = t294 * t295 * t296;
774 t301 = t7 * t255 * t127;
775 t303 = t43 * t285;
776 t304 = t64 * t8;
777 t306 = t7 * t303 * t304;
778 t309 = t147 * t109;
779 t311 = t7 * t108 * t309;
780 t314 = t7 * t114 * t147;
781 t316 = -t275 / 0.6e1 - t278 / 0.6e1 - t282 / 0.2e1 - t289 / 0.2e1 - t195 / 0.6e1 + t213 - t215 / 0.8e1 - t298 / 0.6e1 - t301 / 0.6e1 + t306 / 0.2e1 - t257 / 0.6e1 + t267 - t311 / 0.2e1 - t314 / 0.8e1;
782 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
783 v2rho2[1] = t9 * t316 - t107 - t113 - t125 - t130 - t150 - t182 - t184 - t79;
784
785 t319 = t121 * t121;
786 t320 = t38 * t319;
787 t322 = t7 * t187 * t320;
788 t323 = t322 / 0.6e1;
789 t324 = t278 / 0.3e1;
790 t325 = t72 + t202;
791 t326 = t38 * t325;
792 t328 = t7 * t70 * t326;
793 t329 = t328 / 0.2e1;
794 t330 = t126 * t126;
795 t331 = t64 * t330;
796 t333 = t7 * t249 * t331;
797 t334 = t333 / 0.6e1;
798 t335 = t301 / 0.3e1;
799 t336 = t147 * t126;
800 t338 = t7 * t108 * t336;
801 t339 = -t325;
802 t340 = t64 * t339;
803 t342 = t7 * t108 * t340;
804 t343 = t342 / 0.2e1;
805 t344 = t314 / 0.4e1;
806 t345 = t47 * rho[1];
807 t346 = t48 * t345;
808 t348 = 0.1e1 / t50 / t346;
809 t353 = t138 * t47;
810 t354 = 0.1e1 / t353;
811 t359 = t46 * t46;
812 t360 = t138 * t48;
813 t362 = 0.1e1 / t55 / t360;
814 t364 = t60 * t60;
815 t365 = 0.1e1 / t364;
816 t369 = 0.304e3 / 0.9e1 * t19 * t46 * t348 * t61 - 0.688e3 / 0.9e1 * t93 * t137 * t354 * t143 + 0.128e3 / 0.3e1 * t232 * t359 * t362 * t365;
817 t371 = t7 * t45 * t369;
818 t372 = 0.3e1 / 0.8e1 * t371;
819 t373 = -t323 - t324 - t329 + t213 - t334 - t335 - t338 - t343 + t267 - t344 - t372;
820 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
821 v2rho2[2] = -t124 - t182 - t129 - t184 - 0.3e1 / 0.4e1 * t149 + t9 * t373;
822
823 t376 = t7 * t17 * t165;
824 t378 = t6 * t14;
825 t379 = t165 * t75;
826 t391 = t95 * t217;
827 t393 = 0.1e1 / t29 / t391;
828 t398 = -0.32e2 / 0.3e1 * t19 * sigma[0] * t88 * t35 + 0.80e2 / 0.3e1 * t93 * t20 * t97 * t100 - 0.16e2 * t232 * t94 * t393 * t239;
829 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
830 v2rhosigma[0] = -t376 / 0.2e1 - t155 * t378 * t379 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * t398;
831
832 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
833 v2rhosigma[1] = 0.0e0;
834
835 t403 = t7 * t45 * t178;
836 t404 = t6 * t43;
837 t405 = t178 * t109;
838 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
839 v2rhosigma[2] = -t155 * t404 * t405 / 0.2e1 - t403 / 0.2e1;
840
841 t410 = t378 * t165 * t121;
842 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
843 v2rhosigma[3] = -t155 * t410 / 0.2e1 - t376 / 0.2e1;
844
845 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
846 v2rhosigma[4] = 0.0e0;
847
848 t415 = t404 * t178 * t126;
849 t426 = t138 * t345;
850 t428 = 0.1e1 / t55 / t426;
851 t433 = -0.32e2 / 0.3e1 * t19 * sigma[2] * t133 * t61 + 0.80e2 / 0.3e1 * t93 * t46 * t140 * t143 - 0.16e2 * t232 * t137 * t428 * t365;
852 t434 = t169 * t433;
853 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
854 v2rhosigma[5] = -t403 / 0.2e1 - t155 * t415 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t434;
855
856 t445 = 0.1e1 / t29 / t225;
857 t450 = -0.8e1 * t93 * sigma[0] * t160 * t100 + 0.6e1 * t232 * t20 * t445 * t239 + 0.2e1 * t19 * t26 * t35;
858 t451 = t156 * t450;
859 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
860 v2sigma2[0] = -0.3e1 / 0.8e1 * t155 * t451;
861
862 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
863 v2sigma2[1] = 0.0e0;
864
865 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
866 v2sigma2[2] = 0.0e0;
867
868 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
869 v2sigma2[3] = 0.0e0;
870
871 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
872 v2sigma2[4] = 0.0e0;
873
874 t462 = 0.1e1 / t55 / t353;
875 t467 = -0.8e1 * t93 * sigma[2] * t173 * t143 + 0.6e1 * t232 * t46 * t462 * t365 + 0.2e1 * t19 * t52 * t61;
876 t468 = t169 * t467;
877 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
878 v2sigma2[5] = -0.3e1 / 0.8e1 * t155 * t468;
879
880 #ifndef XC_DONT_COMPILE_KXC
881
882 if(order < 3) return;
883
884
885 t475 = t212 / 0.4e1;
886 t477 = t266 / 0.4e1;
887 t481 = 0.1e1 / t185 / t13;
888 t482 = t481 * t16;
889 t483 = t188 * t75;
890 t484 = t38 * t483;
891 t486 = t7 * t482 * t484;
892 t487 = t486 / 0.9e1;
893 t488 = t186 * t81;
894 t490 = t7 * t488 * t189;
895 t492 = t75 * t203;
896 t494 = t271 * t272 * t492;
897 t495 = t494 / 0.2e1;
898 t496 = t14 * t209;
899 t498 = t7 * t496 * t76;
900 t499 = t498 / 0.3e1;
901 t501 = 0.1e1 / t247 / t42;
902 t502 = t501 * t16;
903 t503 = t250 * t109;
904 t504 = t64 * t503;
905 t506 = t7 * t502 * t504;
906 t507 = t506 / 0.9e1;
907 t508 = t248 * t81;
908 t510 = t7 * t508 * t251;
909 t512 = t109 * t259;
910 t514 = t294 * t295 * t512;
911 t515 = t514 / 0.2e1;
912 t516 = t43 * t209;
913 t518 = t7 * t516 * t110;
914 t519 = t518 / 0.3e1;
915 t520 = t104 * t188;
916 t522 = t7 * t187 * t520;
917 t523 = t522 / 0.2e1;
918 t525 = t7 * t193 * t197;
919 t528 = t7 * t193 * t204;
920 t530 = t243 * t75;
921 t532 = t7 * t70 * t530;
922 t533 = 0.3e1 / 0.2e1 * t532;
923 t534 = t104 * t203;
924 t536 = t7 * t70 * t534;
925 t537 = 0.3e1 / 0.2e1 * t536;
926 t538 = t71 * t71;
927 t539 = 0.1e1 / t538;
928 t540 = t8 * t539;
929 t542 = 0.3e1 * t201 - 0.3e1 * t540;
930 t543 = t38 * t542;
931 t545 = t7 * t70 * t543;
932 t546 = t545 / 0.2e1;
933 t547 = t15 * t285;
934 t549 = t7 * t547 * t38;
935 t550 = 0.5e1 / 0.36e2 * t549;
936 t552 = t7 * t255 * t260;
937 t554 = -t542;
938 t555 = t64 * t554;
939 t557 = t7 * t108 * t555;
940 t558 = t557 / 0.2e1;
941 t559 = t44 * t285;
942 t561 = t7 * t559 * t64;
943 t562 = 0.5e1 / 0.36e2 * t561;
944 t564 = t7 * t210 * t104;
945 t565 = t564 / 0.4e1;
946 t567 = t7 * t82 * t243;
947 t570 = 0.1e1 / t24 / t95;
948 t575 = 0.1e1 / t391;
949 t580 = t95 * t23;
950 t582 = 0.1e1 / t29 / t580;
951 t588 = params->b * t231 * params->gamma;
952 t589 = t233 * sigma[0];
953 t590 = t95 * t95;
954 t592 = 0.1e1 / t24 / t590;
955 t595 = 0.1e1 / t238 / t33;
956 t599 = -0.6688e4 / 0.27e2 * t19 * t20 * t570 * t35 + 0.25504e5 / 0.27e2 * t93 * t94 * t575 * t100 - 0.1152e4 * t232 * t233 * t582 * t239 + 0.4096e4 / 0.9e1 * t588 * t589 * t592 * t595;
957 t601 = t7 * t17 * t599;
958 t602 = 0.3e1 / 0.8e1 * t601;
959 t603 = -t528 / 0.2e1 - t533 - t537 - t546 - t550 - t552 / 0.2e1 - t558 - t562 + t565 - 0.3e1 / 0.8e1 * t567 - t602;
960 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
961 v3rho3[0] = -t195 - t191 / 0.2e1 - t253 / 0.2e1 - t257 - 0.3e1 * t199 - 0.3e1 / 0.2e1 * t206 + t475 - 0.3e1 / 0.2e1 * t262 + t477 - 0.3e1 / 0.4e1 * t215 - 0.9e1 / 0.8e1 * t245 + t9 * (t487 - t490 / 0.6e1 - t495 + t499 + t507 - t510 / 0.6e1 - t515 + t519 - t523 - t525 + t603);
962
963 t608 = t275 / 0.3e1;
964 t609 = t298 / 0.3e1;
965 t612 = t6 * t481;
966 t613 = t4 * t612;
967 t614 = t121 * t188;
968 t616 = t613 * t272 * t614;
969 t618 = t81 * t38;
970 t620 = t271 * t618 * t273;
971 t621 = t620 / 0.9e1;
972 t622 = t285 * t38;
973 t623 = t8 * t75;
974 t625 = t271 * t622 * t623;
975 t627 = t6 * t501;
976 t628 = t4 * t627;
977 t629 = t126 * t250;
978 t631 = t628 * t295 * t629;
979 t633 = t81 * t64;
980 t635 = t294 * t633 * t296;
981 t636 = t635 / 0.9e1;
982 t637 = t16 * t104;
983 t639 = t271 * t637 * t273;
984 t641 = t121 * t203;
985 t643 = t271 * t272 * t641;
986 t645 = t126 * t259;
987 t647 = t294 * t295 * t645;
988 t649 = t285 * t64;
989 t650 = t8 * t109;
990 t652 = t294 * t649 * t650;
991 t661 = t616 / 0.9e1 - t621 - t625 / 0.3e1 + t631 / 0.9e1 - t636 - t639 / 0.3e1 - t643 / 0.6e1 - t647 / 0.6e1 + t652 / 0.3e1 - t490 / 0.18e2 + 0.2e1 / 0.9e1 * t498 - t510 / 0.18e2 + 0.2e1 / 0.9e1 * t518 - t525 / 0.3e1 - t528 / 0.6e1 - t552 / 0.6e1;
992 t663 = t7 * t496 * t122;
993 t666 = 0.1e1 / t80 / t200;
994 t667 = t14 * t666;
995 t669 = t7 * t667 * t287;
996 t671 = t104 * t8;
997 t673 = t7 * t286 * t671;
998 t675 = t7 * t516 * t127;
999 t677 = t147 * t250;
1000 t679 = t7 * t249 * t677;
1001 t682 = t7 * t255 * t309;
1002 t683 = t682 / 0.3e1;
1003 t685 = t7 * t193 * t280;
1004 t686 = t685 / 0.3e1;
1005 t687 = t243 * t121;
1006 t689 = t7 * t70 * t687;
1007 t691 = t43 * t666;
1008 t693 = t7 * t691 * t304;
1009 t695 = t147 * t259;
1010 t697 = t7 * t108 * t695;
1011 t702 = t7 * t286 * t38;
1012 t705 = t7 * t303 * t64;
1013 t708 = t7 * t264 * t147;
1014 t710 = t663 / 0.9e1 + 0.7e1 / 0.6e1 * t669 - t673 + t675 / 0.9e1 - t679 / 0.6e1 - t683 - t686 - t689 / 0.2e1 - 0.7e1 / 0.6e1 * t693 - t697 / 0.2e1 - t550 - t562 + t564 / 0.6e1 - t567 / 0.8e1 - t702 / 0.2e1 + t705 / 0.2e1 + t708 / 0.12e2;
1015 t713 = -t207 + t475 - t263 + t477 - t282 + t306 - t311 - t344 - t215 / 0.2e1 - t246 + t9 * (t661 + t710);
1016 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1017 v3rho3[1] = -0.2e1 / 0.3e1 * t195 - t192 - t254 - 0.2e1 / 0.3e1 * t257 - t608 - t324 - t289 - t609 - t335 - t199 + t713;
1018
1019 t718 = t319 * t75;
1020 t720 = t613 * t272 * t718;
1021 t722 = t325 * t75;
1022 t724 = t271 * t272 * t722;
1023 t726 = t330 * t109;
1024 t728 = t628 * t295 * t726;
1025 t730 = t126 * t8;
1026 t732 = t294 * t649 * t730;
1027 t734 = t16 * t147;
1028 t736 = t294 * t734 * t296;
1029 t738 = t339 * t109;
1030 t740 = t294 * t295 * t738;
1031 t742 = t121 * t8;
1032 t744 = t271 * t622 * t742;
1033 t751 = t720 / 0.9e1 - t724 / 0.6e1 + t728 / 0.9e1 + t732 / 0.3e1 - t736 / 0.3e1 - t740 / 0.6e1 - t744 / 0.3e1 - t621 - t636 + t498 / 0.9e1 + t518 / 0.9e1 + 0.2e1 / 0.9e1 * t663 - t669 / 0.3e1 + 0.2e1 / 0.9e1 * t675 - t683 - t686;
1034 t754 = t7 * t488 * t320;
1035 t757 = t7 * t193 * t326;
1036 t760 = t7 * t508 * t331;
1037 t763 = t7 * t255 * t336;
1038 t766 = t7 * t255 * t340;
1039 t768 = t104 * t319;
1040 t770 = t7 * t187 * t768;
1041 t772 = t104 * t325;
1042 t774 = t7 * t70 * t772;
1043 t777 = -t201 - 0.3e1 * t540;
1044 t778 = t38 * t777;
1045 t780 = t7 * t70 * t778;
1046 t782 = t147 * t8;
1047 t784 = t7 * t303 * t782;
1048 t785 = -t777;
1049 t786 = t64 * t785;
1050 t788 = t7 * t108 * t786;
1051 t790 = t369 * t109;
1052 t792 = t7 * t108 * t790;
1053 t797 = t7 * t114 * t369;
1054 t799 = t693 / 0.3e1 - t754 / 0.18e2 - t757 / 0.6e1 - t760 / 0.18e2 - t763 / 0.3e1 - t766 / 0.6e1 - t770 / 0.6e1 - t774 / 0.2e1 - t780 / 0.2e1 + t784 - t788 / 0.2e1 - t792 / 0.2e1 - t550 - t562 + t564 / 0.12e2 + t708 / 0.6e1 - t797 / 0.8e1;
1055 t802 = t477 - t282 + t306 - t311 - t314 / 0.2e1 - t329 - t338 - t343 - t216 - t372 + t9 * (t751 + t799);
1056 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1057 v3rho3[2] = -t196 - t258 - t608 - 0.2e1 / 0.3e1 * t278 - t289 - t609 - 0.2e1 / 0.3e1 * t301 - t323 - t334 + t475 + t802;
1058
1059 t810 = t663 / 0.3e1;
1060 t811 = t675 / 0.3e1;
1061 t816 = t319 * t121;
1062 t817 = t38 * t816;
1063 t819 = t7 * t482 * t817;
1064 t820 = t819 / 0.9e1;
1065 t821 = t121 * t325;
1066 t823 = t271 * t272 * t821;
1067 t824 = t823 / 0.2e1;
1068 t825 = t330 * t126;
1069 t826 = t64 * t825;
1070 t828 = t7 * t502 * t826;
1071 t829 = t828 / 0.9e1;
1072 t831 = t126 * t339;
1073 t833 = t294 * t295 * t831;
1074 t834 = t833 / 0.2e1;
1075 t835 = t708 / 0.4e1;
1076 t838 = -0.3e1 * t201 - 0.3e1 * t540;
1077 t839 = t38 * t838;
1078 t841 = t7 * t70 * t839;
1079 t842 = t841 / 0.2e1;
1080 t843 = t147 * t330;
1081 t845 = t7 * t249 * t843;
1082 t846 = t845 / 0.2e1;
1083 t847 = t369 * t126;
1084 t849 = t7 * t108 * t847;
1085 t850 = 0.3e1 / 0.2e1 * t849;
1086 t851 = t147 * t339;
1087 t853 = t7 * t108 * t851;
1088 t854 = 0.3e1 / 0.2e1 * t853;
1089 t855 = -t838;
1090 t856 = t64 * t855;
1091 t858 = t7 * t108 * t856;
1092 t859 = t858 / 0.2e1;
1093 t861 = 0.1e1 / t50 / t138;
1094 t866 = 0.1e1 / t426;
1095 t871 = t138 * t49;
1096 t873 = 0.1e1 / t55 / t871;
1097 t878 = t359 * sigma[2];
1098 t879 = t138 * t138;
1099 t881 = 0.1e1 / t50 / t879;
1100 t884 = 0.1e1 / t364 / t59;
1101 t888 = -0.6688e4 / 0.27e2 * t19 * t46 * t861 * t61 + 0.25504e5 / 0.27e2 * t93 * t137 * t866 * t143 - 0.1152e4 * t232 * t359 * t873 * t365 + 0.4096e4 / 0.9e1 * t588 * t878 * t881 * t884;
1102 t890 = t7 * t45 * t888;
1103 t891 = 0.3e1 / 0.8e1 * t890;
1104 t892 = -t834 - t550 - t562 + t835 - 0.3e1 / 0.8e1 * t797 - t842 - t846 - t850 - t854 - t859 - t891;
1105 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1106 v3rho3[3] = -t278 - t301 - t322 / 0.2e1 - t333 / 0.2e1 + t475 + t477 - 0.3e1 / 0.4e1 * t314 - 0.3e1 / 0.2e1 * t328 - 0.3e1 * t338 - 0.3e1 / 0.2e1 * t342 - 0.9e1 / 0.8e1 * t371 + t9 * (t810 + t811 - t754 / 0.6e1 - t757 / 0.2e1 - t760 / 0.6e1 - t763 - t766 / 0.2e1 + t820 - t824 + t829 + t892);
1107
1108 t896 = t7 * t70 * t379;
1109 t900 = t7 * t82 * t165 / 0.6e1;
1110 t902 = t7 * t17 * t398;
1111 t903 = t165 * t188;
1112 t907 = t398 * t75;
1113 t910 = t165 * t203;
1114 t928 = 0.1e1 / t24 / t95 / t218;
1115 t933 = 0.608e3 / 0.9e1 * t19 * sigma[0] * t220 * t35 - 0.2672e4 / 0.9e1 * t93 * t20 * t226 * t100 + 0.400e3 * t232 * t94 * t236 * t239 - 0.512e3 / 0.3e1 * t588 * t233 * t928 * t595;
1116 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1117 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t896 - t900 - t902 - t155 * t270 * t903 / 0.6e1 - t155 * t378 * t907 - t155 * t378 * t910 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * t933;
1118
1119 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1120 v3rho2sigma[1] = 0.0e0;
1121
1122 t938 = t7 * t108 * t405;
1123 t942 = t7 * t114 * t178 / 0.6e1;
1124 t943 = t178 * t250;
1125 t947 = t178 * t259;
1126 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1127 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t938 - t942 - t155 * t293 * t943 / 0.6e1 - t155 * t404 * t947 / 0.2e1;
1128
1129 t953 = t16 * t1;
1130 t954 = t953 * t3;
1131 t955 = t954 * t410;
1132 t957 = t3 * t6;
1133 t958 = t154 * t957;
1134 t959 = t186 * t165;
1135 t960 = t959 * t273;
1136 t964 = t378 * t398 * t121;
1137 t967 = t209 * t1;
1138 t968 = t967 * t3;
1139 t970 = t378 * t165 * t8;
1140 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1141 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t896 - t900 - t902 / 0.2e1 - 0.2e1 / 0.3e1 * t955 - t958 * t960 / 0.6e1 - t155 * t964 / 0.2e1 - t968 * t970 / 0.2e1;
1142
1143 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1144 v3rho2sigma[4] = 0.0e0;
1145
1146 t974 = t954 * t415;
1147 t976 = t248 * t178;
1148 t977 = t976 * t296;
1149 t981 = t404 * t178 * t8;
1150 t984 = t954 * t434;
1151 t987 = t404 * t433 * t109;
1152 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1153 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t938 - t942 - 0.2e1 / 0.3e1 * t974 - t958 * t977 / 0.6e1 + t968 * t981 / 0.2e1 - t984 / 0.2e1 - t155 * t987 / 0.2e1;
1154
1155 t992 = t270 * t165 * t319;
1156 t996 = t378 * t165 * t325;
1157 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1158 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t955 - t900 - t155 * t992 / 0.6e1 - t155 * t996 / 0.2e1;
1159
1160 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1161 v3rho2sigma[7] = 0.0e0;
1162
1163 t1001 = t293 * t178 * t330;
1164 t1005 = t404 * t433 * t126;
1165 t1008 = t404 * t178 * t339;
1166 t1025 = 0.1e1 / t50 / t138 / t346;
1167 t1030 = 0.608e3 / 0.9e1 * t19 * sigma[2] * t348 * t61 - 0.2672e4 / 0.9e1 * t93 * t46 * t354 * t143 + 0.400e3 * t232 * t137 * t362 * t365 - 0.512e3 / 0.3e1 * t588 * t359 * t1025 * t884;
1168 t1031 = t169 * t1030;
1169 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1170 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t974 - t942 - t984 - t155 * t1001 / 0.6e1 - t155 * t1005 - t155 * t1008 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1031;
1171
1172 t1034 = t954 * t451;
1173 t1037 = t378 * t450 * t75;
1174 t1043 = t97 * t100;
1175 t1051 = t95 * t86;
1176 t1053 = 0.1e1 / t24 / t1051;
1177 t1058 = -0.32e2 / 0.3e1 * t19 * t88 * t35 + 0.224e3 / 0.3e1 * t93 * t1043 * sigma[0] - 0.128e3 * t232 * t20 * t393 * t239 + 0.64e2 * t588 * t94 * t1053 * t595;
1178 t1059 = t156 * t1058;
1179 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1180 v3rhosigma2[0] = -t1034 / 0.2e1 - t155 * t1037 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1059;
1181
1182 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1183 v3rhosigma2[1] = 0.0e0;
1184
1185 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1186 v3rhosigma2[2] = 0.0e0;
1187
1188 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1189 v3rhosigma2[3] = 0.0e0;
1190
1191 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1192 v3rhosigma2[4] = 0.0e0;
1193
1194 t1062 = t954 * t468;
1195 t1064 = t404 * t467 * t109;
1196 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1197 v3rhosigma2[5] = -t155 * t1064 / 0.2e1 - t1062 / 0.2e1;
1198
1199 t1068 = t378 * t450 * t121;
1200 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1201 v3rhosigma2[6] = -t155 * t1068 / 0.2e1 - t1034 / 0.2e1;
1202
1203 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1204 v3rhosigma2[7] = 0.0e0;
1205
1206 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1207 v3rhosigma2[8] = 0.0e0;
1208
1209 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1210 v3rhosigma2[9] = 0.0e0;
1211
1212 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1213 v3rhosigma2[10] = 0.0e0;
1214
1215 t1073 = t404 * t467 * t126;
1216 t1079 = t140 * t143;
1217 t1087 = t138 * t131;
1218 t1089 = 0.1e1 / t50 / t1087;
1219 t1094 = -0.32e2 / 0.3e1 * t19 * t133 * t61 + 0.224e3 / 0.3e1 * t93 * t1079 * sigma[2] - 0.128e3 * t232 * t46 * t428 * t365 + 0.64e2 * t588 * t137 * t1089 * t884;
1220 t1095 = t169 * t1094;
1221 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1222 v3rhosigma2[11] = -t1062 / 0.2e1 - t155 * t1073 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t1095;
1223
1224 t1106 = 0.1e1 / t24 / t580;
1225 t1111 = -0.24e2 * t588 * t20 * t1106 * t595 + 0.36e2 * t232 * sigma[0] * t445 * t239 - 0.12e2 * t93 * t160 * t100;
1226 t1112 = t156 * t1111;
1227 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1228 v3sigma3[0] = -0.3e1 / 0.8e1 * t155 * t1112;
1229
1230 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1231 v3sigma3[1] = 0.0e0;
1232
1233 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1234 v3sigma3[2] = 0.0e0;
1235
1236 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1237 v3sigma3[3] = 0.0e0;
1238
1239 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1240 v3sigma3[4] = 0.0e0;
1241
1242 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1243 v3sigma3[5] = 0.0e0;
1244
1245 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1246 v3sigma3[6] = 0.0e0;
1247
1248 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1249 v3sigma3[7] = 0.0e0;
1250
1251 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1252 v3sigma3[8] = 0.0e0;
1253
1254 t1123 = 0.1e1 / t50 / t871;
1255 t1128 = -0.24e2 * t588 * t46 * t1123 * t884 + 0.36e2 * t232 * sigma[2] * t462 * t365 - 0.12e2 * t93 * t173 * t143;
1256 t1129 = t169 * t1128;
1257 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1258 v3sigma3[9] = -0.3e1 / 0.8e1 * t155 * t1129;
1259
1260 #ifndef XC_DONT_COMPILE_LXC
1261
1262 if(order < 4) return;
1263
1264
1265 t1143 = 0.4e1 / 0.9e1 * t486 - 0.2e1 / 0.3e1 * t490 - 0.2e1 * t494 + 0.4e1 / 0.3e1 * t498 + 0.4e1 / 0.9e1 * t506 - 0.2e1 / 0.3e1 * t510 - 0.2e1 * t514 + 0.4e1 / 0.3e1 * t518 - 0.2e1 * t522 - 0.4e1 * t525 - 0.2e1 * t528;
1266 t1147 = 0.5e1 / 0.9e1 * t549;
1267 t1150 = 0.5e1 / 0.9e1 * t561;
1268 t1153 = t259 * t259;
1269 t1159 = t7 * t303 * t110;
1270 t1166 = t7 * t488 * t520;
1271 t1169 = t7 * t496 * t197;
1272 t1172 = t7 * t496 * t204;
1273 t1175 = t7 * t516 * t260;
1274 t1181 = t7 * t193 * t530;
1275 t1184 = t7 * t193 * t534;
1276 t1187 = t7 * t193 * t543;
1277 t1201 = t538 * t9;
1278 t1203 = t8 / t1201;
1279 t1205 = -0.12e2 * t539 + 0.12e2 * t1203;
1280 t1211 = t7 * t255 * t555;
1281 t1218 = t13 * t13;
1282 t1220 = 0.1e1 / t185 / t1218;
1283 t1221 = t1220 * t16;
1284 t1222 = t188 * t188;
1285 t1227 = t481 * t81;
1286 t1229 = t7 * t1227 * t484;
1287 t1231 = -t7 * t249 * t64 * t1153 / 0.2e1 - 0.20e2 / 0.27e2 * t1159 + 0.4e1 / 0.9e1 * t7 * t482 * t104 * t483 - 0.2e1 / 0.3e1 * t1166 + 0.4e1 / 0.3e1 * t1169 + 0.2e1 / 0.3e1 * t1172 + 0.2e1 / 0.3e1 * t1175 - t7 * t187 * t243 * t188 - 0.2e1 * t1181 - 0.2e1 * t1184 - 0.2e1 / 0.3e1 * t1187 - 0.2e1 * t7 * t70 * t599 * t75 - 0.3e1 * t7 * t70 * t243 * t203 - 0.2e1 * t7 * t70 * t104 * t542 - t7 * t70 * t38 * t1205 / 0.2e1 - 0.2e1 / 0.3e1 * t1211 + t7 * t108 * t64 * t1205 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1221 * t38 * t1222 + 0.4e1 / 0.27e2 * t1229;
1288 t1232 = t186 * t209;
1289 t1234 = t7 * t1232 * t189;
1290 t1236 = t203 * t203;
1291 t1242 = t7 * t286 * t76;
1292 t1244 = t42 * t42;
1293 t1246 = 0.1e1 / t247 / t1244;
1294 t1247 = t1246 * t16;
1295 t1248 = t250 * t250;
1296 t1253 = t501 * t81;
1297 t1255 = t7 * t1253 * t504;
1298 t1257 = t248 * t209;
1299 t1259 = t7 * t1257 * t251;
1300 t1266 = t294 * t633 * t512;
1301 t1277 = t271 * t618 * t492;
1302 t1287 = t7 * t547 * t104;
1303 t1290 = t7 * t210 * t243;
1304 t1293 = t7 * t82 * t599;
1305 t1312 = t590 * rho[0];
1306 t1320 = params->b * t231 * t18;
1307 t1326 = 0.1e1 / t238 / t34;
1308 t1337 = 0.10e2 / 0.27e2 * t7 * t15 * t666 * t38;
1309 t1341 = 0.10e2 / 0.27e2 * t7 * t44 * t666 * t64;
1310 t1342 = 0.2e1 / 0.9e1 * t1234 - t7 * t187 * t38 * t1236 / 0.2e1 - 0.20e2 / 0.27e2 * t1242 - 0.5e1 / 0.27e2 * t7 * t1247 * t64 * t1248 + 0.4e1 / 0.27e2 * t1255 + 0.2e1 / 0.9e1 * t1259 + 0.2e1 / 0.3e1 * t628 * t295 * t250 * t259 - 0.2e1 / 0.3e1 * t1266 - 0.2e1 / 0.3e1 * t294 * t295 * t109 * t554 + 0.2e1 / 0.3e1 * t613 * t272 * t188 * t203 - 0.2e1 / 0.3e1 * t1277 - 0.2e1 * t271 * t637 * t492 - 0.2e1 / 0.3e1 * t271 * t272 * t75 * t542 - 0.5e1 / 0.9e1 * t1287 + t1290 / 0.2e1 - t1293 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * (0.167200e6 / 0.81e2 * t19 * t20 / t24 / t96 * t35 - 0.948640e6 / 0.81e2 * t93 * t94 / t234 * t100 + 0.629120e6 / 0.27e2 * t232 * t233 / t29 / t1051 * t239 - 0.532480e6 / 0.27e2 * t588 * t589 / t24 / t1312 * t595 + 0.163840e6 / 0.27e2 * t1320 * t233 * t20 / t590 / t22 * t1326) + t1337 + t1341;
1311 t1345 = -0.6e1 * t532 - 0.6e1 * t536 - 0.2e1 * t545 - t1147 - 0.2e1 * t552 - 0.2e1 * t557 - t1150 + t564 - 0.3e1 / 0.2e1 * t567 - 0.3e1 / 0.2e1 * t601 + t9 * (t1231 + t1342);
1312 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1313 v4rho4[0] = t1143 + t1345;
1314
1315 t1347 = t7 * t667 * t38;
1316 t1353 = t7 * t691 * t64;
1317 t1356 = t7 * t559 * t147;
1318 t1358 = t248 * t285;
1319 t1363 = t7 * t303 * t127;
1320 t1370 = t7 * t508 * t677;
1321 t1373 = t7 * t516 * t309;
1322 t1374 = t1373 / 0.3e1;
1323 t1376 = t7 * t496 * t280;
1324 t1377 = t1376 / 0.3e1;
1325 t1379 = t7 * t667 * t671;
1326 t1386 = 0.1e1 / t80 / t538;
1327 t1389 = t7 * t43 * t1386 * t304;
1328 t1392 = t7 * t255 * t695;
1329 t1395 = t7 * t193 * t687;
1330 t1401 = 0.5e1 / 0.2e1 * t1347 - 0.3e1 / 0.2e1 * t7 * t286 * t104 - 0.5e1 / 0.2e1 * t1353 - 0.5e1 / 0.36e2 * t1356 + t7 * t1358 * t110 / 0.2e1 - 0.5e1 / 0.27e2 * t1363 + t7 * t502 * t147 * t503 / 0.9e1 - t1370 / 0.6e1 + t1374 + t1377 + 0.7e1 / 0.2e1 * t1379 - 0.3e1 / 0.2e1 * t7 * t286 * t243 * t8 + 0.25e2 / 0.6e1 * t1389 - t1392 / 0.2e1 - t1395 / 0.2e1 - t7 * t70 * t599 * t121 / 0.2e1;
1331 t1402 = t186 * t285;
1332 t1411 = t7 * t286 * t122;
1333 t1415 = t7 * t14 * t1386 * t287;
1334 t1429 = -t7 * t1402 * t76 / 0.2e1 - t7 * t108 * t147 * t554 / 0.2e1 - 0.5e1 / 0.27e2 * t1411 - 0.25e2 / 0.6e1 * t1415 - 0.5e1 / 0.9e1 * t1159 - t1166 / 0.6e1 + 0.2e1 / 0.3e1 * t1169 + t1172 / 0.3e1 + t1175 / 0.3e1 - t1181 / 0.2e1 - t1184 / 0.2e1 - t1187 / 0.6e1 - t1211 / 0.6e1 + t1229 / 0.27e2 + t1234 / 0.9e1 - 0.5e1 / 0.9e1 * t1242;
1335 t1449 = t271 * t209 * t38 * t273;
1336 t1450 = t1449 / 0.9e1;
1337 t1451 = t666 * t38;
1338 t1453 = t271 * t1451 * t623;
1339 t1465 = t294 * t209 * t64 * t296;
1340 t1466 = t1465 / 0.9e1;
1341 t1472 = t271 * t81 * t104 * t273;
1342 t1474 = t285 * t104;
1343 t1478 = t271 * t618 * t641;
1344 t1484 = t1255 / 0.27e2 + t1259 / 0.9e1 - t294 * t734 * t512 / 0.2e1 - t271 * t16 * t243 * t273 / 0.2e1 - t271 * t637 * t641 / 0.2e1 - t271 * t272 * t121 * t542 / 0.6e1 + t1450 + 0.7e1 / 0.6e1 * t1453 - t294 * t295 * t126 * t554 / 0.6e1 + t294 * t649 * t8 * t259 / 0.2e1 + t1466 + t613 * t637 * t614 / 0.3e1 - t1472 / 0.3e1 - t271 * t1474 * t623 - t1478 / 0.6e1 - t271 * t622 * t8 * t203 / 0.2e1;
1345 t1486 = t294 * t633 * t645;
1346 t1492 = t666 * t64;
1347 t1494 = t294 * t1492 * t650;
1348 t1497 = t4 * t6 * t1220;
1349 t1503 = t613 * t618 * t614;
1350 t1510 = t4 * t6 * t1246;
1351 t1516 = t628 * t633 * t629;
1352 t1531 = -t1486 / 0.6e1 - t628 * t649 * t8 * t250 / 0.3e1 - 0.7e1 / 0.6e1 * t1494 - 0.5e1 / 0.27e2 * t1497 * t272 * t121 * t483 + t1503 / 0.9e1 + t613 * t622 * t8 * t188 / 0.3e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t126 * t503 + t1516 / 0.9e1 - t1266 / 0.6e1 - t1277 / 0.6e1 + t613 * t272 * t641 * t75 / 0.3e1 + t628 * t295 * t645 * t109 / 0.3e1 - 0.5e1 / 0.12e2 * t1287 + t1290 / 0.4e1 - t1293 / 0.8e1 + t1337 + t1341;
1353 t1536 = t620 / 0.3e1;
1354 t1538 = t635 / 0.3e1;
1355 t1541 = t9 * (t1401 + t1429 + t1484 + t1531) + t616 / 0.3e1 - t1536 - t625 + t631 / 0.3e1 - t1538 - t639 - t643 / 0.2e1 - t647 / 0.2e1 + t652 - t495;
1356 t1545 = -t515 + t487 - t490 / 0.3e1 + t498 + t507 - t510 / 0.3e1 + t518 - t523 - 0.2e1 * t525 - t528 - t533;
1357 t1550 = -t537 - t546 - t552 - t558 + t810 + 0.7e1 / 0.2e1 * t669 - 0.3e1 * t673 + t811 - t679 / 0.2e1 - t682 - t685;
1358 t1558 = -0.3e1 / 0.2e1 * t689 - 0.7e1 / 0.2e1 * t693 - 0.3e1 / 0.2e1 * t697 - t1147 - t1150 + 0.3e1 / 0.4e1 * t564 - 0.3e1 / 0.4e1 * t567 - t602 - 0.3e1 / 0.2e1 * t702 + 0.3e1 / 0.2e1 * t705 + t835;
1359 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1360 v4rho4[1] = t1541 + t1545 + t1550 + t1558;
1361
1362 t1583 = t294 * t1492 * t730;
1363 t1588 = t285 * t147;
1364 t1605 = t294 * t81 * t147 * t296;
1365 t1612 = t294 * t633 * t738;
1366 t1622 = t613 * t618 * t718;
1367 t1625 = t271 * t1451 * t742;
1368 t1631 = 0.2e1 / 0.9e1 * t613 * t637 * t718 + t613 * t272 * t319 * t203 / 0.9e1 - t271 * t637 * t722 / 0.3e1 - t271 * t272 * t777 * t75 / 0.3e1 - t271 * t272 * t325 * t203 / 0.6e1 + t628 * t295 * t330 * t259 / 0.9e1 - 0.7e1 / 0.9e1 * t1583 - t294 * t734 * t645 / 0.3e1 + 0.2e1 / 0.3e1 * t294 * t1588 * t650 - t294 * t295 * t785 * t109 / 0.3e1 - t294 * t295 * t339 * t259 / 0.6e1 + 0.2e1 / 0.9e1 * t628 * t734 * t629 - 0.2e1 / 0.9e1 * t1605 + t628 * t295 * t339 * t250 / 0.9e1 - t1612 / 0.9e1 - 0.2e1 / 0.3e1 * t271 * t1474 * t742 - 0.5e1 / 0.27e2 * t1497 * t272 * t319 * t188 + 0.2e1 / 0.27e2 * t1622 + 0.7e1 / 0.9e1 * t1625 + t613 * t272 * t325 * t188 / 0.9e1;
1369 t1633 = t271 * t618 * t722;
1370 t1640 = t628 * t633 * t726;
1371 t1645 = t7 * t264 * t369;
1372 t1651 = t7 * t1232 * t320;
1373 t1654 = t7 * t496 * t326;
1374 t1657 = t7 * t1257 * t331;
1375 t1663 = t7 * t516 * t336;
1376 t1666 = t7 * t516 * t340;
1377 t1673 = t7 * t255 * t790;
1378 t1679 = 0.1e1 / t80 / t1201;
1379 t1681 = t8 * t8;
1380 t1687 = t7 * t488 * t768;
1381 t1690 = t7 * t193 * t772;
1382 t1692 = -t1633 / 0.9e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t330 * t250 + 0.2e1 / 0.27e2 * t1640 + t7 * t303 * t147 + t1645 / 0.12e2 - t1347 / 0.3e1 + t1353 / 0.3e1 - 0.5e1 / 0.18e2 * t1356 + t1651 / 0.27e2 + t1654 / 0.9e1 + t1657 / 0.27e2 + t7 * t1358 * t127 / 0.3e1 + 0.2e1 / 0.9e1 * t1663 + t1666 / 0.9e1 - t7 * t249 * t369 * t250 / 0.6e1 - t1673 / 0.3e1 - t7 * t1402 * t122 / 0.3e1 - t7 * t186 * t1679 * t38 * t1681 / 0.3e1 - t1687 / 0.9e1 - t1690 / 0.3e1;
1383 t1695 = t7 * t193 * t778;
1384 t1703 = t7 * t691 * t782;
1385 t1706 = t7 * t255 * t786;
1386 t1735 = -t1695 / 0.3e1 - t7 * t248 * t1679 * t64 * t1681 / 0.3e1 - 0.7e1 / 0.3e1 * t1703 - t1706 / 0.3e1 - t7 * t187 * t243 * t319 / 0.6e1 - t7 * t70 * t243 * t325 / 0.2e1 - t7 * t70 * t104 * t777 - t7 * t108 * t369 * t259 / 0.2e1 - 0.10e2 / 0.27e2 * t1363 - t1370 / 0.9e1 + 0.4e1 / 0.9e1 * t1373 + 0.4e1 / 0.9e1 * t1376 - 0.2e1 / 0.3e1 * t1379 + 0.41e2 / 0.9e1 * t1389 - t1392 / 0.3e1 - t1395 / 0.3e1 - 0.10e2 / 0.27e2 * t1411 - 0.41e2 / 0.9e1 * t1415 - 0.10e2 / 0.27e2 * t1159 + 0.2e1 / 0.9e1 * t1169;
1387 t1760 = t1172 / 0.9e1 + t1175 / 0.9e1 + t1234 / 0.27e2 - 0.10e2 / 0.27e2 * t1242 + t1259 / 0.27e2 + 0.4e1 / 0.27e2 * t1449 - 0.2e1 / 0.9e1 * t1453 + 0.4e1 / 0.27e2 * t1465 - 0.2e1 / 0.9e1 * t1472 - t1478 / 0.9e1 - t1486 / 0.9e1 + 0.2e1 / 0.9e1 * t1494 + 0.2e1 / 0.27e2 * t1503 + 0.2e1 / 0.27e2 * t1516 - 0.4e1 / 0.9e1 * t628 * t649 * t296 * t8 + 0.4e1 / 0.9e1 * t613 * t622 * t273 * t8 - 0.5e1 / 0.18e2 * t1287 + t1290 / 0.12e2 + t1337 + t1341;
1388 t1776 = t9 * (t1631 + t1692 + t1735 + t1760) + 0.2e1 / 0.9e1 * t720 - t724 / 0.3e1 + 0.2e1 / 0.9e1 * t728 + 0.2e1 / 0.3e1 * t732 - 0.2e1 / 0.3e1 * t736 - t740 / 0.3e1 - 0.2e1 / 0.3e1 * t744 + 0.2e1 / 0.9e1 * t616 - 0.4e1 / 0.9e1 * t620 - 0.2e1 / 0.3e1 * t625 + 0.2e1 / 0.9e1 * t631 - 0.4e1 / 0.9e1 * t635;
1389 t1790 = -0.2e1 / 0.3e1 * t639 - t643 / 0.3e1 - t647 / 0.3e1 + 0.2e1 / 0.3e1 * t652 - t490 / 0.9e1 + 0.2e1 / 0.3e1 * t498 - t510 / 0.9e1 + 0.2e1 / 0.3e1 * t518 - 0.2e1 / 0.3e1 * t525 - t528 / 0.3e1 - t552 / 0.3e1 + 0.2e1 / 0.3e1 * t663 + 0.5e1 / 0.3e1 * t669;
1390 t1803 = -0.2e1 * t673 + 0.2e1 / 0.3e1 * t675 - t679 / 0.3e1 - 0.4e1 / 0.3e1 * t682 - 0.4e1 / 0.3e1 * t685 - t689 - 0.5e1 / 0.3e1 * t693 - t697 - t754 / 0.9e1 - t757 / 0.3e1 - t760 / 0.9e1 - 0.2e1 / 0.3e1 * t763 - t766 / 0.3e1;
1391 t1810 = -t770 / 0.3e1 - t774 - t780 + 0.2e1 * t784 - t788 - t792 - t1147 - t1150 + t564 / 0.2e1 - t567 / 0.4e1 - t702 + t705 + t708 / 0.2e1 - t797 / 0.4e1;
1392 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1393 v4rho4[2] = t1776 + t1790 + t1803 + t1810;
1394
1395 t1864 = t271 * t618 * t821;
1396 t1870 = -0.5e1 / 0.27e2 * t1497 * t272 * t816 * t75 - t271 * t637 * t821 / 0.2e1 - t271 * t272 * t121 * t777 / 0.2e1 - 0.5e1 / 0.27e2 * t1510 * t295 * t825 * t109 - t628 * t649 * t330 * t8 / 0.3e1 - t294 * t295 * t126 * t785 / 0.2e1 + t294 * t649 * t8 * t339 / 0.2e1 - t271 * t272 * t838 * t75 / 0.6e1 + t628 * t734 * t726 / 0.3e1 + t294 * t1588 * t730 - t294 * t16 * t369 * t296 / 0.2e1 - t294 * t734 * t738 / 0.2e1 - t294 * t295 * t855 * t109 / 0.6e1 + t613 * t622 * t319 * t8 / 0.3e1 - t1864 / 0.6e1 - t271 * t622 * t8 * t325 / 0.2e1;
1397 t1872 = t294 * t633 * t831;
1398 t1882 = t7 * t114 * t888;
1399 t1887 = t7 * t1227 * t817;
1400 t1890 = t7 * t1253 * t826;
1401 t1893 = t7 * t193 * t839;
1402 t1896 = t7 * t508 * t843;
1403 t1899 = t7 * t255 * t847;
1404 t1902 = t7 * t255 * t851;
1405 t1904 = -t1872 / 0.6e1 + t1583 / 0.3e1 - t1605 / 0.3e1 - t1612 / 0.6e1 + t1622 / 0.9e1 - t1625 / 0.3e1 - t1633 / 0.6e1 + t1640 / 0.9e1 - t1882 / 0.8e1 + t1645 / 0.4e1 - 0.5e1 / 0.12e2 * t1356 + t1887 / 0.27e2 + t1890 / 0.27e2 - t1893 / 0.6e1 - t1896 / 0.6e1 - t1899 / 0.2e1 - t1902 / 0.2e1;
1406 t1907 = t7 * t255 * t856;
1407 t1919 = 0.6e1 * t539 + 0.12e2 * t1203;
1408 t1950 = -t1907 / 0.6e1 + t7 * t482 * t104 * t816 / 0.9e1 - t7 * t70 * t104 * t838 / 0.2e1 - t7 * t70 * t38 * t1919 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t303 * t369 * t8 - 0.3e1 / 0.2e1 * t7 * t108 * t147 * t785 + t7 * t108 * t64 * t1919 / 0.2e1 - t7 * t108 * t888 * t109 / 0.2e1 + t1651 / 0.9e1 + t1654 / 0.3e1 + t1657 / 0.9e1 + 0.2e1 / 0.3e1 * t1663 + t1666 / 0.3e1 - t1673 / 0.2e1 - t1687 / 0.6e1 - t1690 / 0.2e1 - t1695 / 0.2e1;
1409 t1967 = t1703 - t1706 / 0.2e1 - 0.5e1 / 0.9e1 * t1363 + t1374 + t1377 - t1389 / 0.3e1 - 0.5e1 / 0.9e1 * t1411 + t1415 / 0.3e1 - 0.5e1 / 0.27e2 * t1159 - 0.5e1 / 0.27e2 * t1242 + t1450 + t1466 + t613 * t272 * t821 * t75 / 0.3e1 + t628 * t295 * t831 * t109 / 0.3e1 - 0.5e1 / 0.36e2 * t1287 + t1337 + t1341;
1410 t1975 = t9 * (t1870 + t1904 + t1950 + t1967) - t824 - t834 + t720 / 0.3e1 - t724 / 0.2e1 + t728 / 0.3e1 + t732 - t736 - t740 / 0.2e1 - t744 - t1536;
1411 t1977 = -t1538 + t499 + t519 + t663 - t669 + t675 - t682 - t685 + t693 - t754 / 0.3e1 - t757;
1412 t1987 = -t760 / 0.3e1 - 0.2e1 * t763 - t766 - t770 / 0.2e1 - 0.3e1 / 0.2e1 * t774 - 0.3e1 / 0.2e1 * t780 + 0.3e1 * t784 - 0.3e1 / 0.2e1 * t788 - 0.3e1 / 0.2e1 * t792 + t820 + t829;
1413 t1990 = -t842 - t846 - t850 - t854 - t859 - t1147 - t1150 + t565 + 0.3e1 / 0.4e1 * t708 - 0.3e1 / 0.4e1 * t797 - t891;
1414 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1415 v4rho4[3] = t1975 + t1977 + t1987 + t1990;
1416
1417 t2003 = 0.4e1 / 0.3e1 * t663 + 0.4e1 / 0.3e1 * t675 - 0.2e1 / 0.3e1 * t754 - 0.2e1 * t757 - 0.2e1 / 0.3e1 * t760 - 0.4e1 * t763 - 0.2e1 * t766 + 0.4e1 / 0.9e1 * t819 - 0.2e1 * t823 + 0.4e1 / 0.9e1 * t828 - 0.2e1 * t833;
1418 t2049 = t879 * rho[1];
1419 t2061 = 0.1e1 / t364 / t60;
1420 t2072 = t339 * t339;
1421 t2082 = 0.12e2 * t539 + 0.12e2 * t1203;
1422 t2107 = 0.2e1 / 0.3e1 * t613 * t272 * t319 * t325 - 0.2e1 / 0.3e1 * t271 * t272 * t121 * t838 + 0.2e1 / 0.3e1 * t628 * t295 * t330 * t339 - 0.2e1 * t294 * t734 * t831 - 0.2e1 / 0.3e1 * t294 * t295 * t126 * t855 - 0.2e1 / 0.3e1 * t1864 - 0.2e1 / 0.3e1 * t1872 - 0.3e1 / 0.8e1 * t7 * t45 * (0.167200e6 / 0.81e2 * t19 * t46 / t50 / t139 * t61 - 0.948640e6 / 0.81e2 * t93 * t137 / t360 * t143 + 0.629120e6 / 0.27e2 * t232 * t359 / t55 / t1087 * t365 - 0.532480e6 / 0.27e2 * t588 * t878 / t50 / t2049 * t884 + 0.163840e6 / 0.27e2 * t1320 * t359 * t46 / t879 / t48 * t2061) - t1882 / 0.2e1 + t1645 / 0.2e1 - 0.5e1 / 0.9e1 * t1356 - t7 * t249 * t64 * t2072 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t502 * t147 * t825 - t7 * t70 * t38 * t2082 / 0.2e1 - t7 * t249 * t369 * t330 - 0.2e1 * t7 * t108 * t888 * t126 - 0.3e1 * t7 * t108 * t369 * t339 - 0.2e1 * t7 * t108 * t147 * t855 + t7 * t108 * t64 * t2082 / 0.2e1;
1423 t2108 = t319 * t319;
1424 t2113 = t325 * t325;
1425 t2118 = t330 * t330;
1426 t2137 = -0.5e1 / 0.27e2 * t7 * t1221 * t38 * t2108 - t7 * t187 * t38 * t2113 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1247 * t64 * t2118 + 0.4e1 / 0.27e2 * t1887 + 0.4e1 / 0.27e2 * t1890 - 0.2e1 / 0.3e1 * t1893 - 0.2e1 / 0.3e1 * t1896 - 0.2e1 * t1899 - 0.2e1 * t1902 - 0.2e1 / 0.3e1 * t1907 + 0.2e1 / 0.9e1 * t1651 + 0.2e1 / 0.3e1 * t1654 + 0.2e1 / 0.9e1 * t1657 + 0.4e1 / 0.3e1 * t1663 + 0.2e1 / 0.3e1 * t1666 - 0.20e2 / 0.27e2 * t1363 - 0.20e2 / 0.27e2 * t1411 + t1337 + t1341;
1427 t2140 = -t1147 - t1150 + t708 - 0.3e1 / 0.2e1 * t797 - 0.2e1 * t841 - 0.2e1 * t845 - 0.6e1 * t849 - 0.6e1 * t853 - 0.2e1 * t858 - 0.3e1 / 0.2e1 * t890 + t9 * (t2107 + t2137);
1428 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1429 v4rho4[4] = t2003 + t2140;
1430
1431 t2142 = t7 * t82 * t398;
1432 t2145 = t7 * t17 * t933;
1433 t2174 = t7 * t70 * t907;
1434 t2177 = t7 * t70 * t910;
1435 t2181 = t7 * t210 * t165 / 0.9e1;
1436 t2199 = t7 * t187 * t903;
1437 t2202 = t7 * t193 * t379;
1438 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1439 v4rho3sigma[0] = -t2142 / 0.2e1 - 0.3e1 / 0.2e1 * t2145 - 0.3e1 / 0.8e1 * t155 * t156 * (-0.13376e5 / 0.27e2 * t19 * sigma[0] * t570 * t35 + 0.89888e5 / 0.27e2 * t93 * t20 * t575 * t100 - 0.66976e5 / 0.9e1 * t232 * t94 * t582 * t239 + 0.61952e5 / 0.9e1 * t588 * t233 * t592 * t595 - 0.20480e5 / 0.9e1 * t1320 * t589 / t590 / t217 * t1326) - 0.4e1 * t2174 - 0.2e1 * t2177 + t2181 - t155 * t270 * t398 * t188 / 0.2e1 - 0.3e1 / 0.2e1 * t155 * t378 * t933 * t75 - 0.3e1 / 0.2e1 * t155 * t378 * t398 * t203 - t155 * t378 * t165 * t542 / 0.2e1 - 0.2e1 / 0.3e1 * t2199 - 0.2e1 / 0.3e1 * t2202 + t155 * t612 * t165 * t483 / 0.9e1 - t958 * t959 * t492 / 0.2e1;
1440
1441 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1442 v4rho3sigma[1] = 0.0e0;
1443
1444 t2212 = t7 * t249 * t943;
1445 t2215 = t7 * t255 * t405;
1446 t2218 = t7 * t108 * t947;
1447 t2222 = t7 * t264 * t178 / 0.9e1;
1448 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1449 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2212 - 0.2e1 / 0.3e1 * t2215 - 0.2e1 * t2218 + t2222 + t155 * t627 * t178 * t503 / 0.9e1 - t958 * t976 * t512 / 0.2e1 - t155 * t404 * t178 * t554 / 0.2e1;
1450
1451 t2239 = 0.4e1 / 0.3e1 * t954 * t964;
1452 t2250 = t81 * t1 * t3;
1453 t2251 = t2250 * t410;
1454 t2253 = t953 * t957;
1455 t2255 = 0.4e1 / 0.9e1 * t2253 * t960;
1456 t2257 = t285 * t1 * t3;
1457 t2258 = t2257 * t970;
1458 t2270 = t481 * t165;
1459 t2274 = t967 * t957;
1460 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1461 v4rho3sigma[3] = -t2142 / 0.3e1 - t2145 / 0.2e1 - 0.4e1 / 0.3e1 * t2174 - 0.2e1 / 0.3e1 * t2177 + t2181 - t2239 - t155 * t378 * t933 * t121 / 0.2e1 - t968 * t378 * t165 / 0.2e1 - 0.2e1 / 0.9e1 * t2199 - 0.4e1 / 0.9e1 * t2202 - 0.2e1 / 0.9e1 * t2251 - t2255 + t2258 / 0.6e1 - t958 * t186 * t398 * t273 / 0.3e1 - t958 * t959 * t641 / 0.6e1 - t968 * t378 * t398 * t8 + t958 * t2270 * t614 / 0.9e1 - t2274 * t959 * t623 / 0.3e1;
1462
1463 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1464 v4rho3sigma[4] = 0.0e0;
1465
1466 t2281 = t2250 * t415;
1467 t2284 = 0.4e1 / 0.9e1 * t2253 * t977;
1468 t2285 = t2257 * t981;
1469 t2287 = t501 * t178;
1470 t2300 = t2250 * t434;
1471 t2303 = 0.4e1 / 0.3e1 * t954 * t987;
1472 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1473 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2212 - 0.4e1 / 0.9e1 * t2215 - 0.2e1 / 0.3e1 * t2218 + t2222 - 0.2e1 / 0.9e1 * t2281 - t2284 - t2285 / 0.6e1 + t958 * t2287 * t629 / 0.9e1 + t2274 * t976 * t650 / 0.3e1 - t958 * t976 * t645 / 0.6e1 + t968 * t404 * t178 / 0.2e1 - t2300 / 0.6e1 - t2303 - t155 * t293 * t433 * t250 / 0.6e1 - t155 * t404 * t433 * t259 / 0.2e1;
1474
1475 t2316 = t954 * t992;
1476 t2328 = t954 * t996;
1477 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1478 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2251 - t2255 - t2239 - 0.4e1 / 0.3e1 * t2258 - 0.2e1 / 0.9e1 * t2202 + t2181 - t2142 / 0.6e1 - 0.2e1 / 0.9e1 * t2316 + t958 * t2270 * t718 / 0.9e1 - t155 * t270 * t398 * t319 / 0.6e1 - t2274 * t959 * t742 / 0.3e1 - 0.2e1 / 0.3e1 * t2328 - t958 * t959 * t722 / 0.6e1 - t155 * t378 * t398 * t325 / 0.2e1 - t155 * t378 * t165 * t777 / 0.2e1;
1479
1480 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1481 v4rho3sigma[7] = 0.0e0;
1482
1483 t2345 = t954 * t1001;
1484 t2353 = t954 * t1005;
1485 t2362 = t954 * t1008;
1486 t2371 = t954 * t1031;
1487 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1488 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2281 - t2284 + 0.4e1 / 0.3e1 * t2285 - 0.2e1 / 0.9e1 * t2215 + t2222 - t2300 / 0.3e1 - t2303 - 0.2e1 / 0.9e1 * t2345 + t958 * t2287 * t726 / 0.9e1 + t2274 * t976 * t730 / 0.3e1 - 0.4e1 / 0.3e1 * t2353 - t958 * t248 * t433 * t296 / 0.3e1 + t968 * t404 * t433 * t8 - 0.2e1 / 0.3e1 * t2362 - t958 * t976 * t738 / 0.6e1 - t155 * t404 * t178 * t785 / 0.2e1 - t2371 / 0.2e1 - t155 * t404 * t1030 * t109 / 0.2e1;
1489
1490 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1491 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2251 - 0.2e1 / 0.3e1 * t2316 - 0.2e1 * t2328 + t2181 + t155 * t612 * t165 * t816 / 0.9e1 - t958 * t959 * t821 / 0.2e1 - t155 * t378 * t165 * t838 / 0.2e1;
1492
1493 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1494 v4rho3sigma[10] = 0.0e0;
1495
1496 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1497 v4rho3sigma[11] = -0.3e1 / 0.8e1 * t155 * t169 * (-0.13376e5 / 0.27e2 * t19 * sigma[2] * t861 * t61 + 0.89888e5 / 0.27e2 * t93 * t46 * t866 * t143 - 0.66976e5 / 0.9e1 * t232 * t137 * t873 * t365 + 0.61952e5 / 0.9e1 * t588 * t359 * t881 * t884 - 0.20480e5 / 0.9e1 * t1320 * t878 / t879 / t345 * t2061) + t2222 - t2300 / 0.2e1 - 0.3e1 / 0.2e1 * t2371 - t155 * t293 * t433 * t330 / 0.2e1 - 0.3e1 / 0.2e1 * t155 * t404 * t1030 * t126 - 0.3e1 / 0.2e1 * t155 * t404 * t433 * t339 - t155 * t404 * t178 * t855 / 0.2e1 - 0.2e1 / 0.3e1 * t2281 - 0.2e1 / 0.3e1 * t2345 - 0.4e1 * t2353 - 0.2e1 * t2362 + t155 * t627 * t178 * t825 / 0.9e1 - t958 * t976 * t831 / 0.2e1;
1498
1499 t2447 = t2250 * t451 / 0.6e1;
1500 t2448 = t954 * t1037;
1501 t2450 = t954 * t1059;
1502 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1503 v4rho2sigma2[0] = -t2447 - 0.4e1 / 0.3e1 * t2448 - t2450 - t155 * t270 * t450 * t188 / 0.6e1 - t155 * t378 * t1058 * t75 - t155 * t378 * t450 * t203 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * (0.608e3 / 0.9e1 * t19 * t220 * t35 - 0.6560e4 / 0.9e1 * t93 * t226 * t100 * sigma[0] + 0.6272e4 / 0.3e1 * t232 * t236 * t239 * t20 - 0.6848e4 / 0.3e1 * t588 * t94 * t928 * t595 + 0.2560e4 / 0.3e1 * t1320 * t233 / t590 / t21 * t1326);
1504
1505 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1506 v4rho2sigma2[1] = 0.0e0;
1507
1508 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1509 v4rho2sigma2[2] = 0.0e0;
1510
1511 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1512 v4rho2sigma2[3] = 0.0e0;
1513
1514 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1515 v4rho2sigma2[4] = 0.0e0;
1516
1517 t2488 = t2250 * t468 / 0.6e1;
1518 t2489 = t954 * t1064;
1519 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1520 v4rho2sigma2[5] = -t2488 - 0.4e1 / 0.3e1 * t2489 - t155 * t293 * t467 * t250 / 0.6e1 - t155 * t404 * t467 * t259 / 0.2e1;
1521
1522 t2501 = t954 * t1068;
1523 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1524 v4rho2sigma2[6] = -t2447 - 0.2e1 / 0.3e1 * t2448 - t2450 / 0.2e1 - 0.2e1 / 0.3e1 * t2501 - t958 * t186 * t450 * t273 / 0.6e1 - t155 * t378 * t1058 * t121 / 0.2e1 - t968 * t378 * t450 * t8 / 0.2e1;
1525
1526 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1527 v4rho2sigma2[7] = 0.0e0;
1528
1529 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1530 v4rho2sigma2[8] = 0.0e0;
1531
1532 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1533 v4rho2sigma2[9] = 0.0e0;
1534
1535 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1536 v4rho2sigma2[10] = 0.0e0;
1537
1538 t2516 = t954 * t1073;
1539 t2526 = t954 * t1095;
1540 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1541 v4rho2sigma2[11] = -t2488 - 0.2e1 / 0.3e1 * t2489 - 0.2e1 / 0.3e1 * t2516 - t958 * t248 * t467 * t296 / 0.6e1 + t968 * t404 * t467 * t8 / 0.2e1 - t2526 / 0.2e1 - t155 * t404 * t1094 * t109 / 0.2e1;
1542
1543 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1544 v4rho2sigma2[12] = -t2447 - 0.4e1 / 0.3e1 * t2501 - t155 * t270 * t450 * t319 / 0.6e1 - t155 * t378 * t450 * t325 / 0.2e1;
1545
1546 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1547 v4rho2sigma2[13] = 0.0e0;
1548
1549 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1550 v4rho2sigma2[14] = 0.0e0;
1551
1552 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1553 v4rho2sigma2[15] = 0.0e0;
1554
1555 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556 v4rho2sigma2[16] = 0.0e0;
1557
1558 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559 v4rho2sigma2[17] = -t2488 - 0.4e1 / 0.3e1 * t2516 - t2526 - t155 * t293 * t467 * t330 / 0.6e1 - t155 * t404 * t1094 * t126 - t155 * t404 * t467 * t339 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t169 * (0.608e3 / 0.9e1 * t19 * t348 * t61 - 0.6560e4 / 0.9e1 * t93 * t354 * t143 * sigma[2] + 0.6272e4 / 0.3e1 * t232 * t362 * t365 * t46 - 0.6848e4 / 0.3e1 * t588 * t137 * t1025 * t884 + 0.2560e4 / 0.3e1 * t1320 * t359 / t879 / t47 * t2061);
1560
1561 t2578 = t954 * t1112;
1562 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1563 v4rhosigma3[0] = -t2578 / 0.2e1 - t155 * t378 * t1111 * t75 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t156 * (0.96e2 * t93 * t1043 - 0.480e3 * t232 * t393 * t239 * sigma[0] + 0.704e3 * t588 * t20 * t1053 * t595 - 0.320e3 * t1320 * t94 / t1312 * t1326);
1564
1565 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1566 v4rhosigma3[1] = 0.0e0;
1567
1568 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1569 v4rhosigma3[2] = 0.0e0;
1570
1571 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1572 v4rhosigma3[3] = 0.0e0;
1573
1574 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1575 v4rhosigma3[4] = 0.0e0;
1576
1577 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1578 v4rhosigma3[5] = 0.0e0;
1579
1580 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1581 v4rhosigma3[6] = 0.0e0;
1582
1583 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1584 v4rhosigma3[7] = 0.0e0;
1585
1586 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1587 v4rhosigma3[8] = 0.0e0;
1588
1589 t2603 = t954 * t1129;
1590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1591 v4rhosigma3[9] = -t155 * t404 * t1128 * t109 / 0.2e1 - t2603 / 0.2e1;
1592
1593 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594 v4rhosigma3[10] = -t155 * t378 * t1111 * t121 / 0.2e1 - t2578 / 0.2e1;
1595
1596 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597 v4rhosigma3[11] = 0.0e0;
1598
1599 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600 v4rhosigma3[12] = 0.0e0;
1601
1602 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603 v4rhosigma3[13] = 0.0e0;
1604
1605 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606 v4rhosigma3[14] = 0.0e0;
1607
1608 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1609 v4rhosigma3[15] = 0.0e0;
1610
1611 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1612 v4rhosigma3[16] = 0.0e0;
1613
1614 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1615 v4rhosigma3[17] = 0.0e0;
1616
1617 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1618 v4rhosigma3[18] = 0.0e0;
1619
1620 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1621 v4rhosigma3[19] = -t2603 / 0.2e1 - t155 * t404 * t1128 * t126 / 0.2e1 - 0.3e1 / 0.8e1 * t155 * t169 * (0.96e2 * t93 * t1079 - 0.480e3 * t232 * t428 * t365 * sigma[2] + 0.704e3 * t588 * t46 * t1089 * t884 - 0.320e3 * t1320 * t137 / t2049 * t2061);
1622
1623 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1624 v4sigma4[0] = -0.3e1 / 0.8e1 * t155 * t156 * (0.72e2 * t232 * t445 * t239 - 0.192e3 * t588 * sigma[0] * t1106 * t595 + 0.120e3 * t1320 * t20 / t590 * t1326);
1625
1626 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1627 v4sigma4[1] = 0.0e0;
1628
1629 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1630 v4sigma4[2] = 0.0e0;
1631
1632 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1633 v4sigma4[3] = 0.0e0;
1634
1635 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1636 v4sigma4[4] = 0.0e0;
1637
1638 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1639 v4sigma4[5] = 0.0e0;
1640
1641 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1642 v4sigma4[6] = 0.0e0;
1643
1644 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645 v4sigma4[7] = 0.0e0;
1646
1647 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1648 v4sigma4[8] = 0.0e0;
1649
1650 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1651 v4sigma4[9] = 0.0e0;
1652
1653 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1654 v4sigma4[10] = 0.0e0;
1655
1656 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1657 v4sigma4[11] = 0.0e0;
1658
1659 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1660 v4sigma4[12] = 0.0e0;
1661
1662 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1663 v4sigma4[13] = 0.0e0;
1664
1665 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1666 v4sigma4[14] = -0.3e1 / 0.8e1 * t155 * t169 * (0.72e2 * t232 * t462 * t365 - 0.192e3 * t588 * sigma[2] * t1123 * t884 + 0.120e3 * t1320 * t46 / t879 * t2061);
1667
1668 #ifndef XC_DONT_COMPILE_MXC
1669
1670 if(order < 5) return;
1671
1672
1673 #endif
1674
1675 #endif
1676
1677 #endif
1678
1679 #endif
1680
1681 #endif
1682
1683
1684 }
1685
1686