38 lduMatrix::preconditioner::
39 addsymMatrixConstructorToTable<DICPreconditioner>
53 rD_(sol.matrix().
diag().size())
70 solveScalar* __restrict__ rDPtr = rD.
begin();
74 const scalar*
const __restrict__ upperPtr = matrix.
upper().
begin();
85 const label nCells = rD.
size();
101 solveScalar* __restrict__ wAPtr = wA.
begin();
102 const solveScalar* __restrict__ rAPtr = rA.
begin();
103 const solveScalar* __restrict__ rDPtr = rD_.begin();
105 const label*
const __restrict__ uPtr =
106 solver_.matrix().lduAddr().upperAddr().begin();
107 const label*
const __restrict__ lPtr =
108 solver_.matrix().lduAddr().lowerAddr().begin();
109 const scalar*
const __restrict__ upperPtr =
110 solver_.matrix().upper().begin();
112 const label nCells = wA.
size();
113 const label
nFaces = solver_.matrix().upper().
size();
114 const label nFacesM1 =
nFaces - 1;
Simplified diagonal-based incomplete Cholesky preconditioner for symmetric matrices (symmetric equiva...
virtual void precondition(solveScalarField &wA, const solveScalarField &rA, const direction cmpt=0) const
Return wA the preconditioned form of residual rA.
static void calcReciprocalD(solveScalarField &, const lduMatrix &)
Calculate the reciprocal of the preconditioned diagonal.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
void size(const label n)
Older name for setAddressableSize.
A cell is defined as a list of faces with extra functionality.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
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.
Abstract base-class for lduMatrix solvers.
const lduMatrix & matrix() const noexcept
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
const lduAddressing & lduAddr() const
Return the LDU addressing.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
const labelList nFaces(UPstream::listGatherValues< label >(aMesh.nFaces()))
lduMatrix::preconditioner::addsymMatrixConstructorToTable< DICPreconditioner > addDICPreconditionerSymMatrixConstructorToTable_
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)