1 | LIB "tst.lib"; tst_init(); |
---|
2 | |
---|
3 | ring R = 0, (x1, x2, x3, x4, x5), lp; |
---|
4 | matrix w[26][11]; |
---|
5 | w = 1, 1, 1, 7, x4 , 12, x3 , 17, x2 , 22, x1, |
---|
6 | 2, 2, 1, 8, x4 , 13, x3 , 18, x2 , 23, x1, |
---|
7 | 3, 3, 1, 9, x4 , 14, x3 , 19, x2 , 24, x1, |
---|
8 | 4, 4, 1, 10, x4 , 15, x3 , 20, x2 , 25, x1, |
---|
9 | 5, 5, 1, 26, 1, 1, 0, 1, 0, 1, 0, |
---|
10 | 6, 2, x5 , 6, 1, 12, x3 , 17, x2 , 22, x1, |
---|
11 | 7, 3, x5 , 7, 1, 13, x3 , 18, x2 , 23, x1, |
---|
12 | 8, 4, x5 , 8, 1, 14, x3 , 19, x2 , 24, x1, |
---|
13 | 9, 5, x5 , 9, 1, 15, x3 , 20, x2 , 25, x1, |
---|
14 | 10, 10, 1, 26, 1, 1, 0, 1, 0, 1, 0, |
---|
15 | 11, 2, x5 , 7, x4 , 11, 1, 17, x2 , 22, x1, |
---|
16 | 12, 3, x5 , 8, x4 , 12, 1, 18, x2 , 23, x1, |
---|
17 | 13, 4, x5 , 9, x4 , 13, 1, 19, x2 , 24, x1, |
---|
18 | 14, 5, x5 , 10, x4 , 14, 1, 20, x2 , 25, x1, |
---|
19 | 15, 15, 1, 26, 1, 1, 0, 1, 0, 1, 0, |
---|
20 | 16, 2, x5 , 7, x4 , 12, x3 , 16, 1, 22, x1, |
---|
21 | 17, 3, x5 , 8, x4 , 13, x3 , 17, 1, 23, x1, |
---|
22 | 18, 4, x5 , 9, x4 , 14, x3 , 18, 1, 24, x1, |
---|
23 | 19, 5, x5 , 10, x4 , 15, x3 , 19, 1, 25, x1, |
---|
24 | 20, 20, 1, 26, 1, 1, 0, 1, 0, 1, 0, |
---|
25 | 21, 2, x5 , 7, x4 , 12, x3 , 17, x2 , 21, 1, |
---|
26 | 22, 3, x5 , 8, x4 , 13, x3 , 18, x2 , 22, 1, |
---|
27 | 23, 4, x5 , 9, x4 , 14, x3 , 19, x2 , 23, 1, |
---|
28 | 24, 5, x5 , 10, x4 , 15, x3 , 20, x2 , 24, 1, |
---|
29 | 25, 25, 1, 26, 1, 1, 0, 1, 0, 1, 0, |
---|
30 | 26, 1, x5 , 6, x4 , 11, x3 , 16, x2 , 21, x1; |
---|
31 | |
---|
32 | matrix m[26][26]; |
---|
33 | int i,j; |
---|
34 | for (i = 1; i <= 26; i++) |
---|
35 | { |
---|
36 | for (j = 1; j <= 5; j++) |
---|
37 | { |
---|
38 | m[i, int(w[i, 2*j])] = w[i, 2*j+1]; |
---|
39 | } |
---|
40 | } |
---|
41 | |
---|
42 | tst_InitTimer(); |
---|
43 | |
---|
44 | poly result = det(m); |
---|
45 | |
---|
46 | tst_ReportTimer("M1"); |
---|
47 | |
---|
48 | size(result); |
---|
49 | lead(result); |
---|
50 | |
---|
51 | map check = R, 2..nvars(R)+1; |
---|
52 | |
---|
53 | det(check(m)) == check(result); |
---|
54 | |
---|
55 | kill R, i, j; |
---|
56 | |
---|
57 | tst_status(1); $ |
---|
58 | |
---|
59 | |
---|