#ifndef ___Sparse2323___
#define ___Sparse2323___
#include <assert.h>
#include <iostream>
#include <stdlib.h>
using namespace std;

class SparseMatrix{
private:
  
  friend void solve(SparseMatrix& M, const Vector& v, Vector& soln);
  friend void mult(SparseMatrix& M, const Vector& v, Vector& prod);
  Matrix A;
public:
  //d = dimension; max = max num of nonzeros
  SparseMatrix(size_t d, size_t max) : A(d,d){
    cout<<"sparsematrix: d = "<<d<<endl;
    for(int i=0; i < d; i++)
      for(int j=0; j < d; j++)
	A(i,j) = 0.0;
  }

  //insert value x at (i,j)
  void insert(size_t i, size_t j, double x){
    A(i,j) = x;
  }

  //insert matrix M at (i,j)
  void insert(size_t i, size_t j, const Matrix& M){
    size_t m = M.getm();
    size_t n = M.getn();
    assert(i < dimension);
    assert(j < dimension);
    assert(i+m-1 < dimension);
    assert(j+n-1 < dimension);
    for(size_t ii=0; ii < m; ii++)
      for(size_t jj=0; jj < n; jj++)
	insert(i+ii, j+jj, M(ii,jj));
  }
  
  void output(const char *s) {A.output(s);}
  void formMatrix(int orderType=0){};
  ~SparseMatrix(){};
};

#endif