From da041f43bc26763c8b1fcaf97ac3acb6abe4c086 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 23 Sep 2019 09:29:38 -0700 Subject: [PATCH] floating point arithmetic refinement Intel 18.4 could not detect 0.0 --- src/prec.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/prec.f90 b/src/prec.f90 index 434ce32e7..99032d8ef 100644 --- a/src/prec.f90 +++ b/src/prec.f90 @@ -127,7 +127,7 @@ logical elemental pure function dEq(a,b,tol) eps = PREAL_EPSILON * maxval(abs([a,b])) endif - dEq = merge(.True.,.False.,abs(a-b) < eps) + dEq = merge(.True.,.False.,abs(a-b) <= eps) end function dEq @@ -170,7 +170,7 @@ logical elemental pure function dEq0(a,tol) eps = PREAL_MIN * 10.0_pReal endif - dEq0 = merge(.True.,.False.,abs(a) < eps) + dEq0 = merge(.True.,.False.,abs(a) <= eps) end function dEq0 @@ -214,7 +214,7 @@ logical elemental pure function cEq(a,b,tol) eps = PREAL_EPSILON * maxval(abs([a,b])) endif - cEq = merge(.True.,.False.,abs(a-b) < eps) + cEq = merge(.True.,.False.,abs(a-b) <= eps) end function cEq