| 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 | |
|---|
