tmatrix.hxx

Go to the documentation of this file.

#ifndef NUMCXX_TMATRIX_H
#define NUMCXX_TMATRIX_H

#include <cmath>
#include <cstdlib>

#include "tarray.hxx"

namespace numcxx
{

  template<typename T> 
  class TMatrix: public TArray<T>, public TLinOperator<T>, public MatrixExpressionBase
  {
  public:
    using TArray<T>::size;
    using TArray<T>::shape;
    using TArray<T>::operator[];
        
        
    TMatrix(index n1, index n2): TArray<T>(n1,n2){_assert_square();};
        

    TMatrix(index n1, index n2, T*data,std::function<void(T*p)> deleter):  TArray<T>(n1,n2,data,deleter){_assert_square();};
        
    TMatrix(index n1, index n2, T*data, std::shared_ptr<void> datamanager): TArray<T>(n1,n2,data,datamanager){_assert_square();};
        

    TMatrix(const  std::initializer_list<std::initializer_list<T>> &il ): TArray<T>(il){_assert_square();};

    TMatrix(const TMatrix<T>& A):TArray<T>(A.shape(0),A.shape(1)){assign(*this,A);}

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

    static std::shared_ptr<TMatrix <T> > create(const  std::initializer_list<std::initializer_list<T>> &il) { return std::make_shared<TMatrix<T>> (il);  }
        
    std::shared_ptr<TMatrix <T> > copy() const  {  return std::make_shared<TMatrix<T>>(*this);  }
        
    std::shared_ptr<TMatrix <T> > clone() const  {  return create(shape(0),shape(1)); } 
        
    T item(index i0,index i1)  { return _data[_idx(i0,i1)];};


    void itemset(index i0, index i1, T x)  { _data[_idx(i0,i1)]=x;};


    const T& xentry(const index i, const index j) const {return _data[_idx(i,j)];}

    std::shared_ptr<TArray1 <T> > const __getitem__(index i0){ return std::shared_ptr<TArray1<T>>(new TArray1<T>(shape(1), &_data[_idx(i0,0)], [](T*p){;}));}


        
    TMatrix<T>&  operator=(const TMatrix<T> &expr) {return static_cast<TMatrix<T>&>(assign(*this,expr));}
        
    TMatrix<T>&  operator=(const T &expr) {return static_cast<TMatrix<T>&>(assign(*this,expr));}
        
    void apply(const TArray<T> &u, TArray<T> &v) const;
        
    std::shared_ptr<TMatrix<T>> calculate_inverse();

    TMatrix(): TArray<T>(){_assert_square();};


    bool is_matrix(){return true;}
        
  private:
    using TArray<T>::_data;
    using TArray<T>::_idx;
    using TArray<T>::_assert_square;
        
  };
    
}

#include "tmatrix.ixx"
#endif