Changeset c2868b in git

Timestamp:

May 9, 2000, 3:26:10 PM (23 years ago)

Author:

Gerhard Pfister <pfister@…>

Branches:

(u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')

Children:

b1926022497677c8d2f91e7c017f1c7a2317bc9a

Parents:

8cd4dfcd0c10c84e83c8f1397608fde5aad6a384

Message:

bug in equidim entfernt


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

File:

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

Singular/LIB/primdec.lib

@@ -1 @@
-// $Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $
+// $Id: primdec.lib,v 1.59 2000-05-09 13:26:10 pfister Exp $
 ///////////////////////////////////////////////////////////////////////////////
 // primdec.lib                                                               //
@@ -11 +11 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-version="$Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $";
+version="$Id: primdec.lib,v 1.59 2000-05-09 13:26:10 pfister Exp $";
 info="
 LIBRARY: primdec.lib   PROCEDURES FOR PRIMARY DECOMPOSITION
@@ -1807 +1807 @@
   ideal te=imap(P,te);
   list equ,equi,indep;
+  int b;
   if(homo==1)
   {
@@ -1817 +1818 @@
   if ((dim(j)==-1)||(size(j)==0)||(nvars(basering)==1)||(dim(j)==0))
   {
+    if((size(reduce(te,j))==0)&&(deg(te[1])>0))
+    {
+      setring P;
+      return(eq);
+    }
     setring P;
+    i=imap(gnir,j);
     eq=i;
     return(eq);
@@ -1855 +1862 @@
   list l=minSat(j,h);
 
-  equ[1]=l[1];
-  attrib(equ[1],"isSB",1);
-
-  if(size(reduce(te,equ[1]))==0)
-  {
-    equ[1]=ideal(1);
-    attrib(equ[1],"isSB",1);
-  } 
+  j=std(j,l[2]);
+  attrib(l[1],"isSB",1);
+  if(size(reduce(te,l[1]))!=0)
+  {
+    te=intersect(te,l[1]);
+    equi=equidim(j,te);  
+    if(dim(l[1])==dim(j))
+    {
+      equi[size(equi)]=intersect(l[1],equi[size(equi)]);
+    }
+    else
+    {
+      equi[size(equi)+1]=l[1];
+    }
+  }
   else
   {
-     te=intersect(te,equ[1]);
-  }
- 
-  j=std(j,l[2]);
-
-  equi=equidim(j,te);
-  attrib(equi[1],"isSB",1);
-
-  if(dim(equ[1])==dim(equi[1]))
-  {
-    equi[1]=intersect(equ[1],equi[1]);
-    equ=equi;
-  }
-  else
-  {
-     for(n=1;n<=size(equi);n++)
-     {
-       if(deg(equi[n][1])>0)
-       {
-         equ[size(equ)+1]=equi[n];
-       }
-    }
-  }
+    equi=equidim(j,te);
+  }
+  b=size(equi);
+
   setring P;
-  eq=imap(gnir,equ);
+  if(b!=0)
+  {
+    eq=imap(gnir,equi);
+  }
   kill gnir;
   return(eq);