OpenFOAM guide/LduMatrix;;faceH

From OpenFOAMWiki

A function built into the matrix to facilitate implementation of the H operator. Used by fvMatrix;;flux.

for (register label face=0; face<l.size(); face++)
{
    faceHpsi[face] = Upper[face]*psi[u[face]]
        - Lower[face]*psi[l[face]];
}

\mathbf{H}_f \left ( \boldsymbol \psi \right ) =
A_n \boldsymbol \psi_n - A_o \boldsymbol \psi_o

Where: