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_b88.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t3, t5, t6, t7, t8, t9, t10;
24 double t11, t12, t13, t14, t15, t16, t17, t18;
25 double t19, t21, t22, t23, t24, t25, t26, t30;
26 double t31, t33, t34, t35, t39, t41;
27
28 #ifndef XC_DONT_COMPILE_VXC
29 double t45, t46, t47, t49, t50, t54, t55, t56;
30 double t60, t62, t63, t65, t66, t67, t68, t71;
31 double t76, t80, t81, t85, t87, t88, t91, t97;
32
33 #ifndef XC_DONT_COMPILE_FXC
34 double t100, t101, t108, t110, t111, t115, t121, t122;
35 double t123, t131, t134, t135, t138, t139, t142, t145;
36 double t146, t149, t154, t159, t165, t174, t175, t176;
37 double t177, t184, t185, t187, t192, t198, t201, t205;
38 double t212, t215, t216, t219, t222, t223, t226, t232;
39
40 #ifndef XC_DONT_COMPILE_KXC
41 double t236, t237, t248, t249, t253, t258, t267, t268;
42 double t269, t270, t275, t276, t277, t283, t287, t290;
43 double t292, t297, t298, t299, t301, t303, t307, t312;
44 double t317, t328, t332, t340, t348, t349, t352, t353;
45 double t354, t357, t358, t369, t371, t375, t380, t381;
46 double t391, t399, t400, t401, t404, t405, t412, t413;
47 double t414, t422, t425, t427, t431, t437, t440, t444;
48 double t448, t453, t454, t459, t463, t468, t472, t478;
49
50 #ifndef XC_DONT_COMPILE_LXC
51 double t493, t498, t503, t512, t525, t526, t527, t537;
52 double t559, t565, t572, t582, t602, t603, t695, t758;
53 double t773, t801, t882, t896, t905;
54 #endif
55
56 #endif
57
58 #endif
59
60 #endif
61
62 #endif
63
64
65 gga_x_b88_params *params;
66
67 assert(p->params != NULL);
68 params = (gga_x_b88_params * )(p->params);
69
70 t1 = M_CBRT3;
71 t3 = POW_1_3(0.1e1 / M_PI);
72 t5 = M_CBRT4;
73 t6 = t5 * t5;
74 t7 = t1 * t3 * t6;
75 t8 = M_CBRT2;
76 t9 = t8 * t8;
77 t10 = POW_1_3(rho[0]);
78 t11 = t9 * t10;
79 t12 = t1 * t1;
80 t13 = params->beta * t12;
81 t14 = 0.1e1 / t3;
82 t15 = t14 * t5;
83 t16 = t13 * t15;
84 t17 = sigma[0] * t9;
85 t18 = rho[0] * rho[0];
86 t19 = t10 * t10;
87 t21 = 0.1e1 / t19 / t18;
88 t22 = params->gamma * params->beta;
89 t23 = sqrt(sigma[0]);
90 t24 = t22 * t23;
91 t25 = t10 * rho[0];
92 t26 = 0.1e1 / t25;
93 t30 = log(t23 * t8 * t26 + sqrt(POW_2(t23 * t8 * t26) + 0.1e1));
94 t31 = t8 * t26 * t30;
95 t33 = 0.10e1 + t24 * t31;
96 t34 = 0.1e1 / t33;
97 t35 = t21 * t34;
98 t39 = 0.10e1 + 0.2e1 / 0.9e1 * t16 * t17 * t35;
99 t41 = t7 * t11 * t39;
100 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
101 zk[0] = -0.3e1 / 0.16e2 * t41;
102
103 #ifndef XC_DONT_COMPILE_VXC
104
105 if(order < 1) return;
106
107
108 t45 = t25 * t1 * t3;
109 t46 = t6 * t9;
110 t47 = t18 * rho[0];
111 t49 = 0.1e1 / t19 / t47;
112 t50 = t49 * t34;
113 t54 = t33 * t33;
114 t55 = 0.1e1 / t54;
115 t56 = t21 * t55;
116 t60 = t8 / t10 / t18 * t30;
117 t62 = t22 * sigma[0];
118 t63 = t9 * t49;
119 t65 = t17 * t21 + 0.1e1;
120 t66 = sqrt(t65);
121 t67 = 0.1e1 / t66;
122 t68 = t63 * t67;
123 t71 = -0.4e1 / 0.3e1 * t24 * t60 - 0.4e1 / 0.3e1 * t62 * t68;
124 t76 = -0.16e2 / 0.27e2 * t16 * t17 * t50 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t71;
125 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
126 vrho[0] = -t41 / 0.4e1 - 0.3e1 / 0.16e2 * t45 * t46 * t76;
127
128 t80 = t13 * t14;
129 t81 = t5 * t9;
130 t85 = t22 / t23;
131 t87 = t9 * t21;
132 t88 = t87 * t67;
133 t91 = t22 * t88 / 0.2e1 + t85 * t31 / 0.2e1;
134 t97 = t46 * (-0.2e1 / 0.9e1 * t16 * t17 * t56 * t91 + 0.2e1 / 0.9e1 * t80 * t81 * t35);
135 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
136 vsigma[0] = -0.3e1 / 0.16e2 * t45 * t97;
137
138 #ifndef XC_DONT_COMPILE_FXC
139
140 if(order < 2) return;
141
142
143 t100 = 0.1e1 / t19;
144 t101 = t9 * t100;
145 t108 = t18 * t18;
146 t110 = 0.1e1 / t19 / t108;
147 t111 = t110 * t34;
148 t115 = t49 * t55;
149 t121 = 0.1e1 / t54 / t33;
150 t122 = t21 * t121;
151 t123 = t71 * t71;
152 t131 = t8 / t10 / t47 * t30;
153 t134 = t9 * t110;
154 t135 = t134 * t67;
155 t138 = sigma[0] * sigma[0];
156 t139 = t22 * t138;
157 t142 = 0.1e1 / t10 / t108 / t47;
158 t145 = 0.1e1 / t66 / t65;
159 t146 = t8 * t142 * t145;
160 t149 = 0.28e2 / 0.9e1 * t24 * t131 + 0.20e2 / 0.3e1 * t62 * t135 - 0.32e2 / 0.9e1 * t139 * t146;
161 t154 = 0.176e3 / 0.81e2 * t16 * t17 * t111 + 0.32e2 / 0.27e2 * t16 * t17 * t115 * t71 + 0.4e1 / 0.9e1 * t16 * t17 * t122 * t123 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t149;
162 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
163 v2rho2[0] = -t7 * t101 * t39 / 0.12e2 - t7 * t11 * t76 / 0.2e1 - 0.3e1 / 0.16e2 * t45 * t46 * t154;
164
165 t159 = t10 * t1 * t3;
166 t165 = t55 * t71;
167 t174 = t13 * t15 * sigma[0];
168 t175 = t121 * t91;
169 t176 = t175 * t71;
170 t177 = t87 * t176;
171 t184 = t22 * t8;
172 t185 = t108 * t18;
173 t187 = 0.1e1 / t10 / t185;
174 t192 = -0.2e1 / 0.3e1 * t85 * t60 - 0.2e1 * t22 * t68 + 0.4e1 / 0.3e1 * t184 * t187 * t145 * sigma[0];
175 t198 = t46 * (-0.16e2 / 0.27e2 * t80 * t81 * t50 - 0.2e1 / 0.9e1 * t16 * t87 * t165 + 0.16e2 / 0.27e2 * t16 * t17 * t115 * t91 + 0.4e1 / 0.9e1 * t174 * t177 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t192);
176 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
177 v2rhosigma[0] = -t159 * t97 / 0.4e1 - 0.3e1 / 0.16e2 * t45 * t198;
178
179 t201 = t55 * t91;
180 t205 = t91 * t91;
181 t212 = t22 / t23 / sigma[0];
182 t215 = 0.1e1 / sigma[0];
183 t216 = t22 * t215;
184 t219 = t108 * rho[0];
185 t222 = t8 / t10 / t219;
186 t223 = t222 * t145;
187 t226 = -t212 * t31 / 0.4e1 + t216 * t88 / 0.4e1 - t22 * t223 / 0.2e1;
188 t232 = t46 * (-0.4e1 / 0.9e1 * t16 * t87 * t201 + 0.4e1 / 0.9e1 * t16 * t17 * t122 * t205 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t226);
189 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
190 v2sigma2[0] = -0.3e1 / 0.16e2 * t45 * t232;
191
192 #ifndef XC_DONT_COMPILE_KXC
193
194 if(order < 3) return;
195
196
197 t236 = 0.1e1 / t19 / rho[0];
198 t237 = t9 * t236;
199 t248 = 0.1e1 / t19 / t219;
200 t249 = t248 * t34;
201 t253 = t110 * t55;
202 t258 = t49 * t121;
203 t267 = t54 * t54;
204 t268 = 0.1e1 / t267;
205 t269 = t21 * t268;
206 t270 = t123 * t71;
207 t275 = t121 * t71;
208 t276 = t275 * t149;
209 t277 = t87 * t276;
210 t283 = t8 / t10 / t108 * t30;
211 t287 = t9 * t248 * t67;
212 t290 = t108 * t108;
213 t292 = 0.1e1 / t10 / t290;
214 t297 = t138 * sigma[0];
215 t298 = t290 * t47;
216 t299 = 0.1e1 / t298;
217 t301 = t65 * t65;
218 t303 = 0.1e1 / t66 / t301;
219 t307 = -0.280e3 / 0.27e2 * t24 * t283 - 0.952e3 / 0.27e2 * t62 * t287 + 0.1184e4 / 0.27e2 * t139 * t8 * t292 * t145 - 0.256e3 / 0.9e1 * t22 * t297 * t299 * t303;
220 t312 = -0.2464e4 / 0.243e3 * t16 * t17 * t249 - 0.176e3 / 0.27e2 * t16 * t17 * t253 * t71 - 0.32e2 / 0.9e1 * t16 * t17 * t258 * t123 + 0.16e2 / 0.9e1 * t16 * t17 * t115 * t149 - 0.4e1 / 0.3e1 * t16 * t17 * t269 * t270 + 0.4e1 / 0.3e1 * t174 * t277 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t307;
221 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
222 v3rho3[0] = t7 * t237 * t39 / 0.18e2 - t7 * t101 * t76 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t154 - 0.3e1 / 0.16e2 * t45 * t46 * t312;
223
224 t317 = t100 * t1 * t3;
225 t328 = t121 * t123;
226 t332 = t55 * t149;
227 t340 = t63 * t176;
228 t348 = t268 * t91 * t123;
229 t349 = t87 * t348;
230 t352 = t121 * t192;
231 t353 = t352 * t71;
232 t354 = t87 * t353;
233 t357 = t175 * t149;
234 t358 = t87 * t357;
235 t369 = t290 * t18;
236 t371 = 0.1e1 / t369 * t303;
237 t375 = 0.14e2 / 0.9e1 * t85 * t131 + 0.74e2 / 0.9e1 * t22 * t135 - 0.124e3 / 0.9e1 * t184 * t142 * t145 * sigma[0] + 0.32e2 / 0.3e1 * t22 * t371 * t138;
238 t380 = 0.176e3 / 0.81e2 * t80 * t81 * t111 + 0.32e2 / 0.27e2 * t16 * t63 * t165 + 0.4e1 / 0.9e1 * t16 * t87 * t328 - 0.2e1 / 0.9e1 * t16 * t87 * t332 - 0.176e3 / 0.81e2 * t16 * t17 * t253 * t91 - 0.64e2 / 0.27e2 * t174 * t340 + 0.32e2 / 0.27e2 * t16 * t17 * t115 * t192 - 0.4e1 / 0.3e1 * t174 * t349 + 0.8e1 / 0.9e1 * t174 * t354 + 0.4e1 / 0.9e1 * t174 * t358 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t375;
239 t381 = t46 * t380;
240 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
241 v3rho2sigma[0] = -t317 * t97 / 0.12e2 - t159 * t198 / 0.2e1 - 0.3e1 / 0.16e2 * t45 * t381;
242
243 t391 = t55 * t192;
244 t399 = t268 * t205;
245 t400 = t399 * t71;
246 t401 = t87 * t400;
247 t404 = t175 * t192;
248 t405 = t87 * t404;
249 t412 = t121 * t226;
250 t413 = t412 * t71;
251 t414 = t87 * t413;
252 t422 = t8 * t187 * t145;
253 t425 = t290 * rho[0];
254 t427 = 0.1e1 / t425 * t303;
255 t431 = t212 * t60 / 0.3e1 - t216 * t68 / 0.3e1 + 0.10e2 / 0.3e1 * t22 * t422 - 0.4e1 * t22 * t427 * sigma[0];
256 t437 = t46 * (0.32e2 / 0.27e2 * t16 * t63 * t201 + 0.8e1 / 0.9e1 * t16 * t177 - 0.4e1 / 0.9e1 * t16 * t87 * t391 - 0.32e2 / 0.27e2 * t16 * t17 * t258 * t205 - 0.4e1 / 0.3e1 * t174 * t401 + 0.8e1 / 0.9e1 * t174 * t405 + 0.16e2 / 0.27e2 * t16 * t17 * t115 * t226 + 0.4e1 / 0.9e1 * t174 * t414 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t431);
257 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
258 v3rhosigma2[0] = -t159 * t232 / 0.4e1 - 0.3e1 / 0.16e2 * t45 * t437;
259
260 t440 = t121 * t205;
261 t444 = t55 * t226;
262 t448 = t205 * t91;
263 t453 = t175 * t226;
264 t454 = t87 * t453;
265 t459 = t22 / t23 / t138;
266 t463 = t22 / t138;
267 t468 = 0.1e1 / t290;
268 t472 = 0.3e1 / 0.8e1 * t459 * t31 - 0.3e1 / 0.8e1 * t463 * t88 - t216 * t223 / 0.4e1 + 0.3e1 / 0.2e1 * t22 * t468 * t303;
269 t478 = t46 * (0.4e1 / 0.3e1 * t16 * t87 * t440 - 0.2e1 / 0.3e1 * t16 * t87 * t444 - 0.4e1 / 0.3e1 * t16 * t17 * t269 * t448 + 0.4e1 / 0.3e1 * t174 * t454 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * t472);
270 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
271 v3sigma3[0] = -0.3e1 / 0.16e2 * t45 * t478;
272
273 #ifndef XC_DONT_COMPILE_LXC
274
275 if(order < 4) return;
276
277
278 t493 = 0.1e1 / t19 / t185;
279 t498 = t248 * t55;
280 t503 = t110 * t121;
281 t512 = t49 * t268;
282 t525 = 0.1e1 / t267 / t33;
283 t526 = t21 * t525;
284 t527 = t123 * t123;
285 t537 = t149 * t149;
286 t559 = t290 * t108;
287 t565 = t138 * t138;
288 t572 = 0.1e1 / t66 / t301 / t65;
289 t582 = 0.41888e5 / 0.729e3 * t16 * t17 * t493 * t34 + 0.9856e4 / 0.243e3 * t16 * t17 * t498 * t71 + 0.704e3 / 0.27e2 * t16 * t17 * t503 * t123 - 0.352e3 / 0.27e2 * t16 * t17 * t253 * t149 + 0.128e3 / 0.9e1 * t16 * t17 * t512 * t270 - 0.128e3 / 0.9e1 * t174 * t63 * t276 + 0.64e2 / 0.27e2 * t16 * t17 * t115 * t307 + 0.16e2 / 0.3e1 * t16 * t17 * t526 * t527 - 0.8e1 * t174 * t87 * t268 * t123 * t149 + 0.4e1 / 0.3e1 * t16 * t17 * t122 * t537 + 0.16e2 / 0.9e1 * t174 * t87 * t275 * t307 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * (0.3640e4 / 0.81e2 * t24 * t222 * t30 + 0.5768e4 / 0.27e2 * t62 * t9 * t493 * t67 - 0.37216e5 / 0.81e2 * t139 * t8 / t10 / t425 * t145 + 0.17920e5 / 0.27e2 * t22 * t297 / t559 * t303 - 0.5120e4 / 0.27e2 * t22 * t565 / t19 / t290 / t185 * t572 * t9);
290 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
291 v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t87 * t39 + 0.2e1 / 0.9e1 * t7 * t237 * t76 - t7 * t101 * t154 / 0.2e1 - t7 * t11 * t312 - 0.3e1 / 0.16e2 * t45 * t46 * t582;
292
293 t602 = t87 * t268;
294 t603 = t91 * t71;
295 t695 = 0.352e3 / 0.27e2 * t174 * t134 * t176 - 0.64e2 / 0.9e1 * t174 * t63 * t353 - 0.32e2 / 0.9e1 * t174 * t63 * t357 - 0.4e1 * t174 * t87 * t268 * t192 * t123 + 0.4e1 / 0.3e1 * t174 * t87 * t121 * t375 * t71 + 0.4e1 / 0.3e1 * t174 * t87 * t352 * t149 + 0.4e1 / 0.9e1 * t174 * t87 * t175 * t307 - 0.2464e4 / 0.243e3 * t80 * t81 * t249 - 0.176e3 / 0.27e2 * t16 * t134 * t165 + 0.16e2 / 0.9e1 * t16 * t63 * t332 - 0.2e1 / 0.9e1 * t16 * t87 * t55 * t307;
296 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
297 v4rho3sigma[0] = t236 * t1 * t3 * t97 / 0.18e2 - t317 * t198 / 0.4e1 - 0.3e1 / 0.4e1 * t159 * t381 - 0.3e1 / 0.16e2 * t45 * t46 * (0.32e2 / 0.3e1 * t174 * t63 * t348 + 0.16e2 / 0.3e1 * t174 * t87 * t525 * t91 * t270 - 0.4e1 * t174 * t602 * t603 * t149 - 0.32e2 / 0.9e1 * t16 * t63 * t328 - 0.4e1 / 0.3e1 * t16 * t87 * t268 * t270 + 0.4e1 / 0.3e1 * t16 * t277 - 0.176e3 / 0.27e2 * t16 * t17 * t253 * t192 + 0.16e2 / 0.9e1 * t16 * t17 * t115 * t375 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * (-0.140e3 / 0.27e2 * t85 * t283 - 0.364e3 / 0.9e1 * t22 * t287 + 0.3320e4 / 0.27e2 * t184 * t292 * t145 * sigma[0] - 0.1952e4 / 0.9e1 * t22 * t299 * t303 * t138 + 0.640e3 / 0.9e1 * t22 / t19 / t290 / t219 * t572 * t297 * t9) + 0.2464e4 / 0.243e3 * t16 * t17 * t498 * t91 + t695);
298
299 t758 = -0.4e1 / 0.9e1 * t16 * t87 * t55 * t375 + 0.16e2 / 0.3e1 * t174 * t87 * t525 * t205 * t123 - 0.16e2 / 0.3e1 * t174 * t602 * t603 * t192 - 0.4e1 / 0.3e1 * t174 * t87 * t268 * t226 * t123 - 0.352e3 / 0.81e2 * t16 * t134 * t201 - 0.128e3 / 0.27e2 * t16 * t340 + 0.16e2 / 0.9e1 * t16 * t354 + 0.8e1 / 0.9e1 * t16 * t358 + 0.32e2 / 0.27e2 * t16 * t17 * t115 * t431 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * (-0.7e1 / 0.9e1 * t212 * t131 + 0.7e1 / 0.9e1 * t216 * t135 - 0.22e2 * t22 * t146 + 0.188e3 / 0.3e1 * t22 * t371 * sigma[0] - 0.80e2 / 0.3e1 * t22 / t19 / t559 * t572 * t138 * t9) - 0.8e1 / 0.3e1 * t16 * t349;
300 t773 = t192 * t192;
301 t801 = 0.352e3 / 0.81e2 * t16 * t17 * t503 * t205 + 0.64e2 / 0.9e1 * t174 * t63 * t400 - 0.128e3 / 0.27e2 * t174 * t63 * t404 - 0.4e1 / 0.3e1 * t174 * t87 * t399 * t149 + 0.8e1 / 0.9e1 * t16 * t17 * t122 * t773 + 0.8e1 / 0.9e1 * t174 * t87 * t175 * t375 - 0.176e3 / 0.81e2 * t16 * t17 * t253 * t226 - 0.64e2 / 0.27e2 * t174 * t63 * t413 + 0.8e1 / 0.9e1 * t174 * t87 * t121 * t431 * t71 + 0.4e1 / 0.9e1 * t174 * t87 * t412 * t149 + 0.64e2 / 0.27e2 * t16 * t63 * t391;
302 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
303 v4rho2sigma2[0] = -t317 * t232 / 0.12e2 - t159 * t437 / 0.2e1 - 0.3e1 / 0.16e2 * t45 * t46 * (t758 + t801);
304
305 t882 = -0.32e2 / 0.9e1 * t16 * t63 * t440 - 0.4e1 * t16 * t401 + 0.8e1 / 0.3e1 * t16 * t405 + 0.16e2 / 0.9e1 * t16 * t63 * t444 + 0.4e1 / 0.3e1 * t16 * t414 - 0.2e1 / 0.3e1 * t16 * t87 * t55 * t431 + 0.32e2 / 0.9e1 * t16 * t17 * t512 * t448 + 0.16e2 / 0.3e1 * t174 * t87 * t525 * t448 * t71 - 0.4e1 * t174 * t87 * t399 * t192 - 0.32e2 / 0.9e1 * t174 * t63 * t453 - 0.4e1 * t174 * t602 * t91 * t226 * t71 + 0.4e1 / 0.3e1 * t174 * t87 * t352 * t226 + 0.4e1 / 0.3e1 * t174 * t87 * t175 * t431 + 0.16e2 / 0.27e2 * t16 * t17 * t115 * t472 + 0.4e1 / 0.9e1 * t174 * t87 * t121 * t472 * t71 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * (-t459 * t60 / 0.2e1 + t463 * t68 / 0.2e1 + t216 * t422 / 0.3e1 - 0.14e2 * t22 * t427 + 0.10e2 * t22 / t19 / t298 * t572 * sigma[0] * t9);
306 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
307 v4rhosigma3[0] = -t159 * t478 / 0.4e1 - 0.3e1 / 0.16e2 * t45 * t46 * t882;
308
309 t896 = t205 * t205;
310 t905 = t226 * t226;
311 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
312 v4sigma4[0] = -0.3e1 / 0.16e2 * t45 * t46 * (-0.16e2 / 0.3e1 * t16 * t87 * t268 * t448 + 0.16e2 / 0.3e1 * t16 * t454 - 0.8e1 / 0.9e1 * t16 * t87 * t55 * t472 + 0.16e2 / 0.3e1 * t16 * t17 * t526 * t896 - 0.8e1 * t174 * t87 * t399 * t226 + 0.4e1 / 0.3e1 * t16 * t17 * t122 * t905 + 0.16e2 / 0.9e1 * t174 * t87 * t175 * t472 - 0.2e1 / 0.9e1 * t16 * t17 * t56 * (-0.15e2 / 0.16e2 * t22 / t23 / t297 * t31 + 0.15e2 / 0.16e2 * t22 / t297 * t88 + 0.5e1 / 0.8e1 * t463 * t223 + 0.3e1 / 0.4e1 * t22 * t215 * t468 * t303 - 0.15e2 / 0.4e1 * t22 / t19 / t369 * t572 * t9));
313
314 #ifndef XC_DONT_COMPILE_MXC
315
316 if(order < 5) return;
317
318
319 #endif
320
321 #endif
322
323 #endif
324
325 #endif
326
327 #endif
328
329
330 }
331
332
333 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 *))334 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
335 {
336
337 #ifndef XC_DONT_COMPILE_EXC
338 double t1, t3, t4, t5, t6, t7, t8, t9;
339 double t10, t11, t12, t13, t14, t15, t17, t18;
340 double t19, t20, t21, t22, t23, t26, t27, t32;
341 double t34;
342
343 #ifndef XC_DONT_COMPILE_VXC
344 double t37, t38, t39, t41, t46, t47, t48, t49;
345 double t50, t52, t56, t57, t58, t59, t63, t64;
346 double t68, t72, t75, t82, t83, t88;
347
348 #ifndef XC_DONT_COMPILE_FXC
349 double t91, t92, t98, t100, t109, t110, t111, t116;
350 double t121, t125, t128, t131, t135, t136, t140, t144;
351 double t147, t157, t158, t169, t171, t172, t176, t177;
352 double t182, t188, t189, t194, t198, t202, t204, t208;
353 double t209, t214;
354
355 #ifndef XC_DONT_COMPILE_KXC
356 double t218, t219, t230, t244, t245, t246, t247, t251;
357 double t252, t257, t262, t266, t268, t273, t274, t275;
358 double t277, t279, t283, t284, t288, t292, t296, t319;
359 double t323, t324, t328, t339, t343, t345, t349, t350;
360 double t354, t355, t363, t364, t374, t375, t379, t386;
361 double t387, t401, t403, t406, t407, t412, t421, t422;
362 double t426, t431, t436, t445, t449, t450, t455;
363
364 #ifndef XC_DONT_COMPILE_LXC
365 double t471, t495, t496, t506, t529, t535, t542, t551;
366 double t602, t603, t664, t689, t713, t776, t819, t865;
367 double t878, t887;
368 #endif
369
370 #endif
371
372 #endif
373
374 #endif
375
376 #endif
377
378
379 gga_x_b88_params *params;
380
381 assert(p->params != NULL);
382 params = (gga_x_b88_params * )(p->params);
383
384 t1 = M_CBRT3;
385 t3 = POW_1_3(0.1e1 / M_PI);
386 t4 = t1 * t3;
387 t5 = M_CBRT4;
388 t6 = t5 * t5;
389 t7 = POW_1_3(rho[0]);
390 t8 = t6 * t7;
391 t9 = t1 * t1;
392 t10 = params->beta * t9;
393 t11 = 0.1e1 / t3;
394 t12 = t10 * t11;
395 t13 = t5 * sigma[0];
396 t14 = rho[0] * rho[0];
397 t15 = t7 * t7;
398 t17 = 0.1e1 / t15 / t14;
399 t18 = params->gamma * params->beta;
400 t19 = sqrt(sigma[0]);
401 t20 = t7 * rho[0];
402 t21 = 0.1e1 / t20;
403 t22 = t19 * t21;
404 t23 = log(t22 + sqrt(t22 * t22 + 0.1e1));
405 t26 = 0.10e1 + t18 * t22 * t23;
406 t27 = 0.1e1 / t26;
407 t32 = 0.10e1 + 0.2e1 / 0.9e1 * t12 * t13 * t17 * t27;
408 t34 = t4 * t8 * t32;
409 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
410 zk[0] = -0.3e1 / 0.8e1 * t34;
411
412 #ifndef XC_DONT_COMPILE_VXC
413
414 if(order < 1) return;
415
416
417 t37 = t20 * t1;
418 t38 = t3 * t6;
419 t39 = t14 * rho[0];
420 t41 = 0.1e1 / t15 / t39;
421 t46 = t11 * t5;
422 t47 = t10 * t46;
423 t48 = sigma[0] * t17;
424 t49 = t26 * t26;
425 t50 = 0.1e1 / t49;
426 t52 = 0.1e1 / t7 / t14;
427 t56 = sigma[0] * t41;
428 t57 = t48 + 0.1e1;
429 t58 = sqrt(t57);
430 t59 = 0.1e1 / t58;
431 t63 = -0.4e1 / 0.3e1 * t18 * t19 * t52 * t23 - 0.4e1 / 0.3e1 * t18 * t56 * t59;
432 t64 = t50 * t63;
433 t68 = -0.16e2 / 0.27e2 * t12 * t13 * t41 * t27 - 0.2e1 / 0.9e1 * t47 * t48 * t64;
434 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
435 vrho[0] = -t34 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * t68;
436
437 t72 = t5 * t17;
438 t75 = 0.1e1 / t19;
439 t82 = t18 * t75 * t21 * t23 / 0.2e1 + t18 * t17 * t59 / 0.2e1;
440 t83 = t50 * t82;
441 t88 = t38 * (0.2e1 / 0.9e1 * t12 * t72 * t27 - 0.2e1 / 0.9e1 * t47 * t48 * t83);
442 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
443 vsigma[0] = -0.3e1 / 0.8e1 * t37 * t88;
444
445 #ifndef XC_DONT_COMPILE_FXC
446
447 if(order < 2) return;
448
449
450 t91 = 0.1e1 / t15;
451 t92 = t6 * t91;
452 t98 = t14 * t14;
453 t100 = 0.1e1 / t15 / t98;
454 t109 = 0.1e1 / t49 / t26;
455 t110 = t63 * t63;
456 t111 = t109 * t110;
457 t116 = 0.1e1 / t7 / t39;
458 t121 = sigma[0] * t100;
459 t125 = sigma[0] * sigma[0];
460 t128 = 0.1e1 / t7 / t98 / t39;
461 t131 = 0.1e1 / t58 / t57;
462 t135 = 0.28e2 / 0.9e1 * t18 * t19 * t116 * t23 + 0.20e2 / 0.3e1 * t18 * t121 * t59 - 0.16e2 / 0.9e1 * t18 * t125 * t128 * t131;
463 t136 = t50 * t135;
464 t140 = 0.176e3 / 0.81e2 * t12 * t13 * t100 * t27 + 0.32e2 / 0.27e2 * t47 * t56 * t64 + 0.4e1 / 0.9e1 * t47 * t48 * t111 - 0.2e1 / 0.9e1 * t47 * t48 * t136;
465 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
466 v2rho2[0] = -t4 * t92 * t32 / 0.6e1 - t4 * t8 * t68 - 0.3e1 / 0.8e1 * t37 * t38 * t140;
467
468 t144 = t7 * t1;
469 t147 = t5 * t41;
470 t157 = t109 * t82;
471 t158 = t157 * t63;
472 t169 = t98 * t14;
473 t171 = 0.1e1 / t7 / t169;
474 t172 = t171 * t131;
475 t176 = -0.2e1 / 0.3e1 * t18 * t75 * t52 * t23 - 0.2e1 * t18 * t41 * t59 + 0.2e1 / 0.3e1 * t18 * t172 * sigma[0];
476 t177 = t50 * t176;
477 t182 = t38 * (-0.16e2 / 0.27e2 * t12 * t147 * t27 - 0.2e1 / 0.9e1 * t12 * t72 * t64 + 0.16e2 / 0.27e2 * t47 * t56 * t83 + 0.4e1 / 0.9e1 * t47 * t48 * t158 - 0.2e1 / 0.9e1 * t47 * t48 * t177);
478 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
479 v2rhosigma[0] = -t144 * t88 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t182;
480
481 t188 = t82 * t82;
482 t189 = t109 * t188;
483 t194 = 0.1e1 / t19 / sigma[0];
484 t198 = 0.1e1 / sigma[0];
485 t202 = t98 * rho[0];
486 t204 = 0.1e1 / t7 / t202;
487 t208 = t18 * t198 * t17 * t59 / 0.4e1 - t18 * t194 * t21 * t23 / 0.4e1 - t18 * t204 * t131 / 0.4e1;
488 t209 = t50 * t208;
489 t214 = t38 * (-0.4e1 / 0.9e1 * t12 * t72 * t83 + 0.4e1 / 0.9e1 * t47 * t48 * t189 - 0.2e1 / 0.9e1 * t47 * t48 * t209);
490 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
491 v2sigma2[0] = -0.3e1 / 0.8e1 * t37 * t214;
492
493 #ifndef XC_DONT_COMPILE_KXC
494
495 if(order < 3) return;
496
497
498 t218 = 0.1e1 / t15 / rho[0];
499 t219 = t6 * t218;
500 t230 = 0.1e1 / t15 / t202;
501 t244 = t49 * t49;
502 t245 = 0.1e1 / t244;
503 t246 = t110 * t63;
504 t247 = t245 * t246;
505 t251 = t109 * t63;
506 t252 = t251 * t135;
507 t257 = 0.1e1 / t7 / t98;
508 t262 = sigma[0] * t230;
509 t266 = t98 * t98;
510 t268 = 0.1e1 / t7 / t266;
511 t273 = t125 * sigma[0];
512 t274 = t266 * t39;
513 t275 = 0.1e1 / t274;
514 t277 = t57 * t57;
515 t279 = 0.1e1 / t58 / t277;
516 t283 = -0.280e3 / 0.27e2 * t18 * t19 * t257 * t23 - 0.952e3 / 0.27e2 * t18 * t262 * t59 + 0.592e3 / 0.27e2 * t18 * t125 * t268 * t131 - 0.64e2 / 0.9e1 * t18 * t273 * t275 * t279;
517 t284 = t50 * t283;
518 t288 = -0.2464e4 / 0.243e3 * t12 * t13 * t230 * t27 - 0.176e3 / 0.27e2 * t47 * t121 * t64 - 0.32e2 / 0.9e1 * t47 * t56 * t111 + 0.16e2 / 0.9e1 * t47 * t56 * t136 - 0.4e1 / 0.3e1 * t47 * t48 * t247 + 0.4e1 / 0.3e1 * t47 * t48 * t252 - 0.2e1 / 0.9e1 * t47 * t48 * t284;
519 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
520 v3rho3[0] = t4 * t219 * t32 / 0.9e1 - t4 * t92 * t68 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t140 - 0.3e1 / 0.8e1 * t37 * t38 * t288;
521
522 t292 = t91 * t1;
523 t296 = t5 * t100;
524 t319 = t245 * t82 * t110;
525 t323 = t109 * t176;
526 t324 = t323 * t63;
527 t328 = t157 * t135;
528 t339 = t128 * t131;
529 t343 = t266 * t14;
530 t345 = 0.1e1 / t343 * t279;
531 t349 = 0.14e2 / 0.9e1 * t18 * t75 * t116 * t23 + 0.74e2 / 0.9e1 * t18 * t100 * t59 - 0.62e2 / 0.9e1 * t18 * t339 * sigma[0] + 0.8e1 / 0.3e1 * t18 * t345 * t125;
532 t350 = t50 * t349;
533 t354 = 0.176e3 / 0.81e2 * t12 * t296 * t27 + 0.32e2 / 0.27e2 * t12 * t147 * t64 + 0.4e1 / 0.9e1 * t12 * t72 * t111 - 0.2e1 / 0.9e1 * t12 * t72 * t136 - 0.176e3 / 0.81e2 * t47 * t121 * t83 - 0.64e2 / 0.27e2 * t47 * t56 * t158 + 0.32e2 / 0.27e2 * t47 * t56 * t177 - 0.4e1 / 0.3e1 * t47 * t48 * t319 + 0.8e1 / 0.9e1 * t47 * t48 * t324 + 0.4e1 / 0.9e1 * t47 * t48 * t328 - 0.2e1 / 0.9e1 * t47 * t48 * t350;
534 t355 = t38 * t354;
535 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
536 v3rho2sigma[0] = -t292 * t88 / 0.6e1 - t144 * t182 - 0.3e1 / 0.8e1 * t37 * t355;
537
538 t363 = t17 * t109;
539 t364 = t82 * t63;
540 t374 = t245 * t188;
541 t375 = t374 * t63;
542 t379 = t157 * t176;
543 t386 = t109 * t208;
544 t387 = t386 * t63;
545 t401 = t266 * rho[0];
546 t403 = 0.1e1 / t401 * t279;
547 t406 = t18 * t194 * t52 * t23 / 0.3e1 - t18 * t198 * t41 * t59 / 0.3e1 + 0.5e1 / 0.3e1 * t18 * t172 - t18 * t403 * sigma[0];
548 t407 = t50 * t406;
549 t412 = t38 * (0.32e2 / 0.27e2 * t12 * t147 * t83 + 0.8e1 / 0.9e1 * t47 * t363 * t364 - 0.4e1 / 0.9e1 * t12 * t72 * t177 - 0.32e2 / 0.27e2 * t47 * t56 * t189 - 0.4e1 / 0.3e1 * t47 * t48 * t375 + 0.8e1 / 0.9e1 * t47 * t48 * t379 + 0.16e2 / 0.27e2 * t47 * t56 * t209 + 0.4e1 / 0.9e1 * t47 * t48 * t387 - 0.2e1 / 0.9e1 * t47 * t48 * t407);
550 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
551 v3rhosigma2[0] = -t144 * t214 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t412;
552
553 t421 = t188 * t82;
554 t422 = t245 * t421;
555 t426 = t157 * t208;
556 t431 = 0.1e1 / t19 / t125;
557 t436 = 0.1e1 / t125;
558 t445 = 0.1e1 / t266;
559 t449 = 0.3e1 / 0.8e1 * t18 * t431 * t21 * t23 - 0.3e1 / 0.8e1 * t18 * t436 * t17 * t59 - t18 * t198 * t204 * t131 / 0.8e1 + 0.3e1 / 0.8e1 * t18 * t445 * t279;
560 t450 = t50 * t449;
561 t455 = t38 * (0.4e1 / 0.3e1 * t12 * t72 * t189 - 0.2e1 / 0.3e1 * t12 * t72 * t209 - 0.4e1 / 0.3e1 * t47 * t48 * t422 + 0.4e1 / 0.3e1 * t47 * t48 * t426 - 0.2e1 / 0.9e1 * t47 * t48 * t450);
562 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
563 v3sigma3[0] = -0.3e1 / 0.8e1 * t37 * t455;
564
565 #ifndef XC_DONT_COMPILE_LXC
566
567 if(order < 4) return;
568
569
570 t471 = 0.1e1 / t15 / t169;
571 t495 = 0.1e1 / t244 / t26;
572 t496 = t110 * t110;
573 t506 = t135 * t135;
574 t529 = t266 * t98;
575 t535 = t125 * t125;
576 t542 = 0.1e1 / t58 / t277 / t57;
577 t551 = 0.41888e5 / 0.729e3 * t12 * t13 * t471 * t27 + 0.9856e4 / 0.243e3 * t47 * t262 * t64 + 0.704e3 / 0.27e2 * t47 * t121 * t111 - 0.352e3 / 0.27e2 * t47 * t121 * t136 + 0.128e3 / 0.9e1 * t47 * t56 * t247 - 0.128e3 / 0.9e1 * t47 * t56 * t252 + 0.64e2 / 0.27e2 * t47 * t56 * t284 + 0.16e2 / 0.3e1 * t47 * t48 * t495 * t496 - 0.8e1 * t47 * t48 * t245 * t110 * t135 + 0.4e1 / 0.3e1 * t47 * t48 * t109 * t506 + 0.16e2 / 0.9e1 * t47 * t48 * t251 * t283 - 0.2e1 / 0.9e1 * t47 * t48 * t50 * (0.3640e4 / 0.81e2 * t18 * t19 * t204 * t23 + 0.5768e4 / 0.27e2 * t18 * sigma[0] * t471 * t59 - 0.18608e5 / 0.81e2 * t18 * t125 / t7 / t401 * t131 + 0.4480e4 / 0.27e2 * t18 * t273 / t529 * t279 - 0.1280e4 / 0.27e2 * t18 * t535 / t15 / t266 / t169 * t542);
578 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
579 v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t17 * t32 + 0.4e1 / 0.9e1 * t4 * t219 * t68 - t4 * t92 * t140 - 0.2e1 * t4 * t8 * t288 - 0.3e1 / 0.8e1 * t37 * t38 * t551;
580
581 t602 = t10 * t46 * sigma[0];
582 t603 = t17 * t245;
583 t664 = -0.4e1 * t602 * t603 * t364 * t135 - 0.32e2 / 0.9e1 * t12 * t147 * t111 - 0.4e1 / 0.3e1 * t12 * t72 * t247 + 0.4e1 / 0.3e1 * t47 * t363 * t63 * t135 - 0.176e3 / 0.27e2 * t47 * t121 * t177 + 0.16e2 / 0.9e1 * t47 * t56 * t350 - 0.2e1 / 0.9e1 * t47 * t48 * t50 * (-0.140e3 / 0.27e2 * t18 * t75 * t257 * t23 - 0.364e3 / 0.9e1 * t18 * t230 * t59 + 0.1660e4 / 0.27e2 * t18 * t268 * t131 * sigma[0] - 0.488e3 / 0.9e1 * t18 * t275 * t279 * t125 + 0.160e3 / 0.9e1 * t18 / t15 / t266 / t202 * t542 * t273) - 0.2464e4 / 0.243e3 * t12 * t5 * t230 * t27 - 0.176e3 / 0.27e2 * t12 * t296 * t64 + 0.16e2 / 0.9e1 * t12 * t147 * t136 - 0.2e1 / 0.9e1 * t12 * t72 * t284;
584 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
585 v4rho3sigma[0] = t218 * t1 * t88 / 0.9e1 - t292 * t182 / 0.2e1 - 0.3e1 / 0.2e1 * t144 * t355 - 0.3e1 / 0.8e1 * t37 * t38 * (0.2464e4 / 0.243e3 * t47 * t262 * t83 + 0.352e3 / 0.27e2 * t47 * t121 * t158 - 0.64e2 / 0.9e1 * t47 * t56 * t324 - 0.32e2 / 0.9e1 * t47 * t56 * t328 - 0.4e1 * t47 * t48 * t245 * t176 * t110 + 0.4e1 / 0.3e1 * t47 * t48 * t109 * t349 * t63 + 0.4e1 / 0.3e1 * t47 * t48 * t323 * t135 + 0.4e1 / 0.9e1 * t47 * t48 * t157 * t283 + 0.32e2 / 0.3e1 * t47 * t56 * t319 + 0.16e2 / 0.3e1 * t47 * t48 * t495 * t82 * t246 + t664);
586
587 t689 = t176 * t176;
588 t713 = -0.8e1 / 0.3e1 * t47 * t603 * t82 * t110 + 0.352e3 / 0.81e2 * t47 * t121 * t189 + 0.64e2 / 0.9e1 * t47 * t56 * t375 - 0.128e3 / 0.27e2 * t47 * t56 * t379 - 0.4e1 / 0.3e1 * t47 * t48 * t374 * t135 + 0.8e1 / 0.9e1 * t47 * t48 * t109 * t689 + 0.8e1 / 0.9e1 * t47 * t48 * t157 * t349 - 0.176e3 / 0.81e2 * t47 * t121 * t209 - 0.64e2 / 0.27e2 * t47 * t56 * t387 + 0.8e1 / 0.9e1 * t47 * t48 * t109 * t406 * t63 + 0.4e1 / 0.9e1 * t47 * t48 * t386 * t135;
589 t776 = 0.16e2 / 0.3e1 * t47 * t48 * t495 * t188 * t110 - 0.16e2 / 0.3e1 * t602 * t603 * t364 * t176 - 0.4e1 / 0.3e1 * t47 * t48 * t245 * t208 * t110 - 0.352e3 / 0.81e2 * t12 * t296 * t83 - 0.128e3 / 0.27e2 * t47 * t41 * t109 * t364 + 0.16e2 / 0.9e1 * t47 * t363 * t176 * t63 + 0.8e1 / 0.9e1 * t47 * t363 * t82 * t135 + 0.32e2 / 0.27e2 * t47 * t56 * t407 - 0.2e1 / 0.9e1 * t47 * t48 * t50 * (-0.7e1 / 0.9e1 * t18 * t194 * t116 * t23 + 0.7e1 / 0.9e1 * t18 * t198 * t100 * t59 - 0.11e2 * t18 * t339 + 0.47e2 / 0.3e1 * t18 * t345 * sigma[0] - 0.20e2 / 0.3e1 * t18 / t15 / t529 * t542 * t125) + 0.64e2 / 0.27e2 * t12 * t147 * t177 - 0.4e1 / 0.9e1 * t12 * t72 * t350;
590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
591 v4rho2sigma2[0] = -t292 * t214 / 0.6e1 - t144 * t412 - 0.3e1 / 0.8e1 * t37 * t38 * (t713 + t776);
592
593 t819 = t82 * t208;
594 t865 = -0.32e2 / 0.9e1 * t12 * t147 * t189 - 0.4e1 * t47 * t603 * t188 * t63 + 0.8e1 / 0.3e1 * t47 * t363 * t82 * t176 + 0.16e2 / 0.9e1 * t12 * t147 * t209 + 0.4e1 / 0.3e1 * t47 * t363 * t208 * t63 - 0.2e1 / 0.3e1 * t12 * t72 * t407 + 0.32e2 / 0.9e1 * t47 * t56 * t422 + 0.16e2 / 0.3e1 * t47 * t48 * t495 * t421 * t63 - 0.4e1 * t47 * t48 * t374 * t176 - 0.32e2 / 0.9e1 * t47 * t56 * t426 - 0.4e1 * t602 * t603 * t819 * t63 + 0.4e1 / 0.3e1 * t47 * t48 * t323 * t208 + 0.4e1 / 0.3e1 * t47 * t48 * t157 * t406 + 0.16e2 / 0.27e2 * t47 * t56 * t450 + 0.4e1 / 0.9e1 * t47 * t48 * t109 * t449 * t63 - 0.2e1 / 0.9e1 * t47 * t48 * t50 * (-t18 * t431 * t52 * t23 / 0.2e1 + t18 * t436 * t41 * t59 / 0.2e1 + t18 * t198 * t171 * t131 / 0.6e1 - 0.7e1 / 0.2e1 * t18 * t403 + 0.5e1 / 0.2e1 * t18 / t15 / t274 * t542 * sigma[0]);
595 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
596 v4rhosigma3[0] = -t144 * t455 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * t865;
597
598 t878 = t188 * t188;
599 t887 = t208 * t208;
600 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
601 v4sigma4[0] = -0.3e1 / 0.8e1 * t37 * t38 * (-0.16e2 / 0.3e1 * t12 * t72 * t422 + 0.16e2 / 0.3e1 * t47 * t363 * t819 - 0.8e1 / 0.9e1 * t12 * t72 * t450 + 0.16e2 / 0.3e1 * t47 * t48 * t495 * t878 - 0.8e1 * t47 * t48 * t374 * t208 + 0.4e1 / 0.3e1 * t47 * t48 * t109 * t887 + 0.16e2 / 0.9e1 * t47 * t48 * t157 * t449 - 0.2e1 / 0.9e1 * t47 * t48 * t50 * (-0.15e2 / 0.16e2 * t18 / t19 / t273 * t21 * t23 + 0.15e2 / 0.16e2 * t18 / t273 * t17 * t59 + 0.5e1 / 0.16e2 * t18 * t436 * t204 * t131 + 0.3e1 / 0.16e2 * t18 * t198 * t445 * t279 - 0.15e2 / 0.16e2 * t18 / t15 / t343 * t542));
602
603 #ifndef XC_DONT_COMPILE_MXC
604
605 if(order < 5) return;
606
607
608 #endif
609
610 #endif
611
612 #endif
613
614 #endif
615
616 #endif
617
618
619 }
620
621
622 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 *))623 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
624 {
625
626 #ifndef XC_DONT_COMPILE_EXC
627 double t1, t3, t4, t5, t6, t7, t8, t9;
628 double t10, t11, t13, t14, t15, t16, t17, t18;
629 double t19, t20, t21, t22, t23, t24, t25, t27;
630 double t28, t29, t31, t32, t33, t36, t37, t42;
631 double t44, t46, t47, t48, t49, t50, t51, t52;
632 double t53, t55, t56, t58, t59, t60, t63, t64;
633 double t69, t71;
634
635 #ifndef XC_DONT_COMPILE_VXC
636 double t73, t74, t75, t76, t77, t78, t80, t81;
637 double t83, t84, t85, t86, t87, t89, t90, t91;
638 double t93, t98, t99, t100, t101, t102, t104, t108;
639 double t109, t110, t111, t115, t116, t120, t122, t123;
640 double t124, t125, t126, t128, t129, t130, t132, t133;
641 double t137, t138, t140, t141, t142, t143, t145, t146;
642 double t147, t149, t154, t155, t156, t158, t162, t163;
643 double t164, t165, t169, t170, t174, t176, t177, t181;
644 double t182, t183, t184, t187, t194, t195, t199, t203;
645 double t204, t207, t214, t215, t219;
646
647 #ifndef XC_DONT_COMPILE_FXC
648 double t223, t225, t226, t227, t228, t229, t230, t232;
649 double t233, t234, t236, t237, t238, t240, t241, t242;
650 double t243, t244, t245, t247, t248, t250, t251, t253;
651 double t254, t256, t257, t258, t260, t269, t270, t271;
652 double t276, t281, t285, t288, t291, t295, t296, t300;
653 double t302, t303, t304, t305, t306, t307, t308, t310;
654 double t311, t312, t314, t315, t316, t317, t319, t320;
655 double t321, t323, t324, t325, t327, t328, t329, t330;
656 double t332, t335, t337, t339, t342, t343, t344, t346;
657 double t350, t351, t352, t353, t355, t358, t360, t361;
658 double t363, t366, t368, t371, t373, t376, t377, t379;
659 double t380, t381, t382, t383, t385, t386, t387, t388;
660 double t390, t391, t392, t393, t395, t396, t397, t399;
661 double t400, t401, t402, t404, t413, t414, t415, t420;
662 double t425, t429, t432, t435, t439, t440, t444, t446;
663 double t447, t448, t451, t453, t454, t458, t468, t469;
664 double t480, t482, t483, t487, t488, t492, t497, t498;
665 double t499, t504, t509, t512, t522, t523, t534, t536;
666 double t537, t541, t542, t546, t547, t553, t554, t559;
667 double t563, t567, t569, t573, t574, t578, t579, t585;
668 double t586, t591, t595, t599, t601, t605, t606, t610;
669 double t611;
670
671 #ifndef XC_DONT_COMPILE_KXC
672 double t618, t620, t623, t625, t626, t628, t630, t632;
673 double t634, t635, t636, t638, t639, t640, t641, t642;
674 double t644, t645, t647, t648, t650, t651, t653, t656;
675 double t670, t671, t672, t673, t677, t678, t683, t688;
676 double t692, t694, t699, t700, t701, t703, t705, t709;
677 double t710, t714, t716, t717, t719, t720, t721, t722;
678 double t724, t725, t727, t729, t731, t733, t734, t735;
679 double t737, t738, t740, t741, t742, t743, t745, t746;
680 double t747, t749, t751, t753, t754, t755, t757, t758;
681 double t759, t761, t762, t764, t766, t767, t769, t770;
682 double t771, t773, t774, t775, t780, t781, t784, t786;
683 double t788, t790, t792, t793, t795, t797, t798, t799;
684 double t801, t803, t805, t806, t807, t808, t810, t812;
685 double t813, t814, t816, t818, t820, t821, t822, t824;
686 double t827, t829, t831, t834, t835, t837, t838, t839;
687 double t841, t844, t845, t847, t849, t851, t853, t855;
688 double t857, t859, t861, t863, t872, t875, t878, t882;
689 double t885, t890, t892, t894, t896, t898, t900, t902;
690 double t904, t906, t908, t910, t912, t914, t916, t919;
691 double t920, t922, t924, t926, t929, t932, t935, t938;
692 double t941, t943, t944, t946, t948, t950, t952, t953;
693 double t955, t957, t959, t967, t971, t974, t982, t983;
694 double t984, t987, t988, t990, t991, t992, t994, t995;
695 double t996, t998, t999, t1000, t1002, t1003, t1005, t1006;
696 double t1008, t1009, t1014, t1015, t1017, t1018, t1019, t1021;
697 double t1022, t1023, t1024, t1026, t1027, t1028, t1030, t1031;
698 double t1033, t1047, t1048, t1049, t1050, t1054, t1055, t1060;
699 double t1065, t1069, t1071, t1076, t1077, t1078, t1080, t1082;
700 double t1086, t1087, t1091, t1093, t1094, t1095, t1099, t1103;
701 double t1105, t1106, t1110, t1113, t1117, t1140, t1144, t1145;
702 double t1149, t1160, t1164, t1166, t1170, t1171, t1175, t1180;
703 double t1184, t1185, t1189, t1195, t1196, t1197, t1199, t1200;
704 double t1201, t1202, t1206, t1209, t1210, t1212, t1216, t1218;
705 double t1219, t1223, t1226, t1229, t1234, t1238, t1243, t1247;
706 double t1250, t1253, t1276, t1280, t1281, t1285, t1296, t1300;
707 double t1302, t1306, t1307, t1311, t1312, t1315, t1318, t1324;
708 double t1325, t1335, t1336, t1340, t1347, t1348, t1362, t1364;
709 double t1367, t1368, t1372, t1373, t1376, t1378, t1382, t1387;
710 double t1393, t1394, t1404, t1405, t1409, t1416, t1417, t1431;
711 double t1433, t1436, t1437, t1441, t1442, t1451, t1452, t1456;
712 double t1461, t1466, t1475, t1479, t1480, t1484, t1485, t1494;
713 double t1495, t1499, t1504, t1509, t1518, t1522, t1523, t1527;
714 double t1528;
715
716 #ifndef XC_DONT_COMPILE_LXC
717 double t1532, t1542, t1566, t1567, t1577, t1600, t1606, t1613;
718 double t1622, t1627, t1630, t1635, t1639, t1641, t1648, t1662;
719 double t1664, t1666, t1672, t1675, t1686, t1688, t1691, t1694;
720 double t1696, t1698, t1699, t1700, t1705, t1707, t1709, t1711;
721 double t1713, t1718, t1720, t1722, t1727, t1729, t1730, t1731;
722 double t1736, t1738, t1744, t1747, t1750, t1761, t1764, t1766;
723 double t1778, t1786, t1789, t1790, t1792, t1793, t1795, t1798;
724 double t1805, t1808, t1809, t1811, t1814, t1817, t1819, t1828;
725 double t1840, t1843, t1845, t1849, t1851, t1857, t1865, t1866;
726 double t1867, t1869, t1872, t1874, t1886, t1888, t1890, t1897;
727 double t1909, t1916, t1922, t1926, t1927, t1938, t1941, t1947;
728 double t1952, t1976, t1983, t1985, t1986, t1993, t1997, t2003;
729 double t2006, t2013, t2020, t2027, t2034, t2037, t2051, t2053;
730 double t2059, t2062, t2065, t2068, t2071, t2074, t2076, t2078;
731 double t2104, t2107, t2110, t2113, t2127, t2145, t2155, t2183;
732 double t2189, t2205, t2221, t2232, t2244, t2255, t2258, t2312;
733 double t2314, t2320, t2337, t2340, t2343, t2355, t2364, t2376;
734 double t2379, t2382, t2385, t2400, t2412, t2420, t2424, t2432;
735 double t2435, t2440, t2445, t2450, t2455, t2477, t2488, t2512;
736 double t2513, t2523, t2546, t2552, t2559, t2568, t2593, t2611;
737 double t2621, t2633, t2661, t2662, t2736, t2742, t2745, t2764;
738 double t2767, t2771, t2780, t2783, t2786, t2789, t2792, t2796;
739 double t2819, t2820, t2822, t2824, t2826, t2827, t2829, t2833;
740 double t2843, t2855, t2858, t2859, t2861, t2874, t2877, t2890;
741 double t2902, t2919, t2927, t2936, t2945, t3082, t3083, t3091;
742 double t3102, t3103, t3105, t3165, t3177, t3221, t3227, t3228;
743 double t3240, t3255, t3265, t3349, t3354, t3396, t3401, t3443;
744 double t3489, t3493, t3543, t3589, t3602, t3611, t3662, t3671;
745 #endif
746
747 #endif
748
749 #endif
750
751 #endif
752
753 #endif
754
755
756 gga_x_b88_params *params;
757
758 assert(p->params != NULL);
759 params = (gga_x_b88_params * )(p->params);
760
761 t1 = M_CBRT3;
762 t3 = POW_1_3(0.1e1 / M_PI);
763 t4 = t1 * t3;
764 t5 = M_CBRT4;
765 t6 = t5 * t5;
766 t7 = t4 * t6;
767 t8 = rho[0] - rho[1];
768 t9 = rho[0] + rho[1];
769 t10 = 0.1e1 / t9;
770 t11 = t8 * t10;
771 t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
772 t14 = POW_1_3(t13);
773 t15 = t14 * t13;
774 t16 = POW_1_3(t9);
775 t17 = t15 * t16;
776 t18 = t1 * t1;
777 t19 = params->beta * t18;
778 t20 = 0.1e1 / t3;
779 t21 = t19 * t20;
780 t22 = t5 * sigma[0];
781 t23 = rho[0] * rho[0];
782 t24 = POW_1_3(rho[0]);
783 t25 = t24 * t24;
784 t27 = 0.1e1 / t25 / t23;
785 t28 = params->gamma * params->beta;
786 t29 = sqrt(sigma[0]);
787 t31 = 0.1e1 / t24 / rho[0];
788 t32 = t29 * t31;
789 t33 = log(t32 + sqrt(t32 * t32 + 0.1e1));
790 t36 = 0.10e1 + t28 * t32 * t33;
791 t37 = 0.1e1 / t36;
792 t42 = 0.10e1 + 0.2e1 / 0.9e1 * t21 * t22 * t27 * t37;
793 t44 = t7 * t17 * t42;
794 t46 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
795 t47 = POW_1_3(t46);
796 t48 = t47 * t46;
797 t49 = t48 * t16;
798 t50 = t5 * sigma[2];
799 t51 = rho[1] * rho[1];
800 t52 = POW_1_3(rho[1]);
801 t53 = t52 * t52;
802 t55 = 0.1e1 / t53 / t51;
803 t56 = sqrt(sigma[2]);
804 t58 = 0.1e1 / t52 / rho[1];
805 t59 = t56 * t58;
806 t60 = log(t59 + sqrt(t59 * t59 + 0.1e1));
807 t63 = 0.10e1 + t28 * t59 * t60;
808 t64 = 0.1e1 / t63;
809 t69 = 0.10e1 + 0.2e1 / 0.9e1 * t21 * t50 * t55 * t64;
810 t71 = t7 * t49 * t69;
811 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
812 zk[0] = -0.3e1 / 0.8e1 * t44 - 0.3e1 / 0.8e1 * t71;
813
814 #ifndef XC_DONT_COMPILE_VXC
815
816 if(order < 1) return;
817
818
819 t73 = 0.3e1 / 0.8e1 * t44;
820 t74 = 0.3e1 / 0.8e1 * t71;
821 t75 = t14 * t16;
822 t76 = t9 * t9;
823 t77 = 0.1e1 / t76;
824 t78 = t8 * t77;
825 t80 = t10 / 0.2e1 - t78 / 0.2e1;
826 t81 = t42 * t80;
827 t83 = t7 * t75 * t81;
828 t84 = t83 / 0.2e1;
829 t85 = t16 * t16;
830 t86 = 0.1e1 / t85;
831 t87 = t15 * t86;
832 t89 = t7 * t87 * t42;
833 t90 = t89 / 0.8e1;
834 t91 = t23 * rho[0];
835 t93 = 0.1e1 / t25 / t91;
836 t98 = t20 * t5;
837 t99 = t19 * t98;
838 t100 = sigma[0] * t27;
839 t101 = t36 * t36;
840 t102 = 0.1e1 / t101;
841 t104 = 0.1e1 / t24 / t23;
842 t108 = sigma[0] * t93;
843 t109 = t100 + 0.1e1;
844 t110 = sqrt(t109);
845 t111 = 0.1e1 / t110;
846 t115 = -0.4e1 / 0.3e1 * t28 * t29 * t104 * t33 - 0.4e1 / 0.3e1 * t28 * t108 * t111;
847 t116 = t102 * t115;
848 t120 = -0.16e2 / 0.27e2 * t21 * t22 * t93 * t37 - 0.2e1 / 0.9e1 * t99 * t100 * t116;
849 t122 = t7 * t17 * t120;
850 t123 = 0.3e1 / 0.8e1 * t122;
851 t124 = t47 * t16;
852 t125 = -t80;
853 t126 = t69 * t125;
854 t128 = t7 * t124 * t126;
855 t129 = t128 / 0.2e1;
856 t130 = t48 * t86;
857 t132 = t7 * t130 * t69;
858 t133 = t132 / 0.8e1;
859 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
860 vrho[0] = -t73 - t74 + t9 * (-t84 - t90 - t123 - t129 - t133);
861
862 t137 = -t10 / 0.2e1 - t78 / 0.2e1;
863 t138 = t42 * t137;
864 t140 = t7 * t75 * t138;
865 t141 = t140 / 0.2e1;
866 t142 = -t137;
867 t143 = t69 * t142;
868 t145 = t7 * t124 * t143;
869 t146 = t145 / 0.2e1;
870 t147 = t51 * rho[1];
871 t149 = 0.1e1 / t53 / t147;
872 t154 = sigma[2] * t55;
873 t155 = t63 * t63;
874 t156 = 0.1e1 / t155;
875 t158 = 0.1e1 / t52 / t51;
876 t162 = sigma[2] * t149;
877 t163 = t154 + 0.1e1;
878 t164 = sqrt(t163);
879 t165 = 0.1e1 / t164;
880 t169 = -0.4e1 / 0.3e1 * t28 * t56 * t158 * t60 - 0.4e1 / 0.3e1 * t28 * t162 * t165;
881 t170 = t156 * t169;
882 t174 = -0.16e2 / 0.27e2 * t21 * t50 * t149 * t64 - 0.2e1 / 0.9e1 * t99 * t154 * t170;
883 t176 = t7 * t49 * t174;
884 t177 = 0.3e1 / 0.8e1 * t176;
885 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
886 vrho[1] = -t73 - t74 + t9 * (-t141 - t90 - t146 - t133 - t177);
887
888 t181 = t16 * t9 * t1;
889 t182 = t181 * t3;
890 t183 = t6 * t15;
891 t184 = t5 * t27;
892 t187 = 0.1e1 / t29;
893 t194 = t28 * t187 * t31 * t33 / 0.2e1 + t28 * t27 * t111 / 0.2e1;
894 t195 = t102 * t194;
895 t199 = -0.2e1 / 0.9e1 * t99 * t100 * t195 + 0.2e1 / 0.9e1 * t21 * t184 * t37;
896 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
897 vsigma[0] = -0.3e1 / 0.8e1 * t182 * t183 * t199;
898
899 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
900 vsigma[1] = 0.0e0;
901
902 t203 = t6 * t48;
903 t204 = t5 * t55;
904 t207 = 0.1e1 / t56;
905 t214 = t28 * t207 * t58 * t60 / 0.2e1 + t28 * t55 * t165 / 0.2e1;
906 t215 = t156 * t214;
907 t219 = -0.2e1 / 0.9e1 * t99 * t154 * t215 + 0.2e1 / 0.9e1 * t21 * t204 * t64;
908 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
909 vsigma[2] = -0.3e1 / 0.8e1 * t182 * t203 * t219;
910
911 #ifndef XC_DONT_COMPILE_FXC
912
913 if(order < 2) return;
914
915
916 t223 = t89 / 0.4e1;
917 t225 = t132 / 0.4e1;
918 t226 = t14 * t14;
919 t227 = 0.1e1 / t226;
920 t228 = t227 * t16;
921 t229 = t80 * t80;
922 t230 = t42 * t229;
923 t232 = t7 * t228 * t230;
924 t233 = t232 / 0.6e1;
925 t234 = t14 * t86;
926 t236 = t7 * t234 * t81;
927 t237 = t236 / 0.3e1;
928 t238 = t120 * t80;
929 t240 = t7 * t75 * t238;
930 t241 = t76 * t9;
931 t242 = 0.1e1 / t241;
932 t243 = t8 * t242;
933 t244 = -t77 + t243;
934 t245 = t42 * t244;
935 t247 = t7 * t75 * t245;
936 t248 = t247 / 0.2e1;
937 t250 = 0.1e1 / t85 / t9;
938 t251 = t15 * t250;
939 t253 = t7 * t251 * t42;
940 t254 = t253 / 0.12e2;
941 t256 = t7 * t87 * t120;
942 t257 = t256 / 0.4e1;
943 t258 = t23 * t23;
944 t260 = 0.1e1 / t25 / t258;
945 t269 = 0.1e1 / t101 / t36;
946 t270 = t115 * t115;
947 t271 = t269 * t270;
948 t276 = 0.1e1 / t24 / t91;
949 t281 = sigma[0] * t260;
950 t285 = sigma[0] * sigma[0];
951 t288 = 0.1e1 / t24 / t258 / t91;
952 t291 = 0.1e1 / t110 / t109;
953 t295 = 0.28e2 / 0.9e1 * t28 * t29 * t276 * t33 + 0.20e2 / 0.3e1 * t28 * t281 * t111 - 0.16e2 / 0.9e1 * t28 * t285 * t288 * t291;
954 t296 = t102 * t295;
955 t300 = 0.176e3 / 0.81e2 * t21 * t22 * t260 * t37 + 0.32e2 / 0.27e2 * t99 * t108 * t116 + 0.4e1 / 0.9e1 * t99 * t100 * t271 - 0.2e1 / 0.9e1 * t99 * t100 * t296;
956 t302 = t7 * t17 * t300;
957 t303 = 0.3e1 / 0.8e1 * t302;
958 t304 = t47 * t47;
959 t305 = 0.1e1 / t304;
960 t306 = t305 * t16;
961 t307 = t125 * t125;
962 t308 = t69 * t307;
963 t310 = t7 * t306 * t308;
964 t311 = t310 / 0.6e1;
965 t312 = t47 * t86;
966 t314 = t7 * t312 * t126;
967 t315 = t314 / 0.3e1;
968 t316 = -t244;
969 t317 = t69 * t316;
970 t319 = t7 * t124 * t317;
971 t320 = t319 / 0.2e1;
972 t321 = t48 * t250;
973 t323 = t7 * t321 * t69;
974 t324 = t323 / 0.12e2;
975 t325 = -t233 - t237 - t240 - t248 + t254 - t257 - t303 - t311 - t315 - t320 + t324;
976 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
977 v2rho2[0] = -t83 - t223 - 0.3e1 / 0.4e1 * t122 - t128 - t225 + t9 * t325;
978
979 t327 = t6 * t227;
980 t328 = t4 * t327;
981 t329 = t16 * t42;
982 t330 = t137 * t80;
983 t332 = t328 * t329 * t330;
984 t335 = t7 * t234 * t138;
985 t337 = t120 * t137;
986 t339 = t7 * t75 * t337;
987 t342 = 0.1e1 / t85 / t76;
988 t343 = t14 * t342;
989 t344 = t42 * t8;
990 t346 = t7 * t343 * t344;
991 t350 = t6 * t305;
992 t351 = t4 * t350;
993 t352 = t16 * t69;
994 t353 = t142 * t125;
995 t355 = t351 * t352 * t353;
996 t358 = t7 * t312 * t143;
997 t360 = t47 * t342;
998 t361 = t69 * t8;
999 t363 = t7 * t360 * t361;
1000 t366 = t174 * t125;
1001 t368 = t7 * t124 * t366;
1002 t371 = t7 * t130 * t174;
1003 t373 = -t332 / 0.6e1 - t335 / 0.6e1 - t339 / 0.2e1 - t346 / 0.2e1 - t236 / 0.6e1 + t254 - t256 / 0.8e1 - t355 / 0.6e1 - t358 / 0.6e1 + t363 / 0.2e1 - t314 / 0.6e1 + t324 - t368 / 0.2e1 - t371 / 0.8e1;
1004 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1005 v2rho2[1] = t9 * t373 - t123 - t129 - t141 - t146 - t177 - t223 - t225 - t84;
1006
1007 t376 = t137 * t137;
1008 t377 = t42 * t376;
1009 t379 = t7 * t228 * t377;
1010 t380 = t379 / 0.6e1;
1011 t381 = t335 / 0.3e1;
1012 t382 = t77 + t243;
1013 t383 = t42 * t382;
1014 t385 = t7 * t75 * t383;
1015 t386 = t385 / 0.2e1;
1016 t387 = t142 * t142;
1017 t388 = t69 * t387;
1018 t390 = t7 * t306 * t388;
1019 t391 = t390 / 0.6e1;
1020 t392 = t358 / 0.3e1;
1021 t393 = t174 * t142;
1022 t395 = t7 * t124 * t393;
1023 t396 = -t382;
1024 t397 = t69 * t396;
1025 t399 = t7 * t124 * t397;
1026 t400 = t399 / 0.2e1;
1027 t401 = t371 / 0.4e1;
1028 t402 = t51 * t51;
1029 t404 = 0.1e1 / t53 / t402;
1030 t413 = 0.1e1 / t155 / t63;
1031 t414 = t169 * t169;
1032 t415 = t413 * t414;
1033 t420 = 0.1e1 / t52 / t147;
1034 t425 = sigma[2] * t404;
1035 t429 = sigma[2] * sigma[2];
1036 t432 = 0.1e1 / t52 / t402 / t147;
1037 t435 = 0.1e1 / t164 / t163;
1038 t439 = 0.28e2 / 0.9e1 * t28 * t56 * t420 * t60 + 0.20e2 / 0.3e1 * t28 * t425 * t165 - 0.16e2 / 0.9e1 * t28 * t429 * t432 * t435;
1039 t440 = t156 * t439;
1040 t444 = 0.176e3 / 0.81e2 * t21 * t50 * t404 * t64 + 0.32e2 / 0.27e2 * t99 * t162 * t170 + 0.4e1 / 0.9e1 * t99 * t154 * t415 - 0.2e1 / 0.9e1 * t99 * t154 * t440;
1041 t446 = t7 * t49 * t444;
1042 t447 = 0.3e1 / 0.8e1 * t446;
1043 t448 = -t380 - t381 - t386 + t254 - t391 - t392 - t395 - t400 + t324 - t401 - t447;
1044 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1045 v2rho2[2] = -t140 - t223 - t145 - t225 - 0.3e1 / 0.4e1 * t176 + t9 * t448;
1046
1047 t451 = t7 * t17 * t199;
1048 t453 = t6 * t14;
1049 t454 = t199 * t80;
1050 t458 = t5 * t93;
1051 t468 = t269 * t194;
1052 t469 = t468 * t115;
1053 t480 = t258 * t23;
1054 t482 = 0.1e1 / t24 / t480;
1055 t483 = t482 * t291;
1056 t487 = -0.2e1 / 0.3e1 * t28 * t187 * t104 * t33 - 0.2e1 * t28 * t93 * t111 + 0.2e1 / 0.3e1 * t28 * t483 * sigma[0];
1057 t488 = t102 * t487;
1058 t492 = -0.16e2 / 0.27e2 * t21 * t458 * t37 - 0.2e1 / 0.9e1 * t21 * t184 * t116 + 0.16e2 / 0.27e2 * t99 * t108 * t195 + 0.4e1 / 0.9e1 * t99 * t100 * t469 - 0.2e1 / 0.9e1 * t99 * t100 * t488;
1059 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1060 v2rhosigma[0] = -t451 / 0.2e1 - t182 * t453 * t454 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t183 * t492;
1061
1062 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1063 v2rhosigma[1] = 0.0e0;
1064
1065 t497 = t7 * t49 * t219;
1066 t498 = t6 * t47;
1067 t499 = t219 * t125;
1068 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1069 v2rhosigma[2] = -t182 * t498 * t499 / 0.2e1 - t497 / 0.2e1;
1070
1071 t504 = t453 * t199 * t137;
1072 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1073 v2rhosigma[3] = -t182 * t504 / 0.2e1 - t451 / 0.2e1;
1074
1075 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1076 v2rhosigma[4] = 0.0e0;
1077
1078 t509 = t498 * t219 * t142;
1079 t512 = t5 * t149;
1080 t522 = t413 * t214;
1081 t523 = t522 * t169;
1082 t534 = t402 * t51;
1083 t536 = 0.1e1 / t52 / t534;
1084 t537 = t536 * t435;
1085 t541 = -0.2e1 / 0.3e1 * t28 * t207 * t158 * t60 - 0.2e1 * t28 * t149 * t165 + 0.2e1 / 0.3e1 * t28 * t537 * sigma[2];
1086 t542 = t156 * t541;
1087 t546 = -0.16e2 / 0.27e2 * t21 * t512 * t64 - 0.2e1 / 0.9e1 * t21 * t204 * t170 + 0.16e2 / 0.27e2 * t99 * t162 * t215 + 0.4e1 / 0.9e1 * t99 * t154 * t523 - 0.2e1 / 0.9e1 * t99 * t154 * t542;
1088 t547 = t203 * t546;
1089 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1090 v2rhosigma[5] = -t497 / 0.2e1 - t182 * t509 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t547;
1091
1092 t553 = t194 * t194;
1093 t554 = t269 * t553;
1094 t559 = 0.1e1 / t29 / sigma[0];
1095 t563 = 0.1e1 / sigma[0];
1096 t567 = t258 * rho[0];
1097 t569 = 0.1e1 / t24 / t567;
1098 t573 = t28 * t563 * t27 * t111 / 0.4e1 - t28 * t559 * t31 * t33 / 0.4e1 - t28 * t569 * t291 / 0.4e1;
1099 t574 = t102 * t573;
1100 t578 = -0.4e1 / 0.9e1 * t21 * t184 * t195 + 0.4e1 / 0.9e1 * t99 * t100 * t554 - 0.2e1 / 0.9e1 * t99 * t100 * t574;
1101 t579 = t183 * t578;
1102 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1103 v2sigma2[0] = -0.3e1 / 0.8e1 * t182 * t579;
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 t585 = t214 * t214;
1118 t586 = t413 * t585;
1119 t591 = 0.1e1 / t56 / sigma[2];
1120 t595 = 0.1e1 / sigma[2];
1121 t599 = t402 * rho[1];
1122 t601 = 0.1e1 / t52 / t599;
1123 t605 = t28 * t595 * t55 * t165 / 0.4e1 - t28 * t591 * t58 * t60 / 0.4e1 - t28 * t601 * t435 / 0.4e1;
1124 t606 = t156 * t605;
1125 t610 = -0.4e1 / 0.9e1 * t21 * t204 * t215 + 0.4e1 / 0.9e1 * t99 * t154 * t586 - 0.2e1 / 0.9e1 * t99 * t154 * t606;
1126 t611 = t203 * t610;
1127 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1128 v2sigma2[5] = -0.3e1 / 0.8e1 * t182 * t611;
1129
1130 #ifndef XC_DONT_COMPILE_KXC
1131
1132 if(order < 3) return;
1133
1134
1135 t618 = t253 / 0.4e1;
1136 t620 = t323 / 0.4e1;
1137 t623 = t120 * t229;
1138 t625 = t7 * t228 * t623;
1139 t626 = t625 / 0.2e1;
1140 t628 = t7 * t234 * t238;
1141 t630 = t7 * t234 * t245;
1142 t632 = t300 * t80;
1143 t634 = t7 * t75 * t632;
1144 t635 = 0.3e1 / 0.2e1 * t634;
1145 t636 = t120 * t244;
1146 t638 = t7 * t75 * t636;
1147 t639 = 0.3e1 / 0.2e1 * t638;
1148 t640 = t76 * t76;
1149 t641 = 0.1e1 / t640;
1150 t642 = t8 * t641;
1151 t644 = 0.3e1 * t242 - 0.3e1 * t642;
1152 t645 = t42 * t644;
1153 t647 = t7 * t75 * t645;
1154 t648 = t647 / 0.2e1;
1155 t650 = t7 * t251 * t120;
1156 t651 = t650 / 0.4e1;
1157 t653 = t7 * t87 * t300;
1158 t656 = 0.1e1 / t25 / t567;
1159 t670 = t101 * t101;
1160 t671 = 0.1e1 / t670;
1161 t672 = t270 * t115;
1162 t673 = t671 * t672;
1163 t677 = t269 * t115;
1164 t678 = t677 * t295;
1165 t683 = 0.1e1 / t24 / t258;
1166 t688 = sigma[0] * t656;
1167 t692 = t258 * t258;
1168 t694 = 0.1e1 / t24 / t692;
1169 t699 = t285 * sigma[0];
1170 t700 = t692 * t91;
1171 t701 = 0.1e1 / t700;
1172 t703 = t109 * t109;
1173 t705 = 0.1e1 / t110 / t703;
1174 t709 = -0.280e3 / 0.27e2 * t28 * t29 * t683 * t33 - 0.952e3 / 0.27e2 * t28 * t688 * t111 + 0.592e3 / 0.27e2 * t28 * t285 * t694 * t291 - 0.64e2 / 0.9e1 * t28 * t699 * t701 * t705;
1175 t710 = t102 * t709;
1176 t714 = -0.2464e4 / 0.243e3 * t21 * t22 * t656 * t37 - 0.176e3 / 0.27e2 * t99 * t281 * t116 - 0.32e2 / 0.9e1 * t99 * t108 * t271 + 0.16e2 / 0.9e1 * t99 * t108 * t296 - 0.4e1 / 0.3e1 * t99 * t100 * t673 + 0.4e1 / 0.3e1 * t99 * t100 * t678 - 0.2e1 / 0.9e1 * t99 * t100 * t710;
1177 t716 = t7 * t17 * t714;
1178 t717 = 0.3e1 / 0.8e1 * t716;
1179 t719 = 0.1e1 / t226 / t13;
1180 t720 = t719 * t16;
1181 t721 = t229 * t80;
1182 t722 = t42 * t721;
1183 t724 = t7 * t720 * t722;
1184 t725 = t724 / 0.9e1;
1185 t727 = t227 * t86;
1186 t729 = t7 * t727 * t230;
1187 t731 = t80 * t244;
1188 t733 = t328 * t329 * t731;
1189 t734 = t733 / 0.2e1;
1190 t735 = t14 * t250;
1191 t737 = t7 * t735 * t81;
1192 t738 = t737 / 0.3e1;
1193 t740 = 0.1e1 / t304 / t46;
1194 t741 = t740 * t16;
1195 t742 = t307 * t125;
1196 t743 = t69 * t742;
1197 t745 = t7 * t741 * t743;
1198 t746 = t745 / 0.9e1;
1199 t747 = t305 * t86;
1200 t749 = t7 * t747 * t308;
1201 t751 = t125 * t316;
1202 t753 = t351 * t352 * t751;
1203 t754 = t753 / 0.2e1;
1204 t755 = t47 * t250;
1205 t757 = t7 * t755 * t126;
1206 t758 = t757 / 0.3e1;
1207 t759 = t15 * t342;
1208 t761 = t7 * t759 * t42;
1209 t762 = 0.5e1 / 0.36e2 * t761;
1210 t764 = t7 * t312 * t317;
1211 t766 = -t644;
1212 t767 = t69 * t766;
1213 t769 = t7 * t124 * t767;
1214 t770 = t769 / 0.2e1;
1215 t771 = t48 * t342;
1216 t773 = t7 * t771 * t69;
1217 t774 = 0.5e1 / 0.36e2 * t773;
1218 t775 = -t729 / 0.6e1 - t734 + t738 + t746 - t749 / 0.6e1 - t754 + t758 - t762 - t764 / 0.2e1 - t770 - t774;
1219 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1220 v3rho3[0] = -t232 / 0.2e1 - t236 - t310 / 0.2e1 - t314 - 0.3e1 * t240 - 0.3e1 / 0.2e1 * t247 + t618 - 0.3e1 / 0.2e1 * t319 + t620 - 0.3e1 / 0.4e1 * t256 - 0.9e1 / 0.8e1 * t302 + t9 * (-t626 - t628 - t630 / 0.2e1 - t635 - t639 - t648 + t651 - 0.3e1 / 0.8e1 * t653 - t717 + t725 + t775);
1221
1222 t780 = t332 / 0.3e1;
1223 t781 = t355 / 0.3e1;
1224 t784 = t137 * t244;
1225 t786 = t328 * t329 * t784;
1226 t788 = t142 * t316;
1227 t790 = t351 * t352 * t788;
1228 t792 = t342 * t69;
1229 t793 = t8 * t125;
1230 t795 = t351 * t792 * t793;
1231 t797 = t6 * t719;
1232 t798 = t4 * t797;
1233 t799 = t137 * t229;
1234 t801 = t798 * t329 * t799;
1235 t803 = t86 * t42;
1236 t805 = t328 * t803 * t330;
1237 t806 = t805 / 0.9e1;
1238 t807 = t342 * t42;
1239 t808 = t8 * t80;
1240 t810 = t328 * t807 * t808;
1241 t812 = t6 * t740;
1242 t813 = t4 * t812;
1243 t814 = t142 * t307;
1244 t816 = t813 * t352 * t814;
1245 t818 = t86 * t69;
1246 t820 = t351 * t818 * t353;
1247 t821 = t820 / 0.9e1;
1248 t822 = t16 * t120;
1249 t824 = t328 * t822 * t330;
1250 t827 = t7 * t755 * t143;
1251 t829 = t174 * t307;
1252 t831 = t7 * t306 * t829;
1253 t834 = t7 * t312 * t366;
1254 t835 = t834 / 0.3e1;
1255 t837 = t7 * t234 * t337;
1256 t838 = t837 / 0.3e1;
1257 t839 = t300 * t137;
1258 t841 = t7 * t75 * t839;
1259 t844 = 0.1e1 / t85 / t241;
1260 t845 = t47 * t844;
1261 t847 = t7 * t845 * t361;
1262 t849 = t174 * t316;
1263 t851 = t7 * t124 * t849;
1264 t853 = -t786 / 0.6e1 - t790 / 0.6e1 + t795 / 0.3e1 + t801 / 0.9e1 - t806 - t810 / 0.3e1 + t816 / 0.9e1 - t821 - t824 / 0.3e1 + t827 / 0.9e1 - t831 / 0.6e1 - t835 - t838 - t841 / 0.2e1 - 0.7e1 / 0.6e1 * t847 - t851 / 0.2e1;
1265 t855 = t7 * t735 * t138;
1266 t857 = t14 * t844;
1267 t859 = t7 * t857 * t344;
1268 t861 = t120 * t8;
1269 t863 = t7 * t343 * t861;
1270 t872 = t7 * t343 * t42;
1271 t875 = t7 * t360 * t69;
1272 t878 = t7 * t321 * t174;
1273 t882 = t855 / 0.9e1 + 0.7e1 / 0.6e1 * t859 - t863 - t729 / 0.18e2 + 0.2e1 / 0.9e1 * t737 - t749 / 0.18e2 + 0.2e1 / 0.9e1 * t757 - t764 / 0.6e1 - t628 / 0.3e1 - t630 / 0.6e1 - t872 / 0.2e1 + t875 / 0.2e1 + t878 / 0.12e2 - t762 - t774 + t650 / 0.6e1 - t653 / 0.8e1;
1274 t885 = -t248 + t618 - t320 + t620 - t339 + t363 - t368 - t401 - t256 / 0.2e1 - t303 + t9 * (t853 + t882);
1275 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1276 v3rho3[1] = -t233 - 0.2e1 / 0.3e1 * t236 - t311 - 0.2e1 / 0.3e1 * t314 - t780 - t381 - t346 - t781 - t392 - t240 + t885;
1277
1278 t890 = t16 * t174;
1279 t892 = t351 * t890 * t353;
1280 t894 = t396 * t125;
1281 t896 = t351 * t352 * t894;
1282 t898 = t376 * t80;
1283 t900 = t798 * t329 * t898;
1284 t902 = t382 * t80;
1285 t904 = t328 * t329 * t902;
1286 t906 = t387 * t125;
1287 t908 = t813 * t352 * t906;
1288 t910 = t142 * t8;
1289 t912 = t351 * t792 * t910;
1290 t914 = t137 * t8;
1291 t916 = t328 * t807 * t914;
1292 t919 = t242 + 0.3e1 * t642;
1293 t920 = t69 * t919;
1294 t922 = t7 * t124 * t920;
1295 t924 = t444 * t125;
1296 t926 = t7 * t124 * t924;
1297 t929 = t7 * t727 * t377;
1298 t932 = t7 * t234 * t383;
1299 t935 = t7 * t747 * t388;
1300 t938 = t7 * t312 * t393;
1301 t941 = t7 * t312 * t397;
1302 t943 = -t892 / 0.3e1 - t896 / 0.6e1 + t900 / 0.9e1 - t904 / 0.6e1 + t908 / 0.9e1 + t912 / 0.3e1 - t916 / 0.3e1 - t806 - t821 - t922 / 0.2e1 - t926 / 0.2e1 - t929 / 0.18e2 - t932 / 0.6e1 - t935 / 0.18e2 - t938 / 0.3e1 - t941 / 0.6e1;
1303 t944 = t120 * t376;
1304 t946 = t7 * t228 * t944;
1305 t948 = t120 * t382;
1306 t950 = t7 * t75 * t948;
1307 t952 = -t919;
1308 t953 = t42 * t952;
1309 t955 = t7 * t75 * t953;
1310 t957 = t174 * t8;
1311 t959 = t7 * t360 * t957;
1312 t967 = t7 * t130 * t444;
1313 t971 = -t946 / 0.6e1 - t950 / 0.2e1 - t955 / 0.2e1 + t959 + 0.2e1 / 0.9e1 * t827 - t835 - t838 + t847 / 0.3e1 + 0.2e1 / 0.9e1 * t855 - t859 / 0.3e1 + t737 / 0.9e1 + t757 / 0.9e1 - t967 / 0.8e1 + t878 / 0.6e1 - t762 - t774 + t650 / 0.12e2;
1314 t974 = t620 - t339 + t363 - t368 - t371 / 0.2e1 - t386 - t395 - t400 - t257 - t447 + t9 * (t943 + t971);
1315 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1316 v3rho3[2] = -t237 - t315 - t780 - 0.2e1 / 0.3e1 * t335 - t346 - t781 - 0.2e1 / 0.3e1 * t358 - t380 - t391 + t618 + t974;
1317
1318 t982 = t855 / 0.3e1;
1319 t983 = t827 / 0.3e1;
1320 t984 = t878 / 0.4e1;
1321 t987 = -0.3e1 * t242 - 0.3e1 * t642;
1322 t988 = t42 * t987;
1323 t990 = t7 * t75 * t988;
1324 t991 = t990 / 0.2e1;
1325 t992 = t174 * t387;
1326 t994 = t7 * t306 * t992;
1327 t995 = t994 / 0.2e1;
1328 t996 = t444 * t142;
1329 t998 = t7 * t124 * t996;
1330 t999 = 0.3e1 / 0.2e1 * t998;
1331 t1000 = t174 * t396;
1332 t1002 = t7 * t124 * t1000;
1333 t1003 = 0.3e1 / 0.2e1 * t1002;
1334 t1005 = -t987;
1335 t1006 = t69 * t1005;
1336 t1008 = t7 * t124 * t1006;
1337 t1009 = t1008 / 0.2e1;
1338 t1014 = t376 * t137;
1339 t1015 = t42 * t1014;
1340 t1017 = t7 * t720 * t1015;
1341 t1018 = t1017 / 0.9e1;
1342 t1019 = t137 * t382;
1343 t1021 = t328 * t329 * t1019;
1344 t1022 = t1021 / 0.2e1;
1345 t1023 = t387 * t142;
1346 t1024 = t69 * t1023;
1347 t1026 = t7 * t741 * t1024;
1348 t1027 = t1026 / 0.9e1;
1349 t1028 = t142 * t396;
1350 t1030 = t351 * t352 * t1028;
1351 t1031 = t1030 / 0.2e1;
1352 t1033 = 0.1e1 / t53 / t599;
1353 t1047 = t155 * t155;
1354 t1048 = 0.1e1 / t1047;
1355 t1049 = t414 * t169;
1356 t1050 = t1048 * t1049;
1357 t1054 = t413 * t169;
1358 t1055 = t1054 * t439;
1359 t1060 = 0.1e1 / t52 / t402;
1360 t1065 = sigma[2] * t1033;
1361 t1069 = t402 * t402;
1362 t1071 = 0.1e1 / t52 / t1069;
1363 t1076 = t429 * sigma[2];
1364 t1077 = t1069 * t147;
1365 t1078 = 0.1e1 / t1077;
1366 t1080 = t163 * t163;
1367 t1082 = 0.1e1 / t164 / t1080;
1368 t1086 = -0.280e3 / 0.27e2 * t28 * t56 * t1060 * t60 - 0.952e3 / 0.27e2 * t28 * t1065 * t165 + 0.592e3 / 0.27e2 * t28 * t429 * t1071 * t435 - 0.64e2 / 0.9e1 * t28 * t1076 * t1078 * t1082;
1369 t1087 = t156 * t1086;
1370 t1091 = -0.2464e4 / 0.243e3 * t21 * t50 * t1033 * t64 - 0.176e3 / 0.27e2 * t99 * t425 * t170 - 0.32e2 / 0.9e1 * t99 * t162 * t415 + 0.16e2 / 0.9e1 * t99 * t162 * t440 - 0.4e1 / 0.3e1 * t99 * t154 * t1050 + 0.4e1 / 0.3e1 * t99 * t154 * t1055 - 0.2e1 / 0.9e1 * t99 * t154 * t1087;
1371 t1093 = t7 * t49 * t1091;
1372 t1094 = 0.3e1 / 0.8e1 * t1093;
1373 t1095 = -t1009 - t929 / 0.6e1 - t932 / 0.2e1 - t935 / 0.6e1 - t938 - t941 / 0.2e1 + t1018 - t1022 + t1027 - t1031 - t1094;
1374 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1375 v3rho3[3] = -t335 - t358 - t379 / 0.2e1 - t390 / 0.2e1 + t618 + t620 - 0.3e1 / 0.4e1 * t371 - 0.3e1 / 0.2e1 * t385 - 0.3e1 * t395 - 0.3e1 / 0.2e1 * t399 - 0.9e1 / 0.8e1 * t446 + t9 * (t982 + t983 - t762 - t774 + t984 - 0.3e1 / 0.8e1 * t967 - t991 - t995 - t999 - t1003 + t1095);
1376
1377 t1099 = t7 * t75 * t454;
1378 t1103 = t7 * t87 * t199 / 0.6e1;
1379 t1105 = t7 * t17 * t492;
1380 t1106 = t199 * t229;
1381 t1110 = t492 * t80;
1382 t1113 = t199 * t244;
1383 t1117 = t5 * t260;
1384 t1140 = t671 * t194 * t270;
1385 t1144 = t269 * t487;
1386 t1145 = t1144 * t115;
1387 t1149 = t468 * t295;
1388 t1160 = t288 * t291;
1389 t1164 = t692 * t23;
1390 t1166 = 0.1e1 / t1164 * t705;
1391 t1170 = 0.14e2 / 0.9e1 * t28 * t187 * t276 * t33 + 0.74e2 / 0.9e1 * t28 * t260 * t111 - 0.62e2 / 0.9e1 * t28 * t1160 * sigma[0] + 0.8e1 / 0.3e1 * t28 * t1166 * t285;
1392 t1171 = t102 * t1170;
1393 t1175 = 0.176e3 / 0.81e2 * t21 * t1117 * t37 + 0.32e2 / 0.27e2 * t21 * t458 * t116 + 0.4e1 / 0.9e1 * t21 * t184 * t271 - 0.2e1 / 0.9e1 * t21 * t184 * t296 - 0.176e3 / 0.81e2 * t99 * t281 * t195 - 0.64e2 / 0.27e2 * t99 * t108 * t469 + 0.32e2 / 0.27e2 * t99 * t108 * t488 - 0.4e1 / 0.3e1 * t99 * t100 * t1140 + 0.8e1 / 0.9e1 * t99 * t100 * t1145 + 0.4e1 / 0.9e1 * t99 * t100 * t1149 - 0.2e1 / 0.9e1 * t99 * t100 * t1171;
1394 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1395 v3rho2sigma[0] = -0.4e1 / 0.3e1 * t1099 - t1103 - t1105 - t182 * t327 * t1106 / 0.6e1 - t182 * t453 * t1110 - t182 * t453 * t1113 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t183 * t1175;
1396
1397 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1398 v3rho2sigma[1] = 0.0e0;
1399
1400 t1180 = t7 * t124 * t499;
1401 t1184 = t7 * t130 * t219 / 0.6e1;
1402 t1185 = t219 * t307;
1403 t1189 = t219 * t316;
1404 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1405 v3rho2sigma[2] = -0.4e1 / 0.3e1 * t1180 - t1184 - t182 * t350 * t1185 / 0.6e1 - t182 * t498 * t1189 / 0.2e1;
1406
1407 t1195 = t16 * t1;
1408 t1196 = t1195 * t3;
1409 t1197 = t1196 * t504;
1410 t1199 = t3 * t6;
1411 t1200 = t181 * t1199;
1412 t1201 = t227 * t199;
1413 t1202 = t1201 * t330;
1414 t1206 = t453 * t492 * t137;
1415 t1209 = t250 * t1;
1416 t1210 = t1209 * t3;
1417 t1212 = t453 * t199 * t8;
1418 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1419 v3rho2sigma[3] = -0.2e1 / 0.3e1 * t1099 - t1103 - t1105 / 0.2e1 - 0.2e1 / 0.3e1 * t1197 - t1200 * t1202 / 0.6e1 - t182 * t1206 / 0.2e1 - t1210 * t1212 / 0.2e1;
1420
1421 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1422 v3rho2sigma[4] = 0.0e0;
1423
1424 t1216 = t1196 * t509;
1425 t1218 = t305 * t219;
1426 t1219 = t1218 * t353;
1427 t1223 = t498 * t219 * t8;
1428 t1226 = t1196 * t547;
1429 t1229 = t498 * t546 * t125;
1430 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1431 v3rho2sigma[5] = -0.2e1 / 0.3e1 * t1180 - t1184 - 0.2e1 / 0.3e1 * t1216 - t1200 * t1219 / 0.6e1 + t1210 * t1223 / 0.2e1 - t1226 / 0.2e1 - t182 * t1229 / 0.2e1;
1432
1433 t1234 = t327 * t199 * t376;
1434 t1238 = t453 * t199 * t382;
1435 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1436 v3rho2sigma[6] = -0.4e1 / 0.3e1 * t1197 - t1103 - t182 * t1234 / 0.6e1 - t182 * t1238 / 0.2e1;
1437
1438 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1439 v3rho2sigma[7] = 0.0e0;
1440
1441 t1243 = t350 * t219 * t387;
1442 t1247 = t498 * t546 * t142;
1443 t1250 = t498 * t219 * t396;
1444 t1253 = t5 * t404;
1445 t1276 = t1048 * t214 * t414;
1446 t1280 = t413 * t541;
1447 t1281 = t1280 * t169;
1448 t1285 = t522 * t439;
1449 t1296 = t432 * t435;
1450 t1300 = t1069 * t51;
1451 t1302 = 0.1e1 / t1300 * t1082;
1452 t1306 = 0.14e2 / 0.9e1 * t28 * t207 * t420 * t60 + 0.74e2 / 0.9e1 * t28 * t404 * t165 - 0.62e2 / 0.9e1 * t28 * t1296 * sigma[2] + 0.8e1 / 0.3e1 * t28 * t1302 * t429;
1453 t1307 = t156 * t1306;
1454 t1311 = 0.176e3 / 0.81e2 * t21 * t1253 * t64 + 0.32e2 / 0.27e2 * t21 * t512 * t170 + 0.4e1 / 0.9e1 * t21 * t204 * t415 - 0.2e1 / 0.9e1 * t21 * t204 * t440 - 0.176e3 / 0.81e2 * t99 * t425 * t215 - 0.64e2 / 0.27e2 * t99 * t162 * t523 + 0.32e2 / 0.27e2 * t99 * t162 * t542 - 0.4e1 / 0.3e1 * t99 * t154 * t1276 + 0.8e1 / 0.9e1 * t99 * t154 * t1281 + 0.4e1 / 0.9e1 * t99 * t154 * t1285 - 0.2e1 / 0.9e1 * t99 * t154 * t1307;
1455 t1312 = t203 * t1311;
1456 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1457 v3rho2sigma[8] = -0.4e1 / 0.3e1 * t1216 - t1184 - t1226 - t182 * t1243 / 0.6e1 - t182 * t1247 - t182 * t1250 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t1312;
1458
1459 t1315 = t1196 * t579;
1460 t1318 = t453 * t578 * t80;
1461 t1324 = t27 * t269;
1462 t1325 = t194 * t115;
1463 t1335 = t671 * t553;
1464 t1336 = t1335 * t115;
1465 t1340 = t468 * t487;
1466 t1347 = t269 * t573;
1467 t1348 = t1347 * t115;
1468 t1362 = t692 * rho[0];
1469 t1364 = 0.1e1 / t1362 * t705;
1470 t1367 = t28 * t559 * t104 * t33 / 0.3e1 - t28 * t563 * t93 * t111 / 0.3e1 + 0.5e1 / 0.3e1 * t28 * t483 - t28 * t1364 * sigma[0];
1471 t1368 = t102 * t1367;
1472 t1372 = 0.32e2 / 0.27e2 * t21 * t458 * t195 + 0.8e1 / 0.9e1 * t99 * t1324 * t1325 - 0.4e1 / 0.9e1 * t21 * t184 * t488 - 0.32e2 / 0.27e2 * t99 * t108 * t554 - 0.4e1 / 0.3e1 * t99 * t100 * t1336 + 0.8e1 / 0.9e1 * t99 * t100 * t1340 + 0.16e2 / 0.27e2 * t99 * t108 * t574 + 0.4e1 / 0.9e1 * t99 * t100 * t1348 - 0.2e1 / 0.9e1 * t99 * t100 * t1368;
1473 t1373 = t183 * t1372;
1474 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1475 v3rhosigma2[0] = -t1315 / 0.2e1 - t182 * t1318 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t1373;
1476
1477 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1478 v3rhosigma2[1] = 0.0e0;
1479
1480 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1481 v3rhosigma2[2] = 0.0e0;
1482
1483 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1484 v3rhosigma2[3] = 0.0e0;
1485
1486 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1487 v3rhosigma2[4] = 0.0e0;
1488
1489 t1376 = t1196 * t611;
1490 t1378 = t498 * t610 * t125;
1491 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1492 v3rhosigma2[5] = -t182 * t1378 / 0.2e1 - t1376 / 0.2e1;
1493
1494 t1382 = t453 * t578 * t137;
1495 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1496 v3rhosigma2[6] = -t182 * t1382 / 0.2e1 - t1315 / 0.2e1;
1497
1498 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1499 v3rhosigma2[7] = 0.0e0;
1500
1501 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1502 v3rhosigma2[8] = 0.0e0;
1503
1504 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1505 v3rhosigma2[9] = 0.0e0;
1506
1507 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1508 v3rhosigma2[10] = 0.0e0;
1509
1510 t1387 = t498 * t610 * t142;
1511 t1393 = t55 * t413;
1512 t1394 = t214 * t169;
1513 t1404 = t1048 * t585;
1514 t1405 = t1404 * t169;
1515 t1409 = t522 * t541;
1516 t1416 = t413 * t605;
1517 t1417 = t1416 * t169;
1518 t1431 = t1069 * rho[1];
1519 t1433 = 0.1e1 / t1431 * t1082;
1520 t1436 = t28 * t591 * t158 * t60 / 0.3e1 - t28 * t595 * t149 * t165 / 0.3e1 + 0.5e1 / 0.3e1 * t28 * t537 - t28 * t1433 * sigma[2];
1521 t1437 = t156 * t1436;
1522 t1441 = 0.32e2 / 0.27e2 * t21 * t512 * t215 + 0.8e1 / 0.9e1 * t99 * t1393 * t1394 - 0.4e1 / 0.9e1 * t21 * t204 * t542 - 0.32e2 / 0.27e2 * t99 * t162 * t586 - 0.4e1 / 0.3e1 * t99 * t154 * t1405 + 0.8e1 / 0.9e1 * t99 * t154 * t1409 + 0.16e2 / 0.27e2 * t99 * t162 * t606 + 0.4e1 / 0.9e1 * t99 * t154 * t1417 - 0.2e1 / 0.9e1 * t99 * t154 * t1437;
1523 t1442 = t203 * t1441;
1524 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1525 v3rhosigma2[11] = -t1376 / 0.2e1 - t182 * t1387 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t1442;
1526
1527 t1451 = t553 * t194;
1528 t1452 = t671 * t1451;
1529 t1456 = t468 * t573;
1530 t1461 = 0.1e1 / t29 / t285;
1531 t1466 = 0.1e1 / t285;
1532 t1475 = 0.1e1 / t692;
1533 t1479 = 0.3e1 / 0.8e1 * t28 * t1461 * t31 * t33 - 0.3e1 / 0.8e1 * t28 * t1466 * t27 * t111 - t28 * t563 * t569 * t291 / 0.8e1 + 0.3e1 / 0.8e1 * t28 * t1475 * t705;
1534 t1480 = t102 * t1479;
1535 t1484 = 0.4e1 / 0.3e1 * t21 * t184 * t554 - 0.2e1 / 0.3e1 * t21 * t184 * t574 - 0.4e1 / 0.3e1 * t99 * t100 * t1452 + 0.4e1 / 0.3e1 * t99 * t100 * t1456 - 0.2e1 / 0.9e1 * t99 * t100 * t1480;
1536 t1485 = t183 * t1484;
1537 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1538 v3sigma3[0] = -0.3e1 / 0.8e1 * t182 * t1485;
1539
1540 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1541 v3sigma3[1] = 0.0e0;
1542
1543 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1544 v3sigma3[2] = 0.0e0;
1545
1546 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1547 v3sigma3[3] = 0.0e0;
1548
1549 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1550 v3sigma3[4] = 0.0e0;
1551
1552 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1553 v3sigma3[5] = 0.0e0;
1554
1555 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1556 v3sigma3[6] = 0.0e0;
1557
1558 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1559 v3sigma3[7] = 0.0e0;
1560
1561 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1562 v3sigma3[8] = 0.0e0;
1563
1564 t1494 = t585 * t214;
1565 t1495 = t1048 * t1494;
1566 t1499 = t522 * t605;
1567 t1504 = 0.1e1 / t56 / t429;
1568 t1509 = 0.1e1 / t429;
1569 t1518 = 0.1e1 / t1069;
1570 t1522 = 0.3e1 / 0.8e1 * t28 * t1504 * t58 * t60 - 0.3e1 / 0.8e1 * t28 * t1509 * t55 * t165 - t28 * t595 * t601 * t435 / 0.8e1 + 0.3e1 / 0.8e1 * t28 * t1518 * t1082;
1571 t1523 = t156 * t1522;
1572 t1527 = 0.4e1 / 0.3e1 * t21 * t204 * t586 - 0.2e1 / 0.3e1 * t21 * t204 * t606 - 0.4e1 / 0.3e1 * t99 * t154 * t1495 + 0.4e1 / 0.3e1 * t99 * t154 * t1499 - 0.2e1 / 0.9e1 * t99 * t154 * t1523;
1573 t1528 = t203 * t1527;
1574 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1575 v3sigma3[9] = -0.3e1 / 0.8e1 * t182 * t1528;
1576
1577 #ifndef XC_DONT_COMPILE_LXC
1578
1579 if(order < 4) return;
1580
1581
1582 t1532 = t328 * t803 * t731;
1583 t1542 = 0.1e1 / t25 / t480;
1584 t1566 = 0.1e1 / t670 / t36;
1585 t1567 = t270 * t270;
1586 t1577 = t295 * t295;
1587 t1600 = t692 * t258;
1588 t1606 = t285 * t285;
1589 t1613 = 0.1e1 / t110 / t703 / t109;
1590 t1622 = 0.41888e5 / 0.729e3 * t21 * t22 * t1542 * t37 + 0.9856e4 / 0.243e3 * t99 * t688 * t116 + 0.704e3 / 0.27e2 * t99 * t281 * t271 - 0.352e3 / 0.27e2 * t99 * t281 * t296 + 0.128e3 / 0.9e1 * t99 * t108 * t673 - 0.128e3 / 0.9e1 * t99 * t108 * t678 + 0.64e2 / 0.27e2 * t99 * t108 * t710 + 0.16e2 / 0.3e1 * t99 * t100 * t1566 * t1567 - 0.8e1 * t99 * t100 * t671 * t270 * t295 + 0.4e1 / 0.3e1 * t99 * t100 * t269 * t1577 + 0.16e2 / 0.9e1 * t99 * t100 * t677 * t709 - 0.2e1 / 0.9e1 * t99 * t100 * t102 * (0.3640e4 / 0.81e2 * t28 * t29 * t569 * t33 + 0.5768e4 / 0.27e2 * t28 * sigma[0] * t1542 * t111 - 0.18608e5 / 0.81e2 * t28 * t285 / t24 / t1362 * t291 + 0.4480e4 / 0.27e2 * t28 * t699 / t1600 * t705 - 0.1280e4 / 0.27e2 * t28 * t1606 / t25 / t692 / t480 * t1613);
1591 t1627 = t7 * t87 * t714;
1592 t1630 = t7 * t759 * t120;
1593 t1635 = 0.10e2 / 0.27e2 * t7 * t15 * t844 * t42;
1594 t1639 = 0.10e2 / 0.27e2 * t7 * t48 * t844 * t69;
1595 t1641 = t7 * t251 * t300;
1596 t1648 = t351 * t818 * t751;
1597 t1662 = t640 * t9;
1598 t1664 = t8 / t1662;
1599 t1666 = -0.12e2 * t641 + 0.12e2 * t1664;
1600 t1672 = t7 * t755 * t317;
1601 t1675 = t7 * t312 * t767;
1602 t1686 = -0.2e1 / 0.3e1 * t1532 - 0.2e1 * t328 * t822 * t731 - 0.2e1 / 0.3e1 * t328 * t329 * t644 * t80 - 0.3e1 / 0.8e1 * t7 * t17 * t1622 - t1627 / 0.2e1 - 0.5e1 / 0.9e1 * t1630 + t1635 + t1639 + t1641 / 0.2e1 + 0.2e1 / 0.3e1 * t813 * t352 * t307 * t316 - 0.2e1 / 0.3e1 * t1648 - 0.2e1 / 0.3e1 * t351 * t352 * t125 * t766 + 0.2e1 / 0.3e1 * t798 * t329 * t229 * t244 - 0.2e1 * t7 * t75 * t120 * t644 - t7 * t75 * t42 * t1666 / 0.2e1 + 0.2e1 / 0.3e1 * t1672 - 0.2e1 / 0.3e1 * t1675 + t7 * t124 * t69 * t1666 / 0.2e1 + 0.4e1 / 0.9e1 * t7 * t720 * t120 * t721;
1603 t1688 = t7 * t727 * t623;
1604 t1691 = t7 * t735 * t238;
1605 t1694 = t7 * t735 * t245;
1606 t1696 = t46 * t46;
1607 t1698 = 0.1e1 / t304 / t1696;
1608 t1699 = t1698 * t16;
1609 t1700 = t307 * t307;
1610 t1705 = t740 * t86;
1611 t1707 = t7 * t1705 * t743;
1612 t1709 = t305 * t250;
1613 t1711 = t7 * t1709 * t308;
1614 t1713 = t316 * t316;
1615 t1718 = t227 * t250;
1616 t1720 = t7 * t1718 * t230;
1617 t1722 = t244 * t244;
1618 t1727 = t13 * t13;
1619 t1729 = 0.1e1 / t226 / t1727;
1620 t1730 = t1729 * t16;
1621 t1731 = t229 * t229;
1622 t1736 = t719 * t86;
1623 t1738 = t7 * t1736 * t722;
1624 t1744 = t7 * t234 * t632;
1625 t1747 = t7 * t234 * t636;
1626 t1750 = t7 * t234 * t645;
1627 t1761 = t7 * t360 * t126;
1628 t1764 = t7 * t343 * t81;
1629 t1766 = -0.2e1 / 0.3e1 * t1688 + 0.4e1 / 0.3e1 * t1691 + 0.2e1 / 0.3e1 * t1694 - 0.5e1 / 0.27e2 * t7 * t1699 * t69 * t1700 + 0.4e1 / 0.27e2 * t1707 + 0.2e1 / 0.9e1 * t1711 - t7 * t306 * t69 * t1713 / 0.2e1 + 0.2e1 / 0.9e1 * t1720 - t7 * t228 * t42 * t1722 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1730 * t42 * t1731 + 0.4e1 / 0.27e2 * t1738 - t7 * t228 * t300 * t229 - 0.2e1 * t1744 - 0.2e1 * t1747 - 0.2e1 / 0.3e1 * t1750 - 0.2e1 * t7 * t75 * t714 * t80 - 0.3e1 * t7 * t75 * t300 * t244 - 0.20e2 / 0.27e2 * t1761 - 0.20e2 / 0.27e2 * t1764;
1630 t1778 = t9 * (t1686 + t1766) - 0.2e1 * t625 - 0.4e1 * t628 - 0.2e1 * t630 - 0.6e1 * t634 - 0.6e1 * t638 - 0.2e1 * t647 + t650 - 0.3e1 / 0.2e1 * t653 - 0.3e1 / 0.2e1 * t716 + 0.4e1 / 0.9e1 * t724;
1631 t1786 = 0.5e1 / 0.9e1 * t761;
1632 t1789 = 0.5e1 / 0.9e1 * t773;
1633 t1790 = -0.2e1 / 0.3e1 * t729 - 0.2e1 * t733 + 0.4e1 / 0.3e1 * t737 + 0.4e1 / 0.9e1 * t745 - 0.2e1 / 0.3e1 * t749 - 0.2e1 * t753 + 0.4e1 / 0.3e1 * t757 - t1786 - 0.2e1 * t764 - 0.2e1 * t769 - t1789;
1634 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1635 v4rho4[0] = t1778 + t1790;
1636
1637 t1792 = t7 * t735 * t337;
1638 t1793 = t1792 / 0.3e1;
1639 t1795 = t7 * t857 * t861;
1640 t1798 = t7 * t360 * t143;
1641 t1805 = t7 * t747 * t829;
1642 t1808 = t7 * t755 * t366;
1643 t1809 = t1808 / 0.3e1;
1644 t1811 = 0.1e1 / t85 / t640;
1645 t1814 = t7 * t14 * t1811 * t344;
1646 t1817 = t7 * t343 * t138;
1647 t1819 = t227 * t342;
1648 t1828 = t7 * t234 * t839;
1649 t1840 = t7 * t47 * t1811 * t361;
1650 t1843 = t7 * t312 * t849;
1651 t1845 = t305 * t342;
1652 t1849 = t1793 + 0.7e1 / 0.2e1 * t1795 - 0.5e1 / 0.27e2 * t1798 + t7 * t741 * t174 * t742 / 0.9e1 - t1805 / 0.6e1 + t1809 - 0.25e2 / 0.6e1 * t1814 - 0.5e1 / 0.27e2 * t1817 - t7 * t1819 * t81 / 0.2e1 - t7 * t124 * t174 * t766 / 0.2e1 - t1828 / 0.2e1 - t7 * t75 * t714 * t137 / 0.2e1 - 0.3e1 / 0.2e1 * t7 * t343 * t300 * t8 + 0.25e2 / 0.6e1 * t1840 - t1843 / 0.2e1 + t7 * t1845 * t126 / 0.2e1;
1653 t1851 = t4 * t6 * t1729;
1654 t1857 = t798 * t803 * t799;
1655 t1865 = t328 * t250 * t42 * t330;
1656 t1866 = t1865 / 0.9e1;
1657 t1867 = t844 * t42;
1658 t1869 = t328 * t1867 * t808;
1659 t1872 = t328 * t803 * t784;
1660 t1874 = t342 * t120;
1661 t1886 = t351 * t818 * t788;
1662 t1888 = t844 * t69;
1663 t1890 = t351 * t1888 * t793;
1664 t1897 = t328 * t86 * t120 * t330;
1665 t1909 = -0.5e1 / 0.27e2 * t1851 * t329 * t137 * t721 + t1857 / 0.9e1 + t798 * t807 * t8 * t229 / 0.3e1 + t1866 + 0.7e1 / 0.6e1 * t1869 - t1872 / 0.6e1 - t328 * t1874 * t808 - t328 * t807 * t8 * t244 / 0.2e1 - t813 * t792 * t8 * t307 / 0.3e1 - t1886 / 0.6e1 - 0.7e1 / 0.6e1 * t1890 + t798 * t822 * t799 / 0.3e1 - t1897 / 0.3e1 - t351 * t890 * t751 / 0.2e1 - t328 * t822 * t784 / 0.2e1 - t328 * t329 * t137 * t644 / 0.6e1;
1666 t1916 = t4 * t6 * t1698;
1667 t1922 = t813 * t818 * t814;
1668 t1926 = t351 * t250 * t69 * t353;
1669 t1927 = t1926 / 0.9e1;
1670 t1938 = t7 * t845 * t69;
1671 t1941 = t7 * t857 * t42;
1672 t1947 = t7 * t771 * t174;
1673 t1952 = -t328 * t16 * t300 * t330 / 0.2e1 - 0.5e1 / 0.27e2 * t1916 * t352 * t142 * t742 + t1922 / 0.9e1 + t1927 - t351 * t352 * t142 * t766 / 0.6e1 + t351 * t792 * t8 * t316 / 0.2e1 - t1532 / 0.6e1 - 0.5e1 / 0.2e1 * t1938 + 0.5e1 / 0.2e1 * t1941 - 0.3e1 / 0.2e1 * t7 * t343 * t120 - 0.5e1 / 0.36e2 * t1947 - t1627 / 0.8e1 - 0.5e1 / 0.12e2 * t1630 + t1635 + t1639 + t1641 / 0.4e1;
1674 t1976 = -t1648 / 0.6e1 + t1672 / 0.3e1 - t1675 / 0.6e1 - t1688 / 0.6e1 + 0.2e1 / 0.3e1 * t1691 + t1694 / 0.3e1 + t1707 / 0.27e2 + t1711 / 0.9e1 + t1720 / 0.9e1 + t1738 / 0.27e2 - t1744 / 0.2e1 - t1747 / 0.2e1 - t1750 / 0.6e1 - 0.5e1 / 0.9e1 * t1761 - 0.5e1 / 0.9e1 * t1764 + t798 * t329 * t330 * t244 / 0.3e1 + t813 * t352 * t353 * t316 / 0.3e1;
1675 t1983 = t805 / 0.3e1;
1676 t1985 = t820 / 0.3e1;
1677 t1986 = t9 * (t1849 + t1909 + t1952 + t1976) - t786 / 0.2e1 - t790 / 0.2e1 + t795 + t801 / 0.3e1 - t1983 - t810 + t816 / 0.3e1 - t1985 - t824 - t754;
1678 t1993 = -t734 + t983 - t831 / 0.2e1 - t834 - t837 - 0.3e1 / 0.2e1 * t841 - 0.7e1 / 0.2e1 * t847 - 0.3e1 / 0.2e1 * t851 + t982 + 0.7e1 / 0.2e1 * t859 - 0.3e1 * t863;
1679 t1997 = -t635 - t639 - t648 + t725 - t729 / 0.3e1 + t737 + t746 - t749 / 0.3e1 + t757 - t764 - t770;
1680 t2003 = -t626 - 0.2e1 * t628 - t630 - 0.3e1 / 0.2e1 * t872 + 0.3e1 / 0.2e1 * t875 + t984 - t1786 - t1789 + 0.3e1 / 0.4e1 * t650 - 0.3e1 / 0.4e1 * t653 - t717;
1681 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1682 v4rho4[1] = t1986 + t1993 + t1997 + t2003;
1683
1684 t2006 = t351 * t818 * t894;
1685 t2013 = t813 * t818 * t906;
1686 t2020 = t351 * t86 * t174 * t353;
1687 t2027 = t798 * t803 * t898;
1688 t2034 = t328 * t803 * t902;
1689 t2037 = t328 * t1867 * t914;
1690 t2051 = 0.1e1 / t85 / t1662;
1691 t2053 = t8 * t8;
1692 t2059 = t7 * t312 * t924;
1693 t2062 = t7 * t1718 * t377;
1694 t2065 = t7 * t234 * t953;
1695 t2068 = t7 * t845 * t957;
1696 t2071 = t7 * t312 * t920;
1697 t2074 = t7 * t727 * t944;
1698 t2076 = -t2006 / 0.9e1 - 0.5e1 / 0.27e2 * t1916 * t352 * t387 * t307 + 0.2e1 / 0.27e2 * t2013 + 0.2e1 / 0.9e1 * t813 * t890 * t814 - 0.2e1 / 0.9e1 * t2020 - 0.5e1 / 0.27e2 * t1851 * t329 * t376 * t229 + 0.2e1 / 0.27e2 * t2027 + t798 * t329 * t382 * t229 / 0.9e1 - t2034 / 0.9e1 + 0.7e1 / 0.9e1 * t2037 + 0.4e1 / 0.9e1 * t798 * t807 * t914 * t80 - 0.4e1 / 0.9e1 * t813 * t792 * t353 * t8 - t7 * t75 * t120 * t952 - t7 * t227 * t2051 * t42 * t2053 / 0.3e1 - t2059 / 0.3e1 + t2062 / 0.27e2 - t2065 / 0.3e1 - 0.7e1 / 0.3e1 * t2068 - t2071 / 0.3e1 - t2074 / 0.9e1;
1699 t2078 = t7 * t234 * t948;
1700 t2104 = t7 * t755 * t393;
1701 t2107 = t7 * t755 * t397;
1702 t2110 = t7 * t735 * t383;
1703 t2113 = t7 * t1709 * t388;
1704 t2127 = -t2078 / 0.3e1 - t7 * t1819 * t138 / 0.3e1 - t7 * t124 * t444 * t316 / 0.2e1 - t7 * t305 * t2051 * t69 * t2053 / 0.3e1 - t7 * t228 * t300 * t376 / 0.6e1 - t7 * t75 * t300 * t382 / 0.2e1 + t7 * t1845 * t143 / 0.3e1 + 0.2e1 / 0.9e1 * t2104 + t2107 / 0.9e1 + t2110 / 0.9e1 + t2113 / 0.27e2 - t7 * t306 * t444 * t307 / 0.6e1 + 0.4e1 / 0.9e1 * t1792 - 0.2e1 / 0.3e1 * t1795 - 0.10e2 / 0.27e2 * t1798 - t1805 / 0.9e1 + 0.4e1 / 0.9e1 * t1808 - 0.41e2 / 0.9e1 * t1814 - 0.10e2 / 0.27e2 * t1817 - t1828 / 0.3e1;
1705 t2145 = t351 * t1888 * t910;
1706 t2155 = t342 * t174;
1707 t2183 = 0.41e2 / 0.9e1 * t1840 - t1843 / 0.3e1 + t813 * t352 * t396 * t307 / 0.9e1 + 0.2e1 / 0.9e1 * t798 * t822 * t898 - t328 * t822 * t902 / 0.3e1 - 0.2e1 / 0.3e1 * t328 * t1874 * t914 - 0.7e1 / 0.9e1 * t2145 + t798 * t329 * t376 * t244 / 0.9e1 - t328 * t329 * t952 * t80 / 0.3e1 + 0.2e1 / 0.3e1 * t351 * t2155 * t793 - t351 * t352 * t919 * t125 / 0.3e1 - t351 * t890 * t788 / 0.3e1 - t351 * t352 * t396 * t316 / 0.6e1 - t328 * t329 * t382 * t244 / 0.6e1 + t813 * t352 * t387 * t316 / 0.9e1 + 0.2e1 / 0.27e2 * t1857 + 0.4e1 / 0.27e2 * t1865 - 0.2e1 / 0.9e1 * t1869 - t1872 / 0.9e1 - t1886 / 0.9e1;
1708 t2189 = t7 * t321 * t444;
1709 t2205 = 0.2e1 / 0.9e1 * t1890 - 0.2e1 / 0.9e1 * t1897 + 0.2e1 / 0.27e2 * t1922 + 0.4e1 / 0.27e2 * t1926 + t2189 / 0.12e2 + t7 * t360 * t174 + t1938 / 0.3e1 - t1941 / 0.3e1 - 0.5e1 / 0.18e2 * t1947 - 0.5e1 / 0.18e2 * t1630 + t1635 + t1639 + t1641 / 0.12e2 + t1672 / 0.9e1 + 0.2e1 / 0.9e1 * t1691 + t1694 / 0.9e1 + t1711 / 0.27e2 + t1720 / 0.27e2 - 0.10e2 / 0.27e2 * t1761 - 0.10e2 / 0.27e2 * t1764;
1710 t2221 = t9 * (t2076 + t2127 + t2183 + t2205) - 0.2e1 / 0.3e1 * t892 - t896 / 0.3e1 + 0.2e1 / 0.9e1 * t900 - t904 / 0.3e1 + 0.2e1 / 0.9e1 * t908 + 0.2e1 / 0.3e1 * t912 - 0.2e1 / 0.3e1 * t916 - t786 / 0.3e1 - t790 / 0.3e1 + 0.2e1 / 0.3e1 * t795 + 0.2e1 / 0.9e1 * t801 - 0.4e1 / 0.9e1 * t805;
1711 t2232 = -0.2e1 / 0.3e1 * t810 + 0.2e1 / 0.9e1 * t816 - 0.4e1 / 0.9e1 * t820 - 0.2e1 / 0.3e1 * t824 - t922 - t926 - t929 / 0.9e1 - t932 / 0.3e1 - t935 / 0.9e1 - 0.2e1 / 0.3e1 * t938 - t941 / 0.3e1 - t946 / 0.3e1 - t950;
1712 t2244 = -t955 + 0.2e1 * t959 + 0.2e1 / 0.3e1 * t827 - t831 / 0.3e1 - 0.4e1 / 0.3e1 * t834 - 0.4e1 / 0.3e1 * t837 - t841 - 0.5e1 / 0.3e1 * t847 - t851 + 0.2e1 / 0.3e1 * t855 + 0.5e1 / 0.3e1 * t859 - 0.2e1 * t863 - t729 / 0.9e1;
1713 t2255 = 0.2e1 / 0.3e1 * t737 - t749 / 0.9e1 + 0.2e1 / 0.3e1 * t757 - t764 / 0.3e1 - 0.2e1 / 0.3e1 * t628 - t630 / 0.3e1 - t967 / 0.4e1 - t872 + t875 + t878 / 0.2e1 - t1786 - t1789 + t650 / 0.2e1 - t653 / 0.4e1;
1714 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1715 v4rho4[2] = t2221 + t2232 + t2244 + t2255;
1716
1717 t2258 = t7 * t130 * t1091;
1718 t2312 = t328 * t803 * t1019;
1719 t2314 = -t2258 / 0.8e1 + t813 * t890 * t906 / 0.3e1 + t351 * t2155 * t910 - t351 * t16 * t444 * t353 / 0.2e1 - t351 * t890 * t894 / 0.2e1 - t351 * t352 * t1005 * t125 / 0.6e1 - t328 * t329 * t987 * t80 / 0.6e1 - t328 * t329 * t137 * t952 / 0.2e1 - 0.5e1 / 0.27e2 * t1916 * t352 * t1023 * t125 - t813 * t792 * t387 * t8 / 0.3e1 - t351 * t352 * t142 * t919 / 0.2e1 + t351 * t792 * t8 * t396 / 0.2e1 - 0.5e1 / 0.27e2 * t1851 * t329 * t1014 * t80 - t328 * t822 * t1019 / 0.2e1 + t798 * t807 * t376 * t8 / 0.3e1 - t2312 / 0.6e1;
1720 t2320 = t351 * t818 * t1028;
1721 t2337 = t7 * t312 * t996;
1722 t2340 = t7 * t312 * t1000;
1723 t2343 = t7 * t312 * t1006;
1724 t2355 = 0.6e1 * t641 + 0.12e2 * t1664;
1725 t2364 = -t328 * t807 * t8 * t382 / 0.2e1 - t2320 / 0.6e1 - t2006 / 0.6e1 + t2013 / 0.9e1 - t2020 / 0.3e1 + t2027 / 0.9e1 - t2034 / 0.6e1 - t2037 / 0.3e1 + t813 * t352 * t1028 * t125 / 0.3e1 + t798 * t329 * t1019 * t80 / 0.3e1 - t2337 / 0.2e1 - t2340 / 0.2e1 - t2343 / 0.6e1 - t7 * t124 * t1091 * t125 / 0.2e1 - t7 * t75 * t120 * t987 / 0.2e1 - t7 * t75 * t42 * t2355 / 0.2e1 + 0.3e1 / 0.2e1 * t7 * t360 * t444 * t8;
1726 t2376 = t7 * t1705 * t1024;
1727 t2379 = t7 * t1736 * t1015;
1728 t2382 = t7 * t234 * t988;
1729 t2385 = t7 * t747 * t992;
1730 t2400 = -0.3e1 / 0.2e1 * t7 * t124 * t174 * t919 + t7 * t124 * t69 * t2355 / 0.2e1 + t2376 / 0.27e2 + t2379 / 0.27e2 - t2382 / 0.6e1 - t2385 / 0.6e1 + t7 * t720 * t120 * t1014 / 0.9e1 - t2059 / 0.2e1 + t2062 / 0.9e1 - t2065 / 0.2e1 + t2068 - t2071 / 0.2e1 - t2074 / 0.6e1 - t2078 / 0.2e1 + 0.2e1 / 0.3e1 * t2104 + t2107 / 0.3e1 + t2110 / 0.3e1;
1731 t2412 = t2113 / 0.9e1 + t1793 - 0.5e1 / 0.9e1 * t1798 + t1809 + t1814 / 0.3e1 - 0.5e1 / 0.9e1 * t1817 - t1840 / 0.3e1 + t2145 / 0.3e1 + t1866 + t1927 + t2189 / 0.4e1 - 0.5e1 / 0.12e2 * t1947 - 0.5e1 / 0.36e2 * t1630 + t1635 + t1639 - 0.5e1 / 0.27e2 * t1761 - 0.5e1 / 0.27e2 * t1764;
1732 t2420 = t9 * (t2314 + t2364 + t2400 + t2412) - t1022 - t1031 - t892 - t896 / 0.2e1 + t900 / 0.3e1 - t904 / 0.2e1 + t908 / 0.3e1 + t912 - t916 - t1983;
1733 t2424 = -t1985 + t1018 + t1027 - t991 - t995 - t999 - t1003 - t1009 - 0.3e1 / 0.2e1 * t922 - 0.3e1 / 0.2e1 * t926 - t929 / 0.3e1;
1734 t2432 = -t932 - t935 / 0.3e1 - 0.2e1 * t938 - t941 - t946 / 0.2e1 - 0.3e1 / 0.2e1 * t950 - 0.3e1 / 0.2e1 * t955 + 0.3e1 * t959 + t827 - t834 - t837;
1735 t2435 = t847 + t855 - t859 + t738 + t758 - t1094 - 0.3e1 / 0.4e1 * t967 + 0.3e1 / 0.4e1 * t878 - t1786 - t1789 + t651;
1736 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1737 v4rho4[3] = t2420 + t2424 + t2432 + t2435;
1738
1739 t2440 = t396 * t396;
1740 t2445 = t376 * t376;
1741 t2450 = t382 * t382;
1742 t2455 = t387 * t387;
1743 t2477 = -0.12e2 * t641 - 0.12e2 * t1664;
1744 t2488 = 0.1e1 / t53 / t534;
1745 t2512 = 0.1e1 / t1047 / t63;
1746 t2513 = t414 * t414;
1747 t2523 = t439 * t439;
1748 t2546 = t1069 * t402;
1749 t2552 = t429 * t429;
1750 t2559 = 0.1e1 / t164 / t1080 / t163;
1751 t2568 = 0.41888e5 / 0.729e3 * t21 * t50 * t2488 * t64 + 0.9856e4 / 0.243e3 * t99 * t1065 * t170 + 0.704e3 / 0.27e2 * t99 * t425 * t415 - 0.352e3 / 0.27e2 * t99 * t425 * t440 + 0.128e3 / 0.9e1 * t99 * t162 * t1050 - 0.128e3 / 0.9e1 * t99 * t162 * t1055 + 0.64e2 / 0.27e2 * t99 * t162 * t1087 + 0.16e2 / 0.3e1 * t99 * t154 * t2512 * t2513 - 0.8e1 * t99 * t154 * t1048 * t414 * t439 + 0.4e1 / 0.3e1 * t99 * t154 * t413 * t2523 + 0.16e2 / 0.9e1 * t99 * t154 * t1054 * t1086 - 0.2e1 / 0.9e1 * t99 * t154 * t156 * (0.3640e4 / 0.81e2 * t28 * t56 * t601 * t60 + 0.5768e4 / 0.27e2 * t28 * sigma[2] * t2488 * t165 - 0.18608e5 / 0.81e2 * t28 * t429 / t52 / t1431 * t435 + 0.4480e4 / 0.27e2 * t28 * t1076 / t2546 * t1082 - 0.1280e4 / 0.27e2 * t28 * t2552 / t53 / t1069 / t534 * t2559);
1752 t2593 = -t7 * t306 * t444 * t387 - t7 * t306 * t69 * t2440 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1730 * t42 * t2445 - t7 * t228 * t42 * t2450 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1699 * t69 * t2455 + 0.4e1 / 0.9e1 * t7 * t741 * t174 * t1023 - 0.2e1 * t7 * t124 * t1091 * t142 - 0.3e1 * t7 * t124 * t444 * t396 - 0.2e1 * t7 * t124 * t174 * t1005 - t7 * t124 * t69 * t2477 / 0.2e1 + t7 * t75 * t42 * t2477 / 0.2e1 - 0.3e1 / 0.8e1 * t7 * t49 * t2568 - t2258 / 0.2e1 + 0.2e1 / 0.3e1 * t798 * t329 * t376 * t382 + 0.2e1 / 0.3e1 * t813 * t352 * t387 * t396 - 0.2e1 / 0.3e1 * t351 * t352 * t1005 * t142 - 0.2e1 * t351 * t890 * t1028 - 0.2e1 / 0.3e1 * t328 * t329 * t987 * t137 - 0.2e1 / 0.3e1 * t2312;
1753 t2611 = -0.2e1 / 0.3e1 * t2320 - 0.2e1 * t2337 - 0.2e1 * t2340 - 0.2e1 / 0.3e1 * t2343 + 0.4e1 / 0.27e2 * t2376 + 0.4e1 / 0.27e2 * t2379 - 0.2e1 / 0.3e1 * t2382 - 0.2e1 / 0.3e1 * t2385 + 0.2e1 / 0.9e1 * t2062 + 0.4e1 / 0.3e1 * t2104 + 0.2e1 / 0.3e1 * t2107 + 0.2e1 / 0.3e1 * t2110 + 0.2e1 / 0.9e1 * t2113 - 0.20e2 / 0.27e2 * t1798 - 0.20e2 / 0.27e2 * t1817 + t2189 / 0.2e1 - 0.5e1 / 0.9e1 * t1947 + t1635 + t1639;
1754 t2621 = t9 * (t2593 + t2611) + 0.4e1 / 0.3e1 * t855 + 0.4e1 / 0.3e1 * t827 - t1786 - t1789 + t878 - 0.3e1 / 0.2e1 * t967 - 0.2e1 * t990 - 0.2e1 * t994 - 0.6e1 * t998 - 0.6e1 * t1002;
1755 t2633 = -0.2e1 * t1008 - 0.2e1 / 0.3e1 * t929 - 0.2e1 * t932 - 0.2e1 / 0.3e1 * t935 - 0.4e1 * t938 - 0.2e1 * t941 + 0.4e1 / 0.9e1 * t1017 - 0.2e1 * t1021 + 0.4e1 / 0.9e1 * t1026 - 0.2e1 * t1030 - 0.3e1 / 0.2e1 * t1093;
1756 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1757 v4rho4[4] = t2621 + t2633;
1758
1759 t2661 = t19 * t98 * sigma[0];
1760 t2662 = t27 * t671;
1761 t2736 = -0.2e1 / 0.9e1 * t21 * t184 * t710 - 0.2e1 / 0.9e1 * t99 * t100 * t102 * (-0.140e3 / 0.27e2 * t28 * t187 * t683 * t33 - 0.364e3 / 0.9e1 * t28 * t656 * t111 + 0.1660e4 / 0.27e2 * t28 * t694 * t291 * sigma[0] - 0.488e3 / 0.9e1 * t28 * t701 * t705 * t285 + 0.160e3 / 0.9e1 * t28 / t25 / t692 / t567 * t1613 * t699) + 0.4e1 / 0.3e1 * t99 * t1324 * t115 * t295 - 0.176e3 / 0.27e2 * t99 * t281 * t488 + 0.16e2 / 0.9e1 * t99 * t108 * t1171 - 0.32e2 / 0.9e1 * t21 * t458 * t271 - 0.4e1 / 0.3e1 * t21 * t184 * t673 + 0.2464e4 / 0.243e3 * t99 * t688 * t195 + 0.352e3 / 0.27e2 * t99 * t281 * t469 - 0.64e2 / 0.9e1 * t99 * t108 * t1145 - 0.32e2 / 0.9e1 * t99 * t108 * t1149;
1762 t2742 = t7 * t17 * t1175;
1763 t2745 = t7 * t87 * t492;
1764 t2764 = t7 * t75 * t1110;
1765 t2767 = t7 * t75 * t1113;
1766 t2771 = t7 * t251 * t199 / 0.9e1;
1767 t2780 = t7 * t228 * t1106;
1768 t2783 = t7 * t234 * t454;
1769 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1770 v4rho3sigma[0] = -0.3e1 / 0.8e1 * t182 * t183 * (-0.4e1 * t99 * t100 * t671 * t487 * t270 + 0.4e1 / 0.3e1 * t99 * t100 * t269 * t1170 * t115 + 0.4e1 / 0.3e1 * t99 * t100 * t1144 * t295 + 0.4e1 / 0.9e1 * t99 * t100 * t468 * t709 + 0.32e2 / 0.3e1 * t99 * t108 * t1140 + 0.16e2 / 0.3e1 * t99 * t100 * t1566 * t194 * t672 - 0.4e1 * t2661 * t2662 * t1325 * t295 - 0.2464e4 / 0.243e3 * t21 * t5 * t656 * t37 - 0.176e3 / 0.27e2 * t21 * t1117 * t116 + 0.16e2 / 0.9e1 * t21 * t458 * t296 + t2736) - 0.3e1 / 0.2e1 * t2742 - t2745 / 0.2e1 - 0.3e1 / 0.2e1 * t182 * t453 * t492 * t244 - t182 * t453 * t199 * t644 / 0.2e1 - t182 * t327 * t492 * t229 / 0.2e1 - 0.3e1 / 0.2e1 * t182 * t453 * t1175 * t80 - 0.4e1 * t2764 - 0.2e1 * t2767 + t2771 + t182 * t797 * t199 * t721 / 0.9e1 - t1200 * t1201 * t731 / 0.2e1 - 0.2e1 / 0.3e1 * t2780 - 0.2e1 / 0.3e1 * t2783;
1771
1772 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1773 v4rho3sigma[1] = 0.0e0;
1774
1775 t2786 = t7 * t306 * t1185;
1776 t2789 = t7 * t312 * t499;
1777 t2792 = t7 * t124 * t1189;
1778 t2796 = t7 * t321 * t219 / 0.9e1;
1779 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1780 v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2786 - 0.2e1 / 0.3e1 * t2789 - 0.2e1 * t2792 + t2796 + t182 * t812 * t219 * t742 / 0.9e1 - t1200 * t1218 * t751 / 0.2e1 - t182 * t498 * t219 * t766 / 0.2e1;
1781
1782 t2819 = t86 * t1 * t3;
1783 t2820 = t2819 * t504;
1784 t2822 = t1195 * t1199;
1785 t2824 = 0.4e1 / 0.9e1 * t2822 * t1202;
1786 t2826 = t342 * t1 * t3;
1787 t2827 = t2826 * t1212;
1788 t2829 = t719 * t199;
1789 t2833 = t1209 * t1199;
1790 t2843 = 0.4e1 / 0.3e1 * t1196 * t1206;
1791 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1792 v4rho3sigma[3] = -t1210 * t453 * t492 * t8 - t1200 * t227 * t492 * t330 / 0.3e1 - t1200 * t1201 * t784 / 0.6e1 - 0.2e1 / 0.9e1 * t2820 - t2824 + t2827 / 0.6e1 + t1200 * t2829 * t799 / 0.9e1 - t2833 * t1201 * t808 / 0.3e1 - t2742 / 0.2e1 - t2745 / 0.3e1 - t1210 * t453 * t199 / 0.2e1 - t2843 - t182 * t453 * t1175 * t137 / 0.2e1 - 0.4e1 / 0.3e1 * t2764 - 0.2e1 / 0.3e1 * t2767 + t2771 - 0.2e1 / 0.9e1 * t2780 - 0.4e1 / 0.9e1 * t2783;
1793
1794 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1795 v4rho3sigma[4] = 0.0e0;
1796
1797 t2855 = t2819 * t509;
1798 t2858 = 0.4e1 / 0.9e1 * t2822 * t1219;
1799 t2859 = t2826 * t1223;
1800 t2861 = t740 * t219;
1801 t2874 = t2819 * t547;
1802 t2877 = 0.4e1 / 0.3e1 * t1196 * t1229;
1803 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1804 v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2786 - 0.4e1 / 0.9e1 * t2789 - 0.2e1 / 0.3e1 * t2792 + t2796 - 0.2e1 / 0.9e1 * t2855 - t2858 - t2859 / 0.6e1 + t1200 * t2861 * t814 / 0.9e1 + t2833 * t1218 * t793 / 0.3e1 - t1200 * t1218 * t788 / 0.6e1 + t1210 * t498 * t219 / 0.2e1 - t2874 / 0.6e1 - t2877 - t182 * t350 * t546 * t307 / 0.6e1 - t182 * t498 * t546 * t316 / 0.2e1;
1805
1806 t2890 = t1196 * t1234;
1807 t2902 = t1196 * t1238;
1808 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1809 v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2820 - t2824 - t2843 - 0.4e1 / 0.3e1 * t2827 - 0.2e1 / 0.9e1 * t2783 + t2771 - t2745 / 0.6e1 - 0.2e1 / 0.9e1 * t2890 + t1200 * t2829 * t898 / 0.9e1 - t182 * t327 * t492 * t376 / 0.6e1 - t2833 * t1201 * t914 / 0.3e1 - 0.2e1 / 0.3e1 * t2902 - t1200 * t1201 * t902 / 0.6e1 - t182 * t453 * t492 * t382 / 0.2e1 - t182 * t453 * t199 * t952 / 0.2e1;
1810
1811 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1812 v4rho3sigma[7] = 0.0e0;
1813
1814 t2919 = t1196 * t1243;
1815 t2927 = t1196 * t1247;
1816 t2936 = t1196 * t1250;
1817 t2945 = t1196 * t1312;
1818 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1819 v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2855 - t2858 + 0.4e1 / 0.3e1 * t2859 - 0.2e1 / 0.9e1 * t2789 + t2796 - t2874 / 0.3e1 - t2877 - 0.2e1 / 0.9e1 * t2919 + t1200 * t2861 * t906 / 0.9e1 + t2833 * t1218 * t910 / 0.3e1 - 0.4e1 / 0.3e1 * t2927 - t1200 * t305 * t546 * t353 / 0.3e1 + t1210 * t498 * t546 * t8 - 0.2e1 / 0.3e1 * t2936 - t1200 * t1218 * t894 / 0.6e1 - t182 * t498 * t219 * t919 / 0.2e1 - t2945 / 0.2e1 - t182 * t498 * t1311 * t125 / 0.2e1;
1820
1821 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1822 v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2820 - 0.2e1 / 0.3e1 * t2890 - 0.2e1 * t2902 + t2771 + t182 * t797 * t199 * t1014 / 0.9e1 - t1200 * t1201 * t1019 / 0.2e1 - t182 * t453 * t199 * t987 / 0.2e1;
1823
1824 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1825 v4rho3sigma[10] = 0.0e0;
1826
1827 t3082 = t19 * t98 * sigma[2];
1828 t3083 = t55 * t1048;
1829 t3091 = 0.352e3 / 0.27e2 * t99 * t425 * t523 - 0.64e2 / 0.9e1 * t99 * t162 * t1281 - 0.32e2 / 0.9e1 * t99 * t162 * t1285 - 0.4e1 * t99 * t154 * t1048 * t541 * t414 + 0.4e1 / 0.3e1 * t99 * t154 * t413 * t1306 * t169 + 0.4e1 / 0.3e1 * t99 * t154 * t1280 * t439 + 0.4e1 / 0.9e1 * t99 * t154 * t522 * t1086 + 0.2464e4 / 0.243e3 * t99 * t1065 * t215 + 0.16e2 / 0.3e1 * t99 * t154 * t2512 * t214 * t1049 - 0.4e1 * t3082 * t3083 * t1394 * t439 + 0.32e2 / 0.3e1 * t99 * t162 * t1276;
1830 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1831 v4rho3sigma[11] = -t1200 * t1218 * t1028 / 0.2e1 + t182 * t812 * t219 * t1023 / 0.9e1 - 0.3e1 / 0.2e1 * t182 * t498 * t546 * t396 - t182 * t498 * t219 * t1005 / 0.2e1 - t182 * t350 * t546 * t387 / 0.2e1 - 0.3e1 / 0.2e1 * t182 * t498 * t1311 * t142 - 0.3e1 / 0.2e1 * t2945 - 0.3e1 / 0.8e1 * t182 * t203 * (0.16e2 / 0.9e1 * t99 * t162 * t1307 - 0.2e1 / 0.9e1 * t99 * t154 * t156 * (-0.140e3 / 0.27e2 * t28 * t207 * t1060 * t60 - 0.364e3 / 0.9e1 * t28 * t1033 * t165 + 0.1660e4 / 0.27e2 * t28 * t1071 * t435 * sigma[2] - 0.488e3 / 0.9e1 * t28 * t1078 * t1082 * t429 + 0.160e3 / 0.9e1 * t28 / t53 / t1069 / t599 * t2559 * t1076) - 0.4e1 / 0.3e1 * t21 * t204 * t1050 + 0.4e1 / 0.3e1 * t99 * t1393 * t169 * t439 - 0.176e3 / 0.27e2 * t99 * t425 * t542 - 0.32e2 / 0.9e1 * t21 * t512 * t415 - 0.2e1 / 0.9e1 * t21 * t204 * t1087 - 0.2464e4 / 0.243e3 * t21 * t5 * t1033 * t64 - 0.176e3 / 0.27e2 * t21 * t1253 * t170 + 0.16e2 / 0.9e1 * t21 * t512 * t440 + t3091) - 0.4e1 * t2927 - 0.2e1 * t2936 - 0.2e1 / 0.3e1 * t2919 - 0.2e1 / 0.3e1 * t2855 - t2874 / 0.2e1 + t2796;
1832
1833 t3102 = t2819 * t579 / 0.6e1;
1834 t3103 = t1196 * t1318;
1835 t3105 = t1196 * t1373;
1836 t3165 = t487 * t487;
1837 t3177 = 0.32e2 / 0.27e2 * t99 * t108 * t1368 - 0.2e1 / 0.9e1 * t99 * t100 * t102 * (-0.7e1 / 0.9e1 * t28 * t559 * t276 * t33 + 0.7e1 / 0.9e1 * t28 * t563 * t260 * t111 - 0.11e2 * t28 * t1160 + 0.47e2 / 0.3e1 * t28 * t1166 * sigma[0] - 0.20e2 / 0.3e1 * t28 / t25 / t1600 * t1613 * t285) - 0.352e3 / 0.81e2 * t21 * t1117 * t195 - 0.128e3 / 0.27e2 * t99 * t93 * t269 * t1325 + 0.16e2 / 0.9e1 * t99 * t1324 * t487 * t115 + 0.8e1 / 0.9e1 * t99 * t1324 * t194 * t295 + 0.64e2 / 0.27e2 * t21 * t458 * t488 - 0.4e1 / 0.9e1 * t21 * t184 * t1171 + 0.8e1 / 0.9e1 * t99 * t100 * t269 * t3165 + 0.8e1 / 0.9e1 * t99 * t100 * t468 * t1170 - 0.176e3 / 0.81e2 * t99 * t281 * t574;
1838 t3221 = -0.64e2 / 0.27e2 * t99 * t108 * t1348 + 0.8e1 / 0.9e1 * t99 * t100 * t269 * t1367 * t115 + 0.4e1 / 0.9e1 * t99 * t100 * t1347 * t295 - 0.8e1 / 0.3e1 * t99 * t2662 * t194 * t270 + 0.352e3 / 0.81e2 * t99 * t281 * t554 + 0.64e2 / 0.9e1 * t99 * t108 * t1336 - 0.128e3 / 0.27e2 * t99 * t108 * t1340 - 0.4e1 / 0.3e1 * t99 * t100 * t1335 * t295 - 0.4e1 / 0.3e1 * t99 * t100 * t671 * t573 * t270 + 0.16e2 / 0.3e1 * t99 * t100 * t1566 * t553 * t270 - 0.16e2 / 0.3e1 * t2661 * t2662 * t1325 * t487;
1839 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1840 v4rho2sigma2[0] = -t3102 - 0.4e1 / 0.3e1 * t3103 - t3105 - t182 * t327 * t578 * t229 / 0.6e1 - t182 * t453 * t1372 * t80 - t182 * t453 * t578 * t244 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t183 * (t3177 + t3221);
1841
1842 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1843 v4rho2sigma2[1] = 0.0e0;
1844
1845 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1846 v4rho2sigma2[2] = 0.0e0;
1847
1848 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1849 v4rho2sigma2[3] = 0.0e0;
1850
1851 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1852 v4rho2sigma2[4] = 0.0e0;
1853
1854 t3227 = t2819 * t611 / 0.6e1;
1855 t3228 = t1196 * t1378;
1856 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1857 v4rho2sigma2[5] = -t3227 - 0.4e1 / 0.3e1 * t3228 - t182 * t350 * t610 * t307 / 0.6e1 - t182 * t498 * t610 * t316 / 0.2e1;
1858
1859 t3240 = t1196 * t1382;
1860 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1861 v4rho2sigma2[6] = -t3102 - 0.2e1 / 0.3e1 * t3103 - t3105 / 0.2e1 - 0.2e1 / 0.3e1 * t3240 - t1200 * t227 * t578 * t330 / 0.6e1 - t182 * t453 * t1372 * t137 / 0.2e1 - t1210 * t453 * t578 * t8 / 0.2e1;
1862
1863 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1864 v4rho2sigma2[7] = 0.0e0;
1865
1866 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1867 v4rho2sigma2[8] = 0.0e0;
1868
1869 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1870 v4rho2sigma2[9] = 0.0e0;
1871
1872 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1873 v4rho2sigma2[10] = 0.0e0;
1874
1875 t3255 = t1196 * t1387;
1876 t3265 = t1196 * t1442;
1877 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1878 v4rho2sigma2[11] = -t3227 - 0.2e1 / 0.3e1 * t3228 - 0.2e1 / 0.3e1 * t3255 - t1200 * t305 * t610 * t353 / 0.6e1 + t1210 * t498 * t610 * t8 / 0.2e1 - t3265 / 0.2e1 - t182 * t498 * t1441 * t125 / 0.2e1;
1879
1880 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1881 v4rho2sigma2[12] = -t3102 - 0.4e1 / 0.3e1 * t3240 - t182 * t327 * t578 * t376 / 0.6e1 - t182 * t453 * t578 * t382 / 0.2e1;
1882
1883 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1884 v4rho2sigma2[13] = 0.0e0;
1885
1886 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1887 v4rho2sigma2[14] = 0.0e0;
1888
1889 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1890 v4rho2sigma2[15] = 0.0e0;
1891
1892 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1893 v4rho2sigma2[16] = 0.0e0;
1894
1895 t3349 = 0.32e2 / 0.27e2 * t99 * t162 * t1437 - 0.2e1 / 0.9e1 * t99 * t154 * t156 * (-0.7e1 / 0.9e1 * t28 * t591 * t420 * t60 + 0.7e1 / 0.9e1 * t28 * t595 * t404 * t165 - 0.11e2 * t28 * t1296 + 0.47e2 / 0.3e1 * t28 * t1302 * sigma[2] - 0.20e2 / 0.3e1 * t28 / t53 / t2546 * t2559 * t429) + 0.16e2 / 0.9e1 * t99 * t1393 * t541 * t169 + 0.8e1 / 0.9e1 * t99 * t1393 * t214 * t439 - 0.352e3 / 0.81e2 * t21 * t1253 * t215 - 0.128e3 / 0.27e2 * t99 * t149 * t413 * t1394 - 0.4e1 / 0.9e1 * t21 * t204 * t1307 + 0.64e2 / 0.27e2 * t21 * t512 * t542 + 0.352e3 / 0.81e2 * t99 * t425 * t586 + 0.64e2 / 0.9e1 * t99 * t162 * t1405 - 0.128e3 / 0.27e2 * t99 * t162 * t1409;
1896 t3354 = t541 * t541;
1897 t3396 = -0.4e1 / 0.3e1 * t99 * t154 * t1404 * t439 + 0.8e1 / 0.9e1 * t99 * t154 * t413 * t3354 + 0.8e1 / 0.9e1 * t99 * t154 * t522 * t1306 - 0.176e3 / 0.81e2 * t99 * t425 * t606 - 0.64e2 / 0.27e2 * t99 * t162 * t1417 + 0.8e1 / 0.9e1 * t99 * t154 * t413 * t1436 * t169 + 0.4e1 / 0.9e1 * t99 * t154 * t1416 * t439 - 0.8e1 / 0.3e1 * t99 * t3083 * t214 * t414 - 0.4e1 / 0.3e1 * t99 * t154 * t1048 * t605 * t414 + 0.16e2 / 0.3e1 * t99 * t154 * t2512 * t585 * t414 - 0.16e2 / 0.3e1 * t3082 * t3083 * t1394 * t541;
1898 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1899 v4rho2sigma2[17] = -t3227 - 0.4e1 / 0.3e1 * t3255 - t3265 - t182 * t350 * t610 * t387 / 0.6e1 - t182 * t498 * t1441 * t142 - t182 * t498 * t610 * t396 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t203 * (t3349 + t3396);
1900
1901 t3401 = t1196 * t1485;
1902 t3443 = t194 * t573;
1903 t3489 = -0.32e2 / 0.9e1 * t21 * t458 * t554 - 0.4e1 * t99 * t2662 * t553 * t115 + 0.8e1 / 0.3e1 * t99 * t1324 * t194 * t487 + 0.16e2 / 0.9e1 * t21 * t458 * t574 + 0.4e1 / 0.3e1 * t99 * t1324 * t573 * t115 - 0.2e1 / 0.3e1 * t21 * t184 * t1368 + 0.32e2 / 0.9e1 * t99 * t108 * t1452 + 0.16e2 / 0.3e1 * t99 * t100 * t1566 * t1451 * t115 - 0.4e1 * t99 * t100 * t1335 * t487 - 0.32e2 / 0.9e1 * t99 * t108 * t1456 - 0.4e1 * t2661 * t2662 * t3443 * t115 + 0.4e1 / 0.3e1 * t99 * t100 * t1144 * t573 + 0.4e1 / 0.3e1 * t99 * t100 * t468 * t1367 + 0.16e2 / 0.27e2 * t99 * t108 * t1480 + 0.4e1 / 0.9e1 * t99 * t100 * t269 * t1479 * t115 - 0.2e1 / 0.9e1 * t99 * t100 * t102 * (-t28 * t1461 * t104 * t33 / 0.2e1 + t28 * t1466 * t93 * t111 / 0.2e1 + t28 * t563 * t482 * t291 / 0.6e1 - 0.7e1 / 0.2e1 * t28 * t1364 + 0.5e1 / 0.2e1 * t28 / t25 / t700 * t1613 * sigma[0]);
1904 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1905 v4rhosigma3[0] = -t3401 / 0.2e1 - t182 * t453 * t1484 * t80 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t183 * t3489;
1906
1907 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1908 v4rhosigma3[1] = 0.0e0;
1909
1910 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1911 v4rhosigma3[2] = 0.0e0;
1912
1913 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1914 v4rhosigma3[3] = 0.0e0;
1915
1916 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1917 v4rhosigma3[4] = 0.0e0;
1918
1919 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1920 v4rhosigma3[5] = 0.0e0;
1921
1922 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1923 v4rhosigma3[6] = 0.0e0;
1924
1925 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1926 v4rhosigma3[7] = 0.0e0;
1927
1928 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1929 v4rhosigma3[8] = 0.0e0;
1930
1931 t3493 = t1196 * t1528;
1932 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1933 v4rhosigma3[9] = -t182 * t498 * t1527 * t125 / 0.2e1 - t3493 / 0.2e1;
1934
1935 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1936 v4rhosigma3[10] = -t182 * t453 * t1484 * t137 / 0.2e1 - t3401 / 0.2e1;
1937
1938 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1939 v4rhosigma3[11] = 0.0e0;
1940
1941 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1942 v4rhosigma3[12] = 0.0e0;
1943
1944 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1945 v4rhosigma3[13] = 0.0e0;
1946
1947 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1948 v4rhosigma3[14] = 0.0e0;
1949
1950 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1951 v4rhosigma3[15] = 0.0e0;
1952
1953 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1954 v4rhosigma3[16] = 0.0e0;
1955
1956 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1957 v4rhosigma3[17] = 0.0e0;
1958
1959 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1960 v4rhosigma3[18] = 0.0e0;
1961
1962 t3543 = t214 * t605;
1963 t3589 = -0.32e2 / 0.9e1 * t21 * t512 * t586 - 0.4e1 * t99 * t3083 * t585 * t169 + 0.8e1 / 0.3e1 * t99 * t1393 * t214 * t541 + 0.16e2 / 0.9e1 * t21 * t512 * t606 + 0.4e1 / 0.3e1 * t99 * t1393 * t605 * t169 - 0.2e1 / 0.3e1 * t21 * t204 * t1437 + 0.32e2 / 0.9e1 * t99 * t162 * t1495 + 0.16e2 / 0.3e1 * t99 * t154 * t2512 * t1494 * t169 - 0.4e1 * t99 * t154 * t1404 * t541 - 0.32e2 / 0.9e1 * t99 * t162 * t1499 - 0.4e1 * t3082 * t3083 * t3543 * t169 + 0.4e1 / 0.3e1 * t99 * t154 * t1280 * t605 + 0.4e1 / 0.3e1 * t99 * t154 * t522 * t1436 + 0.16e2 / 0.27e2 * t99 * t162 * t1523 + 0.4e1 / 0.9e1 * t99 * t154 * t413 * t1522 * t169 - 0.2e1 / 0.9e1 * t99 * t154 * t156 * (-t28 * t1504 * t158 * t60 / 0.2e1 + t28 * t1509 * t149 * t165 / 0.2e1 + t28 * t595 * t536 * t435 / 0.6e1 - 0.7e1 / 0.2e1 * t28 * t1433 + 0.5e1 / 0.2e1 * t28 / t53 / t1077 * t2559 * sigma[2]);
1964 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1965 v4rhosigma3[19] = -t3493 / 0.2e1 - t182 * t498 * t1527 * t142 / 0.2e1 - 0.3e1 / 0.8e1 * t182 * t203 * t3589;
1966
1967 t3602 = t553 * t553;
1968 t3611 = t573 * t573;
1969 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1970 v4sigma4[0] = -0.3e1 / 0.8e1 * t182 * t183 * (-0.16e2 / 0.3e1 * t21 * t184 * t1452 + 0.16e2 / 0.3e1 * t99 * t1324 * t3443 - 0.8e1 / 0.9e1 * t21 * t184 * t1480 + 0.16e2 / 0.3e1 * t99 * t100 * t1566 * t3602 - 0.8e1 * t99 * t100 * t1335 * t573 + 0.4e1 / 0.3e1 * t99 * t100 * t269 * t3611 + 0.16e2 / 0.9e1 * t99 * t100 * t468 * t1479 - 0.2e1 / 0.9e1 * t99 * t100 * t102 * (-0.15e2 / 0.16e2 * t28 / t29 / t699 * t31 * t33 + 0.15e2 / 0.16e2 * t28 / t699 * t27 * t111 + 0.5e1 / 0.16e2 * t28 * t1466 * t569 * t291 + 0.3e1 / 0.16e2 * t28 * t563 * t1475 * t705 - 0.15e2 / 0.16e2 * t28 / t25 / t1164 * t1613));
1971
1972 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1973 v4sigma4[1] = 0.0e0;
1974
1975 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1976 v4sigma4[2] = 0.0e0;
1977
1978 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1979 v4sigma4[3] = 0.0e0;
1980
1981 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1982 v4sigma4[4] = 0.0e0;
1983
1984 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1985 v4sigma4[5] = 0.0e0;
1986
1987 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1988 v4sigma4[6] = 0.0e0;
1989
1990 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1991 v4sigma4[7] = 0.0e0;
1992
1993 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1994 v4sigma4[8] = 0.0e0;
1995
1996 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1997 v4sigma4[9] = 0.0e0;
1998
1999 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2000 v4sigma4[10] = 0.0e0;
2001
2002 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2003 v4sigma4[11] = 0.0e0;
2004
2005 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2006 v4sigma4[12] = 0.0e0;
2007
2008 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2009 v4sigma4[13] = 0.0e0;
2010
2011 t3662 = t585 * t585;
2012 t3671 = t605 * t605;
2013 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2014 v4sigma4[14] = -0.3e1 / 0.8e1 * t182 * t203 * (-0.16e2 / 0.3e1 * t21 * t204 * t1495 + 0.16e2 / 0.3e1 * t99 * t1393 * t3543 - 0.8e1 / 0.9e1 * t21 * t204 * t1523 + 0.16e2 / 0.3e1 * t99 * t154 * t2512 * t3662 - 0.8e1 * t99 * t154 * t1404 * t605 + 0.4e1 / 0.3e1 * t99 * t154 * t413 * t3671 + 0.16e2 / 0.9e1 * t99 * t154 * t522 * t1522 - 0.2e1 / 0.9e1 * t99 * t154 * t156 * (-0.15e2 / 0.16e2 * t28 / t56 / t1076 * t58 * t60 + 0.15e2 / 0.16e2 * t28 / t1076 * t55 * t165 + 0.5e1 / 0.16e2 * t28 * t1509 * t601 * t435 + 0.3e1 / 0.16e2 * t28 * t595 * t1518 * t1082 - 0.15e2 / 0.16e2 * t28 / t53 / t1300 * t2559));
2015
2016 #ifndef XC_DONT_COMPILE_MXC
2017
2018 if(order < 5) return;
2019
2020
2021 #endif
2022
2023 #endif
2024
2025 #endif
2026
2027 #endif
2028
2029 #endif
2030
2031
2032 }
2033
2034