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/lda_exc/lda_c_vwn_1.mpl
11   Type of functional: lda_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,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
20 {
21 
22 #ifndef XC_DONT_COMPILE_EXC
23   double t1, t2, t3, t4, t5, t6, t7, t8;
24   double t9, t10, t12, t14, t15, t19, t20, t21;
25   double t24, t25, t27, t28, t30, t31;
26 
27 #ifndef XC_DONT_COMPILE_VXC
28   double t33, t34, t38, t39, t40, t41, t42, t44;
29   double t45, t46, t50, t55, t57, t58, t59, t60;
30   double t61, t63, t64, t66, t68, t69, t72, t74;
31   double t75, t78, t80, t81, t82, t83;
32 
33 #ifndef XC_DONT_COMPILE_FXC
34   double t90, t92, t93, t95, t97, t102, t103, t104;
35   double t108, t111, t112, t113, t114, t115, t117, t124;
36   double t130, t131, t134, t135, t136, t138, t139, t141;
37   double t143, t144, t146, t149, t152, t156, t158, t160;
38   double t161, t162, t163, t164, t166, t169, t170, t171;
39   double t175, t176, t177, t178, t183, t187, t188, t189;
40   double t192, t193, t194, t195, t196, t198;
41 
42 #ifndef XC_DONT_COMPILE_KXC
43   double t212, t214, t215, t217, t219, t221, t223, t230;
44   double t231, t232, t233, t237, t241, t246, t247, t248;
45   double t249, t250, t253, t260, t266, t267, t269, t271;
46   double t272, t274, t275, t277, t278, t279, t280, t282;
47   double t284, t288, t289, t290, t293, t295, t296, t297;
48   double t298, t299, t301, t304, t305, t306, t307, t309;
49   double t311, t312, t314, t316, t317, t318, t319, t320;
50   double t322, t323, t325, t327, t328, t329, t330, t331;
51   double t333, t334, t336, t340, t341, t345, t359, t360;
52   double t366, t371, t377, t378, t379, t381, t383, t385;
53   double t389, t393, t395, t396, t397, t399, t400, t401;
54   double t403, t404, t405, t407;
55 
56 #ifndef XC_DONT_COMPILE_LXC
57   double t426, t427, t428, t453, t455, t459, t462, t464;
58   double t472, t478, t482, t486, t489, t502, t504, t508;
59   double t512, t516, t520, t526, t527, t528, t530, t534;
60   double t535, t543, t547, t603, t644, t658, t665, t676;
61   double t722, t734, t738, t773, t780, t794;
62 #endif
63 
64 #endif
65 
66 #endif
67 
68 #endif
69 
70 #endif
71 
72 
73 
74   t1 = M_CBRT3;
75   t2 = 0.1e1 / M_PI;
76   t3 = POW_1_3(t2);
77   t4 = t1 * t3;
78   t5 = M_CBRT4;
79   t6 = t5 * t5;
80   t7 = POW_1_3(rho[0]);
81   t8 = 0.1e1 / t7;
82   t9 = t6 * t8;
83   t10 = t4 * t9;
84   t12 = sqrt(t10);
85   t14 = t10 / 0.4e1 + 0.18637200000000000000e1 * t12 + 0.129352e2;
86   t15 = 0.1e1 / t14;
87   t19 = log(t4 * t9 * t15 / 0.4e1);
88   t20 = 0.310907e-1 * t19;
89   t21 = t12 + 0.372744e1;
90   t24 = atan(0.61519908197590802322e1 / t21);
91   t25 = 0.38783294878113014393e-1 * t24;
92   t27 = t12 / 0.2e1 + 0.10498e0;
93   t28 = t27 * t27;
94   t30 = log(t28 * t15);
95   t31 = 0.96902277115443742139e-3 * t30;
96   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
97     zk[0] = t20 + t25 + t31;
98 
99 #ifndef XC_DONT_COMPILE_VXC
100 
101   if(order < 1) return;
102 
103 
104   t33 = 0.1e1 / t7 / rho[0];
105   t34 = t6 * t33;
106   t38 = t4 * t6;
107   t39 = t14 * t14;
108   t40 = 0.1e1 / t39;
109   t41 = t8 * t40;
110   t42 = t4 * t34;
111   t44 = 0.1e1 / t12;
112   t45 = t44 * t1;
113   t46 = t3 * t6;
114   t50 = -t42 / 0.12e2 - 0.31062000000000000000e0 * t45 * t46 * t33;
115   t55 = t1 * t1;
116   t57 = 0.1e1 / t3;
117   t58 = (-t4 * t34 * t15 / 0.12e2 - t38 * t41 * t50 / 0.4e1) * t55 * t57;
118   t59 = t5 * t7;
119   t60 = t59 * t14;
120   t61 = t58 * t60;
121   t63 = t21 * t21;
122   t64 = 0.1e1 / t63;
123   t66 = t64 * t44 * t1;
124   t68 = 0.37846991046400000000e2 * t64 + 0.1e1;
125   t69 = 0.1e1 / t68;
126   t72 = t66 * t46 * t33 * t69;
127   t74 = t27 * t15;
128   t75 = t74 * t44;
129   t78 = t28 * t40;
130   t80 = -t75 * t42 / 0.6e1 - t78 * t50;
131   t81 = 0.1e1 / t28;
132   t82 = t80 * t81;
133   t83 = t82 * t14;
134   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
135     vrho[0] = t20 + t25 + t31 + rho[0] * (0.10363566666666666667e-1 * t61 + 0.39765745675026770179e-1 * t72 + 0.96902277115443742139e-3 * t83);
136 
137 #ifndef XC_DONT_COMPILE_FXC
138 
139   if(order < 2) return;
140 
141 
142   t90 = rho[0] * rho[0];
143   t92 = 0.1e1 / t7 / t90;
144   t93 = t6 * t92;
145   t95 = t4 * t93 * t15;
146   t97 = t33 * t40;
147   t102 = 0.1e1 / t39 / t14;
148   t103 = t8 * t102;
149   t104 = t50 * t50;
150   t108 = t4 * t93;
151   t111 = 0.1e1 / t12 / t10;
152   t112 = t111 * t55;
153   t113 = t3 * t3;
154   t114 = t113 * t5;
155   t115 = t7 * t7;
156   t117 = 0.1e1 / t115 / t90;
157   t124 = t108 / 0.9e1 - 0.20708000000000000000e0 * t112 * t114 * t117 + 0.41416000000000000000e0 * t45 * t46 * t92;
158   t130 = (t95 / 0.9e1 + t38 * t97 * t50 / 0.6e1 + t38 * t103 * t104 / 0.2e1 - t38 * t41 * t124 / 0.4e1) * t55 * t57;
159   t131 = t130 * t60;
160   t134 = t5 / t115;
161   t135 = t134 * t14;
162   t136 = t58 * t135;
163   t138 = t59 * t50;
164   t139 = t58 * t138;
165   t141 = t63 * t21;
166   t143 = 0.1e1 / t141 * t1;
167   t144 = t143 * t3;
168   t146 = t144 * t93 * t69;
169   t149 = t64 * t111 * t55;
170   t152 = t149 * t114 * t117 * t69;
171   t156 = t66 * t46 * t92 * t69;
172   t158 = t63 * t63;
173   t160 = 0.1e1 / t158 / t21;
174   t161 = t160 * t1;
175   t162 = t161 * t3;
176   t163 = t68 * t68;
177   t164 = 0.1e1 / t163;
178   t166 = t162 * t93 * t164;
179   t169 = t27 * t40;
180   t170 = t169 * t45;
181   t171 = t33 * t50;
182   t175 = t74 * t111;
183   t176 = t55 * t113;
184   t177 = t5 * t117;
185   t178 = t176 * t177;
186   t183 = t28 * t102;
187   t187 = t95 / 0.72e2 + t170 * t46 * t171 / 0.3e1 - t175 * t178 / 0.9e1 + 0.2e1 / 0.9e1 * t75 * t108 + 0.2e1 * t183 * t104 - t78 * t124;
188   t188 = t187 * t81;
189   t189 = t188 * t14;
190   t192 = 0.1e1 / t28 / t27;
191   t193 = t80 * t192;
192   t194 = t14 * t44;
193   t195 = t193 * t194;
194   t196 = t195 * t42;
195   t198 = t82 * t50;
196   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
197     v2rho2[0] = 0.20727133333333333334e-1 * t61 + 0.79531491350053540358e-1 * t72 + 0.19380455423088748428e-2 * t83 + rho[0] * (0.10363566666666666667e-1 * t131 + 0.34545222222222222223e-2 * t136 + 0.10363566666666666667e-1 * t139 + 0.13255248558342256726e-1 * t146 + 0.26510497116684513453e-1 * t152 - 0.53020994233369026905e-1 * t156 - 0.50167127350538589834e0 * t166 + 0.96902277115443742139e-3 * t189 + 0.16150379519240623690e-3 * t196 + 0.96902277115443742139e-3 * t198);
198 
199 #ifndef XC_DONT_COMPILE_KXC
200 
201   if(order < 3) return;
202 
203 
204   t212 = t90 * rho[0];
205   t214 = 0.1e1 / t7 / t212;
206   t215 = t6 * t214;
207   t217 = t4 * t215 * t15;
208   t219 = t92 * t40;
209   t221 = t38 * t219 * t50;
210   t223 = t33 * t102;
211   t230 = t39 * t39;
212   t231 = 0.1e1 / t230;
213   t232 = t8 * t231;
214   t233 = t104 * t50;
215   t237 = t50 * t124;
216   t241 = t4 * t215;
217   t246 = 0.1e1 / t12 / t176 / t134 / 0.4e1;
218   t247 = t246 * t2;
219   t248 = t90 * t90;
220   t249 = 0.1e1 / t248;
221   t250 = t247 * t249;
222   t253 = 0.1e1 / t115 / t212;
223   t260 = -0.7e1 / 0.27e2 * t241 - 0.12424800000000000000e1 * t250 + 0.82832000000000000000e0 * t112 * t114 * t253 - 0.96637333333333333333e0 * t45 * t46 * t214;
224   t266 = (-0.7e1 / 0.27e2 * t217 - t221 / 0.3e1 - t38 * t223 * t104 / 0.2e1 + t38 * t97 * t124 / 0.4e1 - 0.3e1 / 0.2e1 * t38 * t232 * t233 + 0.3e1 / 0.2e1 * t38 * t103 * t237 - t38 * t41 * t260 / 0.4e1) * t55 * t57;
225   t267 = t266 * t60;
226   t269 = t130 * t138;
227   t271 = t134 * t50;
228   t272 = t58 * t271;
229   t274 = t59 * t124;
230   t275 = t58 * t274;
231   t277 = 0.1e1 / t113;
232   t278 = t161 * t277;
233   t279 = t215 * t164;
234   t280 = t278 * t279;
235   t282 = t162 * t279;
236   t284 = t130 * t135;
237   t288 = t5 / t115 / rho[0];
238   t289 = t288 * t14;
239   t290 = t58 * t289;
240   t293 = t144 * t215 * t69;
241   t295 = 0.1e1 / t158;
242   t296 = t295 * t55;
243   t297 = t296 * t113;
244   t298 = t5 * t253;
245   t299 = t69 * t44;
246   t301 = t297 * t298 * t299;
247   t304 = 0.1e1 / t158 / t63;
248   t305 = t304 * t55;
249   t306 = t305 * t113;
250   t307 = t164 * t44;
251   t309 = t306 * t298 * t307;
252   t311 = t143 * t277;
253   t312 = t2 * t69;
254   t314 = t311 * t215 * t312;
255   t316 = 0.10363566666666666667e-1 * t267 + 0.20727133333333333334e-1 * t269 + 0.69090444444444444446e-2 * t272 + 0.10363566666666666667e-1 * t275 - 0.26614487661862784322e-1 * t280 + 0.18394613361864149606e1 * t282 + 0.69090444444444444446e-2 * t284 - 0.23030148148148148149e-2 * t290 - 0.48602578047254941329e-1 * t293 + 0.26510497116684513452e-1 * t301 - 0.23411326096918008589e1 * t309 + 0.22092080930570427878e-2 * t314;
256   t317 = t158 * t158;
257   t318 = 0.1e1 / t317;
258   t319 = t318 * t55;
259   t320 = t319 * t113;
260   t322 = 0.1e1 / t163 / t68;
261   t323 = t322 * t44;
262   t325 = t320 * t298 * t323;
263   t327 = t28 * t28;
264   t328 = 0.1e1 / t327;
265   t329 = t80 * t328;
266   t330 = t329 * t14;
267   t331 = t330 * t108;
268   t333 = t64 * t246;
269   t334 = t2 * t249;
270   t336 = t333 * t334 * t69;
271   t340 = t27 * t102;
272   t341 = t340 * t45;
273   t345 = t169 * t112;
274   t359 = t1 * t277 * t6;
275   t360 = t214 * t2;
276   t366 = t176 * t298;
277   t371 = t28 * t231;
278   t377 = -0.11e2 / 0.216e3 * t217 - t221 / 0.24e2 - t341 * t46 * t33 * t104 + t345 * t114 * t117 * t50 / 0.3e1 - 0.2e1 / 0.3e1 * t170 * t46 * t92 * t50 + t170 * t46 * t33 * t124 / 0.2e1 + t359 * t360 * t15 / 0.432e3 - 0.2e1 / 0.3e1 * t74 * t250 + 0.4e1 / 0.9e1 * t175 * t366 - 0.14e2 / 0.27e2 * t75 * t241 - 0.6e1 * t371 * t233 + 0.6e1 * t183 * t237 - t78 * t260;
279   t378 = t377 * t81;
280   t379 = t378 * t14;
281   t381 = t188 * t50;
282   t383 = t82 * t124;
283   t385 = t195 * t108;
284   t389 = t149 * t114 * t253 * t69;
285   t393 = t66 * t46 * t214 * t69;
286   t395 = t187 * t192;
287   t396 = t395 * t194;
288   t397 = t396 * t42;
289   t399 = t50 * t44;
290   t400 = t193 * t399;
291   t401 = t400 * t42;
292   t403 = t14 * t111;
293   t404 = t193 * t403;
294   t405 = t404 * t178;
295   t407 = 0.50631328524251801694e2 * t325 + 0.40375948798101559225e-4 * t331 + 0.15906298270010708072e0 * t336 + 0.96902277115443742139e-3 * t379 + 0.19380455423088748428e-2 * t381 + 0.96902277115443742139e-3 * t383 - 0.21533839358987498253e-3 * t385 - 0.10604198846673805381e0 * t389 + 0.12371565321119439611e0 * t393 + 0.32300759038481247380e-3 * t397 + 0.32300759038481247380e-3 * t401 + 0.10766919679493749127e-3 * t405;
296   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
297     v3rho3[0] = 0.31090700000000000001e-1 * t131 + 0.10363566666666666667e-1 * t136 + 0.31090700000000000001e-1 * t139 + 0.39765745675026770179e-1 * t146 + 0.79531491350053540358e-1 * t152 - 0.15906298270010708072e0 * t156 - 0.15050138205161576950e1 * t166 + 0.29070683134633122642e-2 * t189 + 0.48451138557721871070e-3 * t196 + 0.29070683134633122642e-2 * t198 + rho[0] * (t316 + t407);
298 
299 #ifndef XC_DONT_COMPILE_LXC
300 
301   if(order < 4) return;
302 
303 
304   t426 = 0.1e1 / t7 / t248;
305   t427 = t6 * t426;
306   t428 = t427 * t164;
307   t453 = 0.1e1 / t158 / t141;
308   t455 = 0.1e1 / t115 / t248;
309   t459 = t55 * t57 * t5;
310   t462 = 0.10363566666666666667e-1 * t266 * t135 + 0.20727133333333333334e-1 * t130 * t271 - 0.69090444444444444446e-2 * t58 * t288 * t50 + 0.19517290952032708503e0 * t278 * t428 - 0.77340987998746992661e1 * t162 * t428 - 0.69090444444444444446e-2 * t130 * t289 + 0.38383580246913580248e-2 * t58 * t177 * t14 + 0.20435174860777645786e0 * t144 * t427 * t69 + 0.31090700000000000001e-1 * t130 * t274 + 0.64601518076962494762e-3 * t193 * t14 * t250 + 0.10363566666666666667e-1 * t58 * t134 * t124 + 0.10363566666666666667e-1 * t58 * t59 * t260 - 0.10645795064745113728e0 * t453 * t455 * t164 * t459;
311   t464 = 0.1e1 / t317 / t21;
312   t472 = t163 * t163;
313   t478 = t4 * t427 * t15;
314   t482 = t38 * t214 * t40 * t50;
315   t486 = t38 * t92 * t102 * t104;
316   t489 = t38 * t219 * t124;
317   t502 = 0.1e1 / t230 / t14;
318   t504 = t104 * t104;
319   t508 = t104 * t124;
320   t512 = t124 * t124;
321   t516 = t50 * t260;
322   t520 = t4 * t427;
323   t526 = 0.1e1 / t12 / t2 * rho[0] / 0.48e2;
324   t527 = t526 * t2;
325   t528 = t248 * rho[0];
326   t530 = 0.1e1 / t7 / t528;
327   t534 = 0.1e1 / t528;
328   t535 = t247 * t534;
329   t543 = 0.70e2 / 0.81e2 * t520 - 0.10354000000000000000e1 * t527 * t530 * t38 + 0.99398400000000000000e1 * t535 - 0.36814222222222222222e1 * t112 * t114 * t455 + 0.32212444444444444444e1 * t45 * t46 * t426;
330   t547 = 0.70e2 / 0.81e2 * t478 + 0.28e2 / 0.27e2 * t482 + 0.4e1 / 0.3e1 * t486 - 0.2e1 / 0.3e1 * t489 + 0.2e1 * t38 * t33 * t231 * t233 - 0.2e1 * t38 * t223 * t237 + t38 * t97 * t260 / 0.3e1 + 0.6e1 * t38 * t8 * t502 * t504 - 0.9e1 * t38 * t232 * t508 + 0.3e1 / 0.2e1 * t38 * t103 * t512 + 0.2e1 * t38 * t103 * t516 - t38 * t41 * t543 / 0.4e1;
331   t603 = 0.11e2 / 0.54e2 * t482 - t489 / 0.12e2 + 0.8e1 / 0.3e1 * t169 * t246 * t334 * t50 + t486 / 0.6e1 - 0.11e2 / 0.648e3 * t359 * t426 * t2 * t15 - 0.36e2 * t371 * t508 + 0.8e1 * t183 * t516 + 0.185e3 / 0.864e3 * t478 + 0.16e2 / 0.3e1 * t74 * t535 - t78 * t543 + 0.24e2 * t28 * t502 * t504 + 0.6e1 * t183 * t512;
332   t644 = t5 * t455;
333   t658 = 0.4e1 * t27 * t231 * t45 * t46 * t33 * t233 + 0.8e1 / 0.3e1 * t341 * t46 * t92 * t104 - 0.16e2 / 0.9e1 * t345 * t114 * t253 * t50 + 0.56e2 / 0.27e2 * t170 * t46 * t214 * t50 - 0.4e1 / 0.3e1 * t340 * t112 * t114 * t117 * t104 + 0.2e1 / 0.3e1 * t345 * t114 * t117 * t124 - 0.4e1 / 0.3e1 * t170 * t46 * t92 * t124 + 0.2e1 / 0.3e1 * t170 * t46 * t33 * t260 - 0.5e1 / 0.9e1 * t74 * t527 * t530 * t1 * t46 - 0.160e3 / 0.81e2 * t175 * t176 * t644 + 0.140e3 / 0.81e2 * t75 * t520 - t359 * t360 * t40 * t50 / 0.108e3 - 0.4e1 * t341 * t46 * t171 * t124;
334   t665 = 0.18802527823203489287e2 * t464 * t455 * t322 * t459 - 0.60995923043525943529e3 / t317 / t141 * t455 / t472 * t459 + 0.10363566666666666667e-1 * t547 * t55 * t57 * t60 + 0.31090700000000000001e-1 * t266 * t138 - 0.12725038616008566457e1 * t333 * t2 * t534 * t69 + 0.53020994233369026904e-1 * t295 * t2 * t534 * t69 * t111 - 0.46822652193836017178e1 * t304 * t2 * t534 * t164 * t111 + 0.10126265704850360339e3 * t318 * t2 * t534 * t322 * t111 + 0.29070683134633122642e-2 * t378 * t50 + 0.29070683134633122642e-2 * t188 * t124 + 0.96902277115443742139e-3 * t82 * t260 + 0.96902277115443742139e-3 * (t603 + t658) * t81 * t14 - 0.43067678717974996507e-3 * t404 * t366;
335   t676 = t57 * t5;
336   t722 = 0.32300759038481247380e-3 * t395 * t403 * t178 + 0.32300759038481247380e-3 * t193 * t50 * t111 * t178 + 0.50245625170970829257e-3 * t195 * t241 + 0.44184161861140855756e-2 * t296 * t676 * t455 * t2 * t299 + 0.53834598397468745633e-4 * t80 / t327 / t27 * t14 * t55 * t114 * t253 * t44 + 0.13255248558342256727e0 * t64 * t526 * t2 * t530 * t69 * t1 * t46 + 0.48451138557721871070e-3 * t377 * t192 * t194 * t42 + 0.96902277115443742140e-3 * t395 * t399 * t42 + 0.48451138557721871070e-3 * t193 * t124 * t44 * t42 - 0.64601518076962494760e-3 * t396 * t108 - 0.64601518076962494760e-3 * t400 * t108 + 0.67293247996835932044e-5 * t329 * t14 * t1 * t277 * t6 * t360 - 0.16200859349084980444e-1 * t311 * t427 * t312;
337   t734 = t644 * t307;
338   t738 = t644 * t323;
339   t773 = -0.14804514559303905049e-3 * t330 * t241 + 0.47129772651883579471e0 * t149 * t114 * t455 * t69 - 0.41238551070398132037e0 * t66 * t46 * t426 * t69 - 0.12420094242202632684e0 * t305 * t57 * t734 + 0.26860754033147841838e1 * t319 * t57 * t738 + 0.17168305804406539632e2 * t306 * t734 - 0.37129640917784654576e3 * t320 * t738 - 0.19441031218901976532e0 * t297 * t644 * t299 + 0.17673664744456342301e-1 * t160 * t2 * t455 * t69 * t55 * t676 - 0.23411326096918008589e1 * t453 * t2 * t455 * t164 * t55 * t676 + 0.67508438032335735592e2 * t464 * t2 * t455 * t322 * t55 * t676 + 0.12112784639430467768e-3 * t187 * t328 * t14 * t108 + 0.12112784639430467768e-3 * t329 * t50 * t108;
340   t780 = 0.38760910846177496856e-2 * t379 + 0.77521821692354993712e-2 * t381 + 0.38760910846177496856e-2 * t383 + 0.12920303615392498952e-2 * t397 + 0.12920303615392498952e-2 * t401 + 0.43067678717974996507e-3 * t405 - 0.86135357435949993013e-3 * t385 + 0.63625193080042832287e0 * t336 + rho[0] * (t462 + t665 + t722 + t773) + 0.10604198846673805381e0 * t301 - 0.93645304387672034355e1 * t309 + 0.88368323722281711510e-2 * t314;
341   t794 = 0.20252531409700720677e3 * t325 + 0.16150379519240623690e-3 * t331 - 0.42416795386695221524e0 * t389 + 0.49486261284477758446e0 * t393 + 0.41454266666666666668e-1 * t267 + 0.82908533333333333336e-1 * t269 + 0.27636177777777777779e-1 * t272 + 0.41454266666666666668e-1 * t275 - 0.10645795064745113729e0 * t280 + 0.73578453447456598423e1 * t282 + 0.27636177777777777779e-1 * t284 - 0.92120592592592592596e-2 * t290 - 0.19441031218901976532e0 * t293;
342   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
343     v4rho4[0] = t780 + t794;
344 
345 #ifndef XC_DONT_COMPILE_MXC
346 
347   if(order < 5) return;
348 
349 
350 #endif
351 
352 #endif
353 
354 #endif
355 
356 #endif
357 
358 #endif
359 
360 
361 }
362 
363 
364 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))365 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
366 {
367 
368 #ifndef XC_DONT_COMPILE_EXC
369   double t1, t2, t3, t4, t5, t6, t7, t8;
370   double t9, t10, t12, t14, t15, t19, t20, t21;
371   double t24, t25, t27, t28, t30, t31;
372 
373 #ifndef XC_DONT_COMPILE_VXC
374   double t33, t34, t38, t39, t40, t41, t42, t44;
375   double t45, t46, t50, t55, t57, t58, t59, t60;
376   double t61, t63, t64, t66, t68, t69, t72, t74;
377   double t75, t78, t80, t81, t82, t83;
378 
379 #ifndef XC_DONT_COMPILE_FXC
380   double t90, t92, t93, t95, t97, t102, t103, t104;
381   double t108, t111, t112, t113, t114, t115, t117, t124;
382   double t130, t131, t134, t135, t136, t138, t139, t141;
383   double t143, t144, t146, t149, t152, t156, t158, t160;
384   double t161, t162, t163, t164, t166, t169, t170, t171;
385   double t175, t176, t177, t178, t183, t187, t188, t189;
386   double t192, t193, t194, t195, t196, t198;
387 
388 #ifndef XC_DONT_COMPILE_KXC
389   double t212, t214, t215, t217, t219, t221, t223, t230;
390   double t231, t232, t233, t237, t241, t246, t247, t248;
391   double t249, t250, t253, t260, t266, t267, t269, t271;
392   double t272, t274, t275, t277, t278, t279, t280, t282;
393   double t284, t288, t289, t290, t293, t295, t296, t297;
394   double t298, t299, t301, t304, t305, t306, t307, t309;
395   double t311, t312, t314, t316, t317, t318, t319, t320;
396   double t322, t323, t325, t327, t328, t329, t330, t331;
397   double t333, t334, t336, t340, t341, t345, t359, t360;
398   double t366, t371, t377, t378, t379, t381, t383, t385;
399   double t389, t393, t395, t396, t397, t399, t400, t401;
400   double t403, t404, t405, t407;
401 
402 #ifndef XC_DONT_COMPILE_LXC
403   double t418, t419, t420, t445, t447, t451, t454, t456;
404   double t464, t470, t474, t478, t481, t494, t496, t500;
405   double t504, t508, t512, t518, t519, t520, t522, t526;
406   double t527, t535, t539, t589, t630, t644, t658, t664;
407   double t714, t726, t730, t765, t780, t794;
408 #endif
409 
410 #endif
411 
412 #endif
413 
414 #endif
415 
416 #endif
417 
418 
419 
420   t1 = M_CBRT3;
421   t2 = 0.1e1 / M_PI;
422   t3 = POW_1_3(t2);
423   t4 = t1 * t3;
424   t5 = M_CBRT4;
425   t6 = t5 * t5;
426   t7 = POW_1_3(rho[0]);
427   t8 = 0.1e1 / t7;
428   t9 = t6 * t8;
429   t10 = t4 * t9;
430   t12 = sqrt(t10);
431   t14 = t10 / 0.4e1 + 0.35302100000000000000e1 * t12 + 0.180578e2;
432   t15 = 0.1e1 / t14;
433   t19 = log(t4 * t9 * t15 / 0.4e1);
434   t20 = 0.1554535e-1 * t19;
435   t21 = t12 + 0.706042e1;
436   t24 = atan(0.47309269095601128300e1 / t21);
437   t25 = 0.52491393169780936218e-1 * t24;
438   t27 = t12 / 0.2e1 + 0.32500e0;
439   t28 = t27 * t27;
440   t30 = log(t28 * t15);
441   t31 = 0.22478670955426118383e-2 * t30;
442   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
443     zk[0] = t20 + t25 + t31;
444 
445 #ifndef XC_DONT_COMPILE_VXC
446 
447   if(order < 1) return;
448 
449 
450   t33 = 0.1e1 / t7 / rho[0];
451   t34 = t6 * t33;
452   t38 = t4 * t6;
453   t39 = t14 * t14;
454   t40 = 0.1e1 / t39;
455   t41 = t8 * t40;
456   t42 = t4 * t34;
457   t44 = 0.1e1 / t12;
458   t45 = t44 * t1;
459   t46 = t3 * t6;
460   t50 = -t42 / 0.12e2 - 0.58836833333333333333e0 * t45 * t46 * t33;
461   t55 = t1 * t1;
462   t57 = 0.1e1 / t3;
463   t58 = (-t4 * t34 * t15 / 0.12e2 - t38 * t41 * t50 / 0.4e1) * t55 * t57;
464   t59 = t5 * t7;
465   t60 = t59 * t14;
466   t61 = t58 * t60;
467   t63 = t21 * t21;
468   t64 = 0.1e1 / t63;
469   t66 = t64 * t44 * t1;
470   t68 = 0.22381669423600000000e2 * t64 + 0.1e1;
471   t69 = 0.1e1 / t68;
472   t72 = t66 * t46 * t33 * t69;
473   t74 = t27 * t15;
474   t75 = t74 * t44;
475   t78 = t28 * t40;
476   t80 = -t75 * t42 / 0.6e1 - t78 * t50;
477   t81 = 0.1e1 / t28;
478   t82 = t80 * t81;
479   t83 = t82 * t14;
480   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
481     vrho[0] = t20 + t25 + t31 + rho[0] * (0.51817833333333333333e-2 * t61 + 0.41388824077869423261e-1 * t72 + 0.22478670955426118383e-2 * t83);
482 
483 #ifndef XC_DONT_COMPILE_FXC
484 
485   if(order < 2) return;
486 
487 
488   t90 = rho[0] * rho[0];
489   t92 = 0.1e1 / t7 / t90;
490   t93 = t6 * t92;
491   t95 = t4 * t93 * t15;
492   t97 = t33 * t40;
493   t102 = 0.1e1 / t39 / t14;
494   t103 = t8 * t102;
495   t104 = t50 * t50;
496   t108 = t4 * t93;
497   t111 = 0.1e1 / t12 / t10;
498   t112 = t111 * t55;
499   t113 = t3 * t3;
500   t114 = t113 * t5;
501   t115 = t7 * t7;
502   t117 = 0.1e1 / t115 / t90;
503   t124 = t108 / 0.9e1 - 0.39224555555555555555e0 * t112 * t114 * t117 + 0.78449111111111111110e0 * t45 * t46 * t92;
504   t130 = (t95 / 0.9e1 + t38 * t97 * t50 / 0.6e1 + t38 * t103 * t104 / 0.2e1 - t38 * t41 * t124 / 0.4e1) * t55 * t57;
505   t131 = t130 * t60;
506   t134 = t5 / t115;
507   t135 = t134 * t14;
508   t136 = t58 * t135;
509   t138 = t59 * t50;
510   t139 = t58 * t138;
511   t141 = t63 * t21;
512   t143 = 0.1e1 / t141 * t1;
513   t144 = t143 * t3;
514   t146 = t144 * t93 * t69;
515   t149 = t64 * t111 * t55;
516   t152 = t149 * t114 * t117 * t69;
517   t156 = t66 * t46 * t92 * t69;
518   t158 = t63 * t63;
519   t160 = 0.1e1 / t158 / t21;
520   t161 = t160 * t1;
521   t162 = t161 * t3;
522   t163 = t68 * t68;
523   t164 = 0.1e1 / t163;
524   t166 = t162 * t93 * t164;
525   t169 = t27 * t40;
526   t170 = t169 * t45;
527   t171 = t33 * t50;
528   t175 = t74 * t111;
529   t176 = t55 * t113;
530   t177 = t5 * t117;
531   t178 = t176 * t177;
532   t183 = t28 * t102;
533   t187 = t95 / 0.72e2 + t170 * t46 * t171 / 0.3e1 - t175 * t178 / 0.9e1 + 0.2e1 / 0.9e1 * t75 * t108 + 0.2e1 * t183 * t104 - t78 * t124;
534   t188 = t187 * t81;
535   t189 = t188 * t14;
536   t192 = 0.1e1 / t28 / t27;
537   t193 = t80 * t192;
538   t194 = t14 * t44;
539   t195 = t193 * t194;
540   t196 = t195 * t42;
541   t198 = t82 * t50;
542   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
543     v2rho2[0] = 0.10363566666666666667e-1 * t61 + 0.82777648155738846522e-1 * t72 + 0.44957341910852236766e-2 * t83 + rho[0] * (0.51817833333333333333e-2 * t131 + 0.17272611111111111111e-2 * t136 + 0.51817833333333333333e-2 * t139 + 0.13796274692623141087e-1 * t146 + 0.27592549385246282174e-1 * t152 - 0.55185098770492564348e-1 * t156 - 0.30878365944746984535e0 * t166 + 0.22478670955426118383e-2 * t189 + 0.37464451592376863972e-3 * t196 + 0.22478670955426118383e-2 * t198);
544 
545 #ifndef XC_DONT_COMPILE_KXC
546 
547   if(order < 3) return;
548 
549 
550   t212 = t90 * rho[0];
551   t214 = 0.1e1 / t7 / t212;
552   t215 = t6 * t214;
553   t217 = t4 * t215 * t15;
554   t219 = t92 * t40;
555   t221 = t38 * t219 * t50;
556   t223 = t33 * t102;
557   t230 = t39 * t39;
558   t231 = 0.1e1 / t230;
559   t232 = t8 * t231;
560   t233 = t104 * t50;
561   t237 = t50 * t124;
562   t241 = t4 * t215;
563   t246 = 0.1e1 / t12 / t176 / t134 / 0.4e1;
564   t247 = t246 * t2;
565   t248 = t90 * t90;
566   t249 = 0.1e1 / t248;
567   t250 = t247 * t249;
568   t253 = 0.1e1 / t115 / t212;
569   t260 = -0.7e1 / 0.27e2 * t241 - 0.23534733333333333333e1 * t250 + 0.15689822222222222222e1 * t112 * t114 * t253 - 0.18304792592592592592e1 * t45 * t46 * t214;
570   t266 = (-0.7e1 / 0.27e2 * t217 - t221 / 0.3e1 - t38 * t223 * t104 / 0.2e1 + t38 * t97 * t124 / 0.4e1 - 0.3e1 / 0.2e1 * t38 * t232 * t233 + 0.3e1 / 0.2e1 * t38 * t103 * t237 - t38 * t41 * t260 / 0.4e1) * t55 * t57;
571   t267 = t266 * t60;
572   t269 = t130 * t138;
573   t271 = t134 * t50;
574   t272 = t58 * t271;
575   t274 = t59 * t124;
576   t275 = t58 * t274;
577   t277 = 0.1e1 / t113;
578   t278 = t161 * t277;
579   t279 = t215 * t164;
580   t280 = t278 * t279;
581   t282 = t162 * t279;
582   t284 = t130 * t135;
583   t288 = t5 / t115 / rho[0];
584   t289 = t288 * t14;
585   t290 = t58 * t289;
586   t293 = t144 * t215 * t69;
587   t295 = 0.1e1 / t158;
588   t296 = t295 * t55;
589   t297 = t296 * t113;
590   t298 = t5 * t253;
591   t299 = t69 * t44;
592   t301 = t297 * t298 * t299;
593   t304 = 0.1e1 / t158 / t63;
594   t305 = t304 * t55;
595   t306 = t305 * t113;
596   t307 = t164 * t44;
597   t309 = t306 * t298 * t307;
598   t311 = t143 * t277;
599   t312 = t2 * t69;
600   t314 = t311 * t215 * t312;
601   t316 = 0.51817833333333333333e-2 * t267 + 0.10363566666666666667e-1 * t269 + 0.34545222222222222222e-2 * t272 + 0.51817833333333333333e-2 * t275 - 0.16381481915689750931e-1 * t280 + 0.11322067513073894330e1 * t282 + 0.34545222222222222222e-2 * t284 - 0.11515074074074074074e-2 * t290 - 0.50586340539618183985e-1 * t293 + 0.27592549385246282174e-1 * t301 - 0.14409904107548592783e1 * t309 + 0.22993791154371901812e-2 * t314;
602   t317 = t158 * t158;
603   t318 = 0.1e1 / t317;
604   t319 = t318 * t55;
605   t320 = t319 * t113;
606   t322 = 0.1e1 / t163 / t68;
607   t323 = t322 * t44;
608   t325 = t320 * t298 * t323;
609   t327 = t28 * t28;
610   t328 = 0.1e1 / t327;
611   t329 = t80 * t328;
612   t330 = t329 * t14;
613   t331 = t330 * t108;
614   t333 = t64 * t246;
615   t334 = t2 * t249;
616   t336 = t333 * t334 * t69;
617   t340 = t27 * t102;
618   t341 = t340 * t45;
619   t345 = t169 * t112;
620   t359 = t1 * t277 * t6;
621   t360 = t214 * t2;
622   t366 = t176 * t298;
623   t371 = t28 * t231;
624   t377 = -0.11e2 / 0.216e3 * t217 - t221 / 0.24e2 - t341 * t46 * t33 * t104 + t345 * t114 * t117 * t50 / 0.3e1 - 0.2e1 / 0.3e1 * t170 * t46 * t92 * t50 + t170 * t46 * t33 * t124 / 0.2e1 + t359 * t360 * t15 / 0.432e3 - 0.2e1 / 0.3e1 * t74 * t250 + 0.4e1 / 0.9e1 * t175 * t366 - 0.14e2 / 0.27e2 * t75 * t241 - 0.6e1 * t371 * t233 + 0.6e1 * t183 * t237 - t78 * t260;
625   t378 = t377 * t81;
626   t379 = t378 * t14;
627   t381 = t188 * t50;
628   t383 = t82 * t124;
629   t385 = t195 * t108;
630   t389 = t149 * t114 * t253 * t69;
631   t393 = t66 * t46 * t214 * t69;
632   t395 = t187 * t192;
633   t396 = t395 * t194;
634   t397 = t396 * t42;
635   t399 = t50 * t44;
636   t400 = t193 * t399;
637   t401 = t400 * t42;
638   t403 = t14 * t111;
639   t404 = t193 * t403;
640   t405 = t404 * t178;
641   t407 = 0.18429583437767336288e2 * t325 + 0.93661128980942159930e-4 * t331 + 0.16555529631147769304e0 * t336 + 0.22478670955426118383e-2 * t379 + 0.44957341910852236766e-2 * t381 + 0.22478670955426118383e-2 * t383 - 0.49952602123169151963e-3 * t385 - 0.11037019754098512870e0 * t389 + 0.12876523046448265015e0 * t393 + 0.74928903184753727944e-3 * t397 + 0.74928903184753727944e-3 * t401 + 0.24976301061584575981e-3 * t405;
642   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
643     v3rho3[0] = 0.15545350000000000000e-1 * t131 + 0.51817833333333333334e-2 * t136 + 0.15545350000000000000e-1 * t139 + 0.41388824077869423261e-1 * t146 + 0.82777648155738846522e-1 * t152 - 0.16555529631147769305e0 * t156 - 0.92635097834240953604e0 * t166 + 0.67436012866278355149e-2 * t189 + 0.11239335477713059192e-2 * t196 + 0.67436012866278355149e-2 * t198 + rho[0] * (t316 + t407);
644 
645 #ifndef XC_DONT_COMPILE_LXC
646 
647   if(order < 4) return;
648 
649 
650   t418 = 0.1e1 / t7 / t248;
651   t419 = t6 * t418;
652   t420 = t419 * t164;
653   t445 = 0.1e1 / t158 / t141;
654   t447 = 0.1e1 / t115 / t248;
655   t451 = t55 * t57 * t5;
656   t454 = 0.51817833333333333333e-2 * t266 * t135 + 0.10363566666666666667e-1 * t130 * t271 - 0.34545222222222222222e-2 * t58 * t288 * t50 + 0.12013086738172484016e0 * t278 * t420 - 0.47604147498151601160e1 * t162 * t420 - 0.34545222222222222222e-2 * t130 * t289 + 0.19191790123456790123e-2 * t58 * t177 * t14 + 0.21269256817794009175e0 * t144 * t419 * t69 + 0.15545350000000000000e-1 * t130 * t274 + 0.14985780636950745589e-2 * t193 * t14 * t250 + 0.51817833333333333333e-2 * t58 * t134 * t124 + 0.51817833333333333333e-2 * t58 * t59 * t260 - 0.65525927662759003726e-1 * t445 * t447 * t164 * t451;
657   t456 = 0.1e1 / t317 / t21;
658   t464 = t163 * t163;
659   t470 = t4 * t419 * t15;
660   t474 = t38 * t214 * t40 * t50;
661   t478 = t38 * t92 * t102 * t104;
662   t481 = t38 * t219 * t124;
663   t494 = 0.1e1 / t230 / t14;
664   t496 = t104 * t104;
665   t500 = t104 * t124;
666   t504 = t124 * t124;
667   t508 = t50 * t260;
668   t512 = t4 * t419;
669   t518 = 0.1e1 / t12 / t2 * rho[0] / 0.48e2;
670   t519 = t518 * t2;
671   t520 = t248 * rho[0];
672   t522 = 0.1e1 / t7 / t520;
673   t526 = 0.1e1 / t520;
674   t527 = t247 * t526;
675   t535 = 0.70e2 / 0.81e2 * t512 - 0.19612277777777777777e1 * t519 * t522 * t38 + 0.18827786666666666666e2 * t527 - 0.69732543209876543208e1 * t112 * t114 * t447 + 0.61015975308641975307e1 * t45 * t46 * t418;
676   t539 = 0.70e2 / 0.81e2 * t470 + 0.28e2 / 0.27e2 * t474 + 0.4e1 / 0.3e1 * t478 - 0.2e1 / 0.3e1 * t481 + 0.2e1 * t38 * t33 * t231 * t233 - 0.2e1 * t38 * t223 * t237 + t38 * t97 * t260 / 0.3e1 + 0.6e1 * t38 * t8 * t494 * t496 - 0.9e1 * t38 * t232 * t500 + 0.3e1 / 0.2e1 * t38 * t103 * t504 + 0.2e1 * t38 * t103 * t508 - t38 * t41 * t535 / 0.4e1;
677   t589 = 0.11e2 / 0.54e2 * t474 - t481 / 0.12e2 + 0.8e1 / 0.3e1 * t169 * t246 * t334 * t50 + t478 / 0.6e1 - 0.11e2 / 0.648e3 * t359 * t418 * t2 * t15 - 0.36e2 * t371 * t500 + 0.8e1 * t183 * t508 + 0.185e3 / 0.864e3 * t470 + 0.16e2 / 0.3e1 * t74 * t527 - t78 * t535 + 0.24e2 * t28 * t494 * t496 + 0.6e1 * t183 * t504;
678   t630 = t5 * t447;
679   t644 = 0.4e1 * t27 * t231 * t45 * t46 * t33 * t233 + 0.8e1 / 0.3e1 * t341 * t46 * t92 * t104 - 0.16e2 / 0.9e1 * t345 * t114 * t253 * t50 + 0.56e2 / 0.27e2 * t170 * t46 * t214 * t50 - 0.4e1 / 0.3e1 * t340 * t112 * t114 * t117 * t104 + 0.2e1 / 0.3e1 * t345 * t114 * t117 * t124 - 0.4e1 / 0.3e1 * t170 * t46 * t92 * t124 + 0.2e1 / 0.3e1 * t170 * t46 * t33 * t260 - 0.5e1 / 0.9e1 * t74 * t519 * t522 * t1 * t46 - 0.160e3 / 0.81e2 * t175 * t176 * t630 + 0.140e3 / 0.81e2 * t75 * t512 - t359 * t360 * t40 * t50 / 0.108e3 - 0.4e1 * t341 * t46 * t171 * t124;
680   t658 = 0.68440383742387934953e1 * t456 * t447 * t322 * t451 - 0.13129800378398182140e3 / t317 / t141 * t447 / t464 * t451 + 0.51817833333333333333e-2 * t539 * t55 * t57 * t60 + 0.15545350000000000000e-1 * t266 * t138 - 0.13244423704918215444e1 * t333 * t2 * t526 * t69 + 0.55185098770492564348e-1 * t295 * t2 * t526 * t69 * t111 - 0.28819808215097185566e1 * t304 * t2 * t526 * t164 * t111 + 0.36859166875534672576e2 * t318 * t2 * t526 * t322 * t111 + 0.22478670955426118383e-2 * (t589 + t644) * t81 * t14 + 0.67436012866278355149e-2 * t378 * t50 + 0.67436012866278355149e-2 * t188 * t124 + 0.22478670955426118383e-2 * t82 * t260 + 0.74928903184753727944e-3 * t395 * t403 * t178;
681   t664 = t57 * t5;
682   t714 = 0.74928903184753727944e-3 * t193 * t50 * t111 * t178 + 0.45987582308743803624e-2 * t296 * t664 * t447 * t2 * t299 + 0.12488150530792287991e-3 * t80 / t327 / t27 * t14 * t55 * t114 * t253 * t44 + 0.13796274692623141087e0 * t64 * t518 * t2 * t522 * t69 * t1 * t46 + 0.11239335477713059192e-2 * t377 * t192 * t194 * t42 + 0.22478670955426118383e-2 * t395 * t399 * t42 + 0.11239335477713059192e-2 * t193 * t124 * t44 * t42 - 0.14985780636950745589e-2 * t396 * t108 - 0.14985780636950745589e-2 * t400 * t108 - 0.99905204246338303925e-3 * t404 * t366 + 0.15610188163490359988e-4 * t329 * t14 * t1 * t277 * t6 * t360 + 0.11655607162072802125e-2 * t195 * t241 - 0.16862113513206061329e-1 * t311 * t419 * t312;
683   t726 = t630 * t307;
684   t730 = t630 * t323;
685   t765 = -0.34342413959678791975e-3 * t330 * t241 + 0.49053421129326723867e0 * t149 * t114 * t447 * t69 - 0.42921743488160883383e0 * t66 * t46 * t418 * t69 - 0.76446915606552171013e-1 * t305 * t57 * t726 + 0.97771976774839907070e0 * t319 * t57 * t730 + 0.10567263012202301374e2 * t306 * t726 - 0.13515027854362713278e3 * t320 * t730 - 0.20234536215847273594e0 * t297 * t630 * t299 + 0.18395032923497521449e-1 * t160 * t2 * t447 * t69 * t55 * t664 - 0.14409904107548592783e1 * t445 * t2 * t447 * t164 * t55 * t664 + 0.24572777917023115051e2 * t456 * t2 * t447 * t322 * t55 * t664 + 0.28098338694282647979e-3 * t187 * t328 * t14 * t108 + 0.28098338694282647979e-3 * t329 * t50 * t108;
686   t780 = rho[0] * (t454 + t658 + t714 + t765) + 0.89914683821704473532e-2 * t379 + 0.17982936764340894707e-1 * t381 + 0.89914683821704473532e-2 * t383 + 0.66222118524591077217e0 * t336 + 0.29971561273901491178e-2 * t401 + 0.99905204246338303928e-3 * t405 - 0.19981040849267660785e-2 * t385 + 0.29971561273901491178e-2 * t397 + 0.11037019754098512870e0 * t301 - 0.57639616430194371132e1 * t309 + 0.91975164617487607247e-2 * t314;
687   t794 = 0.73718333751069345151e2 * t325 + 0.37464451592376863973e-3 * t331 - 0.44148079016394051479e0 * t389 + 0.51506092185793060060e0 * t393 + 0.20727133333333333333e-1 * t267 + 0.41454266666666666667e-1 * t269 + 0.13818088888888888889e-1 * t272 + 0.20727133333333333333e-1 * t275 - 0.65525927662759003724e-1 * t280 + 0.45288270052295577318e1 * t282 + 0.13818088888888888889e-1 * t284 - 0.46060296296296296297e-2 * t290 - 0.20234536215847273594e0 * t293;
688   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
689     v4rho4[0] = t780 + t794;
690 
691 #ifndef XC_DONT_COMPILE_MXC
692 
693   if(order < 5) return;
694 
695 
696 #endif
697 
698 #endif
699 
700 #endif
701 
702 #endif
703 
704 #endif
705 
706 
707 }
708 
709 
710 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))711 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
712 {
713 
714 #ifndef XC_DONT_COMPILE_EXC
715   double t1, t2, t3, t4, t5, t6, t7, t8;
716   double t9, t10, t11, t12, t13, t15, t16, t20;
717   double t22, t25, t27, t28, t29, t31, t33, t34;
718   double t35, t36, t37, t38, t40, t41, t43, t44;
719   double t47, t49, t50, t52, t53, t57, t59, t62;
720   double t64, t65, t67, t69, t71;
721 
722 #ifndef XC_DONT_COMPILE_VXC
723   double t73, t74, t78, t79, t80, t81, t82, t83;
724   double t84, t85, t86, t88, t90, t95, t97, t98;
725   double t99, t100, t103, t104, t106, t108, t109, t114;
726   double t115, t118, t120, t121, t122, t125, t126, t127;
727   double t128, t129, t130, t132, t135, t137, t141, t142;
728   double t143, t145, t151, t152, t155, t156, t158, t160;
729   double t161, t166, t167, t170, t172, t173, t174, t177;
730   double t179, t181, t184, t186, t189, t191, t193;
731 
732 #ifndef XC_DONT_COMPILE_FXC
733   double t196, t198, t201, t202, t204, t206, t211, t212;
734   double t213, t217, t218, t220, t221, t222, t223, t224;
735   double t226, t228, t231, t233, t239, t243, t244, t247;
736   double t250, t252, t253, t258, t267, t269, t270, t271;
737   double t272, t273, t278, t279, t280, t284, t285, t286;
738   double t287, t292, t296, t297, t301, t302, t303, t304;
739   double t309, t310, t312, t313, t314, t315, t316, t319;
740   double t320, t321, t323, t326, t327, t328, t331, t334;
741   double t336, t338, t340, t345, t346, t347, t353, t359;
742   double t362, t365, t368, t370, t371, t376, t385, t387;
743   double t388, t389, t390, t391, t396, t397, t398, t402;
744   double t407, t411, t412, t416, t417, t418, t419, t424;
745   double t426, t428, t429, t431, t435, t436, t439, t442;
746   double t445, t448, t450, t452, t454, t459, t460, t464;
747   double t467, t470, t473, t475, t476, t478;
748 
749 #ifndef XC_DONT_COMPILE_KXC
750   double t481, t484, t491, t492, t496, t497, t501, t502;
751   double t503, t505, t506, t507, t512, t513, t514, t515;
752   double t519, t520, t521, t525, t526, t527, t528, t530;
753   double t531, t535, t536, t537, t538, t551, t552, t555;
754   double t556, t557, t560, t561, t565, t567, t569, t571;
755   double t578, t579, t580, t581, t585, t589, t590, t594;
756   double t595, t596, t597, t598, t601, t604, t606, t612;
757   double t617, t620, t623, t624, t631, t632, t636, t650;
758   double t651, t657, t662, t668, t669, t676, t677, t681;
759   double t682, t683, t685, t686, t688, t691, t692, t695;
760   double t698, t700, t704, t705, t708, t711, t714, t716;
761   double t718, t720, t722, t724, t731, t732, t733, t734;
762   double t738, t745, t751, t756, t769, t770, t773, t774;
763   double t777, t778, t781, t784, t785, t790, t793, t794;
764   double t798, t820, t826, t827, t834, t840, t846, t847;
765   double t848, t849, t854, t855, t856, t857, t861, t862;
766   double t866, t867, t868, t869, t871, t872, t876, t877;
767   double t878, t879, t882, t883, t885, t887, t888, t890;
768   double t893, t898, t899, t902, t904, t905, t906, t909;
769   double t919, t922, t932, t934, t937, t939, t940, t942;
770   double t943, t949, t950, t955, t960, t963, t968, t971;
771   double t974, t976, t979, t981, t982, t988, t990, t996;
772   double t999, t1004, t1007, t1009, t1010, t1013;
773 
774 #ifndef XC_DONT_COMPILE_LXC
775   double t1016, t1020, t1038, t1039, t1043, t1049, t1054, t1065;
776   double t1080, t1082, t1087, t1091, t1099, t1100, t1101, t1126;
777   double t1130, t1134, t1137, t1150, t1152, t1156, t1160, t1164;
778   double t1168, t1169, t1173, t1174, t1176, t1178, t1180, t1183;
779   double t1186, t1188, t1192, t1197, t1200, t1207, t1211, t1243;
780   double t1255, t1267, t1271, t1310, t1322, t1328, t1337, t1363;
781   double t1373, t1378, t1381, t1383, t1386, t1389, t1391, t1393;
782   double t1394, t1400, t1405, t1407, t1410, t1412, t1413, t1419;
783   double t1427, t1457, t1461, t1479, t1487, t1495, t1529, t1546;
784   double t1550, t1554, t1557, t1570, t1572, t1576, t1580, t1584;
785   double t1592, t1596, t1621, t1637, t1646, t1704, t1741, t1746;
786   double t1750, t1752, t1755, t1758, t1767, t1769, t1770, t1772;
787   double t1803, t1821, t1824, t1830, t1832, t1833, t1835, t1839;
788   double t1841, t1843, t1854, t1856, t1867, t1870, t1909, t1917;
789   double t1919, t1923, t1949, t1968, t1990, t1997, t2000, t2011;
790   double t2016, t2022, t2025, t2030, t2038;
791 #endif
792 
793 #endif
794 
795 #endif
796 
797 #endif
798 
799 #endif
800 
801 
802 
803   t1 = M_CBRT3;
804   t2 = 0.1e1 / M_PI;
805   t3 = POW_1_3(t2);
806   t4 = t1 * t3;
807   t5 = M_CBRT4;
808   t6 = t5 * t5;
809   t7 = rho[0] + rho[1];
810   t8 = POW_1_3(t7);
811   t9 = 0.1e1 / t8;
812   t10 = t6 * t9;
813   t11 = t4 * t10;
814   t12 = t11 / 0.4e1;
815   t13 = sqrt(t11);
816   t15 = t12 + 0.18637200000000000000e1 * t13 + 0.129352e2;
817   t16 = 0.1e1 / t15;
818   t20 = log(t4 * t10 * t16 / 0.4e1);
819   t22 = t13 + 0.372744e1;
820   t25 = atan(0.61519908197590802322e1 / t22);
821   t27 = t13 / 0.2e1;
822   t28 = t27 + 0.10498e0;
823   t29 = t28 * t28;
824   t31 = log(t29 * t16);
825   t33 = 0.310907e-1 * t20 + 0.38783294878113014393e-1 * t25 + 0.96902277115443742139e-3 * t31;
826   t34 = rho[0] - rho[1];
827   t35 = 0.1e1 / t7;
828   t36 = t34 * t35;
829   t37 = 0.1e1 + t36;
830   t38 = POW_1_3(t37);
831   t40 = 0.1e1 - t36;
832   t41 = POW_1_3(t40);
833   t43 = t38 * t37 + t41 * t40 - 0.2e1;
834   t44 = M_CBRT2;
835   t47 = 0.1e1 / (0.2e1 * t44 - 0.2e1);
836   t49 = -t43 * t47 + 0.1e1;
837   t50 = t33 * t49;
838   t52 = t12 + 0.35302100000000000000e1 * t13 + 0.180578e2;
839   t53 = 0.1e1 / t52;
840   t57 = log(t4 * t10 * t53 / 0.4e1);
841   t59 = t13 + 0.706042e1;
842   t62 = atan(0.47309269095601128300e1 / t59);
843   t64 = t27 + 0.32500e0;
844   t65 = t64 * t64;
845   t67 = log(t65 * t53);
846   t69 = 0.1554535e-1 * t57 + 0.52491393169780936218e-1 * t62 + 0.22478670955426118383e-2 * t67;
847   t71 = t69 * t43 * t47;
848   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
849     zk[0] = t50 + t71;
850 
851 #ifndef XC_DONT_COMPILE_VXC
852 
853   if(order < 1) return;
854 
855 
856   t73 = 0.1e1 / t8 / t7;
857   t74 = t6 * t73;
858   t78 = t4 * t6;
859   t79 = t15 * t15;
860   t80 = 0.1e1 / t79;
861   t81 = t9 * t80;
862   t82 = t4 * t74;
863   t83 = t82 / 0.12e2;
864   t84 = 0.1e1 / t13;
865   t85 = t84 * t1;
866   t86 = t3 * t6;
867   t88 = t85 * t86 * t73;
868   t90 = -t83 - 0.31062000000000000000e0 * t88;
869   t95 = t1 * t1;
870   t97 = 0.1e1 / t3;
871   t98 = (-t4 * t74 * t16 / 0.12e2 - t78 * t81 * t90 / 0.4e1) * t95 * t97;
872   t99 = t5 * t8;
873   t100 = t99 * t15;
874   t103 = t22 * t22;
875   t104 = 0.1e1 / t103;
876   t106 = t104 * t84 * t1;
877   t108 = 0.37846991046400000000e2 * t104 + 0.1e1;
878   t109 = 0.1e1 / t108;
879   t114 = t28 * t16;
880   t115 = t114 * t84;
881   t118 = t29 * t80;
882   t120 = -t115 * t82 / 0.6e1 - t118 * t90;
883   t121 = 0.1e1 / t29;
884   t122 = t120 * t121;
885   t125 = 0.10363566666666666667e-1 * t98 * t100 + 0.39765745675026770179e-1 * t106 * t86 * t73 * t109 + 0.96902277115443742139e-3 * t122 * t15;
886   t126 = t125 * t49;
887   t127 = t7 * t7;
888   t128 = 0.1e1 / t127;
889   t129 = t34 * t128;
890   t130 = t35 - t129;
891   t132 = -t130;
892   t135 = 0.4e1 / 0.3e1 * t38 * t130 + 0.4e1 / 0.3e1 * t41 * t132;
893   t137 = t33 * t135 * t47;
894   t141 = t52 * t52;
895   t142 = 0.1e1 / t141;
896   t143 = t9 * t142;
897   t145 = -t83 - 0.58836833333333333333e0 * t88;
898   t151 = (-t4 * t74 * t53 / 0.12e2 - t78 * t143 * t145 / 0.4e1) * t95 * t97;
899   t152 = t99 * t52;
900   t155 = t59 * t59;
901   t156 = 0.1e1 / t155;
902   t158 = t156 * t84 * t1;
903   t160 = 0.22381669423600000000e2 * t156 + 0.1e1;
904   t161 = 0.1e1 / t160;
905   t166 = t64 * t53;
906   t167 = t166 * t84;
907   t170 = t65 * t142;
908   t172 = -t167 * t82 / 0.6e1 - t170 * t145;
909   t173 = 0.1e1 / t65;
910   t174 = t172 * t173;
911   t177 = 0.51817833333333333333e-2 * t151 * t152 + 0.41388824077869423261e-1 * t158 * t86 * t73 * t161 + 0.22478670955426118383e-2 * t174 * t52;
912   t179 = t177 * t43 * t47;
913   t181 = t69 * t135 * t47;
914   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
915     vrho[0] = t50 + t71 + t7 * (t126 - t137 + t179 + t181);
916 
917   t184 = -t35 - t129;
918   t186 = -t184;
919   t189 = 0.4e1 / 0.3e1 * t38 * t184 + 0.4e1 / 0.3e1 * t41 * t186;
920   t191 = t33 * t189 * t47;
921   t193 = t69 * t189 * t47;
922   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
923     vrho[1] = t50 + t71 + t7 * (t126 - t191 + t179 + t193);
924 
925 #ifndef XC_DONT_COMPILE_FXC
926 
927   if(order < 2) return;
928 
929 
930   t196 = 0.2e1 * t126;
931   t198 = 0.2e1 * t179;
932   t201 = 0.1e1 / t8 / t127;
933   t202 = t6 * t201;
934   t204 = t4 * t202 * t16;
935   t206 = t73 * t80;
936   t211 = 0.1e1 / t79 / t15;
937   t212 = t9 * t211;
938   t213 = t90 * t90;
939   t217 = t4 * t202;
940   t218 = t217 / 0.9e1;
941   t220 = 0.1e1 / t13 / t11;
942   t221 = t220 * t95;
943   t222 = t3 * t3;
944   t223 = t222 * t5;
945   t224 = t8 * t8;
946   t226 = 0.1e1 / t224 / t127;
947   t228 = t221 * t223 * t226;
948   t231 = t85 * t86 * t201;
949   t233 = t218 - 0.20708000000000000000e0 * t228 + 0.41416000000000000000e0 * t231;
950   t239 = (t204 / 0.9e1 + t78 * t206 * t90 / 0.6e1 + t78 * t212 * t213 / 0.2e1 - t78 * t81 * t233 / 0.4e1) * t95 * t97;
951   t243 = t5 / t224;
952   t244 = t243 * t15;
953   t247 = t99 * t90;
954   t250 = t103 * t22;
955   t252 = 0.1e1 / t250 * t1;
956   t253 = t252 * t3;
957   t258 = t104 * t220 * t95;
958   t267 = t103 * t103;
959   t269 = 0.1e1 / t267 / t22;
960   t270 = t269 * t1;
961   t271 = t270 * t3;
962   t272 = t108 * t108;
963   t273 = 0.1e1 / t272;
964   t278 = t28 * t80;
965   t279 = t278 * t85;
966   t280 = t73 * t90;
967   t284 = t114 * t220;
968   t285 = t95 * t222;
969   t286 = t5 * t226;
970   t287 = t285 * t286;
971   t292 = t29 * t211;
972   t296 = t204 / 0.72e2 + t279 * t86 * t280 / 0.3e1 - t284 * t287 / 0.9e1 + 0.2e1 / 0.9e1 * t115 * t217 + 0.2e1 * t292 * t213 - t118 * t233;
973   t297 = t296 * t121;
974   t301 = 0.1e1 / t29 / t28;
975   t302 = t120 * t301;
976   t303 = t15 * t84;
977   t304 = t302 * t303;
978   t309 = 0.10363566666666666667e-1 * t239 * t100 + 0.34545222222222222223e-2 * t98 * t244 + 0.10363566666666666667e-1 * t98 * t247 + 0.13255248558342256726e-1 * t253 * t202 * t109 + 0.26510497116684513453e-1 * t258 * t223 * t226 * t109 - 0.53020994233369026905e-1 * t106 * t86 * t201 * t109 - 0.50167127350538589834e0 * t271 * t202 * t273 + 0.96902277115443742139e-3 * t297 * t15 + 0.16150379519240623690e-3 * t304 * t82 + 0.96902277115443742139e-3 * t122 * t90;
979   t310 = t309 * t49;
980   t312 = t125 * t135 * t47;
981   t313 = 0.2e1 * t312;
982   t314 = t38 * t38;
983   t315 = 0.1e1 / t314;
984   t316 = t130 * t130;
985   t319 = t127 * t7;
986   t320 = 0.1e1 / t319;
987   t321 = t34 * t320;
988   t323 = -0.2e1 * t128 + 0.2e1 * t321;
989   t326 = t41 * t41;
990   t327 = 0.1e1 / t326;
991   t328 = t132 * t132;
992   t331 = -t323;
993   t334 = 0.4e1 / 0.9e1 * t315 * t316 + 0.4e1 / 0.3e1 * t38 * t323 + 0.4e1 / 0.9e1 * t327 * t328 + 0.4e1 / 0.3e1 * t41 * t331;
994   t336 = t33 * t334 * t47;
995   t338 = t4 * t202 * t53;
996   t340 = t73 * t142;
997   t345 = 0.1e1 / t141 / t52;
998   t346 = t9 * t345;
999   t347 = t145 * t145;
1000   t353 = t218 - 0.39224555555555555555e0 * t228 + 0.78449111111111111110e0 * t231;
1001   t359 = (t338 / 0.9e1 + t78 * t340 * t145 / 0.6e1 + t78 * t346 * t347 / 0.2e1 - t78 * t143 * t353 / 0.4e1) * t95 * t97;
1002   t362 = t243 * t52;
1003   t365 = t99 * t145;
1004   t368 = t155 * t59;
1005   t370 = 0.1e1 / t368 * t1;
1006   t371 = t370 * t3;
1007   t376 = t156 * t220 * t95;
1008   t385 = t155 * t155;
1009   t387 = 0.1e1 / t385 / t59;
1010   t388 = t387 * t1;
1011   t389 = t388 * t3;
1012   t390 = t160 * t160;
1013   t391 = 0.1e1 / t390;
1014   t396 = t64 * t142;
1015   t397 = t396 * t85;
1016   t398 = t73 * t145;
1017   t402 = t166 * t220;
1018   t407 = t65 * t345;
1019   t411 = t338 / 0.72e2 + t397 * t86 * t398 / 0.3e1 - t402 * t287 / 0.9e1 + 0.2e1 / 0.9e1 * t167 * t217 + 0.2e1 * t407 * t347 - t170 * t353;
1020   t412 = t411 * t173;
1021   t416 = 0.1e1 / t65 / t64;
1022   t417 = t172 * t416;
1023   t418 = t52 * t84;
1024   t419 = t417 * t418;
1025   t424 = 0.51817833333333333333e-2 * t359 * t152 + 0.17272611111111111111e-2 * t151 * t362 + 0.51817833333333333333e-2 * t151 * t365 + 0.13796274692623141087e-1 * t371 * t202 * t161 + 0.27592549385246282174e-1 * t376 * t223 * t226 * t161 - 0.55185098770492564348e-1 * t158 * t86 * t201 * t161 - 0.30878365944746984535e0 * t389 * t202 * t391 + 0.22478670955426118383e-2 * t412 * t52 + 0.37464451592376863972e-3 * t419 * t82 + 0.22478670955426118383e-2 * t174 * t145;
1026   t426 = t424 * t43 * t47;
1027   t428 = t177 * t135 * t47;
1028   t429 = 0.2e1 * t428;
1029   t431 = t69 * t334 * t47;
1030   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1031     v2rho2[0] = t196 - 0.2e1 * t137 + t198 + 0.2e1 * t181 + t7 * (t310 - t313 - t336 + t426 + t429 + t431);
1032 
1033   t435 = t125 * t189 * t47;
1034   t436 = t315 * t184;
1035   t439 = t38 * t34;
1036   t442 = t327 * t186;
1037   t445 = t41 * t34;
1038   t448 = 0.4e1 / 0.9e1 * t436 * t130 + 0.8e1 / 0.3e1 * t439 * t320 + 0.4e1 / 0.9e1 * t442 * t132 - 0.8e1 / 0.3e1 * t445 * t320;
1039   t450 = t33 * t448 * t47;
1040   t452 = t177 * t189 * t47;
1041   t454 = t69 * t448 * t47;
1042   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1043     v2rho2[1] = t196 - t137 + t198 + t181 - t191 + t193 + t7 * (t310 - t312 - t435 - t450 + t426 + t428 + t452 + t454);
1044 
1045   t459 = 0.2e1 * t435;
1046   t460 = t184 * t184;
1047   t464 = 0.2e1 * t128 + 0.2e1 * t321;
1048   t467 = t186 * t186;
1049   t470 = -t464;
1050   t473 = 0.4e1 / 0.9e1 * t315 * t460 + 0.4e1 / 0.3e1 * t38 * t464 + 0.4e1 / 0.9e1 * t327 * t467 + 0.4e1 / 0.3e1 * t41 * t470;
1051   t475 = t33 * t473 * t47;
1052   t476 = 0.2e1 * t452;
1053   t478 = t69 * t473 * t47;
1054   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1055     v2rho2[2] = t196 - 0.2e1 * t191 + t198 + 0.2e1 * t193 + t7 * (t310 - t459 - t475 + t426 + t476 + t478);
1056 
1057 #ifndef XC_DONT_COMPILE_KXC
1058 
1059   if(order < 3) return;
1060 
1061 
1062   t481 = 0.3e1 * t310;
1063   t484 = 0.3e1 * t426;
1064   t491 = t5 / t224 / t7;
1065   t492 = t491 * t15;
1066   t496 = 0.1e1 / t8 / t319;
1067   t497 = t6 * t496;
1068   t501 = 0.1e1 / t267;
1069   t502 = t501 * t95;
1070   t503 = t502 * t222;
1071   t505 = 0.1e1 / t224 / t319;
1072   t506 = t5 * t505;
1073   t507 = t109 * t84;
1074   t512 = 0.1e1 / t267 / t103;
1075   t513 = t512 * t95;
1076   t514 = t513 * t222;
1077   t515 = t273 * t84;
1078   t519 = 0.1e1 / t222;
1079   t520 = t252 * t519;
1080   t521 = t2 * t109;
1081   t525 = t267 * t267;
1082   t526 = 0.1e1 / t525;
1083   t527 = t526 * t95;
1084   t528 = t527 * t222;
1085   t530 = 0.1e1 / t272 / t108;
1086   t531 = t530 * t84;
1087   t535 = t29 * t29;
1088   t536 = 0.1e1 / t535;
1089   t537 = t120 * t536;
1090   t538 = t537 * t15;
1091   t551 = t296 * t301;
1092   t552 = t551 * t303;
1093   t555 = 0.69090444444444444446e-2 * t239 * t244 - 0.23030148148148148149e-2 * t98 * t492 - 0.48602578047254941329e-1 * t253 * t497 * t109 + 0.26510497116684513452e-1 * t503 * t506 * t507 - 0.23411326096918008589e1 * t514 * t506 * t515 + 0.22092080930570427878e-2 * t520 * t497 * t521 + 0.50631328524251801694e2 * t528 * t506 * t531 + 0.40375948798101559225e-4 * t538 * t217 - 0.21533839358987498253e-3 * t304 * t217 - 0.10604198846673805381e0 * t258 * t223 * t505 * t109 + 0.12371565321119439611e0 * t106 * t86 * t496 * t109 + 0.32300759038481247380e-3 * t552 * t82;
1094   t556 = t90 * t84;
1095   t557 = t302 * t556;
1096   t560 = t15 * t220;
1097   t561 = t302 * t560;
1098   t565 = t4 * t497 * t16;
1099   t567 = t201 * t80;
1100   t569 = t78 * t567 * t90;
1101   t571 = t73 * t211;
1102   t578 = t79 * t79;
1103   t579 = 0.1e1 / t578;
1104   t580 = t9 * t579;
1105   t581 = t213 * t90;
1106   t585 = t90 * t233;
1107   t589 = t4 * t497;
1108   t590 = 0.7e1 / 0.27e2 * t589;
1109   t594 = 0.1e1 / t13 / t285 / t243 / 0.4e1;
1110   t595 = t594 * t2;
1111   t596 = t127 * t127;
1112   t597 = 0.1e1 / t596;
1113   t598 = t595 * t597;
1114   t601 = t221 * t223 * t505;
1115   t604 = t85 * t86 * t496;
1116   t606 = -t590 - 0.12424800000000000000e1 * t598 + 0.82832000000000000000e0 * t601 - 0.96637333333333333333e0 * t604;
1117   t612 = (-0.7e1 / 0.27e2 * t565 - t569 / 0.3e1 - t78 * t571 * t213 / 0.2e1 + t78 * t206 * t233 / 0.4e1 - 0.3e1 / 0.2e1 * t78 * t580 * t581 + 0.3e1 / 0.2e1 * t78 * t212 * t585 - t78 * t81 * t606 / 0.4e1) * t95 * t97;
1118   t617 = t243 * t90;
1119   t620 = t99 * t233;
1120   t623 = t270 * t519;
1121   t624 = t497 * t273;
1122   t631 = t28 * t211;
1123   t632 = t631 * t85;
1124   t636 = t278 * t221;
1125   t650 = t1 * t519 * t6;
1126   t651 = t496 * t2;
1127   t657 = t285 * t506;
1128   t662 = t29 * t579;
1129   t668 = -0.11e2 / 0.216e3 * t565 - t569 / 0.24e2 - t632 * t86 * t73 * t213 + t636 * t223 * t226 * t90 / 0.3e1 - 0.2e1 / 0.3e1 * t279 * t86 * t201 * t90 + t279 * t86 * t73 * t233 / 0.2e1 + t650 * t651 * t16 / 0.432e3 - 0.2e1 / 0.3e1 * t114 * t598 + 0.4e1 / 0.9e1 * t284 * t657 - 0.14e2 / 0.27e2 * t115 * t589 - 0.6e1 * t662 * t581 + 0.6e1 * t292 * t585 - t118 * t606;
1130   t669 = t668 * t121;
1131   t676 = t104 * t594;
1132   t677 = t2 * t597;
1133   t681 = 0.32300759038481247380e-3 * t557 * t82 + 0.10766919679493749127e-3 * t561 * t287 + 0.10363566666666666667e-1 * t612 * t100 + 0.20727133333333333334e-1 * t239 * t247 + 0.69090444444444444446e-2 * t98 * t617 + 0.10363566666666666667e-1 * t98 * t620 - 0.26614487661862784322e-1 * t623 * t624 + 0.18394613361864149606e1 * t271 * t624 + 0.96902277115443742139e-3 * t669 * t15 + 0.19380455423088748428e-2 * t297 * t90 + 0.96902277115443742139e-3 * t122 * t233 + 0.15906298270010708072e0 * t676 * t677 * t109;
1134   t682 = t555 + t681;
1135   t683 = t682 * t49;
1136   t685 = t309 * t135 * t47;
1137   t686 = 0.3e1 * t685;
1138   t688 = t125 * t334 * t47;
1139   t691 = 0.1e1 / t314 / t37;
1140   t692 = t316 * t130;
1141   t695 = t315 * t130;
1142   t698 = t34 * t597;
1143   t700 = 0.6e1 * t320 - 0.6e1 * t698;
1144   t704 = 0.1e1 / t326 / t40;
1145   t705 = t328 * t132;
1146   t708 = t327 * t132;
1147   t711 = -t700;
1148   t714 = -0.8e1 / 0.27e2 * t691 * t692 + 0.4e1 / 0.3e1 * t695 * t323 + 0.4e1 / 0.3e1 * t38 * t700 - 0.8e1 / 0.27e2 * t704 * t705 + 0.4e1 / 0.3e1 * t708 * t331 + 0.4e1 / 0.3e1 * t41 * t711;
1149   t716 = t33 * t714 * t47;
1150   t718 = t4 * t497 * t53;
1151   t720 = t201 * t142;
1152   t722 = t78 * t720 * t145;
1153   t724 = t73 * t345;
1154   t731 = t141 * t141;
1155   t732 = 0.1e1 / t731;
1156   t733 = t9 * t732;
1157   t734 = t347 * t145;
1158   t738 = t145 * t353;
1159   t745 = -t590 - 0.23534733333333333333e1 * t598 + 0.15689822222222222222e1 * t601 - 0.18304792592592592592e1 * t604;
1160   t751 = (-0.7e1 / 0.27e2 * t718 - t722 / 0.3e1 - t78 * t724 * t347 / 0.2e1 + t78 * t340 * t353 / 0.4e1 - 0.3e1 / 0.2e1 * t78 * t733 * t734 + 0.3e1 / 0.2e1 * t78 * t346 * t738 - t78 * t143 * t745 / 0.4e1) * t95 * t97;
1161   t756 = t243 * t145;
1162   t769 = t411 * t416;
1163   t770 = t769 * t418;
1164   t773 = t145 * t84;
1165   t774 = t417 * t773;
1166   t777 = t52 * t220;
1167   t778 = t417 * t777;
1168   t781 = t99 * t353;
1169   t784 = t388 * t519;
1170   t785 = t497 * t391;
1171   t790 = 0.51817833333333333333e-2 * t751 * t152 + 0.10363566666666666667e-1 * t359 * t365 + 0.34545222222222222222e-2 * t151 * t756 - 0.49952602123169151963e-3 * t419 * t217 - 0.11037019754098512870e0 * t376 * t223 * t505 * t161 + 0.12876523046448265015e0 * t158 * t86 * t496 * t161 + 0.74928903184753727944e-3 * t770 * t82 + 0.74928903184753727944e-3 * t774 * t82 + 0.24976301061584575981e-3 * t778 * t287 + 0.51817833333333333333e-2 * t151 * t781 - 0.16381481915689750931e-1 * t784 * t785 + 0.11322067513073894330e1 * t389 * t785;
1172   t793 = t64 * t345;
1173   t794 = t793 * t85;
1174   t798 = t396 * t221;
1175   t820 = t65 * t732;
1176   t826 = -0.11e2 / 0.216e3 * t718 - t722 / 0.24e2 - t794 * t86 * t73 * t347 + t798 * t223 * t226 * t145 / 0.3e1 - 0.2e1 / 0.3e1 * t397 * t86 * t201 * t145 + t397 * t86 * t73 * t353 / 0.2e1 + t650 * t651 * t53 / 0.432e3 - 0.2e1 / 0.3e1 * t166 * t598 + 0.4e1 / 0.9e1 * t402 * t657 - 0.14e2 / 0.27e2 * t167 * t589 - 0.6e1 * t820 * t734 + 0.6e1 * t407 * t738 - t170 * t745;
1177   t827 = t826 * t173;
1178   t834 = t156 * t594;
1179   t840 = t491 * t52;
1180   t846 = 0.1e1 / t385;
1181   t847 = t846 * t95;
1182   t848 = t847 * t222;
1183   t849 = t161 * t84;
1184   t854 = 0.1e1 / t385 / t155;
1185   t855 = t854 * t95;
1186   t856 = t855 * t222;
1187   t857 = t391 * t84;
1188   t861 = t370 * t519;
1189   t862 = t2 * t161;
1190   t866 = t385 * t385;
1191   t867 = 0.1e1 / t866;
1192   t868 = t867 * t95;
1193   t869 = t868 * t222;
1194   t871 = 0.1e1 / t390 / t160;
1195   t872 = t871 * t84;
1196   t876 = t65 * t65;
1197   t877 = 0.1e1 / t876;
1198   t878 = t172 * t877;
1199   t879 = t878 * t52;
1200   t882 = 0.22478670955426118383e-2 * t827 * t52 + 0.44957341910852236766e-2 * t412 * t145 + 0.22478670955426118383e-2 * t174 * t353 + 0.16555529631147769304e0 * t834 * t677 * t161 + 0.34545222222222222222e-2 * t359 * t362 - 0.11515074074074074074e-2 * t151 * t840 - 0.50586340539618183985e-1 * t371 * t497 * t161 + 0.27592549385246282174e-1 * t848 * t506 * t849 - 0.14409904107548592783e1 * t856 * t506 * t857 + 0.22993791154371901812e-2 * t861 * t497 * t862 + 0.18429583437767336288e2 * t869 * t506 * t872 + 0.93661128980942159930e-4 * t879 * t217;
1201   t883 = t790 + t882;
1202   t885 = t883 * t43 * t47;
1203   t887 = t424 * t135 * t47;
1204   t888 = 0.3e1 * t887;
1205   t890 = t177 * t334 * t47;
1206   t893 = t69 * t714 * t47;
1207   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1208     v3rho3[0] = t481 - 0.6e1 * t312 - 0.3e1 * t336 + t484 + 0.6e1 * t428 + 0.3e1 * t431 + t7 * (t683 - t686 - 0.3e1 * t688 - t716 + t885 + t888 + 0.3e1 * t890 + t893);
1209 
1210   t898 = 0.2e1 * t450;
1211   t899 = 0.2e1 * t454;
1212   t902 = t309 * t189 * t47;
1213   t904 = t125 * t448 * t47;
1214   t905 = 0.2e1 * t904;
1215   t906 = t691 * t184;
1216   t909 = t315 * t34;
1217   t919 = t704 * t186;
1218   t922 = t327 * t34;
1219   t932 = -0.8e1 / 0.27e2 * t906 * t316 + 0.16e2 / 0.9e1 * t909 * t320 * t130 + 0.4e1 / 0.9e1 * t436 * t323 + 0.8e1 / 0.3e1 * t38 * t320 - 0.8e1 * t439 * t597 - 0.8e1 / 0.27e2 * t919 * t328 - 0.16e2 / 0.9e1 * t922 * t320 * t132 + 0.4e1 / 0.9e1 * t442 * t331 - 0.8e1 / 0.3e1 * t41 * t320 + 0.8e1 * t445 * t597;
1220   t934 = t33 * t932 * t47;
1221   t937 = t424 * t189 * t47;
1222   t939 = t177 * t448 * t47;
1223   t940 = 0.2e1 * t939;
1224   t942 = t69 * t932 * t47;
1225   t943 = t683 - 0.2e1 * t685 - t688 - t902 - t905 - t934 + t885 + 0.2e1 * t887 + t890 + t937 + t940 + t942;
1226   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1227     v3rho3[1] = t7 * t943 - 0.4e1 * t312 - t336 + 0.4e1 * t428 + t431 - t459 + t476 + t481 + t484 - t898 + t899;
1228 
1229   t949 = t125 * t473 * t47;
1230   t950 = t691 * t460;
1231   t955 = t315 * t464;
1232   t960 = -0.2e1 * t320 - 0.6e1 * t698;
1233   t963 = t704 * t467;
1234   t968 = t327 * t470;
1235   t971 = -t960;
1236   t974 = -0.8e1 / 0.27e2 * t950 * t130 + 0.16e2 / 0.9e1 * t436 * t321 + 0.4e1 / 0.9e1 * t955 * t130 + 0.4e1 / 0.3e1 * t38 * t960 - 0.8e1 / 0.27e2 * t963 * t132 - 0.16e2 / 0.9e1 * t442 * t321 + 0.4e1 / 0.9e1 * t968 * t132 + 0.4e1 / 0.3e1 * t41 * t971;
1237   t976 = t33 * t974 * t47;
1238   t979 = t177 * t473 * t47;
1239   t981 = t69 * t974 * t47;
1240   t982 = t683 - t685 - 0.2e1 * t902 - t905 - t949 - t976 + t885 + t887 + 0.2e1 * t937 + t940 + t979 + t981;
1241   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1242     v3rho3[2] = t7 * t982 - t313 + t429 - 0.4e1 * t435 + 0.4e1 * t452 - t475 + t478 + t481 + t484 - t898 + t899;
1243 
1244   t988 = 0.3e1 * t902;
1245   t990 = t460 * t184;
1246   t996 = -0.6e1 * t320 - 0.6e1 * t698;
1247   t999 = t467 * t186;
1248   t1004 = -t996;
1249   t1007 = -0.8e1 / 0.27e2 * t691 * t990 + 0.4e1 / 0.3e1 * t436 * t464 + 0.4e1 / 0.3e1 * t38 * t996 - 0.8e1 / 0.27e2 * t704 * t999 + 0.4e1 / 0.3e1 * t442 * t470 + 0.4e1 / 0.3e1 * t41 * t1004;
1250   t1009 = t33 * t1007 * t47;
1251   t1010 = 0.3e1 * t937;
1252   t1013 = t69 * t1007 * t47;
1253   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1254     v3rho3[3] = t481 - 0.6e1 * t435 - 0.3e1 * t475 + t484 + 0.6e1 * t452 + 0.3e1 * t478 + t7 * (t683 - t988 - 0.3e1 * t949 - t1009 + t885 + t1010 + 0.3e1 * t979 + t1013);
1255 
1256 #ifndef XC_DONT_COMPILE_LXC
1257 
1258   if(order < 4) return;
1259 
1260 
1261   t1016 = 0.4e1 * t683;
1262   t1020 = 0.4e1 * t885;
1263   t1038 = 0.1e1 / t224 / t596;
1264   t1039 = t5 * t1038;
1265   t1043 = t1039 * t515;
1266   t1049 = t97 * t5;
1267   t1054 = 0.1e1 / t267 / t250;
1268   t1065 = 0.1e1 / t525 / t22;
1269   t1080 = 0.31090700000000000001e-1 * t612 * t247 + 0.31090700000000000001e-1 * t239 * t620 + 0.10363566666666666667e-1 * t98 * t99 * t606 + 0.10363566666666666667e-1 * t98 * t243 * t233 + 0.12112784639430467768e-3 * t537 * t90 * t217 - 0.19441031218901976532e0 * t503 * t1039 * t507 + 0.17168305804406539632e2 * t514 * t1043 + 0.17673664744456342301e-1 * t269 * t2 * t1038 * t109 * t95 * t1049 - 0.23411326096918008589e1 * t1054 * t2 * t1038 * t273 * t95 * t1049 - 0.12420094242202632684e0 * t513 * t97 * t1043 + 0.67508438032335735592e2 * t1065 * t2 * t1038 * t530 * t95 * t1049 + 0.12112784639430467768e-3 * t296 * t536 * t15 * t217 + 0.47129772651883579471e0 * t258 * t223 * t1038 * t109;
1270   t1082 = 0.1e1 / t8 / t596;
1271   t1087 = t6 * t1082;
1272   t1091 = t1039 * t531;
1273   t1099 = t596 * t7;
1274   t1100 = 0.1e1 / t1099;
1275   t1101 = t2 * t1100;
1276   t1126 = t4 * t1087 * t16;
1277   t1130 = t78 * t496 * t80 * t90;
1278   t1134 = t78 * t201 * t211 * t213;
1279   t1137 = t78 * t567 * t233;
1280   t1150 = 0.1e1 / t578 / t15;
1281   t1152 = t213 * t213;
1282   t1156 = t213 * t233;
1283   t1160 = t233 * t233;
1284   t1164 = t90 * t606;
1285   t1168 = t4 * t1087;
1286   t1169 = 0.70e2 / 0.81e2 * t1168;
1287   t1173 = 0.1e1 / t13 / t2 / t35 / 0.48e2;
1288   t1174 = t1173 * t2;
1289   t1176 = 0.1e1 / t8 / t1099;
1290   t1178 = t1174 * t1176 * t78;
1291   t1180 = t595 * t1100;
1292   t1183 = t221 * t223 * t1038;
1293   t1186 = t85 * t86 * t1082;
1294   t1188 = t1169 - 0.10354000000000000000e1 * t1178 + 0.99398400000000000000e1 * t1180 - 0.36814222222222222222e1 * t1183 + 0.32212444444444444444e1 * t1186;
1295   t1192 = 0.70e2 / 0.81e2 * t1126 + 0.28e2 / 0.27e2 * t1130 + 0.4e1 / 0.3e1 * t1134 - 0.2e1 / 0.3e1 * t1137 + 0.2e1 * t78 * t73 * t579 * t581 - 0.2e1 * t78 * t571 * t585 + t78 * t206 * t606 / 0.3e1 + 0.6e1 * t78 * t9 * t1150 * t1152 - 0.9e1 * t78 * t580 * t1156 + 0.3e1 / 0.2e1 * t78 * t212 * t1160 + 0.2e1 * t78 * t212 * t1164 - t78 * t81 * t1188 / 0.4e1;
1296   t1197 = t1087 * t273;
1297   t1200 = -0.41238551070398132037e0 * t106 * t86 * t1082 * t109 - 0.16200859349084980444e-1 * t520 * t1087 * t521 - 0.37129640917784654576e3 * t528 * t1091 - 0.14804514559303905049e-3 * t538 * t589 + 0.26860754033147841838e1 * t527 * t97 * t1091 - 0.12725038616008566457e1 * t676 * t1101 * t109 - 0.46822652193836017178e1 * t512 * t2 * t1100 * t273 * t220 + 0.10126265704850360339e3 * t526 * t2 * t1100 * t530 * t220 + 0.53020994233369026904e-1 * t501 * t2 * t1100 * t109 * t220 + 0.50245625170970829257e-3 * t304 * t589 + 0.64601518076962494762e-3 * t302 * t15 * t598 + 0.10363566666666666667e-1 * t1192 * t95 * t97 * t100 + 0.19517290952032708503e0 * t623 * t1197;
1298   t1207 = t272 * t272;
1299   t1211 = t95 * t97 * t5;
1300   t1243 = -0.77340987998746992661e1 * t271 * t1197 - 0.60995923043525943529e3 / t525 / t250 * t1038 / t1207 * t1211 - 0.10645795064745113728e0 * t1054 * t1038 * t273 * t1211 + 0.18802527823203489287e2 * t1065 * t1038 * t530 * t1211 - 0.69090444444444444446e-2 * t98 * t491 * t90 + 0.20727133333333333334e-1 * t239 * t617 + 0.10363566666666666667e-1 * t612 * t244 + 0.38383580246913580248e-2 * t98 * t286 * t15 + 0.20435174860777645786e0 * t253 * t1087 * t109 - 0.69090444444444444446e-2 * t239 * t492 + 0.29070683134633122642e-2 * t297 * t233 + 0.96902277115443742139e-3 * t122 * t606 + 0.29070683134633122642e-2 * t669 * t90;
1301   t1255 = t285 * t1039;
1302   t1267 = t1082 * t2;
1303   t1271 = -0.4e1 * t632 * t86 * t280 * t233 + 0.16e2 / 0.3e1 * t114 * t1180 + 0.185e3 / 0.864e3 * t1126 - t650 * t651 * t80 * t90 / 0.108e3 - 0.160e3 / 0.81e2 * t284 * t1255 + 0.140e3 / 0.81e2 * t115 * t1168 - 0.36e2 * t662 * t1156 + 0.8e1 * t292 * t1164 + 0.11e2 / 0.54e2 * t1130 + t1134 / 0.6e1 - t1137 / 0.12e2 - 0.11e2 / 0.648e3 * t650 * t1267 * t16;
1304   t1310 = t1176 * t1 * t86;
1305   t1322 = 0.8e1 / 0.3e1 * t278 * t594 * t677 * t90 - t118 * t1188 + 0.24e2 * t29 * t1150 * t1152 + 0.6e1 * t292 * t1160 + 0.4e1 * t28 * t579 * t85 * t86 * t73 * t581 + 0.8e1 / 0.3e1 * t632 * t86 * t201 * t213 - 0.16e2 / 0.9e1 * t636 * t223 * t505 * t90 + 0.56e2 / 0.27e2 * t279 * t86 * t496 * t90 - 0.4e1 / 0.3e1 * t279 * t86 * t201 * t233 + 0.2e1 / 0.3e1 * t279 * t86 * t73 * t606 - 0.5e1 / 0.9e1 * t114 * t1174 * t1310 - 0.4e1 / 0.3e1 * t631 * t221 * t223 * t226 * t213 + 0.2e1 / 0.3e1 * t636 * t223 * t226 * t233;
1306   t1328 = t2 * t1176;
1307   t1337 = t519 * t6 * t651;
1308   t1363 = t1038 * t2;
1309   t1373 = t223 * t505 * t84;
1310   t1378 = 0.96902277115443742139e-3 * (t1271 + t1322) * t121 * t15 + 0.13255248558342256727e0 * t104 * t1173 * t1328 * t109 * t1 * t86 + 0.67293247996835932044e-5 * t537 * t15 * t1 * t1337 + 0.96902277115443742140e-3 * t551 * t556 * t82 + 0.32300759038481247380e-3 * t551 * t560 * t287 + 0.48451138557721871070e-3 * t302 * t233 * t84 * t82 + 0.32300759038481247380e-3 * t302 * t90 * t220 * t287 - 0.64601518076962494760e-3 * t557 * t217 - 0.43067678717974996507e-3 * t561 * t657 + 0.48451138557721871070e-3 * t668 * t301 * t303 * t82 + 0.44184161861140855756e-2 * t502 * t1049 * t1363 * t507 + 0.53834598397468745633e-4 * t120 / t535 / t28 * t15 * t95 * t1373 - 0.64601518076962494760e-3 * t552 * t217;
1311   t1381 = (t1080 + t1200 + t1243 + t1378) * t49;
1312   t1383 = t682 * t135 * t47;
1313   t1386 = t309 * t334 * t47;
1314   t1389 = t125 * t714 * t47;
1315   t1391 = t37 * t37;
1316   t1393 = 0.1e1 / t314 / t1391;
1317   t1394 = t316 * t316;
1318   t1400 = t323 * t323;
1319   t1405 = t34 * t1100;
1320   t1407 = -0.24e2 * t597 + 0.24e2 * t1405;
1321   t1410 = t40 * t40;
1322   t1412 = 0.1e1 / t326 / t1410;
1323   t1413 = t328 * t328;
1324   t1419 = t331 * t331;
1325   t1427 = 0.40e2 / 0.81e2 * t1393 * t1394 - 0.16e2 / 0.9e1 * t691 * t316 * t323 + 0.4e1 / 0.3e1 * t315 * t1400 + 0.16e2 / 0.9e1 * t695 * t700 + 0.4e1 / 0.3e1 * t38 * t1407 + 0.40e2 / 0.81e2 * t1412 * t1413 - 0.16e2 / 0.9e1 * t704 * t328 * t331 + 0.4e1 / 0.3e1 * t327 * t1419 + 0.16e2 / 0.9e1 * t708 * t711 - 0.4e1 / 0.3e1 * t41 * t1407;
1326   t1457 = t1039 * t857;
1327   t1461 = t1039 * t872;
1328   t1479 = -0.13244423704918215444e1 * t834 * t1101 * t161 + 0.36859166875534672576e2 * t867 * t2 * t1100 * t871 * t220 + 0.55185098770492564348e-1 * t846 * t2 * t1100 * t161 * t220 - 0.28819808215097185566e1 * t854 * t2 * t1100 * t391 * t220 - 0.42921743488160883383e0 * t158 * t86 * t1082 * t161 + 0.49053421129326723867e0 * t376 * t223 * t1038 * t161 - 0.76446915606552171013e-1 * t855 * t97 * t1457 + 0.97771976774839907070e0 * t868 * t97 * t1461 + 0.28098338694282647979e-3 * t411 * t877 * t52 * t217 + 0.28098338694282647979e-3 * t878 * t145 * t217 - 0.34342413959678791975e-3 * t879 * t589 - 0.16862113513206061329e-1 * t861 * t1087 * t862 - 0.20234536215847273594e0 * t848 * t1039 * t849;
1329   t1487 = 0.1e1 / t385 / t368;
1330   t1495 = 0.1e1 / t866 / t59;
1331   t1529 = 0.18395032923497521449e-1 * t387 * t2 * t1038 * t161 * t95 * t1049 - 0.14409904107548592783e1 * t1487 * t2 * t1038 * t391 * t95 * t1049 + 0.24572777917023115051e2 * t1495 * t2 * t1038 * t871 * t95 * t1049 + 0.10567263012202301374e2 * t856 * t1457 - 0.13515027854362713278e3 * t869 * t1461 + 0.74928903184753727944e-3 * t417 * t145 * t220 * t287 + 0.15610188163490359988e-4 * t878 * t52 * t1 * t1337 + 0.11239335477713059192e-2 * t826 * t416 * t418 * t82 - 0.14985780636950745589e-2 * t770 * t217 - 0.14985780636950745589e-2 * t774 * t217 - 0.99905204246338303925e-3 * t778 * t657 + 0.11655607162072802125e-2 * t419 * t589 + 0.22478670955426118383e-2 * t769 * t773 * t82;
1332   t1546 = t4 * t1087 * t53;
1333   t1550 = t78 * t496 * t142 * t145;
1334   t1554 = t78 * t201 * t345 * t347;
1335   t1557 = t78 * t720 * t353;
1336   t1570 = 0.1e1 / t731 / t52;
1337   t1572 = t347 * t347;
1338   t1576 = t347 * t353;
1339   t1580 = t353 * t353;
1340   t1584 = t145 * t745;
1341   t1592 = t1169 - 0.19612277777777777777e1 * t1178 + 0.18827786666666666666e2 * t1180 - 0.69732543209876543208e1 * t1183 + 0.61015975308641975307e1 * t1186;
1342   t1596 = 0.70e2 / 0.81e2 * t1546 + 0.28e2 / 0.27e2 * t1550 + 0.4e1 / 0.3e1 * t1554 - 0.2e1 / 0.3e1 * t1557 + 0.2e1 * t78 * t73 * t732 * t734 - 0.2e1 * t78 * t724 * t738 + t78 * t340 * t745 / 0.3e1 + 0.6e1 * t78 * t9 * t1570 * t1572 - 0.9e1 * t78 * t733 * t1576 + 0.3e1 / 0.2e1 * t78 * t346 * t1580 + 0.2e1 * t78 * t346 * t1584 - t78 * t143 * t1592 / 0.4e1;
1343   t1621 = t390 * t390;
1344   t1637 = 0.74928903184753727944e-3 * t769 * t777 * t287 + 0.11239335477713059192e-2 * t417 * t353 * t84 * t82 + 0.51817833333333333333e-2 * t751 * t362 - 0.34545222222222222222e-2 * t151 * t491 * t145 + 0.10363566666666666667e-1 * t359 * t756 + 0.51817833333333333333e-2 * t1596 * t95 * t97 * t152 + 0.12488150530792287991e-3 * t172 / t876 / t64 * t52 * t95 * t1373 + 0.45987582308743803624e-2 * t847 * t1049 * t1363 * t849 + 0.13796274692623141087e0 * t156 * t1173 * t1328 * t161 * t1 * t86 - 0.13129800378398182140e3 / t866 / t368 * t1038 / t1621 * t1211 - 0.65525927662759003726e-1 * t1487 * t1038 * t391 * t1211 + 0.68440383742387934953e1 * t1495 * t1038 * t871 * t1211 + 0.19191790123456790123e-2 * t151 * t286 * t52;
1345   t1646 = t1087 * t391;
1346   t1704 = 0.185e3 / 0.864e3 * t1546 + 0.16e2 / 0.3e1 * t166 * t1180 + 0.8e1 * t407 * t1584 - 0.36e2 * t820 * t1576 - 0.4e1 * t794 * t86 * t398 * t353 - 0.160e3 / 0.81e2 * t402 * t1255 + 0.140e3 / 0.81e2 * t167 * t1168 - t650 * t651 * t142 * t145 / 0.108e3 - 0.16e2 / 0.9e1 * t798 * t223 * t505 * t145 + 0.56e2 / 0.27e2 * t397 * t86 * t496 * t145 + 0.4e1 * t64 * t732 * t85 * t86 * t73 * t734 + 0.8e1 / 0.3e1 * t794 * t86 * t201 * t347;
1347   t1741 = -0.4e1 / 0.3e1 * t793 * t221 * t223 * t226 * t347 + 0.2e1 / 0.3e1 * t798 * t223 * t226 * t353 - 0.4e1 / 0.3e1 * t397 * t86 * t201 * t353 + 0.2e1 / 0.3e1 * t397 * t86 * t73 * t745 - 0.5e1 / 0.9e1 * t166 * t1174 * t1310 - 0.11e2 / 0.648e3 * t650 * t1267 * t53 + 0.8e1 / 0.3e1 * t396 * t594 * t677 * t145 - t1557 / 0.12e2 + 0.11e2 / 0.54e2 * t1550 + t1554 / 0.6e1 + 0.24e2 * t65 * t1570 * t1572 + 0.6e1 * t407 * t1580 - t170 * t1592;
1348   t1746 = -0.34545222222222222222e-2 * t359 * t840 + 0.51817833333333333333e-2 * t151 * t99 * t745 + 0.14985780636950745589e-2 * t417 * t52 * t598 + 0.12013086738172484016e0 * t784 * t1646 + 0.21269256817794009175e0 * t371 * t1087 * t161 - 0.47604147498151601160e1 * t389 * t1646 + 0.51817833333333333333e-2 * t151 * t243 * t353 + 0.15545350000000000000e-1 * t359 * t781 + 0.15545350000000000000e-1 * t751 * t365 + 0.22478670955426118383e-2 * t174 * t745 + 0.67436012866278355149e-2 * t827 * t145 + 0.67436012866278355149e-2 * t412 * t353 + 0.22478670955426118383e-2 * (t1704 + t1741) * t173 * t52;
1349   t1750 = (t1479 + t1529 + t1637 + t1746) * t43 * t47;
1350   t1752 = t883 * t135 * t47;
1351   t1755 = t424 * t334 * t47;
1352   t1758 = t177 * t714 * t47;
1353   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1354     v4rho4[0] = t1016 - 0.12e2 * t685 - 0.12e2 * t688 - 0.4e1 * t716 + t1020 + 0.12e2 * t887 + 0.12e2 * t890 + 0.4e1 * t893 + t7 * (-t33 * t1427 * t47 + t69 * t1427 * t47 + t1381 - 0.4e1 * t1383 - 0.6e1 * t1386 - 0.4e1 * t1389 + t1750 + 0.4e1 * t1752 + 0.6e1 * t1755 + 0.4e1 * t1758);
1355 
1356   t1767 = t682 * t189 * t47;
1357   t1769 = t309 * t448 * t47;
1358   t1770 = 0.3e1 * t1769;
1359   t1772 = t125 * t932 * t47;
1360   t1803 = 0.32e2 * t439 * t1100;
1361   t1821 = 0.32e2 * t445 * t1100;
1362   t1824 = -0.16e2 * t38 * t597 + 0.4e1 / 0.9e1 * t442 * t711 + 0.16e2 * t41 * t597 - 0.8e1 * t909 * t597 * t130 + 0.8e1 * t922 * t597 * t132 + 0.40e2 / 0.81e2 * t1393 * t184 * t692 - 0.8e1 / 0.9e1 * t906 * t130 * t323 - 0.16e2 / 0.9e1 * t691 * t34 * t320 * t316 + 0.8e1 / 0.3e1 * t315 * t320 * t130 + 0.8e1 / 0.3e1 * t909 * t320 * t323 + t1803 + 0.40e2 / 0.81e2 * t1412 * t186 * t705 - 0.8e1 / 0.9e1 * t919 * t132 * t331 + 0.16e2 / 0.9e1 * t704 * t34 * t320 * t328 - 0.8e1 / 0.3e1 * t327 * t320 * t132 - 0.8e1 / 0.3e1 * t922 * t320 * t331 - t1821 + 0.4e1 / 0.9e1 * t436 * t700;
1363   t1830 = t883 * t189 * t47;
1364   t1832 = t424 * t448 * t47;
1365   t1833 = 0.3e1 * t1832;
1366   t1835 = t177 * t932 * t47;
1367   t1839 = -t33 * t1824 * t47 + t69 * t1824 * t47 + t1381 - 0.3e1 * t1383 - 0.3e1 * t1386 - t1389 + t1750 + 0.3e1 * t1752 + 0.3e1 * t1755 + t1758 - t1767 - t1770 - 0.3e1 * t1772 + t1830 + t1833 + 0.3e1 * t1835;
1368   t1841 = 0.6e1 * t904;
1369   t1843 = 0.6e1 * t939;
1370   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1371     v4rho4[1] = t7 * t1839 + t1010 + t1016 + t1020 - t1841 + t1843 - 0.9e1 * t685 - 0.6e1 * t688 - t716 + 0.9e1 * t887 + 0.6e1 * t890 + t893 - 0.3e1 * t934 + 0.3e1 * t942 - t988;
1372 
1373   t1854 = t309 * t473 * t47;
1374   t1856 = t125 * t974 * t47;
1375   t1867 = t34 * t34;
1376   t1870 = 0.1e1 / t596 / t127;
1377   t1909 = 0.40e2 / 0.81e2 * t1393 * t460 * t316 - 0.64e2 / 0.27e2 * t906 * t130 * t34 * t320 - 0.8e1 / 0.27e2 * t950 * t323 + 0.32e2 / 0.9e1 * t315 * t1867 * t1870 + 0.16e2 / 0.9e1 * t436 * t320 - 0.16e2 / 0.3e1 * t436 * t698 - 0.8e1 / 0.27e2 * t691 * t464 * t316 + 0.8e1 / 0.9e1 * t315 * t960 * t130 + 0.4e1 / 0.9e1 * t955 * t323 + t1803 + 0.40e2 / 0.81e2 * t1412 * t467 * t328 + 0.64e2 / 0.27e2 * t919 * t132 * t34 * t320 - 0.8e1 / 0.27e2 * t963 * t331 + 0.32e2 / 0.9e1 * t327 * t1867 * t1870 - 0.16e2 / 0.9e1 * t442 * t320 + 0.16e2 / 0.3e1 * t442 * t698 - 0.8e1 / 0.27e2 * t704 * t470 * t328 + 0.8e1 / 0.9e1 * t327 * t971 * t132 + 0.4e1 / 0.9e1 * t968 * t331 - t1821;
1378   t1917 = t424 * t473 * t47;
1379   t1919 = t177 * t974 * t47;
1380   t1923 = -t33 * t1909 * t47 + t69 * t1909 * t47 + t1381 - 0.2e1 * t1383 - t1386 + t1750 + 0.2e1 * t1752 + t1755 - 0.2e1 * t1767 - 0.4e1 * t1769 - 0.2e1 * t1772 + 0.2e1 * t1830 + 0.4e1 * t1832 + 0.2e1 * t1835 - t1854 - 0.2e1 * t1856 + t1917 + 0.2e1 * t1919;
1381   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1382     v4rho4[2] = t7 * t1923 + t1016 + t1020 - 0.6e1 * t685 - 0.2e1 * t688 + 0.6e1 * t887 + 0.2e1 * t890 - 0.6e1 * t902 - 0.8e1 * t904 - 0.2e1 * t934 + 0.6e1 * t937 + 0.8e1 * t939 + 0.2e1 * t942 - 0.2e1 * t949 - 0.2e1 * t976 + 0.2e1 * t979 + 0.2e1 * t981;
1383 
1384   t1949 = t125 * t1007 * t47;
1385   t1968 = 0.12e2 * t597 + 0.24e2 * t1405;
1386   t1990 = 0.40e2 / 0.81e2 * t1393 * t990 * t130 - 0.16e2 / 0.9e1 * t950 * t321 - 0.8e1 / 0.9e1 * t906 * t464 * t130 + 0.8e1 / 0.3e1 * t909 * t320 * t464 + 0.4e1 / 0.3e1 * t436 * t960 + 0.4e1 / 0.9e1 * t315 * t996 * t130 + 0.4e1 / 0.3e1 * t38 * t1968 + 0.40e2 / 0.81e2 * t1412 * t999 * t132 + 0.16e2 / 0.9e1 * t963 * t321 - 0.8e1 / 0.9e1 * t919 * t470 * t132 - 0.8e1 / 0.3e1 * t922 * t320 * t470 + 0.4e1 / 0.3e1 * t442 * t971 + 0.4e1 / 0.9e1 * t327 * t1004 * t132 - 0.4e1 / 0.3e1 * t41 * t1968;
1387   t1997 = t177 * t1007 * t47;
1388   t2000 = -t33 * t1990 * t47 + t69 * t1990 * t47 + t1381 - t1383 + t1750 + t1752 - 0.3e1 * t1767 - t1770 + 0.3e1 * t1830 + t1833 - 0.3e1 * t1854 - 0.3e1 * t1856 + 0.3e1 * t1917 + 0.3e1 * t1919 - t1949 + t1997;
1389   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1390     v4rho4[3] = t7 * t2000 - t1009 + t1013 + t1016 + t1020 - t1841 + t1843 - t686 + t888 - 0.9e1 * t902 + 0.9e1 * t937 - 0.6e1 * t949 - 0.3e1 * t976 + 0.6e1 * t979 + 0.3e1 * t981;
1391 
1392   t2011 = t460 * t460;
1393   t2016 = t464 * t464;
1394   t2022 = 0.24e2 * t597 + 0.24e2 * t1405;
1395   t2025 = t467 * t467;
1396   t2030 = t470 * t470;
1397   t2038 = 0.40e2 / 0.81e2 * t1393 * t2011 - 0.16e2 / 0.9e1 * t950 * t464 + 0.4e1 / 0.3e1 * t315 * t2016 + 0.16e2 / 0.9e1 * t436 * t996 + 0.4e1 / 0.3e1 * t38 * t2022 + 0.40e2 / 0.81e2 * t1412 * t2025 - 0.16e2 / 0.9e1 * t963 * t470 + 0.4e1 / 0.3e1 * t327 * t2030 + 0.16e2 / 0.9e1 * t442 * t1004 - 0.4e1 / 0.3e1 * t41 * t2022;
1398   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1399     v4rho4[4] = t1016 - 0.12e2 * t902 - 0.12e2 * t949 - 0.4e1 * t1009 + t1020 + 0.12e2 * t937 + 0.12e2 * t979 + 0.4e1 * t1013 + t7 * (-t33 * t2038 * t47 + t69 * t2038 * t47 + t1381 + t1750 - 0.4e1 * t1767 + 0.4e1 * t1830 - 0.6e1 * t1854 + 0.6e1 * t1917 - 0.4e1 * t1949 + 0.4e1 * t1997);
1400 
1401 #ifndef XC_DONT_COMPILE_MXC
1402 
1403   if(order < 5) return;
1404 
1405 
1406 #endif
1407 
1408 #endif
1409 
1410 #endif
1411 
1412 #endif
1413 
1414 #endif
1415 
1416 
1417 }
1418 
1419