correct ref

dG3D/benchmarks/CMakeLists.txt

@@ -122,7 +122,6 @@ add_subdirectory(microBC_reset)
 add_subdirectory(multiscaleCohesiveTest2D_fullDG_changeBCAfterFailure)
 add_subdirectory(elasticFE2)
 add_subdirectory(multiscaleCohesiveTest2D_fullDG_microPathFollowing)
-add_subdirectory(multiscaleCohesiveTest2D_mixedMode_compression)
 add_subdirectory(compositeHyperDamage)
 add_subdirectory(mixedmodeDelaminationExpoLaw)
 add_subdirectory(hyperViscoElastic)

dG3D/benchmarks/multiscaleCohesiveTest2D_mixedMode/model.py

@@ -94,7 +94,7 @@ macromat1.dirichletBC("Face",11,2,0.)
 
 
 lcohNum = 13
-lawCoh = TwoFieldMultiscaleCohesive3DLaw(lcohNum,0.7)
+lawCoh = TwoFieldMultiscaleCohesive3DLaw(lcohNum)
 lawCoh.setCharacteristicLength(1.)
 lawCoh.setLostSolutionUniquenssTolerance(0.05)
 
@@ -161,14 +161,9 @@ mysolver.options("-ksp_type preonly -pc_type lu -pc_factor_shift_type NONZERO -p
 mysolver.displacementBC("Face",11,2,0.0)
 mysolver.displacementBC("Face",12,2,0.0)
 
-#mysolver.displacementBC("Face",11,1,0.0)
-#mysolver.displacementBC("Face",12,1,0.0)
-
 mysolver.periodicBC("Edge",1,3,1,4,0)
 mysolver.periodicBC("Edge",1,3,1,4,1)
 
-#mysolver.displacementBC("Edge",1,1,0.0)
-#mysolver.displacementBC("Edge",3,1,0.01)
 mysolver.displacementBC("Edge",2,0,0.02)
 mysolver.displacementBC("Edge",2,1,0.02)
 mysolver.displacementBC("Edge",4,0,0.0)
@@ -199,5 +194,5 @@ mysolver.archivingForceOnPhysicalGroup("Edge",4,0)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-3.066874e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-3)
+check.equal(-2.816954e+01,mysolver.getArchivedForceOnPhysicalGroup("Edge", 4, 0),1.e-3)