applications/solvers/solidMechanics/elasticIncrAcpSolidFoam/calculateRelativeResidual.H

   1 {
   2     scalar maxDU = gMax(mag(DU.internalField()));
   3
   4     relativeResidual =
   5         gMax
   6         (
   7             mag(DU.internalField() - DU.prevIter().internalField())
   8            /(maxDU + SMALL)
   9         );
  10
  11     if (lduMatrix::debug == 1)
  12     {
  13         Info << "Relative residual = " << relativeResidual << endl;
  14     }
  15 }