Reference documentation for deal.II version Git f6c556c4e9 2020-09-24 06:57:19 +0200
\(\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\}}\)
Public Types | Public Member Functions | Static Public Member Functions | Private Attributes | List of all members
PreconditionUseMatrix< MatrixType, VectorType > Class Template Reference

#include <deal.II/lac/precondition.h>

Inheritance diagram for PreconditionUseMatrix< MatrixType, VectorType >:
[legend]

Public Types

using function_ptr = void(MatrixType::*)(VectorType &, const VectorType &) const
 

Public Member Functions

 PreconditionUseMatrix (const MatrixType &M, const function_ptr method)
 
void vmult (VectorType &dst, const VectorType &src) const
 
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)
 

Static Public Member Functions

static ::ExceptionBaseExcInUse (int arg1, std::string arg2, std::string arg3)
 
static ::ExceptionBaseExcNoSubscriber (std::string arg1, std::string arg2)
 

Private Attributes

const MatrixType & matrix
 
const function_ptr precondition
 

Detailed Description

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
class PreconditionUseMatrix< MatrixType, VectorType >

Preconditioner using a matrix-builtin function. This class forms a preconditioner suitable for the LAC solver classes. Since many preconditioning methods are based on matrix entries, these have to be implemented as member functions of the underlying matrix implementation. This class now is intended to allow easy access to these member functions from LAC solver classes.

It seems that all builtin preconditioners have a relaxation parameter, so please use PreconditionRelaxation for these.

You will usually not want to create a named object of this type, although possible. The most common use is like this:

gmres.solve(
matrix, solution, right_hand_side,
matrix, &SparseMatrix<double>::template precondition_Jacobi<double>));

This creates an unnamed object to be passed as the fourth parameter to the solver function of the SolverGMRES class. It assumes that the SparseMatrix class has a function precondition_Jacobi taking two vectors (source and destination) as parameters (Actually, there is no function like that, the existing function takes a third parameter, denoting the relaxation parameter; this example is therefore only meant to illustrate the general idea).

Note that due to the default template parameters, the above example could be written shorter as follows:

...
gmres.solve(
matrix, solution, right_hand_side,
matrix,&SparseMatrix<double>::template precondition_Jacobi<double>));

Definition at line 357 of file precondition.h.

Member Typedef Documentation

◆ function_ptr

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
using PreconditionUseMatrix< MatrixType, VectorType >::function_ptr = void (MatrixType::*)(VectorType &, const VectorType &) const

Type of the preconditioning function of the matrix.

Definition at line 364 of file precondition.h.

Constructor & Destructor Documentation

◆ PreconditionUseMatrix()

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
PreconditionUseMatrix< MatrixType, VectorType >::PreconditionUseMatrix ( const MatrixType &  M,
const function_ptr  method 
)

Constructor. This constructor stores a reference to the matrix object for later use and selects a preconditioning method, which must be a member function of that matrix.

Member Function Documentation

◆ vmult()

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
void PreconditionUseMatrix< MatrixType, VectorType >::vmult ( VectorType dst,
const VectorType src 
) const

Execute preconditioning. Calls the function passed to the constructor of this object with the two arguments given here.

◆ subscribe()

void Subscriptor::subscribe ( std::atomic< bool > *const  validity,
const std::string &  identifier = "" 
) const
inherited

Subscribes a user of the object by storing the pointer validity. The subscriber may be identified by text supplied as identifier.

Definition at line 136 of file subscriptor.cc.

◆ unsubscribe()

void Subscriptor::unsubscribe ( std::atomic< bool > *const  validity,
const std::string &  identifier = "" 
) const
inherited

Unsubscribes a user from the object.

Note
The identifier and the validity pointer must be the same as the one supplied to subscribe().

Definition at line 156 of file subscriptor.cc.

◆ n_subscriptions()

unsigned int Subscriptor::n_subscriptions ( ) const
inlineinherited

Return the present number of subscriptions to this object. This allows to use this class for reference counted lifetime determination where the last one to unsubscribe also deletes the object.

Definition at line 290 of file subscriptor.h.

◆ list_subscribers() [1/2]

template<typename StreamType >
void Subscriptor::list_subscribers ( StreamType &  stream) const
inlineinherited

List the subscribers to the input stream.

Definition at line 307 of file subscriptor.h.

◆ list_subscribers() [2/2]

void Subscriptor::list_subscribers ( ) const
inherited

List the subscribers to deallog.

Definition at line 204 of file subscriptor.cc.

◆ serialize()

template<class Archive >
void Subscriptor::serialize ( Archive &  ar,
const unsigned int  version 
)
inlineinherited

Read or write the data of this object to or from a stream for the purpose of serialization.

This function does not actually serialize any of the member variables of this class. The reason is that what this class stores is only who subscribes to this object, but who does so at the time of storing the contents of this object does not necessarily have anything to do with who subscribes to the object when it is restored. Consequently, we do not want to overwrite the subscribers at the time of restoring, and then there is no reason to write the subscribers out in the first place.

Definition at line 299 of file subscriptor.h.

Member Data Documentation

◆ matrix

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
const MatrixType& PreconditionUseMatrix< MatrixType, VectorType >::matrix
private

Pointer to the matrix in use.

Definition at line 384 of file precondition.h.

◆ precondition

template<typename MatrixType = SparseMatrix<double>, class VectorType = Vector<double>>
const function_ptr PreconditionUseMatrix< MatrixType, VectorType >::precondition
private

Pointer to the preconditioning function.

Definition at line 389 of file precondition.h.


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