CPPLapack Code

//=============================================================================
//! Samll Real Double-precision General Dence Matrix Class
template<long m, long n> class dgematrix_small
{
public:
  //////// data ////////
  double array[m*n];
  
  //////// constructor ////////
  dgematrix_small(){;}
  dgematrix_small(const double& x){
    for(long k=0; k<m*n; k++){ array[k]=x; }
  }
  ~dgematrix_small(){;}
  
  //////// function ////////
  double& operator()(const long& i, const long& j);
  double  operator()(const long& i, const long& j)const;
  //long M(){ return m; }
  //long N(){ return n; }
  dgematrix_small<m,n>& zero(){
    for(long k=0; k<m*n; k++){ array[k]=0.; }
    return *this;
  }
  dgematrix_small<m,n>& identity(){
    zero();
    for(long k=0; k<std::min(m,n); k++){ (*this)(k,k)=1.; }
    return *this;
  }
  dgematrix_small<m,n>& set(const long& i, const long& j, const double& v){
    (*this)(i,j)=v;
    return *this;
  }  
  
  dcovector_small<m> col(const long& j) const{
    dcovector_small<m> vec;
    for(long i=0; i<m; i++){ vec(i)=(*this)(i,j); }
    return vec;
  }
  
  void write(const char* filename) const;
  void read(const char* filename);
  
  _dgematrix to_dgematrix() const{
    dgematrix mat(m,n);
    for(long i=0; i<m; i++){
      for(long j=0; j<n; j++){
        mat(i,j)=(*this)(i,j);
      }
    }
    return _(mat);
  }
};