Context Navigation

gcdp_s.tst @ 4b81a90

spielwiese

Last change on this file since 4b81a90 was 4b81a90, checked in by Jens Schmidt <schmidt@…>, 26 years ago
* gcd_s.tst: recreated version from directory `tst' * gcd_s.tst: tests added for transcendental extensions and multivariate polynomials git-svn-id: file:///usr/local/Singular/svn/trunk@1414 2c84dea3-7e68-4137-9b89-c4e89433aadc
Property mode set to `100644`
File size: 5.1 KB

Line
1	// $Id: gcdp_s.tst,v 1.2 1998-04-21 15:59:58 schmidt Exp $
2
3	//
4	// gcdp_s.tst - short tests for gcd calculations mod p.
5	//
6	// All univariate non-trivial examples come from gcdUnivP.fex
7	// or gcdUnivPAlpha.fex.
8	//
9	// The example in variables `a', `c', `e', ... came from `untitled2(6/p/1)'
10	// in stdMultivPGcd.in.
11	//
12	// To Do:
13	//
14	// o multivariate gcd calculations with parameters
15	//
16
17	LIB "tst.lib";
18	tst_init();
19	tst_ignore("CVS ID : $Id: gcdp_s.tst,v 1.2 1998-04-21 15:59:58 schmidt Exp $");
20
21	//
22	// - ring r1=32003,x,dp.
23	//
24
25	tst_ignore( "ring r1=32003,x,dp;" );
26	ring r1=32003,x,dp;
27
28	poly f=-9554x^4-12895x^3-10023x^2-6213x;
29	poly g;
30
31	// some trivial examples
32	gcd(0, 0);
33	gcd(0, 3123);
34	gcd(4353, 0);
35
36	gcd(0, f);
37	gcd(f, 0);
38
39	gcd(23123, f);
40	gcd(f, 13123);
41
42	// some less trivial examples
43	f=-9554x^4-12895x^3-10023x^2-6213x;
44	g=-9554x^3-3341x^2+6213*x;
45	gcd(f, g);
46
47	f=-11265x^6+12161x^5+10369x^4-12161x^3+896*x^2;
48	g=10669x^8-10673x^7+5x^6+8x^5-10681x^4+10665x^3+7*x^2;
49	gcd(f, g);
50
51	f=-7918x^9-14406x^8-7256x^7+2092x^6-2198x^5-12539x^4-14631x^3-7150x^2;
52	g=-14833x^13-7011x^12+15121x^11-10864x^10+12943x^9-9871x^8+10354x^7-1437x^6+6604x^5-10394x^4-3231x^3-9348x^2-2092*x;
53	gcd(f, g);
54
55	//
56	// - ring r2=(32003,a),x,dp; minpoly=a^4+8734a^3+a^2+11817a+1.
57	//
58
59	tst_ignore( "ring r2=(32003,a),x,dp;" );
60	ring r2=(32003,a),x,dp;
61	minpoly=a^4+8734a^3+a^2+11817a+1;
62
63	poly f=-9554x^4-12895x^3-10023x^2-6213x;
64	poly g;
65
66	// first, some of the r1 examples
67	gcd(0, 0);
68	gcd(0, 3123);
69	gcd(4353, 0);
70
71	gcd(0, f);
72	gcd(f, 0);
73
74	gcd(23123, f);
75	gcd(f, 13123);
76
77	// some less trivial examples
78	f=-9554x^4-12895x^3-10023x^2-6213x;
79	g=-9554x^3-3341x^2+6213*x;
80	gcd(f, g);
81
82	f=-11265x^6+12161x^5+10369x^4-12161x^3+896*x^2;
83	g=10669x^8-10673x^7+5x^6+8x^5-10681x^4+10665x^3+7*x^2;
84	gcd(f, g);
85
86	// now, examples involving the algebraic variable
87	f=(21147a^3x^5+10147a^3x^4+8142a^2x^4+27671ax^4+29289x^4+29289a^3*x^3);
88	g=(6167a^3x^6+1536a^2x^6+5211ax^6+1536x^6+17534a^3x^5+478a^2x^5+1536ax^5+19085x^5+3203a^3x^4+26578a^2x^4+17052ax^4+26770x^4+26002a^3x^3+24062a^2x^3+192ax^3+6386x^3);
89	gcd(f, g);
90
91	f=(2295a^2x^11+5897a^3x^10+765a^2x^10+27718a^3x^9+1283a^2x^9+23792ax^9+518x^9+22026a^3x^8+11468a^2x^8+16834ax^8+24778x^8+28460a^3x^7+8140a^2x^7+28036ax^7+1613x^7+18196a^3x^6+13274a^2x^6+4638ax^6+31194x^6+12206a^3x^5+10349a^2x^5+30979ax^5+8612*x^5);
92	g=(8192a^3x^10+8957a^2x^10+9659a^3x^9+28672a^2x^9+4553ax^9+12288x^9+5819a^3x^8+320a^2x^8+14380ax^8+12104x^8+12047a^3x^7+19823a^2x^7+6824ax^7+4335x^7+29376a^3x^6+22190a^2x^6+31239ax^6+22647x^6);
93	gcd(f, g);
94
95	//
96	// - ring r3=(32003,t),x,dp.
97	//
98
99	tst_ignore( "ring r3=(32003,t),x,dp;" );
100	ring r3=(32003,t),x,dp;
101
102	poly f=(-9554x^4-12895x^3-10023x^2-6213x);
103	poly g;
104
105	// first, some of the r1 examples (slightly modified)
106	gcd(0, 0);
107	gcd(0, 3123*t);
108	gcd(4353, 0);
109
110	gcd(0, f/t);
111	gcd(f, 0);
112
113	gcd(23123/t, f);
114	gcd(f, 13123);
115
116	// some less trivial examples
117	f=-9554x^4-12895x^3-10023x^2-6213x;
118	g=-9554x^3-3341x^2+6213*x;
119	gcd(f, g);
120
121	// we go on with modified variable names
122	tst_ignore( "ring r3=(32003,a,c,e),(f,h,k),dp;" );
123	kill r3;
124	ring r3=(32003,a,c,e),(f,h,k),dp;
125
126	poly fPoly;
127	poly gPoly;
128
129	fPoly=(16001c^2f^2k^2-16001cefk^3+cf^3hk-ef^2h*k^2);
130	gPoly=(c^3efh-4ce^2f*h^2);
131	gcd(fPoly, gPoly);
132
133	fPoly=(c^4k^2-5c^3fhk+4c^2ehk^2+4c^2f^2h^2-4cefh^2*k);
134	gPoly=(-c^3fk^2+c^2ek^3+2c^2f^2hk-2cefh*k^2);
135	gcd(fPoly, gPoly);
136
137	fPoly=(-8001c^4efk^6+8001c^4f^3hk^4+8001c^3e^2k^7+8001c^3ef^2hk^5+16001c^3f^4h^2k^3+16001c^2e^2fhk^6-16001c^2ef^3h^2k^4);
138	gPoly=(16001c^7ehk^3-15998c^6efh^2k^2-7c^5ef^2h^3k-8c^4e^2fh^3k^2+4c^4ef^3h^4+8c^3e^3h^3k^3+24c^3e^2f^2h^4k-24c^2e^3fh^4k^2-16c^2e^2f^3h^5+16ce^3f^2h^5k);
139	gcd(fPoly, gPoly);
140
141	//
142	// - ring r4=32003,(t,x),dp.
143	//
144	// The examples from r4 are those from r3 with parameters
145	// changed to ring variables.
146	//
147
148	tst_ignore( "ring r4=32003,(t,x),dp;" );
149	ring r4=32003,(t,x),dp;
150
151	poly f=(-9554x^4-12895x^3-10023x^2-6213x);
152	poly g;
153
154	// first, some of the r1 examples (slightly modified)
155	gcd(0, 0);
156	gcd(0, 3123*t);
157	gcd(4353, 0);
158
159	gcd(0, f);
160	gcd(f, 0);
161
162	gcd(23123*t, f);
163	gcd(f, 13123);
164
165	// some less trivial examples
166	f=-9554x^4-12895x^3-10023x^2-6213x;
167	g=-9554x^3-3341x^2+6213*x;
168	gcd(t*f, g);
169
170	// we go on with modified variable names
171	tst_ignore( "ring r4=32003,(a,c,e,f,h,k),dp;" );
172	kill r4;
173	ring r4=32003,(a,c,e,f,h,k),dp;
174
175	poly fPoly;
176	poly gPoly;
177
178	fPoly=(16001c^2f^2k^2-16001cefk^3+cf^3hk-ef^2h*k^2);
179	gPoly=(c^3efh-4ce^2f*h^2);
180	gcd(fPoly, gPoly);
181
182	fPoly=(c^4k^2-5c^3fhk+4c^2ehk^2+4c^2f^2h^2-4cefh^2*k);
183	gPoly=(-c^3fk^2+c^2ek^3+2c^2f^2hk-2cefh*k^2);
184	gcd(fPoly, gPoly);
185
186	fPoly=(-8001c^4efk^6+8001c^4f^3hk^4+8001c^3e^2k^7+8001c^3ef^2hk^5+16001c^3f^4h^2k^3+16001c^2e^2fhk^6-16001c^2ef^3h^2k^4);
187	gPoly=(16001c^7ehk^3-15998c^6efh^2k^2-7c^5ef^2h^3k-8c^4e^2fh^3k^2+4c^4ef^3h^4+8c^3e^3h^3k^3+24c^3e^2f^2h^4k-24c^2e^3fh^4k^2-16c^2e^2f^3h^5+16ce^3f^2h^5k);
188	gcd(fPoly, gPoly);
189	$

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

Download in other formats:

Original Format