Changeset 80b3cd in git


Ignore:
Timestamp:
Aug 3, 1999, 3:01:02 PM (25 years ago)
Author:
Olaf Bachmann <obachman@…>
Branches:
(u'spielwiese', '5b153614cbc72bfa198d75b1e9e33dab2645d9fe')
Children:
46e39d5cb05fa4ee94268859b7b06dfdb5e76678
Parents:
0914249ec3414c61666a68fde6a8c2bb3babd2d2
Message:
* merged changes of gmg manually


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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/primdec.lib

    r091424 r80b3cd  
    1 // $Id: primdec.lib,v 1.43 1999-08-03 11:43:09 obachman Exp $
     1// $Id: primdec.lib,v 1.44 1999-08-03 13:01:02 obachman Exp $
    22///////////////////////////////////////////////////////////////////////////////
    33// primdec.lib                                                               //
     
    1111///////////////////////////////////////////////////////////////////////////////
    1212
    13 version="$Id: primdec.lib,v 1.43 1999-08-03 11:43:09 obachman Exp $";
     13version="$Id: primdec.lib,v 1.44 1999-08-03 13:01:02 obachman Exp $";
    1414info="
    1515LIBRARY: primdec.lib:   PROCEDURES FOR PRIMARY DECOMPOSITION
     
    1717          Wolfram Decker, email: decker@math.uni-sb.de         (SY)
    1818          Hans Schoenemann, email: hannes@mathematik.uni-kl.de (SY)
    19           last change: 99/07/23
    2019
    2120PROCEDURES:
    22   primdecGTZ(I);   complete primary decomposition via Gianni,Trager,Zacharias
    23   primdecSY(I);    complete primary decomposition via Shimoyama-Yokoyama
    24   minAssGTZ(I);    the minimal associated primes via Gianni,Trager,Zacharias
    25   minAssChar(I);  the minimal associated primes using characteristic sets
    26   testPrimary(..); tests the result of the primary decomposition
    27   radical(I);      computes the radical of the ideal I
    28   equiRadical(I);  the radical of the equidimensional part of the ideal I
    29   prepareAss(I);   list of radicals of the equidimensional components of I
     21  primdecGTZ(I);    complete primary decomposition via Gianni,Trager,Zacharias
     22  primdecSY(I...);  complete primary decomposition via Shimoyama-Yokoyama
     23  minAssGTZ(I);     the minimal associated primes via Gianni,Trager,Zacharias
     24  minAssChar(I...); the minimal associated primes using characteristic sets
     25  testPrimary(L,k); tests the result of the primary decomposition
     26  radical(I);       computes the radical of the ideal I
     27  equiRadical(I);   the radical of the equidimensional part of the ideal I
     28  prepareAss(I);    list of radicals of the equidimensional components of I
    3029
    3130REMARK:
    32  These procedures are implemented to be used in characteristic 0.
    33  They work also in positive characteristic >> 0.
    34  In small characteristic primdecGTZ, minAssGTZ, radical and equiRadical may not
    35  terminate and primdecSY and minAssChar may not give a complete decomposition.
     31These procedures are implemented to be used in characteristic 0.
     32They work also in positive characteristic >> 0.
     33In small characteristic primdecGTZ, minAssGTZ, radical and equiRadical may not
     34terminate and primdecSY and minAssChar may not give a complete decomposition.
    3635";
    3736
     
    552551               for(k=2;k<=r;k++)
    553552               {
    554               keepprime[size(l)/2+k-1]=interred(keepprime[i]+ideal(act[1][k]));
     553                  keepprime[size(l)/2+k-1]=interred(keepprime[i]+ideal(act[1][k]));
    555554               }
    556555               keepprime[i]=interred(keepprime[i]+ideal(act[1][1]));
     
    875874                if(lead(primary[2*@k-1][@n])/var(zz)!=0)
    876875                {
    877                 jmap1[zz]=-1/leadcoef(primary[2*@k-1][@n])*primary[2*@k-1][@n]
     876                   jmap1[zz]=-1/leadcoef(primary[2*@k-1][@n])*primary[2*@k-1][@n]
    878877                   +2/leadcoef(primary[2*@k-1][@n])*lead(primary[2*@k-1][@n]);
    879878                    jmap2[zz]=primary[2*@k-1][@n];
     
    15821581   int j,k,odim,ndim,count;
    15831582   attrib(@pr[1],"isSB",1);
    1584  
     1583   if(#[1]==77)
     1584   {
     1585     odim=dim(@pr[1]);
     1586     count=1;
     1587     intvec pos;
     1588     pos[size(@pr)]=0;
     1589     for(j=2;j<=size(@pr);j++)
     1590     {
     1591        attrib(@pr[j],"isSB",1);
     1592        ndim=dim(@pr[j]);
     1593        if(ndim>odim)
     1594        {
     1595           for(k=count;k<=j-1;k++)
     1596           {
     1597              pos[k]=1;
     1598           }
     1599           count=j;
     1600           odim=ndim;
     1601        }
     1602        if(ndim<odim)
     1603        {
     1604           pos[j]=1;
     1605        }
     1606     }
     1607     for(j=1;j<=size(@pr);j++)
     1608     {
     1609        if(pos[j]!=1)
     1610        {
     1611            @res[j]=decomp(@pr[j],2);
     1612        }
     1613        else
     1614        {
     1615           @res[j]=empty;
     1616        }
     1617     }
     1618   }
     1619   else
     1620   {
    15851621     ser=ideal(1);
    15861622     for(j=1;j<=size(@pr);j++)
     
    15951631//       }
    15961632     }
     1633   }
    15971634
    15981635   @res=union(@res);
     
    39934030///////////////////////////////////////////////////////////////////////////////
    39944031
    3995 proc primdecSY(ideal i)
     4032proc primdecSY(ideal i, list #))
    39964033"USAGE:   primdecSY(i); i ideal, c int
    39974034         if c=0, the given ordering of the variables is used.
     
    40484085         Otherwise, the system tries to find an optimal ordering,
    40494086         which in some cases may considerably speed up the algorithm
     4087RETURN:  list = the minimal associated prime ideals of i
    40504088NOTE:    implemented for characteristic 0, works also in char k >> 0,
    40514089         the result may be not compltely decomposed in small characteristic
     
    40714109"USAGE:   equiRadical(i); i ideal
    40724110RETURN:  ideal, intersection of associated primes of i of maximal  dimension
    4073 NOTE:    designed for characteristic 0, works also in char k > 0 if it termi-
    4074          nates, may result in an infinite loop in small characteristic
     4111NOTE:    designed for characteristic 0, works also in char k > 0 if it
     4112         terminates, may result in an infinite loop in small characteristic
    40754113EXAMPLE: example equiRadical; shows an example
    40764114"
Note: See TracChangeset for help on using the changeset viewer.