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