[9542e4] | 1 | // |
---|
| 2 | // solve_s.tst - short tests for |
---|
| 3 | |
---|
| 4 | LIB "tst.lib"; |
---|
| 5 | tst_init(); |
---|
[341696] | 6 | tst_ignore("CVS ID $Id$"); |
---|
[9542e4] | 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 | |
---|
[837fed] | 22 | def A=ures_solve(homog(i,h),1); |
---|
| 23 | setring A; SOL; |
---|
| 24 | kill A; |
---|
| 25 | setring rs1; |
---|
[9542e4] | 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 | |
---|
[837fed] | 38 | def A=ures_solve(i); |
---|
| 39 | setring A; SOL; |
---|
| 40 | kill A; |
---|
| 41 | setring rs2; |
---|
[9542e4] | 42 | mp_res_mat(im); |
---|
| 43 | |
---|
[6ed15c] | 44 | def rinC = fglm_solve(i,30); |
---|
| 45 | setring rinC; |
---|
[9542e4] | 46 | rlist; |
---|
| 47 | |
---|
| 48 | setring rs2; |
---|
[6ed15c] | 49 | kill rinC; |
---|
| 50 | def rinC = triangL_solve(i,30); |
---|
| 51 | setring rinC; |
---|
[9542e4] | 52 | rlist; |
---|
| 53 | |
---|
| 54 | setring rs2; |
---|
[6ed15c] | 55 | kill rinC; |
---|
| 56 | def rinC = triangLf_solve(i,30); |
---|
| 57 | setring rinC; |
---|
[9542e4] | 58 | rlist; |
---|
| 59 | |
---|
| 60 | setring rs2; |
---|
[6ed15c] | 61 | kill rinC; |
---|
| 62 | def rinC = triangM_solve(i,30); |
---|
| 63 | setring rinC; |
---|
[9542e4] | 64 | rlist; |
---|
[6ed15c] | 65 | kill rinC; |
---|
[9542e4] | 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 | |
---|
[5b7df7] | 95 | example simplexOut; |
---|
| 96 | |
---|
[9542e4] | 97 | tst_status(1);$ |
---|
| 98 | |
---|