Changed wrappingLayerEpsilon in lsAdvect.

ViennaTools · Sep 30, 2020 · 6e2fda1 · 6e2fda1
1 parent 3001163
commit 6e2fda1
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/include/lsAdvect.hpp b/include/lsAdvect.hpp
@@ -64,6 +64,7 @@ template <class T, int D> class lsAdvect {
   double advectedTime = 0.;
   unsigned numberOfTimeSteps = 0;
   bool saveAdvectionVelocities = false;
+  static constexpr double wrappingLayerEpsilon = 1e-4;
 
   // SFINAE functions needed for StencilLocalLaxFriedrichs
   template <
@@ -503,7 +504,7 @@ template <class T, int D> class lsAdvect {
 
               // if the lower surface is actually outside, i.e. its LS value is
               // lower or equal
-              if (iterators[lowerLevelSetId].getValue() <= value + 1e-9) {
+              if (iterators[lowerLevelSetId].getValue() <= value + wrappingLayerEpsilon) {
                 velocity =
                     scheme(it.getStartIndices(), velocities, lowerLevelSetId);
                 break;