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

[foam-extend-3.2.git] / applications / solvers / solidMechanics / deprecatedSolvers / contactStressFoam / calculatePatchTraction.H

    Info << "Calculating patch traction" << endl;
    PtrList<vectorField> traction(mesh.boundary().size());

    {
        const surfaceVectorField::GeometricBoundaryField& Apatches =
            mesh.Sf().boundaryField();

        const surfaceScalarField::GeometricBoundaryField& magApatches =
            mesh.magSf().boundaryField();

        forAll(traction, patchI)
        {
            vectorField n = Apatches[patchI]/magApatches[patchI];

            traction.set
            (
                patchI,
                new vectorField
                (
                    patchTraction[patchI]
                  - patchPressure[patchI]*n
                )
            );
        }
    }