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_2d_prm.mpl
11   Type of functional: lda_exc
12 */
13 
14 #define maple2c_order 4
15 #define MAPLE2C_FLAGS (XC_FLAGS_I_HAVE_EXC | XC_FLAGS_I_HAVE_VXC | XC_FLAGS_I_HAVE_FXC | XC_FLAGS_I_HAVE_KXC | XC_FLAGS_I_HAVE_LXC)
16 
17 
18 static inline void
func_unpol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
20 {
21 
22 #ifndef XC_DONT_COMPILE_EXC
23   double t1, t3, t5, t6, t7, t9, t10, t11;
24   double t12, t13, t15, t16, t17, t19, t20, t21;
25   double t23, t25, t26, t27, t28, t30, t31, t33;
26 
27 #ifndef XC_DONT_COMPILE_VXC
28   double t34, t35, t36, t38, t41, t42, t43, t45;
29   double t48, t50, t51, t54, t57, t58, t60, t62;
30   double t64, t66, t68;
31 
32 #ifndef XC_DONT_COMPILE_FXC
33   double t82, t83, t84, t86, t87, t89, t91, t92;
34   double t96, t97, t98, t100, t102, t104, t107, t108;
35   double t111, t113, t114, t117, t120, t122, t123, t125;
36   double t127, t129, t131, t133, t135, t138, t140;
37 
38 #ifndef XC_DONT_COMPILE_KXC
39   double t160, t162, t163, t164, t167, t169, t171, t172;
40   double t174, t175, t177, t178, t181, t182, t186, t187;
41   double t188, t191, t194, t197, t198, t199, t201, t203;
42   double t205, t206, t208, t210, t211, t214, t216, t218;
43   double t220, t222, t224, t226, t229, t232, t233, t235;
44   double t236, t238;
45 
46 #ifndef XC_DONT_COMPILE_LXC
47   double t254, t258, t260, t261, t269, t274, t280, t286;
48   double t289, t292, t297, t301, t303, t304, t319, t362;
49   double t374;
50 #endif
51 
52 #endif
53 
54 #endif
55 
56 #endif
57 
58 #endif
59 
60 
61   lda_c_2d_prm_params *params;
62 
63   assert(p->params != NULL);
64   params = (lda_c_2d_prm_params * )(p->params);
65 
66   assert(params->N > 1);
67 
68   t1 = sqrt(rho[0]);
69   t3 = sqrt(M_PI);
70   t5 = 0.39274000000000000001e1 * t1 + t3 / 0.2e1;
71   t6 = 0.1e1 / t5;
72   t7 = t1 * t6;
73   t9 = 0.39274000000000000001e1 * t7 - 0.1e1;
74   t10 = t1 * t9;
75   t11 = 0.2e1 + params->c;
76   t12 = sqrt(t11);
77   t13 = 0.1e1 / t12;
78   t15 = 0.35445383694248789566e0 * t10 * t13;
79   t16 = 0.1e1 / t11;
80   t17 = t9 * t16;
81   t19 = 0.39995832530297311691e0 * t7 * t17;
82   t20 = t5 * t5;
83   t21 = 0.1e1 / t20;
84   t23 = 0.1e1/POW_3_2(t11);
85   t25 = 0.17722691847124394783e0 * t1 * t21 * t23;
86   t26 = 0.1e1 + params->c;
87   t27 = sqrt(t26);
88   t28 = 0.1e1 / t27;
89   t30 = 0.70890767388497579132e0 * t10 * t28;
90   t31 = 0.1e1 / t26;
91   t33 = 0.39995832530297311691e0 * t7 * t31;
92   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
93     zk[0] = t15 + t19 + t25 + t30 + t33;
94 
95 #ifndef XC_DONT_COMPILE_VXC
96 
97   if(order < 1) return;
98 
99 
100   t34 = 0.1e1 / t1;
101   t35 = t34 * t9;
102   t36 = t35 * t13;
103   t38 = t34 * t6;
104   t41 = 0.19637000000000000000e1 * t38 - 0.77122353800000000002e1 * t21;
105   t42 = t1 * t41;
106   t43 = t42 * t13;
107   t45 = t38 * t17;
108   t48 = t21 * t9 * t16;
109   t50 = t41 * t16;
110   t51 = t7 * t50;
111   t54 = t34 * t21 * t23;
112   t57 = 0.1e1 / t20 / t5;
113   t58 = t57 * t23;
114   t60 = t35 * t28;
115   t62 = t42 * t28;
116   t64 = t38 * t31;
117   t66 = t21 * t31;
118   t68 = 0.17722691847124394783e0 * t36 + 0.35445383694248789566e0 * t43 + 0.19997916265148655846e0 * t45 - 0.78539816339744830968e0 * t48 + 0.39995832530297311691e0 * t51 + 0.88613459235621973915e-1 * t54 - 0.69604099960396348071e0 * t58 + 0.35445383694248789566e0 * t60 + 0.70890767388497579132e0 * t62 + 0.19997916265148655846e0 * t64 - 0.78539816339744830968e0 * t66;
119   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
120     vrho[0] = rho[0] * t68 + t15 + t19 + t25 + t30 + t33;
121 
122 #ifndef XC_DONT_COMPILE_FXC
123 
124   if(order < 2) return;
125 
126 
127   t82 = 0.1e1 / t1 / rho[0];
128   t83 = t82 * t9;
129   t84 = t83 * t13;
130   t86 = t34 * t41;
131   t87 = t86 * t13;
132   t89 = t82 * t6;
133   t91 = 0.1e1 / rho[0];
134   t92 = t91 * t21;
135   t96 = -0.98185000000000000000e0 * t89 - 0.38561176900000000000e1 * t92 + 0.30289033231412000001e2 * t57 * t34;
136   t97 = t1 * t96;
137   t98 = t97 * t13;
138   t100 = t89 * t17;
139   t102 = t92 * t17;
140   t104 = t38 * t50;
141   t107 = t16 * t34;
142   t108 = t57 * t9 * t107;
143   t111 = t21 * t41 * t16;
144   t113 = t96 * t16;
145   t114 = t7 * t113;
146   t117 = t82 * t21 * t23;
147   t120 = t91 * t57 * t23;
148   t122 = t20 * t20;
149   t123 = 0.1e1 / t122;
150   t125 = t123 * t23 * t34;
151   t127 = t83 * t28;
152   t129 = t86 * t28;
153   t131 = t97 * t28;
154   t133 = t89 * t31;
155   t135 = t92 * t31;
156   t138 = t57 * t31 * t34;
157   t140 = -0.88613459235621973915e-1 * t84 + 0.35445383694248789566e0 * t87 + 0.35445383694248789566e0 * t98 - 0.99989581325743279230e-1 * t100 - 0.39269908169872415485e0 * t102 + 0.39995832530297311692e0 * t104 + 0.30845727469271384914e1 * t108 - 0.15707963267948966194e1 * t111 + 0.39995832530297311691e0 * t114 - 0.44306729617810986958e-1 * t117 - 0.34802049980198174035e0 * t120 + 0.41004471327669092612e1 * t125 - 0.17722691847124394783e0 * t127 + 0.70890767388497579132e0 * t129 + 0.70890767388497579132e0 * t131 - 0.99989581325743279230e-1 * t133 - 0.39269908169872415485e0 * t135 + 0.30845727469271384914e1 * t138;
158   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
159     v2rho2[0] = 0.35445383694248789566e0 * t36 + 0.70890767388497579132e0 * t43 + 0.39995832530297311692e0 * t45 - 0.15707963267948966194e1 * t48 + 0.79991665060594623382e0 * t51 + 0.17722691847124394783e0 * t54 - 0.13920819992079269614e1 * t58 + 0.70890767388497579132e0 * t60 + 0.14178153477699515826e1 * t62 + 0.39995832530297311692e0 * t64 - 0.15707963267948966194e1 * t66 + rho[0] * t140;
160 
161 #ifndef XC_DONT_COMPILE_KXC
162 
163   if(order < 3) return;
164 
165 
166   t160 = rho[0] * rho[0];
167   t162 = 0.1e1 / t1 / t160;
168   t163 = t162 * t9;
169   t164 = t163 * t13;
170   t167 = t162 * t21 * t23;
171   t169 = t163 * t28;
172   t171 = t162 * t6;
173   t172 = t171 * t31;
174   t174 = t82 * t41;
175   t175 = t174 * t13;
176   t177 = t34 * t96;
177   t178 = t177 * t13;
178   t181 = 0.1e1 / t160;
179   t182 = t181 * t21;
180   t186 = 0.14727750000000000000e1 * t171 + 0.57841765350000000000e1 * t182 - 0.17843572366957123321e3 * t123 * t91;
181   t187 = t1 * t186;
182   t188 = t187 * t13;
183   t191 = t21 * t96 * t16;
184   t194 = t181 * t57 * t23;
185   t197 = 0.1e1 / t122 / t5;
186   t198 = t197 * t23;
187   t199 = t198 * t91;
188   t201 = t174 * t28;
189   t203 = t177 * t28;
190   t205 = 0.13292018885343296087e0 * t164 + 0.66460094426716480437e-1 * t167 + 0.26584037770686592174e0 * t169 + 0.14998437198861491884e0 * t172 - 0.26584037770686592174e0 * t175 + 0.53168075541373184349e0 * t178 + 0.35445383694248789566e0 * t188 - 0.23561944901923449291e1 * t191 + 0.52203074970297261053e0 * t194 - 0.32208192138457518865e2 * t199 - 0.53168075541373184349e0 * t201 + 0.10633615108274636870e1 * t203;
191   t206 = t187 * t28;
192   t208 = t182 * t31;
193   t210 = t82 * t57;
194   t211 = t210 * t31;
195   t214 = t123 * t31 * t91;
196   t216 = t171 * t17;
197   t218 = t89 * t50;
198   t220 = t182 * t17;
199   t222 = t92 * t50;
200   t224 = t210 * t17;
201   t226 = t38 * t113;
202   t229 = t57 * t41 * t107;
203   t232 = t16 * t91;
204   t233 = t123 * t9 * t232;
205   t235 = t186 * t16;
206   t236 = t7 * t235;
207   t238 = 0.70890767388497579132e0 * t206 + 0.58904862254808623227e0 * t208 + 0.1e-18 * t211 - 0.18171526509422465567e2 * t214 + 0.14998437198861491884e0 * t216 - 0.29996874397722983769e0 * t218 + 0.58904862254808623227e0 * t220 - 0.11780972450961724646e1 * t222 + 0.1e-18 * t224 + 0.59993748795445967538e0 * t226 + 0.92537182407814154744e1 * t229 - 0.18171526509422465567e2 * t233 + 0.39995832530297311691e0 * t236;
208   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
209     v3rho3[0] = -0.29996874397722983769e0 * t100 - 0.26584037770686592174e0 * t84 + 0.11998749759089193508e1 * t104 - 0.11780972450961724646e1 * t102 + 0.92537182407814154744e1 * t108 + 0.11998749759089193507e1 * t114 - 0.13292018885343296087e0 * t117 - 0.53168075541373184349e0 * t127 - 0.29996874397722983769e0 * t133 + 0.10633615108274636870e1 * t87 + 0.10633615108274636870e1 * t98 - 0.47123889803846898582e1 * t111 - 0.10440614994059452211e1 * t120 + 0.12301341398300727784e2 * t125 + 0.21267230216549273739e1 * t129 + 0.21267230216549273739e1 * t131 - 0.11780972450961724646e1 * t135 + 0.92537182407814154744e1 * t138 + rho[0] * (t205 + t238);
210 
211 #ifndef XC_DONT_COMPILE_LXC
212 
213   if(order < 4) return;
214 
215 
216   t254 = 0.53168075541373184348e0 * t164 + 0.26584037770686592174e0 * t167 + 0.10633615108274636870e1 * t169 + 0.59993748795445967538e0 * t172 - 0.10633615108274636870e1 * t175 + 0.21267230216549273740e1 * t178 + 0.14178153477699515827e1 * t188 - 0.94247779607693797163e1 * t191 + 0.20881229988118904421e1 * t194 - 0.12883276855383007546e3 * t199 - 0.21267230216549273740e1 * t201 + 0.42534460433098547479e1 * t203 + 0.28356306955399031652e1 * t206;
217   t258 = t160 * rho[0];
218   t260 = 0.1e1 / t1 / t258;
219   t261 = t260 * t9;
220   t269 = t260 * t6;
221   t274 = t162 * t41;
222   t280 = 0.1e1 / t258;
223   t286 = t280 * t21;
224   t289 = t82 * t96;
225   t292 = t34 * t186;
226   t297 = t162 * t57;
227   t301 = t123 * t181;
228   t303 = -0.36819375000000000000e1 * t269 - 0.14460441337500000000e2 * t286 - 0.22716774923559000000e2 * t297 + 0.14015769222797481226e4 * t197 * t82 + 0.17843572366957123321e3 * t301;
229   t304 = t1 * t303;
230   t319 = -0.33230047213358240218e0 * t261 * t13 - 0.16615023606679120109e0 * t260 * t21 * t23 - 0.66460094426716480435e0 * t261 * t28 - 0.37496092997153729710e0 * t269 * t31 + 0.32208192138457518865e2 * t198 * t181 + 0.53168075541373184348e0 * t274 * t13 - 0.31415926535897932388e1 * t21 * t186 * t16 - 0.13050768742574315264e1 * t280 * t57 * t23 + 0.10633615108274636870e1 * t274 * t28 - 0.14726215563702155806e1 * t286 * t31 - 0.53168075541373184348e0 * t289 * t13 + 0.70890767388497579132e0 * t292 * t13 + 0.35445383694248789566e0 * t304 * t13 - 0.30753353495751819459e1 * t162 * t123 * t23 + 0.31623613451144514898e3 / t122 / t20 * t23 * t82 - 0.10633615108274636870e1 * t289 * t28 + 0.14178153477699515827e1 * t292 * t28;
231   t362 = 0.70890767388497579132e0 * t304 * t28 - 0.23134295601953538688e1 * t297 * t31 + 0.18171526509422465566e2 * t301 * t31 + 0.14273370642621158254e3 * t197 * t31 * t82 - 0.37496092997153729710e0 * t269 * t17 + 0.59993748795445967538e0 * t171 * t50 - 0.14726215563702155806e1 * t286 * t17 - 0.59993748795445967538e0 * t89 * t113 + 0.23561944901923449291e1 * t182 * t50 - 0.23134295601953538688e1 * t297 * t17 - 0.23561944901923449291e1 * t92 * t113 + 0.4e-18 * t210 * t50 + 0.18171526509422465566e2 * t301 * t17 + 0.79991665060594623384e0 * t38 * t235 + 0.18507436481562830949e2 * t57 * t96 * t107 - 0.72686106037689862268e2 * t123 * t41 * t232 + 0.14273370642621158254e3 * t197 * t9 * t16 * t82 + 0.39995832530297311691e0 * t7 * t303 * t16;
232   t374 = 0.23561944901923449292e1 * t208 + 0.4e-18 * t211 - 0.72686106037689862268e2 * t214 + rho[0] * (t319 + t362) + 0.59993748795445967538e0 * t216 - 0.11998749759089193508e1 * t218 + 0.23561944901923449292e1 * t220 - 0.47123889803846898584e1 * t222 + 0.4e-18 * t224 + 0.23997499518178387016e1 * t226 + 0.37014872963125661897e2 * t229 - 0.72686106037689862268e2 * t233 + 0.15998333012118924676e1 * t236;
233   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
234     v4rho4[0] = t254 + t374;
235 
236 #ifndef XC_DONT_COMPILE_MXC
237 
238   if(order < 5) return;
239 
240 
241 #endif
242 
243 #endif
244 
245 #endif
246 
247 #endif
248 
249 #endif
250 
251 
252 }
253 
254 
255 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))256 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
257 {
258 
259 #ifndef XC_DONT_COMPILE_EXC
260   double t1, t3, t5, t6, t7, t9, t10, t11;
261   double t12, t13, t15, t16, t17, t19, t20, t21;
262   double t23, t25, t26, t27, t28, t30, t31, t33;
263 
264 #ifndef XC_DONT_COMPILE_VXC
265   double t34, t35, t36, t38, t41, t42, t43, t45;
266   double t48, t50, t51, t54, t57, t58, t60, t62;
267   double t64, t66, t68;
268 
269 #ifndef XC_DONT_COMPILE_FXC
270   double t82, t83, t84, t86, t87, t89, t91, t92;
271   double t96, t97, t98, t100, t102, t104, t107, t108;
272   double t111, t113, t114, t117, t120, t122, t123, t125;
273   double t127, t129, t131, t133, t135, t138, t140;
274 
275 #ifndef XC_DONT_COMPILE_KXC
276   double t160, t162, t163, t164, t167, t169, t171, t172;
277   double t174, t175, t177, t178, t181, t182, t186, t187;
278   double t188, t191, t194, t197, t198, t199, t201, t203;
279   double t205, t206, t208, t210, t211, t214, t216, t218;
280   double t220, t222, t224, t226, t229, t232, t233, t235;
281   double t236, t238;
282 
283 #ifndef XC_DONT_COMPILE_LXC
284   double t254, t258, t260, t261, t269, t274, t280, t286;
285   double t289, t292, t297, t301, t303, t304, t319, t362;
286   double t374;
287 #endif
288 
289 #endif
290 
291 #endif
292 
293 #endif
294 
295 #endif
296 
297 
298   lda_c_2d_prm_params *params;
299 
300   assert(p->params != NULL);
301   params = (lda_c_2d_prm_params * )(p->params);
302 
303   assert(params->N > 1);
304 
305   t1 = sqrt(rho[0]);
306   t3 = sqrt(M_PI);
307   t5 = 0.39274000000000000001e1 * t1 + t3 / 0.2e1;
308   t6 = 0.1e1 / t5;
309   t7 = t1 * t6;
310   t9 = 0.39274000000000000001e1 * t7 - 0.1e1;
311   t10 = t1 * t9;
312   t11 = 0.2e1 + params->c;
313   t12 = sqrt(t11);
314   t13 = 0.1e1 / t12;
315   t15 = 0.35445383694248789566e0 * t10 * t13;
316   t16 = 0.1e1 / t11;
317   t17 = t9 * t16;
318   t19 = 0.39995832530297311691e0 * t7 * t17;
319   t20 = t5 * t5;
320   t21 = 0.1e1 / t20;
321   t23 = 0.1e1/POW_3_2(t11);
322   t25 = 0.17722691847124394783e0 * t1 * t21 * t23;
323   t26 = 0.1e1 + params->c;
324   t27 = sqrt(t26);
325   t28 = 0.1e1 / t27;
326   t30 = 0.70890767388497579132e0 * t10 * t28;
327   t31 = 0.1e1 / t26;
328   t33 = 0.39995832530297311691e0 * t7 * t31;
329   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
330     zk[0] = t15 + t19 + t25 + t30 + t33;
331 
332 #ifndef XC_DONT_COMPILE_VXC
333 
334   if(order < 1) return;
335 
336 
337   t34 = 0.1e1 / t1;
338   t35 = t34 * t9;
339   t36 = t35 * t13;
340   t38 = t34 * t6;
341   t41 = 0.19637000000000000000e1 * t38 - 0.77122353800000000002e1 * t21;
342   t42 = t1 * t41;
343   t43 = t42 * t13;
344   t45 = t38 * t17;
345   t48 = t21 * t9 * t16;
346   t50 = t41 * t16;
347   t51 = t7 * t50;
348   t54 = t34 * t21 * t23;
349   t57 = 0.1e1 / t20 / t5;
350   t58 = t57 * t23;
351   t60 = t35 * t28;
352   t62 = t42 * t28;
353   t64 = t38 * t31;
354   t66 = t21 * t31;
355   t68 = 0.17722691847124394783e0 * t36 + 0.35445383694248789566e0 * t43 + 0.19997916265148655846e0 * t45 - 0.78539816339744830968e0 * t48 + 0.39995832530297311691e0 * t51 + 0.88613459235621973915e-1 * t54 - 0.69604099960396348071e0 * t58 + 0.35445383694248789566e0 * t60 + 0.70890767388497579132e0 * t62 + 0.19997916265148655846e0 * t64 - 0.78539816339744830968e0 * t66;
356   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
357     vrho[0] = rho[0] * t68 + t15 + t19 + t25 + t30 + t33;
358 
359 #ifndef XC_DONT_COMPILE_FXC
360 
361   if(order < 2) return;
362 
363 
364   t82 = 0.1e1 / t1 / rho[0];
365   t83 = t82 * t9;
366   t84 = t83 * t13;
367   t86 = t34 * t41;
368   t87 = t86 * t13;
369   t89 = t82 * t6;
370   t91 = 0.1e1 / rho[0];
371   t92 = t91 * t21;
372   t96 = -0.98185000000000000000e0 * t89 - 0.38561176900000000000e1 * t92 + 0.30289033231412000001e2 * t57 * t34;
373   t97 = t1 * t96;
374   t98 = t97 * t13;
375   t100 = t89 * t17;
376   t102 = t92 * t17;
377   t104 = t38 * t50;
378   t107 = t16 * t34;
379   t108 = t57 * t9 * t107;
380   t111 = t21 * t41 * t16;
381   t113 = t96 * t16;
382   t114 = t7 * t113;
383   t117 = t82 * t21 * t23;
384   t120 = t91 * t57 * t23;
385   t122 = t20 * t20;
386   t123 = 0.1e1 / t122;
387   t125 = t123 * t23 * t34;
388   t127 = t83 * t28;
389   t129 = t86 * t28;
390   t131 = t97 * t28;
391   t133 = t89 * t31;
392   t135 = t92 * t31;
393   t138 = t57 * t31 * t34;
394   t140 = -0.88613459235621973915e-1 * t84 + 0.35445383694248789566e0 * t87 + 0.35445383694248789566e0 * t98 - 0.99989581325743279230e-1 * t100 - 0.39269908169872415485e0 * t102 + 0.39995832530297311692e0 * t104 + 0.30845727469271384914e1 * t108 - 0.15707963267948966194e1 * t111 + 0.39995832530297311691e0 * t114 - 0.44306729617810986958e-1 * t117 - 0.34802049980198174035e0 * t120 + 0.41004471327669092612e1 * t125 - 0.17722691847124394783e0 * t127 + 0.70890767388497579132e0 * t129 + 0.70890767388497579132e0 * t131 - 0.99989581325743279230e-1 * t133 - 0.39269908169872415485e0 * t135 + 0.30845727469271384914e1 * t138;
395   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
396     v2rho2[0] = 0.35445383694248789566e0 * t36 + 0.70890767388497579132e0 * t43 + 0.39995832530297311692e0 * t45 - 0.15707963267948966194e1 * t48 + 0.79991665060594623382e0 * t51 + 0.17722691847124394783e0 * t54 - 0.13920819992079269614e1 * t58 + 0.70890767388497579132e0 * t60 + 0.14178153477699515826e1 * t62 + 0.39995832530297311692e0 * t64 - 0.15707963267948966194e1 * t66 + rho[0] * t140;
397 
398 #ifndef XC_DONT_COMPILE_KXC
399 
400   if(order < 3) return;
401 
402 
403   t160 = rho[0] * rho[0];
404   t162 = 0.1e1 / t1 / t160;
405   t163 = t162 * t9;
406   t164 = t163 * t13;
407   t167 = t162 * t21 * t23;
408   t169 = t163 * t28;
409   t171 = t162 * t6;
410   t172 = t171 * t31;
411   t174 = t82 * t41;
412   t175 = t174 * t13;
413   t177 = t34 * t96;
414   t178 = t177 * t13;
415   t181 = 0.1e1 / t160;
416   t182 = t181 * t21;
417   t186 = 0.14727750000000000000e1 * t171 + 0.57841765350000000000e1 * t182 - 0.17843572366957123321e3 * t123 * t91;
418   t187 = t1 * t186;
419   t188 = t187 * t13;
420   t191 = t21 * t96 * t16;
421   t194 = t181 * t57 * t23;
422   t197 = 0.1e1 / t122 / t5;
423   t198 = t197 * t23;
424   t199 = t198 * t91;
425   t201 = t174 * t28;
426   t203 = t177 * t28;
427   t205 = 0.13292018885343296087e0 * t164 + 0.66460094426716480437e-1 * t167 + 0.26584037770686592174e0 * t169 + 0.14998437198861491884e0 * t172 - 0.26584037770686592174e0 * t175 + 0.53168075541373184349e0 * t178 + 0.35445383694248789566e0 * t188 - 0.23561944901923449291e1 * t191 + 0.52203074970297261053e0 * t194 - 0.32208192138457518865e2 * t199 - 0.53168075541373184349e0 * t201 + 0.10633615108274636870e1 * t203;
428   t206 = t187 * t28;
429   t208 = t182 * t31;
430   t210 = t82 * t57;
431   t211 = t210 * t31;
432   t214 = t123 * t31 * t91;
433   t216 = t171 * t17;
434   t218 = t89 * t50;
435   t220 = t182 * t17;
436   t222 = t92 * t50;
437   t224 = t210 * t17;
438   t226 = t38 * t113;
439   t229 = t57 * t41 * t107;
440   t232 = t16 * t91;
441   t233 = t123 * t9 * t232;
442   t235 = t186 * t16;
443   t236 = t7 * t235;
444   t238 = 0.70890767388497579132e0 * t206 + 0.58904862254808623227e0 * t208 + 0.1e-18 * t211 - 0.18171526509422465567e2 * t214 + 0.14998437198861491884e0 * t216 - 0.29996874397722983769e0 * t218 + 0.58904862254808623227e0 * t220 - 0.11780972450961724646e1 * t222 + 0.1e-18 * t224 + 0.59993748795445967538e0 * t226 + 0.92537182407814154744e1 * t229 - 0.18171526509422465567e2 * t233 + 0.39995832530297311691e0 * t236;
445   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
446     v3rho3[0] = -0.29996874397722983769e0 * t100 - 0.26584037770686592174e0 * t84 + 0.11998749759089193508e1 * t104 - 0.11780972450961724646e1 * t102 + 0.92537182407814154744e1 * t108 + 0.11998749759089193507e1 * t114 - 0.13292018885343296087e0 * t117 - 0.53168075541373184349e0 * t127 - 0.29996874397722983769e0 * t133 + 0.10633615108274636870e1 * t87 + 0.10633615108274636870e1 * t98 - 0.47123889803846898582e1 * t111 - 0.10440614994059452211e1 * t120 + 0.12301341398300727784e2 * t125 + 0.21267230216549273739e1 * t129 + 0.21267230216549273739e1 * t131 - 0.11780972450961724646e1 * t135 + 0.92537182407814154744e1 * t138 + rho[0] * (t205 + t238);
447 
448 #ifndef XC_DONT_COMPILE_LXC
449 
450   if(order < 4) return;
451 
452 
453   t254 = 0.53168075541373184348e0 * t164 + 0.26584037770686592174e0 * t167 + 0.10633615108274636870e1 * t169 + 0.59993748795445967538e0 * t172 - 0.10633615108274636870e1 * t175 + 0.21267230216549273740e1 * t178 + 0.14178153477699515827e1 * t188 - 0.94247779607693797163e1 * t191 + 0.20881229988118904421e1 * t194 - 0.12883276855383007546e3 * t199 - 0.21267230216549273740e1 * t201 + 0.42534460433098547479e1 * t203 + 0.28356306955399031652e1 * t206;
454   t258 = t160 * rho[0];
455   t260 = 0.1e1 / t1 / t258;
456   t261 = t260 * t9;
457   t269 = t260 * t6;
458   t274 = t162 * t41;
459   t280 = 0.1e1 / t258;
460   t286 = t280 * t21;
461   t289 = t82 * t96;
462   t292 = t34 * t186;
463   t297 = t162 * t57;
464   t301 = t123 * t181;
465   t303 = -0.36819375000000000000e1 * t269 - 0.14460441337500000000e2 * t286 - 0.22716774923559000000e2 * t297 + 0.14015769222797481226e4 * t197 * t82 + 0.17843572366957123321e3 * t301;
466   t304 = t1 * t303;
467   t319 = -0.33230047213358240218e0 * t261 * t13 - 0.16615023606679120109e0 * t260 * t21 * t23 - 0.66460094426716480435e0 * t261 * t28 - 0.37496092997153729710e0 * t269 * t31 + 0.32208192138457518865e2 * t198 * t181 + 0.53168075541373184348e0 * t274 * t13 - 0.31415926535897932388e1 * t21 * t186 * t16 - 0.13050768742574315264e1 * t280 * t57 * t23 + 0.10633615108274636870e1 * t274 * t28 - 0.14726215563702155806e1 * t286 * t31 - 0.53168075541373184348e0 * t289 * t13 + 0.70890767388497579132e0 * t292 * t13 + 0.35445383694248789566e0 * t304 * t13 - 0.30753353495751819459e1 * t162 * t123 * t23 + 0.31623613451144514898e3 / t122 / t20 * t23 * t82 - 0.10633615108274636870e1 * t289 * t28 + 0.14178153477699515827e1 * t292 * t28;
468   t362 = 0.70890767388497579132e0 * t304 * t28 - 0.23134295601953538688e1 * t297 * t31 + 0.18171526509422465566e2 * t301 * t31 + 0.14273370642621158254e3 * t197 * t31 * t82 - 0.37496092997153729710e0 * t269 * t17 + 0.59993748795445967538e0 * t171 * t50 - 0.14726215563702155806e1 * t286 * t17 - 0.59993748795445967538e0 * t89 * t113 + 0.23561944901923449291e1 * t182 * t50 - 0.23134295601953538688e1 * t297 * t17 - 0.23561944901923449291e1 * t92 * t113 + 0.4e-18 * t210 * t50 + 0.18171526509422465566e2 * t301 * t17 + 0.79991665060594623384e0 * t38 * t235 + 0.18507436481562830949e2 * t57 * t96 * t107 - 0.72686106037689862268e2 * t123 * t41 * t232 + 0.14273370642621158254e3 * t197 * t9 * t16 * t82 + 0.39995832530297311691e0 * t7 * t303 * t16;
469   t374 = 0.23561944901923449292e1 * t208 + 0.4e-18 * t211 - 0.72686106037689862268e2 * t214 + rho[0] * (t319 + t362) + 0.59993748795445967538e0 * t216 - 0.11998749759089193508e1 * t218 + 0.23561944901923449292e1 * t220 - 0.47123889803846898584e1 * t222 + 0.4e-18 * t224 + 0.23997499518178387016e1 * t226 + 0.37014872963125661897e2 * t229 - 0.72686106037689862268e2 * t233 + 0.15998333012118924676e1 * t236;
470   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
471     v4rho4[0] = t254 + t374;
472 
473 #ifndef XC_DONT_COMPILE_MXC
474 
475   if(order < 5) return;
476 
477 
478 #endif
479 
480 #endif
481 
482 #endif
483 
484 #endif
485 
486 #endif
487 
488 
489 }
490 
491 
492 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))493 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
494 {
495 
496 #ifndef XC_DONT_COMPILE_EXC
497   double t1, t2, t4, t6, t7, t8, t10, t11;
498   double t12, t13, t14, t16, t17, t18, t20, t21;
499   double t22, t24, t26, t27, t28, t29, t31, t32;
500   double t34;
501 
502 #ifndef XC_DONT_COMPILE_VXC
503   double t35, t36, t37, t39, t42, t43, t44, t46;
504   double t49, t51, t52, t55, t58, t59, t61, t63;
505   double t65, t67, t69;
506 
507 #ifndef XC_DONT_COMPILE_FXC
508   double t83, t84, t85, t87, t88, t90, t92, t93;
509   double t97, t98, t99, t101, t103, t105, t108, t109;
510   double t112, t114, t115, t118, t121, t123, t124, t126;
511   double t128, t130, t132, t134, t136, t139, t141;
512 
513 #ifndef XC_DONT_COMPILE_KXC
514   double t161, t163, t164, t165, t168, t170, t172, t173;
515   double t175, t176, t178, t179, t182, t183, t187, t188;
516   double t189, t192, t195, t198, t199, t200, t202, t204;
517   double t206, t207, t209, t211, t212, t215, t217, t219;
518   double t221, t223, t225, t228, t231, t232, t234, t235;
519   double t237, t239;
520 
521 #ifndef XC_DONT_COMPILE_LXC
522   double t242, t245, t248, t250, t251, t253, t254, t256;
523   double t260, t262, t263, t287, t299, t302, t346, t361;
524   double t375;
525 #endif
526 
527 #endif
528 
529 #endif
530 
531 #endif
532 
533 #endif
534 
535 
536   lda_c_2d_prm_params *params;
537 
538   assert(p->params != NULL);
539   params = (lda_c_2d_prm_params * )(p->params);
540 
541   assert(params->N > 1);
542 
543   t1 = rho[0] + rho[1];
544   t2 = sqrt(t1);
545   t4 = sqrt(M_PI);
546   t6 = 0.39274000000000000001e1 * t2 + t4 / 0.2e1;
547   t7 = 0.1e1 / t6;
548   t8 = t2 * t7;
549   t10 = 0.39274000000000000001e1 * t8 - 0.1e1;
550   t11 = t2 * t10;
551   t12 = 0.2e1 + params->c;
552   t13 = sqrt(t12);
553   t14 = 0.1e1 / t13;
554   t16 = 0.35445383694248789566e0 * t11 * t14;
555   t17 = 0.1e1 / t12;
556   t18 = t10 * t17;
557   t20 = 0.39995832530297311691e0 * t8 * t18;
558   t21 = t6 * t6;
559   t22 = 0.1e1 / t21;
560   t24 = 0.1e1/POW_3_2(t12);
561   t26 = 0.17722691847124394783e0 * t2 * t22 * t24;
562   t27 = 0.1e1 + params->c;
563   t28 = sqrt(t27);
564   t29 = 0.1e1 / t28;
565   t31 = 0.70890767388497579132e0 * t11 * t29;
566   t32 = 0.1e1 / t27;
567   t34 = 0.39995832530297311691e0 * t8 * t32;
568   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
569     zk[0] = t16 + t20 + t26 + t31 + t34;
570 
571 #ifndef XC_DONT_COMPILE_VXC
572 
573   if(order < 1) return;
574 
575 
576   t35 = 0.1e1 / t2;
577   t36 = t35 * t10;
578   t37 = t36 * t14;
579   t39 = t35 * t7;
580   t42 = 0.19637000000000000000e1 * t39 - 0.77122353800000000002e1 * t22;
581   t43 = t2 * t42;
582   t44 = t43 * t14;
583   t46 = t39 * t18;
584   t49 = t22 * t10 * t17;
585   t51 = t42 * t17;
586   t52 = t8 * t51;
587   t55 = t35 * t22 * t24;
588   t58 = 0.1e1 / t21 / t6;
589   t59 = t58 * t24;
590   t61 = t36 * t29;
591   t63 = t43 * t29;
592   t65 = t39 * t32;
593   t67 = t22 * t32;
594   t69 = 0.17722691847124394783e0 * t37 + 0.35445383694248789566e0 * t44 + 0.19997916265148655846e0 * t46 - 0.78539816339744830968e0 * t49 + 0.39995832530297311691e0 * t52 + 0.88613459235621973915e-1 * t55 - 0.69604099960396348071e0 * t59 + 0.35445383694248789566e0 * t61 + 0.70890767388497579132e0 * t63 + 0.19997916265148655846e0 * t65 - 0.78539816339744830968e0 * t67;
595   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
596     vrho[0] = t1 * t69 + t16 + t20 + t26 + t31 + t34;
597 
598   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
599     vrho[1] = vrho[0];
600 
601 #ifndef XC_DONT_COMPILE_FXC
602 
603   if(order < 2) return;
604 
605 
606   t83 = 0.1e1 / t2 / t1;
607   t84 = t83 * t10;
608   t85 = t84 * t14;
609   t87 = t35 * t42;
610   t88 = t87 * t14;
611   t90 = t83 * t7;
612   t92 = 0.1e1 / t1;
613   t93 = t92 * t22;
614   t97 = -0.98185000000000000000e0 * t90 - 0.38561176900000000000e1 * t93 + 0.30289033231412000001e2 * t58 * t35;
615   t98 = t2 * t97;
616   t99 = t98 * t14;
617   t101 = t90 * t18;
618   t103 = t93 * t18;
619   t105 = t39 * t51;
620   t108 = t17 * t35;
621   t109 = t58 * t10 * t108;
622   t112 = t22 * t42 * t17;
623   t114 = t97 * t17;
624   t115 = t8 * t114;
625   t118 = t83 * t22 * t24;
626   t121 = t92 * t58 * t24;
627   t123 = t21 * t21;
628   t124 = 0.1e1 / t123;
629   t126 = t124 * t24 * t35;
630   t128 = t84 * t29;
631   t130 = t87 * t29;
632   t132 = t98 * t29;
633   t134 = t90 * t32;
634   t136 = t93 * t32;
635   t139 = t58 * t32 * t35;
636   t141 = -0.88613459235621973915e-1 * t85 + 0.35445383694248789566e0 * t88 + 0.35445383694248789566e0 * t99 - 0.99989581325743279230e-1 * t101 - 0.39269908169872415485e0 * t103 + 0.39995832530297311692e0 * t105 + 0.30845727469271384914e1 * t109 - 0.15707963267948966194e1 * t112 + 0.39995832530297311691e0 * t115 - 0.44306729617810986958e-1 * t118 - 0.34802049980198174035e0 * t121 + 0.41004471327669092612e1 * t126 - 0.17722691847124394783e0 * t128 + 0.70890767388497579132e0 * t130 + 0.70890767388497579132e0 * t132 - 0.99989581325743279230e-1 * t134 - 0.39269908169872415485e0 * t136 + 0.30845727469271384914e1 * t139;
637   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
638     v2rho2[0] = 0.35445383694248789566e0 * t37 + 0.70890767388497579132e0 * t44 + 0.39995832530297311692e0 * t46 - 0.15707963267948966194e1 * t49 + 0.79991665060594623382e0 * t52 + 0.17722691847124394783e0 * t55 - 0.13920819992079269614e1 * t59 + 0.70890767388497579132e0 * t61 + 0.14178153477699515826e1 * t63 + 0.39995832530297311692e0 * t65 - 0.15707963267948966194e1 * t67 + t1 * t141;
639 
640   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
641     v2rho2[1] = v2rho2[0];
642 
643   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
644     v2rho2[2] = v2rho2[1];
645 
646 #ifndef XC_DONT_COMPILE_KXC
647 
648   if(order < 3) return;
649 
650 
651   t161 = t1 * t1;
652   t163 = 0.1e1 / t2 / t161;
653   t164 = t163 * t10;
654   t165 = t164 * t14;
655   t168 = t163 * t22 * t24;
656   t170 = t164 * t29;
657   t172 = t163 * t7;
658   t173 = t172 * t32;
659   t175 = t83 * t42;
660   t176 = t175 * t14;
661   t178 = t35 * t97;
662   t179 = t178 * t14;
663   t182 = 0.1e1 / t161;
664   t183 = t182 * t22;
665   t187 = 0.14727750000000000000e1 * t172 + 0.57841765350000000000e1 * t183 - 0.17843572366957123321e3 * t124 * t92;
666   t188 = t2 * t187;
667   t189 = t188 * t14;
668   t192 = t22 * t97 * t17;
669   t195 = t182 * t58 * t24;
670   t198 = 0.1e1 / t123 / t6;
671   t199 = t198 * t24;
672   t200 = t199 * t92;
673   t202 = t175 * t29;
674   t204 = t178 * t29;
675   t206 = 0.13292018885343296087e0 * t165 + 0.66460094426716480437e-1 * t168 + 0.26584037770686592174e0 * t170 + 0.14998437198861491884e0 * t173 - 0.26584037770686592174e0 * t176 + 0.53168075541373184349e0 * t179 + 0.35445383694248789566e0 * t189 - 0.23561944901923449291e1 * t192 + 0.52203074970297261053e0 * t195 - 0.32208192138457518865e2 * t200 - 0.53168075541373184349e0 * t202 + 0.10633615108274636870e1 * t204;
676   t207 = t188 * t29;
677   t209 = t183 * t32;
678   t211 = t83 * t58;
679   t212 = t211 * t32;
680   t215 = t124 * t32 * t92;
681   t217 = t90 * t51;
682   t219 = t183 * t18;
683   t221 = t93 * t51;
684   t223 = t211 * t18;
685   t225 = t39 * t114;
686   t228 = t58 * t42 * t108;
687   t231 = t17 * t92;
688   t232 = t124 * t10 * t231;
689   t234 = t187 * t17;
690   t235 = t8 * t234;
691   t237 = t172 * t18;
692   t239 = 0.70890767388497579132e0 * t207 + 0.58904862254808623227e0 * t209 + 0.1e-18 * t212 - 0.18171526509422465567e2 * t215 - 0.29996874397722983769e0 * t217 + 0.58904862254808623227e0 * t219 - 0.11780972450961724646e1 * t221 + 0.1e-18 * t223 + 0.59993748795445967538e0 * t225 + 0.92537182407814154744e1 * t228 - 0.18171526509422465567e2 * t232 + 0.39995832530297311691e0 * t235 + 0.14998437198861491884e0 * t237;
693   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
694     v3rho3[0] = -0.29996874397722983769e0 * t101 - 0.26584037770686592174e0 * t85 + 0.11998749759089193508e1 * t105 - 0.11780972450961724646e1 * t103 + 0.92537182407814154744e1 * t109 + 0.11998749759089193507e1 * t115 - 0.13292018885343296087e0 * t118 - 0.53168075541373184349e0 * t128 - 0.29996874397722983769e0 * t134 + 0.10633615108274636870e1 * t88 + 0.10633615108274636870e1 * t99 - 0.47123889803846898582e1 * t112 - 0.10440614994059452211e1 * t121 + 0.12301341398300727784e2 * t126 + 0.21267230216549273739e1 * t130 + 0.21267230216549273739e1 * t132 - 0.11780972450961724646e1 * t136 + 0.92537182407814154744e1 * t139 + t1 * (t206 + t239);
695 
696   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
697     v3rho3[1] = v3rho3[0];
698 
699   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
700     v3rho3[2] = v3rho3[1];
701 
702   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
703     v3rho3[3] = v3rho3[2];
704 
705 #ifndef XC_DONT_COMPILE_LXC
706 
707   if(order < 4) return;
708 
709 
710   t242 = t83 * t97;
711   t245 = t35 * t187;
712   t248 = t161 * t1;
713   t250 = 0.1e1 / t2 / t248;
714   t251 = t250 * t7;
715   t253 = 0.1e1 / t248;
716   t254 = t253 * t22;
717   t256 = t163 * t58;
718   t260 = t124 * t182;
719   t262 = -0.36819375000000000000e1 * t251 - 0.14460441337500000000e2 * t254 - 0.22716774923559000000e2 * t256 + 0.14015769222797481226e4 * t198 * t83 + 0.17843572366957123321e3 * t260;
720   t263 = t2 * t262;
721   t287 = t250 * t10;
722   t299 = t163 * t42;
723   t302 = -0.53168075541373184348e0 * t242 * t14 + 0.70890767388497579132e0 * t245 * t14 + 0.35445383694248789566e0 * t263 * t14 - 0.30753353495751819459e1 * t163 * t124 * t24 + 0.31623613451144514898e3 / t123 / t21 * t24 * t83 - 0.10633615108274636870e1 * t242 * t29 + 0.14178153477699515827e1 * t245 * t29 + 0.70890767388497579132e0 * t263 * t29 - 0.23134295601953538688e1 * t256 * t32 + 0.18171526509422465566e2 * t260 * t32 + 0.14273370642621158254e3 * t198 * t32 * t83 - 0.33230047213358240218e0 * t287 * t14 - 0.16615023606679120109e0 * t250 * t22 * t24 - 0.66460094426716480435e0 * t287 * t29 - 0.37496092997153729710e0 * t251 * t32 + 0.32208192138457518865e2 * t199 * t182 + 0.53168075541373184348e0 * t299 * t14;
724   t346 = -0.31415926535897932388e1 * t22 * t187 * t17 - 0.13050768742574315264e1 * t253 * t58 * t24 + 0.10633615108274636870e1 * t299 * t29 - 0.14726215563702155806e1 * t254 * t32 + 0.59993748795445967538e0 * t172 * t51 - 0.14726215563702155806e1 * t254 * t18 - 0.59993748795445967538e0 * t90 * t114 + 0.23561944901923449291e1 * t183 * t51 - 0.23134295601953538688e1 * t256 * t18 - 0.23561944901923449291e1 * t93 * t114 + 0.4e-18 * t211 * t51 + 0.18171526509422465566e2 * t260 * t18 + 0.79991665060594623384e0 * t39 * t234 + 0.18507436481562830949e2 * t58 * t97 * t108 - 0.72686106037689862268e2 * t124 * t42 * t231 + 0.14273370642621158254e3 * t198 * t10 * t17 * t83 + 0.39995832530297311691e0 * t8 * t262 * t17 - 0.37496092997153729710e0 * t251 * t18;
725   t361 = t1 * (t302 + t346) + 0.53168075541373184348e0 * t165 + 0.26584037770686592174e0 * t168 + 0.10633615108274636870e1 * t170 + 0.59993748795445967538e0 * t173 - 0.10633615108274636870e1 * t176 + 0.21267230216549273740e1 * t179 + 0.14178153477699515827e1 * t189 - 0.94247779607693797163e1 * t192 + 0.20881229988118904421e1 * t195 - 0.12883276855383007546e3 * t200 - 0.21267230216549273740e1 * t202 + 0.42534460433098547479e1 * t204;
726   t375 = 0.28356306955399031652e1 * t207 + 0.23561944901923449292e1 * t209 + 0.4e-18 * t212 - 0.72686106037689862268e2 * t215 - 0.11998749759089193508e1 * t217 + 0.23561944901923449292e1 * t219 - 0.47123889803846898584e1 * t221 + 0.4e-18 * t223 + 0.23997499518178387016e1 * t225 + 0.37014872963125661897e2 * t228 - 0.72686106037689862268e2 * t232 + 0.15998333012118924676e1 * t235 + 0.59993748795445967538e0 * t237;
727   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
728     v4rho4[0] = t361 + t375;
729 
730   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
731     v4rho4[1] = v4rho4[0];
732 
733   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
734     v4rho4[2] = v4rho4[1];
735 
736   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
737     v4rho4[3] = v4rho4[2];
738 
739   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
740     v4rho4[4] = v4rho4[3];
741 
742 #ifndef XC_DONT_COMPILE_MXC
743 
744   if(order < 5) return;
745 
746 
747 #endif
748 
749 #endif
750 
751 #endif
752 
753 #endif
754 
755 #endif
756 
757 
758 }
759 
760