Possible bug in znprimroot

[Kevin Acres <research@research-systems.com>]

Hi,

Running the following code:

myroots(p) ={
    local(g);
    g=znprimroot(p^2)^p;
    /*vecsort*/(vector(p-1,i,lift(g^i)));
}
{

    v=myroots(1729);
    print([#v,v]);
}

Gives the following message:


         GP/PARI CALCULATOR Version 2.4.2 (development CHANGES-1.1971)
           i686 running cygwin (ix86/GMP-4.2.1 kernel) 32-bit version
 compiled: Dec 23 2007, gcc-3.4.4 (cygming special, gdc 0.12, using dmd 0.125)
                 (readline v5.2 enabled, extended help enabled)

                     Copyright (C) 2000-2006 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.

parisize = 4000000, primelimit = 2000000
   realprecision = 19 significant digits (10 digits displayed)
  *** znprimroot: bug in PARI/GP (Floating Point Exception), please report
Goodbye!


As requested by the message I am reporting it.

Best Regards,

Kevin.