#ifndef _sparsematrix_h_
#define _sparsematrix_h_

#include "matrix.h"
#include <map>
using namespace std;

template <typename ValType>
class  SparseMatrix: public Matrix<ValType>
{
public:
  SparseMatrix(int nrows, int ncols);
  
  ValType & operator()(int i, int j); 
  

private:  
  map<int, ValType> _values; 
};

template <typename ValType>
SparseMatrix<ValType>::SparseMatrix(int nrows, int ncols) : Matrix<ValType>(nrows, ncols)
{     
}

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

#endif