From 671c86c079faa37386f59eb3c814dc837c4df4c5 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Tue, 2 Apr 2013 18:29:33 -0400 Subject: Replace abs() by fabs() Replace calls to abs() with floating arguments by calls to fabs(), to ensure the results are not accidentally truncated to integer. --- Carpet/CarpetInterp2/src/fasterp.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'Carpet/CarpetInterp2/src/fasterp.cc') diff --git a/Carpet/CarpetInterp2/src/fasterp.cc b/Carpet/CarpetInterp2/src/fasterp.cc index fc3ac8a90..3230ba751 100644 --- a/Carpet/CarpetInterp2/src/fasterp.cc +++ b/Carpet/CarpetInterp2/src/fasterp.cc @@ -1190,12 +1190,12 @@ namespace CarpetInterp2 { assert (all (ipos % ext.stride() == ivect(0))); ipos /= ext.stride(); dpos /= rvect(ext.stride()); - if (not (all (abs(dpos) <= rvect(0.5)))) { + if (not (all (fabs(dpos) <= rvect(0.5)))) { cout << "fasterp.cc:659\n" << " dpos=" << dpos << "\n" << " ext=" << ext << "\n"; } - assert (all (abs(dpos) <= rvect(0.5))); + assert (all (fabs(dpos) <= rvect(0.5))); ivect const ind = ipos - ext.lower() / ext.stride(); ivect const ash = pad_shape(ext); -- cgit v1.2.3