Go to the documentation of this file.
38 lduMatrix::solver::addsymMatrixConstructorToTable<PCG>
47 const word& fieldName,
79 lduMatrix::preconditioner::getName(controlDict_) + typeName,
83 label nCells =
psi.size();
85 solveScalar* __restrict__ psiPtr =
psi.begin();
88 solveScalar* __restrict__ pAPtr = pA.begin();
91 solveScalar* __restrict__ wAPtr = wA.begin();
93 solveScalar wArA = solverPerf.
great_;
94 solveScalar wArAold = wArA;
97 matrix_.Amul(wA,
psi, interfaceBouCoeffs_, interfaces_, cmpt);
101 solveScalar* __restrict__ rAPtr = rA.begin();
111 solveScalar normFactor = this->normFactor(
psi, source, wA, pA);
115 Info<<
" Normalisation factor = " << normFactor <<
endl;
146 preconPtr->precondition(wA, rA, cmpt);
160 solveScalar
beta = wArA/wArAold;
170 matrix_.Amul(wA, pA, interfaceBouCoeffs_, interfaces_, cmpt);
172 solveScalar wApA =
gSumProd(wA, pA, matrix().
mesh().comm());
180 solveScalar
alpha = wArA/wApA;
int debug
Static debugging option.
typeOfMag< Type >::type gSumMag(const FieldField< Field, Type > &f)
lduMatrix::solver::addsymMatrixConstructorToTable< PCG > addPCGSymMatrixConstructorToTable_
A class for handling words, derived from Foam::string.
A field of fields is a PtrList of fields with reference counting.
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
lduMatrix is a general matrix class in which the coefficients are stored as three arrays,...
Abstract base-class for lduMatrix solvers.
Ostream & endl(Ostream &os)
Add newline and flush stream.
void setResidualField(const scalarField &residual, const word &fieldName, const bool initial) const
virtual solverPerformance scalarSolve(solveScalarField &psi, const solveScalarField &source, const direction cmpt=0) const
Solve the matrix with this solver.
messageStream Info
Information stream (uses stdout - output is on the master only)
virtual solverPerformance solve(scalarField &psi, const scalarField &source, const direction cmpt=0) const
Solve the matrix with this solver.
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
scalarProduct< Type, Type >::type gSumProd(const UList< Type > &f1, const UList< Type > &f2, const label comm)
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.
Pointer management similar to std::unique_ptr, with some additional methods and type checking.
A Field wrapper with possible data conversion.
FieldType & ref()
Allow modification without const-ref check.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A const Field wrapper with possible data conversion.
const volScalarField & psi
defineTypeNameAndDebug(combustionModel, 0)
A cell is defined as a list of faces with extra functionality.