How do you get an iterator to a row of a boost::numeric::ublas::matrix?

Question

I'm working on matrix multiplication, and I would like an iterator over a single row of a boost matrix? Can this be done?

Currently, I have to get an iterator and advance it. It seems like too much CPU work / non-optimized...

boost::numeric::ublas::matrix<T> aMatrix(2048, 4096);
typename boost::numeric::ublas::unbounded_array<T>::iterator it;
it = aMatrix.data().begin();
offset = row * aMatrix.size2();
advance(it, offset);

Forgo iterators and use the array like operations? – Mikhail Mar 26 '12 at 01:41 — Mikhail, Mar 26 '12 at 01:41

Zak · Accepted Answer · 2012-03-25T16:59:16.630

4

Eureka! Matrix proxies...

boost::numeric::ublas::matrix_row<boost::numeric::ublas::matrix<T> > aRow(aMatrix, row);

edited Mar 25 '12 at 16:59

answered Mar 25 '12 at 16:41

Zak

12,213
21
59
105

How do you get an iterator to a row of a boost::numeric::ublas::matrix?

1 Answers1