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; |
---|
52 | rlist; |
---|
53 | |
---|
54 | setring rs2; |
---|
55 | kill rinC; |
---|
56 | def rinC = triangLf_solve(i,30); |
---|
57 | setring rinC; |
---|
58 | rlist; |
---|
59 | |
---|
60 | setring rs2; |
---|
61 | kill rinC; |
---|
62 | def rinC = triangM_solve(i,30); |
---|
63 | setring rinC; |
---|
64 | rlist; |
---|
65 | kill rinC; |
---|
66 | |
---|
67 | //////////////////////////// |
---|
68 | |
---|
69 | tst_ignore( "ring rs6= 0,(x,y),lp;" ); |
---|
70 | ring rs6= 0,(x,y),lp; |
---|
71 | poly f = 15x5 + x3 + x2 - 10; |
---|
72 | laguerre_solve(f); |
---|
73 | |
---|
74 | laguerre_solve(f,10); |
---|
75 | |
---|
76 | tst_ignore( "ring rsc= (real,20,I),x,lp;" ); |
---|
77 | ring rsc= (real,20,I),x,lp; |
---|
78 | poly f = (15.4+I*5)*x^5 + (25.0e-2+I*2)*x^3 + x2 - 10*I; |
---|
79 | laguerre_solve(f); |
---|
80 | |
---|
81 | //////////////////////////// |
---|
82 | |
---|
83 | ring rs7 = 0,(x),lp; |
---|
84 | ideal v=16,0,11376,1046880,85949136; |
---|
85 | interpolate( 3, v, 4 ); |
---|
86 | |
---|
87 | ring rs8 = 0,(x,y),dp; |
---|
88 | ideal p = 2,3; |
---|
89 | ideal v= 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16; |
---|
90 | poly ip= interpolate( p,v,3 ); |
---|
91 | ip; |
---|
92 | |
---|
93 | //////////////////////////// |
---|
94 | |
---|
95 | example simplexOut; |
---|
96 | |
---|
97 | tst_status(1);$ |
---|
98 | |
---|