Loading...
Searching...
No Matches
PPCG Class Reference

Preconditioned pipelined conjugate gradient solver for symmetric lduMatrices using a run-time selectable preconditioner. More...

#include <PPCG.H>

Inheritance diagram for PPCG:
Inheritance graph
[legend]
Collaboration diagram for PPCG:
Collaboration graph
[legend]

Public Member Functions

 TypeName ("PPCG")
 Runtime type information.
 PPCG (const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls)
 Construct from matrix components and solver controls.
virtual ~PPCG ()=default
 Destructor.
virtual solverPerformance solve (scalarField &psi, const scalarField &source, const direction cmpt=0) const
 Solve the matrix with this solver.
Public Member Functions inherited from lduMatrix::solver
virtual const word & type () const =0
 Runtime type information.
 declareRunTimeSelectionTable (autoPtr, solver, symMatrix,(const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls),(fieldName, matrix, interfaceBouCoeffs, interfaceIntCoeffs, interfaces, solverControls))
 declareRunTimeSelectionTable (autoPtr, solver, asymMatrix,(const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls),(fieldName, matrix, interfaceBouCoeffs, interfaceIntCoeffs, interfaces, solverControls))
 solver (const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls)
 Construct solver for given field name, matrix etc.
virtual ~solver ()=default
 Destructor.
const word & fieldName () const noexcept
const lduMatrix & matrix () const noexcept
const FieldField< Field, scalar > & interfaceBouCoeffs () const noexcept
const FieldField< Field, scalar > & interfaceIntCoeffs () const noexcept
const lduInterfaceFieldPtrsList & interfaces () const noexcept
virtual void read (const dictionary &)
 Read and reset the solver parameters from the given stream.
virtual solverPerformance scalarSolve (solveScalarField &psi, const solveScalarField &source, const direction cmpt=0) const
 Solve with given field and rhs (in solveScalar precision).
solveScalarField::cmptType normFactor (const solveScalarField &psi, const solveScalarField &source, const solveScalarField &Apsi, solveScalarField &tmpField, const lduMatrix::normTypes normType) const
 Return the matrix norm using the specified norm method.
solveScalarField::cmptType normFactor (const solveScalarField &psi, const solveScalarField &source, const solveScalarField &Apsi, solveScalarField &tmpField) const
 Return the matrix norm used to normalise the residual for the stopping criterion.

Protected Member Functions

solverPerformance scalarSolveCG (solveScalarField &psi, const solveScalarField &source, const direction cmpt, const bool cgMode) const
 CG solver. Operates either in conjugate-gradient mode or conjugate residual.
Protected Member Functions inherited from lduMatrix::solver
virtual void readControls ()
 Read the control parameters from controlDict_.

Additional Inherited Members

Static Public Member Functions inherited from lduMatrix::solver
static autoPtr< solver > New (const word &solverName, const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls)
 Return a new solver of given type.
static autoPtr< solver > New (const word &fieldName, const lduMatrix &matrix, const FieldField< Field, scalar > &interfaceBouCoeffs, const FieldField< Field, scalar > &interfaceIntCoeffs, const lduInterfaceFieldPtrsList &interfaces, const dictionary &solverControls)
 Return a new solver given dictionary.
Protected Attributes inherited from lduMatrix::solver
word fieldName_
const lduMatrix & matrix_
const FieldField< Field, scalar > & interfaceBouCoeffs_
const FieldField< Field, scalar > & interfaceIntCoeffs_
lduInterfaceFieldPtrsList interfaces_
dictionary controlDict_
 Dictionary of solution controls.
int log_
 Verbosity level for solver output statements.
label minIter_
 Minimum number of iterations in the solver.
label maxIter_
 Maximum number of iterations in the solver.
lduMatrix::normTypes normType_
 The normalisation type.
scalar tolerance_
 Final convergence tolerance.
scalar relTol_
 Convergence tolerance relative to the initial.
profilingTrigger profiling_
 Profiling instrumentation.

Detailed Description

Preconditioned pipelined conjugate gradient solver for symmetric lduMatrices using a run-time selectable preconditioner.

Reference: 

    P. Ghysels, W. Vanroose.
    "Hiding global synchronization latency in the
     preconditioned Conjugate Gradient algorithm"

and implementation details from 

    Paul Eller, William Gropp
    "Scalable Non-blocking Preconditioned Conjugate Gradient Methods"
Source files

Definition at line 60 of file PPCG.H.

Constructor & Destructor Documentation

◆ PPCG()

PPCG ( const word & fieldName,
const lduMatrix & matrix,
const FieldField< Field, scalar > & interfaceBouCoeffs,
const FieldField< Field, scalar > & interfaceIntCoeffs,
const lduInterfaceFieldPtrsList & interfaces,
const dictionary & solverControls )

Construct from matrix components and solver controls.

Definition at line 264 of file PPCG.C.

References lduMatrix::solver::fieldName(), lduMatrix::solver::interfaceBouCoeffs(), lduMatrix::solver::interfaceIntCoeffs(), lduMatrix::solver::interfaces(), lduMatrix::solver::matrix(), and lduMatrix::solver::solver().

Here is the call graph for this function:

◆ ~PPCG()

virtual ~PPCG ( )
virtualdefault

Destructor.

References psi.

Member Function Documentation

◆ scalarSolveCG()

Foam::solverPerformance scalarSolveCG ( solveScalarField & psi,
const solveScalarField & source,
const direction cmpt,
const bool cgMode ) const
protected

CG solver. Operates either in conjugate-gradient mode or conjugate residual.

Definition at line 74 of file PPCG.C.

References alpha, beta(), SolverPerformance< Type >::checkConvergence(), lduMatrix::solver::controlDict_, delta, Foam::endl(), lduMatrix::solver::fieldName_, SolverPerformance< Type >::finalResidual(), gamma, lduMatrix::preconditioner::getName(), Foam::Info, SolverPerformance< Type >::initialResidual(), lduMatrix::solver::interfaceBouCoeffs_, lduMatrix::solver::interfaces_, lduMatrix::solver::log_, lduMatrix::solver::matrix(), lduMatrix::solver::matrix_, lduMatrix::solver::maxIter_, lduMatrix::solver::minIter_, n, lduMatrix::preconditioner::New(), SolverPerformance< Type >::nIterations(), lduMatrix::solver::normFactor(), p, psi, lduMatrix::solver::relTol_, s(), lduMatrix::solver::tolerance_, Foam::type(), and UPstream::Request::wait().

Referenced by solve(), and PPCR::solve().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ TypeName()

TypeName ( "PPCG" )

Runtime type information.

References lduMatrix::solver::fieldName(), lduMatrix::solver::interfaceBouCoeffs(), lduMatrix::solver::interfaceIntCoeffs(), lduMatrix::solver::interfaces(), and lduMatrix::solver::matrix().

Here is the call graph for this function:

◆ solve()

Foam::solverPerformance solve ( scalarField & psi,
const scalarField & source,
const direction cmpt = 0 ) const
virtual

Solve the matrix with this solver.

Implements lduMatrix::solver.

Reimplemented in PPCR.

Definition at line 288 of file PPCG.C.

References refPtr< T >::ref(), and scalarSolveCG().

Here is the call graph for this function:
The documentation for this class was generated from the following files:
  • src/OpenFOAM/matrices/lduMatrix/solvers/PPCG/PPCG.H
  • src/OpenFOAM/matrices/lduMatrix/solvers/PPCG/PPCG.C