Wie erhalten Sie einen Iterator in eine Zeile einer boost :: numeric :: ublas :: matrix <T>?

Question

Ich arbeite an Matrix-Multiplikation, und ich möchte einen Iterator über eine einzelne Zeile einer Boost-Matrix? Kann das gemacht werden?

Derzeit muss ich einen Iterator bekommen und es voranbringen. Es scheint wie zu viel CPU-Arbeit/nicht-optimierten ...

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); 

Answer 1

Eureka! Matrix proxies ...

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