39 lduMatrix::solver::addsymMatrixConstructorToTable<FPCG>
48 const word& fieldName,
70 void Foam::FPCG::gSumMagProd
72 FixedList<solveScalar, 2>& globalSum,
78 const label nCells = a.
size();
81 for (label cell=0; cell<nCells; ++cell)
83 globalSum[0] += a[cell]*
b[cell];
84 globalSum[1] +=
mag(
b[cell]);
115 solveScalar* __restrict__ psiPtr =
psi.
begin();
118 solveScalar* __restrict__ pAPtr = pA.
begin();
121 solveScalar* __restrict__ wAPtr = wA.
begin();
123 solveScalar wArA = solverPerf.
great_;
124 solveScalar wArAold = wArA;
127 matrix_.Amul(wA,
psi, interfaceBouCoeffs_, interfaces_, cmpt);
131 solveScalar* __restrict__ rAPtr = rA.
begin();
141 solveScalar normFactor = this->normFactor(
psi, source, wA, pA);
145 Info<<
" Normalisation factor = " << normFactor <<
endl;
171 FixedList<solveScalar, 2> globalSum;
180 preconPtr_->precondition(wA, rA, cmpt);
183 gSumMagProd(globalSum, wA, rA, matrix().
mesh().comm());
192 (minIter_ <= 0 || solverPerf.
nIterations() >= minIter_)
201 for (label cell=0; cell<nCells; cell++)
203 pAPtr[cell] = wAPtr[cell];
208 const solveScalar
beta = wArA/wArAold;
210 for (label cell=0; cell<nCells; cell++)
212 pAPtr[cell] = wAPtr[cell] +
beta*pAPtr[cell];
218 matrix_.Amul(wA, pA, interfaceBouCoeffs_, interfaces_, cmpt);
220 solveScalar wApA =
gSumProd(wA, pA, matrix().
mesh().comm());
228 const solveScalar
alpha = wArA/wApA;
230 for (label cell=0; cell<nCells; cell++)
232 psiPtr[cell] +=
alpha*pAPtr[cell];
233 rAPtr[cell] -=
alpha*wAPtr[cell];
244 preconPtr_->setFinished(solverPerf);
249 ConstPrecisionAdaptor<scalar, solveScalar>(rA)(),
void size(const label n)
Older name for setAddressableSize.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
virtual solverPerformance scalarSolve(solveScalarField &psi, const solveScalarField &source, const direction cmpt=0) const
Solve the matrix with this solver.
Field< solveScalar > solveScalarField
A list of keyword definitions, which are a keyword followed by a number of values (eg...
void setResidualField(const scalarField &residual, const word &fieldName, const bool initial) const
Set the residual field using an IOField on the object registry if it exists.
scalarProduct< Type, Type >::type gSumProd(const UList< Type > &f1, const UList< Type > &f2, const label comm)
Ostream & endl(Ostream &os)
Add newline and flush stream.
virtual solverPerformance solve(scalarField &psi, const scalarField &source, const direction cmpt=0) const
Solve the matrix with this solver.
static int & msgType() noexcept
Message tag of standard messages.
A non-const Field/List wrapper with possible data conversion.
A field of fields is a PtrList of fields with reference counting.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
A class for handling words, derived from Foam::string.
static autoPtr< preconditioner > New(const solver &sol, const dictionary &solverControls)
Return a new preconditioner.
A const Field/List wrapper with possible data conversion.
iterator begin() noexcept
Return an iterator to begin traversing the UList.
int debug
Static debugging option.
Container< Type > & ref() const
Return non-const reference to the contents of a non-null managed pointer.
defineTypeNameAndDebug(combustionModel, 0)
typeOfMag< Type >::type gSumMag(const FieldField< Field, Type > &f)
lduMatrix is a general matrix class in which the coefficients are stored as three arrays...
lduMatrix::solver::addsymMatrixConstructorToTable< FPCG > addFPCGSymMatrixConstructorToTable_
messageStream Info
Information stream (stdout output on master, null elsewhere)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
A 'faster' preconditioned conjugate gradient solver for symmetric lduMatrices using a run-time select...
static word getName(const dictionary &)
Find the preconditioner name (directly or from a sub-dictionary)
const volScalarField & psi
const dimensionedScalar alpha
Fine-structure constant: default SI units: [].
void reduce(T &value, const BinaryOp &bop, const int tag=UPstream::msgType(), const label comm=UPstream::worldComm)
Reduce inplace (cf. MPI Allreduce) using linear/tree communication schedule.