34 lduMatrix::preconditioner::
35 addsymMatrixConstructorToTable<FDICPreconditioner>
49 rD_(sol.matrix().
diag()),
50 rDuUpper_(sol.matrix().upper().size()),
51 rDlUpper_(sol.matrix().upper().size())
53 scalar* __restrict__ rDPtr = rD_.
begin();
54 scalar* __restrict__ rDuUpperPtr = rDuUpper_.
begin();
55 scalar* __restrict__ rDlUpperPtr = rDlUpper_.
begin();
57 const label*
const __restrict__ uPtr =
59 const label*
const __restrict__ lPtr =
61 const scalar*
const __restrict__ upperPtr =
95 scalar* __restrict__ wAPtr = wA.
begin();
96 const scalar* __restrict__ rAPtr = rA.
begin();
97 const scalar* __restrict__ rDPtr = rD_.begin();
99 const label*
const __restrict__ uPtr =
100 solver_.matrix().lduAddr().upperAddr().begin();
101 const label*
const __restrict__ lPtr =
102 solver_.matrix().lduAddr().lowerAddr().begin();
104 const scalar*
const __restrict__ rDuUpperPtr = rDuUpper_.begin();
105 const scalar*
const __restrict__ rDlUpperPtr = rDlUpper_.begin();
108 label nFaces = solver_.matrix().upper().size();
109 label nFacesM1 = nFaces - 1;
118 wAPtr[uPtr[
face]] -= rDuUpperPtr[
face]*wAPtr[lPtr[
face]];
123 wAPtr[lPtr[
face]] -= rDlUpperPtr[
face]*wAPtr[uPtr[
face]];
Faster version of the DICPreconditioner diagonal-based incomplete Cholesky preconditioner for symmetr...
virtual void precondition(scalarField &wA, const scalarField &rA, const direction cmpt=0) const
Return wA the preconditioned form of residual rA.
FDICPreconditioner(const lduMatrix::solver &, const dictionary &solverControlsUnused)
Construct from matrix components and preconditioner solver controls.
void size(const label)
Override size to be inconsistent with allocated storage.
iterator begin()
Return an iterator to begin traversing the UList.
A cell is defined as a list of faces with extra functionality.
A list of keyword definitions, which are a keyword followed by any number of values (e....
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.
const solver & solver_
Reference to the base-solver this preconditioner is used with.
Abstract base-class for lduMatrix solvers.
const lduMatrix & matrix() const
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
const lduAddressing & lduAddr() const
Return the LDU addressing.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
lduMatrix::preconditioner::addsymMatrixConstructorToTable< FDICPreconditioner > addFDICPreconditionerSymMatrixConstructorToTable_
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)
defineTypeNameAndDebug(combustionModel, 0)