Reference documentation for deal.II version 8.1.0
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
Classes | Public Member Functions | Private Attributes | List of all members
SolverFGMRES< VECTOR > Class Template Reference

#include <solver_gmres.h>

Inheritance diagram for SolverFGMRES< VECTOR >:
[legend]

Classes

struct  AdditionalData
 

Public Member Functions

 SolverFGMRES (SolverControl &cn, VectorMemory< VECTOR > &mem, const AdditionalData &data=AdditionalData())
 
 SolverFGMRES (SolverControl &cn, const AdditionalData &data=AdditionalData())
 
template<class MATRIX , class PRECONDITIONER >
void solve (const MATRIX &A, VECTOR &x, const VECTOR &b, const PRECONDITIONER &precondition)
 
- Public Member Functions inherited from Solver< VECTOR >
 Solver (SolverControl &solver_control, VectorMemory< VECTOR > &vector_memory)
 
 Solver (SolverControl &solver_control)
 
SolverControlcontrol () const
 
- Public Member Functions inherited from Subscriptor
 Subscriptor ()
 
 Subscriptor (const Subscriptor &)
 
virtual ~Subscriptor ()
 
Subscriptoroperator= (const Subscriptor &)
 
void subscribe (const char *identifier=0) const
 
void unsubscribe (const char *identifier=0) const
 
unsigned int n_subscriptions () const
 
void list_subscribers () const
 
 DeclException3 (ExcInUse, int, char *, std::string &,<< "Object of class "<< arg2<< " is still used by "<< arg1<< " other objects.\n"<< "(Additional information: "<< arg3<< ")\n"<< "Note the entry in the Frequently Asked Questions of "<< "deal.II (linked to from http://www.dealii.org/) for "<< "more information on what this error means.")
 
 DeclException2 (ExcNoSubscriber, char *, char *,<< "No subscriber with identifier \""<< arg2<< "\" did subscribe to this object of class "<< arg1)
 
template<class Archive >
void serialize (Archive &ar, const unsigned int version)
 

Private Attributes

AdditionalData additional_data
 
FullMatrix< doubleH
 
FullMatrix< doubleH1
 

Additional Inherited Members

- Protected Attributes inherited from Solver< VECTOR >
GrowingVectorMemory< VECTOR > static_vector_memory
 
SolverControlcntrl
 
VectorMemory< VECTOR > & memory
 

Detailed Description

template<class VECTOR = Vector<double>>
class SolverFGMRES< VECTOR >

Implementation of the Generalized minimal residual method with flexible preconditioning method.

This version of the GMRES method allows for the use of a different preconditioner in each iteration step. Therefore, it is also more robust with respect to inaccurate evaluation of the preconditioner. An important application is also the use of a Krylov space method inside the preconditioner.

FGMRES needs two vectors in each iteration steps yielding a total of 2 * SolverFGMRESAdditionalData::max_basis_size+1 auxiliary vectors.

Caveat: documentation of this class is not up to date. There are also a few parameters of GMRES we would like to introduce here.

Author
Guido Kanschat, 2003

Definition at line 309 of file solver_gmres.h.

Constructor & Destructor Documentation

template<class VECTOR = Vector<double>>
SolverFGMRES< VECTOR >::SolverFGMRES ( SolverControl cn,
VectorMemory< VECTOR > &  mem,
const AdditionalData data = AdditionalData() 
)

Constructor.

template<class VECTOR = Vector<double>>
SolverFGMRES< VECTOR >::SolverFGMRES ( SolverControl cn,
const AdditionalData data = AdditionalData() 
)

Constructor. Use an object of type GrowingVectorMemory as a default to allocate memory.

Member Function Documentation

template<class VECTOR = Vector<double>>
template<class MATRIX , class PRECONDITIONER >
void SolverFGMRES< VECTOR >::solve ( const MATRIX A,
VECTOR &  x,
const VECTOR &  b,
const PRECONDITIONER &  precondition 
)

Solve the linear system $Ax=b$ for x.

Member Data Documentation

template<class VECTOR = Vector<double>>
AdditionalData SolverFGMRES< VECTOR >::additional_data
private

Additional flags.

Definition at line 362 of file solver_gmres.h.

template<class VECTOR = Vector<double>>
FullMatrix<double> SolverFGMRES< VECTOR >::H
private

Projected system matrix

Definition at line 366 of file solver_gmres.h.

template<class VECTOR = Vector<double>>
FullMatrix<double> SolverFGMRES< VECTOR >::H1
private

Auxiliary matrix for inverting H

Definition at line 370 of file solver_gmres.h.


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