text: /*
    This file was generated automatically with ./scripts/maple2c.pl.
    Do not edit this file directly as it can be overwritten!!
  
    This Source Code Form is subject to the terms of the Mozilla Public
    License, v. 2.0. If a copy of the MPL was not distributed with this
    file, You can obtain one at http://mozilla.org/MPL/2.0/.
  
    Maple version     : Maple 2020 (X86 64 LINUX)
    Maple source      : ./maple/mgga_exc/mgga_xc_zlp.mpl
    Type of functional: mgga_exc
  */
  
  #define maple2c_order 4
  #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)
  
  
  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(XC_COMMA double *, ))
  {
  
  #ifndef XC_DONT_COMPILE_EXC
    double t2, t4, t5, t6, t7, t10, t11, t12;
    double t14, t17, t24, t25, t27, t28, t31, t33;
    double t34, t35, t36, t37, t38;
  
  #ifndef XC_DONT_COMPILE_VXC
    double t41, t42, t44, t49, t50, t53, t55, t58;
    double t61, t64, t67, t68;
  
  #ifndef XC_DONT_COMPILE_FXC
    double t71, t75, t78, t81, t83, t88, t89, t98;
    double t99, t100, t105, t109, t111;
  
  #ifndef XC_DONT_COMPILE_KXC
    double t116, t119, t125, t130, t133, t135, t140, t141;
    double t152, t156, t161, t165, t167, t170;
  
  #ifndef XC_DONT_COMPILE_LXC
    double t218, t224;
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  
  
    t2 = M_CBRT3;
    t4 = POW_1_3(0.1e1 / M_PI);
    t5 = t2 * t4;
    t6 = M_CBRT4;
    t7 = t6 * t6;
    t10 = rho[0] * rho[0];
    t11 = POW_1_3(rho[0]);
    t12 = t11 * t11;
    t14 = 0.1e1 / t12 / t10;
    t17 = 0.1e1 / t12 / rho[0];
    t24 = 0.20710800000000000000e0 * t5 * t7 + 0.53877250000000000000e-2 * t5 * t7 * (-lapl[0] * t17 / 0.8e1 + sigma[0] * t14 / 0.8e1);
    t25 = 0.1e1 / t11;
    t27 = 0.1e1 + 0.48849425066691677572e3 * t25;
    t28 = log(t27);
    t31 = 0.1e1 - 0.20471070000000000000e-2 * t28 * t11;
    t33 = t2 * t2;
    t34 = t24 * t31 * t33;
    t35 = 0.1e1 / t4;
    t36 = t35 * t6;
    t37 = t36 * t11;
    t38 = t34 * t37;
    if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
      zk[0] = -t38 / 0.3e1;
  
  #ifndef XC_DONT_COMPILE_VXC
  
    if(order < 1) return;
  
  
    t41 = t11 * rho[0];
    t42 = t10 * rho[0];
    t44 = 0.1e1 / t12 / t42;
    t49 = -sigma[0] * t44 / 0.3e1 + 0.5e1 / 0.24e2 * lapl[0] * t14;
    t50 = t41 * t49;
    t53 = t41 * t24;
    t55 = 0.1e1 / t27;
    t58 = 0.1e1 / t12;
    t61 = 0.33333333333333333332e0 / rho[0] * t55 - 0.68236900000000000000e-3 * t28 * t58;
    t64 = t33 * t35 * t6;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vrho[0] = -0.4e1 / 0.9e1 * t38 - 0.21550900000000000000e-1 * t50 * t31 - t53 * t61 * t64 / 0.3e1;
  
    t67 = 0.1e1 / t41;
    t68 = t67 * t31;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vsigma[0] = -0.26938625000000000000e-2 * t68;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vlapl[0] = 0.26938625000000000000e-2 * t25 * t31;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vtau[0] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_FXC
  
    if(order < 2) return;
  
  
    t71 = t49 * t31;
    t75 = t24 * t61 * t33;
    t78 = t36 * t58;
    t81 = t10 * t10;
    t83 = 0.1e1 / t12 / t81;
    t88 = 0.11e2 / 0.9e1 * sigma[0] * t83 - 0.5e1 / 0.9e1 * lapl[0] * t44;
    t89 = t41 * t88;
    t98 = 0.1e1 / t11 / t10;
    t99 = t27 * t27;
    t100 = 0.1e1 / t99;
    t105 = -0.22222222222222222221e0 / t10 * t55 + 0.54277138962990752854e2 * t98 * t100 + 0.45491266666666666667e-3 * t28 * t17;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rho2[0] = -0.57469066666666666666e-1 * t71 * t11 - 0.8e1 / 0.9e1 * t75 * t37 - 0.4e1 / 0.27e2 * t34 * t78 - 0.21550900000000000000e-1 * t89 * t31 - 0.43101800000000000000e-1 * t50 * t61 - t53 * t105 * t64 / 0.3e1;
  
    t109 = t98 * t31;
    t111 = t67 * t61;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[0] = 0.35918166666666666667e-2 * t109 - 0.26938625000000000000e-2 * t111;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rholapl[0] = -0.89795416666666666667e-3 * t68 + 0.26938625000000000000e-2 * t25 * t61;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhotau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapl2[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapltau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2tau2[0] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_KXC
  
    if(order < 3) return;
  
  
    t116 = t88 * t31;
    t119 = t49 * t61;
    t125 = t24 * t105 * t33;
    t130 = t36 * t17;
    t133 = t81 * rho[0];
    t135 = 0.1e1 / t12 / t133;
    t140 = -0.154e3 / 0.27e2 * sigma[0] * t135 + 0.55e2 / 0.27e2 * lapl[0] * t83;
    t141 = t41 * t140;
    t152 = 0.1e1 / t11 / t42;
    t156 = 0.1e1 / t99 / t27;
    t161 = 0.37037037037037037035e0 / t42 * t55 - 0.16283141688897225856e3 * t152 * t100 + 0.17676046884046853390e5 * t44 * t156 - 0.75818777777777777778e-3 * t28 * t14;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho3[0] = -0.86203599999999999999e-1 * t116 * t11 - 0.17240720000000000000e0 * t119 * t11 - 0.28734533333333333333e-1 * t71 * t58 - 0.4e1 / 0.3e1 * t125 * t37 - 0.4e1 / 0.9e1 * t75 * t78 + 0.8e1 / 0.81e2 * t34 * t130 - 0.21550900000000000000e-1 * t141 * t31 - 0.64652700000000000000e-1 * t89 * t61 - 0.64652700000000000000e-1 * t50 * t105 - t53 * t161 * t64 / 0.3e1;
  
    t165 = t152 * t31;
    t167 = t98 * t61;
    t170 = 0.26938625000000000000e-2 * t67 * t105;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[0] = -0.83809055555555555556e-2 * t165 + 0.71836333333333333334e-2 * t167 - t170;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[0] = 0.11972722222222222222e-2 * t109 - 0.17959083333333333333e-2 * t111 + 0.26938625000000000000e-2 * t25 * t105;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl3[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3tau3[0] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_LXC
  
    if(order < 4) return;
  
  
    t218 = 0.1e1 / t11 / t81;
    t224 = t99 * t99;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[0] = -0.86203600000000000000e-1 * t50 * t161 - 0.11493813333333333333e0 * t140 * t31 * t11 - 0.34481440000000000000e0 * t88 * t61 * t11 - 0.34481440000000000000e0 * t49 * t105 * t11 - 0.21550900000000000000e-1 * t41 * (0.2618e4 / 0.81e2 * sigma[0] / t12 / t81 / t10 - 0.770e3 / 0.81e2 * lapl[0] * t135) * t31 - 0.86203600000000000000e-1 * t141 * t61 - 0.12930540000000000000e0 * t89 * t105 - 0.57469066666666666666e-1 * t116 * t58 - 0.11493813333333333333e0 * t119 * t58 + 0.25541807407407407407e-1 * t71 * t17 - 0.8e1 / 0.9e1 * t125 * t78 + 0.32e2 / 0.81e2 * t75 * t130 - 0.40e2 / 0.243e3 * t34 * t36 * t14 - t53 * (-0.98765432098765432088e0 / t81 * t55 + 0.60307932181100836503e3 * t218 * t100 - 0.11784031256031235593e6 * t83 * t156 + 0.86346472773757568068e7 / t133 / t224 + 0.20218340740740740741e-2 * t28 * t44) * t64 / 0.3e1 - 0.16e2 / 0.9e1 * t24 * t161 * t33 * t37;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[0] = 0.27936351851851851852e-1 * t218 * t31 - 0.25142716666666666667e-1 * t152 * t61 + 0.10775450000000000000e-1 * t98 * t105 - 0.26938625000000000000e-2 * t67 * t161;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[0] = -0.27936351851851851851e-2 * t165 + 0.35918166666666666666e-2 * t167 - t170 + 0.26938625000000000000e-2 * t25 * t161;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[0] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_MXC
  
    if(order < 5) return;
  
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  
  }
  
  
  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(XC_COMMA double *, ))
  {
  
  #ifndef XC_DONT_COMPILE_EXC
    double t2, t4, t5, t6, t7, t11, t12, t13;
    double t14, t15, t17, t19, t20, t22, t23, t24;
    double t25, t26, t28, t29, t30, t31, t33, t34;
    double t36, t37, t39, t40, t41, t42, t49, t52;
    double t53, t56, t58, t59, t60, t61, t62, t63;
  
  #ifndef XC_DONT_COMPILE_VXC
    double t65, t66, t67, t69, t71, t72, t74, t75;
    double t78, t79, t81, t82, t85, t86, t89, t90;
    double t93, t94, t97, t100, t103, t105, t107, t108;
    double t111, t113, t114, t117, t118, t121, t122, t125;
    double t126, t129, t130, t132, t133;
  
  #ifndef XC_DONT_COMPILE_FXC
    double t135, t136, t139, t141, t142, t144, t145, t147;
    double t149, t152, t153, t158, t159, t160, t163, t164;
    double t165, t166, t169, t170, t171, t174, t175, t178;
    double t179, t182, t187, t188, t189, t193, t196, t199;
    double t201, t202, t206, t207, t210, t211, t216, t217;
    double t220, t221, t224, t225, t228, t232, t233, t236;
    double t237, t242, t243, t248, t249, t252, t253, t256;
    double t257, t261, t263, t267, t269, t270, t273, t274;
    double t277, t279, t280, t282, t283, t284, t287, t289;
    double t290, t293, t296;
  
  #ifndef XC_DONT_COMPILE_KXC
    double t299, t300, t302, t303, t305, t308, t310, t312;
    double t313, t315, t316, t320, t321, t323, t324, t334;
    double t335, t336, t339, t340, t343, t344, t346, t347;
    double t351, t352, t353, t356, t357, t360, t361, t364;
    double t365, t368, t370, t375, t379, t384, t387, t391;
    double t392, t394, t396, t397, t406, t409, t410, t413;
    double t416, t419, t427, t430, t431, t434, t437, t440;
    double t443, t444, t448, t449, t457, t458, t462, t463;
    double t466, t467, t472, t473, t477, t478, t487, t488;
    double t491, t492, t495, t496, t499, t500, t503, t504;
    double t507, t514, t515, t518, t522, t523, t526, t528;
    double t529, t538, t539, t542, t545, t546, t549, t550;
    double t555, t557, t559, t564, t566, t567, t568, t570;
    double t573, t574, t579, t581, t582, t585, t586, t587;
    double t589, t592, t594, t595, t597, t598, t601, t602;
    double t603, t606, t607, t608, t610, t613, t615, t618;
    double t622, t623, t624, t627, t628, t630, t631, t637;
    double t638, t642, t644, t646, t647, t648, t651, t652;
    double t654, t655, t660, t664, t669, t675, t679;
  
  #ifndef XC_DONT_COMPILE_LXC
    double t682, t684, t686, t690, t703, t704, t706, t714;
    double t716, t717, t725, t732, t734, t735, t743, t756;
    double t760, t763, t766, t769, t771, t773, t775, t778;
    double t780, t783, t787, t791, t796, t805, t806, t808;
    double t810, t811, t812, t815, t817, t818, t819, t850;
    double t861, t864, t882, t883, t893, t894, t904, t907;
    double t913, t914, t918, t923, t928, t975, t993, t995;
    double t1015, t1021, t1026, t1028, t1031, t1033, t1036, t1045;
    double t1055, t1057, t1060, t1064, t1097, t1137, t1148, t1155;
    double t1159, t1177, t1184, t1197, t1204, t1211, t1224, t1226;
    double t1228, t1230, t1237, t1239, t1242, t1249, t1251, t1252;
    double t1255, t1256, t1258, t1260, t1262, t1264, t1265, t1267;
    double t1271, t1275, t1277, t1279, t1281, t1284, t1289, t1297;
    double t1301, t1306, t1307, t1309, t1311, t1313, t1315, t1317;
    double t1319, t1322, t1325, t1327, t1330, t1331, t1339, t1355;
    double t1360, t1365, t1369, t1370, t1376, t1380, t1383, t1385;
    double t1389, t1390, t1393, t1395, t1402, t1406, t1408, t1410;
    double t1412, t1416, t1422, t1425, t1431, t1439, t1441, t1444;
    double t1447, t1450, t1454, t1457, t1459, t1478, t1484, t1487;
    double t1492, t1496, t1500, t1505, t1507, t1526, t1534, t1540;
    double t1543, t1548, t1550, t1552, t1554, t1583, t1604;
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  
  
    t2 = M_CBRT3;
    t4 = POW_1_3(0.1e1 / M_PI);
    t5 = t2 * t4;
    t6 = M_CBRT4;
    t7 = t6 * t6;
    t11 = sigma[0] + 0.2e1 * sigma[1] + sigma[2];
    t12 = rho[0] + rho[1];
    t13 = t12 * t12;
    t14 = POW_1_3(t12);
    t15 = t14 * t14;
    t17 = 0.1e1 / t15 / t13;
    t19 = POW_1_3(rho[0]);
    t20 = t19 * t19;
    t22 = 0.1e1 / t20 / rho[0];
    t23 = lapl[0] * t22;
    t24 = rho[0] - rho[1];
    t25 = 0.1e1 / t12;
    t26 = t24 * t25;
    t28 = 0.1e1 / 0.2e1 + t26 / 0.2e1;
    t29 = POW_1_3(t28);
    t30 = t29 * t29;
    t31 = t30 * t28;
    t33 = POW_1_3(rho[1]);
    t34 = t33 * t33;
    t36 = 0.1e1 / t34 / rho[1];
    t37 = lapl[1] * t36;
    t39 = 0.1e1 / 0.2e1 - t26 / 0.2e1;
    t40 = POW_1_3(t39);
    t41 = t40 * t40;
    t42 = t41 * t39;
    t49 = 0.20710800000000000000e0 * t5 * t7 + 0.53877250000000000000e-2 * t5 * t7 * (t11 * t17 / 0.8e1 - t23 * t31 / 0.8e1 - t37 * t42 / 0.8e1);
    t52 = 0.1e1 + 0.48849425066691677572e3 / t14;
    t53 = log(t52);
    t56 = 0.1e1 - 0.20471070000000000000e-2 * t53 * t14;
    t58 = t2 * t2;
    t59 = t49 * t56 * t58;
    t60 = 0.1e1 / t4;
    t61 = t60 * t6;
    t62 = t61 * t14;
    t63 = t59 * t62;
    if(zk != NULL && (p->info->flags & XC_FLAGS_HAVE_EXC))
      zk[0] = -t63 / 0.3e1;
  
  #ifndef XC_DONT_COMPILE_VXC
  
    if(order < 1) return;
  
  
    t65 = 0.4e1 / 0.9e1 * t63;
    t66 = t14 * t12;
    t67 = t13 * t12;
    t69 = 0.1e1 / t15 / t67;
    t71 = t11 * t69 / 0.3e1;
    t72 = rho[0] * rho[0];
    t74 = 0.1e1 / t20 / t72;
    t75 = lapl[0] * t74;
    t78 = 0.1e1 / t13;
    t79 = t24 * t78;
    t81 = t25 / 0.2e1 - t79 / 0.2e1;
    t82 = t30 * t81;
    t85 = -t81;
    t86 = t41 * t85;
    t89 = -t71 + 0.5e1 / 0.24e2 * t75 * t31 - 0.5e1 / 0.24e2 * t23 * t82 - 0.5e1 / 0.24e2 * t37 * t86;
    t90 = t66 * t89;
    t93 = t66 * t49;
    t94 = 0.1e1 / t52;
    t97 = 0.1e1 / t15;
    t100 = 0.33333333333333333332e0 * t25 * t94 - 0.68236900000000000000e-3 * t53 * t97;
    t103 = t58 * t60 * t6;
    t105 = t93 * t100 * t103 / 0.3e1;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vrho[0] = -t65 - 0.21550900000000000000e-1 * t90 * t56 - t105;
  
    t107 = -t25 / 0.2e1 - t79 / 0.2e1;
    t108 = t30 * t107;
    t111 = rho[1] * rho[1];
    t113 = 0.1e1 / t34 / t111;
    t114 = lapl[1] * t113;
    t117 = -t107;
    t118 = t41 * t117;
    t121 = -t71 - 0.5e1 / 0.24e2 * t23 * t108 + 0.5e1 / 0.24e2 * t114 * t42 - 0.5e1 / 0.24e2 * t37 * t118;
    t122 = t66 * t121;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vrho[1] = -t65 - 0.21550900000000000000e-1 * t122 * t56 - t105;
  
    t125 = 0.1e1 / t66;
    t126 = t125 * t56;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vsigma[0] = -0.26938625000000000000e-2 * t126;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vsigma[1] = -0.53877250000000000000e-2 * t126;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vsigma[2] = vsigma[0];
  
    t129 = t66 * t22;
    t130 = t31 * t56;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vlapl[0] = 0.26938625000000000000e-2 * t129 * t130;
  
    t132 = t66 * t36;
    t133 = t42 * t56;
    if(vrho != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vlapl[1] = 0.26938625000000000000e-2 * t132 * t133;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vtau[0] = 0.0e0;
  
    if(vrho != NULL && (p->info->flags & XC_FLAGS_HAVE_VXC))
      vtau[1] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_FXC
  
    if(order < 2) return;
  
  
    t135 = t89 * t56;
    t136 = t135 * t14;
    t139 = t49 * t100 * t58;
    t141 = 0.8e1 / 0.9e1 * t139 * t62;
    t142 = t61 * t97;
    t144 = 0.4e1 / 0.27e2 * t59 * t142;
    t145 = t13 * t13;
    t147 = 0.1e1 / t15 / t145;
    t149 = 0.11e2 / 0.9e1 * t11 * t147;
    t152 = 0.1e1 / t20 / t72 / rho[0];
    t153 = lapl[0] * t152;
    t158 = 0.1e1 / t29;
    t159 = t81 * t81;
    t160 = t158 * t159;
    t163 = 0.1e1 / t67;
    t164 = t24 * t163;
    t165 = -t78 + t164;
    t166 = t30 * t165;
    t169 = 0.1e1 / t40;
    t170 = t85 * t85;
    t171 = t169 * t170;
    t174 = -t165;
    t175 = t41 * t174;
    t178 = t149 - 0.5e1 / 0.9e1 * t153 * t31 + 0.25e2 / 0.36e2 * t75 * t82 - 0.5e1 / 0.36e2 * t23 * t160 - 0.5e1 / 0.24e2 * t23 * t166 - 0.5e1 / 0.36e2 * t37 * t171 - 0.5e1 / 0.24e2 * t37 * t175;
    t179 = t66 * t178;
    t182 = t90 * t100;
    t187 = 0.1e1 / t14 / t13;
    t188 = t52 * t52;
    t189 = 0.1e1 / t188;
    t193 = 0.1e1 / t15 / t12;
    t196 = -0.22222222222222222221e0 * t78 * t94 + 0.54277138962990752854e2 * t187 * t189 + 0.45491266666666666667e-3 * t53 * t193;
    t199 = t93 * t196 * t103 / 0.3e1;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rho2[0] = -0.57469066666666666666e-1 * t136 - t141 - t144 - 0.21550900000000000000e-1 * t179 * t56 - 0.43101800000000000000e-1 * t182 - t199;
  
    t201 = t14 * t121;
    t202 = t201 * t56;
    t206 = t158 * t107;
    t207 = t206 * t81;
    t210 = t30 * t24;
    t211 = t210 * t163;
    t216 = t169 * t117;
    t217 = t216 * t85;
    t220 = t41 * t24;
    t221 = t220 * t163;
    t224 = t149 + 0.25e2 / 0.72e2 * t75 * t108 - 0.5e1 / 0.36e2 * t23 * t207 - 0.5e1 / 0.24e2 * t23 * t211 + 0.25e2 / 0.72e2 * t114 * t86 - 0.5e1 / 0.36e2 * t37 * t217 + 0.5e1 / 0.24e2 * t37 * t221;
    t225 = t66 * t224;
    t228 = t122 * t100;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rho2[1] = -0.28734533333333333333e-1 * t136 - t141 - t144 - 0.28734533333333333333e-1 * t202 - 0.21550900000000000000e-1 * t225 * t56 - 0.21550900000000000000e-1 * t228 - 0.21550900000000000000e-1 * t182 - t199;
  
    t232 = t107 * t107;
    t233 = t158 * t232;
    t236 = t78 + t164;
    t237 = t30 * t236;
    t242 = 0.1e1 / t34 / t111 / rho[1];
    t243 = lapl[1] * t242;
    t248 = t117 * t117;
    t249 = t169 * t248;
    t252 = -t236;
    t253 = t41 * t252;
    t256 = t149 - 0.5e1 / 0.36e2 * t23 * t233 - 0.5e1 / 0.24e2 * t23 * t237 - 0.5e1 / 0.9e1 * t243 * t42 + 0.25e2 / 0.36e2 * t114 * t118 - 0.5e1 / 0.36e2 * t37 * t249 - 0.5e1 / 0.24e2 * t37 * t253;
    t257 = t66 * t256;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rho2[2] = -0.57469066666666666666e-1 * t202 - t141 - t144 - 0.21550900000000000000e-1 * t257 * t56 - 0.43101800000000000000e-1 * t228 - t199;
  
    t261 = t187 * t56;
    t263 = t125 * t100;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[0] = 0.35918166666666666667e-2 * t261 - 0.26938625000000000000e-2 * t263;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[1] = 0.71836333333333333333e-2 * t261 - 0.53877250000000000000e-2 * t263;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[2] = v2rhosigma[0];
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[3] = v2rhosigma[2];
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[4] = v2rhosigma[1];
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhosigma[5] = v2rhosigma[3];
  
    t267 = t14 * t22;
    t269 = 0.35918166666666666667e-2 * t267 * t130;
    t270 = t66 * t74;
    t273 = t30 * t56;
    t274 = t273 * t81;
    t277 = t31 * t100;
    t279 = 0.26938625000000000000e-2 * t129 * t277;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rholapl[0] = t269 - 0.44897708333333333333e-2 * t270 * t130 + 0.44897708333333333333e-2 * t129 * t274 + t279;
  
    t280 = t14 * t36;
    t282 = 0.35918166666666666667e-2 * t280 * t133;
    t283 = t41 * t56;
    t284 = t283 * t85;
    t287 = t42 * t100;
    t289 = 0.26938625000000000000e-2 * t132 * t287;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rholapl[1] = t282 + 0.44897708333333333333e-2 * t132 * t284 + t289;
  
    t290 = t273 * t107;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rholapl[2] = t269 + 0.44897708333333333333e-2 * t129 * t290 + t279;
  
    t293 = t66 * t113;
    t296 = t283 * t117;
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rholapl[3] = t282 - 0.44897708333333333333e-2 * t293 * t133 + 0.44897708333333333333e-2 * t132 * t296 + t289;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhotau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhotau[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhotau[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2rhotau[3] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[3] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[4] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigma2[5] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[3] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[4] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmalapl[5] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[3] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[4] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2sigmatau[5] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapl2[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapl2[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapl2[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapltau[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapltau[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapltau[2] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2lapltau[3] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2tau2[0] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2tau2[1] = 0.0e0;
  
    if(v2rho2 != NULL && (p->info->flags & XC_FLAGS_HAVE_FXC))
      v2tau2[2] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_KXC
  
    if(order < 3) return;
  
  
    t299 = t178 * t56;
    t300 = t299 * t14;
    t302 = t89 * t100;
    t303 = t302 * t14;
    t305 = t135 * t97;
    t308 = t49 * t196 * t58;
    t310 = 0.4e1 / 0.3e1 * t308 * t62;
    t312 = 0.4e1 / 0.9e1 * t139 * t142;
    t313 = t61 * t193;
    t315 = 0.8e1 / 0.81e2 * t59 * t313;
    t316 = t145 * t12;
    t320 = 0.154e3 / 0.27e2 * t11 / t15 / t316;
    t321 = t72 * t72;
    t323 = 0.1e1 / t20 / t321;
    t324 = lapl[0] * t323;
    t334 = 0.1e1 / t29 / t28;
    t335 = t159 * t81;
    t336 = t334 * t335;
    t339 = t158 * t81;
    t340 = t339 * t165;
    t343 = 0.1e1 / t145;
    t344 = t24 * t343;
    t346 = 0.3e1 * t163 - 0.3e1 * t344;
    t347 = t30 * t346;
    t351 = 0.1e1 / t40 / t39;
    t352 = t170 * t85;
    t353 = t351 * t352;
    t356 = t169 * t85;
    t357 = t356 * t174;
    t360 = -t346;
    t361 = t41 * t360;
    t364 = -t320 + 0.55e2 / 0.27e2 * t324 * t31 - 0.25e2 / 0.9e1 * t153 * t82 + 0.25e2 / 0.36e2 * t75 * t160 + 0.25e2 / 0.24e2 * t75 * t166 + 0.5e1 / 0.108e3 * t23 * t336 - 0.5e1 / 0.12e2 * t23 * t340 - 0.5e1 / 0.24e2 * t23 * t347 + 0.5e1 / 0.108e3 * t37 * t353 - 0.5e1 / 0.12e2 * t37 * t357 - 0.5e1 / 0.24e2 * t37 * t361;
    t365 = t66 * t364;
    t368 = t179 * t100;
    t370 = t90 * t196;
    t375 = 0.1e1 / t14 / t67;
    t379 = 0.1e1 / t188 / t52;
    t384 = 0.37037037037037037035e0 * t163 * t94 - 0.16283141688897225856e3 * t375 * t189 + 0.17676046884046853390e5 * t69 * t379 - 0.75818777777777777778e-3 * t53 * t17;
    t387 = t93 * t384 * t103 / 0.3e1;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho3[0] = -0.86203599999999999999e-1 * t300 - 0.17240720000000000000e0 * t303 - 0.28734533333333333333e-1 * t305 - t310 - t312 + t315 - 0.21550900000000000000e-1 * t365 * t56 - 0.64652700000000000000e-1 * t368 - 0.64652700000000000000e-1 * t370 - t387;
  
    t391 = t97 * t121;
    t392 = t391 * t56;
    t394 = t14 * t224;
    t396 = 0.57469066666666666666e-1 * t394 * t56;
    t397 = t201 * t100;
    t406 = t334 * t107 * t159;
    t409 = t23 * t158;
    t410 = t164 * t81;
    t413 = t206 * t165;
    t416 = t30 * t163;
    t419 = t210 * t343;
    t427 = t351 * t117 * t170;
    t430 = t37 * t169;
    t431 = t164 * t85;
    t434 = t216 * t174;
    t437 = t41 * t163;
    t440 = t220 * t343;
    t443 = -t320 - 0.25e2 / 0.27e2 * t153 * t108 + 0.25e2 / 0.54e2 * t75 * t207 + 0.25e2 / 0.36e2 * t75 * t211 + 0.5e1 / 0.108e3 * t23 * t406 - 0.5e1 / 0.18e2 * t409 * t410 - 0.5e1 / 0.36e2 * t23 * t413 - 0.5e1 / 0.24e2 * t23 * t416 + 0.5e1 / 0.8e1 * t23 * t419 + 0.25e2 / 0.108e3 * t114 * t171 + 0.25e2 / 0.72e2 * t114 * t175 + 0.5e1 / 0.108e3 * t37 * t427 + 0.5e1 / 0.18e2 * t430 * t431 - 0.5e1 / 0.36e2 * t37 * t434 + 0.5e1 / 0.24e2 * t37 * t437 - 0.5e1 / 0.8e1 * t37 * t440;
    t444 = t66 * t443;
    t448 = 0.43101800000000000000e-1 * t225 * t100;
    t449 = t122 * t196;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho3[1] = -0.28734533333333333333e-1 * t300 - 0.11493813333333333333e0 * t303 - 0.19156355555555555556e-1 * t305 - t310 - t312 + t315 - 0.95781777777777777777e-2 * t392 - t396 - 0.57469066666666666666e-1 * t397 - 0.21550900000000000000e-1 * t444 * t56 - t448 - 0.21550900000000000000e-1 * t449 - 0.21550900000000000000e-1 * t368 - 0.43101800000000000000e-1 * t370 - t387;
  
    t457 = t14 * t256;
    t458 = t457 * t56;
    t462 = t334 * t232;
    t463 = t462 * t81;
    t466 = t107 * t24;
    t467 = t466 * t163;
    t472 = t158 * t236;
    t473 = t472 * t81;
    t477 = -t163 - 0.3e1 * t344;
    t478 = t30 * t477;
    t487 = t351 * t248;
    t488 = t487 * t85;
    t491 = t117 * t24;
    t492 = t491 * t163;
    t495 = t169 * t252;
    t496 = t495 * t85;
    t499 = -t477;
    t500 = t41 * t499;
    t503 = -t320 + 0.25e2 / 0.108e3 * t75 * t233 + 0.5e1 / 0.108e3 * t23 * t463 - 0.5e1 / 0.18e2 * t409 * t467 + 0.25e2 / 0.72e2 * t75 * t237 - 0.5e1 / 0.36e2 * t23 * t473 - 0.5e1 / 0.24e2 * t23 * t478 - 0.25e2 / 0.27e2 * t243 * t86 + 0.25e2 / 0.54e2 * t114 * t217 - 0.25e2 / 0.36e2 * t114 * t221 + 0.5e1 / 0.108e3 * t37 * t488 + 0.5e1 / 0.18e2 * t430 * t492 - 0.5e1 / 0.36e2 * t37 * t496 - 0.5e1 / 0.24e2 * t37 * t500;
    t504 = t66 * t503;
    t507 = t257 * t100;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho3[2] = -0.19156355555555555555e-1 * t392 - t396 - 0.11493813333333333333e0 * t397 - 0.57469066666666666667e-1 * t303 - t310 - t312 - 0.95781777777777777778e-2 * t305 + t315 - 0.28734533333333333333e-1 * t458 - 0.21550900000000000000e-1 * t504 * t56 - 0.21550900000000000000e-1 * t507 - t448 - 0.43101800000000000000e-1 * t449 - 0.21550900000000000000e-1 * t370 - t387;
  
    t514 = t232 * t107;
    t515 = t334 * t514;
    t518 = t206 * t236;
    t522 = -0.3e1 * t163 - 0.3e1 * t344;
    t523 = t30 * t522;
    t526 = t111 * t111;
    t528 = 0.1e1 / t34 / t526;
    t529 = lapl[1] * t528;
    t538 = t248 * t117;
    t539 = t351 * t538;
    t542 = t216 * t252;
    t545 = -t522;
    t546 = t41 * t545;
    t549 = -t320 + 0.5e1 / 0.108e3 * t23 * t515 - 0.5e1 / 0.12e2 * t23 * t518 - 0.5e1 / 0.24e2 * t23 * t523 + 0.55e2 / 0.27e2 * t529 * t42 - 0.25e2 / 0.9e1 * t243 * t118 + 0.25e2 / 0.36e2 * t114 * t249 + 0.25e2 / 0.24e2 * t114 * t253 + 0.5e1 / 0.108e3 * t37 * t539 - 0.5e1 / 0.12e2 * t37 * t542 - 0.5e1 / 0.24e2 * t37 * t546;
    t550 = t66 * t549;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho3[3] = -0.28734533333333333333e-1 * t392 - 0.86203599999999999999e-1 * t458 - 0.17240720000000000000e0 * t397 - t310 - t312 + t315 - 0.21550900000000000000e-1 * t550 * t56 - 0.64652700000000000000e-1 * t507 - 0.64652700000000000000e-1 * t449 - t387;
  
    t555 = t375 * t56;
    t557 = t187 * t100;
    t559 = t125 * t196;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[0] = -0.83809055555555555556e-2 * t555 + 0.71836333333333333334e-2 * t557 - 0.26938625000000000000e-2 * t559;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[1] = -0.16761811111111111111e-1 * t555 + 0.14367266666666666667e-1 * t557 - 0.53877250000000000000e-2 * t559;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[2] = v3rho2sigma[0];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[3] = v3rho2sigma[2];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[4] = v3rho2sigma[1];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[5] = v3rho2sigma[3];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[6] = v3rho2sigma[5];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[7] = v3rho2sigma[4];
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2sigma[8] = v3rho2sigma[6];
  
    t564 = t97 * t22;
    t566 = 0.11972722222222222222e-2 * t564 * t130;
    t567 = t14 * t74;
    t568 = t567 * t130;
    t570 = t267 * t274;
    t573 = 0.71836333333333333334e-2 * t267 * t277;
    t574 = t66 * t152;
    t579 = t270 * t277;
    t581 = t158 * t56;
    t582 = t581 * t159;
    t585 = t30 * t100;
    t586 = t585 * t81;
    t587 = t129 * t586;
    t589 = t273 * t165;
    t592 = t31 * t196;
    t594 = 0.26938625000000000000e-2 * t129 * t592;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[0] = t566 - 0.11972722222222222222e-1 * t568 + 0.11972722222222222222e-1 * t570 + t573 + 0.11972722222222222222e-1 * t574 * t130 - 0.14965902777777777778e-1 * t270 * t274 - 0.89795416666666666666e-2 * t579 + 0.29931805555555555555e-2 * t129 * t582 + 0.89795416666666666666e-2 * t587 + 0.44897708333333333333e-2 * t129 * t589 + t594;
  
    t595 = t97 * t36;
    t597 = 0.11972722222222222222e-2 * t595 * t133;
    t598 = t280 * t284;
    t601 = 0.71836333333333333334e-2 * t280 * t287;
    t602 = t169 * t56;
    t603 = t602 * t170;
    t606 = t41 * t100;
    t607 = t606 * t85;
    t608 = t132 * t607;
    t610 = t283 * t174;
    t613 = t42 * t196;
    t615 = 0.26938625000000000000e-2 * t132 * t613;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[1] = t597 + 0.11972722222222222222e-1 * t598 + t601 + 0.29931805555555555555e-2 * t132 * t603 + 0.89795416666666666666e-2 * t608 + 0.44897708333333333333e-2 * t132 * t610 + t615;
  
    t618 = t267 * t290;
    t622 = t129 * t158;
    t623 = t56 * t107;
    t624 = t623 * t81;
    t627 = t585 * t107;
    t628 = t129 * t627;
    t630 = t193 * t22;
    t631 = t273 * t24;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[2] = t566 - 0.59863611111111111112e-2 * t568 + 0.59863611111111111112e-2 * t570 + t573 + 0.59863611111111111111e-2 * t618 - 0.74829513888888888888e-2 * t270 * t290 + 0.29931805555555555555e-2 * t622 * t624 + 0.44897708333333333333e-2 * t628 + 0.44897708333333333333e-2 * t630 * t631 - 0.44897708333333333333e-2 * t579 + 0.44897708333333333333e-2 * t587 + t594;
  
    t637 = t14 * t113;
    t638 = t637 * t133;
    t642 = t293 * t287;
    t644 = t280 * t296;
    t646 = t132 * t169;
    t647 = t56 * t117;
    t648 = t647 * t85;
    t651 = t606 * t117;
    t652 = t132 * t651;
    t654 = t193 * t36;
    t655 = t283 * t24;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[3] = t597 + 0.59863611111111111112e-2 * t598 + t601 - 0.59863611111111111111e-2 * t638 - 0.74829513888888888888e-2 * t293 * t284 - 0.44897708333333333333e-2 * t642 + 0.59863611111111111111e-2 * t644 + 0.29931805555555555555e-2 * t646 * t648 + 0.44897708333333333333e-2 * t652 - 0.44897708333333333333e-2 * t654 * t655 + 0.44897708333333333333e-2 * t608 + t615;
  
    t660 = t581 * t232;
    t664 = t273 * t236;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[4] = t566 + 0.11972722222222222222e-1 * t618 + t573 + 0.29931805555555555555e-2 * t129 * t660 + 0.89795416666666666666e-2 * t628 + 0.44897708333333333333e-2 * t129 * t664 + t594;
  
    t669 = t66 * t242;
    t675 = t602 * t248;
    t679 = t283 * t252;
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2lapl[5] = t597 - 0.11972722222222222222e-1 * t638 + 0.11972722222222222222e-1 * t644 + t601 + 0.11972722222222222222e-1 * t669 * t133 - 0.14965902777777777778e-1 * t293 * t296 - 0.89795416666666666666e-2 * t642 + 0.29931805555555555555e-2 * t132 * t675 + 0.89795416666666666666e-2 * t652 + 0.44897708333333333333e-2 * t132 * t679 + t615;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rho2tau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigma2[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmalapl[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhosigmatau[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapl2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rholapltau[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3rhotau2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma3[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2lapl[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigma2tau[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapl2[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[9] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[10] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmalapltau[11] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[6] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[7] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3sigmatau2[8] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl3[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl3[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl3[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl3[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapl2tau[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[3] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[4] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3lapltau2[5] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3tau3[0] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3tau3[1] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3tau3[2] = 0.0e0;
  
    if(v3rho3 != NULL && (p->info->flags & XC_FLAGS_HAVE_KXC))
      v3tau3[3] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_LXC
  
    if(order < 4) return;
  
  
    t682 = t179 * t196;
    t684 = t90 * t384;
    t686 = t145 * t13;
    t690 = 0.2618e4 / 0.81e2 * t11 / t15 / t686;
    t703 = 0.1e1 / t316;
    t704 = t24 * t703;
    t706 = -0.12e2 * t343 + 0.12e2 * t704;
    t714 = t28 * t28;
    t716 = 0.1e1 / t29 / t714;
    t717 = t159 * t159;
    t725 = t165 * t165;
    t732 = t39 * t39;
    t734 = 0.1e1 / t40 / t732;
    t735 = t170 * t170;
    t743 = t174 * t174;
    t756 = t690 - 0.770e3 / 0.81e2 * lapl[0] / t20 / t321 / rho[0] * t31 + 0.1100e4 / 0.81e2 * t324 * t82 - 0.50e2 / 0.9e1 * t153 * t166 + 0.25e2 / 0.18e2 * t75 * t347 - 0.5e1 / 0.24e2 * t23 * t30 * t706 + 0.5e1 / 0.24e2 * t37 * t41 * t706 - 0.5e1 / 0.81e2 * t23 * t716 * t717 + 0.5e1 / 0.18e2 * t23 * t334 * t159 * t165 - 0.5e1 / 0.12e2 * t23 * t158 * t725 - 0.5e1 / 0.9e1 * t23 * t339 * t346 - 0.5e1 / 0.81e2 * t37 * t734 * t735 + 0.5e1 / 0.18e2 * t37 * t351 * t170 * t174 - 0.5e1 / 0.12e2 * t37 * t169 * t743 - 0.5e1 / 0.9e1 * t37 * t356 * t360 - 0.100e3 / 0.27e2 * t153 * t160 - 0.25e2 / 0.81e2 * t75 * t336 + 0.25e2 / 0.9e1 * t75 * t340;
    t760 = t365 * t100;
    t763 = t89 * t196 * t14;
    t766 = t364 * t56 * t14;
    t769 = t178 * t100 * t14;
    t771 = t299 * t97;
    t773 = t302 * t97;
    t775 = t135 * t193;
    t778 = 0.8e1 / 0.9e1 * t308 * t142;
    t780 = 0.32e2 / 0.81e2 * t139 * t313;
    t783 = 0.40e2 / 0.243e3 * t59 * t61 * t17;
    t787 = 0.16e2 / 0.9e1 * t49 * t384 * t58 * t62;
    t791 = 0.1e1 / t14 / t145;
    t796 = t188 * t188;
    t805 = t93 * (-0.98765432098765432088e0 * t343 * t94 + 0.60307932181100836503e3 * t791 * t189 - 0.11784031256031235593e6 * t147 * t379 + 0.86346472773757568068e7 * t703 / t796 + 0.20218340740740740741e-2 * t53 * t69) * t103 / 0.3e1;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[0] = -0.12930540000000000000e0 * t682 - 0.86203600000000000000e-1 * t684 - 0.21550900000000000000e-1 * t66 * t756 * t56 - 0.86203600000000000000e-1 * t760 - 0.34481440000000000000e0 * t763 - 0.11493813333333333333e0 * t766 - 0.34481440000000000000e0 * t769 - 0.57469066666666666666e-1 * t771 - 0.11493813333333333333e0 * t773 + 0.25541807407407407407e-1 * t775 - t778 + t780 - t783 - t787 - t805;
  
    t806 = t122 * t384;
    t808 = t444 * t100;
    t810 = t225 * t196;
    t811 = 0.64652700000000000000e-1 * t810;
    t812 = t391 * t100;
    t815 = t14 * t443 * t56;
    t817 = t394 * t100;
    t818 = 0.17240720000000000000e0 * t817;
    t819 = t201 * t196;
    t850 = 0.5e1 / 0.2e1 * t23 * t210 * t703;
    t861 = 0.5e1 / 0.2e1 * t37 * t220 * t703;
    t864 = -0.5e1 / 0.4e1 * t37 * t41 * t343 + 0.25e2 / 0.72e2 * t114 * t361 + 0.275e3 / 0.81e2 * t324 * t108 + 0.25e2 / 0.24e2 * t75 * t416 + 0.5e1 / 0.4e1 * t23 * t30 * t343 - 0.25e2 / 0.324e3 * t114 * t353 - 0.5e1 / 0.36e2 * t37 * t216 * t360 - 0.25e2 / 0.8e1 * t75 * t419 - 0.5e1 / 0.81e2 * t23 * t716 * t107 * t335 - 0.5e1 / 0.12e2 * t23 * t158 * t163 * t81 - t850 - 0.5e1 / 0.81e2 * t37 * t734 * t117 * t352 + 0.5e1 / 0.12e2 * t37 * t169 * t163 * t85 + t861 - 0.25e2 / 0.108e3 * t75 * t406;
    t882 = t23 * t334;
    t883 = t107 * t81;
    t893 = t37 * t351;
    t894 = t117 * t85;
    t904 = t75 * t158;
    t907 = -0.50e2 / 0.27e2 * t153 * t207 - 0.25e2 / 0.9e1 * t153 * t211 + 0.25e2 / 0.36e2 * t75 * t413 - 0.5e1 / 0.36e2 * t23 * t206 * t346 + 0.25e2 / 0.36e2 * t114 * t357 + t690 + 0.5e1 / 0.4e1 * t409 * t344 * t81 - 0.5e1 / 0.4e1 * t430 * t344 * t85 + 0.5e1 / 0.36e2 * t882 * t883 * t165 + 0.5e1 / 0.36e2 * t882 * t164 * t159 - 0.5e1 / 0.12e2 * t409 * t164 * t165 + 0.5e1 / 0.36e2 * t893 * t894 * t174 - 0.5e1 / 0.36e2 * t893 * t164 * t170 + 0.5e1 / 0.12e2 * t430 * t164 * t174 + 0.25e2 / 0.18e2 * t904 * t410;
    t913 = t97 * t224 * t56;
    t914 = 0.28734533333333333333e-1 * t913;
    t918 = -0.21550900000000000000e-1 * t806 - 0.64652700000000000000e-1 * t808 - t811 - 0.28734533333333333333e-1 * t812 - 0.86203599999999999999e-1 * t815 - t818 - 0.86203599999999999999e-1 * t819 - 0.21550900000000000000e-1 * t66 * (t864 + t907) * t56 - t914 - 0.64652700000000000000e-1 * t682 - 0.64652700000000000000e-1 * t684 - 0.21550900000000000000e-1 * t760;
    t923 = t193 * t121 * t56;
    t928 = -0.25861080000000000000e0 * t763 - 0.28734533333333333333e-1 * t766 - 0.17240719999999999999e0 * t769 + 0.63854518518518518518e-2 * t923 - 0.28734533333333333334e-1 * t771 - 0.86203599999999999999e-1 * t773 + 0.19156355555555555556e-1 * t775 - t778 + t780 - t783 - t787 - t805;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[1] = t918 + t928;
  
    t975 = -0.5e1 / 0.27e2 * t893 * t894 * t164 + 0.5e1 / 0.27e2 * t882 * t883 * t164 + 0.25e2 / 0.36e2 * t75 * t478 - 0.25e2 / 0.27e2 * t243 * t175 - 0.50e2 / 0.81e2 * t153 * t233 - 0.25e2 / 0.27e2 * t153 * t237 - 0.50e2 / 0.81e2 * t243 * t171 - 0.25e2 / 0.36e2 * t114 * t437 + 0.5e1 / 0.108e3 * t37 * t351 * t252 * t170 - 0.5e1 / 0.81e2 * t23 * t716 * t232 * t159 - 0.5e1 / 0.18e2 * t23 * t206 * t163 + 0.5e1 / 0.108e3 * t23 * t334 * t236 * t159 - 0.25e2 / 0.162e3 * t114 * t427 + 0.25e2 / 0.12e2 * t114 * t440 - 0.5e1 / 0.81e2 * t37 * t734 * t248 * t170 + 0.5e1 / 0.18e2 * t37 * t216 * t163 - 0.25e2 / 0.162e3 * t75 * t463;
    t993 = t24 * t24;
    t995 = 0.1e1 / t686;
    t1015 = t114 * t169;
    t1021 = 0.5e1 / 0.108e3 * t23 * t462 * t165 + 0.25e2 / 0.54e2 * t75 * t473 - 0.5e1 / 0.18e2 * t23 * t158 * t477 * t81 - 0.5e1 / 0.36e2 * t23 * t472 * t165 + 0.25e2 / 0.54e2 * t114 * t434 + 0.5e1 / 0.108e3 * t37 * t487 * t174 - 0.5e1 / 0.18e2 * t37 * t169 * t993 * t995 - 0.5e1 / 0.18e2 * t37 * t169 * t499 * t85 - 0.5e1 / 0.36e2 * t37 * t495 * t174 - 0.5e1 / 0.18e2 * t23 * t158 * t993 * t995 - t850 + t861 + t690 + 0.5e1 / 0.6e1 * t409 * t466 * t343 + 0.25e2 / 0.27e2 * t904 * t467 - 0.25e2 / 0.27e2 * t1015 * t431 - 0.5e1 / 0.6e1 * t430 * t491 * t343;
    t1026 = t504 * t100;
    t1028 = t257 * t196;
    t1031 = t14 * t503 * t56;
    t1033 = t457 * t100;
    t1036 = t97 * t256 * t56;
    t1045 = -0.21550900000000000000e-1 * t66 * (t975 + t1021) * t56 - 0.43101800000000000000e-1 * t1026 - 0.21550900000000000000e-1 * t1028 - 0.57469066666666666666e-1 * t1031 - 0.57469066666666666666e-1 * t1033 - 0.95781777777777777777e-2 * t1036 + 0.12770903703703703703e-1 * t923 - 0.38312711111111111110e-1 * t913 - 0.57469066666666666666e-1 * t815 - 0.22987626666666666667e0 * t817 - 0.17240720000000000000e0 * t819 - 0.43101800000000000000e-1 * t808 - 0.86203600000000000000e-1 * t810;
    t1055 = -0.57469066666666666665e-1 * t812 - 0.43101800000000000000e-1 * t806 - 0.57469066666666666667e-1 * t769 - 0.95781777777777777778e-2 * t771 - 0.57469066666666666667e-1 * t773 + 0.12770903703703703704e-1 * t775 - 0.17240720000000000000e0 * t763 - 0.21550900000000000000e-1 * t682 - 0.43101800000000000000e-1 * t684 - t778 + t780 - t783 - t787 - t805;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[2] = t1045 + t1055;
  
    t1057 = t14 * t549 * t56;
    t1060 = t550 * t100;
    t1064 = 0.6e1 * t343 + 0.12e2 * t704;
    t1097 = -0.5e1 / 0.24e2 * t23 * t30 * t1064 + 0.275e3 / 0.81e2 * t529 * t86 + 0.25e2 / 0.24e2 * t114 * t500 + 0.5e1 / 0.24e2 * t37 * t41 * t1064 - 0.25e2 / 0.324e3 * t75 * t515 + 0.25e2 / 0.72e2 * t75 * t523 + 0.25e2 / 0.36e2 * t75 * t518 - 0.5e1 / 0.81e2 * t23 * t716 * t514 * t81 - 0.5e1 / 0.12e2 * t23 * t206 * t477 - 0.5e1 / 0.36e2 * t23 * t158 * t522 * t81 - 0.50e2 / 0.27e2 * t243 * t217 + 0.25e2 / 0.9e1 * t243 * t221;
    t1137 = -0.25e2 / 0.108e3 * t114 * t488 + 0.25e2 / 0.36e2 * t114 * t496 - 0.5e1 / 0.81e2 * t37 * t734 * t538 * t85 - 0.5e1 / 0.12e2 * t37 * t216 * t499 - 0.5e1 / 0.36e2 * t37 * t169 * t545 * t85 + t690 + 0.5e1 / 0.36e2 * t882 * t232 * t24 * t163 + 0.5e1 / 0.36e2 * t882 * t107 * t236 * t81 - 0.5e1 / 0.12e2 * t409 * t164 * t236 - 0.25e2 / 0.18e2 * t1015 * t492 - 0.5e1 / 0.36e2 * t893 * t248 * t24 * t163 + 0.5e1 / 0.36e2 * t893 * t117 * t252 * t85 + 0.5e1 / 0.12e2 * t430 * t164 * t252;
    t1148 = -0.28734533333333333333e-1 * t1057 - 0.28734533333333333333e-1 * t1036 - 0.21550900000000000000e-1 * t1060 - 0.21550900000000000000e-1 * t66 * (t1097 + t1137) * t56 - 0.64652700000000000000e-1 * t1026 - 0.64652700000000000000e-1 * t1028 - 0.86203599999999999999e-1 * t1031 - 0.17240720000000000000e0 * t1033 - 0.64652700000000000000e-1 * t806 - t811 - 0.86203600000000000000e-1 * t812 - t818;
    t1155 = -0.25861080000000000000e0 * t819 - t914 - 0.21550900000000000000e-1 * t684 - 0.86203600000000000000e-1 * t763 + 0.19156355555555555555e-1 * t923 - 0.28734533333333333333e-1 * t773 + 0.63854518518518518519e-2 * t775 - t778 + t780 - t783 - t787 - t805;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[3] = t1148 + t1155;
  
    t1159 = -0.12e2 * t343 - 0.12e2 * t704;
    t1177 = t232 * t232;
    t1184 = t236 * t236;
    t1197 = t248 * t248;
    t1204 = t252 * t252;
    t1211 = 0.25e2 / 0.18e2 * t114 * t546 - 0.5e1 / 0.24e2 * t37 * t41 * t1159 + 0.5e1 / 0.24e2 * t23 * t30 * t1159 - 0.770e3 / 0.81e2 * lapl[1] / t34 / t526 / rho[1] * t42 + 0.1100e4 / 0.81e2 * t529 * t118 - 0.50e2 / 0.9e1 * t243 * t253 - 0.5e1 / 0.81e2 * t23 * t716 * t1177 + 0.5e1 / 0.18e2 * t23 * t462 * t236 - 0.5e1 / 0.12e2 * t23 * t158 * t1184 - 0.5e1 / 0.9e1 * t23 * t206 * t522 - 0.100e3 / 0.27e2 * t243 * t249 - 0.25e2 / 0.81e2 * t114 * t539 + 0.25e2 / 0.9e1 * t114 * t542 - 0.5e1 / 0.81e2 * t37 * t734 * t1197 + 0.5e1 / 0.18e2 * t37 * t487 * t252 - 0.5e1 / 0.12e2 * t37 * t169 * t1204 - 0.5e1 / 0.9e1 * t37 * t216 * t545 + t690;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho4[4] = -0.21550900000000000000e-1 * t66 * t1211 * t56 - 0.11493813333333333333e0 * t1057 - 0.57469066666666666666e-1 * t1036 - 0.86203600000000000000e-1 * t1060 - 0.12930540000000000000e0 * t1028 - 0.34481440000000000000e0 * t1033 - 0.86203600000000000000e-1 * t806 - 0.11493813333333333333e0 * t812 - 0.34481440000000000000e0 * t819 + 0.25541807407407407407e-1 * t923 - t778 + t780 - t783 - t787 - t805;
  
    t1224 = t791 * t56;
    t1226 = t375 * t100;
    t1228 = t187 * t196;
    t1230 = t125 * t384;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[0] = 0.27936351851851851852e-1 * t1224 - 0.25142716666666666667e-1 * t1226 + 0.10775450000000000000e-1 * t1228 - 0.26938625000000000000e-2 * t1230;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[1] = 0.55872703703703703703e-1 * t1224 - 0.50285433333333333334e-1 * t1226 + 0.21550900000000000000e-1 * t1228 - 0.53877250000000000000e-2 * t1230;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[2] = v4rho3sigma[0];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[3] = v4rho3sigma[2];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[4] = v4rho3sigma[1];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[5] = v4rho3sigma[3];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[6] = v4rho3sigma[5];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[7] = v4rho3sigma[4];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[8] = v4rho3sigma[6];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[9] = v4rho3sigma[8];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[10] = v4rho3sigma[7];
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3sigma[11] = v4rho3sigma[9];
  
    t1237 = 0.79818148148148148147e-3 * t630 * t130;
    t1239 = t97 * t74 * t130;
    t1242 = t14 * t152 * t130;
    t1249 = 0.26938625000000000000e-2 * t129 * t31 * t384;
    t1251 = 0.35918166666666666667e-2 * t564 * t277;
    t1252 = t567 * t277;
    t1255 = 0.10775450000000000000e-1 * t267 * t592;
    t1256 = t574 * t277;
    t1258 = t270 * t592;
    t1260 = t564 * t274;
    t1262 = t567 * t274;
    t1264 = -t1237 - 0.59863611111111111110e-2 * t1239 + 0.47890888888888888888e-1 * t1242 - 0.43899981481481481481e-1 * t66 * t323 * t130 + t1249 + t1251 - 0.35918166666666666666e-1 * t1252 + t1255 + 0.35918166666666666666e-1 * t1256 - 0.13469312500000000000e-1 * t1258 + 0.59863611111111111110e-2 * t1260 - 0.59863611111111111111e-1 * t1262;
    t1265 = t267 * t586;
    t1267 = t267 * t589;
    t1271 = t270 * t586;
    t1275 = t158 * t100;
    t1277 = t129 * t1275 * t159;
    t1279 = t30 * t196;
    t1281 = t129 * t1279 * t81;
    t1284 = t129 * t585 * t165;
    t1289 = t267 * t582;
    t1297 = t334 * t56;
    t1301 = 0.35918166666666666666e-1 * t1265 + 0.17959083333333333333e-1 * t1267 + 0.59863611111111111112e-1 * t574 * t274 - 0.44897708333333333334e-1 * t1271 - 0.22448854166666666667e-1 * t270 * t589 + 0.89795416666666666666e-2 * t1277 + 0.13469312500000000000e-1 * t1281 + 0.13469312500000000000e-1 * t1284 + 0.44897708333333333333e-2 * t129 * t273 * t346 + 0.11972722222222222222e-1 * t1289 + 0.89795416666666666665e-2 * t622 * t56 * t81 * t165 - 0.14965902777777777778e-1 * t270 * t582 - 0.99772685185185185183e-3 * t129 * t1297 * t335;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[0] = t1264 + t1301;
  
    t1306 = 0.79818148148148148147e-3 * t654 * t133;
    t1307 = t595 * t284;
    t1309 = t280 * t607;
    t1311 = t280 * t610;
    t1313 = t169 * t100;
    t1315 = t132 * t1313 * t170;
    t1317 = t41 * t196;
    t1319 = t132 * t1317 * t85;
    t1322 = t132 * t606 * t174;
    t1325 = 0.35918166666666666667e-2 * t595 * t287;
    t1327 = 0.10775450000000000000e-1 * t280 * t613;
    t1330 = 0.26938625000000000000e-2 * t132 * t42 * t384;
    t1331 = t351 * t56;
    t1339 = t280 * t603;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[1] = 0.44897708333333333333e-2 * t132 * t283 * t360 - t1306 + 0.59863611111111111110e-2 * t1307 + 0.35918166666666666666e-1 * t1309 + 0.17959083333333333333e-1 * t1311 + 0.89795416666666666666e-2 * t1315 + 0.13469312500000000000e-1 * t1319 + 0.13469312500000000000e-1 * t1322 + t1325 + t1327 + t1330 - 0.99772685185185185183e-3 * t132 * t1331 * t352 + 0.89795416666666666665e-2 * t646 * t56 * t85 * t174 + 0.11972722222222222222e-1 * t1339;
  
    t1355 = 0.44897708333333333333e-2 * t630 * t273 - t1237 - 0.39909074074074074074e-2 * t1239 + 0.15963629629629629630e-1 * t1242 + t1249 + t1251 - 0.23945444444444444444e-1 * t1252 + t1255 + 0.11972722222222222222e-1 * t1256 - 0.89795416666666666666e-2 * t1258 + 0.39909074074074074074e-2 * t1260 - 0.19954537037037037037e-1 * t1262 + 0.23945444444444444444e-1 * t1265 + 0.59863611111111111112e-2 * t1267 - 0.14965902777777777778e-1 * t1271 + 0.29931805555555555555e-2 * t1277 + 0.89795416666666666666e-2 * t1281;
    t1360 = t622 * t100 * t107 * t81;
    t1365 = t129 * t334;
    t1369 = t630 * t158;
    t1370 = t56 * t24;
    t1376 = 0.79818148148148148147e-2 * t267 * t158 * t624;
    t1380 = t267 * t627;
    t1383 = 0.14965902777777777778e-1 * t270 * t627;
    t1385 = t129 * t1279 * t107;
    t1389 = 0.89795416666666666666e-2 * t630 * t585 * t24;
    t1390 = t564 * t290;
    t1393 = 0.19954537037037037037e-1 * t567 * t290;
    t1395 = t17 * t22 * t631;
    t1402 = 0.44897708333333333333e-2 * t1284 + 0.39909074074074074075e-2 * t1289 + 0.59863611111111111110e-2 * t1360 + 0.29931805555555555555e-2 * t622 * t623 * t165 - 0.99772685185185185183e-3 * t1365 * t623 * t159 + 0.59863611111111111110e-2 * t1369 * t1370 * t81 + t1376 - 0.99772685185185185184e-2 * t270 * t158 * t624 + 0.11972722222222222222e-1 * t1380 - t1383 + 0.44897708333333333333e-2 * t1385 + t1389 + 0.19954537037037037037e-2 * t1390 - t1393 - 0.14965902777777777777e-2 * t1395 + 0.19954537037037037037e-1 * t574 * t290 - 0.14965902777777777778e-1 * t193 * t74 * t631;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[2] = t1355 + t1402;
  
    t1406 = t97 * t113 * t133;
    t1408 = t637 * t287;
    t1410 = t293 * t613;
    t1412 = t132 * t351;
    t1416 = t654 * t169;
    t1422 = 0.79818148148148148147e-2 * t280 * t169 * t648;
    t1425 = t646 * t100 * t117 * t85;
    t1431 = t1330 - 0.44897708333333333333e-2 * t654 * t283 - 0.19954537037037037037e-2 * t1406 - 0.11972722222222222222e-1 * t1408 - 0.44897708333333333333e-2 * t1410 - t1306 + t1325 + t1327 - 0.99772685185185185183e-3 * t1412 * t647 * t170 - 0.59863611111111111110e-2 * t1416 * t1370 * t85 + t1422 + 0.59863611111111111110e-2 * t1425 + 0.29931805555555555555e-2 * t646 * t647 * t174 + 0.39909074074074074074e-2 * t1307;
    t1439 = 0.19954537037037037037e-1 * t637 * t284;
    t1441 = 0.14965902777777777778e-1 * t293 * t607;
    t1444 = t280 * t651;
    t1447 = t17 * t36 * t655;
    t1450 = t132 * t1317 * t117;
    t1454 = 0.89795416666666666666e-2 * t654 * t606 * t24;
    t1457 = t595 * t296;
    t1459 = 0.23945444444444444444e-1 * t1309 + 0.59863611111111111112e-2 * t1311 + 0.29931805555555555555e-2 * t1315 + 0.89795416666666666666e-2 * t1319 + 0.44897708333333333333e-2 * t1322 + 0.39909074074074074075e-2 * t1339 - t1439 - t1441 - 0.74829513888888888888e-2 * t293 * t610 + 0.11972722222222222222e-1 * t1444 + 0.14965902777777777777e-2 * t1447 + 0.44897708333333333333e-2 * t1450 - t1454 - 0.49886342592592592592e-2 * t293 * t603 + 0.19954537037037037037e-2 * t1457;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[3] = t1431 + t1459;
  
    t1478 = -t1237 - 0.19954537037037037037e-2 * t1239 + t1249 + t1251 - 0.11972722222222222222e-1 * t1252 + t1255 - 0.44897708333333333333e-2 * t1258 + 0.19954537037037037037e-2 * t1260 + 0.11972722222222222222e-1 * t1265 + 0.44897708333333333333e-2 * t1281 + 0.59863611111111111111e-2 * t1360 + 0.59863611111111111110e-2 * t1369 * t623 * t24 + 0.29931805555555555555e-2 * t622 * t56 * t236 * t81 - 0.99772685185185185183e-3 * t1365 * t56 * t232 * t81;
    t1484 = t129 * t1275 * t232;
    t1487 = t129 * t585 * t236;
    t1492 = t267 * t660;
    t1496 = t267 * t664;
    t1500 = t1376 + 0.23945444444444444444e-1 * t1380 - t1383 + 0.89795416666666666666e-2 * t1385 + t1389 + 0.39909074074074074073e-2 * t1390 - t1393 + 0.11972722222222222222e-1 * t1395 + 0.29931805555555555555e-2 * t1484 + 0.44897708333333333333e-2 * t1487 + 0.44897708333333333333e-2 * t129 * t273 * t477 + 0.39909074074074074073e-2 * t1492 - 0.49886342592592592592e-2 * t270 * t660 + 0.59863611111111111111e-2 * t1496 - 0.74829513888888888888e-2 * t270 * t664;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[4] = t1478 + t1500;
  
    t1505 = t14 * t242 * t133;
    t1507 = t669 * t287;
    t1526 = t1330 - 0.39909074074074074073e-2 * t1406 - 0.23945444444444444444e-1 * t1408 - 0.89795416666666666666e-2 * t1410 + 0.15963629629629629629e-1 * t1505 + 0.11972722222222222222e-1 * t1507 - t1306 + t1325 + t1327 - 0.99772685185185185187e-2 * t293 * t169 * t648 - 0.99772685185185185183e-3 * t1412 * t56 * t248 * t85 - 0.59863611111111111110e-2 * t1416 * t647 * t24 + 0.29931805555555555555e-2 * t646 * t56 * t252 * t85 + t1422 + 0.59863611111111111111e-2 * t1425 + 0.19954537037037037037e-2 * t1307 + 0.11972722222222222222e-1 * t1309;
    t1534 = t293 * t651;
    t1540 = t132 * t1313 * t248;
    t1543 = t132 * t606 * t252;
    t1548 = t637 * t296;
    t1550 = t280 * t675;
    t1552 = t280 * t679;
    t1554 = 0.44897708333333333333e-2 * t1319 - t1439 - t1441 + 0.23945444444444444444e-1 * t1444 - 0.11972722222222222222e-1 * t1447 + 0.89795416666666666666e-2 * t1450 - t1454 + 0.39909074074074074073e-2 * t1457 + 0.19954537037037037037e-1 * t669 * t284 - 0.14965902777777777778e-1 * t1534 + 0.14965902777777777778e-1 * t193 * t113 * t655 + 0.29931805555555555555e-2 * t1540 + 0.44897708333333333333e-2 * t1543 + 0.44897708333333333333e-2 * t132 * t283 * t499 - 0.19954537037037037037e-1 * t1548 + 0.39909074074074074073e-2 * t1550 + 0.59863611111111111111e-2 * t1552;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[5] = t1526 + t1554;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[6] = -0.99772685185185185183e-3 * t129 * t1297 * t514 + 0.89795416666666666665e-2 * t622 * t623 * t236 + 0.44897708333333333333e-2 * t129 * t273 * t522 + 0.89795416666666666666e-2 * t1484 + 0.13469312500000000000e-1 * t1487 + 0.35918166666666666666e-1 * t1380 + 0.13469312500000000000e-1 * t1385 - t1237 + t1249 + t1251 + t1255 + 0.11972722222222222222e-1 * t1492 + 0.17959083333333333333e-1 * t1496 + 0.59863611111111111110e-2 * t1390;
  
    t1583 = t1330 - 0.59863611111111111110e-2 * t1406 - 0.35918166666666666666e-1 * t1408 - 0.13469312500000000000e-1 * t1410 + 0.47890888888888888888e-1 * t1505 + 0.35918166666666666666e-1 * t1507 - 0.43899981481481481481e-1 * t66 * t528 * t133 - t1306 + t1325 + t1327 + 0.89795416666666666665e-2 * t646 * t647 * t252 + 0.35918166666666666666e-1 * t1444;
    t1604 = 0.13469312500000000000e-1 * t1450 + 0.59863611111111111110e-2 * t1457 - 0.44897708333333333334e-1 * t1534 + 0.89795416666666666666e-2 * t1540 + 0.13469312500000000000e-1 * t1543 - 0.59863611111111111111e-1 * t1548 + 0.11972722222222222222e-1 * t1550 + 0.17959083333333333333e-1 * t1552 - 0.14965902777777777778e-1 * t293 * t675 - 0.99772685185185185183e-3 * t132 * t1331 * t538 - 0.22448854166666666667e-1 * t293 * t679 + 0.44897708333333333333e-2 * t132 * t283 * t545 + 0.59863611111111111112e-1 * t669 * t296;
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3lapl[7] = t1583 + t1604;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho3tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigma2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmalapl[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2sigmatau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapl2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2lapltau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rho2tau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma3[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[20] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[21] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[22] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2lapl[23] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[20] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[21] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[22] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigma2tau[23] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapl2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[20] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[21] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[22] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmalapltau[23] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhosigmatau2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapl2tau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rholapltau2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4rhotau3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma4[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3lapl[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma3tau[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapl2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[18] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[19] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[20] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[21] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[22] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2lapltau[23] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigma2tau2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl3[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapl2tau[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[12] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[13] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[14] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[15] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[16] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmalapltau2[17] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[9] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[10] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4sigmatau3[11] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl4[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl3tau[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapl2tau2[8] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[4] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[5] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[6] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_NEEDS_LAPLACIAN) && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4lapltau3[7] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[0] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[1] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[2] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[3] = 0.0e0;
  
    if(v4rho4 != NULL && (p->info->flags & XC_FLAGS_HAVE_LXC))
      v4tau4[4] = 0.0e0;
  
  #ifndef XC_DONT_COMPILE_MXC
  
    if(order < 5) return;
  
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  #endif
  
  
  }
  
  