1 | #include <numeric>
|
---|
2 | #include <boost/numeric/ublas/symmetric.hpp>
|
---|
3 | #include <boost/numeric/ublas/matrix_proxy.hpp>
|
---|
4 | #include <boost/numeric/ublas/io.hpp>
|
---|
5 | using namespace std;
|
---|
6 | using namespace boost::numeric::ublas;
|
---|
7 |
|
---|
8 | int main()
|
---|
9 | {
|
---|
10 | symmetric_matrix<double> m(3, 3);
|
---|
11 | for(unsigned int i = 0; i < 3; i++)
|
---|
12 | for(unsigned int j = 0; j < 3; j++)
|
---|
13 | m(i, j) = 1;
|
---|
14 | for(unsigned int i = 0; i < 3; i++)
|
---|
15 | {
|
---|
16 | matrix_row<symmetric_matrix<double> > row(m, i);
|
---|
17 | cout << accumulate(row.begin(), row.end(), 0.0) << endl; //should be 3 for all rows
|
---|
18 | }
|
---|
19 | return 0;
|
---|
20 | }
|
---|