bug in primdecGTZ() for block ordering

[kroeker@…]

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

[hannes]

fixed by Gerhard with 96a18116acad728d793c95eaa3cf8b553eac13d9