Public Types | |
typedef PreconditionRelaxation < MATRIX > | BaseClass |
Public Member Functions | |
void | initialize (const MATRIX &A, typename BaseClass::AdditionalData parameters=typename BaseClass::AdditionalData()) |
template<class VECTOR > | |
void | vmult (VECTOR &, const VECTOR &) const |
template<class VECTOR > | |
void | Tvmult (VECTOR &, const VECTOR &) const |
Private Attributes | |
std::vector< unsigned int > | pos_right_of_diagonal |
MATRIX
class used is required to have a function precondition_SSOR(VECTOR&, const VECTOR&, double)
// Declare related objects SparseMatrix<double> A; Vector<double> x; Vector<double> b; SolverCG<> solver(...); //...initialize and build A // Define and initialize preconditioner PreconditionSSOR<SparseMatrix<double> > precondition; precondition.initialize (A, .6); solver.solve (A, x, b, precondition);
typedef PreconditionRelaxation<MATRIX> PreconditionSSOR< MATRIX >::BaseClass |
A typedef to the base class.
void PreconditionSSOR< MATRIX >::initialize | ( | const MATRIX & | A, | |
typename BaseClass::AdditionalData | parameters = typename BaseClass::AdditionalData() | |||
) |
Initialize matrix and relaxation parameter. The matrix is just stored in the preconditioner object. The relaxation parameter should be larger than zero and smaller than 2 for numerical reasons. It defaults to 1.
void PreconditionSSOR< MATRIX >::vmult | ( | VECTOR & | , | |
const VECTOR & | ||||
) | const [inline] |
Apply preconditioner.
void PreconditionSSOR< MATRIX >::Tvmult | ( | VECTOR & | , | |
const VECTOR & | ||||
) | const [inline] |
Apply transpose preconditioner. Since this is a symmetric preconditioner, this function is the same as vmult().
std::vector<unsigned int> PreconditionSSOR< MATRIX >::pos_right_of_diagonal [private] |
An array that stores for each matrix row where the first position after the diagonal is located.