const bool every_iteration=false);
protected:
- /**
- * Implementation of the computation of the norm of the residual. This can
- * be replaced by a more problem oriented functional in a derived class.
- */
- virtual double criterion();
-
/**
* Interface for derived class. This function gets the current iteration
* vector, the residual and the update vector in each step. It can be used
VectorType *Vp;
VectorType *Vz;
- /**
- * Within the iteration loop, the square of the residual vector is stored in
- * this variable. The function @p criterion uses this variable to compute
- * the convergence value, which in this class is the norm of the residual
- * vector and thus the square root of the @p res2 value.
- */
- double res2;
-
/**
* Additional parameters.
*/
-template <typename VectorType>
-double
-SolverCG<VectorType>::criterion()
-{
- return std::sqrt(res2);
-}
-
-
-
template <typename VectorType>
void
SolverCG<VectorType>::cleanup()