Changeset d8d3af in git


Ignore:
Timestamp:
May 8, 2000, 5:51:55 PM (24 years ago)
Author:
Gerhard Pfister <pfister@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '3720ae8bfcff4a4649ee98a15552089151d2d59b')
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:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    r33c96a rd8d3af  
    1 // $Id: primdec.lib,v 1.57 2000-05-08 10:07:48 pfister Exp $
     1// $Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $
    22///////////////////////////////////////////////////////////////////////////////
    33// primdec.lib                                                               //
     
    1111///////////////////////////////////////////////////////////////////////////////
    1212
    13 version="$Id: primdec.lib,v 1.57 2000-05-08 10:07:48 pfister Exp $";
     13version="$Id: primdec.lib,v 1.58 2000-05-08 15:51:55 pfister Exp $";
    1414info="
    1515LIBRARY: primdec.lib   PROCEDURES FOR PRIMARY DECOMPOSITION
     
    17641764}
    17651765///////////////////////////////////////////////////////////////////////////////
    1766 proc equidim(ideal i)
     1766proc equidim(ideal i,list #)
    17671767"USAGE:  equidim(i); i ideal           
    17681768 RETURN:  list = list of equidimensional ideals a1,...,as such that
     
    17751775  intvec w;
    17761776  int n;
     1777  ideal te=1;
    17771778  int a=attrib(i,"isSB");
    17781779  int homo=homog(i);
    1779    
     1780  if(size(#)!=0)
     1781  {
     1782    te=#[1];
     1783  }   
    17801784  if(((homo==1)||(a==1))&&(find(ordstr(basering),"l")==0)
    17811785                                &&(find(ordstr(basering),"s")==0))
     
    18001804     ideal j=groebner(i);
    18011805  }
     1806 
     1807  ideal te=imap(P,te);
    18021808  list equ,equi,indep;
    18031809  if(homo==1)
     
    18171823
    18181824  indep=maxIndependSet(j);
     1825
    18191826  string va=string(maxideal(1));
     1827
    18201828  execute "ring gnir1 = ("+charstr(basering)+"),("+indep[1][1]+"),("
    18211829                              +indep[1][2]+");";
    1822   execute "map phi=gnir,"+va+";";
    18231830  if(homo==1)
    18241831  {
    1825      ideal j=std(phi(i),hil,w);
     1832     ideal j=std(imap(gnir,i),hil,w);
    18261833  }
    18271834  else
    18281835  {
    1829      ideal j=groebner(phi(i));
     1836     ideal j=groebner(imap(gnir,i));
    18301837  }
    18311838  string quotring=prepareQuotientring(nvars(basering)-indep[1][3]);
     1839
    18321840  execute quotring;
     1841
    18331842  ideal j=imap(gnir1,j);
     1843
    18341844  kill gnir1;
    18351845  j=clearSB(j);
     
    18431853  ideal h=imap(quring,h);
    18441854  kill quring;
    1845 
    18461855  list l=minSat(j,h);
     1856
    18471857  equ[1]=l[1];
    18481858  attrib(equ[1],"isSB",1);
    1849  
     1859
     1860  if(size(reduce(te,equ[1]))==0)
     1861  {
     1862    equ[1]=ideal(1);
     1863    attrib(equ[1],"isSB",1);
     1864  }
     1865  else
     1866  {
     1867     te=intersect(te,equ[1]);
     1868  }
     1869 
    18501870  j=std(j,l[2]);
    18511871
    1852   equi=equidim(j);
     1872  equi=equidim(j,te);
    18531873  attrib(equi[1],"isSB",1);
    18541874
Note: See TracChangeset for help on using the changeset viewer.