43 if (lowerPtr_ || upperPtr_)
47 Type* __restrict__ HpsiPtr = Hpsi.
begin();
49 const Type* __restrict__ psiPtr =
psi.
begin();
54 const scalar* __restrict__ lowerPtr =
lower().
begin();
55 const scalar* __restrict__ upperPtr =
upper().
begin();
61 HpsiPtr[uPtr[
face]] -= lowerPtr[
face]*psiPtr[lPtr[
face]];
62 HpsiPtr[lPtr[
face]] -= upperPtr[
face]*psiPtr[uPtr[
face]];
83 if (lowerPtr_ || upperPtr_)
105 <<
"Cannot calculate faceH"
106 " the matrix does not have any off-diagonal coefficients."
Generic templated field type.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
void size(const label n)
Older name for setAddressableSize.
tmp< GeometricField< Type, faPatchField, areaMesh > > H() const
Return the H operation source.
A face is a list of labels corresponding to mesh vertices.
virtual const labelUList & upperAddr() const =0
Return upper addressing.
virtual const labelUList & lowerAddr() const =0
Return lower addressing.
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
tmp< Field< Type > > faceH(const Field< Type > &) const
const lduAddressing & lduAddr() const
Return the LDU addressing.
A class for managing temporary objects.
const volScalarField & psi
volScalarField H(IOobject("H", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimLength, Zero))
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
const labelList nFaces(UPstream::listGatherValues< label >(aMesh.nFaces()))
static constexpr const zero Zero
Global zero (0)
errorManipArg< error, int > exit(error &err, const int errNo=1)