Opened 6 years ago

Last modified 6 years ago

#793 new bug

resolve + intersectionDiv fails

Reported by: ren Owned by: anne
Priority: minor Milestone: 4-2-0 and higher
Component: dontKnow Version: 4-1-0
Keywords: Cc:

Description (last modified by ren)

In the following example, 'intersectionDiv' complains because something is not equidimensional (sorry, but I understand too little to determine whether it is an error in 'resolve' or in 'intersectionDiv').

LIB "resolve.lib";
LIB "reszeta.lib";

ring r = 0,(x,y,z),dp;
ideal I = x^3+y^3+z^2*(x+y)+z^5;

list RI = resolve(I);
list DI = intersectionDiv(RI);

And here is another example where it fails albeit for different reasons:

LIB "resolve.lib";
LIB "reszeta.lib";

ring r = 0,(x,y,z),dp;
ideal I = x7+y5+z5+x3yz+x3z2;

list RI = resolve(I);
list DI = intersectionDiv(RI);

Change History (3)

comment:1 Changed 6 years ago by ren

Description: modified (diff)

comment:2 Changed 6 years ago by anne

The bug is in intersectionDiv and not in resolve.

We are already aware of a bug in this case, where algebraic extensions of QQ appear during the identification of exceptional divisors. A rewrite of intersectionDiv has been given out as a master thesis and will hopefully lead to a more robust implementation.

comment:3 Changed 6 years ago by hannes

Owner: changed from somebody to anne
Note: See TracTickets for help on using tickets.