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