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_rc04.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, t3, t4, t5, t6, t7, t12, t13;
24 double t16, t17, t19, t21;
25
26 #ifndef XC_DONT_COMPILE_VXC
27 double t23, t24, t25;
28
29 #ifndef XC_DONT_COMPILE_FXC
30 double t30, t35, t36;
31
32 #ifndef XC_DONT_COMPILE_KXC
33 double t45, t53, t62, t65, t68, t69;
34
35 #ifndef XC_DONT_COMPILE_LXC
36 double t75, t86, t90, t99, t103, t104;
37 #endif
38
39 #endif
40
41 #endif
42
43 #endif
44
45 #endif
46
47
48
49 t1 = M_CBRT3;
50 t3 = POW_1_3(0.1e1 / M_PI);
51 t4 = t1 * t3;
52 t5 = M_CBRT4;
53 t6 = t5 * t5;
54 t7 = POW_1_3(rho[0]);
55 t12 = 0.4888270e1 + 0.79425925000000000000e0 * t4 * t6 / t7;
56 t13 = atan(t12);
57 t16 = t1 * t1;
58 t17 = (-0.655868e0 * t13 + 0.897889e0) * t16;
59 t19 = 0.1e1 / t3 * t5;
60 t21 = t17 * t19 * t7;
61 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
62 zk[0] = t21 / 0.3e1;
63
64 #ifndef XC_DONT_COMPILE_VXC
65
66 if(order < 1) return;
67
68
69 t23 = t12 * t12;
70 t24 = t23 + 0.1e1;
71 t25 = 0.1e1 / t24;
72 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
73 vrho[0] = 0.4e1 / 0.9e1 * t21 + 0.69457230103866666663e0 * t25;
74
75 #ifndef XC_DONT_COMPILE_FXC
76
77 if(order < 2) return;
78
79
80 t30 = t7 * t7;
81 t35 = t24 * t24;
82 t36 = 0.1e1 / t35;
83 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
84 v2rho2[0] = 0.92609640138488888884e0 / rho[0] * t25 + 0.4e1 / 0.27e2 * t17 * t19 / t30 + 0.36778031659583040509e0 * t36 * t12 * t1 * t3 * t6 / t7 / rho[0];
85
86 #ifndef XC_DONT_COMPILE_KXC
87
88 if(order < 3) return;
89
90
91 t45 = rho[0] * rho[0];
92 t53 = t4 * t6;
93 t62 = 0.1e1 / t35 / t24;
94 t65 = t3 * t3;
95 t68 = 0.1e1 / t30 / t45;
96 t69 = t65 * t5 * t68;
97 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
98 v3rho3[0] = -0.61739760092325925923e0 / t45 * t25 - 0.1e-19 / t7 / t45 * t36 * t12 * t53 - 0.8e1 / 0.81e2 * t17 * t19 / t30 / rho[0] + 0.15579355649288896569e1 * t62 * t23 * t16 * t69 - 0.38948389123222241422e0 * t36 * t16 * t69;
99
100 #ifndef XC_DONT_COMPILE_LXC
101
102 if(order < 4) return;
103
104
105 t75 = t45 * rho[0];
106 t86 = 0.1e1 / t30 / t75;
107 t90 = t16 * t65 * t5;
108 t99 = t35 * t35;
109 t103 = t45 * t45;
110 t104 = 0.1e1 / t103;
111 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
112 v4rho4[0] = 0.10289960015387654321e1 / t75 * t25 - 0.32691583697407147117e0 / t7 / t75 * t36 * t12 * t53 - 0.41544948398103724184e1 * t86 * t62 * t23 * t90 + 0.10386237099525931046e1 * t86 * t36 * t90 + 0.40e2 / 0.243e3 * t17 * t19 * t68 + 0.94530758360525579704e1 / t99 * t23 * t12 * t104 - 0.47265379180262789851e1 * t62 * t12 * t104;
113
114 #ifndef XC_DONT_COMPILE_MXC
115
116 if(order < 5) return;
117
118
119 #endif
120
121 #endif
122
123 #endif
124
125 #endif
126
127 #endif
128
129
130 }
131
132
133 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))134 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
135 {
136
137 #ifndef XC_DONT_COMPILE_EXC
138 double t1, t3, t4, t5, t6, t7, t12, t13;
139 double t16, t17, t19, t21;
140
141 #ifndef XC_DONT_COMPILE_VXC
142 double t23, t24, t25;
143
144 #ifndef XC_DONT_COMPILE_FXC
145 double t30, t35, t36;
146
147 #ifndef XC_DONT_COMPILE_KXC
148 double t45, t53, t62, t65, t68, t69;
149
150 #ifndef XC_DONT_COMPILE_LXC
151 double t75, t86, t90, t99, t103, t104;
152 #endif
153
154 #endif
155
156 #endif
157
158 #endif
159
160 #endif
161
162
163
164 t1 = M_CBRT3;
165 t3 = POW_1_3(0.1e1 / M_PI);
166 t4 = t1 * t3;
167 t5 = M_CBRT4;
168 t6 = t5 * t5;
169 t7 = POW_1_3(rho[0]);
170 t12 = 0.4888270e1 + 0.79425925000000000000e0 * t4 * t6 / t7;
171 t13 = atan(t12);
172 t16 = t1 * t1;
173 t17 = (-0.655868e0 * t13 + 0.897889e0) * t16;
174 t19 = 0.1e1 / t3 * t5;
175 t21 = t17 * t19 * t7;
176 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
177 zk[0] = t21 / 0.6e1;
178
179 #ifndef XC_DONT_COMPILE_VXC
180
181 if(order < 1) return;
182
183
184 t23 = t12 * t12;
185 t24 = t23 + 0.1e1;
186 t25 = 0.1e1 / t24;
187 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
188 vrho[0] = 0.2e1 / 0.9e1 * t21 + 0.34728615051933333332e0 * t25;
189
190 #ifndef XC_DONT_COMPILE_FXC
191
192 if(order < 2) return;
193
194
195 t30 = t7 * t7;
196 t35 = t24 * t24;
197 t36 = 0.1e1 / t35;
198 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
199 v2rho2[0] = 0.46304820069244444442e0 / rho[0] * t25 + 0.2e1 / 0.27e2 * t17 * t19 / t30 + 0.18389015829791520255e0 * t36 * t12 * t1 * t3 * t6 / t7 / rho[0];
200
201 #ifndef XC_DONT_COMPILE_KXC
202
203 if(order < 3) return;
204
205
206 t45 = rho[0] * rho[0];
207 t53 = t4 * t6;
208 t62 = 0.1e1 / t35 / t24;
209 t65 = t3 * t3;
210 t68 = 0.1e1 / t30 / t45;
211 t69 = t65 * t5 * t68;
212 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
213 v3rho3[0] = -0.30869880046162962961e0 / t45 * t25 - 0.1e-19 / t7 / t45 * t36 * t12 * t53 - 0.4e1 / 0.81e2 * t17 * t19 / t30 / rho[0] + 0.77896778246444482846e0 * t62 * t23 * t16 * t69 - 0.19474194561611120712e0 * t36 * t16 * t69;
214
215 #ifndef XC_DONT_COMPILE_LXC
216
217 if(order < 4) return;
218
219
220 t75 = t45 * rho[0];
221 t86 = 0.1e1 / t30 / t75;
222 t90 = t16 * t65 * t5;
223 t99 = t35 * t35;
224 t103 = t45 * t45;
225 t104 = 0.1e1 / t103;
226 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
227 v4rho4[0] = 0.51449800076938271602e0 / t75 * t25 - 0.16345791848703573557e0 / t7 / t75 * t36 * t12 * t53 - 0.20772474199051862092e1 * t86 * t62 * t23 * t90 + 0.51931185497629655233e0 * t86 * t36 * t90 + 0.20e2 / 0.243e3 * t17 * t19 * t68 + 0.47265379180262789853e1 / t99 * t23 * t12 * t104 - 0.23632689590131394927e1 * t62 * t12 * t104;
228
229 #ifndef XC_DONT_COMPILE_MXC
230
231 if(order < 5) return;
232
233
234 #endif
235
236 #endif
237
238 #endif
239
240 #endif
241
242 #endif
243
244
245 }
246
247
248 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))249 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
250 {
251
252 #ifndef XC_DONT_COMPILE_EXC
253 double t1, t2, t3, t4, t5, t6, t7, t8;
254 double t9, t10, t12, t13, t14, t15, t17, t18;
255 double t19, t20, t21, t26, t27, t29, t31, t32;
256 double t33, t34, t35, t36;
257
258 #ifndef XC_DONT_COMPILE_VXC
259 double t37, t38, t40, t41, t42, t43, t44, t45;
260 double t46, t48, t49, t52, t56, t57, t58, t60;
261 double t61, t63, t66, t68;
262
263 #ifndef XC_DONT_COMPILE_FXC
264 double t71, t72, t74, t78, t79, t80, t83, t84;
265 double t85, t86, t91, t92, t95, t96, t99, t100;
266 double t101, t103, t107, t108, t111, t114, t116, t118;
267 double t119, t120, t123, t124, t126, t129, t130, t132;
268 double t133, t134, t135, t138, t140, t143, t146, t149;
269 double t152, t154, t158, t160, t164, t168, t171, t174;
270 double t177, t179;
271
272 #ifndef XC_DONT_COMPILE_KXC
273 double t181, t182, t183, t185, t189, t191, t195, t197;
274 double t198, t200, t201, t202, t206, t207, t208, t210;
275 double t214, t217, t222, t224, t226, t229, t230, t233;
276 double t236, t238, t239, t242, t245, t246, t247, t249;
277 double t252, t254, t255, t258, t261, t264, t269, t270;
278 double t272, t273, t275, t277, t279, t281, t284, t289;
279 double t290, t292, t293, t295, t297, t300, t301, t304;
280 double t305, t308, t309, t310, t313, t314, t315, t318;
281 double t322, t325, t326, t329, t339, t342, t352, t354;
282 double t356, t362, t364, t365, t368, t370, t371, t375;
283 double t377, t379, t381, t385, t390, t395, t398, t403;
284 double t406, t409, t411, t413, t421, t426, t432, t435;
285 double t440, t443, t445, t447, t448;
286
287 #ifndef XC_DONT_COMPILE_LXC
288 double t451, t453, t455, t458, t461, t462, t463, t467;
289 double t475, t477, t480, t484, t487, t494, t497, t501;
290 double t506, t509, t510, t516, t522, t523, t525, t530;
291 double t531, t537, t549, t552, t556, t560, t567, t570;
292 double t574, t578, t579, t580, t581, t583, t584, t586;
293 double t589, t591, t597, t600, t601, t602, t604, t608;
294 double t610, t612, t613, t614, t616, t618, t633, t636;
295 double t637, t640, t643, t648, t649, t653, t657, t658;
296 double t662, t664, t683, t701, t716, t722, t725, t730;
297 double t732, t738, t740, t742, t745, t746, t747, t749;
298 double t751, t757, t761, t763, t771, t783, t796, t803;
299 double t805, t820, t823, t862, t867, t869, t874, t879;
300 double t888, t890, t893, t895, t897, t899, t901, t911;
301 double t924, t930, t939, t960, t982, t986, t989, t992;
302 double t1000, t1002, t1019, t1021, t1026, t1032, t1035, t1040;
303 double t1054, t1064;
304 #endif
305
306 #endif
307
308 #endif
309
310 #endif
311
312 #endif
313
314
315
316 t1 = rho[0] - rho[1];
317 t2 = rho[0] + rho[1];
318 t3 = 0.1e1 / t2;
319 t4 = t1 * t3;
320 t5 = 0.1e1 + t4;
321 t6 = POW_1_3(t5);
322 t7 = t6 * t6;
323 t8 = 0.1e1 - t4;
324 t9 = POW_1_3(t8);
325 t10 = t9 * t9;
326 t12 = t7 / 0.2e1 + t10 / 0.2e1;
327 t13 = t12 * t12;
328 t14 = t13 * t12;
329 t15 = M_CBRT3;
330 t17 = POW_1_3(0.1e1 / M_PI);
331 t18 = t15 * t17;
332 t19 = M_CBRT4;
333 t20 = t19 * t19;
334 t21 = POW_1_3(t2);
335 t26 = 0.4888270e1 + 0.79425925000000000000e0 * t18 * t20 / t21;
336 t27 = atan(t26);
337 t29 = -0.655868e0 * t27 + 0.897889e0;
338 t31 = t15 * t15;
339 t32 = t14 * t29 * t31;
340 t33 = 0.1e1 / t17;
341 t34 = t33 * t19;
342 t35 = t34 * t21;
343 t36 = t32 * t35;
344 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
345 zk[0] = t36 / 0.3e1;
346
347 #ifndef XC_DONT_COMPILE_VXC
348
349 if(order < 1) return;
350
351
352 t37 = 0.4e1 / 0.9e1 * t36;
353 t38 = t21 * t2;
354 t40 = t38 * t13 * t29;
355 t41 = t31 * t33;
356 t42 = 0.1e1 / t6;
357 t43 = t2 * t2;
358 t44 = 0.1e1 / t43;
359 t45 = t1 * t44;
360 t46 = t3 - t45;
361 t48 = 0.1e1 / t9;
362 t49 = -t46;
363 t52 = t42 * t46 / 0.3e1 + t48 * t49 / 0.3e1;
364 t56 = t26 * t26;
365 t57 = t56 + 0.1e1;
366 t58 = 0.1e1 / t57;
367 t60 = 0.69457230103866666663e0 * t14 * t58;
368 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
369 vrho[0] = t40 * t41 * t19 * t52 + t37 + t60;
370
371 t61 = -t3 - t45;
372 t63 = -t61;
373 t66 = t42 * t61 / 0.3e1 + t48 * t63 / 0.3e1;
374 t68 = t41 * t19 * t66;
375 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
376 vrho[1] = t40 * t68 + t37 + t60;
377
378 #ifndef XC_DONT_COMPILE_FXC
379
380 if(order < 2) return;
381
382
383 t71 = t13 * t29 * t31;
384 t72 = t21 * t52;
385 t74 = t71 * t34 * t72;
386 t78 = 0.92609640138488888884e0 * t14 * t3 * t58;
387 t79 = t21 * t21;
388 t80 = 0.1e1 / t79;
389 t83 = 0.4e1 / 0.27e2 * t32 * t34 * t80;
390 t84 = t38 * t12;
391 t85 = t84 * t29;
392 t86 = t52 * t52;
393 t91 = t13 * t58;
394 t92 = t91 * t52;
395 t95 = 0.1e1 / t6 / t5;
396 t96 = t46 * t46;
397 t99 = t43 * t2;
398 t100 = 0.1e1 / t99;
399 t101 = t1 * t100;
400 t103 = -0.2e1 * t44 + 0.2e1 * t101;
401 t107 = 0.1e1 / t9 / t8;
402 t108 = t49 * t49;
403 t111 = -t103;
404 t114 = -t95 * t96 / 0.9e1 + t42 * t103 / 0.3e1 - t107 * t108 / 0.9e1 + t48 * t111 / 0.3e1;
405 t116 = t41 * t19 * t114;
406 t118 = t57 * t57;
407 t119 = 0.1e1 / t118;
408 t120 = t14 * t119;
409 t123 = t20 / t38;
410 t124 = t18 * t123;
411 t126 = 0.36778031659583040509e0 * t120 * t26 * t124;
412 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
413 v2rho2[0] = 0.8e1 / 0.3e1 * t74 + t78 + t83 + 0.2e1 * t85 * t41 * t19 * t86 + 0.41674338062319999998e1 * t92 + t40 * t116 + t126;
414
415 t129 = t21 * t13 * t29;
416 t130 = t129 * t68;
417 t132 = t29 * t31;
418 t133 = t84 * t132;
419 t134 = t66 * t52;
420 t135 = t34 * t134;
421 t138 = t91 * t66;
422 t140 = t95 * t61;
423 t143 = t42 * t1;
424 t146 = t107 * t63;
425 t149 = t48 * t1;
426 t152 = -t140 * t46 / 0.9e1 + 0.2e1 / 0.3e1 * t143 * t100 - t146 * t49 / 0.9e1 - 0.2e1 / 0.3e1 * t149 * t100;
427 t154 = t41 * t19 * t152;
428 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
429 v2rho2[1] = 0.4e1 / 0.3e1 * t74 + t78 + t83 + 0.4e1 / 0.3e1 * t130 + 0.2e1 * t133 * t135 + 0.20837169031159999999e1 * t138 + t40 * t154 + 0.20837169031159999999e1 * t92 + t126;
430
431 t158 = t66 * t66;
432 t160 = t41 * t19 * t158;
433 t164 = t61 * t61;
434 t168 = 0.2e1 * t44 + 0.2e1 * t101;
435 t171 = t63 * t63;
436 t174 = -t168;
437 t177 = -t95 * t164 / 0.9e1 + t42 * t168 / 0.3e1 - t107 * t171 / 0.9e1 + t48 * t174 / 0.3e1;
438 t179 = t41 * t19 * t177;
439 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
440 v2rho2[2] = 0.8e1 / 0.3e1 * t130 + t78 + t83 + 0.2e1 * t85 * t160 + 0.41674338062319999998e1 * t138 + t40 * t179 + t126;
441
442 #ifndef XC_DONT_COMPILE_KXC
443
444 if(order < 3) return;
445
446
447 t181 = t12 * t29;
448 t182 = t181 * t31;
449 t183 = t21 * t86;
450 t185 = t182 * t34 * t183;
451 t189 = t71 * t34 * t80 * t52;
452 t191 = t52 * t114;
453 t195 = t13 * t119;
454 t197 = t195 * t52 * t26;
455 t198 = t197 * t124;
456 t200 = t13 * t3;
457 t201 = t58 * t52;
458 t202 = t200 * t201;
459 t206 = 0.61739760092325925923e0 * t14 * t44 * t58;
460 t207 = t12 * t58;
461 t208 = t207 * t86;
462 t210 = t91 * t114;
463 t214 = t71 * t34 * t21 * t114;
464 t217 = 0.1e1 / t21 / t43;
465 t222 = t26 * t15 * t17 * t20;
466 t224 = 0.1e-19 * t14 * t217 * t119 * t222;
467 t226 = 0.1e1 / t79 / t2;
468 t229 = 0.8e1 / 0.81e2 * t32 * t34 * t226;
469 t230 = t86 * t52;
470 t233 = t41 * t19;
471 t236 = t5 * t5;
472 t238 = 0.1e1 / t6 / t236;
473 t239 = t96 * t46;
474 t242 = t95 * t46;
475 t245 = t43 * t43;
476 t246 = 0.1e1 / t245;
477 t247 = t1 * t246;
478 t249 = 0.6e1 * t100 - 0.6e1 * t247;
479 t252 = t8 * t8;
480 t254 = 0.1e1 / t9 / t252;
481 t255 = t108 * t49;
482 t258 = t107 * t49;
483 t261 = -t249;
484 t264 = 0.4e1 / 0.27e2 * t238 * t239 - t242 * t103 / 0.3e1 + t42 * t249 / 0.3e1 + 0.4e1 / 0.27e2 * t254 * t255 - t258 * t111 / 0.3e1 + t48 * t261 / 0.3e1;
485 t269 = 0.1e1 / t118 / t57;
486 t270 = t14 * t269;
487 t272 = t17 * t17;
488 t273 = t31 * t272;
489 t275 = 0.1e1 / t79 / t43;
490 t277 = t273 * t19 * t275;
491 t279 = 0.15579355649288896569e1 * t270 * t56 * t277;
492 t281 = t272 * t19;
493 t284 = 0.38948389123222241422e0 * t120 * t31 * t281 * t275;
494 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
495 v3rho3[0] = 0.8e1 * t185 + 0.4e1 / 0.3e1 * t189 + 0.6e1 * t133 * t34 * t191 + 0.33100228493624736458e1 * t198 + 0.83348676124639999996e1 * t202 - t206 + 0.12502301418695999999e2 * t208 + 0.62511507093479999997e1 * t210 + 0.4e1 * t214 - t224 - t229 + 0.2e1 * t38 * t230 * t29 * t233 + t40 * t41 * t19 * t264 + t279 - t284;
496
497 t289 = t80 * t13 * t29;
498 t290 = t289 * t68;
499 t292 = t21 * t12;
500 t293 = t292 * t132;
501 t295 = 0.16e2 / 0.3e1 * t293 * t135;
502 t297 = t38 * t86 * t29;
503 t300 = t152 * t52;
504 t301 = t34 * t300;
505 t304 = t66 * t114;
506 t305 = t34 * t304;
507 t308 = t66 * t26;
508 t309 = t195 * t308;
509 t310 = t309 * t124;
510 t313 = 0.8e1 / 0.3e1 * t185 + 0.8e1 / 0.9e1 * t189 + 0.22066818995749824306e1 * t198 + 0.4e1 / 0.9e1 * t290 + t295 + 0.2e1 * t297 * t68 + 0.4e1 * t133 * t301 + 0.2e1 * t133 * t305 + 0.11033409497874912153e1 * t310 + 0.55565784083093333330e1 * t202 - t206;
511 t314 = t58 * t66;
512 t315 = t200 * t314;
513 t318 = 0.83348676124639999996e1 * t207 * t134;
514 t322 = 0.41674338062319999998e1 * t91 * t152;
515 t325 = 0.8e1 / 0.3e1 * t129 * t154;
516 t326 = t238 * t61;
517 t329 = t95 * t1;
518 t339 = t254 * t63;
519 t342 = t107 * t1;
520 t352 = 0.4e1 / 0.27e2 * t326 * t96 - 0.4e1 / 0.9e1 * t329 * t100 * t46 - t140 * t103 / 0.9e1 + 0.2e1 / 0.3e1 * t42 * t100 - 0.2e1 * t143 * t246 + 0.4e1 / 0.27e2 * t339 * t108 + 0.4e1 / 0.9e1 * t342 * t100 * t49 - t146 * t111 / 0.9e1 - 0.2e1 / 0.3e1 * t48 * t100 + 0.2e1 * t149 * t246;
521 t354 = t41 * t19 * t352;
522 t356 = 0.27782892041546666665e1 * t315 + t318 + 0.41674338062319999998e1 * t208 + 0.20837169031159999999e1 * t210 + t322 + 0.4e1 / 0.3e1 * t214 - t224 - t229 + t279 + t325 + t40 * t354 - t284;
523 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
524 v3rho3[1] = t313 + t356;
525
526 t362 = t292 * t29 * t160;
527 t364 = t66 * t152;
528 t365 = t34 * t364;
529 t368 = t129 * t179;
530 t370 = t177 * t52;
531 t371 = t34 * t370;
532 t375 = 0.4e1 / 0.9e1 * t189 + 0.11033409497874912153e1 * t198 + 0.8e1 / 0.9e1 * t290 + t295 + 0.22066818995749824305e1 * t310 + 0.8e1 / 0.3e1 * t362 + 0.4e1 * t133 * t365 + 0.4e1 / 0.3e1 * t368 + 0.2e1 * t133 * t371 + 0.27782892041546666665e1 * t202 - t206;
533 t377 = t207 * t158;
534 t379 = t91 * t177;
535 t381 = t38 * t52;
536 t385 = t238 * t164;
537 t390 = t95 * t168;
538 t395 = -0.2e1 * t100 - 0.6e1 * t247;
539 t398 = t254 * t171;
540 t403 = t107 * t174;
541 t406 = -t395;
542 t409 = 0.4e1 / 0.27e2 * t385 * t46 - 0.4e1 / 0.9e1 * t140 * t101 - t390 * t46 / 0.9e1 + t42 * t395 / 0.3e1 + 0.4e1 / 0.27e2 * t398 * t49 + 0.4e1 / 0.9e1 * t146 * t101 - t403 * t49 / 0.9e1 + t48 * t406 / 0.3e1;
543 t411 = t41 * t19 * t409;
544 t413 = 0.55565784083093333331e1 * t315 + t318 + t322 + 0.41674338062319999998e1 * t377 + 0.20837169031159999999e1 * t379 - t224 - t229 + t279 + t325 + 0.2e1 * t381 * t29 * t160 + t40 * t411 - t284;
545 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
546 v3rho3[2] = t375 + t413;
547
548 t421 = t158 * t66;
549 t426 = t164 * t61;
550 t432 = -0.6e1 * t100 - 0.6e1 * t247;
551 t435 = t171 * t63;
552 t440 = -t432;
553 t443 = 0.4e1 / 0.27e2 * t238 * t426 - t140 * t168 / 0.3e1 + t42 * t432 / 0.3e1 + 0.4e1 / 0.27e2 * t254 * t435 - t146 * t174 / 0.3e1 + t48 * t440 / 0.3e1;
554 t445 = t41 * t19 * t443;
555 t447 = t66 * t177;
556 t448 = t34 * t447;
557 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
558 v3rho3[3] = 0.4e1 / 0.3e1 * t290 + 0.33100228493624736458e1 * t310 + 0.8e1 * t362 + 0.4e1 * t368 - t206 + 0.83348676124639999996e1 * t315 + 0.12502301418695999999e2 * t377 + 0.62511507093479999997e1 * t379 - t224 - t229 + t279 + 0.2e1 * t38 * t421 * t29 * t233 + t40 * t445 - t284 + 0.6e1 * t133 * t448;
559
560 #ifndef XC_DONT_COMPILE_LXC
561
562 if(order < 4) return;
563
564
565 t451 = t91 * t264;
566 t453 = t230 * t58;
567 t455 = t12 * t119;
568 t458 = t455 * t86 * t26 * t124;
569 t461 = t18 * t20 * t217;
570 t462 = t197 * t461;
571 t463 = 0.3e-19 * t462;
572 t467 = t181 * t41 * t19 * t21 * t191;
573 t475 = t195 * t114 * t26 * t124;
574 t477 = t13 * t269;
575 t480 = t477 * t52 * t56 * t277;
576 t484 = t71 * t34 * t80 * t114;
577 t487 = t195 * t52 * t277;
578 t494 = 0.32691583697407147117e0 * t14 / t21 / t99 * t119 * t222;
579 t497 = t71 * t34 * t21 * t264;
580 t501 = t14 / t79 / t99;
581 t506 = 0.41544948398103724184e1 * t501 * t269 * t56 * t31 * t281;
582 t509 = 0.1e1 / t6 / t236 / t5;
583 t510 = t96 * t96;
584 t516 = t103 * t103;
585 t522 = 0.1e1 / t245 / t2;
586 t523 = t1 * t522;
587 t525 = -0.24e2 * t246 + 0.24e2 * t523;
588 t530 = 0.1e1 / t9 / t252 / t8;
589 t531 = t108 * t108;
590 t537 = t111 * t111;
591 t549 = 0.83348676124639999996e1 * t451 + 0.16669735224927999999e2 * t453 + 0.13240091397449894583e2 * t458 - t463 + 0.32e2 * t467 + 0.8e1 * t133 * t34 * t52 * t264 + 0.66200456987249472916e1 * t475 + 0.18695226779146675883e2 * t480 + 0.8e1 / 0.3e1 * t484 - 0.46738066947866689707e1 * t487 - t494 + 0.16e2 / 0.3e1 * t497 - t506 + t40 * t41 * t19 * (-0.28e2 / 0.81e2 * t509 * t510 + 0.8e1 / 0.9e1 * t238 * t96 * t103 - t95 * t516 / 0.3e1 - 0.4e1 / 0.9e1 * t242 * t249 + t42 * t525 / 0.3e1 - 0.28e2 / 0.81e2 * t530 * t531 + 0.8e1 / 0.9e1 * t254 * t108 * t111 - t107 * t537 / 0.3e1 - 0.4e1 / 0.9e1 * t258 * t261 - t48 * t525 / 0.3e1);
592 t552 = t182 * t34 * t80 * t86;
593 t556 = t71 * t34 * t226 * t52;
594 t560 = t114 * t114;
595 t567 = 0.10289960015387654321e1 * t14 * t100 * t58;
596 t570 = t230 * t29 * t31 * t35;
597 t574 = 0.40e2 / 0.243e3 * t32 * t34 * t275;
598 t578 = 0.10386237099525931046e1 * t501 * t119 * t273 * t19;
599 t579 = t12 * t3;
600 t580 = t58 * t86;
601 t581 = t579 * t580;
602 t583 = t13 * t44;
603 t584 = t583 * t201;
604 t586 = t207 * t191;
605 t589 = t200 * t58 * t114;
606 t591 = t118 * t118;
607 t597 = 0.94530758360525579704e1 * t14 / t591 * t56 * t26 * t246;
608 t600 = 0.47265379180262789851e1 * t270 * t26 * t246;
609 t601 = 0.16e2 / 0.3e1 * t552 - 0.32e2 / 0.27e2 * t556 + 0.12e2 * t297 * t116 + 0.6e1 * t85 * t41 * t19 * t560 + t567 + 0.32e2 / 0.3e1 * t570 + t574 + t578 + 0.33339470449855999998e2 * t581 - 0.74087712110791111108e1 * t584 + 0.50009205674783999996e2 * t586 + 0.16669735224927999999e2 * t589 + t597 - t600;
610 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
611 v4rho4[0] = t549 + t601;
612
613 t602 = t580 * t66;
614 t604 = t91 * t352;
615 t608 = t309 * t461;
616 t610 = t80 * t12;
617 t612 = t610 * t132 * t135;
618 t613 = 0.8e1 / 0.3e1 * t612;
619 t614 = t293 * t301;
620 t616 = t293 * t305;
621 t618 = t381 * t132;
622 t633 = 0.12502301418695999999e2 * t602 + 0.62511507093479999997e1 * t604 + 0.20837169031159999999e1 * t451 + 0.41674338062319999998e1 * t453 - 0.1e-18 * t608 + t613 + 0.16e2 * t614 + 0.8e1 * t616 + 0.6e1 * t618 * t305 + 0.6e1 * t133 * t34 * t352 * t52 + 0.6e1 * t133 * t34 * t152 * t114 + 0.2e1 * t133 * t34 * t66 * t264;
623 t636 = t195 * t152 * t26 * t124;
624 t637 = 0.33100228493624736458e1 * t636;
625 t640 = t477 * t66 * t56 * t277;
626 t643 = 0.13e-18 * t462;
627 t648 = t579 * t314 * t52;
628 t649 = 0.16669735224927999999e2 * t648;
629 t653 = t637 + 0.46738066947866689708e1 * t640 + 0.66200456987249472917e1 * t458 - t643 + 0.8e1 * t467 + 0.33100228493624736459e1 * t475 + 0.14021420084360006912e2 * t480 + t649 + 0.4e1 / 0.3e1 * t484 - 0.35053550210900017281e1 * t487 - t494 + 0.4e1 / 0.3e1 * t497;
630 t657 = t289 * t154;
631 t658 = 0.4e1 / 0.3e1 * t657;
632 t662 = t195 * t66 * t277;
633 t664 = t129 * t354;
634 t683 = 0.8e1 * t143 * t522;
635 t701 = 0.8e1 * t149 * t522;
636 t716 = -0.28e2 / 0.81e2 * t509 * t61 * t239 + 0.4e1 / 0.9e1 * t326 * t46 * t103 + 0.8e1 / 0.9e1 * t238 * t1 * t100 * t96 - 0.2e1 / 0.3e1 * t95 * t100 * t46 - 0.2e1 / 0.3e1 * t329 * t100 * t103 + t683 - 0.28e2 / 0.81e2 * t530 * t63 * t255 + 0.4e1 / 0.9e1 * t339 * t49 * t111 - 0.8e1 / 0.9e1 * t254 * t1 * t100 * t108 + 0.2e1 / 0.3e1 * t107 * t100 * t49 + 0.2e1 / 0.3e1 * t342 * t100 * t111 - t701 + 0.2e1 * t329 * t246 * t46 - 0.2e1 * t342 * t246 * t49 - t140 * t249 / 0.9e1 - 0.4e1 * t42 * t246 - t146 * t261 / 0.9e1 + 0.4e1 * t48 * t246;
637 t722 = t226 * t13 * t29 * t68;
638 t725 = t183 * t29 * t68;
639 t730 = t455 * t308 * t18 * t123 * t52;
640 t732 = -t506 + 0.8e1 / 0.3e1 * t552 - 0.8e1 / 0.9e1 * t556 + t658 + 0.6e1 * t297 * t154 - 0.11684516736966672427e1 * t662 + 0.4e1 * t664 + t40 * t41 * t19 * t716 - 0.8e1 / 0.27e2 * t722 + 0.8e1 * t725 + 0.66200456987249472917e1 * t730 + t567;
641 t738 = t583 * t314;
642 t740 = t207 * t300;
643 t742 = t207 * t304;
644 t745 = t200 * t58 * t152;
645 t746 = 0.83348676124639999996e1 * t745;
646 t747 = 0.8e1 / 0.3e1 * t570 + t574 + t578 + 0.16669735224927999999e2 * t581 - 0.55565784083093333330e1 * t584 + 0.12502301418695999999e2 * t586 + 0.83348676124639999995e1 * t589 + t597 - t600 - 0.18521928027697777777e1 * t738 + 0.25004602837391999999e2 * t740 + 0.12502301418695999999e2 * t742 + t746;
647 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
648 v4rho4[1] = t633 + t653 + t732 + t747;
649
650 t749 = t201 * t158;
651 t751 = t91 * t409;
652 t757 = t455 * t158 * t26 * t124;
653 t761 = t195 * t177 * t26 * t124;
654 t763 = t293 * t365;
655 t771 = t293 * t371;
656 t783 = 0.83348676124639999996e1 * t749 + 0.41674338062319999998e1 * t751 + 0.83348676124639999996e1 * t602 + 0.41674338062319999998e1 * t604 + 0.22066818995749824305e1 * t757 + 0.11033409497874912153e1 * t761 + 0.32e2 / 0.3e1 * t763 + 0.8e1 * t618 * t365 + 0.4e1 * t133 * t34 * t66 * t352 + 0.16e2 / 0.3e1 * t771 + 0.4e1 * t133 * t34 * t409 * t52 + 0.2e1 * t133 * t34 * t177 * t114 + 0.32e2 / 0.9e1 * t612 + 0.32e2 / 0.3e1 * t614;
657 t796 = 0.16e2 / 0.3e1 * t616 + 0.44133637991499648610e1 * t636 + 0.93476133895733379411e1 * t640 + 0.22066818995749824306e1 * t458 - t643 + 0.11033409497874912153e1 * t475 + 0.93476133895733379415e1 * t480 + 0.22226313633237333332e2 * t648 + 0.4e1 / 0.9e1 * t484 - 0.23369033473933344854e1 * t487 - t494 - t506 + 0.8e1 / 0.9e1 * t552 - 0.16e2 / 0.27e2 * t556 + 0.16e2 / 0.9e1 * t657;
658 t803 = t72 * t29 * t160;
659 t805 = t129 * t411;
660 t820 = t1 * t1;
661 t823 = 0.1e1 / t245 / t43;
662 t862 = -0.28e2 / 0.81e2 * t509 * t164 * t96 + 0.32e2 / 0.27e2 * t326 * t46 * t1 * t100 + 0.4e1 / 0.27e2 * t385 * t103 - 0.8e1 / 0.9e1 * t95 * t820 * t823 - 0.4e1 / 0.9e1 * t140 * t100 + 0.4e1 / 0.3e1 * t140 * t247 + 0.4e1 / 0.27e2 * t238 * t168 * t96 - 0.2e1 / 0.9e1 * t95 * t395 * t46 - t390 * t103 / 0.9e1 + t683 - 0.28e2 / 0.81e2 * t530 * t171 * t108 - 0.32e2 / 0.27e2 * t339 * t49 * t1 * t100 + 0.4e1 / 0.27e2 * t398 * t111 - 0.8e1 / 0.9e1 * t107 * t820 * t823 + 0.4e1 / 0.9e1 * t146 * t100 - 0.4e1 / 0.3e1 * t146 * t247 + 0.4e1 / 0.27e2 * t254 * t174 * t108 - 0.2e1 / 0.9e1 * t107 * t406 * t49 - t403 * t111 / 0.9e1 - t701;
663 t867 = t610 * t29 * t160;
664 t869 = t152 * t152;
665 t874 = t289 * t179;
666 t879 = -0.23369033473933344853e1 * t662 + 0.8e1 / 0.3e1 * t664 - 0.16e2 / 0.27e2 * t722 + 0.16e2 / 0.3e1 * t725 + 0.16e2 / 0.3e1 * t803 + 0.8e1 / 0.3e1 * t805 + 0.2e1 * t38 * t114 * t29 * t160 + t40 * t41 * t19 * t862 + 0.8e1 / 0.9e1 * t867 + 0.4e1 * t85 * t41 * t19 * t869 + 0.4e1 / 0.9e1 * t874 + 0.2e1 * t297 * t179 + 0.88267275982999297221e1 * t730 + t567;
667 t888 = t579 * t58 * t158;
668 t890 = t207 * t364;
669 t893 = t200 * t58 * t177;
670 t895 = t207 * t370;
671 t897 = t574 + t578 + 0.55565784083093333330e1 * t581 - 0.37043856055395555554e1 * t584 + 0.27782892041546666665e1 * t589 + t597 - t600 - 0.37043856055395555554e1 * t738 + 0.16669735224927999999e2 * t740 + 0.83348676124639999996e1 * t742 + 0.11113156816618666666e2 * t745 + 0.55565784083093333331e1 * t888 + 0.16669735224927999999e2 * t890 + 0.27782892041546666665e1 * t893 + 0.83348676124639999996e1 * t895;
672 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
673 v4rho4[2] = t783 + t796 + t879 + t897;
674
675 t899 = t421 * t58;
676 t901 = t91 * t443;
677 t911 = t293 * t448;
678 t924 = 0.41674338062319999998e1 * t899 + 0.20837169031159999999e1 * t901 + 0.12502301418695999999e2 * t749 + 0.62511507093479999997e1 * t751 + 0.6e1 * t133 * t34 * t66 * t409 + 0.66200456987249472914e1 * t757 + 0.33100228493624736458e1 * t761 + 0.8e1 * t911 + 0.2e1 * t133 * t34 * t443 * t52 + 0.6e1 * t618 * t448 + 0.6e1 * t133 * t34 * t152 * t177 + 0.16e2 * t763;
679 t930 = 0.8e1 * t771 + t613 + t637 + 0.14021420084360006912e2 * t640 - t463 + 0.46738066947866689707e1 * t480 + t649 - 0.11684516736966672427e1 * t487 - t494 - t506 - 0.8e1 / 0.27e2 * t556 + t658;
680 t939 = t38 * t158 * t29;
681 t960 = 0.12e2 * t246 + 0.24e2 * t523;
682 t982 = -0.28e2 / 0.81e2 * t509 * t426 * t46 + 0.8e1 / 0.9e1 * t385 * t101 + 0.4e1 / 0.9e1 * t326 * t168 * t46 - 0.2e1 / 0.3e1 * t329 * t100 * t168 - t140 * t395 / 0.3e1 - t95 * t432 * t46 / 0.9e1 + t42 * t960 / 0.3e1 - 0.28e2 / 0.81e2 * t530 * t435 * t49 - 0.8e1 / 0.9e1 * t398 * t101 + 0.4e1 / 0.9e1 * t339 * t174 * t49 + 0.2e1 / 0.3e1 * t342 * t100 * t174 - t146 * t406 / 0.3e1 - t107 * t440 * t49 / 0.9e1 - t48 * t960 / 0.3e1;
683 t986 = t129 * t445;
684 t989 = -0.35053550210900017280e1 * t662 - 0.8e1 / 0.9e1 * t722 + 0.8e1 * t803 + 0.4e1 * t805 + 0.8e1 / 0.3e1 * t867 + 0.4e1 / 0.3e1 * t874 + 0.6e1 * t939 * t154 + t40 * t41 * t19 * t982 + 0.4e1 / 0.3e1 * t986 + 0.66200456987249472916e1 * t730 + t567 + t574;
685 t992 = t21 * t421 * t29 * t233;
686 t1000 = t207 * t447;
687 t1002 = t578 + 0.8e1 / 0.3e1 * t992 - 0.18521928027697777777e1 * t584 + t597 - t600 - 0.55565784083093333331e1 * t738 + t746 + 0.16669735224927999999e2 * t888 + 0.25004602837391999998e2 * t890 + 0.83348676124639999996e1 * t893 + 0.12502301418695999999e2 * t895 + 0.12502301418695999999e2 * t1000;
688 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
689 v4rho4[3] = t924 + t930 + t989 + t1002;
690
691 t1019 = 0.16669735224927999999e2 * t899 + 0.83348676124639999996e1 * t901 + 0.8e1 * t133 * t34 * t443 * t66 + 0.13240091397449894583e2 * t757 + 0.66200456987249472916e1 * t761 + 0.32e2 * t911 - 0.3e-19 * t608 + 0.18695226779146675883e2 * t640 - t494 - t506 - 0.46738066947866689707e1 * t662 - 0.32e2 / 0.27e2 * t722 + 0.16e2 / 0.3e1 * t867 + 0.8e1 / 0.3e1 * t874;
692 t1021 = t164 * t164;
693 t1026 = t168 * t168;
694 t1032 = 0.24e2 * t246 + 0.24e2 * t523;
695 t1035 = t171 * t171;
696 t1040 = t174 * t174;
697 t1054 = t177 * t177;
698 t1064 = 0.16e2 / 0.3e1 * t986 + t40 * t41 * t19 * (-0.28e2 / 0.81e2 * t509 * t1021 + 0.8e1 / 0.9e1 * t385 * t168 - t95 * t1026 / 0.3e1 - 0.4e1 / 0.9e1 * t140 * t432 + t42 * t1032 / 0.3e1 - 0.28e2 / 0.81e2 * t530 * t1035 + 0.8e1 / 0.9e1 * t398 * t174 - t107 * t1040 / 0.3e1 - 0.4e1 / 0.9e1 * t146 * t440 - t48 * t1032 / 0.3e1) + 0.12e2 * t939 * t179 + 0.6e1 * t85 * t41 * t19 * t1054 + t567 + t574 + t578 + 0.32e2 / 0.3e1 * t992 + t597 - t600 - 0.74087712110791111108e1 * t738 + 0.33339470449855999998e2 * t888 + 0.16669735224927999999e2 * t893 + 0.50009205674783999996e2 * t1000;
699 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
700 v4rho4[4] = t1019 + t1064;
701
702 #ifndef XC_DONT_COMPILE_MXC
703
704 if(order < 5) return;
705
706
707 #endif
708
709 #endif
710
711 #endif
712
713 #endif
714
715 #endif
716
717
718 }
719
720