There are still problems with radical: sometimes it gives an empty list instead of an ideal.
/**/ I := ideal(x^2, x*y+1); /**/ radical(I); 
#3 Updated by John Abbott 4 months ago
Here is another case where things go wrong:
/**/ use QQ[x,y,z,t]; /**/ I := ideal(2*x*y*z +2*y*t, -y^3 -x*y*t +x*t^2, x^2*z +2*x*y*z +y^2*t); /**/ radical(I); --> ERROR: I don't know how to evaluate operator + between IDEAL and LIST --> WHERE: at line 1105 (column 19) of radical.cpkg5
#6 Updated by Anna Maria Bigatti 4 months ago
- Status changed from New to Closed
- Assignee set to Anna Maria Bigatti
- % Done changed from 0 to 100
- Estimated time set to 1.51 h
Fixed: all the problem was in ideal(1).
(there is a naughty trick in
I immediately if IsOne(I).