| Re: [eigen] problem computing kernel in eigen |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: Re: [eigen] problem computing kernel in eigen
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Wed, 27 Feb 2013 21:01:57 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=x-received:mime-version:in-reply-to:references:from:date:message-id :subject:to:content-type:content-transfer-encoding; bh=vfp8EEMgrnW8vDkcMqws0lF1LhlUwhbOr/JfX9HJLf0=; b=mcOBpOzcIndDmTPlSkAE1aQwjDfcVvP15ZsVc4RtFQSjv+sPdO0HTw+g1+aT/lNetf JsRNGhUlggBPZlq+1J1bQx2mSCBoUhl0b4nesHZHXNJCZJJxM8xKuWzACieUsrx+rwR8 aq6mveJTZAzYz7MFoQ5g4wy+AzKWz967761ntjYS9s3kHAoticbzvHxlSscteTjCxTl9 nGC18eLyT1Q4yqHYc1wvwAD5AeASIoNCGLcHPM6wOBUCTAChXt6iDkHZLyXlkzMAcuLf 3TykXw4xdSaeAQIKbUClcbJ0fSCIHiC8v8AreBfg/a0rpu48sZjiXcvsFtXF/8K2J7Kp e/cQ==
Good news, thank you for the report.
gael
On Wed, Feb 27, 2013 at 12:45 PM, Vissarion Fisikopoulos
<vfisikop@xxxxxxxxx> wrote:
> Hi,
>
> I tried the latest version and it works fine.
>
> Thanks a lot,
>
> vissarion.
>
>
> 2013/2/26 Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
>>
>> Hi,
>>
>> it seems you are using a quite old version of Eigen. In the 3.1
>> branch, adding internal:: sounds reasonable. I'll apply such a change.
>>
>> In the devel branch, this should work just fine. Since we are going to
>> release soon 3.2-beta1, it would be very helpful for us if you could
>> give it a try. This is also a good way to make sure the next versions
>> will work fine for you! You just have to download the latest version
>> there: http://bitbucket.org/eigen/eigen/get/default.tar.bz2.
>>
>> thanks,
>>
>> gael
>>
>> On Tue, Feb 26, 2013 at 4:09 PM, Vissarion Fisikopoulos
>> <vfisikop@xxxxxxxxx> wrote:
>> > Hello,
>> >
>> > I have a problem with a simple call to the kernel (function call
>> > fullPivLu().kernel() ) of an
>> > Eigen::Matrix<NT,Eigen::Dynamic,Eigen::Dynamic>
>> > . This works for many different instantiations of NT, like mpq_class.
>> > The
>> > problem is when I try to use a rational number type from the CGAL
>> > library,
>> > CGAL::Gmpq. I get the following error from the compiler:
>> >
>> > /usr/include/eigen3/Eigen/src/LU/FullPivLU.h:595:7: error: call of
>> > overloaded ‘abs(const Scalar&)’ is ambiguous
>> >
>> > The error is fixed by replacing
>> >
>> > if(abs(dec().matrixLU().coeff(i,i)) > premultiplied_threshold)
>> >
>> > with
>> >
>> > if(internal::abs(dec().matrixLU().coeff(i,i)) >
>> > premultiplied_threshold)
>> >
>> > in line 595 of /usr/include/eigen3/Eigen/src/LU/FullPivLU.h
>> >
>> > Do you think this is a bug?
>> >
>> > Best,
>> > Vissarion Fisikopoulos.
>> >
>> >
>>
>>
>