GP/PARI CALCULATOR Version 2.8.0 (development 19693-a23853f)

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

    compiled: Oct 17 2016, gcc version 4.8.4 (Ubuntu 4.8.4-2ubuntu1~14.04.3)

                            threading engine: single

                 (readline v6.3 enabled, extended help enabled)

                     Copyright (C) 2000-2016 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 = 8000000, primelimit = 500000

? pol=x^7 + 7*x^3 - 7*x^2 + 7*x + 1;

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

? modu=idealmul(bnf,8,idealprimedec(bnf,7)[1]);

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

? sgl=subgrouplist(bnr,[2],1);

? rnfs = apply(z->rnfkummer(bnr,z), sgl)

  ***   at top-level: rnfs=apply(z->rnfkummer(b

  ***                      ^--------------------

  ***   in function apply: rnfkummer(bnr,z)

  ***                      ^----------------

  *** rnfkummer: bug in PARI/GP (Segmentation Fault), please report.

  ***   Break loop: type 'break' to go back to GP prompt

break>