Boost logo

Ublas :

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