diff --git a/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py b/dG3D/benchmarks/GursonMultipleNonlocal_PathFollowing/Plane_notch.py
index 5772a55d79d6680deaf84c914561a42a940b6828..c27a8d35a960ffe8bcf18924923571ffeca68de1 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(-1.846272e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-2.343186e+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 dfdb80ce703810f8f06714da5fcc31fdac94c463..61beb42713c2451c79b3f7a51a3e3dd76afb2c32 100644
--- a/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
+++ b/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
@@ -194,4 +194,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-8.915829e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
+check.equal(-8.920918e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
diff --git a/dG3D/benchmarks/LocalGurson_PathFollowing/Plane_notch.py b/dG3D/benchmarks/LocalGurson_PathFollowing/Plane_notch.py
index 4a77923288854f65c4318a48b2d55a59024ba48f..3fd9d375e3584bdeb0e12004af31b374fd65a0bc 100644
--- a/dG3D/benchmarks/LocalGurson_PathFollowing/Plane_notch.py
+++ b/dG3D/benchmarks/LocalGurson_PathFollowing/Plane_notch.py
@@ -159,4 +159,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-4.264392e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-4.332014e+02,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 a539e6decd75a363029da5a4d1f1067662754434..bbb3c30f0ffb3e9ce7d46ef99d451fed1ffae857 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.797784e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-4)
+check.equal(-5.765123e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-4)
 
 
 
diff --git a/dG3D/benchmarks/conden_ldbc_SecondOrder/run.py b/dG3D/benchmarks/conden_ldbc_SecondOrder/run.py
index 9a291643138e57d6a419dc9007c0b225abe27a17..2f889ed1ada61f367e1ec065b5db16400dff06ae 100644
--- a/dG3D/benchmarks/conden_ldbc_SecondOrder/run.py
+++ b/dG3D/benchmarks/conden_ldbc_SecondOrder/run.py
@@ -107,7 +107,7 @@ check.equal(1.144128e+03,mysolver.getHomogenizedTangent(0,0,1,1),1.e-4)
 check.equal(9.200879e+02,mysolver.getHomogenizedTangent(0,1,0,1),1.e-4)
 
 check.equal(1.720444e+02,mysolver.getHomogenizedSecondTangent(0,0,0,0,0,0),1.e-4)
-check.equal(83.9775345211,mysolver.getHomogenizedSecondTangent(1,0,1,1,0,1),1.e-4)
+check.equal(8.453953e+01,mysolver.getHomogenizedSecondTangent(1,0,1,1,0,1),1.e-4)
 check.equal(78.2909010986,mysolver.getHomogenizedSecondTangent(1,0,0,1,0,0),1.e-4)
 check.equal(28.491995815,mysolver.getHomogenizedSecondTangent(0,1,0,1,0,0),1.e-4)
 
diff --git a/dG3D/benchmarks/conden_mkbc_SecondOrder/run.py b/dG3D/benchmarks/conden_mkbc_SecondOrder/run.py
index 45d2cd1cf3e4a97a06baadb6dff58b0f723a61d4..9d89b159db0e92c622e7179fa914e5450da3cd1d 100644
--- a/dG3D/benchmarks/conden_mkbc_SecondOrder/run.py
+++ b/dG3D/benchmarks/conden_mkbc_SecondOrder/run.py
@@ -107,4 +107,4 @@ check.equal(1.198088e+03,mysolver.getHomogenizedTangent(0,0,1,1),1.e-4)
 check.equal(8.055230e+02,mysolver.getHomogenizedTangent(0,1,0,1),1.e-4)
 
 check.equal(88.61673754442836,mysolver.getHomogenizedSecondTangent(0,0,0,0,0,0),1.e-4)
-check.equal(47.2887234489,mysolver.getHomogenizedSecondTangent(1,0,0,1,0,0),1.e-4)
+check.equal(4.777566e+01,mysolver.getHomogenizedSecondTangent(1,0,0,1,0,0),1.e-4)
diff --git a/dG3D/benchmarks/electroThermoLaminateFtFv/electroThermoLaminate.py b/dG3D/benchmarks/electroThermoLaminateFtFv/electroThermoLaminate.py
index 1062a15034cb98c66a417b8fd005464c80d0a0d9..59459f6efbb6873d9d967007007dd5dd6738b643 100644
--- a/dG3D/benchmarks/electroThermoLaminateFtFv/electroThermoLaminate.py
+++ b/dG3D/benchmarks/electroThermoLaminateFtFv/electroThermoLaminate.py
@@ -247,7 +247,7 @@ mysolver.solve()
 
 
 check = TestCheck()
-check.equal(-6.476172e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 3),1.e-5)
+check.equal(-6.476339e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 3),1.e-5)
 check.equal(-2.097632e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 85, 4),1.e-5)
 check.equal(3.355705e-03,mysolver.getArchivedNodalValue(88,3,mysolver.displacement),1.e-5)
 check.equal(-3.355705e-02,mysolver.getArchivedNodalValue(88,4,mysolver.displacement),1.e-5)
diff --git a/dG3D/benchmarks/electroThermoSMP/electroThermoSMP.py b/dG3D/benchmarks/electroThermoSMP/electroThermoSMP.py
index 80831b98bac765da455fb60657a72a08e16d00eb..22054131a9f6c0e274a526bd89540ea5a935b179 100644
--- a/dG3D/benchmarks/electroThermoSMP/electroThermoSMP.py
+++ b/dG3D/benchmarks/electroThermoSMP/electroThermoSMP.py
@@ -223,7 +223,7 @@ mysolver.solve()
 check = TestCheck()
 check.equal(0.,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 0),1.e-6)
 check.equal(-2.462150e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 1),1.e-6)
-check.equal(2.193303e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 2),1.e-6)
+check.equal(3.584923e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 2),1.e-6)
 check.equal(4.282709e+04,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 3),1.e-6)
 check.equal(4.429112e+05,mysolver.getArchivedForceOnPhysicalGroup("Face", 247, 4),1.e-6)
 check.equal(6.422382e-09,mysolver.getArchivedNodalValue(93,0,mysolver.displacement),1.e-6)
diff --git a/dG3D/benchmarks/twoHole/twoHole.py b/dG3D/benchmarks/twoHole/twoHole.py
index fbc0fd057f028eb2f9d21d0f178842b8818e119d..93740810077c6986ed2c8db5d260a7435da0c0cd 100644
--- a/dG3D/benchmarks/twoHole/twoHole.py
+++ b/dG3D/benchmarks/twoHole/twoHole.py
@@ -107,6 +107,6 @@ mysolver.archivingNodeDisplacement(57,1,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-6.105636e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 55, 1),6e-2)
+check.equal(-6.629669e-05,mysolver.getArchivedForceOnPhysicalGroup("Face", 55, 1),6e-2)
 check.equal(3e-6,mysolver.getArchivedNodalValue(57,1,mysolver.displacement),1e-5)