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_x.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, t4, t5, t6, t7, t8, t9, t10;
24 double t11, t13;
25
26 #ifndef XC_DONT_COMPILE_FXC
27 double t16;
28
29 #ifndef XC_DONT_COMPILE_LXC
30 double t25;
31 #endif
32
33 #endif
34
35 #endif
36
37
38 lda_x_params *params;
39
40 assert(p->params != NULL);
41 params = (lda_x_params * )(p->params);
42
43 t1 = M_CBRT3;
44 t4 = POW_1_3(0.1e1 / M_PI);
45 t5 = params->alpha * t1 * t4;
46 t6 = M_CBRT4;
47 t7 = t6 * t6;
48 t8 = M_CBRT2;
49 t9 = t8 * t8;
50 t10 = t7 * t9;
51 t11 = POW_1_3(rho[0]);
52 t13 = t5 * t10 * t11;
53 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
54 zk[0] = -0.3e1 / 0.16e2 * t13;
55
56 #ifndef XC_DONT_COMPILE_VXC
57
58 if(order < 1) return;
59
60
61 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
62 vrho[0] = -t13 / 0.4e1;
63
64 #ifndef XC_DONT_COMPILE_FXC
65
66 if(order < 2) return;
67
68
69 t16 = t11 * t11;
70 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
71 v2rho2[0] = -t5 * t10 / t16 / 0.12e2;
72
73 #ifndef XC_DONT_COMPILE_KXC
74
75 if(order < 3) return;
76
77
78 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
79 v3rho3[0] = t5 * t10 / t16 / rho[0] / 0.18e2;
80
81 #ifndef XC_DONT_COMPILE_LXC
82
83 if(order < 4) return;
84
85
86 t25 = rho[0] * rho[0];
87 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
88 v4rho4[0] = -0.5e1 / 0.54e2 * t5 * t10 / t16 / t25;
89
90 #ifndef XC_DONT_COMPILE_MXC
91
92 if(order < 5) return;
93
94
95 #endif
96
97 #endif
98
99 #endif
100
101 #endif
102
103 #endif
104
105
106 }
107
108
109 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))110 func_ferr(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
111 {
112
113 #ifndef XC_DONT_COMPILE_EXC
114 double t1, t2, t4, t5, t6, t7, t8, t10;
115
116 #ifndef XC_DONT_COMPILE_FXC
117 double t13;
118
119 #ifndef XC_DONT_COMPILE_LXC
120 double t22;
121 #endif
122
123 #endif
124
125 #endif
126
127
128 lda_x_params *params;
129
130 assert(p->params != NULL);
131 params = (lda_x_params * )(p->params);
132
133 t1 = M_CBRT3;
134 t2 = params->alpha * t1;
135 t4 = POW_1_3(0.1e1 / M_PI);
136 t5 = M_CBRT4;
137 t6 = t5 * t5;
138 t7 = t4 * t6;
139 t8 = POW_1_3(rho[0]);
140 t10 = t2 * t7 * t8;
141 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
142 zk[0] = -0.3e1 / 0.8e1 * t10;
143
144 #ifndef XC_DONT_COMPILE_VXC
145
146 if(order < 1) return;
147
148
149 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
150 vrho[0] = -t10 / 0.2e1;
151
152 #ifndef XC_DONT_COMPILE_FXC
153
154 if(order < 2) return;
155
156
157 t13 = t8 * t8;
158 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
159 v2rho2[0] = -t2 * t7 / t13 / 0.6e1;
160
161 #ifndef XC_DONT_COMPILE_KXC
162
163 if(order < 3) return;
164
165
166 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
167 v3rho3[0] = t2 * t7 / t13 / rho[0] / 0.9e1;
168
169 #ifndef XC_DONT_COMPILE_LXC
170
171 if(order < 4) return;
172
173
174 t22 = rho[0] * rho[0];
175 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
176 v4rho4[0] = -0.5e1 / 0.27e2 * t2 * t7 / t13 / t22;
177
178 #ifndef XC_DONT_COMPILE_MXC
179
180 if(order < 5) return;
181
182
183 #endif
184
185 #endif
186
187 #endif
188
189 #endif
190
191 #endif
192
193
194 }
195
196
197 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,double * zk,LDA_OUT_PARAMS_NO_EXC (double *))198 func_pol(const xc_func_type *p, int order, const double *rho, double *zk, LDA_OUT_PARAMS_NO_EXC(double *))
199 {
200
201 #ifndef XC_DONT_COMPILE_EXC
202 double t1, t4, t5, t6, t7, t8, t9, t10;
203 double t11, t12, t13, t14, t15, t16, t18, t19;
204 double t21, t22, t25;
205
206 #ifndef XC_DONT_COMPILE_VXC
207 double t27, t30, t31, t32, t33, t34, t35, t37;
208 double t40, t45, t47, t52;
209
210 #ifndef XC_DONT_COMPILE_FXC
211 double t57, t59, t60, t64, t65, t66, t67, t71;
212 double t72, t74, t77, t78, t79, t82, t85, t92;
213 double t93, t95, t98, t101, t104, t109, t113, t117;
214 double t120, t123, t128;
215
216 #ifndef XC_DONT_COMPILE_KXC
217 double t133, t137, t140, t144, t146, t147, t150, t153;
218 double t154, t155, t157, t161, t162, t165, t168, t171;
219 double t179, t180, t183, t184, t187, t197, t200, t212;
220 double t217, t219, t224, t229, t232, t237, t240, t245;
221 double t250, t256, t259, t264, t269;
222
223 #ifndef XC_DONT_COMPILE_LXC
224 double t274, t278, t282, t289, t290, t292, t293, t299;
225 double t305, t306, t308, t311, t313, t314, t320, t338;
226 double t340, t341, t342, t367, t385, t394, t404, t406;
227 double t417, t420, t459, t468, t488, t510, t518, t523;
228 double t529, t532, t537;
229 #endif
230
231 #endif
232
233 #endif
234
235 #endif
236
237 #endif
238
239
240 lda_x_params *params;
241
242 assert(p->params != NULL);
243 params = (lda_x_params * )(p->params);
244
245 t1 = M_CBRT3;
246 t4 = POW_1_3(0.1e1 / M_PI);
247 t5 = params->alpha * t1 * t4;
248 t6 = M_CBRT4;
249 t7 = t6 * t6;
250 t8 = M_CBRT2;
251 t9 = t8 * t8;
252 t10 = t7 * t9;
253 t11 = rho[0] - rho[1];
254 t12 = rho[0] + rho[1];
255 t13 = 0.1e1 / t12;
256 t14 = t11 * t13;
257 t15 = 0.1e1 + t14;
258 t16 = POW_1_3(t15);
259 t18 = 0.1e1 - t14;
260 t19 = POW_1_3(t18);
261 t21 = t16 * t15 + t19 * t18;
262 t22 = POW_1_3(t12);
263 t25 = t5 * t10 * t21 * t22;
264 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
265 zk[0] = -0.3e1 / 0.32e2 * t25;
266
267 #ifndef XC_DONT_COMPILE_VXC
268
269 if(order < 1) return;
270
271
272 t27 = t25 / 0.8e1;
273 t30 = t22 * t12 * params->alpha * t1;
274 t31 = t4 * t7;
275 t32 = t12 * t12;
276 t33 = 0.1e1 / t32;
277 t34 = t11 * t33;
278 t35 = t13 - t34;
279 t37 = -t35;
280 t40 = 0.4e1 / 0.3e1 * t16 * t35 + 0.4e1 / 0.3e1 * t19 * t37;
281 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
282 vrho[0] = -t27 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t40;
283
284 t45 = -t13 - t34;
285 t47 = -t45;
286 t52 = t31 * t9 * (0.4e1 / 0.3e1 * t16 * t45 + 0.4e1 / 0.3e1 * t19 * t47);
287 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
288 vrho[1] = -t27 - 0.3e1 / 0.32e2 * t30 * t52;
289
290 #ifndef XC_DONT_COMPILE_FXC
291
292 if(order < 2) return;
293
294
295 t57 = t5 * t10 * t40 * t22;
296 t59 = t22 * t22;
297 t60 = 0.1e1 / t59;
298 t64 = t5 * t10 * t21 * t60 / 0.24e2;
299 t65 = t16 * t16;
300 t66 = 0.1e1 / t65;
301 t67 = t35 * t35;
302 t71 = 0.1e1 / t32 / t12;
303 t72 = t11 * t71;
304 t74 = -0.2e1 * t33 + 0.2e1 * t72;
305 t77 = t19 * t19;
306 t78 = 0.1e1 / t77;
307 t79 = t37 * t37;
308 t82 = -t74;
309 t85 = 0.4e1 / 0.9e1 * t66 * t67 + 0.4e1 / 0.3e1 * t16 * t74 + 0.4e1 / 0.9e1 * t78 * t79 + 0.4e1 / 0.3e1 * t19 * t82;
310 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
311 v2rho2[0] = -t57 / 0.4e1 - t64 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t85;
312
313 t92 = t22 * params->alpha * t1;
314 t93 = t92 * t52;
315 t95 = t66 * t45;
316 t98 = t16 * t11;
317 t101 = t78 * t47;
318 t104 = t19 * t11;
319 t109 = t31 * t9 * (0.4e1 / 0.9e1 * t95 * t35 + 0.8e1 / 0.3e1 * t98 * t71 + 0.4e1 / 0.9e1 * t101 * t37 - 0.8e1 / 0.3e1 * t104 * t71);
320 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
321 v2rho2[1] = -t57 / 0.8e1 - t64 - t93 / 0.8e1 - 0.3e1 / 0.32e2 * t30 * t109;
322
323 t113 = t45 * t45;
324 t117 = 0.2e1 * t33 + 0.2e1 * t72;
325 t120 = t47 * t47;
326 t123 = -t117;
327 t128 = t31 * t9 * (0.4e1 / 0.9e1 * t66 * t113 + 0.4e1 / 0.3e1 * t16 * t117 + 0.4e1 / 0.9e1 * t78 * t120 + 0.4e1 / 0.3e1 * t19 * t123);
328 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
329 v2rho2[2] = -t93 / 0.4e1 - t64 - 0.3e1 / 0.32e2 * t30 * t128;
330
331 #ifndef XC_DONT_COMPILE_KXC
332
333 if(order < 3) return;
334
335
336 t133 = t5 * t10 * t85 * t22;
337 t137 = t5 * t10 * t40 * t60;
338 t140 = 0.1e1 / t59 / t12;
339 t144 = t5 * t10 * t21 * t140 / 0.36e2;
340 t146 = 0.1e1 / t65 / t15;
341 t147 = t67 * t35;
342 t150 = t66 * t35;
343 t153 = t32 * t32;
344 t154 = 0.1e1 / t153;
345 t155 = t11 * t154;
346 t157 = 0.6e1 * t71 - 0.6e1 * t155;
347 t161 = 0.1e1 / t77 / t18;
348 t162 = t79 * t37;
349 t165 = t78 * t37;
350 t168 = -t157;
351 t171 = -0.8e1 / 0.27e2 * t146 * t147 + 0.4e1 / 0.3e1 * t150 * t74 + 0.4e1 / 0.3e1 * t16 * t157 - 0.8e1 / 0.27e2 * t161 * t162 + 0.4e1 / 0.3e1 * t165 * t82 + 0.4e1 / 0.3e1 * t19 * t168;
352 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
353 v3rho3[0] = -0.3e1 / 0.8e1 * t133 - t137 / 0.8e1 + t144 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t171;
354
355 t179 = t60 * params->alpha * t1;
356 t180 = t179 * t52;
357 t183 = t92 * t109 / 0.4e1;
358 t184 = t146 * t45;
359 t187 = t66 * t11;
360 t197 = t161 * t47;
361 t200 = t78 * t11;
362 t212 = t31 * t9 * (-0.8e1 / 0.27e2 * t184 * t67 + 0.16e2 / 0.9e1 * t187 * t71 * t35 + 0.4e1 / 0.9e1 * t95 * t74 + 0.8e1 / 0.3e1 * t16 * t71 - 0.8e1 * t98 * t154 - 0.8e1 / 0.27e2 * t197 * t79 - 0.16e2 / 0.9e1 * t200 * t71 * t37 + 0.4e1 / 0.9e1 * t101 * t82 - 0.8e1 / 0.3e1 * t19 * t71 + 0.8e1 * t104 * t154);
363 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
364 v3rho3[1] = -t133 / 0.8e1 - t137 / 0.12e2 + t144 - t180 / 0.24e2 - t183 - 0.3e1 / 0.32e2 * t30 * t212;
365
366 t217 = t92 * t128;
367 t219 = t146 * t113;
368 t224 = t66 * t117;
369 t229 = -0.2e1 * t71 - 0.6e1 * t155;
370 t232 = t161 * t120;
371 t237 = t78 * t123;
372 t240 = -t229;
373 t245 = t31 * t9 * (-0.8e1 / 0.27e2 * t219 * t35 + 0.16e2 / 0.9e1 * t95 * t72 + 0.4e1 / 0.9e1 * t224 * t35 + 0.4e1 / 0.3e1 * t16 * t229 - 0.8e1 / 0.27e2 * t232 * t37 - 0.16e2 / 0.9e1 * t101 * t72 + 0.4e1 / 0.9e1 * t237 * t37 + 0.4e1 / 0.3e1 * t19 * t240);
374 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
375 v3rho3[2] = -t180 / 0.12e2 - t183 - t137 / 0.24e2 + t144 - t217 / 0.8e1 - 0.3e1 / 0.32e2 * t30 * t245;
376
377 t250 = t113 * t45;
378 t256 = -0.6e1 * t71 - 0.6e1 * t155;
379 t259 = t120 * t47;
380 t264 = -t256;
381 t269 = t31 * t9 * (-0.8e1 / 0.27e2 * t146 * t250 + 0.4e1 / 0.3e1 * t95 * t117 + 0.4e1 / 0.3e1 * t16 * t256 - 0.8e1 / 0.27e2 * t161 * t259 + 0.4e1 / 0.3e1 * t101 * t123 + 0.4e1 / 0.3e1 * t19 * t264);
382 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
383 v3rho3[3] = -t180 / 0.8e1 - 0.3e1 / 0.8e1 * t217 + t144 - 0.3e1 / 0.32e2 * t30 * t269;
384
385 #ifndef XC_DONT_COMPILE_LXC
386
387 if(order < 4) return;
388
389
390 t274 = t5 * t10 * t171 * t22;
391 t278 = t5 * t10 * t85 * t60;
392 t282 = t5 * t10 * t40 * t140;
393 t289 = 0.5e1 / 0.108e3 * t5 * t10 * t21 / t59 / t32;
394 t290 = t15 * t15;
395 t292 = 0.1e1 / t65 / t290;
396 t293 = t67 * t67;
397 t299 = t74 * t74;
398 t305 = 0.1e1 / t153 / t12;
399 t306 = t11 * t305;
400 t308 = -0.24e2 * t154 + 0.24e2 * t306;
401 t311 = t18 * t18;
402 t313 = 0.1e1 / t77 / t311;
403 t314 = t79 * t79;
404 t320 = t82 * t82;
405 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
406 v4rho4[0] = -t274 / 0.2e1 - t278 / 0.4e1 + t282 / 0.9e1 - t289 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * (0.40e2 / 0.81e2 * t292 * t293 - 0.16e2 / 0.9e1 * t146 * t67 * t74 + 0.4e1 / 0.3e1 * t66 * t299 + 0.16e2 / 0.9e1 * t150 * t157 + 0.4e1 / 0.3e1 * t16 * t308 + 0.40e2 / 0.81e2 * t313 * t314 - 0.16e2 / 0.9e1 * t161 * t79 * t82 + 0.4e1 / 0.3e1 * t78 * t320 + 0.16e2 / 0.9e1 * t165 * t168 - 0.4e1 / 0.3e1 * t19 * t308);
407
408 t338 = t140 * params->alpha * t1 * t52;
409 t340 = t179 * t109;
410 t341 = t340 / 0.8e1;
411 t342 = t92 * t212;
412 t367 = 0.32e2 * t98 * t305;
413 t385 = 0.32e2 * t104 * t305;
414 t394 = -0.8e1 * t187 * t154 * t35 + 0.8e1 * t200 * t154 * t37 + 0.40e2 / 0.81e2 * t292 * t45 * t147 - 0.8e1 / 0.9e1 * t184 * t35 * t74 - 0.16e2 / 0.9e1 * t146 * t11 * t71 * t67 + 0.8e1 / 0.3e1 * t66 * t71 * t35 + 0.8e1 / 0.3e1 * t187 * t71 * t74 + t367 + 0.40e2 / 0.81e2 * t313 * t47 * t162 - 0.8e1 / 0.9e1 * t197 * t37 * t82 + 0.16e2 / 0.9e1 * t161 * t11 * t71 * t79 - 0.8e1 / 0.3e1 * t78 * t71 * t37 - 0.8e1 / 0.3e1 * t200 * t71 * t82 - t385 + 0.4e1 / 0.9e1 * t95 * t157 - 0.16e2 * t16 * t154 + 0.4e1 / 0.9e1 * t101 * t168 + 0.16e2 * t19 * t154;
415 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
416 v4rho4[1] = -t274 / 0.8e1 - t278 / 0.8e1 + t282 / 0.12e2 - t289 + t338 / 0.36e2 - t341 - 0.3e1 / 0.8e1 * t342 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t394;
417
418 t404 = t179 * t128;
419 t406 = t92 * t245;
420 t417 = t11 * t11;
421 t420 = 0.1e1 / t153 / t32;
422 t459 = 0.40e2 / 0.81e2 * t292 * t113 * t67 - 0.64e2 / 0.27e2 * t184 * t35 * t11 * t71 - 0.8e1 / 0.27e2 * t219 * t74 + 0.32e2 / 0.9e1 * t66 * t417 * t420 + 0.16e2 / 0.9e1 * t95 * t71 - 0.16e2 / 0.3e1 * t95 * t155 - 0.8e1 / 0.27e2 * t146 * t117 * t67 + 0.8e1 / 0.9e1 * t66 * t229 * t35 + 0.4e1 / 0.9e1 * t224 * t74 + t367 + 0.40e2 / 0.81e2 * t313 * t120 * t79 + 0.64e2 / 0.27e2 * t197 * t37 * t11 * t71 - 0.8e1 / 0.27e2 * t232 * t82 + 0.32e2 / 0.9e1 * t78 * t417 * t420 - 0.16e2 / 0.9e1 * t101 * t71 + 0.16e2 / 0.3e1 * t101 * t155 - 0.8e1 / 0.27e2 * t161 * t123 * t79 + 0.8e1 / 0.9e1 * t78 * t240 * t37 + 0.4e1 / 0.9e1 * t237 * t82 - t385;
423 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
424 v4rho4[2] = t338 / 0.18e2 - t340 / 0.6e1 - t342 / 0.4e1 - t278 / 0.24e2 + t282 / 0.18e2 - t289 - t404 / 0.24e2 - t406 / 0.4e1 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t459;
425
426 t468 = t92 * t269;
427 t488 = 0.12e2 * t154 + 0.24e2 * t306;
428 t510 = 0.40e2 / 0.81e2 * t292 * t250 * t35 - 0.16e2 / 0.9e1 * t219 * t72 - 0.8e1 / 0.9e1 * t184 * t117 * t35 + 0.8e1 / 0.3e1 * t187 * t71 * t117 + 0.4e1 / 0.3e1 * t95 * t229 + 0.4e1 / 0.9e1 * t66 * t256 * t35 + 0.4e1 / 0.3e1 * t16 * t488 + 0.40e2 / 0.81e2 * t313 * t259 * t37 + 0.16e2 / 0.9e1 * t232 * t72 - 0.8e1 / 0.9e1 * t197 * t123 * t37 - 0.8e1 / 0.3e1 * t200 * t71 * t123 + 0.4e1 / 0.3e1 * t101 * t240 + 0.4e1 / 0.9e1 * t78 * t264 * t37 - 0.4e1 / 0.3e1 * t19 * t488;
429 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
430 v4rho4[3] = t338 / 0.12e2 - t341 - t404 / 0.8e1 - 0.3e1 / 0.8e1 * t406 + t282 / 0.36e2 - t289 - t468 / 0.8e1 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * t510;
431
432 t518 = t113 * t113;
433 t523 = t117 * t117;
434 t529 = 0.24e2 * t154 + 0.24e2 * t306;
435 t532 = t120 * t120;
436 t537 = t123 * t123;
437 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
438 v4rho4[4] = t338 / 0.9e1 - t404 / 0.4e1 - t468 / 0.2e1 - t289 - 0.3e1 / 0.32e2 * t30 * t31 * t9 * (0.40e2 / 0.81e2 * t292 * t518 - 0.16e2 / 0.9e1 * t219 * t117 + 0.4e1 / 0.3e1 * t66 * t523 + 0.16e2 / 0.9e1 * t95 * t256 + 0.4e1 / 0.3e1 * t16 * t529 + 0.40e2 / 0.81e2 * t313 * t532 - 0.16e2 / 0.9e1 * t232 * t123 + 0.4e1 / 0.3e1 * t78 * t537 + 0.16e2 / 0.9e1 * t101 * t264 - 0.4e1 / 0.3e1 * t19 * t529);
439
440 #ifndef XC_DONT_COMPILE_MXC
441
442 if(order < 5) return;
443
444
445 #endif
446
447 #endif
448
449 #endif
450
451 #endif
452
453 #endif
454
455
456 }
457
458