applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/calculateRelativeResidual.H

   1 {
   2     scalarField magDU = mag(DU.internalField());
   3
   4     forAll(magDU, cellI)
   5     {
   6         if (magDU[cellI] < SMALL)
   7         {
   8             magDU[cellI] = SMALL;
   9         }
  10     }
  11
  12     relativeResidual =
  13         gMax
  14         (
  15             mag
  16             (
  17                 DU.internalField()
  18               - DU.prevIter().internalField()
  19             )
  20            /magDU
  21         );
  22
  23 //     resid = (DU - DU.prevIter())/mag(DU);
  24 }