search

jakobkroeker / test_singular

manage a personal issue list for Singular CAS
0 stars 0 forks source link

minAssChar(), minAssGTZ(), primdecSY()... times out and eats mem #223

Open jakobkroeker opened 9 years ago

jakobkroeker commented 9 years ago

taken from /infiniteTests/log/test.minAssGTZ.nomp.minimize.02/bugs/id_3.16.outofmem.bug

LIB("primdec.lib");
ring rng = (0),(x,y,z),(lp,C);
ideal I = x^3+x*y^2*z+y*z^3,-1*x^3*y+1*x*z^2-1;
def result= minAssChar(I); // timeout + outofmem  - finishes with modifications
// in branch `spielwiese.testmode.mod`
result= minAssGTZ(I); // also bad

using 

LIB("primdec.lib");
ring R = (0),(x,y,z),(dp,C);
ideal I = x^3+x*y^2*z+y*z^3,-1*x^3*y+1*x*z^2-1;
def result= minAssChar(I); // finishes

ring S = (0),(x,y,z),(lp,C);
def result = fetch(R,result);
std(result[1]);

I get 

. std(result[1]);
_[1]=y8z3+y7z7-y6z9+2y5z4+2y4z8+3y4z6+y2z5-3y2z3+yz9+1
_[2]=xz22-2xz11+x+y7z16+2y7z5+y6z20-y6z9-y5z22-y5z13-2y5z11+y5z2+3y4z17+3y4z15+3y4z6+y3z21+y3z19-y3z10+5y3z8+y2z23-y2z14-7y2z12+y2z3+2yz18+yz7-3yz5-z20+4z9
_[3]=3xy-xz15+4xz4-y7z9-y6z13+3y6z2+y5z15+y5z6-3y4z10-3y4z8-y3z14-y3z12+3y3z3-y2z16+y2z7+9y2z5-2yz11+z13-6z2
_[4]=x2z11-x2-xy5z2-xy4z6-xy3z8+xy2z12-xy2z3-xyz7+xyz5-xz9-2y4z4-y3z8+y2z12
_[5]=x3+xy2z+yz3

looks ok since Macaulay2 claims that minimalPrimes are

x^3*y-x*z^2+1,
x^3+x*y^2*z+y*z^3,
-x^6+x^5+x*y*z^3+x*y*z-x*z^5-y*z+z^3

with gb

y^8*z^3+y^7*z^7-y^6*z^9+2*y^5*z^4+2*y^4*z^8+3*y^4*z^6+y^2*z^5-3*y^2*z^3+y*z^9+1,
x*z^22-2*x*z^11+x+y^7*z^16+2*y^7*z^5+y^6*z^20-y^6*z^9-y^5*z^22-y^5*z^13-2*y^5*z^11+y^5*z^2+3*y^4*z^17+3*y^4*z^15+3*y^4*z^6+y^3*z^21+y^3*z^19-y^3*z^10+5*y^3*z^8+y^2*z^23-y^2*z^14-7*y^2*z^12+y^2*z^3+2*y*z^18+y*z^7-3*y*z^5-z^20+4*z^9,
3*x*y-x*z^15+4*x*z^4-y^7*z^9-y^6*z^13+3*y^6*z^2+y^5*z^15+y^5*z^6-3*y^4*z^10-3*y^4*z^8-y^3*z^14-y^3*z^12+3*y^3*z^3-y^2*z^16+y^2*z^7+9*y^2*z^5-2*y*z^11+z^13-6*z^2,
3*x^2*z^11-3*x^2+2*x*z^20-2*x*z^9+2*y^7*z^14+2*y^6*z^18-3*y^6*z^7-2*y^5*z^20+y^5*z^11+6*y^4*z^15+3*y^4*z^13-3*y^4*z^4+2*y^3*z^19+2*y^3*z^17-3*y^3*z^8+2*y^2*z^21+4*y^2*z^12-6*y^2*z^10-3*y^2*z+4*y*z^16-3*y*z^5-2*z^18+3*z^7,
3*x^3-2*x*z^20+5*x*z^9-2*y^7*z^14-3*y^7*z^3-2*y^6*z^18+3*y^6*z^7+2*y^5*z^20+2*y^5*z^11+3*y^5*z^9-6*y^4*z^15-6*y^4*z^13-3*y^4*z^4-2*y^3*z^19-2*y^3*z^17+3*y^3*z^8-9*y^3*z^6-2*y^2*z^21+2*y^2*z^12+15*y^2*z^10-4*y*z^16+9*y*z^3+2*z^18-9*z^7-2*y^3*z^17+3*y^3*z^8-9*y^3*z^6-2*y^2*z^21+2*y^2*z^12+15*y^2*z^10-4*y*z^16+9*y*z^3+2*z^18-9*z^7

Macaulay2 input: 

R = QQ[x,y,z, MonomialOrder=>Lex]
I = ideal(x^3+x*y^2*z+y*z^3,-1*x^3*y+1*x*z^2-1)
mp = minimalPrimes I 
gens gb mp#0
jakobkroeker commented 9 years ago

another example: 

LIB("src/tests/test.primdec.mp.timeout");
ring rng = (0),(x,y,z),(lp(1),lp(1),lp(1),C);
ideal I = 10*x^2*z-13*x*y*z-6*x*y-12*y^2,-11*x^3-8*x*y+5*x-11;
printIdealInfo( I ) ;
// numVars     = 3;
// maxAbsCoeff = 13;
// maxDeg      = 3;
// maxTerms    = 4;
// numGens     = 2;
def result= minAssGTZ(I); // finishes in spielwiese.testmode.mod
result;
result = minAssChar(I); // timeout

// simplified: 
LIB("primdec.lib");
ring rng = (0),(x,y,z),lp;
ideal I = x^2*z-x*y*z-x*y-y^2, -x^3-x*y+x-1;
// printIdealInfo( I ) ;
// numVars     = 3;
// maxAbsCoeff = 1;
// maxDeg      = 3;
// maxTerms    = 4;
// numGens     = 2;
def result= minAssGTZ(I);  // finishes in branch `spielwiese.testmode.mod`
result;
result = minAssChar(I); // timeout