Opened 9 years ago

Closed 9 years ago

# bug in primdecGTZ() for block ordering

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

### Description

Example:

```LIB("primdec.lib");

ring rng = 0,(xt, xd, xb ,xy, xg),( dp(2),dp(3) );

ideal I = 39*xt*xy+13*xd*xb, 32*xt^3*xb-35*xt*xy*xg^2;

def L2 = primdecSY(I);
L2;
testPrimary(L2, I);  //ok

def L1 = primdecGTZ(I);
L1;
testPrimary(L1, I); // L1 wrong?
```

output:

```> L2;
[1]:
[1]:
_[1]=xd
_[2]=xt
[2]:
_[1]=xd
_[2]=xt
[2]:
[1]:
_[1]=xy
_[2]=xb
[2]:
_[1]=xy
_[2]=xb
[3]:
[1]:
_[1]=xb
_[2]=xt
[2]:
_[1]=xb
_[2]=xt
[4]:
[1]:
_[1]=3*xt*xy+xd*xb
_[2]=32*xd^2*xb^3-315*xy^3*xg^2
_[3]=32*xt*xd*xb^2+105*xy^2*xg^2
_[4]=32*xt^2*xb-35*xy*xg^2
_[5]=96*xt^3+35*xd*xg^2
[2]:
_[1]=3*xt*xy+xd*xb
_[2]=32*xd^2*xb^3-315*xy^3*xg^2
_[3]=32*xt*xd*xb^2+105*xy^2*xg^2
_[4]=32*xt^2*xb-35*xy*xg^2
_[5]=96*xt^3+35*xd*xg^2

> L1;
[1]:
[1]:
_[1]=xy
_[2]=xb
[2]:
_[1]=xy
_[2]=xb
[2]:
[1]:
_[1]=32*xd^2*xb^3-315*xy^3*xg^2
_[2]=3*xt*xy+xd*xb
_[3]=32*xt*xd*xb^2+105*xy^2*xg^2
_[4]=32*xt^2*xb-35*xy*xg^2
_[5]=96*xt^3+35*xd*xg^2
[2]:
_[1]=32*xd^2*xb^3-315*xy^3*xg^2
_[2]=3*xt*xy+xd*xb
_[3]=32*xt*xd*xb^2+105*xy^2*xg^2
_[4]=32*xt^2*xb-35*xy*xg^2
_[5]=96*xt^3+35*xd*xg^2
[3]:
[1]:
_[1]=xd
_[2]=xt
[2]:
_[1]=xd
_[2]=xt

```

### comment:1 Changed 9 years ago by hannes

Resolution: → fixed new → closed
Note: See TracTickets for help on using tickets.