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_ol2.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, t16, t17, t20, t21;
25 double t22, t23, t28, t29, t30, t32, t34;
26
27 #ifndef XC_DONT_COMPILE_VXC
28 double t38, t39, t40, t42, t43, t49, t52, t53;
29 double t54, t55, t58, t62, t65, t66, t69, t74;
30
31 #ifndef XC_DONT_COMPILE_FXC
32 double t77, t78, t85, t87, t88, t94, t97, t100;
33 double t101, t102, t103, t105, t106, t109, t114, t124;
34 double t125, t127, t131, t134, t135, t138, t139, t140;
35 double t142, t148;
36
37 #ifndef XC_DONT_COMPILE_KXC
38 double t152, t153, t164, t165, t171, t177, t178, t182;
39 double t183, t184, t186, t187, t188, t193, t198, t215;
40 double t217, t221, t230, t235, t236, t240, t243, t245;
41 double t248, t249, t253, t254, t256, t260;
42
43 #ifndef XC_DONT_COMPILE_LXC
44 double t275, t276, t290, t302, t387;
45 #endif
46
47 #endif
48
49 #endif
50
51 #endif
52
53 #endif
54
55
56 gga_x_ol2_params *params;
57
58 assert(p->params != NULL);
59 params = (gga_x_ol2_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->bb * sigma[0];
71 t13 = rho[0] * rho[0];
72 t14 = t10 * t10;
73 t16 = 0.1e1 / t14 / t13;
74 t17 = t9 * t16;
75 t20 = sqrt(sigma[0]);
76 t21 = params->cc * t20;
77 t22 = t10 * rho[0];
78 t23 = 0.1e1 / t22;
79 t28 = 0.4e1 * t20 * t8 * t23 + t8;
80 t29 = 0.1e1 / t28;
81 t30 = t8 * t23 * t29;
82 t32 = params->aa + 0.13888888888888888889e-1 * t12 * t17 + t21 * t30;
83 t34 = t7 * t11 * t32;
84 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
85 zk[0] = -0.3e1 / 0.16e2 * t34;
86
87 #ifndef XC_DONT_COMPILE_VXC
88
89 if(order < 1) return;
90
91
92 t38 = t22 * t1 * t3;
93 t39 = t6 * t9;
94 t40 = t13 * rho[0];
95 t42 = 0.1e1 / t14 / t40;
96 t43 = t9 * t42;
97 t49 = t8 / t10 / t13 * t29;
98 t52 = params->cc * sigma[0];
99 t53 = t28 * t28;
100 t54 = 0.1e1 / t53;
101 t55 = t43 * t54;
102 t58 = -0.37037037037037037037e-1 * t12 * t43 - 0.4e1 / 0.3e1 * t21 * t49 + 0.16e2 / 0.3e1 * t52 * t55;
103 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
104 vrho[0] = -t34 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t39 * t58;
105
106 t62 = params->bb * t9;
107 t65 = 0.1e1 / t20;
108 t66 = params->cc * t65;
109 t69 = params->cc * t9;
110 t74 = t39 * (0.13888888888888888889e-1 * t62 * t16 + t66 * t30 / 0.2e1 - 0.2e1 * t69 * t16 * t54);
111 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
112 vsigma[0] = -0.3e1 / 0.16e2 * t38 * t74;
113
114 #ifndef XC_DONT_COMPILE_FXC
115
116 if(order < 2) return;
117
118
119 t77 = 0.1e1 / t14;
120 t78 = t9 * t77;
121 t85 = t13 * t13;
122 t87 = 0.1e1 / t14 / t85;
123 t88 = t9 * t87;
124 t94 = t8 / t10 / t40 * t29;
125 t97 = t88 * t54;
126 t100 = t20 * sigma[0];
127 t101 = params->cc * t100;
128 t102 = t85 * t13;
129 t103 = 0.1e1 / t102;
130 t105 = 0.1e1 / t53 / t28;
131 t106 = t103 * t105;
132 t109 = 0.13580246913580246914e0 * t12 * t88 + 0.28e2 / 0.9e1 * t21 * t94 - 0.80e2 / 0.3e1 * t52 * t97 + 0.1024e4 / 0.9e1 * t101 * t106;
133 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
134 v2rho2[0] = -t7 * t78 * t32 / 0.12e2 - t7 * t11 * t58 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t39 * t109;
135
136 t114 = t10 * t1 * t3;
137 t124 = t85 * rho[0];
138 t125 = 0.1e1 / t124;
139 t127 = t105 * t20;
140 t131 = t39 * (-0.37037037037037037037e-1 * t62 * t42 - 0.2e1 / 0.3e1 * t66 * t49 + 0.8e1 * t69 * t42 * t54 - 0.128e3 / 0.3e1 * params->cc * t125 * t127);
141 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
142 v2rhosigma[0] = -t114 * t74 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t131;
143
144 t134 = 0.1e1 / t100;
145 t135 = params->cc * t134;
146 t138 = 0.1e1 / sigma[0];
147 t139 = params->cc * t138;
148 t140 = t17 * t54;
149 t142 = 0.1e1 / t85;
150 t148 = t39 * (-t135 * t30 / 0.4e1 - t139 * t140 + 0.16e2 * params->cc * t142 * t105 * t65);
151 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
152 v2sigma2[0] = -0.3e1 / 0.16e2 * t38 * t148;
153
154 #ifndef XC_DONT_COMPILE_KXC
155
156 if(order < 3) return;
157
158
159 t152 = 0.1e1 / t14 / rho[0];
160 t153 = t9 * t152;
161 t164 = 0.1e1 / t14 / t124;
162 t165 = t9 * t164;
163 t171 = t8 / t10 / t85 * t29;
164 t177 = t85 * t40;
165 t178 = 0.1e1 / t177;
166 t182 = sigma[0] * sigma[0];
167 t183 = params->cc * t182;
168 t184 = t85 * t85;
169 t186 = 0.1e1 / t10 / t184;
170 t187 = t53 * t53;
171 t188 = 0.1e1 / t187;
172 t193 = -0.63374485596707818932e0 * t12 * t165 - 0.280e3 / 0.27e2 * t21 * t171 + 0.3808e4 / 0.27e2 * t52 * t165 * t54 - 0.11264e5 / 0.9e1 * t101 * t178 * t105 + 0.16384e5 / 0.9e1 * t183 * t186 * t188 * t8;
173 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
174 v3rho3[0] = t7 * t153 * t32 / 0.18e2 - t7 * t78 * t58 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t109 - 0.3e1 / 0.16e2 * t38 * t39 * t193;
175
176 t198 = t77 * t1 * t3;
177 t215 = params->cc / t10 / t177;
178 t217 = t188 * sigma[0] * t8;
179 t221 = t39 * (0.13580246913580246914e0 * t62 * t87 + 0.14e2 / 0.9e1 * t66 * t94 - 0.296e3 / 0.9e1 * t69 * t87 * t54 + 0.384e3 * params->cc * t103 * t127 - 0.2048e4 / 0.3e1 * t215 * t217);
180 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
181 v3rho2sigma[0] = -t198 * t74 / 0.12e2 - t114 * t131 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t221;
182
183 t230 = t125 * t105;
184 t235 = params->cc / t10 / t102;
185 t236 = t188 * t8;
186 t240 = t39 * (t135 * t49 / 0.3e1 + 0.4e1 / 0.3e1 * t139 * t55 - 0.256e3 / 0.3e1 * t66 * t230 + 0.256e3 * t235 * t236);
187 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
188 v3rhosigma2[0] = -t114 * t148 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t240;
189
190 t243 = t20 * t182;
191 t245 = params->cc / t243;
192 t248 = 0.1e1 / t182;
193 t249 = params->cc * t248;
194 t253 = 0.1e1 / t10 / t124;
195 t254 = params->cc * t253;
196 t256 = t188 * t138 * t8;
197 t260 = t39 * (0.3e1 / 0.8e1 * t245 * t30 + 0.3e1 / 0.2e1 * t249 * t140 - 0.96e2 * t254 * t256);
198 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
199 v3sigma3[0] = -0.3e1 / 0.16e2 * t38 * t260;
200
201 #ifndef XC_DONT_COMPILE_LXC
202
203 if(order < 4) return;
204
205
206 t275 = 0.1e1 / t14 / t102;
207 t276 = t9 * t275;
208 t290 = t184 * rho[0];
209 t302 = 0.1e1 / t187 / t28;
210 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
211 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t17 * t32 + 0.2e1 / 0.9e1 * t7 * t153 * t58 - t7 * t78 * t109 / 0.2e1 - t7 * t11 * t193 - 0.3e1 / 0.16e2 * t38 * t39 * (0.35912208504801097395e1 * t12 * t276 + 0.3640e4 / 0.81e2 * t21 * t8 * t253 * t29 - 0.23072e5 / 0.27e2 * t52 * t276 * t54 + 0.953344e6 / 0.81e2 * t101 / t184 * t105 - 0.950272e6 / 0.27e2 * t183 / t10 / t290 * t188 * t8 + 0.1048576e7 / 0.27e2 * params->cc * t243 / t14 / t184 / t13 * t302 * t9);
212
213 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
214 v4rho3sigma[0] = t152 * t1 * t3 * t74 / 0.18e2 - t198 * t131 / 0.4e1 - 0.3e1 / 0.4e1 * t114 * t221 - 0.3e1 / 0.16e2 * t38 * t39 * (-0.63374485596707818932e0 * t62 * t164 - 0.140e3 / 0.27e2 * t66 * t171 + 0.1456e4 / 0.9e1 * t69 * t164 * t54 - 0.81152e5 / 0.27e2 * params->cc * t178 * t127 + 0.100352e6 / 0.9e1 * params->cc * t186 * t217 - 0.131072e6 / 0.9e1 * params->cc / t14 / t290 * t302 * t100 * t9);
215
216 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
217 v4rho2sigma2[0] = -t198 * t148 / 0.12e2 - t114 * t240 / 0.2e1 - 0.3e1 / 0.16e2 * t38 * t39 * (-0.7e1 / 0.9e1 * t135 * t94 - 0.28e2 / 0.9e1 * t139 * t97 + 0.4096e4 / 0.9e1 * t66 * t106 - 0.8960e4 / 0.3e1 * t215 * t236 + 0.16384e5 / 0.3e1 * params->cc / t14 / t184 * t302 * t9 * t20);
218
219 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
220 v4rhosigma3[0] = -t114 * t260 / 0.4e1 - 0.3e1 / 0.16e2 * t38 * t39 * (-t245 * t49 / 0.2e1 - 0.2e1 * t249 * t55 + 0.32e2 * t135 * t230 + 0.512e3 * t235 * t256 - 0.2048e4 * params->cc / t14 / t177 * t302 * t65 * t9);
221
222 t387 = t182 * sigma[0];
223 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
224 v4sigma4[0] = -0.3e1 / 0.16e2 * t38 * t39 * (-0.15e2 / 0.16e2 * params->cc / t20 / t387 * t30 - 0.15e2 / 0.4e1 * params->cc / t387 * t140 - 0.12e2 * t245 * t142 * t105 + 0.768e3 * params->cc * t275 * t302 * t134 * t9 + 0.96e2 * t254 * t188 * t248 * t8);
225
226 #ifndef XC_DONT_COMPILE_MXC
227
228 if(order < 5) return;
229
230
231 #endif
232
233 #endif
234
235 #endif
236
237 #endif
238
239 #endif
240
241
242 }
243
244
245 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 *))246 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
247 {
248
249 #ifndef XC_DONT_COMPILE_EXC
250 double t1, t3, t4, t5, t6, t7, t8, t9;
251 double t10, t11, t13, t16, t17, t18, t19, t20;
252 double t23, t24, t25, t27, t29;
253
254 #ifndef XC_DONT_COMPILE_VXC
255 double t32, t33, t34, t36, t41, t44, t45, t46;
256 double t47, t50, t56, t57, t64;
257
258 #ifndef XC_DONT_COMPILE_FXC
259 double t67, t68, t74, t76, t81, t84, t87, t88;
260 double t89, t90, t92, t93, t96, t100, t110, t111;
261 double t113, t117, t120, t121, t124, t125, t126, t128;
262 double t134;
263
264 #ifndef XC_DONT_COMPILE_KXC
265 double t138, t139, t150, t155, t161, t162, t166, t167;
266 double t168, t170, t171, t172, t176, t180, t196, t197;
267 double t201, t210, t215, t219, t222, t224, t227, t228;
268 double t232, t233, t234, t238;
269
270 #ifndef XC_DONT_COMPILE_LXC
271 double t254, t267, t278, t357;
272 #endif
273
274 #endif
275
276 #endif
277
278 #endif
279
280 #endif
281
282
283 gga_x_ol2_params *params;
284
285 assert(p->params != NULL);
286 params = (gga_x_ol2_params * )(p->params);
287
288 t1 = M_CBRT3;
289 t3 = POW_1_3(0.1e1 / M_PI);
290 t4 = t1 * t3;
291 t5 = M_CBRT4;
292 t6 = t5 * t5;
293 t7 = POW_1_3(rho[0]);
294 t8 = t6 * t7;
295 t9 = params->bb * sigma[0];
296 t10 = rho[0] * rho[0];
297 t11 = t7 * t7;
298 t13 = 0.1e1 / t11 / t10;
299 t16 = sqrt(sigma[0]);
300 t17 = params->cc * t16;
301 t18 = t7 * rho[0];
302 t19 = 0.1e1 / t18;
303 t20 = M_CBRT2;
304 t23 = 0.4e1 * t16 * t19 + t20;
305 t24 = 0.1e1 / t23;
306 t25 = t19 * t24;
307 t27 = params->aa + 0.13888888888888888889e-1 * t9 * t13 + t17 * t25;
308 t29 = t4 * t8 * t27;
309 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
310 zk[0] = -0.3e1 / 0.8e1 * t29;
311
312 #ifndef XC_DONT_COMPILE_VXC
313
314 if(order < 1) return;
315
316
317 t32 = t18 * t1;
318 t33 = t3 * t6;
319 t34 = t10 * rho[0];
320 t36 = 0.1e1 / t11 / t34;
321 t41 = 0.1e1 / t7 / t10 * t24;
322 t44 = params->cc * sigma[0];
323 t45 = t23 * t23;
324 t46 = 0.1e1 / t45;
325 t47 = t36 * t46;
326 t50 = -0.37037037037037037037e-1 * t9 * t36 - 0.4e1 / 0.3e1 * t17 * t41 + 0.16e2 / 0.3e1 * t44 * t47;
327 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
328 vrho[0] = -t29 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t33 * t50;
329
330 t56 = 0.1e1 / t16;
331 t57 = params->cc * t56;
332 t64 = t33 * (0.13888888888888888889e-1 * params->bb * t13 + t57 * t25 / 0.2e1 - 0.2e1 * params->cc * t13 * t46);
333 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
334 vsigma[0] = -0.3e1 / 0.8e1 * t32 * t64;
335
336 #ifndef XC_DONT_COMPILE_FXC
337
338 if(order < 2) return;
339
340
341 t67 = 0.1e1 / t11;
342 t68 = t6 * t67;
343 t74 = t10 * t10;
344 t76 = 0.1e1 / t11 / t74;
345 t81 = 0.1e1 / t7 / t34 * t24;
346 t84 = t76 * t46;
347 t87 = t16 * sigma[0];
348 t88 = params->cc * t87;
349 t89 = t74 * t10;
350 t90 = 0.1e1 / t89;
351 t92 = 0.1e1 / t45 / t23;
352 t93 = t90 * t92;
353 t96 = 0.13580246913580246914e0 * t9 * t76 + 0.28e2 / 0.9e1 * t17 * t81 - 0.80e2 / 0.3e1 * t44 * t84 + 0.512e3 / 0.9e1 * t88 * t93;
354 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
355 v2rho2[0] = -t4 * t68 * t27 / 0.6e1 - t4 * t8 * t50 - 0.3e1 / 0.8e1 * t32 * t33 * t96;
356
357 t100 = t7 * t1;
358 t110 = t74 * rho[0];
359 t111 = 0.1e1 / t110;
360 t113 = t92 * t16;
361 t117 = t33 * (-0.37037037037037037037e-1 * params->bb * t36 - 0.2e1 / 0.3e1 * t57 * t41 + 0.8e1 * params->cc * t36 * t46 - 0.64e2 / 0.3e1 * params->cc * t111 * t113);
362 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
363 v2rhosigma[0] = -t100 * t64 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t117;
364
365 t120 = 0.1e1 / t87;
366 t121 = params->cc * t120;
367 t124 = 0.1e1 / sigma[0];
368 t125 = params->cc * t124;
369 t126 = t13 * t46;
370 t128 = 0.1e1 / t74;
371 t134 = t33 * (-t121 * t25 / 0.4e1 - t125 * t126 + 0.8e1 * params->cc * t128 * t92 * t56);
372 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
373 v2sigma2[0] = -0.3e1 / 0.8e1 * t32 * t134;
374
375 #ifndef XC_DONT_COMPILE_KXC
376
377 if(order < 3) return;
378
379
380 t138 = 0.1e1 / t11 / rho[0];
381 t139 = t6 * t138;
382 t150 = 0.1e1 / t11 / t110;
383 t155 = 0.1e1 / t7 / t74 * t24;
384 t161 = t74 * t34;
385 t162 = 0.1e1 / t161;
386 t166 = sigma[0] * sigma[0];
387 t167 = params->cc * t166;
388 t168 = t74 * t74;
389 t170 = 0.1e1 / t7 / t168;
390 t171 = t45 * t45;
391 t172 = 0.1e1 / t171;
392 t176 = -0.63374485596707818932e0 * t9 * t150 - 0.280e3 / 0.27e2 * t17 * t155 + 0.3808e4 / 0.27e2 * t44 * t150 * t46 - 0.5632e4 / 0.9e1 * t88 * t162 * t92 + 0.8192e4 / 0.9e1 * t167 * t170 * t172;
393 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
394 v3rho3[0] = t4 * t139 * t27 / 0.9e1 - t4 * t68 * t50 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t96 - 0.3e1 / 0.8e1 * t32 * t33 * t176;
395
396 t180 = t67 * t1;
397 t196 = params->cc / t7 / t161;
398 t197 = t172 * sigma[0];
399 t201 = t33 * (0.13580246913580246914e0 * params->bb * t76 + 0.14e2 / 0.9e1 * t57 * t81 - 0.296e3 / 0.9e1 * params->cc * t76 * t46 + 0.192e3 * params->cc * t90 * t113 - 0.1024e4 / 0.3e1 * t196 * t197);
400 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
401 v3rho2sigma[0] = -t180 * t64 / 0.6e1 - t100 * t117 - 0.3e1 / 0.8e1 * t32 * t201;
402
403 t210 = t111 * t92;
404 t215 = params->cc / t7 / t89;
405 t219 = t33 * (t121 * t41 / 0.3e1 + 0.4e1 / 0.3e1 * t125 * t47 - 0.128e3 / 0.3e1 * t57 * t210 + 0.128e3 * t215 * t172);
406 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
407 v3rhosigma2[0] = -t100 * t134 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t219;
408
409 t222 = t16 * t166;
410 t224 = params->cc / t222;
411 t227 = 0.1e1 / t166;
412 t228 = params->cc * t227;
413 t232 = 0.1e1 / t7 / t110;
414 t233 = params->cc * t232;
415 t234 = t172 * t124;
416 t238 = t33 * (0.3e1 / 0.8e1 * t224 * t25 + 0.3e1 / 0.2e1 * t228 * t126 - 0.48e2 * t233 * t234);
417 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
418 v3sigma3[0] = -0.3e1 / 0.8e1 * t32 * t238;
419
420 #ifndef XC_DONT_COMPILE_LXC
421
422 if(order < 4) return;
423
424
425 t254 = 0.1e1 / t11 / t89;
426 t267 = t168 * rho[0];
427 t278 = 0.1e1 / t171 / t23;
428 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
429 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t13 * t27 + 0.4e1 / 0.9e1 * t4 * t139 * t50 - t4 * t68 * t96 - 0.2e1 * t4 * t8 * t176 - 0.3e1 / 0.8e1 * t32 * t33 * (0.35912208504801097395e1 * t9 * t254 + 0.3640e4 / 0.81e2 * t17 * t232 * t24 - 0.23072e5 / 0.27e2 * t44 * t254 * t46 + 0.476672e6 / 0.81e2 * t88 / t168 * t92 - 0.475136e6 / 0.27e2 * t167 / t7 / t267 * t172 + 0.524288e6 / 0.27e2 * params->cc * t222 / t11 / t168 / t10 * t278);
430
431 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
432 v4rho3sigma[0] = t138 * t1 * t64 / 0.9e1 - t180 * t117 / 0.2e1 - 0.3e1 / 0.2e1 * t100 * t201 - 0.3e1 / 0.8e1 * t32 * t33 * (-0.63374485596707818932e0 * params->bb * t150 - 0.140e3 / 0.27e2 * t57 * t155 + 0.1456e4 / 0.9e1 * params->cc * t150 * t46 - 0.40576e5 / 0.27e2 * params->cc * t162 * t113 + 0.50176e5 / 0.9e1 * params->cc * t170 * t197 - 0.65536e5 / 0.9e1 * params->cc / t11 / t267 * t278 * t87);
433
434 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
435 v4rho2sigma2[0] = -t180 * t134 / 0.6e1 - t100 * t219 - 0.3e1 / 0.8e1 * t32 * t33 * (-0.7e1 / 0.9e1 * t121 * t81 - 0.28e2 / 0.9e1 * t125 * t84 + 0.2048e4 / 0.9e1 * t57 * t93 - 0.4480e4 / 0.3e1 * t196 * t172 + 0.8192e4 / 0.3e1 * params->cc / t11 / t168 * t278 * t16);
436
437 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
438 v4rhosigma3[0] = -t100 * t238 / 0.2e1 - 0.3e1 / 0.8e1 * t32 * t33 * (-t224 * t41 / 0.2e1 - 0.2e1 * t228 * t47 + 0.16e2 * t121 * t210 + 0.256e3 * t215 * t234 - 0.1024e4 * params->cc / t11 / t161 * t278 * t56);
439
440 t357 = t166 * sigma[0];
441 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
442 v4sigma4[0] = -0.3e1 / 0.8e1 * t32 * t33 * (-0.15e2 / 0.16e2 * params->cc / t16 / t357 * t25 - 0.15e2 / 0.4e1 * params->cc / t357 * t126 - 0.6e1 * t224 * t128 * t92 + 0.384e3 * params->cc * t254 * t278 * t120 + 0.48e2 * t233 * t172 * t227);
443
444 #ifndef XC_DONT_COMPILE_MXC
445
446 if(order < 5) return;
447
448
449 #endif
450
451 #endif
452
453 #endif
454
455 #endif
456
457 #endif
458
459
460 }
461
462
463 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 *))464 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
465 {
466
467 #ifndef XC_DONT_COMPILE_EXC
468 double t1, t3, t4, t5, t6, t7, t8, t9;
469 double t10, t11, t13, t14, t15, t16, t17, t18;
470 double t19, t20, t21, t23, t26, t27, t29, t30;
471 double t33, t34, t35, t37, t39, t41, t42, t43;
472 double t44, t45, t46, t47, t48, t50, t53, t54;
473 double t56, t59, t60, t61, t63, t65;
474
475 #ifndef XC_DONT_COMPILE_VXC
476 double t67, t68, t69, t70, t71, t72, t74, t75;
477 double t77, t78, t79, t80, t81, t83, t84, t85;
478 double t87, t92, t95, t96, t97, t98, t101, t103;
479 double t104, t105, t106, t107, t109, t110, t111, t113;
480 double t114, t118, t119, t121, t122, t123, t124, t126;
481 double t127, t128, t130, t135, t138, t139, t140, t141;
482 double t144, t146, t147, t151, t152, t153, t156, t157;
483 double t163, t167, t170, t171, t177;
484
485 #ifndef XC_DONT_COMPILE_FXC
486 double t181, t183, t184, t185, t186, t187, t188, t190;
487 double t191, t192, t194, t195, t196, t198, t199, t200;
488 double t201, t202, t203, t205, t206, t208, t209, t211;
489 double t212, t214, t215, t216, t218, t223, t226, t229;
490 double t230, t231, t232, t234, t235, t238, t240, t241;
491 double t242, t243, t244, t245, t246, t248, t249, t250;
492 double t252, t253, t254, t255, t257, t258, t259, t261;
493 double t262, t263, t265, t266, t267, t268, t270, t273;
494 double t275, t277, t280, t281, t282, t284, t288, t289;
495 double t290, t291, t293, t296, t298, t299, t301, t304;
496 double t306, t309, t311, t314, t315, t317, t318, t319;
497 double t320, t321, t323, t324, t325, t326, t328, t329;
498 double t330, t331, t333, t334, t335, t337, t338, t339;
499 double t340, t342, t347, t350, t353, t354, t355, t356;
500 double t358, t359, t362, t364, t365, t366, t369, t371;
501 double t372, t383, t384, t386, t389, t394, t395, t396;
502 double t401, t406, t416, t417, t419, t422, t423, t426;
503 double t427, t430, t431, t432, t434, t439, t440, t443;
504 double t444, t447, t448, t449, t451, t456, t457;
505
506 #ifndef XC_DONT_COMPILE_KXC
507 double t464, t466, t470, t471, t472, t473, t475, t476;
508 double t477, t479, t481, t483, t484, t485, t487, t488;
509 double t490, t491, t492, t493, t495, t496, t497, t499;
510 double t501, t503, t504, t505, t507, t508, t509, t511;
511 double t512, t514, t517, t519, t521, t522, t523, t525;
512 double t526, t527, t528, t529, t531, t532, t534, t535;
513 double t536, t538, t539, t541, t543, t544, t546, t547;
514 double t548, t550, t551, t553, t554, t556, t559, t564;
515 double t570, t571, t575, t576, t577, t579, t580, t581;
516 double t585, t587, t588, t589, t594, t595, t598, t599;
517 double t600, t602, t604, t606, t607, t608, t609, t611;
518 double t613, t614, t615, t617, t619, t621, t622, t623;
519 double t625, t627, t629, t631, t633, t635, t636, t638;
520 double t643, t646, t649, t651, t653, t656, t657, t659;
521 double t661, t663, t665, t667, t669, t672, t673, t675;
522 double t676, t677, t679, t681, t683, t685, t687, t696;
523 double t699, t704, t706, t708, t710, t712, t714, t716;
524 double t718, t720, t722, t724, t726, t728, t730, t735;
525 double t738, t741, t743, t745, t748, t750, t752, t754;
526 double t756, t759, t760, t762, t764, t766, t767, t768;
527 double t770, t772, t774, t777, t785, t788, t796, t797;
528 double t802, t807, t813, t814, t818, t819, t820, t822;
529 double t823, t824, t828, t830, t831, t832, t836, t837;
530 double t839, t840, t841, t843, t844, t845, t847, t848;
531 double t849, t851, t852, t853, t854, t856, t857, t859;
532 double t860, t862, t863, t864, t866, t867, t868, t869;
533 double t871, t872, t873, t875, t876, t877, t881, t885;
534 double t887, t888, t892, t895, t911, t912, t915, t920;
535 double t924, t925, t929, t935, t936, t937, t939, t940;
536 double t941, t942, t946, t949, t950, t952, t956, t958;
537 double t959, t963, t966, t969, t974, t978, t983, t987;
538 double t990, t1005, t1006, t1009, t1010, t1013, t1016, t1023;
539 double t1028, t1031, t1032, t1035, t1037, t1041, t1046, t1053;
540 double t1058, t1061, t1062, t1065, t1067, t1070, t1071, t1075;
541 double t1076, t1077, t1080, t1081, t1084, t1086, t1089, t1090;
542 double t1094, t1095, t1096, t1099, t1100;
543
544 #ifndef XC_DONT_COMPILE_LXC
545 double t1114, t1118, t1121, t1129, t1143, t1152, t1156, t1158;
546 double t1161, t1164, t1167, t1180, t1191, t1199, t1205, t1207;
547 double t1209, t1210, t1211, t1216, t1218, t1220, t1222, t1224;
548 double t1229, t1231, t1238, t1241, t1244, t1247, t1253, t1256;
549 double t1259, t1273, t1275, t1277, t1283, t1290, t1292, t1293;
550 double t1294, t1299, t1301, t1303, t1305, t1307, t1310, t1319;
551 double t1328, t1331, t1334, t1336, t1341, t1348, t1351, t1352;
552 double t1366, t1367, t1368, t1370, t1372, t1383, t1384, t1393;
553 double t1395, t1399, t1406, t1412, t1414, t1417, t1423, t1430;
554 double t1435, t1438, t1446, t1452, t1455, t1458, t1459, t1461;
555 double t1469, t1471, t1473, t1476, t1496, t1501, t1503, t1506;
556 double t1512, t1520, t1523, t1537, t1540, t1547, t1553, t1556;
557 double t1559, t1561, t1567, t1570, t1573, t1596, t1602, t1605;
558 double t1608, t1615, t1618, t1625, t1632, t1639, t1655, t1666;
559 double t1671, t1694, t1704, t1706, t1725, t1741, t1751, t1762;
560 double t1775, t1790, t1793, t1796, t1799, t1802, t1805, t1808;
561 double t1816, t1838, t1858, t1899, t1906, t1917, t1925, t1932;
562 double t1939, t1943, t1951, t1955, t1986, t1991, t1996, t2001;
563 double t2007, t2035, t2050, t2063, t2074, t2085, t2089, t2101;
564 double t2103, t2106, t2124, t2127, t2131, t2141, t2144, t2170;
565 double t2173, t2176, t2180, t2193, t2194, t2207, t2208, t2210;
566 double t2212, t2223, t2226, t2230, t2239, t2242, t2243, t2245;
567 double t2258, t2261, t2274, t2286, t2303, t2311, t2320, t2329;
568 double t2402, t2403, t2405, t2436, t2437, t2449, t2464, t2474;
569 double t2519, t2543, t2575, t2599;
570 #endif
571
572 #endif
573
574 #endif
575
576 #endif
577
578 #endif
579
580
581 gga_x_ol2_params *params;
582
583 assert(p->params != NULL);
584 params = (gga_x_ol2_params * )(p->params);
585
586 t1 = M_CBRT3;
587 t3 = POW_1_3(0.1e1 / M_PI);
588 t4 = t1 * t3;
589 t5 = M_CBRT4;
590 t6 = t5 * t5;
591 t7 = t4 * t6;
592 t8 = rho[0] - rho[1];
593 t9 = rho[0] + rho[1];
594 t10 = 0.1e1 / t9;
595 t11 = t8 * t10;
596 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
597 t14 = POW_1_3(t13);
598 t15 = t14 * t13;
599 t16 = POW_1_3(t9);
600 t17 = t15 * t16;
601 t18 = params->bb * sigma[0];
602 t19 = rho[0] * rho[0];
603 t20 = POW_1_3(rho[0]);
604 t21 = t20 * t20;
605 t23 = 0.1e1 / t21 / t19;
606 t26 = sqrt(sigma[0]);
607 t27 = params->cc * t26;
608 t29 = 0.1e1 / t20 / rho[0];
609 t30 = M_CBRT2;
610 t33 = 0.4e1 * t26 * t29 + t30;
611 t34 = 0.1e1 / t33;
612 t35 = t29 * t34;
613 t37 = params->aa + 0.13888888888888888889e-1 * t18 * t23 + t27 * t35;
614 t39 = t7 * t17 * t37;
615 t41 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
616 t42 = POW_1_3(t41);
617 t43 = t42 * t41;
618 t44 = t43 * t16;
619 t45 = params->bb * sigma[2];
620 t46 = rho[1] * rho[1];
621 t47 = POW_1_3(rho[1]);
622 t48 = t47 * t47;
623 t50 = 0.1e1 / t48 / t46;
624 t53 = sqrt(sigma[2]);
625 t54 = params->cc * t53;
626 t56 = 0.1e1 / t47 / rho[1];
627 t59 = 0.4e1 * t53 * t56 + t30;
628 t60 = 0.1e1 / t59;
629 t61 = t56 * t60;
630 t63 = params->aa + 0.13888888888888888889e-1 * t45 * t50 + t54 * t61;
631 t65 = t7 * t44 * t63;
632 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
633 zk[0] = -0.3e1 / 0.8e1 * t39 - 0.3e1 / 0.8e1 * t65;
634
635 #ifndef XC_DONT_COMPILE_VXC
636
637 if(order < 1) return;
638
639
640 t67 = 0.3e1 / 0.8e1 * t39;
641 t68 = 0.3e1 / 0.8e1 * t65;
642 t69 = t14 * t16;
643 t70 = t9 * t9;
644 t71 = 0.1e1 / t70;
645 t72 = t8 * t71;
646 t74 = t10 / 0.2e1 - t72 / 0.2e1;
647 t75 = t37 * t74;
648 t77 = t7 * t69 * t75;
649 t78 = t77 / 0.2e1;
650 t79 = t16 * t16;
651 t80 = 0.1e1 / t79;
652 t81 = t15 * t80;
653 t83 = t7 * t81 * t37;
654 t84 = t83 / 0.8e1;
655 t85 = t19 * rho[0];
656 t87 = 0.1e1 / t21 / t85;
657 t92 = 0.1e1 / t20 / t19 * t34;
658 t95 = params->cc * sigma[0];
659 t96 = t33 * t33;
660 t97 = 0.1e1 / t96;
661 t98 = t87 * t97;
662 t101 = -0.37037037037037037037e-1 * t18 * t87 - 0.4e1 / 0.3e1 * t27 * t92 + 0.16e2 / 0.3e1 * t95 * t98;
663 t103 = t7 * t17 * t101;
664 t104 = 0.3e1 / 0.8e1 * t103;
665 t105 = t42 * t16;
666 t106 = -t74;
667 t107 = t63 * t106;
668 t109 = t7 * t105 * t107;
669 t110 = t109 / 0.2e1;
670 t111 = t43 * t80;
671 t113 = t7 * t111 * t63;
672 t114 = t113 / 0.8e1;
673 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
674 vrho[0] = -t67 - t68 + t9 * (-t78 - t84 - t104 - t110 - t114);
675
676 t118 = -t10 / 0.2e1 - t72 / 0.2e1;
677 t119 = t37 * t118;
678 t121 = t7 * t69 * t119;
679 t122 = t121 / 0.2e1;
680 t123 = -t118;
681 t124 = t63 * t123;
682 t126 = t7 * t105 * t124;
683 t127 = t126 / 0.2e1;
684 t128 = t46 * rho[1];
685 t130 = 0.1e1 / t48 / t128;
686 t135 = 0.1e1 / t47 / t46 * t60;
687 t138 = params->cc * sigma[2];
688 t139 = t59 * t59;
689 t140 = 0.1e1 / t139;
690 t141 = t130 * t140;
691 t144 = -0.37037037037037037037e-1 * t45 * t130 - 0.4e1 / 0.3e1 * t54 * t135 + 0.16e2 / 0.3e1 * t138 * t141;
692 t146 = t7 * t44 * t144;
693 t147 = 0.3e1 / 0.8e1 * t146;
694 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
695 vrho[1] = -t67 - t68 + t9 * (-t122 - t84 - t127 - t114 - t147);
696
697 t151 = t16 * t9 * t1;
698 t152 = t151 * t3;
699 t153 = t6 * t15;
700 t156 = 0.1e1 / t26;
701 t157 = params->cc * t156;
702 t163 = 0.13888888888888888889e-1 * params->bb * t23 + t157 * t35 / 0.2e1 - 0.2e1 * params->cc * t23 * t97;
703 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
704 vsigma[0] = -0.3e1 / 0.8e1 * t152 * t153 * t163;
705
706 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
707 vsigma[1] = 0.0e0;
708
709 t167 = t6 * t43;
710 t170 = 0.1e1 / t53;
711 t171 = params->cc * t170;
712 t177 = 0.13888888888888888889e-1 * params->bb * t50 + t171 * t61 / 0.2e1 - 0.2e1 * params->cc * t50 * t140;
713 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
714 vsigma[2] = -0.3e1 / 0.8e1 * t152 * t167 * t177;
715
716 #ifndef XC_DONT_COMPILE_FXC
717
718 if(order < 2) return;
719
720
721 t181 = t83 / 0.4e1;
722 t183 = t113 / 0.4e1;
723 t184 = t14 * t14;
724 t185 = 0.1e1 / t184;
725 t186 = t185 * t16;
726 t187 = t74 * t74;
727 t188 = t37 * t187;
728 t190 = t7 * t186 * t188;
729 t191 = t190 / 0.6e1;
730 t192 = t14 * t80;
731 t194 = t7 * t192 * t75;
732 t195 = t194 / 0.3e1;
733 t196 = t101 * t74;
734 t198 = t7 * t69 * t196;
735 t199 = t70 * t9;
736 t200 = 0.1e1 / t199;
737 t201 = t8 * t200;
738 t202 = -t71 + t201;
739 t203 = t37 * t202;
740 t205 = t7 * t69 * t203;
741 t206 = t205 / 0.2e1;
742 t208 = 0.1e1 / t79 / t9;
743 t209 = t15 * t208;
744 t211 = t7 * t209 * t37;
745 t212 = t211 / 0.12e2;
746 t214 = t7 * t81 * t101;
747 t215 = t214 / 0.4e1;
748 t216 = t19 * t19;
749 t218 = 0.1e1 / t21 / t216;
750 t223 = 0.1e1 / t20 / t85 * t34;
751 t226 = t218 * t97;
752 t229 = t26 * sigma[0];
753 t230 = params->cc * t229;
754 t231 = t216 * t19;
755 t232 = 0.1e1 / t231;
756 t234 = 0.1e1 / t96 / t33;
757 t235 = t232 * t234;
758 t238 = 0.13580246913580246914e0 * t18 * t218 + 0.28e2 / 0.9e1 * t27 * t223 - 0.80e2 / 0.3e1 * t95 * t226 + 0.512e3 / 0.9e1 * t230 * t235;
759 t240 = t7 * t17 * t238;
760 t241 = 0.3e1 / 0.8e1 * t240;
761 t242 = t42 * t42;
762 t243 = 0.1e1 / t242;
763 t244 = t243 * t16;
764 t245 = t106 * t106;
765 t246 = t63 * t245;
766 t248 = t7 * t244 * t246;
767 t249 = t248 / 0.6e1;
768 t250 = t42 * t80;
769 t252 = t7 * t250 * t107;
770 t253 = t252 / 0.3e1;
771 t254 = -t202;
772 t255 = t63 * t254;
773 t257 = t7 * t105 * t255;
774 t258 = t257 / 0.2e1;
775 t259 = t43 * t208;
776 t261 = t7 * t259 * t63;
777 t262 = t261 / 0.12e2;
778 t263 = -t191 - t195 - t198 - t206 + t212 - t215 - t241 - t249 - t253 - t258 + t262;
779 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
780 v2rho2[0] = -t77 - t181 - 0.3e1 / 0.4e1 * t103 - t109 - t183 + t9 * t263;
781
782 t265 = t6 * t185;
783 t266 = t4 * t265;
784 t267 = t16 * t37;
785 t268 = t118 * t74;
786 t270 = t266 * t267 * t268;
787 t273 = t7 * t192 * t119;
788 t275 = t101 * t118;
789 t277 = t7 * t69 * t275;
790 t280 = 0.1e1 / t79 / t70;
791 t281 = t14 * t280;
792 t282 = t37 * t8;
793 t284 = t7 * t281 * t282;
794 t288 = t6 * t243;
795 t289 = t4 * t288;
796 t290 = t16 * t63;
797 t291 = t123 * t106;
798 t293 = t289 * t290 * t291;
799 t296 = t7 * t250 * t124;
800 t298 = t42 * t280;
801 t299 = t63 * t8;
802 t301 = t7 * t298 * t299;
803 t304 = t144 * t106;
804 t306 = t7 * t105 * t304;
805 t309 = t7 * t111 * t144;
806 t311 = -t270 / 0.6e1 - t273 / 0.6e1 - t277 / 0.2e1 - t284 / 0.2e1 - t194 / 0.6e1 + t212 - t214 / 0.8e1 - t293 / 0.6e1 - t296 / 0.6e1 + t301 / 0.2e1 - t252 / 0.6e1 + t262 - t306 / 0.2e1 - t309 / 0.8e1;
807 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
808 v2rho2[1] = t9 * t311 - t104 - t110 - t122 - t127 - t147 - t181 - t183 - t78;
809
810 t314 = t118 * t118;
811 t315 = t37 * t314;
812 t317 = t7 * t186 * t315;
813 t318 = t317 / 0.6e1;
814 t319 = t273 / 0.3e1;
815 t320 = t71 + t201;
816 t321 = t37 * t320;
817 t323 = t7 * t69 * t321;
818 t324 = t323 / 0.2e1;
819 t325 = t123 * t123;
820 t326 = t63 * t325;
821 t328 = t7 * t244 * t326;
822 t329 = t328 / 0.6e1;
823 t330 = t296 / 0.3e1;
824 t331 = t144 * t123;
825 t333 = t7 * t105 * t331;
826 t334 = -t320;
827 t335 = t63 * t334;
828 t337 = t7 * t105 * t335;
829 t338 = t337 / 0.2e1;
830 t339 = t309 / 0.4e1;
831 t340 = t46 * t46;
832 t342 = 0.1e1 / t48 / t340;
833 t347 = 0.1e1 / t47 / t128 * t60;
834 t350 = t342 * t140;
835 t353 = t53 * sigma[2];
836 t354 = params->cc * t353;
837 t355 = t340 * t46;
838 t356 = 0.1e1 / t355;
839 t358 = 0.1e1 / t139 / t59;
840 t359 = t356 * t358;
841 t362 = 0.13580246913580246914e0 * t45 * t342 + 0.28e2 / 0.9e1 * t54 * t347 - 0.80e2 / 0.3e1 * t138 * t350 + 0.512e3 / 0.9e1 * t354 * t359;
842 t364 = t7 * t44 * t362;
843 t365 = 0.3e1 / 0.8e1 * t364;
844 t366 = -t318 - t319 - t324 + t212 - t329 - t330 - t333 - t338 + t262 - t339 - t365;
845 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
846 v2rho2[2] = -t121 - t181 - t126 - t183 - 0.3e1 / 0.4e1 * t146 + t9 * t366;
847
848 t369 = t7 * t17 * t163;
849 t371 = t6 * t14;
850 t372 = t163 * t74;
851 t383 = t216 * rho[0];
852 t384 = 0.1e1 / t383;
853 t386 = t234 * t26;
854 t389 = -0.37037037037037037037e-1 * params->bb * t87 - 0.2e1 / 0.3e1 * t157 * t92 + 0.8e1 * params->cc * t87 * t97 - 0.64e2 / 0.3e1 * params->cc * t384 * t386;
855 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
856 v2rhosigma[0] = -t369 / 0.2e1 - t152 * t371 * t372 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * t389;
857
858 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
859 v2rhosigma[1] = 0.0e0;
860
861 t394 = t7 * t44 * t177;
862 t395 = t6 * t42;
863 t396 = t177 * t106;
864 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
865 v2rhosigma[2] = -t152 * t395 * t396 / 0.2e1 - t394 / 0.2e1;
866
867 t401 = t371 * t163 * t118;
868 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
869 v2rhosigma[3] = -t152 * t401 / 0.2e1 - t369 / 0.2e1;
870
871 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
872 v2rhosigma[4] = 0.0e0;
873
874 t406 = t395 * t177 * t123;
875 t416 = t340 * rho[1];
876 t417 = 0.1e1 / t416;
877 t419 = t358 * t53;
878 t422 = -0.37037037037037037037e-1 * params->bb * t130 - 0.2e1 / 0.3e1 * t171 * t135 + 0.8e1 * params->cc * t130 * t140 - 0.64e2 / 0.3e1 * params->cc * t417 * t419;
879 t423 = t167 * t422;
880 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
881 v2rhosigma[5] = -t394 / 0.2e1 - t152 * t406 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t423;
882
883 t426 = 0.1e1 / t229;
884 t427 = params->cc * t426;
885 t430 = 0.1e1 / sigma[0];
886 t431 = params->cc * t430;
887 t432 = t23 * t97;
888 t434 = 0.1e1 / t216;
889 t439 = -t427 * t35 / 0.4e1 - t431 * t432 + 0.8e1 * params->cc * t434 * t234 * t156;
890 t440 = t153 * t439;
891 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
892 v2sigma2[0] = -0.3e1 / 0.8e1 * t152 * t440;
893
894 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
895 v2sigma2[1] = 0.0e0;
896
897 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
898 v2sigma2[2] = 0.0e0;
899
900 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
901 v2sigma2[3] = 0.0e0;
902
903 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
904 v2sigma2[4] = 0.0e0;
905
906 t443 = 0.1e1 / t353;
907 t444 = params->cc * t443;
908 t447 = 0.1e1 / sigma[2];
909 t448 = params->cc * t447;
910 t449 = t50 * t140;
911 t451 = 0.1e1 / t340;
912 t456 = -t444 * t61 / 0.4e1 - t448 * t449 + 0.8e1 * params->cc * t451 * t358 * t170;
913 t457 = t167 * t456;
914 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
915 v2sigma2[5] = -0.3e1 / 0.8e1 * t152 * t457;
916
917 #ifndef XC_DONT_COMPILE_KXC
918
919 if(order < 3) return;
920
921
922 t464 = t211 / 0.4e1;
923 t466 = t261 / 0.4e1;
924 t470 = 0.1e1 / t184 / t13;
925 t471 = t470 * t16;
926 t472 = t187 * t74;
927 t473 = t37 * t472;
928 t475 = t7 * t471 * t473;
929 t476 = t475 / 0.9e1;
930 t477 = t185 * t80;
931 t479 = t7 * t477 * t188;
932 t481 = t74 * t202;
933 t483 = t266 * t267 * t481;
934 t484 = t483 / 0.2e1;
935 t485 = t14 * t208;
936 t487 = t7 * t485 * t75;
937 t488 = t487 / 0.3e1;
938 t490 = 0.1e1 / t242 / t41;
939 t491 = t490 * t16;
940 t492 = t245 * t106;
941 t493 = t63 * t492;
942 t495 = t7 * t491 * t493;
943 t496 = t495 / 0.9e1;
944 t497 = t243 * t80;
945 t499 = t7 * t497 * t246;
946 t501 = t106 * t254;
947 t503 = t289 * t290 * t501;
948 t504 = t503 / 0.2e1;
949 t505 = t42 * t208;
950 t507 = t7 * t505 * t107;
951 t508 = t507 / 0.3e1;
952 t509 = t101 * t187;
953 t511 = t7 * t186 * t509;
954 t512 = t511 / 0.2e1;
955 t514 = t7 * t192 * t196;
956 t517 = t7 * t192 * t203;
957 t519 = t238 * t74;
958 t521 = t7 * t69 * t519;
959 t522 = 0.3e1 / 0.2e1 * t521;
960 t523 = t101 * t202;
961 t525 = t7 * t69 * t523;
962 t526 = 0.3e1 / 0.2e1 * t525;
963 t527 = t70 * t70;
964 t528 = 0.1e1 / t527;
965 t529 = t8 * t528;
966 t531 = 0.3e1 * t200 - 0.3e1 * t529;
967 t532 = t37 * t531;
968 t534 = t7 * t69 * t532;
969 t535 = t534 / 0.2e1;
970 t536 = t15 * t280;
971 t538 = t7 * t536 * t37;
972 t539 = 0.5e1 / 0.36e2 * t538;
973 t541 = t7 * t250 * t255;
974 t543 = -t531;
975 t544 = t63 * t543;
976 t546 = t7 * t105 * t544;
977 t547 = t546 / 0.2e1;
978 t548 = t43 * t280;
979 t550 = t7 * t548 * t63;
980 t551 = 0.5e1 / 0.36e2 * t550;
981 t553 = t7 * t209 * t101;
982 t554 = t553 / 0.4e1;
983 t556 = t7 * t81 * t238;
984 t559 = 0.1e1 / t21 / t383;
985 t564 = 0.1e1 / t20 / t216 * t34;
986 t570 = t216 * t85;
987 t571 = 0.1e1 / t570;
988 t575 = sigma[0] * sigma[0];
989 t576 = params->cc * t575;
990 t577 = t216 * t216;
991 t579 = 0.1e1 / t20 / t577;
992 t580 = t96 * t96;
993 t581 = 0.1e1 / t580;
994 t585 = -0.63374485596707818932e0 * t18 * t559 - 0.280e3 / 0.27e2 * t27 * t564 + 0.3808e4 / 0.27e2 * t95 * t559 * t97 - 0.5632e4 / 0.9e1 * t230 * t571 * t234 + 0.8192e4 / 0.9e1 * t576 * t579 * t581;
995 t587 = t7 * t17 * t585;
996 t588 = 0.3e1 / 0.8e1 * t587;
997 t589 = -t517 / 0.2e1 - t522 - t526 - t535 - t539 - t541 / 0.2e1 - t547 - t551 + t554 - 0.3e1 / 0.8e1 * t556 - t588;
998 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
999 v3rho3[0] = -t190 / 0.2e1 - t194 - t248 / 0.2e1 - t252 - 0.3e1 * t198 - 0.3e1 / 0.2e1 * t205 + t464 - 0.3e1 / 0.2e1 * t257 + t466 - 0.3e1 / 0.4e1 * t214 - 0.9e1 / 0.8e1 * t240 + t9 * (t476 - t479 / 0.6e1 - t484 + t488 + t496 - t499 / 0.6e1 - t504 + t508 - t512 - t514 + t589);
1000
1001 t594 = t270 / 0.3e1;
1002 t595 = t293 / 0.3e1;
1003 t598 = t6 * t470;
1004 t599 = t4 * t598;
1005 t600 = t118 * t187;
1006 t602 = t599 * t267 * t600;
1007 t604 = t80 * t37;
1008 t606 = t266 * t604 * t268;
1009 t607 = t606 / 0.9e1;
1010 t608 = t280 * t37;
1011 t609 = t8 * t74;
1012 t611 = t266 * t608 * t609;
1013 t613 = t6 * t490;
1014 t614 = t4 * t613;
1015 t615 = t123 * t245;
1016 t617 = t614 * t290 * t615;
1017 t619 = t80 * t63;
1018 t621 = t289 * t619 * t291;
1019 t622 = t621 / 0.9e1;
1020 t623 = t16 * t101;
1021 t625 = t266 * t623 * t268;
1022 t627 = t118 * t202;
1023 t629 = t266 * t267 * t627;
1024 t631 = t123 * t254;
1025 t633 = t289 * t290 * t631;
1026 t635 = t280 * t63;
1027 t636 = t8 * t106;
1028 t638 = t289 * t635 * t636;
1029 t643 = t7 * t281 * t37;
1030 t646 = t7 * t298 * t63;
1031 t649 = t7 * t259 * t144;
1032 t651 = t602 / 0.9e1 - t607 - t611 / 0.3e1 + t617 / 0.9e1 - t622 - t625 / 0.3e1 - t629 / 0.6e1 - t633 / 0.6e1 + t638 / 0.3e1 + t553 / 0.6e1 - t556 / 0.8e1 - t643 / 0.2e1 + t646 / 0.2e1 + t649 / 0.12e2 - t539 - t551;
1033 t653 = t7 * t485 * t119;
1034 t656 = 0.1e1 / t79 / t199;
1035 t657 = t14 * t656;
1036 t659 = t7 * t657 * t282;
1037 t661 = t101 * t8;
1038 t663 = t7 * t281 * t661;
1039 t665 = t7 * t505 * t124;
1040 t667 = t144 * t245;
1041 t669 = t7 * t244 * t667;
1042 t672 = t7 * t250 * t304;
1043 t673 = t672 / 0.3e1;
1044 t675 = t7 * t192 * t275;
1045 t676 = t675 / 0.3e1;
1046 t677 = t238 * t118;
1047 t679 = t7 * t69 * t677;
1048 t681 = t42 * t656;
1049 t683 = t7 * t681 * t299;
1050 t685 = t144 * t254;
1051 t687 = t7 * t105 * t685;
1052 t696 = t653 / 0.9e1 + 0.7e1 / 0.6e1 * t659 - t663 + t665 / 0.9e1 - t669 / 0.6e1 - t673 - t676 - t679 / 0.2e1 - 0.7e1 / 0.6e1 * t683 - t687 / 0.2e1 - t479 / 0.18e2 + 0.2e1 / 0.9e1 * t487 - t499 / 0.18e2 + 0.2e1 / 0.9e1 * t507 - t514 / 0.3e1 - t517 / 0.6e1 - t541 / 0.6e1;
1053 t699 = -t206 + t464 - t258 + t466 - t277 + t301 - t306 - t339 - t214 / 0.2e1 - t241 + t9 * (t651 + t696);
1054 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1055 v3rho3[1] = -t191 - 0.2e1 / 0.3e1 * t194 - t249 - 0.2e1 / 0.3e1 * t252 - t594 - t319 - t284 - t595 - t330 - t198 + t699;
1056
1057 t704 = t320 * t74;
1058 t706 = t266 * t267 * t704;
1059 t708 = t325 * t106;
1060 t710 = t614 * t290 * t708;
1061 t712 = t118 * t8;
1062 t714 = t266 * t608 * t712;
1063 t716 = t123 * t8;
1064 t718 = t289 * t635 * t716;
1065 t720 = t16 * t144;
1066 t722 = t289 * t720 * t291;
1067 t724 = t334 * t106;
1068 t726 = t289 * t290 * t724;
1069 t728 = t314 * t74;
1070 t730 = t599 * t267 * t728;
1071 t735 = t7 * t111 * t362;
1072 t738 = t7 * t250 * t335;
1073 t741 = t7 * t477 * t315;
1074 t743 = -t706 / 0.6e1 + t710 / 0.9e1 - t714 / 0.3e1 + t718 / 0.3e1 - t722 / 0.3e1 - t726 / 0.6e1 - t607 - t622 + t730 / 0.9e1 + t553 / 0.12e2 + t649 / 0.6e1 - t735 / 0.8e1 - t539 - t551 - t738 / 0.6e1 - t741 / 0.18e2;
1075 t745 = t7 * t192 * t321;
1076 t748 = t7 * t497 * t326;
1077 t750 = t101 * t314;
1078 t752 = t7 * t186 * t750;
1079 t754 = t101 * t320;
1080 t756 = t7 * t69 * t754;
1081 t759 = -t200 - 0.3e1 * t529;
1082 t760 = t37 * t759;
1083 t762 = t7 * t69 * t760;
1084 t764 = t144 * t8;
1085 t766 = t7 * t298 * t764;
1086 t767 = -t759;
1087 t768 = t63 * t767;
1088 t770 = t7 * t105 * t768;
1089 t772 = t362 * t106;
1090 t774 = t7 * t105 * t772;
1091 t777 = t7 * t250 * t331;
1092 t785 = -t745 / 0.6e1 - t748 / 0.18e2 - t752 / 0.6e1 - t756 / 0.2e1 - t762 / 0.2e1 + t766 - t770 / 0.2e1 - t774 / 0.2e1 - t777 / 0.3e1 + 0.2e1 / 0.9e1 * t653 - t659 / 0.3e1 + 0.2e1 / 0.9e1 * t665 - t673 - t676 + t683 / 0.3e1 + t487 / 0.9e1 + t507 / 0.9e1;
1093 t788 = t466 - t277 + t301 - t306 - t309 / 0.2e1 - t324 - t333 - t338 - t215 - t365 + t9 * (t743 + t785);
1094 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1095 v3rho3[2] = -t195 - t253 - t594 - 0.2e1 / 0.3e1 * t273 - t284 - t595 - 0.2e1 / 0.3e1 * t296 - t318 - t329 + t464 + t788;
1096
1097 t796 = t653 / 0.3e1;
1098 t797 = t665 / 0.3e1;
1099 t802 = 0.1e1 / t48 / t416;
1100 t807 = 0.1e1 / t47 / t340 * t60;
1101 t813 = t340 * t128;
1102 t814 = 0.1e1 / t813;
1103 t818 = sigma[2] * sigma[2];
1104 t819 = params->cc * t818;
1105 t820 = t340 * t340;
1106 t822 = 0.1e1 / t47 / t820;
1107 t823 = t139 * t139;
1108 t824 = 0.1e1 / t823;
1109 t828 = -0.63374485596707818932e0 * t45 * t802 - 0.280e3 / 0.27e2 * t54 * t807 + 0.3808e4 / 0.27e2 * t138 * t802 * t140 - 0.5632e4 / 0.9e1 * t354 * t814 * t358 + 0.8192e4 / 0.9e1 * t819 * t822 * t824;
1110 t830 = t7 * t44 * t828;
1111 t831 = 0.3e1 / 0.8e1 * t830;
1112 t832 = t649 / 0.4e1;
1113 t836 = -0.3e1 * t200 - 0.3e1 * t529;
1114 t837 = t37 * t836;
1115 t839 = t7 * t69 * t837;
1116 t840 = t839 / 0.2e1;
1117 t841 = t144 * t325;
1118 t843 = t7 * t244 * t841;
1119 t844 = t843 / 0.2e1;
1120 t845 = t362 * t123;
1121 t847 = t7 * t105 * t845;
1122 t848 = 0.3e1 / 0.2e1 * t847;
1123 t849 = t144 * t334;
1124 t851 = t7 * t105 * t849;
1125 t852 = 0.3e1 / 0.2e1 * t851;
1126 t853 = -t836;
1127 t854 = t63 * t853;
1128 t856 = t7 * t105 * t854;
1129 t857 = t856 / 0.2e1;
1130 t859 = t314 * t118;
1131 t860 = t37 * t859;
1132 t862 = t7 * t471 * t860;
1133 t863 = t862 / 0.9e1;
1134 t864 = t118 * t320;
1135 t866 = t266 * t267 * t864;
1136 t867 = t866 / 0.2e1;
1137 t868 = t325 * t123;
1138 t869 = t63 * t868;
1139 t871 = t7 * t491 * t869;
1140 t872 = t871 / 0.9e1;
1141 t873 = t123 * t334;
1142 t875 = t289 * t290 * t873;
1143 t876 = t875 / 0.2e1;
1144 t877 = -t840 - t844 - t848 - t852 - t857 - t777 - t738 / 0.2e1 + t863 - t867 + t872 - t876;
1145 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1146 v3rho3[3] = -t273 - t296 - t317 / 0.2e1 - t328 / 0.2e1 + t464 + t466 - 0.3e1 / 0.4e1 * t309 - 0.3e1 / 0.2e1 * t323 - 0.3e1 * t333 - 0.3e1 / 0.2e1 * t337 - 0.9e1 / 0.8e1 * t364 + t9 * (t796 + t797 - t741 / 0.6e1 - t745 / 0.2e1 - t748 / 0.6e1 - t539 - t551 - t831 + t832 - 0.3e1 / 0.8e1 * t735 + t877);
1147
1148 t881 = t7 * t69 * t372;
1149 t885 = t7 * t81 * t163 / 0.6e1;
1150 t887 = t7 * t17 * t389;
1151 t888 = t163 * t187;
1152 t892 = t389 * t74;
1153 t895 = t163 * t202;
1154 t911 = params->cc / t20 / t570;
1155 t912 = t581 * sigma[0];
1156 t915 = 0.13580246913580246914e0 * params->bb * t218 + 0.14e2 / 0.9e1 * t157 * t223 - 0.296e3 / 0.9e1 * params->cc * t218 * t97 + 0.192e3 * params->cc * t232 * t386 - 0.1024e4 / 0.3e1 * t911 * t912;
1157 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1158 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t881 - t885 - t887 - t152 * t265 * t888 / 0.6e1 - t152 * t371 * t892 - t152 * t371 * t895 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * t915;
1159
1160 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1161 v3rho2sigma[1] = 0.0e0;
1162
1163 t920 = t7 * t105 * t396;
1164 t924 = t7 * t111 * t177 / 0.6e1;
1165 t925 = t177 * t245;
1166 t929 = t177 * t254;
1167 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1168 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t920 - t924 - t152 * t288 * t925 / 0.6e1 - t152 * t395 * t929 / 0.2e1;
1169
1170 t935 = t16 * t1;
1171 t936 = t935 * t3;
1172 t937 = t936 * t401;
1173 t939 = t3 * t6;
1174 t940 = t151 * t939;
1175 t941 = t185 * t163;
1176 t942 = t941 * t268;
1177 t946 = t371 * t389 * t118;
1178 t949 = t208 * t1;
1179 t950 = t949 * t3;
1180 t952 = t371 * t163 * t8;
1181 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1182 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t881 - t885 - t887 / 0.2e1 - 0.2e1 / 0.3e1 * t937 - t940 * t942 / 0.6e1 - t152 * t946 / 0.2e1 - t950 * t952 / 0.2e1;
1183
1184 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1185 v3rho2sigma[4] = 0.0e0;
1186
1187 t956 = t936 * t406;
1188 t958 = t243 * t177;
1189 t959 = t958 * t291;
1190 t963 = t395 * t177 * t8;
1191 t966 = t936 * t423;
1192 t969 = t395 * t422 * t106;
1193 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1194 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t920 - t924 - 0.2e1 / 0.3e1 * t956 - t940 * t959 / 0.6e1 + t950 * t963 / 0.2e1 - t966 / 0.2e1 - t152 * t969 / 0.2e1;
1195
1196 t974 = t265 * t163 * t314;
1197 t978 = t371 * t163 * t320;
1198 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1199 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t937 - t885 - t152 * t974 / 0.6e1 - t152 * t978 / 0.2e1;
1200
1201 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1202 v3rho2sigma[7] = 0.0e0;
1203
1204 t983 = t288 * t177 * t325;
1205 t987 = t395 * t422 * t123;
1206 t990 = t395 * t177 * t334;
1207 t1005 = params->cc / t47 / t813;
1208 t1006 = t824 * sigma[2];
1209 t1009 = 0.13580246913580246914e0 * params->bb * t342 + 0.14e2 / 0.9e1 * t171 * t347 - 0.296e3 / 0.9e1 * params->cc * t342 * t140 + 0.192e3 * params->cc * t356 * t419 - 0.1024e4 / 0.3e1 * t1005 * t1006;
1210 t1010 = t167 * t1009;
1211 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1212 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t956 - t924 - t966 - t152 * t983 / 0.6e1 - t152 * t987 - t152 * t990 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1010;
1213
1214 t1013 = t936 * t440;
1215 t1016 = t371 * t439 * t74;
1216 t1023 = t384 * t234;
1217 t1028 = params->cc / t20 / t231;
1218 t1031 = t427 * t92 / 0.3e1 + 0.4e1 / 0.3e1 * t431 * t98 - 0.128e3 / 0.3e1 * t157 * t1023 + 0.128e3 * t1028 * t581;
1219 t1032 = t153 * t1031;
1220 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1221 v3rhosigma2[0] = -t1013 / 0.2e1 - t152 * t1016 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1032;
1222
1223 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1224 v3rhosigma2[1] = 0.0e0;
1225
1226 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1227 v3rhosigma2[2] = 0.0e0;
1228
1229 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1230 v3rhosigma2[3] = 0.0e0;
1231
1232 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1233 v3rhosigma2[4] = 0.0e0;
1234
1235 t1035 = t936 * t457;
1236 t1037 = t395 * t456 * t106;
1237 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1238 v3rhosigma2[5] = -t152 * t1037 / 0.2e1 - t1035 / 0.2e1;
1239
1240 t1041 = t371 * t439 * t118;
1241 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1242 v3rhosigma2[6] = -t152 * t1041 / 0.2e1 - t1013 / 0.2e1;
1243
1244 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1245 v3rhosigma2[7] = 0.0e0;
1246
1247 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1248 v3rhosigma2[8] = 0.0e0;
1249
1250 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1251 v3rhosigma2[9] = 0.0e0;
1252
1253 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1254 v3rhosigma2[10] = 0.0e0;
1255
1256 t1046 = t395 * t456 * t123;
1257 t1053 = t417 * t358;
1258 t1058 = params->cc / t47 / t355;
1259 t1061 = t444 * t135 / 0.3e1 + 0.4e1 / 0.3e1 * t448 * t141 - 0.128e3 / 0.3e1 * t171 * t1053 + 0.128e3 * t1058 * t824;
1260 t1062 = t167 * t1061;
1261 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1262 v3rhosigma2[11] = -t1035 / 0.2e1 - t152 * t1046 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t1062;
1263
1264 t1065 = t26 * t575;
1265 t1067 = params->cc / t1065;
1266 t1070 = 0.1e1 / t575;
1267 t1071 = params->cc * t1070;
1268 t1075 = 0.1e1 / t20 / t383;
1269 t1076 = params->cc * t1075;
1270 t1077 = t581 * t430;
1271 t1080 = 0.3e1 / 0.8e1 * t1067 * t35 + 0.3e1 / 0.2e1 * t1071 * t432 - 0.48e2 * t1076 * t1077;
1272 t1081 = t153 * t1080;
1273 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1274 v3sigma3[0] = -0.3e1 / 0.8e1 * t152 * t1081;
1275
1276 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1277 v3sigma3[1] = 0.0e0;
1278
1279 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1280 v3sigma3[2] = 0.0e0;
1281
1282 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1283 v3sigma3[3] = 0.0e0;
1284
1285 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1286 v3sigma3[4] = 0.0e0;
1287
1288 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1289 v3sigma3[5] = 0.0e0;
1290
1291 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1292 v3sigma3[6] = 0.0e0;
1293
1294 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1295 v3sigma3[7] = 0.0e0;
1296
1297 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1298 v3sigma3[8] = 0.0e0;
1299
1300 t1084 = t53 * t818;
1301 t1086 = params->cc / t1084;
1302 t1089 = 0.1e1 / t818;
1303 t1090 = params->cc * t1089;
1304 t1094 = 0.1e1 / t47 / t416;
1305 t1095 = params->cc * t1094;
1306 t1096 = t824 * t447;
1307 t1099 = 0.3e1 / 0.8e1 * t1086 * t61 + 0.3e1 / 0.2e1 * t1090 * t449 - 0.48e2 * t1095 * t1096;
1308 t1100 = t167 * t1099;
1309 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1310 v3sigma3[9] = -0.3e1 / 0.8e1 * t152 * t1100;
1311
1312 #ifndef XC_DONT_COMPILE_LXC
1313
1314 if(order < 4) return;
1315
1316
1317 t1114 = 0.4e1 / 0.9e1 * t475 - 0.2e1 / 0.3e1 * t479 - 0.2e1 * t483 + 0.4e1 / 0.3e1 * t487 + 0.4e1 / 0.9e1 * t495 - 0.2e1 / 0.3e1 * t499 - 0.2e1 * t503 + 0.4e1 / 0.3e1 * t507 - 0.2e1 * t511 - 0.4e1 * t514 - 0.2e1 * t517;
1318 t1118 = 0.5e1 / 0.9e1 * t538;
1319 t1121 = 0.5e1 / 0.9e1 * t550;
1320 t1129 = t266 * t604 * t481;
1321 t1143 = t289 * t619 * t501;
1322 t1152 = 0.10e2 / 0.27e2 * t7 * t15 * t656 * t37;
1323 t1156 = 0.10e2 / 0.27e2 * t7 * t43 * t656 * t63;
1324 t1158 = t7 * t536 * t101;
1325 t1161 = t7 * t209 * t238;
1326 t1164 = t7 * t81 * t585;
1327 t1167 = 0.1e1 / t21 / t231;
1328 t1180 = t577 * rho[0];
1329 t1191 = 0.1e1 / t580 / t33;
1330 t1199 = t202 * t202;
1331 t1205 = t7 * t281 * t75;
1332 t1207 = t41 * t41;
1333 t1209 = 0.1e1 / t242 / t1207;
1334 t1210 = t1209 * t16;
1335 t1211 = t245 * t245;
1336 t1216 = t490 * t80;
1337 t1218 = t7 * t1216 * t493;
1338 t1220 = t243 * t208;
1339 t1222 = t7 * t1220 * t246;
1340 t1224 = t254 * t254;
1341 t1229 = 0.2e1 / 0.3e1 * t599 * t267 * t187 * t202 - 0.2e1 / 0.3e1 * t1129 - 0.2e1 * t266 * t623 * t481 - 0.2e1 / 0.3e1 * t266 * t267 * t74 * t531 + 0.2e1 / 0.3e1 * t614 * t290 * t245 * t254 - 0.2e1 / 0.3e1 * t1143 - 0.2e1 / 0.3e1 * t289 * t290 * t106 * t543 + t1152 + t1156 - 0.5e1 / 0.9e1 * t1158 + t1161 / 0.2e1 - t1164 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t17 * (0.35912208504801097395e1 * t18 * t1167 + 0.3640e4 / 0.81e2 * t27 * t1075 * t34 - 0.23072e5 / 0.27e2 * t95 * t1167 * t97 + 0.476672e6 / 0.81e2 * t230 / t577 * t234 - 0.475136e6 / 0.27e2 * t576 / t20 / t1180 * t581 + 0.524288e6 / 0.27e2 * params->cc * t1065 / t21 / t577 / t19 * t1191) - t7 * t186 * t37 * t1199 / 0.2e1 - 0.20e2 / 0.27e2 * t1205 - 0.5e1 / 0.27e2 * t7 * t1210 * t63 * t1211 + 0.4e1 / 0.27e2 * t1218 + 0.2e1 / 0.9e1 * t1222 - t7 * t244 * t63 * t1224 / 0.2e1;
1342 t1231 = t7 * t298 * t107;
1343 t1238 = t7 * t477 * t509;
1344 t1241 = t7 * t485 * t196;
1345 t1244 = t7 * t485 * t203;
1346 t1247 = t7 * t505 * t255;
1347 t1253 = t7 * t192 * t519;
1348 t1256 = t7 * t192 * t523;
1349 t1259 = t7 * t192 * t532;
1350 t1273 = t527 * t9;
1351 t1275 = t8 / t1273;
1352 t1277 = -0.12e2 * t528 + 0.12e2 * t1275;
1353 t1283 = t7 * t250 * t544;
1354 t1290 = t13 * t13;
1355 t1292 = 0.1e1 / t184 / t1290;
1356 t1293 = t1292 * t16;
1357 t1294 = t187 * t187;
1358 t1299 = t470 * t80;
1359 t1301 = t7 * t1299 * t473;
1360 t1303 = t185 * t208;
1361 t1305 = t7 * t1303 * t188;
1362 t1307 = -0.20e2 / 0.27e2 * t1231 + 0.4e1 / 0.9e1 * t7 * t471 * t101 * t472 - 0.2e1 / 0.3e1 * t1238 + 0.4e1 / 0.3e1 * t1241 + 0.2e1 / 0.3e1 * t1244 + 0.2e1 / 0.3e1 * t1247 - t7 * t186 * t238 * t187 - 0.2e1 * t1253 - 0.2e1 * t1256 - 0.2e1 / 0.3e1 * t1259 - 0.2e1 * t7 * t69 * t585 * t74 - 0.3e1 * t7 * t69 * t238 * t202 - 0.2e1 * t7 * t69 * t101 * t531 - t7 * t69 * t37 * t1277 / 0.2e1 - 0.2e1 / 0.3e1 * t1283 + t7 * t105 * t63 * t1277 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1293 * t37 * t1294 + 0.4e1 / 0.27e2 * t1301 + 0.2e1 / 0.9e1 * t1305;
1363 t1310 = -0.6e1 * t521 - 0.6e1 * t525 - 0.2e1 * t534 - t1118 - 0.2e1 * t541 - 0.2e1 * t546 - t1121 + t553 - 0.3e1 / 0.2e1 * t556 - 0.3e1 / 0.2e1 * t587 + t9 * (t1229 + t1307);
1364 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1365 v4rho4[0] = t1114 + t1310;
1366
1367 t1319 = t185 * t280;
1368 t1328 = t7 * t281 * t119;
1369 t1331 = 0.1e1 / t79 / t527;
1370 t1334 = t7 * t14 * t1331 * t282;
1371 t1336 = t243 * t280;
1372 t1341 = t7 * t298 * t124;
1373 t1348 = t7 * t497 * t667;
1374 t1351 = t7 * t505 * t304;
1375 t1352 = t1351 / 0.3e1;
1376 t1366 = t266 * t208 * t37 * t268;
1377 t1367 = t1366 / 0.9e1;
1378 t1368 = t656 * t37;
1379 t1370 = t266 * t1368 * t609;
1380 t1372 = t614 * t290 * t631 * t106 / 0.3e1 - t7 * t69 * t585 * t118 / 0.2e1 - t7 * t1319 * t75 / 0.2e1 - t7 * t105 * t144 * t543 / 0.2e1 - 0.5e1 / 0.27e2 * t1328 - 0.25e2 / 0.6e1 * t1334 + t7 * t1336 * t107 / 0.2e1 - 0.5e1 / 0.27e2 * t1341 + t7 * t491 * t144 * t492 / 0.9e1 - t1348 / 0.6e1 + t1352 - t266 * t16 * t238 * t268 / 0.2e1 - t266 * t623 * t627 / 0.2e1 - t266 * t267 * t118 * t531 / 0.6e1 + t1367 + 0.7e1 / 0.6e1 * t1370;
1381 t1383 = t289 * t208 * t63 * t291;
1382 t1384 = t1383 / 0.9e1;
1383 t1393 = t266 * t80 * t101 * t268;
1384 t1395 = t280 * t101;
1385 t1399 = t266 * t604 * t627;
1386 t1406 = t289 * t619 * t631;
1387 t1412 = t656 * t63;
1388 t1414 = t289 * t1412 * t636;
1389 t1417 = t4 * t6 * t1292;
1390 t1423 = t599 * t604 * t600;
1391 t1430 = t4 * t6 * t1209;
1392 t1435 = -t289 * t290 * t123 * t543 / 0.6e1 + t289 * t635 * t8 * t254 / 0.2e1 + t1384 - t289 * t720 * t501 / 0.2e1 + t599 * t623 * t600 / 0.3e1 - t1393 / 0.3e1 - t266 * t1395 * t609 - t1399 / 0.6e1 - t266 * t608 * t8 * t202 / 0.2e1 - t1406 / 0.6e1 - t614 * t635 * t8 * t245 / 0.3e1 - 0.7e1 / 0.6e1 * t1414 - 0.5e1 / 0.27e2 * t1417 * t267 * t118 * t472 + t1423 / 0.9e1 + t599 * t608 * t8 * t187 / 0.3e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t123 * t492;
1393 t1438 = t614 * t619 * t615;
1394 t1446 = t7 * t657 * t37;
1395 t1452 = t7 * t681 * t63;
1396 t1455 = t7 * t548 * t144;
1397 t1458 = t7 * t485 * t275;
1398 t1459 = t1458 / 0.3e1;
1399 t1461 = t7 * t657 * t661;
1400 t1469 = t7 * t42 * t1331 * t299;
1401 t1471 = t1438 / 0.9e1 - t1129 / 0.6e1 - t1143 / 0.6e1 + t1152 + t1156 - 0.5e1 / 0.12e2 * t1158 + t1161 / 0.4e1 - t1164 / 0.8e1 + 0.5e1 / 0.2e1 * t1446 - 0.3e1 / 0.2e1 * t7 * t281 * t101 - 0.5e1 / 0.2e1 * t1452 - 0.5e1 / 0.36e2 * t1455 + t1459 + 0.7e1 / 0.2e1 * t1461 - 0.3e1 / 0.2e1 * t7 * t281 * t238 * t8 + 0.25e2 / 0.6e1 * t1469;
1402 t1473 = t7 * t250 * t685;
1403 t1476 = t7 * t192 * t677;
1404 t1496 = -t1473 / 0.2e1 - t1476 / 0.2e1 - 0.5e1 / 0.9e1 * t1205 + t1218 / 0.27e2 + t1222 / 0.9e1 - 0.5e1 / 0.9e1 * t1231 - t1238 / 0.6e1 + 0.2e1 / 0.3e1 * t1241 + t1244 / 0.3e1 + t1247 / 0.3e1 - t1253 / 0.2e1 - t1256 / 0.2e1 - t1259 / 0.6e1 - t1283 / 0.6e1 + t1301 / 0.27e2 + t1305 / 0.9e1 + t599 * t267 * t627 * t74 / 0.3e1;
1405 t1501 = t606 / 0.3e1;
1406 t1503 = t621 / 0.3e1;
1407 t1506 = t9 * (t1372 + t1435 + t1471 + t1496) + t602 / 0.3e1 - t1501 - t611 + t617 / 0.3e1 - t1503 - t625 - t629 / 0.2e1 - t633 / 0.2e1 + t638 - t484;
1408 t1512 = -t504 + 0.3e1 / 0.4e1 * t553 - 0.3e1 / 0.4e1 * t556 - t588 - 0.3e1 / 0.2e1 * t643 + 0.3e1 / 0.2e1 * t646 + t832 - t1118 - t1121 + t796 + 0.7e1 / 0.2e1 * t659;
1409 t1520 = -0.3e1 * t663 + t797 - t669 / 0.2e1 - t672 - t675 - 0.3e1 / 0.2e1 * t679 - 0.7e1 / 0.2e1 * t683 - 0.3e1 / 0.2e1 * t687 + t476 - t479 / 0.3e1 + t487;
1410 t1523 = t496 - t499 / 0.3e1 + t507 - t512 - 0.2e1 * t514 - t517 - t522 - t526 - t535 - t541 - t547;
1411 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1412 v4rho4[1] = t1506 + t1512 + t1520 + t1523;
1413
1414 t1537 = t7 * t681 * t764;
1415 t1540 = t7 * t250 * t768;
1416 t1547 = t7 * t250 * t772;
1417 t1553 = t7 * t505 * t331;
1418 t1556 = t7 * t505 * t335;
1419 t1559 = 0.1e1 / t79 / t1273;
1420 t1561 = t8 * t8;
1421 t1567 = t7 * t477 * t750;
1422 t1570 = t7 * t192 * t754;
1423 t1573 = t7 * t192 * t760;
1424 t1596 = -0.4e1 / 0.9e1 * t614 * t635 * t291 * t8 + 0.4e1 / 0.9e1 * t599 * t608 * t268 * t8 - t7 * t1319 * t119 / 0.3e1 - 0.7e1 / 0.3e1 * t1537 - t1540 / 0.3e1 - t7 * t244 * t362 * t245 / 0.6e1 - t1547 / 0.3e1 + t7 * t1336 * t124 / 0.3e1 + 0.2e1 / 0.9e1 * t1553 + t1556 / 0.9e1 - t7 * t185 * t1559 * t37 * t1561 / 0.3e1 - t1567 / 0.9e1 - t1570 / 0.3e1 - t1573 / 0.3e1 - t7 * t186 * t238 * t314 / 0.6e1 - t7 * t69 * t238 * t320 / 0.2e1 - t7 * t69 * t101 * t759 - t7 * t105 * t362 * t254 / 0.2e1 - t7 * t243 * t1559 * t63 * t1561 / 0.3e1 - 0.10e2 / 0.27e2 * t1328;
1425 t1602 = t7 * t1303 * t315;
1426 t1605 = t7 * t485 * t321;
1427 t1608 = t7 * t1220 * t326;
1428 t1615 = t599 * t604 * t728;
1429 t1618 = t266 * t1368 * t712;
1430 t1625 = t266 * t604 * t704;
1431 t1632 = t614 * t619 * t708;
1432 t1639 = t289 * t619 * t724;
1433 t1655 = -0.41e2 / 0.9e1 * t1334 - 0.10e2 / 0.27e2 * t1341 - t1348 / 0.9e1 + 0.4e1 / 0.9e1 * t1351 + t1602 / 0.27e2 + t1605 / 0.9e1 + t1608 / 0.27e2 - 0.5e1 / 0.27e2 * t1417 * t267 * t314 * t187 + 0.2e1 / 0.27e2 * t1615 + 0.7e1 / 0.9e1 * t1618 + t599 * t267 * t320 * t187 / 0.9e1 - t1625 / 0.9e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t325 * t245 + 0.2e1 / 0.27e2 * t1632 + t614 * t290 * t334 * t245 / 0.9e1 - t1639 / 0.9e1 + 0.2e1 / 0.9e1 * t599 * t623 * t728 + t599 * t267 * t314 * t202 / 0.9e1 - t266 * t623 * t704 / 0.3e1 - t266 * t267 * t759 * t74 / 0.3e1;
1434 t1666 = t289 * t1412 * t716;
1435 t1671 = t280 * t144;
1436 t1694 = t289 * t80 * t144 * t291;
1437 t1704 = -t266 * t267 * t320 * t202 / 0.6e1 + t614 * t290 * t325 * t254 / 0.9e1 - 0.7e1 / 0.9e1 * t1666 - 0.2e1 / 0.3e1 * t266 * t1395 * t712 + 0.2e1 / 0.3e1 * t289 * t1671 * t636 - t289 * t290 * t767 * t106 / 0.3e1 - t289 * t720 * t631 / 0.3e1 - t289 * t290 * t334 * t254 / 0.6e1 + 0.4e1 / 0.27e2 * t1366 - 0.2e1 / 0.9e1 * t1370 + 0.4e1 / 0.27e2 * t1383 + 0.2e1 / 0.9e1 * t614 * t720 * t615 - 0.2e1 / 0.9e1 * t1694 - 0.2e1 / 0.9e1 * t1393 - t1399 / 0.9e1 - t1406 / 0.9e1 + 0.2e1 / 0.9e1 * t1414 + 0.2e1 / 0.27e2 * t1423 + 0.2e1 / 0.27e2 * t1438 + t7 * t298 * t144;
1438 t1706 = t7 * t259 * t362;
1439 t1725 = t1706 / 0.12e2 + t1152 + t1156 - 0.5e1 / 0.18e2 * t1158 + t1161 / 0.12e2 - t1446 / 0.3e1 + t1452 / 0.3e1 - 0.5e1 / 0.18e2 * t1455 + 0.4e1 / 0.9e1 * t1458 - 0.2e1 / 0.3e1 * t1461 + 0.41e2 / 0.9e1 * t1469 - t1473 / 0.3e1 - t1476 / 0.3e1 - 0.10e2 / 0.27e2 * t1205 + t1222 / 0.27e2 - 0.10e2 / 0.27e2 * t1231 + 0.2e1 / 0.9e1 * t1241 + t1244 / 0.9e1 + t1247 / 0.9e1 + t1305 / 0.27e2;
1440 t1741 = t9 * (t1596 + t1655 + t1704 + t1725) - t706 / 0.3e1 + 0.2e1 / 0.9e1 * t710 - 0.2e1 / 0.3e1 * t714 + 0.2e1 / 0.3e1 * t718 - 0.2e1 / 0.3e1 * t722 - t726 / 0.3e1 + 0.2e1 / 0.9e1 * t602 - 0.4e1 / 0.9e1 * t606 - 0.2e1 / 0.3e1 * t611 + 0.2e1 / 0.9e1 * t617 - 0.4e1 / 0.9e1 * t621 + 0.2e1 / 0.9e1 * t730;
1441 t1751 = -0.2e1 / 0.3e1 * t625 - t629 / 0.3e1 - t633 / 0.3e1 + 0.2e1 / 0.3e1 * t638 + t553 / 0.2e1 - t556 / 0.4e1 - t643 + t646 + t649 / 0.2e1 - t735 / 0.4e1 - t1118 - t1121 - t738 / 0.3e1;
1442 t1762 = -t741 / 0.9e1 - t745 / 0.3e1 - t748 / 0.9e1 - t752 / 0.3e1 - t756 - t762 + 0.2e1 * t766 - t770 - t774 - 0.2e1 / 0.3e1 * t777 + 0.2e1 / 0.3e1 * t653 + 0.5e1 / 0.3e1 * t659 - 0.2e1 * t663;
1443 t1775 = 0.2e1 / 0.3e1 * t665 - t669 / 0.3e1 - 0.4e1 / 0.3e1 * t672 - 0.4e1 / 0.3e1 * t675 - t679 - 0.5e1 / 0.3e1 * t683 - t687 - t479 / 0.9e1 + 0.2e1 / 0.3e1 * t487 - t499 / 0.9e1 + 0.2e1 / 0.3e1 * t507 - 0.2e1 / 0.3e1 * t514 - t517 / 0.3e1 - t541 / 0.3e1;
1444 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1445 v4rho4[2] = t1741 + t1751 + t1762 + t1775;
1446
1447 t1790 = t7 * t192 * t837;
1448 t1793 = t7 * t497 * t841;
1449 t1796 = t7 * t250 * t845;
1450 t1799 = t7 * t250 * t849;
1451 t1802 = t7 * t250 * t854;
1452 t1805 = t7 * t1299 * t860;
1453 t1808 = t7 * t1216 * t869;
1454 t1816 = 0.6e1 * t528 + 0.12e2 * t1275;
1455 t1838 = t599 * t267 * t864 * t74 / 0.3e1 + t614 * t290 * t873 * t106 / 0.3e1 - t7 * t105 * t828 * t106 / 0.2e1 - t1790 / 0.6e1 - t1793 / 0.6e1 - t1796 / 0.2e1 - t1799 / 0.2e1 - t1802 / 0.6e1 + t1805 / 0.27e2 + t1808 / 0.27e2 - t7 * t69 * t101 * t836 / 0.2e1 - t7 * t69 * t37 * t1816 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t298 * t362 * t8 - 0.3e1 / 0.2e1 * t7 * t105 * t144 * t767 + t7 * t105 * t63 * t1816 / 0.2e1 + t7 * t471 * t101 * t859 / 0.9e1;
1456 t1858 = t1537 - t1540 / 0.2e1 - t1547 / 0.2e1 + 0.2e1 / 0.3e1 * t1553 + t1556 / 0.3e1 - t1567 / 0.6e1 - t1570 / 0.2e1 - t1573 / 0.2e1 - 0.5e1 / 0.9e1 * t1328 + t1334 / 0.3e1 - 0.5e1 / 0.9e1 * t1341 + t1352 + t1602 / 0.9e1 + t1605 / 0.3e1 + t1608 / 0.9e1 + t289 * t1671 * t716 - t289 * t16 * t362 * t291 / 0.2e1;
1457 t1899 = t289 * t619 * t873;
1458 t1906 = t266 * t604 * t864;
1459 t1917 = -t289 * t720 * t724 / 0.2e1 - t289 * t290 * t853 * t106 / 0.6e1 - 0.5e1 / 0.27e2 * t1417 * t267 * t859 * t74 - t266 * t623 * t864 / 0.2e1 - t266 * t267 * t118 * t759 / 0.2e1 - 0.5e1 / 0.27e2 * t1430 * t290 * t868 * t106 - t614 * t635 * t325 * t8 / 0.3e1 - t289 * t290 * t123 * t767 / 0.2e1 + t289 * t635 * t8 * t334 / 0.2e1 - t266 * t608 * t8 * t320 / 0.2e1 - t1899 / 0.6e1 + t599 * t608 * t314 * t8 / 0.3e1 - t1906 / 0.6e1 - t266 * t267 * t836 * t74 / 0.6e1 + t614 * t720 * t708 / 0.3e1 + t1615 / 0.9e1 - t1618 / 0.3e1;
1460 t1925 = t7 * t111 * t828;
1461 t1932 = -t1625 / 0.6e1 + t1632 / 0.9e1 - t1639 / 0.6e1 + t1666 / 0.3e1 + t1367 + t1384 - t1694 / 0.3e1 + t1706 / 0.4e1 - t1925 / 0.8e1 + t1152 + t1156 - 0.5e1 / 0.36e2 * t1158 - 0.5e1 / 0.12e2 * t1455 + t1459 - t1469 / 0.3e1 - 0.5e1 / 0.27e2 * t1205 - 0.5e1 / 0.27e2 * t1231;
1462 t1939 = t9 * (t1838 + t1858 + t1917 + t1932) - t867 - t876 - t706 / 0.2e1 + t710 / 0.3e1 - t714 + t718 - t722 - t726 / 0.2e1 - t1501 - t1503;
1463 t1943 = t730 / 0.3e1 + t554 + 0.3e1 / 0.4e1 * t649 - 0.3e1 / 0.4e1 * t735 - t831 - t1118 - t1121 - t738 - t840 - t844 - t848;
1464 t1951 = -t852 - t857 + t863 + t872 - t741 / 0.3e1 - t745 - t748 / 0.3e1 - t752 / 0.2e1 - 0.3e1 / 0.2e1 * t756 - 0.3e1 / 0.2e1 * t762 + 0.3e1 * t766;
1465 t1955 = -0.3e1 / 0.2e1 * t770 - 0.3e1 / 0.2e1 * t774 - 0.2e1 * t777 + t653 - t659 + t665 - t672 - t675 + t683 + t488 + t508;
1466 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1467 v4rho4[3] = t1939 + t1943 + t1951 + t1955;
1468
1469 t1986 = t320 * t320;
1470 t1991 = t325 * t325;
1471 t1996 = t334 * t334;
1472 t2001 = t314 * t314;
1473 t2007 = 0.12e2 * t528 + 0.12e2 * t1275;
1474 t2035 = 0.2e1 / 0.3e1 * t614 * t290 * t325 * t334 - 0.2e1 / 0.3e1 * t289 * t290 * t853 * t123 + 0.2e1 / 0.3e1 * t599 * t267 * t314 * t320 - 0.2e1 / 0.3e1 * t266 * t267 * t836 * t118 - 0.2e1 * t289 * t720 * t873 + 0.4e1 / 0.9e1 * t7 * t491 * t144 * t868 - t7 * t186 * t37 * t1986 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1210 * t63 * t1991 - t7 * t244 * t63 * t1996 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1293 * t37 * t2001 - t7 * t69 * t37 * t2007 / 0.2e1 - t7 * t244 * t362 * t325 - 0.3e1 * t7 * t105 * t362 * t334 - 0.2e1 * t7 * t105 * t144 * t853 + t7 * t105 * t63 * t2007 / 0.2e1 - 0.2e1 * t7 * t105 * t828 * t123 - 0.2e1 / 0.3e1 * t1790 - 0.2e1 / 0.3e1 * t1793 - 0.2e1 * t1796;
1475 t2050 = 0.1e1 / t48 / t355;
1476 t2063 = t820 * rho[1];
1477 t2074 = 0.1e1 / t823 / t59;
1478 t2085 = -0.2e1 * t1799 - 0.2e1 / 0.3e1 * t1802 + 0.4e1 / 0.27e2 * t1805 + 0.4e1 / 0.27e2 * t1808 + 0.4e1 / 0.3e1 * t1553 + 0.2e1 / 0.3e1 * t1556 - 0.20e2 / 0.27e2 * t1328 - 0.20e2 / 0.27e2 * t1341 + 0.2e1 / 0.9e1 * t1602 + 0.2e1 / 0.3e1 * t1605 + 0.2e1 / 0.9e1 * t1608 - 0.2e1 / 0.3e1 * t1899 - 0.2e1 / 0.3e1 * t1906 - 0.3e1 / 0.8e1 * t7 * t44 * (0.35912208504801097395e1 * t45 * t2050 + 0.3640e4 / 0.81e2 * t54 * t1094 * t60 - 0.23072e5 / 0.27e2 * t138 * t2050 * t140 + 0.476672e6 / 0.81e2 * t354 / t820 * t358 - 0.475136e6 / 0.27e2 * t819 / t47 / t2063 * t824 + 0.524288e6 / 0.27e2 * params->cc * t1084 / t48 / t820 / t46 * t2074) + t1706 / 0.2e1 - t1925 / 0.2e1 + t1152 + t1156 - 0.5e1 / 0.9e1 * t1455;
1479 t2089 = 0.4e1 / 0.3e1 * t653 + 0.4e1 / 0.3e1 * t665 - 0.2e1 / 0.3e1 * t741 - 0.2e1 * t745 - 0.2e1 / 0.3e1 * t748 - t1118 - t1121 - 0.3e1 / 0.2e1 * t830 + t9 * (t2035 + t2085) + t649 - 0.3e1 / 0.2e1 * t735;
1480 t2101 = -0.2e1 * t839 - 0.2e1 * t843 - 0.6e1 * t847 - 0.6e1 * t851 - 0.2e1 * t856 - 0.4e1 * t777 - 0.2e1 * t738 + 0.4e1 / 0.9e1 * t862 - 0.2e1 * t866 + 0.4e1 / 0.9e1 * t871 - 0.2e1 * t875;
1481 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1482 v4rho4[4] = t2089 + t2101;
1483
1484 t2103 = t7 * t186 * t888;
1485 t2106 = t7 * t192 * t372;
1486 t2124 = t7 * t69 * t892;
1487 t2127 = t7 * t69 * t895;
1488 t2131 = t7 * t209 * t163 / 0.9e1;
1489 t2141 = t7 * t81 * t389;
1490 t2144 = t7 * t17 * t915;
1491 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1492 v4rho3sigma[0] = -0.2e1 / 0.3e1 * t2103 - 0.2e1 / 0.3e1 * t2106 + t152 * t598 * t163 * t472 / 0.9e1 - t940 * t941 * t481 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t371 * t389 * t202 - t152 * t371 * t163 * t531 / 0.2e1 - 0.4e1 * t2124 - 0.2e1 * t2127 + t2131 - t152 * t265 * t389 * t187 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t371 * t915 * t74 - t2141 / 0.2e1 - 0.3e1 / 0.2e1 * t2144 - 0.3e1 / 0.8e1 * t152 * t153 * (-0.63374485596707818932e0 * params->bb * t559 - 0.140e3 / 0.27e2 * t157 * t564 + 0.1456e4 / 0.9e1 * params->cc * t559 * t97 - 0.40576e5 / 0.27e2 * params->cc * t571 * t386 + 0.50176e5 / 0.9e1 * params->cc * t579 * t912 - 0.65536e5 / 0.9e1 * params->cc / t21 / t1180 * t1191 * t229);
1493
1494 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1495 v4rho3sigma[1] = 0.0e0;
1496
1497 t2170 = t7 * t244 * t925;
1498 t2173 = t7 * t250 * t396;
1499 t2176 = t7 * t105 * t929;
1500 t2180 = t7 * t259 * t177 / 0.9e1;
1501 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1502 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2170 - 0.2e1 / 0.3e1 * t2173 - 0.2e1 * t2176 + t2180 + t152 * t613 * t177 * t492 / 0.9e1 - t940 * t958 * t501 / 0.2e1 - t152 * t395 * t177 * t543 / 0.2e1;
1503
1504 t2193 = t280 * t1 * t3;
1505 t2194 = t2193 * t952;
1506 t2207 = t80 * t1 * t3;
1507 t2208 = t2207 * t401;
1508 t2210 = t935 * t939;
1509 t2212 = 0.4e1 / 0.9e1 * t2210 * t942;
1510 t2223 = 0.4e1 / 0.3e1 * t936 * t946;
1511 t2226 = t470 * t163;
1512 t2230 = t949 * t939;
1513 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1514 v4rho3sigma[3] = t2194 / 0.6e1 - t940 * t185 * t389 * t268 / 0.3e1 - t940 * t941 * t627 / 0.6e1 - t950 * t371 * t389 * t8 - 0.2e1 / 0.9e1 * t2208 - t2212 - 0.2e1 / 0.9e1 * t2103 - 0.4e1 / 0.9e1 * t2106 - t152 * t371 * t915 * t118 / 0.2e1 - t950 * t371 * t163 / 0.2e1 - t2223 - 0.4e1 / 0.3e1 * t2124 - 0.2e1 / 0.3e1 * t2127 + t2131 + t940 * t2226 * t600 / 0.9e1 - t2230 * t941 * t609 / 0.3e1 - t2141 / 0.3e1 - t2144 / 0.2e1;
1515
1516 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1517 v4rho3sigma[4] = 0.0e0;
1518
1519 t2239 = t2207 * t406;
1520 t2242 = 0.4e1 / 0.9e1 * t2210 * t959;
1521 t2243 = t2193 * t963;
1522 t2245 = t490 * t177;
1523 t2258 = t2207 * t423;
1524 t2261 = 0.4e1 / 0.3e1 * t936 * t969;
1525 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1526 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2170 - 0.4e1 / 0.9e1 * t2173 - 0.2e1 / 0.3e1 * t2176 + t2180 - 0.2e1 / 0.9e1 * t2239 - t2242 - t2243 / 0.6e1 + t940 * t2245 * t615 / 0.9e1 + t2230 * t958 * t636 / 0.3e1 - t940 * t958 * t631 / 0.6e1 + t950 * t395 * t177 / 0.2e1 - t2258 / 0.6e1 - t2261 - t152 * t288 * t422 * t245 / 0.6e1 - t152 * t395 * t422 * t254 / 0.2e1;
1527
1528 t2274 = t936 * t974;
1529 t2286 = t936 * t978;
1530 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1531 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2208 - t2212 - t2223 - 0.4e1 / 0.3e1 * t2194 - 0.2e1 / 0.9e1 * t2106 + t2131 - t2141 / 0.6e1 - 0.2e1 / 0.9e1 * t2274 + t940 * t2226 * t728 / 0.9e1 - t152 * t265 * t389 * t314 / 0.6e1 - t2230 * t941 * t712 / 0.3e1 - 0.2e1 / 0.3e1 * t2286 - t940 * t941 * t704 / 0.6e1 - t152 * t371 * t389 * t320 / 0.2e1 - t152 * t371 * t163 * t759 / 0.2e1;
1532
1533 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1534 v4rho3sigma[7] = 0.0e0;
1535
1536 t2303 = t936 * t983;
1537 t2311 = t936 * t987;
1538 t2320 = t936 * t990;
1539 t2329 = t936 * t1010;
1540 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1541 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2239 - t2242 + 0.4e1 / 0.3e1 * t2243 - 0.2e1 / 0.9e1 * t2173 + t2180 - t2258 / 0.3e1 - t2261 - 0.2e1 / 0.9e1 * t2303 + t940 * t2245 * t708 / 0.9e1 + t2230 * t958 * t716 / 0.3e1 - 0.4e1 / 0.3e1 * t2311 - t940 * t243 * t422 * t291 / 0.3e1 + t950 * t395 * t422 * t8 - 0.2e1 / 0.3e1 * t2320 - t940 * t958 * t724 / 0.6e1 - t152 * t395 * t177 * t767 / 0.2e1 - t2329 / 0.2e1 - t152 * t395 * t1009 * t106 / 0.2e1;
1542
1543 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1544 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2208 - 0.2e1 / 0.3e1 * t2274 - 0.2e1 * t2286 + t2131 + t152 * t598 * t163 * t859 / 0.9e1 - t940 * t941 * t864 / 0.2e1 - t152 * t371 * t163 * t836 / 0.2e1;
1545
1546 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1547 v4rho3sigma[10] = 0.0e0;
1548
1549 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1550 v4rho3sigma[11] = -t152 * t288 * t422 * t325 / 0.2e1 - 0.3e1 / 0.2e1 * t152 * t395 * t1009 * t123 - 0.3e1 / 0.2e1 * t152 * t395 * t422 * t334 - t152 * t395 * t177 * t853 / 0.2e1 - t940 * t958 * t873 / 0.2e1 + t152 * t613 * t177 * t868 / 0.9e1 - 0.2e1 / 0.3e1 * t2303 - 0.4e1 * t2311 - 0.2e1 * t2320 - 0.2e1 / 0.3e1 * t2239 - 0.3e1 / 0.2e1 * t2329 - t2258 / 0.2e1 + t2180 - 0.3e1 / 0.8e1 * t152 * t167 * (-0.63374485596707818932e0 * params->bb * t802 - 0.140e3 / 0.27e2 * t171 * t807 + 0.1456e4 / 0.9e1 * params->cc * t802 * t140 - 0.40576e5 / 0.27e2 * params->cc * t814 * t419 + 0.50176e5 / 0.9e1 * params->cc * t822 * t1006 - 0.65536e5 / 0.9e1 * params->cc / t48 / t2063 * t2074 * t353);
1551
1552 t2402 = t2207 * t440 / 0.6e1;
1553 t2403 = t936 * t1016;
1554 t2405 = t936 * t1032;
1555 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556 v4rho2sigma2[0] = -t2402 - 0.4e1 / 0.3e1 * t2403 - t2405 - t152 * t265 * t439 * t187 / 0.6e1 - t152 * t371 * t1031 * t74 - t152 * t371 * t439 * t202 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * (-0.7e1 / 0.9e1 * t427 * t223 - 0.28e2 / 0.9e1 * t431 * t226 + 0.2048e4 / 0.9e1 * t157 * t235 - 0.4480e4 / 0.3e1 * t911 * t581 + 0.8192e4 / 0.3e1 * params->cc / t21 / t577 * t1191 * t26);
1557
1558 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559 v4rho2sigma2[1] = 0.0e0;
1560
1561 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1562 v4rho2sigma2[2] = 0.0e0;
1563
1564 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1565 v4rho2sigma2[3] = 0.0e0;
1566
1567 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1568 v4rho2sigma2[4] = 0.0e0;
1569
1570 t2436 = t2207 * t457 / 0.6e1;
1571 t2437 = t936 * t1037;
1572 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1573 v4rho2sigma2[5] = -t2436 - 0.4e1 / 0.3e1 * t2437 - t152 * t288 * t456 * t245 / 0.6e1 - t152 * t395 * t456 * t254 / 0.2e1;
1574
1575 t2449 = t936 * t1041;
1576 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1577 v4rho2sigma2[6] = -t2402 - 0.2e1 / 0.3e1 * t2403 - t2405 / 0.2e1 - 0.2e1 / 0.3e1 * t2449 - t940 * t185 * t439 * t268 / 0.6e1 - t152 * t371 * t1031 * t118 / 0.2e1 - t950 * t371 * t439 * t8 / 0.2e1;
1578
1579 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1580 v4rho2sigma2[7] = 0.0e0;
1581
1582 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1583 v4rho2sigma2[8] = 0.0e0;
1584
1585 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1586 v4rho2sigma2[9] = 0.0e0;
1587
1588 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1589 v4rho2sigma2[10] = 0.0e0;
1590
1591 t2464 = t936 * t1046;
1592 t2474 = t936 * t1062;
1593 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594 v4rho2sigma2[11] = -t2436 - 0.2e1 / 0.3e1 * t2437 - 0.2e1 / 0.3e1 * t2464 - t940 * t243 * t456 * t291 / 0.6e1 + t950 * t395 * t456 * t8 / 0.2e1 - t2474 / 0.2e1 - t152 * t395 * t1061 * t106 / 0.2e1;
1595
1596 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597 v4rho2sigma2[12] = -t2402 - 0.4e1 / 0.3e1 * t2449 - t152 * t265 * t439 * t314 / 0.6e1 - t152 * t371 * t439 * t320 / 0.2e1;
1598
1599 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1600 v4rho2sigma2[13] = 0.0e0;
1601
1602 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1603 v4rho2sigma2[14] = 0.0e0;
1604
1605 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1606 v4rho2sigma2[15] = 0.0e0;
1607
1608 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1609 v4rho2sigma2[16] = 0.0e0;
1610
1611 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1612 v4rho2sigma2[17] = -t2436 - 0.4e1 / 0.3e1 * t2464 - t2474 - t152 * t288 * t456 * t325 / 0.6e1 - t152 * t395 * t1061 * t123 - t152 * t395 * t456 * t334 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t167 * (-0.7e1 / 0.9e1 * t444 * t347 - 0.28e2 / 0.9e1 * t448 * t350 + 0.2048e4 / 0.9e1 * t171 * t359 - 0.4480e4 / 0.3e1 * t1005 * t824 + 0.8192e4 / 0.3e1 * params->cc / t48 / t820 * t2074 * t53);
1613
1614 t2519 = t936 * t1081;
1615 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1616 v4rhosigma3[0] = -t2519 / 0.2e1 - t152 * t371 * t1080 * t74 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t153 * (-t1067 * t92 / 0.2e1 - 0.2e1 * t1071 * t98 + 0.16e2 * t427 * t1023 + 0.256e3 * t1028 * t1077 - 0.1024e4 * params->cc / t21 / t570 * t1191 * t156);
1617
1618 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1619 v4rhosigma3[1] = 0.0e0;
1620
1621 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1622 v4rhosigma3[2] = 0.0e0;
1623
1624 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1625 v4rhosigma3[3] = 0.0e0;
1626
1627 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1628 v4rhosigma3[4] = 0.0e0;
1629
1630 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1631 v4rhosigma3[5] = 0.0e0;
1632
1633 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1634 v4rhosigma3[6] = 0.0e0;
1635
1636 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1637 v4rhosigma3[7] = 0.0e0;
1638
1639 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1640 v4rhosigma3[8] = 0.0e0;
1641
1642 t2543 = t936 * t1100;
1643 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1644 v4rhosigma3[9] = -t152 * t395 * t1099 * t106 / 0.2e1 - t2543 / 0.2e1;
1645
1646 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1647 v4rhosigma3[10] = -t152 * t371 * t1080 * t118 / 0.2e1 - t2519 / 0.2e1;
1648
1649 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1650 v4rhosigma3[11] = 0.0e0;
1651
1652 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1653 v4rhosigma3[12] = 0.0e0;
1654
1655 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1656 v4rhosigma3[13] = 0.0e0;
1657
1658 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1659 v4rhosigma3[14] = 0.0e0;
1660
1661 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1662 v4rhosigma3[15] = 0.0e0;
1663
1664 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1665 v4rhosigma3[16] = 0.0e0;
1666
1667 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1668 v4rhosigma3[17] = 0.0e0;
1669
1670 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1671 v4rhosigma3[18] = 0.0e0;
1672
1673 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1674 v4rhosigma3[19] = -t2543 / 0.2e1 - t152 * t395 * t1099 * t123 / 0.2e1 - 0.3e1 / 0.8e1 * t152 * t167 * (-t1086 * t135 / 0.2e1 - 0.2e1 * t1090 * t141 + 0.16e2 * t444 * t1053 + 0.256e3 * t1058 * t1096 - 0.1024e4 * params->cc / t48 / t813 * t2074 * t170);
1675
1676 t2575 = t575 * sigma[0];
1677 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1678 v4sigma4[0] = -0.3e1 / 0.8e1 * t152 * t153 * (-0.15e2 / 0.16e2 * params->cc / t26 / t2575 * t35 - 0.15e2 / 0.4e1 * params->cc / t2575 * t432 - 0.6e1 * t1067 * t434 * t234 + 0.384e3 * params->cc * t1167 * t1191 * t426 + 0.48e2 * t1076 * t581 * t1070);
1679
1680 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1681 v4sigma4[1] = 0.0e0;
1682
1683 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1684 v4sigma4[2] = 0.0e0;
1685
1686 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1687 v4sigma4[3] = 0.0e0;
1688
1689 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1690 v4sigma4[4] = 0.0e0;
1691
1692 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1693 v4sigma4[5] = 0.0e0;
1694
1695 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1696 v4sigma4[6] = 0.0e0;
1697
1698 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1699 v4sigma4[7] = 0.0e0;
1700
1701 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1702 v4sigma4[8] = 0.0e0;
1703
1704 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1705 v4sigma4[9] = 0.0e0;
1706
1707 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1708 v4sigma4[10] = 0.0e0;
1709
1710 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1711 v4sigma4[11] = 0.0e0;
1712
1713 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1714 v4sigma4[12] = 0.0e0;
1715
1716 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1717 v4sigma4[13] = 0.0e0;
1718
1719 t2599 = t818 * sigma[2];
1720 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1721 v4sigma4[14] = -0.3e1 / 0.8e1 * t152 * t167 * (-0.15e2 / 0.16e2 * params->cc / t53 / t2599 * t61 - 0.15e2 / 0.4e1 * params->cc / t2599 * t449 - 0.6e1 * t1086 * t451 * t358 + 0.384e3 * params->cc * t2050 * t2074 * t443 + 0.48e2 * t1095 * t824 * t1089);
1722
1723 #ifndef XC_DONT_COMPILE_MXC
1724
1725 if(order < 5) return;
1726
1727
1728 #endif
1729
1730 #endif
1731
1732 #endif
1733
1734 #endif
1735
1736 #endif
1737
1738
1739 }
1740
1741