bug in radicalZ or primdecZ

[kroeker@…]

Hello,

there is either a bug in radicalZ or primdecZ.

Example:

LIB("primdecint.lib");
ring rng = integer,(x, y, z),dp;
ideal j = 2*x, y*z+4, x*z, 2*y^2-6*z^2, 3*z^3+4*y;

def rJ = radicalZ(j);

primdecZ(rJ);

primdecZ output:

[1]:
   [1]:
      _[1]=x
      _[2]=y*z+4
      _[3]=y^2-3*z^2
      _[4]=3*z^3+4*y
   [2]:
      _[1]=x
      _[2]=y*z+4
      _[3]=y^2-3*z^2
      _[4]=3*z^3+4*y
[2]:
   [1]:
      _[1]=2
      _[2]=z
   [2]:
      _[1]=2
      _[2]=z
[3]:
   [1]:
      _[1]=2
      _[2]=y
      _[3]=x
      _[4]=z^3
   [2]:
      _[1]=2
      _[2]=z
      _[3]=y
      _[4]=x

[hannes]

fixed by adi (a33be22758b3f3bf40bc27d6fc504cb2487a4322)