Re: [eigen] How to obtain the matrices Q and R of a QR decomposition

[ Thread Index | Date Index | More Archives ]

this time I confirm the issue. In the meantime use -msse2 only.


On Tue, Nov 16, 2010 at 3:40 PM, Sebastian Birk
<birk@xxxxxxxxxxxxxxxxxxxxx> wrote:
> Hi,
> when I compile that QR decomposition code using SSE extensions like with
> this optimization parameters
>  -O3 -mtune=core2 -mfpmath=sse -msse4.1
> the computation of the QR decomposition causes an segmentation fault,
> which I tracked back to the source line
>  return pmul(a, pconj(b));
> in the function
>  template<> struct conj_helper<Packet1cd, Packet1cd, false,true>
> in SSE/Complex.h.
> Is this an error in Eigen or is it me again doing weird stuff? ;-)
> Sebastian

Mail converted by MHonArc 2.6.19+