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

Return to minimal.k CVS log

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

Diff for /OpenXM/src/k097/lib/minimal/minimal.k between version 1.11 and 1.13

-version 1.11, 2000/05/19 11:16:51
+version 1.13, 2000/06/08 08:37:53
 Line 1
 Line 1
 Line 1
- /* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.10 2000/05/07 02:10:44 takayama Exp $ */
+ /* $OpenXM: OpenXM/src/k097/lib/minimal/minimal.k,v 1.12 2000/05/24 15:24:54 takayama Exp $ */
  #define DEBUG 1
  /* #define ORDINARY 1 */
  /* If you run this program on openxm version 1.1.2 (FreeBSD),
 Line 1641  def sm1_gb(f,v) {
 Line 1641  def sm1_gb(f,v) {
 Line 1641  def sm1_gb(f,v) {
    sm1(" [f v] gb /FunctionValue set ");
  }
- def test5() {
+ def SisExact_h(c,v) {
-   local a,b,c,cc,v;
+   local a;
-   a = Sannfs3_laScala2("x^3-y^2*z^2");
+   v = ToString_array(v);
-   b = a[0];
+   a = [c,v];
-   v = [x,y,z];
+   sm1(a," isExact /FunctionValue set ");
-   c = Skernel(b[0],v);
-   c = c[0];
-   sm1_pmat([c,b[1],v]);
-   Println("-----------------------------------");
-   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(a);
  }
- def test6() {
-   local a,b,c,cc,v;
+ def SisComplex(a) {
-   a = Sannfs3("x^3-y^2*z^2");
+   local n,i,j,k,b,p,q;
-   b = a[0];
+   n = Length(a);
-   v = [x,y,z];
+   for (i=0; i<n-1; i++) {
-   c = Skernel(b[0],v);
+     if (Length(a[i+1]) != 0) {
-   c = c[0];
+       b = a[i+1]*a[i];
-   sm1_pmat([c,b[1],v]);
+       p = Length(b); q = Length(b[0]);
-   Println("-------ker = im for minimal ?---------------------");
+       for (j=0; j<p; j++) {
-   cc = sm1_res_div(c,b[1],v);
+         for (k=0; k<q; k++) {
-   sm1_pmat(sm1_gb(cc,v));
+           if (!IsZero(b[j,k])) {
-   c = Skernel(b[1],v);
+              Print("Is is not complex at ");
-   c = c[0];
+              Println([i,j,k]);
-   cc = sm1_res_div(c,b[2],v);
+              return(false);
-   sm1_pmat(sm1_gb(cc,v));
+           }
-   Println("------ ker=im for Schreyer ?------------------");
+         }
-   b = a[3];
+       }
-   c = Skernel(b[0],v);
+     }
-   c = c[0];
+   }
-   sm1_pmat([c,b[1],v]);
+   return(true);
-   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(a);
- }

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>