fixed poisson solver, thus virial errors