Changeset a23a7fd in git

Timestamp:

Dec 6, 2000, 3:56:26 PM (23 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

ac31058626d6128287fd40e307e34fc22c7c0beb

Parents:

a79bc624a96888a1bce8c52f3c8f6a33640de235

Message:

Singular/LIB/gaussman.lib


git-svn-id: file:///usr/local/Singular/svn/trunk@4824 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• Singular/LIB/gaussman.lib (modified) (8 diffs)

Singular/LIB/gaussman.lib

@@ -1 +1 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-version="$Id: gaussman.lib,v 1.13 2000-12-06 12:34:50 mschulze Exp $";
+version="$Id: gaussman.lib,v 1.14 2000-12-06 14:56:26 mschulze Exp $";
 info="
 LIBRARY:  gaussman.lib  GAUSS-MANIN CONNECTION OF A SINGULARITY
@@ -820 +820 @@
   }
 
+  dbprint(printlevel-voice+2,
+    "//gaussman::vjacob: compute multiplication of Jacobian algebra");
   list M;
   for(i=ncols(m);i>=1;i--)
@@ -836 +838 @@
   while(r0<number(s[ncols(s)]-s[1]))
   {
+    dbprint(printlevel-voice+2,
+      "//gaussman::vjacob: compute next difference r0 of spectral numbers");
     r1=number(s[ncols(s)]-s[1]);
     for(j=ncols(s);j>=1;j--)
@@ -856 +860 @@
     }
     r0=r1;
+    dbprint(printlevel-voice+2,"//gaussman::vjacob: r0="+string(r1));
 
     l=ideal();
@@ -863 +868 @@
     }
 
-    j0=mu-d[ncols(s)]+1;
-    for(j=ncols(s);j>=1;j--)
+    dbprint(printlevel-voice+2,"//gaussman::vjacob: compute conditions");
+    j=ncols(s);
+    j0=mu;
+    while(j>=1)
     {
       i0=1;
@@ -883 +890 @@
         if(i0>1)
         {
-          l[k]=l[k],M0[1..i0-1,j0..j0+d[j]-1];
+          l[k]=l[k],M0[1..i0-1,j0-d[j]+1..j0];
         }
       }
       j0=j0-d[j];
-    }
-
+      j--;
+    }
+
+    dbprint(printlevel-voice+2,"//gaussman::vjacob: compute condition matrix");
     L=transpose(module(l[1]));
     for(k=2;k<=ncols(m);k++)
@@ -895 +904 @@
     }
 
+    dbprint(printlevel-voice+2,
+      "//gaussman::vjacob: compute kernel of condition matrix");
     v0=v0+list(syz(L));
     s0=s0,r0;
   }
 
+  dbprint(printlevel-voice+2,"//gaussman::vjacob: compute graded parts");
   option(redSB);
   for(i=1;i<size(v0);i++)
@@ -906 +918 @@
   }
 
+  dbprint(printlevel-voice+2,
+    "//gaussman::vjacob: remove trivial graded parts");
   i=1;
   while(size(v0[i])==0)