| 1 | // |
|---|
| 2 | // test script for bareiss command |
|---|
| 3 | // |
|---|
| 4 | pagelength = 10000; |
|---|
| 5 | LIB "lib0"; |
|---|
| 6 | ring r = 32003,(x,y,z),dp; |
|---|
| 7 | r; |
|---|
| 8 | matrix m1[3][3]=1,2,3,2,3,4,1,3,2; |
|---|
| 9 | m1; |
|---|
| 10 | matrix mm1[3][3]=bareiss(m1); |
|---|
| 11 | mm1; |
|---|
| 12 | mm1[3,3]-det(m1); |
|---|
| 13 | "------------------------------------"; |
|---|
| 14 | // |
|---|
| 15 | 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; |
|---|
| 16 | m2; |
|---|
| 17 | matrix mm2[5][6]=bareiss(m2); |
|---|
| 18 | mm2; |
|---|
| 19 | "------------------------------------"; |
|---|
| 20 | // |
|---|
| 21 | matrix m3[3][4]=xy,0,x2z,1,z2y,y2z21,0,1,zx3,1,0,z3y3; |
|---|
| 22 | m3; |
|---|
| 23 | matrix mm3[3][4]=bareiss(m3); |
|---|
| 24 | mm3; |
|---|
| 25 | "-------------------------------------"; |
|---|
| 26 | // |
|---|
| 27 | poly s1 = 2x3y2+23x5y+104xy6+56z2x; |
|---|
| 28 | poly s2 = 3/2x3y2-2/3z6y4+12x2z4; |
|---|
| 29 | poly s3 = -2x4y5+x4z5; |
|---|
| 30 | poly s4 = 7/8x7y3z-1; |
|---|
| 31 | poly s5 = 4x2y3-2x2y2; |
|---|
| 32 | poly s6 = 1/2y2x3+4x2z2; |
|---|
| 33 | poly s7 = x+y+z; |
|---|
| 34 | poly s8 = x2+y2+z2; |
|---|
| 35 | poly s9 = xz+xy+yz; |
|---|
| 36 | matrix m4[3][3]=s1,0,s3,1,s5,0,s7,s8,s9; |
|---|
| 37 | m4; |
|---|
| 38 | matrix mm4[3][3]=bareiss(m4); |
|---|
| 39 | mm4; |
|---|
| 40 | mm4[3,2]-det(m4); |
|---|
| 41 | "---------------------------------------"; |
|---|
| 42 | // |
|---|
| 43 | matrix m5[3][2]=s8,s7,s9,s6,s5,0; |
|---|
| 44 | m5; |
|---|
| 45 | matrix mm5[3][2]=bareiss(m5); |
|---|
| 46 | mm5; |
|---|
| 47 | "-------------------------------------"; |
|---|
| 48 | ring r1=0,(x(1..9)),(dp); |
|---|
| 49 | matrix m6[3][3]=maxideal(1); |
|---|
| 50 | pmat(m6); |
|---|
| 51 | matrix mm6=bareiss(m6); |
|---|
| 52 | pmat(mm6); |
|---|
| 53 | mm6[3,3]-det(m6); |
|---|
| 54 | "--------------------------------------"; |
|---|
| 55 | ring r2=0,(x(1..12)),ds; |
|---|
| 56 | matrix m7[3][4]= maxideal(1); |
|---|
| 57 | pmat(m7); |
|---|
| 58 | m7=bareiss(m7); |
|---|
| 59 | pmat(m7,8); |
|---|
| 60 | matrix m77=m7; |
|---|
| 61 | m77[1,1]=x(1)^2; |
|---|
| 62 | m77; |
|---|
| 63 | m77=bareiss(m77); |
|---|
| 64 | m77; |
|---|
| 65 | "--------------------------------------"; |
|---|
| 66 | matrix m8[4][3]=maxideal(1); |
|---|
| 67 | pmat(m8); |
|---|
| 68 | m8=bareiss(m8); |
|---|
| 69 | pmat(m8,8); |
|---|
| 70 | "--------------------------------------"; |
|---|
| 71 | matrix m9[4][4]=maxideal(2); |
|---|
| 72 | pmat(m9); |
|---|
| 73 | m9=bareiss(m9); |
|---|
| 74 | pmat(m9,8); |
|---|
| 75 | "------------------------------------"; |
|---|
| 76 | matrix m10[4][4]=maxideal(1); |
|---|
| 77 | pmat(m10); |
|---|
| 78 | m10=bareiss(m10); |
|---|
| 79 | pmat(m10,8); |
|---|
| 80 | "------------------------------------"; |
|---|
| 81 | matrix m11[5][5]=maxideal(1),maxideal(1); |
|---|
| 82 | pmat(m11); |
|---|
| 83 | matrix mm11=bareiss(m11); |
|---|
| 84 | mm11; |
|---|
| 85 | pmat(mm11,8); |
|---|
| 86 | mm11[2,2]-det(m11); |
|---|
| 87 | "-------------------------------------"; |
|---|
| 88 | listvar(all); |
|---|
| 89 | kill r,r1,r2; |
|---|
| 90 | $; |
|---|
