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