Changeset 3c23b65 in git for Singular

Timestamp:

Dec 16, 1999, 6:06:23 PM (24 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')

Children:

56025355795b496dc24087a340263f72e5c24ad5

Parents:

7fe30d0984ffc5114d6a954e3abc953620514049

Message:

mschulze: modified invunit


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

File:

• Singular/LIB/mondromy.lib (modified) (4 diffs)

Singular/LIB/mondromy.lib

@@ -1 +1 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-version="$Id: mondromy.lib,v 1.12 1999-12-15 22:29:51 greuel Exp $";
+version="$Id: mondromy.lib,v 1.13 1999-12-16 17:06:23 mschulze Exp $";
 info="
 LIBRARY:  mondromy.lib  PROCEDURES TO COMPUTE THE MONODROMY OF A SINGULARITY
@@ -24 +24 @@
 LIB "sing.lib";
 LIB "jordan.lib";
-LIB "poly.lib";
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -163 +162 @@
   if(pcvmindeg(u)==0)
   {
-    poly c = jet(u,0);
-    poly h = 1-u/c;
-    int d = mindeg(h);
-    int m = n div d;
-    if( m*d < n )
-    {
-      m=m+1;
-    }
-    poly g = 1;
-    poly v = 1;
     dbprint(printlevel-voice+2,"//computing inverse...");
     int t=timer;
-    int ii;
-    for( ii=1; ii<=m; ii++ )
-    {
-      g = jet(g*h,n);
-      v = v+g;
-    }
-    v = jet(v,n)/c;    
+    poly u0=jet(u,0);
+    u=jet(1-u/u0,n);
+    poly ui=u;
+    poly v=1+u;
+    int m=n div pcvmindeg(u);
+    int i;
+    for(i=1;i<m;i++)
+    {
+      ui=jet(ui*u,n);
+      v=v+ui;
+    }
+    v=jet(v,n)/u0;    
     dbprint(printlevel-voice+2,"//...inverse computed ["+string(timer-t)+
       " secs, "+string((memory(1)+1023)/1024)+" K]");
-
     return(v);
   }
@@ -196 +189 @@
 { "EXAMPLE:"; echo=2;
   ring R=0,(x,y),dp;
-  invunit(3+x+x2+y3,10);
+  invunit(2+x3+xy4,10);
 }