Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm > Class Template Reference

#include <dune/pdelab/backend/istl/ovlp_amg_dg_backend.hh>

Public Types
typedef Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type >	domain_type
	The domain type of the preconditioner.
typedef Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type >	range_type
	The range type of the preconditioner.
typedef X::field_type	field_type
	The field type of the preconditioner.

Public Member Functions
SolverCategory::Category	category () const override
	Category of the preconditioner (see SolverCategory::Category)
	OvlpDGAMGPrec (const DGGFS &dggfs_, DGMatrix &dgmatrix_, DGPrec &dgprec_, const DGCC &dgcc_, const CGGFS &cggfs_, CGPrec &cgprec_, const CGCC &cgcc_, P &p_, const DGHelper &dghelper_, const Comm &comm_, int n1_, int n2_)
	Constructor. More...
virtual void	pre (V &x, W &b) override
	Prepare the preconditioner. More...
virtual void	apply (V &x, const W &b) override
	Apply the precondioner. More...
virtual void	post (V &x) override
	Clean up. More...
virtual void	pre (Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &x, Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > &b)=0
	Prepare the preconditioner. More...
virtual void	apply (Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &v, const Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > &d)=0
	Apply one step of the preconditioner to the system A(v)=d. More...
virtual void	post (Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &x)=0
	Clean up. More...

Detailed Description

template<class DGGFS, class DGMatrix, class DGPrec, class DGCC, class CGGFS, class CGPrec, class CGCC, class P, class DGHelper, class Comm>
class Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm >

An ISTL preconditioner for DG based on AMG applied to CG subspace

The template parameters are: DGGFS DG space DGMatrix BCRSMatrix assembled with DG DGPrec preconditioner to be used for DG CGPrec preconditioner to be used on CG subspace P BCRSMatrix for grid transfer

Constructor & Destructor Documentation

OvlpDGAMGPrec()

template<class DGGFS , class DGMatrix , class DGPrec , class DGCC , class CGGFS , class CGPrec , class CGCC , class P , class DGHelper , class Comm >

Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm >::OvlpDGAMGPrec ( const DGGFS &  dggfs_, DGMatrix &  dgmatrix_, DGPrec &  dgprec_, const DGCC &  dgcc_, const CGGFS &  cggfs_, CGPrec &  cgprec_, const CGCC &  cgcc_, P &  p_, const DGHelper &  dghelper_, const Comm &  comm_, int  n1_, int  n2_  )

inline

Constructor.

Constructor gets all parameters to operate the prec.

Parameters

A	The matrix to operate on.
n	The number of iterations to perform.
w	The relaxation factor.

Member Function Documentation

apply() [1/2]

virtual void Dune::Preconditioner< Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > , Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > >::apply ( Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &  v, const Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > &  d  )

pure virtualinherited

Apply one step of the preconditioner to the system A(v)=d.

On entry v=0 and d=b-A(x) (although this might not be computed in that way. On exit v contains the update, i.e one step computes \( v = M^{-1} d \) where \( M \) is the approximate inverse of the operator \( A \) characterizing the preconditioner.

Parameters

[out]	v	The update to be computed
	d	The current defect.

apply() [2/2]

template<class DGGFS , class DGMatrix , class DGPrec , class DGCC , class CGGFS , class CGPrec , class CGCC , class P , class DGHelper , class Comm >

virtual void Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm >::apply ( V &  x, const W &  b  )

inlineoverridevirtual

Apply the precondioner.

Apply one step of the preconditioner to the system A(v)=d.

On entry v=0 and d=b-A(x) (although this might not be computed in that way. On exit v contains the update, i.e one step computes \( v = M^{-1} d \) where \( M \) is the approximate inverse of the operator \( A \) characterizing the preconditioner.

Parameters

[out]	v	The update to be computed
	d	The current defect.

References Dune::All_All_Interface, Dune::ForwardCommunication, and Dune::PDELab::set_constrained_dofs().

post() [1/2]

virtual void Dune::Preconditioner< Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > , Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > >::post ( Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &  x )

pure virtualinherited

Clean up.

This method is called after the last apply call for the linear system to be solved. Memory may be deallocated safely here. x is the solution of the linear equation.

Parameters

x	The right hand side of the equation.

post() [2/2]

template<class DGGFS , class DGMatrix , class DGPrec , class DGCC , class CGGFS , class CGPrec , class CGCC , class P , class DGHelper , class Comm >

virtual void Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm >::post ( V &  x )

inlineoverridevirtual

Clean up.

Clean up.

This method is called after the last apply call for the linear system to be solved. Memory may be deallocated safely here. x is the solution of the linear equation.

Parameters

x	The right hand side of the equation.

pre() [1/2]

virtual void Dune::Preconditioner< Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > , Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > >::pre ( Dune::PDELab::Backend::Vector< DGGFS, DGPrec::domain_type::field_type > &  x, Dune::PDELab::Backend::Vector< DGGFS, DGPrec::range_type::field_type > &  b  )

pure virtualinherited

Prepare the preconditioner.

A solver solves a linear operator equation A(x)=b by applying one or several steps of the preconditioner. The method pre() is called before the first apply operation. b and x are right hand side and solution vector of the linear system respectively. It may. e.g., scale the system, allocate memory or compute a (I)LU decomposition. Note: The ILU decomposition could also be computed in the constructor or with a separate method of the derived method if several linear systems with the same matrix are to be solved.

Note

if a preconditioner is copied (e.g. for a second thread) again the pre() method has to be called to ensure proper memory management.

X x(0.0);
Y b = ...; // rhs
Preconditioner<X,Y> prec(...);
prec.pre(x,b);   // prepare the preconditioner
prec.apply(x,b); // can be called multiple times now...
prec.post(x);    // cleanup internal state

Parameters

x	The left hand side of the equation.
b	The right hand side of the equation.

pre() [2/2]

template<class DGGFS , class DGMatrix , class DGPrec , class DGCC , class CGGFS , class CGPrec , class CGCC , class P , class DGHelper , class Comm >

virtual void Dune::PDELab::OvlpDGAMGPrec< DGGFS, DGMatrix, DGPrec, DGCC, CGGFS, CGPrec, CGCC, P, DGHelper, Comm >::pre ( V &  x, W &  b  )

inlineoverridevirtual

Prepare the preconditioner.

Prepare the preconditioner.

A solver solves a linear operator equation A(x)=b by applying one or several steps of the preconditioner. The method pre() is called before the first apply operation. b and x are right hand side and solution vector of the linear system respectively. It may. e.g., scale the system, allocate memory or compute a (I)LU decomposition. Note: The ILU decomposition could also be computed in the constructor or with a separate method of the derived method if several linear systems with the same matrix are to be solved.

Note

if a preconditioner is copied (e.g. for a second thread) again the pre() method has to be called to ensure proper memory management.

X x(0.0);
Y b = ...; // rhs
Preconditioner<X,Y> prec(...);
prec.pre(x,b);   // prepare the preconditioner
prec.apply(x,b); // can be called multiple times now...
prec.post(x);    // cleanup internal state

Parameters

x	The left hand side of the equation.
b	The right hand side of the equation.

---

The documentation for this class was generated from the following file:

• dune/pdelab/backend/istl/ovlp_amg_dg_backend.hh