diff --git a/dG3D/benchmarks/GursonThomasonShearTransition/model.py b/dG3D/benchmarks/GursonThomasonShearTransition/model.py
index 3587c22ab2de93754a3d4b9a2950530f930478c0..d8e95bd726aa1582081b9c75185dbb504e5f14c4 100755
--- a/dG3D/benchmarks/GursonThomasonShearTransition/model.py
+++ b/dG3D/benchmarks/GursonThomasonShearTransition/model.py
@@ -194,8 +194,8 @@ mysolver.solve()
 
 
 check = TestCheck()
-check.equal(-1.915989e-04,mysolver.getArchivedNodalValue(2,0,mysolver.displacement),1.e-4)
-check.equal(1.660548e-04,mysolver.getArchivedNodalValue(3,0,mysolver.displacement),1.e-4)
+check.equal(-1.870109e-04,mysolver.getArchivedNodalValue(2,0,mysolver.displacement),1.e-4)
+check.equal(1.685898e-04,mysolver.getArchivedNodalValue(3,0,mysolver.displacement),1.e-4)
 
 
 
diff --git a/dG3D/benchmarks/HOwithConstantBodyForce/run.py b/dG3D/benchmarks/HOwithConstantBodyForce/run.py
index 3e03c4eac58f8bb76fbd42df75144a330ab64cb9..a148dde4ffa06dac9a7747f98145029597a570af 100644
--- a/dG3D/benchmarks/HOwithConstantBodyForce/run.py
+++ b/dG3D/benchmarks/HOwithConstantBodyForce/run.py
@@ -182,7 +182,7 @@ mysolver.solve()
 import linecache
 total_line_number = sum(1 for line in open('E_0_GP_0_stress.csv'))
 homoStress = linecache.getline('E_0_GP_0_stress.csv',total_line_number)
-print(homoStress)
+#print(homoStress)
 val = float(homoStress.split(';')[1])
 check = TestCheck()
 check.equal(-0.014840731,val,1.e-4)
diff --git a/dG3D/benchmarks/nanoIndentation/ni.py b/dG3D/benchmarks/nanoIndentation/ni.py
index 221c298661066079534c45d7a4206d0b485f2903..0c74b577ebaf7468f9268658777326c85cf8bece 100644
--- a/dG3D/benchmarks/nanoIndentation/ni.py
+++ b/dG3D/benchmarks/nanoIndentation/ni.py
@@ -166,7 +166,7 @@ mysolver.archivingRigidContactForce(22, 2, 1)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-2.820332e-04,mysolver.getArchivedNodalValue(1001,2,mysolver.displacement),1.e-4)
+check.equal(-2.711273e-04,mysolver.getArchivedNodalValue(1001,2,mysolver.displacement),1.e-4)