Package: pari
Verstion: latest from git

The following works in 2.7.3 (producing [x^6 - 3*x^5 + 37*x^4 - 10*x^3 - 312*x^2 + 1737*x - 12101] instead of 0), but fails in the current version in git.

John


          GP/PARI CALCULATOR Version 2.8.0 (development 18144-cd853e9)

          amd64 running linux (x86-64/GMP-5.1.3 kernel) 64-bit version

     compiled: Oct 25 2015, gcc version 4.8.4 (Ubuntu 4.8.4-2ubuntu1~14.04)

                            threading engine: single

                 (readline v6.3 enabled, extended help enabled)


                     Copyright (C) 2000-2015 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 ?15 for how to get moral (and possibly technical) support.


parisize = 128000000, primelimit = 500000

? default(realprecision, 1000)

? modulus=[[385, 0, 305; 0, 385, 265; 0, 0, 5], [0]]

%2 = [[385, 0, 305; 0, 385, 265; 0, 0, 5], [0]]

? pol=x^3 - x^2 + 1

%3 = x^3 - x^2 + 1

? bnf=subst(bnfinit(pol),x,y);

? bnfcertify(bnf)

%5 = 1

? bnr=bnrinit(bnf,modulus,1);

? sublis=subgrouplist(bnr,[2])

%7 = [[1, 0, 0; 0, 2, 1; 0, 0, 1]]

? rnfkummer(bnr,%[1])

%8 = 0