
Ublas : 
Subject: Re: [ublas] symmetric matrices, row and columnmajor
From: Gunter Winkler (guwi17_at_[hidden])
Date: 20090225 14:18:24
Tiago Requeijo schrieb:
> Looking at functional.hpp, I think the problem is indeed the
> transposed_structure template. As far as I understand, the template
> simply takes an upper matrix and regards it as the transpose of a
> lower matrix. While this is (of course) correct, it doesn't take into
> account the row/column major storage:
>
> If we take the upper matrix U
>
> 1 2 3
>  4 5
>   6
>
> and regard it as the transpose of a lower matrix L
>
> 1
> 2 4
> 3 5 6
>
> without a row/column major swap, then we're storing the row major
> version of the lower matrix or, in other words, the column major
> version of the upper matrix:
>
> a) row_major for U: 1 2 3 4 5 6
> b) row_major for L: 1 2 4 3 5 6
>
> Since the transposed_structure template simply maps the elements to
> the transpose and doesn't change the storage, we are in fact storing
> b) instead of a)
Ok. Thus it is a bug. Could you please open a ticket at
https://svn.boost.org/trac/boost/newticket and attach you examples. Then
we can easily keep track of the problem an I'll spent some time to fix it.
mfg
Gunter