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_1d_csc.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, t8, t9, t13;
24 double t14, t15, t16, t17, t18, t21, t24, t25;
25 double t26, t29, t30, t31, t33, t34;
26
27 #ifndef XC_DONT_COMPILE_VXC
28 double t37, t40, t41, t42, t44, t49, t50, t52;
29 double t54, t55, t56, t64, t65;
30
31 #ifndef XC_DONT_COMPILE_FXC
32 double t67, t70, t71, t74, t77, t78, t81, t82;
33 double t89, t94, t96, t98, t99, t100, t102, t104;
34 double t105, t106, t110, t111, t112, t113, t117, t122;
35 double t128, t129;
36
37 #ifndef XC_DONT_COMPILE_KXC
38 double t131, t134, t135, t138, t143, t146, t149, t154;
39 double t157, t158, t160, t163, t164, t169, t170, t171;
40 double t174, t177, t178, t179, t180, t183, t184, t191;
41 double t193, t195, t198, t199, t201, t202, t217, t226;
42 double t228, t232, t242, t252, t253, t255;
43
44 #ifndef XC_DONT_COMPILE_LXC
45 double t261, t274, t290, t297, t298, t311, t317, t329;
46 double t332, t377, t390, t423, t454, t470;
47 #endif
48
49 #endif
50
51 #endif
52
53 #endif
54
55 #endif
56
57
58 lda_c_1d_csc_params *params;
59
60 assert(p->params != NULL);
61 params = (lda_c_1d_csc_params * )(p->params);
62
63 t1 = 0.1e1 / rho[0];
64 t2 = t1 / 0.2e1;
65 t3 = params->para[4];
66 t4 = rho[0] * rho[0];
67 t5 = 0.1e1 / t4;
68 t8 = t2 + t3 * t5 / 0.4e1;
69 t9 = params->para[7];
70 t13 = params->para[9];
71 t14 = pow(t2, t13);
72 t15 = params->para[8] * t14;
73 t16 = 0.1e1 + t9 * t1 / 0.2e1 + t15;
74 t17 = log(t16);
75 t18 = t8 * t17;
76 t21 = params->para[1];
77 t24 = params->para[5];
78 t25 = pow(t2, t24);
79 t26 = params->para[2] * t25;
80 t29 = params->para[6];
81 t30 = pow(t2, t29);
82 t31 = params->para[3] * t30;
83 t33 = t21 * t1 + 0.2e1 * t26 + 0.2e1 * t31 + 0.2e1 * params->para[0];
84 t34 = 0.1e1 / t33;
85 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
86 zk[0] = -t18 * t34;
87
88 #ifndef XC_DONT_COMPILE_VXC
89
90 if(order < 1) return;
91
92
93 t37 = 0.1e1 / t4 / rho[0];
94 t40 = -t3 * t37 / 0.2e1 - t5 / 0.2e1;
95 t41 = rho[0] * t40;
96 t42 = t17 * t34;
97 t44 = rho[0] * t8;
98 t49 = -t9 * t5 / 0.2e1 - t15 * t13 * t1;
99 t50 = 0.1e1 / t16;
100 t52 = t49 * t50 * t34;
101 t54 = t33 * t33;
102 t55 = 0.1e1 / t54;
103 t56 = t17 * t55;
104 t64 = -0.2e1 * t26 * t24 * t1 - 0.2e1 * t31 * t29 * t1 - t21 * t5;
105 t65 = t56 * t64;
106 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
107 vrho[0] = -t41 * t42 - t44 * t52 + t44 * t65 + (-t18 * t34);
108
109 #ifndef XC_DONT_COMPILE_FXC
110
111 if(order < 2) return;
112
113
114 t67 = t40 * t17;
115 t70 = t8 * t49;
116 t71 = t50 * t34;
117 t74 = t55 * t64;
118 t77 = t4 * t4;
119 t78 = 0.1e1 / t77;
120 t81 = t37 + 0.3e1 / 0.2e1 * t3 * t78;
121 t82 = rho[0] * t81;
122 t89 = t13 * t13;
123 t94 = t15 * t13 * t5 + t15 * t89 * t5 + t9 * t37;
124 t96 = t94 * t50 * t34;
125 t98 = t49 * t49;
126 t99 = t16 * t16;
127 t100 = 0.1e1 / t99;
128 t102 = t98 * t100 * t34;
129 t104 = t44 * t49;
130 t105 = t50 * t55;
131 t106 = t105 * t64;
132 t110 = 0.1e1 / t54 / t33;
133 t111 = t17 * t110;
134 t112 = t64 * t64;
135 t113 = t111 * t112;
136 t117 = t24 * t24;
137 t122 = t29 * t29;
138 t128 = 0.2e1 * t26 * t117 * t5 + 0.2e1 * t31 * t122 * t5 + 0.2e1 * t26 * t24 * t5 + 0.2e1 * t31 * t29 * t5 + 0.2e1 * t21 * t37;
139 t129 = t56 * t128;
140 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
141 v2rho2[0] = t44 * t102 + 0.2e1 * t104 * t106 - 0.2e1 * t44 * t113 + t44 * t129 + 0.2e1 * t18 * t74 - 0.2e1 * t67 * t34 - 0.2e1 * t41 * t52 + 0.2e1 * t41 * t65 - t82 * t42 - t44 * t96 - 0.2e1 * t70 * t71;
142
143 #ifndef XC_DONT_COMPILE_KXC
144
145 if(order < 3) return;
146
147
148 t131 = t81 * t17;
149 t134 = t8 * t98;
150 t135 = t100 * t34;
151 t138 = t110 * t112;
152 t143 = t40 * t49;
153 t146 = t8 * t94;
154 t149 = t55 * t128;
155 t154 = 0.1e1 / t77 / rho[0];
156 t157 = -0.6e1 * t3 * t154 - 0.3e1 * t78;
157 t158 = rho[0] * t157;
158 t160 = t41 * t49;
159 t163 = t44 * t94;
160 t164 = t135 * t49;
161 t169 = t44 * t98;
162 t170 = t100 * t55;
163 t171 = t170 * t64;
164 t174 = t105 * t128;
165 t177 = 0.3e1 * t104 * t174 + 0.6e1 * t160 * t106 + 0.3e1 * t163 * t106 - 0.3e1 * t131 * t34 + 0.3e1 * t134 * t135 - 0.6e1 * t18 * t138 - 0.6e1 * t143 * t71 - 0.3e1 * t146 * t71 + 0.3e1 * t18 * t149 - t158 * t42 + 0.3e1 * t163 * t164 - 0.3e1 * t169 * t171 + 0.6e1 * t67 * t74;
166 t178 = t44 * t17;
167 t179 = t110 * t64;
168 t180 = t179 * t128;
169 t183 = t50 * t110;
170 t184 = t183 * t112;
171 t191 = t98 * t49;
172 t193 = 0.1e1 / t99 / t16;
173 t195 = t191 * t193 * t34;
174 t198 = t54 * t54;
175 t199 = 0.1e1 / t198;
176 t201 = t112 * t64;
177 t202 = t17 * t199 * t201;
178 t217 = t89 * t13;
179 t226 = -0.2e1 * t15 * t13 * t37 - t15 * t217 * t37 - 0.3e1 * t15 * t89 * t37 - 0.3e1 * t9 * t78;
180 t228 = t226 * t50 * t34;
181 t232 = t117 * t24;
182 t242 = t122 * t29;
183 t252 = -0.6e1 * t26 * t117 * t37 - 0.6e1 * t31 * t122 * t37 - 0.2e1 * t26 * t232 * t37 - 0.4e1 * t26 * t24 * t37 - 0.2e1 * t31 * t242 * t37 - 0.4e1 * t31 * t29 * t37 - 0.6e1 * t21 * t78;
184 t253 = t56 * t252;
185 t255 = 0.3e1 * t41 * t102 - 0.6e1 * t104 * t184 + 0.6e1 * t70 * t106 - 0.6e1 * t41 * t113 + 0.3e1 * t41 * t129 - 0.6e1 * t178 * t180 - 0.2e1 * t44 * t195 + 0.6e1 * t44 * t202 - t44 * t228 + t44 * t253 - 0.3e1 * t41 * t96 - 0.3e1 * t82 * t52 + 0.3e1 * t82 * t65;
186 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
187 v3rho3[0] = t177 + t255;
188
189 #ifndef XC_DONT_COMPILE_LXC
190
191 if(order < 4) return;
192
193
194 t261 = t117 * t117;
195 t274 = t122 * t122;
196 t290 = t94 * t94;
197 t297 = t98 * t98;
198 t298 = t99 * t99;
199 t311 = t112 * t112;
200 t317 = t128 * t128;
201 t329 = t193 * t34;
202 t332 = -0.4e1 * t157 * t17 * t34 + t44 * t56 * (0.22e2 * t26 * t117 * t78 + 0.22e2 * t31 * t122 * t78 + 0.12e2 * t26 * t232 * t78 + 0.12e2 * t26 * t24 * t78 + 0.12e2 * t31 * t242 * t78 + 0.2e1 * t26 * t261 * t78 + 0.2e1 * t31 * t274 * t78 + 0.12e2 * t31 * t29 * t78 + 0.24e2 * t21 * t154) + 0.3e1 * t44 * t290 * t100 * t34 + 0.12e2 * t146 * t106 + 0.6e1 * t44 * t297 / t298 * t34 - 0.12e2 * t134 * t171 + 0.12e2 * t70 * t174 - 0.24e2 * t44 * t17 / t198 / t33 * t311 - 0.24e2 * t18 * t180 - 0.6e1 * t44 * t111 * t317 - 0.24e2 * t70 * t184 + 0.12e2 * t40 * t98 * t135 - 0.24e2 * t67 * t138 - 0.8e1 * t8 * t191 * t329;
203 t377 = 0.24e2 * t18 * t199 * t201 + 0.12e2 * t131 * t74 - 0.12e2 * t81 * t49 * t71 - 0.12e2 * t40 * t94 * t71 + 0.12e2 * t67 * t149 - 0.4e1 * t8 * t226 * t71 + 0.4e1 * t18 * t55 * t252 - rho[0] * (0.12e2 * t154 + 0.30e2 * t3 / t77 / t4) * t42 + 0.12e2 * t160 * t135 * t94 - 0.12e2 * t41 * t98 * t171 + 0.12e2 * t82 * t49 * t106 + 0.12e2 * t41 * t94 * t106 + 0.12e2 * t160 * t174 - 0.24e2 * t41 * t17 * t180;
204 t390 = t44 * t226;
205 t423 = 0.8e1 * t44 * t191 * t193 * t55 * t64 + 0.12e2 * t169 * t100 * t110 * t112 + 0.24e2 * t104 * t50 * t199 * t201 + 0.36e2 * t178 * t199 * t112 * t128 - 0.12e2 * t163 * t170 * t49 * t64 + 0.4e1 * t104 * t105 * t252 - 0.6e1 * t169 * t170 * t128 - 0.12e2 * t163 * t329 * t98 - 0.8e1 * t178 * t179 * t252 + 0.4e1 * t390 * t106 - 0.24e2 * t160 * t184 + 0.6e1 * t163 * t174 - 0.12e2 * t163 * t184 + 0.4e1 * t390 * t164;
206 t454 = t89 * t89;
207 t470 = -0.24e2 * t104 * t183 * t128 * t64 - 0.8e1 * t41 * t195 + 0.24e2 * t143 * t106 + 0.24e2 * t41 * t202 + 0.12e2 * t146 * t164 + 0.6e1 * t82 * t102 - 0.12e2 * t82 * t113 + 0.4e1 * t158 * t65 + 0.6e1 * t82 * t129 - 0.4e1 * t158 * t52 - 0.6e1 * t82 * t96 - 0.4e1 * t41 * t228 + 0.4e1 * t41 * t253 - t44 * (0.6e1 * t15 * t13 * t78 + 0.6e1 * t15 * t217 * t78 + t15 * t454 * t78 + 0.11e2 * t15 * t89 * t78 + 0.12e2 * t9 * t154) * t50 * t34;
208 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
209 v4rho4[0] = t332 + t377 + t423 + t470;
210
211 #ifndef XC_DONT_COMPILE_MXC
212
213 if(order < 5) return;
214
215
216 #endif
217
218 #endif
219
220 #endif
221
222 #endif
223
224 #endif
225
226
227 }
228
229
230 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))231 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
232 {
233
234 #ifndef XC_DONT_COMPILE_EXC
235 double t1, t2, t3, t4, t5, t8, t9, t13;
236 double t14, t15, t16, t17, t18, t21, t24, t25;
237 double t26, t29, t30, t31, t33, t34;
238
239 #ifndef XC_DONT_COMPILE_VXC
240 double t37, t40, t41, t42, t44, t49, t50, t52;
241 double t54, t55, t56, t64, t65;
242
243 #ifndef XC_DONT_COMPILE_FXC
244 double t67, t70, t71, t74, t77, t78, t81, t82;
245 double t89, t94, t96, t98, t99, t100, t102, t104;
246 double t105, t106, t110, t111, t112, t113, t117, t122;
247 double t128, t129;
248
249 #ifndef XC_DONT_COMPILE_KXC
250 double t131, t134, t135, t138, t143, t146, t149, t154;
251 double t157, t158, t160, t163, t164, t169, t170, t171;
252 double t174, t177, t178, t179, t180, t183, t184, t191;
253 double t193, t195, t198, t199, t201, t202, t217, t226;
254 double t228, t232, t242, t252, t253, t255;
255
256 #ifndef XC_DONT_COMPILE_LXC
257 double t263, t269, t275, t288, t304, t311, t312, t329;
258 double t332, t377, t390, t423, t454, t470;
259 #endif
260
261 #endif
262
263 #endif
264
265 #endif
266
267 #endif
268
269
270 lda_c_1d_csc_params *params;
271
272 assert(p->params != NULL);
273 params = (lda_c_1d_csc_params * )(p->params);
274
275 t1 = 0.1e1 / rho[0];
276 t2 = t1 / 0.2e1;
277 t3 = params->ferro[4];
278 t4 = rho[0] * rho[0];
279 t5 = 0.1e1 / t4;
280 t8 = t2 + t3 * t5 / 0.4e1;
281 t9 = params->ferro[7];
282 t13 = params->ferro[9];
283 t14 = pow(t2, t13);
284 t15 = params->ferro[8] * t14;
285 t16 = 0.1e1 + t9 * t1 / 0.2e1 + t15;
286 t17 = log(t16);
287 t18 = t8 * t17;
288 t21 = params->ferro[1];
289 t24 = params->ferro[5];
290 t25 = pow(t2, t24);
291 t26 = params->ferro[2] * t25;
292 t29 = params->ferro[6];
293 t30 = pow(t2, t29);
294 t31 = params->ferro[3] * t30;
295 t33 = t21 * t1 + 0.2e1 * t26 + 0.2e1 * t31 + 0.2e1 * params->ferro[0];
296 t34 = 0.1e1 / t33;
297 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
298 zk[0] = -t18 * t34;
299
300 #ifndef XC_DONT_COMPILE_VXC
301
302 if(order < 1) return;
303
304
305 t37 = 0.1e1 / t4 / rho[0];
306 t40 = -t3 * t37 / 0.2e1 - t5 / 0.2e1;
307 t41 = rho[0] * t40;
308 t42 = t17 * t34;
309 t44 = rho[0] * t8;
310 t49 = -t9 * t5 / 0.2e1 - t15 * t13 * t1;
311 t50 = 0.1e1 / t16;
312 t52 = t49 * t50 * t34;
313 t54 = t33 * t33;
314 t55 = 0.1e1 / t54;
315 t56 = t17 * t55;
316 t64 = -0.2e1 * t26 * t24 * t1 - 0.2e1 * t31 * t29 * t1 - t21 * t5;
317 t65 = t56 * t64;
318 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
319 vrho[0] = -t41 * t42 - t44 * t52 + t44 * t65 + (-t18 * t34);
320
321 #ifndef XC_DONT_COMPILE_FXC
322
323 if(order < 2) return;
324
325
326 t67 = t40 * t17;
327 t70 = t8 * t49;
328 t71 = t50 * t34;
329 t74 = t55 * t64;
330 t77 = t4 * t4;
331 t78 = 0.1e1 / t77;
332 t81 = t37 + 0.3e1 / 0.2e1 * t3 * t78;
333 t82 = rho[0] * t81;
334 t89 = t13 * t13;
335 t94 = t15 * t13 * t5 + t15 * t89 * t5 + t9 * t37;
336 t96 = t94 * t50 * t34;
337 t98 = t49 * t49;
338 t99 = t16 * t16;
339 t100 = 0.1e1 / t99;
340 t102 = t98 * t100 * t34;
341 t104 = t44 * t49;
342 t105 = t50 * t55;
343 t106 = t105 * t64;
344 t110 = 0.1e1 / t54 / t33;
345 t111 = t17 * t110;
346 t112 = t64 * t64;
347 t113 = t111 * t112;
348 t117 = t24 * t24;
349 t122 = t29 * t29;
350 t128 = 0.2e1 * t26 * t117 * t5 + 0.2e1 * t31 * t122 * t5 + 0.2e1 * t26 * t24 * t5 + 0.2e1 * t31 * t29 * t5 + 0.2e1 * t21 * t37;
351 t129 = t56 * t128;
352 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
353 v2rho2[0] = t44 * t102 + 0.2e1 * t104 * t106 - 0.2e1 * t44 * t113 + t44 * t129 + 0.2e1 * t18 * t74 - 0.2e1 * t67 * t34 - 0.2e1 * t41 * t52 + 0.2e1 * t41 * t65 - t82 * t42 - t44 * t96 - 0.2e1 * t70 * t71;
354
355 #ifndef XC_DONT_COMPILE_KXC
356
357 if(order < 3) return;
358
359
360 t131 = t81 * t17;
361 t134 = t8 * t98;
362 t135 = t100 * t34;
363 t138 = t110 * t112;
364 t143 = t40 * t49;
365 t146 = t8 * t94;
366 t149 = t55 * t128;
367 t154 = 0.1e1 / t77 / rho[0];
368 t157 = -0.6e1 * t3 * t154 - 0.3e1 * t78;
369 t158 = rho[0] * t157;
370 t160 = t41 * t49;
371 t163 = t44 * t94;
372 t164 = t135 * t49;
373 t169 = t44 * t98;
374 t170 = t100 * t55;
375 t171 = t170 * t64;
376 t174 = t105 * t128;
377 t177 = 0.3e1 * t104 * t174 + 0.6e1 * t160 * t106 + 0.3e1 * t163 * t106 - 0.3e1 * t131 * t34 + 0.3e1 * t134 * t135 - 0.6e1 * t18 * t138 - 0.6e1 * t143 * t71 - 0.3e1 * t146 * t71 + 0.3e1 * t18 * t149 - t158 * t42 + 0.3e1 * t163 * t164 - 0.3e1 * t169 * t171 + 0.6e1 * t67 * t74;
378 t178 = t44 * t17;
379 t179 = t110 * t64;
380 t180 = t179 * t128;
381 t183 = t50 * t110;
382 t184 = t183 * t112;
383 t191 = t98 * t49;
384 t193 = 0.1e1 / t99 / t16;
385 t195 = t191 * t193 * t34;
386 t198 = t54 * t54;
387 t199 = 0.1e1 / t198;
388 t201 = t112 * t64;
389 t202 = t17 * t199 * t201;
390 t217 = t89 * t13;
391 t226 = -0.2e1 * t15 * t13 * t37 - t15 * t217 * t37 - 0.3e1 * t15 * t89 * t37 - 0.3e1 * t9 * t78;
392 t228 = t226 * t50 * t34;
393 t232 = t117 * t24;
394 t242 = t122 * t29;
395 t252 = -0.6e1 * t26 * t117 * t37 - 0.6e1 * t31 * t122 * t37 - 0.2e1 * t26 * t232 * t37 - 0.4e1 * t26 * t24 * t37 - 0.2e1 * t31 * t242 * t37 - 0.4e1 * t31 * t29 * t37 - 0.6e1 * t21 * t78;
396 t253 = t56 * t252;
397 t255 = 0.3e1 * t41 * t102 - 0.6e1 * t104 * t184 + 0.6e1 * t70 * t106 - 0.6e1 * t41 * t113 + 0.3e1 * t41 * t129 - 0.6e1 * t178 * t180 - 0.2e1 * t44 * t195 + 0.6e1 * t44 * t202 - t44 * t228 + t44 * t253 - 0.3e1 * t41 * t96 - 0.3e1 * t82 * t52 + 0.3e1 * t82 * t65;
398 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
399 v3rho3[0] = t177 + t255;
400
401 #ifndef XC_DONT_COMPILE_LXC
402
403 if(order < 4) return;
404
405
406 t263 = t112 * t112;
407 t269 = t128 * t128;
408 t275 = t117 * t117;
409 t288 = t122 * t122;
410 t304 = t94 * t94;
411 t311 = t98 * t98;
412 t312 = t99 * t99;
413 t329 = t193 * t34;
414 t332 = -0.24e2 * t70 * t184 + 0.12e2 * t70 * t174 - 0.24e2 * t44 * t17 / t198 / t33 * t263 - 0.24e2 * t18 * t180 - 0.6e1 * t44 * t111 * t269 + t44 * t56 * (0.22e2 * t26 * t117 * t78 + 0.22e2 * t31 * t122 * t78 + 0.12e2 * t26 * t232 * t78 + 0.12e2 * t26 * t24 * t78 + 0.12e2 * t31 * t242 * t78 + 0.2e1 * t26 * t275 * t78 + 0.2e1 * t31 * t288 * t78 + 0.12e2 * t31 * t29 * t78 + 0.24e2 * t21 * t154) + 0.3e1 * t44 * t304 * t100 * t34 + 0.12e2 * t146 * t106 + 0.6e1 * t44 * t311 / t312 * t34 - 0.12e2 * t134 * t171 - 0.4e1 * t157 * t17 * t34 + 0.12e2 * t40 * t98 * t135 - 0.24e2 * t67 * t138 - 0.8e1 * t8 * t191 * t329;
415 t377 = 0.24e2 * t18 * t199 * t201 + 0.12e2 * t131 * t74 - 0.12e2 * t81 * t49 * t71 - 0.12e2 * t40 * t94 * t71 + 0.12e2 * t67 * t149 - 0.4e1 * t8 * t226 * t71 + 0.4e1 * t18 * t55 * t252 - rho[0] * (0.12e2 * t154 + 0.30e2 * t3 / t77 / t4) * t42 + 0.12e2 * t160 * t135 * t94 - 0.12e2 * t41 * t98 * t171 + 0.12e2 * t82 * t49 * t106 + 0.12e2 * t41 * t94 * t106 + 0.12e2 * t160 * t174 - 0.24e2 * t41 * t17 * t180;
416 t390 = t44 * t226;
417 t423 = 0.8e1 * t44 * t191 * t193 * t55 * t64 + 0.12e2 * t169 * t100 * t110 * t112 + 0.24e2 * t104 * t50 * t199 * t201 + 0.36e2 * t178 * t199 * t112 * t128 - 0.12e2 * t163 * t170 * t49 * t64 + 0.4e1 * t104 * t105 * t252 - 0.6e1 * t169 * t170 * t128 - 0.12e2 * t163 * t329 * t98 - 0.8e1 * t178 * t179 * t252 + 0.4e1 * t390 * t106 - 0.24e2 * t160 * t184 + 0.6e1 * t163 * t174 - 0.12e2 * t163 * t184 + 0.4e1 * t390 * t164;
418 t454 = t89 * t89;
419 t470 = -0.24e2 * t104 * t183 * t128 * t64 - 0.8e1 * t41 * t195 + 0.24e2 * t143 * t106 + 0.24e2 * t41 * t202 + 0.12e2 * t146 * t164 + 0.6e1 * t82 * t102 - 0.12e2 * t82 * t113 + 0.4e1 * t158 * t65 + 0.6e1 * t82 * t129 - 0.4e1 * t158 * t52 - 0.6e1 * t82 * t96 - 0.4e1 * t41 * t228 + 0.4e1 * t41 * t253 - t44 * (0.6e1 * t15 * t13 * t78 + 0.6e1 * t15 * t217 * t78 + t15 * t454 * t78 + 0.11e2 * t15 * t89 * t78 + 0.12e2 * t9 * t154) * t50 * t34;
420 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
421 v4rho4[0] = t332 + t377 + t423 + t470;
422
423 #ifndef XC_DONT_COMPILE_MXC
424
425 if(order < 5) return;
426
427
428 #endif
429
430 #endif
431
432 #endif
433
434 #endif
435
436 #endif
437
438
439 }
440
441
442 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))443 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
444 {
445
446 #ifndef XC_DONT_COMPILE_EXC
447 double t1, t2, t3, t4, t5, t6, t9, t10;
448 double t14, t15, t16, t17, t18, t19, t22, t25;
449 double t26, t27, t30, t31, t32, t34, t35, t36;
450 double t37, t40, t41, t45, t46, t47, t48, t49;
451 double t50, t53, t56, t57, t58, t61, t62, t63;
452 double t65, t66, t68, t69, t70, t71, t72;
453
454 #ifndef XC_DONT_COMPILE_VXC
455 double t74, t77, t78, t79, t84, t85, t86, t87;
456 double t88, t89, t90, t98, t99, t100, t103, t104;
457 double t110, t111, t112, t113, t115, t116, t124, t125;
458 double t127, t128, t129, t130, t131, t132, t133, t134;
459
460 #ifndef XC_DONT_COMPILE_FXC
461 double t139, t140, t141, t142, t143, t144, t145, t146;
462 double t149, t150, t151, t152, t153, t154, t155, t156;
463 double t158, t163, t164, t165, t166, t167, t168, t169;
464 double t170, t171, t172, t173, t174, t175, t177, t178;
465 double t179, t180, t181, t183, t188, t194, t195, t196;
466 double t199, t200, t202, t208, t213, t214, t216, t217;
467 double t218, t219, t220, t222, t223, t227, t228, t229;
468 double t233, t238, t244, t245, t247, t248, t249, t250;
469 double t251, t252, t253, t254, t255, t256, t257, t258;
470 double t259, t260, t261, t263, t265;
471
472 #ifndef XC_DONT_COMPILE_KXC
473 double t267, t268, t269, t270, t271, t272, t273, t274;
474 double t276, t277, t278, t279, t282, t283, t284, t287;
475 double t290, t297, t306, t307, t311, t321, t331, t332;
476 double t334, t339, t340, t342, t343, t346, t347, t348;
477 double t349, t352, t357, t358, t361, t364, t365, t371;
478 double t376, t377, t379, t380, t381, t382, t383, t384;
479 double t385, t386, t388, t389, t390, t391, t392, t393;
480 double t394, t395, t396, t397, t398, t399, t400, t401;
481 double t402, t403, t404, t405, t406, t407, t408, t410;
482 double t411, t412, t413, t414, t415, t416, t417, t418;
483 double t419, t420, t421, t422, t423, t424, t425, t426;
484 double t427, t428, t429, t432, t441, t442, t443, t446;
485 double t456, t466, t467, t468, t469, t470, t471, t472;
486 double t473, t474, t475, t476, t477, t478, t479, t480;
487 double t481, t484, t485, t486, t487, t488, t489, t490;
488 double t491, t492, t493, t494, t495, t496, t497, t498;
489 double t499, t500, t501, t502, t503, t506, t509;
490
491 #ifndef XC_DONT_COMPILE_LXC
492 double t512, t514, t516, t521, t524, t526, t528, t530;
493 double t532, t536, t539, t542, t543, t545, t546, t548;
494 double t553, t580, t597, t610, t635, t639, t645, t649;
495 double t651, t659, t705, t709, t712, t715, t717, t720;
496 double t734, t737, t750, t765, t768, t771, t772, t775;
497 double t778, t780, t782, t783, t785, t787, t788, t791;
498 double t792, t794, t798, t801, t804, t806, t808, t811;
499 double t813, t815, t817, t819, t823, t827, t830, t835;
500 double t844, t848, t850, t851, t853, t854, t855, t857;
501 double t858, t859, t861, t862, t863, t864, t866, t867;
502 double t868, t871, t872, t873, t874, t875, t876, t877;
503 double t878, t879, t880, t881, t882, t883, t884, t885;
504 double t886, t887, t888, t889, t890, t891, t892, t893;
505 double t894, t895, t896, t897, t898, t899, t900, t901;
506 double t902, t903, t904, t905, t907, t908, t909, t912;
507 double t913, t914, t915, t916, t922, t923, t925, t929;
508 double t930, t934, t935, t939;
509 #endif
510
511 #endif
512
513 #endif
514
515 #endif
516
517 #endif
518
519
520 lda_c_1d_csc_params *params;
521
522 assert(p->params != NULL);
523 params = (lda_c_1d_csc_params * )(p->params);
524
525 t1 = rho[0] + rho[1];
526 t2 = 0.1e1 / t1;
527 t3 = t2 / 0.2e1;
528 t4 = params->para[4];
529 t5 = t1 * t1;
530 t6 = 0.1e1 / t5;
531 t9 = t3 + t4 * t6 / 0.4e1;
532 t10 = params->para[7];
533 t14 = params->para[9];
534 t15 = pow(t3, t14);
535 t16 = params->para[8] * t15;
536 t17 = 0.1e1 + t10 * t2 / 0.2e1 + t16;
537 t18 = log(t17);
538 t19 = t9 * t18;
539 t22 = params->para[1];
540 t25 = params->para[5];
541 t26 = pow(t3, t25);
542 t27 = params->para[2] * t26;
543 t30 = params->para[6];
544 t31 = pow(t3, t30);
545 t32 = params->para[3] * t31;
546 t34 = t22 * t2 + 0.2e1 * t27 + 0.2e1 * t32 + 0.2e1 * params->para[0];
547 t35 = 0.1e1 / t34;
548 t36 = t19 * t35;
549 t37 = params->ferro[4];
550 t40 = t3 + t37 * t6 / 0.4e1;
551 t41 = params->ferro[7];
552 t45 = params->ferro[9];
553 t46 = pow(t3, t45);
554 t47 = params->ferro[8] * t46;
555 t48 = 0.1e1 + t41 * t2 / 0.2e1 + t47;
556 t49 = log(t48);
557 t50 = t40 * t49;
558 t53 = params->ferro[1];
559 t56 = params->ferro[5];
560 t57 = pow(t3, t56);
561 t58 = params->ferro[2] * t57;
562 t61 = params->ferro[6];
563 t62 = pow(t3, t61);
564 t63 = params->ferro[3] * t62;
565 t65 = t53 * t2 + 0.2e1 * t58 + 0.2e1 * t63 + 0.2e1 * params->ferro[0];
566 t66 = 0.1e1 / t65;
567 t68 = -t50 * t66 + t36;
568 t69 = rho[0] - rho[1];
569 t70 = t69 * t69;
570 t71 = t68 * t70;
571 t72 = t71 * t6;
572 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
573 zk[0] = -t36 + t72;
574
575 #ifndef XC_DONT_COMPILE_VXC
576
577 if(order < 1) return;
578
579
580 t74 = 0.1e1 / t5 / t1;
581 t77 = -t4 * t74 / 0.2e1 - t6 / 0.2e1;
582 t78 = t77 * t18;
583 t79 = t78 * t35;
584 t84 = -t10 * t6 / 0.2e1 - t16 * t14 * t2;
585 t85 = t9 * t84;
586 t86 = 0.1e1 / t17;
587 t87 = t86 * t35;
588 t88 = t85 * t87;
589 t89 = t34 * t34;
590 t90 = 0.1e1 / t89;
591 t98 = -0.2e1 * t27 * t25 * t2 - 0.2e1 * t32 * t30 * t2 - t22 * t6;
592 t99 = t90 * t98;
593 t100 = t19 * t99;
594 t103 = -t37 * t74 / 0.2e1 - t6 / 0.2e1;
595 t104 = t103 * t49;
596 t110 = -t41 * t6 / 0.2e1 - t47 * t45 * t2;
597 t111 = t40 * t110;
598 t112 = 0.1e1 / t48;
599 t113 = t112 * t66;
600 t115 = t65 * t65;
601 t116 = 0.1e1 / t115;
602 t124 = -0.2e1 * t58 * t56 * t2 - 0.2e1 * t63 * t61 * t2 - t53 * t6;
603 t125 = t116 * t124;
604 t127 = -t104 * t66 - t111 * t113 + t50 * t125 - t100 + t79 + t88;
605 t128 = t127 * t70;
606 t129 = t128 * t6;
607 t130 = t68 * t69;
608 t131 = t130 * t6;
609 t132 = 0.2e1 * t131;
610 t133 = t71 * t74;
611 t134 = 0.2e1 * t133;
612 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
613 vrho[0] = -t36 + t72 + t1 * (-t79 - t88 + t100 + t129 + t132 - t134);
614
615 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
616 vrho[1] = -t36 + t72 + t1 * (-t79 - t88 + t100 + t129 - t132 - t134);
617
618 #ifndef XC_DONT_COMPILE_FXC
619
620 if(order < 2) return;
621
622
623 t139 = 0.2e1 * t79;
624 t140 = 0.2e1 * t88;
625 t141 = 0.2e1 * t100;
626 t142 = 0.2e1 * t129;
627 t143 = 0.4e1 * t131;
628 t144 = 0.4e1 * t133;
629 t145 = t5 * t5;
630 t146 = 0.1e1 / t145;
631 t149 = t74 + 0.3e1 / 0.2e1 * t4 * t146;
632 t150 = t149 * t18;
633 t151 = t150 * t35;
634 t152 = t77 * t84;
635 t153 = t152 * t87;
636 t154 = 0.2e1 * t153;
637 t155 = t78 * t99;
638 t156 = 0.2e1 * t155;
639 t158 = t14 * t14;
640 t163 = t16 * t14 * t6 + t16 * t158 * t6 + t10 * t74;
641 t164 = t9 * t163;
642 t165 = t164 * t87;
643 t166 = t84 * t84;
644 t167 = t9 * t166;
645 t168 = t17 * t17;
646 t169 = 0.1e1 / t168;
647 t170 = t169 * t35;
648 t171 = t167 * t170;
649 t172 = t86 * t90;
650 t173 = t172 * t98;
651 t174 = t85 * t173;
652 t175 = 0.2e1 * t174;
653 t177 = 0.1e1 / t89 / t34;
654 t178 = t98 * t98;
655 t179 = t177 * t178;
656 t180 = t19 * t179;
657 t181 = 0.2e1 * t180;
658 t183 = t25 * t25;
659 t188 = t30 * t30;
660 t194 = 0.2e1 * t27 * t183 * t6 + 0.2e1 * t32 * t188 * t6 + 0.2e1 * t27 * t25 * t6 + 0.2e1 * t32 * t30 * t6 + 0.2e1 * t22 * t74;
661 t195 = t90 * t194;
662 t196 = t19 * t195;
663 t199 = t74 + 0.3e1 / 0.2e1 * t37 * t146;
664 t200 = t199 * t49;
665 t202 = t103 * t110;
666 t208 = t45 * t45;
667 t213 = t47 * t208 * t6 + t47 * t45 * t6 + t41 * t74;
668 t214 = t40 * t213;
669 t216 = t110 * t110;
670 t217 = t40 * t216;
671 t218 = t48 * t48;
672 t219 = 0.1e1 / t218;
673 t220 = t219 * t66;
674 t222 = t112 * t116;
675 t223 = t222 * t124;
676 t227 = 0.1e1 / t115 / t65;
677 t228 = t124 * t124;
678 t229 = t227 * t228;
679 t233 = t56 * t56;
680 t238 = t61 * t61;
681 t244 = 0.2e1 * t58 * t233 * t6 + 0.2e1 * t63 * t238 * t6 + 0.2e1 * t58 * t56 * t6 + 0.2e1 * t63 * t61 * t6 + 0.2e1 * t53 * t74;
682 t245 = t116 * t244;
683 t247 = 0.2e1 * t104 * t125 + 0.2e1 * t111 * t223 - 0.2e1 * t202 * t113 - t214 * t113 - t200 * t66 + t217 * t220 - 0.2e1 * t50 * t229 + t50 * t245 + t151 + t154 - t156 + t165 - t171 - t175 + t181 - t196;
684 t248 = t247 * t70;
685 t249 = t248 * t6;
686 t250 = t127 * t69;
687 t251 = t250 * t6;
688 t252 = 0.4e1 * t251;
689 t253 = t128 * t74;
690 t254 = 0.4e1 * t253;
691 t255 = t68 * t6;
692 t256 = 0.2e1 * t255;
693 t257 = t130 * t74;
694 t258 = 0.8e1 * t257;
695 t259 = t71 * t146;
696 t260 = 0.6e1 * t259;
697 t261 = -t151 - t154 + t156 - t165 + t171 + t175 - t181 + t196 + t249 + t252 - t254 + t256 - t258 + t260;
698 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
699 v2rho2[0] = t1 * t261 - t139 - t140 + t141 + t142 + t143 - t144;
700
701 t263 = -t151 - t154 + t156 - t165 + t171 + t175 - t181 + t196 + t249 - t254 - t256 + t260;
702 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
703 v2rho2[1] = t1 * t263 - t139 - t140 + t141 + t142 - t144;
704
705 t265 = -t151 - t154 + t156 - t165 + t171 + t175 - t181 + t196 + t249 - t252 - t254 + t256 + t258 + t260;
706 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
707 v2rho2[2] = t1 * t265 - t139 - t140 + t141 + t142 - t143 - t144;
708
709 #ifndef XC_DONT_COMPILE_KXC
710
711 if(order < 3) return;
712
713
714 t267 = t127 * t6;
715 t268 = 0.6e1 * t267;
716 t269 = t68 * t74;
717 t270 = 0.12e2 * t269;
718 t271 = t128 * t146;
719 t272 = 0.18e2 * t271;
720 t273 = t130 * t146;
721 t274 = 0.36e2 * t273;
722 t276 = 0.1e1 / t145 / t1;
723 t277 = t71 * t276;
724 t278 = 0.24e2 * t277;
725 t279 = 0.3e1 * t146;
726 t282 = -0.6e1 * t4 * t276 - t279;
727 t283 = t282 * t18;
728 t284 = t283 * t35;
729 t287 = t199 * t110;
730 t290 = t103 * t213;
731 t297 = t208 * t45;
732 t306 = -0.3e1 * t47 * t208 * t74 - t47 * t297 * t74 - 0.2e1 * t47 * t45 * t74 - 0.3e1 * t41 * t146;
733 t307 = t40 * t306;
734 t311 = t233 * t56;
735 t321 = t238 * t61;
736 t331 = -0.6e1 * t58 * t233 * t74 - 0.6e1 * t63 * t238 * t74 - 0.2e1 * t58 * t311 * t74 - 0.2e1 * t63 * t321 * t74 - 0.4e1 * t58 * t56 * t74 - 0.4e1 * t63 * t61 * t74 - 0.6e1 * t53 * t146;
737 t332 = t116 * t331;
738 t334 = t103 * t216;
739 t339 = t216 * t110;
740 t340 = t40 * t339;
741 t342 = 0.1e1 / t218 / t48;
742 t343 = t342 * t66;
743 t346 = t115 * t115;
744 t347 = 0.1e1 / t346;
745 t348 = t228 * t124;
746 t349 = t347 * t348;
747 t352 = t220 * t110;
748 t357 = t219 * t116;
749 t358 = t357 * t124;
750 t361 = t222 * t244;
751 t364 = t227 * t124;
752 t365 = t364 * t244;
753 t371 = t112 * t227 * t228;
754 t376 = -0.6e1 * t37 * t276 - t279;
755 t377 = t376 * t49;
756 t379 = -0.6e1 * t104 * t229 + 0.3e1 * t104 * t245 + 0.3e1 * t111 * t361 - 0.6e1 * t111 * t371 - 0.3e1 * t287 * t113 - 0.3e1 * t290 * t113 - t307 * t113 + 0.3e1 * t200 * t125 + 0.6e1 * t202 * t223 + 0.3e1 * t214 * t223 + 0.3e1 * t214 * t352 - 0.3e1 * t217 * t358 + 0.3e1 * t334 * t220 + t50 * t332 - 0.2e1 * t340 * t343 + 0.6e1 * t50 * t349 - 0.6e1 * t50 * t365 - t377 * t66;
757 t380 = t172 * t194;
758 t381 = t85 * t380;
759 t382 = 0.3e1 * t381;
760 t383 = t177 * t98;
761 t384 = t383 * t194;
762 t385 = t19 * t384;
763 t386 = 0.6e1 * t385;
764 t388 = t86 * t177 * t178;
765 t389 = t85 * t388;
766 t390 = 0.6e1 * t389;
767 t391 = t152 * t173;
768 t392 = 0.6e1 * t391;
769 t393 = t170 * t84;
770 t394 = t164 * t393;
771 t395 = 0.3e1 * t394;
772 t396 = t164 * t173;
773 t397 = 0.3e1 * t396;
774 t398 = t169 * t90;
775 t399 = t398 * t98;
776 t400 = t167 * t399;
777 t401 = 0.3e1 * t400;
778 t402 = t77 * t166;
779 t403 = t402 * t170;
780 t404 = 0.3e1 * t403;
781 t405 = t78 * t179;
782 t406 = 0.6e1 * t405;
783 t407 = t166 * t84;
784 t408 = t9 * t407;
785 t410 = 0.1e1 / t168 / t17;
786 t411 = t410 * t35;
787 t412 = t408 * t411;
788 t413 = 0.2e1 * t412;
789 t414 = t89 * t89;
790 t415 = 0.1e1 / t414;
791 t416 = t178 * t98;
792 t417 = t415 * t416;
793 t418 = t19 * t417;
794 t419 = 0.6e1 * t418;
795 t420 = t150 * t99;
796 t421 = 0.3e1 * t420;
797 t422 = t149 * t84;
798 t423 = t422 * t87;
799 t424 = 0.3e1 * t423;
800 t425 = t77 * t163;
801 t426 = t425 * t87;
802 t427 = 0.3e1 * t426;
803 t428 = t78 * t195;
804 t429 = 0.3e1 * t428;
805 t432 = t158 * t14;
806 t441 = -0.2e1 * t16 * t14 * t74 - 0.3e1 * t16 * t158 * t74 - t16 * t432 * t74 - 0.3e1 * t10 * t146;
807 t442 = t9 * t441;
808 t443 = t442 * t87;
809 t446 = t183 * t25;
810 t456 = t188 * t30;
811 t466 = -0.6e1 * t27 * t183 * t74 - 0.6e1 * t32 * t188 * t74 - 0.4e1 * t27 * t25 * t74 - 0.2e1 * t27 * t446 * t74 - 0.4e1 * t32 * t30 * t74 - 0.2e1 * t32 * t456 * t74 - 0.6e1 * t22 * t146;
812 t467 = t90 * t466;
813 t468 = t19 * t467;
814 t469 = t284 - t382 + t386 + t390 - t392 - t395 - t397 + t401 - t404 + t406 + t413 - t419 - t421 + t424 + t427 - t429 + t443 - t468;
815 t470 = t379 + t469;
816 t471 = t470 * t70;
817 t472 = t471 * t6;
818 t473 = t248 * t74;
819 t474 = 0.6e1 * t473;
820 t475 = t250 * t74;
821 t476 = 0.24e2 * t475;
822 t477 = t247 * t69;
823 t478 = t477 * t6;
824 t479 = 0.6e1 * t478;
825 t480 = t268 - t270 + t272 + t274 - t278 - t284 + t472 - t474 - t476 + t479 + t382 - t386 - t390;
826 t481 = t392 + t395 + t397 - t401 + t404 - t406 - t413 + t419 + t421 - t424 - t427 + t429 - t443 + t468;
827 t484 = 0.3e1 * t171;
828 t485 = 0.6e1 * t174;
829 t486 = 0.6e1 * t180;
830 t487 = 0.6e1 * t155;
831 t488 = 0.6e1 * t153;
832 t489 = 0.3e1 * t165;
833 t490 = 0.3e1 * t196;
834 t491 = 0.12e2 * t251;
835 t492 = 0.12e2 * t253;
836 t493 = 0.24e2 * t257;
837 t494 = 0.18e2 * t259;
838 t495 = 0.3e1 * t151;
839 t496 = 0.3e1 * t249;
840 t497 = 0.6e1 * t255;
841 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
842 v3rho3[0] = t1 * (t480 + t481) + t484 + t485 - t486 + t487 - t488 - t489 + t490 + t491 - t492 - t493 + t494 - t495 + t496 + t497;
843
844 t498 = 0.2e1 * t267;
845 t499 = 0.4e1 * t269;
846 t500 = 0.12e2 * t273;
847 t501 = 0.8e1 * t475;
848 t502 = 0.2e1 * t478;
849 t503 = -t498 + t499 + t272 + t500 - t278 - t284 + t472 - t474 - t501 + t502 + t382 - t386 - t390;
850 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
851 v3rho3[1] = t1 * (t503 + t481) + t484 + t485 - t486 + t487 - t488 - t489 + t490 + t252 - t492 - t258 + t494 - t495 + t496 - t256;
852
853 t506 = -t498 + t499 + t272 - t500 - t278 - t284 + t472 - t474 + t501 - t502 + t382 - t386 - t390;
854 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
855 v3rho3[2] = t1 * (t506 + t481) + t484 + t485 - t486 + t487 - t488 - t489 + t490 - t252 - t492 + t258 + t494 - t495 + t496 - t256;
856
857 t509 = t268 - t270 + t272 - t274 - t278 - t284 + t472 - t474 + t476 - t479 + t382 - t386 - t390;
858 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
859 v3rho3[3] = t1 * (t509 + t481) + t484 + t485 - t486 + t487 - t488 - t489 + t490 - t491 - t492 + t493 + t494 - t495 + t496 + t497;
860
861 #ifndef XC_DONT_COMPILE_LXC
862
863 if(order < 4) return;
864
865
866 t512 = t127 * t74;
867 t514 = t68 * t146;
868 t516 = t247 * t6;
869 t521 = 0.36e2 * t19 * t415 * t178 * t194;
870 t524 = 0.8e1 * t19 * t383 * t466;
871 t526 = 0.24e2 * t78 * t384;
872 t528 = 0.4e1 * t442 * t393;
873 t530 = 0.4e1 * t442 * t173;
874 t532 = 0.6e1 * t164 * t380;
875 t536 = 0.8e1 * t408 * t410 * t90 * t98;
876 t539 = 0.6e1 * t167 * t398 * t194;
877 t542 = 0.4e1 * t85 * t172 * t466;
878 t543 = -0.48e2 * t512 + 0.72e2 * t514 + 0.12e2 * t516 + t521 - t524 - t526 + t528 + t530 + t532 + t536 - t539 + t542;
879 t545 = 0.8e1 * t471 * t74;
880 t546 = 0.12e2 * t276;
881 t548 = 0.1e1 / t145 / t5;
882 t553 = (0.30e2 * t4 * t548 + t546) * t18 * t35;
883 t580 = t208 * t208;
884 t597 = t233 * t233;
885 t610 = t238 * t238;
886 t635 = t228 * t228;
887 t639 = t244 * t244;
888 t645 = t213 * t213;
889 t649 = t216 * t216;
890 t651 = t218 * t218;
891 t659 = -0.12e2 * t214 * t343 * t216 - 0.8e1 * t50 * t364 * t331 - (0.30e2 * t37 * t548 + t546) * t49 * t66 + 0.6e1 * t200 * t245 - 0.4e1 * t376 * t110 * t113 - 0.6e1 * t199 * t213 * t113 - 0.4e1 * t103 * t306 * t113 + 0.4e1 * t104 * t332 - t40 * (0.11e2 * t47 * t208 * t146 + 0.6e1 * t47 * t297 * t146 + 0.6e1 * t47 * t45 * t146 + t47 * t580 * t146 + 0.12e2 * t41 * t276) * t113 + t50 * t116 * (0.22e2 * t58 * t233 * t146 + 0.22e2 * t63 * t238 * t146 + 0.12e2 * t58 * t311 * t146 + 0.12e2 * t63 * t321 * t146 + 0.12e2 * t58 * t56 * t146 + 0.2e1 * t58 * t597 * t146 + 0.12e2 * t63 * t61 * t146 + 0.2e1 * t63 * t610 * t146 + 0.24e2 * t53 * t276) + 0.6e1 * t199 * t216 * t220 - 0.12e2 * t200 * t229 + 0.4e1 * t377 * t125 - 0.24e2 * t50 / t346 / t65 * t635 - 0.6e1 * t50 * t227 * t639 + 0.24e2 * t104 * t349 + 0.3e1 * t40 * t645 * t220 + 0.6e1 * t40 * t649 / t651 * t66 - 0.8e1 * t103 * t339 * t343;
892 t705 = 0.4e1 * t307 * t352 + 0.4e1 * t307 * t223 + 0.6e1 * t214 * t361 + 0.8e1 * t340 * t342 * t116 * t124 - 0.6e1 * t217 * t357 * t244 + 0.4e1 * t111 * t222 * t331 + 0.36e2 * t50 * t347 * t228 * t244 - 0.24e2 * t202 * t371 - 0.12e2 * t214 * t371 + 0.12e2 * t202 * t220 * t213 - 0.12e2 * t334 * t358 + 0.12e2 * t287 * t223 + 0.12e2 * t290 * t223 + 0.12e2 * t202 * t361 - 0.24e2 * t104 * t365 + 0.12e2 * t217 * t219 * t227 * t228 + 0.24e2 * t111 * t112 * t347 * t348 - t521 + t524;
893 t709 = 0.4e1 * t282 * t84 * t87;
894 t712 = 0.6e1 * t149 * t163 * t87;
895 t715 = 0.4e1 * t77 * t441 * t87;
896 t717 = 0.4e1 * t78 * t467;
897 t720 = t158 * t158;
898 t734 = t9 * (0.6e1 * t16 * t14 * t146 + 0.11e2 * t16 * t158 * t146 + 0.6e1 * t16 * t432 * t146 + t16 * t720 * t146 + 0.12e2 * t10 * t276) * t87;
899 t737 = t183 * t183;
900 t750 = t188 * t188;
901 t765 = t19 * t90 * (0.22e2 * t27 * t183 * t146 + 0.22e2 * t32 * t188 * t146 + 0.12e2 * t27 * t25 * t146 + 0.12e2 * t27 * t446 * t146 + 0.2e1 * t27 * t737 * t146 + 0.12e2 * t32 * t30 * t146 + 0.12e2 * t32 * t456 * t146 + 0.2e1 * t32 * t750 * t146 + 0.24e2 * t22 * t276);
902 t768 = t178 * t178;
903 t771 = 0.24e2 * t19 / t414 / t34 * t768;
904 t772 = t194 * t194;
905 t775 = 0.6e1 * t19 * t177 * t772;
906 t778 = 0.6e1 * t149 * t166 * t170;
907 t780 = 0.12e2 * t150 * t179;
908 t782 = 0.4e1 * t283 * t99;
909 t783 = t526 - t528 - t530 - t532 - t536 + t539 - t542 + t553 + t709 + t712 + t715 - t717 + t734 - t765 + t771 + t775 - t778 + t780 - t782;
910 t785 = 0.6e1 * t150 * t195;
911 t787 = 0.24e2 * t78 * t417;
912 t788 = t163 * t163;
913 t791 = 0.3e1 * t9 * t788 * t170;
914 t792 = t166 * t166;
915 t794 = t168 * t168;
916 t798 = 0.6e1 * t9 * t792 / t794 * t35;
917 t801 = 0.8e1 * t77 * t407 * t411;
918 t804 = 0.12e2 * t164 * t411 * t166;
919 t806 = 0.24e2 * t152 * t388;
920 t808 = 0.12e2 * t164 * t388;
921 t811 = 0.12e2 * t152 * t170 * t163;
922 t813 = 0.12e2 * t402 * t399;
923 t815 = 0.12e2 * t422 * t173;
924 t817 = 0.12e2 * t425 * t173;
925 t819 = 0.12e2 * t152 * t380;
926 t823 = 0.12e2 * t167 * t169 * t177 * t178;
927 t827 = 0.24e2 * t85 * t86 * t415 * t416;
928 t830 = 0.24e2 * t85 * t86 * t384;
929 t835 = 0.12e2 * t164 * t169 * t90 * t84 * t98;
930 t844 = -0.12e2 * t214 * t219 * t116 * t110 * t124 - 0.24e2 * t111 * t112 * t365 - t785 - t787 - t791 - t798 + t801 + t804 + t806 + t808 - t811 + t813 - t815 - t817 - t819 - t823 - t827 + t830 + t835;
931 t848 = (t659 + t705 + t783 + t844) * t70 * t6;
932 t850 = t470 * t69 * t6;
933 t851 = 0.8e1 * t850;
934 t853 = 0.120e3 * t71 * t548;
935 t854 = t250 * t146;
936 t855 = 0.144e3 * t854;
937 t857 = 0.96e2 * t128 * t276;
938 t858 = t130 * t276;
939 t859 = 0.192e3 * t858;
940 t861 = 0.36e2 * t248 * t146;
941 t862 = t477 * t74;
942 t863 = 0.48e2 * t862;
943 t864 = -t545 - t553 + t848 + t851 + t853 + t855 - t857 - t859 + t861 - t863 - t709 - t712 - t715;
944 t866 = t717 - t734 + t765 - t771 - t775 + t778 - t780 + t782 + t785 + t787 + t791 + t798;
945 t867 = -t801 - t804 - t806 - t808 + t811 - t813 + t815 + t817 + t819 + t823 + t827 - t830 - t835;
946 t868 = t866 + t867;
947 t871 = 0.24e2 * t267;
948 t872 = 0.48e2 * t269;
949 t873 = 0.72e2 * t271;
950 t874 = 0.144e3 * t273;
951 t875 = 0.96e2 * t277;
952 t876 = 0.4e1 * t284;
953 t877 = 0.4e1 * t472;
954 t878 = 0.24e2 * t473;
955 t879 = 0.96e2 * t475;
956 t880 = 0.24e2 * t478;
957 t881 = 0.12e2 * t381;
958 t882 = 0.24e2 * t385;
959 t883 = 0.24e2 * t389;
960 t884 = t1 * (t543 + t864 + t868) + t871 - t872 + t873 + t874 - t875 - t876 + t877 - t878 - t879 + t880 + t881 - t882 - t883;
961 t885 = 0.24e2 * t391;
962 t886 = 0.12e2 * t394;
963 t887 = 0.12e2 * t396;
964 t888 = 0.12e2 * t400;
965 t889 = 0.12e2 * t403;
966 t890 = 0.24e2 * t405;
967 t891 = 0.8e1 * t412;
968 t892 = 0.24e2 * t418;
969 t893 = 0.12e2 * t420;
970 t894 = 0.12e2 * t423;
971 t895 = 0.12e2 * t426;
972 t896 = 0.12e2 * t428;
973 t897 = 0.4e1 * t443;
974 t898 = 0.4e1 * t468;
975 t899 = t885 + t886 + t887 - t888 + t889 - t890 - t891 + t892 + t893 - t894 - t895 + t896 - t897 + t898;
976 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
977 v4rho4[0] = t884 + t899;
978
979 t900 = t521 - t524 - t526 + t528 + t530 + t532 + t536 - t539 + t542 - t545 - t553;
980 t901 = 0.4e1 * t850;
981 t902 = 0.72e2 * t854;
982 t903 = 0.96e2 * t858;
983 t904 = 0.24e2 * t862;
984 t905 = t848 + t901 + t853 + t902 - t857 - t903 + t861 - t904 - t709 - t712 - t715 + t717;
985 t907 = -t734 + t765 - t771 - t775 + t778 - t780 + t782 + t785 + t787 + t791 + t798 - t801;
986 t908 = -t804 - t806 - t808 + t811 - t813 + t815 + t817 + t819 + t823 + t827 - t830 - t835;
987 t909 = t907 + t908;
988 t912 = 0.72e2 * t273;
989 t913 = 0.48e2 * t475;
990 t914 = 0.12e2 * t478;
991 t915 = t1 * (t900 + t905 + t909) + t873 + t912 - t875 - t876 + t877 - t878 - t913 + t914 + t881 - t882 - t883 + t885;
992 t916 = t886 + t887 - t888 + t889 - t890 - t891 + t892 + t893 - t894 - t895 + t896 - t897 + t898;
993 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
994 v4rho4[1] = t915 + t916;
995
996 t922 = 0.16e2 * t512 - 0.24e2 * t514 - 0.4e1 * t516 + t521 - t524 - t526 + t528 + t530 + t532 + t536 - t539;
997 t923 = t542 - t545 - t553 + t848 + t853 - t857 + t861 - t709 - t712 - t715 + t717 - t734;
998 t925 = t765 - t771 - t775 + t778 - t780 + t782 + t785 + t787 + t791 + t798 - t801;
999 t929 = -0.8e1 * t267 + 0.16e2 * t269 + t1 * (t922 + t923 + t925 + t908) + t873 - t875 - t876 + t877 - t878 + t881 - t882 - t883 + t885;
1000 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1001 v4rho4[2] = t929 + t916;
1002
1003 t930 = t848 - t901 + t853 - t902 - t857 + t903 + t861 + t904 - t709 - t712 - t715 + t717;
1004 t934 = t1 * (t900 + t930 + t909) + t873 - t912 - t875 - t876 + t877 - t878 + t913 - t914 + t881 - t882 - t883 + t885;
1005 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1006 v4rho4[3] = t934 + t916;
1007
1008 t935 = -t545 - t553 + t848 - t851 + t853 - t855 - t857 + t859 + t861 + t863 - t709 - t712 - t715;
1009 t939 = t1 * (t543 + t935 + t868) + t871 - t872 + t873 - t874 - t875 - t876 + t877 - t878 + t879 - t880 + t881 - t882 - t883;
1010 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
1011 v4rho4[4] = t939 + t899;
1012
1013 #ifndef XC_DONT_COMPILE_MXC
1014
1015 if(order < 5) return;
1016
1017
1018 #endif
1019
1020 #endif
1021
1022 #endif
1023
1024 #endif
1025
1026 #endif
1027
1028
1029 }
1030
1031