//
// test script for bareiss command
//
pagelength = 10000;
ring r = 0,(x,y,z),dp;
r;
matrix m1[3][3]=x,2,3,2,3x,4,1,3,2x;
m1;
list mm1=bareiss(m1);
mm1;
vector vdet=mm1[1][3];
poly pdet=det(m1);
if(leadcoef(pdet)!=leadcoef(vdet))
{
  vdet+pdet*gen(3);
}
else
{
  "error";
  vdet-pdet*gen(3);
}
"------------------------------------";
//
matrix m2[5][6]=0,2,1,-3,12,-9,1,2,5,3,6,-1,2,4,10,6,12,-2,-3,2,1;
list mm2=bareiss(m2);
print(mm2[1]);
"------------------------------------";
//
matrix m3[3][4]=xy,0,x2z,1,z2y,y2z21,0,1,zx3,1,0,z3y3;
list mm3=bareiss(m3);
mm3;
print(mm3[1]);
"-------------------------------------";
 //
poly s1 = 2x3y2+23x5y+104xy6+56z2x;
poly s2 = 3/2x3y2-2/3z6y4+12x2z4;
poly s3 = -2x4y5+x4z5;
poly s4 = 7/8x7y3z-1;
poly s5 = 4x2y3-2x2y2;
poly s6 = 1/2y2x3+4x2z2;
poly s7 = x+y+z;
poly s8 = x2+y2+z2;
poly s9 = xz+xy+yz;
matrix m4[3][3]=s1,0,s3,1,s5,0,s7,s8,s9;
m4;
list mm4=bareiss(m4);
mm4;
vector vdet=mm4[1][3];
poly pdet=det(m4);
if(leadcoef(pdet)!=leadcoef(vdet))
{
  vdet+pdet*gen(3);
}
else
{
  "error";
  vdet-pdet*gen(3);
}
"---------------------------------------";
//
matrix m5[3][2]=s8,s7,s9,s6,s5,0;
module mm=m5;
list mm5=bareiss(mm);
print(mm5[1]);
kill r;
"-------------------------------------";
ring r1=0,(x(1..9)),(dp);
matrix m6[3][3]=maxideal(1);
list mm6=bareiss(m6);
print(mm6[1]);
vector p1=mm6[1][3];
poly p2=det(m6);
if(leadcoef(p1)!=leadcoef(p2))
{
  p1+p2*gen(3);
}
else
{
  p1-p2*gen(3);
}
"--------------------------------------";
ring r2=0,(x(1..12)),ds;
matrix m7[3][4]= maxideal(1);
print(m7);
list mm7=bareiss(m7);
print(mm7[1]);
"--------------------------------------";
matrix m8[4][3]=maxideal(1);
print(m8);
list mm8=bareiss(m8);
print(mm8[1]);
"--------------------------------------";
matrix m9[4][4]=maxideal(2);
print(m9);
list mm9=bareiss(m9);
mm9;
"------------------------------------";
matrix m11[5][5]=maxideal(1),maxideal(1);
print(m11);
list mm11=bareiss(m11);
mm11;
"-------------------------------------";
listvar(all);
LIB "tst.lib";tst_status(1);$