Bill Allombert on Tue, 07 Sep 2004 18:10:37 +0200


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

Re: 64bit regression


On Tue, Sep 07, 2004 at 11:33:14AM -0400, Igor Schein wrote:
> Hi,
> 
>           GP/PARI CALCULATOR Version 2.2.8 (development CHANGES-1.998)
>    UltraSparc running solaris (portable C kernel) 64-bit version -- debugging
>                         compiled: Sep  7 2004, gcc-3.4.1
>             (readline not compiled in, extended help not available)
>    echo = 1 (on)
> ? intnum(x=0,Pi,sin(x),1)
> 
> Never converges.
> 
> This comes right out of sumiter bench.

Could you compile the program below on that box and send me the output ?

Cheers,
Bill.

#include <math.h>
typedef unsigned long ulong;
main()
{
    ulong a=18446744073709551431UL;
    double beta=sqrt((double)a);
    double p=beta*(1UL << 32);
    ulong u = (ulong) p;
    long l= (long) p;
    printf("beta=%g\np=%g\nu=%lu\nl=%ld\n",beta,p,u,l);
    return 0;
}