diff --git a/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py b/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py
index ab55dfe8275854688f324e9546b3125a5fa4dde1..7c0f828e621e7757517831a8a6e4e709ec77f480 100755
--- a/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py
+++ b/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py
@@ -205,7 +205,7 @@ mysolver.archivingNodeDisplacement(73,2,nstepArchForce)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-6.666707e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 61, 1),1.e-4)
+check.equal(-7.725560e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 61, 1),1.e-4)
 
 
 
diff --git a/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py b/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py
index 294bc3b553e033387900908d0e238b987a68a916..0f685b6e06b247f7f02b8d6868f6656130db003c 100644
--- a/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py
+++ b/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py
@@ -180,4 +180,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-2.084549e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-2.084728e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
diff --git a/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py b/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
index 88b81a0a9dbe6ba814381242ef8fc962f4984420..64722f9441dc431d91d398a447fedbde2a79b600 100644
--- a/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
+++ b/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
@@ -195,4 +195,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-9.131203e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
+check.equal(-9.176614e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
diff --git a/dG3D/benchmarks/GursonThomasonThermo_NonlocalLogPorosity/Plane_notch.py b/dG3D/benchmarks/GursonThomasonThermo_NonlocalLogPorosity/Plane_notch.py
index a3c1fb718f13386de49c7737be9ecf36d248bdca..7d79d2d44560adc352b3f2c6c8101dc3122e17c5 100644
--- a/dG3D/benchmarks/GursonThomasonThermo_NonlocalLogPorosity/Plane_notch.py
+++ b/dG3D/benchmarks/GursonThomasonThermo_NonlocalLogPorosity/Plane_notch.py
@@ -196,4 +196,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-1.912441e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-1.912437e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
diff --git a/dG3D/benchmarks/Thomason_planeStrain_withLodeEnhanced/model.py b/dG3D/benchmarks/Thomason_planeStrain_withLodeEnhanced/model.py
index 9f179dc9b3b015cf1d19af511404ed51e5a32078..2ff7930f107e0623c1e97819860736fb315f3812 100644
--- a/dG3D/benchmarks/Thomason_planeStrain_withLodeEnhanced/model.py
+++ b/dG3D/benchmarks/Thomason_planeStrain_withLodeEnhanced/model.py
@@ -164,7 +164,7 @@ mysolver.solve()
 
 
 check = TestCheck()
-check.equal(-5.738027e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-4)
+check.equal(-5.774703e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-4)
 
 
 
diff --git a/dG3D/benchmarks/pathFollowing/arch3D.py b/dG3D/benchmarks/pathFollowing/arch3D.py
index 91ead88698a3aea4bca7e468951e5358646bb5e8..28a3f8d1e95eb163b810ff7e680ef7e6bf32e840 100644
--- a/dG3D/benchmarks/pathFollowing/arch3D.py
+++ b/dG3D/benchmarks/pathFollowing/arch3D.py
@@ -85,4 +85,4 @@ mysolver.archivingForceOnPhysicalGroup('Face',57,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(1.340224e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", 57, 1),1.e-4)
+check.equal(1.481245e+01,mysolver.getArchivedForceOnPhysicalGroup("Face", 57, 1),1.e-4)
diff --git a/dG3D/benchmarks/pathFollowing_cohesive_selectiveUpdate/model.py b/dG3D/benchmarks/pathFollowing_cohesive_selectiveUpdate/model.py
index 86a8f23259314153eb81f739c695003d2a4e78e6..da9c9a22b7b22f0d5299aa52b611d63a862f00d0 100644
--- a/dG3D/benchmarks/pathFollowing_cohesive_selectiveUpdate/model.py
+++ b/dG3D/benchmarks/pathFollowing_cohesive_selectiveUpdate/model.py
@@ -132,4 +132,4 @@ mysolver.archivingForceOnPhysicalGroup('Face',1,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(2.227505e-11,mysolver.getArchivedForceOnPhysicalGroup("Face", 1, 1),1.e-2)
+check.equal(-7.680655e-11,mysolver.getArchivedForceOnPhysicalGroup("Face", 1, 1),1.e1)