I have the impression that fix setforce operates after the virial is computed in the pair routine.
So, in spite of the fix setforce 0 the group to which it applies still contributes to the pressure.
Is this right?
Is there any way -besides setting the pair interaction to zero- to make the contribution of a group
to the pressure zero in a meaningful way?
thanks for your help,