Go to the documentation of this file.
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();
57 const label nFaces =
upper().size();
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."
const lduAddressing & lduAddr() const
Return the LDU addressing.
A class for managing temporary objects.
static constexpr const zero Zero
Global zero (0)
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
volScalarField H(IOobject("H", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), mesh, dimensionedScalar(dimLength, Zero))
virtual const labelUList & upperAddr() const =0
Return upper addressing.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
Generic templated field type.
tmp< Field< Type > > H(const Field< Type > &) const
errorManipArg< error, int > exit(error &err, const int errNo=1)
#define FatalErrorInFunction
Report an error message using Foam::FatalError.
virtual const labelUList & lowerAddr() const =0
Return lower addressing.
tmp< Field< Type > > faceH(const Field< Type > &) const
A face is a list of labels corresponding to mesh vertices.
void size(const label n)
Older name for setAddressableSize.
const volScalarField & psi