1 | // |
---|
2 | // solve_s.tst - short tests for |
---|
3 | |
---|
4 | LIB "tst.lib"; |
---|
5 | tst_init(); |
---|
6 | tst_ignore("CVS ID $Id$"); |
---|
7 | |
---|
8 | LIB "solve.lib"; |
---|
9 | |
---|
10 | //////////////////////////// |
---|
11 | |
---|
12 | //tst_ignore( "ring rs1 = 0,(x,y,z),lp;" ); |
---|
13 | //ring rs1 = 0,(h,x,y,z),lp; |
---|
14 | // |
---|
15 | //poly f0= x+y+z+1; |
---|
16 | //poly f1= x^2 + 6*x + 3*y + 6*z - 4; |
---|
17 | //poly f2= y^2 + 2*x -7*y + 5 + 2*z; |
---|
18 | //poly f3= x^2 + y^2 + z^2 - 1; |
---|
19 | //ideal i=f1,f2,f3; |
---|
20 | //ideal im=f0,f1,f2,f3; |
---|
21 | |
---|
22 | //def A=ures_solve(homog(i,h),1); |
---|
23 | //setring A; SOL; |
---|
24 | //kill A; |
---|
25 | //setring rs1; |
---|
26 | //mp_res_mat(homog(im,h),1); |
---|
27 | |
---|
28 | //tst_ignore( "ring rs2 = 0,(x,y,z),lp;" ); |
---|
29 | //ring rs2 = 0,(x,y,z),lp; |
---|
30 | // |
---|
31 | //poly f0= x+y+z; |
---|
32 | //poly f1= x^2 + 6*x + 3*y + 6*z - 4; |
---|
33 | //poly f2= y^2 + 2*x -7*y + 5 + 2*z; |
---|
34 | //poly f3= x^2 + y^2 + z^2 - 1; |
---|
35 | //ideal i=f1,f2,f3; |
---|
36 | //ideal im=f0,f1,f2,f3; |
---|
37 | // |
---|
38 | //def A=ures_solve(i); |
---|
39 | //setring A; SOL; |
---|
40 | //kill A; |
---|
41 | //setring rs2; |
---|
42 | //mp_res_mat(im); |
---|
43 | // |
---|
44 | //def rinC=fglm_solve(i,30); |
---|
45 | //setring rinC; |
---|
46 | //rlist; |
---|
47 | // |
---|
48 | //setring rs2; |
---|
49 | //kill rinC; |
---|
50 | //def rinC=triangL_solve(i,30); |
---|
51 | //setring rinC; rlist; |
---|
52 | // |
---|
53 | //setring rs2; |
---|
54 | //kill rinC; |
---|
55 | //def rinC=triangLf_solve(i,30); |
---|
56 | //setring rinC; rlist; |
---|
57 | // |
---|
58 | //setring rs2; |
---|
59 | //kill rinC; |
---|
60 | //def rinC=triangM_solve(i,30); |
---|
61 | //setring rinC; rlist; |
---|
62 | //kill rinC; |
---|
63 | |
---|
64 | //////////////////////////// |
---|
65 | |
---|
66 | tst_ignore( "ring rs3 = 0,(x1,x2,x3),lp;" ); |
---|
67 | ring rs3 = 0,(x0,x1,x2,x3),lp; |
---|
68 | |
---|
69 | poly f0= x1+x2+x3+1; |
---|
70 | poly f1= -1 + 6*x1^2 - 4*x1*x2 + 3*x2^2 - 8*x1*x3 + 10*x2*x3 + 13*x3^2; |
---|
71 | poly f2= -2*x1 - 2*x1^2 + x2 - 2*x2^2 + 3*x1 - 4*x1*x3 - 4*x2*x3 - 4*x3^2; |
---|
72 | poly f3= x1 - 5*x1^2 - x2 + 2*x1*x2 - 2*x2^2 + 8*x1*x3 - 10*x2*x3 - 13*x3^2; |
---|
73 | ideal i=f1,f2,f3; |
---|
74 | ideal im=f0,f1,f2,f3; |
---|
75 | |
---|
76 | def A=ures_solve(homog(i,x0),1); |
---|
77 | setring A; SOL; |
---|
78 | kill A; |
---|
79 | setring rs3; |
---|
80 | mp_res_mat(homog(im,x0),1); |
---|
81 | |
---|
82 | tst_ignore( "ring rs4 = 0,(x1,x2,x3),lp;" ); |
---|
83 | ring rs4 = 0,(x1,x2,x3),lp; |
---|
84 | |
---|
85 | poly f0= x1+x2+x3; |
---|
86 | poly f1= -1 + 6*x1^2 - 4*x1*x2 + 3*x2^2 - 8*x1*x3 + 10*x2*x3 + 13*x3^2; |
---|
87 | poly f2= -2*x1 - 2*x1^2 + x2 - 2*x2^2 + 3*x1 - 4*x1*x3 - 4*x2*x3 - 4*x3^2; |
---|
88 | poly f3= x1 - 5*x1^2 - x2 + 2*x1*x2 - 2*x2^2 + 8*x1*x3 - 10*x2*x3 - 13*x3^2; |
---|
89 | ideal i=f1,f2,f3; |
---|
90 | ideal im=f0,f1,f2,f3; |
---|
91 | |
---|
92 | def A=ures_solve(i); |
---|
93 | setring A; SOL; |
---|
94 | kill A; |
---|
95 | setring rs4; |
---|
96 | mp_res_mat(im); |
---|
97 | |
---|
98 | def rinC=fglm_solve(i,30); |
---|
99 | setring rinC; rlist; |
---|
100 | |
---|
101 | setring rs4; |
---|
102 | kill rinC; |
---|
103 | def rinC=triangL_solve(i,30); |
---|
104 | setring rinC;rlist; |
---|
105 | |
---|
106 | setring rs4; |
---|
107 | kill rinC; |
---|
108 | def rinC=triangLf_solve(i,30); |
---|
109 | setring rinC;rlist; |
---|
110 | |
---|
111 | setring rs4; |
---|
112 | kill rinC; |
---|
113 | def rinC=triangM_solve(i,30); |
---|
114 | setring rinC;rlist; |
---|
115 | |
---|
116 | //////////////////////////// |
---|
117 | |
---|
118 | //tst_ignore( "ring rs6= 0,(x,y),lp;" ); |
---|
119 | //ring rs6= 0,(x,y),lp; |
---|
120 | //poly f = 15x5 + x3 + x2 - 10; |
---|
121 | //laguerre_solve(f); |
---|
122 | // |
---|
123 | //laguerre_solve(f,10); |
---|
124 | // |
---|
125 | //tst_ignore( "ring rsc= (real,20,I),x,lp;" ); |
---|
126 | //ring rsc= (real,20,I),x,lp; |
---|
127 | //poly f = (15.4+I*5)*x^5 + (25.0e-2+I*2)*x^3 + x2 - 10*I; |
---|
128 | //laguerre_solve(f); |
---|
129 | // |
---|
130 | //////////////////////////// |
---|
131 | |
---|
132 | //ring rs7 = 0,(x),lp; |
---|
133 | //ideal v=16,0,11376,1046880,85949136; |
---|
134 | //interpolate( 3, v, 4 ); |
---|
135 | // |
---|
136 | //ring rs8 = 0,(x,y),dp; |
---|
137 | //ideal p = 2,3; |
---|
138 | //ideal v= 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16; |
---|
139 | //poly ip= interpolate( p,v,3 ); |
---|
140 | //ip; |
---|
141 | |
---|
142 | //////////////////////////// |
---|
143 | |
---|
144 | tst_status(1);$ |
---|
145 | |
---|