Changeset 9526639 in git
- Timestamp:
- Feb 12, 2002, 3:35:24 PM (21 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'f875bbaccd0831e36aaed09ff6adeb3eb45aeb94')
- Children:
- 6d07e39974898a3c343d0298c4eb25a179334ab7
- Parents:
- bd59178f61b0d020cc0cf4b16a448d78a242dfce
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/gaussman.lib
rbd5917 r9526639 1 1 /////////////////////////////////////////////////////////////////////////////// 2 version="$Id: gaussman.lib,v 1.6 5 2002-01-24 13:32:33mschulze Exp $";2 version="$Id: gaussman.lib,v 1.66 2002-02-12 14:35:24 mschulze Exp $"; 3 3 category="Singularities"; 4 4 … … 118 118 matrix gmsmatrix; matrix(gmsjacob)*gmsmatrix==matrix(gmsstd) 119 119 ideal gmsbasis; monomial vector space basis of Jacobian algebra 120 int gmsmax weight; maximal weight of variables120 int gmsmaxdeg; maximal weight of variables 121 121 @end format 122 122 KEYWORDS: singularities; Gauss-Manin connection; Brieskorn lattice … … 151 151 } 152 152 } 153 matrix a=l[2];153 matrix B=l[2]; 154 154 ideal m=kbase(g); 155 155 156 int gmsmax weight;156 int gmsmaxdeg; 157 157 for(i=nvars(R);i>=1;i--) 158 158 { 159 if(deg(var(i))>gmsmax weight)160 { 161 gmsmax weight=deg(var(i));159 if(deg(var(i))>gmsmaxdeg) 160 { 161 gmsmaxdeg=deg(var(i)); 162 162 } 163 163 } … … 183 183 184 184 execute("ring G="+string(charstr(R))+",("+s+","+varstr(R)+"),(ws("+ 185 string(deg(highcorner(g))+2*gmsmax weight)+"),"+os+",c);");185 string(deg(highcorner(g))+2*gmsmaxdeg)+"),"+os+",c);"); 186 186 187 187 poly gmspoly=imap(R,t); 188 188 ideal gmsjacob=imap(R,dt); 189 189 ideal gmsstd=imap(R,g); 190 matrix gmsmatrix=imap(R, a);190 matrix gmsmatrix=imap(R,B); 191 191 ideal gmsbasis=imap(R,m); 192 192 193 193 attrib(gmsstd,"isSB",1); 194 export gmspoly,gmsjacob,gmsstd,gmsmatrix,gmsbasis,gmsmax weight;194 export gmspoly,gmsjacob,gmsstd,gmsmatrix,gmsbasis,gmsmaxdeg; 195 195 196 196 return(G); … … 207 207 print(gmsmatrix); 208 208 print(gmsbasis); 209 gmsmax weight;209 gmsmaxdeg; 210 210 } 211 211 /////////////////////////////////////////////////////////////////////////////// … … 243 243 if(Kmax>=0) 244 244 { 245 p=jet(jet(p,K,v),(Kmax+1)*deg(var(1))-2*gmsmax weight);245 p=jet(jet(p,K,v),(Kmax+1)*deg(var(1))-2*gmsmaxdeg); 246 246 } 247 247 … … 272 272 p[k]=p[k]+ 273 273 jet(jet(diff(s*gmsmatrix[j,i],var(j+1))*var(1),Kmax,v), 274 (Kmax+1)*deg(var(1))-2*gmsmax weight);274 (Kmax+1)*deg(var(1))-2*gmsmaxdeg); 275 275 } 276 276 else … … 285 285 p[k]=p[k]-lead(p[k]); 286 286 } 287 while(deg(lead(p[k]))>(K+1)*deg(var(1))-2*gmsmax weight&&287 while(deg(lead(p[k]))>(K+1)*deg(var(1))-2*gmsmaxdeg&& 288 288 deg(lead(p[k]),v)<=K) 289 289 { … … 577 577 setring(R); 578 578 list l=imap(G,l); 579 kill G,gmsmax weight;579 kill G,gmsmaxdeg; 580 580 581 581 return(l); … … 909 909 setring(R); 910 910 list l=imap(G,l); 911 kill G,gmsmax weight;911 kill G,gmsmaxdeg; 912 912 attrib(l[5],"isSB",1); 913 913 … … 973 973 module H0=imap(G,H0); 974 974 ideal e=imap(G,e); 975 kill G,gmsmax weight;975 kill G,gmsmaxdeg; 976 976 977 977 dbprint(printlevel-voice+2,"// transform to Jordan basis");
Note: See TracChangeset
for help on using the changeset viewer.