Changeset d8d3af in git

Timestamp:

May 8, 2000, 5:51:55 PM (23 years ago)

Author:

Gerhard Pfister <pfister@…>

Branches:

(u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')

Children:

62e2d6009f561be67e002fef8503d15b30c4f8c4

Parents:

33c96aa0ca7b10d8e8b1b12a0b0bb5abaad3da9d

Message:

bug in equidim gefixt


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

File:

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

Singular/LIB/primdec.lib

@@ -1 @@
-// $Id: primdec.lib,v 1.57 2000-05-08 10:07:48 pfister Exp $
+// $Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $
 ///////////////////////////////////////////////////////////////////////////////
 // primdec.lib                                                               //
@@ -11 +11 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-version="$Id: primdec.lib,v 1.57 2000-05-08 10:07:48 pfister Exp $";
+version="$Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $";
 info="
 LIBRARY: primdec.lib   PROCEDURES FOR PRIMARY DECOMPOSITION
@@ -1764 +1764 @@
 }
 ///////////////////////////////////////////////////////////////////////////////
-proc equidim(ideal i)
+proc equidim(ideal i,list #)
 "USAGE:  equidim(i); i ideal           
  RETURN:  list = list of equidimensional ideals a1,...,as such that
@@ -1775 +1775 @@
   intvec w;
   int n;
+  ideal te=1;
   int a=attrib(i,"isSB");
   int homo=homog(i);
-   
+  if(size(#)!=0)
+  {
+    te=#[1];
+  }   
   if(((homo==1)||(a==1))&&(find(ordstr(basering),"l")==0)
                                 &&(find(ordstr(basering),"s")==0))
@@ -1800 +1804 @@
      ideal j=groebner(i);
   }
+  
+  ideal te=imap(P,te);
   list equ,equi,indep;
   if(homo==1)
@@ -1817 +1823 @@
 
   indep=maxIndependSet(j);
+
   string va=string(maxideal(1));
+
   execute "ring gnir1 = ("+charstr(basering)+"),("+indep[1][1]+"),("
                               +indep[1][2]+");";
-  execute "map phi=gnir,"+va+";";
   if(homo==1)
   {
-     ideal j=std(phi(i),hil,w);
+     ideal j=std(imap(gnir,i),hil,w);
   }
   else
   {
-     ideal j=groebner(phi(i));
+     ideal j=groebner(imap(gnir,i));
   }
   string quotring=prepareQuotientring(nvars(basering)-indep[1][3]);
+
   execute quotring;
+
   ideal j=imap(gnir1,j);
+
   kill gnir1;
   j=clearSB(j);
@@ -1843 +1853 @@
   ideal h=imap(quring,h);
   kill quring;
-
   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);
+  } 
+  else
+  {
+     te=intersect(te,equ[1]);
+  }
+ 
   j=std(j,l[2]);
 
-  equi=equidim(j);
+  equi=equidim(j,te);
   attrib(equi[1],"isSB",1);