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]]; }
Where:
- A are matrix coefficients;
- o and n are the owner and neighbour indices; and
- f is a face index.