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