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