\r packageaurel.gp

\\cyclotomic examples
abbnf = abelianbnfinit(polcyclo(216)); \\degree 72, 5 seconds
getcyc(abbnf)
abbnf = abelianbnfinit(polcyclo(504)); \\degree 144, 15 seconds
getcyc(abbnf)
abbnf = abelianbnfinit(polcyclo(1260)); \\degree 288, 3 minutes
getcyc(abbnf)

pol = multiquad([-1,2,3,5]);
abbnf = abelianbnfinit(pol);
getcyc(abbnf)
bnf = bnfinit(pol);
bnf.cyc

pol = multiquad([-1,2,3,5,7]);
abbnf = abelianbnfinit(pol);
getcyc(abbnf)
bnf = bnfinit(pol);
bnf.cyc

pol = multiquad([-1,2,3,5,7,11]);
abbnf = abelianbnfinit(pol);
getcyc(abbnf)
\\bnf = bnfinit(pol); \\very long...

pol = kummer(3,[2,3,5]);
abbnf = abelianbnfinit(pol);
getcyc(abbnf)
\\bnf = bnfinit(pol); bnf.cyc \\1 min 30s

pol = galoisgetpol(3^4,3)[1];
abbnf = abelianbnfinit(pol);
getcyc(abbnf)

pol = galoisgetpol(3^4,7)[1];
abbnf = abelianbnfinit(pol);
getcyc(abbnf)

pol = polcyclo(6552);
abbnf = abelianbnfinit(pol); \\degree 1728, 4h