idealintersect(nfinit(x^2+1),0,1)
Q = nfinit(y);
idealaddtoone(Q,[1,[;]])
idealaddtoone(Q,[1,[;],Mat(1/2)])
ideallog(Q,2,idealstar(Q,4,1))
idealstar(nfinit(y^2+1),matdiagonal([6,12]))

