=================================================================== RCS file: /home/cvs/OpenXM/src/k097/lib/minimal/minimal-test.k,v retrieving revision 1.7 retrieving revision 1.8 diff -u -p -r1.7 -r1.8 --- OpenXM/src/k097/lib/minimal/minimal-test.k 2000/07/30 02:26:25 1.7 +++ OpenXM/src/k097/lib/minimal/minimal-test.k 2000/07/31 01:21:41 1.8 @@ -1,4 +1,4 @@ -/* $OpenXM: OpenXM/src/k097/lib/minimal/minimal-test.k,v 1.6 2000/07/26 02:21:31 takayama Exp $ */ +/* $OpenXM: OpenXM/src/k097/lib/minimal/minimal-test.k,v 1.7 2000/07/30 02:26:25 takayama Exp $ */ load["minimal.k"]; def sm1_resol1(p) { sm1(" p resol1 /FunctionValue set "); @@ -41,96 +41,6 @@ def test8() { SisComplex(a): */ -def test8a() { - local p,pp,ans,b,c,cc,ww, ans_all; - f = "x^3-y^2*z^2"; - p = Sannfs(f,"x,y,z"); - sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set "); - ww = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - /* Removed "x",1, ... ===> It causes an error. I do not know the reason.*/ - Sweyl("x,y,z",ww); - pp = Map(p,"Spoly"); - /* return(pp); */ - /* pp = - [y*Dy-z*Dz , -2*x*Dx-3*y*Dy+1 , 2*x*Dy*Dz^2-3*y*Dx^2 , - 2*x*Dy^2*Dz-3*z*Dx^2 , 2*x*z*Dz^3-3*y^2*Dx^2+4*x*Dz^2 ] - */ - ans_all = Sschreyer(pp); - ans = ans_all[0]; - /* ans = sm1_resol1([pp,"x,y,z",ww]); */ - /* Schreyer is in ans. */ - - v = [x,y,z]; - b = ans; - Println("------ ker=im for Schreyer ?----- wrong method!!!-----------"); - c = Skernel(b[0],v); - c = c[0]; - sm1_pmat([c,b[1],v]); - cc = sm1_res_div(c,b[1],v); - sm1_pmat(sm1_gb(cc,v)); - c = Skernel(b[1],v); - c = c[0]; - cc = sm1_res_div(c,b[2],v); - sm1_pmat(sm1_gb(cc,v)); - return(ans); -} - -/* Comparing two constructions */ -def test9() { - local p,pp,ans,b,c,cc,ww,ww2,ans_all,ans2; - f = "x^3-y^2*z^2"; - p = Sannfs(f,"x,y,z"); - ww2 = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set "); - Sweyl("x,y,z",ww2); - pp = Map(p,"Spoly"); - ans = sm1_resol1([pp,"x,y,z",ww2]); - - f = "x^3-y^2*z^2"; - p = Sannfs(f,"x,y,z"); - sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set "); - ww = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - Sweyl("x,y,z",ww); - pp = Map(p,"Spoly"); - ans_all = Sschreyer(pp); - ans2 = ans_all[0]; - - return([ans,ans2]); - -} - -/* Check if the complex by Sschreyer() is exact or not in our example? */ -def test10() { - local p,pp,ans,b,c,cc,ww,ww2,ans_all,ans2, r; - f = "x^3-y^2*z^2"; - p = Sannfs(f,"x,y,z"); - ww2 = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set "); - Sweyl("x,y,z",ww2); - pp = Map(p,"Spoly"); - ans = sm1_resol1([pp,"x,y,z",ww2]); - - f = "x^3-y^2*z^2"; - p = Sannfs(f,"x,y,z"); - sm1(" p 0 get { [(x) (y) (z) (Dx) (Dy) (Dz)] laplace0 } map /p set "); - ww = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - Sweyl("x,y,z",ww); - pp = Map(p,"Spoly"); - ans_all = Sschreyer(pp); /* Schreyer by LaScala-Stillman */ - ans2 = ans_all[0]; - - sm1(" /gb.verbose 1 def "); - - ww2 = [["x",-1,"y",-1,"z",-1,"Dx",1,"Dy",1,"Dz",1]]; - Sweyl("x,y,z",ww2); - ans2 = ReParse(ans2); - r= IsExact_h(ans2,[x,y,z]); - Print(r); - - return([r,[ans,ans2]]); - -} - def test11() { local a; a = test_ann3("x^3-y^2*z^2"); @@ -263,7 +173,7 @@ def test14a() { ans2 = [x1*Dx1+2*x2*Dx2+3*x3*Dx3 , Dx1^2-Dx2*h , -Dx1*Dx2+Dx3*h , Dx2^2-Dx1*Dx3 ]; ans2 = ReParse(ans2); - return(Sminimal(ans2,"homogenized")); + return(Sminimal(ans2,["homogenized"])); } def test15() { @@ -289,7 +199,7 @@ def test15b() { Sweyl("x1,x2,x3",ww2); ans2 = [Dx1^2-Dx2*h , -Dx1*Dx2+Dx3*h , Dx2^2-Dx1*Dx3 ]; ans2 = ReParse(ans2); - return(Sminimal(ans2,"homogenized")); + return(Sminimal(ans2,["homogenized"])); } def test15c() {