Changeset a3432c in git for Singular/LIB/primdec.lib

Timestamp:

Sep 24, 1998, 2:09:26 PM (25 years ago)

Author:

Hans Schönemann <hannes@…>

Branches:

(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')

Children:

2386d98357886538fbc1e8f5e3e458fcb811ff1f

Parents:

56c52a7879fbc6de62cefba1da86b2edf2aadd4c

Message:

minor optimizations


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

File:

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

Singular/LIB/primdec.lib

@@ -1 @@
-// $Id: primdec.lib,v 1.26 1998-09-24 09:59:51 Singular Exp $
+// $Id: primdec.lib,v 1.27 1998-09-24 12:09:26 Singular Exp $
 ////////////////////////////////////////////////////////////////////////////////
 // primdec.lib                                                                //
@@ -11 +11 @@
 ////////////////////////////////////////////////////////////////////////////////
 
-version="$Id: primdec.lib,v 1.26 1998-09-24 09:59:51 Singular Exp $";
+version="$Id: primdec.lib,v 1.27 1998-09-24 12:09:26 Singular Exp $";
 info="
 LIBRARY: primdec.lib: PROCEDURE FOR PRIMARY DECOMPOSITION
@@ -383 +383 @@
    if(size(k1)==size(k2))
    {
-      for(j=size(k1);j>0;j--)
+      for(j=ncols(k1);j>0;j--)
       {
          if(leadexp(k1[j])!=leadexp(k2[j]))
@@ -401 +401 @@
    if(size(k1)==size(k2))
    {
-      for(j=size(k1);j>0;j--)
+      for(j=ncols(k1);j>0;j--)
       {
          if(leadexp(k1[j])!=leadexp(k2[j]))
@@ -431 +431 @@
 
   ideal i=maxideal(1);
-  int k=size(i);
+  int k=ncols(i);
   i[k]=0;
-  i=randomid(i,size(i),b);
+  i=randomid(i,ncols(i),b);
   ideal ires=maxideal(1);
   ires[k]=i[1]+var(k);
@@ -472 +472 @@
       else
       {
-         while (lead(i[m])/var(n-1)==0)
-        {
-            m=m+1;
-         }
-         m=m-1;
+        while (lead(i[m])/var(n-1)==0)
+        {
+          m=m+1;
+        }
+        m=m-1;
       }
       //check whether i[m] =(c*var(n)+h)^e modulo prm for some
@@ -482 +482 @@
       //if not (0) is returned, else var(n)+h is added to prm
 
-         e=deg(lead(i[m]));
-         t=leadcoef(i[m])*e*var(n)+(i[m]-lead(i[m]))/var(n)^(e-1);
-
-         i[m]=poly(e)^e*leadcoef(i[m])^(e-1)*i[m];
-
-         if (reduce(i[m]-t^e,prm,1) !=0)
-         {
-           return(ideal(0));
-         }
-         h=interred(t);
-         t=h[1];
+      e=deg(lead(i[m]));
+      t=leadcoef(i[m])*e*var(n)+(i[m]-lead(i[m]))/var(n)^(e-1);
+
+      i[m]=poly(e)^e*leadcoef(i[m])^(e-1)*i[m];
+
+      if (reduce(i[m]-t^e,prm,1) !=0)
+      {
+        return(ideal(0));
+      }
+      h=interred(t);
+      t=h[1];
 
       prm = prm,t;
@@ -530 +530 @@
          continue;
       }
-
 
       if(size(l[2*i])==0)