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_pbepow.mpl
11   Type of functional: gga_exc
12 */
13 
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16 
17 
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,const double * sigma,double * zk,GGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
20 {
21 
22 #ifndef XC_DONT_COMPILE_EXC
23   double t1, t3, t5, t6, t7, t8, t9, t10;
24   double t11, t12, t13, t14, t15, t16, t17, t18;
25   double t19, t20, t22, t26, t27, t28, t29, t30;
26   double t31, t32, t34, t37, t39;
27 
28 #ifndef XC_DONT_COMPILE_VXC
29   double t44, t45, t46, t48, t49, t50, t53, t54;
30   double t55, t58, t61, t62, t63, t64, t65, t67;
31   double t69, t70, t71, t74, t77, t86, t90, t93;
32   double t97;
33 
34 #ifndef XC_DONT_COMPILE_FXC
35   double t100, t101, t109, t110, t111, t114, t118, t119;
36   double t120, t123, t126, t128, t130, t133, t134, t135;
37   double t136, t139, t142, t145, t150, t155, t156, t160;
38   double t164, t169, t174, t175, t179, t183, t189, t194;
39   double t198, t202;
40 
41 #ifndef XC_DONT_COMPILE_KXC
42   double t206, t207, t219, t220, t226, t230, t234, t235;
43   double t236, t239, t242, t246, t251, t252, t256, t257;
44   double t259, t261, t262, t266, t269, t274, t281, t285;
45   double t295, t296, t297, t300, t304, t307, t310, t319;
46   double t322, t324, t329, t332, t333, t347, t351, t354;
47   double t357, t361, t366, t370, t375, t379, t388, t391;
48   double t397, t402, t406;
49 
50 #ifndef XC_DONT_COMPILE_LXC
51   double t422, t446, t447, t448, t454, t473, t481, t487;
52   double t494, t526, t579, t589, t593, t635, t654, t684;
53   double t754, t767, t773;
54 #endif
55 
56 #endif
57 
58 #endif
59 
60 #endif
61 
62 #endif
63 
64 
65 
66   t1 = M_CBRT3;
67   t3 = POW_1_3(0.1e1 / M_PI);
68   t5 = M_CBRT4;
69   t6 = t5 * t5;
70   t7 = t1 * t3 * t6;
71   t8 = M_CBRT2;
72   t9 = t8 * t8;
73   t10 = POW_1_3(rho[0]);
74   t11 = t9 * t10;
75   t12 = M_CBRT6;
76   t13 = M_PI * M_PI;
77   t14 = POW_1_3(t13);
78   t15 = t14 * t14;
79   t16 = 0.1e1 / t15;
80   t17 = t12 * t16;
81   t18 = sigma[0] * t9;
82   t19 = rho[0] * rho[0];
83   t20 = t10 * t10;
84   t22 = 0.1e1 / t20 / t19;
85   t26 = 0.10000000000000000000e-1 + 0.11376190545424808251e-1 * t17 * t18 * t22;
86   t27 = t17 * sigma[0];
87   t28 = t9 * t22;
88   t29 = 0.1e1 / t26;
89   t30 = t28 * t29;
90   t31 = t27 * t30;
91   t32 = pow(t31, 0.1010e3);
92   t34 = 0.1e1 - 0.45256672948257153200e-196 * t32;
93   t37 = 0.99196000000000000000e0 + 0.8040e0 * t26 * t34;
94   t39 = t7 * t11 * t37;
95   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
96     zk[0] = -0.3e1 / 0.16e2 * t39;
97 
98 #ifndef XC_DONT_COMPILE_VXC
99 
100   if(order < 1) return;
101 
102 
103   t44 = t10 * rho[0] * t1 * t3;
104   t45 = t6 * t9;
105   t46 = t19 * rho[0];
106   t48 = 0.1e1 / t20 / t46;
107   t49 = t9 * t48;
108   t50 = t49 * t34;
109   t53 = pow(t31, 0.1000e3);
110   t54 = t26 * t53;
111   t55 = t49 * t29;
112   t58 = t12 * t12;
113   t61 = t58 / t14 / t13;
114   t62 = sigma[0] * sigma[0];
115   t63 = t61 * t62;
116   t64 = t19 * t19;
117   t65 = t64 * t19;
118   t67 = 0.1e1 / t10 / t65;
119   t69 = t26 * t26;
120   t70 = 0.1e1 / t69;
121   t71 = t8 * t67 * t70;
122   t74 = -0.8e1 / 0.3e1 * t27 * t55 + 0.60673016242265644006e-1 * t63 * t71;
123   t77 = -0.24390552529390788890e-1 * t27 * t50 - 0.36750228700902738685e-194 * t54 * t74;
124   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
125     vrho[0] = -t39 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t45 * t77;
126 
127   t86 = t64 * rho[0];
128   t90 = t8 / t10 / t86 * t70;
129   t93 = t17 * t30 - 0.22752381090849616502e-1 * t61 * sigma[0] * t90;
130   t97 = t45 * (0.91464571985215458338e-2 * t17 * t28 * t34 - 0.36750228700902738685e-194 * t54 * t93);
131   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
132     vsigma[0] = -0.3e1 / 0.16e2 * t44 * t97;
133 
134 #ifndef XC_DONT_COMPILE_FXC
135 
136   if(order < 2) return;
137 
138 
139   t100 = 0.1e1 / t20;
140   t101 = t9 * t100;
141   t109 = 0.1e1 / t20 / t64;
142   t110 = t9 * t109;
143   t111 = t110 * t34;
144   t114 = t53 * t74;
145   t118 = pow(t31, 0.990e2);
146   t119 = t26 * t118;
147   t120 = t74 * t74;
148   t123 = t110 * t29;
149   t126 = t64 * t46;
150   t128 = 0.1e1 / t10 / t126;
151   t130 = t8 * t128 * t70;
152   t133 = t62 * sigma[0];
153   t134 = t64 * t64;
154   t135 = t134 * t19;
155   t136 = 0.1e1 / t135;
156   t139 = 0.1e1 / t69 / t26;
157   t142 = 0.88e2 / 0.9e1 * t27 * t123 - 0.54605714618039079605e0 * t63 * t130 + 0.45349544206085259001e-3 * t133 * t136 * t139;
158   t145 = 0.89432025941099559263e-1 * t27 * t111 + 0.22297472228768489022e-195 * t27 * t49 * t114 - 0.36750228700902738685e-192 * t119 * t120 - 0.36750228700902738685e-194 * t54 * t142;
159   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
160     v2rho2[0] = -t7 * t101 * t37 / 0.12e2 - t7 * t11 * t77 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t45 * t145;
161 
162   t150 = t10 * t1 * t3;
163   t155 = t17 * t9;
164   t156 = t22 * t53;
165   t160 = t53 * t93;
166   t164 = t93 * t74;
167   t169 = t61 * t8;
168   t174 = t134 * rho[0];
169   t175 = 0.1e1 / t174;
170   t179 = -0.8e1 / 0.3e1 * t17 * t55 + 0.18201904872679693202e0 * t169 * t67 * t70 * sigma[0] - 0.17006079077281972125e-3 * t62 * t175 * t139;
171   t183 = t45 * (-0.24390552529390788890e-1 * t17 * t50 - 0.41807760428940916917e-196 * t155 * t156 * t74 + 0.11148736114384244511e-195 * t27 * t49 * t160 - 0.36750228700902738685e-192 * t119 * t164 - 0.36750228700902738685e-194 * t54 * t179);
172   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
173     v2rhosigma[0] = -t150 * t97 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t183;
174 
175   t189 = t93 * t93;
176   t194 = 0.1e1 / t134;
177   t198 = -0.45504762181699233004e-1 * t61 * t90 + 0.63772796539807395468e-4 * sigma[0] * t194 * t139;
178   t202 = t45 * (-0.83615520857881833834e-196 * t155 * t156 * t93 - 0.36750228700902738685e-192 * t119 * t189 - 0.36750228700902738685e-194 * t54 * t198);
179   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
180     v2sigma2[0] = -0.3e1 / 0.16e2 * t44 * t202;
181 
182 #ifndef XC_DONT_COMPILE_KXC
183 
184   if(order < 3) return;
185 
186 
187   t206 = 0.1e1 / t20 / rho[0];
188   t207 = t9 * t206;
189   t219 = t9 / t20 / t86;
190   t220 = t219 * t34;
191   t226 = t118 * t120;
192   t230 = t53 * t142;
193   t234 = pow(t31, 0.980e2);
194   t235 = t26 * t234;
195   t236 = t120 * t74;
196   t239 = t74 * t142;
197   t242 = t219 * t29;
198   t246 = 0.1e1 / t10 / t134;
199   t251 = t134 * t46;
200   t252 = 0.1e1 / t251;
201   t256 = t62 * t62;
202   t257 = t134 * t86;
203   t259 = 0.1e1 / t20 / t257;
204   t261 = t69 * t69;
205   t262 = 0.1e1 / t261;
206   t266 = -0.1232e4 / 0.27e2 * t27 * t242 + 0.45976663419139076902e1 * t63 * t8 * t246 * t70 - 0.86164133991561992102e-2 * t133 * t252 * t139 + 0.41272404482927321238e-4 * t256 * t259 * t262 * t155;
207   t269 = -0.41734945439179794323e0 * t27 * t220 - 0.12263609725822668962e-194 * t27 * t110 * t114 + 0.33446208343152733533e-193 * t27 * t49 * t226 + 0.33446208343152733533e-195 * t27 * t49 * t230 - 0.36382726413893711298e-190 * t235 * t236 - 0.11025068610270821606e-191 * t119 * t239 - 0.36750228700902738685e-194 * t54 * t266;
208   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
209     v3rho3[0] = t7 * t207 * t37 / 0.18e2 - t7 * t101 * t77 / 0.4e1 - 0.3e1 / 0.4e1 * t7 * t11 * t145 - 0.3e1 / 0.16e2 * t44 * t45 * t269;
210 
211   t274 = t100 * t1 * t3;
212   t281 = t48 * t53;
213   t285 = t22 * t118;
214   t295 = t17 * t18;
215   t296 = t48 * t118;
216   t297 = t296 * t164;
217   t300 = t53 * t179;
218   t304 = t93 * t120;
219   t307 = t179 * t74;
220   t310 = t93 * t142;
221   t319 = t136 * t139;
222   t322 = t134 * t64;
223   t324 = 0.1e1 / t20 / t322;
224   t329 = 0.88e2 / 0.9e1 * t17 * t123 - 0.13145820185824222868e1 * t169 * t128 * t70 * sigma[0] + 0.28910334431379352613e-2 * t319 * t62 - 0.15477151681097745464e-4 * t133 * t324 * t262 * t155;
225   t332 = 0.89432025941099559263e-1 * t17 * t111 + 0.22297472228768489022e-195 * t155 * t281 * t74 - 0.41807760428940916917e-194 * t155 * t285 * t120 - 0.41807760428940916917e-196 * t155 * t156 * t142 - 0.40878699086075563207e-195 * t27 * t110 * t160 + 0.22297472228768489022e-193 * t295 * t297 + 0.22297472228768489022e-195 * t27 * t49 * t300 - 0.36382726413893711298e-190 * t235 * t304 - 0.73500457401805477370e-192 * t119 * t307 - 0.36750228700902738685e-192 * t119 * t310 - 0.36750228700902738685e-194 * t54 * t329;
226   t333 = t45 * t332;
227   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
228     v3rho2sigma[0] = -t274 * t97 / 0.12e2 - t150 * t183 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t333;
229 
230   t347 = t118 * t189;
231   t351 = t189 * t74;
232   t354 = t93 * t179;
233   t357 = t53 * t198;
234   t361 = t198 * t74;
235   t366 = t175 * t139;
236   t370 = 0.1e1 / t20 / t251;
237   t375 = 0.24269206496906257602e0 * t61 * t71 - 0.85030395386409860625e-3 * t366 * sigma[0] + 0.58039318804116545489e-5 * t62 * t370 * t262 * t155;
238   t379 = t45 * (0.22297472228768489022e-195 * t155 * t281 * t93 - 0.83615520857881833834e-194 * t155 * t285 * t164 - 0.83615520857881833834e-196 * t155 * t156 * t179 + 0.11148736114384244511e-193 * t27 * t49 * t347 - 0.36382726413893711298e-190 * t235 * t351 - 0.73500457401805477370e-192 * t119 * t354 + 0.11148736114384244511e-195 * t27 * t49 * t357 - 0.36750228700902738685e-192 * t119 * t361 - 0.36750228700902738685e-194 * t54 * t375);
239   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
240     v3rhosigma2[0] = -t150 * t202 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t379;
241 
242   t388 = t189 * t93;
243   t391 = t93 * t198;
244   t397 = 0.1e1 / t20 / t135;
245   t402 = 0.19131838961942218641e-3 * t194 * t139 - 0.21764744551543704558e-5 * sigma[0] * t397 * t262 * t155;
246   t406 = t45 * (-0.12542328128682275075e-193 * t155 * t285 * t189 - 0.12542328128682275075e-195 * t155 * t156 * t198 - 0.36382726413893711298e-190 * t235 * t388 - 0.11025068610270821606e-191 * t119 * t391 - 0.36750228700902738685e-194 * t54 * t402);
247   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
248     v3sigma3[0] = -0.3e1 / 0.16e2 * t44 * t406;
249 
250 #ifndef XC_DONT_COMPILE_LXC
251 
252   if(order < 4) return;
253 
254 
255   t422 = t9 / t20 / t65;
256   t446 = pow(t31, 0.970e2);
257   t447 = t26 * t446;
258   t448 = t120 * t120;
259   t454 = t142 * t142;
260   t473 = t134 * t65;
261   t481 = t134 * t134;
262   t487 = 0.1e1 / t261 / t26;
263   t494 = 0.23649802415535216783e1 * t27 * t422 * t34 + 0.76306904960674384652e-194 * t27 * t219 * t114 - 0.24527219451645337924e-192 * t27 * t110 * t226 - 0.24527219451645337924e-194 * t27 * t110 * t230 + 0.44148995012961608264e-191 * t27 * t49 * t234 * t236 + 0.13378483337261093413e-192 * t295 * t296 * t239 + 0.44594944457536978044e-195 * t27 * t49 * t53 * t266 - 0.35655071885615837072e-188 * t447 * t448 - 0.21829635848336226779e-189 * t235 * t120 * t142 - 0.11025068610270821606e-191 * t119 * t454 - 0.14700091480361095474e-191 * t119 * t74 * t266 - 0.36750228700902738685e-194 * t54 * (0.20944e5 / 0.81e2 * t27 * t422 * t29 - 0.41082373442262981619e2 * t63 * t8 / t10 / t174 * t70 + 0.12914542422244057647e0 * t133 / t322 * t139 - 0.13482318797756258271e-2 * t256 / t20 / t473 * t262 * t155 + 0.10016485070200026960e-4 * t256 * sigma[0] / t10 / t481 / rho[0] * t487 * t169);
264   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
265     v4rho4[0] = -0.5e1 / 0.54e2 * t7 * t28 * t37 + 0.2e1 / 0.9e1 * t7 * t207 * t77 - t7 * t101 * t145 / 0.2e1 - t7 * t11 * t269 - 0.3e1 / 0.16e2 * t44 * t45 * t494;
266 
267   t526 = t22 * t234;
268   t579 = t109 * t53;
269   t589 = t48 * t234;
270   t593 = -0.35655071885615837072e-188 * t447 * t93 * t236 - 0.10914817924168113389e-189 * t235 * t164 * t142 - 0.10914817924168113389e-189 * t235 * t179 * t120 - 0.11025068610270821606e-191 * t119 * t329 * t74 - 0.11025068610270821606e-191 * t119 * t179 * t142 - 0.36750228700902738685e-192 * t119 * t93 * t266 - 0.41734945439179794323e0 * t17 * t220 - 0.12263609725822668962e-194 * t155 * t579 * t74 + 0.33446208343152733533e-195 * t155 * t281 * t142 - 0.41807760428940916917e-196 * t155 * t156 * t266 + 0.33111746259721206198e-191 * t295 * t589 * t304;
271   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
272     v4rho3sigma[0] = t206 * t1 * t3 * t97 / 0.18e2 - t274 * t183 / 0.4e1 - 0.3e1 / 0.4e1 * t150 * t333 - 0.3e1 / 0.16e2 * t44 * t45 * (0.66892416686305467067e-193 * t295 * t296 * t307 + 0.33446208343152733533e-193 * t295 * t296 * t310 + 0.19076726240168596163e-194 * t27 * t219 * t160 - 0.12263609725822668962e-192 * t295 * t109 * t118 * t164 + 0.33446208343152733533e-195 * t27 * t49 * t53 * t329 - 0.12263609725822668962e-194 * t27 * t110 * t300 - 0.41389682824651507748e-192 * t155 * t526 * t236 - 0.12542328128682275075e-193 * t155 * t285 * t239 + 0.33446208343152733533e-193 * t155 * t296 * t120 - 0.36750228700902738685e-194 * t54 * (-0.1232e4 / 0.27e2 * t17 * t242 + 0.10233515406195471956e2 * t169 * t246 * t70 * sigma[0] - 0.38736069009364492063e-1 * t252 * t139 * t62 + 0.45915549987256644877e-3 * t259 * t262 * t133 * t155 - 0.37561819013250101099e-5 * t256 / t10 / t481 * t487 * t169) + t593);
273 
274   t635 = 0.22297472228768489022e-195 * t27 * t49 * t53 * t375 - 0.16723104171576366767e-193 * t155 * t285 * t307 - 0.83615520857881833834e-194 * t155 * t285 * t310 - 0.81757398172151126414e-195 * t155 * t579 * t93 + 0.44594944457536978044e-193 * t155 * t297 - 0.83615520857881833834e-196 * t155 * t156 * t329 + 0.44594944457536978044e-195 * t155 * t281 * t179 - 0.36382726413893711298e-190 * t235 * t198 * t120 - 0.35655071885615837072e-188 * t447 * t189 * t120 - 0.14553090565557484519e-189 * t235 * t164 * t179 - 0.40878699086075563207e-195 * t27 * t110 * t357;
275   t654 = t179 * t179;
276   t684 = 0.22297472228768489022e-193 * t295 * t296 * t361 + 0.44594944457536978045e-193 * t295 * t296 * t354 - 0.40878699086075563207e-193 * t27 * t110 * t347 + 0.22074497506480804132e-191 * t295 * t589 * t351 - 0.82779365649303015496e-192 * t155 * t526 * t304 - 0.36382726413893711298e-190 * t235 * t189 * t142 - 0.73500457401805477370e-192 * t119 * t654 - 0.73500457401805477370e-192 * t119 * t93 * t329 - 0.73500457401805477370e-192 * t119 * t375 * t74 - 0.36750228700902738685e-192 * t119 * t198 * t142 - 0.36750228700902738685e-194 * t54 * (-0.15370497448040629815e1 * t61 * t130 + 0.94667173530202978162e-2 * t319 * sigma[0] - 0.14509829701029136372e-3 * t324 * t262 * t62 * t155 + 0.14085682129968787912e-5 * t133 / t10 / t134 / t126 * t487 * t169);
277   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
278     v4rho2sigma2[0] = -t274 * t202 / 0.12e2 - t150 * t379 / 0.2e1 - 0.3e1 / 0.16e2 * t44 * t45 * (t635 + t684);
279 
280   t754 = 0.33446208343152733533e-193 * t155 * t296 * t189 - 0.12416904847395452324e-191 * t155 * t526 * t351 - 0.25084656257364550150e-193 * t155 * t285 * t354 + 0.33446208343152733533e-195 * t155 * t281 * t198 - 0.12542328128682275075e-193 * t155 * t285 * t361 - 0.12542328128682275075e-195 * t155 * t156 * t375 + 0.11037248753240402066e-191 * t27 * t49 * t234 * t388 - 0.35655071885615837072e-188 * t447 * t388 * t74 - 0.10914817924168113389e-189 * t235 * t189 * t179 + 0.33446208343152733536e-193 * t295 * t296 * t391 - 0.10914817924168113390e-189 * t235 * t391 * t74 - 0.11025068610270821606e-191 * t119 * t179 * t198 - 0.11025068610270821606e-191 * t119 * t93 * t375 + 0.11148736114384244511e-195 * t27 * t49 * t53 * t402 - 0.36750228700902738685e-192 * t119 * t402 * t74 - 0.36750228700902738685e-194 * t54 * (-0.15305471169553774913e-2 * t366 + 0.40627523162881581842e-4 * t370 * t262 * t12 * t16 * sigma[0] * t9 - 0.52821307987382954669e-6 * t62 / t10 / t473 * t487 * t169);
281   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
282     v4rhosigma3[0] = -t150 * t406 / 0.4e1 - 0.3e1 / 0.16e2 * t44 * t45 * t754;
283 
284   t767 = t189 * t189;
285   t773 = t198 * t198;
286   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
287     v4sigma4[0] = -0.3e1 / 0.16e2 * t44 * t45 * (-0.16555873129860603099e-191 * t155 * t526 * t388 - 0.50169312514729100301e-193 * t155 * t285 * t391 - 0.16723104171576366767e-195 * t155 * t156 * t402 - 0.35655071885615837072e-188 * t447 * t767 - 0.21829635848336226779e-189 * t235 * t189 * t198 - 0.11025068610270821606e-191 * t119 * t773 - 0.14700091480361095474e-191 * t119 * t93 * t402 - 0.36750228700902738685e-194 * t54 * (-0.87058978206174818234e-5 * t397 * t262 * t155 + 0.19807990495268608001e-6 * sigma[0] / t10 / t257 * t487 * t169));
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, t18;
315   double t19, t22, t23, t25, t26, t28, t31, t33;
316 
317 #ifndef XC_DONT_COMPILE_VXC
318   double t37, t38, t39, t41, t42, t46, t47, t51;
319   double t53, t54, t55, t56, t57, t59, t61, t62;
320   double t66, t69, t78, t80, t85, t89;
321 
322 #ifndef XC_DONT_COMPILE_FXC
323   double t92, t93, t100, t101, t105, t106, t107, t110;
324   double t111, t112, t118, t120, t125, t126, t127, t128;
325   double t131, t134, t137, t141, t147, t151, t154, t160;
326   double t164, t165, t169, t173, t179, t185, t189, t193;
327 
328 #ifndef XC_DONT_COMPILE_KXC
329   double t197, t198, t209, t210, t214, t215, t218, t219;
330   double t222, t225, t226, t227, t230, t237, t242, t243;
331   double t247, t248, t250, t252, t253, t255, t258, t261;
332   double t265, t274, t281, t287, t290, t293, t296, t302;
333   double t306, t309, t311, t315, t318, t319, t326, t327;
334   double t328, t334, t337, t340, t343, t346, t351, t355;
335   double t359, t363, t372, t375, t381, t385, t389;
336 
337 #ifndef XC_DONT_COMPILE_LXC
338   double t406, t410, t414, t421, t431, t432, t433, t439;
339   double t458, t465, t473, t479, t559, t573, t621, t630;
340   double t663, t727, t740, t746;
341 #endif
342 
343 #endif
344 
345 #endif
346 
347 #endif
348 
349 #endif
350 
351 
352 
353   t1 = M_CBRT3;
354   t3 = POW_1_3(0.1e1 / M_PI);
355   t4 = t1 * t3;
356   t5 = M_CBRT4;
357   t6 = t5 * t5;
358   t7 = POW_1_3(rho[0]);
359   t8 = t6 * t7;
360   t9 = M_CBRT6;
361   t10 = M_PI * M_PI;
362   t11 = POW_1_3(t10);
363   t12 = t11 * t11;
364   t13 = 0.1e1 / t12;
365   t14 = t9 * t13;
366   t15 = rho[0] * rho[0];
367   t16 = t7 * t7;
368   t18 = 0.1e1 / t16 / t15;
369   t19 = sigma[0] * t18;
370   t22 = 0.10000000000000000000e-1 + 0.11376190545424808251e-1 * t14 * t19;
371   t23 = 0.1e1 / t22;
372   t25 = t14 * t19 * t23;
373   t26 = pow(t25, 0.1010e3);
374   t28 = 0.1e1 - 0.45256672948257153200e-196 * t26;
375   t31 = 0.99196000000000000000e0 + 0.8040e0 * t22 * t28;
376   t33 = t4 * t8 * t31;
377   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
378     zk[0] = -0.3e1 / 0.8e1 * t33;
379 
380 #ifndef XC_DONT_COMPILE_VXC
381 
382   if(order < 1) return;
383 
384 
385   t37 = t7 * rho[0] * t1;
386   t38 = t3 * t6;
387   t39 = t15 * rho[0];
388   t41 = 0.1e1 / t16 / t39;
389   t42 = sigma[0] * t41;
390   t46 = pow(t25, 0.1000e3);
391   t47 = t22 * t46;
392   t51 = t9 * t9;
393   t53 = 0.1e1 / t11 / t10;
394   t54 = t51 * t53;
395   t55 = sigma[0] * sigma[0];
396   t56 = t15 * t15;
397   t57 = t56 * t15;
398   t59 = 0.1e1 / t7 / t57;
399   t61 = t22 * t22;
400   t62 = 0.1e1 / t61;
401   t66 = -0.8e1 / 0.3e1 * t14 * t42 * t23 + 0.30336508121132822003e-1 * t54 * t55 * t59 * t62;
402   t69 = -0.24390552529390788890e-1 * t14 * t42 * t28 - 0.36750228700902738685e-194 * t47 * t66;
403   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
404     vrho[0] = -t33 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * t69;
405 
406   t78 = t56 * rho[0];
407   t80 = 0.1e1 / t7 / t78;
408   t85 = t14 * t18 * t23 - 0.11376190545424808251e-1 * t54 * sigma[0] * t80 * t62;
409   t89 = t38 * (0.91464571985215458338e-2 * t14 * t18 * t28 - 0.36750228700902738685e-194 * t47 * t85);
410   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
411     vsigma[0] = -0.3e1 / 0.8e1 * t37 * t89;
412 
413 #ifndef XC_DONT_COMPILE_FXC
414 
415   if(order < 2) return;
416 
417 
418   t92 = 0.1e1 / t16;
419   t93 = t6 * t92;
420   t100 = 0.1e1 / t16 / t56;
421   t101 = sigma[0] * t100;
422   t105 = t14 * sigma[0];
423   t106 = t41 * t46;
424   t107 = t106 * t66;
425   t110 = pow(t25, 0.990e2);
426   t111 = t22 * t110;
427   t112 = t66 * t66;
428   t118 = t56 * t39;
429   t120 = 0.1e1 / t7 / t118;
430   t125 = t55 * sigma[0];
431   t126 = t56 * t56;
432   t127 = t126 * t15;
433   t128 = 0.1e1 / t127;
434   t131 = 0.1e1 / t61 / t22;
435   t134 = 0.88e2 / 0.9e1 * t14 * t101 * t23 - 0.27302857309019539803e0 * t54 * t55 * t120 * t62 + 0.11337386051521314750e-3 * t125 * t128 * t131;
436   t137 = 0.89432025941099559263e-1 * t14 * t101 * t28 + 0.22297472228768489022e-195 * t105 * t107 - 0.36750228700902738685e-192 * t111 * t112 - 0.36750228700902738685e-194 * t47 * t134;
437   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
438     v2rho2[0] = -t4 * t93 * t31 / 0.6e1 - t4 * t8 * t69 - 0.3e1 / 0.8e1 * t37 * t38 * t137;
439 
440   t141 = t7 * t1;
441   t147 = t18 * t46;
442   t151 = t106 * t85;
443   t154 = t85 * t66;
444   t160 = t59 * t62;
445   t164 = t126 * rho[0];
446   t165 = 0.1e1 / t164;
447   t169 = -0.8e1 / 0.3e1 * t14 * t41 * t23 + 0.91009524363398466008e-1 * t54 * t160 * sigma[0] - 0.42515197693204930312e-4 * t55 * t165 * t131;
448   t173 = t38 * (-0.24390552529390788890e-1 * t14 * t41 * t28 - 0.41807760428940916917e-196 * t14 * t147 * t66 + 0.11148736114384244511e-195 * t105 * t151 - 0.36750228700902738685e-192 * t111 * t154 - 0.36750228700902738685e-194 * t47 * t169);
449   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
450     v2rhosigma[0] = -t141 * t89 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t173;
451 
452   t179 = t85 * t85;
453   t185 = 0.1e1 / t126;
454   t189 = -0.22752381090849616502e-1 * t54 * t80 * t62 + 0.15943199134951848867e-4 * sigma[0] * t185 * t131;
455   t193 = t38 * (-0.83615520857881833834e-196 * t14 * t147 * t85 - 0.36750228700902738685e-192 * t111 * t179 - 0.36750228700902738685e-194 * t47 * t189);
456   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
457     v2sigma2[0] = -0.3e1 / 0.8e1 * t37 * t193;
458 
459 #ifndef XC_DONT_COMPILE_KXC
460 
461   if(order < 3) return;
462 
463 
464   t197 = 0.1e1 / t16 / rho[0];
465   t198 = t6 * t197;
466   t209 = 0.1e1 / t16 / t78;
467   t210 = sigma[0] * t209;
468   t214 = t100 * t46;
469   t215 = t214 * t66;
470   t218 = t41 * t110;
471   t219 = t218 * t112;
472   t222 = t106 * t134;
473   t225 = pow(t25, 0.980e2);
474   t226 = t22 * t225;
475   t227 = t112 * t66;
476   t230 = t66 * t134;
477   t237 = 0.1e1 / t7 / t126;
478   t242 = t126 * t39;
479   t243 = 0.1e1 / t242;
480   t247 = t55 * t55;
481   t248 = t126 * t78;
482   t250 = 0.1e1 / t16 / t248;
483   t252 = t61 * t61;
484   t253 = 0.1e1 / t252;
485   t255 = t253 * t9 * t13;
486   t258 = -0.1232e4 / 0.27e2 * t14 * t210 * t23 + 0.22988331709569538451e1 * t54 * t55 * t237 * t62 - 0.21541033497890498025e-2 * t125 * t243 * t131 + 0.10318101120731830309e-4 * t247 * t250 * t255;
487   t261 = -0.41734945439179794323e0 * t14 * t210 * t28 - 0.12263609725822668962e-194 * t105 * t215 + 0.33446208343152733533e-193 * t105 * t219 + 0.33446208343152733533e-195 * t105 * t222 - 0.36382726413893711298e-190 * t226 * t227 - 0.11025068610270821606e-191 * t111 * t230 - 0.36750228700902738685e-194 * t47 * t258;
488   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
489     v3rho3[0] = t4 * t198 * t31 / 0.9e1 - t4 * t93 * t69 / 0.2e1 - 0.3e1 / 0.2e1 * t4 * t8 * t137 - 0.3e1 / 0.8e1 * t37 * t38 * t261;
490 
491   t265 = t92 * t1;
492   t274 = t18 * t110;
493   t281 = t214 * t85;
494   t287 = t106 * t169;
495   t290 = t85 * t112;
496   t293 = t169 * t66;
497   t296 = t85 * t134;
498   t302 = t120 * t62;
499   t306 = t128 * t131;
500   t309 = t126 * t56;
501   t311 = 0.1e1 / t16 / t309;
502   t315 = 0.88e2 / 0.9e1 * t14 * t100 * t23 - 0.65729100929121114339e0 * t54 * t302 * sigma[0] + 0.72275836078448381531e-3 * t306 * t55 - 0.38692879202744363659e-5 * t125 * t311 * t255;
503   t318 = 0.89432025941099559263e-1 * t14 * t100 * t28 + 0.22297472228768489022e-195 * t14 * t107 - 0.41807760428940916917e-194 * t14 * t274 * t112 - 0.41807760428940916917e-196 * t14 * t147 * t134 - 0.40878699086075563207e-195 * t105 * t281 + 0.22297472228768489022e-193 * t105 * t218 * t154 + 0.22297472228768489022e-195 * t105 * t287 - 0.36382726413893711298e-190 * t226 * t290 - 0.73500457401805477370e-192 * t111 * t293 - 0.36750228700902738685e-192 * t111 * t296 - 0.36750228700902738685e-194 * t47 * t315;
504   t319 = t38 * t318;
505   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
506     v3rho2sigma[0] = -t265 * t89 / 0.6e1 - t141 * t173 - 0.3e1 / 0.8e1 * t37 * t319;
507 
508   t326 = t14 * t18;
509   t327 = t110 * t85;
510   t328 = t327 * t66;
511   t334 = t218 * t179;
512   t337 = t179 * t66;
513   t340 = t85 * t169;
514   t343 = t106 * t189;
515   t346 = t189 * t66;
516   t351 = t165 * t131;
517   t355 = 0.1e1 / t16 / t242;
518   t359 = 0.12134603248453128801e0 * t54 * t160 - 0.21257598846602465156e-3 * t351 * sigma[0] + 0.14509829701029136372e-5 * t55 * t355 * t255;
519   t363 = t38 * (0.22297472228768489022e-195 * t14 * t151 - 0.83615520857881833834e-194 * t326 * t328 - 0.83615520857881833834e-196 * t14 * t147 * t169 + 0.11148736114384244511e-193 * t105 * t334 - 0.36382726413893711298e-190 * t226 * t337 - 0.73500457401805477370e-192 * t111 * t340 + 0.11148736114384244511e-195 * t105 * t343 - 0.36750228700902738685e-192 * t111 * t346 - 0.36750228700902738685e-194 * t47 * t359);
520   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
521     v3rhosigma2[0] = -t141 * t193 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t363;
522 
523   t372 = t179 * t85;
524   t375 = t85 * t189;
525   t381 = 0.1e1 / t16 / t127;
526   t385 = 0.47829597404855546601e-4 * t185 * t131 - 0.54411861378859261395e-6 * sigma[0] * t381 * t255;
527   t389 = t38 * (-0.12542328128682275075e-193 * t14 * t274 * t179 - 0.12542328128682275075e-195 * t14 * t147 * t189 - 0.36382726413893711298e-190 * t226 * t372 - 0.11025068610270821606e-191 * t111 * t375 - 0.36750228700902738685e-194 * t47 * t385);
528   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
529     v3sigma3[0] = -0.3e1 / 0.8e1 * t37 * t389;
530 
531 #ifndef XC_DONT_COMPILE_LXC
532 
533   if(order < 4) return;
534 
535 
536   t406 = sigma[0] / t16 / t57;
537   t410 = t209 * t46;
538   t414 = t100 * t110;
539   t421 = t41 * t225;
540   t431 = pow(t25, 0.970e2);
541   t432 = t22 * t431;
542   t433 = t112 * t112;
543   t439 = t134 * t134;
544   t458 = t126 * t57;
545   t465 = t126 * t126;
546   t473 = 0.1e1 / t252 / t22 * t51 * t53;
547   t479 = 0.23649802415535216783e1 * t14 * t406 * t28 + 0.76306904960674384652e-194 * t105 * t410 * t66 - 0.24527219451645337924e-192 * t105 * t414 * t112 - 0.24527219451645337924e-194 * t105 * t214 * t134 + 0.44148995012961608264e-191 * t105 * t421 * t227 + 0.13378483337261093413e-192 * t105 * t218 * t230 + 0.44594944457536978044e-195 * t105 * t106 * t258 - 0.35655071885615837072e-188 * t432 * t433 - 0.21829635848336226779e-189 * t226 * t112 * t134 - 0.11025068610270821606e-191 * t111 * t439 - 0.14700091480361095474e-191 * t111 * t66 * t258 - 0.36750228700902738685e-194 * t47 * (0.20944e5 / 0.81e2 * t14 * t406 * t23 - 0.20541186721131490809e2 * t54 * t55 / t7 / t164 * t62 + 0.32286356055610144116e-1 * t125 / t309 * t131 - 0.33705796994390645677e-3 * t247 / t16 / t458 * t255 + 0.12520606337750033699e-5 * t247 * sigma[0] / t7 / t465 / rho[0] * t473);
548   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
549     v4rho4[0] = -0.5e1 / 0.27e2 * t4 * t6 * t18 * t31 + 0.4e1 / 0.9e1 * t4 * t198 * t69 - t4 * t93 * t137 - 0.2e1 * t4 * t8 * t261 - 0.3e1 / 0.8e1 * t37 * t38 * t479;
550 
551   t559 = t18 * t225;
552   t573 = -0.12263609725822668962e-192 * t105 * t414 * t154 + 0.66892416686305467067e-193 * t105 * t218 * t293 + 0.33446208343152733533e-193 * t105 * t218 * t296 - 0.10914817924168113389e-189 * t226 * t169 * t112 - 0.11025068610270821606e-191 * t111 * t315 * t66 - 0.11025068610270821606e-191 * t111 * t169 * t134 - 0.36750228700902738685e-192 * t111 * t85 * t258 - 0.41389682824651507748e-192 * t14 * t559 * t227 - 0.12542328128682275075e-193 * t326 * t110 * t66 * t134 - 0.12263609725822668962e-194 * t105 * t214 * t169 + 0.33446208343152733533e-195 * t105 * t106 * t315;
553   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
554     v4rho3sigma[0] = t197 * t1 * t89 / 0.9e1 - t265 * t173 / 0.2e1 - 0.3e1 / 0.2e1 * t141 * t319 - 0.3e1 / 0.8e1 * t37 * t38 * (0.33111746259721206198e-191 * t105 * t421 * t290 - 0.36750228700902738685e-194 * t47 * (-0.1232e4 / 0.27e2 * t14 * t209 * t23 + 0.51167577030977359778e1 * t54 * t237 * t62 * sigma[0] - 0.96840172523411230156e-2 * t243 * t131 * t55 + 0.11478887496814161219e-3 * t250 * t253 * t125 * t9 * t13 - 0.46952273766562626373e-6 * t247 / t7 / t465 * t473) + 0.33446208343152733533e-193 * t14 * t219 - 0.41734945439179794323e0 * t14 * t209 * t28 - 0.12263609725822668962e-194 * t14 * t215 + 0.33446208343152733533e-195 * t14 * t222 - 0.41807760428940916917e-196 * t14 * t147 * t258 - 0.35655071885615837072e-188 * t432 * t85 * t227 - 0.10914817924168113389e-189 * t226 * t154 * t134 + 0.19076726240168596163e-194 * t105 * t410 * t85 + t573);
555 
556   t621 = t169 * t169;
557   t630 = -0.36750228700902738685e-194 * t47 * (-0.76852487240203149073e0 * t54 * t302 + 0.23666793382550744540e-2 * t306 * sigma[0] - 0.36274574252572840931e-4 * t311 * t253 * t55 * t9 * t13 + 0.17607102662460984890e-6 * t125 / t7 / t126 / t118 * t473) - 0.82779365649303015496e-192 * t326 * t225 * t85 * t112 - 0.40878699086075563207e-193 * t105 * t414 * t179 + 0.22074497506480804132e-191 * t105 * t421 * t337 + 0.44594944457536978045e-193 * t105 * t218 * t340 - 0.40878699086075563207e-195 * t105 * t214 * t189 + 0.22297472228768489022e-193 * t105 * t218 * t346 - 0.36382726413893711298e-190 * t226 * t179 * t134 - 0.73500457401805477370e-192 * t111 * t621 - 0.73500457401805477370e-192 * t111 * t85 * t315 - 0.73500457401805477370e-192 * t111 * t359 * t66;
558   t663 = -0.36750228700902738685e-192 * t111 * t189 * t134 - 0.81757398172151126414e-195 * t14 * t281 + 0.44594944457536978044e-193 * t14 * t41 * t328 - 0.16723104171576366767e-193 * t326 * t110 * t169 * t66 - 0.83615520857881833834e-194 * t326 * t327 * t134 + 0.22297472228768489022e-195 * t105 * t106 * t359 + 0.44594944457536978044e-195 * t14 * t287 - 0.83615520857881833834e-196 * t14 * t147 * t315 - 0.35655071885615837072e-188 * t432 * t179 * t112 - 0.14553090565557484519e-189 * t226 * t154 * t169 - 0.36382726413893711298e-190 * t226 * t189 * t112;
559   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
560     v4rho2sigma2[0] = -t265 * t193 / 0.6e1 - t141 * t363 - 0.3e1 / 0.8e1 * t37 * t38 * (t630 + t663);
561 
562   t727 = 0.33446208343152733533e-193 * t14 * t334 - 0.12416904847395452324e-191 * t326 * t225 * t179 * t66 - 0.25084656257364550150e-193 * t326 * t327 * t169 + 0.33446208343152733533e-195 * t14 * t343 - 0.12542328128682275075e-193 * t326 * t110 * t189 * t66 - 0.12542328128682275075e-195 * t14 * t147 * t359 + 0.11037248753240402066e-191 * t105 * t421 * t372 - 0.35655071885615837072e-188 * t432 * t372 * t66 - 0.10914817924168113389e-189 * t226 * t179 * t169 + 0.33446208343152733536e-193 * t105 * t218 * t375 - 0.10914817924168113390e-189 * t226 * t375 * t66 - 0.11025068610270821606e-191 * t111 * t169 * t189 - 0.11025068610270821606e-191 * t111 * t85 * t359 + 0.11148736114384244511e-195 * t105 * t106 * t385 - 0.36750228700902738685e-192 * t111 * t385 * t66 - 0.36750228700902738685e-194 * t47 * (-0.38263677923884437281e-3 * t351 + 0.10156880790720395460e-4 * t355 * t253 * t105 - 0.66026634984228693337e-7 * t55 / t7 / t458 * t473);
563   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
564     v4rhosigma3[0] = -t141 * t389 / 0.2e1 - 0.3e1 / 0.8e1 * t37 * t38 * t727;
565 
566   t740 = t179 * t179;
567   t746 = t189 * t189;
568   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
569     v4sigma4[0] = -0.3e1 / 0.8e1 * t37 * t38 * (-0.16555873129860603099e-191 * t14 * t559 * t372 - 0.50169312514729100301e-193 * t326 * t327 * t189 - 0.16723104171576366767e-195 * t14 * t147 * t385 - 0.35655071885615837072e-188 * t432 * t740 - 0.21829635848336226779e-189 * t226 * t179 * t189 - 0.11025068610270821606e-191 * t111 * t746 - 0.14700091480361095474e-191 * t111 * t85 * t385 - 0.36750228700902738685e-194 * t47 * (-0.21764744551543704558e-5 * t381 * t253 * t14 + 0.24759988119085760001e-7 * sigma[0] / t7 / t248 * t473));
570 
571 #ifndef XC_DONT_COMPILE_MXC
572 
573   if(order < 5) return;
574 
575 
576 #endif
577 
578 #endif
579 
580 #endif
581 
582 #endif
583 
584 #endif
585 
586 
587 }
588 
589 
590 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 *))591 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, double *zk, GGA_OUT_PARAMS_NO_EXC(double *))
592 {
593 
594 #ifndef XC_DONT_COMPILE_EXC
595   double t1, t3, t4, t5, t6, t7, t8, t9;
596   double t10, t11, t13, t14, t15, t16, t17, t18;
597   double t19, t20, t21, t22, t23, t24, t25, t26;
598   double t28, t29, t32, t33, t35, t36, t38, t41;
599   double t43, t45, t46, t47, t48, t49, t50, t51;
600   double t53, t54, t57, t58, t60, t61, t63, t66;
601   double t68;
602 
603 #ifndef XC_DONT_COMPILE_VXC
604   double t70, t71, t72, t73, t74, t75, t77, t78;
605   double t80, t81, t82, t83, t84, t86, t87, t88;
606   double t90, t91, t95, t96, t100, t102, t103, t104;
607   double t105, t106, t108, t110, t111, t115, t118, t120;
608   double t121, t122, t123, t124, t126, t127, t128, t130;
609   double t131, t135, t136, t138, t139, t140, t141, t143;
610   double t144, t145, t147, t148, t152, t153, t157, t158;
611   double t159, t161, t163, t164, t168, t171, t173, t174;
612   double t178, t179, t180, t186, t188, t193, t196, t200;
613   double t206, t208, t213, t216;
614 
615 #ifndef XC_DONT_COMPILE_FXC
616   double t220, t222, t223, t224, t225, t226, t227, t229;
617   double t230, t231, t233, t234, t235, t237, t238, t239;
618   double t240, t241, t242, t244, t245, t247, t248, t250;
619   double t251, t253, t254, t256, t257, t261, t262, t263;
620   double t266, t267, t268, t274, t276, t281, t282, t283;
621   double t284, t287, t290, t293, t295, t296, t297, t298;
622   double t299, t300, t301, t303, t304, t305, t307, t308;
623   double t309, t310, t312, t313, t314, t316, t317, t318;
624   double t320, t321, t322, t323, t325, t328, t330, t332;
625   double t335, t336, t337, t339, t343, t344, t345, t346;
626   double t348, t351, t353, t354, t356, t359, t361, t364;
627   double t366, t369, t370, t372, t373, t374, t375, t376;
628   double t378, t379, t380, t381, t383, t384, t385, t386;
629   double t388, t389, t390, t392, t393, t394, t396, t397;
630   double t401, t402, t403, t406, t407, t408, t414, t416;
631   double t421, t422, t423, t424, t427, t430, t433, t435;
632   double t436, t437, t440, t442, t443, t450, t454, t457;
633   double t463, t467, t468, t472, t475, t480, t481, t482;
634   double t487, t492, t498, t502, t505, t511, t515, t516;
635   double t520, t523, t524, t530, t536, t540, t543, t544;
636   double t550, t556, t560, t563, t564;
637 
638 #ifndef XC_DONT_COMPILE_KXC
639   double t571, t573, t577, t578, t579, t580, t582, t583;
640   double t584, t586, t588, t590, t591, t592, t594, t595;
641   double t597, t598, t599, t600, t602, t603, t604, t606;
642   double t608, t610, t611, t612, t614, t615, t616, t618;
643   double t619, t621, t624, t626, t628, t629, t630, t632;
644   double t633, t634, t635, t636, t638, t639, t641, t642;
645   double t643, t645, t646, t648, t649, t651, t654, t655;
646   double t659, t660, t663, t664, t667, t670, t671, t672;
647   double t675, t682, t687, t688, t692, t693, t695, t697;
648   double t698, t700, t703, t706, t708, t709, t711, t713;
649   double t714, t716, t717, t718, t720, t721, t722, t727;
650   double t728, t731, t732, t733, t735, t737, t739, t740;
651   double t741, t742, t744, t746, t747, t748, t750, t752;
652   double t754, t755, t756, t758, t760, t762, t764, t766;
653   double t768, t769, t771, t774, t777, t780, t784, t786;
654   double t787, t789, t791, t793, t796, t798, t800, t802;
655   double t804, t806, t808, t810, t813, t814, t816, t817;
656   double t818, t820, t829, t832, t837, t839, t841, t843;
657   double t845, t847, t849, t851, t853, t855, t857, t859;
658   double t861, t863, t866, t870, t872, t874, t876, t878;
659   double t880, t881, t883, t885, t887, t888, t889, t891;
660   double t893, t895, t898, t901, t904, t907, t911, t918;
661   double t921, t929, t930, t932, t935, t936, t938, t939;
662   double t940, t942, t943, t944, t946, t947, t949, t951;
663   double t952, t953, t954, t956, t957, t959, t960, t964;
664   double t965, t968, t969, t972, t975, t976, t977, t980;
665   double t987, t992, t993, t997, t998, t1000, t1002, t1003;
666   double t1005, t1008, t1011, t1013, t1014, t1018, t1019, t1021;
667   double t1022, t1023, t1025, t1026, t1027, t1028, t1030, t1031;
668   double t1032, t1034, t1035, t1036, t1040, t1044, t1046, t1047;
669   double t1051, t1054, t1063, t1070, t1076, t1079, t1082, t1085;
670   double t1091, t1095, t1098, t1100, t1104, t1107, t1112, t1116;
671   double t1117, t1121, t1127, t1128, t1129, t1131, t1132, t1133;
672   double t1134, t1138, t1141, t1142, t1144, t1148, t1150, t1151;
673   double t1155, t1158, t1161, t1166, t1170, t1175, t1179, t1182;
674   double t1190, t1197, t1203, t1206, t1209, t1212, t1218, t1222;
675   double t1225, t1227, t1231, t1234, t1235, t1238, t1241, t1246;
676   double t1247, t1248, t1254, t1257, t1260, t1263, t1266, t1271;
677   double t1275, t1279, t1282, t1283, t1286, t1288, t1292, t1297;
678   double t1302, t1303, t1304, t1310, t1313, t1316, t1319, t1322;
679   double t1327, t1331, t1335, t1338, t1339, t1348, t1351, t1357;
680   double t1361, t1364, t1365, t1374, t1377, t1383, t1387, t1390;
681   double t1391;
682 
683 #ifndef XC_DONT_COMPILE_LXC
684   double t1394, t1396, t1398, t1406, t1410, t1413, t1417, t1421;
685   double t1428, t1438, t1439, t1440, t1446, t1465, t1472, t1480;
686   double t1486, t1491, t1494, t1497, t1508, t1522, t1524, t1526;
687   double t1527, t1528, t1533, t1535, t1537, t1539, t1541, t1547;
688   double t1549, t1550, t1552, t1553, t1554, t1559, t1561, t1563;
689   double t1565, t1567, t1573, t1580, t1583, t1586, t1589, t1595;
690   double t1598, t1601, t1621, t1623, t1636, t1641, t1646, t1647;
691   double t1648, t1650, t1662, t1663, t1671, t1673, t1676, t1682;
692   double t1689, t1695, t1702, t1705, t1708, t1711, t1712, t1713;
693   double t1715, t1722, t1726, t1728, t1733, t1740, t1743, t1744;
694   double t1745, t1754, t1757, t1763, t1766, t1768, t1778, t1780;
695   double t1784, t1791, t1806, t1807, t1809, t1833, t1838, t1840;
696   double t1843, t1850, t1857, t1860, t1875, t1884, t1895, t1911;
697   double t1918, t1925, t1932, t1939, t1946, t1949, t1956, t1969;
698   double t1977, t2016, t2019, t2021, t2027, t2030, t2033, t2038;
699   double t2041, t2044, t2047, t2050, t2057, t2063, t2066, t2083;
700   double t2108, t2112, t2118, t2121, t2131, t2134, t2138, t2145;
701   double t2187, t2216, t2222, t2225, t2228, t2230, t2233, t2236;
702   double t2239, t2242, t2269, t2288, t2292, t2296, t2300, t2304;
703   double t2311, t2321, t2322, t2323, t2329, t2348, t2355, t2363;
704   double t2369, t2396, t2401, t2406, t2411, t2429, t2438, t2462;
705   double t2472, t2484, t2514, t2568, t2590, t2593, t2597, t2599;
706   double t2602, t2612, t2615, t2618, t2621, t2624, t2628, t2644;
707   double t2657, t2658, t2660, t2662, t2664, t2665, t2676, t2680;
708   double t2687, t2690, t2691, t2693, t2706, t2709, t2722, t2734;
709   double t2751, t2759, t2768, t2777, t2821, t2880, t2915, t2916;
710   double t2918, t2964, t2982, t3012, t3018, t3019, t3031, t3046;
711   double t3056, t3104, t3116, t3165, t3170, t3233, t3237, t3308;
712   double t3321, t3327, t3357, t3363;
713 #endif
714 
715 #endif
716 
717 #endif
718 
719 #endif
720 
721 #endif
722 
723 
724 
725   t1 = M_CBRT3;
726   t3 = POW_1_3(0.1e1 / M_PI);
727   t4 = t1 * t3;
728   t5 = M_CBRT4;
729   t6 = t5 * t5;
730   t7 = t4 * t6;
731   t8 = rho[0] - rho[1];
732   t9 = rho[0] + rho[1];
733   t10 = 0.1e1 / t9;
734   t11 = t8 * t10;
735   t13 = 0.1e1 / 0.2e1 + t11 / 0.2e1;
736   t14 = POW_1_3(t13);
737   t15 = t14 * t13;
738   t16 = POW_1_3(t9);
739   t17 = t15 * t16;
740   t18 = M_CBRT6;
741   t19 = M_PI * M_PI;
742   t20 = POW_1_3(t19);
743   t21 = t20 * t20;
744   t22 = 0.1e1 / t21;
745   t23 = t18 * t22;
746   t24 = rho[0] * rho[0];
747   t25 = POW_1_3(rho[0]);
748   t26 = t25 * t25;
749   t28 = 0.1e1 / t26 / t24;
750   t29 = sigma[0] * t28;
751   t32 = 0.10000000000000000000e-1 + 0.11376190545424808251e-1 * t23 * t29;
752   t33 = 0.1e1 / t32;
753   t35 = t23 * t29 * t33;
754   t36 = pow(t35, 0.1010e3);
755   t38 = 0.1e1 - 0.45256672948257153200e-196 * t36;
756   t41 = 0.99196000000000000000e0 + 0.8040e0 * t32 * t38;
757   t43 = t7 * t17 * t41;
758   t45 = 0.1e1 / 0.2e1 - t11 / 0.2e1;
759   t46 = POW_1_3(t45);
760   t47 = t46 * t45;
761   t48 = t47 * t16;
762   t49 = rho[1] * rho[1];
763   t50 = POW_1_3(rho[1]);
764   t51 = t50 * t50;
765   t53 = 0.1e1 / t51 / t49;
766   t54 = sigma[2] * t53;
767   t57 = 0.10000000000000000000e-1 + 0.11376190545424808251e-1 * t23 * t54;
768   t58 = 0.1e1 / t57;
769   t60 = t23 * t54 * t58;
770   t61 = pow(t60, 0.1010e3);
771   t63 = 0.1e1 - 0.45256672948257153200e-196 * t61;
772   t66 = 0.99196000000000000000e0 + 0.8040e0 * t57 * t63;
773   t68 = t7 * t48 * t66;
774   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
775     zk[0] = -0.3e1 / 0.8e1 * t43 - 0.3e1 / 0.8e1 * t68;
776 
777 #ifndef XC_DONT_COMPILE_VXC
778 
779   if(order < 1) return;
780 
781 
782   t70 = 0.3e1 / 0.8e1 * t43;
783   t71 = 0.3e1 / 0.8e1 * t68;
784   t72 = t14 * t16;
785   t73 = t9 * t9;
786   t74 = 0.1e1 / t73;
787   t75 = t8 * t74;
788   t77 = t10 / 0.2e1 - t75 / 0.2e1;
789   t78 = t41 * t77;
790   t80 = t7 * t72 * t78;
791   t81 = t80 / 0.2e1;
792   t82 = t16 * t16;
793   t83 = 0.1e1 / t82;
794   t84 = t15 * t83;
795   t86 = t7 * t84 * t41;
796   t87 = t86 / 0.8e1;
797   t88 = t24 * rho[0];
798   t90 = 0.1e1 / t26 / t88;
799   t91 = sigma[0] * t90;
800   t95 = pow(t35, 0.1000e3);
801   t96 = t32 * t95;
802   t100 = t18 * t18;
803   t102 = 0.1e1 / t20 / t19;
804   t103 = t100 * t102;
805   t104 = sigma[0] * sigma[0];
806   t105 = t24 * t24;
807   t106 = t105 * t24;
808   t108 = 0.1e1 / t25 / t106;
809   t110 = t32 * t32;
810   t111 = 0.1e1 / t110;
811   t115 = -0.8e1 / 0.3e1 * t23 * t91 * t33 + 0.30336508121132822003e-1 * t103 * t104 * t108 * t111;
812   t118 = -0.24390552529390788890e-1 * t23 * t91 * t38 - 0.36750228700902738685e-194 * t96 * t115;
813   t120 = t7 * t17 * t118;
814   t121 = 0.3e1 / 0.8e1 * t120;
815   t122 = t46 * t16;
816   t123 = -t77;
817   t124 = t66 * t123;
818   t126 = t7 * t122 * t124;
819   t127 = t126 / 0.2e1;
820   t128 = t47 * t83;
821   t130 = t7 * t128 * t66;
822   t131 = t130 / 0.8e1;
823   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
824     vrho[0] = -t70 - t71 + t9 * (-t81 - t87 - t121 - t127 - t131);
825 
826   t135 = -t10 / 0.2e1 - t75 / 0.2e1;
827   t136 = t41 * t135;
828   t138 = t7 * t72 * t136;
829   t139 = t138 / 0.2e1;
830   t140 = -t135;
831   t141 = t66 * t140;
832   t143 = t7 * t122 * t141;
833   t144 = t143 / 0.2e1;
834   t145 = t49 * rho[1];
835   t147 = 0.1e1 / t51 / t145;
836   t148 = sigma[2] * t147;
837   t152 = pow(t60, 0.1000e3);
838   t153 = t57 * t152;
839   t157 = sigma[2] * sigma[2];
840   t158 = t49 * t49;
841   t159 = t158 * t49;
842   t161 = 0.1e1 / t50 / t159;
843   t163 = t57 * t57;
844   t164 = 0.1e1 / t163;
845   t168 = -0.8e1 / 0.3e1 * t23 * t148 * t58 + 0.30336508121132822003e-1 * t103 * t157 * t161 * t164;
846   t171 = -0.24390552529390788890e-1 * t23 * t148 * t63 - 0.36750228700902738685e-194 * t153 * t168;
847   t173 = t7 * t48 * t171;
848   t174 = 0.3e1 / 0.8e1 * t173;
849   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
850     vrho[1] = -t70 - t71 + t9 * (-t139 - t87 - t144 - t131 - t174);
851 
852   t178 = t16 * t9 * t1;
853   t179 = t178 * t3;
854   t180 = t6 * t15;
855   t186 = t105 * rho[0];
856   t188 = 0.1e1 / t25 / t186;
857   t193 = t23 * t28 * t33 - 0.11376190545424808251e-1 * t103 * sigma[0] * t188 * t111;
858   t196 = 0.91464571985215458338e-2 * t23 * t28 * t38 - 0.36750228700902738685e-194 * t96 * t193;
859   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
860     vsigma[0] = -0.3e1 / 0.8e1 * t179 * t180 * t196;
861 
862   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
863     vsigma[1] = 0.0e0;
864 
865   t200 = t6 * t47;
866   t206 = t158 * rho[1];
867   t208 = 0.1e1 / t50 / t206;
868   t213 = t23 * t53 * t58 - 0.11376190545424808251e-1 * t103 * sigma[2] * t208 * t164;
869   t216 = 0.91464571985215458338e-2 * t23 * t53 * t63 - 0.36750228700902738685e-194 * t153 * t213;
870   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
871     vsigma[2] = -0.3e1 / 0.8e1 * t179 * t200 * t216;
872 
873 #ifndef XC_DONT_COMPILE_FXC
874 
875   if(order < 2) return;
876 
877 
878   t220 = t86 / 0.4e1;
879   t222 = t130 / 0.4e1;
880   t223 = t14 * t14;
881   t224 = 0.1e1 / t223;
882   t225 = t224 * t16;
883   t226 = t77 * t77;
884   t227 = t41 * t226;
885   t229 = t7 * t225 * t227;
886   t230 = t229 / 0.6e1;
887   t231 = t14 * t83;
888   t233 = t7 * t231 * t78;
889   t234 = t233 / 0.3e1;
890   t235 = t118 * t77;
891   t237 = t7 * t72 * t235;
892   t238 = t73 * t9;
893   t239 = 0.1e1 / t238;
894   t240 = t8 * t239;
895   t241 = -t74 + t240;
896   t242 = t41 * t241;
897   t244 = t7 * t72 * t242;
898   t245 = t244 / 0.2e1;
899   t247 = 0.1e1 / t82 / t9;
900   t248 = t15 * t247;
901   t250 = t7 * t248 * t41;
902   t251 = t250 / 0.12e2;
903   t253 = t7 * t84 * t118;
904   t254 = t253 / 0.4e1;
905   t256 = 0.1e1 / t26 / t105;
906   t257 = sigma[0] * t256;
907   t261 = t23 * sigma[0];
908   t262 = t90 * t95;
909   t263 = t262 * t115;
910   t266 = pow(t35, 0.990e2);
911   t267 = t32 * t266;
912   t268 = t115 * t115;
913   t274 = t105 * t88;
914   t276 = 0.1e1 / t25 / t274;
915   t281 = t104 * sigma[0];
916   t282 = t105 * t105;
917   t283 = t282 * t24;
918   t284 = 0.1e1 / t283;
919   t287 = 0.1e1 / t110 / t32;
920   t290 = 0.88e2 / 0.9e1 * t23 * t257 * t33 - 0.27302857309019539803e0 * t103 * t104 * t276 * t111 + 0.11337386051521314750e-3 * t281 * t284 * t287;
921   t293 = 0.89432025941099559263e-1 * t23 * t257 * t38 + 0.22297472228768489022e-195 * t261 * t263 - 0.36750228700902738685e-192 * t267 * t268 - 0.36750228700902738685e-194 * t96 * t290;
922   t295 = t7 * t17 * t293;
923   t296 = 0.3e1 / 0.8e1 * t295;
924   t297 = t46 * t46;
925   t298 = 0.1e1 / t297;
926   t299 = t298 * t16;
927   t300 = t123 * t123;
928   t301 = t66 * t300;
929   t303 = t7 * t299 * t301;
930   t304 = t303 / 0.6e1;
931   t305 = t46 * t83;
932   t307 = t7 * t305 * t124;
933   t308 = t307 / 0.3e1;
934   t309 = -t241;
935   t310 = t66 * t309;
936   t312 = t7 * t122 * t310;
937   t313 = t312 / 0.2e1;
938   t314 = t47 * t247;
939   t316 = t7 * t314 * t66;
940   t317 = t316 / 0.12e2;
941   t318 = -t230 - t234 - t237 - t245 + t251 - t254 - t296 - t304 - t308 - t313 + t317;
942   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
943     v2rho2[0] = -t80 - t220 - 0.3e1 / 0.4e1 * t120 - t126 - t222 + t9 * t318;
944 
945   t320 = t6 * t224;
946   t321 = t4 * t320;
947   t322 = t16 * t41;
948   t323 = t135 * t77;
949   t325 = t321 * t322 * t323;
950   t328 = t7 * t231 * t136;
951   t330 = t118 * t135;
952   t332 = t7 * t72 * t330;
953   t335 = 0.1e1 / t82 / t73;
954   t336 = t14 * t335;
955   t337 = t41 * t8;
956   t339 = t7 * t336 * t337;
957   t343 = t6 * t298;
958   t344 = t4 * t343;
959   t345 = t16 * t66;
960   t346 = t140 * t123;
961   t348 = t344 * t345 * t346;
962   t351 = t7 * t305 * t141;
963   t353 = t46 * t335;
964   t354 = t66 * t8;
965   t356 = t7 * t353 * t354;
966   t359 = t171 * t123;
967   t361 = t7 * t122 * t359;
968   t364 = t7 * t128 * t171;
969   t366 = -t325 / 0.6e1 - t328 / 0.6e1 - t332 / 0.2e1 - t339 / 0.2e1 - t233 / 0.6e1 + t251 - t253 / 0.8e1 - t348 / 0.6e1 - t351 / 0.6e1 + t356 / 0.2e1 - t307 / 0.6e1 + t317 - t361 / 0.2e1 - t364 / 0.8e1;
970   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
971     v2rho2[1] = t9 * t366 - t121 - t127 - t139 - t144 - t174 - t220 - t222 - t81;
972 
973   t369 = t135 * t135;
974   t370 = t41 * t369;
975   t372 = t7 * t225 * t370;
976   t373 = t372 / 0.6e1;
977   t374 = t328 / 0.3e1;
978   t375 = t74 + t240;
979   t376 = t41 * t375;
980   t378 = t7 * t72 * t376;
981   t379 = t378 / 0.2e1;
982   t380 = t140 * t140;
983   t381 = t66 * t380;
984   t383 = t7 * t299 * t381;
985   t384 = t383 / 0.6e1;
986   t385 = t351 / 0.3e1;
987   t386 = t171 * t140;
988   t388 = t7 * t122 * t386;
989   t389 = -t375;
990   t390 = t66 * t389;
991   t392 = t7 * t122 * t390;
992   t393 = t392 / 0.2e1;
993   t394 = t364 / 0.4e1;
994   t396 = 0.1e1 / t51 / t158;
995   t397 = sigma[2] * t396;
996   t401 = t23 * sigma[2];
997   t402 = t147 * t152;
998   t403 = t402 * t168;
999   t406 = pow(t60, 0.990e2);
1000   t407 = t57 * t406;
1001   t408 = t168 * t168;
1002   t414 = t158 * t145;
1003   t416 = 0.1e1 / t50 / t414;
1004   t421 = t157 * sigma[2];
1005   t422 = t158 * t158;
1006   t423 = t422 * t49;
1007   t424 = 0.1e1 / t423;
1008   t427 = 0.1e1 / t163 / t57;
1009   t430 = 0.88e2 / 0.9e1 * t23 * t397 * t58 - 0.27302857309019539803e0 * t103 * t157 * t416 * t164 + 0.11337386051521314750e-3 * t421 * t424 * t427;
1010   t433 = 0.89432025941099559263e-1 * t23 * t397 * t63 + 0.22297472228768489022e-195 * t401 * t403 - 0.36750228700902738685e-192 * t407 * t408 - 0.36750228700902738685e-194 * t153 * t430;
1011   t435 = t7 * t48 * t433;
1012   t436 = 0.3e1 / 0.8e1 * t435;
1013   t437 = -t373 - t374 - t379 + t251 - t384 - t385 - t388 - t393 + t317 - t394 - t436;
1014   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1015     v2rho2[2] = -t138 - t220 - t143 - t222 - 0.3e1 / 0.4e1 * t173 + t9 * t437;
1016 
1017   t440 = t7 * t17 * t196;
1018   t442 = t6 * t14;
1019   t443 = t196 * t77;
1020   t450 = t28 * t95;
1021   t454 = t262 * t193;
1022   t457 = t193 * t115;
1023   t463 = t108 * t111;
1024   t467 = t282 * rho[0];
1025   t468 = 0.1e1 / t467;
1026   t472 = -0.8e1 / 0.3e1 * t23 * t90 * t33 + 0.91009524363398466008e-1 * t103 * t463 * sigma[0] - 0.42515197693204930312e-4 * t104 * t468 * t287;
1027   t475 = -0.24390552529390788890e-1 * t23 * t90 * t38 - 0.41807760428940916917e-196 * t23 * t450 * t115 + 0.11148736114384244511e-195 * t261 * t454 - 0.36750228700902738685e-192 * t267 * t457 - 0.36750228700902738685e-194 * t96 * t472;
1028   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1029     v2rhosigma[0] = -t440 / 0.2e1 - t179 * t442 * t443 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t180 * t475;
1030 
1031   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1032     v2rhosigma[1] = 0.0e0;
1033 
1034   t480 = t7 * t48 * t216;
1035   t481 = t6 * t46;
1036   t482 = t216 * t123;
1037   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1038     v2rhosigma[2] = -t179 * t481 * t482 / 0.2e1 - t480 / 0.2e1;
1039 
1040   t487 = t442 * t196 * t135;
1041   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1042     v2rhosigma[3] = -t179 * t487 / 0.2e1 - t440 / 0.2e1;
1043 
1044   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1045     v2rhosigma[4] = 0.0e0;
1046 
1047   t492 = t481 * t216 * t140;
1048   t498 = t53 * t152;
1049   t502 = t402 * t213;
1050   t505 = t213 * t168;
1051   t511 = t161 * t164;
1052   t515 = t422 * rho[1];
1053   t516 = 0.1e1 / t515;
1054   t520 = -0.8e1 / 0.3e1 * t23 * t147 * t58 + 0.91009524363398466008e-1 * t103 * t511 * sigma[2] - 0.42515197693204930312e-4 * t157 * t516 * t427;
1055   t523 = -0.24390552529390788890e-1 * t23 * t147 * t63 - 0.41807760428940916917e-196 * t23 * t498 * t168 + 0.11148736114384244511e-195 * t401 * t502 - 0.36750228700902738685e-192 * t407 * t505 - 0.36750228700902738685e-194 * t153 * t520;
1056   t524 = t200 * t523;
1057   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1058     v2rhosigma[5] = -t480 / 0.2e1 - t179 * t492 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t524;
1059 
1060   t530 = t193 * t193;
1061   t536 = 0.1e1 / t282;
1062   t540 = -0.22752381090849616502e-1 * t103 * t188 * t111 + 0.15943199134951848867e-4 * sigma[0] * t536 * t287;
1063   t543 = -0.83615520857881833834e-196 * t23 * t450 * t193 - 0.36750228700902738685e-192 * t267 * t530 - 0.36750228700902738685e-194 * t96 * t540;
1064   t544 = t180 * t543;
1065   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1066     v2sigma2[0] = -0.3e1 / 0.8e1 * t179 * t544;
1067 
1068   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1069     v2sigma2[1] = 0.0e0;
1070 
1071   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1072     v2sigma2[2] = 0.0e0;
1073 
1074   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1075     v2sigma2[3] = 0.0e0;
1076 
1077   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1078     v2sigma2[4] = 0.0e0;
1079 
1080   t550 = t213 * t213;
1081   t556 = 0.1e1 / t422;
1082   t560 = -0.22752381090849616502e-1 * t103 * t208 * t164 + 0.15943199134951848867e-4 * sigma[2] * t556 * t427;
1083   t563 = -0.83615520857881833834e-196 * t23 * t498 * t213 - 0.36750228700902738685e-192 * t407 * t550 - 0.36750228700902738685e-194 * t153 * t560;
1084   t564 = t200 * t563;
1085   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1086     v2sigma2[5] = -0.3e1 / 0.8e1 * t179 * t564;
1087 
1088 #ifndef XC_DONT_COMPILE_KXC
1089 
1090   if(order < 3) return;
1091 
1092 
1093   t571 = t250 / 0.4e1;
1094   t573 = t316 / 0.4e1;
1095   t577 = 0.1e1 / t223 / t13;
1096   t578 = t577 * t16;
1097   t579 = t226 * t77;
1098   t580 = t41 * t579;
1099   t582 = t7 * t578 * t580;
1100   t583 = t582 / 0.9e1;
1101   t584 = t224 * t83;
1102   t586 = t7 * t584 * t227;
1103   t588 = t77 * t241;
1104   t590 = t321 * t322 * t588;
1105   t591 = t590 / 0.2e1;
1106   t592 = t14 * t247;
1107   t594 = t7 * t592 * t78;
1108   t595 = t594 / 0.3e1;
1109   t597 = 0.1e1 / t297 / t45;
1110   t598 = t597 * t16;
1111   t599 = t300 * t123;
1112   t600 = t66 * t599;
1113   t602 = t7 * t598 * t600;
1114   t603 = t602 / 0.9e1;
1115   t604 = t298 * t83;
1116   t606 = t7 * t604 * t301;
1117   t608 = t123 * t309;
1118   t610 = t344 * t345 * t608;
1119   t611 = t610 / 0.2e1;
1120   t612 = t46 * t247;
1121   t614 = t7 * t612 * t124;
1122   t615 = t614 / 0.3e1;
1123   t616 = t118 * t226;
1124   t618 = t7 * t225 * t616;
1125   t619 = t618 / 0.2e1;
1126   t621 = t7 * t231 * t235;
1127   t624 = t7 * t231 * t242;
1128   t626 = t293 * t77;
1129   t628 = t7 * t72 * t626;
1130   t629 = 0.3e1 / 0.2e1 * t628;
1131   t630 = t118 * t241;
1132   t632 = t7 * t72 * t630;
1133   t633 = 0.3e1 / 0.2e1 * t632;
1134   t634 = t73 * t73;
1135   t635 = 0.1e1 / t634;
1136   t636 = t8 * t635;
1137   t638 = 0.3e1 * t239 - 0.3e1 * t636;
1138   t639 = t41 * t638;
1139   t641 = t7 * t72 * t639;
1140   t642 = t641 / 0.2e1;
1141   t643 = t15 * t335;
1142   t645 = t7 * t643 * t41;
1143   t646 = 0.5e1 / 0.36e2 * t645;
1144   t648 = t7 * t248 * t118;
1145   t649 = t648 / 0.4e1;
1146   t651 = t7 * t84 * t293;
1147   t654 = 0.1e1 / t26 / t186;
1148   t655 = sigma[0] * t654;
1149   t659 = t256 * t95;
1150   t660 = t659 * t115;
1151   t663 = t90 * t266;
1152   t664 = t663 * t268;
1153   t667 = t262 * t290;
1154   t670 = pow(t35, 0.980e2);
1155   t671 = t32 * t670;
1156   t672 = t268 * t115;
1157   t675 = t115 * t290;
1158   t682 = 0.1e1 / t25 / t282;
1159   t687 = t282 * t88;
1160   t688 = 0.1e1 / t687;
1161   t692 = t104 * t104;
1162   t693 = t282 * t186;
1163   t695 = 0.1e1 / t26 / t693;
1164   t697 = t110 * t110;
1165   t698 = 0.1e1 / t697;
1166   t700 = t698 * t18 * t22;
1167   t703 = -0.1232e4 / 0.27e2 * t23 * t655 * t33 + 0.22988331709569538451e1 * t103 * t104 * t682 * t111 - 0.21541033497890498025e-2 * t281 * t688 * t287 + 0.10318101120731830309e-4 * t692 * t695 * t700;
1168   t706 = -0.41734945439179794323e0 * t23 * t655 * t38 - 0.12263609725822668962e-194 * t261 * t660 + 0.33446208343152733533e-193 * t261 * t664 + 0.33446208343152733533e-195 * t261 * t667 - 0.36382726413893711298e-190 * t671 * t672 - 0.11025068610270821606e-191 * t267 * t675 - 0.36750228700902738685e-194 * t96 * t703;
1169   t708 = t7 * t17 * t706;
1170   t709 = 0.3e1 / 0.8e1 * t708;
1171   t711 = t7 * t305 * t310;
1172   t713 = -t638;
1173   t714 = t66 * t713;
1174   t716 = t7 * t122 * t714;
1175   t717 = t716 / 0.2e1;
1176   t718 = t47 * t335;
1177   t720 = t7 * t718 * t66;
1178   t721 = 0.5e1 / 0.36e2 * t720;
1179   t722 = -t624 / 0.2e1 - t629 - t633 - t642 - t646 + t649 - 0.3e1 / 0.8e1 * t651 - t709 - t711 / 0.2e1 - t717 - t721;
1180   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1181     v3rho3[0] = -t229 / 0.2e1 - t233 - t303 / 0.2e1 - t307 - 0.3e1 * t237 - 0.3e1 / 0.2e1 * t244 + t571 - 0.3e1 / 0.2e1 * t312 + t573 - 0.3e1 / 0.4e1 * t253 - 0.9e1 / 0.8e1 * t295 + t9 * (t583 - t586 / 0.6e1 - t591 + t595 + t603 - t606 / 0.6e1 - t611 + t615 - t619 - t621 + t722);
1182 
1183   t727 = t325 / 0.3e1;
1184   t728 = t348 / 0.3e1;
1185   t731 = t6 * t577;
1186   t732 = t4 * t731;
1187   t733 = t135 * t226;
1188   t735 = t732 * t322 * t733;
1189   t737 = t83 * t41;
1190   t739 = t321 * t737 * t323;
1191   t740 = t739 / 0.9e1;
1192   t741 = t335 * t41;
1193   t742 = t8 * t77;
1194   t744 = t321 * t741 * t742;
1195   t746 = t6 * t597;
1196   t747 = t4 * t746;
1197   t748 = t140 * t300;
1198   t750 = t747 * t345 * t748;
1199   t752 = t83 * t66;
1200   t754 = t344 * t752 * t346;
1201   t755 = t754 / 0.9e1;
1202   t756 = t16 * t118;
1203   t758 = t321 * t756 * t323;
1204   t760 = t135 * t241;
1205   t762 = t321 * t322 * t760;
1206   t764 = t140 * t309;
1207   t766 = t344 * t345 * t764;
1208   t768 = t335 * t66;
1209   t769 = t8 * t123;
1210   t771 = t344 * t768 * t769;
1211   t774 = t7 * t314 * t171;
1212   t777 = t7 * t336 * t41;
1213   t780 = t7 * t353 * t66;
1214   t784 = t735 / 0.9e1 - t740 - t744 / 0.3e1 + t750 / 0.9e1 - t755 - t758 / 0.3e1 - t762 / 0.6e1 - t766 / 0.6e1 + t771 / 0.3e1 + t774 / 0.12e2 - t777 / 0.2e1 + t780 / 0.2e1 - t721 - t646 + t648 / 0.6e1 - t651 / 0.8e1;
1215   t786 = 0.1e1 / t82 / t238;
1216   t787 = t46 * t786;
1217   t789 = t7 * t787 * t354;
1218   t791 = t171 * t309;
1219   t793 = t7 * t122 * t791;
1220   t796 = t7 * t592 * t136;
1221   t798 = t14 * t786;
1222   t800 = t7 * t798 * t337;
1223   t802 = t118 * t8;
1224   t804 = t7 * t336 * t802;
1225   t806 = t7 * t612 * t141;
1226   t808 = t171 * t300;
1227   t810 = t7 * t299 * t808;
1228   t813 = t7 * t305 * t359;
1229   t814 = t813 / 0.3e1;
1230   t816 = t7 * t231 * t330;
1231   t817 = t816 / 0.3e1;
1232   t818 = t293 * t135;
1233   t820 = t7 * t72 * t818;
1234   t829 = -0.7e1 / 0.6e1 * t789 - t793 / 0.2e1 + t796 / 0.9e1 + 0.7e1 / 0.6e1 * t800 - t804 + t806 / 0.9e1 - t810 / 0.6e1 - t814 - t817 - t820 / 0.2e1 - t586 / 0.18e2 + 0.2e1 / 0.9e1 * t594 - t606 / 0.18e2 + 0.2e1 / 0.9e1 * t614 - t621 / 0.3e1 - t624 / 0.6e1 - t711 / 0.6e1;
1235   t832 = -t245 + t571 - t313 + t573 - t332 + t356 - t361 - t394 - t253 / 0.2e1 - t296 + t9 * (t784 + t829);
1236   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1237     v3rho3[1] = -t230 - 0.2e1 / 0.3e1 * t233 - t304 - 0.2e1 / 0.3e1 * t307 - t727 - t374 - t339 - t728 - t385 - t237 + t832;
1238 
1239   t837 = t369 * t77;
1240   t839 = t732 * t322 * t837;
1241   t841 = t135 * t8;
1242   t843 = t321 * t741 * t841;
1243   t845 = t375 * t77;
1244   t847 = t321 * t322 * t845;
1245   t849 = t380 * t123;
1246   t851 = t747 * t345 * t849;
1247   t853 = t140 * t8;
1248   t855 = t344 * t768 * t853;
1249   t857 = t16 * t171;
1250   t859 = t344 * t857 * t346;
1251   t861 = t389 * t123;
1252   t863 = t344 * t345 * t861;
1253   t866 = t7 * t128 * t433;
1254   t870 = t118 * t369;
1255   t872 = t7 * t225 * t870;
1256   t874 = t118 * t375;
1257   t876 = t7 * t72 * t874;
1258   t878 = t839 / 0.9e1 - t843 / 0.3e1 - t847 / 0.6e1 + t851 / 0.9e1 + t855 / 0.3e1 - t859 / 0.3e1 - t863 / 0.6e1 - t740 - t755 - t866 / 0.8e1 + t774 / 0.6e1 - t721 - t646 + t648 / 0.12e2 - t872 / 0.6e1 - t876 / 0.2e1;
1259   t880 = -t239 - 0.3e1 * t636;
1260   t881 = t41 * t880;
1261   t883 = t7 * t72 * t881;
1262   t885 = t171 * t8;
1263   t887 = t7 * t353 * t885;
1264   t888 = -t880;
1265   t889 = t66 * t888;
1266   t891 = t7 * t122 * t889;
1267   t893 = t433 * t123;
1268   t895 = t7 * t122 * t893;
1269   t898 = t7 * t231 * t376;
1270   t901 = t7 * t604 * t381;
1271   t904 = t7 * t305 * t386;
1272   t907 = t7 * t305 * t390;
1273   t911 = t7 * t584 * t370;
1274   t918 = -t883 / 0.2e1 + t887 - t891 / 0.2e1 - t895 / 0.2e1 - t898 / 0.6e1 - t901 / 0.18e2 - t904 / 0.3e1 - t907 / 0.6e1 + t789 / 0.3e1 - t911 / 0.18e2 + 0.2e1 / 0.9e1 * t796 - t800 / 0.3e1 + 0.2e1 / 0.9e1 * t806 - t814 - t817 + t594 / 0.9e1 + t614 / 0.9e1;
1275   t921 = t573 - t332 + t356 - t361 - t364 / 0.2e1 - t379 - t388 - t393 - t254 - t436 + t9 * (t878 + t918);
1276   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1277     v3rho3[2] = -t234 - t308 - t727 - 0.2e1 / 0.3e1 * t328 - t339 - t728 - 0.2e1 / 0.3e1 * t351 - t373 - t384 + t571 + t921;
1278 
1279   t929 = t796 / 0.3e1;
1280   t930 = t806 / 0.3e1;
1281   t932 = t774 / 0.4e1;
1282   t935 = -0.3e1 * t239 - 0.3e1 * t636;
1283   t936 = t41 * t935;
1284   t938 = t7 * t72 * t936;
1285   t939 = t938 / 0.2e1;
1286   t940 = t171 * t380;
1287   t942 = t7 * t299 * t940;
1288   t943 = t942 / 0.2e1;
1289   t944 = t433 * t140;
1290   t946 = t7 * t122 * t944;
1291   t947 = 0.3e1 / 0.2e1 * t946;
1292   t949 = t171 * t389;
1293   t951 = t7 * t122 * t949;
1294   t952 = 0.3e1 / 0.2e1 * t951;
1295   t953 = -t935;
1296   t954 = t66 * t953;
1297   t956 = t7 * t122 * t954;
1298   t957 = t956 / 0.2e1;
1299   t959 = 0.1e1 / t51 / t206;
1300   t960 = sigma[2] * t959;
1301   t964 = t396 * t152;
1302   t965 = t964 * t168;
1303   t968 = t147 * t406;
1304   t969 = t968 * t408;
1305   t972 = t402 * t430;
1306   t975 = pow(t60, 0.980e2);
1307   t976 = t57 * t975;
1308   t977 = t408 * t168;
1309   t980 = t168 * t430;
1310   t987 = 0.1e1 / t50 / t422;
1311   t992 = t422 * t145;
1312   t993 = 0.1e1 / t992;
1313   t997 = t157 * t157;
1314   t998 = t422 * t206;
1315   t1000 = 0.1e1 / t51 / t998;
1316   t1002 = t163 * t163;
1317   t1003 = 0.1e1 / t1002;
1318   t1005 = t1003 * t18 * t22;
1319   t1008 = -0.1232e4 / 0.27e2 * t23 * t960 * t58 + 0.22988331709569538451e1 * t103 * t157 * t987 * t164 - 0.21541033497890498025e-2 * t421 * t993 * t427 + 0.10318101120731830309e-4 * t997 * t1000 * t1005;
1320   t1011 = -0.41734945439179794323e0 * t23 * t960 * t63 - 0.12263609725822668962e-194 * t401 * t965 + 0.33446208343152733533e-193 * t401 * t969 + 0.33446208343152733533e-195 * t401 * t972 - 0.36382726413893711298e-190 * t976 * t977 - 0.11025068610270821606e-191 * t407 * t980 - 0.36750228700902738685e-194 * t153 * t1008;
1321   t1013 = t7 * t48 * t1011;
1322   t1014 = 0.3e1 / 0.8e1 * t1013;
1323   t1018 = t369 * t135;
1324   t1019 = t41 * t1018;
1325   t1021 = t7 * t578 * t1019;
1326   t1022 = t1021 / 0.9e1;
1327   t1023 = t135 * t375;
1328   t1025 = t321 * t322 * t1023;
1329   t1026 = t1025 / 0.2e1;
1330   t1027 = t380 * t140;
1331   t1028 = t66 * t1027;
1332   t1030 = t7 * t598 * t1028;
1333   t1031 = t1030 / 0.9e1;
1334   t1032 = t140 * t389;
1335   t1034 = t344 * t345 * t1032;
1336   t1035 = t1034 / 0.2e1;
1337   t1036 = -t952 - t957 - t1014 - t898 / 0.2e1 - t901 / 0.6e1 - t904 - t907 / 0.2e1 + t1022 - t1026 + t1031 - t1035;
1338   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1339     v3rho3[3] = -t328 - t351 - t372 / 0.2e1 - t383 / 0.2e1 + t571 + t573 - 0.3e1 / 0.4e1 * t364 - 0.3e1 / 0.2e1 * t378 - 0.3e1 * t388 - 0.3e1 / 0.2e1 * t392 - 0.9e1 / 0.8e1 * t435 + t9 * (t929 + t930 - t911 / 0.6e1 - t646 - t721 + t932 - 0.3e1 / 0.8e1 * t866 - t939 - t943 - t947 + t1036);
1340 
1341   t1040 = t7 * t72 * t443;
1342   t1044 = t7 * t84 * t196 / 0.6e1;
1343   t1046 = t7 * t17 * t475;
1344   t1047 = t196 * t226;
1345   t1051 = t475 * t77;
1346   t1054 = t196 * t241;
1347   t1063 = t28 * t266;
1348   t1070 = t659 * t193;
1349   t1076 = t262 * t472;
1350   t1079 = t193 * t268;
1351   t1082 = t472 * t115;
1352   t1085 = t193 * t290;
1353   t1091 = t276 * t111;
1354   t1095 = t284 * t287;
1355   t1098 = t282 * t105;
1356   t1100 = 0.1e1 / t26 / t1098;
1357   t1104 = 0.88e2 / 0.9e1 * t23 * t256 * t33 - 0.65729100929121114339e0 * t103 * t1091 * sigma[0] + 0.72275836078448381531e-3 * t1095 * t104 - 0.38692879202744363659e-5 * t281 * t1100 * t700;
1358   t1107 = 0.89432025941099559263e-1 * t23 * t256 * t38 + 0.22297472228768489022e-195 * t23 * t263 - 0.41807760428940916917e-194 * t23 * t1063 * t268 - 0.41807760428940916917e-196 * t23 * t450 * t290 - 0.40878699086075563207e-195 * t261 * t1070 + 0.22297472228768489022e-193 * t261 * t663 * t457 + 0.22297472228768489022e-195 * t261 * t1076 - 0.36382726413893711298e-190 * t671 * t1079 - 0.73500457401805477370e-192 * t267 * t1082 - 0.36750228700902738685e-192 * t267 * t1085 - 0.36750228700902738685e-194 * t96 * t1104;
1359   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1360     v3rho2sigma[0] = -0.4e1 / 0.3e1 * t1040 - t1044 - t1046 - t179 * t320 * t1047 / 0.6e1 - t179 * t442 * t1051 - t179 * t442 * t1054 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t180 * t1107;
1361 
1362   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1363     v3rho2sigma[1] = 0.0e0;
1364 
1365   t1112 = t7 * t122 * t482;
1366   t1116 = t7 * t128 * t216 / 0.6e1;
1367   t1117 = t216 * t300;
1368   t1121 = t216 * t309;
1369   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1370     v3rho2sigma[2] = -0.4e1 / 0.3e1 * t1112 - t1116 - t179 * t343 * t1117 / 0.6e1 - t179 * t481 * t1121 / 0.2e1;
1371 
1372   t1127 = t16 * t1;
1373   t1128 = t1127 * t3;
1374   t1129 = t1128 * t487;
1375   t1131 = t3 * t6;
1376   t1132 = t178 * t1131;
1377   t1133 = t224 * t196;
1378   t1134 = t1133 * t323;
1379   t1138 = t442 * t475 * t135;
1380   t1141 = t247 * t1;
1381   t1142 = t1141 * t3;
1382   t1144 = t442 * t196 * t8;
1383   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1384     v3rho2sigma[3] = -0.2e1 / 0.3e1 * t1040 - t1044 - t1046 / 0.2e1 - 0.2e1 / 0.3e1 * t1129 - t1132 * t1134 / 0.6e1 - t179 * t1138 / 0.2e1 - t1142 * t1144 / 0.2e1;
1385 
1386   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1387     v3rho2sigma[4] = 0.0e0;
1388 
1389   t1148 = t1128 * t492;
1390   t1150 = t298 * t216;
1391   t1151 = t1150 * t346;
1392   t1155 = t481 * t216 * t8;
1393   t1158 = t1128 * t524;
1394   t1161 = t481 * t523 * t123;
1395   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1396     v3rho2sigma[5] = -0.2e1 / 0.3e1 * t1112 - t1116 - 0.2e1 / 0.3e1 * t1148 - t1132 * t1151 / 0.6e1 + t1142 * t1155 / 0.2e1 - t1158 / 0.2e1 - t179 * t1161 / 0.2e1;
1397 
1398   t1166 = t320 * t196 * t369;
1399   t1170 = t442 * t196 * t375;
1400   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1401     v3rho2sigma[6] = -0.4e1 / 0.3e1 * t1129 - t1044 - t179 * t1166 / 0.6e1 - t179 * t1170 / 0.2e1;
1402 
1403   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1404     v3rho2sigma[7] = 0.0e0;
1405 
1406   t1175 = t343 * t216 * t380;
1407   t1179 = t481 * t523 * t140;
1408   t1182 = t481 * t216 * t389;
1409   t1190 = t53 * t406;
1410   t1197 = t964 * t213;
1411   t1203 = t402 * t520;
1412   t1206 = t213 * t408;
1413   t1209 = t520 * t168;
1414   t1212 = t213 * t430;
1415   t1218 = t416 * t164;
1416   t1222 = t424 * t427;
1417   t1225 = t422 * t158;
1418   t1227 = 0.1e1 / t51 / t1225;
1419   t1231 = 0.88e2 / 0.9e1 * t23 * t396 * t58 - 0.65729100929121114339e0 * t103 * t1218 * sigma[2] + 0.72275836078448381531e-3 * t1222 * t157 - 0.38692879202744363659e-5 * t421 * t1227 * t1005;
1420   t1234 = 0.89432025941099559263e-1 * t23 * t396 * t63 + 0.22297472228768489022e-195 * t23 * t403 - 0.41807760428940916917e-194 * t23 * t1190 * t408 - 0.41807760428940916917e-196 * t23 * t498 * t430 - 0.40878699086075563207e-195 * t401 * t1197 + 0.22297472228768489022e-193 * t401 * t968 * t505 + 0.22297472228768489022e-195 * t401 * t1203 - 0.36382726413893711298e-190 * t976 * t1206 - 0.73500457401805477370e-192 * t407 * t1209 - 0.36750228700902738685e-192 * t407 * t1212 - 0.36750228700902738685e-194 * t153 * t1231;
1421   t1235 = t200 * t1234;
1422   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1423     v3rho2sigma[8] = -0.4e1 / 0.3e1 * t1148 - t1116 - t1158 - t179 * t1175 / 0.6e1 - t179 * t1179 - t179 * t1182 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t1235;
1424 
1425   t1238 = t1128 * t544;
1426   t1241 = t442 * t543 * t77;
1427   t1246 = t23 * t28;
1428   t1247 = t266 * t193;
1429   t1248 = t1247 * t115;
1430   t1254 = t663 * t530;
1431   t1257 = t530 * t115;
1432   t1260 = t193 * t472;
1433   t1263 = t262 * t540;
1434   t1266 = t540 * t115;
1435   t1271 = t468 * t287;
1436   t1275 = 0.1e1 / t26 / t687;
1437   t1279 = 0.12134603248453128801e0 * t103 * t463 - 0.21257598846602465156e-3 * t1271 * sigma[0] + 0.14509829701029136372e-5 * t104 * t1275 * t700;
1438   t1282 = 0.22297472228768489022e-195 * t23 * t454 - 0.83615520857881833834e-194 * t1246 * t1248 - 0.83615520857881833834e-196 * t23 * t450 * t472 + 0.11148736114384244511e-193 * t261 * t1254 - 0.36382726413893711298e-190 * t671 * t1257 - 0.73500457401805477370e-192 * t267 * t1260 + 0.11148736114384244511e-195 * t261 * t1263 - 0.36750228700902738685e-192 * t267 * t1266 - 0.36750228700902738685e-194 * t96 * t1279;
1439   t1283 = t180 * t1282;
1440   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1441     v3rhosigma2[0] = -t1238 / 0.2e1 - t179 * t1241 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t1283;
1442 
1443   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1444     v3rhosigma2[1] = 0.0e0;
1445 
1446   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1447     v3rhosigma2[2] = 0.0e0;
1448 
1449   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1450     v3rhosigma2[3] = 0.0e0;
1451 
1452   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1453     v3rhosigma2[4] = 0.0e0;
1454 
1455   t1286 = t1128 * t564;
1456   t1288 = t481 * t563 * t123;
1457   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1458     v3rhosigma2[5] = -t179 * t1288 / 0.2e1 - t1286 / 0.2e1;
1459 
1460   t1292 = t442 * t543 * t135;
1461   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1462     v3rhosigma2[6] = -t179 * t1292 / 0.2e1 - t1238 / 0.2e1;
1463 
1464   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1465     v3rhosigma2[7] = 0.0e0;
1466 
1467   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1468     v3rhosigma2[8] = 0.0e0;
1469 
1470   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1471     v3rhosigma2[9] = 0.0e0;
1472 
1473   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1474     v3rhosigma2[10] = 0.0e0;
1475 
1476   t1297 = t481 * t563 * t140;
1477   t1302 = t23 * t53;
1478   t1303 = t406 * t213;
1479   t1304 = t1303 * t168;
1480   t1310 = t968 * t550;
1481   t1313 = t550 * t168;
1482   t1316 = t213 * t520;
1483   t1319 = t402 * t560;
1484   t1322 = t560 * t168;
1485   t1327 = t516 * t427;
1486   t1331 = 0.1e1 / t51 / t992;
1487   t1335 = 0.12134603248453128801e0 * t103 * t511 - 0.21257598846602465156e-3 * t1327 * sigma[2] + 0.14509829701029136372e-5 * t157 * t1331 * t1005;
1488   t1338 = 0.22297472228768489022e-195 * t23 * t502 - 0.83615520857881833834e-194 * t1302 * t1304 - 0.83615520857881833834e-196 * t23 * t498 * t520 + 0.11148736114384244511e-193 * t401 * t1310 - 0.36382726413893711298e-190 * t976 * t1313 - 0.73500457401805477370e-192 * t407 * t1316 + 0.11148736114384244511e-195 * t401 * t1319 - 0.36750228700902738685e-192 * t407 * t1322 - 0.36750228700902738685e-194 * t153 * t1335;
1489   t1339 = t200 * t1338;
1490   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1491     v3rhosigma2[11] = -t1286 / 0.2e1 - t179 * t1297 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t1339;
1492 
1493   t1348 = t530 * t193;
1494   t1351 = t193 * t540;
1495   t1357 = 0.1e1 / t26 / t283;
1496   t1361 = 0.47829597404855546601e-4 * t536 * t287 - 0.54411861378859261395e-6 * sigma[0] * t1357 * t700;
1497   t1364 = -0.12542328128682275075e-193 * t23 * t1063 * t530 - 0.12542328128682275075e-195 * t23 * t450 * t540 - 0.36382726413893711298e-190 * t671 * t1348 - 0.11025068610270821606e-191 * t267 * t1351 - 0.36750228700902738685e-194 * t96 * t1361;
1498   t1365 = t180 * t1364;
1499   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1500     v3sigma3[0] = -0.3e1 / 0.8e1 * t179 * t1365;
1501 
1502   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1503     v3sigma3[1] = 0.0e0;
1504 
1505   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1506     v3sigma3[2] = 0.0e0;
1507 
1508   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1509     v3sigma3[3] = 0.0e0;
1510 
1511   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1512     v3sigma3[4] = 0.0e0;
1513 
1514   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1515     v3sigma3[5] = 0.0e0;
1516 
1517   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1518     v3sigma3[6] = 0.0e0;
1519 
1520   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1521     v3sigma3[7] = 0.0e0;
1522 
1523   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1524     v3sigma3[8] = 0.0e0;
1525 
1526   t1374 = t550 * t213;
1527   t1377 = t213 * t560;
1528   t1383 = 0.1e1 / t51 / t423;
1529   t1387 = 0.47829597404855546601e-4 * t556 * t427 - 0.54411861378859261395e-6 * sigma[2] * t1383 * t1005;
1530   t1390 = -0.12542328128682275075e-193 * t23 * t1190 * t550 - 0.12542328128682275075e-195 * t23 * t498 * t560 - 0.36382726413893711298e-190 * t976 * t1374 - 0.11025068610270821606e-191 * t407 * t1377 - 0.36750228700902738685e-194 * t153 * t1387;
1531   t1391 = t200 * t1390;
1532   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1533     v3sigma3[9] = -0.3e1 / 0.8e1 * t179 * t1391;
1534 
1535 #ifndef XC_DONT_COMPILE_LXC
1536 
1537   if(order < 4) return;
1538 
1539 
1540   t1394 = t634 * t9;
1541   t1396 = t8 / t1394;
1542   t1398 = 0.12e2 * t635 - 0.12e2 * t1396;
1543   t1406 = 0.10e2 / 0.27e2 * t7 * t15 * t786 * t41;
1544   t1410 = 0.10e2 / 0.27e2 * t7 * t47 * t786 * t66;
1545   t1413 = sigma[0] / t26 / t106;
1546   t1417 = t654 * t95;
1547   t1421 = t256 * t266;
1548   t1428 = t90 * t670;
1549   t1438 = pow(t35, 0.970e2);
1550   t1439 = t32 * t1438;
1551   t1440 = t268 * t268;
1552   t1446 = t290 * t290;
1553   t1465 = t282 * t106;
1554   t1472 = t282 * t282;
1555   t1480 = 0.1e1 / t697 / t32 * t100 * t102;
1556   t1486 = 0.23649802415535216783e1 * t23 * t1413 * t38 + 0.76306904960674384652e-194 * t261 * t1417 * t115 - 0.24527219451645337924e-192 * t261 * t1421 * t268 - 0.24527219451645337924e-194 * t261 * t659 * t290 + 0.44148995012961608264e-191 * t261 * t1428 * t672 + 0.13378483337261093413e-192 * t261 * t663 * t675 + 0.44594944457536978044e-195 * t261 * t262 * t703 - 0.35655071885615837072e-188 * t1439 * t1440 - 0.21829635848336226779e-189 * t671 * t268 * t290 - 0.11025068610270821606e-191 * t267 * t1446 - 0.14700091480361095474e-191 * t267 * t115 * t703 - 0.36750228700902738685e-194 * t96 * (0.20944e5 / 0.81e2 * t23 * t1413 * t33 - 0.20541186721131490809e2 * t103 * t104 / t25 / t467 * t111 + 0.32286356055610144116e-1 * t281 / t1098 * t287 - 0.33705796994390645677e-3 * t692 / t26 / t1465 * t700 + 0.12520606337750033699e-5 * t692 * sigma[0] / t25 / t1472 / rho[0] * t1480);
1557   t1491 = t7 * t84 * t706;
1558   t1494 = t7 * t643 * t118;
1559   t1497 = t7 * t248 * t293;
1560   t1508 = t321 * t737 * t588;
1561   t1522 = t344 * t752 * t608;
1562   t1524 = t13 * t13;
1563   t1526 = 0.1e1 / t223 / t1524;
1564   t1527 = t1526 * t16;
1565   t1528 = t226 * t226;
1566   t1533 = t577 * t83;
1567   t1535 = t7 * t1533 * t580;
1568   t1537 = t224 * t247;
1569   t1539 = t7 * t1537 * t227;
1570   t1541 = t241 * t241;
1571   t1547 = t7 * t336 * t78;
1572   t1549 = -t7 * t122 * t66 * t1398 / 0.2e1 + t1406 + t1410 - 0.3e1 / 0.8e1 * t7 * t17 * t1486 - t1491 / 0.2e1 - 0.5e1 / 0.9e1 * t1494 + t1497 / 0.2e1 - 0.2e1 / 0.3e1 * t344 * t345 * t123 * t713 + 0.2e1 / 0.3e1 * t732 * t322 * t226 * t241 - 0.2e1 / 0.3e1 * t1508 - 0.2e1 * t321 * t756 * t588 - 0.2e1 / 0.3e1 * t321 * t322 * t77 * t638 + 0.2e1 / 0.3e1 * t747 * t345 * t300 * t309 - 0.2e1 / 0.3e1 * t1522 - 0.5e1 / 0.27e2 * t7 * t1527 * t41 * t1528 + 0.4e1 / 0.27e2 * t1535 + 0.2e1 / 0.9e1 * t1539 - t7 * t225 * t41 * t1541 / 0.2e1 - 0.20e2 / 0.27e2 * t1547;
1573   t1550 = t45 * t45;
1574   t1552 = 0.1e1 / t297 / t1550;
1575   t1553 = t1552 * t16;
1576   t1554 = t300 * t300;
1577   t1559 = t597 * t83;
1578   t1561 = t7 * t1559 * t600;
1579   t1563 = t298 * t247;
1580   t1565 = t7 * t1563 * t301;
1581   t1567 = t309 * t309;
1582   t1573 = t7 * t353 * t124;
1583   t1580 = t7 * t584 * t616;
1584   t1583 = t7 * t592 * t235;
1585   t1586 = t7 * t592 * t242;
1586   t1589 = t7 * t612 * t310;
1587   t1595 = t7 * t231 * t626;
1588   t1598 = t7 * t231 * t630;
1589   t1601 = t7 * t231 * t639;
1590   t1621 = t7 * t305 * t714;
1591   t1623 = -0.5e1 / 0.27e2 * t7 * t1553 * t66 * t1554 + 0.4e1 / 0.27e2 * t1561 + 0.2e1 / 0.9e1 * t1565 - t7 * t299 * t66 * t1567 / 0.2e1 - 0.20e2 / 0.27e2 * t1573 + 0.4e1 / 0.9e1 * t7 * t578 * t118 * t579 - 0.2e1 / 0.3e1 * t1580 + 0.4e1 / 0.3e1 * t1583 + 0.2e1 / 0.3e1 * t1586 + 0.2e1 / 0.3e1 * t1589 - t7 * t225 * t293 * t226 - 0.2e1 * t1595 - 0.2e1 * t1598 - 0.2e1 / 0.3e1 * t1601 - 0.2e1 * t7 * t72 * t706 * t77 - 0.3e1 * t7 * t72 * t293 * t241 - 0.2e1 * t7 * t72 * t118 * t638 + t7 * t72 * t41 * t1398 / 0.2e1 - 0.2e1 / 0.3e1 * t1621;
1592   t1636 = t9 * (t1549 + t1623) + 0.4e1 / 0.9e1 * t582 - 0.2e1 / 0.3e1 * t586 - 0.2e1 * t590 + 0.4e1 / 0.3e1 * t594 + 0.4e1 / 0.9e1 * t602 - 0.2e1 / 0.3e1 * t606 - 0.2e1 * t610 + 0.4e1 / 0.3e1 * t614 - 0.2e1 * t618 - 0.4e1 * t621;
1593   t1641 = 0.5e1 / 0.9e1 * t645;
1594   t1646 = 0.5e1 / 0.9e1 * t720;
1595   t1647 = -0.2e1 * t624 - 0.6e1 * t628 - 0.6e1 * t632 - 0.2e1 * t641 - t1641 + t648 - 0.3e1 / 0.2e1 * t651 - 0.3e1 / 0.2e1 * t708 - 0.2e1 * t711 - 0.2e1 * t716 - t1646;
1596   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1597     v4rho4[0] = t1636 + t1647;
1598 
1599   t1648 = t786 * t41;
1600   t1650 = t321 * t1648 * t742;
1601   t1662 = t344 * t247 * t66 * t346;
1602   t1663 = t1662 / 0.9e1;
1603   t1671 = t786 * t66;
1604   t1673 = t344 * t1671 * t769;
1605   t1676 = t4 * t6 * t1526;
1606   t1682 = t732 * t737 * t733;
1607   t1689 = t4 * t6 * t1552;
1608   t1695 = t747 * t752 * t748;
1609   t1702 = 0.1e1 / t82 / t634;
1610   t1705 = t7 * t46 * t1702 * t354;
1611   t1708 = t7 * t305 * t791;
1612   t1711 = t7 * t592 * t330;
1613   t1712 = t1711 / 0.3e1;
1614   t1713 = 0.7e1 / 0.6e1 * t1650 - t344 * t345 * t140 * t713 / 0.6e1 + t344 * t768 * t8 * t309 / 0.2e1 + t1663 - t344 * t857 * t608 / 0.2e1 - t747 * t768 * t8 * t300 / 0.3e1 - 0.7e1 / 0.6e1 * t1673 - 0.5e1 / 0.27e2 * t1676 * t322 * t135 * t579 + t1682 / 0.9e1 + t732 * t741 * t8 * t226 / 0.3e1 - 0.5e1 / 0.27e2 * t1689 * t345 * t140 * t599 + t1695 / 0.9e1 - t7 * t72 * t706 * t135 / 0.2e1 + 0.25e2 / 0.6e1 * t1705 - t1708 / 0.2e1 + t1712;
1615   t1715 = t7 * t798 * t802;
1616   t1722 = t7 * t336 * t136;
1617   t1726 = t7 * t14 * t1702 * t337;
1618   t1728 = t298 * t335;
1619   t1733 = t7 * t353 * t141;
1620   t1740 = t7 * t604 * t808;
1621   t1743 = t7 * t612 * t359;
1622   t1744 = t1743 / 0.3e1;
1623   t1745 = t224 * t335;
1624   t1754 = t7 * t231 * t818;
1625   t1757 = t7 * t798 * t41;
1626   t1763 = t7 * t718 * t171;
1627   t1766 = t7 * t787 * t66;
1628   t1768 = 0.7e1 / 0.2e1 * t1715 - 0.3e1 / 0.2e1 * t7 * t336 * t293 * t8 - 0.5e1 / 0.27e2 * t1722 - 0.25e2 / 0.6e1 * t1726 + t7 * t1728 * t124 / 0.2e1 - 0.5e1 / 0.27e2 * t1733 + t7 * t598 * t171 * t599 / 0.9e1 - t1740 / 0.6e1 + t1744 - t7 * t1745 * t78 / 0.2e1 - t7 * t122 * t171 * t713 / 0.2e1 - t1754 / 0.2e1 + 0.5e1 / 0.2e1 * t1757 - 0.3e1 / 0.2e1 * t7 * t336 * t118 - 0.5e1 / 0.36e2 * t1763 - 0.5e1 / 0.2e1 * t1766;
1629   t1778 = t321 * t83 * t118 * t323;
1630   t1780 = t335 * t118;
1631   t1784 = t321 * t737 * t760;
1632   t1791 = t344 * t752 * t764;
1633   t1806 = t321 * t247 * t41 * t323;
1634   t1807 = t1806 / 0.9e1;
1635   t1809 = t1406 + t1410 - t1491 / 0.8e1 - 0.5e1 / 0.12e2 * t1494 + t1497 / 0.4e1 + t732 * t756 * t733 / 0.3e1 - t1778 / 0.3e1 - t321 * t1780 * t742 - t1784 / 0.6e1 - t321 * t741 * t8 * t241 / 0.2e1 - t1791 / 0.6e1 - t321 * t16 * t293 * t323 / 0.2e1 - t321 * t756 * t760 / 0.2e1 - t321 * t322 * t135 * t638 / 0.6e1 + t1807 - t1508 / 0.6e1;
1636   t1833 = -t1522 / 0.6e1 + t1535 / 0.27e2 + t1539 / 0.9e1 - 0.5e1 / 0.9e1 * t1547 + t1561 / 0.27e2 + t1565 / 0.9e1 - 0.5e1 / 0.9e1 * t1573 - t1580 / 0.6e1 + 0.2e1 / 0.3e1 * t1583 + t1586 / 0.3e1 + t1589 / 0.3e1 - t1595 / 0.2e1 - t1598 / 0.2e1 - t1601 / 0.6e1 - t1621 / 0.6e1 + t732 * t322 * t760 * t77 / 0.3e1 + t747 * t345 * t764 * t123 / 0.3e1;
1637   t1838 = t739 / 0.3e1;
1638   t1840 = t754 / 0.3e1;
1639   t1843 = t9 * (t1713 + t1768 + t1809 + t1833) + t735 / 0.3e1 - t1838 - t744 + t750 / 0.3e1 - t1840 - t758 - t762 / 0.2e1 - t766 / 0.2e1 + t771 - t591;
1640   t1850 = -t611 + t932 - 0.3e1 / 0.2e1 * t777 + 0.3e1 / 0.2e1 * t780 - t1646 - t1641 + 0.3e1 / 0.4e1 * t648 - 0.3e1 / 0.4e1 * t651 - t709 - 0.7e1 / 0.2e1 * t789 - 0.3e1 / 0.2e1 * t793;
1641   t1857 = t929 + 0.7e1 / 0.2e1 * t800 - 0.3e1 * t804 + t930 - t810 / 0.2e1 - t813 - t816 - 0.3e1 / 0.2e1 * t820 + t583 - t586 / 0.3e1 + t594;
1642   t1860 = t603 - t606 / 0.3e1 + t614 - t619 - 0.2e1 * t621 - t624 - t629 - t633 - t642 - t711 - t717;
1643   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1644     v4rho4[1] = t1843 + t1850 + t1857 + t1860;
1645 
1646   t1875 = 0.2e1 / 0.9e1 * t839 - 0.2e1 / 0.3e1 * t843 - t847 / 0.3e1 + 0.2e1 / 0.9e1 * t851 + 0.2e1 / 0.3e1 * t855 - 0.2e1 / 0.3e1 * t859 - t863 / 0.3e1 + 0.2e1 / 0.9e1 * t735 - 0.4e1 / 0.9e1 * t739 - 0.2e1 / 0.3e1 * t744 + 0.2e1 / 0.9e1 * t750 - 0.4e1 / 0.9e1 * t754 - 0.2e1 / 0.3e1 * t758;
1647   t1884 = -t762 / 0.3e1 - t766 / 0.3e1 + 0.2e1 / 0.3e1 * t771 - t866 / 0.4e1 + t774 / 0.2e1 - t777 + t780 - t1646 - t1641 + t648 / 0.2e1 - t651 / 0.4e1 - t872 / 0.3e1 - t876;
1648   t1895 = -t883 + 0.2e1 * t887 - t891 - t895 - t898 / 0.3e1 - t901 / 0.9e1 - 0.2e1 / 0.3e1 * t904 - t907 / 0.3e1 - 0.5e1 / 0.3e1 * t789 - t793 - t911 / 0.9e1 + 0.2e1 / 0.3e1 * t796 + 0.5e1 / 0.3e1 * t800;
1649   t1911 = t7 * t314 * t433;
1650   t1918 = t344 * t83 * t171 * t346;
1651   t1925 = t344 * t752 * t861;
1652   t1932 = t321 * t737 * t845;
1653   t1939 = t747 * t752 * t849;
1654   t1946 = t732 * t737 * t837;
1655   t1949 = t321 * t1648 * t841;
1656   t1956 = t344 * t1671 * t853;
1657   t1969 = t335 * t171;
1658   t1977 = t7 * t353 * t171 + t1911 / 0.12e2 + 0.2e1 / 0.9e1 * t747 * t857 * t748 - 0.2e1 / 0.9e1 * t1918 + t747 * t345 * t389 * t300 / 0.9e1 - t1925 / 0.9e1 + t732 * t322 * t375 * t226 / 0.9e1 - t1932 / 0.9e1 - 0.5e1 / 0.27e2 * t1689 * t345 * t380 * t300 + 0.2e1 / 0.27e2 * t1939 - 0.5e1 / 0.27e2 * t1676 * t322 * t369 * t226 + 0.2e1 / 0.27e2 * t1946 + 0.7e1 / 0.9e1 * t1949 + t747 * t345 * t380 * t309 / 0.9e1 - 0.7e1 / 0.9e1 * t1956 - t344 * t857 * t764 / 0.3e1 - t344 * t345 * t389 * t309 / 0.6e1 - t321 * t322 * t880 * t77 / 0.3e1 + 0.2e1 / 0.3e1 * t344 * t1969 * t769 - t344 * t345 * t888 * t123 / 0.3e1;
1659   t2016 = t7 * t584 * t870;
1660   t2019 = 0.1e1 / t82 / t1394;
1661   t2021 = t8 * t8;
1662   t2027 = t7 * t612 * t390;
1663   t2030 = t7 * t592 * t376;
1664   t2033 = t7 * t1563 * t381;
1665   t2038 = -t321 * t322 * t375 * t241 / 0.6e1 - t321 * t756 * t845 / 0.3e1 + 0.2e1 / 0.9e1 * t732 * t756 * t837 + t732 * t322 * t369 * t241 / 0.9e1 - 0.2e1 / 0.3e1 * t321 * t1780 * t841 - 0.2e1 / 0.9e1 * t1650 + 0.4e1 / 0.27e2 * t1662 + 0.2e1 / 0.9e1 * t1673 + 0.2e1 / 0.27e2 * t1682 + 0.2e1 / 0.27e2 * t1695 - t7 * t72 * t293 * t375 / 0.2e1 - t7 * t72 * t118 * t880 - t7 * t122 * t433 * t309 / 0.2e1 - t7 * t225 * t293 * t369 / 0.6e1 - t2016 / 0.9e1 - t7 * t224 * t2019 * t41 * t2021 / 0.3e1 + t2027 / 0.9e1 + t2030 / 0.9e1 + t2033 / 0.27e2 + t7 * t1728 * t141 / 0.3e1;
1666   t2041 = t7 * t612 * t386;
1667   t2044 = t7 * t231 * t881;
1668   t2047 = t7 * t787 * t885;
1669   t2050 = t7 * t305 * t889;
1670   t2057 = t7 * t305 * t893;
1671   t2063 = t7 * t231 * t874;
1672   t2066 = t7 * t1537 * t370;
1673   t2083 = 0.2e1 / 0.9e1 * t2041 - t2044 / 0.3e1 - 0.7e1 / 0.3e1 * t2047 - t2050 / 0.3e1 - t7 * t299 * t433 * t300 / 0.6e1 - t2057 / 0.3e1 - t7 * t1745 * t136 / 0.3e1 - t2063 / 0.3e1 + t2066 / 0.27e2 - t7 * t298 * t2019 * t66 * t2021 / 0.3e1 + 0.41e2 / 0.9e1 * t1705 - t1708 / 0.3e1 + 0.4e1 / 0.9e1 * t1711 - 0.2e1 / 0.3e1 * t1715 - 0.10e2 / 0.27e2 * t1722 - 0.41e2 / 0.9e1 * t1726 - 0.10e2 / 0.27e2 * t1733 - t1740 / 0.9e1 + 0.4e1 / 0.9e1 * t1743 - t1754 / 0.3e1;
1674   t2108 = -t1757 / 0.3e1 - 0.5e1 / 0.18e2 * t1763 + t1766 / 0.3e1 + t1406 + t1410 - 0.5e1 / 0.18e2 * t1494 + t1497 / 0.12e2 - 0.2e1 / 0.9e1 * t1778 - t1784 / 0.9e1 - t1791 / 0.9e1 + 0.4e1 / 0.27e2 * t1806 + t1539 / 0.27e2 - 0.10e2 / 0.27e2 * t1547 + t1565 / 0.27e2 - 0.10e2 / 0.27e2 * t1573 + 0.2e1 / 0.9e1 * t1583 + t1586 / 0.9e1 + t1589 / 0.9e1 - 0.4e1 / 0.9e1 * t747 * t768 * t346 * t8 + 0.4e1 / 0.9e1 * t732 * t741 * t323 * t8;
1675   t2112 = -0.2e1 * t804 + 0.2e1 / 0.3e1 * t806 - t810 / 0.3e1 - 0.4e1 / 0.3e1 * t813 - 0.4e1 / 0.3e1 * t816 - t820 - t586 / 0.9e1 + 0.2e1 / 0.3e1 * t594 - t606 / 0.9e1 + 0.2e1 / 0.3e1 * t614 - 0.2e1 / 0.3e1 * t621 - t624 / 0.3e1 - t711 / 0.3e1 + t9 * (t1977 + t2038 + t2083 + t2108);
1676   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1677     v4rho4[2] = t1875 + t1884 + t1895 + t2112;
1678 
1679   t2118 = -t1014 - t1026 - t1035 + t839 / 0.3e1 - t843 - t847 / 0.2e1 + t851 / 0.3e1 + t855 - t859 - t863 / 0.2e1 - t1838;
1680   t2121 = -t1840 - 0.3e1 / 0.4e1 * t866 + 0.3e1 / 0.4e1 * t774 - t1646 - t1641 + t649 - t939 - t943 - t947 - t952 - t957;
1681   t2131 = t1022 + t1031 - t872 / 0.2e1 - 0.3e1 / 0.2e1 * t876 - 0.3e1 / 0.2e1 * t883 + 0.3e1 * t887 - 0.3e1 / 0.2e1 * t891 - 0.3e1 / 0.2e1 * t895 - t898 - t901 / 0.3e1 - 0.2e1 * t904;
1682   t2134 = t7 * t128 * t1011;
1683   t2138 = t344 * t752 * t1032;
1684   t2145 = t321 * t737 * t1023;
1685   t2187 = -t2134 / 0.8e1 + t1911 / 0.4e1 - t2138 / 0.6e1 + t732 * t741 * t369 * t8 / 0.3e1 - t2145 / 0.6e1 - t321 * t741 * t8 * t375 / 0.2e1 - 0.5e1 / 0.27e2 * t1689 * t345 * t1027 * t123 - t747 * t768 * t380 * t8 / 0.3e1 + t344 * t1969 * t853 - t344 * t16 * t433 * t346 / 0.2e1 - t344 * t857 * t861 / 0.2e1 - t344 * t345 * t953 * t123 / 0.6e1 - t321 * t322 * t935 * t77 / 0.6e1 + t747 * t857 * t849 / 0.3e1 - t344 * t345 * t140 * t888 / 0.2e1 + t344 * t768 * t8 * t389 / 0.2e1;
1686   t2216 = 0.6e1 * t635 + 0.12e2 * t1396;
1687   t2222 = t7 * t1533 * t1019;
1688   t2225 = t7 * t1559 * t1028;
1689   t2228 = t7 * t305 * t944;
1690   t2230 = -0.5e1 / 0.27e2 * t1676 * t322 * t1018 * t77 - t321 * t756 * t1023 / 0.2e1 - t321 * t322 * t135 * t880 / 0.2e1 - t1918 / 0.3e1 - t1925 / 0.6e1 - t1932 / 0.6e1 + t1939 / 0.9e1 + t1946 / 0.9e1 - t1949 / 0.3e1 + t1956 / 0.3e1 + t1663 + t7 * t578 * t118 * t1018 / 0.9e1 - t7 * t72 * t118 * t935 / 0.2e1 - t7 * t72 * t41 * t2216 / 0.2e1 + t2222 / 0.27e2 + t2225 / 0.27e2 - t2228 / 0.2e1;
1691   t2233 = t7 * t305 * t949;
1692   t2236 = t7 * t305 * t954;
1693   t2239 = t7 * t231 * t936;
1694   t2242 = t7 * t604 * t940;
1695   t2269 = -t2233 / 0.2e1 - t2236 / 0.6e1 - t2239 / 0.6e1 - t2242 / 0.6e1 + 0.3e1 / 0.2e1 * t7 * t353 * t433 * t8 - 0.3e1 / 0.2e1 * t7 * t122 * t171 * t888 + t7 * t122 * t66 * t2216 / 0.2e1 - t7 * t122 * t1011 * t123 / 0.2e1 - t2016 / 0.6e1 + t2027 / 0.3e1 + t2030 / 0.3e1 + t2033 / 0.9e1 + 0.2e1 / 0.3e1 * t2041 - t2044 / 0.2e1 + t2047 - t2050 / 0.2e1 - t2057 / 0.2e1;
1696   t2288 = -t2063 / 0.2e1 + t2066 / 0.9e1 - t1705 / 0.3e1 + t1712 - 0.5e1 / 0.9e1 * t1722 + t1726 / 0.3e1 - 0.5e1 / 0.9e1 * t1733 + t1744 - 0.5e1 / 0.12e2 * t1763 + t1406 + t1410 - 0.5e1 / 0.36e2 * t1494 + t1807 - 0.5e1 / 0.27e2 * t1547 - 0.5e1 / 0.27e2 * t1573 + t747 * t345 * t1032 * t123 / 0.3e1 + t732 * t322 * t1023 * t77 / 0.3e1;
1697   t2292 = -t907 + t789 - t911 / 0.3e1 + t796 - t800 + t806 - t813 - t816 + t595 + t615 + t9 * (t2187 + t2230 + t2269 + t2288);
1698   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1699     v4rho4[3] = t2118 + t2121 + t2131 + t2292;
1700 
1701   t2296 = sigma[2] / t51 / t159;
1702   t2300 = t959 * t152;
1703   t2304 = t396 * t406;
1704   t2311 = t147 * t975;
1705   t2321 = pow(t60, 0.970e2);
1706   t2322 = t57 * t2321;
1707   t2323 = t408 * t408;
1708   t2329 = t430 * t430;
1709   t2348 = t422 * t159;
1710   t2355 = t422 * t422;
1711   t2363 = 0.1e1 / t1002 / t57 * t100 * t102;
1712   t2369 = 0.23649802415535216783e1 * t23 * t2296 * t63 + 0.76306904960674384652e-194 * t401 * t2300 * t168 - 0.24527219451645337924e-192 * t401 * t2304 * t408 - 0.24527219451645337924e-194 * t401 * t964 * t430 + 0.44148995012961608264e-191 * t401 * t2311 * t977 + 0.13378483337261093413e-192 * t401 * t968 * t980 + 0.44594944457536978044e-195 * t401 * t402 * t1008 - 0.35655071885615837072e-188 * t2322 * t2323 - 0.21829635848336226779e-189 * t976 * t408 * t430 - 0.11025068610270821606e-191 * t407 * t2329 - 0.14700091480361095474e-191 * t407 * t168 * t1008 - 0.36750228700902738685e-194 * t153 * (0.20944e5 / 0.81e2 * t23 * t2296 * t58 - 0.20541186721131490809e2 * t103 * t157 / t50 / t515 * t164 + 0.32286356055610144116e-1 * t421 / t1225 * t427 - 0.33705796994390645677e-3 * t997 / t51 / t2348 * t1005 + 0.12520606337750033699e-5 * t997 * sigma[2] / t50 / t2355 / rho[1] * t2363);
1713   t2396 = t389 * t389;
1714   t2401 = t369 * t369;
1715   t2406 = t375 * t375;
1716   t2411 = t380 * t380;
1717   t2429 = -0.12e2 * t635 - 0.12e2 * t1396;
1718   t2438 = -0.3e1 / 0.8e1 * t7 * t48 * t2369 - t2134 / 0.2e1 + t1911 / 0.2e1 + 0.2e1 / 0.3e1 * t732 * t322 * t369 * t375 - 0.2e1 / 0.3e1 * t344 * t345 * t953 * t140 - 0.2e1 / 0.3e1 * t321 * t322 * t935 * t135 - 0.2e1 * t344 * t857 * t1032 + 0.2e1 / 0.3e1 * t747 * t345 * t380 * t389 - 0.2e1 / 0.3e1 * t2138 - 0.2e1 / 0.3e1 * t2145 - t7 * t299 * t66 * t2396 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1527 * t41 * t2401 - t7 * t225 * t41 * t2406 / 0.2e1 - 0.5e1 / 0.27e2 * t7 * t1553 * t66 * t2411 + 0.4e1 / 0.9e1 * t7 * t598 * t171 * t1027 - 0.3e1 * t7 * t122 * t433 * t389 - 0.2e1 * t7 * t122 * t171 * t953 - t7 * t122 * t66 * t2429 / 0.2e1 - 0.2e1 * t7 * t122 * t1011 * t140;
1719   t2462 = t7 * t72 * t41 * t2429 / 0.2e1 - t7 * t299 * t433 * t380 + 0.4e1 / 0.27e2 * t2222 + 0.4e1 / 0.27e2 * t2225 - 0.2e1 * t2228 - 0.2e1 * t2233 - 0.2e1 / 0.3e1 * t2236 - 0.2e1 / 0.3e1 * t2239 - 0.2e1 / 0.3e1 * t2242 + 0.2e1 / 0.3e1 * t2027 + 0.2e1 / 0.3e1 * t2030 + 0.2e1 / 0.9e1 * t2033 + 0.4e1 / 0.3e1 * t2041 + 0.2e1 / 0.9e1 * t2066 - 0.20e2 / 0.27e2 * t1722 - 0.20e2 / 0.27e2 * t1733 - 0.5e1 / 0.9e1 * t1763 + t1406 + t1410;
1720   t2472 = t9 * (t2438 + t2462) + 0.4e1 / 0.3e1 * t796 + 0.4e1 / 0.3e1 * t806 - 0.2e1 / 0.3e1 * t911 - t1641 - t1646 + t774 - 0.3e1 / 0.2e1 * t866 - 0.2e1 * t938 - 0.2e1 * t942 - 0.6e1 * t946;
1721   t2484 = -0.6e1 * t951 - 0.2e1 * t956 - 0.3e1 / 0.2e1 * t1013 - 0.2e1 * t898 - 0.2e1 / 0.3e1 * t901 - 0.4e1 * t904 - 0.2e1 * t907 + 0.4e1 / 0.9e1 * t1021 - 0.2e1 * t1025 + 0.4e1 / 0.9e1 * t1030 - 0.2e1 * t1034;
1722   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1723     v4rho4[4] = t2472 + t2484;
1724 
1725   t2514 = t28 * t670;
1726   t2568 = -0.10914817924168113389e-189 * t671 * t457 * t290 - 0.35655071885615837072e-188 * t1439 * t193 * t672 - 0.41734945439179794323e0 * t23 * t654 * t38 - 0.12263609725822668962e-194 * t23 * t660 + 0.33446208343152733533e-195 * t23 * t667 - 0.41807760428940916917e-196 * t23 * t450 * t703 + 0.33111746259721206198e-191 * t261 * t1428 * t1079 - 0.10914817924168113389e-189 * t671 * t472 * t268 - 0.11025068610270821606e-191 * t267 * t1104 * t115 - 0.11025068610270821606e-191 * t267 * t472 * t290 - 0.36750228700902738685e-192 * t267 * t193 * t703;
1727   t2590 = t7 * t72 * t1051;
1728   t2593 = t7 * t72 * t1054;
1729   t2597 = t7 * t248 * t196 / 0.9e1;
1730   t2599 = t7 * t17 * t1107;
1731   t2602 = t7 * t84 * t475;
1732   t2612 = t7 * t225 * t1047;
1733   t2615 = t7 * t231 * t443;
1734   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1735     v4rho3sigma[0] = -0.3e1 / 0.8e1 * t179 * t180 * (-0.36750228700902738685e-194 * t96 * (-0.1232e4 / 0.27e2 * t23 * t654 * t33 + 0.51167577030977359778e1 * t103 * t682 * t111 * sigma[0] - 0.96840172523411230156e-2 * t688 * t287 * t104 + 0.11478887496814161219e-3 * t695 * t698 * t281 * t18 * t22 - 0.46952273766562626373e-6 * t692 / t25 / t1472 * t1480) + 0.33446208343152733533e-195 * t261 * t262 * t1104 - 0.12263609725822668962e-194 * t261 * t659 * t472 - 0.41389682824651507748e-192 * t23 * t2514 * t672 - 0.12542328128682275075e-193 * t1246 * t266 * t115 * t290 + 0.33446208343152733533e-193 * t23 * t664 + 0.19076726240168596163e-194 * t261 * t1417 * t193 - 0.12263609725822668962e-192 * t261 * t1421 * t457 + 0.66892416686305467067e-193 * t261 * t663 * t1082 + 0.33446208343152733533e-193 * t261 * t663 * t1085 + t2568) - 0.3e1 / 0.2e1 * t179 * t442 * t475 * t241 - t179 * t442 * t196 * t638 / 0.2e1 - 0.3e1 / 0.2e1 * t179 * t442 * t1107 * t77 - t179 * t320 * t475 * t226 / 0.2e1 - 0.4e1 * t2590 - 0.2e1 * t2593 + t2597 - 0.3e1 / 0.2e1 * t2599 - t2602 / 0.2e1 + t179 * t731 * t196 * t579 / 0.9e1 - t1132 * t1133 * t588 / 0.2e1 - 0.2e1 / 0.3e1 * t2612 - 0.2e1 / 0.3e1 * t2615;
1736 
1737   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1738     v4rho3sigma[1] = 0.0e0;
1739 
1740   t2618 = t7 * t299 * t1117;
1741   t2621 = t7 * t305 * t482;
1742   t2624 = t7 * t122 * t1121;
1743   t2628 = t7 * t314 * t216 / 0.9e1;
1744   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1745     v4rho3sigma[2] = -0.2e1 / 0.3e1 * t2618 - 0.2e1 / 0.3e1 * t2621 - 0.2e1 * t2624 + t2628 + t179 * t746 * t216 * t599 / 0.9e1 - t1132 * t1150 * t608 / 0.2e1 - t179 * t481 * t216 * t713 / 0.2e1;
1746 
1747   t2644 = 0.4e1 / 0.3e1 * t1128 * t1138;
1748   t2657 = t83 * t1 * t3;
1749   t2658 = t2657 * t487;
1750   t2660 = t1127 * t1131;
1751   t2662 = 0.4e1 / 0.9e1 * t2660 * t1134;
1752   t2664 = t335 * t1 * t3;
1753   t2665 = t2664 * t1144;
1754   t2676 = t577 * t196;
1755   t2680 = t1141 * t1131;
1756   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1757     v4rho3sigma[3] = -t1142 * t442 * t196 / 0.2e1 - t2644 - t179 * t442 * t1107 * t135 / 0.2e1 - 0.4e1 / 0.3e1 * t2590 - 0.2e1 / 0.3e1 * t2593 + t2597 - t2599 / 0.2e1 - t2602 / 0.3e1 - t1142 * t442 * t475 * t8 - 0.2e1 / 0.9e1 * t2658 - t2662 + t2665 / 0.6e1 - t1132 * t224 * t475 * t323 / 0.3e1 - t1132 * t1133 * t760 / 0.6e1 - 0.2e1 / 0.9e1 * t2612 - 0.4e1 / 0.9e1 * t2615 + t1132 * t2676 * t733 / 0.9e1 - t2680 * t1133 * t742 / 0.3e1;
1758 
1759   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1760     v4rho3sigma[4] = 0.0e0;
1761 
1762   t2687 = t2657 * t492;
1763   t2690 = 0.4e1 / 0.9e1 * t2660 * t1151;
1764   t2691 = t2664 * t1155;
1765   t2693 = t597 * t216;
1766   t2706 = t2657 * t524;
1767   t2709 = 0.4e1 / 0.3e1 * t1128 * t1161;
1768   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1769     v4rho3sigma[5] = -0.2e1 / 0.9e1 * t2618 - 0.4e1 / 0.9e1 * t2621 - 0.2e1 / 0.3e1 * t2624 + t2628 - 0.2e1 / 0.9e1 * t2687 - t2690 - t2691 / 0.6e1 + t1132 * t2693 * t748 / 0.9e1 + t2680 * t1150 * t769 / 0.3e1 - t1132 * t1150 * t764 / 0.6e1 + t1142 * t481 * t216 / 0.2e1 - t2706 / 0.6e1 - t2709 - t179 * t343 * t523 * t300 / 0.6e1 - t179 * t481 * t523 * t309 / 0.2e1;
1770 
1771   t2722 = t1128 * t1166;
1772   t2734 = t1128 * t1170;
1773   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1774     v4rho3sigma[6] = -0.4e1 / 0.9e1 * t2658 - t2662 - t2644 - 0.4e1 / 0.3e1 * t2665 - 0.2e1 / 0.9e1 * t2615 + t2597 - t2602 / 0.6e1 - 0.2e1 / 0.9e1 * t2722 + t1132 * t2676 * t837 / 0.9e1 - t179 * t320 * t475 * t369 / 0.6e1 - t2680 * t1133 * t841 / 0.3e1 - 0.2e1 / 0.3e1 * t2734 - t1132 * t1133 * t845 / 0.6e1 - t179 * t442 * t475 * t375 / 0.2e1 - t179 * t442 * t196 * t880 / 0.2e1;
1775 
1776   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1777     v4rho3sigma[7] = 0.0e0;
1778 
1779   t2751 = t1128 * t1175;
1780   t2759 = t1128 * t1179;
1781   t2768 = t1128 * t1182;
1782   t2777 = t1128 * t1235;
1783   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1784     v4rho3sigma[8] = -0.4e1 / 0.9e1 * t2687 - t2690 + 0.4e1 / 0.3e1 * t2691 - 0.2e1 / 0.9e1 * t2621 + t2628 - t2706 / 0.3e1 - t2709 - 0.2e1 / 0.9e1 * t2751 + t1132 * t2693 * t849 / 0.9e1 + t2680 * t1150 * t853 / 0.3e1 - 0.4e1 / 0.3e1 * t2759 - t1132 * t298 * t523 * t346 / 0.3e1 + t1142 * t481 * t523 * t8 - 0.2e1 / 0.3e1 * t2768 - t1132 * t1150 * t861 / 0.6e1 - t179 * t481 * t216 * t888 / 0.2e1 - t2777 / 0.2e1 - t179 * t481 * t1234 * t123 / 0.2e1;
1785 
1786   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1787     v4rho3sigma[9] = -0.2e1 / 0.3e1 * t2658 - 0.2e1 / 0.3e1 * t2722 - 0.2e1 * t2734 + t2597 + t179 * t731 * t196 * t1018 / 0.9e1 - t1132 * t1133 * t1023 / 0.2e1 - t179 * t442 * t196 * t935 / 0.2e1;
1788 
1789   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1790     v4rho3sigma[10] = 0.0e0;
1791 
1792   t2821 = t53 * t975;
1793   t2880 = -0.12542328128682275075e-193 * t1302 * t406 * t168 * t430 + 0.66892416686305467067e-193 * t401 * t968 * t1209 + 0.33446208343152733533e-193 * t401 * t968 * t1212 + 0.19076726240168596163e-194 * t401 * t2300 * t213 - 0.12263609725822668962e-192 * t401 * t2304 * t505 - 0.11025068610270821606e-191 * t407 * t520 * t430 - 0.36750228700902738685e-192 * t407 * t213 * t1008 - 0.10914817924168113389e-189 * t976 * t520 * t408 - 0.11025068610270821606e-191 * t407 * t1231 * t168 - 0.36750228700902738685e-194 * t153 * (-0.1232e4 / 0.27e2 * t23 * t959 * t58 + 0.51167577030977359778e1 * t103 * t987 * t164 * sigma[2] - 0.96840172523411230156e-2 * t993 * t427 * t157 + 0.11478887496814161219e-3 * t1000 * t1003 * t421 * t18 * t22 - 0.46952273766562626373e-6 * t997 / t50 / t2355 * t2363) + 0.33111746259721206198e-191 * t401 * t2311 * t1206;
1794   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1795     v4rho3sigma[11] = -0.3e1 / 0.8e1 * t179 * t200 * (0.33446208343152733533e-195 * t401 * t402 * t1231 - 0.35655071885615837072e-188 * t2322 * t213 * t977 - 0.10914817924168113389e-189 * t976 * t505 * t430 - 0.41734945439179794323e0 * t23 * t959 * t63 - 0.12263609725822668962e-194 * t23 * t965 + 0.33446208343152733533e-195 * t23 * t972 - 0.41807760428940916917e-196 * t23 * t498 * t1008 - 0.12263609725822668962e-194 * t401 * t964 * t520 + 0.33446208343152733533e-193 * t23 * t969 - 0.41389682824651507748e-192 * t23 * t2821 * t977 + t2880) - t179 * t481 * t216 * t953 / 0.2e1 - 0.3e1 / 0.2e1 * t179 * t481 * t523 * t389 - t179 * t343 * t523 * t380 / 0.2e1 - 0.3e1 / 0.2e1 * t179 * t481 * t1234 * t140 - 0.3e1 / 0.2e1 * t2777 - t2706 / 0.2e1 + t2628 - t1132 * t1150 * t1032 / 0.2e1 + t179 * t746 * t216 * t1027 / 0.9e1 - 0.2e1 / 0.3e1 * t2751 - 0.4e1 * t2759 - 0.2e1 * t2768 - 0.2e1 / 0.3e1 * t2687;
1796 
1797   t2915 = t2657 * t544 / 0.6e1;
1798   t2916 = t1128 * t1241;
1799   t2918 = t1128 * t1283;
1800   t2964 = -0.40878699086075563207e-195 * t261 * t659 * t540 + 0.22297472228768489022e-193 * t261 * t663 * t1266 - 0.82779365649303015496e-192 * t1246 * t670 * t193 * t268 - 0.40878699086075563207e-193 * t261 * t1421 * t530 + 0.22074497506480804132e-191 * t261 * t1428 * t1257 + 0.44594944457536978045e-193 * t261 * t663 * t1260 + 0.22297472228768489022e-195 * t261 * t262 * t1279 - 0.81757398172151126414e-195 * t23 * t1070 + 0.44594944457536978044e-193 * t23 * t90 * t1248 - 0.16723104171576366767e-193 * t1246 * t266 * t472 * t115 - 0.83615520857881833834e-194 * t1246 * t1247 * t290;
1801   t2982 = t472 * t472;
1802   t3012 = -0.36382726413893711298e-190 * t671 * t540 * t268 - 0.35655071885615837072e-188 * t1439 * t530 * t268 - 0.14553090565557484519e-189 * t671 * t457 * t472 + 0.44594944457536978044e-195 * t23 * t1076 - 0.83615520857881833834e-196 * t23 * t450 * t1104 - 0.36382726413893711298e-190 * t671 * t530 * t290 - 0.73500457401805477370e-192 * t267 * t2982 - 0.73500457401805477370e-192 * t267 * t193 * t1104 - 0.73500457401805477370e-192 * t267 * t1279 * t115 - 0.36750228700902738685e-192 * t267 * t540 * t290 - 0.36750228700902738685e-194 * t96 * (-0.76852487240203149073e0 * t103 * t1091 + 0.23666793382550744540e-2 * t1095 * sigma[0] - 0.36274574252572840931e-4 * t1100 * t698 * t104 * t18 * t22 + 0.17607102662460984890e-6 * t281 / t25 / t282 / t274 * t1480);
1803   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1804     v4rho2sigma2[0] = -t2915 - 0.4e1 / 0.3e1 * t2916 - t2918 - t179 * t320 * t543 * t226 / 0.6e1 - t179 * t442 * t1282 * t77 - t179 * t442 * t543 * t241 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t180 * (t2964 + t3012);
1805 
1806   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1807     v4rho2sigma2[1] = 0.0e0;
1808 
1809   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1810     v4rho2sigma2[2] = 0.0e0;
1811 
1812   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1813     v4rho2sigma2[3] = 0.0e0;
1814 
1815   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1816     v4rho2sigma2[4] = 0.0e0;
1817 
1818   t3018 = t2657 * t564 / 0.6e1;
1819   t3019 = t1128 * t1288;
1820   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1821     v4rho2sigma2[5] = -t3018 - 0.4e1 / 0.3e1 * t3019 - t179 * t343 * t563 * t300 / 0.6e1 - t179 * t481 * t563 * t309 / 0.2e1;
1822 
1823   t3031 = t1128 * t1292;
1824   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1825     v4rho2sigma2[6] = -t2915 - 0.2e1 / 0.3e1 * t2916 - t2918 / 0.2e1 - 0.2e1 / 0.3e1 * t3031 - t1132 * t224 * t543 * t323 / 0.6e1 - t179 * t442 * t1282 * t135 / 0.2e1 - t1142 * t442 * t543 * t8 / 0.2e1;
1826 
1827   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1828     v4rho2sigma2[7] = 0.0e0;
1829 
1830   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1831     v4rho2sigma2[8] = 0.0e0;
1832 
1833   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1834     v4rho2sigma2[9] = 0.0e0;
1835 
1836   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1837     v4rho2sigma2[10] = 0.0e0;
1838 
1839   t3046 = t1128 * t1297;
1840   t3056 = t1128 * t1339;
1841   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1842     v4rho2sigma2[11] = -t3018 - 0.2e1 / 0.3e1 * t3019 - 0.2e1 / 0.3e1 * t3046 - t1132 * t298 * t563 * t346 / 0.6e1 + t1142 * t481 * t563 * t8 / 0.2e1 - t3056 / 0.2e1 - t179 * t481 * t1338 * t123 / 0.2e1;
1843 
1844   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1845     v4rho2sigma2[12] = -t2915 - 0.4e1 / 0.3e1 * t3031 - t179 * t320 * t543 * t369 / 0.6e1 - t179 * t442 * t543 * t375 / 0.2e1;
1846 
1847   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1848     v4rho2sigma2[13] = 0.0e0;
1849 
1850   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1851     v4rho2sigma2[14] = 0.0e0;
1852 
1853   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1854     v4rho2sigma2[15] = 0.0e0;
1855 
1856   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1857     v4rho2sigma2[16] = 0.0e0;
1858 
1859   t3104 = t520 * t520;
1860   t3116 = 0.22297472228768489022e-195 * t401 * t402 * t1335 - 0.81757398172151126414e-195 * t23 * t1197 + 0.44594944457536978044e-193 * t23 * t147 * t1304 - 0.16723104171576366767e-193 * t1302 * t406 * t520 * t168 - 0.83615520857881833834e-194 * t1302 * t1303 * t430 - 0.36750228700902738685e-192 * t407 * t560 * t430 - 0.36382726413893711298e-190 * t976 * t550 * t430 - 0.73500457401805477370e-192 * t407 * t3104 - 0.73500457401805477370e-192 * t407 * t213 * t1231 - 0.73500457401805477370e-192 * t407 * t1335 * t168 - 0.36382726413893711298e-190 * t976 * t560 * t408;
1861   t3165 = -0.83615520857881833834e-196 * t23 * t498 * t1231 - 0.35655071885615837072e-188 * t2322 * t550 * t408 - 0.14553090565557484519e-189 * t976 * t505 * t520 + 0.44594944457536978044e-195 * t23 * t1203 - 0.36750228700902738685e-194 * t153 * (-0.76852487240203149073e0 * t103 * t1218 + 0.23666793382550744540e-2 * t1222 * sigma[2] - 0.36274574252572840931e-4 * t1227 * t1003 * t157 * t18 * t22 + 0.17607102662460984890e-6 * t421 / t50 / t422 / t414 * t2363) - 0.40878699086075563207e-193 * t401 * t2304 * t550 + 0.22074497506480804132e-191 * t401 * t2311 * t1313 + 0.44594944457536978045e-193 * t401 * t968 * t1316 - 0.40878699086075563207e-195 * t401 * t964 * t560 + 0.22297472228768489022e-193 * t401 * t968 * t1322 - 0.82779365649303015496e-192 * t1302 * t975 * t213 * t408;
1862   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1863     v4rho2sigma2[17] = -t3018 - 0.4e1 / 0.3e1 * t3046 - t3056 - t179 * t343 * t563 * t380 / 0.6e1 - t179 * t481 * t1338 * t140 - t179 * t481 * t563 * t389 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t200 * (t3116 + t3165);
1864 
1865   t3170 = t1128 * t1365;
1866   t3233 = 0.33446208343152733533e-193 * t23 * t1254 - 0.12416904847395452324e-191 * t1246 * t670 * t530 * t115 - 0.25084656257364550150e-193 * t1246 * t1247 * t472 + 0.33446208343152733533e-195 * t23 * t1263 - 0.12542328128682275075e-193 * t1246 * t266 * t540 * t115 - 0.12542328128682275075e-195 * t23 * t450 * t1279 + 0.11037248753240402066e-191 * t261 * t1428 * t1348 - 0.35655071885615837072e-188 * t1439 * t1348 * t115 - 0.10914817924168113389e-189 * t671 * t530 * t472 + 0.33446208343152733536e-193 * t261 * t663 * t1351 - 0.10914817924168113390e-189 * t671 * t1351 * t115 - 0.11025068610270821606e-191 * t267 * t472 * t540 - 0.11025068610270821606e-191 * t267 * t193 * t1279 + 0.11148736114384244511e-195 * t261 * t262 * t1361 - 0.36750228700902738685e-192 * t267 * t1361 * t115 - 0.36750228700902738685e-194 * t96 * (-0.38263677923884437281e-3 * t1271 + 0.10156880790720395460e-4 * t1275 * t698 * t261 - 0.66026634984228693337e-7 * t104 / t25 / t1465 * t1480);
1867   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1868     v4rhosigma3[0] = -t3170 / 0.2e1 - t179 * t442 * t1364 * t77 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t180 * t3233;
1869 
1870   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1871     v4rhosigma3[1] = 0.0e0;
1872 
1873   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1874     v4rhosigma3[2] = 0.0e0;
1875 
1876   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1877     v4rhosigma3[3] = 0.0e0;
1878 
1879   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1880     v4rhosigma3[4] = 0.0e0;
1881 
1882   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1883     v4rhosigma3[5] = 0.0e0;
1884 
1885   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1886     v4rhosigma3[6] = 0.0e0;
1887 
1888   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1889     v4rhosigma3[7] = 0.0e0;
1890 
1891   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1892     v4rhosigma3[8] = 0.0e0;
1893 
1894   t3237 = t1128 * t1391;
1895   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1896     v4rhosigma3[9] = -t179 * t481 * t1390 * t123 / 0.2e1 - t3237 / 0.2e1;
1897 
1898   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1899     v4rhosigma3[10] = -t179 * t442 * t1364 * t135 / 0.2e1 - t3170 / 0.2e1;
1900 
1901   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1902     v4rhosigma3[11] = 0.0e0;
1903 
1904   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1905     v4rhosigma3[12] = 0.0e0;
1906 
1907   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1908     v4rhosigma3[13] = 0.0e0;
1909 
1910   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1911     v4rhosigma3[14] = 0.0e0;
1912 
1913   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1914     v4rhosigma3[15] = 0.0e0;
1915 
1916   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1917     v4rhosigma3[16] = 0.0e0;
1918 
1919   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1920     v4rhosigma3[17] = 0.0e0;
1921 
1922   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1923     v4rhosigma3[18] = 0.0e0;
1924 
1925   t3308 = 0.33446208343152733533e-193 * t23 * t1310 - 0.12416904847395452324e-191 * t1302 * t975 * t550 * t168 - 0.25084656257364550150e-193 * t1302 * t1303 * t520 + 0.33446208343152733533e-195 * t23 * t1319 - 0.12542328128682275075e-193 * t1302 * t406 * t560 * t168 - 0.12542328128682275075e-195 * t23 * t498 * t1335 + 0.11037248753240402066e-191 * t401 * t2311 * t1374 - 0.35655071885615837072e-188 * t2322 * t1374 * t168 - 0.10914817924168113389e-189 * t976 * t550 * t520 + 0.33446208343152733536e-193 * t401 * t968 * t1377 - 0.10914817924168113390e-189 * t976 * t1377 * t168 - 0.11025068610270821606e-191 * t407 * t520 * t560 - 0.11025068610270821606e-191 * t407 * t213 * t1335 + 0.11148736114384244511e-195 * t401 * t402 * t1387 - 0.36750228700902738685e-192 * t407 * t1387 * t168 - 0.36750228700902738685e-194 * t153 * (-0.38263677923884437281e-3 * t1327 + 0.10156880790720395460e-4 * t1331 * t1003 * t401 - 0.66026634984228693337e-7 * t157 / t50 / t2348 * t2363);
1926   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1927     v4rhosigma3[19] = -t3237 / 0.2e1 - t179 * t481 * t1390 * t140 / 0.2e1 - 0.3e1 / 0.8e1 * t179 * t200 * t3308;
1928 
1929   t3321 = t530 * t530;
1930   t3327 = t540 * t540;
1931   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1932     v4sigma4[0] = -0.3e1 / 0.8e1 * t179 * t180 * (-0.16555873129860603099e-191 * t23 * t2514 * t1348 - 0.50169312514729100301e-193 * t1246 * t1247 * t540 - 0.16723104171576366767e-195 * t23 * t450 * t1361 - 0.35655071885615837072e-188 * t1439 * t3321 - 0.21829635848336226779e-189 * t671 * t530 * t540 - 0.11025068610270821606e-191 * t267 * t3327 - 0.14700091480361095474e-191 * t267 * t193 * t1361 - 0.36750228700902738685e-194 * t96 * (-0.21764744551543704558e-5 * t1357 * t698 * t23 + 0.24759988119085760001e-7 * sigma[0] / t25 / t693 * t1480));
1933 
1934   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1935     v4sigma4[1] = 0.0e0;
1936 
1937   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1938     v4sigma4[2] = 0.0e0;
1939 
1940   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1941     v4sigma4[3] = 0.0e0;
1942 
1943   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1944     v4sigma4[4] = 0.0e0;
1945 
1946   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1947     v4sigma4[5] = 0.0e0;
1948 
1949   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1950     v4sigma4[6] = 0.0e0;
1951 
1952   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1953     v4sigma4[7] = 0.0e0;
1954 
1955   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1956     v4sigma4[8] = 0.0e0;
1957 
1958   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1959     v4sigma4[9] = 0.0e0;
1960 
1961   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1962     v4sigma4[10] = 0.0e0;
1963 
1964   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1965     v4sigma4[11] = 0.0e0;
1966 
1967   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1968     v4sigma4[12] = 0.0e0;
1969 
1970   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1971     v4sigma4[13] = 0.0e0;
1972 
1973   t3357 = t550 * t550;
1974   t3363 = t560 * t560;
1975   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1976     v4sigma4[14] = -0.3e1 / 0.8e1 * t179 * t200 * (-0.16555873129860603099e-191 * t23 * t2821 * t1374 - 0.50169312514729100301e-193 * t1302 * t1303 * t560 - 0.16723104171576366767e-195 * t23 * t498 * t1387 - 0.35655071885615837072e-188 * t2322 * t3357 - 0.21829635848336226779e-189 * t976 * t550 * t560 - 0.11025068610270821606e-191 * t407 * t3363 - 0.14700091480361095474e-191 * t407 * t213 * t1387 - 0.36750228700902738685e-194 * t153 * (-0.21764744551543704558e-5 * t1383 * t1003 * t23 + 0.24759988119085760001e-7 * sigma[2] / t50 / t998 * t2363));
1977 
1978 #ifndef XC_DONT_COMPILE_MXC
1979 
1980   if(order < 5) return;
1981 
1982 
1983 #endif
1984 
1985 #endif
1986 
1987 #endif
1988 
1989 #endif
1990 
1991 #endif
1992 
1993 
1994 }
1995 
1996