dune-pdelab 2.7-git
|
Nonoverlapping parallel CG solver with Jacobi preconditioner. More...
#include <dune/pdelab/backend/istl/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 Dune::template FieldTraits< typename V::ElementType >::real_type reduction) |
solve the given linear system More... | |
const LinearSolverResult< double > & | result () const |
Return access to result data. More... | |
Nonoverlapping parallel CG solver with Jacobi preconditioner.
|
inlineexplicit |
make a linear solver object
gfs_ | A grid function space |
maxiter_ | Maximum number of iterations to do. |
verbose_ | Verbosity level, directly handed to the CGSolver. |
|
inline |
solve the given linear system
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.
|
inline |
compute global norm of a vector
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). |
|
inline |
Return access to result data.