floating point arithmetic refinement
Intel 18.4 could not detect 0.0
This commit is contained in:
parent
96811f77ca
commit
da041f43bc
|
@ -127,7 +127,7 @@ logical elemental pure function dEq(a,b,tol)
|
||||||
eps = PREAL_EPSILON * maxval(abs([a,b]))
|
eps = PREAL_EPSILON * maxval(abs([a,b]))
|
||||||
endif
|
endif
|
||||||
|
|
||||||
dEq = merge(.True.,.False.,abs(a-b) < eps)
|
dEq = merge(.True.,.False.,abs(a-b) <= eps)
|
||||||
|
|
||||||
end function dEq
|
end function dEq
|
||||||
|
|
||||||
|
@ -170,7 +170,7 @@ logical elemental pure function dEq0(a,tol)
|
||||||
eps = PREAL_MIN * 10.0_pReal
|
eps = PREAL_MIN * 10.0_pReal
|
||||||
endif
|
endif
|
||||||
|
|
||||||
dEq0 = merge(.True.,.False.,abs(a) < eps)
|
dEq0 = merge(.True.,.False.,abs(a) <= eps)
|
||||||
|
|
||||||
end function dEq0
|
end function dEq0
|
||||||
|
|
||||||
|
@ -214,7 +214,7 @@ logical elemental pure function cEq(a,b,tol)
|
||||||
eps = PREAL_EPSILON * maxval(abs([a,b]))
|
eps = PREAL_EPSILON * maxval(abs([a,b]))
|
||||||
endif
|
endif
|
||||||
|
|
||||||
cEq = merge(.True.,.False.,abs(a-b) < eps)
|
cEq = merge(.True.,.False.,abs(a-b) <= eps)
|
||||||
|
|
||||||
end function cEq
|
end function cEq
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue