Go to the documentation of this file.
40 lduMatrix::smoother::addasymMatrixConstructorToTable<DILUSmoother>
49 const word& fieldName,
64 rD_(matrix_.diag().size())
67 std::copy(
diag.begin(),
diag.end(), rD_.begin());
69 DILUPreconditioner::calcReciprocalD(rD_, matrix_);
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();
117 const label u = uPtr[
face];
118 rAPtr[u] -= rDPtr[u]*lowerPtr[
face]*rAPtr[lPtr[
face]];
121 const label nFacesM1 = nFaces - 1;
124 const label l = lPtr[
face];
125 rAPtr[l] -= rDPtr[l]*upperPtr[
face]*rAPtr[uPtr[
face]];
A class for handling words, derived from Foam::string.
A field of fields is a PtrList of fields with reference counting.
lduMatrix::smoother::addasymMatrixConstructorToTable< DILUSmoother > addDILUSmootherAsymMatrixConstructorToTable_
void diag(pointPatchField< vector > &, const pointPatchField< tensor > &)
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
#define forAll(list, i)
Loop across all elements in list.
Abstract base-class for lduMatrix smoothers.
void smooth(solveScalarField &psi, const scalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
void scalarSmooth(solveScalarField &psi, const solveScalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
A const Field wrapper with possible data conversion.
A face is a list of labels corresponding to mesh vertices.
DILUSmoother(const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces)
Construct from matrix components.
const volScalarField & psi
defineTypeNameAndDebug(combustionModel, 0)
void sweep(volScalarField &field, const volScalarField &alpha, const label nLayers, const scalar alphaDiff=0.2)
void smooth(volScalarField &field, const scalar coeff)