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