? rnfkummer(bnrinit(bnfinit(y^6-2*y^5+8*y^4-18*y^3+24*y^2-18*y+6),168,1),matdiagonal([2,1,1,1,1,1,1,1]))

  ***   at top-level: rnfkummer(bnrinit(bn

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

  *** rnfkummer: elements not coprime in zlog_prk1:

    [1/3, 0, 0, 0, 1/3, -1/3]~

    [4, 0, 0, 2, 0, 0; 0, 4, 2, 0, 0, 0; 0, 0, 2, 0, 0, 0; 0, 0, 0, 2, 0, 0; 0, 0, 0, 0, 2, 0; 0, 0, 0, 0, 0, 2]