MatrixFreeOperators::MGInterfaceOperator< OperatorType > Class Template Reference

#include <deal.II/matrix_free/operators.h>

Inheritance diagram for MatrixFreeOperators::MGInterfaceOperator< OperatorType >:

Public Types
typedef OperatorType::value_type	value_type

Public Member Functions
	MGInterfaceOperator ()
void	clear ()
void	initialize (const OperatorType &operator_in)
template<typename VectorType >
void	vmult (VectorType &dst, const VectorType &src) const
template<typename VectorType >
void	Tvmult (VectorType &dst, const VectorType &src) const
Public Member Functions inherited from Subscriptor
	Subscriptor ()
	Subscriptor (const Subscriptor &)
	Subscriptor (Subscriptor &&)
virtual	~Subscriptor ()
Subscriptor &	operator= (const Subscriptor &)
Subscriptor &	operator= (Subscriptor &&)
void	subscribe (const char *identifier=0) const
void	unsubscribe (const char *identifier=0) const
unsigned int	n_subscriptions () const
void	list_subscribers () const
template<class Archive >
void	serialize (Archive &ar, const unsigned int version)

Private Attributes
SmartPointer< const OperatorType >	mf_base_operator

Additional Inherited Members
Static Public Member Functions inherited from Subscriptor
static ::ExceptionBase &	ExcInUse (int arg1, char *arg2, std::string &arg3)
static ::ExceptionBase &	ExcNoSubscriber (char *arg1, char *arg2)

Detailed Description

template<typename OperatorType>
class MatrixFreeOperators::MGInterfaceOperator< OperatorType >

Auxiliary class to provide interface vmult/Tvmult methods required in adaptive geometric multgrids. OperatorType class should be derived from MatrixFreeOperators::Base class.

The adaptive multigrid realization in deal.II implements an approach called local smoothing. This means that the smoothing on the finest level only covers the local part of the mesh defined by the fixed (finest) grid level and ignores parts of the computational domain where the terminal cells are coarser than this level. As the method progresses to coarser levels, more and more of the global mesh will be covered. At some coarser level, the whole mesh will be covered. Since all level matrices in the multigrid method cover a single level in the mesh, no hanging nodes appear on the level matrices. At the interface between multigrid levels, homogeneous Dirichlet boundary conditions are set while smoothing. When the residual is transferred to the next coarser level, however, the coupling over the multigrid interface needs to be taken into account. This is done by the so-called interface (or edge) matrices that compute the part of the residual that is missed by the level matrix with homogeneous Dirichlet conditions. We refer to the Multigrid paper by Janssen and Kanschat for more details.

For the implementation of those interface matrices, most infrastructure is already in place and provided by MatrixFreeOperators::Base through the two multiplication routines vmult_interface_down() and vmult_interface_up(). The only thing MGInterfaceOperator does is wrapping those operations and make them accessible via vmult() and Tvmult interface as expected by the multigrid routines (that were originally written for matrices, hence expecting those names). Note that the vmult_interface_down is used during the restriction phase of the multigrid V-cycle, whereas vmult_interface_up is used during the prolongation phase.

Author

Martin Kronbichler, 2016

Definition at line 465 of file operators.h.

Member Typedef Documentation

value_type

template<typename OperatorType >

typedef OperatorType::value_type MatrixFreeOperators::MGInterfaceOperator< OperatorType >::value_type

Number typedef.

Definition at line 471 of file operators.h.

Constructor & Destructor Documentation

MGInterfaceOperator()

template<typename OperatorType >

MatrixFreeOperators::MGInterfaceOperator< OperatorType >::MGInterfaceOperator

(

)

Default constructor.

Definition at line 1403 of file operators.h.

Member Function Documentation

clear()

template<typename OperatorType >

void MatrixFreeOperators::MGInterfaceOperator< OperatorType >::clear

(

)

Clear the pointer to the OperatorType object.

Definition at line 1414 of file operators.h.

initialize()

template<typename OperatorType >

void MatrixFreeOperators::MGInterfaceOperator< OperatorType >::initialize

(

const OperatorType &

operator_in

)

Initialize this class with an operator operator_in.

Definition at line 1423 of file operators.h.

vmult()

template<typename OperatorType >

template<typename VectorType >

void MatrixFreeOperators::MGInterfaceOperator< OperatorType >::vmult	(	VectorType &	dst,
		const VectorType &	src
	)		const

vmult operator, see class description for more info.

Definition at line 1433 of file operators.h.

Tvmult()

template<typename OperatorType >

template<typename VectorType >

void MatrixFreeOperators::MGInterfaceOperator< OperatorType >::Tvmult	(	VectorType &	dst,
		const VectorType &	src
	)		const

Tvmult operator, see class description for more info.

Definition at line 1455 of file operators.h.

Member Data Documentation

mf_base_operator

template<typename OperatorType >

SmartPointer<const OperatorType> MatrixFreeOperators::MGInterfaceOperator< OperatorType >::mf_base_operator

private

Const pointer to the operator class.

Definition at line 506 of file operators.h.

---

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

• deal.II/matrix_free/operators.h