primdecGTZ() example with incorrect radical entry

[kroeker@…]

LIB("primdec.lib");
ring rng = (2),(x,y,z),dp;
ideal I = x^2*y+x*z^2+1, x^2*z+z^2, z^3+1;
system("random",1426357199);
def L2 = primdecGTZ(I);

//here is the problem
idealsEqual( radical( L2[2][1] ),  L2[2][2]  ); // 0 
L2[2][2] ;
radical( L2[2][1] );

output

> ring rng = (2),(x,y,z),dp;
> ideal I = x^2*y+x*z^2+1, x^2*z+z^2, z^3+1;
> system("random",1426357199);
> def L2 = primdecGTZ(I);

> 
. //here is the problem
. idealsEqual( radical( L2[2][1] ),  L2[2][2]  ); // 0 
0
> L2[2][2] ;
_[1]=z2+z+1
_[2]=y2+z+1
_[3]=yz2+x+z
> radical( L2[2][1] );
_[1]=y+z
_[2]=x+z+1
_[3]=z2+z+1

[hannes]

fixed with ad75673d46dfdfa95bfef0dfdad5d5e3c7f7de3b, c513741db768d8779fd1b824a7737684ae3abad9