Opened 8 years ago

Closed 8 years ago

# bug in primdecGTZ() with minpoly, char=5

Reported by: Owned by: kroeker@… pfister minor 4-1-0 and higher singular-libs 4-0-2

### Description

Example:

```LIB("primdec.lib");
system("random",1580702082);
ring rng = (5,vv),(x,y,z),dp;
minpoly = (vv^2-2);
ideal I = (-2*vv)*x*z+(-vv)*y^3,x*y*z-2*x*y+2*y^3,-x+z^2+(-2*vv);
list  PL,ML;
list L1 =  primdecGTZ (I);
list L2 = primdecSY (I);
L1;
L2;

testPrimary(L1 , I) ;  // fails
testPrimary(L2 , I) ;  //ok

idealsEqual(L1[1][2],L2[3][2]); //ok
idealsEqual(L1[1][1],L2[3][1]); //ok
idealsEqual(L1[3][2],L2[1][2]); //ok
idealsEqual(L1[3][1],L2[1][1]); //ok
idealsEqual(L1[2][2],L2[2][2]); //ok
idealsEqual(L1[2][1],L2[2][1]); //ok

// here is the issue:
idealsEqual(L1[4][2],L2[4][2]); //ok
idealsEqual(L1[4][1],L2[4][1]); //=0
```

output

```> L1[4][1]; //wrong
_[1]=y
_[2]=x
_[3]=z^2-x+(-2*vv)

> L2[4][1];
_[1]=x
_[2]=z^2+(-2*vv)
_[3]=y^3
```

### comment:1 Changed 8 years ago by Oleksandr

Owner: changed from somebody to pfister

### comment:2 Changed 8 years ago by hannes

Component: dontKnow → singular-libs → fixed new → closed
Note: See TracTickets for help on using tickets.