Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS > Class Template Reference

Nonoverlapping parallel CG solver with Jacobi preconditioner. More...

#include <dune/pdelab/backend/novlpistlsolverbackend.hh>

Public Member Functions
	ISTLBackend_NOVLP_CG_Jacobi (const GFS &gfs_, unsigned maxiter_=5000, int verbose_=1)
	make a linear solver object More...
template<class V >
V::ElementType	norm (const V &v) const
	compute global norm of a vector More...
template<class M , class V , class W >
void	apply (M &A, V &z, W &r, typename V::ElementType reduction)
	solve the given linear system More...
const LinearSolverResult < double > &	result () const
	Return access to result data. More...

Detailed Description

template<class GFS>
class Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS >

Nonoverlapping parallel CG solver with Jacobi preconditioner.

Constructor & Destructor Documentation

template<class GFS >

Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS >::ISTLBackend_NOVLP_CG_Jacobi ( const GFS &  gfs_, unsigned  maxiter_ = 5000, int  verbose_ = 1  )

inlineexplicit

make a linear solver object

Parameters

gfs_	A grid function space
maxiter_	Maximum number of iterations to do.
verbose_	Verbosity level, directly handed to the CGSolver.

Member Function Documentation

template<class GFS >

template<class M , class V , class W >

void Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS >::apply ( M &  A, V &  z, W &  r, typename V::ElementType  reduction  )

inline

solve the given linear system

Parameters

A	The matrix to solve. Should be a matrix from one of PDELabs ISTL backends (only ISTLBCRSMatrixBackend at the moment).
z	The solution vector to be computed
r	Right hand side
reduction	to be achieved

Solve the linear system A*z=r such that norm(A*z0-r)/norm(A*z-r) < reduction where z0 is the initial value of z.

References Dune::PDELab::LinearSolverResult< RFType >::conv_rate, Dune::PDELab::LinearSolverResult< RFType >::converged, Dune::PDELab::LinearSolverResult< RFType >::elapsed, Dune::PDELab::LinearSolverResult< RFType >::iterations, Dune::PDELab::istl::raw(), and Dune::PDELab::LinearSolverResult< RFType >::reduction.

template<class GFS >

template<class V >

V::ElementType Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS >::norm ( const V &  v ) const

inline

compute global norm of a vector

Parameters

v	The vector to compute the norm of. Should be an inconsistent vector (i.e. the entries corresponding a DoF on the border should only contain the summand of this process).

template<class GFS >

const LinearSolverResult<double>& Dune::PDELab::ISTLBackend_NOVLP_CG_Jacobi< GFS >::result ( ) const

inline

Return access to result data.

---

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

• novlpistlsolverbackend.hh