Changeset b87095 in git for Singular/LIB/gaussman.lib

Timestamp:

Feb 8, 2001, 3:20:29 PM (23 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'spielwiese', '6e5adcba05493683b94648c659a729c189812c77')

Children:

b6b6153e2b2b50dc5dc315b69fdb649c1d799fa9

Parents:

ab26fe772f0657938f6600e9c5b9cbdeeab7d32f

Message:

*mschulze: update for new jet and reduce commands


git-svn-id: file:///usr/local/Singular/svn/trunk@5212 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.32 2001-02-08 08:09:33 mschulze Exp $";
+version="$Id: gaussman.lib,v 1.33 2001-02-08 14:20:29 mschulze Exp $";
 category="Singularities";
 
@@ -145 +145 @@
 
     dbprint(printlevel-voice+2,"//gaussman::monodromy: compute C");
-    C=coeffs(system("rednf",sJ,w,U),m);
+    C=coeffs(reduce(w,sJ,U),m);
     A0=A0+C*var(1)^k;
 
@@ -177 +177 @@
         "//gaussman::monodromy: compute A on saturation");
       l=division(H*var(1),A0*H+var(1)^2*diff(matrix(H),var(1)));
-      A=system("series",N-1,module(l[1]),l[2]);
+      A=jet(l[1],N-1,l[2]);
       if(mide==0)
       {
@@ -397 +397 @@
 
     dbprint(printlevel-voice+2,"//gaussman::vfiltration: compute C");
-    C=coeffs(system("rednf",sJ,w,U),m);
+    C=coeffs(reduce(w,sJ,U),m);
     A=A+C*var(1)^k;
 
@@ -455 +455 @@
     "//gaussman::vfiltration: transform H0 to saturation");
   l=division(H,H0);
-  H0=system("series",N-1,module(l[1]),l[2]);
+  H0=jet(l[1],N-1,l[2]);
 
   dbprint(printlevel-voice+2,
@@ -490 +490 @@
     "//gaussman::vfiltration: compute A on saturation");
   l=division(H*var(1),A*H+var(1)^2*diff(matrix(H),var(1)));
-  A=system("series",N-1,module(l[1]),l[2]);
+  A=jet(l[1],N-1,l[2]);
 
   dbprint(printlevel-voice+2,"//gaussman::vfiltration: compute matrix M of A");
@@ -728 +728 @@
     "//gaussman::vfiltjacalg: compute multiplication in Jacobian algebra");
   list M;
+  matrix U=freemodule(ncols(m));
   for(i=ncols(m);i>=1;i--)
   {
-    M[i]=lift(V,coeffs(system("rednf",sJ,m[i]*m),m)*V);
+    M[i]=lift(V,coeffs(reduce(m[i]*m,sJ,U),m)*V);
   }
 
@@ -884 +885 @@
 }
 ///////////////////////////////////////////////////////////////////////////////
-
-proc tst_gaussm(poly f)
-{
-  echo=2;
-  basering;
-  f;
-  print(monodromy(f));
-  list l=vfiltration(f);
-  l;
-  vfiltjacalg(l);
-  gamma(l);
-}
-///////////////////////////////////////////////////////////////////////////////