+ if (EVDW_PME(ic->vdwtype))
+ {
+ /* We have PME for both Coulomb and VdW, set rvdw equal to rcoulomb */
+ ic->rvdw = set->rcut_coulomb;
+ ic->ewaldcoeff_lj = set->ewaldcoeff_lj;
+ if (ic->vdw_modifier == eintmodPOTSHIFT)
+ {
+ real crc2;
+
+ ic->dispersion_shift.cpot = -pow(ic->rvdw, -6.0);
+ ic->repulsion_shift.cpot = -pow(ic->rvdw, -12.0);
+ ic->sh_invrc6 = -ic->dispersion_shift.cpot;
+ crc2 = sqr(ic->ewaldcoeff_lj*ic->rvdw);
+ ic->sh_lj_ewald = (exp(-crc2)*(1 + crc2 + 0.5*crc2*crc2) - 1)*pow(ic->rvdw, -6.0);
+ }
+ }