| 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:52:56 +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=aUZDFqWjcvzlzZxZccj6apo/VDwA9WBAdcChBO166fY=; b=KUYF0PXHfwWYXfMdV2RkkdmbZAL7Mr5YhRXD7kQGi04XMqhb5YGLDSjmH7ku6sHWa6 TXc0rJUWgWibaWW5b1S2YYyh4cTiNAHeLb2VOVgkBsoInJBuvdmvTm+o5hkluI76NW61 X2/jJJvCDHbGvPppYk6JykvIRGPho6MLgAvHQ=
- 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=nd+9F0E5OZZNr8Ts5WoFTWrr7PoNrC3EREplyKSY5t08R7kfO6I9PFLxYzppzqUfBq v8cAJMe3MtYrWiVspZzfplznd3mhDIXNNfy8PRkH/sP9W5BwHa/qBFGXOohrLCThCkIJ t9/GZx6kZney7zvObyKP8B9FmV+us3Cg1rebk=
ok, that was a renaming issue "only". Problem fixed.
gael
On Tue, Nov 16, 2010 at 3:47 PM, Gael Guennebaud
<gael.guennebaud@xxxxxxxxx> wrote:
> 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
>>
>>
>