40 lduMatrix::smoother::addasymMatrixConstructorToTable<DILUSmoother>
49 const word& fieldName,
64 rD_(matrix_.
diag().size())
83 const solveScalar*
const __restrict__ rDPtr = rD_.begin();
85 const label*
const __restrict__ uPtr =
86 matrix_.lduAddr().upperAddr().begin();
87 const label*
const __restrict__ lPtr =
88 matrix_.lduAddr().lowerAddr().begin();
90 const scalar*
const __restrict__ upperPtr = matrix_.upper().begin();
91 const scalar*
const __restrict__ lowerPtr = matrix_.lower().begin();
95 solveScalar* __restrict__ rAPtr = rA.
begin();
114 const label nFaces = matrix_.upper().size();
115 for (label face=0; face<nFaces; face++)
117 const label u = uPtr[face];
118 rAPtr[u] -= rDPtr[u]*lowerPtr[face]*rAPtr[lPtr[face]];
121 const label nFacesM1 = nFaces - 1;
122 for (label face=nFacesM1; face>=0; face--)
124 const label l = lPtr[face];
125 rAPtr[l] -= rDPtr[l]*upperPtr[face]*rAPtr[uPtr[face]];
144 ConstPrecisionAdaptor<scalar, solveScalar>(source),
const scalarField & diag() const
Simplified diagonal-based incomplete LU smoother for asymmetric matrices.
static void calcReciprocalD(solveScalarField &, const lduMatrix &)
Calculate the reciprocal of the preconditioned diagonal.
Field< solveScalar > solveScalarField
DILUSmoother(const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces)
Construct from matrix components.
A field of fields is a PtrList of fields with reference counting.
void smooth(volScalarField &field, const scalar coeff)
#define forAll(list, i)
Loop across all elements in list.
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)
const lduMatrix & matrix_
void sweep(volScalarField &field, const volScalarField &alpha, const label nLayers, const scalar alphaDiff=0.2)
A class for handling words, derived from Foam::string.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
void smooth(solveScalarField &psi, const scalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
defineTypeNameAndDebug(combustionModel, 0)
lduMatrix::smoother::addasymMatrixConstructorToTable< DILUSmoother > addDILUSmootherAsymMatrixConstructorToTable_
lduMatrix is a general matrix class in which the coefficients are stored as three arrays...
void scalarSmooth(solveScalarField &psi, const solveScalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
const volScalarField & psi