#ifndef _densematrix_h_
#define _densematrix_h_

#include <vector>
using namespace std;

#include "matrix.h"

template <typename ValType>
class  DenseMatrix: public Matrix<ValType>
{
public:
  DenseMatrix(int nrows, int ncols);  
   
  ValType & operator()(int i, int j);
 
private:  
  vector<ValType> _values;
};

template <typename ValType>
DenseMatrix<ValType>::DenseMatrix(int nrows, int ncols) : Matrix<ValType>(nrows, ncols)
{
  _values.resize(nrows * ncols);
}

template <typename ValType>
ValType & DenseMatrix<ValType>::operator()(int i, int j)
{
  return _values[j + i * this->_ncols];
}


#endif