Bug fix in sr_smonotonic: avoid a floating point rounding error

edf-hpc · Oct 3, 2024 · 9340766 · 9340766
1 parent 7080042
commit 9340766
Showing 1 changed file with 4 additions and 3 deletions.
diff --git a/interflop_backends/interflop_verrou/vr_roundingOp.hxx b/interflop_backends/interflop_verrou/vr_roundingOp.hxx
@@ -421,10 +421,11 @@ public:
 	}else{
 	  resm=nextAwayFromZero<RealType>(res);
 	  const RealType resHash(-resp);
-	  const RealType limit=1.-RAND::randRatioFromResult(&vr_rand, &resHash);
+	  //	  const RealType limit=1.-RAND::randRatioFromResult(&vr_rand, &resHash);
+	  const RealType x=RAND::randRatioFromResult(&vr_rand, &resHash);
 	  const RealType u(resp-resm);
-	  const RealType errorTh(u+error);
-	  down =( errorTh <= (limit * u));
+	  //	  const RealType errorTh(u+error);
+	  down =( error <= (-x) * u);
 	}
 
       }