#include <numeric>
#include <boost/numeric/ublas/symmetric.hpp>
#include <boost/numeric/ublas/matrix_proxy.hpp>
#include <boost/numeric/ublas/io.hpp>
using namespace std;
using namespace boost::numeric::ublas;

int main()
{
	symmetric_matrix<double> m(3, 3);
	for(unsigned int i = 0; i < 3; i++)
		for(unsigned int j = 0; j < 3; j++)
			m(i, j) = 1;
	for(unsigned int i = 0; i < 3; i++)
	{
		matrix_row<symmetric_matrix<double> > row(m, i);
		cout << accumulate(row.begin(), row.end(), 0.0) << endl; //should be 3 for all rows
	}
	return 0;
}