here is finally an example where 'idealsEqual()' fails

[jakobkroeker]

// test: test.std.correct.random.03

LIB "primdec.lib";
ring rng = (26053),(x,y,z),(Ws(1,-3,-5),C);
short = 0 ;
ideal I = 66*x*y^2*z^2+99*x^2*z^3+100*y+62*x,
          24*x*y^2*z^2+101*y^3*z-68*x*y^2*z+20*x^2*z^2;

ideal gI =  std(I);
ideal ggI =  std(gI);
lead(gI);
lead(ggI);
gI;
ggI;
idealsEqual(I,gI );
reduce(I,gI);
reduce(gI,std(gI));

output

> lead(gI);
_[1]=y^6
_[2]=x*y^2*z^2
_[3]=y^3*z
_[4]=x^3*z^3
> lead(ggI);
_[1]=x*y^5*z
_[2]=y^6
_[3]=x*y^2*z^2
_[4]=y^3*z
_[5]=x^3*z^3
> gI;
gI[1]=y^6+3353*x^2*y^3*z^2-4160*x^2*y^4*z-1945*x*y^5-6114*x^3*y^2*z^2+1032*x^2*y^3*z-4753*x^4*z^3-3426*x^2*y^4-6655*x^3*y^2*z-9239*x*y^3-4165*x^2*y^2
gI[2]=x*y^2*z^2+7603*y^3*z-4345*x*y^2*z+4343*x^2*z^2
gI[3]=y^3*z+3353*x^2*z^3+12123*x*y^2*z+1032*x^2*z^2+5229*y+12100*x
gI[4]=x^3*z^3-12683*y^4+3955*x^2*y^2*z-3756*x*y^3+5789*x^3*z^2-10250*x^2*y*z-4446*x^2*y^2-6355*x^3*z-8157*x*y-10789*x^2
> ggI;
ggI[1]=x*y^5*z+7603*y^6+11985*x^2*y^4*z+10269*x*y^5+4343*x^2*y^3*z+5122*x^2*y^4-3039*x^3*y^2*z
ggI[2]=y^6+3353*x^2*y^3*z^2-4160*x^2*y^4*z-1945*x*y^5-6114*x^3*y^2*z^2+1032*x^2*y^3*z-4753*x^4*z^3-3426*x^2*y^4-6655*x^3*y^2*z-9239*x*y^3-4165*x^2*y^2
ggI[3]=x*y^2*z^2+7603*y^3*z-4345*x*y^2*z+4343*x^2*z^2
ggI[4]=y^3*z+3353*x^2*z^3+12123*x*y^2*z+1032*x^2*z^2+5229*y+12100*x
ggI[5]=x^3*z^3-12683*y^4+3955*x^2*y^2*z-3756*x*y^3+5789*x^3*z^2-10250*x^2*y*z-4446*x^2*y^2-6355*x^3*z-8157*x*y-10789*x^2
> idealsEqual(I,gI );
0
> reduce(I,gI);
_[1]=0
_[2]=0
> reduce(gI,std(gI));
_[1]=0
_[2]=0
_[3]=0
_[4]=0

[jakobkroeker]

update: seems fixed now (which commit?) size(ggI) is 4, and ggi is

ggI[1]=y^6+12287*x^4*z^4-4160*x^2*y^4*z-1945*x*y^5-11853*x^3*y^2*z^2+1032*x^2*y^3*z-3426*x^2*y^4-6655*x^3*y^2*z-9239*x*y^3+832*x^2*y*z-4165*x^2*y^2-6779*x^3*z
ggI[2]=x*y^2*z^2-13025*x^2*z^3+791*y-3157*x
ggI[3]=y^3*z+3353*x^2*z^3+12123*x*y^2*z+1032*x^2*z^2+5229*y+12100*x
ggI[4]=x^3*z^3-12683*y^4+3955*x^2*y^2*z-3756*x*y^3+5789*x^3*z^2-10250*x^2*y*z-4446*x^2*y^2-6355*x^3*z-8157*x*y-10789*x^2

lead(ggI):

> lead(ggI);
_[1]=y^6
_[2]=x*y^2*z^2
_[3]=y^3*z
_[4]=x^3*z^3

[jakobkroeker]

bugfix was merged