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