Re: [Boost-Users] Re: uBlas - diagonal_matrix

Hi Patrick, you wrote:

Answering myself,

That happens sometimes ;-)

without a symmetric type I can just check for a diagonal matrix inside the pinv function and branch, see code below.

As far as I've understood you're looking for a type safe diagonal variant of banded_matrix. Would it be helpful, if we'd go to add something like // Diagonal matrix class template class diagonal_matrix: public banded_matrix { public: BOOST_UBLAS_USING banded_matrix::operator =; typedef banded_matrix matrix_type; // Construction and destruction BOOST_UBLAS_INLINE diagonal_matrix (): matrix_type () {} BOOST_UBLAS_INLINE diagonal_matrix (std::size_t size1, std::size_t size2): matrix_type (size1, size2) {} BOOST_UBLAS_INLINE ~diagonal_matrix () {} // Assignment BOOST_UBLAS_INLINE diagonal_matrix &operator = (const diagonal_matrix &m) { matrix_type::operator = (m); return *this; } }; and the corresponding diagonal_adaptor to the code base?

Are there some recommendations for namespaces for functions like "pinv" ? Inside the loop, shall I use size_type?

I currently haven't enough time to look into a SVD implementation. So I'll cc this reply to http://groups.yahoo.com/group/ublas-dev. HTH, Joerg [snip]