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

-version 1.11, 2000/05/19 11:16:51
+version 1.12, 2000/05/24 15:24:54
 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.11 2000/05/19 11:16:51 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 SisComplex(a) {
-   local a,b,c,cc,v;
+   local n,i,j,k,b,p,q;
-   a = Sannfs3_laScala2("x^3-y^2*z^2");
+   n = Length(a);
-   b = a[0];
+   for (i=0; i<n-1; i++) {
-   v = [x,y,z];
+     if (Length(a[i+1]) != 0) {
-   c = Skernel(b[0],v);
+       b = a[i+1]*a[i];
-   c = c[0];
+       p = Length(b); q = Length(b[0]);
-   sm1_pmat([c,b[1],v]);
+       for (j=0; j<p; j++) {
-   Println("-----------------------------------");
+         for (k=0; k<q; k++) {
-   cc = sm1_res_div(c,b[1],v);
+           if (!IsZero(b[j,k])) {
-   sm1_pmat(sm1_gb(cc,v));
+              Print("Is is not complex at ");
-   c = Skernel(b[1],v);
+              Println([i,j,k]);
-   c = c[0];
+              return(false);
-   cc = sm1_res_div(c,b[2],v);
+           }
-   sm1_pmat(sm1_gb(cc,v));
+         }
-   return(a);
+       }
- }
+     }
- def test6() {
+   }
-   local a,b,c,cc,v;
+   return(true);
-   a = Sannfs3("x^3-y^2*z^2");
+ }
-   b = a[0];
-   v = [x,y,z];
-   c = Skernel(b[0],v);
-   c = c[0];
-   sm1_pmat([c,b[1],v]);
-   Println("-------ker = im for minimal ?---------------------");
-   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));
-   Println("------ ker=im for Schreyer ?------------------");
-   b = a[3];
-   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(a);
- }

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