 LIB "toric.lib";
ring r=0,(x,y,z),wp(3,2,1);
// call with toric ideal (of the matrix A=(1,1,1))
ideal I=xy,xz;
ideal J=toric_std(I);
J;
==> J[1]=yz
==> J[2]=xz
// call with the same ideal, but badly chosen generators:
// 1) not only binomials
I=xy,2xyz;
J=toric_std(I);
==> ? Generator 2 of the input ideal is no difference of monomials.
==> ? leaving toric.lib::toric_std
// 2) binomials whose monomials are not relatively prime
I=xy,xyyz,yz;
J=toric_std(I);
==> Warning: The monomials of generator 2 of the input ideal are not relative\
ly prime.
J;
==> J[1]=yz
==> J[2]=xz
// call with a nontoric ideal that seems to be toric
I=xyz,xyz;
J=toric_std(I);
J;
==> J[1]=y21
==> J[2]=xyz
// comparison with real standard basis and saturation
ideal H=std(I);
H;
==> H[1]=xyz
==> H[2]=y2zz
LIB "elim.lib";
sat(H,xyz);
==> [1]:
==> _[1]=xyz
==> _[2]=y21
==> [2]:
==> 1
