From: Giuseppe Bilotta (giuseppe.bilotta_at_[hidden])
Date: 2008-02-23 17:41:53
I was wondering if it was possible to 'decompose' a matrix into its diagonal
and lower and upper triangles, in uBLAS.
The idea would be, given some matrix M, to have some methods like
to return a diagonal_adaptor with the diagonal of M,
to return triangular matrices of the appropriate type with the lower and upper
triangles of the matrix M (thus excluding the diagonal of the matrix), in such
a way that
M = M.diagonal() + M.lower_triangle() + M.upper_triangle()
Is there anything like this in uBLAS?
-- Giuseppe "Oblomov" Bilotta