39 lduMatrix::smoother::addsymMatrixConstructorToTable<FDICSmoother>
48 const word& fieldName,
63 rD_(matrix_.
diag().size()),
64 rDuUpper_(matrix_.upper().size()),
65 rDlUpper_(matrix_.upper().size())
67 solveScalar* __restrict__ rDPtr = rD_.
begin();
68 solveScalar* __restrict__ rDuUpperPtr = rDuUpper_.
begin();
69 solveScalar* __restrict__ rDlUpperPtr = rDlUpper_.
begin();
71 const label*
const __restrict__ uPtr =
73 const label*
const __restrict__ lPtr =
75 const scalar*
const __restrict__ upperPtr =
103 const solveScalar*
const __restrict__ rDuUpperPtr = rDuUpper_.begin();
104 const solveScalar*
const __restrict__ rDlUpperPtr = rDlUpper_.begin();
106 const label*
const __restrict__ uPtr =
107 matrix_.lduAddr().upperAddr().begin();
108 const label*
const __restrict__ lPtr =
109 matrix_.lduAddr().lowerAddr().begin();
113 solveScalar* __restrict__ rAPtr = rA.
begin();
115 for (label sweep=0; sweep<nSweeps; sweep++)
135 rAPtr[uPtr[
face]] -= rDuUpperPtr[
face]*rAPtr[lPtr[
face]];
138 const label nFacesM1 =
nFaces - 1;
141 rAPtr[lPtr[
face]] -= rDlUpperPtr[
face]*rAPtr[uPtr[
face]];
A const Field/List wrapper with possible data conversion.
static void calcReciprocalD(solveScalarField &, const lduMatrix &)
Calculate the reciprocal of the preconditioned diagonal.
Simplified diagonal-based incomplete Cholesky smoother for symmetric matrices.
void smooth(solveScalarField &psi, const scalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
virtual void scalarSmooth(solveScalarField &psi, const solveScalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
A field of fields is a PtrList of fields with reference counting.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
void size(const label n)
Older name for setAddressableSize.
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 lduMatrix & matrix_
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
const lduAddressing & lduAddr() const
Return the LDU addressing.
A class for handling words, derived from Foam::string.
#define defineTypeNameAndDebug(Type, DebugSwitch)
Define the typeName and debug information.
const volScalarField & psi
const labelList nFaces(UPstream::listGatherValues< label >(aMesh.nFaces()))
lduMatrix::smoother::addsymMatrixConstructorToTable< FDICSmoother > addFDICSmootherSymMatrixConstructorToTable_
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)
#define forAll(list, i)
Loop across all elements in list.