577     real delta, ep_p[6], ep_f[6], ep_i[6], d_ef0, d_ei0;  in vwn_fourth()  local
634     d_ef0 = ep_f[5] - ep_p[5];  in vwn_fourth()
714     d_ef2bi = ef3 + (-ef0 + ef1)/rho2 - (-d_ef0 + ef2)/rho;  in vwn_fourth()
722     d_delta_indep = 4*(d_ef0 - d_ei0)*f_zeta*zeta3 +  in vwn_fourth()
723         d_ei0*f_zet1*(1 - zeta4) + d_ef0*f_zet1*zeta4;  in vwn_fourth()
753     d_ctrm2_indep = d_bterm_indep*f_zet2 + (-d_ef0 + d_ei0 + ef2 - ei2)*  in vwn_fourth()
755         d_ei0*f_zet2*(-1 + zeta4) - d_ef0*f_zet2*zeta4;  in vwn_fourth()
772     d_dterm_indep = 12*(ef0 - ei0)*f_zeta*zeta2 + 12*d_ef0*f_zeta*rho*  in vwn_fourth()
774         8*d_ef0*f_zet1*rho*zeta3 + d_ef0*f_zet2*rho*zeta4 +  in vwn_fourth()
787         (-d_ef0 + d_ei0 + ef2 - ei2)*(12*f_zeta*zeta2 +  in vwn_fourth()
788                                       8*f_zet1*zeta3) - d_ef0*f_zet2*zeta4 +  in vwn_fourth()
802         12*(d_ef0 - d_ei0)*rho*(2*f_zeta*zeta + 3*f_zet1*zeta2 +  in vwn_fourth()
803                                 f_zet2*zeta3) + d_ef0*f_zet3*rho*zeta4 -  in vwn_fourth()