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_k_pw86.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t2, t3, t4, t5, t6, t7, t8;
24 double t9, t12, t13, t14, t15, t17, t23, t24;
25 double t25, t26, t27, t29, t33, t34, t35, t38;
26 double t39, t41;
27
28 #ifndef XC_DONT_COMPILE_VXC
29 double t44, t45, t46, t47, t49, t50, t51, t52;
30 double t53, t55, t59, t61, t65, t66, t69, t76;
31 double t82, t83;
32
33 #ifndef XC_DONT_COMPILE_FXC
34 double t85, t90, t91, t96, t97, t98, t103, t109;
35 double t114, t117, t122, t126, t137, t138, t141, t142;
36 double t150, t151;
37
38 #ifndef XC_DONT_COMPILE_KXC
39 double t155, t160, t164, t171, t173, t174, t175, t179;
40 double t184, t189, t194, t197, t202, t210, t213, t214;
41 double t218, t229, t230, t236, t240, t245, t246, t254;
42 double t255, t258, t259, t263, t266, t267, t268;
43
44 #ifndef XC_DONT_COMPILE_LXC
45 double t301, t302, t310, t342, t343, t406, t467, t471;
46 double t479;
47 #endif
48
49 #endif
50
51 #endif
52
53 #endif
54
55 #endif
56
57
58
59 t1 = M_CBRT6;
60 t2 = t1 * t1;
61 t3 = M_PI * M_PI;
62 t4 = POW_1_3(t3);
63 t5 = t4 * t4;
64 t6 = t2 * t5;
65 t7 = M_CBRT2;
66 t8 = POW_1_3(rho[0]);
67 t9 = t8 * t8;
68 t12 = t1 / t5;
69 t13 = t7 * t7;
70 t14 = sigma[0] * t13;
71 t15 = rho[0] * rho[0];
72 t17 = 0.1e1 / t9 / t15;
73 t23 = t2 / t4 / t3;
74 t24 = sigma[0] * sigma[0];
75 t25 = t24 * t7;
76 t26 = t15 * t15;
77 t27 = t26 * rho[0];
78 t29 = 0.1e1 / t8 / t27;
79 t33 = t24 * sigma[0];
80 t34 = t26 * t26;
81 t35 = 0.1e1 / t34;
82 t38 = 0.1e1 + 0.91999999999999999998e-1 * t12 * t14 * t17 + 0.32187500000000000000e-1 * t23 * t25 * t29 + 0.35645771717653941627e-5 * t33 * t35;
83 t39 = pow(t38, 0.1e1 / 0.15e2);
84 t41 = t6 * t7 * t9 * t39;
85 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
86 zk[0] = 0.3e1 / 0.20e2 * t41;
87
88 #ifndef XC_DONT_COMPILE_VXC
89
90 if(order < 1) return;
91
92
93 t44 = t9 * rho[0] * t2;
94 t45 = t44 * t5;
95 t46 = t39 * t39;
96 t47 = t46 * t46;
97 t49 = t47 * t47;
98 t50 = t49 * t47 * t46;
99 t51 = 0.1e1 / t50;
100 t52 = t7 * t51;
101 t53 = t15 * rho[0];
102 t55 = 0.1e1 / t9 / t53;
103 t59 = t26 * t15;
104 t61 = 0.1e1 / t8 / t59;
105 t65 = t34 * rho[0];
106 t66 = 0.1e1 / t65;
107 t69 = -0.24533333333333333333e0 * t12 * t14 * t55 - 0.17166666666666666667e0 * t23 * t25 * t61 - 0.28516617374123153302e-4 * t33 * t66;
108 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
109 vrho[0] = t41 / 0.4e1 + t45 * t52 * t69 / 0.100e3;
110
111 t76 = sigma[0] * t7;
112 t82 = 0.91999999999999999998e-1 * t12 * t13 * t17 + 0.64375000000000000000e-1 * t23 * t76 * t29 + 0.10693731515296182488e-4 * t24 * t35;
113 t83 = t52 * t82;
114 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
115 vsigma[0] = t45 * t83 / 0.100e3;
116
117 #ifndef XC_DONT_COMPILE_FXC
118
119 if(order < 2) return;
120
121
122 t85 = 0.1e1 / t8;
123 t90 = t6 * t7;
124 t91 = t9 * t51;
125 t96 = 0.1e1 / t50 / t38;
126 t97 = t7 * t96;
127 t98 = t69 * t69;
128 t103 = 0.1e1 / t9 / t26;
129 t109 = 0.1e1 / t8 / t26 / t53;
130 t114 = 0.1e1 / t34 / t15;
131 t117 = 0.89955555555555555554e0 * t12 * t14 * t103 + 0.10872222222222222222e1 * t23 * t25 * t109 + 0.25664955636710837972e-3 * t33 * t114;
132 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
133 v2rho2[0] = t6 * t7 * t85 * t39 / 0.6e1 + t90 * t91 * t69 / 0.30e2 - 0.7e1 / 0.750e3 * t45 * t97 * t98 + t45 * t52 * t117 / 0.100e3;
134
135 t122 = t9 * t2 * t5;
136 t126 = t97 * t82 * t69;
137 t137 = -0.24533333333333333333e0 * t12 * t13 * t55 - 0.34333333333333333333e0 * t23 * t76 * t61 - 0.85549852122369459904e-4 * t24 * t66;
138 t138 = t52 * t137;
139 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
140 v2rhosigma[0] = t122 * t83 / 0.60e2 - 0.7e1 / 0.750e3 * t45 * t126 + t45 * t138 / 0.100e3;
141
142 t141 = t82 * t82;
143 t142 = t97 * t141;
144 t150 = 0.64375000000000000000e-1 * t23 * t7 * t29 + 0.21387463030592364976e-4 * sigma[0] * t35;
145 t151 = t52 * t150;
146 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
147 v2sigma2[0] = -0.7e1 / 0.750e3 * t45 * t142 + t45 * t151 / 0.100e3;
148
149 #ifndef XC_DONT_COMPILE_KXC
150
151 if(order < 3) return;
152
153
154 t155 = 0.1e1 / t8 / rho[0];
155 t160 = t85 * t51;
156 t164 = t9 * t96;
157 t171 = t38 * t38;
158 t173 = 0.1e1 / t50 / t171;
159 t174 = t7 * t173;
160 t175 = t98 * t69;
161 t179 = t69 * t117;
162 t184 = 0.1e1 / t9 / t27;
163 t189 = 0.1e1 / t8 / t34;
164 t194 = 0.1e1 / t34 / t53;
165 t197 = -0.41979259259259259259e1 * t12 * t14 * t184 - 0.79729629629629629628e1 * t23 * t25 * t189 - 0.25664955636710837972e-2 * t33 * t194;
166 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
167 v3rho3[0] = -t6 * t7 * t155 * t39 / 0.18e2 + t90 * t160 * t69 / 0.30e2 - 0.7e1 / 0.150e3 * t90 * t164 * t98 + t90 * t91 * t117 / 0.20e2 + 0.203e3 / 0.11250e5 * t45 * t174 * t175 - 0.7e1 / 0.250e3 * t45 * t97 * t179 + t45 * t52 * t197 / 0.100e3;
168
169 t202 = t85 * t2 * t5;
170 t210 = t174 * t82 * t98;
171 t213 = t137 * t69;
172 t214 = t97 * t213;
173 t218 = t97 * t82 * t117;
174 t229 = 0.89955555555555555554e0 * t12 * t13 * t103 + 0.21744444444444444444e1 * t23 * t76 * t109 + 0.76994866910132513914e-3 * t24 * t114;
175 t230 = t52 * t229;
176 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
177 v3rho2sigma[0] = t202 * t83 / 0.90e2 - 0.7e1 / 0.225e3 * t122 * t126 + t122 * t138 / 0.30e2 + 0.203e3 / 0.11250e5 * t45 * t210 - 0.7e1 / 0.375e3 * t45 * t214 - 0.7e1 / 0.750e3 * t45 * t218 + t45 * t230 / 0.100e3;
178
179 t236 = t174 * t141 * t69;
180 t240 = t97 * t82 * t137;
181 t245 = t150 * t69;
182 t246 = t97 * t245;
183 t254 = -0.34333333333333333333e0 * t23 * t7 * t61 - 0.17109970424473891981e-3 * sigma[0] * t66;
184 t255 = t52 * t254;
185 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
186 v3rhosigma2[0] = -0.7e1 / 0.450e3 * t122 * t142 + 0.203e3 / 0.11250e5 * t45 * t236 - 0.7e1 / 0.375e3 * t45 * t240 + t122 * t151 / 0.60e2 - 0.7e1 / 0.750e3 * t45 * t246 + t45 * t255 / 0.100e3;
187
188 t258 = t141 * t82;
189 t259 = t174 * t258;
190 t263 = t97 * t82 * t150;
191 t266 = t61 * t2;
192 t267 = t5 * t7;
193 t268 = t267 * t51;
194 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
195 v3sigma3[0] = 0.203e3 / 0.11250e5 * t45 * t259 - 0.7e1 / 0.250e3 * t45 * t263 + 0.21387463030592364976e-6 * t266 * t268;
196
197 #ifndef XC_DONT_COMPILE_LXC
198
199 if(order < 4) return;
200
201
202 t301 = t7 / t50 / t171 / t38;
203 t302 = t98 * t98;
204 t310 = t117 * t117;
205 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
206 v4rho4[0] = 0.2e1 / 0.27e2 * t6 * t7 / t8 / t15 * t39 - 0.2e1 / 0.135e3 * t90 * t155 * t51 * t69 - 0.14e2 / 0.225e3 * t90 * t85 * t96 * t98 + t90 * t160 * t117 / 0.15e2 + 0.406e3 / 0.3375e4 * t90 * t9 * t173 * t175 - 0.14e2 / 0.75e2 * t90 * t164 * t179 + t90 * t91 * t197 / 0.15e2 - 0.4466e4 / 0.84375e5 * t45 * t301 * t302 + 0.203e3 / 0.1875e4 * t45 * t174 * t98 * t117 - 0.7e1 / 0.250e3 * t45 * t97 * t310 - 0.14e2 / 0.375e3 * t45 * t97 * t69 * t197 + t45 * t52 * (0.23788246913580246913e2 * t12 * t14 / t9 / t59 + 0.66441358024691358023e2 * t23 * t25 / t8 / t65 + 0.28231451200381921769e-1 * t33 / t34 / t26) / 0.100e3;
207
208 t342 = t44 * t267;
209 t343 = t173 * t82;
210 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
211 v4rho3sigma[0] = 0.203e3 / 0.2250e4 * t122 * t210 - 0.4466e4 / 0.84375e5 * t45 * t301 * t82 * t175 + 0.203e3 / 0.3750e4 * t342 * t343 * t179 - t155 * t2 * t5 * t83 / 0.270e3 - 0.7e1 / 0.225e3 * t202 * t126 - 0.7e1 / 0.75e2 * t122 * t214 - 0.7e1 / 0.150e3 * t122 * t218 + 0.203e3 / 0.3750e4 * t45 * t174 * t137 * t98 - 0.7e1 / 0.250e3 * t45 * t97 * t229 * t69 - 0.7e1 / 0.250e3 * t45 * t97 * t137 * t117 - 0.7e1 / 0.750e3 * t45 * t97 * t82 * t197 + t202 * t138 / 0.30e2 + t122 * t230 / 0.20e2 + t45 * t52 * (-0.41979259259259259259e1 * t12 * t13 * t184 - 0.15945925925925925926e2 * t23 * t76 * t189 - 0.76994866910132513914e-2 * t24 * t194) / 0.100e3;
212
213 t406 = t137 * t137;
214 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
215 v4rho2sigma2[0] = -0.7e1 / 0.675e3 * t202 * t142 + 0.203e3 / 0.3375e4 * t122 * t236 - 0.14e2 / 0.225e3 * t122 * t240 - 0.4466e4 / 0.84375e5 * t45 * t301 * t141 * t98 + 0.406e3 / 0.5625e4 * t342 * t343 * t213 + 0.203e3 / 0.11250e5 * t45 * t174 * t141 * t117 - 0.7e1 / 0.375e3 * t45 * t97 * t406 - 0.7e1 / 0.375e3 * t45 * t97 * t82 * t229 + t202 * t151 / 0.90e2 - 0.7e1 / 0.225e3 * t122 * t246 + t122 * t255 / 0.30e2 + 0.203e3 / 0.11250e5 * t45 * t174 * t150 * t98 - 0.7e1 / 0.375e3 * t45 * t97 * t254 * t69 - 0.7e1 / 0.750e3 * t45 * t97 * t150 * t117 + t45 * t52 * (0.21744444444444444444e1 * t23 * t7 * t109 + 0.15398973382026502783e-2 * sigma[0] * t114) / 0.100e3;
216
217 t467 = t266 * t5;
218 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
219 v4rhosigma3[0] = 0.203e3 / 0.6750e4 * t122 * t259 - 0.4466e4 / 0.84375e5 * t45 * t301 * t258 * t69 + 0.203e3 / 0.3750e4 * t45 * t174 * t141 * t137 - 0.7e1 / 0.150e3 * t122 * t263 + 0.203e3 / 0.3750e4 * t342 * t343 * t245 - 0.7e1 / 0.250e3 * t45 * t97 * t137 * t150 - 0.7e1 / 0.250e3 * t45 * t97 * t82 * t254 - 0.13545393252708497818e-5 * t109 * t2 * t268 - 0.19961632161886207311e-6 * t467 * t97 * t69;
220
221 t471 = t141 * t141;
222 t479 = t150 * t150;
223 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
224 v4sigma4[0] = -0.4466e4 / 0.84375e5 * t45 * t301 * t471 + 0.203e3 / 0.1875e4 * t45 * t174 * t141 * t150 - 0.7e1 / 0.250e3 * t45 * t97 * t479 - 0.79846528647544829244e-6 * t467 * t97 * t82;
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, t2, t3, t4, t5, t6, t7, t8;
251 double t10, t11, t13, t19, t20, t21, t22, t24;
252 double t28, t29, t30, t33, t34, t36;
253
254 #ifndef XC_DONT_COMPILE_VXC
255 double t39, t40, t41, t43, t44, t45, t46, t47;
256 double t49, t53, t55, t59, t60, t63, t74, t75;
257
258 #ifndef XC_DONT_COMPILE_FXC
259 double t77, t81, t86, t87, t88, t93, t99, t104;
260 double t107, t111, t114, t115, t116, t126, t127, t130;
261 double t131, t138, t139;
262
263 #ifndef XC_DONT_COMPILE_KXC
264 double t143, t147, t158, t160, t161, t162, t166, t167;
265 double t171, t176, t181, t184, t188, t191, t196, t197;
266 double t200, t201, t204, t214, t215, t220, t221, t224;
267 double t229, t230, t237, t238, t241, t242, t245, t248;
268
269 #ifndef XC_DONT_COMPILE_LXC
270 double t278, t279, t280, t288, t379, t438, t445;
271 #endif
272
273 #endif
274
275 #endif
276
277 #endif
278
279 #endif
280
281
282
283 t1 = M_CBRT6;
284 t2 = t1 * t1;
285 t3 = M_PI * M_PI;
286 t4 = POW_1_3(t3);
287 t5 = t4 * t4;
288 t6 = t2 * t5;
289 t7 = POW_1_3(rho[0]);
290 t8 = t7 * t7;
291 t10 = t1 / t5;
292 t11 = rho[0] * rho[0];
293 t13 = 0.1e1 / t8 / t11;
294 t19 = t2 / t4 / t3;
295 t20 = sigma[0] * sigma[0];
296 t21 = t11 * t11;
297 t22 = t21 * rho[0];
298 t24 = 0.1e1 / t7 / t22;
299 t28 = t20 * sigma[0];
300 t29 = t21 * t21;
301 t30 = 0.1e1 / t29;
302 t33 = 0.1e1 + 0.91999999999999999998e-1 * t10 * sigma[0] * t13 + 0.16093750000000000000e-1 * t19 * t20 * t24 + 0.89114429294134854068e-6 * t28 * t30;
303 t34 = pow(t33, 0.1e1 / 0.15e2);
304 t36 = t6 * t8 * t34;
305 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
306 zk[0] = 0.3e1 / 0.10e2 * t36;
307
308 #ifndef XC_DONT_COMPILE_VXC
309
310 if(order < 1) return;
311
312
313 t39 = t8 * rho[0] * t2;
314 t40 = t34 * t34;
315 t41 = t40 * t40;
316 t43 = t41 * t41;
317 t44 = t43 * t41 * t40;
318 t45 = 0.1e1 / t44;
319 t46 = t5 * t45;
320 t47 = t11 * rho[0];
321 t49 = 0.1e1 / t8 / t47;
322 t53 = t21 * t11;
323 t55 = 0.1e1 / t7 / t53;
324 t59 = t29 * rho[0];
325 t60 = 0.1e1 / t59;
326 t63 = -0.24533333333333333333e0 * t10 * sigma[0] * t49 - 0.85833333333333333333e-1 * t19 * t20 * t55 - 0.71291543435307883254e-5 * t28 * t60;
327 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
328 vrho[0] = t36 / 0.2e1 + t39 * t46 * t63 / 0.50e2;
329
330 t74 = 0.91999999999999999998e-1 * t10 * t13 + 0.32187500000000000000e-1 * t19 * sigma[0] * t24 + 0.26734328788240456220e-5 * t20 * t30;
331 t75 = t46 * t74;
332 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
333 vsigma[0] = t39 * t75 / 0.50e2;
334
335 #ifndef XC_DONT_COMPILE_FXC
336
337 if(order < 2) return;
338
339
340 t77 = 0.1e1 / t7;
341 t81 = t8 * t45;
342 t86 = 0.1e1 / t44 / t33;
343 t87 = t5 * t86;
344 t88 = t63 * t63;
345 t93 = 0.1e1 / t8 / t21;
346 t99 = 0.1e1 / t7 / t21 / t47;
347 t104 = 0.1e1 / t29 / t11;
348 t107 = 0.89955555555555555554e0 * t10 * sigma[0] * t93 + 0.54361111111111111111e0 * t19 * t20 * t99 + 0.64162389091777094929e-4 * t28 * t104;
349 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
350 v2rho2[0] = t6 * t77 * t34 / 0.3e1 + t6 * t81 * t63 / 0.15e2 - 0.7e1 / 0.375e3 * t39 * t87 * t88 + t39 * t46 * t107 / 0.50e2;
351
352 t111 = t8 * t2;
353 t114 = t39 * t5;
354 t115 = t86 * t74;
355 t116 = t115 * t63;
356 t126 = -0.24533333333333333333e0 * t10 * t49 - 0.17166666666666666667e0 * t19 * sigma[0] * t55 - 0.21387463030592364976e-4 * t20 * t60;
357 t127 = t46 * t126;
358 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
359 v2rhosigma[0] = t111 * t75 / 0.30e2 - 0.7e1 / 0.375e3 * t114 * t116 + t39 * t127 / 0.50e2;
360
361 t130 = t74 * t74;
362 t131 = t87 * t130;
363 t138 = 0.32187500000000000000e-1 * t19 * t24 + 0.53468657576480912440e-5 * sigma[0] * t30;
364 t139 = t46 * t138;
365 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
366 v2sigma2[0] = -0.7e1 / 0.375e3 * t39 * t131 + t39 * t139 / 0.50e2;
367
368 #ifndef XC_DONT_COMPILE_KXC
369
370 if(order < 3) return;
371
372
373 t143 = 0.1e1 / t7 / rho[0];
374 t147 = t77 * t45;
375 t158 = t33 * t33;
376 t160 = 0.1e1 / t44 / t158;
377 t161 = t5 * t160;
378 t162 = t88 * t63;
379 t166 = t86 * t63;
380 t167 = t166 * t107;
381 t171 = 0.1e1 / t8 / t22;
382 t176 = 0.1e1 / t7 / t29;
383 t181 = 0.1e1 / t29 / t47;
384 t184 = -0.41979259259259259259e1 * t10 * sigma[0] * t171 - 0.39864814814814814815e1 * t19 * t20 * t176 - 0.64162389091777094929e-3 * t28 * t181;
385 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
386 v3rho3[0] = -t6 * t143 * t34 / 0.9e1 + t6 * t147 * t63 / 0.15e2 - 0.7e1 / 0.75e2 * t6 * t8 * t86 * t88 + t6 * t81 * t107 / 0.10e2 + 0.203e3 / 0.5625e4 * t39 * t161 * t162 - 0.7e1 / 0.125e3 * t114 * t167 + t39 * t46 * t184 / 0.50e2;
387
388 t188 = t77 * t2;
389 t191 = t111 * t5;
390 t196 = t160 * t74;
391 t197 = t196 * t88;
392 t200 = t86 * t126;
393 t201 = t200 * t63;
394 t204 = t115 * t107;
395 t214 = 0.89955555555555555554e0 * t10 * t93 + 0.10872222222222222222e1 * t19 * sigma[0] * t99 + 0.19248716727533128478e-3 * t20 * t104;
396 t215 = t46 * t214;
397 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
398 v3rho2sigma[0] = t188 * t75 / 0.45e2 - 0.14e2 / 0.225e3 * t191 * t116 + t111 * t127 / 0.15e2 + 0.203e3 / 0.5625e4 * t114 * t197 - 0.14e2 / 0.375e3 * t114 * t201 - 0.7e1 / 0.375e3 * t114 * t204 + t39 * t215 / 0.50e2;
399
400 t220 = t160 * t130;
401 t221 = t220 * t63;
402 t224 = t115 * t126;
403 t229 = t86 * t138;
404 t230 = t229 * t63;
405 t237 = -0.17166666666666666667e0 * t19 * t55 - 0.42774926061184729952e-4 * sigma[0] * t60;
406 t238 = t46 * t237;
407 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
408 v3rhosigma2[0] = -0.7e1 / 0.225e3 * t111 * t131 + 0.203e3 / 0.5625e4 * t114 * t221 - 0.14e2 / 0.375e3 * t114 * t224 + t111 * t139 / 0.30e2 - 0.7e1 / 0.375e3 * t114 * t230 + t39 * t238 / 0.50e2;
409
410 t241 = t130 * t74;
411 t242 = t161 * t241;
412 t245 = t115 * t138;
413 t248 = t55 * t2;
414 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
415 v3sigma3[0] = 0.203e3 / 0.5625e4 * t39 * t242 - 0.7e1 / 0.125e3 * t114 * t245 + 0.10693731515296182488e-6 * t248 * t46;
416
417 #ifndef XC_DONT_COMPILE_LXC
418
419 if(order < 4) return;
420
421
422 t278 = 0.1e1 / t44 / t158 / t33;
423 t279 = t5 * t278;
424 t280 = t88 * t88;
425 t288 = t107 * t107;
426 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
427 v4rho4[0] = 0.4e1 / 0.27e2 * t6 / t7 / t11 * t34 - 0.4e1 / 0.135e3 * t6 * t143 * t45 * t63 - 0.28e2 / 0.225e3 * t6 * t77 * t86 * t88 + 0.2e1 / 0.15e2 * t6 * t147 * t107 + 0.812e3 / 0.3375e4 * t6 * t8 * t160 * t162 - 0.28e2 / 0.75e2 * t191 * t167 + 0.2e1 / 0.15e2 * t6 * t81 * t184 - 0.8932e4 / 0.84375e5 * t39 * t279 * t280 + 0.406e3 / 0.1875e4 * t114 * t160 * t88 * t107 - 0.7e1 / 0.125e3 * t39 * t87 * t288 - 0.28e2 / 0.375e3 * t114 * t166 * t184 + t39 * t46 * (0.23788246913580246913e2 * t10 * sigma[0] / t8 / t53 + 0.33220679012345679012e2 * t19 * t20 / t7 / t59 + 0.70578628000954804422e-2 * t28 / t29 / t21) / 0.50e2;
428
429 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
430 v4rho3sigma[0] = -t143 * t2 * t75 / 0.135e3 - 0.14e2 / 0.225e3 * t188 * t5 * t116 - 0.14e2 / 0.75e2 * t191 * t201 - 0.7e1 / 0.75e2 * t191 * t204 + 0.203e3 / 0.1875e4 * t114 * t160 * t126 * t88 - 0.7e1 / 0.125e3 * t114 * t86 * t214 * t63 - 0.7e1 / 0.125e3 * t114 * t200 * t107 - 0.7e1 / 0.375e3 * t114 * t115 * t184 + 0.203e3 / 0.1125e4 * t191 * t197 - 0.8932e4 / 0.84375e5 * t114 * t278 * t74 * t162 + 0.203e3 / 0.1875e4 * t114 * t196 * t63 * t107 + t188 * t127 / 0.15e2 + t111 * t215 / 0.10e2 + t39 * t46 * (-0.41979259259259259259e1 * t10 * t171 - 0.79729629629629629628e1 * t19 * sigma[0] * t176 - 0.19248716727533128478e-2 * t20 * t181) / 0.50e2;
431
432 t379 = t126 * t126;
433 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
434 v4rho2sigma2[0] = -0.14e2 / 0.675e3 * t188 * t131 + 0.406e3 / 0.3375e4 * t191 * t221 - 0.28e2 / 0.225e3 * t191 * t224 - 0.8932e4 / 0.84375e5 * t114 * t278 * t130 * t88 + 0.812e3 / 0.5625e4 * t114 * t196 * t63 * t126 + 0.203e3 / 0.5625e4 * t114 * t220 * t107 - 0.14e2 / 0.375e3 * t39 * t87 * t379 - 0.14e2 / 0.375e3 * t114 * t115 * t214 + t188 * t139 / 0.45e2 - 0.14e2 / 0.225e3 * t191 * t230 + t111 * t238 / 0.15e2 + 0.203e3 / 0.5625e4 * t114 * t160 * t138 * t88 - 0.14e2 / 0.375e3 * t114 * t86 * t237 * t63 - 0.7e1 / 0.375e3 * t114 * t229 * t107 + t39 * t46 * (0.10872222222222222222e1 * t19 * t99 + 0.38497433455066256957e-3 * sigma[0] * t104) / 0.50e2;
435
436 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
437 v4rhosigma3[0] = 0.203e3 / 0.3375e4 * t111 * t242 - 0.8932e4 / 0.84375e5 * t114 * t278 * t241 * t63 + 0.203e3 / 0.1875e4 * t114 * t220 * t126 - 0.7e1 / 0.75e2 * t191 * t245 + 0.203e3 / 0.1875e4 * t114 * t196 * t138 * t63 - 0.7e1 / 0.125e3 * t114 * t200 * t138 - 0.7e1 / 0.125e3 * t114 * t115 * t237 - 0.67726966263542489091e-6 * t99 * t2 * t46 - 0.99808160809431036555e-7 * t248 * t87 * t63;
438
439 t438 = t130 * t130;
440 t445 = t138 * t138;
441 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
442 v4sigma4[0] = -0.8932e4 / 0.84375e5 * t39 * t279 * t438 + 0.406e3 / 0.1875e4 * t114 * t220 * t138 - 0.7e1 / 0.125e3 * t39 * t87 * t445 - 0.39923264323772414622e-6 * t248 * t87 * t74;
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, t2, t3, t4, t5, t6, t7, t8;
469 double t9, t10, t12, t13, t14, t15, t16, t17;
470 double t20, t21, t22, t23, t25, t31, t32, t33;
471 double t34, t36, t40, t41, t42, t45, t46, t48;
472 double t50, t51, t52, t53, t55, t56, t57, t59;
473 double t63, t64, t65, t67, t71, t72, t73, t76;
474 double t77, t79;
475
476 #ifndef XC_DONT_COMPILE_VXC
477 double t81, t82, t83, t84, t85, t86, t87, t89;
478 double t91, t92, t93, t96, t97, t98, t99, t100;
479 double t102, t103, t104, t105, t106, t108, t112, t114;
480 double t118, t119, t122, t124, t125, t126, t127, t128;
481 double t130, t131, t134, t135, t139, t141, t142, t143;
482 double t145, t146, t147, t148, t149, t151, t152, t153;
483 double t154, t155, t157, t161, t163, t167, t168, t171;
484 double t173, t174, t178, t179, t180, t188, t191, t199;
485
486 #ifndef XC_DONT_COMPILE_FXC
487 double t202, t204, t205, t206, t207, t209, t210, t211;
488 double t213, t214, t215, t217, t218, t219, t220, t221;
489 double t222, t224, t225, t227, t230, t231, t232, t234;
490 double t235, t237, t238, t239, t241, t242, t244, t250;
491 double t255, t258, t260, t261, t262, t263, t264, t266;
492 double t267, t268, t270, t271, t272, t274, t275, t278;
493 double t279, t280, t282, t284, t287, t289, t291, t294;
494 double t295, t297, t301, t303, t306, t308, t310, t313;
495 double t315, t317, t319, t321, t324, t326, t327, t328;
496 double t329, t331, t332, t333, t335, t336, t337, t338;
497 double t340, t341, t342, t344, t345, t346, t348, t349;
498 double t350, t352, t353, t355, t361, t366, t369, t371;
499 double t372, t373, t376, t378, t379, t383, t384, t395;
500 double t400, t401, t402, t407, t412, t415, t417, t427;
501 double t428, t431, t432, t439, t440, t443, t444, t451;
502 double t452;
503
504 #ifndef XC_DONT_COMPILE_KXC
505 double t458, t462, t467, t468, t471, t472, t473, t475;
506 double t476, t477, t479, t480, t481, t483, t484, t486;
507 double t487, t488, t490, t491, t492, t494, t496, t498;
508 double t500, t501, t503, t504, t505, t506, t508, t510;
509 double t511, t512, t513, t515, t516, t518, t521, t526;
510 double t531, t534, t536, t537, t539, t540, t542, t543;
511 double t545, t546, t547, t549, t550, t552, t554, t556;
512 double t557, t558, t560, t562, t564, t565, t566, t568;
513 double t569, t571, t572, t573, t575, t576, t578, t580;
514 double t582, t583, t584, t589, t591, t592, t593, t594;
515 double t596, t597, t599, t602, t604, t606, t609, t611;
516 double t613, t617, t621, t624, t625, t626, t627, t630;
517 double t631, t632, t634, t637, t638, t639, t641, t643;
518 double t645, t647, t649, t651, t653, t655, t656, t658;
519 double t661, t662, t663, t665, t667, t669, t671, t672;
520 double t674, t676, t678, t681, t682, t683, t685, t697;
521 double t700, t702, t705, t707, t709, t711, t714, t716;
522 double t718, t720, t723, t726, t729, t732, t734, t736;
523 double t739, t746, t748, t749, t750, t751, t753, t755;
524 double t758, t760, t762, t764, t766, t768, t770, t772;
525 double t774, t776, t778, t780, t782, t784, t785, t787;
526 double t789, t791, t796, t799, t808, t810, t811, t812;
527 double t814, t815, t817, t818, t819, t821, t822, t824;
528 double t825, t826, t828, t829, t830, t832, t833, t835;
529 double t837, t838, t839, t841, t842, t844, t849, t854;
530 double t857, t859, t860, t861, t863, t864, t865, t867;
531 double t868, t869, t870, t872, t873, t874, t876, t877;
532 double t882, t883, t884, t886, t890, t894, t896, t899;
533 double t901, t902, t906, t907, t908, t909, t912, t916;
534 double t920, t921, t925, t929, t940, t945, t949, t950;
535 double t951, t955, t963, t964, t966, t967, t968, t969;
536 double t972, t975, t976, t979, t980, t982, t986, t988;
537 double t989, t990, t991, t995, t998, t1000, t1001, t1002;
538 double t1003, t1006, t1008, t1009, t1014, t1018, t1025, t1028;
539 double t1031, t1032, t1036, t1039, t1041, t1044, t1045, t1049;
540 double t1059, t1060, t1064, t1065, t1067, t1071, t1075, t1079;
541 double t1080, t1081, t1084, t1085, t1092, t1093, t1097, t1098;
542 double t1100, t1104, t1105, t1106, t1110, t1113, t1114, t1118;
543 double t1122, t1126, t1129, t1130, t1133, t1134, t1141, t1142;
544 double t1145, t1146, t1150, t1153, t1156, t1157, t1161, t1164;
545
546 #ifndef XC_DONT_COMPILE_LXC
547 double t1170, t1174, t1176, t1187, t1190, t1193, t1204, t1210;
548 double t1220, t1243, t1254, t1257, t1259, t1262, t1283, t1285;
549 double t1287, t1289, t1293, t1294, t1297, t1301, t1304, t1306;
550 double t1309, t1310, t1315, t1318, t1320, t1323, t1326, t1330;
551 double t1332, t1336, t1340, t1343, t1344, t1349, t1352, t1354;
552 double t1358, t1364, t1365, t1367, t1368, t1371, t1380, t1384;
553 double t1385, t1396, t1410, t1413, t1417, t1424, t1431, t1432;
554 double t1434, t1437, t1450, t1457, t1460, t1461, t1467, t1469;
555 double t1471, t1472, t1474, t1488, t1491, t1494, t1501, t1502;
556 double t1504, t1525, t1539, t1542, t1544, t1545, t1549, t1551;
557 double t1552, t1554, t1555, t1560, t1572, t1580, t1585, t1586;
558 double t1598, t1604, t1608, t1612, t1621, t1622, t1623, t1628;
559 double t1630, t1635, t1639, t1646, t1651, t1654, t1661, t1672;
560 double t1674, t1684, t1691, t1698, t1705, t1715, t1722, t1733;
561 double t1744, t1758, t1761, t1776, t1787, t1789, t1796, t1799;
562 double t1801, t1803, t1806, t1809, t1816, t1820, t1823, t1831;
563 double t1844, t1858, t1881, t1884, t1900, t1911, t1927, t1943;
564 double t1969, t1975, t1982, t1988, t1991, t1994, t1997, t2000;
565 double t2002, t2065, t2077, t2080, t2083, t2086, t2089, t2092;
566 double t2096, t2100, t2119, t2124, t2130, t2139, t2149, t2154;
567 double t2194, t2226, t2229, t2234, t2236, t2240, t2244, t2248;
568 double t2251, t2256, t2284, t2298, t2314, t2326, t2328, t2331;
569 double t2334, t2337, t2352, t2355, t2358, t2360, t2365, t2382;
570 double t2392, t2394, t2402, t2413, t2414, t2423, t2425, t2428;
571 double t2441, t2443, t2446, t2449, t2453, t2454, t2467, t2468;
572 double t2471, t2472, t2486, t2493, t2494, t2496, t2511, t2515;
573 double t2536, t2541, t2546, t2549, t2553, t2556, t2560, t2564;
574 double t2567, t2576, t2588, t2603, t2623, t2625, t2627, t2629;
575 double t2631, t2633, t2650, t2659, t2663, t2670, t2683, t2685;
576 double t2704, t2743, t2766, t2787, t2788, t2805, t2807, t2809;
577 double t2823, t2825, t2827, t2844, t2854, t2859, t2860, t2865;
578 double t2866, t2876, t2878, t2879, t2881, t2891, t2892, t2905;
579 double t2924, t2942, t2952, t2954, t2958, t2964, t2974, t2976;
580 double t2980, t3049, t3075, t3079, t3081, t3082, t3096, t3112;
581 double t3125, t3126, t3132, t3137, t3190, t3198, t3206, t3214;
582 #endif
583
584 #endif
585
586 #endif
587
588 #endif
589
590 #endif
591
592
593
594 t1 = M_CBRT6;
595 t2 = t1 * t1;
596 t3 = M_PI * M_PI;
597 t4 = POW_1_3(t3);
598 t5 = t4 * t4;
599 t6 = t2 * t5;
600 t7 = rho[0] - rho[1];
601 t8 = rho[0] + rho[1];
602 t9 = 0.1e1 / t8;
603 t10 = t7 * t9;
604 t12 = 0.1e1 / 0.2e1 + t10 / 0.2e1;
605 t13 = POW_1_3(t12);
606 t14 = t13 * t13;
607 t15 = t14 * t12;
608 t16 = POW_1_3(t8);
609 t17 = t16 * t16;
610 t20 = t1 / t5;
611 t21 = rho[0] * rho[0];
612 t22 = POW_1_3(rho[0]);
613 t23 = t22 * t22;
614 t25 = 0.1e1 / t23 / t21;
615 t31 = t2 / t4 / t3;
616 t32 = sigma[0] * sigma[0];
617 t33 = t21 * t21;
618 t34 = t33 * rho[0];
619 t36 = 0.1e1 / t22 / t34;
620 t40 = t32 * sigma[0];
621 t41 = t33 * t33;
622 t42 = 0.1e1 / t41;
623 t45 = 0.1e1 + 0.91999999999999999998e-1 * t20 * sigma[0] * t25 + 0.16093750000000000000e-1 * t31 * t32 * t36 + 0.89114429294134854068e-6 * t40 * t42;
624 t46 = pow(t45, 0.1e1 / 0.15e2);
625 t48 = t6 * t15 * t17 * t46;
626 t50 = 0.1e1 / 0.2e1 - t10 / 0.2e1;
627 t51 = POW_1_3(t50);
628 t52 = t51 * t51;
629 t53 = t52 * t50;
630 t55 = rho[1] * rho[1];
631 t56 = POW_1_3(rho[1]);
632 t57 = t56 * t56;
633 t59 = 0.1e1 / t57 / t55;
634 t63 = sigma[2] * sigma[2];
635 t64 = t55 * t55;
636 t65 = t64 * rho[1];
637 t67 = 0.1e1 / t56 / t65;
638 t71 = t63 * sigma[2];
639 t72 = t64 * t64;
640 t73 = 0.1e1 / t72;
641 t76 = 0.1e1 + 0.91999999999999999998e-1 * t20 * sigma[2] * t59 + 0.16093750000000000000e-1 * t31 * t63 * t67 + 0.89114429294134854068e-6 * t71 * t73;
642 t77 = pow(t76, 0.1e1 / 0.15e2);
643 t79 = t6 * t53 * t17 * t77;
644 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
645 zk[0] = 0.3e1 / 0.10e2 * t48 + 0.3e1 / 0.10e2 * t79;
646
647 #ifndef XC_DONT_COMPILE_VXC
648
649 if(order < 1) return;
650
651
652 t81 = 0.3e1 / 0.10e2 * t48;
653 t82 = 0.3e1 / 0.10e2 * t79;
654 t83 = t6 * t14;
655 t84 = t17 * t46;
656 t85 = t8 * t8;
657 t86 = 0.1e1 / t85;
658 t87 = t7 * t86;
659 t89 = t9 / 0.2e1 - t87 / 0.2e1;
660 t91 = t83 * t84 * t89;
661 t92 = t91 / 0.2e1;
662 t93 = 0.1e1 / t16;
663 t96 = t6 * t15 * t93 * t46;
664 t97 = t96 / 0.5e1;
665 t98 = t6 * t15;
666 t99 = t46 * t46;
667 t100 = t99 * t99;
668 t102 = t100 * t100;
669 t103 = t102 * t100 * t99;
670 t104 = 0.1e1 / t103;
671 t105 = t17 * t104;
672 t106 = t21 * rho[0];
673 t108 = 0.1e1 / t23 / t106;
674 t112 = t33 * t21;
675 t114 = 0.1e1 / t22 / t112;
676 t118 = t41 * rho[0];
677 t119 = 0.1e1 / t118;
678 t122 = -0.24533333333333333333e0 * t20 * sigma[0] * t108 - 0.85833333333333333333e-1 * t31 * t32 * t114 - 0.71291543435307883254e-5 * t40 * t119;
679 t124 = t98 * t105 * t122;
680 t125 = t124 / 0.50e2;
681 t126 = t6 * t52;
682 t127 = t17 * t77;
683 t128 = -t89;
684 t130 = t126 * t127 * t128;
685 t131 = t130 / 0.2e1;
686 t134 = t6 * t53 * t93 * t77;
687 t135 = t134 / 0.5e1;
688 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
689 vrho[0] = t81 + t82 + t8 * (t92 + t97 + t125 + t131 + t135);
690
691 t139 = -t9 / 0.2e1 - t87 / 0.2e1;
692 t141 = t83 * t84 * t139;
693 t142 = t141 / 0.2e1;
694 t143 = -t139;
695 t145 = t126 * t127 * t143;
696 t146 = t145 / 0.2e1;
697 t147 = t6 * t53;
698 t148 = t77 * t77;
699 t149 = t148 * t148;
700 t151 = t149 * t149;
701 t152 = t151 * t149 * t148;
702 t153 = 0.1e1 / t152;
703 t154 = t17 * t153;
704 t155 = t55 * rho[1];
705 t157 = 0.1e1 / t57 / t155;
706 t161 = t64 * t55;
707 t163 = 0.1e1 / t56 / t161;
708 t167 = t72 * rho[1];
709 t168 = 0.1e1 / t167;
710 t171 = -0.24533333333333333333e0 * t20 * sigma[2] * t157 - 0.85833333333333333333e-1 * t31 * t63 * t163 - 0.71291543435307883254e-5 * t71 * t168;
711 t173 = t147 * t154 * t171;
712 t174 = t173 / 0.50e2;
713 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
714 vrho[1] = t81 + t82 + t8 * (t142 + t97 + t146 + t135 + t174);
715
716 t178 = t17 * t8 * t2;
717 t179 = t178 * t5;
718 t180 = t15 * t104;
719 t188 = 0.91999999999999999998e-1 * t20 * t25 + 0.32187500000000000000e-1 * t31 * sigma[0] * t36 + 0.26734328788240456220e-5 * t32 * t42;
720 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
721 vsigma[0] = t179 * t180 * t188 / 0.50e2;
722
723 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
724 vsigma[1] = 0.0e0;
725
726 t191 = t53 * t153;
727 t199 = 0.91999999999999999998e-1 * t20 * t59 + 0.32187500000000000000e-1 * t31 * sigma[2] * t67 + 0.26734328788240456220e-5 * t63 * t73;
728 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
729 vsigma[2] = t179 * t191 * t199 / 0.50e2;
730
731 #ifndef XC_DONT_COMPILE_FXC
732
733 if(order < 2) return;
734
735
736 t202 = 0.2e1 / 0.5e1 * t96;
737 t204 = 0.2e1 / 0.5e1 * t134;
738 t205 = 0.1e1 / t13;
739 t206 = t6 * t205;
740 t207 = t89 * t89;
741 t209 = t206 * t84 * t207;
742 t210 = t209 / 0.3e1;
743 t211 = t93 * t46;
744 t213 = t83 * t211 * t89;
745 t214 = 0.2e1 / 0.3e1 * t213;
746 t215 = t89 * t122;
747 t217 = t83 * t105 * t215;
748 t218 = t217 / 0.15e2;
749 t219 = t85 * t8;
750 t220 = 0.1e1 / t219;
751 t221 = t7 * t220;
752 t222 = -t86 + t221;
753 t224 = t83 * t84 * t222;
754 t225 = t224 / 0.2e1;
755 t227 = 0.1e1 / t16 / t8;
756 t230 = t6 * t15 * t227 * t46;
757 t231 = t230 / 0.15e2;
758 t232 = t93 * t104;
759 t234 = t98 * t232 * t122;
760 t235 = 0.2e1 / 0.75e2 * t234;
761 t237 = 0.1e1 / t103 / t45;
762 t238 = t17 * t237;
763 t239 = t122 * t122;
764 t241 = t98 * t238 * t239;
765 t242 = 0.7e1 / 0.375e3 * t241;
766 t244 = 0.1e1 / t23 / t33;
767 t250 = 0.1e1 / t22 / t33 / t106;
768 t255 = 0.1e1 / t41 / t21;
769 t258 = 0.89955555555555555554e0 * t20 * sigma[0] * t244 + 0.54361111111111111111e0 * t31 * t32 * t250 + 0.64162389091777094929e-4 * t40 * t255;
770 t260 = t98 * t105 * t258;
771 t261 = t260 / 0.50e2;
772 t262 = 0.1e1 / t51;
773 t263 = t6 * t262;
774 t264 = t128 * t128;
775 t266 = t263 * t127 * t264;
776 t267 = t266 / 0.3e1;
777 t268 = t93 * t77;
778 t270 = t126 * t268 * t128;
779 t271 = 0.2e1 / 0.3e1 * t270;
780 t272 = -t222;
781 t274 = t126 * t127 * t272;
782 t275 = t274 / 0.2e1;
783 t278 = t6 * t53 * t227 * t77;
784 t279 = t278 / 0.15e2;
785 t280 = t210 + t214 + t218 + t225 - t231 + t235 - t242 + t261 + t267 + t271 + t275 - t279;
786 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
787 v2rho2[0] = t91 + t202 + t124 / 0.25e2 + t130 + t204 + t8 * t280;
788
789 t282 = t139 * t89;
790 t284 = t206 * t84 * t282;
791 t287 = t83 * t211 * t139;
792 t289 = t139 * t122;
793 t291 = t83 * t105 * t289;
794 t294 = 0.1e1 / t16 / t85;
795 t295 = t294 * t46;
796 t297 = t83 * t295 * t7;
797 t301 = t143 * t128;
798 t303 = t263 * t127 * t301;
799 t306 = t126 * t268 * t143;
800 t308 = t294 * t77;
801 t310 = t126 * t308 * t7;
802 t313 = t171 * t128;
803 t315 = t126 * t154 * t313;
804 t317 = t93 * t153;
805 t319 = t147 * t317 * t171;
806 t321 = t284 / 0.3e1 + t287 / 0.3e1 + t291 / 0.30e2 + t297 / 0.2e1 + t213 / 0.3e1 - t231 + t234 / 0.75e2 + t303 / 0.3e1 + t306 / 0.3e1 - t310 / 0.2e1 + t270 / 0.3e1 - t279 + t315 / 0.30e2 + t319 / 0.75e2;
807 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
808 v2rho2[1] = t8 * t321 + t125 + t131 + t142 + t146 + t174 + t202 + t204 + t92;
809
810 t324 = t139 * t139;
811 t326 = t206 * t84 * t324;
812 t327 = t326 / 0.3e1;
813 t328 = 0.2e1 / 0.3e1 * t287;
814 t329 = t86 + t221;
815 t331 = t83 * t84 * t329;
816 t332 = t331 / 0.2e1;
817 t333 = t143 * t143;
818 t335 = t263 * t127 * t333;
819 t336 = t335 / 0.3e1;
820 t337 = 0.2e1 / 0.3e1 * t306;
821 t338 = t143 * t171;
822 t340 = t126 * t154 * t338;
823 t341 = t340 / 0.15e2;
824 t342 = -t329;
825 t344 = t126 * t127 * t342;
826 t345 = t344 / 0.2e1;
827 t346 = 0.2e1 / 0.75e2 * t319;
828 t348 = 0.1e1 / t152 / t76;
829 t349 = t17 * t348;
830 t350 = t171 * t171;
831 t352 = t147 * t349 * t350;
832 t353 = 0.7e1 / 0.375e3 * t352;
833 t355 = 0.1e1 / t57 / t64;
834 t361 = 0.1e1 / t56 / t64 / t155;
835 t366 = 0.1e1 / t72 / t55;
836 t369 = 0.89955555555555555554e0 * t20 * sigma[2] * t355 + 0.54361111111111111111e0 * t31 * t63 * t361 + 0.64162389091777094929e-4 * t71 * t366;
837 t371 = t147 * t154 * t369;
838 t372 = t371 / 0.50e2;
839 t373 = t327 + t328 + t332 - t231 + t336 + t337 + t341 + t345 - t279 + t346 - t353 + t372;
840 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
841 v2rho2[2] = t141 + t202 + t145 + t204 + t173 / 0.25e2 + t8 * t373;
842
843 t376 = t98 * t105 * t188;
844 t378 = t14 * t104;
845 t379 = t188 * t89;
846 t383 = t15 * t237;
847 t384 = t188 * t122;
848 t395 = -0.24533333333333333333e0 * t20 * t108 - 0.17166666666666666667e0 * t31 * sigma[0] * t114 - 0.21387463030592364976e-4 * t32 * t119;
849 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
850 v2rhosigma[0] = t376 / 0.30e2 + t179 * t378 * t379 / 0.30e2 - 0.7e1 / 0.375e3 * t179 * t383 * t384 + t179 * t180 * t395 / 0.50e2;
851
852 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
853 v2rhosigma[1] = 0.0e0;
854
855 t400 = t147 * t154 * t199;
856 t401 = t52 * t153;
857 t402 = t199 * t128;
858 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
859 v2rhosigma[2] = t179 * t401 * t402 / 0.30e2 + t400 / 0.30e2;
860
861 t407 = t378 * t188 * t139;
862 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
863 v2rhosigma[3] = t179 * t407 / 0.30e2 + t376 / 0.30e2;
864
865 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
866 v2rhosigma[4] = 0.0e0;
867
868 t412 = t401 * t199 * t143;
869 t415 = t53 * t348;
870 t417 = t415 * t199 * t171;
871 t427 = -0.24533333333333333333e0 * t20 * t157 - 0.17166666666666666667e0 * t31 * sigma[2] * t163 - 0.21387463030592364976e-4 * t63 * t168;
872 t428 = t191 * t427;
873 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
874 v2rhosigma[5] = t400 / 0.30e2 + t179 * t412 / 0.30e2 - 0.7e1 / 0.375e3 * t179 * t417 + t179 * t428 / 0.50e2;
875
876 t431 = t188 * t188;
877 t432 = t383 * t431;
878 t439 = 0.32187500000000000000e-1 * t31 * t36 + 0.53468657576480912440e-5 * sigma[0] * t42;
879 t440 = t180 * t439;
880 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
881 v2sigma2[0] = -0.7e1 / 0.375e3 * t179 * t432 + t179 * t440 / 0.50e2;
882
883 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
884 v2sigma2[1] = 0.0e0;
885
886 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
887 v2sigma2[2] = 0.0e0;
888
889 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
890 v2sigma2[3] = 0.0e0;
891
892 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
893 v2sigma2[4] = 0.0e0;
894
895 t443 = t199 * t199;
896 t444 = t415 * t443;
897 t451 = 0.32187500000000000000e-1 * t31 * t67 + 0.53468657576480912440e-5 * sigma[2] * t73;
898 t452 = t191 * t451;
899 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
900 v2sigma2[5] = -0.7e1 / 0.375e3 * t179 * t444 + t179 * t452 / 0.50e2;
901
902 #ifndef XC_DONT_COMPILE_KXC
903
904 if(order < 3) return;
905
906
907 t458 = t230 / 0.5e1;
908 t462 = t278 / 0.5e1;
909 t467 = t6 * t15 * t294 * t46;
910 t468 = 0.4e1 / 0.45e2 * t467;
911 t471 = t6 * t53 * t294 * t77;
912 t472 = 0.4e1 / 0.45e2 * t471;
913 t473 = t89 * t258;
914 t475 = t83 * t105 * t473;
915 t476 = t475 / 0.10e2;
916 t477 = t122 * t258;
917 t479 = t98 * t238 * t477;
918 t480 = 0.7e1 / 0.125e3 * t479;
919 t481 = t128 * t272;
920 t483 = t263 * t127 * t481;
921 t484 = t89 * t239;
922 t486 = t83 * t238 * t484;
923 t487 = 0.7e1 / 0.75e2 * t486;
924 t488 = t207 * t122;
925 t490 = t206 * t105 * t488;
926 t491 = t490 / 0.15e2;
927 t492 = t89 * t222;
928 t494 = t206 * t84 * t492;
929 t496 = t83 * t232 * t215;
930 t498 = t222 * t122;
931 t500 = t83 * t105 * t498;
932 t501 = t500 / 0.10e2;
933 t503 = t83 * t211 * t222;
934 t504 = t85 * t85;
935 t505 = 0.1e1 / t504;
936 t506 = t7 * t505;
937 t508 = 0.3e1 * t220 - 0.3e1 * t506;
938 t510 = t83 * t84 * t508;
939 t511 = t510 / 0.2e1;
940 t512 = t468 + t472 + t476 - t480 + t483 - t487 + t491 + t494 + 0.2e1 / 0.15e2 * t496 + t501 + t503 + t511;
941 t513 = t227 * t104;
942 t515 = t98 * t513 * t122;
943 t516 = t515 / 0.75e2;
944 t518 = t98 * t232 * t258;
945 t521 = 0.1e1 / t23 / t34;
946 t526 = 0.1e1 / t22 / t41;
947 t531 = 0.1e1 / t41 / t106;
948 t534 = -0.41979259259259259259e1 * t20 * sigma[0] * t521 - 0.39864814814814814815e1 * t31 * t32 * t526 - 0.64162389091777094929e-3 * t40 * t531;
949 t536 = t98 * t105 * t534;
950 t537 = t536 / 0.50e2;
951 t539 = t126 * t268 * t272;
952 t540 = -t508;
953 t542 = t126 * t127 * t540;
954 t543 = t542 / 0.2e1;
955 t545 = 0.1e1 / t13 / t12;
956 t546 = t6 * t545;
957 t547 = t207 * t89;
958 t549 = t546 * t84 * t547;
959 t550 = t549 / 0.9e1;
960 t552 = t206 * t211 * t207;
961 t554 = t227 * t46;
962 t556 = t83 * t554 * t89;
963 t557 = t556 / 0.3e1;
964 t558 = t93 * t237;
965 t560 = t98 * t558 * t239;
966 t562 = t45 * t45;
967 t564 = 0.1e1 / t103 / t562;
968 t565 = t17 * t564;
969 t566 = t239 * t122;
970 t568 = t98 * t565 * t566;
971 t569 = 0.203e3 / 0.5625e4 * t568;
972 t571 = 0.1e1 / t51 / t50;
973 t572 = t6 * t571;
974 t573 = t264 * t128;
975 t575 = t572 * t127 * t573;
976 t576 = t575 / 0.9e1;
977 t578 = t263 * t268 * t264;
978 t580 = t227 * t77;
979 t582 = t126 * t580 * t128;
980 t583 = t582 / 0.3e1;
981 t584 = -t516 + t518 / 0.25e2 + t537 + t539 + t543 - t550 + 0.2e1 / 0.3e1 * t552 - t557 - 0.14e2 / 0.375e3 * t560 + t569 - t576 + 0.2e1 / 0.3e1 * t578 - t583;
982 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
983 v3rho3[0] = t209 + 0.2e1 * t213 + t217 / 0.5e1 + 0.3e1 / 0.2e1 * t224 - t458 + 0.2e1 / 0.25e2 * t234 + 0.3e1 / 0.50e2 * t260 + 0.3e1 / 0.2e1 * t274 - t462 - 0.7e1 / 0.125e3 * t241 + t266 + 0.2e1 * t270 + t8 * (t512 + t584);
984
985 t589 = t210 + 0.4e1 / 0.3e1 * t213 + t218 + t225 - t458 + 0.4e1 / 0.75e2 * t234 + t261 + t275 - t462 - t310 - t242;
986 t591 = 0.2e1 / 0.3e1 * t284;
987 t592 = t291 / 0.15e2;
988 t593 = 0.2e1 / 0.3e1 * t303;
989 t594 = t315 / 0.15e2;
990 t596 = 0.1e1 / t16 / t219;
991 t597 = t596 * t77;
992 t599 = t126 * t597 * t7;
993 t602 = t83 * t554 * t139;
994 t604 = t596 * t46;
995 t606 = t83 * t604 * t7;
996 t609 = t126 * t580 * t143;
997 t611 = t227 * t153;
998 t613 = t147 * t611 * t171;
999 t617 = t6 * t14 * t294 * t46;
1000 t621 = t6 * t52 * t294 * t77;
1001 t624 = t6 * t205 * t17;
1002 t625 = t104 * t139;
1003 t626 = t625 * t215;
1004 t627 = t624 * t626;
1005 t630 = t126 * t317 * t313;
1006 t631 = 0.2e1 / 0.45e2 * t630;
1007 t632 = t139 * t222;
1008 t634 = t206 * t84 * t632;
1009 t637 = t83 * t232 * t289;
1010 t638 = 0.2e1 / 0.45e2 * t637;
1011 t639 = t139 * t258;
1012 t641 = t83 * t105 * t639;
1013 t643 = t143 * t272;
1014 t645 = t263 * t127 * t643;
1015 t647 = t7 * t128;
1016 t649 = t263 * t308 * t647;
1017 t651 = t171 * t272;
1018 t653 = t126 * t154 * t651;
1019 t655 = 0.5e1 / 0.6e1 * t599 - t602 / 0.9e1 - 0.5e1 / 0.6e1 * t606 - t609 / 0.9e1 - t613 / 0.225e3 + t617 / 0.2e1 - t621 / 0.2e1 + t468 + t472 + 0.2e1 / 0.45e2 * t627 + t631 + t634 / 0.3e1 + t638 + t641 / 0.30e2 + t645 / 0.3e1 - 0.2e1 / 0.3e1 * t649 + t653 / 0.30e2;
1020 t656 = t139 * t207;
1021 t658 = t546 * t84 * t656;
1022 t661 = t206 * t211 * t282;
1023 t662 = 0.4e1 / 0.9e1 * t661;
1024 t663 = t7 * t89;
1025 t665 = t206 * t295 * t663;
1026 t667 = t139 * t239;
1027 t669 = t83 * t238 * t667;
1028 t671 = t294 * t104;
1029 t672 = t7 * t122;
1030 t674 = t83 * t671 * t672;
1031 t676 = t143 * t264;
1032 t678 = t572 * t127 * t676;
1033 t681 = t263 * t268 * t301;
1034 t682 = 0.4e1 / 0.9e1 * t681;
1035 t683 = t171 * t264;
1036 t685 = t263 * t154 * t683;
1037 t697 = -t658 / 0.9e1 + t662 + 0.2e1 / 0.3e1 * t665 - 0.7e1 / 0.225e3 * t669 + t674 / 0.15e2 - t678 / 0.9e1 + t682 + t685 / 0.45e2 + 0.2e1 / 0.45e2 * t496 + t503 / 0.3e1 - 0.2e1 / 0.225e3 * t515 + t518 / 0.75e2 + t539 / 0.3e1 + 0.2e1 / 0.9e1 * t552 - 0.2e1 / 0.9e1 * t556 - 0.14e2 / 0.1125e4 * t560 + 0.2e1 / 0.9e1 * t578 - 0.2e1 / 0.9e1 * t582;
1038 t700 = t267 + 0.4e1 / 0.3e1 * t270 + t591 + t328 + t592 + t297 + t593 + t337 + t594 + t346 + t8 * (t655 + t697);
1039 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1040 v3rho3[1] = t589 + t700;
1041
1042 t702 = t214 - t458 + t235 - t462 - t310 + t332 + t345 + t372 + t271 + t591 + 0.4e1 / 0.3e1 * t287;
1043 t705 = t350 * t128;
1044 t707 = t126 * t349 * t705;
1045 t709 = t369 * t128;
1046 t711 = t126 * t154 * t709;
1047 t714 = -t220 - 0.3e1 * t506;
1048 t716 = t83 * t84 * t714;
1049 t718 = -t714;
1050 t720 = t126 * t127 * t718;
1051 t723 = t206 * t211 * t324;
1052 t726 = t83 * t211 * t329;
1053 t729 = t263 * t268 * t333;
1054 t732 = t126 * t268 * t342;
1055 t734 = t93 * t348;
1056 t736 = t147 * t734 * t350;
1057 t739 = t147 * t317 * t369;
1058 t746 = -0.7e1 / 0.225e3 * t707 + t711 / 0.30e2 + t716 / 0.2e1 + t720 / 0.2e1 + 0.2e1 / 0.9e1 * t723 + t726 / 0.3e1 + 0.2e1 / 0.9e1 * t729 + t732 / 0.3e1 - 0.14e2 / 0.1125e4 * t736 + t739 / 0.75e2 - 0.2e1 / 0.3e1 * t599 - 0.2e1 / 0.9e1 * t602 + 0.2e1 / 0.3e1 * t606 - 0.2e1 / 0.9e1 * t609 - 0.2e1 / 0.225e3 * t613 + t468 + t472;
1059 t748 = t6 * t262 * t17;
1060 t749 = t153 * t143;
1061 t750 = t749 * t313;
1062 t751 = t748 * t750;
1063 t753 = t139 * t7;
1064 t755 = t206 * t295 * t753;
1065 t758 = t126 * t317 * t338;
1066 t760 = t324 * t89;
1067 t762 = t546 * t84 * t760;
1068 t764 = t324 * t122;
1069 t766 = t206 * t105 * t764;
1070 t768 = t329 * t89;
1071 t770 = t206 * t84 * t768;
1072 t772 = t329 * t122;
1073 t774 = t83 * t105 * t772;
1074 t776 = t333 * t128;
1075 t778 = t572 * t127 * t776;
1076 t780 = t143 * t7;
1077 t782 = t263 * t308 * t780;
1078 t784 = t294 * t153;
1079 t785 = t7 * t171;
1080 t787 = t126 * t784 * t785;
1081 t789 = t342 * t128;
1082 t791 = t263 * t127 * t789;
1083 t796 = 0.2e1 / 0.45e2 * t751 + 0.2e1 / 0.3e1 * t755 + 0.2e1 / 0.45e2 * t758 - t762 / 0.9e1 + t766 / 0.45e2 + t770 / 0.3e1 + t774 / 0.30e2 - t778 / 0.9e1 - 0.2e1 / 0.3e1 * t782 - t787 / 0.15e2 + t791 / 0.3e1 + t631 + t638 + t662 + t682 - t515 / 0.225e3 - t556 / 0.9e1 - t582 / 0.9e1;
1084 t799 = t592 + t297 + t593 + 0.4e1 / 0.3e1 * t306 + t594 + 0.4e1 / 0.75e2 * t319 + t327 + t336 + t341 - t353 + t8 * (t746 + t796);
1085 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1086 v3rho3[2] = t702 + t799;
1087
1088 t808 = t143 * t350;
1089 t810 = t126 * t349 * t808;
1090 t811 = 0.7e1 / 0.75e2 * t810;
1091 t812 = t139 * t329;
1092 t814 = t206 * t84 * t812;
1093 t815 = t333 * t171;
1094 t817 = t263 * t154 * t815;
1095 t818 = t817 / 0.15e2;
1096 t819 = t143 * t342;
1097 t821 = t263 * t127 * t819;
1098 t822 = t342 * t171;
1099 t824 = t126 * t154 * t822;
1100 t825 = t824 / 0.10e2;
1101 t826 = t143 * t369;
1102 t828 = t126 * t154 * t826;
1103 t829 = t828 / 0.10e2;
1104 t830 = t171 * t369;
1105 t832 = t147 * t349 * t830;
1106 t833 = 0.7e1 / 0.125e3 * t832;
1107 t835 = -0.3e1 * t220 - 0.3e1 * t506;
1108 t837 = t83 * t84 * t835;
1109 t838 = t837 / 0.2e1;
1110 t839 = -t835;
1111 t841 = t126 * t127 * t839;
1112 t842 = t841 / 0.2e1;
1113 t844 = 0.1e1 / t57 / t65;
1114 t849 = 0.1e1 / t56 / t72;
1115 t854 = 0.1e1 / t72 / t155;
1116 t857 = -0.41979259259259259259e1 * t20 * sigma[2] * t844 - 0.39864814814814814815e1 * t31 * t63 * t849 - 0.64162389091777094929e-3 * t71 * t854;
1117 t859 = t147 * t154 * t857;
1118 t860 = t859 / 0.50e2;
1119 t861 = t324 * t139;
1120 t863 = t546 * t84 * t861;
1121 t864 = t863 / 0.9e1;
1122 t865 = t333 * t143;
1123 t867 = t572 * t127 * t865;
1124 t868 = t867 / 0.9e1;
1125 t869 = -t811 + t814 + t818 + t821 + t825 + t829 - t833 + t838 + t842 + t860 - t864 - t868;
1126 t870 = t76 * t76;
1127 t872 = 0.1e1 / t152 / t870;
1128 t873 = t17 * t872;
1129 t874 = t350 * t171;
1130 t876 = t147 * t873 * t874;
1131 t877 = 0.203e3 / 0.5625e4 * t876;
1132 t882 = t602 / 0.3e1;
1133 t883 = t609 / 0.3e1;
1134 t884 = t613 / 0.75e2;
1135 t886 = t877 + 0.2e1 / 0.3e1 * t723 + t726 + 0.2e1 / 0.3e1 * t729 + t732 - 0.14e2 / 0.375e3 * t736 + t739 / 0.25e2 - t882 - t883 - t884 + t468 + t472 + 0.2e1 / 0.15e2 * t758;
1136 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1137 v3rho3[3] = -t458 - t462 + 0.3e1 / 0.2e1 * t331 + 0.3e1 / 0.2e1 * t344 + 0.3e1 / 0.50e2 * t371 + 0.2e1 * t287 + 0.2e1 * t306 + 0.2e1 / 0.25e2 * t319 + t326 + t335 + t340 / 0.5e1 - 0.7e1 / 0.125e3 * t352 + t8 * (t869 + t886);
1138
1139 t890 = t83 * t105 * t379;
1140 t894 = t98 * t232 * t188 / 0.45e2;
1141 t896 = t98 * t238 * t384;
1142 t899 = t98 * t105 * t395;
1143 t901 = t205 * t104;
1144 t902 = t188 * t207;
1145 t906 = t5 * t14;
1146 t907 = t178 * t906;
1147 t908 = t237 * t188;
1148 t909 = t908 * t215;
1149 t912 = t395 * t89;
1150 t916 = t188 * t222;
1151 t920 = t15 * t564;
1152 t921 = t188 * t239;
1153 t925 = t395 * t122;
1154 t929 = t188 * t258;
1155 t940 = 0.89955555555555555554e0 * t20 * t244 + 0.10872222222222222222e1 * t31 * sigma[0] * t250 + 0.19248716727533128478e-3 * t32 * t255;
1156 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1157 v3rho2sigma[0] = t890 / 0.9e1 + t894 - 0.14e2 / 0.225e3 * t896 + t899 / 0.15e2 + t179 * t901 * t902 / 0.45e2 - 0.14e2 / 0.225e3 * t907 * t909 + t179 * t378 * t912 / 0.15e2 + t179 * t378 * t916 / 0.30e2 + 0.203e3 / 0.5625e4 * t179 * t920 * t921 - 0.14e2 / 0.375e3 * t179 * t383 * t925 - 0.7e1 / 0.375e3 * t179 * t383 * t929 + t179 * t180 * t940 / 0.50e2;
1158
1159 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1160 v3rho2sigma[1] = 0.0e0;
1161
1162 t945 = t126 * t154 * t402;
1163 t949 = t147 * t317 * t199 / 0.45e2;
1164 t950 = t262 * t153;
1165 t951 = t199 * t264;
1166 t955 = t199 * t272;
1167 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1168 v3rho2sigma[2] = t945 / 0.9e1 + t949 + t179 * t950 * t951 / 0.45e2 + t179 * t401 * t955 / 0.30e2;
1169
1170 t963 = t17 * t2 * t5;
1171 t964 = t963 * t407;
1172 t966 = t5 * t205;
1173 t967 = t178 * t966;
1174 t968 = t104 * t188;
1175 t969 = t968 * t282;
1176 t972 = t908 * t289;
1177 t975 = t395 * t139;
1178 t976 = t378 * t975;
1179 t979 = t227 * t2;
1180 t980 = t979 * t5;
1181 t982 = t378 * t188 * t7;
1182 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1183 v3rho2sigma[3] = t890 / 0.18e2 + t894 - 0.7e1 / 0.225e3 * t896 + t899 / 0.30e2 + t964 / 0.18e2 + t967 * t969 / 0.45e2 - 0.7e1 / 0.225e3 * t907 * t972 + t179 * t976 / 0.30e2 + t980 * t982 / 0.30e2;
1184
1185 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1186 v3rho2sigma[4] = 0.0e0;
1187
1188 t986 = t963 * t412;
1189 t988 = t5 * t262;
1190 t989 = t178 * t988;
1191 t990 = t153 * t199;
1192 t991 = t990 * t301;
1193 t995 = t401 * t199 * t7;
1194 t998 = t963 * t417;
1195 t1000 = t5 * t52;
1196 t1001 = t178 * t1000;
1197 t1002 = t348 * t199;
1198 t1003 = t1002 * t313;
1199 t1006 = t963 * t428;
1200 t1008 = t427 * t128;
1201 t1009 = t401 * t1008;
1202 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1203 v3rho2sigma[5] = t945 / 0.18e2 + t949 + t986 / 0.18e2 + t989 * t991 / 0.45e2 - t980 * t995 / 0.30e2 - 0.7e1 / 0.225e3 * t998 - 0.7e1 / 0.225e3 * t1001 * t1003 + t1006 / 0.30e2 + t179 * t1009 / 0.30e2;
1204
1205 t1014 = t901 * t188 * t324;
1206 t1018 = t378 * t188 * t329;
1207 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1208 v3rho2sigma[6] = t964 / 0.9e1 + t894 + t179 * t1014 / 0.45e2 + t179 * t1018 / 0.30e2;
1209
1210 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1211 v3rho2sigma[7] = 0.0e0;
1212
1213 t1025 = t950 * t199 * t333;
1214 t1028 = t1002 * t338;
1215 t1031 = t427 * t143;
1216 t1032 = t401 * t1031;
1217 t1036 = t401 * t199 * t342;
1218 t1039 = t53 * t872;
1219 t1041 = t1039 * t199 * t350;
1220 t1044 = t427 * t171;
1221 t1045 = t415 * t1044;
1222 t1049 = t415 * t199 * t369;
1223 t1059 = 0.89955555555555555554e0 * t20 * t355 + 0.10872222222222222222e1 * t31 * sigma[2] * t361 + 0.19248716727533128478e-3 * t63 * t366;
1224 t1060 = t191 * t1059;
1225 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1226 v3rho2sigma[8] = t986 / 0.9e1 + t949 - 0.14e2 / 0.225e3 * t998 + t1006 / 0.15e2 + t179 * t1025 / 0.45e2 - 0.14e2 / 0.225e3 * t1001 * t1028 + t179 * t1032 / 0.15e2 + t179 * t1036 / 0.30e2 + 0.203e3 / 0.5625e4 * t179 * t1041 - 0.14e2 / 0.375e3 * t179 * t1045 - 0.7e1 / 0.375e3 * t179 * t1049 + t179 * t1060 / 0.50e2;
1227
1228 t1064 = 0.7e1 / 0.225e3 * t963 * t432;
1229 t1065 = t14 * t237;
1230 t1067 = t1065 * t431 * t89;
1231 t1071 = t920 * t431 * t122;
1232 t1075 = t383 * t188 * t395;
1233 t1079 = t963 * t440 / 0.30e2;
1234 t1080 = t439 * t89;
1235 t1081 = t378 * t1080;
1236 t1084 = t439 * t122;
1237 t1085 = t383 * t1084;
1238 t1092 = -0.17166666666666666667e0 * t31 * t114 - 0.42774926061184729952e-4 * sigma[0] * t119;
1239 t1093 = t180 * t1092;
1240 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1241 v3rhosigma2[0] = -t1064 - 0.7e1 / 0.225e3 * t179 * t1067 + 0.203e3 / 0.5625e4 * t179 * t1071 - 0.14e2 / 0.375e3 * t179 * t1075 + t1079 + t179 * t1081 / 0.30e2 - 0.7e1 / 0.375e3 * t179 * t1085 + t179 * t1093 / 0.50e2;
1242
1243 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1244 v3rhosigma2[1] = 0.0e0;
1245
1246 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1247 v3rhosigma2[2] = 0.0e0;
1248
1249 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1250 v3rhosigma2[3] = 0.0e0;
1251
1252 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1253 v3rhosigma2[4] = 0.0e0;
1254
1255 t1097 = 0.7e1 / 0.225e3 * t963 * t444;
1256 t1098 = t52 * t348;
1257 t1100 = t1098 * t443 * t128;
1258 t1104 = t963 * t452 / 0.30e2;
1259 t1105 = t451 * t128;
1260 t1106 = t401 * t1105;
1261 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1262 v3rhosigma2[5] = -t1097 - 0.7e1 / 0.225e3 * t179 * t1100 + t1104 + t179 * t1106 / 0.30e2;
1263
1264 t1110 = t1065 * t431 * t139;
1265 t1113 = t439 * t139;
1266 t1114 = t378 * t1113;
1267 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1268 v3rhosigma2[6] = -t1064 - 0.7e1 / 0.225e3 * t179 * t1110 + t1079 + t179 * t1114 / 0.30e2;
1269
1270 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1271 v3rhosigma2[7] = 0.0e0;
1272
1273 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1274 v3rhosigma2[8] = 0.0e0;
1275
1276 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1277 v3rhosigma2[9] = 0.0e0;
1278
1279 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1280 v3rhosigma2[10] = 0.0e0;
1281
1282 t1118 = t1098 * t443 * t143;
1283 t1122 = t1039 * t443 * t171;
1284 t1126 = t415 * t199 * t427;
1285 t1129 = t451 * t143;
1286 t1130 = t401 * t1129;
1287 t1133 = t451 * t171;
1288 t1134 = t415 * t1133;
1289 t1141 = -0.17166666666666666667e0 * t31 * t163 - 0.42774926061184729952e-4 * sigma[2] * t168;
1290 t1142 = t191 * t1141;
1291 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1292 v3rhosigma2[11] = -t1097 - 0.7e1 / 0.225e3 * t179 * t1118 + 0.203e3 / 0.5625e4 * t179 * t1122 - 0.14e2 / 0.375e3 * t179 * t1126 + t1104 + t179 * t1130 / 0.30e2 - 0.7e1 / 0.375e3 * t179 * t1134 + t179 * t1142 / 0.50e2;
1293
1294 t1145 = t431 * t188;
1295 t1146 = t920 * t1145;
1296 t1150 = t383 * t188 * t439;
1297 t1153 = t180 * t42;
1298 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1299 v3sigma3[0] = 0.203e3 / 0.5625e4 * t179 * t1146 - 0.7e1 / 0.125e3 * t179 * t1150 + 0.10693731515296182488e-6 * t179 * t1153;
1300
1301 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1302 v3sigma3[1] = 0.0e0;
1303
1304 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1305 v3sigma3[2] = 0.0e0;
1306
1307 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1308 v3sigma3[3] = 0.0e0;
1309
1310 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1311 v3sigma3[4] = 0.0e0;
1312
1313 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1314 v3sigma3[5] = 0.0e0;
1315
1316 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1317 v3sigma3[6] = 0.0e0;
1318
1319 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1320 v3sigma3[7] = 0.0e0;
1321
1322 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1323 v3sigma3[8] = 0.0e0;
1324
1325 t1156 = t443 * t199;
1326 t1157 = t1039 * t1156;
1327 t1161 = t415 * t199 * t451;
1328 t1164 = t191 * t73;
1329 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1330 v3sigma3[9] = 0.203e3 / 0.5625e4 * t179 * t1157 - 0.7e1 / 0.125e3 * t179 * t1161 + 0.10693731515296182488e-6 * t179 * t1164;
1331
1332 #ifndef XC_DONT_COMPILE_LXC
1333
1334 if(order < 4) return;
1335
1336
1337 t1170 = 0.28e2 / 0.135e3 * t6 * t53 * t596 * t77;
1338 t1174 = 0.28e2 / 0.135e3 * t6 * t15 * t596 * t46;
1339 t1176 = t83 * t232 * t498;
1340 t1187 = t83 * t513 * t215;
1341 t1190 = t83 * t232 * t473;
1342 t1193 = t98 * t558 * t477;
1343 t1204 = t263 * t268 * t481;
1344 t1210 = -t1170 - t1174 + 0.4e1 / 0.15e2 * t1176 + 0.4e1 / 0.3e1 * t206 * t84 * t508 * t89 + 0.2e1 / 0.15e2 * t83 * t105 * t508 * t122 - 0.4e1 / 0.45e2 * t1187 + 0.4e1 / 0.15e2 * t1190 - 0.56e2 / 0.375e3 * t1193 + 0.2e1 / 0.15e2 * t83 * t105 * t534 * t89 - 0.28e2 / 0.375e3 * t98 * t238 * t534 * t122 + 0.8e1 / 0.3e1 * t1204 + 0.4e1 / 0.3e1 * t263 * t127 * t540 * t128;
1345 t1220 = t206 * t232 * t488;
1346 t1243 = t83 * t558 * t484;
1347 t1254 = t206 * t211 * t492;
1348 t1257 = t98 * t513 * t258;
1349 t1259 = -0.4e1 / 0.135e3 * t546 * t105 * t547 * t122 - 0.2e1 / 0.3e1 * t546 * t84 * t207 * t222 + 0.8e1 / 0.45e2 * t1220 + 0.2e1 / 0.15e2 * t206 * t105 * t207 * t258 + t83 * t105 * t222 * t258 / 0.5e1 + 0.406e3 / 0.1875e4 * t98 * t565 * t239 * t258 - 0.2e1 / 0.3e1 * t572 * t127 * t264 * t272 - 0.28e2 / 0.225e3 * t206 * t238 * t207 * t239 - 0.56e2 / 0.225e3 * t1243 + 0.812e3 / 0.3375e4 * t83 * t565 * t89 * t566 - 0.14e2 / 0.75e2 * t83 * t238 * t222 * t239 + 0.8e1 / 0.3e1 * t1254 - 0.2e1 / 0.75e2 * t1257;
1350 t1262 = t98 * t232 * t534;
1351 t1283 = t126 * t268 * t540;
1352 t1285 = t504 * t8;
1353 t1287 = t7 / t1285;
1354 t1289 = 0.12e2 * t505 - 0.12e2 * t1287;
1355 t1293 = t308 * t128;
1356 t1294 = t126 * t1293;
1357 t1297 = t83 * t554 * t222;
1358 t1301 = t98 * t227 * t237 * t239;
1359 t1304 = t126 * t580 * t272;
1360 t1306 = t12 * t12;
1361 t1309 = t6 / t13 / t1306;
1362 t1310 = t207 * t207;
1363 t1315 = t546 * t211 * t547;
1364 t1318 = t206 * t554 * t207;
1365 t1320 = t222 * t222;
1366 t1323 = 0.4e1 / 0.75e2 * t1262 + t98 * t105 * (0.23788246913580246913e2 * t20 * sigma[0] / t23 / t112 + 0.33220679012345679012e2 * t31 * t32 / t22 / t118 + 0.70578628000954804422e-2 * t40 / t41 / t33) / 0.50e2 + 0.4e1 / 0.3e1 * t1283 + t126 * t127 * t1289 / 0.2e1 + 0.16e2 / 0.27e2 * t1294 - 0.2e1 / 0.3e1 * t1297 + 0.28e2 / 0.1125e4 * t1301 - 0.2e1 / 0.3e1 * t1304 + 0.4e1 / 0.27e2 * t1309 * t84 * t1310 - 0.8e1 / 0.27e2 * t1315 - 0.4e1 / 0.9e1 * t1318 + t206 * t84 * t1320;
1367 t1326 = t98 * t93 * t564 * t566;
1368 t1330 = 0.1e1 / t103 / t562 / t45;
1369 t1332 = t239 * t239;
1370 t1336 = t258 * t258;
1371 t1340 = t50 * t50;
1372 t1343 = t6 / t51 / t1340;
1373 t1344 = t264 * t264;
1374 t1349 = t572 * t268 * t573;
1375 t1352 = t263 * t580 * t264;
1376 t1354 = t272 * t272;
1377 t1358 = t83 * t211 * t508;
1378 t1364 = t295 * t89;
1379 t1365 = t83 * t1364;
1380 t1367 = t671 * t122;
1381 t1368 = t98 * t1367;
1382 t1371 = t6 * t14 * t17;
1383 t1380 = 0.1624e4 / 0.16875e5 * t1326 - 0.8932e4 / 0.84375e5 * t98 * t17 * t1330 * t1332 - 0.7e1 / 0.125e3 * t98 * t238 * t1336 + 0.4e1 / 0.27e2 * t1343 * t127 * t1344 - 0.8e1 / 0.27e2 * t1349 - 0.4e1 / 0.9e1 * t1352 + t263 * t127 * t1354 + 0.4e1 / 0.3e1 * t1358 - t83 * t84 * t1289 / 0.2e1 + 0.16e2 / 0.27e2 * t1365 + 0.16e2 / 0.675e3 * t1368 - 0.28e2 / 0.75e2 * t1371 * t237 * t89 * t477 + 0.4e1 / 0.15e2 * t624 * t104 * t89 * t498;
1384 t1384 = 0.16e2 / 0.45e2 * t467;
1385 t1385 = 0.16e2 / 0.45e2 * t471;
1386 t1396 = t8 * (t1210 + t1259 + t1323 + t1380) + t1384 + t1385 + 0.2e1 / 0.5e1 * t475 - 0.28e2 / 0.125e3 * t479 + 0.4e1 * t483 - 0.28e2 / 0.75e2 * t486 + 0.4e1 / 0.15e2 * t490 + 0.4e1 * t494 + 0.8e1 / 0.15e2 * t496 + 0.2e1 / 0.5e1 * t500 + 0.4e1 * t503 + 0.2e1 * t510;
1387 t1410 = -0.4e1 / 0.75e2 * t515 + 0.4e1 / 0.25e2 * t518 + 0.2e1 / 0.25e2 * t536 + 0.4e1 * t539 + 0.2e1 * t542 - 0.4e1 / 0.9e1 * t549 + 0.8e1 / 0.3e1 * t552 - 0.4e1 / 0.3e1 * t556 - 0.56e2 / 0.375e3 * t560 + 0.812e3 / 0.5625e4 * t568 - 0.4e1 / 0.9e1 * t575 + 0.8e1 / 0.3e1 * t578 - 0.4e1 / 0.3e1 * t582;
1388 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1389 v4rho4[0] = t1396 + t1410;
1390
1391 t1413 = t6 * t14 * t596 * t46;
1392 t1417 = t6 * t52 * t596 * t77;
1393 t1424 = t546 * t211 * t656;
1394 t1431 = t206 * t554 * t282;
1395 t1432 = 0.2e1 / 0.9e1 * t1431;
1396 t1434 = t206 * t604 * t663;
1397 t1437 = t83 * t558 * t667;
1398 t1450 = t83 * t596 * t104 * t672;
1399 t1457 = t572 * t268 * t676;
1400 t1460 = t263 * t580 * t301;
1401 t1461 = 0.2e1 / 0.9e1 * t1460;
1402 t1467 = t263 * t317 * t683;
1403 t1469 = -0.2e1 * t1413 + 0.2e1 * t1417 - t1170 - t1174 + 0.4e1 / 0.27e2 * t1309 * t84 * t139 * t547 - 0.2e1 / 0.9e1 * t1424 - t546 * t295 * t7 * t207 / 0.3e1 - t1432 - 0.5e1 / 0.3e1 * t1434 - 0.14e2 / 0.225e3 * t1437 + 0.203e3 / 0.3375e4 * t83 * t565 * t139 * t566 - 0.7e1 / 0.75e2 * t83 * t294 * t237 * t7 * t239 - t1450 / 0.6e1 + 0.4e1 / 0.27e2 * t1343 * t127 * t143 * t573 - 0.2e1 / 0.9e1 * t1457 - t1461 - t572 * t154 * t171 * t573 / 0.135e3 + 0.2e1 / 0.45e2 * t1467;
1404 t1471 = t126 * t611 * t313;
1405 t1472 = t1471 / 0.45e2;
1406 t1474 = t206 * t211 * t632;
1407 t1488 = t263 * t268 * t643;
1408 t1491 = t263 * t597 * t647;
1409 t1494 = t126 * t317 * t651;
1410 t1501 = t83 * t513 * t289;
1411 t1502 = t1501 / 0.45e2;
1412 t1504 = t83 * t232 * t639;
1413 t1525 = -t1472 + 0.2e1 / 0.3e1 * t1474 + t206 * t295 * t7 * t222 + t83 * t671 * t7 * t258 / 0.10e2 + t572 * t308 * t7 * t264 / 0.3e1 + 0.2e1 / 0.3e1 * t1488 + 0.5e1 / 0.3e1 * t1491 + t1494 / 0.15e2 + t206 * t84 * t139 * t508 / 0.3e1 - t1502 + t1504 / 0.15e2 + t83 * t105 * t139 * t534 / 0.30e2 + t263 * t127 * t143 * t540 / 0.3e1 - t263 * t308 * t7 * t272 + t126 * t154 * t171 * t540 / 0.30e2 + t1176 / 0.15e2 - 0.2e1 / 0.45e2 * t1187 + t1190 / 0.15e2 - 0.14e2 / 0.375e3 * t1193;
1414 t1539 = 0.1e1 / t16 / t504;
1415 t1542 = t126 * t1539 * t77 * t7;
1416 t1544 = t295 * t139;
1417 t1545 = t83 * t1544;
1418 t1549 = t83 * t1539 * t46 * t7;
1419 t1551 = t308 * t143;
1420 t1552 = t126 * t1551;
1421 t1554 = t784 * t171;
1422 t1555 = t147 * t1554;
1423 t1560 = 0.2e1 / 0.3e1 * t1204 + 0.2e1 / 0.45e2 * t1220 - 0.14e2 / 0.225e3 * t1243 + 0.2e1 / 0.3e1 * t1254 - t1257 / 0.75e2 + t1262 / 0.75e2 + t1283 / 0.3e1 + 0.4e1 / 0.9e1 * t1294 + t206 * t1364 + t83 * t1367 / 0.10e2 - 0.8e1 / 0.3e1 * t1542 + 0.4e1 / 0.27e2 * t1545 + 0.8e1 / 0.3e1 * t1549 + 0.4e1 / 0.27e2 * t1552 + 0.4e1 / 0.675e3 * t1555 - t263 * t1293 - t1297 / 0.3e1 + 0.14e2 / 0.1125e4 * t1301;
1424 t1572 = t6 * t205 * t93 * t626;
1425 t1580 = t237 * t139;
1426 t1585 = t6 * t571 * t17;
1427 t1586 = t77 * t143;
1428 t1598 = t6 * t205 * t294;
1429 t1604 = t6 * t545 * t17;
1430 t1608 = t46 * t139;
1431 t1612 = -t1304 / 0.3e1 - 0.2e1 / 0.27e2 * t1315 - 0.2e1 / 0.9e1 * t1318 + 0.406e3 / 0.16875e5 * t1326 - 0.2e1 / 0.27e2 * t1349 - 0.2e1 / 0.9e1 * t1352 + t1358 / 0.3e1 + 0.4e1 / 0.9e1 * t1365 + 0.4e1 / 0.225e3 * t1368 + 0.4e1 / 0.45e2 * t1572 + t624 * t625 * t498 / 0.15e2 + t624 * t625 * t473 / 0.15e2 - 0.7e1 / 0.75e2 * t1371 * t1580 * t477 - t1585 * t1586 * t481 / 0.3e1 + t748 * t153 * t171 * t481 / 0.15e2 - 0.14e2 / 0.225e3 * t624 * t1580 * t484 + 0.2e1 / 0.15e2 * t1598 * t104 * t7 * t215 - t1604 * t625 * t488 / 0.45e2 - t1604 * t1608 * t492 / 0.3e1;
1432 t1621 = 0.2e1 / 0.15e2 * t630;
1433 t1622 = t8 * (t1469 + t1525 + t1560 + t1612) + 0.5e1 / 0.2e1 * t599 - t882 - 0.5e1 / 0.2e1 * t606 - t883 - t884 + 0.3e1 / 0.2e1 * t617 - 0.3e1 / 0.2e1 * t621 + t1384 + t1385 + 0.2e1 / 0.15e2 * t627 + t1621;
1434 t1623 = 0.2e1 / 0.15e2 * t637;
1435 t1628 = 0.4e1 / 0.3e1 * t661;
1436 t1630 = t634 + t1623 + t641 / 0.10e2 + t645 - 0.2e1 * t649 + t653 / 0.10e2 + t476 - t480 + t483 - t658 / 0.3e1 + t1628 + 0.2e1 * t665;
1437 t1635 = 0.4e1 / 0.3e1 * t681;
1438 t1639 = -0.7e1 / 0.75e2 * t669 + t674 / 0.5e1 - t678 / 0.3e1 + t1635 + t685 / 0.15e2 - t487 + t491 + t494 + 0.4e1 / 0.15e2 * t496 + t501 + 0.2e1 * t503 + t511;
1439 t1646 = -t515 / 0.25e2 + 0.2e1 / 0.25e2 * t518 + t537 + 0.2e1 * t539 + t543 - t550 + 0.4e1 / 0.3e1 * t552 - t556 - 0.28e2 / 0.375e3 * t560 + t569 - t576 + 0.4e1 / 0.3e1 * t578 - t582;
1440 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1441 v4rho4[1] = t1622 + t1630 + t1639 + t1646;
1442
1443 t1651 = t206 * t604 * t753;
1444 t1654 = t126 * t611 * t338;
1445 t1661 = t546 * t211 * t760;
1446 t1672 = t206 * t211 * t768;
1447 t1674 = 0.2e1 / 0.3e1 * t1413 - 0.2e1 / 0.3e1 * t1417 - t1170 - t1174 - 0.10e2 / 0.9e1 * t1651 - 0.2e1 / 0.135e3 * t1654 + 0.4e1 / 0.27e2 * t1309 * t84 * t324 * t207 - 0.4e1 / 0.27e2 * t1661 - 0.14e2 / 0.675e3 * t206 * t238 * t324 * t239 - t546 * t84 * t329 * t207 / 0.9e1 + 0.4e1 / 0.9e1 * t1672;
1448 t1684 = t572 * t268 * t776;
1449 t1691 = t263 * t268 * t789;
1450 t1698 = t126 * t734 * t705;
1451 t1705 = t126 * t317 * t709;
1452 t1715 = -0.7e1 / 0.225e3 * t83 * t238 * t329 * t239 + 0.4e1 / 0.27e2 * t1343 * t127 * t333 * t264 - 0.4e1 / 0.27e2 * t1684 - t572 * t127 * t342 * t264 / 0.9e1 + 0.4e1 / 0.9e1 * t1691 - 0.14e2 / 0.675e3 * t263 * t349 * t350 * t264 - 0.28e2 / 0.675e3 * t1698 + t263 * t154 * t369 * t264 / 0.45e2 + 0.2e1 / 0.45e2 * t1705 + t83 * t105 * t714 * t122 / 0.15e2 + 0.2e1 / 0.3e1 * t263 * t127 * t718 * t128;
1453 t1722 = t206 * t232 * t764;
1454 t1733 = t83 * t232 * t772;
1455 t1744 = t263 * t597 * t780;
1456 t1758 = -t546 * t84 * t324 * t222 / 0.9e1 + 0.4e1 / 0.135e3 * t1722 + t206 * t105 * t324 * t258 / 0.45e2 + t206 * t84 * t329 * t222 / 0.3e1 + 0.2e1 / 0.45e2 * t1733 + t83 * t105 * t329 * t258 / 0.30e2 - t572 * t127 * t333 * t272 / 0.9e1 + 0.10e2 / 0.9e1 * t1744 + t263 * t127 * t342 * t272 / 0.3e1 - 0.7e1 / 0.225e3 * t126 * t349 * t350 * t272 + t126 * t154 * t369 * t272 / 0.30e2;
1457 t1761 = t126 * t596 * t153 * t785;
1458 t1776 = t1761 / 0.9e1 + 0.2e1 / 0.3e1 * t206 * t84 * t714 * t89 - 0.4e1 / 0.27e2 * t1424 - 0.8e1 / 0.27e2 * t1431 + 0.8e1 / 0.9e1 * t1434 - 0.28e2 / 0.675e3 * t1437 + 0.4e1 / 0.45e2 * t1450 - 0.4e1 / 0.27e2 * t1457 - 0.8e1 / 0.27e2 * t1460 + 0.4e1 / 0.135e3 * t1467 - 0.4e1 / 0.135e3 * t1471;
1459 t1787 = 0.1e1 / t16 / t1285;
1460 t1789 = t7 * t7;
1461 t1796 = t83 * t211 * t714;
1462 t1799 = t126 * t268 * t718;
1463 t1801 = 0.4e1 / 0.9e1 * t1474 + 0.4e1 / 0.9e1 * t1488 - 0.8e1 / 0.9e1 * t1491 + 0.2e1 / 0.45e2 * t1494 - 0.4e1 / 0.135e3 * t1501 + 0.2e1 / 0.45e2 * t1504 - 0.2e1 / 0.135e3 * t1187 + 0.2e1 / 0.3e1 * t263 * t1787 * t77 * t1789 + 0.2e1 / 0.3e1 * t206 * t1544 + 0.2e1 / 0.3e1 * t1796 + 0.2e1 / 0.3e1 * t1799;
1464 t1803 = t206 * t554 * t324;
1465 t1806 = t83 * t554 * t329;
1466 t1809 = t263 * t580 * t333;
1467 t1816 = t126 * t580 * t342;
1468 t1820 = t147 * t227 * t348 * t350;
1469 t1823 = t147 * t611 * t369;
1470 t1831 = -0.2e1 / 0.27e2 * t1803 - t1806 / 0.9e1 - 0.2e1 / 0.27e2 * t1809 - 0.2e1 / 0.3e1 * t263 * t1551 - t126 * t1554 / 0.15e2 - t1816 / 0.9e1 + 0.14e2 / 0.3375e4 * t1820 - t1823 / 0.225e3 + 0.2e1 / 0.3e1 * t206 * t1787 * t46 * t1789 - t1257 / 0.225e3 + 0.8e1 / 0.27e2 * t1294;
1471 t1844 = -0.32e2 / 0.9e1 * t1542 + 0.8e1 / 0.27e2 * t1545 + 0.32e2 / 0.9e1 * t1549 + 0.8e1 / 0.27e2 * t1552 + 0.8e1 / 0.675e3 * t1555 - t1297 / 0.9e1 + 0.14e2 / 0.3375e4 * t1301 - t1304 / 0.9e1 - 0.2e1 / 0.27e2 * t1318 - 0.2e1 / 0.27e2 * t1352 + 0.8e1 / 0.27e2 * t1365;
1472 t1858 = t6 * t262 * t294;
1473 t1881 = t6 * t262 * t93 * t750;
1474 t1884 = 0.8e1 / 0.675e3 * t1368 - 0.4e1 / 0.9e1 * t6 * t545 * t294 * t1608 * t663 + 0.4e1 / 0.45e2 * t1598 * t625 * t672 + 0.2e1 / 0.45e2 * t748 * t749 * t651 - 0.4e1 / 0.45e2 * t1858 * t153 * t7 * t313 - 0.2e1 / 0.135e3 * t1604 * t104 * t324 * t215 + 0.2e1 / 0.45e2 * t624 * t104 * t329 * t215 + 0.4e1 / 0.9e1 * t6 * t571 * t294 * t1586 * t647 - 0.2e1 / 0.135e3 * t1585 * t749 * t683 + 0.8e1 / 0.135e3 * t1881 + 0.8e1 / 0.135e3 * t1572;
1475 t1900 = t8 * (t1674 + t1715 + t1758 + t1776 + t1801 + t1831 + t1844 + t1884) - 0.14e2 / 0.225e3 * t707 + t711 / 0.15e2 + t716 + t720 + 0.4e1 / 0.9e1 * t723 + 0.2e1 / 0.3e1 * t726 + 0.4e1 / 0.9e1 * t729 + 0.2e1 / 0.3e1 * t732 - 0.28e2 / 0.1125e4 * t736 + 0.2e1 / 0.75e2 * t739 + t599 / 0.3e1 - 0.2e1 / 0.3e1 * t602 - t606 / 0.3e1;
1476 t1911 = -0.2e1 / 0.3e1 * t609 - 0.2e1 / 0.75e2 * t613 + t617 - t621 + t1384 + t1385 + 0.4e1 / 0.45e2 * t751 + 0.4e1 / 0.45e2 * t627 + 0.4e1 / 0.3e1 * t755 + 0.4e1 / 0.45e2 * t758 - 0.2e1 / 0.9e1 * t762 + 0.2e1 / 0.45e2 * t766 + 0.2e1 / 0.3e1 * t770 + t774 / 0.15e2;
1477 t1927 = -0.2e1 / 0.9e1 * t778 - 0.4e1 / 0.3e1 * t782 - 0.2e1 / 0.15e2 * t787 + 0.2e1 / 0.3e1 * t791 + 0.8e1 / 0.45e2 * t630 + 0.2e1 / 0.3e1 * t634 + 0.8e1 / 0.45e2 * t637 + t641 / 0.15e2 + 0.2e1 / 0.3e1 * t645 - 0.4e1 / 0.3e1 * t649 + t653 / 0.15e2 - 0.2e1 / 0.9e1 * t658 + 0.16e2 / 0.9e1 * t661 + 0.4e1 / 0.3e1 * t665;
1478 t1943 = -0.14e2 / 0.225e3 * t669 + 0.2e1 / 0.15e2 * t674 - 0.2e1 / 0.9e1 * t678 + 0.16e2 / 0.9e1 * t681 + 0.2e1 / 0.45e2 * t685 + 0.4e1 / 0.45e2 * t496 + 0.2e1 / 0.3e1 * t503 - 0.2e1 / 0.75e2 * t515 + 0.2e1 / 0.75e2 * t518 + 0.2e1 / 0.3e1 * t539 + 0.4e1 / 0.9e1 * t552 - 0.2e1 / 0.3e1 * t556 - 0.28e2 / 0.1125e4 * t560 + 0.4e1 / 0.9e1 * t578 - 0.2e1 / 0.3e1 * t582;
1479 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1480 v4rho4[2] = t1900 + t1911 + t1927 + t1943;
1481
1482 t1969 = t6 * t52 * t17;
1483 t1975 = t126 * t734 * t808;
1484 t1982 = t206 * t211 * t812;
1485 t1988 = t263 * t317 * t815;
1486 t1991 = t263 * t268 * t819;
1487 t1994 = t126 * t317 * t822;
1488 t1997 = t126 * t317 * t826;
1489 t2000 = t147 * t734 * t830;
1490 t2002 = -t1604 * t1608 * t768 / 0.3e1 + t624 * t625 * t772 / 0.15e2 - t1585 * t153 * t333 * t313 / 0.45e2 - 0.2e1 / 0.15e2 * t1858 * t749 * t785 - t1585 * t1586 * t789 / 0.3e1 + t748 * t153 * t342 * t313 / 0.15e2 + t748 * t749 * t709 / 0.15e2 - 0.7e1 / 0.75e2 * t1969 * t348 * t171 * t709 - t1170 - t1174 - 0.14e2 / 0.225e3 * t1975 - t546 * t295 * t324 * t7 / 0.3e1 + 0.2e1 / 0.3e1 * t1982 + t206 * t295 * t7 * t329 + 0.2e1 / 0.45e2 * t1988 + 0.2e1 / 0.3e1 * t1991 + t1994 / 0.15e2 + t1997 / 0.15e2 - 0.14e2 / 0.375e3 * t2000;
1491 t2065 = 0.203e3 / 0.3375e4 * t126 * t873 * t874 * t128 + 0.4e1 / 0.27e2 * t1309 * t84 * t861 * t89 - t546 * t105 * t861 * t122 / 0.135e3 + t206 * t84 * t139 * t714 + 0.4e1 / 0.27e2 * t1343 * t127 * t865 * t128 + t572 * t308 * t333 * t7 / 0.3e1 + t263 * t127 * t143 * t718 - t263 * t308 * t7 * t342 + t126 * t154 * t718 * t171 / 0.10e2 - t126 * t784 * t7 * t369 / 0.10e2 + 0.7e1 / 0.75e2 * t126 * t294 * t348 * t7 * t350 + t206 * t84 * t835 * t89 / 0.3e1 + t83 * t105 * t835 * t122 / 0.30e2 + t263 * t127 * t839 * t128 / 0.3e1 + t126 * t154 * t857 * t128 / 0.30e2 + 0.4e1 / 0.3e1 * t1651 - 0.2e1 / 0.45e2 * t1654 - 0.2e1 / 0.9e1 * t1661 + 0.2e1 / 0.3e1 * t1672;
1492 t2077 = t147 * t93 * t872 * t874;
1493 t2080 = t546 * t211 * t861;
1494 t2083 = t572 * t268 * t865;
1495 t2086 = t83 * t211 * t835;
1496 t2089 = t126 * t268 * t839;
1497 t2092 = t147 * t317 * t857;
1498 t2096 = 0.6e1 * t505 + 0.12e2 * t1287;
1499 t2100 = -0.2e1 / 0.9e1 * t1684 + 0.2e1 / 0.3e1 * t1691 - 0.14e2 / 0.225e3 * t1698 + t1705 / 0.15e2 + 0.2e1 / 0.45e2 * t1722 + t1733 / 0.15e2 - 0.4e1 / 0.3e1 * t1744 - 0.2e1 / 0.15e2 * t1761 - t1432 - t1461 - t1472 - t1502 + 0.406e3 / 0.16875e5 * t2077 - 0.2e1 / 0.27e2 * t2080 - 0.2e1 / 0.27e2 * t2083 + t2086 / 0.3e1 + t2089 / 0.3e1 + t2092 / 0.75e2 + t83 * t84 * t2096 / 0.2e1;
1500 t2119 = t348 * t143;
1501 t2124 = -t126 * t127 * t2096 / 0.2e1 + t1796 + t1799 - 0.2e1 / 0.9e1 * t1803 - t1806 / 0.3e1 - 0.2e1 / 0.9e1 * t1809 - t1816 / 0.3e1 + 0.14e2 / 0.1125e4 * t1820 - t1823 / 0.75e2 + 0.4e1 / 0.27e2 * t1294 + t1542 / 0.3e1 + 0.4e1 / 0.9e1 * t1545 - t1549 / 0.3e1 + 0.4e1 / 0.9e1 * t1552 + 0.4e1 / 0.225e3 * t1555 + 0.4e1 / 0.27e2 * t1365 + 0.4e1 / 0.675e3 * t1368 - 0.14e2 / 0.225e3 * t748 * t2119 * t705 + 0.4e1 / 0.45e2 * t1881;
1502 t2130 = t8 * (t2002 + t2065 + t2100 + t2124) - 0.7e1 / 0.75e2 * t707 + t711 / 0.10e2 - t811 + t814 + t818 + t821 + t825 + t829 - t833 + t838 + t842;
1503 t2139 = t860 - t864 - t868 + t877 + 0.3e1 / 0.2e1 * t716 + 0.3e1 / 0.2e1 * t720 + 0.4e1 / 0.3e1 * t723 + 0.2e1 * t726 + 0.4e1 / 0.3e1 * t729 + 0.2e1 * t732 - 0.28e2 / 0.375e3 * t736 + 0.2e1 / 0.25e2 * t739;
1504 t2149 = -0.2e1 * t599 - t602 + 0.2e1 * t606 - t609 - t613 / 0.25e2 + t1384 + t1385 + 0.2e1 / 0.15e2 * t751 + 0.2e1 * t755 + 0.4e1 / 0.15e2 * t758 - t762 / 0.3e1 + t766 / 0.15e2;
1505 t2154 = t770 + t774 / 0.10e2 - t778 / 0.3e1 - 0.2e1 * t782 - t787 / 0.5e1 + t791 + t1621 + t1623 + t1628 + t1635 - t516 - t557 - t583;
1506 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1507 v4rho4[3] = t2130 + t2139 + t2149 + t2154;
1508
1509 t2194 = 0.4e1 / 0.15e2 * t748 * t749 * t822 - 0.28e2 / 0.75e2 * t1969 * t2119 * t830 - t1170 - t1174 + 0.406e3 / 0.1875e4 * t147 * t873 * t350 * t369 - 0.2e1 / 0.3e1 * t546 * t84 * t324 * t329 + 0.2e1 / 0.15e2 * t126 * t154 * t839 * t171 + 0.2e1 / 0.15e2 * t126 * t154 * t857 * t143 - 0.28e2 / 0.375e3 * t147 * t349 * t857 * t171 - 0.14e2 / 0.75e2 * t126 * t349 * t342 * t350 + 0.4e1 / 0.3e1 * t206 * t84 * t835 * t139 + 0.4e1 / 0.3e1 * t263 * t127 * t839 * t143;
1510 t2226 = -0.28e2 / 0.225e3 * t263 * t349 * t333 * t350 + 0.812e3 / 0.3375e4 * t126 * t873 * t143 * t874 - 0.56e2 / 0.225e3 * t1975 + 0.8e1 / 0.3e1 * t1982 + 0.8e1 / 0.45e2 * t1988 + 0.8e1 / 0.3e1 * t1991 + 0.4e1 / 0.15e2 * t1994 + 0.4e1 / 0.15e2 * t1997 - 0.56e2 / 0.375e3 * t2000 - 0.4e1 / 0.135e3 * t572 * t154 * t865 * t171 - 0.2e1 / 0.3e1 * t572 * t127 * t333 * t342 + 0.2e1 / 0.15e2 * t263 * t154 * t333 * t369 + t126 * t154 * t342 * t369 / 0.5e1;
1511 t2229 = t342 * t342;
1512 t2234 = 0.1e1 / t152 / t870 / t76;
1513 t2236 = t350 * t350;
1514 t2240 = t369 * t369;
1515 t2244 = t324 * t324;
1516 t2248 = t329 * t329;
1517 t2251 = t333 * t333;
1518 t2256 = 0.12e2 * t505 + 0.12e2 * t1287;
1519 t2284 = -0.4e1 / 0.45e2 * t1654 + t263 * t127 * t2229 - 0.8932e4 / 0.84375e5 * t147 * t17 * t2234 * t2236 - 0.7e1 / 0.125e3 * t147 * t349 * t2240 + 0.4e1 / 0.27e2 * t1309 * t84 * t2244 + t206 * t84 * t2248 + 0.4e1 / 0.27e2 * t1343 * t127 * t2251 + t83 * t84 * t2256 / 0.2e1 - t126 * t127 * t2256 / 0.2e1 + t147 * t154 * (0.23788246913580246913e2 * t20 * sigma[2] / t57 / t161 + 0.33220679012345679012e2 * t31 * t63 / t56 / t167 + 0.70578628000954804422e-2 * t71 / t72 / t64) / 0.50e2 + 0.1624e4 / 0.16875e5 * t2077 - 0.8e1 / 0.27e2 * t2080;
1520 t2298 = -0.8e1 / 0.27e2 * t2083 + 0.4e1 / 0.3e1 * t2086 + 0.4e1 / 0.3e1 * t2089 + 0.4e1 / 0.75e2 * t2092 - 0.4e1 / 0.9e1 * t1803 - 0.2e1 / 0.3e1 * t1806 - 0.4e1 / 0.9e1 * t1809 - 0.2e1 / 0.3e1 * t1816 + 0.28e2 / 0.1125e4 * t1820 - 0.2e1 / 0.75e2 * t1823 + 0.16e2 / 0.27e2 * t1545 + 0.16e2 / 0.27e2 * t1552 + 0.16e2 / 0.675e3 * t1555;
1521 t2314 = t8 * (t2194 + t2226 + t2284 + t2298) - 0.28e2 / 0.75e2 * t810 + 0.4e1 * t814 + 0.4e1 / 0.15e2 * t817 + 0.4e1 * t821 + 0.2e1 / 0.5e1 * t824 + 0.2e1 / 0.5e1 * t828 - 0.28e2 / 0.125e3 * t832 + 0.2e1 * t837 + 0.2e1 * t841 + 0.2e1 / 0.25e2 * t859 - 0.4e1 / 0.9e1 * t863 - 0.4e1 / 0.9e1 * t867;
1522 t2326 = 0.812e3 / 0.5625e4 * t876 + 0.8e1 / 0.3e1 * t723 + 0.4e1 * t726 + 0.8e1 / 0.3e1 * t729 + 0.4e1 * t732 - 0.56e2 / 0.375e3 * t736 + 0.4e1 / 0.25e2 * t739 - 0.4e1 / 0.3e1 * t602 - 0.4e1 / 0.3e1 * t609 - 0.4e1 / 0.75e2 * t613 + t1384 + t1385 + 0.8e1 / 0.15e2 * t758;
1523 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1524 v4rho4[4] = t2314 + t2326;
1525
1526 t2328 = t83 * t105 * t916;
1527 t2331 = t98 * t558 * t384;
1528 t2334 = t98 * t238 * t925;
1529 t2337 = t98 * t238 * t929;
1530 t2352 = t206 * t105 * t902;
1531 t2355 = t83 * t232 * t379;
1532 t2358 = t98 * t565 * t921;
1533 t2360 = t545 * t104;
1534 t2365 = t15 * t1330;
1535 t2382 = t2328 / 0.6e1 - 0.14e2 / 0.225e3 * t2331 - 0.14e2 / 0.75e2 * t2334 - 0.7e1 / 0.75e2 * t2337 + t179 * t901 * t395 * t207 / 0.15e2 + t179 * t378 * t940 * t89 / 0.10e2 + t179 * t378 * t395 * t222 / 0.10e2 + t2352 / 0.9e1 + t2355 / 0.9e1 + 0.203e3 / 0.1125e4 * t2358 - t179 * t2360 * t188 * t547 / 0.135e3 - 0.8932e4 / 0.84375e5 * t179 * t2365 * t188 * t566 + t179 * t378 * t188 * t508 / 0.30e2 + 0.203e3 / 0.1875e4 * t179 * t920 * t395 * t239 - 0.7e1 / 0.125e3 * t179 * t383 * t940 * t122;
1536 t2392 = t83 * t105 * t912;
1537 t2394 = t1371 * t909;
1538 t2402 = t237 * t395;
1539 t2413 = t178 * t5 * t15;
1540 t2414 = t564 * t188;
1541 t2423 = t98 * t513 * t188 / 0.135e3;
1542 t2425 = t98 * t232 * t395;
1543 t2428 = t98 * t105 * t940;
1544 t2441 = -0.7e1 / 0.125e3 * t179 * t383 * t395 * t258 - 0.7e1 / 0.375e3 * t179 * t383 * t188 * t534 + t2392 / 0.3e1 - 0.14e2 / 0.45e2 * t2394 - 0.14e2 / 0.225e3 * t967 * t908 * t488 + t967 * t968 * t492 / 0.15e2 - 0.14e2 / 0.75e2 * t907 * t2402 * t215 - 0.7e1 / 0.75e2 * t907 * t908 * t498 - 0.7e1 / 0.75e2 * t907 * t908 * t473 + 0.203e3 / 0.1875e4 * t2413 * t2414 * t477 + 0.203e3 / 0.1125e4 * t907 * t2414 * t484 - t2423 + t2425 / 0.15e2 + t2428 / 0.10e2 + t179 * t180 * (-0.41979259259259259259e1 * t20 * t521 - 0.79729629629629629628e1 * t31 * sigma[0] * t526 - 0.19248716727533128478e-2 * t32 * t531) / 0.50e2;
1545 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1546 v4rho3sigma[0] = t2382 + t2441;
1547
1548 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1549 v4rho3sigma[1] = 0.0e0;
1550
1551 t2443 = t263 * t154 * t951;
1552 t2446 = t126 * t317 * t402;
1553 t2449 = t126 * t154 * t955;
1554 t2453 = t147 * t611 * t199 / 0.135e3;
1555 t2454 = t571 * t153;
1556 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1557 v4rho3sigma[2] = t2443 / 0.9e1 + t2446 / 0.9e1 + t2449 / 0.6e1 - t2453 - t179 * t2454 * t199 * t573 / 0.135e3 + t989 * t990 * t481 / 0.15e2 + t179 * t401 * t199 * t540 / 0.30e2;
1558
1559 t2467 = t93 * t2 * t5;
1560 t2468 = t2467 * t407;
1561 t2471 = t294 * t2 * t5;
1562 t2472 = t2471 * t982;
1563 t2486 = t963 * t976 / 0.9e1;
1564 t2493 = 0.2e1 / 0.27e2 * t624 * t969;
1565 t2494 = t2468 / 0.27e2 + t2472 / 0.90e2 + t980 * t378 * t395 * t7 / 0.15e2 + t2328 / 0.18e2 - 0.28e2 / 0.675e3 * t2331 - 0.14e2 / 0.225e3 * t2334 - 0.7e1 / 0.225e3 * t2337 + t2352 / 0.27e2 + 0.2e1 / 0.27e2 * t2355 + 0.203e3 / 0.3375e4 * t2358 + t2486 + t179 * t378 * t940 * t139 / 0.30e2 + t2392 / 0.9e1 + t2493;
1566 t2496 = 0.14e2 / 0.135e3 * t1371 * t972;
1567 t2511 = t178 * t5 * t545;
1568 t2515 = t979 * t966;
1569 t2536 = -t2496 + 0.2e1 / 0.45e2 * t967 * t104 * t395 * t282 + t967 * t968 * t632 / 0.45e2 - 0.14e2 / 0.225e3 * t907 * t2402 * t289 - 0.7e1 / 0.225e3 * t907 * t908 * t639 - t2511 * t968 * t656 / 0.135e3 + 0.2e1 / 0.45e2 * t2515 * t968 * t663 + 0.203e3 / 0.3375e4 * t907 * t2414 * t667 - 0.14e2 / 0.225e3 * t979 * t906 * t908 * t672 - 0.14e2 / 0.135e3 * t2394 - t2423 + 0.2e1 / 0.45e2 * t2425 + t2428 / 0.30e2 + t980 * t378 * t188 / 0.30e2 - 0.28e2 / 0.675e3 * t967 * t908 * t282 * t122;
1570 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1571 v4rho3sigma[3] = t2494 + t2536;
1572
1573 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1574 v4rho3sigma[4] = 0.0e0;
1575
1576 t2541 = t2471 * t995;
1577 t2546 = t2467 * t428;
1578 t2549 = t963 * t1009 / 0.9e1;
1579 t2553 = t2467 * t412;
1580 t2556 = 0.2e1 / 0.27e2 * t748 * t991;
1581 t2560 = t979 * t988;
1582 t2564 = t2467 * t417;
1583 t2567 = 0.14e2 / 0.135e3 * t1969 * t1003;
1584 t2576 = t178 * t5 * t571;
1585 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1586 v4rho3sigma[5] = t179 * t401 * t427 * t272 / 0.30e2 - t2541 / 0.90e2 - t980 * t401 * t199 / 0.30e2 + t2546 / 0.45e2 + t2549 + t2449 / 0.18e2 - t2453 + t2443 / 0.27e2 + 0.2e1 / 0.27e2 * t2446 + t2553 / 0.27e2 + t2556 + t989 * t990 * t643 / 0.45e2 - 0.2e1 / 0.45e2 * t2560 * t990 * t647 - 0.14e2 / 0.675e3 * t2564 - t2567 - 0.7e1 / 0.225e3 * t1001 * t1002 * t651 + t179 * t950 * t427 * t264 / 0.45e2 - t2576 * t990 * t676 / 0.135e3 - 0.14e2 / 0.675e3 * t989 * t1002 * t683;
1587
1588 t2588 = t963 * t1014;
1589 t2603 = t963 * t1018;
1590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1591 v4rho3sigma[6] = 0.2e1 / 0.27e2 * t2468 + t2493 - t2496 + t2486 + t2472 / 0.9e1 + t2355 / 0.27e2 - t2423 - 0.14e2 / 0.675e3 * t2331 + t2425 / 0.45e2 + t2588 / 0.27e2 - t2511 * t968 * t760 / 0.135e3 - 0.14e2 / 0.675e3 * t967 * t908 * t764 + t179 * t901 * t395 * t324 / 0.45e2 + 0.2e1 / 0.45e2 * t2515 * t968 * t753 + t2603 / 0.18e2 + t967 * t968 * t768 / 0.45e2 - 0.7e1 / 0.225e3 * t907 * t908 * t772 + t179 * t378 * t395 * t329 / 0.30e2 + t179 * t378 * t188 * t714 / 0.30e2;
1592
1593 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1594 v4rho3sigma[7] = 0.0e0;
1595
1596 t2623 = t963 * t1025;
1597 t2625 = t963 * t1032;
1598 t2627 = t963 * t1036;
1599 t2629 = t963 * t1041;
1600 t2631 = t963 * t1045;
1601 t2633 = t963 * t1049;
1602 t2650 = -0.28e2 / 0.675e3 * t989 * t1002 * t338 * t128 + t2623 / 0.27e2 + t2625 / 0.9e1 + t2627 / 0.18e2 + 0.203e3 / 0.3375e4 * t2629 - 0.14e2 / 0.225e3 * t2631 - 0.7e1 / 0.225e3 * t2633 - t2541 / 0.9e1 + t2549 + 0.2e1 / 0.27e2 * t2553 - 0.28e2 / 0.675e3 * t2564 - t980 * t401 * t427 * t7 / 0.15e2 + t179 * t401 * t199 * t718 / 0.30e2 + t179 * t401 * t1059 * t128 / 0.30e2;
1603 t2659 = t872 * t199;
1604 t2663 = t348 * t427;
1605 t2670 = t1969 * t1028;
1606 t2683 = t963 * t1060;
1607 t2685 = t2446 / 0.27e2 + 0.2e1 / 0.45e2 * t989 * t153 * t427 * t301 + t989 * t990 * t789 / 0.45e2 + 0.203e3 / 0.3375e4 * t1001 * t2659 * t705 - 0.14e2 / 0.225e3 * t1001 * t2663 * t313 - 0.7e1 / 0.225e3 * t1001 * t1002 * t709 - 0.14e2 / 0.135e3 * t2670 - t2576 * t990 * t776 / 0.135e3 - 0.2e1 / 0.45e2 * t2560 * t990 * t780 + 0.14e2 / 0.225e3 * t979 * t1000 * t1002 * t785 + t2556 - t2567 - t2453 + 0.2e1 / 0.45e2 * t2546 + t2683 / 0.30e2;
1608 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1609 v4rho3sigma[8] = t2650 + t2685;
1610
1611 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1612 v4rho3sigma[9] = t2468 / 0.9e1 + t2588 / 0.9e1 + t2603 / 0.6e1 - t2423 - t179 * t2360 * t188 * t861 / 0.135e3 + t967 * t968 * t812 / 0.15e2 + t179 * t378 * t188 * t835 / 0.30e2;
1613
1614 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1615 v4rho3sigma[10] = 0.0e0;
1616
1617 t2704 = t53 * t2234;
1618 t2743 = -t179 * t2454 * t199 * t865 / 0.135e3 - 0.8932e4 / 0.84375e5 * t179 * t2704 * t199 * t874 + t2623 / 0.9e1 + t2625 / 0.3e1 + t2627 / 0.6e1 + 0.203e3 / 0.1125e4 * t2629 - 0.14e2 / 0.75e2 * t2631 - 0.7e1 / 0.75e2 * t2633 + t179 * t950 * t427 * t333 / 0.15e2 + t179 * t401 * t1059 * t143 / 0.10e2 + t179 * t401 * t427 * t342 / 0.10e2 + t179 * t401 * t199 * t839 / 0.30e2 + 0.203e3 / 0.1875e4 * t179 * t1039 * t427 * t350 - 0.7e1 / 0.125e3 * t179 * t415 * t1059 * t171 - 0.7e1 / 0.125e3 * t179 * t415 * t427 * t369;
1619 t2766 = t178 * t5 * t53;
1620 t2787 = -0.7e1 / 0.375e3 * t179 * t415 * t199 * t857 + t2553 / 0.9e1 - 0.14e2 / 0.225e3 * t2564 - 0.14e2 / 0.225e3 * t989 * t1002 * t815 + t989 * t990 * t819 / 0.15e2 - 0.14e2 / 0.75e2 * t1001 * t2663 * t338 - 0.7e1 / 0.75e2 * t1001 * t1002 * t822 - 0.7e1 / 0.75e2 * t1001 * t1002 * t826 + 0.203e3 / 0.1875e4 * t2766 * t2659 * t830 + 0.203e3 / 0.1125e4 * t1001 * t2659 * t808 - 0.14e2 / 0.45e2 * t2670 - t2453 + t2546 / 0.15e2 + t2683 / 0.10e2 + t179 * t191 * (-0.41979259259259259259e1 * t20 * t844 - 0.79729629629629629628e1 * t31 * sigma[2] * t849 - 0.19248716727533128478e-2 * t63 * t854) / 0.50e2;
1621 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1622 v4rho3sigma[11] = t2743 + t2787;
1623
1624 t2788 = t205 * t237;
1625 t2805 = t963 * t1067;
1626 t2807 = t963 * t1071;
1627 t2809 = t963 * t1075;
1628 t2823 = t963 * t1081;
1629 t2825 = t963 * t1085;
1630 t2827 = -0.14e2 / 0.675e3 * t179 * t2788 * t431 * t207 - 0.8932e4 / 0.84375e5 * t179 * t2365 * t431 * t239 + t179 * t901 * t439 * t207 / 0.45e2 + 0.203e3 / 0.5625e4 * t179 * t920 * t439 * t239 - 0.14e2 / 0.135e3 * t2805 + 0.406e3 / 0.3375e4 * t2807 - 0.28e2 / 0.225e3 * t2809 - 0.7e1 / 0.225e3 * t179 * t1065 * t431 * t222 + 0.203e3 / 0.5625e4 * t179 * t920 * t431 * t258 - 0.14e2 / 0.375e3 * t179 * t383 * t188 * t940 + t2823 / 0.9e1 - 0.14e2 / 0.225e3 * t2825;
1631 t2844 = t564 * t431;
1632 t2854 = t237 * t439;
1633 t2859 = 0.14e2 / 0.675e3 * t2467 * t432;
1634 t2860 = t395 * t395;
1635 t2865 = t2467 * t440 / 0.45e2;
1636 t2866 = t963 * t1093;
1637 t2876 = t179 * t378 * t1092 * t89 / 0.15e2 + t179 * t378 * t439 * t222 / 0.30e2 - 0.14e2 / 0.375e3 * t179 * t383 * t1092 * t122 - 0.7e1 / 0.375e3 * t179 * t383 * t439 * t258 + 0.406e3 / 0.3375e4 * t907 * t2844 * t215 - 0.28e2 / 0.225e3 * t907 * t908 * t912 + 0.812e3 / 0.5625e4 * t2413 * t2414 * t925 - 0.14e2 / 0.225e3 * t907 * t2854 * t215 - t2859 - 0.14e2 / 0.375e3 * t179 * t383 * t2860 + t2865 + t2866 / 0.15e2 + t179 * t180 * (0.10872222222222222222e1 * t31 * t250 + 0.38497433455066256957e-3 * sigma[0] * t255) / 0.50e2;
1638 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1639 v4rho2sigma2[0] = t2827 + t2876;
1640
1641 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1642 v4rho2sigma2[1] = 0.0e0;
1643
1644 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1645 v4rho2sigma2[2] = 0.0e0;
1646
1647 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1648 v4rho2sigma2[3] = 0.0e0;
1649
1650 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1651 v4rho2sigma2[4] = 0.0e0;
1652
1653 t2878 = 0.14e2 / 0.675e3 * t2467 * t444;
1654 t2879 = t963 * t1100;
1655 t2881 = t262 * t348;
1656 t2891 = t2467 * t452 / 0.45e2;
1657 t2892 = t963 * t1106;
1658 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1659 v4rho2sigma2[5] = -t2878 - 0.14e2 / 0.135e3 * t2879 - 0.14e2 / 0.675e3 * t179 * t2881 * t443 * t264 - 0.7e1 / 0.225e3 * t179 * t1098 * t443 * t272 + t2891 + t2892 / 0.9e1 + t179 * t950 * t451 * t264 / 0.45e2 + t179 * t401 * t451 * t272 / 0.30e2;
1660
1661 t2905 = t963 * t1110;
1662 t2924 = t963 * t1114;
1663 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1664 v4rho2sigma2[6] = -t2859 - 0.7e1 / 0.135e3 * t2805 + 0.203e3 / 0.3375e4 * t2807 - 0.14e2 / 0.225e3 * t2809 - 0.7e1 / 0.135e3 * t2905 - 0.14e2 / 0.675e3 * t967 * t237 * t431 * t282 + 0.203e3 / 0.3375e4 * t907 * t2844 * t289 - 0.14e2 / 0.225e3 * t907 * t908 * t975 - 0.7e1 / 0.225e3 * t980 * t1065 * t431 * t7 + t2865 + t2823 / 0.18e2 - 0.7e1 / 0.225e3 * t2825 + t2866 / 0.30e2 + t2924 / 0.18e2 + t967 * t104 * t439 * t282 / 0.45e2 - 0.7e1 / 0.225e3 * t907 * t2854 * t289 + t179 * t378 * t1092 * t139 / 0.30e2 + t980 * t378 * t439 * t7 / 0.30e2;
1665
1666 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1667 v4rho2sigma2[7] = 0.0e0;
1668
1669 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1670 v4rho2sigma2[8] = 0.0e0;
1671
1672 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1673 v4rho2sigma2[9] = 0.0e0;
1674
1675 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1676 v4rho2sigma2[10] = 0.0e0;
1677
1678 t2942 = t963 * t1118;
1679 t2952 = t963 * t1122;
1680 t2954 = t872 * t443;
1681 t2958 = t963 * t1126;
1682 t2964 = t963 * t1130;
1683 t2974 = t963 * t1134;
1684 t2976 = t348 * t451;
1685 t2980 = t963 * t1142;
1686 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1687 v4rho2sigma2[11] = -t2878 - 0.7e1 / 0.135e3 * t2879 - 0.7e1 / 0.135e3 * t2942 - 0.14e2 / 0.675e3 * t989 * t348 * t443 * t301 + 0.7e1 / 0.225e3 * t980 * t1098 * t443 * t7 + 0.203e3 / 0.3375e4 * t2952 + 0.203e3 / 0.3375e4 * t1001 * t2954 * t313 - 0.14e2 / 0.225e3 * t2958 - 0.14e2 / 0.225e3 * t1001 * t1002 * t1008 + t2891 + t2892 / 0.18e2 + t2964 / 0.18e2 + t989 * t153 * t451 * t301 / 0.45e2 - t980 * t401 * t451 * t7 / 0.30e2 - 0.7e1 / 0.225e3 * t2974 - 0.7e1 / 0.225e3 * t1001 * t2976 * t313 + t2980 / 0.30e2 + t179 * t401 * t1141 * t128 / 0.30e2;
1688
1689 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1690 v4rho2sigma2[12] = -t2859 - 0.14e2 / 0.135e3 * t2905 - 0.14e2 / 0.675e3 * t179 * t2788 * t431 * t324 - 0.7e1 / 0.225e3 * t179 * t1065 * t431 * t329 + t2865 + t2924 / 0.9e1 + t179 * t901 * t439 * t324 / 0.45e2 + t179 * t378 * t439 * t329 / 0.30e2;
1691
1692 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1693 v4rho2sigma2[13] = 0.0e0;
1694
1695 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1696 v4rho2sigma2[14] = 0.0e0;
1697
1698 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1699 v4rho2sigma2[15] = 0.0e0;
1700
1701 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1702 v4rho2sigma2[16] = 0.0e0;
1703
1704 t3049 = t179 * t950 * t451 * t333 / 0.45e2 + 0.203e3 / 0.5625e4 * t179 * t1039 * t451 * t350 - 0.14e2 / 0.675e3 * t179 * t2881 * t443 * t333 - 0.8932e4 / 0.84375e5 * t179 * t2704 * t443 * t350 - 0.7e1 / 0.225e3 * t179 * t1098 * t443 * t342 + 0.203e3 / 0.5625e4 * t179 * t1039 * t443 * t369 - 0.14e2 / 0.375e3 * t179 * t415 * t199 * t1059 + t179 * t401 * t1141 * t143 / 0.15e2 + t179 * t401 * t451 * t342 / 0.30e2 - 0.14e2 / 0.375e3 * t179 * t415 * t1141 * t171 - 0.7e1 / 0.375e3 * t179 * t415 * t451 * t369 - 0.14e2 / 0.135e3 * t2942;
1705 t3075 = t427 * t427;
1706 t3079 = 0.406e3 / 0.3375e4 * t2952 - 0.28e2 / 0.225e3 * t2958 + t2964 / 0.9e1 - 0.14e2 / 0.225e3 * t2974 - 0.14e2 / 0.225e3 * t1001 * t2976 * t338 + 0.406e3 / 0.3375e4 * t1001 * t2954 * t338 - 0.28e2 / 0.225e3 * t1001 * t1002 * t1031 + 0.812e3 / 0.5625e4 * t2766 * t2659 * t1044 - t2878 + t2891 + t2980 / 0.15e2 + t179 * t191 * (0.10872222222222222222e1 * t31 * t361 + 0.38497433455066256957e-3 * sigma[2] * t366) / 0.50e2 - 0.14e2 / 0.375e3 * t179 * t415 * t3075;
1707 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1708 v4rho2sigma2[17] = t3049 + t3079;
1709
1710 t3081 = 0.203e3 / 0.3375e4 * t963 * t1146;
1711 t3082 = t14 * t564;
1712 t3096 = 0.7e1 / 0.75e2 * t963 * t1150;
1713 t3112 = 0.17822885858826970813e-6 * t963 * t1153;
1714 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1715 v4rhosigma3[0] = t3081 + 0.203e3 / 0.3375e4 * t179 * t3082 * t1145 * t89 - 0.8932e4 / 0.84375e5 * t179 * t2365 * t1145 * t122 + 0.203e3 / 0.1875e4 * t179 * t920 * t431 * t395 - t3096 - 0.7e1 / 0.75e2 * t907 * t908 * t1080 + 0.203e3 / 0.1875e4 * t2413 * t2414 * t1084 - 0.7e1 / 0.125e3 * t179 * t383 * t395 * t439 - 0.7e1 / 0.125e3 * t179 * t383 * t188 * t1092 + t3112 + 0.17822885858826970813e-6 * t179 * t378 * t42 * t89 - 0.99808160809431036555e-7 * t179 * t383 * t42 * t122 - 0.85549852122369459904e-6 * t179 * t180 * t119;
1716
1717 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1718 v4rhosigma3[1] = 0.0e0;
1719
1720 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1721 v4rhosigma3[2] = 0.0e0;
1722
1723 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1724 v4rhosigma3[3] = 0.0e0;
1725
1726 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1727 v4rhosigma3[4] = 0.0e0;
1728
1729 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1730 v4rhosigma3[5] = 0.0e0;
1731
1732 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1733 v4rhosigma3[6] = 0.0e0;
1734
1735 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1736 v4rhosigma3[7] = 0.0e0;
1737
1738 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1739 v4rhosigma3[8] = 0.0e0;
1740
1741 t3125 = 0.203e3 / 0.3375e4 * t963 * t1157;
1742 t3126 = t52 * t872;
1743 t3132 = 0.7e1 / 0.75e2 * t963 * t1161;
1744 t3137 = 0.17822885858826970813e-6 * t963 * t1164;
1745 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1746 v4rhosigma3[9] = t3125 + 0.203e3 / 0.3375e4 * t179 * t3126 * t1156 * t128 - t3132 - 0.7e1 / 0.75e2 * t1001 * t1002 * t1105 + t3137 + 0.17822885858826970813e-6 * t179 * t401 * t73 * t128;
1747
1748 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1749 v4rhosigma3[10] = t3081 + 0.203e3 / 0.3375e4 * t179 * t3082 * t1145 * t139 - t3096 - 0.7e1 / 0.75e2 * t907 * t908 * t1113 + t3112 + 0.17822885858826970813e-6 * t179 * t378 * t42 * t139;
1750
1751 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1752 v4rhosigma3[11] = 0.0e0;
1753
1754 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1755 v4rhosigma3[12] = 0.0e0;
1756
1757 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1758 v4rhosigma3[13] = 0.0e0;
1759
1760 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1761 v4rhosigma3[14] = 0.0e0;
1762
1763 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1764 v4rhosigma3[15] = 0.0e0;
1765
1766 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1767 v4rhosigma3[16] = 0.0e0;
1768
1769 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1770 v4rhosigma3[17] = 0.0e0;
1771
1772 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1773 v4rhosigma3[18] = 0.0e0;
1774
1775 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1776 v4rhosigma3[19] = t3125 + 0.203e3 / 0.3375e4 * t179 * t3126 * t1156 * t143 - 0.8932e4 / 0.84375e5 * t179 * t2704 * t1156 * t171 + 0.203e3 / 0.1875e4 * t179 * t1039 * t443 * t427 - t3132 - 0.7e1 / 0.75e2 * t1001 * t1002 * t1129 + 0.203e3 / 0.1875e4 * t2766 * t2659 * t1133 - 0.7e1 / 0.125e3 * t179 * t415 * t427 * t451 - 0.7e1 / 0.125e3 * t179 * t415 * t199 * t1141 + t3137 + 0.17822885858826970813e-6 * t179 * t401 * t73 * t143 - 0.99808160809431036555e-7 * t179 * t415 * t73 * t171 - 0.85549852122369459904e-6 * t179 * t191 * t168;
1777
1778 t3190 = t431 * t431;
1779 t3198 = t439 * t439;
1780 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1781 v4sigma4[0] = -0.8932e4 / 0.84375e5 * t179 * t2365 * t3190 + 0.406e3 / 0.1875e4 * t179 * t920 * t431 * t439 - 0.7e1 / 0.125e3 * t179 * t383 * t3198 - 0.39923264323772414622e-6 * t179 * t383 * t188 * t42;
1782
1783 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1784 v4sigma4[1] = 0.0e0;
1785
1786 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1787 v4sigma4[2] = 0.0e0;
1788
1789 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1790 v4sigma4[3] = 0.0e0;
1791
1792 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1793 v4sigma4[4] = 0.0e0;
1794
1795 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1796 v4sigma4[5] = 0.0e0;
1797
1798 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1799 v4sigma4[6] = 0.0e0;
1800
1801 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1802 v4sigma4[7] = 0.0e0;
1803
1804 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1805 v4sigma4[8] = 0.0e0;
1806
1807 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1808 v4sigma4[9] = 0.0e0;
1809
1810 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1811 v4sigma4[10] = 0.0e0;
1812
1813 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1814 v4sigma4[11] = 0.0e0;
1815
1816 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1817 v4sigma4[12] = 0.0e0;
1818
1819 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1820 v4sigma4[13] = 0.0e0;
1821
1822 t3206 = t443 * t443;
1823 t3214 = t451 * t451;
1824 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1825 v4sigma4[14] = -0.8932e4 / 0.84375e5 * t179 * t2704 * t3206 + 0.406e3 / 0.1875e4 * t179 * t1039 * t443 * t451 - 0.7e1 / 0.125e3 * t179 * t415 * t3214 - 0.39923264323772414622e-6 * t179 * t415 * t199 * t73;
1826
1827 #ifndef XC_DONT_COMPILE_MXC
1828
1829 if(order < 5) return;
1830
1831
1832 #endif
1833
1834 #endif
1835
1836 #endif
1837
1838 #endif
1839
1840 #endif
1841
1842
1843 }
1844
1845