dune-istl 2.9.0
|
Sequential overlapping Schwarz preconditioner. More...
#include <dune/istl/overlappingschwarz.hh>
Public Types | |
typedef M | matrix_type |
The type of the matrix to precondition. More... | |
typedef X | domain_type |
The domain type of the preconditioner. More... | |
typedef X | range_type |
The range type of the preconditioner. More... | |
typedef TM | Mode |
The mode (additive or multiplicative) of the Schwarz method. More... | |
typedef X::field_type | field_type |
The field type of the preconditioner. More... | |
typedef matrix_type::size_type | size_type |
The return type of the size method. More... | |
typedef TA | allocator |
The allocator to use. More... | |
typedef std::set< size_type, std::less< size_type >, typename std::allocator_traits< TA >::template rebind_alloc< size_type > > | subdomain_type |
The type for the subdomain to row index mapping. More... | |
typedef std::vector< subdomain_type, typename std::allocator_traits< TA >::template rebind_alloc< subdomain_type > > | subdomain_vector |
The vector type containing the subdomain to row index mapping. More... | |
typedef SLList< size_type, typename std::allocator_traits< TA >::template rebind_alloc< size_type > > | subdomain_list |
The type for the row to subdomain mapping. More... | |
typedef std::vector< subdomain_list, typename std::allocator_traits< TA >::template rebind_alloc< subdomain_list > > | rowtodomain_vector |
The vector type containing the row index to subdomain mapping. More... | |
typedef TD | slu |
The type for the subdomain solver in use. More... | |
typedef std::vector< slu, typename std::allocator_traits< TA >::template rebind_alloc< slu > > | slu_vector |
The vector type containing subdomain solvers. More... | |
Public Member Functions | |
SeqOverlappingSchwarz (const matrix_type &mat, const subdomain_vector &subDomains, field_type relaxationFactor=1, bool onTheFly_=true) | |
Construct the overlapping Schwarz method. More... | |
SeqOverlappingSchwarz (const matrix_type &mat, const rowtodomain_vector &rowToDomain, field_type relaxationFactor=1, bool onTheFly_=true) | |
virtual void | pre (X &x, X &b) |
Prepare the preconditioner. More... | |
virtual void | apply (X &v, const X &d) |
Apply the precondtioner. More... | |
virtual void | post (X &x) |
Postprocess the preconditioner. More... | |
template<bool forward> | |
void | apply (X &v, const X &d) |
Apply one step of the preconditioner to the system A(v)=d. More... | |
virtual SolverCategory::Category | category () const |
Category of the preconditioner (see SolverCategory::Category) More... | |
Sequential overlapping Schwarz preconditioner.
M | The matrix type. |
X | The range and domain type. |
TM | The Schwarz mode. Currently supported modes are AdditiveSchwarzMode, MultiplicativeSchwarzMode, and SymmetricMultiplicativeSchwarzMode. (Default values is AdditiveSchwarzMode) |
TD | The type of the local subdomain solver to be used. |
TA | The type of the allocator to use. |
typedef TA Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::allocator |
The allocator to use.
typedef X Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::domain_type |
The domain type of the preconditioner.
typedef X::field_type Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::field_type |
The field type of the preconditioner.
typedef M Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::matrix_type |
The type of the matrix to precondition.
typedef TM Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::Mode |
The mode (additive or multiplicative) of the Schwarz method.
Either AdditiveSchwarzMode or MultiplicativeSchwarzMode
typedef X Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::range_type |
The range type of the preconditioner.
typedef std::vector<subdomain_list, typename std::allocator_traits<TA>::template rebind_alloc<subdomain_list> > Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::rowtodomain_vector |
The vector type containing the row index to subdomain mapping.
typedef matrix_type::size_type Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::size_type |
The return type of the size method.
typedef TD Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::slu |
The type for the subdomain solver in use.
typedef std::vector<slu, typename std::allocator_traits<TA>::template rebind_alloc<slu> > Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::slu_vector |
The vector type containing subdomain solvers.
typedef SLList<size_type, typename std::allocator_traits<TA>::template rebind_alloc<size_type> > Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::subdomain_list |
The type for the row to subdomain mapping.
typedef std::set<size_type, std::less<size_type>, typename std::allocator_traits<TA>::template rebind_alloc<size_type> > Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::subdomain_type |
The type for the subdomain to row index mapping.
typedef std::vector<subdomain_type, typename std::allocator_traits<TA>::template rebind_alloc<subdomain_type> > Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >::subdomain_vector |
The vector type containing the subdomain to row index mapping.
|
virtual |
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 where
is the approximate inverse of the operator
characterizing the preconditioner.
[out] | v | The update to be computed |
d | The current defect. |
Implements Dune::Preconditioner< X, X >.
|
inlinevirtual |
Category of the preconditioner (see SolverCategory::Category)
Implements Dune::Preconditioner< X, X >.
|
inlinevirtual |
Postprocess the preconditioner.
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.
x | The right hand side of the equation. |
Implements Dune::Preconditioner< X, X >.
|
inlinevirtual |
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.
x | The left hand side of the equation. |
b | The right hand side of the equation. |
Implements Dune::Preconditioner< X, X >.