solver.hh

Go to the documentation of this file.

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
 
#ifndef DUNE_ISTL_SOLVER_HH
#define DUNE_ISTL_SOLVER_HH
 
#include <iomanip>
#include <ostream>
 
namespace Dune
{
  struct InverseOperatorResult
  {
    InverseOperatorResult ()
    {
      clear();
    }
 
    void clear ()
    {
      iterations = 0;
      reduction = 0;
      converged = false;
      conv_rate = 1;
      elapsed = 0;
    }
 
    int iterations;
 
    double reduction;
 
    bool converged;
 
    double conv_rate;
 
    double elapsed;
  };
 
 
  //=====================================================================
  template<class X, class Y>
  class InverseOperator {
  public:
    typedef X domain_type;
 
    typedef Y range_type;
 
    typedef typename X::field_type field_type;
 
    virtual void apply (X& x, Y& b, InverseOperatorResult& res) = 0;
 
    virtual void apply (X& x, Y& b, double reduction, InverseOperatorResult& res) = 0;
 
    virtual ~InverseOperator () {}
 
  protected:
    // spacing values
    enum { iterationSpacing = 5 , normSpacing = 16 };
 
    void printHeader(std::ostream& s) const
    {
      s << std::setw(iterationSpacing)  << " Iter";
      s << std::setw(normSpacing) << "Defect";
      s << std::setw(normSpacing) << "Rate" << std::endl;
    }
 
    template <class DataType>
    void printOutput(std::ostream& s,
                     const double iter,
                     const DataType& norm,
                     const DataType& norm_old) const
    {
      const DataType rate = norm/norm_old;
      s << std::setw(iterationSpacing)  << iter << " ";
      s << std::setw(normSpacing) << norm << " ";
      s << std::setw(normSpacing) << rate << std::endl;
    }
 
    template <class DataType>
    void printOutput(std::ostream& s,
                     const double iter,
                     const DataType& norm) const
    {
      s << std::setw(iterationSpacing)  << iter << " ";
      s << std::setw(normSpacing) << norm << std::endl;
    }
  };
 
}
 
#endif