Change the density tolerance

cb-geo · Nov 3, 2020 · ad851bf · ad851bf
1 parent 1b34a27
commit ad851bf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/include/node.tcc b/include/node.tcc
@@ -595,7 +595,7 @@ void mpm::Node<Tdim, Tdof, Tnphases>::compute_density() {
         density_(phase) = mass_(phase) / volume_(phase);
 
       // Check to see if value is below threshold
-      if (std::abs(density_(phase)) < 1.E-15) density_(phase) = 0.;
+      if (std::abs(density_(phase)) < tolerance) density_(phase) = 0.;
     }
   }
 }