Merge branch 'master' of ssh://git.code.sf.net/p/foam-extend/foam-extend-3.2

[foam-extend-3.2.git] / applications / solvers / solidMechanics / elasticIncrSolidFoam / calculateRelativeResidual.H

{
    scalarField magDU = mag(DU.internalField());

    forAll(magDU, cellI)
    {
        if (magDU[cellI] < SMALL)
        {
            magDU[cellI] = SMALL;
        }
    }

    relativeResidual =
        gMax
        (
            mag
            (
                DU.internalField()
              - DU.prevIter().internalField()
            )/magDU
        );
}