Re: 64bit regression

[Igor Schein <igor@txc.com>]

On Tue, Sep 07, 2004 at 06:03:50PM +0200, Bill Allombert wrote:
> 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;
> }

beta=4.29497e+09
p=1.84467e+19
u=9223372036854775807
l=9223372036854775807

Igor