ir->nkx = ir->nky = ir->nkz = 0;
calc_grid(NULL, state.box, fourierspacing*fac, &ir->nkx, &ir->nky, &ir->nkz);
- /* Adjust other radii since various conditions neet to be fulfilled */
+ /* Adjust other radii since various conditions need to be fulfilled */
if (eelPME == ir->coulombtype)
{
/* plain PME, rcoulomb must be equal to rlist */
if (bScaleRvdw && evdwCUT == ir->vdwtype)
{
- /* For vdw cutoff, rvdw >= rlist */
- ir->rvdw = max(info->rvdw[0], ir->rlist);
+ if ( ecutsVERLET == ir->cutoff_scheme)
+ {
+ /* With Verlet, the van der Waals radius must always equal the Coulomb radius */
+ ir->rvdw = ir->rcoulomb;
+ }
+ else
+ {
+ /* For vdw cutoff, rvdw >= rlist */
+ ir->rvdw = max(info->rvdw[0], ir->rlist);
+ }
}
ir->rlistlong = max_cutoff(ir->rlist, max_cutoff(ir->rvdw, ir->rcoulomb));