pcg.hxx

Go to the documentation of this file.

#ifndef CH_MATRIX_CLASSES__PCG_HXX
#define CH_MATRIX_CLASSES__PCG_HXX

#include "IterativeSystemObject.hxx"

namespace CH_Matrix_Classes {


  class PCG: public IterativeSolverObject
  {
  private: 
    Integer maxit; 
    Real resnorm;   
    Real old_resnorm; 
    Real avg_reduction; 
    Real termprec;     
    Integer nmult;  
    Integer err;    
    
    Matrix Ap;   
    Matrix r;    
    Matrix z;    
    Matrix p;    

    Integer stall; 


    std::ostream* myout;   
    int print_level;   

  public:
    PCG(std::ostream* out=0,int pril=-1); 
    ~PCG(){}
    
  
    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