rnfkummer(bnrinit(bnfinit(y^3+y-1,1),[[429, 234, 78; 0, 39, 0; 0, 0, 39], [0]], 1), [2, 1, 1; 0, 1, 0; 0, 0, 1])
used to work, and currently produces 0.
John
* John Jones [2015-10-25 20:18]:
[...]
> 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.
> ? 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
Now fixed in master.
Thanks for your report !
K.B.
--
Karim Belabas, IMB (UMR 5251) Tel: (+33) (0)5 40 00 26 17
Universite de Bordeaux Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation http://www.math.u-bordeaux.fr/~kbelabas/
F-33405 Talence (France) http://pari.math.u-bordeaux.fr/ [PARI/GP]
`