Karim BELABAS on Wed, 27 Nov 2002 15:56:15 +0100 (MET)


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: bug on Itanium


On Fri, 22 Nov 2002, Igor Schein wrote:
> On Thu, Nov 21, 2002 at 07:46:23PM +0100, Karim BELABAS wrote:
>> On Thu, 21 Nov 2002, Bill Allombert wrote:
>>
>>> On Thu, Nov 21, 2002 at 12:32:43PM +0100, Bill Allombert wrote:
>>>> On the other hand, on ia64 with gcc 3.2 and 2.96, I get
>>>> qflll(1/mathilbert(7)) --> oo loop
>>>
>>> It also happen on alpha, so I suppose it is 64 bits specific.
>>
>> I very recently modified the floating point LLL code for exact inputs, and
>> did not pay attention to [ pre-existing ] dependance on sizeof(long) [ my
>> main machines all have 32-bit architectures ].
>>
>> I have fixed the problem.
>
> There's still cases resulting in infinite loop on Itanium,
> here's the smallest one I was able to produce:
>
> nffactor(nfinit(polsubcyclo(19,6,y)),x^64-x^52-5712805568332331698655486483884351420000000*x^20-1)

I have slightly changed the LLL strategy, fixing the above.

Anything more ?

    Karim.
-- 
Karim Belabas                    Tel: (+33) (0)1 69 15 57 48
Dép. de Mathématiques, Bât. 425  Fax: (+33) (0)1 69 15 60 19
Université Paris-Sud             Email: Karim.Belabas@math.u-psud.fr
F-91405 Orsay (France)           http://www.math.u-psud.fr/~belabas/
--
PARI/GP Home Page: http://www.parigp-home.de/