1 | // |
---|
2 | // solve_s.tst - short tests for |
---|
3 | |
---|
4 | LIB "tst.lib"; |
---|
5 | tst_init(); |
---|
6 | |
---|
7 | LIB "solve.lib"; |
---|
8 | |
---|
9 | //////////////////////////// |
---|
10 | |
---|
11 | tst_ignore( "ring rs1 = 0,(x,y,z),lp;" ); |
---|
12 | ring rs1 = 0,(h,x,y,z),lp; |
---|
13 | |
---|
14 | poly f0= x+y+z+1; |
---|
15 | poly f1= x^2 + 6*x + 3*y + 6*z - 4; |
---|
16 | poly f2= y^2 + 2*x -7*y + 5 + 2*z; |
---|
17 | poly f3= x^2 + y^2 + z^2 - 1; |
---|
18 | ideal i=f1,f2,f3; |
---|
19 | ideal im=f0,f1,f2,f3; |
---|
20 | |
---|
21 | def A=ures_solve(homog(i,h),1); |
---|
22 | setring A; SOL; |
---|
23 | kill A; |
---|
24 | setring rs1; |
---|
25 | mp_res_mat(homog(im,h),1); |
---|
26 | |
---|
27 | tst_ignore( "ring rs2 = 0,(x,y,z),lp;" ); |
---|
28 | ring rs2 = 0,(x,y,z),lp; |
---|
29 | |
---|
30 | poly f0= x+y+z; |
---|
31 | poly f1= x^2 + 6*x + 3*y + 6*z - 4; |
---|
32 | poly f2= y^2 + 2*x -7*y + 5 + 2*z; |
---|
33 | poly f3= x^2 + y^2 + z^2 - 1; |
---|
34 | ideal i=f1,f2,f3; |
---|
35 | ideal im=f0,f1,f2,f3; |
---|
36 | |
---|
37 | def A=ures_solve(i); |
---|
38 | setring A; SOL; |
---|
39 | kill A; |
---|
40 | setring rs2; |
---|
41 | mp_res_mat(im); |
---|
42 | |
---|
43 | def rinC = fglm_solve(i,30); |
---|
44 | setring rinC; |
---|
45 | rlist; |
---|
46 | |
---|
47 | setring rs2; |
---|
48 | kill rinC; |
---|
49 | def rinC = triangL_solve(i,30); |
---|
50 | setring rinC; |
---|
51 | rlist; |
---|
52 | |
---|
53 | setring rs2; |
---|
54 | kill rinC; |
---|
55 | def rinC = triangLf_solve(i,30); |
---|
56 | setring rinC; |
---|
57 | rlist; |
---|
58 | |
---|
59 | setring rs2; |
---|
60 | kill rinC; |
---|
61 | def rinC = triangM_solve(i,30); |
---|
62 | setring rinC; |
---|
63 | rlist; |
---|
64 | kill rinC; |
---|
65 | |
---|
66 | //////////////////////////// |
---|
67 | |
---|
68 | tst_ignore( "ring rs6= 0,(x,y),lp;" ); |
---|
69 | ring rs6= 0,(x,y),lp; |
---|
70 | poly f = 15x5 + x3 + x2 - 10; |
---|
71 | laguerre_solve(f); |
---|
72 | |
---|
73 | laguerre_solve(f,10); |
---|
74 | |
---|
75 | tst_ignore( "ring rsc= (real,20,I),x,lp;" ); |
---|
76 | ring rsc= (real,20,I),x,lp; |
---|
77 | poly f = (15.4+I*5)*x^5 + (25.0e-2+I*2)*x^3 + x2 - 10*I; |
---|
78 | laguerre_solve(f); |
---|
79 | |
---|
80 | //////////////////////////// |
---|
81 | |
---|
82 | ring rs7 = 0,(x),lp; |
---|
83 | ideal v=16,0,11376,1046880,85949136; |
---|
84 | interpolate( 3, v, 4 ); |
---|
85 | |
---|
86 | ring rs8 = 0,(x,y),dp; |
---|
87 | ideal p = 2,3; |
---|
88 | ideal v= 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16; |
---|
89 | poly ip= interpolate( p,v,3 ); |
---|
90 | ip; |
---|
91 | |
---|
92 | //////////////////////////// |
---|
93 | |
---|
94 | example simplexOut; |
---|
95 | |
---|
96 | tst_status(1);$ |
---|
97 | |
---|