Changeset 9526639 in git for Singular/LIB


Ignore:
Timestamp:
Feb 12, 2002, 3:35:24 PM (22 years ago)
Author:
Mathias Schulze <mschulze@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38dfc5131670d387a89455159ed1e071997eec94')
Children:
6d07e39974898a3c343d0298c4eb25a179334ab7
Parents:
bd59178f61b0d020cc0cf4b16a448d78a242dfce
Message:
*** empty log message ***


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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/gaussman.lib

    rbd5917 r9526639  
    11///////////////////////////////////////////////////////////////////////////////
    2 version="$Id: gaussman.lib,v 1.65 2002-01-24 13:32:33 mschulze Exp $";
     2version="$Id: gaussman.lib,v 1.66 2002-02-12 14:35:24 mschulze Exp $";
    33category="Singularities";
    44
     
    118118  matrix gmsmatrix;  matrix(gmsjacob)*gmsmatrix==matrix(gmsstd)
    119119  ideal gmsbasis;  monomial vector space basis of Jacobian algebra
    120   int gmsmaxweight;  maximal weight of variables
     120  int gmsmaxdeg;  maximal weight of variables
    121121@end format
    122122KEYWORDS: singularities; Gauss-Manin connection; Brieskorn lattice
     
    151151    }
    152152  } 
    153   matrix a=l[2];
     153  matrix B=l[2];
    154154  ideal m=kbase(g);
    155155
    156   int gmsmaxweight;
     156  int gmsmaxdeg;
    157157  for(i=nvars(R);i>=1;i--)
    158158  {
    159     if(deg(var(i))>gmsmaxweight)
    160     {
    161       gmsmaxweight=deg(var(i));
     159    if(deg(var(i))>gmsmaxdeg)
     160    {
     161      gmsmaxdeg=deg(var(i));
    162162    }
    163163  }
     
    183183
    184184  execute("ring G="+string(charstr(R))+",("+s+","+varstr(R)+"),(ws("+
    185     string(deg(highcorner(g))+2*gmsmaxweight)+"),"+os+",c);");
     185    string(deg(highcorner(g))+2*gmsmaxdeg)+"),"+os+",c);");
    186186
    187187  poly gmspoly=imap(R,t);
    188188  ideal gmsjacob=imap(R,dt);
    189189  ideal gmsstd=imap(R,g);
    190   matrix gmsmatrix=imap(R,a);
     190  matrix gmsmatrix=imap(R,B);
    191191  ideal gmsbasis=imap(R,m);
    192192
    193193  attrib(gmsstd,"isSB",1);
    194   export gmspoly,gmsjacob,gmsstd,gmsmatrix,gmsbasis,gmsmaxweight;
     194  export gmspoly,gmsjacob,gmsstd,gmsmatrix,gmsbasis,gmsmaxdeg;
    195195
    196196  return(G);
     
    207207  print(gmsmatrix);
    208208  print(gmsbasis);
    209   gmsmaxweight;
     209  gmsmaxdeg;
    210210}
    211211///////////////////////////////////////////////////////////////////////////////
     
    243243  if(Kmax>=0)
    244244  {
    245     p=jet(jet(p,K,v),(Kmax+1)*deg(var(1))-2*gmsmaxweight);
     245    p=jet(jet(p,K,v),(Kmax+1)*deg(var(1))-2*gmsmaxdeg);
    246246  }
    247247
     
    272272            p[k]=p[k]+
    273273              jet(jet(diff(s*gmsmatrix[j,i],var(j+1))*var(1),Kmax,v),
    274                 (Kmax+1)*deg(var(1))-2*gmsmaxweight);
     274                (Kmax+1)*deg(var(1))-2*gmsmaxdeg);
    275275          }
    276276          else
     
    285285        p[k]=p[k]-lead(p[k]);
    286286      }
    287       while(deg(lead(p[k]))>(K+1)*deg(var(1))-2*gmsmaxweight&&
     287      while(deg(lead(p[k]))>(K+1)*deg(var(1))-2*gmsmaxdeg&&
    288288        deg(lead(p[k]),v)<=K)
    289289      {
     
    577577  setring(R);
    578578  list l=imap(G,l);
    579   kill G,gmsmaxweight;
     579  kill G,gmsmaxdeg;
    580580
    581581  return(l);
     
    909909  setring(R);
    910910  list l=imap(G,l);
    911   kill G,gmsmaxweight;
     911  kill G,gmsmaxdeg;
    912912  attrib(l[5],"isSB",1);
    913913
     
    973973  module H0=imap(G,H0);
    974974  ideal e=imap(G,e);
    975   kill G,gmsmaxweight;
     975  kill G,gmsmaxdeg;
    976976
    977977  dbprint(printlevel-voice+2,"// transform to Jordan basis");
Note: See TracChangeset for help on using the changeset viewer.