checkpoint test using all the norm instead of normal part

NonLinearSolver/materialLaw/mlawNonLocalPorousCoupled.cpp

@@ -413,14 +413,17 @@ void mlawNonLocalPorousCoupledLaw::checkCoalescence(IPNonLocalPorosity* q1, cons
         // Get normal traction on the interface from tau stress
         double sigmax = 0.;
         static SVector3 Norm;
+        static SVector3 tractionForce;
         Norm = q1->getConstRefToReferenceOutwardNormal();
         if (Norm.norm() > 0){ // is it usefulll ???
           Norm.normalize();
           for (int i=0; i<3; i++){
+            tractionForce(i) = 0.;
             for (int j=0; j<3; j++){
-              sigmax += Norm(i)*kCor(i,j)*Norm(j);
+              tractionForce(i) += kCor(i,j)*Norm(j);
             }
           }
+          sigmax = tractionForce.norm();
         }
         else{
           Msg::Fatal("mlawNonLocalPorousCoupledLaw::updatePlasticState: zero normal");