From ad851bfede223c8bcd616babd4b4f6f72ba7106c Mon Sep 17 00:00:00 2001
From: tianchiTJ <149181511@qq.com>
Date: Tue, 3 Nov 2020 11:29:47 -0800
Subject: [PATCH] Change the density tolerance

---
 include/node.tcc | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/include/node.tcc b/include/node.tcc
index 7d26cc6d2..e4840f1a0 100644
--- 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.;
     }
   }
 }