void Foam::GAMGPreconditioner::precondition
(
scalarField& wA,
const scalarField& rA,
const direction cmpt
) const
{
wA = 0.0;
scalarField AwA(wA.size());
scalarField finestCorrection(wA.size());
scalarField finestResidual(rA);
// Create coarse grid correction fields
PtrList<scalarField> coarseCorrX;
// Create coarse grid sources
PtrList<scalarField> coarseB;
// Create the smoothers for all levels
PtrList<lduSmoother> smoothers;
// Initialise the above data structures
GAMG_.initVcycle(coarseCorrX, coarseB, smoothers);
for (label cycle=0; cycle<nVcycles_; cycle++)
{
GAMG_.Vcycle
(
smoothers,
wA,
rA,
AwA,
finestCorrection,
finestResidual,
coarseCorrX,
coarseB,
cmpt
);
if (cycle < nVcycles_ - 1)
{
// Calculate finest level residual field
matrix_.Amul(AwA, wA, coupleBouCoeffs_, interfaces_, cmpt);
finestResidual = rA;
finestResidual -= AwA;
}
}
}
void Foam::GAMGPreconditioner::precondition
(
solveScalarField& wA,
const solveScalarField& rA_ss,
const direction cmpt
) const
{
wA = 0.0;
solveScalarField AwA(wA.size());
solveScalarField finestCorrection(wA.size());
solveScalarField finestResidual(rA_ss);
// Create coarse grid correction fields
PtrList<solveScalarField> coarseCorrFields;
// Create coarse grid sources
PtrList<solveScalarField> coarseSources;
// Create the smoothers for all levels
PtrList<lduMatrix::smoother> smoothers;
// Scratch fields if processor-agglomerated coarse level meshes
// are bigger than original. Usually not needed
solveScalarField ApsiScratch;
solveScalarField finestCorrectionScratch;
// Initialise the above data structures
initVcycle
(
coarseCorrFields,
coarseSources,
smoothers,
ApsiScratch,
finestCorrectionScratch
);
scalarField rA_s;
for (label cycle=0; cycle<nVcycles_; cycle++)
{
const scalarField& rA =
ConstFieldWrapper<scalar, solveScalar>::get(rA_ss, rA_s);
Vcycle
(
smoothers,
wA,
rA,
AwA,
finestCorrection,
finestResidual,
(ApsiScratch.size() ? ApsiScratch : AwA),
(
finestCorrectionScratch.size()
? finestCorrectionScratch
: finestCorrection
),
coarseCorrFields,
coarseSources,
cmpt
);
if (cycle < nVcycles_-1)
{
// Calculate finest level residual field
matrix_.Amul(AwA, wA, interfaceBouCoeffs_, interfaces_, cmpt);
finestResidual = rA_ss;
finestResidual -= AwA;
}
}
}
