Use numerical limits.

erf-model · Nov 6, 2024 · a0aa615 · a0aa615
1 parent 48e7860
commit a0aa615
Show file tree

Hide file tree

Showing 3 changed files with 3 additions and 3 deletions.
diff --git a/Source/Microphysics/Kessler/ERF_Kessler.cpp b/Source/Microphysics/Kessler/ERF_Kessler.cpp
@@ -177,7 +177,7 @@ void Kessler::AdvanceKessler (const SolverChoice &solverChoice)
                                        {
                                            return { ma_fz_arr[box_no](i,j,k) };
                                        });
-        wt_max = get<0>(max) + 1.0e-16;
+        wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
         n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
         AMREX_ALWAYS_ASSERT(n_substep >= 1);
         coef /= Real(n_substep);

diff --git a/Source/Microphysics/SAM/ERF_IceFall.cpp b/Source/Microphysics/SAM/ERF_IceFall.cpp
@@ -80,7 +80,7 @@ void SAM::IceFall (const SolverChoice& sc) {
                          {
                              return { ma_fz_arr[box_no](i,j,k) };
                          });
-    wt_max = get<0>(max) + 1.0e-16;
+    wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
     n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
     AMREX_ALWAYS_ASSERT(n_substep >= 1);
     coef /= Real(n_substep);

diff --git a/Source/Microphysics/SAM/ERF_PrecipFall.cpp b/Source/Microphysics/SAM/ERF_PrecipFall.cpp
@@ -124,7 +124,7 @@ SAM::PrecipFall (const SolverChoice& sc)
                          {
                              return { ma_fz_arr[box_no](i,j,k) };
                          });
-    wt_max = get<0>(max) + 1.0e-16;
+    wt_max = get<0>(max) + std::numeric_limits<Real>::epsilon();
     n_substep = int( std::ceil(wt_max * coef / CFL_MAX) );
     AMREX_ALWAYS_ASSERT(n_substep >= 1);
     coef /= Real(n_substep);