Context Navigation

P.tst @ 2e49a5

spielwiese

Last change on this file since 2e49a5 was 2e49a5, checked in by Olaf Bachmann <obachman@…>, 23 years ago
* initial commit git-svn-id: file:///usr/local/Singular/svn/trunk@5001 2c84dea3-7e68-4137-9b89-c4e89433aadc
Property mode set to `100644`
File size: 9.7 KB

Line
1	LIB "tst.lib"; tst_init();
2	option(prot);
3
4	ring R = 0, x, lp;
5
6	matrix w[101][21];
7	w = 1, 1, 1, 12, 9, 22, 8, 32, 7, 42, 6, 52, 5, 62, 4, 72, 3, 82, 2, 92, 1 ,
8	2, 2, 1, 13, 9, 23, 8, 33, 7, 43, 6, 53, 5, 63, 4, 73, 3, 83, 2, 93, 1 ,
9	3, 3, 1, 14, 9, 24, 8, 34, 7, 44, 6, 54, 5, 64, 4, 74, 3, 84, 2, 94, 1 ,
10	4, 4, 1, 15, 9, 25, 8, 35, 7, 45, 6, 55, 5, 65, 4, 75, 3, 85, 2, 95, 1 ,
11	5, 5, 1, 16, 9, 26, 8, 36, 7, 46, 6, 56, 5, 66, 4, 76, 3, 86, 2, 96, 1 ,
12	6, 6, 1, 17, 9, 27, 8, 37, 7, 47, 6, 57, 5, 67, 4, 77, 3, 87, 2, 97, 1 ,
13	7, 7, 1, 18, 9, 28, 8, 38, 7, 48, 6, 58, 5, 68, 4, 78, 3, 88, 2, 98, 1 ,
14	8, 8, 1, 19, 9, 29, 8, 39, 7, 49, 6, 59, 5, 69, 4, 79, 3, 89, 2, 99, 1 ,
15	9, 9, 1, 20, 9, 30, 8, 40, 7, 50, 6, 60, 5, 70, 4, 80, 3, 90, 2, 100, 1 ,
16	10, 10, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
17	11, 2, 10, 11, 1, 22, 8, 32, 7, 42, 6, 52, 5, 62, 4, 72, 3, 82, 2, 92, 1 ,
18	12, 3, 10, 12, 1, 23, 8, 33, 7, 43, 6, 53, 5, 63, 4, 73, 3, 83, 2, 93, 1 ,
19	13, 4, 10, 13, 1, 24, 8, 34, 7, 44, 6, 54, 5, 64, 4, 74, 3, 84, 2, 94, 1 ,
20	14, 5, 10, 14, 1, 25, 8, 35, 7, 45, 6, 55, 5, 65, 4, 75, 3, 85, 2, 95, 1 ,
21	15, 6, 10, 15, 1, 26, 8, 36, 7, 46, 6, 56, 5, 66, 4, 76, 3, 86, 2, 96, 1 ,
22	16, 7, 10, 16, 1, 27, 8, 37, 7, 47, 6, 57, 5, 67, 4, 77, 3, 87, 2, 97, 1 ,
23	17, 8, 10, 17, 1, 28, 8, 38, 7, 48, 6, 58, 5, 68, 4, 78, 3, 88, 2, 98, 1 ,
24	18, 9, 10, 18, 1, 29, 8, 39, 7, 49, 6, 59, 5, 69, 4, 79, 3, 89, 2, 99, 1 ,
25	19, 10, 10, 19, 1, 30, 8, 40, 7, 50, 6, 60, 5, 70, 4, 80, 3, 90, 2, 100, 1 ,
26	20, 20, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
27	21, 2, 10, 12, 9, 21, 1, 32, 7, 42, 6, 52, 5, 62, 4, 72, 3, 82, 2, 92, 1 ,
28	22, 3, 10, 13, 9, 22, 1, 33, 7, 43, 6, 53, 5, 63, 4, 73, 3, 83, 2, 93, 1 ,
29	23, 4, 10, 14, 9, 23, 1, 34, 7, 44, 6, 54, 5, 64, 4, 74, 3, 84, 2, 94, 1 ,
30	24, 5, 10, 15, 9, 24, 1, 35, 7, 45, 6, 55, 5, 65, 4, 75, 3, 85, 2, 95, 1 ,
31	25, 6, 10, 16, 9, 25, 1, 36, 7, 46, 6, 56, 5, 66, 4, 76, 3, 86, 2, 96, 1 ,
32	26, 7, 10, 17, 9, 26, 1, 37, 7, 47, 6, 57, 5, 67, 4, 77, 3, 87, 2, 97, 1 ,
33	27, 8, 10, 18, 9, 27, 1, 38, 7, 48, 6, 58, 5, 68, 4, 78, 3, 88, 2, 98, 1 ,
34	28, 9, 10, 19, 9, 28, 1, 39, 7, 49, 6, 59, 5, 69, 4, 79, 3, 89, 2, 99, 1 ,
35	29, 10, 10, 20, 9, 29, 1, 40, 7, 50, 6, 60, 5, 70, 4, 80, 3, 90, 2, 100, 1 ,
36	30, 30, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
37	31, 2, 10, 12, 9, 22, 8, 31, 1, 42, 6, 52, 5, 62, 4, 72, 3, 82, 2, 92, 1 ,
38	32, 3, 10, 13, 9, 23, 8, 32, 1, 43, 6, 53, 5, 63, 4, 73, 3, 83, 2, 93, 1 ,
39	33, 4, 10, 14, 9, 24, 8, 33, 1, 44, 6, 54, 5, 64, 4, 74, 3, 84, 2, 94, 1 ,
40	34, 5, 10, 15, 9, 25, 8, 34, 1, 45, 6, 55, 5, 65, 4, 75, 3, 85, 2, 95, 1 ,
41	35, 6, 10, 16, 9, 26, 8, 35, 1, 46, 6, 56, 5, 66, 4, 76, 3, 86, 2, 96, 1 ,
42	36, 7, 10, 17, 9, 27, 8, 36, 1, 47, 6, 57, 5, 67, 4, 77, 3, 87, 2, 97, 1 ,
43	37, 8, 10, 18, 9, 28, 8, 37, 1, 48, 6, 58, 5, 68, 4, 78, 3, 88, 2, 98, 1 ,
44	38, 9, 10, 19, 9, 29, 8, 38, 1, 49, 6, 59, 5, 69, 4, 79, 3, 89, 2, 99, 1 ,
45	39, 10, 10, 20, 9, 30, 8, 39, 1, 50, 6, 60, 5, 70, 4, 80, 3, 90, 2, 100, 1 ,
46	40, 40, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
47	41, 2, 10, 12, 9, 22, 8, 32, 7, 41, 1, 52, 5, 62, 4, 72, 3, 82, 2, 92, 1 ,
48	42, 3, 10, 13, 9, 23, 8, 33, 7, 42, 1, 53, 5, 63, 4, 73, 3, 83, 2, 93, 1 ,
49	43, 4, 10, 14, 9, 24, 8, 34, 7, 43, 1, 54, 5, 64, 4, 74, 3, 84, 2, 94, 1 ,
50	44, 5, 10, 15, 9, 25, 8, 35, 7, 44, 1, 55, 5, 65, 4, 75, 3, 85, 2, 95, 1 ,
51	45, 6, 10, 16, 9, 26, 8, 36, 7, 45, 1, 56, 5, 66, 4, 76, 3, 86, 2, 96, 1 ,
52	46, 7, 10, 17, 9, 27, 8, 37, 7, 46, 1, 57, 5, 67, 4, 77, 3, 87, 2, 97, 1 ,
53	47, 8, 10, 18, 9, 28, 8, 38, 7, 47, 1, 58, 5, 68, 4, 78, 3, 88, 2, 98, 1 ,
54	48, 9, 10, 19, 9, 29, 8, 39, 7, 48, 1, 59, 5, 69, 4, 79, 3, 89, 2, 99, 1 ,
55	49, 10, 10, 20, 9, 30, 8, 40, 7, 49, 1, 60, 5, 70, 4, 80, 3, 90, 2, 100, 1 ,
56	50, 50, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
57	51, 2, 10, 12, 9, 22, 8, 32, 7, 42, 6, 51, 1, 62, 4, 72, 3, 82, 2, 92, 1 ,
58	52, 3, 10, 13, 9, 23, 8, 33, 7, 43, 6, 52, 1, 63, 4, 73, 3, 83, 2, 93, 1 ,
59	53, 4, 10, 14, 9, 24, 8, 34, 7, 44, 6, 53, 1, 64, 4, 74, 3, 84, 2, 94, 1 ,
60	54, 5, 10, 15, 9, 25, 8, 35, 7, 45, 6, 54, 1, 65, 4, 75, 3, 85, 2, 95, 1 ,
61	55, 6, 10, 16, 9, 26, 8, 36, 7, 46, 6, 55, 1, 66, 4, 76, 3, 86, 2, 96, 1 ,
62	56, 7, 10, 17, 9, 27, 8, 37, 7, 47, 6, 56, 1, 67, 4, 77, 3, 87, 2, 97, 1 ,
63	57, 8, 10, 18, 9, 28, 8, 38, 7, 48, 6, 57, 1, 68, 4, 78, 3, 88, 2, 98, 1 ,
64	58, 9, 10, 19, 9, 29, 8, 39, 7, 49, 6, 58, 1, 69, 4, 79, 3, 89, 2, 99, 1 ,
65	59, 10, 10, 20, 9, 30, 8, 40, 7, 50, 6, 59, 1, 70, 4, 80, 3, 90, 2, 100, 1 ,
66	60, 60, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
67	61, 2, 10, 12, 9, 22, 8, 32, 7, 42, 6, 52, 5, 61, 1, 72, 3, 82, 2, 92, 1 ,
68	62, 3, 10, 13, 9, 23, 8, 33, 7, 43, 6, 53, 5, 62, 1, 73, 3, 83, 2, 93, 1 ,
69	63, 4, 10, 14, 9, 24, 8, 34, 7, 44, 6, 54, 5, 63, 1, 74, 3, 84, 2, 94, 1 ,
70	64, 5, 10, 15, 9, 25, 8, 35, 7, 45, 6, 55, 5, 64, 1, 75, 3, 85, 2, 95, 1 ,
71	65, 6, 10, 16, 9, 26, 8, 36, 7, 46, 6, 56, 5, 65, 1, 76, 3, 86, 2, 96, 1 ,
72	66, 7, 10, 17, 9, 27, 8, 37, 7, 47, 6, 57, 5, 66, 1, 77, 3, 87, 2, 97, 1 ,
73	67, 8, 10, 18, 9, 28, 8, 38, 7, 48, 6, 58, 5, 67, 1, 78, 3, 88, 2, 98, 1 ,
74	68, 9, 10, 19, 9, 29, 8, 39, 7, 49, 6, 59, 5, 68, 1, 79, 3, 89, 2, 99, 1 ,
75	69, 10, 10, 20, 9, 30, 8, 40, 7, 50, 6, 60, 5, 69, 1, 80, 3, 90, 2, 100, 1 ,
76	70, 70, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
77	71, 2, 10, 12, 9, 22, 8, 32, 7, 42, 6, 52, 5, 62, 4, 71, 1, 82, 2, 92, 1 ,
78	72, 3, 10, 13, 9, 23, 8, 33, 7, 43, 6, 53, 5, 63, 4, 72, 1, 83, 2, 93, 1 ,
79	73, 4, 10, 14, 9, 24, 8, 34, 7, 44, 6, 54, 5, 64, 4, 73, 1, 84, 2, 94, 1 ,
80	74, 5, 10, 15, 9, 25, 8, 35, 7, 45, 6, 55, 5, 65, 4, 74, 1, 85, 2, 95, 1 ,
81	75, 6, 10, 16, 9, 26, 8, 36, 7, 46, 6, 56, 5, 66, 4, 75, 1, 86, 2, 96, 1 ,
82	76, 7, 10, 17, 9, 27, 8, 37, 7, 47, 6, 57, 5, 67, 4, 76, 1, 87, 2, 97, 1 ,
83	77, 8, 10, 18, 9, 28, 8, 38, 7, 48, 6, 58, 5, 68, 4, 77, 1, 88, 2, 98, 1 ,
84	78, 9, 10, 19, 9, 29, 8, 39, 7, 49, 6, 59, 5, 69, 4, 78, 1, 89, 2, 99, 1 ,
85	79, 10, 10, 20, 9, 30, 8, 40, 7, 50, 6, 60, 5, 70, 4, 79, 1, 90, 2, 100, 1 ,
86	80, 80, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
87	81, 2, 10, 12, 9, 22, 8, 32, 7, 42, 6, 52, 5, 62, 4, 72, 3, 81, 1, 92, 1 ,
88	82, 3, 10, 13, 9, 23, 8, 33, 7, 43, 6, 53, 5, 63, 4, 73, 3, 82, 1, 93, 1 ,
89	83, 4, 10, 14, 9, 24, 8, 34, 7, 44, 6, 54, 5, 64, 4, 74, 3, 83, 1, 94, 1 ,
90	84, 5, 10, 15, 9, 25, 8, 35, 7, 45, 6, 55, 5, 65, 4, 75, 3, 84, 1, 95, 1 ,
91	85, 6, 10, 16, 9, 26, 8, 36, 7, 46, 6, 56, 5, 66, 4, 76, 3, 85, 1, 96, 1 ,
92	86, 7, 10, 17, 9, 27, 8, 37, 7, 47, 6, 57, 5, 67, 4, 77, 3, 86, 1, 97, 1 ,
93	87, 8, 10, 18, 9, 28, 8, 38, 7, 48, 6, 58, 5, 68, 4, 78, 3, 87, 1, 98, 1 ,
94	88, 9, 10, 19, 9, 29, 8, 39, 7, 49, 6, 59, 5, 69, 4, 79, 3, 88, 1, 99, 1 ,
95	89, 10, 10, 20, 9, 30, 8, 40, 7, 50, 6, 60, 5, 70, 4, 80, 3, 89, 1, 100, 1 ,
96	90, 90, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
97	91, 2, 10, 12, 9, 22, 8, 32, 7, 42, 6, 52, 5, 62, 4, 72, 3, 82, 2, 91, 1 ,
98	92, 3, 10, 13, 9, 23, 8, 33, 7, 43, 6, 53, 5, 63, 4, 73, 3, 83, 2, 92, 1 ,
99	93, 4, 10, 14, 9, 24, 8, 34, 7, 44, 6, 54, 5, 64, 4, 74, 3, 84, 2, 93, 1 ,
100	94, 5, 10, 15, 9, 25, 8, 35, 7, 45, 6, 55, 5, 65, 4, 75, 3, 85, 2, 94, 1 ,
101	95, 6, 10, 16, 9, 26, 8, 36, 7, 46, 6, 56, 5, 66, 4, 76, 3, 86, 2, 95, 1 ,
102	96, 7, 10, 17, 9, 27, 8, 37, 7, 47, 6, 57, 5, 67, 4, 77, 3, 87, 2, 96, 1 ,
103	97, 8, 10, 18, 9, 28, 8, 38, 7, 48, 6, 58, 5, 68, 4, 78, 3, 88, 2, 97, 1 ,
104	98, 9, 10, 19, 9, 29, 8, 39, 7, 49, 6, 59, 5, 69, 4, 79, 3, 89, 2, 98, 1 ,
105	99, 10, 10, 20, 9, 30, 8, 40, 7, 50, 6, 60, 5, 70, 4, 80, 3, 90, 2, 99, 1 ,
106	100, 100, 1, 101, 1, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0, 1,0 ,
107	101, 1, 10, 11, 9, 21, 8, 31, 7, 41, 6, 51, 5, 61, 4, 71, 3, 81, 2, 91, 1 ;
108
109	matrix m[101][101];
110
111	int i, j;
112	for (i = 1; i <= 101; i++)
113	{
114	for (j = 1; j <= 10; j++)
115	{
116	m[i, int(w[i, 2j])] = w[i, 2j+1];
117	}
118	}
119
120	poly result, a;
121
122	tst_InitTimer();
123	result = det(module(m));
124	tst_ReportTimer("P");
125	result;
126
127	tst_InitTimer();
128	result = det(m - x*freemodule(101));
129	tst_ReportTimer("Q");
130
131	size(result);
132	lead(result);
133	subst(result, x, 2);
134
135	matrix m2[101][101];
136
137	for (i = 1; i <= 101; i++)
138	{
139	for (j = 1; j <= 101; j++)
140	{
141	m2[i, j] = m[i, j];
142	}
143	}
144	for (i = 1; i <= 101; i++)
145	{
146	a = m2[i, 1];
147	for (j = 1; j <= 100; j++)
148	{
149	m2[i, j] = m2[i, j+1];
150	}
151	m2[i, 101] = a;
152	}
153	for (i = 1; i <= 101; i++)
154	{
155	for (j = 1; j <= 101; j++)
156	{
157	if (m[i, j] <> 0)
158	{
159	m2[i, j] = m[i, j];
160	}
161	}
162	}
163
164
165	tst_InitTimer();
166	result = det(module(m2));
167	tst_ReportTimer("P'");
168	result;
169
170	tst_InitTimer();
171	/* this goes nowhere !
172	result = det(m2 - x*freemodule(101));
173	tst_ReportTimer("Q'");
174
175	size(result);
176	lead(result);
177	subst(result, x, 2);
178	*/
179	tst_OutTimer(100000, "Q': killed");
180
181	kill i,j, R;
182
183	tst_status(1); $

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format