diff --git a/srcs/FEM/solverFEM.cpp b/srcs/FEM/solverFEM.cpp
index 97bcda367a82fab7ef26824712b9e056fac5069f..7000141bdb6d86985fb81d1639fedbe0c5c8674d 100644
--- a/srcs/FEM/solverFEM.cpp
+++ b/srcs/FEM/solverFEM.cpp
@@ -190,7 +190,7 @@ void solverFEMnonLinear(std::map<int, double> &electrostaticPressure,
         previousTotalNodalForces[i] = 0;
         currentTotalNodalForces[i] = 0;
     }
-    double residualTolerance = 1e-12; // empirical
+    double residualTolerance = 1e-10; // empirical
 
     double start5 = omp_get_wtime();
     step_name = "further initialization of global variables for the iterative algorithm";