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/mgga_exc/mgga_x_gdme.mpl
11 Type of functional: mgga_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,const double * sigma,const double * lapl,const double * tau,double * zk,MGGA_OUT_PARAMS_NO_EXC (double *))19 func_unpol(const xc_func_type *p, int order, const double *rho, const double *sigma, const double *lapl, const double *tau, double *zk, MGGA_OUT_PARAMS_NO_EXC(double *))
20 {
21
22 #ifndef XC_DONT_COMPILE_EXC
23 double t2, t4, t6, t7, t8, t9, t10, t11;
24 double t16, t18, t19, t20, t29, t30, t31, t32;
25 double t33, t35, t36, t38, t46, t48;
26
27 #ifndef XC_DONT_COMPILE_VXC
28 double t51, t52, t54, t55, t56, t58, t59, t65;
29 double t69, t70;
30
31 #ifndef XC_DONT_COMPILE_FXC
32 double t74, t79, t80, t83, t85, t92, t96;
33
34 #ifndef XC_DONT_COMPILE_KXC
35 double t103, t108, t110, t117, t122;
36
37 #ifndef XC_DONT_COMPILE_LXC
38 double t138, t150;
39 #endif
40
41 #endif
42
43 #endif
44
45 #endif
46
47 #endif
48
49
50 mgga_x_gdme_params *params;
51
52 assert(p->params != NULL);
53 params = (mgga_x_gdme_params * )(p->params);
54
55 t2 = M_CBRT3;
56 t4 = POW_1_3(0.1e1 / M_PI);
57 t6 = M_CBRT4;
58 t7 = t6 * t6;
59 t8 = t2 * t4 * t7;
60 t9 = M_CBRT2;
61 t10 = t9 * t9;
62 t11 = POW_1_3(rho[0]);
63 t16 = 0.1e1 / t4;
64 t18 = M_PI * M_PI;
65 t19 = POW_1_3(t18);
66 t20 = t19 * t19;
67 t29 = 0.1e1 / t19 / t18;
68 t30 = params->a * params->a;
69 t31 = t30 - params->a + 0.1e1 / 0.2e1;
70 t32 = t31 * lapl[0];
71 t33 = t11 * t11;
72 t35 = 0.1e1 / t33 / rho[0];
73 t36 = t10 * t35;
74 t38 = tau[0] * t10;
75 t46 = 0.2e1 / 0.9e1 * (params->AA + 0.3e1 / 0.5e1 * params->BB) * t9 * t16 * t6 / t20 + params->BB * t2 * t16 * t6 * t10 * t29 * (t32 * t36 - 0.2e1 * t38 * t35) / 0.27e2;
76 t48 = t8 * t10 * t11 * t46;
77 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
78 zk[0] = -0.3e1 / 0.16e2 * t48;
79
80 #ifndef XC_DONT_COMPILE_VXC
81
82 if(order < 1) return;
83
84
85 t51 = t11 * rho[0];
86 t52 = t2 * t2;
87 t54 = t51 * t52 * t9;
88 t55 = params->BB * t29;
89 t56 = rho[0] * rho[0];
90 t58 = 0.1e1 / t33 / t56;
91 t59 = t10 * t58;
92 t65 = t55 * (-0.5e1 / 0.3e1 * t32 * t59 + 0.10e2 / 0.3e1 * t38 * t58);
93 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
94 vrho[0] = -t48 / 0.4e1 - t54 * t65 / 0.18e2;
95
96 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
97 vsigma[0] = 0.0e0;
98
99 t69 = 0.1e1 / t11 * t52;
100 t70 = t55 * t31;
101 if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
102 vlapl[0] = -t69 * t70 / 0.9e1;
103
104 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
105 vtau[0] = 0.2e1 / 0.9e1 * t69 * t55;
106
107 #ifndef XC_DONT_COMPILE_FXC
108
109 if(order < 2) return;
110
111
112 t74 = 0.1e1 / t33;
113 t79 = t52 * t9;
114 t80 = t79 * t11;
115 t83 = t56 * rho[0];
116 t85 = 0.1e1 / t33 / t83;
117 t92 = t55 * (0.40e2 / 0.9e1 * t32 * t10 * t85 - 0.80e2 / 0.9e1 * t38 * t85);
118 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
119 v2rho2[0] = -t8 * t10 * t74 * t46 / 0.12e2 - 0.4e1 / 0.27e2 * t80 * t65 - t54 * t92 / 0.18e2;
120
121 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
122 v2rhosigma[0] = 0.0e0;
123
124 t96 = 0.1e1 / t51 * t52;
125 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
126 v2rholapl[0] = t96 * t70 / 0.27e2;
127
128 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
129 v2rhotau[0] = -0.2e1 / 0.27e2 * t96 * t55;
130
131 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
132 v2sigma2[0] = 0.0e0;
133
134 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
135 v2sigmalapl[0] = 0.0e0;
136
137 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
138 v2sigmatau[0] = 0.0e0;
139
140 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
141 v2lapl2[0] = 0.0e0;
142
143 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
144 v2lapltau[0] = 0.0e0;
145
146 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
147 v2tau2[0] = 0.0e0;
148
149 #ifndef XC_DONT_COMPILE_KXC
150
151 if(order < 3) return;
152
153
154 t103 = t79 * t74;
155 t108 = t56 * t56;
156 t110 = 0.1e1 / t33 / t108;
157 t117 = t55 * (-0.440e3 / 0.27e2 * t32 * t10 * t110 + 0.880e3 / 0.27e2 * t38 * t110);
158 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
159 v3rho3[0] = t8 * t36 * t46 / 0.18e2 - 0.2e1 / 0.27e2 * t103 * t65 - 0.2e1 / 0.9e1 * t80 * t92 - t54 * t117 / 0.18e2;
160
161 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
162 v3rho2sigma[0] = 0.0e0;
163
164 t122 = 0.1e1 / t11 / t56 * t52;
165 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
166 v3rho2lapl[0] = -0.4e1 / 0.81e2 * t122 * t70;
167
168 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
169 v3rho2tau[0] = 0.8e1 / 0.81e2 * t122 * t55;
170
171 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
172 v3rhosigma2[0] = 0.0e0;
173
174 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
175 v3rhosigmalapl[0] = 0.0e0;
176
177 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
178 v3rhosigmatau[0] = 0.0e0;
179
180 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
181 v3rholapl2[0] = 0.0e0;
182
183 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
184 v3rholapltau[0] = 0.0e0;
185
186 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
187 v3rhotau2[0] = 0.0e0;
188
189 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
190 v3sigma3[0] = 0.0e0;
191
192 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
193 v3sigma2lapl[0] = 0.0e0;
194
195 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
196 v3sigma2tau[0] = 0.0e0;
197
198 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
199 v3sigmalapl2[0] = 0.0e0;
200
201 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
202 v3sigmalapltau[0] = 0.0e0;
203
204 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
205 v3sigmatau2[0] = 0.0e0;
206
207 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
208 v3lapl3[0] = 0.0e0;
209
210 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
211 v3lapl2tau[0] = 0.0e0;
212
213 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
214 v3lapltau2[0] = 0.0e0;
215
216 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
217 v3tau3[0] = 0.0e0;
218
219 #ifndef XC_DONT_COMPILE_LXC
220
221 if(order < 4) return;
222
223
224 t138 = 0.1e1 / t33 / t108 / rho[0];
225 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
226 v4rho4[0] = -0.5e1 / 0.54e2 * t8 * t59 * t46 + 0.16e2 / 0.243e3 * t79 * t35 * t65 - 0.4e1 / 0.27e2 * t103 * t92 - 0.8e1 / 0.27e2 * t80 * t117 - t54 * t55 * (0.6160e4 / 0.81e2 * t32 * t10 * t138 - 0.12320e5 / 0.81e2 * t38 * t138) / 0.18e2;
227
228 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
229 v4rho3sigma[0] = 0.0e0;
230
231 t150 = 0.1e1 / t11 / t83 * t52;
232 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
233 v4rho3lapl[0] = 0.28e2 / 0.243e3 * t150 * t70;
234
235 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
236 v4rho3tau[0] = -0.56e2 / 0.243e3 * t150 * t55;
237
238 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
239 v4rho2sigma2[0] = 0.0e0;
240
241 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
242 v4rho2sigmalapl[0] = 0.0e0;
243
244 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
245 v4rho2sigmatau[0] = 0.0e0;
246
247 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
248 v4rho2lapl2[0] = 0.0e0;
249
250 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
251 v4rho2lapltau[0] = 0.0e0;
252
253 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
254 v4rho2tau2[0] = 0.0e0;
255
256 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
257 v4rhosigma3[0] = 0.0e0;
258
259 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
260 v4rhosigma2lapl[0] = 0.0e0;
261
262 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
263 v4rhosigma2tau[0] = 0.0e0;
264
265 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
266 v4rhosigmalapl2[0] = 0.0e0;
267
268 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
269 v4rhosigmalapltau[0] = 0.0e0;
270
271 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
272 v4rhosigmatau2[0] = 0.0e0;
273
274 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
275 v4rholapl3[0] = 0.0e0;
276
277 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
278 v4rholapl2tau[0] = 0.0e0;
279
280 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
281 v4rholapltau2[0] = 0.0e0;
282
283 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
284 v4rhotau3[0] = 0.0e0;
285
286 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
287 v4sigma4[0] = 0.0e0;
288
289 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
290 v4sigma3lapl[0] = 0.0e0;
291
292 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
293 v4sigma3tau[0] = 0.0e0;
294
295 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
296 v4sigma2lapl2[0] = 0.0e0;
297
298 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
299 v4sigma2lapltau[0] = 0.0e0;
300
301 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
302 v4sigma2tau2[0] = 0.0e0;
303
304 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
305 v4sigmalapl3[0] = 0.0e0;
306
307 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
308 v4sigmalapl2tau[0] = 0.0e0;
309
310 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
311 v4sigmalapltau2[0] = 0.0e0;
312
313 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
314 v4sigmatau3[0] = 0.0e0;
315
316 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
317 v4lapl4[0] = 0.0e0;
318
319 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
320 v4lapl3tau[0] = 0.0e0;
321
322 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
323 v4lapl2tau2[0] = 0.0e0;
324
325 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
326 v4lapltau3[0] = 0.0e0;
327
328 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
329 v4tau4[0] = 0.0e0;
330
331 #ifndef XC_DONT_COMPILE_MXC
332
333 if(order < 5) return;
334
335
336 #endif
337
338 #endif
339
340 #endif
341
342 #endif
343
344 #endif
345
346
347 }
348
349
350 static inline void
func_ferr(const xc_func_type * p,int order,const double * rho,const double * sigma,const double * lapl,const double * tau,double * zk,MGGA_OUT_PARAMS_NO_EXC (double *))351 func_ferr(const xc_func_type *p, int order, const double *rho, const double *sigma, const double *lapl, const double *tau, double *zk, MGGA_OUT_PARAMS_NO_EXC(double *))
352 {
353
354 #ifndef XC_DONT_COMPILE_EXC
355 double t2, t4, t5, t6, t7, t8, t12, t14;
356 double t16, t17, t18, t25, t28, t29, t30, t31;
357 double t32, t34, t43, t45;
358
359 #ifndef XC_DONT_COMPILE_VXC
360 double t48, t49, t51, t52, t53, t55, t61, t65;
361 double t67, t71;
362
363 #ifndef XC_DONT_COMPILE_FXC
364 double t73, t79, t82, t84, t90, t94;
365
366 #ifndef XC_DONT_COMPILE_KXC
367 double t104, t109, t111, t117, t122;
368
369 #ifndef XC_DONT_COMPILE_LXC
370 double t141, t152;
371 #endif
372
373 #endif
374
375 #endif
376
377 #endif
378
379 #endif
380
381
382 mgga_x_gdme_params *params;
383
384 assert(p->params != NULL);
385 params = (mgga_x_gdme_params * )(p->params);
386
387 t2 = M_CBRT3;
388 t4 = POW_1_3(0.1e1 / M_PI);
389 t5 = t2 * t4;
390 t6 = M_CBRT4;
391 t7 = t6 * t6;
392 t8 = POW_1_3(rho[0]);
393 t12 = M_CBRT2;
394 t14 = 0.1e1 / t4;
395 t16 = M_PI * M_PI;
396 t17 = POW_1_3(t16);
397 t18 = t17 * t17;
398 t25 = t12 * t12;
399 t28 = 0.1e1 / t17 / t16;
400 t29 = params->a * params->a;
401 t30 = t29 - params->a + 0.1e1 / 0.2e1;
402 t31 = t30 * lapl[0];
403 t32 = t8 * t8;
404 t34 = 0.1e1 / t32 / rho[0];
405 t43 = 0.2e1 / 0.9e1 * (params->AA + 0.3e1 / 0.5e1 * params->BB) * t12 * t14 * t6 / t18 + params->BB * t2 * t14 * t6 * t25 * t28 * (t31 * t34 - 0.2e1 * tau[0] * t34) / 0.27e2;
406 t45 = t5 * t7 * t8 * t43;
407 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
408 zk[0] = -0.3e1 / 0.8e1 * t45;
409
410 #ifndef XC_DONT_COMPILE_VXC
411
412 if(order < 1) return;
413
414
415 t48 = t8 * rho[0];
416 t49 = t2 * t2;
417 t51 = t48 * t49 * params->BB;
418 t52 = t25 * t28;
419 t53 = rho[0] * rho[0];
420 t55 = 0.1e1 / t32 / t53;
421 t61 = t52 * (-0.5e1 / 0.3e1 * t31 * t55 + 0.10e2 / 0.3e1 * tau[0] * t55);
422 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
423 vrho[0] = -t45 / 0.2e1 - t51 * t61 / 0.18e2;
424
425 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
426 vsigma[0] = 0.0e0;
427
428 t65 = 0.1e1 / t8 * t49;
429 t67 = t52 * t30;
430 if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
431 vlapl[0] = -t65 * params->BB * t67 / 0.18e2;
432
433 t71 = params->BB * t25 * t28;
434 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
435 vtau[0] = t65 * t71 / 0.9e1;
436
437 #ifndef XC_DONT_COMPILE_FXC
438
439 if(order < 2) return;
440
441
442 t73 = 0.1e1 / t32;
443 t79 = t49 * t8 * params->BB;
444 t82 = t53 * rho[0];
445 t84 = 0.1e1 / t32 / t82;
446 t90 = t52 * (0.40e2 / 0.9e1 * t31 * t84 - 0.80e2 / 0.9e1 * tau[0] * t84);
447 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
448 v2rho2[0] = -t5 * t7 * t73 * t43 / 0.6e1 - 0.4e1 / 0.27e2 * t79 * t61 - t51 * t90 / 0.18e2;
449
450 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
451 v2rhosigma[0] = 0.0e0;
452
453 t94 = 0.1e1 / t48 * t49;
454 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
455 v2rholapl[0] = t94 * params->BB * t67 / 0.54e2;
456
457 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
458 v2rhotau[0] = -t94 * t71 / 0.27e2;
459
460 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
461 v2sigma2[0] = 0.0e0;
462
463 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
464 v2sigmalapl[0] = 0.0e0;
465
466 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
467 v2sigmatau[0] = 0.0e0;
468
469 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
470 v2lapl2[0] = 0.0e0;
471
472 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
473 v2lapltau[0] = 0.0e0;
474
475 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
476 v2tau2[0] = 0.0e0;
477
478 #ifndef XC_DONT_COMPILE_KXC
479
480 if(order < 3) return;
481
482
483 t104 = t49 * t73 * params->BB;
484 t109 = t53 * t53;
485 t111 = 0.1e1 / t32 / t109;
486 t117 = t52 * (-0.440e3 / 0.27e2 * t31 * t111 + 0.880e3 / 0.27e2 * tau[0] * t111);
487 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
488 v3rho3[0] = t5 * t7 * t34 * t43 / 0.9e1 - 0.2e1 / 0.27e2 * t104 * t61 - 0.2e1 / 0.9e1 * t79 * t90 - t51 * t117 / 0.18e2;
489
490 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
491 v3rho2sigma[0] = 0.0e0;
492
493 t122 = 0.1e1 / t8 / t53 * t49;
494 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
495 v3rho2lapl[0] = -0.2e1 / 0.81e2 * t122 * params->BB * t67;
496
497 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
498 v3rho2tau[0] = 0.4e1 / 0.81e2 * t122 * t71;
499
500 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
501 v3rhosigma2[0] = 0.0e0;
502
503 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
504 v3rhosigmalapl[0] = 0.0e0;
505
506 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
507 v3rhosigmatau[0] = 0.0e0;
508
509 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
510 v3rholapl2[0] = 0.0e0;
511
512 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
513 v3rholapltau[0] = 0.0e0;
514
515 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
516 v3rhotau2[0] = 0.0e0;
517
518 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
519 v3sigma3[0] = 0.0e0;
520
521 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
522 v3sigma2lapl[0] = 0.0e0;
523
524 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
525 v3sigma2tau[0] = 0.0e0;
526
527 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
528 v3sigmalapl2[0] = 0.0e0;
529
530 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
531 v3sigmalapltau[0] = 0.0e0;
532
533 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
534 v3sigmatau2[0] = 0.0e0;
535
536 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
537 v3lapl3[0] = 0.0e0;
538
539 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
540 v3lapl2tau[0] = 0.0e0;
541
542 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
543 v3lapltau2[0] = 0.0e0;
544
545 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
546 v3tau3[0] = 0.0e0;
547
548 #ifndef XC_DONT_COMPILE_LXC
549
550 if(order < 4) return;
551
552
553 t141 = 0.1e1 / t32 / t109 / rho[0];
554 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
555 v4rho4[0] = -0.5e1 / 0.27e2 * t5 * t7 * t55 * t43 + 0.16e2 / 0.243e3 * t49 * t34 * params->BB * t61 - 0.4e1 / 0.27e2 * t104 * t90 - 0.8e1 / 0.27e2 * t79 * t117 - t51 * t52 * (0.6160e4 / 0.81e2 * t31 * t141 - 0.12320e5 / 0.81e2 * tau[0] * t141) / 0.18e2;
556
557 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
558 v4rho3sigma[0] = 0.0e0;
559
560 t152 = 0.1e1 / t8 / t82 * t49;
561 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
562 v4rho3lapl[0] = 0.14e2 / 0.243e3 * t152 * params->BB * t67;
563
564 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
565 v4rho3tau[0] = -0.28e2 / 0.243e3 * t152 * t71;
566
567 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
568 v4rho2sigma2[0] = 0.0e0;
569
570 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
571 v4rho2sigmalapl[0] = 0.0e0;
572
573 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
574 v4rho2sigmatau[0] = 0.0e0;
575
576 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
577 v4rho2lapl2[0] = 0.0e0;
578
579 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
580 v4rho2lapltau[0] = 0.0e0;
581
582 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
583 v4rho2tau2[0] = 0.0e0;
584
585 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
586 v4rhosigma3[0] = 0.0e0;
587
588 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
589 v4rhosigma2lapl[0] = 0.0e0;
590
591 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
592 v4rhosigma2tau[0] = 0.0e0;
593
594 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
595 v4rhosigmalapl2[0] = 0.0e0;
596
597 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
598 v4rhosigmalapltau[0] = 0.0e0;
599
600 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
601 v4rhosigmatau2[0] = 0.0e0;
602
603 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
604 v4rholapl3[0] = 0.0e0;
605
606 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
607 v4rholapl2tau[0] = 0.0e0;
608
609 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
610 v4rholapltau2[0] = 0.0e0;
611
612 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
613 v4rhotau3[0] = 0.0e0;
614
615 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
616 v4sigma4[0] = 0.0e0;
617
618 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
619 v4sigma3lapl[0] = 0.0e0;
620
621 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
622 v4sigma3tau[0] = 0.0e0;
623
624 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
625 v4sigma2lapl2[0] = 0.0e0;
626
627 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
628 v4sigma2lapltau[0] = 0.0e0;
629
630 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
631 v4sigma2tau2[0] = 0.0e0;
632
633 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
634 v4sigmalapl3[0] = 0.0e0;
635
636 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
637 v4sigmalapl2tau[0] = 0.0e0;
638
639 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
640 v4sigmalapltau2[0] = 0.0e0;
641
642 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
643 v4sigmatau3[0] = 0.0e0;
644
645 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
646 v4lapl4[0] = 0.0e0;
647
648 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
649 v4lapl3tau[0] = 0.0e0;
650
651 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
652 v4lapl2tau2[0] = 0.0e0;
653
654 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
655 v4lapltau3[0] = 0.0e0;
656
657 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
658 v4tau4[0] = 0.0e0;
659
660 #ifndef XC_DONT_COMPILE_MXC
661
662 if(order < 5) return;
663
664
665 #endif
666
667 #endif
668
669 #endif
670
671 #endif
672
673 #endif
674
675
676 }
677
678
679 static inline void
func_pol(const xc_func_type * p,int order,const double * rho,const double * sigma,const double * lapl,const double * tau,double * zk,MGGA_OUT_PARAMS_NO_EXC (double *))680 func_pol(const xc_func_type *p, int order, const double *rho, const double *sigma, const double *lapl, const double *tau, double *zk, MGGA_OUT_PARAMS_NO_EXC(double *))
681 {
682
683 #ifndef XC_DONT_COMPILE_EXC
684 double t2, t4, t5, t6, t7, t8, t9, t10;
685 double t11, t12, t14, t15, t16, t17, t21, t23;
686 double t25, t26, t27, t31, t33, t34, t35, t37;
687 double t38, t39, t40, t41, t42, t44, t53, t55;
688 double t57, t58, t59, t61, t62, t63, t65, t74;
689 double t76;
690
691 #ifndef XC_DONT_COMPILE_VXC
692 double t78, t79, t80, t81, t82, t83, t85, t86;
693 double t88, t89, t90, t91, t94, t95, t96, t97;
694 double t98, t99, t100, t102, t107, t109, t110, t111;
695 double t112, t113, t114, t116, t117, t120, t121, t125;
696 double t126, t128, t129, t130, t131, t133, t134, t135;
697 double t136, t137, t139, t144, t146, t147, t148, t152;
698 double t154, t155, t156, t157, t161, t162, t163, t166;
699 double t168, t170, t172;
700
701 #ifndef XC_DONT_COMPILE_FXC
702 double t174, t176, t177, t178, t179, t180, t181, t183;
703 double t184, t185, t187, t188, t189, t190, t191, t192;
704 double t193, t194, t195, t196, t197, t198, t199, t200;
705 double t202, t203, t205, t208, t209, t210, t211, t212;
706 double t215, t220, t222, t223, t224, t225, t226, t227;
707 double t228, t229, t231, t232, t233, t235, t236, t237;
708 double t238, t240, t241, t244, t245, t246, t249, t250;
709 double t251, t253, t256, t258, t259, t260, t263, t264;
710 double t265, t267, t272, t273, t274, t276, t279, t281;
711 double t282, t284, t287, t288, t290, t291, t293, t294;
712 double t296, t299, t300, t302, t303, t304, t305, t306;
713 double t308, t309, t310, t311, t313, t314, t315, t316;
714 double t317, t318, t319, t320, t321, t323, t324, t325;
715 double t328, t333, t335, t336, t337, t338, t340, t341;
716 double t343, t344, t346, t349, t350, t353, t354, t355;
717 double t356, t358, t362, t367, t370, t371, t374, t376;
718 double t377, t381, t384, t385, t386, t389, t390, t394;
719 double t395, t399;
720
721 #ifndef XC_DONT_COMPILE_KXC
722 double t403, t407, t412, t413, t414, t415, t417, t418;
723 double t419, t421, t423, t425, t426, t427, t429, t430;
724 double t432, t433, t434, t435, t437, t438, t439, t441;
725 double t443, t445, t446, t447, t448, t450, t451, t452;
726 double t453, t455, t456, t457, t458, t459, t463, t464;
727 double t465, t467, t468, t469, t470, t471, t472, t473;
728 double t475, t477, t482, t484, t485, t486, t488, t490;
729 double t491, t492, t494, t495, t497, t498, t501, t502;
730 double t504, t506, t507, t509, t510, t513, t514, t515;
731 double t519, t520, t523, t524, t526, t528, t530, t531;
732 double t533, t535, t537, t539, t540, t541, t542, t543;
733 double t546, t547, t549, t551, t553, t554, t555, t556;
734 double t558, t559, t561, t562, t563, t565, t566, t568;
735 double t570, t572, t573, t574, t575, t577, t579, t580;
736 double t581, t583, t584, t585, t586, t587, t591, t593;
737 double t595, t598, t599, t601, t603, t604, t607, t609;
738 double t611, t614, t625, t626, t628, t631, t636, t638;
739 double t641, t642, t644, t646, t648, t650, t652, t654;
740 double t655, t656, t658, t660, t662, t663, t664, t666;
741 double t668, t670, t672, t674, t675, t676, t677, t679;
742 double t680, t683, t684, t686, t688, t689, t691, t694;
743 double t697, t700, t703, t714, t715, t717, t720, t728;
744 double t734, t735, t737, t738, t739, t741, t742, t743;
745 double t744, t746, t747, t748, t750, t751, t753, t754;
746 double t755, t756, t757, t758, t759, t760, t761, t762;
747 double t763, t764, t765, t766, t768, t770, t775, t777;
748 double t778, t779, t781, t782, t784, t785, t786, t787;
749 double t789, t790, t791, t794, t796, t798, t799, t801;
750 double t802, t804, t808, t812, t815, t816, t819, t821;
751 double t823, t824, t825, t827, t831, t836, t838, t839;
752 double t840, t841, t842, t846, t849, t850, t852, t856;
753 double t858, t859, t860, t861, t864, t866, t869, t872;
754 double t877, t881, t887, t891, t895, t898, t899, t902;
755 double t905, t906, t909, t913, t917, t921, t924, t927;
756 double t929, t933, t938, t940, t943, t944, t947, t948;
757 double t952, t954, t957, t958, t961, t964, t968, t969;
758 double t972, t973, t978, t979, t982, t983, t986, t987;
759 double t991;
760
761 #ifndef XC_DONT_COMPILE_LXC
762 double t1007, t1009, t1011, t1013, t1036, t1042, t1044, t1045;
763 double t1052, t1054, t1056, t1058, t1063, t1064, t1066, t1067;
764 double t1068, t1073, t1075, t1077, t1079, t1081, t1086, t1088;
765 double t1090, t1096, t1099, t1101, t1103, t1107, t1118, t1126;
766 double t1129, t1132, t1134, t1136, t1137, t1138, t1143, t1145;
767 double t1150, t1154, t1155, t1168, t1176, t1179, t1180, t1182;
768 double t1188, t1190, t1192, t1195, t1200, t1201, t1204, t1210;
769 double t1217, t1224, t1230, t1232, t1234, t1242, t1243, t1244;
770 double t1247, t1248, t1253, t1254, t1255, t1257, t1259, t1277;
771 double t1278, t1285, t1287, t1292, t1319, t1322, t1325, t1331;
772 double t1335, t1337, t1342, t1345, t1346, t1363, t1366, t1368;
773 double t1376, t1378, t1380, t1381, t1383, t1387, t1393, t1398;
774 double t1404, t1406, t1413, t1416, t1423, t1426, t1433, t1440;
775 double t1454, t1468, t1469, t1488, t1504, t1515, t1536, t1544;
776 double t1546, t1552, t1560, t1563, t1565, t1571, t1576, t1581;
777 double t1584, t1601, t1617, t1631, t1640, t1655, t1712, t1741;
778 double t1754, t1765, t1768, t1771, t1773, t1776, t1785, t1795;
779 double t1797, t1799, t1801, t1808, t1814, t1818, t1828, t1834;
780 double t1838, t1845, t1851, t1863, t1873, t1878, t1883, t1898;
781 double t1943, t1956, t1966, t1967, t1981, t1983, t1986, t1987;
782 double t1989, t1990, t1999, t2007, t2009, t2011, t2013, t2015;
783 double t2017, t2021, t2022, t2023, t2028, t2042, t2043, t2047;
784 double t2051, t2052, t2057, t2058, t2060, t2072, t2074, t2077;
785 double t2084, t2088, t2089, t2092, t2093, t2097, t2098, t2102;
786 double t2108, t2111, t2124, t2126, t2134, t2138, t2140, t2156;
787 double t2158, t2162, t2166, t2168, t2176, t2178, t2186, t2199;
788 double t2229, t2253, t2256, t2257, t2259, t2261, t2263, t2265;
789 double t2282, t2284, t2286, t2289, t2307, t2310, t2312, t2313;
790 double t2335, t2342, t2345, t2346, t2351, t2361, t2364, t2377;
791 double t2389, t2406, t2414, t2423, t2432;
792 #endif
793
794 #endif
795
796 #endif
797
798 #endif
799
800 #endif
801
802
803 mgga_x_gdme_params *params;
804
805 assert(p->params != NULL);
806 params = (mgga_x_gdme_params * )(p->params);
807
808 t2 = M_CBRT3;
809 t4 = POW_1_3(0.1e1 / M_PI);
810 t5 = t2 * t4;
811 t6 = M_CBRT4;
812 t7 = t6 * t6;
813 t8 = t5 * t7;
814 t9 = rho[0] - rho[1];
815 t10 = rho[0] + rho[1];
816 t11 = 0.1e1 / t10;
817 t12 = t9 * t11;
818 t14 = 0.1e1 / 0.2e1 + t12 / 0.2e1;
819 t15 = POW_1_3(t14);
820 t16 = t15 * t14;
821 t17 = POW_1_3(t10);
822 t21 = M_CBRT2;
823 t23 = 0.1e1 / t4;
824 t25 = M_PI * M_PI;
825 t26 = POW_1_3(t25);
826 t27 = t26 * t26;
827 t31 = 0.2e1 / 0.9e1 * (params->AA + 0.3e1 / 0.5e1 * params->BB) * t21 * t23 * t6 / t27;
828 t33 = params->BB * t2 * t23;
829 t34 = t21 * t21;
830 t35 = t6 * t34;
831 t37 = 0.1e1 / t26 / t25;
832 t38 = params->a * params->a;
833 t39 = t38 - params->a + 0.1e1 / 0.2e1;
834 t40 = t39 * lapl[0];
835 t41 = POW_1_3(rho[0]);
836 t42 = t41 * t41;
837 t44 = 0.1e1 / t42 / rho[0];
838 t53 = t31 + t33 * t35 * t37 * (t40 * t44 - 0.2e1 * tau[0] * t44) / 0.27e2;
839 t55 = t8 * t16 * t17 * t53;
840 t57 = 0.1e1 / 0.2e1 - t12 / 0.2e1;
841 t58 = POW_1_3(t57);
842 t59 = t58 * t57;
843 t61 = t39 * lapl[1];
844 t62 = POW_1_3(rho[1]);
845 t63 = t62 * t62;
846 t65 = 0.1e1 / t63 / rho[1];
847 t74 = t31 + t33 * t35 * t37 * (t61 * t65 - 0.2e1 * tau[1] * t65) / 0.27e2;
848 t76 = t8 * t59 * t17 * t74;
849 if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
850 zk[0] = -0.3e1 / 0.8e1 * t55 - 0.3e1 / 0.8e1 * t76;
851
852 #ifndef XC_DONT_COMPILE_VXC
853
854 if(order < 1) return;
855
856
857 t78 = 0.3e1 / 0.8e1 * t55;
858 t79 = 0.3e1 / 0.8e1 * t76;
859 t80 = t15 * t17;
860 t81 = t10 * t10;
861 t82 = 0.1e1 / t81;
862 t83 = t9 * t82;
863 t85 = t11 / 0.2e1 - t83 / 0.2e1;
864 t86 = t53 * t85;
865 t88 = t8 * t80 * t86;
866 t89 = t88 / 0.2e1;
867 t90 = t17 * t17;
868 t91 = 0.1e1 / t90;
869 t94 = t8 * t16 * t91 * t53;
870 t95 = t94 / 0.8e1;
871 t96 = t2 * t2;
872 t97 = t96 * t16;
873 t98 = t97 * t17;
874 t99 = params->BB * t34;
875 t100 = rho[0] * rho[0];
876 t102 = 0.1e1 / t42 / t100;
877 t107 = -0.5e1 / 0.3e1 * t40 * t102 + 0.10e2 / 0.3e1 * tau[0] * t102;
878 t109 = t99 * t37 * t107;
879 t110 = t98 * t109;
880 t111 = t110 / 0.18e2;
881 t112 = t58 * t17;
882 t113 = -t85;
883 t114 = t74 * t113;
884 t116 = t8 * t112 * t114;
885 t117 = t116 / 0.2e1;
886 t120 = t8 * t59 * t91 * t74;
887 t121 = t120 / 0.8e1;
888 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
889 vrho[0] = -t78 - t79 + t10 * (-t89 - t95 - t111 - t117 - t121);
890
891 t125 = -t11 / 0.2e1 - t83 / 0.2e1;
892 t126 = t53 * t125;
893 t128 = t8 * t80 * t126;
894 t129 = t128 / 0.2e1;
895 t130 = -t125;
896 t131 = t74 * t130;
897 t133 = t8 * t112 * t131;
898 t134 = t133 / 0.2e1;
899 t135 = t96 * t59;
900 t136 = t135 * t17;
901 t137 = rho[1] * rho[1];
902 t139 = 0.1e1 / t63 / t137;
903 t144 = -0.5e1 / 0.3e1 * t61 * t139 + 0.10e2 / 0.3e1 * tau[1] * t139;
904 t146 = t99 * t37 * t144;
905 t147 = t136 * t146;
906 t148 = t147 / 0.18e2;
907 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
908 vrho[1] = -t78 - t79 + t10 * (-t129 - t95 - t134 - t121 - t148);
909
910 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
911 vsigma[0] = 0.0e0;
912
913 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
914 vsigma[1] = 0.0e0;
915
916 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
917 vsigma[2] = 0.0e0;
918
919 t152 = t17 * t10 * t96;
920 t154 = t152 * t16 * params->BB;
921 t155 = t34 * t37;
922 t156 = t39 * t44;
923 t157 = t155 * t156;
924 if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
925 vlapl[0] = -t154 * t157 / 0.18e2;
926
927 t161 = t152 * t59 * params->BB;
928 t162 = t39 * t65;
929 t163 = t155 * t162;
930 if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
931 vlapl[1] = -t161 * t163 / 0.18e2;
932
933 t166 = t152 * t16;
934 t168 = t99 * t37 * t44;
935 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
936 vtau[0] = t166 * t168 / 0.9e1;
937
938 t170 = t152 * t59;
939 t172 = t99 * t37 * t65;
940 if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
941 vtau[1] = t170 * t172 / 0.9e1;
942
943 #ifndef XC_DONT_COMPILE_FXC
944
945 if(order < 2) return;
946
947
948 t174 = t94 / 0.4e1;
949 t176 = t120 / 0.4e1;
950 t177 = t15 * t15;
951 t178 = 0.1e1 / t177;
952 t179 = t178 * t17;
953 t180 = t85 * t85;
954 t181 = t53 * t180;
955 t183 = t8 * t179 * t181;
956 t184 = t183 / 0.6e1;
957 t185 = t15 * t91;
958 t187 = t8 * t185 * t86;
959 t188 = t187 / 0.3e1;
960 t189 = t96 * t15;
961 t190 = t17 * params->BB;
962 t191 = t189 * t190;
963 t192 = t107 * t85;
964 t193 = t155 * t192;
965 t194 = t191 * t193;
966 t195 = 0.4e1 / 0.27e2 * t194;
967 t196 = t81 * t10;
968 t197 = 0.1e1 / t196;
969 t198 = t9 * t197;
970 t199 = -t82 + t198;
971 t200 = t53 * t199;
972 t202 = t8 * t80 * t200;
973 t203 = t202 / 0.2e1;
974 t205 = 0.1e1 / t90 / t10;
975 t208 = t8 * t16 * t205 * t53;
976 t209 = t208 / 0.12e2;
977 t210 = t97 * t91;
978 t211 = t210 * t109;
979 t212 = t211 / 0.27e2;
980 t215 = 0.1e1 / t42 / t100 / rho[0];
981 t220 = 0.40e2 / 0.9e1 * t40 * t215 - 0.80e2 / 0.9e1 * tau[0] * t215;
982 t222 = t99 * t37 * t220;
983 t223 = t98 * t222;
984 t224 = t223 / 0.18e2;
985 t225 = t58 * t58;
986 t226 = 0.1e1 / t225;
987 t227 = t226 * t17;
988 t228 = t113 * t113;
989 t229 = t74 * t228;
990 t231 = t8 * t227 * t229;
991 t232 = t231 / 0.6e1;
992 t233 = t58 * t91;
993 t235 = t8 * t233 * t114;
994 t236 = t235 / 0.3e1;
995 t237 = -t199;
996 t238 = t74 * t237;
997 t240 = t8 * t112 * t238;
998 t241 = t240 / 0.2e1;
999 t244 = t8 * t59 * t205 * t74;
1000 t245 = t244 / 0.12e2;
1001 t246 = -t184 - t188 - t195 - t203 + t209 - t212 - t224 - t232 - t236 - t241 + t245;
1002 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1003 v2rho2[0] = -t88 - t174 - t110 / 0.9e1 - t116 - t176 + t10 * t246;
1004
1005 t249 = t5 * t7 * t178;
1006 t250 = t17 * t53;
1007 t251 = t125 * t85;
1008 t253 = t249 * t250 * t251;
1009 t256 = t8 * t185 * t126;
1010 t258 = t107 * t125;
1011 t259 = t155 * t258;
1012 t260 = t191 * t259;
1013 t263 = 0.1e1 / t90 / t81;
1014 t264 = t15 * t263;
1015 t265 = t53 * t9;
1016 t267 = t8 * t264 * t265;
1017 t272 = t5 * t7 * t226;
1018 t273 = t17 * t74;
1019 t274 = t130 * t113;
1020 t276 = t272 * t273 * t274;
1021 t279 = t8 * t233 * t131;
1022 t281 = t58 * t263;
1023 t282 = t74 * t9;
1024 t284 = t8 * t281 * t282;
1025 t287 = t96 * t58;
1026 t288 = t287 * t190;
1027 t290 = t155 * t144 * t113;
1028 t291 = t288 * t290;
1029 t293 = t135 * t91;
1030 t294 = t293 * t146;
1031 t296 = -t253 / 0.6e1 - t256 / 0.6e1 - 0.2e1 / 0.27e2 * t260 - t267 / 0.2e1 - t187 / 0.6e1 + t209 - t211 / 0.54e2 - t276 / 0.6e1 - t279 / 0.6e1 + t284 / 0.2e1 - t235 / 0.6e1 + t245 - 0.2e1 / 0.27e2 * t291 - t294 / 0.54e2;
1032 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1033 v2rho2[1] = t10 * t296 - t111 - t117 - t129 - t134 - t148 - t174 - t176 - t89;
1034
1035 t299 = t125 * t125;
1036 t300 = t53 * t299;
1037 t302 = t8 * t179 * t300;
1038 t303 = t302 / 0.6e1;
1039 t304 = t256 / 0.3e1;
1040 t305 = t82 + t198;
1041 t306 = t53 * t305;
1042 t308 = t8 * t80 * t306;
1043 t309 = t308 / 0.2e1;
1044 t310 = t130 * t130;
1045 t311 = t74 * t310;
1046 t313 = t8 * t227 * t311;
1047 t314 = t313 / 0.6e1;
1048 t315 = t279 / 0.3e1;
1049 t316 = t144 * t130;
1050 t317 = t155 * t316;
1051 t318 = t288 * t317;
1052 t319 = 0.4e1 / 0.27e2 * t318;
1053 t320 = -t305;
1054 t321 = t74 * t320;
1055 t323 = t8 * t112 * t321;
1056 t324 = t323 / 0.2e1;
1057 t325 = t294 / 0.27e2;
1058 t328 = 0.1e1 / t63 / t137 / rho[1];
1059 t333 = 0.40e2 / 0.9e1 * t61 * t328 - 0.80e2 / 0.9e1 * tau[1] * t328;
1060 t335 = t99 * t37 * t333;
1061 t336 = t136 * t335;
1062 t337 = t336 / 0.18e2;
1063 t338 = -t303 - t304 - t309 + t209 - t314 - t315 - t319 - t324 + t245 - t325 - t337;
1064 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1065 v2rho2[2] = -t128 - t174 - t133 - t176 - t147 / 0.9e1 + t10 * t338;
1066
1067 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1068 v2rhosigma[0] = 0.0e0;
1069
1070 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1071 v2rhosigma[1] = 0.0e0;
1072
1073 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1074 v2rhosigma[2] = 0.0e0;
1075
1076 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1077 v2rhosigma[3] = 0.0e0;
1078
1079 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1080 v2rhosigma[4] = 0.0e0;
1081
1082 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1083 v2rhosigma[5] = 0.0e0;
1084
1085 t340 = t97 * t190;
1086 t341 = t340 * t157;
1087 t343 = t15 * params->BB;
1088 t344 = t152 * t343;
1089 t346 = t155 * t156 * t85;
1090 t349 = t39 * t102;
1091 t350 = t155 * t349;
1092 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1093 v2rholapl[0] = -0.2e1 / 0.27e2 * t341 - 0.2e1 / 0.27e2 * t344 * t346 + 0.5e1 / 0.54e2 * t154 * t350;
1094
1095 t353 = t135 * t190;
1096 t354 = t353 * t163;
1097 t355 = t58 * params->BB;
1098 t356 = t152 * t355;
1099 t358 = t155 * t162 * t113;
1100 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1101 v2rholapl[1] = -0.2e1 / 0.27e2 * t356 * t358 - 0.2e1 / 0.27e2 * t354;
1102
1103 t362 = t155 * t156 * t125;
1104 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1105 v2rholapl[2] = -0.2e1 / 0.27e2 * t344 * t362 - 0.2e1 / 0.27e2 * t341;
1106
1107 t367 = t155 * t162 * t130;
1108 t370 = t39 * t139;
1109 t371 = t155 * t370;
1110 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1111 v2rholapl[3] = -0.2e1 / 0.27e2 * t354 - 0.2e1 / 0.27e2 * t356 * t367 + 0.5e1 / 0.54e2 * t161 * t371;
1112
1113 t374 = t98 * t168;
1114 t376 = t44 * t85;
1115 t377 = t155 * t376;
1116 t381 = t99 * t37 * t102;
1117 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1118 v2rhotau[0] = 0.4e1 / 0.27e2 * t374 + 0.4e1 / 0.27e2 * t344 * t377 - 0.5e1 / 0.27e2 * t166 * t381;
1119
1120 t384 = t136 * t172;
1121 t385 = t65 * t113;
1122 t386 = t155 * t385;
1123 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1124 v2rhotau[1] = 0.4e1 / 0.27e2 * t356 * t386 + 0.4e1 / 0.27e2 * t384;
1125
1126 t389 = t44 * t125;
1127 t390 = t155 * t389;
1128 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1129 v2rhotau[2] = 0.4e1 / 0.27e2 * t344 * t390 + 0.4e1 / 0.27e2 * t374;
1130
1131 t394 = t65 * t130;
1132 t395 = t155 * t394;
1133 t399 = t99 * t37 * t139;
1134 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1135 v2rhotau[3] = 0.4e1 / 0.27e2 * t384 + 0.4e1 / 0.27e2 * t356 * t395 - 0.5e1 / 0.27e2 * t170 * t399;
1136
1137 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1138 v2sigma2[0] = 0.0e0;
1139
1140 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1141 v2sigma2[1] = 0.0e0;
1142
1143 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1144 v2sigma2[2] = 0.0e0;
1145
1146 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1147 v2sigma2[3] = 0.0e0;
1148
1149 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1150 v2sigma2[4] = 0.0e0;
1151
1152 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1153 v2sigma2[5] = 0.0e0;
1154
1155 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1156 v2sigmalapl[0] = 0.0e0;
1157
1158 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1159 v2sigmalapl[1] = 0.0e0;
1160
1161 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1162 v2sigmalapl[2] = 0.0e0;
1163
1164 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1165 v2sigmalapl[3] = 0.0e0;
1166
1167 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1168 v2sigmalapl[4] = 0.0e0;
1169
1170 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1171 v2sigmalapl[5] = 0.0e0;
1172
1173 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1174 v2sigmatau[0] = 0.0e0;
1175
1176 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1177 v2sigmatau[1] = 0.0e0;
1178
1179 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1180 v2sigmatau[2] = 0.0e0;
1181
1182 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1183 v2sigmatau[3] = 0.0e0;
1184
1185 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1186 v2sigmatau[4] = 0.0e0;
1187
1188 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1189 v2sigmatau[5] = 0.0e0;
1190
1191 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1192 v2lapl2[0] = 0.0e0;
1193
1194 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1195 v2lapl2[1] = 0.0e0;
1196
1197 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1198 v2lapl2[2] = 0.0e0;
1199
1200 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1201 v2lapltau[0] = 0.0e0;
1202
1203 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1204 v2lapltau[1] = 0.0e0;
1205
1206 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1207 v2lapltau[2] = 0.0e0;
1208
1209 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
1210 v2lapltau[3] = 0.0e0;
1211
1212 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1213 v2tau2[0] = 0.0e0;
1214
1215 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1216 v2tau2[1] = 0.0e0;
1217
1218 if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
1219 v2tau2[2] = 0.0e0;
1220
1221 #ifndef XC_DONT_COMPILE_KXC
1222
1223 if(order < 3) return;
1224
1225
1226 t403 = t208 / 0.4e1;
1227 t407 = t244 / 0.4e1;
1228 t412 = 0.1e1 / t225 / t57;
1229 t413 = t412 * t17;
1230 t414 = t228 * t113;
1231 t415 = t74 * t414;
1232 t417 = t8 * t413 * t415;
1233 t418 = t417 / 0.9e1;
1234 t419 = t226 * t91;
1235 t421 = t8 * t419 * t229;
1236 t423 = t113 * t237;
1237 t425 = t272 * t273 * t423;
1238 t426 = t425 / 0.2e1;
1239 t427 = t58 * t205;
1240 t429 = t8 * t427 * t114;
1241 t430 = t429 / 0.3e1;
1242 t432 = 0.1e1 / t177 / t14;
1243 t433 = t432 * t17;
1244 t434 = t180 * t85;
1245 t435 = t53 * t434;
1246 t437 = t8 * t433 * t435;
1247 t438 = t437 / 0.9e1;
1248 t439 = t178 * t91;
1249 t441 = t8 * t439 * t181;
1250 t443 = t85 * t199;
1251 t445 = t249 * t250 * t443;
1252 t446 = t445 / 0.2e1;
1253 t447 = t96 * t178;
1254 t448 = t447 * t190;
1255 t450 = t155 * t107 * t180;
1256 t451 = t448 * t450;
1257 t452 = 0.2e1 / 0.27e2 * t451;
1258 t453 = t15 * t205;
1259 t455 = t8 * t453 * t86;
1260 t456 = t455 / 0.3e1;
1261 t457 = t91 * params->BB;
1262 t458 = t189 * t457;
1263 t459 = t458 * t193;
1264 t463 = t155 * t220 * t85;
1265 t464 = t191 * t463;
1266 t465 = 0.2e1 / 0.9e1 * t464;
1267 t467 = t155 * t107 * t199;
1268 t468 = t191 * t467;
1269 t469 = 0.2e1 / 0.9e1 * t468;
1270 t470 = t97 * t205;
1271 t471 = t470 * t109;
1272 t472 = t471 / 0.27e2;
1273 t473 = t210 * t222;
1274 t475 = t100 * t100;
1275 t477 = 0.1e1 / t42 / t475;
1276 t482 = -0.440e3 / 0.27e2 * t40 * t477 + 0.880e3 / 0.27e2 * tau[0] * t477;
1277 t484 = t99 * t37 * t482;
1278 t485 = t98 * t484;
1279 t486 = t485 / 0.18e2;
1280 t488 = t8 * t233 * t238;
1281 t490 = t81 * t81;
1282 t491 = 0.1e1 / t490;
1283 t492 = t9 * t491;
1284 t494 = -0.3e1 * t197 + 0.3e1 * t492;
1285 t495 = t74 * t494;
1286 t497 = t8 * t112 * t495;
1287 t498 = t497 / 0.2e1;
1288 t501 = t8 * t59 * t263 * t74;
1289 t502 = 0.5e1 / 0.36e2 * t501;
1290 t504 = t8 * t185 * t200;
1291 t506 = -t494;
1292 t507 = t53 * t506;
1293 t509 = t8 * t80 * t507;
1294 t510 = t509 / 0.2e1;
1295 t513 = t8 * t16 * t263 * t53;
1296 t514 = 0.5e1 / 0.36e2 * t513;
1297 t515 = -t465 - t469 + t472 - t473 / 0.18e2 - t486 - t488 / 0.2e1 - t498 - t502 - t504 / 0.2e1 - t510 - t514;
1298 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1299 v3rho3[0] = -0.3e1 / 0.2e1 * t202 + t403 - t211 / 0.9e1 - t223 / 0.6e1 - 0.3e1 / 0.2e1 * t240 + t407 - t231 / 0.2e1 - t235 - t183 / 0.2e1 - t187 - 0.4e1 / 0.9e1 * t194 + t10 * (t418 - t421 / 0.6e1 - t426 + t430 + t438 - t441 / 0.6e1 - t446 - t452 + t456 - 0.4e1 / 0.27e2 * t459 + t515);
1300
1301 t519 = t276 / 0.3e1;
1302 t520 = 0.4e1 / 0.27e2 * t291;
1303 t523 = t253 / 0.3e1;
1304 t524 = 0.4e1 / 0.27e2 * t260;
1305 t526 = t130 * t237;
1306 t528 = t272 * t273 * t526;
1307 t530 = t263 * t74;
1308 t531 = t9 * t113;
1309 t533 = t272 * t530 * t531;
1310 t535 = t125 * t199;
1311 t537 = t249 * t250 * t535;
1312 t539 = t263 * params->BB;
1313 t540 = t189 * t539;
1314 t541 = t107 * t9;
1315 t542 = t155 * t541;
1316 t543 = t540 * t542;
1317 t546 = t5 * t7 * t412;
1318 t547 = t130 * t228;
1319 t549 = t546 * t273 * t547;
1320 t551 = t91 * t74;
1321 t553 = t272 * t551 * t274;
1322 t554 = t553 / 0.9e1;
1323 t555 = t96 * t226;
1324 t556 = t555 * t190;
1325 t558 = t155 * t144 * t228;
1326 t559 = t556 * t558;
1327 t561 = t287 * t457;
1328 t562 = t561 * t290;
1329 t563 = 0.4e1 / 0.81e2 * t562;
1330 t565 = t5 * t7 * t432;
1331 t566 = t125 * t180;
1332 t568 = t565 * t250 * t566;
1333 t570 = t91 * t53;
1334 t572 = t249 * t570 * t251;
1335 t573 = t572 / 0.9e1;
1336 t574 = t263 * t53;
1337 t575 = t9 * t85;
1338 t577 = t249 * t574 * t575;
1339 t579 = t144 * t237;
1340 t580 = t155 * t579;
1341 t581 = t288 * t580;
1342 t583 = t458 * t259;
1343 t584 = 0.4e1 / 0.81e2 * t583;
1344 t585 = t220 * t125;
1345 t586 = t155 * t585;
1346 t587 = t191 * t586;
1347 t591 = t8 * t264 * t53;
1348 t593 = -t528 / 0.6e1 + t533 / 0.3e1 - t537 / 0.6e1 - 0.4e1 / 0.27e2 * t543 + t549 / 0.9e1 - t554 - 0.2e1 / 0.81e2 * t559 - t563 + t568 / 0.9e1 - t573 - t577 / 0.3e1 - 0.2e1 / 0.27e2 * t581 - t584 - 0.2e1 / 0.27e2 * t587 - 0.4e1 / 0.81e2 * t459 - t591 / 0.2e1;
1349 t595 = t8 * t281 * t74;
1350 t598 = 0.1e1 / t90 / t196;
1351 t599 = t58 * t598;
1352 t601 = t8 * t599 * t282;
1353 t603 = t135 * t205;
1354 t604 = t603 * t146;
1355 t607 = t8 * t453 * t126;
1356 t609 = t15 * t598;
1357 t611 = t8 * t609 * t265;
1358 t614 = t8 * t427 * t131;
1359 t625 = t155 * t258 * t85;
1360 t626 = t448 * t625;
1361 t628 = t595 / 0.2e1 - t502 - t514 - 0.7e1 / 0.6e1 * t601 + t604 / 0.81e2 + t607 / 0.9e1 + 0.7e1 / 0.6e1 * t611 + t614 / 0.9e1 + 0.2e1 / 0.9e1 * t429 - t441 / 0.18e2 + 0.2e1 / 0.9e1 * t455 + 0.2e1 / 0.81e2 * t471 - t473 / 0.54e2 - t488 / 0.6e1 - t504 / 0.6e1 - t421 / 0.18e2 - 0.4e1 / 0.81e2 * t626;
1362 t631 = -t325 - t232 - 0.2e1 / 0.3e1 * t235 - t523 - t304 - t524 - t267 - t184 - 0.2e1 / 0.3e1 * t187 - t195 + t10 * (t593 + t628);
1363 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1364 v3rho3[1] = t284 - t203 + t403 - 0.2e1 / 0.27e2 * t211 - t224 - t241 + t407 - t519 - t315 - t520 + t631;
1365
1366 t636 = t320 * t113;
1367 t638 = t272 * t273 * t636;
1368 t641 = t155 * t333 * t113;
1369 t642 = t288 * t641;
1370 t644 = t125 * t9;
1371 t646 = t249 * t574 * t644;
1372 t648 = t561 * t317;
1373 t650 = t299 * t85;
1374 t652 = t565 * t250 * t650;
1375 t654 = t107 * t299;
1376 t655 = t155 * t654;
1377 t656 = t448 * t655;
1378 t658 = t305 * t85;
1379 t660 = t249 * t250 * t658;
1380 t662 = t107 * t305;
1381 t663 = t155 * t662;
1382 t664 = t191 * t663;
1383 t666 = t310 * t113;
1384 t668 = t546 * t273 * t666;
1385 t670 = t130 * t9;
1386 t672 = t272 * t530 * t670;
1387 t674 = t287 * t539;
1388 t675 = t144 * t9;
1389 t676 = t155 * t675;
1390 t677 = t674 * t676;
1391 t679 = -t638 / 0.6e1 - 0.2e1 / 0.27e2 * t642 - t646 / 0.3e1 - 0.4e1 / 0.81e2 * t648 + t652 / 0.9e1 - 0.2e1 / 0.81e2 * t656 - t660 / 0.6e1 - 0.2e1 / 0.27e2 * t664 + t668 / 0.9e1 + t672 / 0.3e1 + 0.4e1 / 0.27e2 * t677 - t554 - t563 - t573 - t584 - t502;
1392 t680 = t293 * t335;
1393 t683 = -t197 - 0.3e1 * t492;
1394 t684 = t53 * t683;
1395 t686 = t8 * t80 * t684;
1396 t688 = -t683;
1397 t689 = t74 * t688;
1398 t691 = t8 * t112 * t689;
1399 t694 = t8 * t439 * t300;
1400 t697 = t8 * t185 * t306;
1401 t700 = t8 * t419 * t311;
1402 t703 = t8 * t233 * t321;
1403 t714 = t155 * t316 * t113;
1404 t715 = t556 * t714;
1405 t717 = -t514 - t680 / 0.54e2 - t686 / 0.2e1 - t691 / 0.2e1 - t694 / 0.18e2 - t697 / 0.6e1 - t700 / 0.18e2 - t703 / 0.6e1 + t601 / 0.3e1 + 0.2e1 / 0.81e2 * t604 + 0.2e1 / 0.9e1 * t607 - t611 / 0.3e1 + 0.2e1 / 0.9e1 * t614 + t429 / 0.9e1 + t455 / 0.9e1 + t471 / 0.81e2 - 0.4e1 / 0.81e2 * t715;
1406 t720 = -t520 - 0.2e1 / 0.27e2 * t294 - t303 - t314 - t236 - t523 - 0.2e1 / 0.3e1 * t256 - t524 - t267 - t188 + t10 * (t679 + t717);
1407 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1408 v3rho3[2] = -t337 + t284 - t309 - t324 + t403 - t212 + t407 - t319 - t519 - 0.2e1 / 0.3e1 * t279 + t720;
1409
1410 t728 = t604 / 0.27e2;
1411 t734 = t299 * t125;
1412 t735 = t53 * t734;
1413 t737 = t8 * t433 * t735;
1414 t738 = t737 / 0.9e1;
1415 t739 = t125 * t305;
1416 t741 = t249 * t250 * t739;
1417 t742 = t741 / 0.2e1;
1418 t743 = t310 * t130;
1419 t744 = t74 * t743;
1420 t746 = t8 * t413 * t744;
1421 t747 = t746 / 0.9e1;
1422 t748 = t130 * t320;
1423 t750 = t272 * t273 * t748;
1424 t751 = t750 / 0.2e1;
1425 t753 = t144 * t310;
1426 t754 = t155 * t753;
1427 t755 = t556 * t754;
1428 t756 = 0.2e1 / 0.27e2 * t755;
1429 t757 = t333 * t130;
1430 t758 = t155 * t757;
1431 t759 = t288 * t758;
1432 t760 = 0.2e1 / 0.9e1 * t759;
1433 t761 = t144 * t320;
1434 t762 = t155 * t761;
1435 t763 = t288 * t762;
1436 t764 = 0.2e1 / 0.9e1 * t763;
1437 t765 = t607 / 0.3e1;
1438 t766 = t614 / 0.3e1;
1439 t768 = t137 * t137;
1440 t770 = 0.1e1 / t63 / t768;
1441 t775 = -0.440e3 / 0.27e2 * t61 * t770 + 0.880e3 / 0.27e2 * tau[1] * t770;
1442 t777 = t99 * t37 * t775;
1443 t778 = t136 * t777;
1444 t779 = t778 / 0.18e2;
1445 t781 = -0.3e1 * t197 - 0.3e1 * t492;
1446 t782 = t53 * t781;
1447 t784 = t8 * t80 * t782;
1448 t785 = t784 / 0.2e1;
1449 t786 = -t781;
1450 t787 = t74 * t786;
1451 t789 = t8 * t112 * t787;
1452 t790 = t789 / 0.2e1;
1453 t791 = -t756 - t760 - t764 + t765 + t766 - 0.4e1 / 0.27e2 * t648 - t502 - t514 - t779 - t785 - t790;
1454 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1455 v3rho3[3] = -t336 / 0.6e1 - 0.3e1 / 0.2e1 * t308 - 0.3e1 / 0.2e1 * t323 + t403 + t407 - 0.4e1 / 0.9e1 * t318 - t279 - t294 / 0.9e1 - t302 / 0.2e1 - t313 / 0.2e1 - t256 + t10 * (t728 - t694 / 0.6e1 - t697 / 0.2e1 - t700 / 0.6e1 - t703 / 0.2e1 - t680 / 0.18e2 + t738 - t742 + t747 - t751 + t791);
1456
1457 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1458 v3rho2sigma[0] = 0.0e0;
1459
1460 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1461 v3rho2sigma[1] = 0.0e0;
1462
1463 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1464 v3rho2sigma[2] = 0.0e0;
1465
1466 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1467 v3rho2sigma[3] = 0.0e0;
1468
1469 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1470 v3rho2sigma[4] = 0.0e0;
1471
1472 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1473 v3rho2sigma[5] = 0.0e0;
1474
1475 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1476 v3rho2sigma[6] = 0.0e0;
1477
1478 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1479 v3rho2sigma[7] = 0.0e0;
1480
1481 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1482 v3rho2sigma[8] = 0.0e0;
1483
1484 t794 = t191 * t346;
1485 t796 = t97 * t457;
1486 t798 = 0.2e1 / 0.81e2 * t796 * t157;
1487 t799 = t340 * t350;
1488 t801 = t178 * params->BB;
1489 t802 = t152 * t801;
1490 t804 = t155 * t156 * t180;
1491 t808 = t155 * t349 * t85;
1492 t812 = t155 * t156 * t199;
1493 t815 = t39 * t215;
1494 t816 = t155 * t815;
1495 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1496 v3rho2lapl[0] = -0.16e2 / 0.81e2 * t794 - t798 + 0.20e2 / 0.81e2 * t799 - 0.2e1 / 0.81e2 * t802 * t804 + 0.20e2 / 0.81e2 * t344 * t808 - 0.2e1 / 0.27e2 * t344 * t812 - 0.20e2 / 0.81e2 * t154 * t816;
1497
1498 t819 = t288 * t358;
1499 t821 = t135 * t457;
1500 t823 = 0.2e1 / 0.81e2 * t821 * t163;
1501 t824 = t226 * params->BB;
1502 t825 = t152 * t824;
1503 t827 = t155 * t162 * t228;
1504 t831 = t155 * t162 * t237;
1505 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1506 v3rho2lapl[1] = -0.16e2 / 0.81e2 * t819 - t823 - 0.2e1 / 0.81e2 * t825 * t827 - 0.2e1 / 0.27e2 * t356 * t831;
1507
1508 t836 = t191 * t362;
1509 t838 = t801 * t34;
1510 t839 = t152 * t838;
1511 t840 = t37 * t39;
1512 t841 = t389 * t85;
1513 t842 = t840 * t841;
1514 t846 = t155 * t349 * t125;
1515 t849 = t205 * t96;
1516 t850 = t849 * t343;
1517 t852 = t155 * t156 * t9;
1518 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1519 v3rho2lapl[2] = -0.8e1 / 0.81e2 * t794 - t798 + 0.10e2 / 0.81e2 * t799 - 0.8e1 / 0.81e2 * t836 - 0.2e1 / 0.81e2 * t839 * t842 + 0.10e2 / 0.81e2 * t344 * t846 - 0.2e1 / 0.27e2 * t850 * t852;
1520
1521 t856 = t288 * t367;
1522 t858 = t824 * t34;
1523 t859 = t152 * t858;
1524 t860 = t394 * t113;
1525 t861 = t840 * t860;
1526 t864 = t849 * t355;
1527 t866 = t155 * t162 * t9;
1528 t869 = t353 * t371;
1529 t872 = t155 * t370 * t113;
1530 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1531 v3rho2lapl[3] = -0.8e1 / 0.81e2 * t819 - t823 - 0.8e1 / 0.81e2 * t856 - 0.2e1 / 0.81e2 * t859 * t861 + 0.2e1 / 0.27e2 * t864 * t866 + 0.10e2 / 0.81e2 * t869 + 0.10e2 / 0.81e2 * t356 * t872;
1532
1533 t877 = t155 * t156 * t299;
1534 t881 = t155 * t156 * t305;
1535 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1536 v3rho2lapl[4] = -0.16e2 / 0.81e2 * t836 - t798 - 0.2e1 / 0.81e2 * t802 * t877 - 0.2e1 / 0.27e2 * t344 * t881;
1537
1538 t887 = t155 * t162 * t310;
1539 t891 = t155 * t370 * t130;
1540 t895 = t155 * t162 * t320;
1541 t898 = t39 * t328;
1542 t899 = t155 * t898;
1543 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1544 v3rho2lapl[5] = -0.16e2 / 0.81e2 * t856 - t823 + 0.20e2 / 0.81e2 * t869 - 0.2e1 / 0.81e2 * t825 * t887 + 0.20e2 / 0.81e2 * t356 * t891 - 0.2e1 / 0.27e2 * t356 * t895 - 0.20e2 / 0.81e2 * t161 * t899;
1545
1546 t902 = t191 * t377;
1547 t905 = 0.4e1 / 0.81e2 * t210 * t168;
1548 t906 = t98 * t381;
1549 t909 = t155 * t44 * t180;
1550 t913 = t155 * t102 * t85;
1551 t917 = t155 * t44 * t199;
1552 t921 = t99 * t37 * t215;
1553 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1554 v3rho2tau[0] = 0.32e2 / 0.81e2 * t902 + t905 - 0.40e2 / 0.81e2 * t906 + 0.4e1 / 0.81e2 * t802 * t909 - 0.40e2 / 0.81e2 * t344 * t913 + 0.4e1 / 0.27e2 * t344 * t917 + 0.40e2 / 0.81e2 * t166 * t921;
1555
1556 t924 = t288 * t386;
1557 t927 = 0.4e1 / 0.81e2 * t293 * t172;
1558 t929 = t155 * t65 * t228;
1559 t933 = t155 * t65 * t237;
1560 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1561 v3rho2tau[1] = 0.32e2 / 0.81e2 * t924 + t927 + 0.4e1 / 0.81e2 * t825 * t929 + 0.4e1 / 0.27e2 * t356 * t933;
1562
1563 t938 = t191 * t390;
1564 t940 = t155 * t841;
1565 t943 = t102 * t125;
1566 t944 = t155 * t943;
1567 t947 = t44 * t9;
1568 t948 = t155 * t947;
1569 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1570 v3rho2tau[2] = 0.16e2 / 0.81e2 * t902 + t905 - 0.20e2 / 0.81e2 * t906 + 0.16e2 / 0.81e2 * t938 + 0.4e1 / 0.81e2 * t802 * t940 - 0.20e2 / 0.81e2 * t344 * t944 + 0.4e1 / 0.27e2 * t850 * t948;
1571
1572 t952 = t288 * t395;
1573 t954 = t155 * t860;
1574 t957 = t65 * t9;
1575 t958 = t155 * t957;
1576 t961 = t136 * t399;
1577 t964 = t155 * t139 * t113;
1578 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1579 v3rho2tau[3] = 0.16e2 / 0.81e2 * t924 + t927 + 0.16e2 / 0.81e2 * t952 + 0.4e1 / 0.81e2 * t825 * t954 - 0.4e1 / 0.27e2 * t864 * t958 - 0.20e2 / 0.81e2 * t961 - 0.20e2 / 0.81e2 * t356 * t964;
1580
1581 t968 = t44 * t299;
1582 t969 = t155 * t968;
1583 t972 = t44 * t305;
1584 t973 = t155 * t972;
1585 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1586 v3rho2tau[4] = 0.32e2 / 0.81e2 * t938 + t905 + 0.4e1 / 0.81e2 * t802 * t969 + 0.4e1 / 0.27e2 * t344 * t973;
1587
1588 t978 = t65 * t310;
1589 t979 = t155 * t978;
1590 t982 = t139 * t130;
1591 t983 = t155 * t982;
1592 t986 = t65 * t320;
1593 t987 = t155 * t986;
1594 t991 = t99 * t37 * t328;
1595 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1596 v3rho2tau[5] = 0.32e2 / 0.81e2 * t952 + t927 - 0.40e2 / 0.81e2 * t961 + 0.4e1 / 0.81e2 * t825 * t979 - 0.40e2 / 0.81e2 * t356 * t983 + 0.4e1 / 0.27e2 * t356 * t987 + 0.40e2 / 0.81e2 * t170 * t991;
1597
1598 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1599 v3rhosigma2[0] = 0.0e0;
1600
1601 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1602 v3rhosigma2[1] = 0.0e0;
1603
1604 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1605 v3rhosigma2[2] = 0.0e0;
1606
1607 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1608 v3rhosigma2[3] = 0.0e0;
1609
1610 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1611 v3rhosigma2[4] = 0.0e0;
1612
1613 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1614 v3rhosigma2[5] = 0.0e0;
1615
1616 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1617 v3rhosigma2[6] = 0.0e0;
1618
1619 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1620 v3rhosigma2[7] = 0.0e0;
1621
1622 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1623 v3rhosigma2[8] = 0.0e0;
1624
1625 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1626 v3rhosigma2[9] = 0.0e0;
1627
1628 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1629 v3rhosigma2[10] = 0.0e0;
1630
1631 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1632 v3rhosigma2[11] = 0.0e0;
1633
1634 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1635 v3rhosigmalapl[0] = 0.0e0;
1636
1637 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1638 v3rhosigmalapl[1] = 0.0e0;
1639
1640 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1641 v3rhosigmalapl[2] = 0.0e0;
1642
1643 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1644 v3rhosigmalapl[3] = 0.0e0;
1645
1646 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1647 v3rhosigmalapl[4] = 0.0e0;
1648
1649 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1650 v3rhosigmalapl[5] = 0.0e0;
1651
1652 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1653 v3rhosigmalapl[6] = 0.0e0;
1654
1655 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1656 v3rhosigmalapl[7] = 0.0e0;
1657
1658 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1659 v3rhosigmalapl[8] = 0.0e0;
1660
1661 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1662 v3rhosigmalapl[9] = 0.0e0;
1663
1664 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1665 v3rhosigmalapl[10] = 0.0e0;
1666
1667 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1668 v3rhosigmalapl[11] = 0.0e0;
1669
1670 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1671 v3rhosigmatau[0] = 0.0e0;
1672
1673 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1674 v3rhosigmatau[1] = 0.0e0;
1675
1676 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1677 v3rhosigmatau[2] = 0.0e0;
1678
1679 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1680 v3rhosigmatau[3] = 0.0e0;
1681
1682 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1683 v3rhosigmatau[4] = 0.0e0;
1684
1685 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1686 v3rhosigmatau[5] = 0.0e0;
1687
1688 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1689 v3rhosigmatau[6] = 0.0e0;
1690
1691 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1692 v3rhosigmatau[7] = 0.0e0;
1693
1694 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1695 v3rhosigmatau[8] = 0.0e0;
1696
1697 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1698 v3rhosigmatau[9] = 0.0e0;
1699
1700 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1701 v3rhosigmatau[10] = 0.0e0;
1702
1703 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1704 v3rhosigmatau[11] = 0.0e0;
1705
1706 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1707 v3rholapl2[0] = 0.0e0;
1708
1709 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1710 v3rholapl2[1] = 0.0e0;
1711
1712 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1713 v3rholapl2[2] = 0.0e0;
1714
1715 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1716 v3rholapl2[3] = 0.0e0;
1717
1718 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1719 v3rholapl2[4] = 0.0e0;
1720
1721 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1722 v3rholapl2[5] = 0.0e0;
1723
1724 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1725 v3rholapltau[0] = 0.0e0;
1726
1727 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1728 v3rholapltau[1] = 0.0e0;
1729
1730 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1731 v3rholapltau[2] = 0.0e0;
1732
1733 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1734 v3rholapltau[3] = 0.0e0;
1735
1736 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1737 v3rholapltau[4] = 0.0e0;
1738
1739 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1740 v3rholapltau[5] = 0.0e0;
1741
1742 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1743 v3rholapltau[6] = 0.0e0;
1744
1745 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1746 v3rholapltau[7] = 0.0e0;
1747
1748 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1749 v3rhotau2[0] = 0.0e0;
1750
1751 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1752 v3rhotau2[1] = 0.0e0;
1753
1754 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1755 v3rhotau2[2] = 0.0e0;
1756
1757 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1758 v3rhotau2[3] = 0.0e0;
1759
1760 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1761 v3rhotau2[4] = 0.0e0;
1762
1763 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1764 v3rhotau2[5] = 0.0e0;
1765
1766 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1767 v3sigma3[0] = 0.0e0;
1768
1769 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1770 v3sigma3[1] = 0.0e0;
1771
1772 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1773 v3sigma3[2] = 0.0e0;
1774
1775 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1776 v3sigma3[3] = 0.0e0;
1777
1778 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1779 v3sigma3[4] = 0.0e0;
1780
1781 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1782 v3sigma3[5] = 0.0e0;
1783
1784 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1785 v3sigma3[6] = 0.0e0;
1786
1787 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1788 v3sigma3[7] = 0.0e0;
1789
1790 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1791 v3sigma3[8] = 0.0e0;
1792
1793 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1794 v3sigma3[9] = 0.0e0;
1795
1796 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1797 v3sigma2lapl[0] = 0.0e0;
1798
1799 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1800 v3sigma2lapl[1] = 0.0e0;
1801
1802 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1803 v3sigma2lapl[2] = 0.0e0;
1804
1805 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1806 v3sigma2lapl[3] = 0.0e0;
1807
1808 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1809 v3sigma2lapl[4] = 0.0e0;
1810
1811 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1812 v3sigma2lapl[5] = 0.0e0;
1813
1814 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1815 v3sigma2lapl[6] = 0.0e0;
1816
1817 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1818 v3sigma2lapl[7] = 0.0e0;
1819
1820 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1821 v3sigma2lapl[8] = 0.0e0;
1822
1823 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1824 v3sigma2lapl[9] = 0.0e0;
1825
1826 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1827 v3sigma2lapl[10] = 0.0e0;
1828
1829 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1830 v3sigma2lapl[11] = 0.0e0;
1831
1832 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1833 v3sigma2tau[0] = 0.0e0;
1834
1835 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1836 v3sigma2tau[1] = 0.0e0;
1837
1838 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1839 v3sigma2tau[2] = 0.0e0;
1840
1841 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1842 v3sigma2tau[3] = 0.0e0;
1843
1844 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1845 v3sigma2tau[4] = 0.0e0;
1846
1847 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1848 v3sigma2tau[5] = 0.0e0;
1849
1850 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1851 v3sigma2tau[6] = 0.0e0;
1852
1853 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1854 v3sigma2tau[7] = 0.0e0;
1855
1856 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1857 v3sigma2tau[8] = 0.0e0;
1858
1859 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1860 v3sigma2tau[9] = 0.0e0;
1861
1862 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1863 v3sigma2tau[10] = 0.0e0;
1864
1865 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1866 v3sigma2tau[11] = 0.0e0;
1867
1868 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1869 v3sigmalapl2[0] = 0.0e0;
1870
1871 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1872 v3sigmalapl2[1] = 0.0e0;
1873
1874 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1875 v3sigmalapl2[2] = 0.0e0;
1876
1877 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1878 v3sigmalapl2[3] = 0.0e0;
1879
1880 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1881 v3sigmalapl2[4] = 0.0e0;
1882
1883 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1884 v3sigmalapl2[5] = 0.0e0;
1885
1886 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1887 v3sigmalapl2[6] = 0.0e0;
1888
1889 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1890 v3sigmalapl2[7] = 0.0e0;
1891
1892 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1893 v3sigmalapl2[8] = 0.0e0;
1894
1895 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1896 v3sigmalapltau[0] = 0.0e0;
1897
1898 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1899 v3sigmalapltau[1] = 0.0e0;
1900
1901 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1902 v3sigmalapltau[2] = 0.0e0;
1903
1904 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1905 v3sigmalapltau[3] = 0.0e0;
1906
1907 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1908 v3sigmalapltau[4] = 0.0e0;
1909
1910 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1911 v3sigmalapltau[5] = 0.0e0;
1912
1913 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1914 v3sigmalapltau[6] = 0.0e0;
1915
1916 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1917 v3sigmalapltau[7] = 0.0e0;
1918
1919 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1920 v3sigmalapltau[8] = 0.0e0;
1921
1922 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1923 v3sigmalapltau[9] = 0.0e0;
1924
1925 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1926 v3sigmalapltau[10] = 0.0e0;
1927
1928 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1929 v3sigmalapltau[11] = 0.0e0;
1930
1931 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1932 v3sigmatau2[0] = 0.0e0;
1933
1934 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1935 v3sigmatau2[1] = 0.0e0;
1936
1937 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1938 v3sigmatau2[2] = 0.0e0;
1939
1940 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1941 v3sigmatau2[3] = 0.0e0;
1942
1943 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1944 v3sigmatau2[4] = 0.0e0;
1945
1946 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1947 v3sigmatau2[5] = 0.0e0;
1948
1949 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1950 v3sigmatau2[6] = 0.0e0;
1951
1952 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1953 v3sigmatau2[7] = 0.0e0;
1954
1955 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
1956 v3sigmatau2[8] = 0.0e0;
1957
1958 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1959 v3lapl3[0] = 0.0e0;
1960
1961 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1962 v3lapl3[1] = 0.0e0;
1963
1964 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1965 v3lapl3[2] = 0.0e0;
1966
1967 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1968 v3lapl3[3] = 0.0e0;
1969
1970 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1971 v3lapl2tau[0] = 0.0e0;
1972
1973 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1974 v3lapl2tau[1] = 0.0e0;
1975
1976 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1977 v3lapl2tau[2] = 0.0e0;
1978
1979 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1980 v3lapl2tau[3] = 0.0e0;
1981
1982 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1983 v3lapl2tau[4] = 0.0e0;
1984
1985 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1986 v3lapl2tau[5] = 0.0e0;
1987
1988 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1989 v3lapltau2[0] = 0.0e0;
1990
1991 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1992 v3lapltau2[1] = 0.0e0;
1993
1994 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1995 v3lapltau2[2] = 0.0e0;
1996
1997 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
1998 v3lapltau2[3] = 0.0e0;
1999
2000 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
2001 v3lapltau2[4] = 0.0e0;
2002
2003 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
2004 v3lapltau2[5] = 0.0e0;
2005
2006 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
2007 v3tau3[0] = 0.0e0;
2008
2009 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
2010 v3tau3[1] = 0.0e0;
2011
2012 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
2013 v3tau3[2] = 0.0e0;
2014
2015 if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
2016 v3tau3[3] = 0.0e0;
2017
2018 #ifndef XC_DONT_COMPILE_LXC
2019
2020 if(order < 4) return;
2021
2022
2023 t1007 = t272 * t551 * t423;
2024 t1009 = t850 * t193;
2025 t1011 = t458 * t463;
2026 t1013 = t458 * t467;
2027 t1036 = t96 * t432 * t190;
2028 t1042 = t249 * t570 * t443;
2029 t1044 = t447 * t457;
2030 t1045 = t1044 * t450;
2031 t1052 = t8 * t185 * t507;
2032 t1054 = t178 * t205;
2033 t1056 = t8 * t1054 * t181;
2034 t1058 = t199 * t199;
2035 t1063 = -0.2e1 / 0.3e1 * t249 * t250 * t85 * t506 - 0.4e1 / 0.27e2 * t448 * t155 * t220 * t180 + 0.2e1 / 0.3e1 * t546 * t273 * t228 * t237 - 0.2e1 / 0.3e1 * t1007 + 0.16e2 / 0.81e2 * t1009 - 0.8e1 / 0.27e2 * t1011 - 0.8e1 / 0.27e2 * t1013 - 0.8e1 / 0.27e2 * t191 * t155 * t482 * t85 - 0.4e1 / 0.9e1 * t191 * t155 * t220 * t199 - 0.8e1 / 0.27e2 * t191 * t155 * t107 * t506 - 0.2e1 / 0.3e1 * t272 * t273 * t113 * t494 + 0.2e1 / 0.3e1 * t565 * t250 * t180 * t199 + 0.16e2 / 0.243e3 * t1036 * t155 * t107 * t434 - 0.2e1 / 0.3e1 * t1042 - 0.8e1 / 0.81e2 * t1045 - 0.8e1 / 0.27e2 * t448 * t155 * t192 * t199 - 0.2e1 / 0.3e1 * t1052 + 0.2e1 / 0.9e1 * t1056 - t8 * t179 * t53 * t1058 / 0.2e1;
2036 t1064 = t57 * t57;
2037 t1066 = 0.1e1 / t225 / t1064;
2038 t1067 = t1066 * t17;
2039 t1068 = t228 * t228;
2040 t1073 = t412 * t91;
2041 t1075 = t8 * t1073 * t415;
2042 t1077 = t226 * t205;
2043 t1079 = t8 * t1077 * t229;
2044 t1081 = t237 * t237;
2045 t1086 = t490 * t10;
2046 t1088 = t9 / t1086;
2047 t1090 = -0.12e2 * t491 + 0.12e2 * t1088;
2048 t1096 = t8 * t427 * t238;
2049 t1099 = t8 * t453 * t200;
2050 t1101 = t470 * t222;
2051 t1103 = t210 * t484;
2052 t1107 = 0.1e1 / t42 / t475 / rho[0];
2053 t1118 = t8 * t233 * t495;
2054 t1126 = t8 * t264 * t86;
2055 t1129 = t97 * t263 * t109;
2056 t1132 = t8 * t281 * t114;
2057 t1134 = t14 * t14;
2058 t1136 = 0.1e1 / t177 / t1134;
2059 t1137 = t1136 * t17;
2060 t1138 = t180 * t180;
2061 t1143 = t432 * t91;
2062 t1145 = t8 * t1143 * t435;
2063 t1150 = 0.10e2 / 0.27e2 * t8 * t16 * t598 * t53;
2064 t1154 = 0.10e2 / 0.27e2 * t8 * t59 * t598 * t74;
2065 t1155 = -0.5e1 / 0.27e2 * t8 * t1067 * t74 * t1068 + 0.4e1 / 0.27e2 * t1075 + 0.2e1 / 0.9e1 * t1079 - t8 * t227 * t74 * t1081 / 0.2e1 - t8 * t80 * t53 * t1090 / 0.2e1 + 0.2e1 / 0.3e1 * t1096 + 0.2e1 / 0.3e1 * t1099 + 0.2e1 / 0.27e2 * t1101 - 0.2e1 / 0.27e2 * t1103 - t98 * t99 * t37 * (0.6160e4 / 0.81e2 * t40 * t1107 - 0.12320e5 / 0.81e2 * tau[0] * t1107) / 0.18e2 - 0.2e1 / 0.3e1 * t1118 + t8 * t112 * t74 * t1090 / 0.2e1 - 0.20e2 / 0.27e2 * t1126 - 0.20e2 / 0.243e3 * t1129 - 0.20e2 / 0.27e2 * t1132 - 0.5e1 / 0.27e2 * t8 * t1137 * t53 * t1138 + 0.4e1 / 0.27e2 * t1145 + t1150 + t1154;
2066 t1168 = t10 * (t1063 + t1155) + 0.4e1 / 0.9e1 * t417 - 0.2e1 / 0.3e1 * t421 - 0.2e1 * t425 + 0.4e1 / 0.3e1 * t429 + 0.4e1 / 0.9e1 * t437 - 0.2e1 / 0.3e1 * t441 - 0.2e1 * t445 - 0.8e1 / 0.27e2 * t451 + 0.4e1 / 0.3e1 * t455 - 0.16e2 / 0.27e2 * t459;
2067 t1176 = 0.5e1 / 0.9e1 * t501;
2068 t1179 = 0.5e1 / 0.9e1 * t513;
2069 t1180 = -0.8e1 / 0.9e1 * t464 - 0.8e1 / 0.9e1 * t468 + 0.4e1 / 0.27e2 * t471 - 0.2e1 / 0.9e1 * t473 - 0.2e1 / 0.9e1 * t485 - 0.2e1 * t488 - 0.2e1 * t497 - t1176 - 0.2e1 * t504 - 0.2e1 * t509 - t1179;
2070 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2071 v4rho4[0] = t1168 + t1180;
2072
2073 t1182 = t272 * t551 * t526;
2074 t1188 = t598 * t74;
2075 t1190 = t272 * t1188 * t531;
2076 t1192 = t561 * t580;
2077 t1195 = t96 * t412 * t190;
2078 t1200 = t555 * t457;
2079 t1201 = t1200 * t558;
2080 t1204 = t5 * t7 * t1136;
2081 t1210 = t565 * t570 * t566;
2082 t1217 = t249 * t570 * t535;
2083 t1224 = t5 * t7 * t1066;
2084 t1230 = t546 * t551 * t547;
2085 t1232 = t598 * params->BB;
2086 t1234 = t189 * t1232 * t542;
2087 t1242 = t249 * t205 * t53 * t251;
2088 t1243 = t1242 / 0.9e1;
2089 t1244 = -t1182 / 0.6e1 - t546 * t530 * t9 * t228 / 0.3e1 - 0.7e1 / 0.6e1 * t1190 - 0.2e1 / 0.27e2 * t1192 + 0.4e1 / 0.243e3 * t1195 * t155 * t144 * t414 - 0.2e1 / 0.81e2 * t1201 - 0.5e1 / 0.27e2 * t1204 * t250 * t125 * t434 + t1210 / 0.9e1 + t565 * t574 * t9 * t180 / 0.3e1 - t1217 / 0.6e1 - t249 * t574 * t9 * t199 / 0.2e1 - 0.5e1 / 0.27e2 * t1224 * t273 * t130 * t414 + t1230 / 0.9e1 + 0.14e2 / 0.27e2 * t1234 - 0.2e1 / 0.9e1 * t540 * t155 * t220 * t9 + t1243;
2090 t1247 = t272 * t205 * t74 * t274;
2091 t1248 = t1247 / 0.9e1;
2092 t1253 = t850 * t259;
2093 t1254 = 0.4e1 / 0.81e2 * t1253;
2094 t1255 = t458 * t586;
2095 t1257 = t598 * t53;
2096 t1259 = t249 * t1257 * t575;
2097 t1277 = t864 * t290;
2098 t1278 = 0.4e1 / 0.81e2 * t1277;
2099 t1285 = t1248 - t249 * t250 * t125 * t506 / 0.6e1 + t1254 - 0.2e1 / 0.27e2 * t1255 + 0.7e1 / 0.6e1 * t1259 - 0.2e1 / 0.27e2 * t191 * t155 * t482 * t125 - t272 * t273 * t130 * t494 / 0.6e1 + t272 * t530 * t9 * t237 / 0.2e1 - 0.2e1 / 0.27e2 * t288 * t155 * t144 * t494 + t1278 - t1007 / 0.6e1 + 0.8e1 / 0.81e2 * t1009 - 0.2e1 / 0.27e2 * t1011 - 0.2e1 / 0.27e2 * t1013 - t1042 / 0.6e1 - 0.2e1 / 0.81e2 * t1045;
2100 t1287 = t447 * t539;
2101 t1292 = t1044 * t625;
2102 t1319 = 0.1e1 / t90 / t490;
2103 t1322 = t8 * t58 * t1319 * t282;
2104 t1325 = t8 * t281 * t131;
2105 t1331 = t8 * t264 * t126;
2106 t1335 = t8 * t15 * t1319 * t265;
2107 t1337 = t226 * t263;
2108 t1342 = t135 * t263 * t146;
2109 t1345 = -0.4e1 / 0.27e2 * t1287 * t155 * t541 * t85 - 0.4e1 / 0.81e2 * t1292 - 0.2e1 / 0.27e2 * t448 * t155 * t585 * t85 + t546 * t273 * t526 * t113 / 0.3e1 - 0.2e1 / 0.27e2 * t556 * t155 * t579 * t113 + 0.4e1 / 0.81e2 * t1036 * t155 * t258 * t180 + t565 * t250 * t535 * t85 / 0.3e1 - 0.2e1 / 0.27e2 * t448 * t155 * t258 * t199 + 0.25e2 / 0.6e1 * t1322 - 0.5e1 / 0.27e2 * t1325 - 0.2e1 / 0.9e1 * t189 * t263 * t109 - 0.5e1 / 0.27e2 * t1331 - 0.25e2 / 0.6e1 * t1335 + t8 * t1337 * t114 / 0.2e1 - 0.5e1 / 0.243e3 * t1342 - t1052 / 0.6e1;
2110 t1346 = t178 * t263;
2111 t1363 = t8 * t599 * t74;
2112 t1366 = t8 * t609 * t53;
2113 t1368 = -t8 * t1346 * t86 / 0.2e1 + t1056 / 0.9e1 + t1075 / 0.27e2 + t1079 / 0.9e1 + t1096 / 0.3e1 + t1099 / 0.3e1 + t1101 / 0.27e2 - t1103 / 0.54e2 - t1118 / 0.6e1 - 0.5e1 / 0.9e1 * t1126 - 0.5e1 / 0.81e2 * t1129 - 0.5e1 / 0.9e1 * t1132 + t1145 / 0.27e2 - 0.5e1 / 0.2e1 * t1363 + 0.5e1 / 0.2e1 * t1366 + t1150 + t1154;
2114 t1376 = t553 / 0.3e1;
2115 t1378 = 0.4e1 / 0.27e2 * t562;
2116 t1380 = t572 / 0.3e1;
2117 t1381 = t10 * (t1244 + t1285 + t1345 + t1368) - t528 / 0.2e1 + t533 - t537 / 0.2e1 - 0.4e1 / 0.9e1 * t543 + t549 / 0.3e1 - t1376 - 0.2e1 / 0.27e2 * t559 - t1378 + t568 / 0.3e1 - t1380;
2118 t1383 = 0.4e1 / 0.27e2 * t583;
2119 t1387 = -t577 - 0.2e1 / 0.9e1 * t581 - t1383 - 0.2e1 / 0.9e1 * t587 - t446 - t452 - 0.8e1 / 0.27e2 * t459 - t465 - t469 - t426 - 0.3e1 / 0.2e1 * t591;
2120 t1393 = 0.3e1 / 0.2e1 * t595 - t1176 - t1179 - 0.7e1 / 0.2e1 * t601 + t728 + t765 + 0.7e1 / 0.2e1 * t611 + t766 + t429 + t438 - t441 / 0.3e1;
2121 t1398 = t455 + t471 / 0.9e1 - t473 / 0.9e1 - t486 - t488 - t498 - t504 - t510 + t418 - t421 / 0.3e1 - 0.4e1 / 0.27e2 * t626;
2122 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2123 v4rho4[1] = t1381 + t1387 + t1393 + t1398;
2124
2125 t1404 = t561 * t641;
2126 t1406 = t864 * t317;
2127 t1413 = t546 * t551 * t666;
2128 t1416 = t287 * t1232 * t676;
2129 t1423 = t272 * t551 * t636;
2130 t1426 = t249 * t1257 * t644;
2131 t1433 = t249 * t570 * t658;
2132 t1440 = t565 * t570 * t650;
2133 t1454 = t1044 * t655;
2134 t1468 = -0.2e1 / 0.81e2 * t556 * t155 * t333 * t228 - 0.4e1 / 0.81e2 * t1404 + 0.8e1 / 0.243e3 * t1406 - 0.5e1 / 0.27e2 * t1224 * t273 * t310 * t228 + 0.2e1 / 0.27e2 * t1413 - 0.28e2 / 0.81e2 * t1416 + t546 * t273 * t320 * t228 / 0.9e1 - t1423 / 0.9e1 + 0.7e1 / 0.9e1 * t1426 + t565 * t250 * t305 * t180 / 0.9e1 - t1433 / 0.9e1 - 0.5e1 / 0.27e2 * t1204 * t250 * t299 * t180 + 0.2e1 / 0.27e2 * t1440 - t272 * t273 * t688 * t113 / 0.3e1 - t272 * t273 * t320 * t237 / 0.6e1 + t565 * t250 * t299 * t199 / 0.9e1 - 0.4e1 / 0.243e3 * t1454 - 0.2e1 / 0.81e2 * t448 * t155 * t220 * t299 - t249 * t250 * t683 * t85 / 0.3e1 - t249 * t250 * t305 * t199 / 0.6e1;
2135 t1469 = t458 * t663;
2136 t1488 = t272 * t1188 * t670;
2137 t1504 = -0.4e1 / 0.81e2 * t1469 - 0.2e1 / 0.27e2 * t288 * t155 * t333 * t237 - 0.2e1 / 0.27e2 * t191 * t155 * t220 * t305 - 0.4e1 / 0.27e2 * t191 * t155 * t107 * t683 + t546 * t273 * t310 * t237 / 0.9e1 - 0.7e1 / 0.9e1 * t1488 - t1182 / 0.9e1 + 0.2e1 / 0.9e1 * t1190 - 0.4e1 / 0.81e2 * t1192 - 0.4e1 / 0.243e3 * t1201 + 0.2e1 / 0.27e2 * t1210 - t1217 / 0.9e1 + 0.2e1 / 0.27e2 * t1230 - 0.8e1 / 0.81e2 * t1234 + 0.4e1 / 0.27e2 * t1242 + 0.4e1 / 0.27e2 * t1247 + 0.16e2 / 0.243e3 * t1253 - 0.4e1 / 0.81e2 * t1255 - 0.2e1 / 0.9e1 * t1259 + 0.16e2 / 0.243e3 * t1277;
2138 t1515 = t555 * t539;
2139 t1536 = t1200 * t714;
2140 t1544 = 0.1e1 / t90 / t1086;
2141 t1546 = t9 * t9;
2142 t1552 = t8 * t185 * t684;
2143 t1560 = t8 * t233 * t689;
2144 t1563 = t8 * t427 * t321;
2145 t1565 = t603 * t335;
2146 t1571 = t8 * t1077 * t311;
2147 t1576 = 0.8e1 / 0.243e3 * t1009 - 0.4e1 / 0.81e2 * t556 * t155 * t316 * t237 - 0.4e1 / 0.9e1 * t546 * t530 * t274 * t9 + 0.8e1 / 0.81e2 * t1515 * t155 * t675 * t113 + 0.8e1 / 0.243e3 * t1036 * t155 * t654 * t85 - 0.4e1 / 0.81e2 * t448 * t155 * t662 * t85 - 0.8e1 / 0.81e2 * t1287 * t155 * t258 * t9 + 0.8e1 / 0.243e3 * t1195 * t155 * t316 * t228 - 0.8e1 / 0.243e3 * t1536 + 0.4e1 / 0.9e1 * t565 * t574 * t251 * t9 - 0.8e1 / 0.243e3 * t1292 - t8 * t178 * t1544 * t53 * t1546 / 0.3e1 - t1552 / 0.3e1 - t8 * t226 * t1544 * t74 * t1546 / 0.3e1 - t1560 / 0.3e1 + t1563 / 0.9e1 + t1565 / 0.81e2 - t8 * t1346 * t126 / 0.3e1 + t1571 / 0.27e2 + t8 * t1337 * t131 / 0.3e1;
2148 t1581 = t8 * t1054 * t300;
2149 t1584 = t8 * t453 * t306;
2150 t1601 = 0.4e1 / 0.27e2 * t287 * t263 * t146 + t1581 / 0.27e2 + t1584 / 0.9e1 + 0.41e2 / 0.9e1 * t1322 - 0.10e2 / 0.27e2 * t1325 - 0.10e2 / 0.27e2 * t1331 - 0.41e2 / 0.9e1 * t1335 - 0.10e2 / 0.243e3 * t1342 + t1056 / 0.27e2 + t1079 / 0.27e2 + t1096 / 0.9e1 + t1099 / 0.9e1 + t1101 / 0.81e2 - 0.10e2 / 0.27e2 * t1126 - 0.10e2 / 0.243e3 * t1129 - 0.10e2 / 0.27e2 * t1132 + t1363 / 0.3e1 - t1366 / 0.3e1 + t1150 + t1154;
2151 t1617 = t10 * (t1468 + t1504 + t1576 + t1601) - t638 / 0.3e1 - 0.4e1 / 0.27e2 * t642 - 0.2e1 / 0.3e1 * t646 - 0.8e1 / 0.81e2 * t648 + 0.2e1 / 0.9e1 * t652 - 0.4e1 / 0.81e2 * t656 - t660 / 0.3e1 - 0.4e1 / 0.27e2 * t664 + 0.2e1 / 0.9e1 * t668 + 0.2e1 / 0.3e1 * t672 + 0.8e1 / 0.27e2 * t677 - t528 / 0.3e1;
2152 t1631 = 0.2e1 / 0.3e1 * t533 - t537 / 0.3e1 - 0.8e1 / 0.27e2 * t543 + 0.2e1 / 0.9e1 * t549 - 0.4e1 / 0.9e1 * t553 - 0.4e1 / 0.81e2 * t559 - 0.16e2 / 0.81e2 * t562 + 0.2e1 / 0.9e1 * t568 - 0.4e1 / 0.9e1 * t572 - 0.2e1 / 0.3e1 * t577 - 0.4e1 / 0.27e2 * t581 - 0.16e2 / 0.81e2 * t583 - 0.4e1 / 0.27e2 * t587;
2153 t1640 = -0.8e1 / 0.81e2 * t459 - t591 + t595 - t1176 - t1179 - t680 / 0.27e2 - t686 - t691 - t694 / 0.9e1 - t697 / 0.3e1 - t700 / 0.9e1 - t703 / 0.3e1 - 0.5e1 / 0.3e1 * t601;
2154 t1655 = 0.2e1 / 0.27e2 * t604 + 0.2e1 / 0.3e1 * t607 + 0.5e1 / 0.3e1 * t611 + 0.2e1 / 0.3e1 * t614 + 0.2e1 / 0.3e1 * t429 - t441 / 0.9e1 + 0.2e1 / 0.3e1 * t455 + 0.2e1 / 0.27e2 * t471 - t473 / 0.27e2 - t488 / 0.3e1 - t504 / 0.3e1 - t421 / 0.9e1 - 0.8e1 / 0.81e2 * t715 - 0.8e1 / 0.81e2 * t626;
2155 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2156 v4rho4[2] = t1617 + t1631 + t1640 + t1655;
2157
2158 t1712 = -0.2e1 / 0.27e2 * t191 * t155 * t107 * t781 - t272 * t273 * t786 * t113 / 0.6e1 - 0.5e1 / 0.27e2 * t1224 * t273 * t743 * t113 - t546 * t530 * t310 * t9 / 0.3e1 - t272 * t273 * t130 * t688 / 0.2e1 + t272 * t530 * t9 * t320 / 0.2e1 + 0.2e1 / 0.9e1 * t674 * t155 * t333 * t9 - 0.2e1 / 0.9e1 * t288 * t155 * t144 * t688 - 0.2e1 / 0.27e2 * t288 * t155 * t775 * t113 - 0.5e1 / 0.27e2 * t1204 * t250 * t734 * t85 + 0.4e1 / 0.243e3 * t1036 * t155 * t107 * t734 - t249 * t250 * t125 * t683 / 0.2e1 - t249 * t250 * t781 * t85 / 0.6e1 - 0.2e1 / 0.27e2 * t1404 + 0.8e1 / 0.81e2 * t1406 + t1413 / 0.9e1;
2159 t1741 = 0.4e1 / 0.27e2 * t1416 - t1423 / 0.6e1 - t1426 / 0.3e1 - t1433 / 0.6e1 + t1440 / 0.9e1 - 0.2e1 / 0.81e2 * t1454 - 0.2e1 / 0.27e2 * t1469 + t1488 / 0.3e1 + t1243 + t1248 + t1254 + t1278 + t546 * t273 * t748 * t113 / 0.3e1 + 0.4e1 / 0.81e2 * t1195 * t155 * t753 * t113 + 0.4e1 / 0.27e2 * t1515 * t155 * t316 * t9 - 0.2e1 / 0.27e2 * t556 * t155 * t757 * t113 - 0.2e1 / 0.27e2 * t556 * t155 * t761 * t113;
2160 t1754 = 0.6e1 * t491 + 0.12e2 * t1088;
2161 t1765 = t8 * t185 * t782;
2162 t1768 = t8 * t233 * t787;
2163 t1771 = t8 * t1073 * t744;
2164 t1773 = t293 * t777;
2165 t1776 = t8 * t1143 * t735;
2166 t1785 = t565 * t250 * t739 * t85 / 0.3e1 - 0.2e1 / 0.27e2 * t448 * t155 * t258 * t305 - 0.4e1 / 0.81e2 * t1536 - t8 * t80 * t53 * t1754 / 0.2e1 + t8 * t112 * t74 * t1754 / 0.2e1 - t1765 / 0.6e1 - t1768 / 0.6e1 + t1771 / 0.27e2 - t1773 / 0.54e2 + t1776 / 0.27e2 - t1552 / 0.2e1 - t1560 / 0.2e1 + t1563 / 0.3e1 + t1565 / 0.27e2 + t1571 / 0.9e1 + t1581 / 0.9e1 + t1584 / 0.3e1;
2167 t1795 = t272 * t551 * t748;
2168 t1797 = t1200 * t754;
2169 t1799 = t561 * t758;
2170 t1801 = t561 * t762;
2171 t1808 = t249 * t570 * t739;
2172 t1814 = -t1322 / 0.3e1 - 0.5e1 / 0.9e1 * t1325 - 0.5e1 / 0.9e1 * t1331 + t1335 / 0.3e1 - 0.5e1 / 0.81e2 * t1342 - 0.5e1 / 0.27e2 * t1126 - 0.5e1 / 0.243e3 * t1129 - 0.5e1 / 0.27e2 * t1132 + t1150 + t1154 - t1795 / 0.6e1 - 0.2e1 / 0.81e2 * t1797 - 0.2e1 / 0.27e2 * t1799 - 0.2e1 / 0.27e2 * t1801 + t565 * t574 * t299 * t9 / 0.3e1 - t1808 / 0.6e1 - t249 * t574 * t9 * t305 / 0.2e1;
2173 t1818 = t10 * (t1712 + t1741 + t1785 + t1814) - t785 - t790 + t738 + t747 - t779 - t742 - t751 - t756 - t760 - t764;
2174 t1828 = -t638 / 0.2e1 - 0.2e1 / 0.9e1 * t642 - t646 - 0.8e1 / 0.27e2 * t648 + t652 / 0.3e1 - 0.2e1 / 0.27e2 * t656 - t660 / 0.2e1 - 0.2e1 / 0.9e1 * t664 + t668 / 0.3e1 + t672 + 0.4e1 / 0.9e1 * t677;
2175 t1834 = -t1376 - t1378 - t1380 - t1383 - t1176 - t1179 - t680 / 0.9e1 - 0.3e1 / 0.2e1 * t686 - 0.3e1 / 0.2e1 * t691 - t694 / 0.3e1 - t697;
2176 t1838 = -t700 / 0.3e1 - t703 + t601 + t604 / 0.9e1 + t607 - t611 + t614 + t430 + t456 + t472 - 0.4e1 / 0.27e2 * t715;
2177 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2178 v4rho4[3] = t1818 + t1828 + t1834 + t1838;
2179
2180 t1845 = t320 * t320;
2181 t1851 = 0.12e2 * t491 + 0.12e2 * t1088;
2182 t1863 = 0.1e1 / t63 / t768 / rho[1];
2183 t1873 = t299 * t299;
2184 t1878 = t305 * t305;
2185 t1883 = t310 * t310;
2186 t1898 = 0.16e2 / 0.81e2 * t1406 - 0.8e1 / 0.27e2 * t556 * t155 * t316 * t320 - t8 * t227 * t74 * t1845 / 0.2e1 - t8 * t80 * t53 * t1851 / 0.2e1 + t8 * t112 * t74 * t1851 / 0.2e1 - t136 * t99 * t37 * (0.6160e4 / 0.81e2 * t61 * t1863 - 0.12320e5 / 0.81e2 * tau[1] * t1863) / 0.18e2 - 0.5e1 / 0.27e2 * t8 * t1137 * t53 * t1873 - t8 * t179 * t53 * t1878 / 0.2e1 - 0.5e1 / 0.27e2 * t8 * t1067 * t74 * t1883 - 0.2e1 / 0.3e1 * t1765 - 0.2e1 / 0.3e1 * t1768 + 0.4e1 / 0.27e2 * t1771 - 0.2e1 / 0.27e2 * t1773 + 0.4e1 / 0.27e2 * t1776 + 0.2e1 / 0.3e1 * t1563 + 0.2e1 / 0.27e2 * t1565 + 0.2e1 / 0.9e1 * t1571 + 0.2e1 / 0.9e1 * t1581 + 0.2e1 / 0.3e1 * t1584;
2187 t1943 = -0.20e2 / 0.27e2 * t1325 - 0.20e2 / 0.27e2 * t1331 - 0.20e2 / 0.243e3 * t1342 + t1150 + t1154 - 0.2e1 / 0.3e1 * t272 * t273 * t130 * t786 - 0.4e1 / 0.27e2 * t556 * t155 * t333 * t310 - 0.8e1 / 0.27e2 * t288 * t155 * t775 * t130 - 0.4e1 / 0.9e1 * t288 * t155 * t333 * t320 - 0.8e1 / 0.27e2 * t288 * t155 * t144 * t786 + 0.2e1 / 0.3e1 * t565 * t250 * t299 * t305 - 0.2e1 / 0.3e1 * t249 * t250 * t125 * t781 + 0.2e1 / 0.3e1 * t546 * t273 * t310 * t320 + 0.16e2 / 0.243e3 * t1195 * t155 * t144 * t743 - 0.2e1 / 0.3e1 * t1795 - 0.8e1 / 0.81e2 * t1797 - 0.8e1 / 0.27e2 * t1799 - 0.8e1 / 0.27e2 * t1801 - 0.2e1 / 0.3e1 * t1808;
2188 t1956 = t10 * (t1898 + t1943) + 0.4e1 / 0.27e2 * t604 - 0.2e1 / 0.3e1 * t694 - 0.2e1 * t697 - 0.2e1 / 0.3e1 * t700 - 0.2e1 * t703 - 0.2e1 / 0.9e1 * t680 + 0.4e1 / 0.9e1 * t737 - 0.2e1 * t741 + 0.4e1 / 0.9e1 * t746 - 0.2e1 * t750;
2189 t1966 = -0.8e1 / 0.27e2 * t755 - 0.8e1 / 0.9e1 * t759 - 0.8e1 / 0.9e1 * t763 + 0.4e1 / 0.3e1 * t607 + 0.4e1 / 0.3e1 * t614 - 0.16e2 / 0.27e2 * t648 - t1176 - t1179 - 0.2e1 / 0.9e1 * t778 - 0.2e1 * t784 - 0.2e1 * t789;
2190 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2191 v4rho4[4] = t1956 + t1966;
2192
2193 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2194 v4rho3sigma[0] = 0.0e0;
2195
2196 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2197 v4rho3sigma[1] = 0.0e0;
2198
2199 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2200 v4rho3sigma[2] = 0.0e0;
2201
2202 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2203 v4rho3sigma[3] = 0.0e0;
2204
2205 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2206 v4rho3sigma[4] = 0.0e0;
2207
2208 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2209 v4rho3sigma[5] = 0.0e0;
2210
2211 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2212 v4rho3sigma[6] = 0.0e0;
2213
2214 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2215 v4rho3sigma[7] = 0.0e0;
2216
2217 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2218 v4rho3sigma[8] = 0.0e0;
2219
2220 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2221 v4rho3sigma[9] = 0.0e0;
2222
2223 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2224 v4rho3sigma[10] = 0.0e0;
2225
2226 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2227 v4rho3sigma[11] = 0.0e0;
2228
2229 t1967 = t340 * t816;
2230 t1981 = t191 * t812;
2231 t1983 = t205 * params->BB;
2232 t1986 = 0.4e1 / 0.243e3 * t97 * t1983 * t157;
2233 t1987 = t796 * t350;
2234 t1989 = t432 * params->BB;
2235 t1990 = t152 * t1989;
2236 t1999 = t376 * t199;
2237 t2007 = t448 * t804;
2238 t2009 = t458 * t346;
2239 t2011 = t191 * t808;
2240 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2241 v4rho3lapl[0] = -0.80e2 / 0.81e2 * t1967 + 0.10e2 / 0.27e2 * t344 * t155 * t349 * t199 - 0.2e1 / 0.27e2 * t344 * t155 * t156 * t506 + 0.220e3 / 0.243e3 * t154 * t155 * t39 * t477 - 0.8e1 / 0.27e2 * t1981 + t1986 + 0.10e2 / 0.81e2 * t1987 + 0.4e1 / 0.243e3 * t1990 * t155 * t156 * t434 + 0.10e2 / 0.81e2 * t802 * t155 * t349 * t180 - 0.2e1 / 0.27e2 * t839 * t840 * t1999 - 0.80e2 / 0.81e2 * t344 * t155 * t815 * t85 - 0.8e1 / 0.81e2 * t2007 - 0.8e1 / 0.81e2 * t2009 + 0.80e2 / 0.81e2 * t2011;
2242
2243 t2013 = t556 * t827;
2244 t2015 = t561 * t358;
2245 t2017 = t288 * t831;
2246 t2021 = 0.4e1 / 0.243e3 * t135 * t1983 * t163;
2247 t2022 = t412 * params->BB;
2248 t2023 = t152 * t2022;
2249 t2028 = t385 * t237;
2250 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2251 v4rho3lapl[1] = -0.8e1 / 0.81e2 * t2013 - 0.8e1 / 0.81e2 * t2015 - 0.8e1 / 0.27e2 * t2017 + t2021 + 0.4e1 / 0.243e3 * t2023 * t155 * t162 * t414 - 0.2e1 / 0.27e2 * t859 * t840 * t2028 - 0.2e1 / 0.27e2 * t356 * t155 * t162 * t494;
2252
2253 t2042 = t152 * t1989 * t34;
2254 t2043 = t389 * t180;
2255 t2047 = t943 * t85;
2256 t2051 = t849 * t838;
2257 t2052 = t947 * t85;
2258 t2057 = 0.80e2 / 0.243e3 * t191 * t846;
2259 t2058 = t540 * t852;
2260 t2060 = t389 * t199;
2261 t2072 = t458 * t362;
2262 t2074 = t17 * t96;
2263 t2077 = 0.16e2 / 0.243e3 * t2074 * t838 * t842;
2264 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2265 v4rho3lapl[2] = -0.2e1 / 0.27e2 * t850 * t157 - 0.80e2 / 0.243e3 * t1967 - 0.8e1 / 0.81e2 * t1981 + t1986 + 0.20e2 / 0.243e3 * t1987 + 0.4e1 / 0.243e3 * t2042 * t840 * t2043 + 0.20e2 / 0.243e3 * t839 * t840 * t2047 - 0.4e1 / 0.81e2 * t2051 * t840 * t2052 + t2057 + 0.2e1 / 0.81e2 * t2058 - 0.2e1 / 0.81e2 * t839 * t840 * t2060 - 0.80e2 / 0.243e3 * t344 * t155 * t815 * t125 + 0.20e2 / 0.81e2 * t850 * t155 * t349 * t9 - 0.8e1 / 0.243e3 * t2072 - t2077 - 0.8e1 / 0.243e3 * t2007 - 0.16e2 / 0.243e3 * t2009 + 0.80e2 / 0.243e3 * t2011;
2266
2267 t2084 = t561 * t367;
2268 t2088 = 0.16e2 / 0.243e3 * t2074 * t858 * t861;
2269 t2089 = t674 * t866;
2270 t2092 = t152 * t2022 * t34;
2271 t2093 = t394 * t228;
2272 t2097 = t849 * t858;
2273 t2098 = t957 * t113;
2274 t2102 = t394 * t237;
2275 t2108 = t821 * t371;
2276 t2111 = 0.80e2 / 0.243e3 * t288 * t872;
2277 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2278 v4rho3lapl[3] = -0.8e1 / 0.243e3 * t2013 - 0.16e2 / 0.243e3 * t2015 - 0.8e1 / 0.81e2 * t2017 + t2021 - 0.8e1 / 0.243e3 * t2084 - t2088 - 0.2e1 / 0.81e2 * t2089 + 0.4e1 / 0.243e3 * t2092 * t840 * t2093 + 0.4e1 / 0.81e2 * t2097 * t840 * t2098 - 0.2e1 / 0.81e2 * t859 * t840 * t2102 + 0.2e1 / 0.27e2 * t864 * t163 + 0.10e2 / 0.243e3 * t2108 + t2111 + 0.10e2 / 0.243e3 * t825 * t155 * t370 * t228 + 0.10e2 / 0.81e2 * t356 * t155 * t370 * t237;
2279
2280 t2124 = t448 * t877;
2281 t2126 = t968 * t85;
2282 t2134 = t389 * t9;
2283 t2138 = t191 * t881;
2284 t2140 = t972 * t85;
2285 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2286 v4rho3lapl[4] = -0.16e2 / 0.243e3 * t2072 - t2077 + t2057 - 0.16e2 / 0.81e2 * t2058 - 0.8e1 / 0.243e3 * t2009 + t1986 + 0.10e2 / 0.243e3 * t1987 - 0.8e1 / 0.243e3 * t2124 + 0.4e1 / 0.243e3 * t2042 * t840 * t2126 + 0.10e2 / 0.243e3 * t802 * t155 * t349 * t299 - 0.4e1 / 0.81e2 * t2051 * t840 * t2134 - 0.8e1 / 0.81e2 * t2138 - 0.2e1 / 0.81e2 * t839 * t840 * t2140 + 0.10e2 / 0.81e2 * t344 * t155 * t349 * t305 - 0.2e1 / 0.27e2 * t344 * t155 * t156 * t683;
2287
2288 t2156 = t556 * t887;
2289 t2158 = t978 * t113;
2290 t2162 = t394 * t9;
2291 t2166 = t288 * t891;
2292 t2168 = t982 * t113;
2293 t2176 = t288 * t895;
2294 t2178 = t986 * t113;
2295 t2186 = t353 * t899;
2296 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2297 v4rho3lapl[5] = -0.16e2 / 0.243e3 * t2084 - t2088 + 0.16e2 / 0.81e2 * t2089 - 0.8e1 / 0.243e3 * t2015 + t2021 + 0.20e2 / 0.243e3 * t2108 + t2111 - 0.8e1 / 0.243e3 * t2156 + 0.4e1 / 0.243e3 * t2092 * t840 * t2158 + 0.4e1 / 0.81e2 * t2097 * t840 * t2162 + 0.80e2 / 0.243e3 * t2166 + 0.20e2 / 0.243e3 * t859 * t840 * t2168 - 0.20e2 / 0.81e2 * t864 * t155 * t370 * t9 - 0.8e1 / 0.81e2 * t2176 - 0.2e1 / 0.81e2 * t859 * t840 * t2178 - 0.2e1 / 0.27e2 * t356 * t155 * t162 * t688 - 0.80e2 / 0.243e3 * t2186 - 0.80e2 / 0.243e3 * t356 * t155 * t898 * t113;
2298
2299 t2199 = t389 * t305;
2300 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2301 v4rho3lapl[6] = -0.8e1 / 0.81e2 * t2072 - 0.8e1 / 0.81e2 * t2124 - 0.8e1 / 0.27e2 * t2138 + t1986 + 0.4e1 / 0.243e3 * t1990 * t155 * t156 * t734 - 0.2e1 / 0.27e2 * t839 * t840 * t2199 - 0.2e1 / 0.27e2 * t344 * t155 * t156 * t781;
2302
2303 t2229 = t394 * t320;
2304 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2305 v4rho3lapl[7] = 0.10e2 / 0.27e2 * t356 * t155 * t370 * t320 - 0.2e1 / 0.27e2 * t356 * t155 * t162 * t786 + 0.220e3 / 0.243e3 * t161 * t155 * t39 * t770 - 0.80e2 / 0.81e2 * t2186 + 0.10e2 / 0.81e2 * t2108 + t2021 + 0.4e1 / 0.243e3 * t2023 * t155 * t162 * t743 + 0.10e2 / 0.81e2 * t825 * t155 * t370 * t310 - 0.2e1 / 0.27e2 * t859 * t840 * t2229 - 0.80e2 / 0.81e2 * t356 * t155 * t898 * t130 + 0.80e2 / 0.81e2 * t2166 - 0.8e1 / 0.27e2 * t2176 - 0.8e1 / 0.81e2 * t2156 - 0.8e1 / 0.81e2 * t2084;
2306
2307 t2253 = t191 * t917;
2308 t2256 = 0.8e1 / 0.243e3 * t470 * t168;
2309 t2257 = t210 * t381;
2310 t2259 = t98 * t921;
2311 t2261 = t448 * t909;
2312 t2263 = t458 * t377;
2313 t2265 = t191 * t913;
2314 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2315 v4rho3tau[0] = -0.20e2 / 0.27e2 * t344 * t155 * t102 * t199 + 0.4e1 / 0.27e2 * t344 * t155 * t44 * t506 - 0.440e3 / 0.243e3 * t166 * t99 * t37 * t477 + 0.16e2 / 0.27e2 * t2253 - t2256 - 0.20e2 / 0.81e2 * t2257 + 0.160e3 / 0.81e2 * t2259 + 0.16e2 / 0.81e2 * t2261 + 0.16e2 / 0.81e2 * t2263 - 0.160e3 / 0.81e2 * t2265 - 0.8e1 / 0.243e3 * t1990 * t155 * t44 * t434 - 0.20e2 / 0.81e2 * t802 * t155 * t102 * t180 + 0.4e1 / 0.27e2 * t802 * t155 * t1999 + 0.160e3 / 0.81e2 * t344 * t155 * t215 * t85;
2316
2317 t2282 = t556 * t929;
2318 t2284 = t561 * t386;
2319 t2286 = t288 * t933;
2320 t2289 = 0.8e1 / 0.243e3 * t603 * t172;
2321 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2322 v4rho3tau[1] = 0.16e2 / 0.81e2 * t2282 + 0.16e2 / 0.81e2 * t2284 + 0.16e2 / 0.27e2 * t2286 - t2289 - 0.8e1 / 0.243e3 * t2023 * t155 * t65 * t414 + 0.4e1 / 0.27e2 * t825 * t155 * t2028 + 0.4e1 / 0.27e2 * t356 * t155 * t65 * t494;
2323
2324 t2307 = t458 * t390;
2325 t2310 = 0.32e2 / 0.243e3 * t448 * t940;
2326 t2312 = 0.160e3 / 0.243e3 * t191 * t944;
2327 t2313 = t540 * t948;
2328 t2335 = t849 * t801;
2329 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2330 v4rho3tau[2] = 0.4e1 / 0.27e2 * t849 * t15 * t168 + 0.16e2 / 0.81e2 * t2253 - t2256 - 0.40e2 / 0.243e3 * t2257 + 0.160e3 / 0.243e3 * t2259 + 0.16e2 / 0.243e3 * t2307 + t2310 - t2312 - 0.4e1 / 0.81e2 * t2313 + 0.4e1 / 0.81e2 * t802 * t155 * t2060 + 0.160e3 / 0.243e3 * t344 * t155 * t215 * t125 - 0.40e2 / 0.81e2 * t850 * t155 * t102 * t9 + 0.16e2 / 0.243e3 * t2261 + 0.32e2 / 0.243e3 * t2263 - 0.160e3 / 0.243e3 * t2265 - 0.8e1 / 0.243e3 * t1990 * t155 * t2043 - 0.40e2 / 0.243e3 * t802 * t155 * t2047 + 0.8e1 / 0.81e2 * t2335 * t155 * t2052;
2331
2332 t2342 = t561 * t395;
2333 t2345 = 0.32e2 / 0.243e3 * t556 * t954;
2334 t2346 = t674 * t958;
2335 t2351 = t849 * t824;
2336 t2361 = t293 * t399;
2337 t2364 = 0.160e3 / 0.243e3 * t288 * t964;
2338 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2339 v4rho3tau[3] = 0.16e2 / 0.243e3 * t2282 + 0.32e2 / 0.243e3 * t2284 + 0.16e2 / 0.81e2 * t2286 - t2289 + 0.16e2 / 0.243e3 * t2342 + t2345 + 0.4e1 / 0.81e2 * t2346 - 0.8e1 / 0.243e3 * t2023 * t155 * t2093 - 0.8e1 / 0.81e2 * t2351 * t155 * t2098 + 0.4e1 / 0.81e2 * t825 * t155 * t2102 - 0.4e1 / 0.27e2 * t849 * t58 * t172 - 0.20e2 / 0.243e3 * t2361 - t2364 - 0.20e2 / 0.243e3 * t825 * t155 * t139 * t228 - 0.20e2 / 0.81e2 * t356 * t155 * t139 * t237;
2340
2341 t2377 = t448 * t969;
2342 t2389 = t191 * t973;
2343 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2344 v4rho3tau[4] = 0.32e2 / 0.243e3 * t2307 + t2310 - t2312 + 0.32e2 / 0.81e2 * t2313 + 0.16e2 / 0.243e3 * t2263 - t2256 - 0.20e2 / 0.243e3 * t2257 + 0.16e2 / 0.243e3 * t2377 - 0.8e1 / 0.243e3 * t1990 * t155 * t2126 - 0.20e2 / 0.243e3 * t802 * t155 * t102 * t299 + 0.8e1 / 0.81e2 * t2335 * t155 * t2134 + 0.16e2 / 0.81e2 * t2389 + 0.4e1 / 0.81e2 * t802 * t155 * t2140 - 0.20e2 / 0.81e2 * t344 * t155 * t102 * t305 + 0.4e1 / 0.27e2 * t344 * t155 * t44 * t683;
2345
2346 t2406 = t556 * t979;
2347 t2414 = t288 * t983;
2348 t2423 = t288 * t987;
2349 t2432 = t136 * t991;
2350 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2351 v4rho3tau[5] = 0.32e2 / 0.243e3 * t2342 + t2345 - 0.32e2 / 0.81e2 * t2346 + 0.16e2 / 0.243e3 * t2284 - t2289 - 0.40e2 / 0.243e3 * t2361 - t2364 + 0.16e2 / 0.243e3 * t2406 - 0.8e1 / 0.243e3 * t2023 * t155 * t2158 - 0.8e1 / 0.81e2 * t2351 * t155 * t2162 - 0.160e3 / 0.243e3 * t2414 - 0.40e2 / 0.243e3 * t825 * t155 * t2168 + 0.40e2 / 0.81e2 * t864 * t155 * t139 * t9 + 0.16e2 / 0.81e2 * t2423 + 0.4e1 / 0.81e2 * t825 * t155 * t2178 + 0.4e1 / 0.27e2 * t356 * t155 * t65 * t688 + 0.160e3 / 0.243e3 * t2432 + 0.160e3 / 0.243e3 * t356 * t155 * t328 * t113;
2352
2353 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2354 v4rho3tau[6] = 0.16e2 / 0.81e2 * t2307 + 0.16e2 / 0.81e2 * t2377 + 0.16e2 / 0.27e2 * t2389 - t2256 - 0.8e1 / 0.243e3 * t1990 * t155 * t44 * t734 + 0.4e1 / 0.27e2 * t802 * t155 * t2199 + 0.4e1 / 0.27e2 * t344 * t155 * t44 * t781;
2355
2356 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2357 v4rho3tau[7] = -0.20e2 / 0.27e2 * t356 * t155 * t139 * t320 + 0.4e1 / 0.27e2 * t356 * t155 * t65 * t786 - 0.440e3 / 0.243e3 * t170 * t99 * t37 * t770 + 0.160e3 / 0.81e2 * t2432 - 0.20e2 / 0.81e2 * t2361 - t2289 - 0.8e1 / 0.243e3 * t2023 * t155 * t65 * t743 - 0.20e2 / 0.81e2 * t825 * t155 * t139 * t310 + 0.4e1 / 0.27e2 * t825 * t155 * t2229 + 0.160e3 / 0.81e2 * t356 * t155 * t328 * t130 - 0.160e3 / 0.81e2 * t2414 + 0.16e2 / 0.27e2 * t2423 + 0.16e2 / 0.81e2 * t2406 + 0.16e2 / 0.81e2 * t2342;
2358
2359 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2360 v4rho2sigma2[0] = 0.0e0;
2361
2362 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2363 v4rho2sigma2[1] = 0.0e0;
2364
2365 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2366 v4rho2sigma2[2] = 0.0e0;
2367
2368 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2369 v4rho2sigma2[3] = 0.0e0;
2370
2371 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2372 v4rho2sigma2[4] = 0.0e0;
2373
2374 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2375 v4rho2sigma2[5] = 0.0e0;
2376
2377 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2378 v4rho2sigma2[6] = 0.0e0;
2379
2380 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2381 v4rho2sigma2[7] = 0.0e0;
2382
2383 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2384 v4rho2sigma2[8] = 0.0e0;
2385
2386 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2387 v4rho2sigma2[9] = 0.0e0;
2388
2389 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2390 v4rho2sigma2[10] = 0.0e0;
2391
2392 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2393 v4rho2sigma2[11] = 0.0e0;
2394
2395 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2396 v4rho2sigma2[12] = 0.0e0;
2397
2398 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2399 v4rho2sigma2[13] = 0.0e0;
2400
2401 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2402 v4rho2sigma2[14] = 0.0e0;
2403
2404 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2405 v4rho2sigma2[15] = 0.0e0;
2406
2407 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2408 v4rho2sigma2[16] = 0.0e0;
2409
2410 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2411 v4rho2sigma2[17] = 0.0e0;
2412
2413 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2414 v4rho2sigmalapl[0] = 0.0e0;
2415
2416 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2417 v4rho2sigmalapl[1] = 0.0e0;
2418
2419 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2420 v4rho2sigmalapl[2] = 0.0e0;
2421
2422 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2423 v4rho2sigmalapl[3] = 0.0e0;
2424
2425 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2426 v4rho2sigmalapl[4] = 0.0e0;
2427
2428 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2429 v4rho2sigmalapl[5] = 0.0e0;
2430
2431 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2432 v4rho2sigmalapl[6] = 0.0e0;
2433
2434 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2435 v4rho2sigmalapl[7] = 0.0e0;
2436
2437 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2438 v4rho2sigmalapl[8] = 0.0e0;
2439
2440 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2441 v4rho2sigmalapl[9] = 0.0e0;
2442
2443 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2444 v4rho2sigmalapl[10] = 0.0e0;
2445
2446 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2447 v4rho2sigmalapl[11] = 0.0e0;
2448
2449 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2450 v4rho2sigmalapl[12] = 0.0e0;
2451
2452 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2453 v4rho2sigmalapl[13] = 0.0e0;
2454
2455 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2456 v4rho2sigmalapl[14] = 0.0e0;
2457
2458 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2459 v4rho2sigmalapl[15] = 0.0e0;
2460
2461 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2462 v4rho2sigmalapl[16] = 0.0e0;
2463
2464 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2465 v4rho2sigmalapl[17] = 0.0e0;
2466
2467 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2468 v4rho2sigmatau[0] = 0.0e0;
2469
2470 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2471 v4rho2sigmatau[1] = 0.0e0;
2472
2473 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2474 v4rho2sigmatau[2] = 0.0e0;
2475
2476 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2477 v4rho2sigmatau[3] = 0.0e0;
2478
2479 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2480 v4rho2sigmatau[4] = 0.0e0;
2481
2482 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2483 v4rho2sigmatau[5] = 0.0e0;
2484
2485 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2486 v4rho2sigmatau[6] = 0.0e0;
2487
2488 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2489 v4rho2sigmatau[7] = 0.0e0;
2490
2491 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2492 v4rho2sigmatau[8] = 0.0e0;
2493
2494 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2495 v4rho2sigmatau[9] = 0.0e0;
2496
2497 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2498 v4rho2sigmatau[10] = 0.0e0;
2499
2500 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2501 v4rho2sigmatau[11] = 0.0e0;
2502
2503 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2504 v4rho2sigmatau[12] = 0.0e0;
2505
2506 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2507 v4rho2sigmatau[13] = 0.0e0;
2508
2509 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2510 v4rho2sigmatau[14] = 0.0e0;
2511
2512 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2513 v4rho2sigmatau[15] = 0.0e0;
2514
2515 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2516 v4rho2sigmatau[16] = 0.0e0;
2517
2518 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2519 v4rho2sigmatau[17] = 0.0e0;
2520
2521 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2522 v4rho2lapl2[0] = 0.0e0;
2523
2524 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2525 v4rho2lapl2[1] = 0.0e0;
2526
2527 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2528 v4rho2lapl2[2] = 0.0e0;
2529
2530 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2531 v4rho2lapl2[3] = 0.0e0;
2532
2533 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2534 v4rho2lapl2[4] = 0.0e0;
2535
2536 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2537 v4rho2lapl2[5] = 0.0e0;
2538
2539 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2540 v4rho2lapl2[6] = 0.0e0;
2541
2542 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2543 v4rho2lapl2[7] = 0.0e0;
2544
2545 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2546 v4rho2lapl2[8] = 0.0e0;
2547
2548 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2549 v4rho2lapltau[0] = 0.0e0;
2550
2551 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2552 v4rho2lapltau[1] = 0.0e0;
2553
2554 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2555 v4rho2lapltau[2] = 0.0e0;
2556
2557 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2558 v4rho2lapltau[3] = 0.0e0;
2559
2560 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2561 v4rho2lapltau[4] = 0.0e0;
2562
2563 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2564 v4rho2lapltau[5] = 0.0e0;
2565
2566 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2567 v4rho2lapltau[6] = 0.0e0;
2568
2569 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2570 v4rho2lapltau[7] = 0.0e0;
2571
2572 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2573 v4rho2lapltau[8] = 0.0e0;
2574
2575 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2576 v4rho2lapltau[9] = 0.0e0;
2577
2578 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2579 v4rho2lapltau[10] = 0.0e0;
2580
2581 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2582 v4rho2lapltau[11] = 0.0e0;
2583
2584 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2585 v4rho2tau2[0] = 0.0e0;
2586
2587 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2588 v4rho2tau2[1] = 0.0e0;
2589
2590 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2591 v4rho2tau2[2] = 0.0e0;
2592
2593 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2594 v4rho2tau2[3] = 0.0e0;
2595
2596 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2597 v4rho2tau2[4] = 0.0e0;
2598
2599 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2600 v4rho2tau2[5] = 0.0e0;
2601
2602 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2603 v4rho2tau2[6] = 0.0e0;
2604
2605 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2606 v4rho2tau2[7] = 0.0e0;
2607
2608 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2609 v4rho2tau2[8] = 0.0e0;
2610
2611 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2612 v4rhosigma3[0] = 0.0e0;
2613
2614 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2615 v4rhosigma3[1] = 0.0e0;
2616
2617 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2618 v4rhosigma3[2] = 0.0e0;
2619
2620 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2621 v4rhosigma3[3] = 0.0e0;
2622
2623 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2624 v4rhosigma3[4] = 0.0e0;
2625
2626 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2627 v4rhosigma3[5] = 0.0e0;
2628
2629 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2630 v4rhosigma3[6] = 0.0e0;
2631
2632 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2633 v4rhosigma3[7] = 0.0e0;
2634
2635 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2636 v4rhosigma3[8] = 0.0e0;
2637
2638 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2639 v4rhosigma3[9] = 0.0e0;
2640
2641 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2642 v4rhosigma3[10] = 0.0e0;
2643
2644 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2645 v4rhosigma3[11] = 0.0e0;
2646
2647 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2648 v4rhosigma3[12] = 0.0e0;
2649
2650 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2651 v4rhosigma3[13] = 0.0e0;
2652
2653 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2654 v4rhosigma3[14] = 0.0e0;
2655
2656 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2657 v4rhosigma3[15] = 0.0e0;
2658
2659 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2660 v4rhosigma3[16] = 0.0e0;
2661
2662 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2663 v4rhosigma3[17] = 0.0e0;
2664
2665 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2666 v4rhosigma3[18] = 0.0e0;
2667
2668 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2669 v4rhosigma3[19] = 0.0e0;
2670
2671 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2672 v4rhosigma2lapl[0] = 0.0e0;
2673
2674 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2675 v4rhosigma2lapl[1] = 0.0e0;
2676
2677 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2678 v4rhosigma2lapl[2] = 0.0e0;
2679
2680 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2681 v4rhosigma2lapl[3] = 0.0e0;
2682
2683 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2684 v4rhosigma2lapl[4] = 0.0e0;
2685
2686 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2687 v4rhosigma2lapl[5] = 0.0e0;
2688
2689 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2690 v4rhosigma2lapl[6] = 0.0e0;
2691
2692 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2693 v4rhosigma2lapl[7] = 0.0e0;
2694
2695 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2696 v4rhosigma2lapl[8] = 0.0e0;
2697
2698 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2699 v4rhosigma2lapl[9] = 0.0e0;
2700
2701 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2702 v4rhosigma2lapl[10] = 0.0e0;
2703
2704 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2705 v4rhosigma2lapl[11] = 0.0e0;
2706
2707 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2708 v4rhosigma2lapl[12] = 0.0e0;
2709
2710 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2711 v4rhosigma2lapl[13] = 0.0e0;
2712
2713 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2714 v4rhosigma2lapl[14] = 0.0e0;
2715
2716 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2717 v4rhosigma2lapl[15] = 0.0e0;
2718
2719 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2720 v4rhosigma2lapl[16] = 0.0e0;
2721
2722 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2723 v4rhosigma2lapl[17] = 0.0e0;
2724
2725 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2726 v4rhosigma2lapl[18] = 0.0e0;
2727
2728 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2729 v4rhosigma2lapl[19] = 0.0e0;
2730
2731 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2732 v4rhosigma2lapl[20] = 0.0e0;
2733
2734 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2735 v4rhosigma2lapl[21] = 0.0e0;
2736
2737 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2738 v4rhosigma2lapl[22] = 0.0e0;
2739
2740 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2741 v4rhosigma2lapl[23] = 0.0e0;
2742
2743 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2744 v4rhosigma2tau[0] = 0.0e0;
2745
2746 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2747 v4rhosigma2tau[1] = 0.0e0;
2748
2749 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2750 v4rhosigma2tau[2] = 0.0e0;
2751
2752 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2753 v4rhosigma2tau[3] = 0.0e0;
2754
2755 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2756 v4rhosigma2tau[4] = 0.0e0;
2757
2758 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2759 v4rhosigma2tau[5] = 0.0e0;
2760
2761 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2762 v4rhosigma2tau[6] = 0.0e0;
2763
2764 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2765 v4rhosigma2tau[7] = 0.0e0;
2766
2767 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2768 v4rhosigma2tau[8] = 0.0e0;
2769
2770 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2771 v4rhosigma2tau[9] = 0.0e0;
2772
2773 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2774 v4rhosigma2tau[10] = 0.0e0;
2775
2776 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2777 v4rhosigma2tau[11] = 0.0e0;
2778
2779 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2780 v4rhosigma2tau[12] = 0.0e0;
2781
2782 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2783 v4rhosigma2tau[13] = 0.0e0;
2784
2785 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2786 v4rhosigma2tau[14] = 0.0e0;
2787
2788 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2789 v4rhosigma2tau[15] = 0.0e0;
2790
2791 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2792 v4rhosigma2tau[16] = 0.0e0;
2793
2794 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2795 v4rhosigma2tau[17] = 0.0e0;
2796
2797 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2798 v4rhosigma2tau[18] = 0.0e0;
2799
2800 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2801 v4rhosigma2tau[19] = 0.0e0;
2802
2803 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2804 v4rhosigma2tau[20] = 0.0e0;
2805
2806 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2807 v4rhosigma2tau[21] = 0.0e0;
2808
2809 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2810 v4rhosigma2tau[22] = 0.0e0;
2811
2812 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2813 v4rhosigma2tau[23] = 0.0e0;
2814
2815 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2816 v4rhosigmalapl2[0] = 0.0e0;
2817
2818 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2819 v4rhosigmalapl2[1] = 0.0e0;
2820
2821 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2822 v4rhosigmalapl2[2] = 0.0e0;
2823
2824 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2825 v4rhosigmalapl2[3] = 0.0e0;
2826
2827 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2828 v4rhosigmalapl2[4] = 0.0e0;
2829
2830 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2831 v4rhosigmalapl2[5] = 0.0e0;
2832
2833 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2834 v4rhosigmalapl2[6] = 0.0e0;
2835
2836 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2837 v4rhosigmalapl2[7] = 0.0e0;
2838
2839 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2840 v4rhosigmalapl2[8] = 0.0e0;
2841
2842 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2843 v4rhosigmalapl2[9] = 0.0e0;
2844
2845 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2846 v4rhosigmalapl2[10] = 0.0e0;
2847
2848 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2849 v4rhosigmalapl2[11] = 0.0e0;
2850
2851 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2852 v4rhosigmalapl2[12] = 0.0e0;
2853
2854 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2855 v4rhosigmalapl2[13] = 0.0e0;
2856
2857 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2858 v4rhosigmalapl2[14] = 0.0e0;
2859
2860 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2861 v4rhosigmalapl2[15] = 0.0e0;
2862
2863 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2864 v4rhosigmalapl2[16] = 0.0e0;
2865
2866 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2867 v4rhosigmalapl2[17] = 0.0e0;
2868
2869 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2870 v4rhosigmalapltau[0] = 0.0e0;
2871
2872 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2873 v4rhosigmalapltau[1] = 0.0e0;
2874
2875 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2876 v4rhosigmalapltau[2] = 0.0e0;
2877
2878 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2879 v4rhosigmalapltau[3] = 0.0e0;
2880
2881 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2882 v4rhosigmalapltau[4] = 0.0e0;
2883
2884 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2885 v4rhosigmalapltau[5] = 0.0e0;
2886
2887 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2888 v4rhosigmalapltau[6] = 0.0e0;
2889
2890 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2891 v4rhosigmalapltau[7] = 0.0e0;
2892
2893 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2894 v4rhosigmalapltau[8] = 0.0e0;
2895
2896 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2897 v4rhosigmalapltau[9] = 0.0e0;
2898
2899 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2900 v4rhosigmalapltau[10] = 0.0e0;
2901
2902 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2903 v4rhosigmalapltau[11] = 0.0e0;
2904
2905 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2906 v4rhosigmalapltau[12] = 0.0e0;
2907
2908 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2909 v4rhosigmalapltau[13] = 0.0e0;
2910
2911 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2912 v4rhosigmalapltau[14] = 0.0e0;
2913
2914 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2915 v4rhosigmalapltau[15] = 0.0e0;
2916
2917 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2918 v4rhosigmalapltau[16] = 0.0e0;
2919
2920 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2921 v4rhosigmalapltau[17] = 0.0e0;
2922
2923 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2924 v4rhosigmalapltau[18] = 0.0e0;
2925
2926 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2927 v4rhosigmalapltau[19] = 0.0e0;
2928
2929 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2930 v4rhosigmalapltau[20] = 0.0e0;
2931
2932 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2933 v4rhosigmalapltau[21] = 0.0e0;
2934
2935 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2936 v4rhosigmalapltau[22] = 0.0e0;
2937
2938 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2939 v4rhosigmalapltau[23] = 0.0e0;
2940
2941 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2942 v4rhosigmatau2[0] = 0.0e0;
2943
2944 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2945 v4rhosigmatau2[1] = 0.0e0;
2946
2947 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2948 v4rhosigmatau2[2] = 0.0e0;
2949
2950 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2951 v4rhosigmatau2[3] = 0.0e0;
2952
2953 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2954 v4rhosigmatau2[4] = 0.0e0;
2955
2956 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2957 v4rhosigmatau2[5] = 0.0e0;
2958
2959 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2960 v4rhosigmatau2[6] = 0.0e0;
2961
2962 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2963 v4rhosigmatau2[7] = 0.0e0;
2964
2965 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2966 v4rhosigmatau2[8] = 0.0e0;
2967
2968 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2969 v4rhosigmatau2[9] = 0.0e0;
2970
2971 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2972 v4rhosigmatau2[10] = 0.0e0;
2973
2974 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2975 v4rhosigmatau2[11] = 0.0e0;
2976
2977 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2978 v4rhosigmatau2[12] = 0.0e0;
2979
2980 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2981 v4rhosigmatau2[13] = 0.0e0;
2982
2983 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2984 v4rhosigmatau2[14] = 0.0e0;
2985
2986 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2987 v4rhosigmatau2[15] = 0.0e0;
2988
2989 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2990 v4rhosigmatau2[16] = 0.0e0;
2991
2992 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
2993 v4rhosigmatau2[17] = 0.0e0;
2994
2995 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2996 v4rholapl3[0] = 0.0e0;
2997
2998 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
2999 v4rholapl3[1] = 0.0e0;
3000
3001 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3002 v4rholapl3[2] = 0.0e0;
3003
3004 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3005 v4rholapl3[3] = 0.0e0;
3006
3007 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3008 v4rholapl3[4] = 0.0e0;
3009
3010 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3011 v4rholapl3[5] = 0.0e0;
3012
3013 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3014 v4rholapl3[6] = 0.0e0;
3015
3016 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3017 v4rholapl3[7] = 0.0e0;
3018
3019 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3020 v4rholapl2tau[0] = 0.0e0;
3021
3022 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3023 v4rholapl2tau[1] = 0.0e0;
3024
3025 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3026 v4rholapl2tau[2] = 0.0e0;
3027
3028 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3029 v4rholapl2tau[3] = 0.0e0;
3030
3031 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3032 v4rholapl2tau[4] = 0.0e0;
3033
3034 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3035 v4rholapl2tau[5] = 0.0e0;
3036
3037 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3038 v4rholapl2tau[6] = 0.0e0;
3039
3040 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3041 v4rholapl2tau[7] = 0.0e0;
3042
3043 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3044 v4rholapl2tau[8] = 0.0e0;
3045
3046 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3047 v4rholapl2tau[9] = 0.0e0;
3048
3049 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3050 v4rholapl2tau[10] = 0.0e0;
3051
3052 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3053 v4rholapl2tau[11] = 0.0e0;
3054
3055 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3056 v4rholapltau2[0] = 0.0e0;
3057
3058 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3059 v4rholapltau2[1] = 0.0e0;
3060
3061 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3062 v4rholapltau2[2] = 0.0e0;
3063
3064 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3065 v4rholapltau2[3] = 0.0e0;
3066
3067 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3068 v4rholapltau2[4] = 0.0e0;
3069
3070 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3071 v4rholapltau2[5] = 0.0e0;
3072
3073 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3074 v4rholapltau2[6] = 0.0e0;
3075
3076 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3077 v4rholapltau2[7] = 0.0e0;
3078
3079 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3080 v4rholapltau2[8] = 0.0e0;
3081
3082 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3083 v4rholapltau2[9] = 0.0e0;
3084
3085 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3086 v4rholapltau2[10] = 0.0e0;
3087
3088 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3089 v4rholapltau2[11] = 0.0e0;
3090
3091 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3092 v4rhotau3[0] = 0.0e0;
3093
3094 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3095 v4rhotau3[1] = 0.0e0;
3096
3097 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3098 v4rhotau3[2] = 0.0e0;
3099
3100 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3101 v4rhotau3[3] = 0.0e0;
3102
3103 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3104 v4rhotau3[4] = 0.0e0;
3105
3106 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3107 v4rhotau3[5] = 0.0e0;
3108
3109 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3110 v4rhotau3[6] = 0.0e0;
3111
3112 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3113 v4rhotau3[7] = 0.0e0;
3114
3115 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3116 v4sigma4[0] = 0.0e0;
3117
3118 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3119 v4sigma4[1] = 0.0e0;
3120
3121 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3122 v4sigma4[2] = 0.0e0;
3123
3124 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3125 v4sigma4[3] = 0.0e0;
3126
3127 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3128 v4sigma4[4] = 0.0e0;
3129
3130 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3131 v4sigma4[5] = 0.0e0;
3132
3133 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3134 v4sigma4[6] = 0.0e0;
3135
3136 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3137 v4sigma4[7] = 0.0e0;
3138
3139 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3140 v4sigma4[8] = 0.0e0;
3141
3142 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3143 v4sigma4[9] = 0.0e0;
3144
3145 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3146 v4sigma4[10] = 0.0e0;
3147
3148 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3149 v4sigma4[11] = 0.0e0;
3150
3151 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3152 v4sigma4[12] = 0.0e0;
3153
3154 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3155 v4sigma4[13] = 0.0e0;
3156
3157 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3158 v4sigma4[14] = 0.0e0;
3159
3160 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3161 v4sigma3lapl[0] = 0.0e0;
3162
3163 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3164 v4sigma3lapl[1] = 0.0e0;
3165
3166 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3167 v4sigma3lapl[2] = 0.0e0;
3168
3169 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3170 v4sigma3lapl[3] = 0.0e0;
3171
3172 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3173 v4sigma3lapl[4] = 0.0e0;
3174
3175 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3176 v4sigma3lapl[5] = 0.0e0;
3177
3178 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3179 v4sigma3lapl[6] = 0.0e0;
3180
3181 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3182 v4sigma3lapl[7] = 0.0e0;
3183
3184 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3185 v4sigma3lapl[8] = 0.0e0;
3186
3187 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3188 v4sigma3lapl[9] = 0.0e0;
3189
3190 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3191 v4sigma3lapl[10] = 0.0e0;
3192
3193 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3194 v4sigma3lapl[11] = 0.0e0;
3195
3196 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3197 v4sigma3lapl[12] = 0.0e0;
3198
3199 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3200 v4sigma3lapl[13] = 0.0e0;
3201
3202 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3203 v4sigma3lapl[14] = 0.0e0;
3204
3205 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3206 v4sigma3lapl[15] = 0.0e0;
3207
3208 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3209 v4sigma3lapl[16] = 0.0e0;
3210
3211 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3212 v4sigma3lapl[17] = 0.0e0;
3213
3214 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3215 v4sigma3lapl[18] = 0.0e0;
3216
3217 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3218 v4sigma3lapl[19] = 0.0e0;
3219
3220 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3221 v4sigma3tau[0] = 0.0e0;
3222
3223 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3224 v4sigma3tau[1] = 0.0e0;
3225
3226 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3227 v4sigma3tau[2] = 0.0e0;
3228
3229 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3230 v4sigma3tau[3] = 0.0e0;
3231
3232 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3233 v4sigma3tau[4] = 0.0e0;
3234
3235 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3236 v4sigma3tau[5] = 0.0e0;
3237
3238 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3239 v4sigma3tau[6] = 0.0e0;
3240
3241 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3242 v4sigma3tau[7] = 0.0e0;
3243
3244 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3245 v4sigma3tau[8] = 0.0e0;
3246
3247 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3248 v4sigma3tau[9] = 0.0e0;
3249
3250 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3251 v4sigma3tau[10] = 0.0e0;
3252
3253 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3254 v4sigma3tau[11] = 0.0e0;
3255
3256 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3257 v4sigma3tau[12] = 0.0e0;
3258
3259 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3260 v4sigma3tau[13] = 0.0e0;
3261
3262 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3263 v4sigma3tau[14] = 0.0e0;
3264
3265 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3266 v4sigma3tau[15] = 0.0e0;
3267
3268 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3269 v4sigma3tau[16] = 0.0e0;
3270
3271 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3272 v4sigma3tau[17] = 0.0e0;
3273
3274 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3275 v4sigma3tau[18] = 0.0e0;
3276
3277 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3278 v4sigma3tau[19] = 0.0e0;
3279
3280 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3281 v4sigma2lapl2[0] = 0.0e0;
3282
3283 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3284 v4sigma2lapl2[1] = 0.0e0;
3285
3286 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3287 v4sigma2lapl2[2] = 0.0e0;
3288
3289 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3290 v4sigma2lapl2[3] = 0.0e0;
3291
3292 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3293 v4sigma2lapl2[4] = 0.0e0;
3294
3295 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3296 v4sigma2lapl2[5] = 0.0e0;
3297
3298 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3299 v4sigma2lapl2[6] = 0.0e0;
3300
3301 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3302 v4sigma2lapl2[7] = 0.0e0;
3303
3304 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3305 v4sigma2lapl2[8] = 0.0e0;
3306
3307 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3308 v4sigma2lapl2[9] = 0.0e0;
3309
3310 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3311 v4sigma2lapl2[10] = 0.0e0;
3312
3313 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3314 v4sigma2lapl2[11] = 0.0e0;
3315
3316 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3317 v4sigma2lapl2[12] = 0.0e0;
3318
3319 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3320 v4sigma2lapl2[13] = 0.0e0;
3321
3322 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3323 v4sigma2lapl2[14] = 0.0e0;
3324
3325 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3326 v4sigma2lapl2[15] = 0.0e0;
3327
3328 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3329 v4sigma2lapl2[16] = 0.0e0;
3330
3331 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3332 v4sigma2lapl2[17] = 0.0e0;
3333
3334 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3335 v4sigma2lapltau[0] = 0.0e0;
3336
3337 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3338 v4sigma2lapltau[1] = 0.0e0;
3339
3340 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3341 v4sigma2lapltau[2] = 0.0e0;
3342
3343 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3344 v4sigma2lapltau[3] = 0.0e0;
3345
3346 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3347 v4sigma2lapltau[4] = 0.0e0;
3348
3349 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3350 v4sigma2lapltau[5] = 0.0e0;
3351
3352 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3353 v4sigma2lapltau[6] = 0.0e0;
3354
3355 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3356 v4sigma2lapltau[7] = 0.0e0;
3357
3358 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3359 v4sigma2lapltau[8] = 0.0e0;
3360
3361 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3362 v4sigma2lapltau[9] = 0.0e0;
3363
3364 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3365 v4sigma2lapltau[10] = 0.0e0;
3366
3367 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3368 v4sigma2lapltau[11] = 0.0e0;
3369
3370 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3371 v4sigma2lapltau[12] = 0.0e0;
3372
3373 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3374 v4sigma2lapltau[13] = 0.0e0;
3375
3376 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3377 v4sigma2lapltau[14] = 0.0e0;
3378
3379 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3380 v4sigma2lapltau[15] = 0.0e0;
3381
3382 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3383 v4sigma2lapltau[16] = 0.0e0;
3384
3385 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3386 v4sigma2lapltau[17] = 0.0e0;
3387
3388 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3389 v4sigma2lapltau[18] = 0.0e0;
3390
3391 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3392 v4sigma2lapltau[19] = 0.0e0;
3393
3394 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3395 v4sigma2lapltau[20] = 0.0e0;
3396
3397 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3398 v4sigma2lapltau[21] = 0.0e0;
3399
3400 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3401 v4sigma2lapltau[22] = 0.0e0;
3402
3403 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3404 v4sigma2lapltau[23] = 0.0e0;
3405
3406 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3407 v4sigma2tau2[0] = 0.0e0;
3408
3409 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3410 v4sigma2tau2[1] = 0.0e0;
3411
3412 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3413 v4sigma2tau2[2] = 0.0e0;
3414
3415 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3416 v4sigma2tau2[3] = 0.0e0;
3417
3418 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3419 v4sigma2tau2[4] = 0.0e0;
3420
3421 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3422 v4sigma2tau2[5] = 0.0e0;
3423
3424 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3425 v4sigma2tau2[6] = 0.0e0;
3426
3427 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3428 v4sigma2tau2[7] = 0.0e0;
3429
3430 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3431 v4sigma2tau2[8] = 0.0e0;
3432
3433 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3434 v4sigma2tau2[9] = 0.0e0;
3435
3436 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3437 v4sigma2tau2[10] = 0.0e0;
3438
3439 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3440 v4sigma2tau2[11] = 0.0e0;
3441
3442 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3443 v4sigma2tau2[12] = 0.0e0;
3444
3445 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3446 v4sigma2tau2[13] = 0.0e0;
3447
3448 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3449 v4sigma2tau2[14] = 0.0e0;
3450
3451 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3452 v4sigma2tau2[15] = 0.0e0;
3453
3454 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3455 v4sigma2tau2[16] = 0.0e0;
3456
3457 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3458 v4sigma2tau2[17] = 0.0e0;
3459
3460 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3461 v4sigmalapl3[0] = 0.0e0;
3462
3463 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3464 v4sigmalapl3[1] = 0.0e0;
3465
3466 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3467 v4sigmalapl3[2] = 0.0e0;
3468
3469 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3470 v4sigmalapl3[3] = 0.0e0;
3471
3472 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3473 v4sigmalapl3[4] = 0.0e0;
3474
3475 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3476 v4sigmalapl3[5] = 0.0e0;
3477
3478 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3479 v4sigmalapl3[6] = 0.0e0;
3480
3481 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3482 v4sigmalapl3[7] = 0.0e0;
3483
3484 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3485 v4sigmalapl3[8] = 0.0e0;
3486
3487 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3488 v4sigmalapl3[9] = 0.0e0;
3489
3490 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3491 v4sigmalapl3[10] = 0.0e0;
3492
3493 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3494 v4sigmalapl3[11] = 0.0e0;
3495
3496 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3497 v4sigmalapl2tau[0] = 0.0e0;
3498
3499 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3500 v4sigmalapl2tau[1] = 0.0e0;
3501
3502 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3503 v4sigmalapl2tau[2] = 0.0e0;
3504
3505 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3506 v4sigmalapl2tau[3] = 0.0e0;
3507
3508 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3509 v4sigmalapl2tau[4] = 0.0e0;
3510
3511 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3512 v4sigmalapl2tau[5] = 0.0e0;
3513
3514 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3515 v4sigmalapl2tau[6] = 0.0e0;
3516
3517 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3518 v4sigmalapl2tau[7] = 0.0e0;
3519
3520 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3521 v4sigmalapl2tau[8] = 0.0e0;
3522
3523 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3524 v4sigmalapl2tau[9] = 0.0e0;
3525
3526 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3527 v4sigmalapl2tau[10] = 0.0e0;
3528
3529 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3530 v4sigmalapl2tau[11] = 0.0e0;
3531
3532 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3533 v4sigmalapl2tau[12] = 0.0e0;
3534
3535 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3536 v4sigmalapl2tau[13] = 0.0e0;
3537
3538 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3539 v4sigmalapl2tau[14] = 0.0e0;
3540
3541 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3542 v4sigmalapl2tau[15] = 0.0e0;
3543
3544 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3545 v4sigmalapl2tau[16] = 0.0e0;
3546
3547 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3548 v4sigmalapl2tau[17] = 0.0e0;
3549
3550 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3551 v4sigmalapltau2[0] = 0.0e0;
3552
3553 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3554 v4sigmalapltau2[1] = 0.0e0;
3555
3556 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3557 v4sigmalapltau2[2] = 0.0e0;
3558
3559 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3560 v4sigmalapltau2[3] = 0.0e0;
3561
3562 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3563 v4sigmalapltau2[4] = 0.0e0;
3564
3565 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3566 v4sigmalapltau2[5] = 0.0e0;
3567
3568 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3569 v4sigmalapltau2[6] = 0.0e0;
3570
3571 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3572 v4sigmalapltau2[7] = 0.0e0;
3573
3574 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3575 v4sigmalapltau2[8] = 0.0e0;
3576
3577 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3578 v4sigmalapltau2[9] = 0.0e0;
3579
3580 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3581 v4sigmalapltau2[10] = 0.0e0;
3582
3583 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3584 v4sigmalapltau2[11] = 0.0e0;
3585
3586 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3587 v4sigmalapltau2[12] = 0.0e0;
3588
3589 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3590 v4sigmalapltau2[13] = 0.0e0;
3591
3592 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3593 v4sigmalapltau2[14] = 0.0e0;
3594
3595 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3596 v4sigmalapltau2[15] = 0.0e0;
3597
3598 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3599 v4sigmalapltau2[16] = 0.0e0;
3600
3601 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3602 v4sigmalapltau2[17] = 0.0e0;
3603
3604 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3605 v4sigmatau3[0] = 0.0e0;
3606
3607 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3608 v4sigmatau3[1] = 0.0e0;
3609
3610 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3611 v4sigmatau3[2] = 0.0e0;
3612
3613 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3614 v4sigmatau3[3] = 0.0e0;
3615
3616 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3617 v4sigmatau3[4] = 0.0e0;
3618
3619 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3620 v4sigmatau3[5] = 0.0e0;
3621
3622 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3623 v4sigmatau3[6] = 0.0e0;
3624
3625 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3626 v4sigmatau3[7] = 0.0e0;
3627
3628 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3629 v4sigmatau3[8] = 0.0e0;
3630
3631 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3632 v4sigmatau3[9] = 0.0e0;
3633
3634 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3635 v4sigmatau3[10] = 0.0e0;
3636
3637 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3638 v4sigmatau3[11] = 0.0e0;
3639
3640 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3641 v4lapl4[0] = 0.0e0;
3642
3643 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3644 v4lapl4[1] = 0.0e0;
3645
3646 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3647 v4lapl4[2] = 0.0e0;
3648
3649 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3650 v4lapl4[3] = 0.0e0;
3651
3652 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3653 v4lapl4[4] = 0.0e0;
3654
3655 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3656 v4lapl3tau[0] = 0.0e0;
3657
3658 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3659 v4lapl3tau[1] = 0.0e0;
3660
3661 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3662 v4lapl3tau[2] = 0.0e0;
3663
3664 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3665 v4lapl3tau[3] = 0.0e0;
3666
3667 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3668 v4lapl3tau[4] = 0.0e0;
3669
3670 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3671 v4lapl3tau[5] = 0.0e0;
3672
3673 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3674 v4lapl3tau[6] = 0.0e0;
3675
3676 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3677 v4lapl3tau[7] = 0.0e0;
3678
3679 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3680 v4lapl2tau2[0] = 0.0e0;
3681
3682 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3683 v4lapl2tau2[1] = 0.0e0;
3684
3685 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3686 v4lapl2tau2[2] = 0.0e0;
3687
3688 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3689 v4lapl2tau2[3] = 0.0e0;
3690
3691 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3692 v4lapl2tau2[4] = 0.0e0;
3693
3694 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3695 v4lapl2tau2[5] = 0.0e0;
3696
3697 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3698 v4lapl2tau2[6] = 0.0e0;
3699
3700 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3701 v4lapl2tau2[7] = 0.0e0;
3702
3703 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3704 v4lapl2tau2[8] = 0.0e0;
3705
3706 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3707 v4lapltau3[0] = 0.0e0;
3708
3709 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3710 v4lapltau3[1] = 0.0e0;
3711
3712 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3713 v4lapltau3[2] = 0.0e0;
3714
3715 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3716 v4lapltau3[3] = 0.0e0;
3717
3718 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3719 v4lapltau3[4] = 0.0e0;
3720
3721 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3722 v4lapltau3[5] = 0.0e0;
3723
3724 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3725 v4lapltau3[6] = 0.0e0;
3726
3727 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
3728 v4lapltau3[7] = 0.0e0;
3729
3730 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3731 v4tau4[0] = 0.0e0;
3732
3733 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3734 v4tau4[1] = 0.0e0;
3735
3736 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3737 v4tau4[2] = 0.0e0;
3738
3739 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3740 v4tau4[3] = 0.0e0;
3741
3742 if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
3743 v4tau4[4] = 0.0e0;
3744
3745 #ifndef XC_DONT_COMPILE_MXC
3746
3747 if(order < 5) return;
3748
3749
3750 #endif
3751
3752 #endif
3753
3754 #endif
3755
3756 #endif
3757
3758 #endif
3759
3760
3761 }
3762
3763