Changeset 962b1f in git

Timestamp:

Feb 25, 2002, 7:14:54 PM (21 years ago)

Author:

Mathias Schulze <mschulze@…>

Branches:

(u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')

Children:

97403d7b5d41e537394fcbd7c261d7b70302488c

Parents:

8fe22ccc9de914e68da65f1ac8874c938ddd2809

Message:

*mschulze: liftbound replaced by new division


git-svn-id: file:///usr/local/Singular/svn/trunk@5925 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.69 2002-02-21 17:52:23 mschulze Exp $";
+version="$Id: gaussman.lib,v 1.70 2002-02-25 18:14:39 mschulze Exp $";
 category="Singularities";
 
@@ -42 +42 @@
 ";
 
-LIB "linalg.lib";
+LIB "linalg2.lib";
 
 ///////////////////////////////////////////////////////////////////////////////
@@ -312 +312 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-static proc liftbound(module U,module H,int k,int K)
-{
-  def R=basering;
-  ring T=0,s,ds;
-  module U=imap(R,U);
-  module H=imap(R,H);
-  degBound=k+K+1;
-  H=jet(lift(U,H),K);
-  degBound=0;
-  setring(R);
-  return(imap(T,H));
-}
-///////////////////////////////////////////////////////////////////////////////
-
 static proc saturate()
 {
@@ -356 +342 @@
 
   dbprint(printlevel-voice+2,"// transform H'' to saturation of H''");
-  H0=liftbound(H,freemodule(mu)*s^k,k,k);
+  H0=division(freemodule(mu)*s^k,H,k,intvec(1))[1];
 
   return(A0,r,H,H0,k);
@@ -372 +358 @@
 
   dbprint(printlevel-voice+2,"// transform A to saturation of H''");
-  matrix A=liftbound(H*s,A0*H+s^2*diff(matrix(H),s),k0+1,K);
+  matrix A=division(module(A0*H+s^2*diff(matrix(H),s)),H,K+1,intvec(1))[1]/s;
 
   return(A,A0,r);
@@ -381 +367 @@
 {
   dbprint(printlevel-voice+2,
-    "// compute eigenvalues e with multiplicities m of A");
+    "// compute eigenvalues e with multiplicities m of A0");
   matrix A;
   A,A0,r=basisrep(A0,r,H,k0,0);
@@ -1099 +1085 @@
     U=matrix(U)+s^k*V;
   }
+  attrib(U,"isSB",1);
 
   dbprint(printlevel-voice+2,"// transform to V-splitting basis");
   A=jet(A,0);
-  H=std(liftbound(U,H,0,k0+k1));
+  H=std(division(H,U,k0+k1,intvec(1))[1]);
 
   dbprint(printlevel-voice+2,"// compute V-leading terms of H''");
@@ -1210 +1197 @@
   degBound=0;
   H=simplify(jet(H,k0+k1),1);
+  attrib(H,"isSB",1);
 
   dbprint(printlevel-voice+2,"// compute matrix A0+sA1 of t");
-  A=liftbound(H,s*A*H+s^2*diff(matrix(H),s),k0+k1,1);
+  A=division(module(s*A*H+s^2*diff(matrix(H),s)),H,1,intvec(1))[1];
   A0=jet(A,0);
   A=jet(A,1)/s;