//
// solve_s.tst - short tests for 

LIB "tst.lib";
tst_init();

LIB "solve.lib";

////////////////////////////

tst_ignore( "ring rs1 = 0,(x,y,z),lp;" );
ring rs1 = 0,(h,x,y,z),lp;

poly f0= x+y+z+1;
poly f1= x^2 + 6*x + 3*y + 6*z - 4;
poly f2= y^2 + 2*x -7*y + 5 + 2*z;
poly f3= x^2 + y^2 + z^2 - 1;
ideal i=f1,f2,f3;
ideal im=f0,f1,f2,f3;

def A=ures_solve(homog(i,h),1);
setring A; SOL;
kill A;
setring rs1;
mp_res_mat(homog(im,h),1);

tst_ignore( "ring rs2 = 0,(x,y,z),lp;" );
ring rs2 = 0,(x,y,z),lp;

poly f0= x+y+z;
poly f1= x^2 + 6*x + 3*y + 6*z - 4;
poly f2= y^2 + 2*x -7*y + 5 + 2*z;
poly f3= x^2 + y^2 + z^2 - 1;
ideal i=f1,f2,f3;
ideal im=f0,f1,f2,f3;

def A=ures_solve(i);
setring A; SOL;
kill A;
setring rs2;
mp_res_mat(im);

def rinC = fglm_solve(i,30);
setring rinC;
rlist;

setring rs2;
kill rinC;
def rinC = triangL_solve(i,30);
setring rinC;
rlist;

setring rs2;
kill rinC;
def rinC = triangLf_solve(i,30);
setring rinC;
rlist;

setring rs2;
kill rinC;
def rinC = triangM_solve(i,30);
setring rinC;
rlist;
kill rinC;

////////////////////////////

tst_ignore( "ring rs6= 0,(x,y),lp;" );
ring rs6= 0,(x,y),lp;
poly f = 15x5 + x3 + x2 - 10;
laguerre_solve(f);

laguerre_solve(f,10);

tst_ignore( "ring rsc= (real,20,I),x,lp;" );
ring rsc= (real,20,I),x,lp;
poly f = (15.4+I*5)*x^5 + (25.0e-2+I*2)*x^3 + x2 - 10*I;
laguerre_solve(f);

////////////////////////////

ring rs7 = 0,(x),lp;
ideal v=16,0,11376,1046880,85949136;
interpolate( 3, v, 4 );

ring rs8 = 0,(x,y),dp;
ideal p = 2,3;
ideal v= 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16;
poly ip= interpolate( p,v,3 );
ip;

////////////////////////////

example simplexOut;

tst_status(1);$