diff --git a/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py b/dG3D/benchmarks/CohesiveBand_Plate/notchedPlate.py
index 7c0f828e621e7757517831a8a6e4e709ec77f480..aae9cd42b8ce09f5665e153132c3fd6119cb5260 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(-7.725560e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 61, 1),1.e-4)
+check.equal(-7.525179e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 61, 1),1.e-4)
 
 
 
diff --git a/dG3D/benchmarks/GursonThomasonMultipleNonlocalVar_PathFollowing/Plane_notch.py b/dG3D/benchmarks/GursonThomasonMultipleNonlocalVar_PathFollowing/Plane_notch.py
index 617ccb838f2748ce770fa96cd52879b16942e1c3..e6dadd6a33728babc34a214f5c9ca2217f64e39e 100644
--- a/dG3D/benchmarks/GursonThomasonMultipleNonlocalVar_PathFollowing/Plane_notch.py
+++ b/dG3D/benchmarks/GursonThomasonMultipleNonlocalVar_PathFollowing/Plane_notch.py
@@ -189,4 +189,4 @@ mysolver.archivingIPOnPhysicalGroup("Face", 700, IPField.LOCAL_POROSITY,IPField.
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-1.889379e+03,3,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
+check.equal(-2.558598e+03,mysolver.getArchivedForceOnPhysicalGroup("Edge", 1400, 0),1.e-6)
diff --git a/dG3D/benchmarks/GursonThomason_cube_withLodeEnhanced/cube.py b/dG3D/benchmarks/GursonThomason_cube_withLodeEnhanced/cube.py
index 109f885467221055e20d6c5499d7f6711159e471..57c2b2aeac362f9b16a0c7244f11f6b5be4aa624 100644
--- a/dG3D/benchmarks/GursonThomason_cube_withLodeEnhanced/cube.py
+++ b/dG3D/benchmarks/GursonThomason_cube_withLodeEnhanced/cube.py
@@ -198,7 +198,7 @@ check = TestCheck()
 import csv
 data = csv.reader(open('IPVolume11val_LOCAL_POROSITYMax.csv'), delimiter=';')
 porosity = list(data)
-check.equal(5.856826e-01,float(porosity[-1][1]),1e-6)
+check.equal(5.473613e-01,float(porosity[-1][1]),1e-6)
 
 
 
diff --git a/dG3D/benchmarks/Gurson_Cube_Transition/cube.py b/dG3D/benchmarks/Gurson_Cube_Transition/cube.py
index 1e5c4a79a0ddfbc26d8f2ee2aa3aac403d9e53f4..50741c868a7783a4d1729f591d1c6b1ec3c82ebe 100755
--- a/dG3D/benchmarks/Gurson_Cube_Transition/cube.py
+++ b/dG3D/benchmarks/Gurson_Cube_Transition/cube.py
@@ -169,7 +169,7 @@ check = TestCheck()
 import csv
 data = csv.reader(open('IPVolume29val_DAMAGEMean.csv'), delimiter=';')
 porosity = list(data)
-check.equal(2.759742e-02,float(porosity[-1][1]),1e-6)
+check.equal(2.758505e-02,float(porosity[-1][1]),1e-6)
 
 
 
diff --git a/dG3D/benchmarks/taylorMPI/taylorCG.py b/dG3D/benchmarks/taylorMPI/taylorCG.py
index e1da15d9a40b2f12b66221c017877779c85500ed..c54ce21acc174db3c373ba7c556b0d95a2314975 100644
--- a/dG3D/benchmarks/taylorMPI/taylorCG.py
+++ b/dG3D/benchmarks/taylorMPI/taylorCG.py
@@ -73,7 +73,7 @@ mysolver.solve()
 # perform some check (pay attention that only rank 0 really perform the check)
 try:
   import linecache
-  linedisp = linecache.getline('NodalDisplacement1comp0_part3.csv',26050)
+  linedisp = linecache.getline('NodalDisplacementPhysical5Num1comp0_part3.csv',26050)
 except:
   print('Cannot read the results in the files')
   import os