almostEqual for arbitrary floating point types

ajnonaka · Sep 16, 2024 · 1ede927 · 1ede927
1 parent 0c79af8
commit 1ede927
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/Src/Base/AMReX_RKIntegrator.H b/Src/Base/AMReX_RKIntegrator.H
@@ -252,7 +252,7 @@ public:
         {
             // Adjust step size to reach output time
             // protect against roundoff
-            if (std::abs(time_out - time_current - dt) <= 1.e5 * std::numeric_limits<amrex::Real>::epsilon() * std::abs(dt)) {
+            if (almostEqual(std::abs(time_out-time_current), std::abs(dt)), 10.) {
                 dt = time_out - time_current;
                 stop = true;
             }