Changeset 51534b6 in git

Timestamp:

Jan 22, 2002, 12:18:28 PM (22 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'spielwiese', 'e7cc1ebecb61be8b9ca6c18016352af89940b21a')

Children:

ddb18ba0330724c7ea767a4c646a66e5b1e8399f

Parents:

4d97ce4808e8a46aa628bb52456558bf2b489d33

Message:

*** empty log message ***


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

File:

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

Singular/LIB/gaussman.lib

@@ -1 +1 @@
 ///////////////////////////////////////////////////////////////////////////////
-version="$Id: gaussman.lib,v 1.60 2002-01-16 13:16:17 mschulze Exp $";
+version="$Id: gaussman.lib,v 1.61 2002-01-22 11:18:28 mschulze Exp $";
 category="Singularities";
 
@@ -385 +385 @@
 {
   def R=basering;
-  ring S=0,s,ds;
+  ring T=0,s,ds;
   module U=imap(R,U);
   module H=imap(R,H);
@@ -392 +392 @@
   degBound=0;
   setring(R);
-  return(imap(S,H));
+  return(imap(T,H));
 }
 ///////////////////////////////////////////////////////////////////////////////
@@ -574 +574 @@
   A,A0,r,H0,U0,e,m=transf(A,A0,r,H,H0,e,m,k0,k1,0,0);
 
+  list l=jordan(A,e,m);
   setring(R);
-  return(jordan(imap(G,A),imap(G,e),m));
+  list l=imap(G,l);
+  kill G,gmsmaxweight;
+
+  return(l);
 }
 example
@@ -901 +905 @@
   module v=simplify(jet(H*U0*H0,2*k0)/s^(2*k0),1);
 
+  kill l;
+  list l=sppnf(a,w,v)+list(gmsbasis,gmsstd);
   setring(R);
-  ideal g=imap(G,gmsstd);
-  attrib(g,"isSB",1);
-  return(sppnf(imap(G,a),w,imap(G,v))+list(imap(G,gmsbasis),g));
+  list l=imap(G,l);
+  kill G,gmsmaxweight;
+  attrib(l[5],"isSB",1);
+
+  return(l);
 }
 example
@@ -961 +969 @@
   A,A0,r,H0,U0,e,m=transf(A,A0,r,H,H0,e,m,k0,k1,k0+k1,k0+k1);
 
+  ring S=0,s,(ds,c);
+  matrix A=imap(G,A);
+  module H0=imap(G,H0);
+  ideal e=imap(G,e);
+  kill G,gmsmaxweight;
+
   dbprint(printlevel-voice+2,"// transform to Jordan basis");
   module U=jordanbasis(A,e,m)[1];
   matrix V=inverse(U);
   A=V*A*U;
-  H=V*H0;
+  module H=V*H0;
 
   dbprint(printlevel-voice+2,"// compute splitting of V-filtration");
@@ -972 +986 @@
   V=matrix(0,mu,mu);
   matrix v[mu^2][1];
-  A0=commutator(jet(A,0));
+  matrix A0=commutator(jet(A,0));
   for(k=1;k<=k0+k1;k++)
   {
@@ -995 +1009 @@
   {
     i0=leadexp(H1[k])[nvars(basering)+1];
-    j0=ord(H1[k])/deg(s);
+    j0=ord(H1[k]);//deg(s);
     H0[k]=lead(H1[k]);
     H1[k]=H1[k]-lead(H1[k]);
@@ -1001 +1015 @@
     {
       i=leadexp(H1[k])[nvars(basering)+1];
-      j=ord(H1[k])/deg(s);
+      j=ord(H1[k]);//deg(s);
       while(A[i,i]+j==A[i0,i0]+j0)
       {
@@ -1007 +1021 @@
         H1[k]=H1[k]-lead(H1[k]);
         i=leadexp(H1[k])[nvars(basering)+1];
-        j=ord(H1[k])/deg(s);
+        j=ord(H1[k]);//deg(s);
       }
     }
@@ -1059 +1073 @@
     I=list();
     I0=module();
-    U0=std(0);
+    module U0=std(0);
     for(i=size(F);i>=1;i--)
     {
@@ -1089 +1103 @@
 
   dbprint(printlevel-voice+2,"// compute reduced standard basis of H''");
-  ring S=0,s,ds;
-  module H=imap(G,H);
-  degBound=k0+k1+1;
+  degBound=k0+k1+2;
   option("redSB");
-
   H=std(H);
-
+  option("noredSB");
   degBound=0;
   H=simplify(jet(H,k0+k1),1);
-  setring(G);
-  H=imap(S,H);
 
   dbprint(printlevel-voice+2,"// compute matrix A0+sA1 of t");
@@ -1107 +1116 @@
 
   setring(R);
-  return(list(imap(G,A0),imap(G,A)));
+  matrix A0=imap(S,A0);
+  matrix A1=imap(S,A);
+  kill S;
+  return(list(A0,A1));
 }
 example