9 Mar
2011
9 Mar
'11
9:09 a.m.
On 09/03/2011 09:59, sguazt wrote:
namespace ublas = ::boost::numeric::ublas;
template<typename MatrixT> inline ublas::matrix_range<MatrixT> topleft(MatrixT& M) { return ublas::subrange(M, 0, 3, 0, 3); }
and use this function just like the case here below:
matrix<double> A(9,9); topleft(A) = zero_matrix<double>(3,3);
You might want to provide a const overload as well.