minres.hxx

Go to the documentation of this file.

#ifndef CH_MATRIX_CLASSES__MINRES_HXX
#define CH_MATRIX_CLASSES__MINRES_HXX

#include "IterativeSystemObject.hxx"

namespace CH_Matrix_Classes {


  class MinRes: public IterativeSolverObject
  {
  private: 
    Integer maxit; 
    Real resnorm;   
    Real old_resnorm; 
    Real avg_reduction; 
    Real termprec;     
    Integer nmult;  
    Integer err;    

    
    Matrix vecvcurr;    
    Matrix vecvprev;    
    Matrix vecwcurr;    
    Matrix vecwprev;    
    Matrix veczcurr;    
    Matrix veczprev;    
    Matrix vecAv;       
    
    Integer stall;     

    std::ostream* myout;   
    int print_level;   

  public:
    MinRes(std::ostream* out=0,int pril=-1); 
    virtual ~MinRes(){}
    
  
    void set_maxit(Integer in_maxit) {maxit=in_maxit;}
    Integer get_maxit() const {return maxit;}

    int get_err() const {return err;}
    Integer get_nmult() const {return nmult;}
    Real get_residual_norm() const {return resnorm;}
    virtual Real get_avg_reduction() const {return avg_reduction;}
    virtual Real get_termprec() const {return termprec;}

    
    int compute(IterativeSystemObject& system, 
                Matrix& x,        
                Real termprec,  
                Matrix* storex=0, 
                Integer storestep=0 
                );

  
    virtual void set_out(std::ostream* in_out=0,int in_print_level=1)
    {myout=in_out;print_level=in_print_level;}

  };

  
}

#endif