diff --git a/dG3D/benchmarks/defoDefoContactSlaveEdge2DQuadPert/squares.py b/dG3D/benchmarks/defoDefoContactSlaveEdge2DQuadPert/squares.py
index 0a12e0e09347f72078e49cad824670ce50a99196..612e3583cb931b2e425db33ba3dc6ffe0a003b25 100644
--- a/dG3D/benchmarks/defoDefoContactSlaveEdge2DQuadPert/squares.py
+++ b/dG3D/benchmarks/defoDefoContactSlaveEdge2DQuadPert/squares.py
@@ -32,6 +32,7 @@ space1 = 0 # function space (Lagrange=0)
 # creation of law
 
 law1 = dG3DLinearElasticMaterialLaw(lawnum,rho,E,nu)
+law1.setUseBarF(True)
 
 # creation of ElasticField
 nfield1 = 1265 # number of the field (physical number of surface)
@@ -44,7 +45,7 @@ myfield2 = dG3DDomain(1100,nfield2,space1,lawnum,fullDg,2)
 
 # creation of Solver
 mysolver = nonLinearMechSolver(1200)
-mysolver.createModel(geofile,meshfile,2,2)
+mysolver.createModel(geofile,meshfile,2,1)
 #mysolver.loadModel(meshfile)
 mysolver.addDomain(myfield1)
 mysolver.addDomain(myfield2)
@@ -94,6 +95,6 @@ mysolver.archivingForceOnPhysicalGroup("Edge", 156, 2)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(2.418994e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 2),1.e-6)
+check.equal(2.451079e+07,mysolver.getArchivedForceOnPhysicalGroup("Edge", 12, 2),1.e-6)