CPPLapack Code

//=============================================================================
//! Base class for Small Real Double-precision General Dence Matrix Class
class dgematrix_small_base
{
public:
  //////// constructor ////////
  dgematrix_small_base(){;}
  virtual ~dgematrix_small_base(){;}
};

//=============================================================================
//! Samll Real Double-precision General Dence Matrix Class
template<long m, long n> class dgematrix_small : public dgematrix_small_base
{
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){ return array[n*i+j]; }
  double  operator()(const long& i, const long& j)const{ return array[n*i+j]; }
  //long get_m(){ return m; }
  //long get_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<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;
  }  
};

///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////

//=============================================================================
template<long m, long n>
inline std::ostream& operator<<(std::ostream& s, const dgematrix_small<m,n>& A)
{
  s << std::setiosflags(std::ios::showpos);
  for(long i=0; i<m; i++){
    for(long j=0; j<n; j++){
      s << " " << A(i,j);
    }
    s << std::endl;
  }
  return s;
}