Default iterative solver for boundary integral equations. More...

#include </home/wojtek/Projects/BEM/bempp-sven/bempp/lib/linalg/default_iterative_solver.hpp>

Inheritance diagram for Bempp::DefaultIterativeSolver< BasisFunctionType, ResultType >:

Public Types
typedef Solver < BasisFunctionType, ResultType >	Base

Public Member Functions
	DefaultIterativeSolver (const BoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, ConvergenceTestMode::Mode mode=ConvergenceTestMode::TEST_CONVERGENCE_IN_DUAL_TO_RANGE)
	Constructor of the `DefaultIterativeSolver` class. More...

	DefaultIterativeSolver (const BlockedBoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, ConvergenceTestMode::Mode mode=ConvergenceTestMode::TEST_CONVERGENCE_IN_DUAL_TO_RANGE)
	Constructor of the `DefaultIterativeSolver` class. More...

BEMPP_DEPRECATED void	setPreconditioner (const Preconditioner< ResultType > &preconditioner)
	Define a preconditioner. More...

void	initializeSolver (const Teuchos::RCP< Teuchos::ParameterList > &paramList)
	Initialize a Belos iterative solver. More...

void	initializeSolver (const Teuchos::RCP< Teuchos::ParameterList > &paramList, const Preconditioner< ResultType > &preconditioner)
	Initialize a preconditioned Belos iterative solver. More...

Public Member Functions inherited from Bempp::Solver< BasisFunctionType, ResultType >
Solution< BasisFunctionType, ResultType >	solve (const GridFunction< BasisFunctionType, ResultType > &rhs) const
	Solve a standard (non-blocked) boundary integral equation. More...

BlockedSolution < BasisFunctionType, ResultType >	solve (const std::vector< GridFunction< BasisFunctionType, ResultType > > &rhs) const
	Solve a block-operator system of boundary integral equations. More...

Private Member Functions
virtual Solution < BasisFunctionType, ResultType >	solveImplNonblocked (const GridFunction< BasisFunctionType, ResultType > &rhs) const

virtual BlockedSolution < BasisFunctionType, ResultType >	solveImplBlocked (const std::vector< GridFunction< BasisFunctionType, ResultType > > &rhs) const

Private Attributes
boost::scoped_ptr< Impl >	m_impl

Additional Inherited Members
Static Protected Member Functions inherited from Bempp::Solver< BasisFunctionType, ResultType >
static void	checkConsistency (const BoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, const GridFunction< BasisFunctionType, ResultType > &rhs, ConvergenceTestMode::Mode mode)

static void	checkConsistency (const BlockedBoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, const std::vector< GridFunction< BasisFunctionType, ResultType > > &rhs, ConvergenceTestMode::Mode mode)

static std::vector < GridFunction < BasisFunctionType, ResultType > >	canonicalizeBlockedRhs (const BlockedBoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, const std::vector< GridFunction< BasisFunctionType, ResultType > > &rhs, ConvergenceTestMode::Mode mode)

static void	constructBlockedGridFunction (const arma::Col< ResultType > &solution, const BlockedBoundaryOperator< BasisFunctionType, ResultType > &boundaryOp, std::vector< GridFunction< BasisFunctionType, ResultType > > &solutionFunctions)

Detailed Description

template<typename BasisFunctionType, typename ResultType>
class Bempp::DefaultIterativeSolver< BasisFunctionType, ResultType >

Default iterative solver for boundary integral equations.

This class provides an interface to various iterative solvers available via the Stratimikos interface to the Belos solver family from Trilinos (see Stratimikos documentation). Convergence can be tested either in range space or in the dual space to the range space. A standard Galerkin discretisation of the form $Ax=b$ , maps into the dual space of the range of the operator. By choosing to test in the range space the equation $M^\dagger Ax=M^\dagger b$ is solved, where $M$ is the mass matrix, mapping from the range space into its dual and $M^\dagger$ is its pseudoinverse.