41 if (lowerPtr_ || upperPtr_)
45 Type* __restrict__ HpsiPtr = Hpsi.
begin();
47 const Type* __restrict__ psiPtr = psi.
begin();
52 const scalar* __restrict__ lowerPtr =
lower().
begin();
53 const scalar* __restrict__ upperPtr =
upper().
begin();
55 register const label nFaces =
upper().
size();
59 HpsiPtr[uPtr[
face]] -= lowerPtr[
face]*psiPtr[lPtr[
face]];
60 HpsiPtr[lPtr[
face]] -= upperPtr[
face]*psiPtr[uPtr[
face]];
81 if (lowerPtr_ || upperPtr_)
101 FatalErrorIn(
"lduMatrix::faceH(const Field<Type>& psi) const")
102 <<
"Cannot calculate faceH"
103 " the matrix does not have any off-diagonal coefficients."