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_chachiyo_mod.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, t5, t7, t8, t9, t13;
24 double t14, t16, t17, t18, t22, t23;
25
26 #ifndef XC_DONT_COMPILE_VXC
27 double t24, t33, t34;
28
29 #ifndef XC_DONT_COMPILE_FXC
30 double t49, t52, t53, t54;
31
32 #ifndef XC_DONT_COMPILE_KXC
33 double t57, t63, t74, t81, t83;
34
35 #ifndef XC_DONT_COMPILE_LXC
36 double t96, t118, t122, t123;
37 #endif
38
39 #endif
40
41 #endif
42
43 #endif
44
45 #endif
46
47
48 lda_c_chachiyo_mod_params *params;
49
50 assert(p->params != NULL);
51 params = (lda_c_chachiyo_mod_params * )(p->params);
52
53 t1 = M_CBRT3;
54 t2 = t1 * t1;
55 t3 = params->bp * t2;
56 t5 = POW_1_3(0.1e1 / M_PI);
57 t7 = M_CBRT4;
58 t8 = 0.1e1 / t5 * t7;
59 t9 = POW_1_3(rho[0]);
60 t13 = params->bp * t1;
61 t14 = t5 * t5;
62 t16 = t7 * t7;
63 t17 = 0.1e1 / t14 * t16;
64 t18 = t9 * t9;
65 t22 = 0.1e1 + t3 * t8 * t9 / 0.3e1 + t13 * t17 * t18 / 0.3e1;
66 t23 = log(t22);
67 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
68 zk[0] = params->ap * t23;
69
70 #ifndef XC_DONT_COMPILE_VXC
71
72 if(order < 1) return;
73
74
75 t24 = rho[0] * params->ap;
76 t33 = t3 * t8 / t18 / 0.9e1 + 0.2e1 / 0.9e1 * t13 * t17 / t9;
77 t34 = 0.1e1 / t22;
78 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
79 vrho[0] = t24 * t33 * t34 + (params->ap * t23);
80
81 #ifndef XC_DONT_COMPILE_FXC
82
83 if(order < 2) return;
84
85
86 t49 = -0.2e1 / 0.27e2 * t3 * t8 / t18 / rho[0] - 0.2e1 / 0.27e2 * t13 * t17 / t9 / rho[0];
87 t52 = t33 * t33;
88 t53 = t22 * t22;
89 t54 = 0.1e1 / t53;
90 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
91 v2rho2[0] = t24 * t49 * t34 - t24 * t52 * t54 + 0.2e1 * params->ap * t33 * t34;
92
93 #ifndef XC_DONT_COMPILE_KXC
94
95 if(order < 3) return;
96
97
98 t57 = params->ap * t49;
99 t63 = rho[0] * rho[0];
100 t74 = 0.10e2 / 0.81e2 * t3 * t8 / t18 / t63 + 0.8e1 / 0.81e2 * t13 * t17 / t9 / t63;
101 t81 = t52 * t33;
102 t83 = 0.1e1 / t53 / t22;
103 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
104 v3rho3[0] = -0.3e1 * t24 * t49 * t54 * t33 + t24 * t74 * t34 + 0.2e1 * t24 * t81 * t83 - 0.3e1 * params->ap * t52 * t54 + 0.3e1 * t57 * t34;
105
106 #ifndef XC_DONT_COMPILE_LXC
107
108 if(order < 4) return;
109
110
111 t96 = t63 * rho[0];
112 t118 = t49 * t49;
113 t122 = t52 * t52;
114 t123 = t53 * t53;
115 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
116 v4rho4[0] = 0.4e1 * params->ap * t74 * t34 - 0.12e2 * t57 * t54 * t33 + 0.8e1 * params->ap * t81 * t83 + t24 * (-0.80e2 / 0.243e3 * t3 * t8 / t18 / t96 - 0.56e2 / 0.243e3 * t13 * t17 / t9 / t96) * t34 - 0.4e1 * t24 * t74 * t54 * t33 + 0.12e2 * t24 * t49 * t83 * t52 - 0.3e1 * t24 * t118 * t54 - 0.6e1 * t24 * t122 / t123;
117
118 #ifndef XC_DONT_COMPILE_MXC
119
120 if(order < 5) return;
121
122
123 #endif
124
125 #endif
126
127 #endif
128
129 #endif
130
131 #endif
132
133
134 }
135
136
137 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))138 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
139 {
140
141 #ifndef XC_DONT_COMPILE_EXC
142 double t1, t2, t3, t5, t7, t8, t9, t13;
143 double t14, t16, t17, t18, t22, t23;
144
145 #ifndef XC_DONT_COMPILE_VXC
146 double t24, t33, t34;
147
148 #ifndef XC_DONT_COMPILE_FXC
149 double t49, t52, t53, t54;
150
151 #ifndef XC_DONT_COMPILE_KXC
152 double t57, t63, t74, t81, t83;
153
154 #ifndef XC_DONT_COMPILE_LXC
155 double t96, t118, t122, t123;
156 #endif
157
158 #endif
159
160 #endif
161
162 #endif
163
164 #endif
165
166
167 lda_c_chachiyo_mod_params *params;
168
169 assert(p->params != NULL);
170 params = (lda_c_chachiyo_mod_params * )(p->params);
171
172 t1 = M_CBRT3;
173 t2 = t1 * t1;
174 t3 = params->bf * t2;
175 t5 = POW_1_3(0.1e1 / M_PI);
176 t7 = M_CBRT4;
177 t8 = 0.1e1 / t5 * t7;
178 t9 = POW_1_3(rho[0]);
179 t13 = params->bf * t1;
180 t14 = t5 * t5;
181 t16 = t7 * t7;
182 t17 = 0.1e1 / t14 * t16;
183 t18 = t9 * t9;
184 t22 = 0.1e1 + t3 * t8 * t9 / 0.3e1 + t13 * t17 * t18 / 0.3e1;
185 t23 = log(t22);
186 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
187 zk[0] = params->af * t23;
188
189 #ifndef XC_DONT_COMPILE_VXC
190
191 if(order < 1) return;
192
193
194 t24 = rho[0] * params->af;
195 t33 = t3 * t8 / t18 / 0.9e1 + 0.2e1 / 0.9e1 * t13 * t17 / t9;
196 t34 = 0.1e1 / t22;
197 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
198 vrho[0] = t24 * t33 * t34 + (params->af * t23);
199
200 #ifndef XC_DONT_COMPILE_FXC
201
202 if(order < 2) return;
203
204
205 t49 = -0.2e1 / 0.27e2 * t3 * t8 / t18 / rho[0] - 0.2e1 / 0.27e2 * t13 * t17 / t9 / rho[0];
206 t52 = t33 * t33;
207 t53 = t22 * t22;
208 t54 = 0.1e1 / t53;
209 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
210 v2rho2[0] = t24 * t49 * t34 - t24 * t52 * t54 + 0.2e1 * params->af * t33 * t34;
211
212 #ifndef XC_DONT_COMPILE_KXC
213
214 if(order < 3) return;
215
216
217 t57 = params->af * t49;
218 t63 = rho[0] * rho[0];
219 t74 = 0.10e2 / 0.81e2 * t3 * t8 / t18 / t63 + 0.8e1 / 0.81e2 * t13 * t17 / t9 / t63;
220 t81 = t52 * t33;
221 t83 = 0.1e1 / t53 / t22;
222 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
223 v3rho3[0] = -0.3e1 * t24 * t49 * t54 * t33 + t24 * t74 * t34 + 0.2e1 * t24 * t81 * t83 - 0.3e1 * params->af * t52 * t54 + 0.3e1 * t57 * t34;
224
225 #ifndef XC_DONT_COMPILE_LXC
226
227 if(order < 4) return;
228
229
230 t96 = t63 * rho[0];
231 t118 = t49 * t49;
232 t122 = t52 * t52;
233 t123 = t53 * t53;
234 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
235 v4rho4[0] = 0.4e1 * params->af * t74 * t34 - 0.12e2 * t57 * t54 * t33 + 0.8e1 * params->af * t81 * t83 + t24 * (-0.80e2 / 0.243e3 * t3 * t8 / t18 / t96 - 0.56e2 / 0.243e3 * t13 * t17 / t9 / t96) * t34 - 0.4e1 * t24 * t74 * t54 * t33 + 0.12e2 * t24 * t49 * t83 * t52 - 0.3e1 * t24 * t118 * t54 - 0.6e1 * t24 * t122 / t123;
236
237 #ifndef XC_DONT_COMPILE_MXC
238
239 if(order < 5) return;
240
241
242 #endif
243
244 #endif
245
246 #endif
247
248 #endif
249
250 #endif
251
252
253 }
254
255
256 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))257 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
258 {
259
260 #ifndef XC_DONT_COMPILE_EXC
261 double t1, t2, t3, t5, t7, t8, t9, t10;
262 double t11, t14, t15, t17, t18, t19, t20, t23;
263 double t24, t25, t26, t29, t32, t33, t35, t36;
264 double t37, t38, t39, t40, t41, t42, t43, t44;
265 double t46, t47, t50, t51;
266
267 #ifndef XC_DONT_COMPILE_VXC
268 double t53, t57, t60, t62, t63, t68, t70, t72;
269 double t73, t74, t75, t76, t77, t78, t79, t81;
270 double t82, t85, t86, t87, t90, t92, t95, t96;
271 double t97;
272
273 #ifndef XC_DONT_COMPILE_FXC
274 double t100, t101, t105, t109, t112, t113, t114, t115;
275 double t117, t118, t119, t123, t124, t126, t128, t129;
276 double t131, t132, t133, t134, t135, t136, t137, t138;
277 double t139, t141, t142, t145, t146, t147, t149, t153;
278 double t154, t157, t160, t161, t162, t166, t168, t169;
279 double t171, t174, t177, t180, t183, t184, t189, t190;
280 double t191, t192, t193, t197, t200, t203, t206, t207;
281 double t208;
282
283 #ifndef XC_DONT_COMPILE_KXC
284 double t211, t212, t213, t219, t224, t228, t229, t230;
285 double t231, t232, t236, t237, t238, t244, t246, t252;
286 double t255, t256, t257, t258, t259, t260, t261, t263;
287 double t265, t266, t267, t268, t269, t270, t271, t273;
288 double t274, t277, t280, t281, t282, t284, t287, t289;
289 double t290, t293, t296, t299, t300, t301, t305, t306;
290 double t310, t312, t313, t314, t315, t316, t317, t319;
291 double t320, t322, t323, t325, t328, t338, t341, t351;
292 double t352, t354, t359, t361, t362, t364, t365, t367;
293 double t369, t370, t372, t377, t382, t385, t390, t393;
294 double t396, t397, t399, t404, t407, t408, t409, t410;
295 double t411, t412, t413, t419, t422, t427, t430, t431;
296 double t432;
297
298 #ifndef XC_DONT_COMPILE_LXC
299 double t435, t438, t443, t454, t458, t460, t470, t472;
300 double t475, t476, t479, t480, t482, t485, t487, t488;
301 double t489, t491, t493, t499, t500, t502, t504, t506;
302 double t512, t513, t519, t525, t526, t528, t533, t534;
303 double t540, t551, t553, t554, t557, t564, t565, t566;
304 double t568, t581, t587, t589, t591, t592, t593, t620;
305 double t638, t647, t650, t654, t660, t671, t674, t713;
306 double t716, t723, t734, t742, t747, t750, t755, t757;
307 double t759, t779, t780, t782, t802, t824, t827, t830;
308 double t845, t852, t864, t876, t881, t887, t890, t895;
309 double t906;
310 #endif
311
312 #endif
313
314 #endif
315
316 #endif
317
318 #endif
319
320
321 lda_c_chachiyo_mod_params *params;
322
323 assert(p->params != NULL);
324 params = (lda_c_chachiyo_mod_params * )(p->params);
325
326 t1 = M_CBRT3;
327 t2 = t1 * t1;
328 t3 = params->bp * t2;
329 t5 = POW_1_3(0.1e1 / M_PI);
330 t7 = M_CBRT4;
331 t8 = 0.1e1 / t5 * t7;
332 t9 = rho[0] + rho[1];
333 t10 = POW_1_3(t9);
334 t11 = t8 * t10;
335 t14 = params->bp * t1;
336 t15 = t5 * t5;
337 t17 = t7 * t7;
338 t18 = 0.1e1 / t15 * t17;
339 t19 = t10 * t10;
340 t20 = t18 * t19;
341 t23 = 0.1e1 + t3 * t11 / 0.3e1 + t14 * t20 / 0.3e1;
342 t24 = log(t23);
343 t25 = params->ap * t24;
344 t26 = params->bf * t2;
345 t29 = params->bf * t1;
346 t32 = 0.1e1 + t26 * t11 / 0.3e1 + t29 * t20 / 0.3e1;
347 t33 = log(t32);
348 t35 = params->af * t33 - t25;
349 t36 = rho[0] - rho[1];
350 t37 = 0.1e1 / t9;
351 t38 = t36 * t37;
352 t39 = 0.1e1 + t38;
353 t40 = POW_1_3(t39);
354 t41 = t40 * t40;
355 t42 = 0.1e1 - t38;
356 t43 = POW_1_3(t42);
357 t44 = t43 * t43;
358 t46 = t41 / 0.2e1 + t44 / 0.2e1;
359 t47 = t46 * t46;
360 t50 = -0.2e1 * t47 * t46 + 0.2e1;
361 t51 = t35 * t50;
362 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
363 zk[0] = t25 + t51;
364
365 #ifndef XC_DONT_COMPILE_VXC
366
367 if(order < 1) return;
368
369
370 t53 = t8 / t19;
371 t57 = t18 / t10;
372 t60 = t3 * t53 / 0.9e1 + 0.2e1 / 0.9e1 * t14 * t57;
373 t62 = 0.1e1 / t23;
374 t63 = params->ap * t60 * t62;
375 t68 = t26 * t53 / 0.9e1 + 0.2e1 / 0.9e1 * t29 * t57;
376 t70 = 0.1e1 / t32;
377 t72 = params->af * t68 * t70 - t63;
378 t73 = t72 * t50;
379 t74 = t35 * t47;
380 t75 = 0.1e1 / t40;
381 t76 = t9 * t9;
382 t77 = 0.1e1 / t76;
383 t78 = t36 * t77;
384 t79 = t37 - t78;
385 t81 = 0.1e1 / t43;
386 t82 = -t79;
387 t85 = t75 * t79 / 0.3e1 + t81 * t82 / 0.3e1;
388 t86 = t74 * t85;
389 t87 = 0.6e1 * t86;
390 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
391 vrho[0] = t25 + t51 + t9 * (t63 + t73 - t87);
392
393 t90 = -t37 - t78;
394 t92 = -t90;
395 t95 = t75 * t90 / 0.3e1 + t81 * t92 / 0.3e1;
396 t96 = t74 * t95;
397 t97 = 0.6e1 * t96;
398 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
399 vrho[1] = t25 + t51 + t9 * (t63 + t73 - t97);
400
401 #ifndef XC_DONT_COMPILE_FXC
402
403 if(order < 2) return;
404
405
406 t100 = 0.2e1 * t63;
407 t101 = 0.2e1 * t73;
408 t105 = t8 / t19 / t9;
409 t109 = t18 / t10 / t9;
410 t112 = -0.2e1 / 0.27e2 * t3 * t105 - 0.2e1 / 0.27e2 * t14 * t109;
411 t113 = params->ap * t112;
412 t114 = t113 * t62;
413 t115 = t60 * t60;
414 t117 = t23 * t23;
415 t118 = 0.1e1 / t117;
416 t119 = params->ap * t115 * t118;
417 t123 = -0.2e1 / 0.27e2 * t26 * t105 - 0.2e1 / 0.27e2 * t29 * t109;
418 t124 = params->af * t123;
419 t126 = t68 * t68;
420 t128 = t32 * t32;
421 t129 = 0.1e1 / t128;
422 t131 = -params->af * t126 * t129 + t124 * t70 - t114 + t119;
423 t132 = t131 * t50;
424 t133 = t72 * t47;
425 t134 = t133 * t85;
426 t135 = 0.12e2 * t134;
427 t136 = t35 * t46;
428 t137 = t85 * t85;
429 t138 = t136 * t137;
430 t139 = 0.12e2 * t138;
431 t141 = 0.1e1 / t40 / t39;
432 t142 = t79 * t79;
433 t145 = t76 * t9;
434 t146 = 0.1e1 / t145;
435 t147 = t36 * t146;
436 t149 = -0.2e1 * t77 + 0.2e1 * t147;
437 t153 = 0.1e1 / t43 / t42;
438 t154 = t82 * t82;
439 t157 = -t149;
440 t160 = -t141 * t142 / 0.9e1 + t75 * t149 / 0.3e1 - t153 * t154 / 0.9e1 + t81 * t157 / 0.3e1;
441 t161 = t74 * t160;
442 t162 = 0.6e1 * t161;
443 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
444 v2rho2[0] = t100 + t101 - 0.12e2 * t86 + t9 * (t114 - t119 + t132 - t135 - t139 - t162);
445
446 t166 = t133 * t95;
447 t168 = t95 * t85;
448 t169 = t136 * t168;
449 t171 = t141 * t90;
450 t174 = t75 * t36;
451 t177 = t153 * t92;
452 t180 = t81 * t36;
453 t183 = -t171 * t79 / 0.9e1 + 0.2e1 / 0.3e1 * t174 * t146 - t177 * t82 / 0.9e1 - 0.2e1 / 0.3e1 * t180 * t146;
454 t184 = t74 * t183;
455 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
456 v2rho2[1] = t100 + t101 - t87 - t97 + t9 * (t114 - t119 + t132 - 0.6e1 * t134 - 0.6e1 * t166 - 0.12e2 * t169 - 0.6e1 * t184);
457
458 t189 = 0.12e2 * t166;
459 t190 = t95 * t95;
460 t191 = t136 * t190;
461 t192 = 0.12e2 * t191;
462 t193 = t90 * t90;
463 t197 = 0.2e1 * t77 + 0.2e1 * t147;
464 t200 = t92 * t92;
465 t203 = -t197;
466 t206 = -t141 * t193 / 0.9e1 + t75 * t197 / 0.3e1 - t153 * t200 / 0.9e1 + t81 * t203 / 0.3e1;
467 t207 = t74 * t206;
468 t208 = 0.6e1 * t207;
469 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
470 v2rho2[2] = t100 + t101 - 0.12e2 * t96 + t9 * (t114 - t119 + t132 - t189 - t192 - t208);
471
472 #ifndef XC_DONT_COMPILE_KXC
473
474 if(order < 3) return;
475
476
477 t211 = 0.3e1 * t114;
478 t212 = 0.3e1 * t119;
479 t213 = 0.3e1 * t132;
480 t219 = t8 / t19 / t76;
481 t224 = t18 / t10 / t76;
482 t228 = params->ap * (0.10e2 / 0.81e2 * t3 * t219 + 0.8e1 / 0.81e2 * t14 * t224);
483 t229 = t228 * t62;
484 t230 = t118 * t60;
485 t231 = t113 * t230;
486 t232 = 0.3e1 * t231;
487 t236 = 0.1e1 / t117 / t23;
488 t237 = params->ap * t115 * t60 * t236;
489 t238 = 0.2e1 * t237;
490 t244 = params->af * (0.10e2 / 0.81e2 * t26 * t219 + 0.8e1 / 0.81e2 * t29 * t224);
491 t246 = t129 * t68;
492 t252 = 0.1e1 / t128 / t32;
493 t255 = 0.2e1 * params->af * t126 * t68 * t252 - 0.3e1 * t124 * t246 + t244 * t70 - t229 + t232 - t238;
494 t256 = t255 * t50;
495 t257 = t131 * t47;
496 t258 = t257 * t85;
497 t259 = 0.18e2 * t258;
498 t260 = t72 * t46;
499 t261 = t260 * t137;
500 t263 = t133 * t160;
501 t265 = t137 * t85;
502 t266 = t35 * t265;
503 t267 = 0.12e2 * t266;
504 t268 = t85 * t160;
505 t269 = t136 * t268;
506 t270 = 0.36e2 * t269;
507 t271 = t39 * t39;
508 t273 = 0.1e1 / t40 / t271;
509 t274 = t142 * t79;
510 t277 = t141 * t79;
511 t280 = t76 * t76;
512 t281 = 0.1e1 / t280;
513 t282 = t36 * t281;
514 t284 = 0.6e1 * t146 - 0.6e1 * t282;
515 t287 = t42 * t42;
516 t289 = 0.1e1 / t43 / t287;
517 t290 = t154 * t82;
518 t293 = t153 * t82;
519 t296 = -t284;
520 t299 = 0.4e1 / 0.27e2 * t273 * t274 - t277 * t149 / 0.3e1 + t75 * t284 / 0.3e1 + 0.4e1 / 0.27e2 * t289 * t290 - t293 * t157 / 0.3e1 + t81 * t296 / 0.3e1;
521 t300 = t74 * t299;
522 t301 = 0.6e1 * t300;
523 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
524 v3rho3[0] = t211 - t212 + t213 - 0.36e2 * t134 - 0.36e2 * t138 - 0.18e2 * t161 + t9 * (t229 - t232 + t238 + t256 - t259 - 0.36e2 * t261 - 0.18e2 * t263 - t267 - t270 - t301);
525
526 t305 = 0.24e2 * t169;
527 t306 = 0.12e2 * t184;
528 t310 = t257 * t95;
529 t312 = t260 * t168;
530 t313 = 0.24e2 * t312;
531 t314 = t133 * t183;
532 t315 = 0.12e2 * t314;
533 t316 = t35 * t137;
534 t317 = t316 * t95;
535 t319 = t183 * t85;
536 t320 = t136 * t319;
537 t322 = t95 * t160;
538 t323 = t136 * t322;
539 t325 = t273 * t90;
540 t328 = t141 * t36;
541 t338 = t289 * t92;
542 t341 = t153 * t36;
543 t351 = 0.4e1 / 0.27e2 * t325 * t142 - 0.4e1 / 0.9e1 * t328 * t146 * t79 - t171 * t149 / 0.9e1 + 0.2e1 / 0.3e1 * t75 * t146 - 0.2e1 * t174 * t281 + 0.4e1 / 0.27e2 * t338 * t154 + 0.4e1 / 0.9e1 * t341 * t146 * t82 - t177 * t157 / 0.9e1 - 0.2e1 / 0.3e1 * t81 * t146 + 0.2e1 * t180 * t281;
544 t352 = t74 * t351;
545 t354 = t229 - t232 + t238 + t256 - 0.12e2 * t258 - 0.12e2 * t261 - 0.6e1 * t263 - 0.6e1 * t310 - t313 - t315 - 0.12e2 * t317 - 0.24e2 * t320 - 0.12e2 * t323 - 0.6e1 * t352;
546 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
547 v3rho3[1] = t9 * t354 - 0.24e2 * t134 - t139 - t162 - t189 + t211 - t212 + t213 - t305 - t306;
548
549 t359 = t260 * t190;
550 t361 = t35 * t85;
551 t362 = t361 * t190;
552 t364 = t95 * t183;
553 t365 = t136 * t364;
554 t367 = t133 * t206;
555 t369 = t206 * t85;
556 t370 = t136 * t369;
557 t372 = t273 * t193;
558 t377 = t141 * t197;
559 t382 = -0.2e1 * t146 - 0.6e1 * t282;
560 t385 = t289 * t200;
561 t390 = t153 * t203;
562 t393 = -t382;
563 t396 = 0.4e1 / 0.27e2 * t372 * t79 - 0.4e1 / 0.9e1 * t171 * t147 - t377 * t79 / 0.9e1 + t75 * t382 / 0.3e1 + 0.4e1 / 0.27e2 * t385 * t82 + 0.4e1 / 0.9e1 * t177 * t147 - t390 * t82 / 0.9e1 + t81 * t393 / 0.3e1;
564 t397 = t74 * t396;
565 t399 = t229 - t232 + t238 + t256 - 0.6e1 * t258 - 0.12e2 * t310 - t313 - t315 - 0.12e2 * t359 - 0.12e2 * t362 - 0.24e2 * t365 - 0.6e1 * t367 - 0.12e2 * t370 - 0.6e1 * t397;
566 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
567 v3rho3[2] = t9 * t399 - t135 - 0.24e2 * t166 - t192 - t208 + t211 - t212 + t213 - t305 - t306;
568
569 t404 = 0.18e2 * t310;
570 t407 = t190 * t95;
571 t408 = t35 * t407;
572 t409 = 0.12e2 * t408;
573 t410 = t95 * t206;
574 t411 = t136 * t410;
575 t412 = 0.36e2 * t411;
576 t413 = t193 * t90;
577 t419 = -0.6e1 * t146 - 0.6e1 * t282;
578 t422 = t200 * t92;
579 t427 = -t419;
580 t430 = 0.4e1 / 0.27e2 * t273 * t413 - t171 * t197 / 0.3e1 + t75 * t419 / 0.3e1 + 0.4e1 / 0.27e2 * t289 * t422 - t177 * t203 / 0.3e1 + t81 * t427 / 0.3e1;
581 t431 = t74 * t430;
582 t432 = 0.6e1 * t431;
583 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
584 v3rho3[3] = t211 - t212 + t213 - 0.36e2 * t166 - 0.36e2 * t191 - 0.18e2 * t207 + t9 * (t229 - t232 + t238 + t256 - t404 - 0.36e2 * t359 - 0.18e2 * t367 - t409 - t412 - t432);
585
586 #ifndef XC_DONT_COMPILE_LXC
587
588 if(order < 4) return;
589
590
591 t435 = 0.4e1 * t256;
592 t438 = t8 / t19 / t145;
593 t443 = t18 / t10 / t145;
594 t454 = t123 * t123;
595 t458 = t126 * t126;
596 t460 = t128 * t128;
597 t470 = params->ap * (-0.80e2 / 0.243e3 * t3 * t438 - 0.56e2 / 0.243e3 * t14 * t443) * t62;
598 t472 = 0.4e1 * t228 * t230;
599 t475 = 0.12e2 * t113 * t236 * t115;
600 t476 = t112 * t112;
601 t479 = 0.3e1 * params->ap * t476 * t118;
602 t480 = t115 * t115;
603 t482 = t117 * t117;
604 t485 = 0.6e1 * params->ap * t480 / t482;
605 t487 = (params->af * (-0.80e2 / 0.243e3 * t26 * t438 - 0.56e2 / 0.243e3 * t29 * t443) * t70 - 0.4e1 * t244 * t246 + 0.12e2 * t124 * t252 * t126 - 0.3e1 * params->af * t454 * t129 - 0.6e1 * params->af * t458 / t460 - t470 + t472 - t475 + t479 + t485) * t50;
606 t488 = t131 * t46;
607 t489 = t488 * t137;
608 t491 = t260 * t268;
609 t493 = t160 * t160;
610 t499 = t255 * t47;
611 t500 = t499 * t85;
612 t502 = t257 * t160;
613 t504 = t72 * t265;
614 t506 = t133 * t299;
615 t512 = 0.1e1 / t40 / t271 / t39;
616 t513 = t142 * t142;
617 t519 = t149 * t149;
618 t525 = 0.1e1 / t280 / t9;
619 t526 = t36 * t525;
620 t528 = -0.24e2 * t281 + 0.24e2 * t526;
621 t533 = 0.1e1 / t43 / t287 / t42;
622 t534 = t154 * t154;
623 t540 = t157 * t157;
624 t551 = t487 + t475 - t472 - t479 - t485 - 0.72e2 * t489 - 0.144e3 * t491 - 0.36e2 * t136 * t493 - 0.48e2 * t136 * t85 * t299 + t470 - 0.24e2 * t500 - 0.36e2 * t502 - 0.48e2 * t504 - 0.24e2 * t506 - 0.72e2 * t316 * t160 - 0.6e1 * t74 * (-0.28e2 / 0.81e2 * t512 * t513 + 0.8e1 / 0.9e1 * t273 * t142 * t149 - t141 * t519 / 0.3e1 - 0.4e1 / 0.9e1 * t277 * t284 + t75 * t528 / 0.3e1 - 0.28e2 / 0.81e2 * t533 * t534 + 0.8e1 / 0.9e1 * t289 * t154 * t157 - t153 * t540 / 0.3e1 - 0.4e1 / 0.9e1 * t293 * t296 - t81 * t528 / 0.3e1);
625 t553 = 0.12e2 * t231;
626 t554 = 0.8e1 * t237;
627 t557 = 0.4e1 * t229;
628 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
629 v4rho4[0] = t9 * t551 - 0.72e2 * t258 - 0.144e3 * t261 - 0.72e2 * t263 - 0.48e2 * t266 - 0.144e3 * t269 - 0.24e2 * t300 + t435 - t553 + t554 + t557;
630
631 t564 = t488 * t168;
632 t565 = 0.36e2 * t564;
633 t566 = t260 * t319;
634 t568 = t260 * t322;
635 t581 = -0.36e2 * t136 * t183 * t160 - 0.12e2 * t136 * t95 * t299 - 0.36e2 * t136 * t351 * t85 - 0.36e2 * t361 * t322 - t472 + t475 + t487 - 0.36e2 * t491 - 0.12e2 * t504 - t565 - 0.72e2 * t566 - 0.36e2 * t568;
636 t587 = t72 * t137 * t95;
637 t589 = t499 * t95;
638 t591 = t257 * t183;
639 t592 = 0.18e2 * t591;
640 t593 = t133 * t351;
641 t620 = 0.8e1 * t174 * t525;
642 t638 = 0.8e1 * t180 * t525;
643 t647 = 0.2e1 * t328 * t281 * t79 - 0.2e1 * t341 * t281 * t82 - 0.28e2 / 0.81e2 * t512 * t90 * t274 + 0.4e1 / 0.9e1 * t325 * t79 * t149 + 0.8e1 / 0.9e1 * t273 * t36 * t146 * t142 - 0.2e1 / 0.3e1 * t141 * t146 * t79 - 0.2e1 / 0.3e1 * t328 * t146 * t149 + t620 - 0.28e2 / 0.81e2 * t533 * t92 * t290 + 0.4e1 / 0.9e1 * t338 * t82 * t157 - 0.8e1 / 0.9e1 * t289 * t36 * t146 * t154 + 0.2e1 / 0.3e1 * t153 * t146 * t82 + 0.2e1 / 0.3e1 * t341 * t146 * t157 - t638 - t171 * t284 / 0.9e1 - 0.4e1 * t75 * t281 - t177 * t296 / 0.9e1 + 0.4e1 * t81 * t281;
644 t650 = -0.36e2 * t316 * t183 - 0.6e1 * t74 * t647 + t470 - t479 - t485 - 0.36e2 * t489 - 0.18e2 * t500 - 0.18e2 * t502 - 0.6e1 * t506 - 0.36e2 * t587 - 0.6e1 * t589 - t592 - 0.18e2 * t593;
645 t654 = 0.72e2 * t312;
646 t660 = 0.36e2 * t314;
647 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
648 v4rho4[1] = t435 + t9 * (t581 + t650) - t553 + t554 - 0.72e2 * t261 - t270 - t654 - 0.36e2 * t317 - 0.72e2 * t320 - 0.36e2 * t323 + t557 - 0.54e2 * t258 - 0.36e2 * t263 - t267 - t301 - t404 - t660 - 0.18e2 * t352;
649
650 t671 = t36 * t36;
651 t674 = 0.1e1 / t280 / t76;
652 t713 = -0.28e2 / 0.81e2 * t512 * t193 * t142 + 0.32e2 / 0.27e2 * t325 * t79 * t36 * t146 + 0.4e1 / 0.27e2 * t372 * t149 - 0.8e1 / 0.9e1 * t141 * t671 * t674 - 0.4e1 / 0.9e1 * t171 * t146 + 0.4e1 / 0.3e1 * t171 * t282 + 0.4e1 / 0.27e2 * t273 * t197 * t142 - 0.2e1 / 0.9e1 * t141 * t382 * t79 - t377 * t149 / 0.9e1 + t620 - 0.28e2 / 0.81e2 * t533 * t200 * t154 - 0.32e2 / 0.27e2 * t338 * t82 * t36 * t146 + 0.4e1 / 0.27e2 * t385 * t157 - 0.8e1 / 0.9e1 * t153 * t671 * t674 + 0.4e1 / 0.9e1 * t177 * t146 - 0.4e1 / 0.3e1 * t177 * t282 + 0.4e1 / 0.27e2 * t289 * t203 * t154 - 0.2e1 / 0.9e1 * t153 * t393 * t82 - t390 * t157 / 0.9e1 - t638;
653 t716 = t260 * t364;
654 t723 = t260 * t369;
655 t734 = -0.12e2 * t136 * t206 * t160 - 0.24e2 * t136 * t95 * t351 - 0.24e2 * t136 * t396 * t85 - 0.48e2 * t361 * t364 - 0.6e1 * t74 * t713 + t470 - t472 + t475 - t479 + t487 - 0.48e2 * t564 - 0.48e2 * t566 - 0.24e2 * t568 - 0.48e2 * t716 - 0.24e2 * t723;
656 t742 = t183 * t183;
657 t747 = t488 * t190;
658 t750 = t72 * t85 * t190;
659 t755 = t257 * t206;
660 t757 = t133 * t396;
661 t759 = -0.12e2 * t35 * t160 * t190 - 0.24e2 * t136 * t742 - 0.12e2 * t316 * t206 - t485 - 0.12e2 * t489 - 0.12e2 * t500 - 0.6e1 * t502 - 0.24e2 * t587 - 0.12e2 * t589 - 0.24e2 * t591 - 0.12e2 * t593 - 0.12e2 * t747 - 0.24e2 * t750 - 0.6e1 * t755 - 0.12e2 * t757;
662 t779 = -0.24e2 * t370 + t557 - 0.36e2 * t258 - 0.12e2 * t263 - 0.36e2 * t310 - 0.48e2 * t314 - 0.12e2 * t352 - 0.24e2 * t359 - 0.24e2 * t362 - 0.12e2 * t367 - 0.12e2 * t397;
663 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
664 v4rho4[2] = t435 + t9 * (t734 + t759) - t553 + t554 - 0.24e2 * t261 - 0.96e2 * t312 - 0.24e2 * t317 - 0.48e2 * t320 - 0.24e2 * t323 - 0.48e2 * t365 + t779;
665
666 t780 = t72 * t407;
667 t782 = t133 * t430;
668 t802 = 0.12e2 * t281 + 0.24e2 * t526;
669 t824 = -0.28e2 / 0.81e2 * t512 * t413 * t79 + 0.8e1 / 0.9e1 * t372 * t147 + 0.4e1 / 0.9e1 * t325 * t197 * t79 - 0.2e1 / 0.3e1 * t328 * t146 * t197 - t171 * t382 / 0.3e1 - t141 * t419 * t79 / 0.9e1 + t75 * t802 / 0.3e1 - 0.28e2 / 0.81e2 * t533 * t422 * t82 - 0.8e1 / 0.9e1 * t385 * t147 + 0.4e1 / 0.9e1 * t338 * t203 * t82 + 0.2e1 / 0.3e1 * t341 * t146 * t203 - t177 * t393 / 0.3e1 - t153 * t427 * t82 / 0.9e1 - t81 * t802 / 0.3e1;
670 t827 = t35 * t190;
671 t830 = t260 * t410;
672 t845 = -0.36e2 * t136 * t183 * t206 - 0.36e2 * t136 * t95 * t396 - 0.12e2 * t136 * t430 * t85 - 0.36e2 * t827 * t183 - 0.36e2 * t361 * t410 - 0.6e1 * t74 * t824 + t487 - 0.72e2 * t716 - 0.36e2 * t723 - 0.12e2 * t780 - 0.6e1 * t782 - 0.36e2 * t830;
673 t852 = -t472 + t475 - t565 + t470 - t479 - t485 - 0.6e1 * t500 - 0.18e2 * t589 - t592 - 0.36e2 * t747 - 0.36e2 * t750 - 0.18e2 * t755 - 0.18e2 * t757;
674 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
675 v4rho4[3] = t435 + t9 * (t845 + t852) - t553 + t554 - t654 - 0.72e2 * t365 - 0.36e2 * t370 - t412 + t557 - t259 - 0.54e2 * t310 - t660 - 0.72e2 * t359 - 0.36e2 * t362 - 0.36e2 * t367 - 0.18e2 * t397 - t409 - t432;
676
677 t864 = t206 * t206;
678 t876 = t193 * t193;
679 t881 = t197 * t197;
680 t887 = 0.24e2 * t281 + 0.24e2 * t526;
681 t890 = t200 * t200;
682 t895 = t203 * t203;
683 t906 = t487 - 0.48e2 * t780 + t475 - t472 - t479 - t485 - 0.144e3 * t830 - 0.36e2 * t136 * t864 - 0.48e2 * t136 * t95 * t430 + t470 - 0.24e2 * t589 - 0.72e2 * t747 - 0.36e2 * t755 - 0.24e2 * t782 - 0.72e2 * t827 * t206 - 0.6e1 * t74 * (-0.28e2 / 0.81e2 * t512 * t876 + 0.8e1 / 0.9e1 * t372 * t197 - t141 * t881 / 0.3e1 - 0.4e1 / 0.9e1 * t171 * t419 + t75 * t887 / 0.3e1 - 0.28e2 / 0.81e2 * t533 * t890 + 0.8e1 / 0.9e1 * t385 * t203 - t153 * t895 / 0.3e1 - 0.4e1 / 0.9e1 * t177 * t427 - t81 * t887 / 0.3e1);
684 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
685 v4rho4[4] = t9 * t906 - 0.72e2 * t310 - 0.144e3 * t359 - 0.72e2 * t367 - 0.48e2 * t408 - 0.144e3 * t411 - 0.24e2 * t431 + t435 - t553 + t554 + t557;
686
687 #ifndef XC_DONT_COMPILE_MXC
688
689 if(order < 5) return;
690
691
692 #endif
693
694 #endif
695
696 #endif
697
698 #endif
699
700 #endif
701
702
703 }
704
705