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_k_tf.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, t4, t5, t7, t8, t9, t10;
24   double t11, t13;
25 
26 #ifndef XC_DONT_COMPILE_LXC
27   double t22;
28 #endif
29 
30 #endif
31 
32 
33   lda_k_tf_params *params;
34 
35   assert(p->params != NULL);
36   params = (lda_k_tf_params * )(p->params);
37 
38   t1 = M_CBRT3;
39   t2 = params->ax * t1;
40   t4 = POW_1_3(0.1e1 / M_PI);
41   t5 = t4 * t4;
42   t7 = M_CBRT4;
43   t8 = t7 * t7;
44   t9 = 0.1e1 / t5 * t8;
45   t10 = POW_1_3(rho[0]);
46   t11 = t10 * t10;
47   t13 = t2 * t9 * t11;
48   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
49     zk[0] = t13 / 0.3e1;
50 
51 #ifndef XC_DONT_COMPILE_VXC
52 
53   if(order < 1) return;
54 
55 
56   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
57     vrho[0] = 0.5e1 / 0.9e1 * t13;
58 
59 #ifndef XC_DONT_COMPILE_FXC
60 
61   if(order < 2) return;
62 
63 
64   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
65     v2rho2[0] = 0.10e2 / 0.27e2 * t2 * t9 / t10;
66 
67 #ifndef XC_DONT_COMPILE_KXC
68 
69   if(order < 3) return;
70 
71 
72   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
73     v3rho3[0] = -0.10e2 / 0.81e2 * t2 * t9 / t10 / rho[0];
74 
75 #ifndef XC_DONT_COMPILE_LXC
76 
77   if(order < 4) return;
78 
79 
80   t22 = rho[0] * rho[0];
81   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
82     v4rho4[0] = 0.40e2 / 0.243e3 * t2 * t9 / t10 / t22;
83 
84 #ifndef XC_DONT_COMPILE_MXC
85 
86   if(order < 5) return;
87 
88 
89 #endif
90 
91 #endif
92 
93 #endif
94 
95 #endif
96 
97 #endif
98 
99 
100 }
101 
102 
103 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))104 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
105 {
106 
107 #ifndef XC_DONT_COMPILE_EXC
108   double t1, t2, t4, t5, t7, t8, t10, t11;
109   double t12, t13, t14, t16;
110 
111 #ifndef XC_DONT_COMPILE_LXC
112   double t25;
113 #endif
114 
115 #endif
116 
117 
118   lda_k_tf_params *params;
119 
120   assert(p->params != NULL);
121   params = (lda_k_tf_params * )(p->params);
122 
123   t1 = M_CBRT2;
124   t2 = t1 * t1;
125   t4 = M_CBRT3;
126   t5 = params->ax * t2 * t4;
127   t7 = POW_1_3(0.1e1 / M_PI);
128   t8 = t7 * t7;
129   t10 = M_CBRT4;
130   t11 = t10 * t10;
131   t12 = 0.1e1 / t8 * t11;
132   t13 = POW_1_3(rho[0]);
133   t14 = t13 * t13;
134   t16 = t5 * t12 * t14;
135   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
136     zk[0] = t16 / 0.3e1;
137 
138 #ifndef XC_DONT_COMPILE_VXC
139 
140   if(order < 1) return;
141 
142 
143   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
144     vrho[0] = 0.5e1 / 0.9e1 * t16;
145 
146 #ifndef XC_DONT_COMPILE_FXC
147 
148   if(order < 2) return;
149 
150 
151   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
152     v2rho2[0] = 0.10e2 / 0.27e2 * t5 * t12 / t13;
153 
154 #ifndef XC_DONT_COMPILE_KXC
155 
156   if(order < 3) return;
157 
158 
159   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
160     v3rho3[0] = -0.10e2 / 0.81e2 * t5 * t12 / t13 / rho[0];
161 
162 #ifndef XC_DONT_COMPILE_LXC
163 
164   if(order < 4) return;
165 
166 
167   t25 = rho[0] * rho[0];
168   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
169     v4rho4[0] = 0.40e2 / 0.243e3 * t5 * t12 / t13 / t25;
170 
171 #ifndef XC_DONT_COMPILE_MXC
172 
173   if(order < 5) return;
174 
175 
176 #endif
177 
178 #endif
179 
180 #endif
181 
182 #endif
183 
184 #endif
185 
186 
187 }
188 
189 
190 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))191 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
192 {
193 
194 #ifndef XC_DONT_COMPILE_EXC
195   double t1, t2, t3, t4, t5, t6, t7, t9;
196   double t10, t11, t16, t17, t19, t20, t21, t22;
197   double t23, t24, t25, t26, t27, t28;
198 
199 #ifndef XC_DONT_COMPILE_VXC
200   double t29, t31, t32, t33, t34, t35, t37, t40;
201   double t43, t46, t48, t51;
202 
203 #ifndef XC_DONT_COMPILE_FXC
204   double t56, t57, t59, t60, t62, t63, t64, t68;
205   double t69, t71, t74, t75, t78, t81, t86, t88;
206   double t90, t93, t96, t99, t102, t107, t111, t114;
207   double t117, t120;
208 
209 #ifndef XC_DONT_COMPILE_KXC
210   double t125, t126, t128, t131, t132, t134, t136, t137;
211   double t140, t143, t144, t145, t147, t151, t152, t155;
212   double t158, t161, t167, t169, t173, t174, t177, t187;
213   double t190, t200, t207, t209, t214, t219, t222, t227;
214   double t230, t233, t239, t245, t248, t253, t256;
215 
216 #ifndef XC_DONT_COMPILE_LXC
217   double t262, t264, t266, t272, t273, t275, t276, t282;
218   double t288, t289, t291, t294, t296, t297, t303, t320;
219   double t323, t324, t326, t351, t369, t378, t388, t391;
220   double t402, t405, t444, t453, t473, t495, t502, t507;
221   double t513, t516, t521;
222 #endif
223 
224 #endif
225 
226 #endif
227 
228 #endif
229 
230 #endif
231 
232 
233   lda_k_tf_params *params;
234 
235   assert(p->params != NULL);
236   params = (lda_k_tf_params * )(p->params);
237 
238   t1 = rho[0] - rho[1];
239   t2 = rho[0] + rho[1];
240   t3 = 0.1e1 / t2;
241   t4 = t1 * t3;
242   t5 = 0.1e1 + t4;
243   t6 = POW_1_3(t5);
244   t7 = t6 * t6;
245   t9 = 0.1e1 - t4;
246   t10 = POW_1_3(t9);
247   t11 = t10 * t10;
248   t16 = M_CBRT3;
249   t17 = params->ax * (t11 * t9 / 0.2e1 + t7 * t5 / 0.2e1) * t16;
250   t19 = POW_1_3(0.1e1 / M_PI);
251   t20 = t19 * t19;
252   t21 = 0.1e1 / t20;
253   t22 = M_CBRT4;
254   t23 = t22 * t22;
255   t24 = t21 * t23;
256   t25 = POW_1_3(t2);
257   t26 = t25 * t25;
258   t27 = t24 * t26;
259   t28 = t17 * t27;
260   if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
261     zk[0] = t28 / 0.3e1;
262 
263 #ifndef XC_DONT_COMPILE_VXC
264 
265   if(order < 1) return;
266 
267 
268   t29 = 0.5e1 / 0.9e1 * t28;
269   t31 = t26 * t2 * params->ax;
270   t32 = t2 * t2;
271   t33 = 0.1e1 / t32;
272   t34 = t1 * t33;
273   t35 = t3 - t34;
274   t37 = -t35;
275   t40 = 0.5e1 / 0.6e1 * t11 * t37 + 0.5e1 / 0.6e1 * t7 * t35;
276   t43 = t16 * t21 * t23;
277   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
278     vrho[0] = t29 + t31 * t40 * t43 / 0.3e1;
279 
280   t46 = -t3 - t34;
281   t48 = -t46;
282   t51 = 0.5e1 / 0.6e1 * t11 * t48 + 0.5e1 / 0.6e1 * t7 * t46;
283   if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
284     vrho[1] = t29 + t31 * t51 * t43 / 0.3e1;
285 
286 #ifndef XC_DONT_COMPILE_FXC
287 
288   if(order < 2) return;
289 
290 
291   t56 = params->ax * t40 * t16;
292   t57 = t56 * t27;
293   t59 = 0.1e1 / t25;
294   t60 = t24 * t59;
295   t62 = 0.10e2 / 0.27e2 * t17 * t60;
296   t63 = 0.1e1 / t6;
297   t64 = t35 * t35;
298   t68 = 0.1e1 / t32 / t2;
299   t69 = t1 * t68;
300   t71 = -0.2e1 * t33 + 0.2e1 * t69;
301   t74 = 0.1e1 / t10;
302   t75 = t37 * t37;
303   t78 = -t71;
304   t81 = 0.5e1 / 0.9e1 * t63 * t64 + 0.5e1 / 0.6e1 * t7 * t71 + 0.5e1 / 0.9e1 * t74 * t75 + 0.5e1 / 0.6e1 * t11 * t78;
305   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
306     v2rho2[0] = 0.10e2 / 0.9e1 * t57 + t62 + t31 * t81 * t43 / 0.3e1;
307 
308   t86 = t26 * params->ax;
309   t88 = t86 * t51 * t43;
310   t90 = t63 * t46;
311   t93 = t7 * t1;
312   t96 = t74 * t48;
313   t99 = t11 * t1;
314   t102 = 0.5e1 / 0.9e1 * t90 * t35 + 0.5e1 / 0.3e1 * t93 * t68 + 0.5e1 / 0.9e1 * t96 * t37 - 0.5e1 / 0.3e1 * t99 * t68;
315   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
316     v2rho2[1] = 0.5e1 / 0.9e1 * t57 + t62 + 0.5e1 / 0.9e1 * t88 + t31 * t102 * t43 / 0.3e1;
317 
318   t107 = t46 * t46;
319   t111 = 0.2e1 * t33 + 0.2e1 * t69;
320   t114 = t48 * t48;
321   t117 = -t111;
322   t120 = 0.5e1 / 0.9e1 * t63 * t107 + 0.5e1 / 0.6e1 * t7 * t111 + 0.5e1 / 0.9e1 * t74 * t114 + 0.5e1 / 0.6e1 * t11 * t117;
323   if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
324     v2rho2[2] = 0.10e2 / 0.9e1 * t88 + t62 + t31 * t120 * t43 / 0.3e1;
325 
326 #ifndef XC_DONT_COMPILE_KXC
327 
328   if(order < 3) return;
329 
330 
331   t125 = params->ax * t81 * t16;
332   t126 = t125 * t27;
333   t128 = t56 * t60;
334   t131 = 0.1e1 / t25 / t2;
335   t132 = t24 * t131;
336   t134 = 0.10e2 / 0.81e2 * t17 * t132;
337   t136 = 0.1e1 / t6 / t5;
338   t137 = t64 * t35;
339   t140 = t63 * t35;
340   t143 = t32 * t32;
341   t144 = 0.1e1 / t143;
342   t145 = t1 * t144;
343   t147 = 0.6e1 * t68 - 0.6e1 * t145;
344   t151 = 0.1e1 / t10 / t9;
345   t152 = t75 * t37;
346   t155 = t74 * t37;
347   t158 = -t147;
348   t161 = -0.5e1 / 0.27e2 * t136 * t137 + 0.5e1 / 0.3e1 * t140 * t71 + 0.5e1 / 0.6e1 * t7 * t147 - 0.5e1 / 0.27e2 * t151 * t152 + 0.5e1 / 0.3e1 * t155 * t78 + 0.5e1 / 0.6e1 * t11 * t158;
349   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
350     v3rho3[0] = 0.5e1 / 0.3e1 * t126 + 0.10e2 / 0.9e1 * t128 - t134 + t31 * t161 * t43 / 0.3e1;
351 
352   t167 = t59 * params->ax;
353   t169 = t167 * t51 * t43;
354   t173 = 0.10e2 / 0.9e1 * t86 * t102 * t43;
355   t174 = t136 * t46;
356   t177 = t63 * t1;
357   t187 = t151 * t48;
358   t190 = t74 * t1;
359   t200 = -0.5e1 / 0.27e2 * t174 * t64 + 0.20e2 / 0.9e1 * t177 * t68 * t35 + 0.5e1 / 0.9e1 * t90 * t71 + 0.5e1 / 0.3e1 * t7 * t68 - 0.5e1 * t93 * t144 - 0.5e1 / 0.27e2 * t187 * t75 - 0.20e2 / 0.9e1 * t190 * t68 * t37 + 0.5e1 / 0.9e1 * t96 * t78 - 0.5e1 / 0.3e1 * t11 * t68 + 0.5e1 * t99 * t144;
360   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
361     v3rho3[1] = 0.5e1 / 0.9e1 * t126 + 0.20e2 / 0.27e2 * t128 - t134 + 0.10e2 / 0.27e2 * t169 + t173 + t31 * t200 * t43 / 0.3e1;
362 
363   t207 = t86 * t120 * t43;
364   t209 = t136 * t107;
365   t214 = t63 * t111;
366   t219 = -0.2e1 * t68 - 0.6e1 * t145;
367   t222 = t151 * t114;
368   t227 = t74 * t117;
369   t230 = -t219;
370   t233 = -0.5e1 / 0.27e2 * t209 * t35 + 0.20e2 / 0.9e1 * t90 * t69 + 0.5e1 / 0.9e1 * t214 * t35 + 0.5e1 / 0.6e1 * t7 * t219 - 0.5e1 / 0.27e2 * t222 * t37 - 0.20e2 / 0.9e1 * t96 * t69 + 0.5e1 / 0.9e1 * t227 * t37 + 0.5e1 / 0.6e1 * t11 * t230;
371   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
372     v3rho3[2] = 0.20e2 / 0.27e2 * t169 + t173 + 0.10e2 / 0.27e2 * t128 - t134 + 0.5e1 / 0.9e1 * t207 + t31 * t233 * t43 / 0.3e1;
373 
374   t239 = t107 * t46;
375   t245 = -0.6e1 * t68 - 0.6e1 * t145;
376   t248 = t114 * t48;
377   t253 = -t245;
378   t256 = -0.5e1 / 0.27e2 * t136 * t239 + 0.5e1 / 0.3e1 * t90 * t111 + 0.5e1 / 0.6e1 * t7 * t245 - 0.5e1 / 0.27e2 * t151 * t248 + 0.5e1 / 0.3e1 * t96 * t117 + 0.5e1 / 0.6e1 * t11 * t253;
379   if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
380     v3rho3[3] = 0.10e2 / 0.9e1 * t169 + 0.5e1 / 0.3e1 * t207 - t134 + t31 * t256 * t43 / 0.3e1;
381 
382 #ifndef XC_DONT_COMPILE_LXC
383 
384   if(order < 4) return;
385 
386 
387   t262 = params->ax * t161 * t16 * t27;
388   t264 = t125 * t60;
389   t266 = t56 * t132;
390   t272 = 0.40e2 / 0.243e3 * t17 * t24 / t25 / t32;
391   t273 = t5 * t5;
392   t275 = 0.1e1 / t6 / t273;
393   t276 = t64 * t64;
394   t282 = t71 * t71;
395   t288 = 0.1e1 / t143 / t2;
396   t289 = t1 * t288;
397   t291 = -0.24e2 * t144 + 0.24e2 * t289;
398   t294 = t9 * t9;
399   t296 = 0.1e1 / t10 / t294;
400   t297 = t75 * t75;
401   t303 = t78 * t78;
402   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
403     v4rho4[0] = 0.20e2 / 0.9e1 * t262 + 0.20e2 / 0.9e1 * t264 - 0.40e2 / 0.81e2 * t266 + t272 + t31 * (0.20e2 / 0.81e2 * t275 * t276 - 0.10e2 / 0.9e1 * t136 * t64 * t71 + 0.5e1 / 0.3e1 * t63 * t282 + 0.20e2 / 0.9e1 * t140 * t147 + 0.5e1 / 0.6e1 * t7 * t291 + 0.20e2 / 0.81e2 * t296 * t297 - 0.10e2 / 0.9e1 * t151 * t75 * t78 + 0.5e1 / 0.3e1 * t74 * t303 + 0.20e2 / 0.9e1 * t155 * t158 - 0.5e1 / 0.6e1 * t11 * t291) * t43 / 0.3e1;
404 
405   t320 = t131 * params->ax * t51 * t43;
406   t323 = t167 * t102 * t43;
407   t324 = 0.10e2 / 0.9e1 * t323;
408   t326 = t86 * t200 * t43;
409   t351 = 0.20e2 * t93 * t288;
410   t369 = 0.20e2 * t99 * t288;
411   t378 = -0.10e2 * t177 * t144 * t35 + 0.10e2 * t190 * t144 * t37 + 0.20e2 / 0.81e2 * t275 * t46 * t137 - 0.5e1 / 0.9e1 * t174 * t35 * t71 - 0.10e2 / 0.9e1 * t136 * t1 * t68 * t64 + 0.10e2 / 0.3e1 * t63 * t68 * t35 + 0.10e2 / 0.3e1 * t177 * t68 * t71 + t351 + 0.20e2 / 0.81e2 * t296 * t48 * t152 - 0.5e1 / 0.9e1 * t187 * t37 * t78 + 0.10e2 / 0.9e1 * t151 * t1 * t68 * t75 - 0.10e2 / 0.3e1 * t74 * t68 * t37 - 0.10e2 / 0.3e1 * t190 * t68 * t78 - t369 + 0.5e1 / 0.9e1 * t90 * t147 - 0.10e2 * t7 * t144 + 0.5e1 / 0.9e1 * t96 * t158 + 0.10e2 * t11 * t144;
412   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
413     v4rho4[1] = 0.5e1 / 0.9e1 * t262 + 0.10e2 / 0.9e1 * t264 - 0.10e2 / 0.27e2 * t266 + t272 - 0.10e2 / 0.81e2 * t320 + t324 + 0.5e1 / 0.3e1 * t326 + t31 * t378 * t43 / 0.3e1;
414 
415   t388 = t167 * t120 * t43;
416   t391 = t86 * t233 * t43;
417   t402 = t1 * t1;
418   t405 = 0.1e1 / t143 / t32;
419   t444 = 0.20e2 / 0.81e2 * t275 * t107 * t64 - 0.40e2 / 0.27e2 * t174 * t35 * t1 * t68 - 0.5e1 / 0.27e2 * t209 * t71 + 0.40e2 / 0.9e1 * t63 * t402 * t405 + 0.20e2 / 0.9e1 * t90 * t68 - 0.20e2 / 0.3e1 * t90 * t145 - 0.5e1 / 0.27e2 * t136 * t111 * t64 + 0.10e2 / 0.9e1 * t63 * t219 * t35 + 0.5e1 / 0.9e1 * t214 * t71 + t351 + 0.20e2 / 0.81e2 * t296 * t114 * t75 + 0.40e2 / 0.27e2 * t187 * t37 * t1 * t68 - 0.5e1 / 0.27e2 * t222 * t78 + 0.40e2 / 0.9e1 * t74 * t402 * t405 - 0.20e2 / 0.9e1 * t96 * t68 + 0.20e2 / 0.3e1 * t96 * t145 - 0.5e1 / 0.27e2 * t151 * t117 * t75 + 0.10e2 / 0.9e1 * t74 * t230 * t37 + 0.5e1 / 0.9e1 * t227 * t78 - t369;
420   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
421     v4rho4[2] = -0.20e2 / 0.81e2 * t320 + 0.40e2 / 0.27e2 * t323 + 0.10e2 / 0.9e1 * t326 + 0.10e2 / 0.27e2 * t264 - 0.20e2 / 0.81e2 * t266 + t272 + 0.10e2 / 0.27e2 * t388 + 0.10e2 / 0.9e1 * t391 + t31 * t444 * t43 / 0.3e1;
422 
423   t453 = t86 * t256 * t43;
424   t473 = 0.12e2 * t144 + 0.24e2 * t289;
425   t495 = 0.20e2 / 0.81e2 * t275 * t239 * t35 - 0.10e2 / 0.9e1 * t209 * t69 - 0.5e1 / 0.9e1 * t174 * t111 * t35 + 0.10e2 / 0.3e1 * t177 * t68 * t111 + 0.5e1 / 0.3e1 * t90 * t219 + 0.5e1 / 0.9e1 * t63 * t245 * t35 + 0.5e1 / 0.6e1 * t7 * t473 + 0.20e2 / 0.81e2 * t296 * t248 * t37 + 0.10e2 / 0.9e1 * t222 * t69 - 0.5e1 / 0.9e1 * t187 * t117 * t37 - 0.10e2 / 0.3e1 * t190 * t68 * t117 + 0.5e1 / 0.3e1 * t96 * t230 + 0.5e1 / 0.9e1 * t74 * t253 * t37 - 0.5e1 / 0.6e1 * t11 * t473;
426   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
427     v4rho4[3] = -0.10e2 / 0.27e2 * t320 + t324 + 0.10e2 / 0.9e1 * t388 + 0.5e1 / 0.3e1 * t391 - 0.10e2 / 0.81e2 * t266 + t272 + 0.5e1 / 0.9e1 * t453 + t31 * t495 * t43 / 0.3e1;
428 
429   t502 = t107 * t107;
430   t507 = t111 * t111;
431   t513 = 0.24e2 * t144 + 0.24e2 * t289;
432   t516 = t114 * t114;
433   t521 = t117 * t117;
434   if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
435     v4rho4[4] = -0.40e2 / 0.81e2 * t320 + 0.20e2 / 0.9e1 * t388 + 0.20e2 / 0.9e1 * t453 + t272 + t31 * (0.20e2 / 0.81e2 * t275 * t502 - 0.10e2 / 0.9e1 * t209 * t111 + 0.5e1 / 0.3e1 * t63 * t507 + 0.20e2 / 0.9e1 * t90 * t245 + 0.5e1 / 0.6e1 * t7 * t513 + 0.20e2 / 0.81e2 * t296 * t516 - 0.10e2 / 0.9e1 * t222 * t117 + 0.5e1 / 0.3e1 * t74 * t521 + 0.20e2 / 0.9e1 * t96 * t253 - 0.5e1 / 0.6e1 * t11 * t513) * t43 / 0.3e1;
436 
437 #ifndef XC_DONT_COMPILE_MXC
438 
439   if(order < 5) return;
440 
441 
442 #endif
443 
444 #endif
445 
446 #endif
447 
448 #endif
449 
450 #endif
451 
452 
453 }
454 
455