1 /*
2 This file was generated automatically with ./scripts/maple2c.pl.
3 Do not edit this file directly as it can be overwritten!!
4
5 This Source Code Form is subject to the terms of the Mozilla Public
6 License, v. 2.0. If a copy of the MPL was not distributed with this
7 file, You can obtain one at http://mozilla.org/MPL/2.0/.
8
9 Maple version : Maple 2016 (X86 64 LINUX)
10 Maple source : ./maple/gga_exc/gga_x_ssb_sw.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t3, t5, t6, t7, t8, t9, t10;
24 double t11, t12, t14, t15, t16, t17, t18, t19;
25 double t20, t21, t23, t29, t30, t36, t37, t40;
26 double t42, t44, t45, t47, t51, t52, t57, t59;
27
28 #ifndef XC_DONT_COMPILE_VXC
29 double t64, t65, t66, t68, t73, t75, t76, t78;
30 double t80, t81, t82, t83, t90, t91, t92, t93;
31 double t94, t95, t96, t97, t98, t99, t101, t104;
32 double t108, t115, t122, t124, t128;
33
34 #ifndef XC_DONT_COMPILE_FXC
35 double t131, t132, t140, t145, t147, t149, t152, t153;
36 double t154, t155, t157, t159, t171, t172, t175, t179;
37 double t180, t183, t184, t187, t192, t195, t200, t202;
38 double t208, t214, t215, t216, t219, t221, t227, t230;
39 double t235, t239, t242, t244, t250;
40
41 #ifndef XC_DONT_COMPILE_KXC
42 double t254, t255, t266, t272, t277, t278, t284, t287;
43 double t288, t289, t292, t293, t304, t306, t311, t314;
44 double t317, t318, t319, t324, t325, t328, t333, t338;
45 double t346, t347, t348, t351, t353, t361, t365, t367;
46 double t370, t373, t379, t387, t391, t399, t401, t406;
47 double t412, t418, t424, t426, t433, t439;
48
49 #ifndef XC_DONT_COMPILE_LXC
50 double t454, t465, t479, t483, t485, t506, t509, t510;
51 double t516, t518, t534, t627;
52 #endif
53
54 #endif
55
56 #endif
57
58 #endif
59
60 #endif
61
62
63 gga_x_ssb_sw_params *params;
64
65 assert(p->params != NULL);
66 params = (gga_x_ssb_sw_params * )(p->params);
67
68 t1 = M_CBRT3;
69 t3 = POW_1_3(0.1e1 / M_PI);
70 t5 = M_CBRT4;
71 t6 = t5 * t5;
72 t7 = t1 * t3 * t6;
73 t8 = M_CBRT2;
74 t9 = t8 * t8;
75 t10 = POW_1_3(rho[0]);
76 t11 = t9 * t10;
77 t12 = M_CBRT6;
78 t14 = M_PI * M_PI;
79 t15 = POW_1_3(t14);
80 t16 = t15 * t15;
81 t17 = 0.1e1 / t16;
82 t18 = params->B * t12 * t17;
83 t19 = sigma[0] * t9;
84 t20 = rho[0] * rho[0];
85 t21 = t10 * t10;
86 t23 = 0.1e1 / t21 / t20;
87 t29 = 0.1e1 + params->C * t12 * t17 * t19 * t23 / 0.24e2;
88 t30 = 0.1e1 / t29;
89 t36 = params->D * t12 * t17;
90 t37 = t12 * t12;
91 t40 = 0.1e1 / t15 / t14;
92 t42 = sigma[0] * sigma[0];
93 t44 = t20 * t20;
94 t45 = t44 * rho[0];
95 t47 = 0.1e1 / t10 / t45;
96 t51 = 0.1e1 + params->E * t37 * t40 * t42 * t8 * t47 / 0.288e3;
97 t52 = 0.1e1 / t51;
98 t57 = params->A + t18 * t19 * t23 * t30 / 0.24e2 - t36 * t19 * t23 * t52 / 0.24e2;
99 t59 = t7 * t11 * t57;
100 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
101 zk[0] = -0.3e1 / 0.16e2 * t59;
102
103 #ifndef XC_DONT_COMPILE_VXC
104
105 if(order < 1) return;
106
107
108 t64 = t10 * rho[0] * t1 * t3;
109 t65 = t6 * t9;
110 t66 = t20 * rho[0];
111 t68 = 0.1e1 / t21 / t66;
112 t73 = params->B * t37;
113 t75 = t73 * t40 * t42;
114 t76 = t44 * t20;
115 t78 = 0.1e1 / t10 / t76;
116 t80 = t29 * t29;
117 t81 = 0.1e1 / t80;
118 t82 = t81 * params->C;
119 t83 = t8 * t78 * t82;
120 t90 = t14 * t14;
121 t91 = 0.1e1 / t90;
122 t92 = params->D * t91;
123 t93 = t42 * sigma[0];
124 t94 = t92 * t93;
125 t95 = t44 * t44;
126 t96 = t95 * rho[0];
127 t97 = 0.1e1 / t96;
128 t98 = t51 * t51;
129 t99 = 0.1e1 / t98;
130 t101 = t97 * t99 * params->E;
131 t104 = -t18 * t19 * t68 * t30 / 0.9e1 + t75 * t83 / 0.108e3 + t36 * t19 * t68 * t52 / 0.9e1 - t94 * t101 / 0.108e3;
132 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
133 vrho[0] = -t59 / 0.4e1 - 0.3e1 / 0.16e2 * t64 * t65 * t104;
134
135 t108 = t9 * t23;
136 t115 = t8 * t47 * t82;
137 t122 = 0.1e1 / t95;
138 t124 = t122 * t99 * params->E;
139 t128 = t65 * (t18 * t108 * t30 / 0.24e2 - t73 * t40 * sigma[0] * t115 / 0.288e3 - t36 * t108 * t52 / 0.24e2 + t92 * t42 * t124 / 0.288e3);
140 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
141 vsigma[0] = -0.3e1 / 0.16e2 * t64 * t128;
142
143 #ifndef XC_DONT_COMPILE_FXC
144
145 if(order < 2) return;
146
147
148 t131 = 0.1e1 / t21;
149 t132 = t9 * t131;
150 t140 = 0.1e1 / t21 / t44;
151 t145 = t44 * t66;
152 t147 = 0.1e1 / t10 / t145;
153 t149 = t8 * t147 * t82;
154 t152 = params->B * t91;
155 t153 = t152 * t93;
156 t154 = t95 * t20;
157 t155 = 0.1e1 / t154;
158 t157 = 0.1e1 / t80 / t29;
159 t159 = params->C * params->C;
160 t171 = t42 * t42;
161 t172 = t171 * sigma[0];
162 t175 = 0.1e1 / t10 / t95 / t145;
163 t179 = 0.1e1 / t98 / t51;
164 t180 = params->E * params->E;
165 t183 = t37 * t40 * t8;
166 t184 = t179 * t180 * t183;
167 t187 = 0.11e2 / 0.27e2 * t18 * t19 * t140 * t30 - t75 * t149 / 0.12e2 + 0.2e1 / 0.81e2 * t153 * t155 * t157 * t159 - 0.11e2 / 0.27e2 * t36 * t19 * t140 * t52 + 0.35e2 / 0.324e3 * t94 * t155 * t99 * params->E - t92 * t172 * t175 * t184 / 0.2916e4;
168 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
169 v2rho2[0] = -t7 * t132 * t57 / 0.12e2 - t7 * t11 * t104 / 0.2e1 - 0.3e1 / 0.16e2 * t64 * t65 * t187;
170
171 t192 = t10 * t1 * t3;
172 t195 = t9 * t68;
173 t200 = t73 * t40 * t8;
174 t202 = params->C * sigma[0];
175 t208 = t97 * t157 * t159;
176 t214 = t92 * t97;
177 t215 = t99 * params->E;
178 t216 = t215 * t42;
179 t219 = t95 * t76;
180 t221 = 0.1e1 / t10 / t219;
181 t227 = t65 * (-t18 * t195 * t30 / 0.9e1 + t200 * t78 * t81 * t202 / 0.36e2 - t152 * t42 * t208 / 0.108e3 + t36 * t195 * t52 / 0.9e1 - t214 * t216 / 0.27e2 + t92 * t171 * t221 * t184 / 0.7776e4);
182 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
183 v2rhosigma[0] = -t192 * t128 / 0.4e1 - 0.3e1 / 0.16e2 * t64 * t227;
184
185 t230 = t73 * t40;
186 t235 = t122 * t157 * t159;
187 t239 = t215 * sigma[0];
188 t242 = t95 * t45;
189 t244 = 0.1e1 / t10 / t242;
190 t250 = t65 * (-t230 * t115 / 0.144e3 + t152 * sigma[0] * t235 / 0.288e3 + t92 * t122 * t239 / 0.96e2 - t92 * t93 * t244 * t184 / 0.20736e5);
191 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
192 v2sigma2[0] = -0.3e1 / 0.16e2 * t64 * t250;
193
194 #ifndef XC_DONT_COMPILE_KXC
195
196 if(order < 3) return;
197
198
199 t254 = 0.1e1 / t21 / rho[0];
200 t255 = t9 * t254;
201 t266 = 0.1e1 / t21 / t45;
202 t272 = 0.1e1 / t10 / t95;
203 t277 = t95 * t66;
204 t278 = 0.1e1 / t277;
205 t284 = 0.1e1 / t21 / t242;
206 t287 = t80 * t80;
207 t288 = 0.1e1 / t287;
208 t289 = t159 * params->C;
209 t292 = t12 * t17 * t9;
210 t293 = t288 * t289 * t292;
211 t304 = t95 * t95;
212 t306 = 0.1e1 / t10 / t304;
213 t311 = t171 * t93;
214 t314 = 0.1e1 / t21 / t304 / t45;
215 t317 = t98 * t98;
216 t318 = 0.1e1 / t317;
217 t319 = t180 * params->E;
218 t324 = t12 / t16 / t90 * t9;
219 t325 = t318 * t319 * t324;
220 t328 = -0.154e3 / 0.81e2 * t18 * t19 * t266 * t30 + 0.341e3 / 0.486e3 * t75 * t8 * t272 * t82 - 0.38e2 / 0.81e2 * t153 * t278 * t157 * t159 + 0.2e1 / 0.243e3 * t152 * t171 * t284 * t293 + 0.154e3 / 0.81e2 * t36 * t19 * t266 * t52 - 0.569e3 / 0.486e3 * t94 * t278 * t99 * params->E + t92 * t172 * t306 * t184 / 0.108e3 - t92 * t311 * t314 * t325 / 0.8748e4;
221 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
222 v3rho3[0] = t7 * t255 * t57 / 0.18e2 - t7 * t132 * t104 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t187 - 0.3e1 / 0.16e2 * t64 * t65 * t328;
223
224 t333 = t131 * t1 * t3;
225 t338 = t9 * t140;
226 t346 = t152 * t155;
227 t347 = t157 * t159;
228 t348 = t347 * t42;
229 t351 = t95 * t44;
230 t353 = 0.1e1 / t21 / t351;
231 t361 = t92 * t155;
232 t365 = t92 * t175 * t179;
233 t367 = t180 * t171 * t183;
234 t370 = t171 * t42;
235 t373 = 0.1e1 / t21 / t304 / t44;
236 t379 = t65 * (0.11e2 / 0.27e2 * t18 * t338 * t30 - 0.65e2 / 0.324e3 * t200 * t147 * t81 * t202 + 0.17e2 / 0.108e3 * t346 * t348 - t152 * t93 * t353 * t293 / 0.324e3 - 0.11e2 / 0.27e2 * t36 * t338 * t52 + 0.29e2 / 0.81e2 * t361 * t216 - 0.25e2 / 0.7776e4 * t365 * t367 + t92 * t370 * t373 * t325 / 0.23328e5);
237 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
238 v3rho2sigma[0] = -t333 * t128 / 0.12e2 - t192 * t227 / 0.2e1 - 0.3e1 / 0.16e2 * t64 * t379;
239
240 t387 = t347 * sigma[0];
241 t391 = 0.1e1 / t21 / t277;
242 t399 = t92 * t221 * t179;
243 t401 = t180 * t93 * t183;
244 t406 = 0.1e1 / t21 / t304 / t66;
245 t412 = t65 * (t230 * t83 / 0.27e2 - 0.5e1 / 0.108e3 * t152 * t97 * t387 + t152 * t42 * t391 * t293 / 0.864e3 - t214 * t239 / 0.12e2 + t399 * t401 / 0.972e3 - t92 * t172 * t406 * t325 / 0.62208e5);
246 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
247 v3rhosigma2[0] = -t192 * t250 / 0.4e1 - 0.3e1 / 0.16e2 * t64 * t412;
248
249 t418 = 0.1e1 / t21 / t154;
250 t424 = t92 * t244 * t179;
251 t426 = t180 * t42 * t183;
252 t433 = 0.1e1 / t21 / t304 / t20;
253 t439 = t65 * (t152 * t235 / 0.96e2 - t152 * sigma[0] * t418 * t293 / 0.2304e4 - t424 * t426 / 0.3456e4 + t92 * t124 / 0.96e2 + t92 * t171 * t433 * t325 / 0.165888e6);
254 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
255 v3sigma3[0] = -0.3e1 / 0.16e2 * t64 * t439;
256
257 #ifndef XC_DONT_COMPILE_LXC
258
259 if(order < 4) return;
260
261
262 t454 = 0.1e1 / t21 / t76;
263 t465 = 0.1e1 / t351;
264 t479 = t172 / t10 / t304 / rho[0];
265 t483 = t159 * t159;
266 t485 = 0.1e1 / t287 / t29 * t483 * t183;
267 t506 = t90 * t90;
268 t509 = params->D / t506 / t90;
269 t510 = t171 * t171;
270 t516 = 0.1e1 / t317 / t51;
271 t518 = t180 * t180;
272 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
273 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t108 * t57 + 0.2e1 / 0.9e1 * t7 * t255 * t104 - t7 * t132 * t187 / 0.2e1 - t7 * t11 * t328 - 0.3e1 / 0.16e2 * t64 * t65 * (0.2618e4 / 0.243e3 * t18 * t19 * t454 * t30 - 0.3047e4 / 0.486e3 * t75 * t8 / t10 / t96 * t82 + 0.5126e4 / 0.729e3 * t153 * t465 * t157 * t159 - 0.196e3 / 0.729e3 * t152 * t171 / t21 / t219 * t293 + 0.16e2 / 0.2187e4 * t152 * t479 * t485 - 0.2618e4 / 0.243e3 * t36 * t19 * t454 * t52 + 0.19393e5 / 0.1458e4 * t94 * t465 * t99 * params->E - 0.5107e4 / 0.26244e5 * t92 * t479 * t184 + 0.73e2 / 0.13122e5 * t92 * t311 / t21 / t304 / t76 * t325 - 0.2e1 / 0.19683e5 * t509 * t510 * sigma[0] / t304 / t351 * t516 * t518);
274
275 t534 = t9 * t266;
276 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
277 v4rho3sigma[0] = t254 * t1 * t3 * t128 / 0.18e2 - t333 * t227 / 0.4e1 - 0.3e1 / 0.4e1 * t192 * t379 - 0.3e1 / 0.16e2 * t64 * t65 * (-0.154e3 / 0.81e2 * t18 * t534 * t30 + 0.253e3 / 0.162e3 * t200 * t272 * t81 * t202 - 0.1025e4 / 0.486e3 * t152 * t278 * t348 + 0.89e2 / 0.972e3 * t152 * t284 * t288 * t289 * t93 * t292 - 0.2e1 / 0.729e3 * t152 * t171 * t306 * t485 + 0.154e3 / 0.81e2 * t36 * t534 * t52 - 0.892e3 / 0.243e3 * t92 * t278 * t216 + 0.2189e4 / 0.34992e5 * t92 * t306 * t179 * t367 - 0.137e3 / 0.69984e5 * t92 * t314 * t318 * t319 * t370 * t324 + t509 * t510 / t304 / t277 * t516 * t518 / 0.26244e5);
278
279 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
280 v4rho2sigma2[0] = -t333 * t250 / 0.12e2 - t192 * t412 / 0.2e1 - 0.3e1 / 0.16e2 * t64 * t65 * (-0.19e2 / 0.81e2 * t230 * t149 + 0.167e3 / 0.324e3 * t346 * t387 - 0.25e2 / 0.864e3 * t152 * t353 * t288 * t289 * t42 * t292 + t152 * t93 * t175 * t485 / 0.972e3 + 0.3e1 / 0.4e1 * t361 * t239 - 0.13e2 / 0.729e3 * t365 * t401 + 0.41e2 / 0.62208e5 * t92 * t373 * t318 * t319 * t172 * t324 - t509 * t311 / t304 / t154 * t516 * t518 / 0.69984e5);
281
282 t627 = t289 * t12;
283 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
284 v4rhosigma3[0] = -t192 * t439 / 0.4e1 - 0.3e1 / 0.16e2 * t64 * t65 * (-t152 * t208 / 0.12e2 + 0.7e1 / 0.864e3 * t152 * t391 * t288 * t627 * t17 * sigma[0] * t9 - t152 * t42 * t221 * t485 / 0.2592e4 + 0.11e2 / 0.2592e4 * t399 * t426 - 0.13e2 / 0.62208e5 * t92 * t406 * t318 * t319 * t171 * t324 - t92 * t101 / 0.12e2 + t509 * t370 / t304 / t96 * t516 * t518 / 0.186624e6);
285
286 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
287 v4sigma4[0] = -0.3e1 / 0.16e2 * t64 * t65 * (-t152 * t418 * t288 * t627 * t17 * t9 / 0.576e3 + t152 * sigma[0] * t244 * t485 / 0.6912e4 + 0.5e1 / 0.82944e5 * t92 * t433 * t318 * t319 * t93 * t324 - 0.5e1 / 0.6912e4 * t424 * t180 * sigma[0] * t183 - t509 * t172 / t304 / t95 * t516 * t518 / 0.497664e6);
288
289 #ifndef XC_DONT_COMPILE_MXC
290
291 if(order < 5) return;
292
293
294 #endif
295
296 #endif
297
298 #endif
299
300 #endif
301
302 #endif
303
304
305 }
306
307
308 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 *))309 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
310 {
311
312 #ifndef XC_DONT_COMPILE_EXC
313 double t1, t3, t4, t5, t6, t7, t8, t9;
314 double t10, t11, t12, t13, t14, t15, t16, t17;
315 double t19, t20, t22, t26, t27, t31, t32, t33;
316 double t36, t37, t39, t40, t42, t46, t47, t51;
317 double t53;
318
319 #ifndef XC_DONT_COMPILE_VXC
320 double t57, t58, t59, t61, t62, t67, t68, t70;
321 double t72, t73, t74, t81, t82, t83, t84, t85;
322 double t86, t87, t88, t89, t90, t92, t95, t99;
323 double t111, t113, t117;
324
325 #ifndef XC_DONT_COMPILE_FXC
326 double t120, t121, t128, t129, t133, t135, t140, t141;
327 double t142, t143, t145, t147, t158, t159, t162, t166;
328 double t167, t169, t170, t173, t177, t180, t184, t185;
329 double t191, t197, t198, t199, t202, t204, t210, t219;
330 double t223, t226, t228, t234;
331
332 #ifndef XC_DONT_COMPILE_KXC
333 double t238, t239, t250, t251, t256, t261, t262, t268;
334 double t271, t272, t273, t275, t276, t286, t288, t293;
335 double t296, t299, t300, t301, t305, t306, t309, t313;
336 double t317, t321, t325, t326, t327, t330, t332, t340;
337 double t344, t346, t349, t352, t358, t367, t371, t379;
338 double t381, t386, t392, t398, t404, t406, t413, t419;
339
340 #ifndef XC_DONT_COMPILE_LXC
341 double t436, t446, t460, t464, t466, t486, t489, t490;
342 double t496, t498, t513;
343 #endif
344
345 #endif
346
347 #endif
348
349 #endif
350
351 #endif
352
353
354 gga_x_ssb_sw_params *params;
355
356 assert(p->params != NULL);
357 params = (gga_x_ssb_sw_params * )(p->params);
358
359 t1 = M_CBRT3;
360 t3 = POW_1_3(0.1e1 / M_PI);
361 t4 = t1 * t3;
362 t5 = M_CBRT4;
363 t6 = t5 * t5;
364 t7 = POW_1_3(rho[0]);
365 t8 = t6 * t7;
366 t9 = M_CBRT6;
367 t10 = params->B * t9;
368 t11 = M_PI * M_PI;
369 t12 = POW_1_3(t11);
370 t13 = t12 * t12;
371 t14 = 0.1e1 / t13;
372 t15 = t10 * t14;
373 t16 = rho[0] * rho[0];
374 t17 = t7 * t7;
375 t19 = 0.1e1 / t17 / t16;
376 t20 = sigma[0] * t19;
377 t22 = t14 * sigma[0];
378 t26 = 0.1e1 + params->C * t9 * t22 * t19 / 0.24e2;
379 t27 = 0.1e1 / t26;
380 t31 = params->D * t9;
381 t32 = t31 * t14;
382 t33 = t9 * t9;
383 t36 = 0.1e1 / t12 / t11;
384 t37 = sigma[0] * sigma[0];
385 t39 = t16 * t16;
386 t40 = t39 * rho[0];
387 t42 = 0.1e1 / t7 / t40;
388 t46 = 0.1e1 + params->E * t33 * t36 * t37 * t42 / 0.576e3;
389 t47 = 0.1e1 / t46;
390 t51 = params->A + t15 * t20 * t27 / 0.24e2 - t32 * t20 * t47 / 0.24e2;
391 t53 = t4 * t8 * t51;
392 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
393 zk[0] = -0.3e1 / 0.8e1 * t53;
394
395 #ifndef XC_DONT_COMPILE_VXC
396
397 if(order < 1) return;
398
399
400 t57 = t7 * rho[0] * t1;
401 t58 = t3 * t6;
402 t59 = t16 * rho[0];
403 t61 = 0.1e1 / t17 / t59;
404 t62 = sigma[0] * t61;
405 t67 = params->B * t33 * t36;
406 t68 = t39 * t16;
407 t70 = 0.1e1 / t7 / t68;
408 t72 = t26 * t26;
409 t73 = 0.1e1 / t72;
410 t74 = t73 * params->C;
411 t81 = t11 * t11;
412 t82 = 0.1e1 / t81;
413 t83 = params->D * t82;
414 t84 = t37 * sigma[0];
415 t85 = t83 * t84;
416 t86 = t39 * t39;
417 t87 = t86 * rho[0];
418 t88 = 0.1e1 / t87;
419 t89 = t46 * t46;
420 t90 = 0.1e1 / t89;
421 t92 = t88 * t90 * params->E;
422 t95 = -t15 * t62 * t27 / 0.9e1 + t67 * t37 * t70 * t74 / 0.216e3 + t32 * t62 * t47 / 0.9e1 - t85 * t92 / 0.432e3;
423 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
424 vrho[0] = -t53 / 0.2e1 - 0.3e1 / 0.8e1 * t57 * t58 * t95;
425
426 t99 = t14 * t19;
427 t111 = 0.1e1 / t86;
428 t113 = t111 * t90 * params->E;
429 t117 = t58 * (t10 * t99 * t27 / 0.24e2 - t67 * sigma[0] * t42 * t74 / 0.576e3 - t31 * t99 * t47 / 0.24e2 + t83 * t37 * t113 / 0.1152e4);
430 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
431 vsigma[0] = -0.3e1 / 0.8e1 * t57 * t117;
432
433 #ifndef XC_DONT_COMPILE_FXC
434
435 if(order < 2) return;
436
437
438 t120 = 0.1e1 / t17;
439 t121 = t6 * t120;
440 t128 = 0.1e1 / t17 / t39;
441 t129 = sigma[0] * t128;
442 t133 = t39 * t59;
443 t135 = 0.1e1 / t7 / t133;
444 t140 = params->B * t82;
445 t141 = t140 * t84;
446 t142 = t86 * t16;
447 t143 = 0.1e1 / t142;
448 t145 = 0.1e1 / t72 / t26;
449 t147 = params->C * params->C;
450 t158 = t37 * t37;
451 t159 = t158 * sigma[0];
452 t162 = 0.1e1 / t7 / t86 / t133;
453 t166 = 0.1e1 / t89 / t46;
454 t167 = params->E * params->E;
455 t169 = t33 * t36;
456 t170 = t166 * t167 * t169;
457 t173 = 0.11e2 / 0.27e2 * t15 * t129 * t27 - t67 * t37 * t135 * t74 / 0.24e2 + t141 * t143 * t145 * t147 / 0.162e3 - 0.11e2 / 0.27e2 * t32 * t129 * t47 + 0.35e2 / 0.1296e4 * t85 * t143 * t90 * params->E - t83 * t159 * t162 * t170 / 0.23328e5;
458 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
459 v2rho2[0] = -t4 * t121 * t51 / 0.6e1 - t4 * t8 * t95 - 0.3e1 / 0.8e1 * t57 * t58 * t173;
460
461 t177 = t7 * t1;
462 t180 = t14 * t61;
463 t184 = t70 * t73;
464 t185 = params->C * sigma[0];
465 t191 = t88 * t145 * t147;
466 t197 = t83 * t88;
467 t198 = t90 * params->E;
468 t199 = t198 * t37;
469 t202 = t86 * t68;
470 t204 = 0.1e1 / t7 / t202;
471 t210 = t58 * (-t10 * t180 * t27 / 0.9e1 + t67 * t184 * t185 / 0.72e2 - t140 * t37 * t191 / 0.432e3 + t31 * t180 * t47 / 0.9e1 - t197 * t199 / 0.108e3 + t83 * t158 * t204 * t170 / 0.62208e5);
472 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
473 v2rhosigma[0] = -t177 * t117 / 0.2e1 - 0.3e1 / 0.8e1 * t57 * t210;
474
475 t219 = t111 * t145 * t147;
476 t223 = t198 * sigma[0];
477 t226 = t86 * t40;
478 t228 = 0.1e1 / t7 / t226;
479 t234 = t58 * (-t67 * t42 * t73 * params->C / 0.288e3 + t140 * sigma[0] * t219 / 0.1152e4 + t83 * t111 * t223 / 0.384e3 - t83 * t84 * t228 * t170 / 0.165888e6);
480 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
481 v2sigma2[0] = -0.3e1 / 0.8e1 * t57 * t234;
482
483 #ifndef XC_DONT_COMPILE_KXC
484
485 if(order < 3) return;
486
487
488 t238 = 0.1e1 / t17 / rho[0];
489 t239 = t6 * t238;
490 t250 = 0.1e1 / t17 / t40;
491 t251 = sigma[0] * t250;
492 t256 = 0.1e1 / t7 / t86;
493 t261 = t86 * t59;
494 t262 = 0.1e1 / t261;
495 t268 = 0.1e1 / t17 / t226;
496 t271 = t72 * t72;
497 t272 = 0.1e1 / t271;
498 t273 = t147 * params->C;
499 t275 = t9 * t14;
500 t276 = t272 * t273 * t275;
501 t286 = t86 * t86;
502 t288 = 0.1e1 / t7 / t286;
503 t293 = t158 * t84;
504 t296 = 0.1e1 / t17 / t286 / t40;
505 t299 = t89 * t89;
506 t300 = 0.1e1 / t299;
507 t301 = t167 * params->E;
508 t305 = t9 / t13 / t81;
509 t306 = t300 * t301 * t305;
510 t309 = -0.154e3 / 0.81e2 * t15 * t251 * t27 + 0.341e3 / 0.972e3 * t67 * t37 * t256 * t74 - 0.19e2 / 0.162e3 * t141 * t262 * t145 * t147 + t140 * t158 * t268 * t276 / 0.486e3 + 0.154e3 / 0.81e2 * t32 * t251 * t47 - 0.569e3 / 0.1944e4 * t85 * t262 * t90 * params->E + t83 * t159 * t288 * t170 / 0.864e3 - t83 * t293 * t296 * t306 / 0.139968e6;
511 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
512 v3rho3[0] = t4 * t239 * t51 / 0.9e1 - t4 * t121 * t95 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t173 - 0.3e1 / 0.8e1 * t57 * t58 * t309;
513
514 t313 = t120 * t1;
515 t317 = t14 * t128;
516 t321 = t135 * t73;
517 t325 = t140 * t143;
518 t326 = t145 * t147;
519 t327 = t326 * t37;
520 t330 = t86 * t39;
521 t332 = 0.1e1 / t17 / t330;
522 t340 = t83 * t143;
523 t344 = t83 * t162 * t166;
524 t346 = t167 * t158 * t169;
525 t349 = t158 * t37;
526 t352 = 0.1e1 / t17 / t286 / t39;
527 t358 = t58 * (0.11e2 / 0.27e2 * t10 * t317 * t27 - 0.65e2 / 0.648e3 * t67 * t321 * t185 + 0.17e2 / 0.432e3 * t325 * t327 - t140 * t84 * t332 * t276 / 0.1296e4 - 0.11e2 / 0.27e2 * t31 * t317 * t47 + 0.29e2 / 0.324e3 * t340 * t199 - 0.25e2 / 0.62208e5 * t344 * t346 + t83 * t349 * t352 * t306 / 0.373248e6);
528 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
529 v3rho2sigma[0] = -t313 * t117 / 0.6e1 - t177 * t210 - 0.3e1 / 0.8e1 * t57 * t358;
530
531 t367 = t326 * sigma[0];
532 t371 = 0.1e1 / t17 / t261;
533 t379 = t83 * t204 * t166;
534 t381 = t167 * t84 * t169;
535 t386 = 0.1e1 / t17 / t286 / t59;
536 t392 = t58 * (t67 * t184 * params->C / 0.54e2 - 0.5e1 / 0.432e3 * t140 * t88 * t367 + t140 * t37 * t371 * t276 / 0.3456e4 - t197 * t223 / 0.48e2 + t379 * t381 / 0.7776e4 - t83 * t159 * t386 * t306 / 0.995328e6);
537 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
538 v3rhosigma2[0] = -t177 * t234 / 0.2e1 - 0.3e1 / 0.8e1 * t57 * t392;
539
540 t398 = 0.1e1 / t17 / t142;
541 t404 = t83 * t228 * t166;
542 t406 = t167 * t37 * t169;
543 t413 = 0.1e1 / t17 / t286 / t16;
544 t419 = t58 * (t140 * t219 / 0.384e3 - t140 * sigma[0] * t398 * t276 / 0.9216e4 - t404 * t406 / 0.27648e5 + t83 * t113 / 0.384e3 + t83 * t158 * t413 * t306 / 0.2654208e7);
545 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
546 v3sigma3[0] = -0.3e1 / 0.8e1 * t57 * t419;
547
548 #ifndef XC_DONT_COMPILE_LXC
549
550 if(order < 4) return;
551
552
553 t436 = sigma[0] / t17 / t68;
554 t446 = 0.1e1 / t330;
555 t460 = t159 / t7 / t286 / rho[0];
556 t464 = t147 * t147;
557 t466 = 0.1e1 / t271 / t26 * t464 * t169;
558 t486 = t81 * t81;
559 t489 = params->D / t486 / t81;
560 t490 = t158 * t158;
561 t496 = 0.1e1 / t299 / t46;
562 t498 = t167 * t167;
563 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
564 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t19 * t51 + 0.4e1 / 0.9e1 * t4 * t239 * t95 - t4 * t121 * t173 - 0.2e1 * t4 * t8 * t309 - 0.3e1 / 0.8e1 * t57 * t58 * (0.2618e4 / 0.243e3 * t15 * t436 * t27 - 0.3047e4 / 0.972e3 * t67 * t37 / t7 / t87 * t74 + 0.2563e4 / 0.1458e4 * t141 * t446 * t145 * t147 - 0.49e2 / 0.729e3 * t140 * t158 / t17 / t202 * t276 + 0.2e1 / 0.2187e4 * t140 * t460 * t466 - 0.2618e4 / 0.243e3 * t32 * t436 * t47 + 0.19393e5 / 0.5832e4 * t85 * t446 * t90 * params->E - 0.5107e4 / 0.209952e6 * t83 * t460 * t170 + 0.73e2 / 0.209952e6 * t83 * t293 / t17 / t286 / t68 * t306 - t489 * t490 * sigma[0] / t286 / t330 * t496 * t498 / 0.629856e6);
565
566 t513 = t14 * t250;
567 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
568 v4rho3sigma[0] = t238 * t1 * t117 / 0.9e1 - t313 * t210 / 0.2e1 - 0.3e1 / 0.2e1 * t177 * t358 - 0.3e1 / 0.8e1 * t57 * t58 * (-0.154e3 / 0.81e2 * t10 * t513 * t27 + 0.253e3 / 0.324e3 * t67 * t256 * t73 * t185 - 0.1025e4 / 0.1944e4 * t140 * t262 * t327 + 0.89e2 / 0.3888e4 * t140 * t268 * t272 * t273 * t84 * t275 - t140 * t158 * t288 * t466 / 0.2916e4 + 0.154e3 / 0.81e2 * t31 * t513 * t47 - 0.223e3 / 0.243e3 * t83 * t262 * t199 + 0.2189e4 / 0.279936e6 * t83 * t288 * t166 * t346 - 0.137e3 / 0.1119744e7 * t83 * t296 * t300 * t301 * t349 * t305 + t489 * t490 / t286 / t261 * t496 * t498 / 0.1679616e7);
569
570 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
571 v4rho2sigma2[0] = -t313 * t234 / 0.6e1 - t177 * t392 - 0.3e1 / 0.8e1 * t57 * t58 * (-0.19e2 / 0.162e3 * t67 * t321 * params->C + 0.167e3 / 0.1296e4 * t325 * t367 - 0.25e2 / 0.3456e4 * t140 * t332 * t272 * t273 * t37 * t275 + t140 * t84 * t162 * t466 / 0.7776e4 + 0.3e1 / 0.16e2 * t340 * t223 - 0.13e2 / 0.5832e4 * t344 * t381 + 0.41e2 / 0.995328e6 * t83 * t352 * t300 * t301 * t159 * t305 - t489 * t293 / t286 / t142 * t496 * t498 / 0.4478976e7);
572
573 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
574 v4rhosigma3[0] = -t177 * t419 / 0.2e1 - 0.3e1 / 0.8e1 * t57 * t58 * (-t140 * t191 / 0.48e2 + 0.7e1 / 0.3456e4 * t140 * t371 * t272 * t273 * t9 * t22 - t140 * t37 * t204 * t466 / 0.20736e5 + 0.11e2 / 0.20736e5 * t379 * t406 - 0.13e2 / 0.995328e6 * t83 * t386 * t300 * t301 * t158 * t305 - t83 * t92 / 0.48e2 + t489 * t349 / t286 / t87 * t496 * t498 / 0.11943936e8);
575
576 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
577 v4sigma4[0] = -0.3e1 / 0.8e1 * t57 * t58 * (-t140 * t398 * t276 / 0.2304e4 + t140 * sigma[0] * t228 * t466 / 0.55296e5 + 0.5e1 / 0.1327104e7 * t83 * t413 * t300 * t301 * t84 * t305 - 0.5e1 / 0.55296e5 * t404 * t167 * sigma[0] * t169 - t489 * t159 / t286 / t86 * t496 * t498 / 0.31850496e8);
578
579 #ifndef XC_DONT_COMPILE_MXC
580
581 if(order < 5) return;
582
583
584 #endif
585
586 #endif
587
588 #endif
589
590 #endif
591
592 #endif
593
594
595 }
596
597
598 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 *))599 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
600 {
601
602 #ifndef XC_DONT_COMPILE_EXC
603 double t1, t3, t4, t5, t6, t7, t8, t9;
604 double t10, t11, t13, t14, t15, t16, t17, t18;
605 double t19, t20, t21, t22, t23, t24, t25, t26;
606 double t27, t29, t30, t31, t32, t36, t37, t41;
607 double t42, t43, t44, t46, t47, t49, t50, t52;
608 double t56, t57, t61, t63, t65, t66, t67, t68;
609 double t69, t70, t71, t73, t74, t75, t79, t80;
610 double t84, t86, t87, t89, t93, t94, t98, t100;
611
612 #ifndef XC_DONT_COMPILE_VXC
613 double t102, t103, t104, t105, t106, t107, t109, t110;
614 double t112, t113, t114, t115, t116, t118, t119, t120;
615 double t122, t123, t128, t129, t131, t133, t134, t135;
616 double t142, t143, t144, t145, t146, t147, t148, t149;
617 double t150, t151, t153, t156, t158, t159, t160, t161;
618 double t162, t164, t165, t166, t168, t169, t173, t174;
619 double t176, t177, t178, t179, t181, t182, t183, t185;
620 double t186, t190, t192, t194, t195, t196, t203, t204;
621 double t205, t206, t207, t208, t209, t211, t214, t216;
622 double t217, t221, t222, t223, t224, t236, t238, t241;
623 double t245, t246, t258, t260, t263;
624
625 #ifndef XC_DONT_COMPILE_FXC
626 double t267, t269, t270, t271, t272, t273, t274, t276;
627 double t277, t278, t280, t281, t282, t284, t285, t286;
628 double t287, t288, t289, t291, t292, t294, t295, t297;
629 double t298, t300, t301, t303, t304, t308, t310, t315;
630 double t316, t317, t318, t320, t322, t333, t334, t337;
631 double t341, t342, t344, t345, t348, t350, t351, t352;
632 double t353, t354, t355, t356, t358, t359, t360, t362;
633 double t363, t364, t365, t367, t368, t369, t371, t372;
634 double t373, t375, t376, t377, t378, t380, t383, t385;
635 double t387, t390, t391, t392, t394, t398, t399, t400;
636 double t401, t403, t406, t408, t409, t411, t414, t416;
637 double t419, t421, t424, t425, t427, t428, t429, t430;
638 double t431, t433, t434, t435, t436, t438, t439, t440;
639 double t441, t443, t444, t445, t447, t448, t449, t451;
640 double t452, t456, t458, t463, t464, t465, t467, t479;
641 double t480, t483, t487, t489, t492, t494, t495, t496;
642 double t499, t501, t502, t506, t510, t511, t517, t523;
643 double t524, t525, t528, t530, t535, t540, t541, t542;
644 double t547, t552, t555, t559, t560, t566, t572, t573;
645 double t574, t577, t579, t584, t585, t594, t598, t601;
646 double t603, t608, t609, t618, t622, t625, t627, t632;
647 double t633;
648
649 #ifndef XC_DONT_COMPILE_KXC
650 double t640, t642, t646, t647, t648, t649, t651, t652;
651 double t653, t655, t657, t659, t660, t661, t663, t664;
652 double t666, t667, t668, t669, t671, t672, t673, t675;
653 double t677, t679, t680, t681, t683, t684, t685, t687;
654 double t688, t690, t693, t695, t697, t698, t699, t701;
655 double t702, t703, t704, t705, t707, t708, t710, t711;
656 double t712, t714, t715, t717, t719, t720, t722, t723;
657 double t724, t726, t727, t729, t730, t732, t735, t736;
658 double t741, t746, t747, t753, t756, t757, t758, t760;
659 double t761, t771, t773, t778, t781, t784, t785, t786;
660 double t790, t791, t794, t796, t797, t798, t803, t804;
661 double t808, t811, t814, t826, t829, t830, t832, t834;
662 double t835, t837, t839, t841, t843, t846, t847, t849;
663 double t850, t851, t853, t855, t857, t859, t861, t863;
664 double t865, t867, t869, t871, t873, t875, t876, t878;
665 double t880, t881, t882, t884, t886, t888, t889, t890;
666 double t891, t893, t895, t896, t897, t899, t901, t903;
667 double t904, t905, t908, t914, t925, t928, t931, t933;
668 double t935, t938, t940, t942, t944, t946, t949, t950;
669 double t952, t954, t956, t957, t958, t960, t962, t964;
670 double t966, t968, t970, t972, t974, t976, t978, t980;
671 double t982, t984, t986, t988, t990, t992, t994, t997;
672 double t1005, t1006, t1011, t1012, t1014, t1015, t1016, t1018;
673 double t1019, t1020, t1021, t1023, t1024, t1026, t1028, t1029;
674 double t1030, t1033, t1034, t1036, t1037, t1038, t1040, t1041;
675 double t1042, t1044, t1045, t1046, t1048, t1049, t1050, t1051;
676 double t1053, t1054, t1056, t1057, t1062, t1067, t1068, t1074;
677 double t1077, t1078, t1080, t1090, t1092, t1097, t1100, t1103;
678 double t1104, t1106, t1109, t1111, t1112, t1113, t1117, t1121;
679 double t1123, t1124, t1128, t1131, t1135, t1139, t1143, t1144;
680 double t1145, t1148, t1150, t1158, t1162, t1164, t1167, t1170;
681 double t1175, t1180, t1184, t1185, t1189, t1195, t1196, t1197;
682 double t1199, t1200, t1201, t1202, t1206, t1209, t1210, t1212;
683 double t1216, t1218, t1219, t1223, t1226, t1229, t1234, t1238;
684 double t1243, t1247, t1250, t1253, t1257, t1261, t1262, t1263;
685 double t1266, t1268, t1276, t1280, t1282, t1285, t1288, t1293;
686 double t1294, t1297, t1300, t1307, t1311, t1319, t1321, t1326;
687 double t1331, t1332, t1335, t1337, t1341, t1346, t1353, t1357;
688 double t1365, t1367, t1372, t1377, t1378, t1384, t1390, t1392;
689 double t1399, t1404, t1405, t1411, t1417, t1419, t1426, t1431;
690 double t1432;
691
692 #ifndef XC_DONT_COMPILE_LXC
693 double t1446, t1450, t1453, t1458, t1468, t1482, t1486, t1488;
694 double t1508, t1511, t1512, t1518, t1520, t1529, t1532, t1535;
695 double t1542, t1553, t1565, t1569, t1570, t1576, t1578, t1580;
696 double t1581, t1582, t1587, t1589, t1591, t1593, t1595, t1600;
697 double t1602, t1612, t1614, t1616, t1622, t1629, t1631, t1632;
698 double t1633, t1638, t1640, t1642, t1644, t1647, t1650, t1656;
699 double t1659, t1662, t1673, t1676, t1678, t1681, t1683, t1686;
700 double t1692, t1698, t1707, t1708, t1710, t1713, t1716, t1719;
701 double t1723, t1725, t1730, t1737, t1739, t1741, t1742, t1744;
702 double t1751, t1769, t1770, t1771, t1773, t1785, t1786, t1790;
703 double t1797, t1799, t1803, t1810, t1816, t1818, t1821, t1827;
704 double t1834, t1840, t1844, t1867, t1876, t1879, t1887, t1891;
705 double t1893, t1894, t1899, t1902, t1905, t1923, t1926, t1929;
706 double t1932, t1934, t1940, t1947, t1953, t1956, t1962, t1964;
707 double t1966, t1998, t2013, t2018, t2031, t2034, t2041, t2048;
708 double t2055, t2062, t2071, t2096, t2108, t2120, t2131, t2146;
709 double t2149, t2152, t2155, t2158, t2161, t2164, t2204, t2219;
710 double t2226, t2235, t2245, t2248, t2267, t2284, t2303, t2309;
711 double t2310, t2321, t2326, t2339, t2347, t2352, t2357, t2362;
712 double t2376, t2404, t2414, t2428, t2433, t2453, t2459, t2473;
713 double t2502, t2505, t2507, t2510, t2512, t2561, t2564, t2568;
714 double t2586, t2589, t2599, t2602, t2605, t2609, t2626, t2637;
715 double t2638, t2640, t2642, t2644, t2645, t2657, t2661, t2668;
716 double t2671, t2672, t2674, t2687, t2690, t2703, t2715, t2732;
717 double t2740, t2749, t2758, t2778, t2856, t2857, t2859, t2908;
718 double t2909, t2921, t2936, t2946, t3009, t3019, t3048;
719 #endif
720
721 #endif
722
723 #endif
724
725 #endif
726
727 #endif
728
729
730 gga_x_ssb_sw_params *params;
731
732 assert(p->params != NULL);
733 params = (gga_x_ssb_sw_params * )(p->params);
734
735 t1 = M_CBRT3;
736 t3 = POW_1_3(0.1e1 / M_PI);
737 t4 = t1 * t3;
738 t5 = M_CBRT4;
739 t6 = t5 * t5;
740 t7 = t4 * t6;
741 t8 = rho[0] - rho[1];
742 t9 = rho[0] + rho[1];
743 t10 = 0.1e1 / t9;
744 t11 = t8 * t10;
745 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
746 t14 = POW_1_3(t13);
747 t15 = t14 * t13;
748 t16 = POW_1_3(t9);
749 t17 = t15 * t16;
750 t18 = M_CBRT6;
751 t19 = params->B * t18;
752 t20 = M_PI * M_PI;
753 t21 = POW_1_3(t20);
754 t22 = t21 * t21;
755 t23 = 0.1e1 / t22;
756 t24 = t19 * t23;
757 t25 = rho[0] * rho[0];
758 t26 = POW_1_3(rho[0]);
759 t27 = t26 * t26;
760 t29 = 0.1e1 / t27 / t25;
761 t30 = sigma[0] * t29;
762 t31 = params->C * t18;
763 t32 = t23 * sigma[0];
764 t36 = 0.1e1 + t31 * t32 * t29 / 0.24e2;
765 t37 = 0.1e1 / t36;
766 t41 = params->D * t18;
767 t42 = t41 * t23;
768 t43 = t18 * t18;
769 t44 = params->E * t43;
770 t46 = 0.1e1 / t21 / t20;
771 t47 = sigma[0] * sigma[0];
772 t49 = t25 * t25;
773 t50 = t49 * rho[0];
774 t52 = 0.1e1 / t26 / t50;
775 t56 = 0.1e1 + t44 * t46 * t47 * t52 / 0.576e3;
776 t57 = 0.1e1 / t56;
777 t61 = params->A + t24 * t30 * t37 / 0.24e2 - t42 * t30 * t57 / 0.24e2;
778 t63 = t7 * t17 * t61;
779 t65 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
780 t66 = POW_1_3(t65);
781 t67 = t66 * t65;
782 t68 = t67 * t16;
783 t69 = rho[1] * rho[1];
784 t70 = POW_1_3(rho[1]);
785 t71 = t70 * t70;
786 t73 = 0.1e1 / t71 / t69;
787 t74 = sigma[2] * t73;
788 t75 = t23 * sigma[2];
789 t79 = 0.1e1 + t31 * t75 * t73 / 0.24e2;
790 t80 = 0.1e1 / t79;
791 t84 = sigma[2] * sigma[2];
792 t86 = t69 * t69;
793 t87 = t86 * rho[1];
794 t89 = 0.1e1 / t70 / t87;
795 t93 = 0.1e1 + t44 * t46 * t84 * t89 / 0.576e3;
796 t94 = 0.1e1 / t93;
797 t98 = params->A + t24 * t74 * t80 / 0.24e2 - t42 * t74 * t94 / 0.24e2;
798 t100 = t7 * t68 * t98;
799 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
800 zk[0] = -0.3e1 / 0.8e1 * t63 - 0.3e1 / 0.8e1 * t100;
801
802 #ifndef XC_DONT_COMPILE_VXC
803
804 if(order < 1) return;
805
806
807 t102 = 0.3e1 / 0.8e1 * t63;
808 t103 = 0.3e1 / 0.8e1 * t100;
809 t104 = t14 * t16;
810 t105 = t9 * t9;
811 t106 = 0.1e1 / t105;
812 t107 = t8 * t106;
813 t109 = t10 / 0.2e1 - t107 / 0.2e1;
814 t110 = t61 * t109;
815 t112 = t7 * t104 * t110;
816 t113 = t112 / 0.2e1;
817 t114 = t16 * t16;
818 t115 = 0.1e1 / t114;
819 t116 = t15 * t115;
820 t118 = t7 * t116 * t61;
821 t119 = t118 / 0.8e1;
822 t120 = t25 * rho[0];
823 t122 = 0.1e1 / t27 / t120;
824 t123 = sigma[0] * t122;
825 t128 = params->B * t43 * t46;
826 t129 = t49 * t25;
827 t131 = 0.1e1 / t26 / t129;
828 t133 = t36 * t36;
829 t134 = 0.1e1 / t133;
830 t135 = t134 * params->C;
831 t142 = t20 * t20;
832 t143 = 0.1e1 / t142;
833 t144 = params->D * t143;
834 t145 = t47 * sigma[0];
835 t146 = t144 * t145;
836 t147 = t49 * t49;
837 t148 = t147 * rho[0];
838 t149 = 0.1e1 / t148;
839 t150 = t56 * t56;
840 t151 = 0.1e1 / t150;
841 t153 = t149 * t151 * params->E;
842 t156 = -t24 * t123 * t37 / 0.9e1 + t128 * t47 * t131 * t135 / 0.216e3 + t42 * t123 * t57 / 0.9e1 - t146 * t153 / 0.432e3;
843 t158 = t7 * t17 * t156;
844 t159 = 0.3e1 / 0.8e1 * t158;
845 t160 = t66 * t16;
846 t161 = -t109;
847 t162 = t98 * t161;
848 t164 = t7 * t160 * t162;
849 t165 = t164 / 0.2e1;
850 t166 = t67 * t115;
851 t168 = t7 * t166 * t98;
852 t169 = t168 / 0.8e1;
853 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
854 vrho[0] = -t102 - t103 + t9 * (-t113 - t119 - t159 - t165 - t169);
855
856 t173 = -t10 / 0.2e1 - t107 / 0.2e1;
857 t174 = t61 * t173;
858 t176 = t7 * t104 * t174;
859 t177 = t176 / 0.2e1;
860 t178 = -t173;
861 t179 = t98 * t178;
862 t181 = t7 * t160 * t179;
863 t182 = t181 / 0.2e1;
864 t183 = t69 * rho[1];
865 t185 = 0.1e1 / t71 / t183;
866 t186 = sigma[2] * t185;
867 t190 = t86 * t69;
868 t192 = 0.1e1 / t70 / t190;
869 t194 = t79 * t79;
870 t195 = 0.1e1 / t194;
871 t196 = t195 * params->C;
872 t203 = t84 * sigma[2];
873 t204 = t144 * t203;
874 t205 = t86 * t86;
875 t206 = t205 * rho[1];
876 t207 = 0.1e1 / t206;
877 t208 = t93 * t93;
878 t209 = 0.1e1 / t208;
879 t211 = t207 * t209 * params->E;
880 t214 = -t24 * t186 * t80 / 0.9e1 + t128 * t84 * t192 * t196 / 0.216e3 + t42 * t186 * t94 / 0.9e1 - t204 * t211 / 0.432e3;
881 t216 = t7 * t68 * t214;
882 t217 = 0.3e1 / 0.8e1 * t216;
883 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
884 vrho[1] = -t102 - t103 + t9 * (-t177 - t119 - t182 - t169 - t217);
885
886 t221 = t16 * t9 * t1;
887 t222 = t221 * t3;
888 t223 = t6 * t15;
889 t224 = t23 * t29;
890 t236 = 0.1e1 / t147;
891 t238 = t236 * t151 * params->E;
892 t241 = t19 * t224 * t37 / 0.24e2 - t128 * sigma[0] * t52 * t135 / 0.576e3 - t41 * t224 * t57 / 0.24e2 + t144 * t47 * t238 / 0.1152e4;
893 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
894 vsigma[0] = -0.3e1 / 0.8e1 * t222 * t223 * t241;
895
896 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
897 vsigma[1] = 0.0e0;
898
899 t245 = t6 * t67;
900 t246 = t23 * t73;
901 t258 = 0.1e1 / t205;
902 t260 = t258 * t209 * params->E;
903 t263 = t19 * t246 * t80 / 0.24e2 - t128 * sigma[2] * t89 * t196 / 0.576e3 - t41 * t246 * t94 / 0.24e2 + t144 * t84 * t260 / 0.1152e4;
904 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
905 vsigma[2] = -0.3e1 / 0.8e1 * t222 * t245 * t263;
906
907 #ifndef XC_DONT_COMPILE_FXC
908
909 if(order < 2) return;
910
911
912 t267 = t118 / 0.4e1;
913 t269 = t168 / 0.4e1;
914 t270 = t14 * t14;
915 t271 = 0.1e1 / t270;
916 t272 = t271 * t16;
917 t273 = t109 * t109;
918 t274 = t61 * t273;
919 t276 = t7 * t272 * t274;
920 t277 = t276 / 0.6e1;
921 t278 = t14 * t115;
922 t280 = t7 * t278 * t110;
923 t281 = t280 / 0.3e1;
924 t282 = t156 * t109;
925 t284 = t7 * t104 * t282;
926 t285 = t105 * t9;
927 t286 = 0.1e1 / t285;
928 t287 = t8 * t286;
929 t288 = -t106 + t287;
930 t289 = t61 * t288;
931 t291 = t7 * t104 * t289;
932 t292 = t291 / 0.2e1;
933 t294 = 0.1e1 / t114 / t9;
934 t295 = t15 * t294;
935 t297 = t7 * t295 * t61;
936 t298 = t297 / 0.12e2;
937 t300 = t7 * t116 * t156;
938 t301 = t300 / 0.4e1;
939 t303 = 0.1e1 / t27 / t49;
940 t304 = sigma[0] * t303;
941 t308 = t49 * t120;
942 t310 = 0.1e1 / t26 / t308;
943 t315 = params->B * t143;
944 t316 = t315 * t145;
945 t317 = t147 * t25;
946 t318 = 0.1e1 / t317;
947 t320 = 0.1e1 / t133 / t36;
948 t322 = params->C * params->C;
949 t333 = t47 * t47;
950 t334 = t333 * sigma[0];
951 t337 = 0.1e1 / t26 / t147 / t308;
952 t341 = 0.1e1 / t150 / t56;
953 t342 = params->E * params->E;
954 t344 = t43 * t46;
955 t345 = t341 * t342 * t344;
956 t348 = 0.11e2 / 0.27e2 * t24 * t304 * t37 - t128 * t47 * t310 * t135 / 0.24e2 + t316 * t318 * t320 * t322 / 0.162e3 - 0.11e2 / 0.27e2 * t42 * t304 * t57 + 0.35e2 / 0.1296e4 * t146 * t318 * t151 * params->E - t144 * t334 * t337 * t345 / 0.23328e5;
957 t350 = t7 * t17 * t348;
958 t351 = 0.3e1 / 0.8e1 * t350;
959 t352 = t66 * t66;
960 t353 = 0.1e1 / t352;
961 t354 = t353 * t16;
962 t355 = t161 * t161;
963 t356 = t98 * t355;
964 t358 = t7 * t354 * t356;
965 t359 = t358 / 0.6e1;
966 t360 = t66 * t115;
967 t362 = t7 * t360 * t162;
968 t363 = t362 / 0.3e1;
969 t364 = -t288;
970 t365 = t98 * t364;
971 t367 = t7 * t160 * t365;
972 t368 = t367 / 0.2e1;
973 t369 = t67 * t294;
974 t371 = t7 * t369 * t98;
975 t372 = t371 / 0.12e2;
976 t373 = -t277 - t281 - t284 - t292 + t298 - t301 - t351 - t359 - t363 - t368 + t372;
977 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
978 v2rho2[0] = -t112 - t267 - 0.3e1 / 0.4e1 * t158 - t164 - t269 + t9 * t373;
979
980 t375 = t6 * t271;
981 t376 = t4 * t375;
982 t377 = t16 * t61;
983 t378 = t173 * t109;
984 t380 = t376 * t377 * t378;
985 t383 = t7 * t278 * t174;
986 t385 = t156 * t173;
987 t387 = t7 * t104 * t385;
988 t390 = 0.1e1 / t114 / t105;
989 t391 = t14 * t390;
990 t392 = t61 * t8;
991 t394 = t7 * t391 * t392;
992 t398 = t6 * t353;
993 t399 = t4 * t398;
994 t400 = t16 * t98;
995 t401 = t178 * t161;
996 t403 = t399 * t400 * t401;
997 t406 = t7 * t360 * t179;
998 t408 = t66 * t390;
999 t409 = t98 * t8;
1000 t411 = t7 * t408 * t409;
1001 t414 = t214 * t161;
1002 t416 = t7 * t160 * t414;
1003 t419 = t7 * t166 * t214;
1004 t421 = -t380 / 0.6e1 - t383 / 0.6e1 - t387 / 0.2e1 - t394 / 0.2e1 - t280 / 0.6e1 + t298 - t300 / 0.8e1 - t403 / 0.6e1 - t406 / 0.6e1 + t411 / 0.2e1 - t362 / 0.6e1 + t372 - t416 / 0.2e1 - t419 / 0.8e1;
1005 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1006 v2rho2[1] = t9 * t421 - t113 - t159 - t165 - t177 - t182 - t217 - t267 - t269;
1007
1008 t424 = t173 * t173;
1009 t425 = t61 * t424;
1010 t427 = t7 * t272 * t425;
1011 t428 = t427 / 0.6e1;
1012 t429 = t383 / 0.3e1;
1013 t430 = t106 + t287;
1014 t431 = t61 * t430;
1015 t433 = t7 * t104 * t431;
1016 t434 = t433 / 0.2e1;
1017 t435 = t178 * t178;
1018 t436 = t98 * t435;
1019 t438 = t7 * t354 * t436;
1020 t439 = t438 / 0.6e1;
1021 t440 = t406 / 0.3e1;
1022 t441 = t214 * t178;
1023 t443 = t7 * t160 * t441;
1024 t444 = -t430;
1025 t445 = t98 * t444;
1026 t447 = t7 * t160 * t445;
1027 t448 = t447 / 0.2e1;
1028 t449 = t419 / 0.4e1;
1029 t451 = 0.1e1 / t71 / t86;
1030 t452 = sigma[2] * t451;
1031 t456 = t86 * t183;
1032 t458 = 0.1e1 / t70 / t456;
1033 t463 = t315 * t203;
1034 t464 = t205 * t69;
1035 t465 = 0.1e1 / t464;
1036 t467 = 0.1e1 / t194 / t79;
1037 t479 = t84 * t84;
1038 t480 = t479 * sigma[2];
1039 t483 = 0.1e1 / t70 / t205 / t456;
1040 t487 = 0.1e1 / t208 / t93;
1041 t489 = t487 * t342 * t344;
1042 t492 = 0.11e2 / 0.27e2 * t24 * t452 * t80 - t128 * t84 * t458 * t196 / 0.24e2 + t463 * t465 * t467 * t322 / 0.162e3 - 0.11e2 / 0.27e2 * t42 * t452 * t94 + 0.35e2 / 0.1296e4 * t204 * t465 * t209 * params->E - t144 * t480 * t483 * t489 / 0.23328e5;
1043 t494 = t7 * t68 * t492;
1044 t495 = 0.3e1 / 0.8e1 * t494;
1045 t496 = -t428 - t429 - t434 + t298 - t439 - t440 - t443 - t448 + t372 - t449 - t495;
1046 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1047 v2rho2[2] = -t176 - t267 - t181 - t269 - 0.3e1 / 0.4e1 * t216 + t9 * t496;
1048
1049 t499 = t7 * t17 * t241;
1050 t501 = t6 * t14;
1051 t502 = t241 * t109;
1052 t506 = t23 * t122;
1053 t510 = t131 * t134;
1054 t511 = params->C * sigma[0];
1055 t517 = t149 * t320 * t322;
1056 t523 = t144 * t149;
1057 t524 = t151 * params->E;
1058 t525 = t524 * t47;
1059 t528 = t147 * t129;
1060 t530 = 0.1e1 / t26 / t528;
1061 t535 = -t19 * t506 * t37 / 0.9e1 + t128 * t510 * t511 / 0.72e2 - t315 * t47 * t517 / 0.432e3 + t41 * t506 * t57 / 0.9e1 - t523 * t525 / 0.108e3 + t144 * t333 * t530 * t345 / 0.62208e5;
1062 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1063 v2rhosigma[0] = -t499 / 0.2e1 - t222 * t501 * t502 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t223 * t535;
1064
1065 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1066 v2rhosigma[1] = 0.0e0;
1067
1068 t540 = t7 * t68 * t263;
1069 t541 = t6 * t66;
1070 t542 = t263 * t161;
1071 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1072 v2rhosigma[2] = -t222 * t541 * t542 / 0.2e1 - t540 / 0.2e1;
1073
1074 t547 = t501 * t241 * t173;
1075 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1076 v2rhosigma[3] = -t222 * t547 / 0.2e1 - t499 / 0.2e1;
1077
1078 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1079 v2rhosigma[4] = 0.0e0;
1080
1081 t552 = t541 * t263 * t178;
1082 t555 = t23 * t185;
1083 t559 = t192 * t195;
1084 t560 = params->C * sigma[2];
1085 t566 = t207 * t467 * t322;
1086 t572 = t144 * t207;
1087 t573 = t209 * params->E;
1088 t574 = t573 * t84;
1089 t577 = t205 * t190;
1090 t579 = 0.1e1 / t70 / t577;
1091 t584 = -t19 * t555 * t80 / 0.9e1 + t128 * t559 * t560 / 0.72e2 - t315 * t84 * t566 / 0.432e3 + t41 * t555 * t94 / 0.9e1 - t572 * t574 / 0.108e3 + t144 * t479 * t579 * t489 / 0.62208e5;
1092 t585 = t245 * t584;
1093 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1094 v2rhosigma[5] = -t540 / 0.2e1 - t222 * t552 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t585;
1095
1096 t594 = t236 * t320 * t322;
1097 t598 = t524 * sigma[0];
1098 t601 = t147 * t50;
1099 t603 = 0.1e1 / t26 / t601;
1100 t608 = -t128 * t52 * t134 * params->C / 0.288e3 + t315 * sigma[0] * t594 / 0.1152e4 + t144 * t236 * t598 / 0.384e3 - t144 * t145 * t603 * t345 / 0.165888e6;
1101 t609 = t223 * t608;
1102 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1103 v2sigma2[0] = -0.3e1 / 0.8e1 * t222 * t609;
1104
1105 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1106 v2sigma2[1] = 0.0e0;
1107
1108 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1109 v2sigma2[2] = 0.0e0;
1110
1111 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1112 v2sigma2[3] = 0.0e0;
1113
1114 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1115 v2sigma2[4] = 0.0e0;
1116
1117 t618 = t258 * t467 * t322;
1118 t622 = t573 * sigma[2];
1119 t625 = t205 * t87;
1120 t627 = 0.1e1 / t70 / t625;
1121 t632 = -t128 * t89 * t195 * params->C / 0.288e3 + t315 * sigma[2] * t618 / 0.1152e4 + t144 * t258 * t622 / 0.384e3 - t144 * t203 * t627 * t489 / 0.165888e6;
1122 t633 = t245 * t632;
1123 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1124 v2sigma2[5] = -0.3e1 / 0.8e1 * t222 * t633;
1125
1126 #ifndef XC_DONT_COMPILE_KXC
1127
1128 if(order < 3) return;
1129
1130
1131 t640 = t297 / 0.4e1;
1132 t642 = t371 / 0.4e1;
1133 t646 = 0.1e1 / t270 / t13;
1134 t647 = t646 * t16;
1135 t648 = t273 * t109;
1136 t649 = t61 * t648;
1137 t651 = t7 * t647 * t649;
1138 t652 = t651 / 0.9e1;
1139 t653 = t271 * t115;
1140 t655 = t7 * t653 * t274;
1141 t657 = t109 * t288;
1142 t659 = t376 * t377 * t657;
1143 t660 = t659 / 0.2e1;
1144 t661 = t14 * t294;
1145 t663 = t7 * t661 * t110;
1146 t664 = t663 / 0.3e1;
1147 t666 = 0.1e1 / t352 / t65;
1148 t667 = t666 * t16;
1149 t668 = t355 * t161;
1150 t669 = t98 * t668;
1151 t671 = t7 * t667 * t669;
1152 t672 = t671 / 0.9e1;
1153 t673 = t353 * t115;
1154 t675 = t7 * t673 * t356;
1155 t677 = t161 * t364;
1156 t679 = t399 * t400 * t677;
1157 t680 = t679 / 0.2e1;
1158 t681 = t66 * t294;
1159 t683 = t7 * t681 * t162;
1160 t684 = t683 / 0.3e1;
1161 t685 = t156 * t273;
1162 t687 = t7 * t272 * t685;
1163 t688 = t687 / 0.2e1;
1164 t690 = t7 * t278 * t282;
1165 t693 = t7 * t278 * t289;
1166 t695 = t348 * t109;
1167 t697 = t7 * t104 * t695;
1168 t698 = 0.3e1 / 0.2e1 * t697;
1169 t699 = t156 * t288;
1170 t701 = t7 * t104 * t699;
1171 t702 = 0.3e1 / 0.2e1 * t701;
1172 t703 = t105 * t105;
1173 t704 = 0.1e1 / t703;
1174 t705 = t8 * t704;
1175 t707 = 0.3e1 * t286 - 0.3e1 * t705;
1176 t708 = t61 * t707;
1177 t710 = t7 * t104 * t708;
1178 t711 = t710 / 0.2e1;
1179 t712 = t15 * t390;
1180 t714 = t7 * t712 * t61;
1181 t715 = 0.5e1 / 0.36e2 * t714;
1182 t717 = t7 * t360 * t365;
1183 t719 = -t707;
1184 t720 = t98 * t719;
1185 t722 = t7 * t160 * t720;
1186 t723 = t722 / 0.2e1;
1187 t724 = t67 * t390;
1188 t726 = t7 * t724 * t98;
1189 t727 = 0.5e1 / 0.36e2 * t726;
1190 t729 = t7 * t295 * t156;
1191 t730 = t729 / 0.4e1;
1192 t732 = t7 * t116 * t348;
1193 t735 = 0.1e1 / t27 / t50;
1194 t736 = sigma[0] * t735;
1195 t741 = 0.1e1 / t26 / t147;
1196 t746 = t147 * t120;
1197 t747 = 0.1e1 / t746;
1198 t753 = 0.1e1 / t27 / t601;
1199 t756 = t133 * t133;
1200 t757 = 0.1e1 / t756;
1201 t758 = t322 * params->C;
1202 t760 = t18 * t23;
1203 t761 = t757 * t758 * t760;
1204 t771 = t147 * t147;
1205 t773 = 0.1e1 / t26 / t771;
1206 t778 = t333 * t145;
1207 t781 = 0.1e1 / t27 / t771 / t50;
1208 t784 = t150 * t150;
1209 t785 = 0.1e1 / t784;
1210 t786 = t342 * params->E;
1211 t790 = t18 / t22 / t142;
1212 t791 = t785 * t786 * t790;
1213 t794 = -0.154e3 / 0.81e2 * t24 * t736 * t37 + 0.341e3 / 0.972e3 * t128 * t47 * t741 * t135 - 0.19e2 / 0.162e3 * t316 * t747 * t320 * t322 + t315 * t333 * t753 * t761 / 0.486e3 + 0.154e3 / 0.81e2 * t42 * t736 * t57 - 0.569e3 / 0.1944e4 * t146 * t747 * t151 * params->E + t144 * t334 * t773 * t345 / 0.864e3 - t144 * t778 * t781 * t791 / 0.139968e6;
1214 t796 = t7 * t17 * t794;
1215 t797 = 0.3e1 / 0.8e1 * t796;
1216 t798 = -t693 / 0.2e1 - t698 - t702 - t711 - t715 - t717 / 0.2e1 - t723 - t727 + t730 - 0.3e1 / 0.8e1 * t732 - t797;
1217 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1218 v3rho3[0] = -t276 / 0.2e1 - t280 - t358 / 0.2e1 - t362 - 0.3e1 * t284 - 0.3e1 / 0.2e1 * t291 + t640 - 0.3e1 / 0.2e1 * t367 + t642 - 0.3e1 / 0.4e1 * t300 - 0.9e1 / 0.8e1 * t350 + t9 * (t652 - t655 / 0.6e1 - t660 + t664 + t672 - t675 / 0.6e1 - t680 + t684 - t688 - t690 + t798);
1219
1220 t803 = t380 / 0.3e1;
1221 t804 = t403 / 0.3e1;
1222 t808 = t7 * t408 * t98;
1223 t811 = t7 * t369 * t214;
1224 t814 = t7 * t391 * t61;
1225 t826 = t7 * t661 * t174;
1226 t829 = 0.1e1 / t114 / t285;
1227 t830 = t14 * t829;
1228 t832 = t7 * t830 * t392;
1229 t834 = t808 / 0.2e1 + t811 / 0.12e2 - t814 / 0.2e1 - t715 - t727 + t729 / 0.6e1 - t732 / 0.8e1 - t655 / 0.18e2 + 0.2e1 / 0.9e1 * t663 - t675 / 0.18e2 + 0.2e1 / 0.9e1 * t683 - t690 / 0.3e1 - t693 / 0.6e1 - t717 / 0.6e1 + t826 / 0.9e1 + 0.7e1 / 0.6e1 * t832;
1230 t835 = t156 * t8;
1231 t837 = t7 * t391 * t835;
1232 t839 = t7 * t681 * t179;
1233 t841 = t214 * t355;
1234 t843 = t7 * t354 * t841;
1235 t846 = t7 * t360 * t414;
1236 t847 = t846 / 0.3e1;
1237 t849 = t7 * t278 * t385;
1238 t850 = t849 / 0.3e1;
1239 t851 = t348 * t173;
1240 t853 = t7 * t104 * t851;
1241 t855 = t66 * t829;
1242 t857 = t7 * t855 * t409;
1243 t859 = t214 * t364;
1244 t861 = t7 * t160 * t859;
1245 t863 = t16 * t156;
1246 t865 = t376 * t863 * t378;
1247 t867 = t173 * t288;
1248 t869 = t376 * t377 * t867;
1249 t871 = t178 * t364;
1250 t873 = t399 * t400 * t871;
1251 t875 = t390 * t98;
1252 t876 = t8 * t161;
1253 t878 = t399 * t875 * t876;
1254 t880 = t6 * t646;
1255 t881 = t4 * t880;
1256 t882 = t173 * t273;
1257 t884 = t881 * t377 * t882;
1258 t886 = t115 * t61;
1259 t888 = t376 * t886 * t378;
1260 t889 = t888 / 0.9e1;
1261 t890 = t390 * t61;
1262 t891 = t8 * t109;
1263 t893 = t376 * t890 * t891;
1264 t895 = t6 * t666;
1265 t896 = t4 * t895;
1266 t897 = t178 * t355;
1267 t899 = t896 * t400 * t897;
1268 t901 = t115 * t98;
1269 t903 = t399 * t901 * t401;
1270 t904 = t903 / 0.9e1;
1271 t905 = -t837 + t839 / 0.9e1 - t843 / 0.6e1 - t847 - t850 - t853 / 0.2e1 - 0.7e1 / 0.6e1 * t857 - t861 / 0.2e1 - t865 / 0.3e1 - t869 / 0.6e1 - t873 / 0.6e1 + t878 / 0.3e1 + t884 / 0.9e1 - t889 - t893 / 0.3e1 + t899 / 0.9e1 - t904;
1272 t908 = -t292 + t640 - t368 + t642 - t387 + t411 - t416 - t449 - t300 / 0.2e1 - t351 + t9 * (t834 + t905);
1273 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1274 v3rho3[1] = -t277 - 0.2e1 / 0.3e1 * t280 - t359 - 0.2e1 / 0.3e1 * t362 - t803 - t429 - t394 - t804 - t440 - t284 + t908;
1275
1276 t914 = t7 * t166 * t492;
1277 t925 = t7 * t653 * t425;
1278 t928 = t7 * t278 * t431;
1279 t931 = t7 * t673 * t436;
1280 t933 = -t914 / 0.8e1 + t811 / 0.6e1 - t715 - t727 + t729 / 0.12e2 + t663 / 0.9e1 + t683 / 0.9e1 + 0.2e1 / 0.9e1 * t826 - t832 / 0.3e1 + 0.2e1 / 0.9e1 * t839 - t847 - t850 + t857 / 0.3e1 - t925 / 0.18e2 - t928 / 0.6e1 - t931 / 0.18e2;
1281 t935 = t7 * t360 * t441;
1282 t938 = t7 * t360 * t445;
1283 t940 = t156 * t424;
1284 t942 = t7 * t272 * t940;
1285 t944 = t156 * t430;
1286 t946 = t7 * t104 * t944;
1287 t949 = -t286 - 0.3e1 * t705;
1288 t950 = t61 * t949;
1289 t952 = t7 * t104 * t950;
1290 t954 = t214 * t8;
1291 t956 = t7 * t408 * t954;
1292 t957 = -t949;
1293 t958 = t98 * t957;
1294 t960 = t7 * t160 * t958;
1295 t962 = t492 * t161;
1296 t964 = t7 * t160 * t962;
1297 t966 = t424 * t109;
1298 t968 = t881 * t377 * t966;
1299 t970 = t430 * t109;
1300 t972 = t376 * t377 * t970;
1301 t974 = t435 * t161;
1302 t976 = t896 * t400 * t974;
1303 t978 = t178 * t8;
1304 t980 = t399 * t875 * t978;
1305 t982 = t16 * t214;
1306 t984 = t399 * t982 * t401;
1307 t986 = t444 * t161;
1308 t988 = t399 * t400 * t986;
1309 t990 = t173 * t8;
1310 t992 = t376 * t890 * t990;
1311 t994 = -t935 / 0.3e1 - t938 / 0.6e1 - t942 / 0.6e1 - t946 / 0.2e1 - t952 / 0.2e1 + t956 - t960 / 0.2e1 - t964 / 0.2e1 - t889 - t904 + t968 / 0.9e1 - t972 / 0.6e1 + t976 / 0.9e1 + t980 / 0.3e1 - t984 / 0.3e1 - t988 / 0.6e1 - t992 / 0.3e1;
1312 t997 = t642 - t387 + t411 - t416 - t419 / 0.2e1 - t434 - t443 - t448 - t301 - t495 + t9 * (t933 + t994);
1313 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1314 v3rho3[2] = -t281 - t363 - t803 - 0.2e1 / 0.3e1 * t383 - t394 - t804 - 0.2e1 / 0.3e1 * t406 - t428 - t439 + t640 + t997;
1315
1316 t1005 = t826 / 0.3e1;
1317 t1006 = t839 / 0.3e1;
1318 t1011 = t424 * t173;
1319 t1012 = t61 * t1011;
1320 t1014 = t7 * t647 * t1012;
1321 t1015 = t1014 / 0.9e1;
1322 t1016 = t173 * t430;
1323 t1018 = t376 * t377 * t1016;
1324 t1019 = t1018 / 0.2e1;
1325 t1020 = t435 * t178;
1326 t1021 = t98 * t1020;
1327 t1023 = t7 * t667 * t1021;
1328 t1024 = t1023 / 0.9e1;
1329 t1026 = t178 * t444;
1330 t1028 = t399 * t400 * t1026;
1331 t1029 = t1028 / 0.2e1;
1332 t1030 = t811 / 0.4e1;
1333 t1033 = -0.3e1 * t286 - 0.3e1 * t705;
1334 t1034 = t61 * t1033;
1335 t1036 = t7 * t104 * t1034;
1336 t1037 = t1036 / 0.2e1;
1337 t1038 = t214 * t435;
1338 t1040 = t7 * t354 * t1038;
1339 t1041 = t1040 / 0.2e1;
1340 t1042 = t492 * t178;
1341 t1044 = t7 * t160 * t1042;
1342 t1045 = 0.3e1 / 0.2e1 * t1044;
1343 t1046 = t214 * t444;
1344 t1048 = t7 * t160 * t1046;
1345 t1049 = 0.3e1 / 0.2e1 * t1048;
1346 t1050 = -t1033;
1347 t1051 = t98 * t1050;
1348 t1053 = t7 * t160 * t1051;
1349 t1054 = t1053 / 0.2e1;
1350 t1056 = 0.1e1 / t71 / t87;
1351 t1057 = sigma[2] * t1056;
1352 t1062 = 0.1e1 / t70 / t205;
1353 t1067 = t205 * t183;
1354 t1068 = 0.1e1 / t1067;
1355 t1074 = 0.1e1 / t71 / t625;
1356 t1077 = t194 * t194;
1357 t1078 = 0.1e1 / t1077;
1358 t1080 = t1078 * t758 * t760;
1359 t1090 = t205 * t205;
1360 t1092 = 0.1e1 / t70 / t1090;
1361 t1097 = t479 * t203;
1362 t1100 = 0.1e1 / t71 / t1090 / t87;
1363 t1103 = t208 * t208;
1364 t1104 = 0.1e1 / t1103;
1365 t1106 = t1104 * t786 * t790;
1366 t1109 = -0.154e3 / 0.81e2 * t24 * t1057 * t80 + 0.341e3 / 0.972e3 * t128 * t84 * t1062 * t196 - 0.19e2 / 0.162e3 * t463 * t1068 * t467 * t322 + t315 * t479 * t1074 * t1080 / 0.486e3 + 0.154e3 / 0.81e2 * t42 * t1057 * t94 - 0.569e3 / 0.1944e4 * t204 * t1068 * t209 * params->E + t144 * t480 * t1092 * t489 / 0.864e3 - t144 * t1097 * t1100 * t1106 / 0.139968e6;
1367 t1111 = t7 * t68 * t1109;
1368 t1112 = 0.3e1 / 0.8e1 * t1111;
1369 t1113 = -t1029 - t715 - t727 + t1030 - 0.3e1 / 0.8e1 * t914 - t1037 - t1041 - t1045 - t1049 - t1054 - t1112;
1370 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1371 v3rho3[3] = -t383 - t406 - t427 / 0.2e1 - t438 / 0.2e1 + t640 + t642 - 0.3e1 / 0.4e1 * t419 - 0.3e1 / 0.2e1 * t433 - 0.3e1 * t443 - 0.3e1 / 0.2e1 * t447 - 0.9e1 / 0.8e1 * t494 + t9 * (t1005 + t1006 - t925 / 0.6e1 - t928 / 0.2e1 - t931 / 0.6e1 - t935 - t938 / 0.2e1 + t1015 - t1019 + t1024 + t1113);
1372
1373 t1117 = t7 * t104 * t502;
1374 t1121 = t7 * t116 * t241 / 0.6e1;
1375 t1123 = t7 * t17 * t535;
1376 t1124 = t241 * t273;
1377 t1128 = t535 * t109;
1378 t1131 = t241 * t288;
1379 t1135 = t23 * t303;
1380 t1139 = t310 * t134;
1381 t1143 = t315 * t318;
1382 t1144 = t320 * t322;
1383 t1145 = t1144 * t47;
1384 t1148 = t147 * t49;
1385 t1150 = 0.1e1 / t27 / t1148;
1386 t1158 = t144 * t318;
1387 t1162 = t144 * t337 * t341;
1388 t1164 = t342 * t333 * t344;
1389 t1167 = t333 * t47;
1390 t1170 = 0.1e1 / t27 / t771 / t49;
1391 t1175 = 0.11e2 / 0.27e2 * t19 * t1135 * t37 - 0.65e2 / 0.648e3 * t128 * t1139 * t511 + 0.17e2 / 0.432e3 * t1143 * t1145 - t315 * t145 * t1150 * t761 / 0.1296e4 - 0.11e2 / 0.27e2 * t41 * t1135 * t57 + 0.29e2 / 0.324e3 * t1158 * t525 - 0.25e2 / 0.62208e5 * t1162 * t1164 + t144 * t1167 * t1170 * t791 / 0.373248e6;
1392 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1393 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t1117 - t1121 - t1123 - t222 * t375 * t1124 / 0.6e1 - t222 * t501 * t1128 - t222 * t501 * t1131 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t223 * t1175;
1394
1395 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1396 v3rho2sigma[1] = 0.0e0;
1397
1398 t1180 = t7 * t160 * t542;
1399 t1184 = t7 * t166 * t263 / 0.6e1;
1400 t1185 = t263 * t355;
1401 t1189 = t263 * t364;
1402 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1403 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t1180 - t1184 - t222 * t398 * t1185 / 0.6e1 - t222 * t541 * t1189 / 0.2e1;
1404
1405 t1195 = t16 * t1;
1406 t1196 = t1195 * t3;
1407 t1197 = t1196 * t547;
1408 t1199 = t3 * t6;
1409 t1200 = t221 * t1199;
1410 t1201 = t271 * t241;
1411 t1202 = t1201 * t378;
1412 t1206 = t501 * t535 * t173;
1413 t1209 = t294 * t1;
1414 t1210 = t1209 * t3;
1415 t1212 = t501 * t241 * t8;
1416 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1417 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t1117 - t1121 - t1123 / 0.2e1 - 0.2e1 / 0.3e1 * t1197 - t1200 * t1202 / 0.6e1 - t222 * t1206 / 0.2e1 - t1210 * t1212 / 0.2e1;
1418
1419 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1420 v3rho2sigma[4] = 0.0e0;
1421
1422 t1216 = t1196 * t552;
1423 t1218 = t353 * t263;
1424 t1219 = t1218 * t401;
1425 t1223 = t541 * t263 * t8;
1426 t1226 = t1196 * t585;
1427 t1229 = t541 * t584 * t161;
1428 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1429 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t1180 - t1184 - 0.2e1 / 0.3e1 * t1216 - t1200 * t1219 / 0.6e1 + t1210 * t1223 / 0.2e1 - t1226 / 0.2e1 - t222 * t1229 / 0.2e1;
1430
1431 t1234 = t375 * t241 * t424;
1432 t1238 = t501 * t241 * t430;
1433 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1434 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t1197 - t1121 - t222 * t1234 / 0.6e1 - t222 * t1238 / 0.2e1;
1435
1436 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1437 v3rho2sigma[7] = 0.0e0;
1438
1439 t1243 = t398 * t263 * t435;
1440 t1247 = t541 * t584 * t178;
1441 t1250 = t541 * t263 * t444;
1442 t1253 = t23 * t451;
1443 t1257 = t458 * t195;
1444 t1261 = t315 * t465;
1445 t1262 = t467 * t322;
1446 t1263 = t1262 * t84;
1447 t1266 = t205 * t86;
1448 t1268 = 0.1e1 / t71 / t1266;
1449 t1276 = t144 * t465;
1450 t1280 = t144 * t483 * t487;
1451 t1282 = t342 * t479 * t344;
1452 t1285 = t479 * t84;
1453 t1288 = 0.1e1 / t71 / t1090 / t86;
1454 t1293 = 0.11e2 / 0.27e2 * t19 * t1253 * t80 - 0.65e2 / 0.648e3 * t128 * t1257 * t560 + 0.17e2 / 0.432e3 * t1261 * t1263 - t315 * t203 * t1268 * t1080 / 0.1296e4 - 0.11e2 / 0.27e2 * t41 * t1253 * t94 + 0.29e2 / 0.324e3 * t1276 * t574 - 0.25e2 / 0.62208e5 * t1280 * t1282 + t144 * t1285 * t1288 * t1106 / 0.373248e6;
1455 t1294 = t245 * t1293;
1456 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1457 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t1216 - t1184 - t1226 - t222 * t1243 / 0.6e1 - t222 * t1247 - t222 * t1250 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t1294;
1458
1459 t1297 = t1196 * t609;
1460 t1300 = t501 * t608 * t109;
1461 t1307 = t1144 * sigma[0];
1462 t1311 = 0.1e1 / t27 / t746;
1463 t1319 = t144 * t530 * t341;
1464 t1321 = t342 * t145 * t344;
1465 t1326 = 0.1e1 / t27 / t771 / t120;
1466 t1331 = t128 * t510 * params->C / 0.54e2 - 0.5e1 / 0.432e3 * t315 * t149 * t1307 + t315 * t47 * t1311 * t761 / 0.3456e4 - t523 * t598 / 0.48e2 + t1319 * t1321 / 0.7776e4 - t144 * t334 * t1326 * t791 / 0.995328e6;
1467 t1332 = t223 * t1331;
1468 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1469 v3rhosigma2[0] = -t1297 / 0.2e1 - t222 * t1300 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t1332;
1470
1471 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1472 v3rhosigma2[1] = 0.0e0;
1473
1474 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1475 v3rhosigma2[2] = 0.0e0;
1476
1477 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1478 v3rhosigma2[3] = 0.0e0;
1479
1480 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1481 v3rhosigma2[4] = 0.0e0;
1482
1483 t1335 = t1196 * t633;
1484 t1337 = t541 * t632 * t161;
1485 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1486 v3rhosigma2[5] = -t222 * t1337 / 0.2e1 - t1335 / 0.2e1;
1487
1488 t1341 = t501 * t608 * t173;
1489 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1490 v3rhosigma2[6] = -t222 * t1341 / 0.2e1 - t1297 / 0.2e1;
1491
1492 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1493 v3rhosigma2[7] = 0.0e0;
1494
1495 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1496 v3rhosigma2[8] = 0.0e0;
1497
1498 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1499 v3rhosigma2[9] = 0.0e0;
1500
1501 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1502 v3rhosigma2[10] = 0.0e0;
1503
1504 t1346 = t541 * t632 * t178;
1505 t1353 = t1262 * sigma[2];
1506 t1357 = 0.1e1 / t71 / t1067;
1507 t1365 = t144 * t579 * t487;
1508 t1367 = t342 * t203 * t344;
1509 t1372 = 0.1e1 / t71 / t1090 / t183;
1510 t1377 = t128 * t559 * params->C / 0.54e2 - 0.5e1 / 0.432e3 * t315 * t207 * t1353 + t315 * t84 * t1357 * t1080 / 0.3456e4 - t572 * t622 / 0.48e2 + t1365 * t1367 / 0.7776e4 - t144 * t480 * t1372 * t1106 / 0.995328e6;
1511 t1378 = t245 * t1377;
1512 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1513 v3rhosigma2[11] = -t1335 / 0.2e1 - t222 * t1346 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t1378;
1514
1515 t1384 = 0.1e1 / t27 / t317;
1516 t1390 = t144 * t603 * t341;
1517 t1392 = t342 * t47 * t344;
1518 t1399 = 0.1e1 / t27 / t771 / t25;
1519 t1404 = t315 * t594 / 0.384e3 - t315 * sigma[0] * t1384 * t761 / 0.9216e4 - t1390 * t1392 / 0.27648e5 + t144 * t238 / 0.384e3 + t144 * t333 * t1399 * t791 / 0.2654208e7;
1520 t1405 = t223 * t1404;
1521 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1522 v3sigma3[0] = -0.3e1 / 0.8e1 * t222 * t1405;
1523
1524 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1525 v3sigma3[1] = 0.0e0;
1526
1527 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1528 v3sigma3[2] = 0.0e0;
1529
1530 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1531 v3sigma3[3] = 0.0e0;
1532
1533 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1534 v3sigma3[4] = 0.0e0;
1535
1536 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1537 v3sigma3[5] = 0.0e0;
1538
1539 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1540 v3sigma3[6] = 0.0e0;
1541
1542 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1543 v3sigma3[7] = 0.0e0;
1544
1545 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1546 v3sigma3[8] = 0.0e0;
1547
1548 t1411 = 0.1e1 / t71 / t464;
1549 t1417 = t144 * t627 * t487;
1550 t1419 = t342 * t84 * t344;
1551 t1426 = 0.1e1 / t71 / t1090 / t69;
1552 t1431 = t315 * t618 / 0.384e3 - t315 * sigma[2] * t1411 * t1080 / 0.9216e4 - t1417 * t1419 / 0.27648e5 + t144 * t260 / 0.384e3 + t144 * t479 * t1426 * t1106 / 0.2654208e7;
1553 t1432 = t245 * t1431;
1554 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1555 v3sigma3[9] = -0.3e1 / 0.8e1 * t222 * t1432;
1556
1557 #ifndef XC_DONT_COMPILE_LXC
1558
1559 if(order < 4) return;
1560
1561
1562 t1446 = 0.4e1 / 0.9e1 * t651 - 0.2e1 / 0.3e1 * t655 - 0.2e1 * t659 + 0.4e1 / 0.3e1 * t663 + 0.4e1 / 0.9e1 * t671 - 0.2e1 / 0.3e1 * t675 - 0.2e1 * t679 + 0.4e1 / 0.3e1 * t683 - 0.2e1 * t687 - 0.4e1 * t690 - 0.2e1 * t693;
1563 t1450 = 0.5e1 / 0.9e1 * t714;
1564 t1453 = 0.5e1 / 0.9e1 * t726;
1565 t1458 = sigma[0] / t27 / t129;
1566 t1468 = 0.1e1 / t1148;
1567 t1482 = t334 / t26 / t771 / rho[0];
1568 t1486 = t322 * t322;
1569 t1488 = 0.1e1 / t756 / t36 * t1486 * t344;
1570 t1508 = t142 * t142;
1571 t1511 = params->D / t1508 / t142;
1572 t1512 = t333 * t333;
1573 t1518 = 0.1e1 / t784 / t56;
1574 t1520 = t342 * t342;
1575 t1529 = t7 * t712 * t156;
1576 t1532 = t7 * t295 * t348;
1577 t1535 = t7 * t116 * t794;
1578 t1542 = t399 * t901 * t677;
1579 t1553 = t376 * t886 * t657;
1580 t1565 = 0.10e2 / 0.27e2 * t7 * t67 * t829 * t98;
1581 t1569 = 0.10e2 / 0.27e2 * t7 * t15 * t829 * t61;
1582 t1570 = t288 * t288;
1583 t1576 = t7 * t391 * t110;
1584 t1578 = t65 * t65;
1585 t1580 = 0.1e1 / t352 / t1578;
1586 t1581 = t1580 * t16;
1587 t1582 = t355 * t355;
1588 t1587 = t666 * t115;
1589 t1589 = t7 * t1587 * t669;
1590 t1591 = t353 * t294;
1591 t1593 = t7 * t1591 * t356;
1592 t1595 = t364 * t364;
1593 t1600 = -0.3e1 / 0.8e1 * t7 * t17 * (0.2618e4 / 0.243e3 * t24 * t1458 * t37 - 0.3047e4 / 0.972e3 * t128 * t47 / t26 / t148 * t135 + 0.2563e4 / 0.1458e4 * t316 * t1468 * t320 * t322 - 0.49e2 / 0.729e3 * t315 * t333 / t27 / t528 * t761 + 0.2e1 / 0.2187e4 * t315 * t1482 * t1488 - 0.2618e4 / 0.243e3 * t42 * t1458 * t57 + 0.19393e5 / 0.5832e4 * t146 * t1468 * t151 * params->E - 0.5107e4 / 0.209952e6 * t144 * t1482 * t345 + 0.73e2 / 0.209952e6 * t144 * t778 / t27 / t771 / t129 * t791 - t1511 * t1512 * sigma[0] / t771 / t1148 * t1518 * t1520 / 0.629856e6) - 0.5e1 / 0.9e1 * t1529 + t1532 / 0.2e1 - t1535 / 0.2e1 + 0.2e1 / 0.3e1 * t896 * t400 * t355 * t364 - 0.2e1 / 0.3e1 * t1542 - 0.2e1 / 0.3e1 * t399 * t400 * t161 * t719 + 0.2e1 / 0.3e1 * t881 * t377 * t273 * t288 - 0.2e1 / 0.3e1 * t1553 - 0.2e1 * t376 * t863 * t657 - 0.2e1 / 0.3e1 * t376 * t377 * t109 * t707 + t1565 + t1569 - t7 * t272 * t61 * t1570 / 0.2e1 - 0.20e2 / 0.27e2 * t1576 - 0.5e1 / 0.27e2 * t7 * t1581 * t98 * t1582 + 0.4e1 / 0.27e2 * t1589 + 0.2e1 / 0.9e1 * t1593 - t7 * t354 * t98 * t1595 / 0.2e1;
1594 t1602 = t7 * t408 * t162;
1595 t1612 = t703 * t9;
1596 t1614 = t8 / t1612;
1597 t1616 = -0.12e2 * t704 + 0.12e2 * t1614;
1598 t1622 = t7 * t360 * t720;
1599 t1629 = t13 * t13;
1600 t1631 = 0.1e1 / t270 / t1629;
1601 t1632 = t1631 * t16;
1602 t1633 = t273 * t273;
1603 t1638 = t646 * t115;
1604 t1640 = t7 * t1638 * t649;
1605 t1642 = t271 * t294;
1606 t1644 = t7 * t1642 * t274;
1607 t1647 = t7 * t661 * t289;
1608 t1650 = t7 * t681 * t365;
1609 t1656 = t7 * t278 * t695;
1610 t1659 = t7 * t278 * t699;
1611 t1662 = t7 * t278 * t708;
1612 t1673 = t7 * t653 * t685;
1613 t1676 = t7 * t661 * t282;
1614 t1678 = -0.20e2 / 0.27e2 * t1602 - 0.3e1 * t7 * t104 * t348 * t288 - 0.2e1 * t7 * t104 * t156 * t707 - t7 * t104 * t61 * t1616 / 0.2e1 - 0.2e1 / 0.3e1 * t1622 + t7 * t160 * t98 * t1616 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1632 * t61 * t1633 + 0.4e1 / 0.27e2 * t1640 + 0.2e1 / 0.9e1 * t1644 + 0.2e1 / 0.3e1 * t1647 + 0.2e1 / 0.3e1 * t1650 - t7 * t272 * t348 * t273 - 0.2e1 * t1656 - 0.2e1 * t1659 - 0.2e1 / 0.3e1 * t1662 - 0.2e1 * t7 * t104 * t794 * t109 + 0.4e1 / 0.9e1 * t7 * t647 * t156 * t648 - 0.2e1 / 0.3e1 * t1673 + 0.4e1 / 0.3e1 * t1676;
1615 t1681 = -0.6e1 * t697 - 0.6e1 * t701 - 0.2e1 * t710 - t1450 - 0.2e1 * t717 - 0.2e1 * t722 - t1453 + t729 - 0.3e1 / 0.2e1 * t732 - 0.3e1 / 0.2e1 * t796 + t9 * (t1600 + t1678);
1616 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1617 v4rho4[0] = t1446 + t1681;
1618
1619 t1683 = t7 * t855 * t98;
1620 t1686 = t7 * t724 * t214;
1621 t1692 = t7 * t278 * t851;
1622 t1698 = t271 * t390;
1623 t1707 = t7 * t681 * t414;
1624 t1708 = t1707 / 0.3e1;
1625 t1710 = 0.1e1 / t114 / t703;
1626 t1713 = t7 * t66 * t1710 * t409;
1627 t1716 = t7 * t360 * t859;
1628 t1719 = t7 * t391 * t174;
1629 t1723 = t7 * t14 * t1710 * t392;
1630 t1725 = t353 * t390;
1631 t1730 = t7 * t408 * t179;
1632 t1737 = t7 * t673 * t841;
1633 t1739 = -0.5e1 / 0.2e1 * t1683 - 0.5e1 / 0.36e2 * t1686 - 0.3e1 / 0.2e1 * t7 * t391 * t156 - t1692 / 0.2e1 - t7 * t104 * t794 * t173 / 0.2e1 - t7 * t1698 * t110 / 0.2e1 - t7 * t160 * t214 * t719 / 0.2e1 + t1708 + 0.25e2 / 0.6e1 * t1713 - t1716 / 0.2e1 - 0.5e1 / 0.27e2 * t1719 - 0.25e2 / 0.6e1 * t1723 + t7 * t1725 * t162 / 0.2e1 - 0.5e1 / 0.27e2 * t1730 + t7 * t667 * t214 * t668 / 0.9e1 - t1737 / 0.6e1;
1634 t1741 = t7 * t661 * t385;
1635 t1742 = t1741 / 0.3e1;
1636 t1744 = t7 * t830 * t835;
1637 t1751 = t7 * t830 * t61;
1638 t1769 = t376 * t294 * t61 * t378;
1639 t1770 = t1769 / 0.9e1;
1640 t1771 = t829 * t61;
1641 t1773 = t376 * t1771 * t891;
1642 t1785 = t399 * t294 * t98 * t401;
1643 t1786 = t1785 / 0.9e1;
1644 t1790 = t1742 + 0.7e1 / 0.2e1 * t1744 - 0.3e1 / 0.2e1 * t7 * t391 * t348 * t8 + 0.5e1 / 0.2e1 * t1751 - 0.5e1 / 0.12e2 * t1529 + t1532 / 0.4e1 - t1535 / 0.8e1 - t376 * t16 * t348 * t378 / 0.2e1 - t376 * t863 * t867 / 0.2e1 - t376 * t377 * t173 * t707 / 0.6e1 + t1770 + 0.7e1 / 0.6e1 * t1773 - t399 * t400 * t178 * t719 / 0.6e1 + t399 * t875 * t8 * t364 / 0.2e1 + t1786 - t399 * t982 * t677 / 0.2e1;
1645 t1797 = t376 * t115 * t156 * t378;
1646 t1799 = t390 * t156;
1647 t1803 = t376 * t886 * t867;
1648 t1810 = t399 * t901 * t871;
1649 t1816 = t829 * t98;
1650 t1818 = t399 * t1816 * t876;
1651 t1821 = t4 * t6 * t1631;
1652 t1827 = t881 * t886 * t882;
1653 t1834 = t4 * t6 * t1580;
1654 t1840 = t896 * t901 * t897;
1655 t1844 = t881 * t863 * t882 / 0.3e1 - t1797 / 0.3e1 - t376 * t1799 * t891 - t1803 / 0.6e1 - t376 * t890 * t8 * t288 / 0.2e1 - t1810 / 0.6e1 - t896 * t875 * t8 * t355 / 0.3e1 - 0.7e1 / 0.6e1 * t1818 - 0.5e1 / 0.27e2 * t1821 * t377 * t173 * t648 + t1827 / 0.9e1 + t881 * t890 * t8 * t273 / 0.3e1 - 0.5e1 / 0.27e2 * t1834 * t400 * t178 * t668 + t1840 / 0.9e1 - t1542 / 0.6e1 - t1553 / 0.6e1 + t1565;
1656 t1867 = t1569 - 0.5e1 / 0.9e1 * t1576 + t1589 / 0.27e2 + t1593 / 0.9e1 - 0.5e1 / 0.9e1 * t1602 - t1622 / 0.6e1 + t1640 / 0.27e2 + t1644 / 0.9e1 + t1647 / 0.3e1 + t1650 / 0.3e1 - t1656 / 0.2e1 - t1659 / 0.2e1 - t1662 / 0.6e1 - t1673 / 0.6e1 + 0.2e1 / 0.3e1 * t1676 + t881 * t377 * t867 * t109 / 0.3e1 + t896 * t400 * t871 * t161 / 0.3e1;
1657 t1876 = t9 * (t1739 + t1790 + t1844 + t1867) + 0.3e1 / 0.2e1 * t808 + t1030 - 0.3e1 / 0.2e1 * t814 - t797 - t1450 - t1453 + 0.3e1 / 0.4e1 * t729 - 0.3e1 / 0.4e1 * t732 + t652 - t655 / 0.3e1;
1658 t1879 = t663 + t672 - t675 / 0.3e1 + t683 - t688 - 0.2e1 * t690 - t693 - t698 - t702 - t711 - t717;
1659 t1887 = -t723 + t1005 + 0.7e1 / 0.2e1 * t832 - 0.3e1 * t837 + t1006 - t843 / 0.2e1 - t846 - t849 - 0.3e1 / 0.2e1 * t853 - 0.7e1 / 0.2e1 * t857 - 0.3e1 / 0.2e1 * t861;
1660 t1891 = t888 / 0.3e1;
1661 t1893 = t903 / 0.3e1;
1662 t1894 = -t660 - t680 - t865 - t869 / 0.2e1 - t873 / 0.2e1 + t878 + t884 / 0.3e1 - t1891 - t893 + t899 / 0.3e1 - t1893;
1663 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1664 v4rho4[1] = t1876 + t1879 + t1887 + t1894;
1665
1666 t1899 = t7 * t369 * t492;
1667 t1902 = t7 * t855 * t954;
1668 t1905 = t7 * t360 * t958;
1669 t1923 = t7 * t653 * t940;
1670 t1926 = t7 * t278 * t944;
1671 t1929 = t7 * t278 * t950;
1672 t1932 = 0.1e1 / t114 / t1612;
1673 t1934 = t8 * t8;
1674 t1940 = t7 * t681 * t445;
1675 t1947 = t7 * t360 * t962;
1676 t1953 = t7 * t661 * t431;
1677 t1956 = t7 * t1591 * t436;
1678 t1962 = t7 * t681 * t441;
1679 t1964 = t7 * t408 * t214 + t1899 / 0.12e2 - 0.7e1 / 0.3e1 * t1902 - t1905 / 0.3e1 - t7 * t272 * t348 * t424 / 0.6e1 - t7 * t104 * t348 * t430 / 0.2e1 - t7 * t104 * t156 * t949 - t7 * t160 * t492 * t364 / 0.2e1 - t1923 / 0.9e1 - t1926 / 0.3e1 - t1929 / 0.3e1 - t7 * t353 * t1932 * t98 * t1934 / 0.3e1 + t1940 / 0.9e1 - t7 * t354 * t492 * t355 / 0.6e1 - t1947 / 0.3e1 - t7 * t1698 * t174 / 0.3e1 + t1953 / 0.9e1 + t1956 / 0.27e2 + t7 * t1725 * t179 / 0.3e1 + 0.2e1 / 0.9e1 * t1962;
1680 t1966 = t7 * t1642 * t425;
1681 t1998 = t1966 / 0.27e2 + t1683 / 0.3e1 - 0.5e1 / 0.18e2 * t1686 - t7 * t271 * t1932 * t61 * t1934 / 0.3e1 - t1692 / 0.3e1 + 0.4e1 / 0.9e1 * t1707 + 0.41e2 / 0.9e1 * t1713 - t1716 / 0.3e1 - 0.10e2 / 0.27e2 * t1719 - 0.41e2 / 0.9e1 * t1723 - 0.10e2 / 0.27e2 * t1730 - t1737 / 0.9e1 + 0.4e1 / 0.9e1 * t1741 - 0.2e1 / 0.3e1 * t1744 - t1751 / 0.3e1 - 0.5e1 / 0.18e2 * t1529 + t1532 / 0.12e2 + 0.2e1 / 0.9e1 * t881 * t863 * t966 + t881 * t377 * t424 * t288 / 0.9e1 - t376 * t863 * t970 / 0.3e1;
1682 t2013 = t399 * t1816 * t978;
1683 t2018 = t390 * t214;
1684 t2031 = t881 * t886 * t966;
1685 t2034 = t376 * t1771 * t990;
1686 t2041 = t376 * t886 * t970;
1687 t2048 = t896 * t901 * t974;
1688 t2055 = t399 * t115 * t214 * t401;
1689 t2062 = t399 * t901 * t986;
1690 t2071 = -t376 * t377 * t949 * t109 / 0.3e1 - t376 * t377 * t430 * t288 / 0.6e1 + t896 * t400 * t435 * t364 / 0.9e1 - 0.7e1 / 0.9e1 * t2013 - t399 * t982 * t871 / 0.3e1 + 0.2e1 / 0.3e1 * t399 * t2018 * t876 - t399 * t400 * t957 * t161 / 0.3e1 - t399 * t400 * t444 * t364 / 0.6e1 + 0.2e1 / 0.27e2 * t2031 + 0.7e1 / 0.9e1 * t2034 + t881 * t377 * t430 * t273 / 0.9e1 - t2041 / 0.9e1 - 0.5e1 / 0.27e2 * t1834 * t400 * t435 * t355 + 0.2e1 / 0.27e2 * t2048 + 0.2e1 / 0.9e1 * t896 * t982 * t897 - 0.2e1 / 0.9e1 * t2055 + t896 * t400 * t444 * t355 / 0.9e1 - t2062 / 0.9e1 - 0.2e1 / 0.3e1 * t376 * t1799 * t990 - 0.5e1 / 0.27e2 * t1821 * t377 * t424 * t273;
1691 t2096 = 0.4e1 / 0.27e2 * t1769 - 0.2e1 / 0.9e1 * t1773 + 0.4e1 / 0.27e2 * t1785 - 0.2e1 / 0.9e1 * t1797 - t1803 / 0.9e1 - t1810 / 0.9e1 + 0.2e1 / 0.9e1 * t1818 + 0.2e1 / 0.27e2 * t1827 + 0.2e1 / 0.27e2 * t1840 + t1565 + t1569 - 0.10e2 / 0.27e2 * t1576 + t1593 / 0.27e2 - 0.10e2 / 0.27e2 * t1602 + t1644 / 0.27e2 + t1647 / 0.9e1 + t1650 / 0.9e1 + 0.2e1 / 0.9e1 * t1676 - 0.4e1 / 0.9e1 * t896 * t875 * t401 * t8 + 0.4e1 / 0.9e1 * t881 * t890 * t378 * t8;
1692 t2108 = t9 * (t1964 + t1998 + t2071 + t2096) - t914 / 0.4e1 + t808 + t811 / 0.2e1 - t814 - t1450 - t1453 + t729 / 0.2e1 - t732 / 0.4e1 - t655 / 0.9e1 + 0.2e1 / 0.3e1 * t663 - t675 / 0.9e1 + 0.2e1 / 0.3e1 * t683;
1693 t2120 = -0.2e1 / 0.3e1 * t690 - t693 / 0.3e1 - t717 / 0.3e1 + 0.2e1 / 0.3e1 * t826 + 0.5e1 / 0.3e1 * t832 - 0.2e1 * t837 + 0.2e1 / 0.3e1 * t839 - t843 / 0.3e1 - 0.4e1 / 0.3e1 * t846 - 0.4e1 / 0.3e1 * t849 - t853 - 0.5e1 / 0.3e1 * t857 - t861;
1694 t2131 = -t925 / 0.9e1 - t928 / 0.3e1 - t931 / 0.9e1 - 0.2e1 / 0.3e1 * t935 - t938 / 0.3e1 - t942 / 0.3e1 - t946 - t952 + 0.2e1 * t956 - t960 - t964 - 0.2e1 / 0.3e1 * t865 - t869 / 0.3e1;
1695 t2146 = -t873 / 0.3e1 + 0.2e1 / 0.3e1 * t878 + 0.2e1 / 0.9e1 * t884 - 0.4e1 / 0.9e1 * t888 - 0.2e1 / 0.3e1 * t893 + 0.2e1 / 0.9e1 * t899 - 0.4e1 / 0.9e1 * t903 + 0.2e1 / 0.9e1 * t968 - t972 / 0.3e1 + 0.2e1 / 0.9e1 * t976 + 0.2e1 / 0.3e1 * t980 - 0.2e1 / 0.3e1 * t984 - t988 / 0.3e1 - 0.2e1 / 0.3e1 * t992;
1696 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1697 v4rho4[2] = t2108 + t2120 + t2131 + t2146;
1698
1699 t2149 = t7 * t360 * t1046;
1700 t2152 = t7 * t360 * t1051;
1701 t2155 = t7 * t278 * t1034;
1702 t2158 = t7 * t673 * t1038;
1703 t2161 = t7 * t360 * t1042;
1704 t2164 = t7 * t166 * t1109;
1705 t2204 = -t2149 / 0.2e1 - t2152 / 0.6e1 - t2155 / 0.6e1 - t2158 / 0.6e1 - t2161 / 0.2e1 - t2164 / 0.8e1 - t399 * t16 * t492 * t401 / 0.2e1 - t399 * t982 * t986 / 0.2e1 - t399 * t400 * t1050 * t161 / 0.6e1 - t376 * t863 * t1016 / 0.2e1 - t376 * t377 * t173 * t949 / 0.2e1 - 0.5e1 / 0.27e2 * t1834 * t400 * t1020 * t161 - t896 * t875 * t435 * t8 / 0.3e1 - t399 * t400 * t178 * t957 / 0.2e1 + t399 * t875 * t8 * t444 / 0.2e1 - t376 * t377 * t1033 * t109 / 0.6e1;
1706 t2219 = t376 * t886 * t1016;
1707 t2226 = t399 * t901 * t1026;
1708 t2235 = -0.6e1 * t704 - 0.12e2 * t1614;
1709 t2245 = t7 * t1638 * t1012;
1710 t2248 = t7 * t1587 * t1021;
1711 t2267 = t896 * t982 * t974 / 0.3e1 + t399 * t2018 * t978 - 0.5e1 / 0.27e2 * t1821 * t377 * t1011 * t109 + t881 * t890 * t424 * t8 / 0.3e1 - t2219 / 0.6e1 - t376 * t890 * t8 * t430 / 0.2e1 - t2226 / 0.6e1 + t1899 / 0.4e1 - 0.3e1 / 0.2e1 * t7 * t160 * t214 * t957 - t7 * t160 * t98 * t2235 / 0.2e1 - t7 * t160 * t1109 * t161 / 0.2e1 + t2245 / 0.27e2 + t2248 / 0.27e2 + t7 * t647 * t156 * t1011 / 0.9e1 - t7 * t104 * t156 * t1033 / 0.2e1 + t7 * t104 * t61 * t2235 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t408 * t492 * t8;
1712 t2284 = t1902 - t1905 / 0.2e1 - t1923 / 0.6e1 - t1926 / 0.2e1 - t1929 / 0.2e1 + t1940 / 0.3e1 - t1947 / 0.2e1 + t1953 / 0.3e1 + t1956 / 0.9e1 + 0.2e1 / 0.3e1 * t1962 + t1966 / 0.9e1 - 0.5e1 / 0.12e2 * t1686 + t1708 - t1713 / 0.3e1 - 0.5e1 / 0.9e1 * t1719 + t1723 / 0.3e1 - 0.5e1 / 0.9e1 * t1730;
1713 t2303 = t1742 - 0.5e1 / 0.36e2 * t1529 + t2013 / 0.3e1 + t2031 / 0.9e1 - t2034 / 0.3e1 - t2041 / 0.6e1 + t2048 / 0.9e1 - t2055 / 0.3e1 - t2062 / 0.6e1 + t1770 + t1786 + t1565 + t1569 - 0.5e1 / 0.27e2 * t1576 - 0.5e1 / 0.27e2 * t1602 + t881 * t377 * t1016 * t109 / 0.3e1 + t896 * t400 * t1026 * t161 / 0.3e1;
1714 t2309 = t9 * (t2204 + t2267 + t2284 + t2303) - t1045 - t1049 - t1054 + t1024 - t1037 - t1041 + t1015 - t1112 - 0.3e1 / 0.4e1 * t914 + 0.3e1 / 0.4e1 * t811;
1715 t2310 = -t1450 - t1453 + t730 + t664 + t684 + t826 - t832 + t839 - t846 - t849 + t857;
1716 t2321 = -t925 / 0.3e1 - t928 - t931 / 0.3e1 - 0.2e1 * t935 - t938 - t942 / 0.2e1 - 0.3e1 / 0.2e1 * t946 - 0.3e1 / 0.2e1 * t952 + 0.3e1 * t956 - 0.3e1 / 0.2e1 * t960 - 0.3e1 / 0.2e1 * t964;
1717 t2326 = -t1891 - t1893 + t968 / 0.3e1 - t972 / 0.2e1 + t976 / 0.3e1 + t980 - t984 - t988 / 0.2e1 - t992 - t1019 - t1029;
1718 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1719 v4rho4[3] = t2309 + t2310 + t2321 + t2326;
1720
1721 t2339 = 0.4e1 / 0.3e1 * t826 + 0.4e1 / 0.3e1 * t839 - 0.2e1 / 0.3e1 * t925 - 0.2e1 * t928 - 0.2e1 / 0.3e1 * t931 - 0.4e1 * t935 - 0.2e1 * t938 + 0.4e1 / 0.9e1 * t1014 - 0.2e1 * t1018 + 0.4e1 / 0.9e1 * t1023 - 0.2e1 * t1028;
1722 t2347 = t424 * t424;
1723 t2352 = t430 * t430;
1724 t2357 = t435 * t435;
1725 t2362 = t444 * t444;
1726 t2376 = -0.12e2 * t704 - 0.12e2 * t1614;
1727 t2404 = sigma[2] / t71 / t190;
1728 t2414 = 0.1e1 / t1266;
1729 t2428 = t480 / t70 / t1090 / rho[1];
1730 t2433 = 0.1e1 / t1077 / t79 * t1486 * t344;
1731 t2453 = t479 * t479;
1732 t2459 = 0.1e1 / t1103 / t93;
1733 t2473 = -0.5e1 / 0.27e2 * t7 * t1632 * t61 * t2347 - t7 * t272 * t61 * t2352 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1581 * t98 * t2357 - t7 * t354 * t98 * t2362 / 0.2e1 - 0.3e1 * t7 * t160 * t492 * t444 - 0.2e1 * t7 * t160 * t214 * t1050 - t7 * t160 * t98 * t2376 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t667 * t214 * t1020 + t7 * t104 * t61 * t2376 / 0.2e1 - t7 * t354 * t492 * t435 - 0.2e1 * t7 * t160 * t1109 * t178 - 0.2e1 * t2149 - 0.2e1 / 0.3e1 * t2152 - 0.2e1 / 0.3e1 * t2155 - 0.2e1 / 0.3e1 * t2158 - 0.2e1 * t2161 - 0.3e1 / 0.8e1 * t7 * t68 * (0.2618e4 / 0.243e3 * t24 * t2404 * t80 - 0.3047e4 / 0.972e3 * t128 * t84 / t70 / t206 * t196 + 0.2563e4 / 0.1458e4 * t463 * t2414 * t467 * t322 - 0.49e2 / 0.729e3 * t315 * t479 / t71 / t577 * t1080 + 0.2e1 / 0.2187e4 * t315 * t2428 * t2433 - 0.2618e4 / 0.243e3 * t42 * t2404 * t94 + 0.19393e5 / 0.5832e4 * t204 * t2414 * t209 * params->E - 0.5107e4 / 0.209952e6 * t144 * t2428 * t489 + 0.73e2 / 0.209952e6 * t144 * t1097 / t71 / t1090 / t190 * t1106 - t1511 * t2453 * sigma[2] / t1090 / t1266 * t2459 * t1520 / 0.629856e6) - t2164 / 0.2e1 + 0.2e1 / 0.3e1 * t881 * t377 * t424 * t430;
1734 t2502 = -0.2e1 / 0.3e1 * t376 * t377 * t173 * t1033 + 0.2e1 / 0.3e1 * t896 * t400 * t435 * t444 - 0.2e1 * t399 * t982 * t1026 - 0.2e1 / 0.3e1 * t399 * t400 * t178 * t1050 - 0.2e1 / 0.3e1 * t2219 - 0.2e1 / 0.3e1 * t2226 + t1899 / 0.2e1 + 0.4e1 / 0.27e2 * t2245 + 0.4e1 / 0.27e2 * t2248 + 0.2e1 / 0.3e1 * t1940 + 0.2e1 / 0.3e1 * t1953 + 0.2e1 / 0.9e1 * t1956 + 0.4e1 / 0.3e1 * t1962 + 0.2e1 / 0.9e1 * t1966 - 0.5e1 / 0.9e1 * t1686 - 0.20e2 / 0.27e2 * t1719 - 0.20e2 / 0.27e2 * t1730 + t1565 + t1569;
1735 t2505 = -t1450 - t1453 + t811 - 0.3e1 / 0.2e1 * t914 - 0.2e1 * t1036 - 0.2e1 * t1040 - 0.6e1 * t1044 - 0.6e1 * t1048 - 0.2e1 * t1053 - 0.3e1 / 0.2e1 * t1111 + t9 * (t2473 + t2502);
1736 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1737 v4rho4[4] = t2339 + t2505;
1738
1739 t2507 = t7 * t116 * t535;
1740 t2510 = t7 * t17 * t1175;
1741 t2512 = t23 * t735;
1742 t2561 = t7 * t104 * t1128;
1743 t2564 = t7 * t104 * t1131;
1744 t2568 = t7 * t295 * t241 / 0.9e1;
1745 t2586 = t7 * t272 * t1124;
1746 t2589 = t7 * t278 * t502;
1747 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1748 v4rho3sigma[0] = -t2507 / 0.2e1 - 0.3e1 / 0.2e1 * t2510 - 0.3e1 / 0.8e1 * t222 * t223 * (-0.154e3 / 0.81e2 * t19 * t2512 * t37 + 0.253e3 / 0.324e3 * t128 * t741 * t134 * t511 - 0.1025e4 / 0.1944e4 * t315 * t747 * t1145 + 0.89e2 / 0.3888e4 * t315 * t753 * t757 * t758 * t145 * t760 - t315 * t333 * t773 * t1488 / 0.2916e4 + 0.154e3 / 0.81e2 * t41 * t2512 * t57 - 0.223e3 / 0.243e3 * t144 * t747 * t525 + 0.2189e4 / 0.279936e6 * t144 * t773 * t341 * t1164 - 0.137e3 / 0.1119744e7 * t144 * t781 * t785 * t786 * t1167 * t790 + t1511 * t1512 / t771 / t746 * t1518 * t1520 / 0.1679616e7) - 0.4e1 * t2561 - 0.2e1 * t2564 + t2568 - t222 * t375 * t535 * t273 / 0.2e1 - 0.3e1 / 0.2e1 * t222 * t501 * t1175 * t109 - 0.3e1 / 0.2e1 * t222 * t501 * t535 * t288 - t222 * t501 * t241 * t707 / 0.2e1 - 0.2e1 / 0.3e1 * t2586 - 0.2e1 / 0.3e1 * t2589 + t222 * t880 * t241 * t648 / 0.9e1 - t1200 * t1201 * t657 / 0.2e1;
1749
1750 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1751 v4rho3sigma[1] = 0.0e0;
1752
1753 t2599 = t7 * t354 * t1185;
1754 t2602 = t7 * t360 * t542;
1755 t2605 = t7 * t160 * t1189;
1756 t2609 = t7 * t369 * t263 / 0.9e1;
1757 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1758 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2599 - 0.2e1 / 0.3e1 * t2602 - 0.2e1 * t2605 + t2609 + t222 * t895 * t263 * t668 / 0.9e1 - t1200 * t1218 * t677 / 0.2e1 - t222 * t541 * t263 * t719 / 0.2e1;
1759
1760 t2626 = 0.4e1 / 0.3e1 * t1196 * t1206;
1761 t2637 = t115 * t1 * t3;
1762 t2638 = t2637 * t547;
1763 t2640 = t1195 * t1199;
1764 t2642 = 0.4e1 / 0.9e1 * t2640 * t1202;
1765 t2644 = t390 * t1 * t3;
1766 t2645 = t2644 * t1212;
1767 t2657 = t646 * t241;
1768 t2661 = t1209 * t1199;
1769 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1770 v4rho3sigma[3] = -t2507 / 0.3e1 - t2510 / 0.2e1 - 0.4e1 / 0.3e1 * t2561 - 0.2e1 / 0.3e1 * t2564 + t2568 - t2626 - t222 * t501 * t1175 * t173 / 0.2e1 - t1210 * t501 * t241 / 0.2e1 - 0.2e1 / 0.9e1 * t2586 - 0.4e1 / 0.9e1 * t2589 - 0.2e1 / 0.9e1 * t2638 - t2642 + t2645 / 0.6e1 - t1200 * t271 * t535 * t378 / 0.3e1 - t1200 * t1201 * t867 / 0.6e1 - t1210 * t501 * t535 * t8 + t1200 * t2657 * t882 / 0.9e1 - t2661 * t1201 * t891 / 0.3e1;
1771
1772 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1773 v4rho3sigma[4] = 0.0e0;
1774
1775 t2668 = t2637 * t552;
1776 t2671 = 0.4e1 / 0.9e1 * t2640 * t1219;
1777 t2672 = t2644 * t1223;
1778 t2674 = t666 * t263;
1779 t2687 = t2637 * t585;
1780 t2690 = 0.4e1 / 0.3e1 * t1196 * t1229;
1781 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1782 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2599 - 0.4e1 / 0.9e1 * t2602 - 0.2e1 / 0.3e1 * t2605 + t2609 - 0.2e1 / 0.9e1 * t2668 - t2671 - t2672 / 0.6e1 + t1200 * t2674 * t897 / 0.9e1 + t2661 * t1218 * t876 / 0.3e1 - t1200 * t1218 * t871 / 0.6e1 + t1210 * t541 * t263 / 0.2e1 - t2687 / 0.6e1 - t2690 - t222 * t398 * t584 * t355 / 0.6e1 - t222 * t541 * t584 * t364 / 0.2e1;
1783
1784 t2703 = t1196 * t1234;
1785 t2715 = t1196 * t1238;
1786 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1787 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2638 - t2642 - t2626 - 0.4e1 / 0.3e1 * t2645 - 0.2e1 / 0.9e1 * t2589 + t2568 - t2507 / 0.6e1 - 0.2e1 / 0.9e1 * t2703 + t1200 * t2657 * t966 / 0.9e1 - t222 * t375 * t535 * t424 / 0.6e1 - t2661 * t1201 * t990 / 0.3e1 - 0.2e1 / 0.3e1 * t2715 - t1200 * t1201 * t970 / 0.6e1 - t222 * t501 * t535 * t430 / 0.2e1 - t222 * t501 * t241 * t949 / 0.2e1;
1788
1789 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1790 v4rho3sigma[7] = 0.0e0;
1791
1792 t2732 = t1196 * t1243;
1793 t2740 = t1196 * t1247;
1794 t2749 = t1196 * t1250;
1795 t2758 = t1196 * t1294;
1796 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1797 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2668 - t2671 + 0.4e1 / 0.3e1 * t2672 - 0.2e1 / 0.9e1 * t2602 + t2609 - t2687 / 0.3e1 - t2690 - 0.2e1 / 0.9e1 * t2732 + t1200 * t2674 * t974 / 0.9e1 + t2661 * t1218 * t978 / 0.3e1 - 0.4e1 / 0.3e1 * t2740 - t1200 * t353 * t584 * t401 / 0.3e1 + t1210 * t541 * t584 * t8 - 0.2e1 / 0.3e1 * t2749 - t1200 * t1218 * t986 / 0.6e1 - t222 * t541 * t263 * t957 / 0.2e1 - t2758 / 0.2e1 - t222 * t541 * t1293 * t161 / 0.2e1;
1798
1799 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1800 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2638 - 0.2e1 / 0.3e1 * t2703 - 0.2e1 * t2715 + t2568 + t222 * t880 * t241 * t1011 / 0.9e1 - t1200 * t1201 * t1016 / 0.2e1 - t222 * t501 * t241 * t1033 / 0.2e1;
1801
1802 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1803 v4rho3sigma[10] = 0.0e0;
1804
1805 t2778 = t23 * t1056;
1806 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1807 v4rho3sigma[11] = -0.3e1 / 0.8e1 * t222 * t245 * (-0.154e3 / 0.81e2 * t19 * t2778 * t80 + 0.253e3 / 0.324e3 * t128 * t1062 * t195 * t560 - 0.1025e4 / 0.1944e4 * t315 * t1068 * t1263 + 0.89e2 / 0.3888e4 * t315 * t1074 * t1078 * t758 * t203 * t760 - t315 * t479 * t1092 * t2433 / 0.2916e4 + 0.154e3 / 0.81e2 * t41 * t2778 * t94 - 0.223e3 / 0.243e3 * t144 * t1068 * t574 + 0.2189e4 / 0.279936e6 * t144 * t1092 * t487 * t1282 - 0.137e3 / 0.1119744e7 * t144 * t1100 * t1104 * t786 * t1285 * t790 + t1511 * t2453 / t1090 / t1067 * t2459 * t1520 / 0.1679616e7) + t2609 - t2687 / 0.2e1 - 0.3e1 / 0.2e1 * t2758 - t222 * t398 * t584 * t435 / 0.2e1 - 0.3e1 / 0.2e1 * t222 * t541 * t1293 * t178 - 0.3e1 / 0.2e1 * t222 * t541 * t584 * t444 - t222 * t541 * t263 * t1050 / 0.2e1 - 0.2e1 / 0.3e1 * t2668 - 0.2e1 / 0.3e1 * t2732 - 0.4e1 * t2740 - 0.2e1 * t2749 + t222 * t895 * t263 * t1020 / 0.9e1 - t1200 * t1218 * t1026 / 0.2e1;
1808
1809 t2856 = t2637 * t609 / 0.6e1;
1810 t2857 = t1196 * t1300;
1811 t2859 = t1196 * t1332;
1812 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1813 v4rho2sigma2[0] = -t2856 - 0.4e1 / 0.3e1 * t2857 - t2859 - t222 * t375 * t608 * t273 / 0.6e1 - t222 * t501 * t1331 * t109 - t222 * t501 * t608 * t288 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t223 * (-0.19e2 / 0.162e3 * t128 * t1139 * params->C + 0.167e3 / 0.1296e4 * t1143 * t1307 - 0.25e2 / 0.3456e4 * t315 * t1150 * t757 * t758 * t47 * t760 + t315 * t145 * t337 * t1488 / 0.7776e4 + 0.3e1 / 0.16e2 * t1158 * t598 - 0.13e2 / 0.5832e4 * t1162 * t1321 + 0.41e2 / 0.995328e6 * t144 * t1170 * t785 * t786 * t334 * t790 - t1511 * t778 / t771 / t317 * t1518 * t1520 / 0.4478976e7);
1814
1815 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1816 v4rho2sigma2[1] = 0.0e0;
1817
1818 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1819 v4rho2sigma2[2] = 0.0e0;
1820
1821 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1822 v4rho2sigma2[3] = 0.0e0;
1823
1824 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1825 v4rho2sigma2[4] = 0.0e0;
1826
1827 t2908 = t2637 * t633 / 0.6e1;
1828 t2909 = t1196 * t1337;
1829 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1830 v4rho2sigma2[5] = -t2908 - 0.4e1 / 0.3e1 * t2909 - t222 * t398 * t632 * t355 / 0.6e1 - t222 * t541 * t632 * t364 / 0.2e1;
1831
1832 t2921 = t1196 * t1341;
1833 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1834 v4rho2sigma2[6] = -t2856 - 0.2e1 / 0.3e1 * t2857 - t2859 / 0.2e1 - 0.2e1 / 0.3e1 * t2921 - t1200 * t271 * t608 * t378 / 0.6e1 - t222 * t501 * t1331 * t173 / 0.2e1 - t1210 * t501 * t608 * t8 / 0.2e1;
1835
1836 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1837 v4rho2sigma2[7] = 0.0e0;
1838
1839 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1840 v4rho2sigma2[8] = 0.0e0;
1841
1842 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1843 v4rho2sigma2[9] = 0.0e0;
1844
1845 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1846 v4rho2sigma2[10] = 0.0e0;
1847
1848 t2936 = t1196 * t1346;
1849 t2946 = t1196 * t1378;
1850 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1851 v4rho2sigma2[11] = -t2908 - 0.2e1 / 0.3e1 * t2909 - 0.2e1 / 0.3e1 * t2936 - t1200 * t353 * t632 * t401 / 0.6e1 + t1210 * t541 * t632 * t8 / 0.2e1 - t2946 / 0.2e1 - t222 * t541 * t1377 * t161 / 0.2e1;
1852
1853 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1854 v4rho2sigma2[12] = -t2856 - 0.4e1 / 0.3e1 * t2921 - t222 * t375 * t608 * t424 / 0.6e1 - t222 * t501 * t608 * t430 / 0.2e1;
1855
1856 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1857 v4rho2sigma2[13] = 0.0e0;
1858
1859 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1860 v4rho2sigma2[14] = 0.0e0;
1861
1862 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1863 v4rho2sigma2[15] = 0.0e0;
1864
1865 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1866 v4rho2sigma2[16] = 0.0e0;
1867
1868 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1869 v4rho2sigma2[17] = -t2908 - 0.4e1 / 0.3e1 * t2936 - t2946 - t222 * t398 * t632 * t435 / 0.6e1 - t222 * t541 * t1377 * t178 - t222 * t541 * t632 * t444 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t245 * (-0.19e2 / 0.162e3 * t128 * t1257 * params->C + 0.167e3 / 0.1296e4 * t1261 * t1353 - 0.25e2 / 0.3456e4 * t315 * t1268 * t1078 * t758 * t84 * t760 + t315 * t203 * t483 * t2433 / 0.7776e4 + 0.3e1 / 0.16e2 * t1276 * t622 - 0.13e2 / 0.5832e4 * t1280 * t1367 + 0.41e2 / 0.995328e6 * t144 * t1288 * t1104 * t786 * t480 * t790 - t1511 * t1097 / t1090 / t464 * t2459 * t1520 / 0.4478976e7);
1870
1871 t3009 = t1196 * t1405;
1872 t3019 = t758 * t18;
1873 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1874 v4rhosigma3[0] = -t3009 / 0.2e1 - t222 * t501 * t1404 * t109 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t223 * (-t315 * t517 / 0.48e2 + 0.7e1 / 0.3456e4 * t315 * t1311 * t757 * t3019 * t32 - t315 * t47 * t530 * t1488 / 0.20736e5 + 0.11e2 / 0.20736e5 * t1319 * t1392 - 0.13e2 / 0.995328e6 * t144 * t1326 * t785 * t786 * t333 * t790 - t144 * t153 / 0.48e2 + t1511 * t1167 / t771 / t148 * t1518 * t1520 / 0.11943936e8);
1875
1876 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1877 v4rhosigma3[1] = 0.0e0;
1878
1879 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1880 v4rhosigma3[2] = 0.0e0;
1881
1882 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1883 v4rhosigma3[3] = 0.0e0;
1884
1885 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1886 v4rhosigma3[4] = 0.0e0;
1887
1888 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1889 v4rhosigma3[5] = 0.0e0;
1890
1891 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1892 v4rhosigma3[6] = 0.0e0;
1893
1894 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1895 v4rhosigma3[7] = 0.0e0;
1896
1897 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1898 v4rhosigma3[8] = 0.0e0;
1899
1900 t3048 = t1196 * t1432;
1901 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1902 v4rhosigma3[9] = -t222 * t541 * t1431 * t161 / 0.2e1 - t3048 / 0.2e1;
1903
1904 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1905 v4rhosigma3[10] = -t222 * t501 * t1404 * t173 / 0.2e1 - t3009 / 0.2e1;
1906
1907 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1908 v4rhosigma3[11] = 0.0e0;
1909
1910 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1911 v4rhosigma3[12] = 0.0e0;
1912
1913 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1914 v4rhosigma3[13] = 0.0e0;
1915
1916 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1917 v4rhosigma3[14] = 0.0e0;
1918
1919 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1920 v4rhosigma3[15] = 0.0e0;
1921
1922 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1923 v4rhosigma3[16] = 0.0e0;
1924
1925 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1926 v4rhosigma3[17] = 0.0e0;
1927
1928 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1929 v4rhosigma3[18] = 0.0e0;
1930
1931 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1932 v4rhosigma3[19] = -t3048 / 0.2e1 - t222 * t541 * t1431 * t178 / 0.2e1 - 0.3e1 / 0.8e1 * t222 * t245 * (-t315 * t566 / 0.48e2 + 0.7e1 / 0.3456e4 * t315 * t1357 * t1078 * t3019 * t75 - t315 * t84 * t579 * t2433 / 0.20736e5 + 0.11e2 / 0.20736e5 * t1365 * t1419 - 0.13e2 / 0.995328e6 * t144 * t1372 * t1104 * t786 * t479 * t790 - t144 * t211 / 0.48e2 + t1511 * t1285 / t1090 / t206 * t2459 * t1520 / 0.11943936e8);
1933
1934 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1935 v4sigma4[0] = -0.3e1 / 0.8e1 * t222 * t223 * (-t315 * t1384 * t761 / 0.2304e4 + t315 * sigma[0] * t603 * t1488 / 0.55296e5 + 0.5e1 / 0.1327104e7 * t144 * t1399 * t785 * t786 * t145 * t790 - 0.5e1 / 0.55296e5 * t1390 * t342 * sigma[0] * t344 - t1511 * t334 / t771 / t147 * t1518 * t1520 / 0.31850496e8);
1936
1937 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1938 v4sigma4[1] = 0.0e0;
1939
1940 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1941 v4sigma4[2] = 0.0e0;
1942
1943 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1944 v4sigma4[3] = 0.0e0;
1945
1946 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1947 v4sigma4[4] = 0.0e0;
1948
1949 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1950 v4sigma4[5] = 0.0e0;
1951
1952 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1953 v4sigma4[6] = 0.0e0;
1954
1955 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1956 v4sigma4[7] = 0.0e0;
1957
1958 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1959 v4sigma4[8] = 0.0e0;
1960
1961 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1962 v4sigma4[9] = 0.0e0;
1963
1964 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1965 v4sigma4[10] = 0.0e0;
1966
1967 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1968 v4sigma4[11] = 0.0e0;
1969
1970 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1971 v4sigma4[12] = 0.0e0;
1972
1973 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1974 v4sigma4[13] = 0.0e0;
1975
1976 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1977 v4sigma4[14] = -0.3e1 / 0.8e1 * t222 * t245 * (-t315 * t1411 * t1080 / 0.2304e4 + t315 * sigma[2] * t627 * t2433 / 0.55296e5 + 0.5e1 / 0.1327104e7 * t144 * t1426 * t1104 * t786 * t203 * t790 - 0.5e1 / 0.55296e5 * t1417 * t342 * sigma[2] * t344 - t1511 * t480 / t1090 / t205 * t2459 * t1520 / 0.31850496e8);
1978
1979 #ifndef XC_DONT_COMPILE_MXC
1980
1981 if(order < 5) return;
1982
1983
1984 #endif
1985
1986 #endif
1987
1988 #endif
1989
1990 #endif
1991
1992 #endif
1993
1994
1995 }
1996
1997