/* posted from kyriales at ublas-dev@yahoogroups.com */
#include <boost/numeric/ublas/matrix_sparse.hpp>
#include <boost/numeric/ublas/io.hpp>

using std::cout;
using std::endl;

typedef boost::numeric::ublas::compressed_matrix<double> spm_t;

int main() {
  {
    boost::numeric::ublas::compressed_matrix<double> mat(2,1);
    mat(1,0) = 1.0;
    cout<<mat<<endl;
  }
  {
    boost::numeric::ublas::compressed_matrix<double> mat(10,3);
    mat(6,0) = 1.;
    mat(7,1) = 1.;
    mat(8,2) = 1.;
    cout<<mat<<endl;
    mat(9,0) += 42;
    cout<<mat<<endl;
  }
  return EXIT_SUCCESS;
}

/*
Output is:
[2,1]((0),(1))

[10,3]((0,0,0),(0,0,0),(0,0,0),(0,0,0),(0,0,0),(0,0,0),(1,0,0),(0,1,0),(0,0,1),(0,0,0))
[10,3]((0,0,0),(0,0,0),(0,0,0),(0,0,0),(0,0,0),(0,0,0),(1,0,0),(0,1,0),(0,0,1),(42,0,0))
*/