
Ublas : 
From: Giuseppe Bilotta (giuseppe.bilotta_at_[hidden])
Date: 20080223 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