lanczos.hxx

Go to the documentation of this file.

#ifndef CH_MATRIX_CLASSES__LANCZOS_HXX
#define CH_MATRIX_CLASSES__LANCZOS_HXX

#ifndef CH_MATRIX_CLASSES__MATRIX_HXX
#include "matrix.hxx"
#endif

namespace CH_Matrix_Classes {


  class Lanczosmatrix
  {
  public:
    virtual ~Lanczosmatrix(){}

    virtual Integer lanczosdim() const=0;
  
    virtual Integer lanczosflops() const=0;

    virtual int lanczosmult(const Matrix& A,Matrix& B) const=0;
  };


  class Lanczos
  {
  public:
    virtual ~Lanczos(){}
    
  
    virtual void set_mineig(Real ie)=0;
    virtual void set_maxmult(Integer mop)=0;
    virtual void set_maxiter(Integer mi)=0;
    virtual void set_relprec(Real relprec)=0; 
    virtual void set_nblockmult(Integer nb)=0;
    virtual void set_nchebit(Integer nc)=0;
    virtual void enable_stop_above(Real ub)=0;
    virtual void disable_stop_above()=0;

    virtual void set_retlanvecs(Integer nl)=0; 

    virtual int get_lanczosvecs(Matrix& val,Matrix& vecs) const=0;
    
    virtual Real get_relprec(void)=0;

    virtual int get_err() const=0;
    
    virtual Integer get_iter() const=0;
    
    virtual Integer get_nmult() const=0;

    
    virtual int compute(const Lanczosmatrix* bigmat, 
                        Matrix& eigval,        
                        Matrix& eigvec,        
                        Integer nreig,         
                        Integer in_blocksz=0,  
                        Integer maxcol=0       
                        )=0;

  
    virtual void set_out(std::ostream* out=0,int print_level=1)=0;

    virtual std::ostream& save(std::ostream& out) const =0;

    virtual std::istream& restore(std::istream& in) =0;

  };

  
}

#endif