From c24e46c547db11342d3f4d20a3d5fafe54411cc9 Mon Sep 17 00:00:00 2001 From: Axel Kohlmeyer Date: Sat, 5 Mar 2022 12:30:40 -0500 Subject: [PATCH] don't apply modulus to doubles but integers --- src/EXTRA-COMPUTE/compute_pressure_cartesian.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/EXTRA-COMPUTE/compute_pressure_cartesian.cpp b/src/EXTRA-COMPUTE/compute_pressure_cartesian.cpp index 999eeba254..3df0e82e3c 100644 --- a/src/EXTRA-COMPUTE/compute_pressure_cartesian.cpp +++ b/src/EXTRA-COMPUTE/compute_pressure_cartesian.cpp @@ -372,13 +372,13 @@ void ComputePressureCartesian::compute_pressure_1d(double fpair, double xi, doub xj -= (domain->boxhi[dir1] - domain->boxlo[dir1]); // Integrating contour from bin_s to bin_e - bin_s = (int) lround(((xi - domain->boxlo[dir1]) / bin_width1)) % nbins1; - bin_e = (int) lround(((xj - domain->boxlo[dir1]) / bin_width1)) % nbins1; + bin_s = ((int) lround((xi - domain->boxlo[dir1]) / bin_width1)) % nbins1; + bin_e = ((int) lround((xj - domain->boxlo[dir1]) / bin_width1)) % nbins1; // If not periodic in dir1 if (domain->periodicity[dir1] == 0) { - bin_s = (int) lround(((xi - domain->boxlo[dir1]) / bin_width1)); - bin_e = (int) lround(((xj - domain->boxlo[dir1]) / bin_width1)); + bin_s = ((int) lround((xi - domain->boxlo[dir1]) / bin_width1)); + bin_e = ((int) lround((xj - domain->boxlo[dir1]) / bin_width1)); if (bin_e == nbins1) bin_e--; if (bin_s == nbins1) bin_s--;