tsparsematrix.hxx

Go to the documentation of this file.

#ifndef NUMCXX_TSPARSEMATRIX_H
#define NUMCXX_TSPARSEMATRIX_H

#include <vector>
#include "tarray1.hxx"
#include "tmatrix.hxx"

namespace numcxx
{
    
  template<typename T>   class TSolverUMFPACK;
  template<typename T>   class TPreconJacobi;
  template<typename T>   class TPreconILU;

  template<typename T> 
  class TSparseMatrix: public TLinOperator<T>, public SparseMatrixExpressionBase
  {
    friend class TSolverUMFPACK<T>;
    friend class TPreconJacobi<T>;
    friend class TPreconILU<T>;
  public:
    typedef T value_type;

    TSparseMatrix():n(0) {};

    TSparseMatrix(index n1, index n2);


    static std::shared_ptr<TSparseMatrix <T> > create(index n1, index n2) {return std::make_shared<TSparseMatrix<T>>(n1,n2);}

    TSparseMatrix(const  std::initializer_list<std::initializer_list<T>> &il);

    static std::shared_ptr<TSparseMatrix <T> > create(const  std::initializer_list<std::initializer_list<T>> &il);


    T& operator()(int i, int j);
        
    void  flush();


    void  clear() {  flush(); (*pIA)=0.0;}

    void apply(const TArray<T> &U,TArray<T> &V ) const;

    std::shared_ptr<TMatrix<T>> copy_as_dense();

    index shape(int idim)  const  {return n;}
        
    TSparseMatrix(const TSparseMatrix<T>& A)=delete;


    // std::shared_ptr<TSparseMatrix <T> > copy() const;
    // std::shared_ptr<TSparseMatrix <T> > clone() const;
    // T xentry(const index i, const index j) const {return _data[_idx(i,j)];}
    // template <typename VAL>
    // TSparseMatrix<T>&  operator=(const VAL  &expr)  { assign(*this,expr); return *this;}
    // TSparseMatrix<T>&  operator=(const TSparseMatrix<T> &expr) { assign(*this,expr); return *this;}
    // void apply(const TArray1<T> &u, TArray1<T> &v) const;

    std::shared_ptr<TArray1<int>> pIA;
        
    std::shared_ptr<TArray1<int>> pJA;  
        
    std::shared_ptr<TArray1 <T> > pA; 

    std::shared_ptr<TMatrix<T>> calculate_inverse();

    bool pattern_changed()  const  {return _pattern_changed;}

  private:




    void pattern_changed(bool chg) {_pattern_changed=chg;};

        
    bool _pattern_changed=false;
    const index n;
    int maxrow=0;

    bool _first_flush_done=false;
    bool empty()  const { return !_first_flush_done;};
    class RowEntry;

    class Extension;

    std::shared_ptr<Extension>  pExt=nullptr;
        
  };
    
}

#include "tsparsematrix.ixx"
#endif