Go to the documentation of this file.
32 template<
class Type,
class DType,
class LUType>
35 const word& fieldName,
51 template<
class Type,
class DType,
class LUType>
55 const word preconditionerName(this->controlDict_.getWord(
"preconditioner"));
60 preconditionerName + typeName,
66 label nCells =
psi.size();
68 Type* __restrict__ psiPtr =
psi.begin();
71 Type* __restrict__ pAPtr = pA.begin();
74 Type* __restrict__ pTPtr = pT.begin();
77 Type* __restrict__ wAPtr = wA.begin();
80 Type* __restrict__ wTPtr = wT.begin();
86 this->matrix_.Amul(wA,
psi);
87 this->matrix_.Tmul(wT,
psi);
92 Type* __restrict__ rAPtr = rA.begin();
93 Type* __restrict__ rTPtr = rT.begin();
96 Type normFactor = this->normFactor(
psi, wA, pA);
100 Info<<
" Normalisation factor = " << normFactor <<
endl;
105 solverPerf.finalResidual() = solverPerf.initialResidual();
108 if (!solverPerf.checkConvergence(this->tolerance_, this->relTol_))
125 preconPtr->precondition(wA, rA);
126 preconPtr->preconditionT(wT, rT);
156 this->matrix_.Amul(wA, pA);
157 this->matrix_.Tmul(wT, pT);
164 solverPerf.checkSingularity
189 solverPerf.finalResidual() =
194 nIter++ < this->maxIter_
195 && !(solverPerf.checkConvergence(this->tolerance_, this->relTol_))
199 solverPerf.nIterations() =
A class for handling words, derived from Foam::string.
dimensioned< Type > cmptMultiply(const dimensioned< Type > &, const dimensioned< Type > &)
static constexpr const zero Zero
Global zero (0)
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
LduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
A class representing the concept of 1 (one) that can be used to avoid manipulating objects known to b...
Ostream & endl(Ostream &os)
Add newline and flush stream.
Abstract base-class for LduMatrix solvers.
Type gSumCmptMag(const UList< Type > &f, const label comm)
void cmptMag(FieldField< Field, Type > &cf, const FieldField< Field, Type > &f)
Generic templated field type.
tmp< DimensionedField< scalar, GeoMesh > > stabilise(const DimensionedField< scalar, GeoMesh > &dsf, const dimensioned< scalar > &ds)
messageStream Info
Information stream (uses stdout - output is on the master only)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
virtual SolverPerformance< Type > solve(Field< Type > &psi) const
Solve the matrix with this solver.
A list of keyword definitions, which are a keyword followed by a number of values (eg,...
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Global function forwards to reuseTmpDimensionedField::New.
dimensioned< Type > cmptDivide(const dimensioned< Type > &, const dimensioned< Type > &)
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
Traits class for primitives.
const volScalarField & psi
A cell is defined as a list of faces with extra functionality.
Type gSumCmptProd(const UList< Type > &f1, const UList< Type > &f2, const label comm)
Preconditioned bi-conjugate gradient solver for asymmetric lduMatrices using a run-time selectable pr...