OpenXM/src/k097/lib/minimal/minimal-test.k - diff

Return to minimal-test.k CVS log

Up to [local] / OpenXM / src / k097 / lib / minimal

-version 1.1, 2000/05/24 15:31:28
+version 1.2, 2000/06/08 08:37:53
 Line 1
 Line 1
 Line 1
- /* $OpenXM$ */
+ /* $OpenXM: OpenXM/src/k097/lib/minimal/minimal-test.k,v 1.1 2000/05/24 15:31:28 takayama Exp $ */
  load["minimal.k"];
  def test5() {
    local a,b,c,cc,v;
 Line 173  def test9() {
 Line 173  def test9() {
 Line 173  def test9() {
    ans2 = ans_all[0];
    return([ans,ans2]);
+ }
+ /* Check if the complex is exact or not? */
+ 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];
+   r= SisExact_h(ans2,[x,y,z]);
+   Print(r);
+   return([r,[ans,ans2]]);
  }