Bigmatrix.hxx

Go to the documentation of this file.

#ifndef CONICBUNDLE_BIGMATRIX_HXX
#define CONICBUNDLE_BIGMATRIX_HXX

#include "lanczos.hxx"
#include "SparseCoeffmatMatrix.hxx"

namespace ConicBundle {

  class Bigmatrix:
    public CH_Matrix_Classes::Lanczosmatrix,
    protected CH_Matrix_Classes::Memarrayuser
  {
  private:
    typedef std::map<CH_Matrix_Classes::Integer,CH_Matrix_Classes::Real> Indexhash;

    CH_Matrix_Classes::Real    tol; 

    CH_Matrix_Classes::Integer dim; 

    CH_Matrix_Classes::Integer nz;  
    CH_Matrix_Classes::Indexmatrix colnz; 
    std::vector<CH_Matrix_Classes::Indexmatrix> rowind; 
    std::vector<CH_Matrix_Classes::Matrix> rowval; 

    CH_Matrix_Classes::Matrix di;   
    std::vector<Indexhash> rowhash;  
    
    std::vector<CoeffmatPointer> mcp; 

    CH_Matrix_Classes::Matrix mcv; 
    
    mutable CH_Matrix_Classes::Integer nmult; 
    
    //---
    bool use_dense; 
    mutable CH_Matrix_Classes::Symmatrix symrep;  
    bool symrep_init; 
    
    //--- temporary variables for lanczosmult
    mutable CH_Matrix_Classes::Matrix At; 
    mutable CH_Matrix_Classes::Matrix Bt; 
    CH_Matrix_Classes::Indexmatrix collecti; 
    CH_Matrix_Classes::Indexmatrix collectj; 
    CH_Matrix_Classes::Matrix collectval; 
    
 public:
    Bigmatrix();
    ~Bigmatrix();
    void clear();
    CH_Matrix_Classes::Integer get_nmult() const {return nmult;}
    void reset_nmult() {nmult=0;}
    int init(const CH_Matrix_Classes::Matrix& yi,CH_Matrix_Classes::Integer indim,const CoeffmatPointer C,
             const SparseCoeffmatVector* Ai,const bool dense=false);

    void set_tol(CH_Matrix_Classes::Real t){tol=t;}

    int get_dense() const {return use_dense;}

    virtual CH_Matrix_Classes::Integer lanczosdim() const;

    virtual CH_Matrix_Classes::Integer lanczosflops() const;

    virtual int lanczosmult(const CH_Matrix_Classes::Matrix& A,CH_Matrix_Classes::Matrix& B) const;

    int make_symmatrix(CH_Matrix_Classes::Symmatrix& S) const;

    const CH_Matrix_Classes::Symmatrix& get_symrep() const
    { make_symmatrix(symrep); return symrep; }
    
    friend std::ostream& operator<<(std::ostream& out,const Bigmatrix&);
};


} //end namespace ConicBundle

#endif