Re: bnrinit() bug

[Phil Carmody <thefatphil@yahoo.co.uk>]

--- Igor Schein <igor@txc.com> wrote:
> Package: pari
> Version: 2.2.8 (development CHANGES-1.944)
> Severity: normal
> 
> \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\
>    realprecision = 28 significant digits
>    seriesprecision = 16 significant terms
>    format = g0.28
> 
> parisize = 8000000, primelimit = 500000
> ?
>
bnrinit(bnfinit(y^16-232*y^14+17564*y^12-592696*y^10+10090294*y^8-87937112*y^6+356253116*y^4-482477960*y^2+143400625),60,1);
>   *** bnrinit: impossible inverse modulo: Mod(5, 300).
> \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

I assume that's 32-bit x86? 2 other arch's - one good (straight 2.2.8), one bad
(2.2.8 patched to 1.944).

Athlon64 with GMP:
<<<
                 GP/PARI CALCULATOR Version 2.2.8 (development)
          x86_64 running linux (portable C/GMP kernel) 64-bit version
compiled: Jan  2 2004, gcc-3.3.2 20031218 (Gentoo Linux 3.3.2-r5,
propolice-3.3-7)
                (readline v4.3 enabled, extended help available)

                       Copyright (C) 2003 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and 
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 38 significant digits
   seriesprecision = 16 significant terms
   format = g0.38

parisize = 8000000, primelimit = 500000

[...]
?
bnrinit(bnfinit(y^16-232*y^14+17564*y^12-592696*y^10+10090294*y^8-87937112*y^6+356253116*y^4-482477960*y^2+143400625),60,1);
  ***   the PARI stack overflows !
  current stack size: 64000000 (61.035 Mbytes)
  [hint] you can increase GP stack with allocatemem()

? allocatemem()                                                                
  ***   Warning: doubling stack size; new stack = 128000000 (122.070 Mbytes).
? allocatemem()
  ***   Warning: doubling stack size; new stack = 256000000 (244.141 Mbytes).
?
bnrinit(bnfinit(y^16-232*y^14+17564*y^12-592696*y^10+10090294*y^8-87937112*y^6+356253116*y^4-482477960*y^2+143400625),60,1);
>>>

Alpha (21164)
<<<
phil@megaspaz:Oalpha-gmp-gcc$ ./gp-sta 
Reading GPRC: /home/phil/.gprc ...Done.

          GP/PARI CALCULATOR Version 2.2.8 (development CHANGES-1.944)
          Alpha running linux (Alpha/GMP-4.1.2 kernel) 64-bit version
        compiled: May 31 2004, gcc-3.3.3 20040125 (prerelease) (Debian)
              (readline v4.3 enabled, extended help not available)

                       Copyright (C) 2003 The PARI Group

PARI/GP is free software, covered by the GNU General Public License, and 
comes WITHOUT ANY WARRANTY WHATSOEVER.

Type ? for help, \q to quit.
Type ?12 for how to get moral (and possibly technical) support.

   realprecision = 38 significant digits
   seriesprecision = 16 significant terms
   format = g0.38
[...]
(09:38) gp > allocatemem()
  *** allocatemem: Warning: doubling stack size; new stack = 256000000 (244.141
Mbytes).
(09:38) gp >
bnrinit(bnfinit(y^16-232*y^14+17564*y^12-592696*y^10+10090294*y^8-87937112*y^6+356253116*y^4-482477960*y^2+143400625),60,1);
  *** bnrinit: impossible inverse modulo: Mod(5, 300).
>>>



=====
When inserting a CD, hold down shift to stop the AutoRun feature
In the Device Manager, disable the SbcpHid device.
http://www.cs.princeton.edu/~jhalderm/cd3/


	
		
__________________________________
Do you Yahoo!?
Friends.  Fun.  Try the all-new Yahoo! Messenger.
http://messenger.yahoo.com/