applications/solvers/solidMechanics/elasticOrthoAcpSolidFoam/checkCohesivePatchOrdering.H

   1 {
   2     label sizeByTwo =  mesh.boundary()[cohesivePatchID].size()/2;
   3     vectorField Cczm = mesh.boundary()[cohesivePatchID].Cf();
   4
   5     for(label i = 0; i < sizeByTwo; i++)
   6     {
   7         Cczm[i] -= Cczm[sizeByTwo + i];
   8         Cczm[sizeByTwo + i] = Cczm[i];
   9     }
  10
  11     scalar error = sum(mag(Cczm));
  12
  13     if(error < SMALL)
  14     {
  15         Info << "Cohesive patch face ordering: OK" << endl;
  16     }
  17     else
  18     {
  19         Info << "Cohesive patch face ordering: error (" << error << ")"
  20             << endl;
  21     }
  22 }