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

From: Giuseppe Bilotta (giuseppe.bilotta_at_[hidden])
Date: 2008-03-02 07:30:55

On Saturday 01 March 2008 23:33, Gunter Winkler wrote:

> Am Samstag, 23. Februar 2008 23:41 schrieb Giuseppe Bilotta:
>> M = M.diagonal() + M.lower_triangle() + M.upper_triangle()
>>
>> Is there anything like this in uBLAS?
> 
> there are 
> 
> triangular_adaptor and banded_adaptor
> 
> Thus it should be posible to write
> 
> M = matrix<...>(n,n);
> L = triangular_adaptor<lower>(M);
> D = banded_adaptor(M,0,0);
> U = triangular_adaptor<upper>(M);
> 
> which satifies  M == L + D + U;
> 
> (see triangular.hpp and banded.hpp)

Thanks you very much.

-- 
Giuseppe "Oblomov" Bilotta