[ublas] Splitting matrix into diagonal, lower, upper triangles?

From: Giuseppe Bilotta (giuseppe.bilotta_at_[hidden])
Date: 2008-02-23 17:41:53

Hello all,

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

M.diagonal()

to return a diagonal_adaptor with the diagonal of M,

M.lower_triangle()
M.upper_triangle()

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