Another numerical edge case added for the virial simulations
Trying to increase the numerical accuracy for virial integration as its sensitive for the constant volume stepping
Fixed a numerical instability in the virial stepping