#include<boost/numeric/ublas/matrix_sparse.hpp>
#include<boost/numeric/ublas/io.hpp>

int main() 
{ 
	boost::numeric::ublas::coordinate_matrix<double> matrix_mask(3, 3, 4);
	matrix_mask(1,1) = 2.1;
	matrix_mask(0,1) = 1.1;
	matrix_mask(2,2) = 4.1;
	matrix_mask(0,0) = 0.1;

	std::cout << matrix_mask << std::endl;
	std::copy( matrix_mask.index1_data().begin(),  matrix_mask.index1_data().end(),    
	   std::ostream_iterator<int>(std::cout, ","));
	std::cout << std::endl;
	std::copy(matrix_mask.index2_data().begin(), matrix_mask.index2_data().end(),    
	   std::ostream_iterator<int>(std::cout, ","));
	std::cout << std::endl;
	std::copy(matrix_mask.value_data().begin(), matrix_mask.value_data().end(),
	   std::ostream_iterator<double>(std::cout, ","));
		std::cout << std::endl;
}