From 9671a632b55110992eb6878a61c2488b05abce15 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Sun, 26 Apr 2020 14:58:07 +0200
Subject: [PATCH] compare floats with tolerances

---
 src/rotations.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/rotations.f90 b/src/rotations.f90
index 98c529933..8ff5aca4e 100644
--- a/src/rotations.f90
+++ b/src/rotations.f90
@@ -527,7 +527,7 @@ pure function qu2ho(qu) result(ho)
 
   omega = 2.0 * acos(math_clip(qu(1),-1.0_pReal,1.0_pReal))
 
-  if (dEq0(omega)) then
+  if (dEq0(omega,tol=1.e-5_pReal)) then
     ho = [ 0.0_pReal, 0.0_pReal, 0.0_pReal ]
   else
     ho = qu(2:4)
@@ -1298,7 +1298,7 @@ pure function Lambert_BallToCube(xyz) result(cube)
   integer,                 dimension(3,2) :: p
 
   rs = norm2(xyz)
-  if (rs > R1) then
+  if (rs > R1+1.e-6_pReal) then
     cube = IEEE_value(cube,IEEE_positive_inf)
     return
   endif