Go to the documentation of this file.
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 =
72 matrix_.lduAddr().upperAddr().begin();
73 const label*
const __restrict__ lPtr =
74 matrix_.lduAddr().lowerAddr().begin();
75 const scalar*
const __restrict__ upperPtr =
76 matrix_.upper().begin();
78 const label nFaces = matrix_.upper().size();
81 std::copy(
diag.begin(),
diag.end(), rD_.begin());
83 DICPreconditioner::calcReciprocalD(rD_, matrix_);
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();
132 const label nFaces = matrix_.upper().size();
135 rAPtr[uPtr[
face]] -= rDuUpperPtr[
face]*rAPtr[lPtr[
face]];
138 const label nFacesM1 = nFaces - 1;
141 rAPtr[lPtr[
face]] -= rDlUpperPtr[
face]*rAPtr[uPtr[
face]];
void smooth(solveScalarField &psi, const scalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
A class for handling words, derived from Foam::string.
A field of fields is a PtrList of fields with reference counting.
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.
lduMatrix::smoother::addsymMatrixConstructorToTable< FDICSmoother > addFDICSmootherSymMatrixConstructorToTable_
Abstract base-class for lduMatrix smoothers.
virtual void scalarSmooth(solveScalarField &psi, const solveScalarField &source, const direction cmpt, const label nSweeps) const
Smooth the solution for a given number of sweeps.
FDICSmoother(const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces)
Construct from matrix components.
A const Field/List wrapper with possible data conversion.
A face is a list of labels corresponding to mesh vertices.
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)