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