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_2d_b86_mgc.mpl
11 Type of functional: gga_exc
12 */
13
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16
17
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t1, t2, t3, t4, t5, t6, t7, t8;
24 double t9, t11, t12, t13, t14, t15, t18, t20;
25
26 #ifndef XC_DONT_COMPILE_VXC
27 double t23, t24, t25, t26, t30, t31, t32, t35;
28 double t38, t44, t45, t46, t50;
29
30 #ifndef XC_DONT_COMPILE_FXC
31 double t53, t59, t60, t64, t65, t69, t71, t73;
32 double t75, t78, t82, t87, t91, t92, t96, t102;
33 double t107;
34
35 #ifndef XC_DONT_COMPILE_KXC
36 double t110, t126, t130, t132, t136, t139, t143, t149;
37 double t152, t156, t161, t167, t171, t176, t185;
38
39 #ifndef XC_DONT_COMPILE_LXC
40 double t210, t211, t219, t221;
41 #endif
42
43 #endif
44
45 #endif
46
47 #endif
48
49 #endif
50
51
52
53 t1 = sqrt(M_PI);
54 t2 = 0.1e1 / t1;
55 t3 = M_SQRT2;
56 t4 = t2 * t3;
57 t5 = sqrt(rho[0]);
58 t6 = rho[0] * rho[0];
59 t7 = t6 * rho[0];
60 t8 = 0.1e1 / t7;
61 t9 = sigma[0] * t8;
62 t11 = 0.1e1 + 0.16646e-1 * t9;
63 t12 = POW_1_4(t11);
64 t13 = t12 * t12;
65 t14 = t13 * t12;
66 t15 = 0.1e1 / t14;
67 t18 = 0.1e1 + 0.44094220675901974970e-2 * t9 * t15;
68 t20 = t4 * t5 * t18;
69 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
70 zk[0] = -0.4e1 / 0.3e1 * t20;
71
72 #ifndef XC_DONT_COMPILE_VXC
73
74 if(order < 1) return;
75
76
77 t23 = t5 * rho[0];
78 t24 = t23 * t2;
79 t25 = t6 * t6;
80 t26 = 0.1e1 / t25;
81 t30 = sigma[0] * sigma[0];
82 t31 = t25 * t7;
83 t32 = 0.1e1 / t31;
84 t35 = 0.1e1 / t14 / t11;
85 t38 = -0.13228266202770592491e-1 * sigma[0] * t26 * t15 + 0.16514828940848946195e-3 * t30 * t32 * t35;
86 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
87 vrho[0] = -0.2e1 * t20 - 0.4e1 / 0.3e1 * t24 * t3 * t38;
88
89 t44 = t25 * t6;
90 t45 = 0.1e1 / t44;
91 t46 = sigma[0] * t45;
92 t50 = t3 * (0.44094220675901974970e-2 * t8 * t15 - 0.55049429802829820651e-4 * t46 * t35);
93 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
94 vsigma[0] = -0.4e1 / 0.3e1 * t24 * t50;
95
96 #ifndef XC_DONT_COMPILE_FXC
97
98 if(order < 2) return;
99
100
101 t53 = 0.1e1 / t5;
102 t59 = t25 * rho[0];
103 t60 = 0.1e1 / t59;
104 t64 = t25 * t25;
105 t65 = 0.1e1 / t64;
106 t69 = t30 * sigma[0];
107 t71 = 0.1e1 / t64 / t7;
108 t73 = t11 * t11;
109 t75 = 0.1e1 / t14 / t73;
110 t78 = 0.52913064811082369964e-1 * sigma[0] * t60 * t15 - 0.16514828940848946195e-2 * t30 * t65 * t35 + 0.14432556733842006814e-4 * t69 * t71 * t75;
111 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
112 v2rho2[0] = -t4 * t53 * t18 - 0.4e1 * t4 * t5 * t38 - 0.4e1 / 0.3e1 * t24 * t3 * t78;
113
114 t82 = t5 * t2;
115 t87 = t32 * t35;
116 t91 = 0.1e1 / t64 / t6;
117 t92 = t30 * t91;
118 t96 = t3 * (-0.13228266202770592491e-1 * t26 * t15 + 0.49544486822546838586e-3 * t87 * sigma[0] - 0.48108522446140022714e-5 * t92 * t75);
119 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
120 v2rhosigma[0] = -0.2e1 * t82 * t50 - 0.4e1 / 0.3e1 * t24 * t96;
121
122 t102 = 0.1e1 / t64 / rho[0];
123 t107 = t3 * (-0.11009885960565964130e-3 * t45 * t35 + 0.16036174148713340905e-5 * sigma[0] * t102 * t75);
124 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
125 v2sigma2[0] = -0.4e1 / 0.3e1 * t24 * t107;
126
127 #ifndef XC_DONT_COMPILE_KXC
128
129 if(order < 3) return;
130
131
132 t110 = 0.1e1 / t23;
133 t126 = 0.1e1 / t64 / t25;
134 t130 = t30 * t30;
135 t132 = 0.1e1 / t64 / t31;
136 t136 = 0.1e1 / t14 / t73 / t11;
137 t139 = -0.26456532405541184982e0 * t46 * t15 + 0.15193642625581030499e-1 * t30 * t102 * t35 - 0.30308369141068214309e-3 * t69 * t126 * t75 + 0.19820157999801558748e-5 * t130 * t132 * t136;
138 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
139 v3rho3[0] = t4 * t110 * t18 / 0.2e1 - 0.3e1 * t4 * t53 * t38 - 0.6e1 * t4 * t5 * t78 - 0.4e1 / 0.3e1 * t24 * t3 * t139;
140
141 t143 = t53 * t2;
142 t149 = t65 * t35;
143 t152 = t71 * t75;
144 t156 = 0.1e1 / t64 / t44;
145 t161 = t3 * (0.52913064811082369964e-1 * t60 * t15 - 0.39635589458037470869e-2 * t149 * sigma[0] + 0.91406192647666043157e-4 * t152 * t30 - 0.66067193332671862493e-6 * t69 * t156 * t136);
146 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
147 v3rho2sigma[0] = -t143 * t50 - 0.4e1 * t82 * t96 - 0.4e1 / 0.3e1 * t24 * t161;
148
149 t167 = t91 * t75;
150 t171 = 0.1e1 / t64 / t59;
151 t176 = t3 * (0.66059315763395784780e-3 * t87 - 0.24054261223070011357e-4 * t167 * sigma[0] + 0.22022397777557287498e-6 * t30 * t171 * t136);
152 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
153 v3rhosigma2[0] = -0.2e1 * t82 * t107 - 0.4e1 / 0.3e1 * t24 * t176;
154
155 t185 = t3 * (0.48108522446140022714e-5 * t102 * t75 - 0.73407992591857624994e-7 * sigma[0] * t126 * t136);
156 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
157 v3sigma3[0] = -0.4e1 / 0.3e1 * t24 * t185;
158
159 #ifndef XC_DONT_COMPILE_LXC
160
161 if(order < 4) return;
162
163
164 t210 = t64 * t64;
165 t211 = 0.1e1 / t210;
166 t219 = t73 * t73;
167 t221 = 0.1e1 / t14 / t219;
168 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
169 v4rho4[0] = -0.3e1 / 0.4e1 * t4 / t5 / t6 * t18 + 0.2e1 * t4 * t110 * t38 - 0.6e1 * t4 * t53 * t78 - 0.8e1 * t4 * t5 * t139 - 0.4e1 / 0.3e1 * t24 * t3 * (0.15873919443324710989e1 * sigma[0] * t32 * t15 - 0.14665168099473864221e0 * t92 * t35 + 0.49647995164416503440e-2 * t69 * t171 * t75 - 0.71352568799285611491e-4 * t130 * t211 * t136 + 0.37116714382278384028e-6 * t130 * sigma[0] / t210 / t7 * t221);
170
171 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
172 v4rho3sigma[0] = t110 * t2 * t50 / 0.2e1 - 0.3e1 * t143 * t96 - 0.6e1 * t82 * t161 - 0.4e1 / 0.3e1 * t24 * t3 * (-0.26456532405541184982e0 * t45 * t15 + 0.33690251039331850238e-1 * t102 * t35 * sigma[0] - 0.13518494807365346382e-2 * t126 * t75 * t30 + 0.21802173799781714623e-4 * t132 * t136 * t69 - 0.12372238127426128009e-6 * t130 / t210 / t6 * t221);
173
174 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
175 v4rho2sigma2[0] = -t143 * t107 - 0.4e1 * t82 * t176 - 0.4e1 / 0.3e1 * t24 * t3 * (-0.46241521034377049346e-2 * t149 + 0.29827283916606814083e-3 * t152 * sigma[0] - 0.61662713777160404994e-5 * t156 * t136 * t30 + 0.41240793758087093365e-7 * t69 / t210 / rho[0] * t221);
176
177 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
178 v4rhosigma3[0] = -0.2e1 * t82 * t185 - 0.4e1 / 0.3e1 * t24 * t3 * (-0.43297670201526020443e-4 * t167 + 0.15415678444290101249e-5 * t171 * t136 * sigma[0] - 0.13746931252695697789e-7 * t30 * t211 * t221);
179
180 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
181 v4sigma4[0] = -0.4e1 / 0.3e1 * t24 * t3 * (-0.29363197036743049997e-6 * t126 * t136 + 0.45823104175652325962e-8 * sigma[0] * t132 * t221);
182
183 #ifndef XC_DONT_COMPILE_MXC
184
185 if(order < 5) return;
186
187
188 #endif
189
190 #endif
191
192 #endif
193
194 #endif
195
196 #endif
197
198
199 }
200
201
202 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 *))203 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
204 {
205
206 #ifndef XC_DONT_COMPILE_EXC
207 double t1, t2, t3, t4, t5, t6, t7, t8;
208 double t10, t11, t12, t13, t14, t17, t18;
209
210 #ifndef XC_DONT_COMPILE_VXC
211 double t21, t22, t23, t24, t28, t29, t30, t33;
212 double t36, t41, t42, t43, t46;
213
214 #ifndef XC_DONT_COMPILE_FXC
215 double t50, t55, t56, t60, t61, t65, t67, t69;
216 double t71, t74, t81, t85, t86, t89, t95, t99;
217
218 #ifndef XC_DONT_COMPILE_KXC
219 double t103, t115, t119, t121, t125, t128, t137, t140;
220 double t144, t148, t154, t158, t162, t170;
221
222 #ifndef XC_DONT_COMPILE_LXC
223 double t192, t193, t201, t203;
224 #endif
225
226 #endif
227
228 #endif
229
230 #endif
231
232 #endif
233
234
235
236 t1 = sqrt(M_PI);
237 t2 = 0.1e1 / t1;
238 t3 = sqrt(rho[0]);
239 t4 = t2 * t3;
240 t5 = rho[0] * rho[0];
241 t6 = t5 * rho[0];
242 t7 = 0.1e1 / t6;
243 t8 = sigma[0] * t7;
244 t10 = 0.1e1 + 0.8323e-2 * t8;
245 t11 = POW_1_4(t10);
246 t12 = t11 * t11;
247 t13 = t12 * t11;
248 t14 = 0.1e1 / t13;
249 t17 = 0.1e1 + 0.22047110337950987485e-2 * t8 * t14;
250 t18 = t4 * t17;
251 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
252 zk[0] = -0.8e1 / 0.3e1 * t18;
253
254 #ifndef XC_DONT_COMPILE_VXC
255
256 if(order < 1) return;
257
258
259 t21 = t3 * rho[0];
260 t22 = t21 * t2;
261 t23 = t5 * t5;
262 t24 = 0.1e1 / t23;
263 t28 = sigma[0] * sigma[0];
264 t29 = t23 * t6;
265 t30 = 0.1e1 / t29;
266 t33 = 0.1e1 / t13 / t10;
267 t36 = -0.66141331013852962455e-2 * sigma[0] * t24 * t14 + 0.41287072352122365488e-4 * t28 * t30 * t33;
268 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
269 vrho[0] = -0.4e1 * t18 - 0.8e1 / 0.3e1 * t22 * t36;
270
271 t41 = t23 * t5;
272 t42 = 0.1e1 / t41;
273 t43 = sigma[0] * t42;
274 t46 = 0.22047110337950987485e-2 * t7 * t14 - 0.13762357450707455163e-4 * t43 * t33;
275 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
276 vsigma[0] = -0.8e1 / 0.3e1 * t22 * t46;
277
278 #ifndef XC_DONT_COMPILE_FXC
279
280 if(order < 2) return;
281
282
283 t50 = t2 / t3;
284 t55 = t23 * rho[0];
285 t56 = 0.1e1 / t55;
286 t60 = t23 * t23;
287 t61 = 0.1e1 / t60;
288 t65 = t28 * sigma[0];
289 t67 = 0.1e1 / t60 / t6;
290 t69 = t10 * t10;
291 t71 = 0.1e1 / t13 / t69;
292 t74 = 0.26456532405541184982e-1 * sigma[0] * t56 * t14 - 0.41287072352122365489e-3 * t28 * t61 * t33 + 0.18040695917302508518e-5 * t65 * t67 * t71;
293 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
294 v2rho2[0] = -0.2e1 * t50 * t17 - 0.8e1 * t4 * t36 - 0.8e1 / 0.3e1 * t22 * t74;
295
296 t81 = t30 * t33;
297 t85 = 0.1e1 / t60 / t5;
298 t86 = t28 * t85;
299 t89 = -0.66141331013852962455e-2 * t24 * t14 + 0.12386121705636709647e-3 * t81 * sigma[0] - 0.60135653057675028394e-6 * t86 * t71;
300 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
301 v2rhosigma[0] = -0.4e1 * t4 * t46 - 0.8e1 / 0.3e1 * t22 * t89;
302
303 t95 = 0.1e1 / t60 / rho[0];
304 t99 = -0.27524714901414910326e-4 * t42 * t33 + 0.20045217685891676131e-6 * sigma[0] * t95 * t71;
305 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
306 v2sigma2[0] = -0.8e1 / 0.3e1 * t22 * t99;
307
308 #ifndef XC_DONT_COMPILE_KXC
309
310 if(order < 3) return;
311
312
313 t103 = t2 / t21;
314 t115 = 0.1e1 / t60 / t23;
315 t119 = t28 * t28;
316 t121 = 0.1e1 / t60 / t29;
317 t125 = 0.1e1 / t13 / t69 / t10;
318 t128 = -0.13228266202770592491e0 * t43 * t14 + 0.37984106563952576250e-2 * t28 * t95 * t33 - 0.37885461426335267888e-4 * t65 * t115 * t71 + 0.12387598749875974218e-6 * t119 * t121 * t125;
319 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
320 v3rho3[0] = t103 * t17 - 0.6e1 * t50 * t36 - 0.12e2 * t4 * t74 - 0.8e1 / 0.3e1 * t22 * t128;
321
322 t137 = t61 * t33;
323 t140 = t67 * t71;
324 t144 = 0.1e1 / t60 / t41;
325 t148 = 0.26456532405541184982e-1 * t56 * t14 - 0.99088973645093677176e-3 * t137 * sigma[0] + 0.11425774080958255395e-4 * t140 * t28 - 0.41291995832919914059e-7 * t65 * t144 * t125;
326 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
327 v3rho2sigma[0] = -0.2e1 * t50 * t46 - 0.8e1 * t4 * t89 - 0.8e1 / 0.3e1 * t22 * t148;
328
329 t154 = t85 * t71;
330 t158 = 0.1e1 / t60 / t55;
331 t162 = 0.16514828940848946196e-3 * t81 - 0.30067826528837514197e-5 * t154 * sigma[0] + 0.13763998610973304686e-7 * t28 * t158 * t125;
332 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
333 v3rhosigma2[0] = -0.4e1 * t4 * t99 - 0.8e1 / 0.3e1 * t22 * t162;
334
335 t170 = 0.60135653057675028394e-6 * t95 * t71 - 0.45879995369911015621e-8 * sigma[0] * t115 * t125;
336 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
337 v3sigma3[0] = -0.8e1 / 0.3e1 * t22 * t170;
338
339 #ifndef XC_DONT_COMPILE_LXC
340
341 if(order < 4) return;
342
343
344 t192 = t60 * t60;
345 t193 = 0.1e1 / t192;
346 t201 = t69 * t69;
347 t203 = 0.1e1 / t13 / t201;
348 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
349 v4rho4[0] = -0.3e1 / 0.2e1 * t2 / t3 / t5 * t17 + 0.4e1 * t103 * t36 - 0.12e2 * t50 * t74 - 0.16e2 * t4 * t128 - 0.8e1 / 0.3e1 * t22 * (0.79369597216623554946e0 * sigma[0] * t30 * t14 - 0.36662920248684660554e-1 * t86 * t33 + 0.62059993955520629303e-3 * t65 * t158 * t71 - 0.44595355499553507184e-5 * t119 * t193 * t125 + 0.11598973244461995009e-7 * t119 * sigma[0] / t192 / t6 * t203);
350
351 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
352 v4rho3sigma[0] = t103 * t46 - 0.6e1 * t50 * t89 - 0.12e2 * t4 * t148 - 0.8e1 / 0.3e1 * t22 * (-0.13228266202770592491e0 * t42 * t14 + 0.84225627598329625600e-2 * t95 * t33 * sigma[0] - 0.16898118509206682978e-3 * t115 * t71 * t28 + 0.13626358624863571640e-5 * t121 * t125 * t65 - 0.38663244148206650030e-8 * t119 / t192 / t5 * t203);
353
354 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
355 v4rho2sigma2[0] = -0.2e1 * t50 * t99 - 0.8e1 * t4 * t162 - 0.8e1 / 0.3e1 * t22 * (-0.11560380258594262337e-2 * t137 + 0.37284104895758517604e-4 * t140 * sigma[0] - 0.38539196110725253122e-6 * t144 * t125 * t28 + 0.12887748049402216676e-8 * t65 / t192 / rho[0] * t203);
356
357 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
358 v4rhosigma3[0] = -0.4e1 * t4 * t170 - 0.8e1 / 0.3e1 * t22 * (-0.54122087751907525555e-5 * t154 + 0.96347990276813132804e-7 * t158 * t125 * sigma[0] - 0.42959160164674055589e-9 * t28 * t193 * t203);
359
360 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
361 v4sigma4[0] = -0.8e1 / 0.3e1 * t22 * (-0.18351998147964406248e-7 * t115 * t125 + 0.14319720054891351863e-9 * sigma[0] * t121 * t203);
362
363 #ifndef XC_DONT_COMPILE_MXC
364
365 if(order < 5) return;
366
367
368 #endif
369
370 #endif
371
372 #endif
373
374 #endif
375
376 #endif
377
378
379 }
380
381
382 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 *))383 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
384 {
385
386 #ifndef XC_DONT_COMPILE_EXC
387 double t1, t2, t3, t4, t5, t6, t8, t9;
388 double t10, t11, t12, t13, t14, t15, t16, t18;
389 double t19, t20, t21, t22, t25, t26, t27, t29;
390 double t30, t31, t32, t33, t34, t35, t36, t38;
391 double t39, t40, t41, t42, t45, t46, t47;
392
393 #ifndef XC_DONT_COMPILE_VXC
394 double t49, t50, t51, t52, t53, t54, t56, t58;
395 double t59, t60, t61, t62, t63, t64, t65, t69;
396 double t70, t71, t74, t77, t78, t79, t80, t81;
397 double t82, t84, t85, t86, t87, t88, t92, t94;
398 double t95, t96, t98, t99, t100, t101, t105, t106;
399 double t107, t110, t113, t114, t115, t116, t119, t120;
400 double t123, t124, t125, t128, t134, t135, t136, t139;
401
402 #ifndef XC_DONT_COMPILE_FXC
403 double t144, t147, t148, t149, t150, t152, t153, t155;
404 double t156, t158, t159, t160, t161, t162, t163, t165;
405 double t166, t167, t168, t169, t170, t171, t172, t173;
406 double t174, t175, t179, t180, t184, t186, t188, t190;
407 double t193, t194, t195, t196, t197, t198, t199, t201;
408 double t202, t204, t205, t206, t208, t209, t210, t211;
409 double t212, t213, t215, t216, t217, t218, t221, t224;
410 double t227, t228, t230, t234, t235, t236, t237, t240;
411 double t242, t244, t248, t250, t251, t253, t258, t260;
412 double t261, t262, t263, t265, t266, t267, t269, t270;
413 double t271, t273, t274, t275, t277, t278, t279, t280;
414 double t281, t285, t286, t290, t292, t294, t296, t299;
415 double t300, t301, t302, t303, t305, t306, t308, t314;
416 double t318, t319, t322, t326, t327, t328, t332, t336;
417 double t341, t345, t346, t349, t350, t356, t360, t361;
418 double t367, t371, t372;
419
420 #ifndef XC_DONT_COMPILE_KXC
421 double t381, t383, t386, t387, t388, t389, t390, t398;
422 double t402, t404, t408, t411, t412, t413, t414, t415;
423 double t416, t417, t419, t421, t423, t424, t425, t426;
424 double t428, t429, t430, t431, t432, t434, t436, t438;
425 double t439, t440, t441, t444, t445, t447, t448, t450;
426 double t453, t456, t457, t459, t460, t461, t462, t463;
427 double t465, t467, t468, t469, t471, t473, t475, t476;
428 double t477, t478, t483, t485, t486, t488, t489, t490;
429 double t491, t493, t495, t496, t501, t504, t505, t507;
430 double t508, t510, t516, t517, t519, t520, t521, t523;
431 double t525, t527, t529, t532, t533, t535, t536, t539;
432 double t540, t541, t542, t543, t544, t545, t546, t547;
433 double t548, t550, t551, t552, t553, t554, t555, t556;
434 double t557, t558, t560, t561, t563, t564, t566, t567;
435 double t568, t569, t571, t574, t579, t580, t585, t587;
436 double t590, t594, t596, t599, t601, t604, t607, t610;
437 double t612, t614, t616, t617, t619, t624, t625, t627;
438 double t628, t629, t630, t631, t632, t634, t635, t636;
439 double t637, t638, t640, t641, t642, t644, t645, t646;
440 double t648, t651, t667, t671, t673, t677, t680, t681;
441 double t682, t683, t684, t685, t691, t693, t694, t695;
442 double t696, t698, t700, t701, t702, t703, t704, t707;
443 double t709, t710, t712, t713, t715, t716, t718, t719;
444 double t720, t722, t723, t724, t728, t730, t732, t733;
445 double t734, t736, t740, t749, t752, t756, t760, t765;
446 double t767, t769, t770, t779, t780, t782, t783, t784;
447 double t787, t790, t791, t795, t797, t798, t799, t802;
448 double t805, t807, t808, t812, t815, t820, t823, t826;
449 double t831, t834, t838, t842, t843, t846, t848, t849;
450 double t853, t857, t861, t862, t865, t866, t867, t870;
451 double t874, t878, t882, t886, t887, t895, t896, t904;
452 double t905;
453
454 #ifndef XC_DONT_COMPILE_LXC
455 double t919, t926, t929, t931, t933, t935, t938, t948;
456 double t949, t957, t959, t966, t969, t971, t975, t976;
457 double t979, t982, t991, t994, t997, t1003, t1006, t1008;
458 double t1015, t1017, t1019, t1024, t1030, t1033, t1034, t1039;
459 double t1042, t1044, t1048, t1049, t1051, t1054, t1055, t1060;
460 double t1066, t1079, t1084, t1087, t1088, t1089, t1093, t1099;
461 double t1109, t1110, t1113, t1116, t1120, t1121, t1122, t1127;
462 double t1130, t1131, t1133, t1134, t1137, t1145, t1148, t1157;
463 double t1173, t1177, t1178, t1180, t1184, t1185, t1186, t1187;
464 double t1188, t1190, t1194, t1195, t1197, t1198, t1199, t1200;
465 double t1204, t1206, t1215, t1219, t1221, t1222, t1247, t1257;
466 double t1264, t1267, t1268, t1273, t1275, t1278, t1283, t1288;
467 double t1298, t1301, t1303, t1308, t1311, t1316, t1319, t1321;
468 double t1326, t1328, t1330, t1333, t1341, t1344, t1369, t1397;
469 double t1413, t1418, t1419, t1422, t1427, t1431, t1435, t1452;
470 double t1465, t1479, t1493, t1508, t1511, t1516, t1521, t1561;
471 double t1573, t1576, t1579, t1582, t1585, t1588, t1600, t1604;
472 double t1629, t1633, t1653, t1661, t1668, t1677, t1685, t1698;
473 double t1715, t1716, t1724, t1726, t1754, t1771, t1775, t1776;
474 double t1780, t1784, t1808, t1811, t1813, t1816, t1819, t1820;
475 double t1824, t1835, t1838, t1861, t1864, t1866, t1874, t1877;
476 double t1880, t1883, t1884, t1900, t1901, t1904, t1916, t1922;
477 double t1932, t1935, t1936, t1948, t1951, t1952, t1960, t1965;
478 double t1976, t1993, t2001, t2010, t2019, t2021, t2080, t2081;
479 double t2083, t2085, t2089, t2112, t2113, t2115, t2124, t2137;
480 double t2146, t2148, t2185, t2187, t2202, t2203;
481 #endif
482
483 #endif
484
485 #endif
486
487 #endif
488
489 #endif
490
491
492
493 t1 = sqrt(M_PI);
494 t2 = 0.1e1 / t1;
495 t3 = rho[0] - rho[1];
496 t4 = rho[0] + rho[1];
497 t5 = 0.1e1 / t4;
498 t6 = t3 * t5;
499 t8 = 0.1e1 / 0.2e1 + t6 / 0.2e1;
500 t9 = sqrt(t8);
501 t10 = t9 * t8;
502 t11 = t2 * t10;
503 t12 = sqrt(t4);
504 t13 = rho[0] * rho[0];
505 t14 = t13 * rho[0];
506 t15 = 0.1e1 / t14;
507 t16 = sigma[0] * t15;
508 t18 = 0.1e1 + 0.8323e-2 * t16;
509 t19 = POW_1_4(t18);
510 t20 = t19 * t19;
511 t21 = t20 * t19;
512 t22 = 0.1e1 / t21;
513 t25 = 0.1e1 + 0.22047110337950987485e-2 * t16 * t22;
514 t26 = t12 * t25;
515 t27 = t11 * t26;
516 t29 = 0.1e1 / 0.2e1 - t6 / 0.2e1;
517 t30 = sqrt(t29);
518 t31 = t30 * t29;
519 t32 = t2 * t31;
520 t33 = rho[1] * rho[1];
521 t34 = t33 * rho[1];
522 t35 = 0.1e1 / t34;
523 t36 = sigma[2] * t35;
524 t38 = 0.1e1 + 0.8323e-2 * t36;
525 t39 = POW_1_4(t38);
526 t40 = t39 * t39;
527 t41 = t40 * t39;
528 t42 = 0.1e1 / t41;
529 t45 = 0.1e1 + 0.22047110337950987485e-2 * t36 * t42;
530 t46 = t12 * t45;
531 t47 = t32 * t46;
532 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
533 zk[0] = -0.8e1 / 0.3e1 * t27 - 0.8e1 / 0.3e1 * t47;
534
535 #ifndef XC_DONT_COMPILE_VXC
536
537 if(order < 1) return;
538
539
540 t49 = 0.8e1 / 0.3e1 * t27;
541 t50 = 0.8e1 / 0.3e1 * t47;
542 t51 = t2 * t9;
543 t52 = t4 * t4;
544 t53 = 0.1e1 / t52;
545 t54 = t3 * t53;
546 t56 = t5 / 0.2e1 - t54 / 0.2e1;
547 t58 = t51 * t26 * t56;
548 t59 = 0.4e1 * t58;
549 t60 = 0.1e1 / t12;
550 t61 = t60 * t25;
551 t62 = t11 * t61;
552 t63 = 0.4e1 / 0.3e1 * t62;
553 t64 = t13 * t13;
554 t65 = 0.1e1 / t64;
555 t69 = sigma[0] * sigma[0];
556 t70 = t64 * t14;
557 t71 = 0.1e1 / t70;
558 t74 = 0.1e1 / t21 / t18;
559 t77 = -0.66141331013852962455e-2 * sigma[0] * t65 * t22 + 0.41287072352122365488e-4 * t69 * t71 * t74;
560 t78 = t12 * t77;
561 t79 = t11 * t78;
562 t80 = 0.8e1 / 0.3e1 * t79;
563 t81 = t2 * t30;
564 t82 = -t56;
565 t84 = t81 * t46 * t82;
566 t85 = 0.4e1 * t84;
567 t86 = t60 * t45;
568 t87 = t32 * t86;
569 t88 = 0.4e1 / 0.3e1 * t87;
570 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
571 vrho[0] = -t49 - t50 + t4 * (-t59 - t63 - t80 - t85 - t88);
572
573 t92 = -t5 / 0.2e1 - t54 / 0.2e1;
574 t94 = t51 * t26 * t92;
575 t95 = 0.4e1 * t94;
576 t96 = -t92;
577 t98 = t81 * t46 * t96;
578 t99 = 0.4e1 * t98;
579 t100 = t33 * t33;
580 t101 = 0.1e1 / t100;
581 t105 = sigma[2] * sigma[2];
582 t106 = t100 * t34;
583 t107 = 0.1e1 / t106;
584 t110 = 0.1e1 / t41 / t38;
585 t113 = -0.66141331013852962455e-2 * sigma[2] * t101 * t42 + 0.41287072352122365488e-4 * t105 * t107 * t110;
586 t114 = t12 * t113;
587 t115 = t32 * t114;
588 t116 = 0.8e1 / 0.3e1 * t115;
589 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
590 vrho[1] = -t49 - t50 + t4 * (-t95 - t63 - t99 - t88 - t116);
591
592 t119 = t12 * t4;
593 t120 = t119 * t2;
594 t123 = t64 * t13;
595 t124 = 0.1e1 / t123;
596 t125 = sigma[0] * t124;
597 t128 = 0.22047110337950987485e-2 * t15 * t22 - 0.13762357450707455163e-4 * t125 * t74;
598 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
599 vsigma[0] = -0.8e1 / 0.3e1 * t120 * t10 * t128;
600
601 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
602 vsigma[1] = 0.0e0;
603
604 t134 = t100 * t33;
605 t135 = 0.1e1 / t134;
606 t136 = sigma[2] * t135;
607 t139 = 0.22047110337950987485e-2 * t35 * t42 - 0.13762357450707455163e-4 * t136 * t110;
608 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
609 vsigma[2] = -0.8e1 / 0.3e1 * t120 * t31 * t139;
610
611 #ifndef XC_DONT_COMPILE_FXC
612
613 if(order < 2) return;
614
615
616 t144 = 0.8e1 / 0.3e1 * t62;
617 t147 = 0.8e1 / 0.3e1 * t87;
618 t148 = 0.1e1 / t9;
619 t149 = t2 * t148;
620 t150 = t56 * t56;
621 t152 = t149 * t26 * t150;
622 t153 = 0.2e1 * t152;
623 t155 = t51 * t61 * t56;
624 t156 = 0.4e1 * t155;
625 t158 = t51 * t78 * t56;
626 t159 = 0.8e1 * t158;
627 t160 = t52 * t4;
628 t161 = 0.1e1 / t160;
629 t162 = t3 * t161;
630 t163 = -t53 + t162;
631 t165 = t51 * t26 * t163;
632 t166 = 0.4e1 * t165;
633 t167 = 0.1e1 / t119;
634 t168 = t167 * t25;
635 t169 = t11 * t168;
636 t170 = 0.2e1 / 0.3e1 * t169;
637 t171 = t60 * t77;
638 t172 = t11 * t171;
639 t173 = 0.8e1 / 0.3e1 * t172;
640 t174 = t64 * rho[0];
641 t175 = 0.1e1 / t174;
642 t179 = t64 * t64;
643 t180 = 0.1e1 / t179;
644 t184 = t69 * sigma[0];
645 t186 = 0.1e1 / t179 / t14;
646 t188 = t18 * t18;
647 t190 = 0.1e1 / t21 / t188;
648 t193 = 0.26456532405541184982e-1 * sigma[0] * t175 * t22 - 0.41287072352122365489e-3 * t69 * t180 * t74 + 0.18040695917302508518e-5 * t184 * t186 * t190;
649 t194 = t12 * t193;
650 t195 = t11 * t194;
651 t196 = 0.8e1 / 0.3e1 * t195;
652 t197 = 0.1e1 / t30;
653 t198 = t2 * t197;
654 t199 = t82 * t82;
655 t201 = t198 * t46 * t199;
656 t202 = 0.2e1 * t201;
657 t204 = t81 * t86 * t82;
658 t205 = 0.4e1 * t204;
659 t206 = -t163;
660 t208 = t81 * t46 * t206;
661 t209 = 0.4e1 * t208;
662 t210 = t167 * t45;
663 t211 = t32 * t210;
664 t212 = 0.2e1 / 0.3e1 * t211;
665 t213 = -t153 - t156 - t159 - t166 + t170 - t173 - t196 - t202 - t205 - t209 + t212;
666 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
667 v2rho2[0] = -0.8e1 * t58 - t144 - 0.16e2 / 0.3e1 * t79 - 0.8e1 * t84 - t147 + t4 * t213;
668
669 t215 = t149 * t12;
670 t216 = t25 * t92;
671 t217 = t216 * t56;
672 t218 = t215 * t217;
673 t221 = t51 * t61 * t92;
674 t224 = t51 * t78 * t92;
675 t227 = 0.1e1 / t12 / t52;
676 t228 = t227 * t25;
677 t230 = t51 * t228 * t3;
678 t234 = t198 * t12;
679 t235 = t45 * t96;
680 t236 = t235 * t82;
681 t237 = t234 * t236;
682 t240 = t81 * t86 * t96;
683 t242 = t227 * t45;
684 t244 = t81 * t242 * t3;
685 t248 = t81 * t114 * t82;
686 t250 = t60 * t113;
687 t251 = t32 * t250;
688 t253 = -0.2e1 * t218 - 0.2e1 * t221 - 0.4e1 * t224 - 0.4e1 * t230 - 0.2e1 * t155 + t170 - 0.4e1 / 0.3e1 * t172 - 0.2e1 * t237 - 0.2e1 * t240 + 0.4e1 * t244 - 0.2e1 * t204 + t212 - 0.4e1 * t248 - 0.4e1 / 0.3e1 * t251;
689 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
690 v2rho2[1] = t4 * t253 - t116 - t144 - t147 - t59 - t80 - t85 - t95 - t99;
691
692 t258 = t92 * t92;
693 t260 = t149 * t26 * t258;
694 t261 = 0.2e1 * t260;
695 t262 = 0.4e1 * t221;
696 t263 = t53 + t162;
697 t265 = t51 * t26 * t263;
698 t266 = 0.4e1 * t265;
699 t267 = t96 * t96;
700 t269 = t198 * t46 * t267;
701 t270 = 0.2e1 * t269;
702 t271 = 0.4e1 * t240;
703 t273 = t81 * t114 * t96;
704 t274 = 0.8e1 * t273;
705 t275 = -t263;
706 t277 = t81 * t46 * t275;
707 t278 = 0.4e1 * t277;
708 t279 = 0.8e1 / 0.3e1 * t251;
709 t280 = t100 * rho[1];
710 t281 = 0.1e1 / t280;
711 t285 = t100 * t100;
712 t286 = 0.1e1 / t285;
713 t290 = t105 * sigma[2];
714 t292 = 0.1e1 / t285 / t34;
715 t294 = t38 * t38;
716 t296 = 0.1e1 / t41 / t294;
717 t299 = 0.26456532405541184982e-1 * sigma[2] * t281 * t42 - 0.41287072352122365489e-3 * t105 * t286 * t110 + 0.18040695917302508518e-5 * t290 * t292 * t296;
718 t300 = t12 * t299;
719 t301 = t32 * t300;
720 t302 = 0.8e1 / 0.3e1 * t301;
721 t303 = -t261 - t262 - t266 + t170 - t270 - t271 - t274 - t278 + t212 - t279 - t302;
722 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
723 v2rho2[2] = -0.8e1 * t94 - t144 - 0.8e1 * t98 - t147 - 0.16e2 / 0.3e1 * t115 + t4 * t303;
724
725 t305 = t12 * t128;
726 t306 = t11 * t305;
727 t308 = t9 * t128;
728 t314 = t71 * t74;
729 t318 = 0.1e1 / t179 / t13;
730 t319 = t69 * t318;
731 t322 = -0.66141331013852962455e-2 * t65 * t22 + 0.12386121705636709647e-3 * t314 * sigma[0] - 0.60135653057675028394e-6 * t319 * t190;
732 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
733 v2rhosigma[0] = -0.4e1 * t306 - 0.4e1 * t120 * t308 * t56 - 0.8e1 / 0.3e1 * t120 * t10 * t322;
734
735 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
736 v2rhosigma[1] = 0.0e0;
737
738 t326 = t12 * t139;
739 t327 = t32 * t326;
740 t328 = t30 * t139;
741 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
742 v2rhosigma[2] = -0.4e1 * t120 * t328 * t82 - 0.4e1 * t327;
743
744 t332 = t308 * t92;
745 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
746 v2rhosigma[3] = -0.4e1 * t120 * t332 - 0.4e1 * t306;
747
748 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
749 v2rhosigma[4] = 0.0e0;
750
751 t336 = t328 * t96;
752 t341 = t107 * t110;
753 t345 = 0.1e1 / t285 / t33;
754 t346 = t105 * t345;
755 t349 = -0.66141331013852962455e-2 * t101 * t42 + 0.12386121705636709647e-3 * t341 * sigma[2] - 0.60135653057675028394e-6 * t346 * t296;
756 t350 = t31 * t349;
757 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
758 v2rhosigma[5] = -0.4e1 * t327 - 0.4e1 * t120 * t336 - 0.8e1 / 0.3e1 * t120 * t350;
759
760 t356 = 0.1e1 / t179 / rho[0];
761 t360 = -0.27524714901414910326e-4 * t124 * t74 + 0.20045217685891676131e-6 * sigma[0] * t356 * t190;
762 t361 = t10 * t360;
763 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
764 v2sigma2[0] = -0.8e1 / 0.3e1 * t120 * t361;
765
766 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
767 v2sigma2[1] = 0.0e0;
768
769 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
770 v2sigma2[2] = 0.0e0;
771
772 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
773 v2sigma2[3] = 0.0e0;
774
775 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
776 v2sigma2[4] = 0.0e0;
777
778 t367 = 0.1e1 / t285 / rho[1];
779 t371 = -0.27524714901414910326e-4 * t135 * t110 + 0.20045217685891676131e-6 * sigma[2] * t367 * t296;
780 t372 = t31 * t371;
781 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
782 v2sigma2[5] = -0.8e1 / 0.3e1 * t120 * t372;
783
784 #ifndef XC_DONT_COMPILE_KXC
785
786 if(order < 3) return;
787
788
789 t381 = 0.2e1 * t169;
790 t383 = 0.2e1 * t211;
791 t386 = t167 * t77;
792 t387 = t11 * t386;
793 t388 = 0.2e1 * t387;
794 t389 = t60 * t193;
795 t390 = t11 * t389;
796 t398 = 0.1e1 / t179 / t64;
797 t402 = t69 * t69;
798 t404 = 0.1e1 / t179 / t70;
799 t408 = 0.1e1 / t21 / t188 / t18;
800 t411 = -0.13228266202770592491e0 * t125 * t22 + 0.37984106563952576250e-2 * t69 * t356 * t74 - 0.37885461426335267888e-4 * t184 * t398 * t190 + 0.12387598749875974218e-6 * t402 * t404 * t408;
801 t412 = t12 * t411;
802 t413 = t11 * t412;
803 t414 = 0.8e1 / 0.3e1 * t413;
804 t415 = 0.1e1 / t10;
805 t416 = t2 * t415;
806 t417 = t150 * t56;
807 t419 = t416 * t26 * t417;
808 t421 = t149 * t61 * t150;
809 t423 = t25 * t56;
810 t424 = t423 * t163;
811 t425 = t215 * t424;
812 t426 = 0.6e1 * t425;
813 t428 = t51 * t168 * t56;
814 t429 = 0.3e1 * t428;
815 t430 = 0.1e1 / t31;
816 t431 = t2 * t430;
817 t432 = t199 * t82;
818 t434 = t431 * t46 * t432;
819 t436 = t198 * t86 * t199;
820 t438 = t45 * t82;
821 t439 = t438 * t206;
822 t440 = t234 * t439;
823 t441 = 0.6e1 * t440;
824 t444 = t81 * t210 * t82;
825 t445 = 0.3e1 * t444;
826 t447 = t149 * t78 * t150;
827 t448 = 0.6e1 * t447;
828 t450 = t51 * t171 * t56;
829 t453 = t51 * t61 * t163;
830 t456 = t51 * t194 * t56;
831 t457 = 0.12e2 * t456;
832 t459 = t51 * t78 * t163;
833 t460 = 0.12e2 * t459;
834 t461 = t52 * t52;
835 t462 = 0.1e1 / t461;
836 t463 = t3 * t462;
837 t465 = 0.3e1 * t161 - 0.3e1 * t463;
838 t467 = t51 * t26 * t465;
839 t468 = 0.4e1 * t467;
840 t469 = t11 * t228;
841 t471 = t81 * t86 * t206;
842 t473 = -t465;
843 t475 = t81 * t46 * t473;
844 t476 = 0.4e1 * t475;
845 t477 = t32 * t242;
846 t478 = t445 - t448 - 0.12e2 * t450 - 0.6e1 * t453 - t457 - t460 - t468 - t469 - 0.6e1 * t471 - t476 - t477;
847 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
848 v3rho3[0] = -0.6e1 * t152 - 0.12e2 * t155 - 0.6e1 * t201 - 0.12e2 * t204 - 0.24e2 * t158 - 0.12e2 * t165 + t381 - 0.12e2 * t208 + t383 - 0.8e1 * t172 - 0.8e1 * t195 + t4 * (t388 - 0.4e1 * t390 - t414 + t419 - 0.3e1 * t421 - t426 + t429 + t434 - 0.3e1 * t436 - t441 + t478);
849
850 t483 = 0.8e1 * t224;
851 t485 = 0.8e1 * t244;
852 t486 = 0.8e1 * t248;
853 t488 = 0.4e1 * t218;
854 t489 = 0.8e1 * t230;
855 t490 = 0.4e1 * t237;
856 t491 = t51 * t228;
857 t493 = t81 * t242;
858 t495 = t167 * t113;
859 t496 = t32 * t495;
860 t501 = t51 * t194 * t92;
861 t504 = 0.1e1 / t12 / t160;
862 t505 = t504 * t45;
863 t507 = t81 * t505 * t3;
864 t508 = 0.8e1 * t507;
865 t510 = t81 * t114 * t206;
866 t516 = t51 * t168 * t92;
867 t517 = t504 * t25;
868 t519 = t51 * t517 * t3;
869 t520 = 0.8e1 * t519;
870 t521 = t227 * t77;
871 t523 = t51 * t521 * t3;
872 t525 = -0.4e1 * t491 + 0.4e1 * t493 + 0.2e1 / 0.3e1 * t496 - t469 - t477 + 0.4e1 / 0.3e1 * t387 - 0.4e1 / 0.3e1 * t390 - 0.4e1 * t501 - t508 - 0.4e1 * t510 - 0.4e1 * t450 - 0.2e1 * t453 - 0.2e1 * t471 + t516 + t520 - 0.8e1 * t523;
873 t527 = t81 * t210 * t96;
874 t529 = t198 * t114 * t199;
875 t532 = t81 * t250 * t82;
876 t533 = 0.4e1 * t532;
877 t535 = t51 * t171 * t92;
878 t536 = 0.4e1 * t535;
879 t539 = t416 * t12;
880 t540 = t216 * t150;
881 t541 = t539 * t540;
882 t542 = t149 * t60;
883 t543 = t542 * t217;
884 t544 = 0.2e1 * t543;
885 t545 = t149 * t227;
886 t546 = t25 * t3;
887 t547 = t546 * t56;
888 t548 = t545 * t547;
889 t550 = t431 * t12;
890 t551 = t235 * t199;
891 t552 = t550 * t551;
892 t553 = t198 * t60;
893 t554 = t553 * t236;
894 t555 = 0.2e1 * t554;
895 t556 = t77 * t92;
896 t557 = t556 * t56;
897 t558 = t215 * t557;
898 t560 = t216 * t163;
899 t561 = t215 * t560;
900 t563 = t235 * t206;
901 t564 = t234 * t563;
902 t566 = t198 * t227;
903 t567 = t45 * t3;
904 t568 = t567 * t82;
905 t569 = t566 * t568;
906 t571 = t527 - 0.2e1 * t529 - t533 - t536 - t421 + 0.2e1 * t428 - t436 + 0.2e1 * t444 + t541 - t544 - 0.4e1 * t548 + t552 - t555 - 0.4e1 * t558 - 0.2e1 * t561 - 0.2e1 * t564 + 0.4e1 * t569;
907 t574 = t485 - t486 - t279 - 0.16e2 / 0.3e1 * t172 - t196 - t488 - t262 - t489 - t490 - t271 + t4 * (t525 + t571);
908 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
909 v3rho3[1] = -t153 - 0.8e1 * t155 - t202 - 0.8e1 * t204 - t159 - t166 + t381 - t209 + t383 - t483 + t574;
910
911 t579 = t60 * t299;
912 t580 = t32 * t579;
913 t585 = t161 + 0.3e1 * t463;
914 t587 = t81 * t46 * t585;
915 t590 = t81 * t300 * t82;
916 t594 = t149 * t61 * t258;
917 t596 = t51 * t61 * t263;
918 t599 = t198 * t86 * t267;
919 t601 = t81 * t250 * t96;
920 t604 = t81 * t86 * t275;
921 t607 = t149 * t78 * t258;
922 t610 = t51 * t78 * t263;
923 t612 = -t585;
924 t614 = t51 * t26 * t612;
925 t616 = -0.4e1 / 0.3e1 * t580 + 0.4e1 / 0.3e1 * t496 - t469 - t477 + 0.2e1 / 0.3e1 * t387 - 0.4e1 * t587 - 0.4e1 * t590 + 0.4e1 * t507 - t594 - 0.2e1 * t596 - t599 - 0.4e1 * t601 - 0.2e1 * t604 - 0.2e1 * t607 - 0.4e1 * t610 - 0.4e1 * t614;
926 t617 = t227 * t113;
927 t619 = t81 * t617 * t3;
928 t624 = t216 * t3;
929 t625 = t545 * t624;
930 t627 = t25 * t258;
931 t628 = t627 * t56;
932 t629 = t539 * t628;
933 t630 = t25 * t263;
934 t631 = t630 * t56;
935 t632 = t215 * t631;
936 t634 = t45 * t267;
937 t635 = t634 * t82;
938 t636 = t550 * t635;
939 t637 = t235 * t3;
940 t638 = t566 * t637;
941 t640 = t113 * t96;
942 t641 = t640 * t82;
943 t642 = t234 * t641;
944 t644 = t45 * t275;
945 t645 = t644 * t82;
946 t646 = t234 * t645;
947 t648 = 0.8e1 * t619 + 0.2e1 * t516 - 0.4e1 * t519 + 0.2e1 * t527 - t533 - t536 + t428 + t444 - 0.4e1 * t625 + t629 - 0.2e1 * t632 + t636 + 0.4e1 * t638 - 0.4e1 * t642 - 0.2e1 * t646 - t544 - t555;
948 t651 = -t278 - t173 - t302 - t488 - 0.8e1 * t221 - t489 - t490 - 0.8e1 * t240 - t261 - t270 + t4 * (t616 + t648);
949 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
950 v3rho3[2] = -t156 - t205 + t381 + t383 - t483 + t485 - t486 - 0.16e2 / 0.3e1 * t251 - t266 - t274 + t651;
951
952 t667 = 0.1e1 / t285 / t100;
953 t671 = t105 * t105;
954 t673 = 0.1e1 / t285 / t106;
955 t677 = 0.1e1 / t41 / t294 / t38;
956 t680 = -0.13228266202770592491e0 * t136 * t42 + 0.37984106563952576250e-2 * t105 * t367 * t110 - 0.37885461426335267888e-4 * t290 * t667 * t296 + 0.12387598749875974218e-6 * t671 * t673 * t677;
957 t681 = t12 * t680;
958 t682 = t32 * t681;
959 t683 = 0.8e1 / 0.3e1 * t682;
960 t684 = 0.3e1 * t516;
961 t685 = 0.3e1 * t527;
962 t691 = t258 * t92;
963 t693 = t416 * t26 * t691;
964 t694 = t216 * t263;
965 t695 = t215 * t694;
966 t696 = 0.6e1 * t695;
967 t698 = t267 * t96;
968 t700 = t431 * t46 * t698;
969 t701 = t235 * t275;
970 t702 = t234 * t701;
971 t703 = 0.6e1 * t702;
972 t704 = 0.2e1 * t496;
973 t707 = -0.3e1 * t161 - 0.3e1 * t463;
974 t709 = t51 * t26 * t707;
975 t710 = 0.4e1 * t709;
976 t712 = t198 * t114 * t267;
977 t713 = 0.6e1 * t712;
978 t715 = t81 * t300 * t96;
979 t716 = 0.12e2 * t715;
980 t718 = t81 * t114 * t275;
981 t719 = 0.12e2 * t718;
982 t720 = -t707;
983 t722 = t81 * t46 * t720;
984 t723 = 0.4e1 * t722;
985 t724 = t700 - t703 - t469 - t477 + t704 - 0.4e1 * t580 - t710 - t713 - t716 - t719 - t723;
986 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
987 v3rho3[3] = t381 + t383 - 0.8e1 * t251 - 0.12e2 * t265 - 0.24e2 * t273 - 0.12e2 * t277 - 0.8e1 * t301 - 0.12e2 * t221 - 0.12e2 * t240 - 0.6e1 * t260 - 0.6e1 * t269 + t4 * (-t683 + t684 + t685 - 0.3e1 * t594 - 0.6e1 * t596 - 0.3e1 * t599 - 0.12e2 * t601 - 0.6e1 * t604 + t693 - t696 + t724);
988
989 t728 = t51 * t305 * t56;
990 t730 = t60 * t128;
991 t732 = 0.2e1 * t11 * t730;
992 t733 = t12 * t322;
993 t734 = t11 * t733;
994 t736 = t148 * t128;
995 t740 = t9 * t322;
996 t749 = t180 * t74;
997 t752 = t186 * t190;
998 t756 = 0.1e1 / t179 / t123;
999 t760 = 0.26456532405541184982e-1 * t175 * t22 - 0.99088973645093677176e-3 * t749 * sigma[0] + 0.11425774080958255395e-4 * t752 * t69 - 0.41291995832919914059e-7 * t184 * t756 * t408;
1000 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1001 v3rho2sigma[0] = -0.12e2 * t728 - t732 - 0.8e1 * t734 - 0.2e1 * t120 * t736 * t150 - 0.8e1 * t120 * t740 * t56 - 0.4e1 * t120 * t308 * t163 - 0.8e1 / 0.3e1 * t120 * t10 * t760;
1002
1003 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1004 v3rho2sigma[1] = 0.0e0;
1005
1006 t765 = t81 * t326 * t82;
1007 t767 = t60 * t139;
1008 t769 = 0.2e1 * t32 * t767;
1009 t770 = t197 * t139;
1010 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1011 v3rho2sigma[2] = -0.2e1 * t120 * t770 * t199 - 0.4e1 * t120 * t328 * t206 - 0.12e2 * t765 - t769;
1012
1013 t779 = t12 * t2;
1014 t780 = t779 * t332;
1015 t782 = t120 * t148;
1016 t783 = t128 * t92;
1017 t784 = t783 * t56;
1018 t787 = t740 * t92;
1019 t790 = t167 * t2;
1020 t791 = t308 * t3;
1021 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1022 v3rho2sigma[3] = -0.4e1 * t120 * t787 - 0.2e1 * t782 * t784 - 0.4e1 * t790 * t791 - 0.6e1 * t728 - t732 - 0.4e1 * t734 - 0.6e1 * t780;
1023
1024 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1025 v3rho2sigma[4] = 0.0e0;
1026
1027 t795 = t779 * t336;
1028 t797 = t120 * t197;
1029 t798 = t139 * t96;
1030 t799 = t798 * t82;
1031 t802 = t328 * t3;
1032 t805 = t779 * t350;
1033 t807 = t30 * t349;
1034 t808 = t807 * t82;
1035 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1036 v3rho2sigma[5] = -0.4e1 * t120 * t808 + 0.4e1 * t790 * t802 - 0.2e1 * t797 * t799 - 0.6e1 * t765 - t769 - 0.6e1 * t795 - 0.4e1 * t805;
1037
1038 t812 = t736 * t258;
1039 t815 = t308 * t263;
1040 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1041 v3rho2sigma[6] = -0.2e1 * t120 * t812 - 0.4e1 * t120 * t815 - t732 - 0.12e2 * t780;
1042
1043 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1044 v3rho2sigma[7] = 0.0e0;
1045
1046 t820 = t770 * t267;
1047 t823 = t807 * t96;
1048 t826 = t328 * t275;
1049 t831 = t286 * t110;
1050 t834 = t292 * t296;
1051 t838 = 0.1e1 / t285 / t134;
1052 t842 = 0.26456532405541184982e-1 * t281 * t42 - 0.99088973645093677176e-3 * t831 * sigma[2] + 0.11425774080958255395e-4 * t834 * t105 - 0.41291995832919914059e-7 * t290 * t838 * t677;
1053 t843 = t31 * t842;
1054 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1055 v3rho2sigma[8] = -0.12e2 * t795 - t769 - 0.8e1 * t805 - 0.2e1 * t120 * t820 - 0.8e1 * t120 * t823 - 0.4e1 * t120 * t826 - 0.8e1 / 0.3e1 * t120 * t843;
1056
1057 t846 = t779 * t361;
1058 t848 = t9 * t360;
1059 t849 = t848 * t56;
1060 t853 = t318 * t190;
1061 t857 = 0.1e1 / t179 / t174;
1062 t861 = 0.16514828940848946196e-3 * t314 - 0.30067826528837514197e-5 * t853 * sigma[0] + 0.13763998610973304686e-7 * t69 * t857 * t408;
1063 t862 = t10 * t861;
1064 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1065 v3rhosigma2[0] = -0.4e1 * t846 - 0.4e1 * t120 * t849 - 0.8e1 / 0.3e1 * t120 * t862;
1066
1067 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1068 v3rhosigma2[1] = 0.0e0;
1069
1070 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1071 v3rhosigma2[2] = 0.0e0;
1072
1073 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1074 v3rhosigma2[3] = 0.0e0;
1075
1076 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1077 v3rhosigma2[4] = 0.0e0;
1078
1079 t865 = t779 * t372;
1080 t866 = t30 * t371;
1081 t867 = t866 * t82;
1082 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1083 v3rhosigma2[5] = -0.4e1 * t120 * t867 - 0.4e1 * t865;
1084
1085 t870 = t848 * t92;
1086 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1087 v3rhosigma2[6] = -0.4e1 * t120 * t870 - 0.4e1 * t846;
1088
1089 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1090 v3rhosigma2[7] = 0.0e0;
1091
1092 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1093 v3rhosigma2[8] = 0.0e0;
1094
1095 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1096 v3rhosigma2[9] = 0.0e0;
1097
1098 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1099 v3rhosigma2[10] = 0.0e0;
1100
1101 t874 = t866 * t96;
1102 t878 = t345 * t296;
1103 t882 = 0.1e1 / t285 / t280;
1104 t886 = 0.16514828940848946196e-3 * t341 - 0.30067826528837514197e-5 * t878 * sigma[2] + 0.13763998610973304686e-7 * t105 * t882 * t677;
1105 t887 = t31 * t886;
1106 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1107 v3rhosigma2[11] = -0.4e1 * t865 - 0.4e1 * t120 * t874 - 0.8e1 / 0.3e1 * t120 * t887;
1108
1109 t895 = 0.60135653057675028394e-6 * t356 * t190 - 0.45879995369911015621e-8 * sigma[0] * t398 * t408;
1110 t896 = t10 * t895;
1111 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1112 v3sigma3[0] = -0.8e1 / 0.3e1 * t120 * t896;
1113
1114 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1115 v3sigma3[1] = 0.0e0;
1116
1117 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1118 v3sigma3[2] = 0.0e0;
1119
1120 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1121 v3sigma3[3] = 0.0e0;
1122
1123 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1124 v3sigma3[4] = 0.0e0;
1125
1126 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1127 v3sigma3[5] = 0.0e0;
1128
1129 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1130 v3sigma3[6] = 0.0e0;
1131
1132 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1133 v3sigma3[7] = 0.0e0;
1134
1135 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1136 v3sigma3[8] = 0.0e0;
1137
1138 t904 = 0.60135653057675028394e-6 * t367 * t296 - 0.45879995369911015621e-8 * sigma[2] * t667 * t677;
1139 t905 = t31 * t904;
1140 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1141 v3sigma3[9] = -0.8e1 / 0.3e1 * t120 * t905;
1142
1143 #ifndef XC_DONT_COMPILE_LXC
1144
1145 if(order < 4) return;
1146
1147
1148 t919 = 0.8e1 * t387 - 0.16e2 * t390 - 0.32e2 / 0.3e1 * t413 + 0.4e1 * t419 - 0.12e2 * t421 - 0.24e2 * t425 + 0.12e2 * t428 + 0.4e1 * t434 - 0.12e2 * t436 - 0.24e2 * t440 + 0.12e2 * t444;
1149 t926 = 0.4e1 * t469;
1150 t929 = 0.4e1 * t477;
1151 t931 = 0.5e1 / 0.2e1 * t11 * t517;
1152 t933 = 0.5e1 / 0.2e1 * t32 * t505;
1153 t935 = t11 * t167 * t193;
1154 t938 = t11 * t60 * t411;
1155 t948 = t179 * t179;
1156 t949 = 0.1e1 / t948;
1157 t957 = t188 * t188;
1158 t959 = 0.1e1 / t21 / t957;
1159 t966 = t11 * t521;
1160 t969 = t198 * t210 * t199;
1161 t971 = t206 * t206;
1162 t975 = t242 * t82;
1163 t976 = t81 * t975;
1164 t979 = t51 * t386 * t56;
1165 t982 = t51 * t389 * t56;
1166 t991 = t149 * t171 * t150;
1167 t994 = t51 * t168 * t163;
1168 t997 = t81 * t210 * t206;
1169 t1003 = t51 * t171 * t163;
1170 t1006 = t51 * t61 * t465;
1171 t1008 = t931 + t933 + 0.4e1 * t935 - 0.16e2 / 0.3e1 * t938 - 0.8e1 / 0.3e1 * t11 * t12 * (0.79369597216623554946e0 * sigma[0] * t71 * t22 - 0.36662920248684660554e-1 * t319 * t74 + 0.62059993955520629303e-3 * t184 * t857 * t190 - 0.44595355499553507184e-5 * t402 * t949 * t408 + 0.11598973244461995009e-7 * t402 * sigma[0] / t948 / t14 * t959) - 0.4e1 * t966 + 0.3e1 * t969 - 0.6e1 * t198 * t46 * t971 - 0.6e1 * t976 + 0.12e2 * t979 - 0.24e2 * t982 - 0.16e2 * t51 * t412 * t56 + 0.4e1 * t416 * t78 * t417 - 0.12e2 * t991 + 0.6e1 * t994 + 0.6e1 * t997 - 0.12e2 * t149 * t194 * t150 - 0.24e2 * t1003 - 0.8e1 * t1006;
1172 t1015 = t461 * t4;
1173 t1017 = t3 / t1015;
1174 t1019 = -0.12e2 * t462 + 0.12e2 * t1017;
1175 t1024 = t81 * t86 * t473;
1176 t1030 = t8 * t8;
1177 t1033 = t2 / t9 / t1030;
1178 t1034 = t150 * t150;
1179 t1039 = t416 * t61 * t417;
1180 t1042 = t149 * t168 * t150;
1181 t1044 = t163 * t163;
1182 t1048 = t228 * t56;
1183 t1049 = t51 * t1048;
1184 t1051 = t29 * t29;
1185 t1054 = t2 / t30 / t1051;
1186 t1055 = t199 * t199;
1187 t1060 = t431 * t86 * t432;
1188 t1066 = t542 * t424;
1189 t1079 = t553 * t439;
1190 t1084 = -0.24e2 * t51 * t194 * t163 - 0.16e2 * t51 * t78 * t465 - 0.4e1 * t51 * t26 * t1019 - 0.8e1 * t1024 + 0.4e1 * t81 * t46 * t1019 - 0.3e1 / 0.2e1 * t1033 * t26 * t1034 + 0.2e1 * t1039 + 0.3e1 * t1042 - 0.6e1 * t149 * t26 * t1044 - 0.6e1 * t1049 - 0.3e1 / 0.2e1 * t1054 * t46 * t1055 + 0.2e1 * t1060 + 0.6e1 * t539 * t25 * t150 * t163 - 0.12e2 * t1066 - 0.24e2 * t215 * t77 * t56 * t163 - 0.8e1 * t215 * t423 * t465 + 0.6e1 * t550 * t45 * t199 * t206 - 0.12e2 * t1079 - 0.8e1 * t234 * t438 * t473;
1191 t1087 = -0.24e2 * t447 - 0.48e2 * t450 - 0.24e2 * t453 - 0.48e2 * t456 - 0.48e2 * t459 - 0.16e2 * t467 - t926 - 0.24e2 * t471 - 0.16e2 * t475 - t929 + t4 * (t1008 + t1084);
1192 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1193 v4rho4[0] = t919 + t1087;
1194
1195 t1088 = t32 * t617;
1196 t1089 = t51 * t517;
1197 t1093 = t81 * t505;
1198 t1099 = t51 * t389 * t92;
1199 t1109 = t228 * t92;
1200 t1110 = t51 * t1109;
1201 t1113 = 0.1e1 / t12 / t461;
1202 t1116 = t51 * t1113 * t25 * t3;
1203 t1120 = -t1088 + 0.18e2 * t1089 - 0.12e2 * t51 * t521 - 0.18e2 * t1093 + t931 + t933 + 0.2e1 * t935 - 0.4e1 / 0.3e1 * t938 - 0.3e1 * t966 - 0.6e1 * t1099 - 0.4e1 * t51 * t412 * t92 - 0.4e1 * t81 * t114 * t473 - 0.6e1 * t149 * t1048 - 0.3e1 / 0.2e1 * t1110 - 0.27e2 * t1116 + 0.6e1 * t198 * t975;
1204 t1121 = t242 * t96;
1205 t1122 = t81 * t1121;
1206 t1127 = t198 * t250 * t199;
1207 t1130 = t81 * t495 * t82;
1208 t1131 = 0.3e1 * t1130;
1209 t1133 = t51 * t386 * t92;
1210 t1134 = 0.3e1 * t1133;
1211 t1137 = t51 * t504 * t77 * t3;
1212 t1145 = t81 * t1113 * t45 * t3;
1213 t1148 = t81 * t250 * t206;
1214 t1157 = -0.3e1 / 0.2e1 * t1122 + t431 * t114 * t432 - 0.3e1 * t1127 + t1131 + t1134 + 0.24e2 * t1137 - 0.12e2 * t51 * t227 * t193 * t3 + 0.27e2 * t1145 - 0.6e1 * t1148 + 0.3e1 / 0.2e1 * t969 - 0.9e1 / 0.2e1 * t976 + 0.6e1 * t979 - 0.6e1 * t982 - 0.3e1 * t991 + 0.3e1 * t994 + 0.3e1 * t997;
1215 t1173 = t1033 * t12;
1216 t1177 = t416 * t60;
1217 t1178 = t1177 * t540;
1218 t1180 = t416 * t227;
1219 t1184 = t149 * t167;
1220 t1185 = t1184 * t217;
1221 t1186 = 0.3e1 / 0.2e1 * t1185;
1222 t1187 = t149 * t504;
1223 t1188 = t1187 * t547;
1224 t1190 = t1054 * t12;
1225 t1194 = t431 * t60;
1226 t1195 = t1194 * t551;
1227 t1197 = -0.6e1 * t1003 - 0.2e1 * t1006 - 0.2e1 * t1024 + t1039 / 0.2e1 + 0.3e1 / 0.2e1 * t1042 - 0.9e1 / 0.2e1 * t1049 + t1060 / 0.2e1 + 0.6e1 * t566 * t567 * t206 - 0.6e1 * t234 * t113 * t82 * t206 - 0.3e1 / 0.2e1 * t1173 * t216 * t417 + 0.3e1 / 0.2e1 * t1178 + 0.3e1 * t1180 * t546 * t150 + t1186 + 0.12e2 * t1188 - 0.3e1 / 0.2e1 * t1190 * t235 * t432 + 0.3e1 / 0.2e1 * t1195;
1228 t1198 = t198 * t167;
1229 t1199 = t1198 * t236;
1230 t1200 = 0.3e1 / 0.2e1 * t1199;
1231 t1204 = t542 * t557;
1232 t1206 = t542 * t560;
1233 t1215 = t431 * t227;
1234 t1219 = t553 * t563;
1235 t1221 = t198 * t504;
1236 t1222 = t1221 * t568;
1237 t1247 = 0.3e1 * t539 * t216 * t56 * t163 - 0.6e1 * t215 * t193 * t92 * t56 + 0.3e1 * t550 * t235 * t82 * t206 - 0.12e2 * t545 * t77 * t3 * t56 - 0.3e1 * t1215 * t567 * t199 + 0.3e1 * t539 * t556 * t150 - 0.6e1 * t215 * t556 * t163 - 0.6e1 * t545 * t546 * t163 - 0.2e1 * t215 * t216 * t465 - 0.2e1 * t234 * t235 * t473 - 0.3e1 * t1066 - 0.3e1 * t1079 + t1200 - 0.6e1 * t1204 - 0.3e1 * t1206 - 0.3e1 * t1219 - 0.12e2 * t1222;
1238 t1257 = t4 * (t1120 + t1157 + t1197 + t1247) - 0.12e2 * t491 + 0.12e2 * t493 + t704 - t926 - t929 + 0.6e1 * t387 - 0.8e1 * t390 - t414 - 0.12e2 * t501 - 0.24e2 * t507;
1239 t1264 = -0.12e2 * t510 - 0.24e2 * t450 - 0.12e2 * t453 - t457 - t460 - t468 - 0.12e2 * t471 - t476 + t684 + 0.24e2 * t519 - 0.24e2 * t523;
1240 t1267 = 0.12e2 * t532;
1241 t1268 = 0.12e2 * t535;
1242 t1273 = t685 - 0.6e1 * t529 - t1267 - t1268 + t419 - 0.6e1 * t421 + 0.9e1 * t428 + t434 - 0.6e1 * t436 + 0.9e1 * t444 - t448;
1243 t1275 = 0.6e1 * t543;
1244 t1278 = 0.6e1 * t554;
1245 t1283 = 0.3e1 * t541 - t1275 - 0.12e2 * t548 + 0.3e1 * t552 - t1278 - 0.12e2 * t558 - 0.6e1 * t561 - 0.6e1 * t564 + 0.12e2 * t569 - t426 - t441;
1246 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1247 v4rho4[1] = t1257 + t1264 + t1273 + t1283;
1248
1249 t1288 = t32 * t167 * t299;
1250 t1298 = t149 * t168 * t258;
1251 t1301 = t51 * t168 * t263;
1252 t1303 = t198 * t210 * t267;
1253 t1308 = t81 * t495 * t96;
1254 t1311 = t81 * t210 * t275;
1255 t1316 = t81 * t579 * t82;
1256 t1319 = 0.1e1 / t12 / t1015;
1257 t1321 = t3 * t3;
1258 t1326 = t149 * t171 * t258;
1259 t1328 = 0.8e1 * t81 * t617 + 0.2e1 / 0.3e1 * t1288 - 0.2e1 * t1088 - 0.4e1 * t1089 + 0.4e1 * t1093 + t931 + t933 + 0.2e1 / 0.3e1 * t935 - 0.2e1 * t966 - 0.4e1 * t149 * t1109 + t1298 / 0.2e1 + t1301 + t1303 / 0.2e1 + 0.4e1 * t198 * t1121 + 0.2e1 * t1308 + t1311 - 0.2e1 * t198 * t300 * t199 - 0.4e1 * t1316 - 0.4e1 * t149 * t1319 * t25 * t1321 - 0.2e1 * t1326;
1260 t1330 = t51 * t171 * t263;
1261 t1333 = t51 * t61 * t612;
1262 t1341 = t81 * t504 * t113 * t3;
1263 t1344 = t81 * t86 * t585;
1264 t1369 = -0.4e1 * t1330 - 0.4e1 * t1333 - 0.4e1 * t198 * t1319 * t45 * t1321 - 0.16e2 * t1341 - 0.4e1 * t1344 - 0.2e1 * t149 * t194 * t258 - 0.4e1 * t51 * t194 * t263 - 0.8e1 * t51 * t78 * t612 - 0.4e1 * t81 * t300 * t206 - 0.4e1 * t1099 - 0.3e1 * t1110 - 0.32e2 * t1116 - 0.3e1 * t1122 - 0.2e1 * t1127 + 0.4e1 * t1130 + 0.4e1 * t1133 - 0.8e1 * t1137 + 0.32e2 * t1145 - 0.4e1 * t1148 + t969 / 0.2e1;
1265 t1397 = t1221 * t637;
1266 t1413 = t1187 * t624;
1267 t1418 = -0.3e1 * t976 + 0.2e1 * t979 + t994 + t997 + t1042 / 0.2e1 - 0.3e1 * t1049 - 0.8e1 * t545 * t556 * t3 + 0.2e1 * t539 * t77 * t258 * t56 + t539 * t627 * t163 - 0.4e1 * t215 * t77 * t263 * t56 - 0.4e1 * t215 * t25 * t612 * t56 - 0.2e1 * t215 * t630 * t163 + t550 * t634 * t206 - 0.8e1 * t1397 - 0.4e1 * t234 * t640 * t206 + 0.8e1 * t566 * t113 * t3 * t82 - 0.4e1 * t234 * t45 * t585 * t82 - 0.2e1 * t234 * t644 * t206 + 0.8e1 * t1413 - 0.3e1 / 0.2e1 * t1173 * t627 * t150;
1268 t1419 = t1177 * t628;
1269 t1422 = t542 * t631;
1270 t1427 = t1194 * t635;
1271 t1431 = t553 * t641;
1272 t1435 = t553 * t645;
1273 t1452 = t1419 + t539 * t630 * t150 - 0.2e1 * t1422 - 0.3e1 / 0.2e1 * t1190 * t634 * t199 + t1427 + 0.2e1 * t550 * t640 * t199 - 0.4e1 * t1431 + t550 * t644 * t199 - 0.2e1 * t1435 + t1178 + 0.2e1 * t1185 - 0.4e1 * t1188 + t1195 + 0.2e1 * t1199 - 0.4e1 * t1204 - 0.2e1 * t1206 - 0.2e1 * t1219 + 0.4e1 * t1222 + 0.4e1 * t1180 * t216 * t3 * t56 - 0.4e1 * t1215 * t235 * t82 * t3;
1274 t1465 = t4 * (t1328 + t1369 + t1418 + t1452) - 0.8e1 / 0.3e1 * t580 - 0.8e1 * t491 + 0.8e1 * t493 + 0.4e1 * t496 - t926 - t929 + 0.4e1 * t387 - 0.8e1 / 0.3e1 * t390 - 0.8e1 * t587 - 0.8e1 * t590 - 0.8e1 * t501 - t508;
1275 t1479 = -0.8e1 * t510 - 0.2e1 * t594 - 0.4e1 * t596 - 0.2e1 * t599 - 0.8e1 * t601 - 0.4e1 * t604 - 0.4e1 * t607 - 0.8e1 * t610 - 0.8e1 * t614 + 0.16e2 * t619 - 0.8e1 * t450 - 0.4e1 * t453 - 0.4e1 * t471;
1276 t1493 = 0.6e1 * t516 + t520 - 0.16e2 * t523 + 0.6e1 * t527 - 0.4e1 * t529 - 0.16e2 * t532 - 0.16e2 * t535 - 0.2e1 * t421 + 0.6e1 * t428 - 0.2e1 * t436 + 0.6e1 * t444 - 0.8e1 * t625 + 0.2e1 * t629;
1277 t1508 = -0.4e1 * t632 + 0.2e1 * t636 + 0.8e1 * t638 - 0.8e1 * t642 - 0.4e1 * t646 + 0.2e1 * t541 - 0.8e1 * t543 - 0.8e1 * t548 + 0.2e1 * t552 - 0.8e1 * t554 - 0.8e1 * t558 - 0.4e1 * t561 - 0.4e1 * t564 + 0.8e1 * t569;
1278 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1279 v4rho4[2] = t1465 + t1479 + t1493 + t1508;
1280
1281 t1511 = t32 * t60 * t680;
1282 t1516 = t542 * t694;
1283 t1521 = t553 * t701;
1284 t1561 = -0.4e1 / 0.3e1 * t1511 + 0.3e1 * t1180 * t627 * t3 - 0.3e1 * t1516 - 0.6e1 * t545 * t546 * t263 - 0.3e1 * t1521 - 0.3e1 / 0.2e1 * t1173 * t25 * t691 * t56 - 0.6e1 * t215 * t556 * t263 - 0.6e1 * t215 * t216 * t612 - 0.3e1 / 0.2e1 * t1190 * t45 * t698 * t82 - 0.3e1 * t1215 * t634 * t3 - 0.6e1 * t234 * t235 * t585 + 0.6e1 * t566 * t567 * t275 - 0.2e1 * t215 * t25 * t707 * t56 + 0.3e1 * t550 * t113 * t267 * t82 + 0.12e2 * t566 * t640 * t3 - 0.6e1 * t234 * t299 * t96 * t82;
1285 t1573 = t431 * t86 * t698;
1286 t1576 = t51 * t61 * t707;
1287 t1579 = t198 * t250 * t267;
1288 t1582 = t81 * t579 * t96;
1289 t1585 = t81 * t250 * t275;
1290 t1588 = t81 * t86 * t720;
1291 t1600 = 0.6e1 * t462 + 0.12e2 * t1017;
1292 t1604 = -0.6e1 * t234 * t113 * t275 * t82 - 0.2e1 * t234 * t45 * t720 * t82 + 0.2e1 * t1288 - 0.3e1 * t1088 + t931 + t933 - t966 + t1573 / 0.2e1 - 0.2e1 * t1576 - 0.3e1 * t1579 - 0.6e1 * t1582 - 0.6e1 * t1585 - 0.2e1 * t1588 - 0.4e1 * t81 * t681 * t82 + t416 * t78 * t691 - 0.4e1 * t51 * t78 * t707 - 0.4e1 * t51 * t26 * t1600;
1293 t1629 = t416 * t61 * t691;
1294 t1633 = 0.12e2 * t81 * t227 * t299 * t3 - 0.12e2 * t81 * t114 * t585 + 0.4e1 * t81 * t46 * t1600 + 0.3e1 / 0.2e1 * t1298 + 0.3e1 * t1301 + 0.3e1 / 0.2e1 * t1303 + 0.6e1 * t1308 + 0.3e1 * t1311 - 0.6e1 * t1316 - 0.3e1 * t1326 - 0.6e1 * t1330 - 0.6e1 * t1333 + 0.12e2 * t1341 - 0.6e1 * t1344 + t1629 / 0.2e1 - 0.9e1 / 0.2e1 * t1110 + 0.3e1 * t1116;
1295 t1653 = -0.9e1 / 0.2e1 * t1122 + t1131 + t1134 - 0.3e1 * t1145 - 0.3e1 / 0.2e1 * t976 - 0.3e1 / 0.2e1 * t1049 + 0.6e1 * t1397 - 0.6e1 * t1413 + 0.3e1 / 0.2e1 * t1419 - 0.3e1 * t1422 + 0.3e1 / 0.2e1 * t1427 - 0.6e1 * t1431 - 0.3e1 * t1435 + t1186 + t1200 + 0.3e1 * t539 * t216 * t263 * t56 + 0.3e1 * t550 * t235 * t275 * t82;
1296 t1661 = t4 * (t1561 + t1604 + t1633 + t1653) - 0.8e1 * t580 - t683 + 0.6e1 * t496 - t926 - t929 + t388 - 0.12e2 * t587 - 0.12e2 * t590 + t693 + t700;
1297 t1668 = -t710 - t713 - t716 - t719 - t723 + 0.12e2 * t507 - 0.6e1 * t594 - 0.12e2 * t596 - 0.6e1 * t599 - 0.24e2 * t601 - 0.12e2 * t604;
1298 t1677 = -0.6e1 * t607 - 0.12e2 * t610 - 0.12e2 * t614 + 0.24e2 * t619 + 0.9e1 * t516 - 0.12e2 * t519 + 0.9e1 * t527 - t1267 - t1268 + t429 + t445;
1299 t1685 = -t696 - t703 - 0.12e2 * t625 + 0.3e1 * t629 - 0.6e1 * t632 + 0.3e1 * t636 + 0.12e2 * t638 - 0.12e2 * t642 - 0.6e1 * t646 - t1275 - t1278;
1300 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1301 v4rho4[3] = t1661 + t1668 + t1677 + t1685;
1302
1303 t1698 = -0.32e2 / 0.3e1 * t682 + 0.12e2 * t516 + 0.12e2 * t527 - 0.12e2 * t594 - 0.24e2 * t596 - 0.12e2 * t599 - 0.48e2 * t601 - 0.24e2 * t604 + 0.4e1 * t693 - 0.24e2 * t695 + 0.4e1 * t700;
1304 t1715 = t285 * t285;
1305 t1716 = 0.1e1 / t1715;
1306 t1724 = t294 * t294;
1307 t1726 = 0.1e1 / t41 / t1724;
1308 t1754 = 0.12e2 * t462 + 0.12e2 * t1017;
1309 t1771 = t258 * t258;
1310 t1775 = -0.8e1 / 0.3e1 * t32 * t12 * (0.79369597216623554946e0 * sigma[2] * t107 * t42 - 0.36662920248684660554e-1 * t346 * t110 + 0.62059993955520629303e-3 * t290 * t882 * t296 - 0.44595355499553507184e-5 * t671 * t1716 * t677 + 0.11598973244461995009e-7 * t671 * sigma[2] / t1715 / t34 * t1726) - 0.16e2 / 0.3e1 * t1511 + 0.6e1 * t539 * t627 * t263 - 0.8e1 * t215 * t216 * t707 + 0.6e1 * t550 * t634 * t275 - 0.24e2 * t234 * t640 * t275 - 0.8e1 * t234 * t235 * t720 - 0.12e2 * t1516 - 0.12e2 * t1521 + 0.4e1 * t1288 - 0.4e1 * t1088 + t931 + t933 - 0.4e1 * t51 * t26 * t1754 - 0.12e2 * t198 * t300 * t267 - 0.24e2 * t81 * t300 * t275 - 0.16e2 * t81 * t114 * t720 + 0.4e1 * t81 * t46 * t1754 - 0.3e1 / 0.2e1 * t1033 * t26 * t1771;
1311 t1776 = t263 * t263;
1312 t1780 = t267 * t267;
1313 t1784 = t275 * t275;
1314 t1808 = -0.6e1 * t149 * t26 * t1776 - 0.3e1 / 0.2e1 * t1054 * t46 * t1780 - 0.6e1 * t198 * t46 * t1784 - 0.16e2 * t81 * t681 * t96 + 0.4e1 * t431 * t114 * t698 + 0.2e1 * t1573 - 0.8e1 * t1576 - 0.12e2 * t1579 - 0.24e2 * t1582 - 0.24e2 * t1585 - 0.8e1 * t1588 + 0.3e1 * t1298 + 0.6e1 * t1301 + 0.3e1 * t1303 + 0.12e2 * t1308 + 0.6e1 * t1311 + 0.2e1 * t1629 - 0.6e1 * t1110 - 0.6e1 * t1122;
1315 t1811 = -0.24e2 * t702 - t926 - t929 + 0.8e1 * t496 - 0.16e2 * t580 - 0.16e2 * t709 - 0.24e2 * t712 - 0.48e2 * t715 - 0.48e2 * t718 - 0.16e2 * t722 + t4 * (t1775 + t1808);
1316 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1317 v4rho4[4] = t1698 + t1811;
1318
1319 t1813 = t51 * t733 * t56;
1320 t1816 = t51 * t305 * t163;
1321 t1819 = t11 * t167 * t128;
1322 t1820 = t148 * t322;
1323 t1824 = t9 * t760;
1324 t1835 = t11 * t60 * t322;
1325 t1838 = t11 * t12 * t760;
1326 t1861 = t149 * t305 * t150;
1327 t1864 = t51 * t730 * t56;
1328 t1866 = t415 * t128;
1329 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1330 v4rho3sigma[0] = -0.36e2 * t1813 - 0.18e2 * t1816 + t1819 - 0.6e1 * t120 * t1820 * t150 - 0.12e2 * t120 * t1824 * t56 - 0.12e2 * t120 * t740 * t163 - 0.4e1 * t120 * t308 * t465 - 0.6e1 * t1835 - 0.12e2 * t1838 - 0.8e1 / 0.3e1 * t120 * t10 * (-0.13228266202770592491e0 * t124 * t22 + 0.84225627598329625600e-2 * t356 * t74 * sigma[0] - 0.16898118509206682978e-3 * t398 * t190 * t69 + 0.13626358624863571640e-5 * t404 * t408 * t184 - 0.38663244148206650030e-8 * t402 / t948 / t13 * t959) - 0.9e1 * t1861 - 0.9e1 * t1864 + t120 * t1866 * t417 - 0.6e1 * t782 * t128 * t56 * t163;
1331
1332 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1333 v4rho3sigma[1] = 0.0e0;
1334
1335 t1874 = t198 * t326 * t199;
1336 t1877 = t81 * t767 * t82;
1337 t1880 = t81 * t326 * t206;
1338 t1883 = t32 * t167 * t139;
1339 t1884 = t430 * t139;
1340 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1341 v4rho3sigma[2] = -0.6e1 * t797 * t139 * t82 * t206 + t120 * t1884 * t432 - 0.4e1 * t120 * t328 * t473 - 0.9e1 * t1874 - 0.9e1 * t1877 - 0.18e2 * t1880 + t1883;
1342
1343 t1900 = t60 * t2;
1344 t1901 = t1900 * t332;
1345 t1904 = 0.6e1 * t215 * t784;
1346 t1916 = 0.12e2 * t779 * t787;
1347 t1922 = t120 * t415;
1348 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1349 v4rho3sigma[3] = -0.4e1 * t1184 * t128 * t3 * t56 - 0.4e1 * t782 * t322 * t92 * t56 - 0.4e1 * t120 * t1824 * t92 + t1922 * t783 * t150 - 0.2e1 * t782 * t783 * t163 - 0.8e1 * t790 * t740 * t3 - 0.4e1 * t790 * t308 - 0.12e2 * t1813 - 0.6e1 * t1816 + t1819 - 0.4e1 * t1835 - 0.4e1 * t1838 - 0.3e1 * t1861 - 0.6e1 * t1864 - 0.3e1 * t1901 - t1904 - t1916;
1350
1351 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1352 v4rho3sigma[4] = 0.0e0;
1353
1354 t1932 = t1900 * t336;
1355 t1935 = 0.6e1 * t234 * t799;
1356 t1936 = t120 * t430;
1357 t1948 = t1900 * t350;
1358 t1951 = 0.12e2 * t779 * t808;
1359 t1952 = t197 * t349;
1360 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1361 v4rho3sigma[5] = 0.4e1 * t1198 * t139 * t3 * t82 - 0.2e1 * t120 * t1952 * t199 - 0.4e1 * t120 * t807 * t206 + t1936 * t798 * t199 - 0.2e1 * t797 * t798 * t206 + 0.4e1 * t790 * t328 - 0.3e1 * t1874 - 0.6e1 * t1877 - 0.6e1 * t1880 + t1883 - 0.3e1 * t1932 - t1935 - 0.2e1 * t1948 - t1951;
1362
1363 t1960 = t227 * t2;
1364 t1965 = t779 * t812;
1365 t1976 = t779 * t815;
1366 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1367 v4rho3sigma[6] = t1922 * t128 * t258 * t56 - 0.2e1 * t782 * t128 * t263 * t56 - 0.4e1 * t1184 * t783 * t3 - 0.2e1 * t120 * t1820 * t258 - 0.4e1 * t120 * t740 * t263 - 0.4e1 * t120 * t308 * t612 - 0.12e2 * t1960 * t791 + t1819 - 0.2e1 * t1835 - 0.3e1 * t1864 - 0.6e1 * t1901 - t1904 - t1916 - 0.3e1 * t1965 - 0.6e1 * t1976;
1368
1369 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1370 v4rho3sigma[7] = 0.0e0;
1371
1372 t1993 = t779 * t820;
1373 t2001 = t779 * t823;
1374 t2010 = t779 * t826;
1375 t2019 = t779 * t843;
1376 t2021 = t30 * t842;
1377 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1378 v4rho3sigma[8] = t1936 * t139 * t267 * t82 - 0.2e1 * t797 * t139 * t275 * t82 - 0.4e1 * t797 * t349 * t96 * t82 + 0.4e1 * t1198 * t798 * t3 - 0.4e1 * t120 * t2021 * t82 - 0.4e1 * t120 * t328 * t585 + 0.8e1 * t790 * t807 * t3 + 0.12e2 * t1960 * t802 - 0.3e1 * t1877 + t1883 - 0.6e1 * t1932 - t1935 - 0.4e1 * t1948 - t1951 - 0.3e1 * t1993 - 0.12e2 * t2001 - 0.6e1 * t2010 - 0.4e1 * t2019;
1379
1380 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1381 v4rho3sigma[9] = t120 * t1866 * t691 - 0.4e1 * t120 * t308 * t707 - 0.6e1 * t782 * t783 * t263 + t1819 - 0.9e1 * t1901 - 0.9e1 * t1965 - 0.18e2 * t1976;
1382
1383 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1384 v4rho3sigma[10] = 0.0e0;
1385
1386 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1387 v4rho3sigma[11] = t1883 - 0.9e1 * t1932 - 0.9e1 * t1993 - 0.6e1 * t1948 - 0.12e2 * t2019 - 0.6e1 * t120 * t1952 * t267 - 0.12e2 * t120 * t2021 * t96 - 0.12e2 * t120 * t807 * t275 - 0.4e1 * t120 * t328 * t720 - 0.8e1 / 0.3e1 * t120 * t31 * (-0.13228266202770592491e0 * t135 * t42 + 0.84225627598329625600e-2 * t367 * t110 * sigma[2] - 0.16898118509206682978e-3 * t667 * t296 * t105 + 0.13626358624863571640e-5 * t673 * t677 * t290 - 0.38663244148206650030e-8 * t671 / t1715 / t33 * t1726) - 0.36e2 * t2001 - 0.18e2 * t2010 + t120 * t1884 * t698 - 0.6e1 * t797 * t798 * t275;
1388
1389 t2080 = 0.2e1 * t1900 * t361;
1390 t2081 = t779 * t849;
1391 t2083 = t779 * t862;
1392 t2085 = t148 * t360;
1393 t2089 = t9 * t861;
1394 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1395 v4rho2sigma2[0] = -t2080 - 0.12e2 * t2081 - 0.8e1 * t2083 - 0.2e1 * t120 * t2085 * t150 - 0.8e1 * t120 * t2089 * t56 - 0.4e1 * t120 * t848 * t163 - 0.8e1 / 0.3e1 * t120 * t10 * (-0.11560380258594262337e-2 * t749 + 0.37284104895758517604e-4 * t752 * sigma[0] - 0.38539196110725253122e-6 * t756 * t408 * t69 + 0.12887748049402216676e-8 * t184 / t948 / rho[0] * t959);
1396
1397 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1398 v4rho2sigma2[1] = 0.0e0;
1399
1400 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1401 v4rho2sigma2[2] = 0.0e0;
1402
1403 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1404 v4rho2sigma2[3] = 0.0e0;
1405
1406 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1407 v4rho2sigma2[4] = 0.0e0;
1408
1409 t2112 = 0.2e1 * t1900 * t372;
1410 t2113 = t779 * t867;
1411 t2115 = t197 * t371;
1412 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1413 v4rho2sigma2[5] = -0.2e1 * t120 * t2115 * t199 - 0.4e1 * t120 * t866 * t206 - t2112 - 0.12e2 * t2113;
1414
1415 t2124 = t779 * t870;
1416 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1417 v4rho2sigma2[6] = -0.2e1 * t782 * t360 * t92 * t56 - 0.4e1 * t120 * t2089 * t92 - 0.4e1 * t790 * t848 * t3 - t2080 - 0.6e1 * t2081 - 0.4e1 * t2083 - 0.6e1 * t2124;
1418
1419 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1420 v4rho2sigma2[7] = 0.0e0;
1421
1422 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1423 v4rho2sigma2[8] = 0.0e0;
1424
1425 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1426 v4rho2sigma2[9] = 0.0e0;
1427
1428 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1429 v4rho2sigma2[10] = 0.0e0;
1430
1431 t2137 = t779 * t874;
1432 t2146 = t779 * t887;
1433 t2148 = t30 * t886;
1434 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1435 v4rho2sigma2[11] = -0.2e1 * t797 * t371 * t96 * t82 - 0.4e1 * t120 * t2148 * t82 + 0.4e1 * t790 * t866 * t3 - t2112 - 0.6e1 * t2113 - 0.6e1 * t2137 - 0.4e1 * t2146;
1436
1437 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1438 v4rho2sigma2[12] = -0.2e1 * t120 * t2085 * t258 - 0.4e1 * t120 * t848 * t263 - t2080 - 0.12e2 * t2124;
1439
1440 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1441 v4rho2sigma2[13] = 0.0e0;
1442
1443 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1444 v4rho2sigma2[14] = 0.0e0;
1445
1446 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1447 v4rho2sigma2[15] = 0.0e0;
1448
1449 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1450 v4rho2sigma2[16] = 0.0e0;
1451
1452 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1453 v4rho2sigma2[17] = -t2112 - 0.12e2 * t2137 - 0.8e1 * t2146 - 0.2e1 * t120 * t2115 * t267 - 0.8e1 * t120 * t2148 * t96 - 0.4e1 * t120 * t866 * t275 - 0.8e1 / 0.3e1 * t120 * t31 * (-0.11560380258594262337e-2 * t831 + 0.37284104895758517604e-4 * t834 * sigma[2] - 0.38539196110725253122e-6 * t838 * t677 * t105 + 0.12887748049402216676e-8 * t290 / t1715 / rho[1] * t1726);
1454
1455 t2185 = t779 * t896;
1456 t2187 = t9 * t895;
1457 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1458 v4rhosigma3[0] = -0.4e1 * t2185 - 0.4e1 * t120 * t2187 * t56 - 0.8e1 / 0.3e1 * t120 * t10 * (-0.54122087751907525555e-5 * t853 + 0.96347990276813132804e-7 * t857 * t408 * sigma[0] - 0.42959160164674055589e-9 * t69 * t949 * t959);
1459
1460 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1461 v4rhosigma3[1] = 0.0e0;
1462
1463 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1464 v4rhosigma3[2] = 0.0e0;
1465
1466 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1467 v4rhosigma3[3] = 0.0e0;
1468
1469 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1470 v4rhosigma3[4] = 0.0e0;
1471
1472 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1473 v4rhosigma3[5] = 0.0e0;
1474
1475 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1476 v4rhosigma3[6] = 0.0e0;
1477
1478 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1479 v4rhosigma3[7] = 0.0e0;
1480
1481 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1482 v4rhosigma3[8] = 0.0e0;
1483
1484 t2202 = t779 * t905;
1485 t2203 = t30 * t904;
1486 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1487 v4rhosigma3[9] = -0.4e1 * t120 * t2203 * t82 - 0.4e1 * t2202;
1488
1489 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1490 v4rhosigma3[10] = -0.4e1 * t120 * t2187 * t92 - 0.4e1 * t2185;
1491
1492 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1493 v4rhosigma3[11] = 0.0e0;
1494
1495 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1496 v4rhosigma3[12] = 0.0e0;
1497
1498 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1499 v4rhosigma3[13] = 0.0e0;
1500
1501 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1502 v4rhosigma3[14] = 0.0e0;
1503
1504 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1505 v4rhosigma3[15] = 0.0e0;
1506
1507 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1508 v4rhosigma3[16] = 0.0e0;
1509
1510 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1511 v4rhosigma3[17] = 0.0e0;
1512
1513 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1514 v4rhosigma3[18] = 0.0e0;
1515
1516 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1517 v4rhosigma3[19] = -0.4e1 * t2202 - 0.4e1 * t120 * t2203 * t96 - 0.8e1 / 0.3e1 * t120 * t31 * (-0.54122087751907525555e-5 * t878 + 0.96347990276813132804e-7 * t882 * t677 * sigma[2] - 0.42959160164674055589e-9 * t105 * t1716 * t1726);
1518
1519 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1520 v4sigma4[0] = -0.8e1 / 0.3e1 * t120 * t10 * (-0.18351998147964406248e-7 * t398 * t408 + 0.14319720054891351863e-9 * sigma[0] * t404 * t959);
1521
1522 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1523 v4sigma4[1] = 0.0e0;
1524
1525 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1526 v4sigma4[2] = 0.0e0;
1527
1528 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1529 v4sigma4[3] = 0.0e0;
1530
1531 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1532 v4sigma4[4] = 0.0e0;
1533
1534 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1535 v4sigma4[5] = 0.0e0;
1536
1537 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1538 v4sigma4[6] = 0.0e0;
1539
1540 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1541 v4sigma4[7] = 0.0e0;
1542
1543 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1544 v4sigma4[8] = 0.0e0;
1545
1546 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1547 v4sigma4[9] = 0.0e0;
1548
1549 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1550 v4sigma4[10] = 0.0e0;
1551
1552 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1553 v4sigma4[11] = 0.0e0;
1554
1555 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1556 v4sigma4[12] = 0.0e0;
1557
1558 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1559 v4sigma4[13] = 0.0e0;
1560
1561 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1562 v4sigma4[14] = -0.8e1 / 0.3e1 * t120 * t31 * (-0.18351998147964406248e-7 * t667 * t677 + 0.14319720054891351863e-9 * sigma[2] * t673 * t1726);
1563
1564 #ifndef XC_DONT_COMPILE_MXC
1565
1566 if(order < 5) return;
1567
1568
1569 #endif
1570
1571 #endif
1572
1573 #endif
1574
1575 #endif
1576
1577 #endif
1578
1579
1580 }
1581
1582