I have 2 matrices Qex,R of type float which are returned from a
function (qr decomposition). I then do


Float_Matrix M=Qex*R; //wrong result

Float_Matrix Qex_copy(Qex);//copy constructor
Float_Matrix R_copy(R);//copy constructor

Float_Matrix Mcopy=Qex_copy*R_copy;//correct result

For some reason M is incorrect!

I do a cout<<Qex<<R; and cout<<Qex_copy<<R_copy;
to visually compare Qex with Qex_copy and R with R_copy and they are
equal!!!
However I think they are not really equal but cant figure out where.

I tried running the same program with type complex (Complex_Matrix) and
everything works fine.

Any idea? Can it be a misalignment issue???

bluekite2000@gmail.com wrote:[color=blue]
> I have 2 matrices Qex,R of type float which are returned from a
> function (qr decomposition). I then do
>
>
> Float_Matrix M=Qex*R; //wrong result
>
> Float_Matrix Qex_copy(Qex);//copy constructor
> Float_Matrix R_copy(R);//copy constructor
>
> Float_Matrix Mcopy=Qex_copy*R_copy;//correct result
>
> For some reason M is incorrect!
>
> I do a cout<<Qex<<R; and cout<<Qex_copy<<R_copy;
> to visually compare Qex with Qex_copy and R with R_copy and they are
> equal!!!
> However I think they are not really equal but cant figure out where.
>
> I tried running the same program with type complex (Complex_Matrix)
> and everything works fine.
>
> Any idea? Can it be a misalignment issue???[/color]

No idea. But I don't think it's "a misalignment issue". But do
check out FAQ 5.8.

V