diff --git a/dG3D/benchmarks/GursonPlasticInstabilities_PathFollowing/Plane_notch.py b/dG3D/benchmarks/GursonPlasticInstabilities_PathFollowing/Plane_notch.py
index 834633c08e7080c7b5dfe7a7dd45cf82eb6592b4..9db524010e423a289abde705c91577f838127fa3 100644
--- a/dG3D/benchmarks/GursonPlasticInstabilities_PathFollowing/Plane_notch.py
+++ b/dG3D/benchmarks/GursonPlasticInstabilities_PathFollowing/Plane_notch.py
@@ -190,4 +190,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-1.957893e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-1.957854e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
diff --git a/dG3D/benchmarks/MFH_FiniteStrainDPS_20/cube.py b/dG3D/benchmarks/MFH_FiniteStrainDPS_20/cube.py
index 26a7c53f8f8e7d2368ea719df71dd1c0b88bf6e6..2afc6b1a1f8b4a221892eab57656b9b364675f5b 100644
--- a/dG3D/benchmarks/MFH_FiniteStrainDPS_20/cube.py
+++ b/dG3D/benchmarks/MFH_FiniteStrainDPS_20/cube.py
@@ -119,5 +119,5 @@ else:
   check.equal(1.999675e-01,float(linesVF.split(';')[1]))
   check.equal(9.882490e-01,float(linesAR1.split(';')[1]))
   check.equal(9.882490e-01,float(linesAR2.split(';')[1]))
-  check.equal(2.570723e-01,float(linesRXX.split(';')[1]))
+  check.equal(2.489432e-01,float(linesRXX.split(';')[1]))
 
diff --git a/dG3D/benchmarks/MFH_FiniteStrainDPS_El20/cube.py b/dG3D/benchmarks/MFH_FiniteStrainDPS_El20/cube.py
index e866a8b58eba92cbdd39d4549c4e302e3f745497..e0b041af1ee548c253808a6893f9143b60715643 100644
--- a/dG3D/benchmarks/MFH_FiniteStrainDPS_El20/cube.py
+++ b/dG3D/benchmarks/MFH_FiniteStrainDPS_El20/cube.py
@@ -132,7 +132,7 @@ else:
   check.equal(2.002383e-01,float(linesVF.split(';')[1]))
   check.equal(3.033909e+00,float(linesAR1.split(';')[1]))
   check.equal(3.033909e+00,float(linesAR2.split(';')[1]))
-  check.equal(-1.092195e-01,float(linesRXX.split(';')[1]))
+  check.equal(4.385424e-01,float(linesRXX.split(';')[1]))
   check.equal(1.723483e+09,float(linesISIG.split(';')[1]))
   check.equal(1.012394e+09,float(linesMSIG.split(';')[1]))
 
diff --git a/dG3D/benchmarks/MFH_FiniteStrainMMC1/cube.py b/dG3D/benchmarks/MFH_FiniteStrainMMC1/cube.py
index 6bac407c3ebde804edc50b179e2ce0df87780f70..fb7fa94c961d165a6763719eccdd2639f80b018a 100644
--- a/dG3D/benchmarks/MFH_FiniteStrainMMC1/cube.py
+++ b/dG3D/benchmarks/MFH_FiniteStrainMMC1/cube.py
@@ -127,7 +127,7 @@ else:
   check.equal(1.998084e-01,float(linesVF.split(';')[1]))
   check.equal(9.986400e-01,float(linesAR1.split(';')[1]))
   check.equal(9.986400e-01,float(linesAR2.split(';')[1]))
-  check.equal(6.897763e-01,float(linesRXX.split(';')[1]))
+  check.equal(6.863363e-01,float(linesRXX.split(';')[1]))
   check.equal(3.600329e+08,float(linesSigZZ.split(';')[1]))
   check.equal(5.248672e+08,float(linesISIG.split(';')[1]))
   check.equal(3.188737e+08,float(linesMSIG.split(';')[1]))
diff --git a/dG3D/benchmarks/torchOpenHole/model.py b/dG3D/benchmarks/torchOpenHole/model.py
index e7433df4c84ffc81fb44ea0a0ede9800635a47b1..9e415d4eb34e9528991a21e58a6d55c9624333ad 100644
--- a/dG3D/benchmarks/torchOpenHole/model.py
+++ b/dG3D/benchmarks/torchOpenHole/model.py
@@ -152,7 +152,7 @@ mysolver.archivingNodeDisplacement(6,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-5.408420e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 2, 1),1.e-4)
+check.equal(-5.412217e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 2, 1),1.e-4)