Reference documentation for deal.II version 9.4.1
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Private Attributes | List of all members
MatrixFreeOperators::MGInterfaceOperator< OperatorType > Class Template Reference

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

Inheritance diagram for MatrixFreeOperators::MGInterfaceOperator< OperatorType >:
[legend]

Public Types

using value_type = typename OperatorType::value_type
 
using size_type = typename OperatorType::size_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
 
template<typename VectorType >
void initialize_dof_vector (VectorType &vec) const
 

Private Attributes

SmartPointer< const OperatorType > mf_base_operator
 

Subscriptor functionality

Classes derived from Subscriptor provide a facility to subscribe to this object. This is mostly used by the SmartPointer class.

std::atomic< unsigned intcounter
 
std::map< std::string, unsigned intcounter_map
 
std::vector< std::atomic< bool > * > validity_pointers
 
const std::type_info * object_info
 
void subscribe (std::atomic< bool > *const validity, const std::string &identifier="") const
 
void unsubscribe (std::atomic< bool > *const validity, const std::string &identifier="") const
 
unsigned int n_subscriptions () const
 
template<typename StreamType >
void list_subscribers (StreamType &stream) const
 
void list_subscribers () const
 
template<class Archive >
void serialize (Archive &ar, const unsigned int version)
 
using map_value_type = decltype(counter_map)::value_type
 
using map_iterator = decltype(counter_map)::iterator
 
static std::mutex mutex
 
static ::ExceptionBaseExcInUse (int arg1, std::string arg2, std::string arg3)
 
static ::ExceptionBaseExcNoSubscriber (std::string arg1, std::string arg2)
 
void check_no_subscribers () const noexcept
 

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.

Definition at line 540 of file operators.h.

Member Typedef Documentation

◆ value_type

template<typename OperatorType >
using MatrixFreeOperators::MGInterfaceOperator< OperatorType >::value_type = typename OperatorType::value_type

Number alias.

Definition at line 546 of file operators.h.

◆ size_type

template<typename OperatorType >
using MatrixFreeOperators::MGInterfaceOperator< OperatorType >::size_type = typename OperatorType::size_type

Size type.

Definition at line 551 of file operators.h.

Constructor & Destructor Documentation

◆ MGInterfaceOperator()

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

Default constructor.

Definition at line 1741 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 1750 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 1759 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 1769 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 1789 of file operators.h.

◆ initialize_dof_vector()

template<typename OperatorType >
template<typename VectorType >
void MatrixFreeOperators::MGInterfaceOperator< OperatorType >::initialize_dof_vector ( VectorType &  vec) const

A wrapper for initialize_dof_vector() of OperatorType object.

Definition at line 1809 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 596 of file operators.h.


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