From e069dcb119317e27a3237e0408207d79f30237fe Mon Sep 17 00:00:00 2001 From: Luke Shingles Date: Fri, 7 Mar 2025 21:37:46 +0000 Subject: [PATCH] Update grid.cc --- grid.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/grid.cc b/grid.cc index 1f6649b14..0292429a0 100644 --- a/grid.cc +++ b/grid.cc @@ -2514,7 +2514,7 @@ auto get_totmassradionuclide(const int z, const int a) -> double { std::tie(d_coordminboundary, d_coordmaxboundary) = get_coordboundary_distances_cylindrical2d( pos, dir, pktposgridcoord, pktvelgridcoord, cellindex, tstart, cellcoordmax); - if (pktvelgridcoord[1] > (pktposgridcoord[1] / tstart)) { + if ((pktvelgridcoord[1] * tstart) > pktposgridcoord[1]) { d_coordminboundary[1] = -1; } else { d_coordmaxboundary[1] = -1; @@ -2534,7 +2534,7 @@ auto get_totmassradionuclide(const int z, const int a) -> double { // boundary, regardless of direction. for (int d = 0; d < 3; d++) { - if (pktvelgridcoord[d] > (pktposgridcoord[d] / tstart)) { + if ((pktvelgridcoord[d] * tstart) > pktposgridcoord[d]) { d_coordminboundary[d] = -1; const double t_coordmaxboundary = ((pktposgridcoord[d] - (pktvelgridcoord[d] * tstart)) / (cellcoordmax[d] - (pktvelgridcoord[d] * globals::tmin)) * globals::tmin) -