diff --git a/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py b/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
index 74d6bd453f3e401bcaaf6d6e78ecf9dbfb38d86d..097a6c65758c4bd80743b9f2d69fe2560c67c51e 100644
--- a/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
+++ b/dG3D/benchmarks/GursonThomasonFullCoalescence/Plane_notch.py
@@ -190,4 +190,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-9.076623e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
+check.equal(-8.747734e+02,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-4)
diff --git a/dG3D/benchmarks/GursonThomason_cube/cube.py b/dG3D/benchmarks/GursonThomason_cube/cube.py
index af3e3334e39c10b45ca24d2efec24783e32727c8..d9b943f3cc2c76f638595e03f1db1fc7e3947d86 100644
--- a/dG3D/benchmarks/GursonThomason_cube/cube.py
+++ b/dG3D/benchmarks/GursonThomason_cube/cube.py
@@ -186,7 +186,7 @@ check = TestCheck()
 import csv
 data = csv.reader(open('IPVolume11val_LOCAL_POROSITYMax.csv'), delimiter=';')
 porosity = list(data)
-check.equal(3.800967e-01,float(porosity[-1][1]),1e-6)
+check.equal(3.733499e-01,float(porosity[-1][1]),1e-6)
 
 
 
diff --git a/dG3D/benchmarks/mixedmodeDelamination/delam.py b/dG3D/benchmarks/mixedmodeDelamination/delam.py
index ee5217265ebd9a1f91cec24b8ba8bb67ba2d6853..8a0a561cacdc976f0abf036124b4afec4e66b702 100644
--- a/dG3D/benchmarks/mixedmodeDelamination/delam.py
+++ b/dG3D/benchmarks/mixedmodeDelamination/delam.py
@@ -108,6 +108,6 @@ mysolver.archivingNodeDisplacement(58,2,1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(1.194041e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", 56, 2),1.e-3)
+check.equal(1.196025e+00,mysolver.getArchivedForceOnPhysicalGroup("Face", 56, 2),1.e-3)