Re: [eigen] How to obtain the matrices Q and R of a QR decomposition |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] How to obtain the matrices Q and R of a QR decomposition
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Tue, 16 Nov 2010 15:47:28 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:mime-version:received:in-reply-to :references:from:date:message-id:subject:to:content-type :content-transfer-encoding; bh=GfK3TyXZAEDyFFbEoRczPad275kCBJ+k1geidCubTwc=; b=ioiwJBvlU6UA6KbU6t79nh2EafkhGGMWkhoc6FJQZI2iyUcZMomqCrwSp5GfIK/NY7 TUmLB2hFdSgqNawLJsOE6J7cWB5VfP1uI0NcK7QZ1iTypTbeCBwPsnFFc+A66yzY2C1v P9Lxo5HzxAnG2NY+uOuMDioBx2d0k9IdGWYoU=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :content-type:content-transfer-encoding; b=ad62idc0cnZS0/Hkwbcbs1TAnDRsmcpBLjireKbqhAW5/W4vYO5w3Gh21m2SkBanOx HSKoMJDf8DQolsfoejNWWIBHBxkcR9wCbWh/QsJWWxDlJmI4z3e3u5/CbAh9QHbOP5rj bolxKyPeEXOguomqzvVPhf8+rvSO74kWLqRrU=
this time I confirm the issue. In the meantime use -msse2 only.
gael
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
>
>