Changeset 712167 in git
- Timestamp:
- Jul 13, 2007, 2:08:34 PM (16 years ago)
- Branches:
- (u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
- Children:
- 043cba9d4d180ed9cd19fd2212f6b7e506874565
- Parents:
- be4db46303a610786bdd624ae324eb267ae77931
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/hnoether.lib
rbe4db4 r712167 1 version="$Id: hnoether.lib,v 1.5 3 2006-08-02 15:40:51Singular Exp $";1 version="$Id: hnoether.lib,v 1.54 2007-07-13 12:08:34 Singular Exp $"; 2 2 category="Singularities"; 3 3 info=" … … 892 892 { "EXAMPLE:"; echo = 2; 893 893 ring exring = 7,(x,y),ds; 894 list hne=develop(4x98+2x49y7+x11y14+2y14);895 print( hne[1]);894 list Hne=develop(4x98+2x49y7+x11y14+2y14); 895 print(Hne[1]); 896 896 // therefore the HNE is: 897 897 // z(-1)= 3*z(0)^7 + z(0)^7*z(1), … … 902 902 // (the missing x in the last line indicates that it is not complete.) 903 903 hne[2]; 904 param( hne);904 param(Hne); 905 905 // parametrization: x(t)= -t^14+O(t^21), y(t)= -3t^98+O(t^105) 906 906 // (the term -t^109 in y may have a wrong coefficient) 907 displayHNE( hne);907 displayHNE(Hne); 908 908 } 909 909 … … 1059 1059 ring exring=0,(x,y,t),ds; 1060 1060 poly f=x3+2xy2+y2; 1061 list hne=develop(f);1062 list hne_extended=extdevelop( hne,10);1061 list Hne=develop(f); 1062 list hne_extended=extdevelop(Hne,10); 1063 1063 // compare the HNE matrices ... 1064 print( hne[1]);1065 print( hne_extended[1]);1064 print(Hne[1]); 1065 print(Hne_extended[1]); 1066 1066 // ... and the resulting parametrizations: 1067 param( hne);1067 param(Hne); 1068 1068 param(hne_extended); 1069 1069 param(hne_extended,0); … … 1257 1257 { "EXAMPLE:"; echo = 2; 1258 1258 ring exring=0,(x,y),dp; 1259 list hne=develop(y4+2x3y2+x6+x5y);1260 list INV=invariants( hne);1259 list Hne=develop(y4+2x3y2+x6+x5y); 1260 list INV=invariants(Hne); 1261 1261 INV[1]; // the characteristic exponents 1262 1262 INV[2]; // the generators of the semigroup of values … … 1265 1265 INV[6]; // the sequence of multiplicities 1266 1266 // To display the invariants more 'nicely': 1267 displayInvariants( hne);1267 displayInvariants(Hne); 1268 1268 ///////////////////////////// 1269 1269 INV=invariants((x2-y3)*(x3-y5)); … … 1422 1422 { "EXAMPLE:"; echo = 2; 1423 1423 ring exring=0,(x,y),dp; 1424 list hne=develop(y4+2x3y2+x6+x5y);1425 displayInvariants( hne);1424 list Hne=develop(y4+2x3y2+x6+x5y); 1425 displayInvariants(Hne); 1426 1426 } 1427 1427 /////////////////////////////////////////////////////////////////////////////// … … 1605 1605 "EXAMPLE:"; echo = 2; 1606 1606 ring r=0,(x,y),dp; 1607 list hne=hnexpansion((x2-y3)*(x2+y3));1608 intersection( hne[1],hne[2]);1607 list Hne=hnexpansion((x2-y3)*(x2+y3)); 1608 intersection(Hne[1],Hne[2]); 1609 1609 } 1610 1610 /////////////////////////////////////////////////////////////////////////////// … … 1761 1761 "EXAMPLE:"; echo = 2; 1762 1762 ring r=0,(x,y),dp; 1763 list hne=hnexpansion((x6-y10)*(x+y2-y3)*(x+y2+y3));1764 multsequence( hne[1])," | ",multsequence(hne[2])," | ",1765 multsequence( hne[3])," | ",multsequence(hne[4]);1766 multsequence( hne);1763 list Hne=hnexpansion((x6-y10)*(x+y2-y3)*(x+y2+y3)); 1764 multsequence(Hne[1])," | ",multsequence(Hne[2])," | ", 1765 multsequence(Hne[3])," | ",multsequence(Hne[4]); 1766 multsequence(Hne); 1767 1767 // The meaning of the entries of the 2nd matrix is as follows: 1768 displayMultsequence( hne);1768 displayMultsequence(Hne); 1769 1769 } 1770 1770 /////////////////////////////////////////////////////////////////////////////// … … 2542 2542 "EXAMPLE:"; echo = 2; 2543 2543 ring exring=0,(x,y),dp; 2544 list hne=hnexpansion(x14-3y2x11-y3x10-y2x9+3y4x8+y5x7+3y4x6+x5*(-y6+y5)2544 list Hne=hnexpansion(x14-3y2x11-y3x10-y2x9+3y4x8+y5x7+3y4x6+x5*(-y6+y5) 2545 2545 -3y6x3-y7x2+y8); 2546 displayHNE( hne); // HNE of 1st,3rd branch is finite2547 print(extdevelop( hne[1],5)[1]);2548 list ehne=extdevelop( hne[2],5);2546 displayHNE(Hne); // HNE of 1st,3rd branch is finite 2547 print(extdevelop(Hne[1],5)[1]); 2548 list ehne=extdevelop(Hne[2],5); 2549 2549 displayHNE(ehne); 2550 param( hne[2]);2550 param(Hne[2]); 2551 2551 param(ehne); 2552 2552 … … 2585 2585 "EXAMPLE:"; echo = 2; 2586 2586 ring r=0,(x,y),dp; 2587 list hne=develop(x2+y3+y4);2588 hne;2589 stripHNE( hne);2587 list Hne=develop(x2+y3+y4); 2588 Hne; 2589 stripHNE(Hne); 2590 2590 } 2591 2591 /////////////////////////////////////////////////////////////////////////////// … … 2854 2854 ring r=0,(x,y),dp; 2855 2855 // First, an example which requires no field extension: 2856 list hne=hnexpansion(x4-y6);2857 size( hne); // number of branches2858 displayHNE( hne); // HN expansion of branches2859 param( hne[1]); // parametrization of 1st branch2860 param( hne[2]); // parametrization of 2nd branch2856 list Hne=hnexpansion(x4-y6); 2857 size(Hne); // number of branches 2858 displayHNE(Hne); // HN expansion of branches 2859 param(Hne[1]); // parametrization of 1st branch 2860 param(Hne[2]); // parametrization of 2nd branch 2861 2861 2862 2862 // An example which requires a field extension:
Note: See TracChangeset
for help on using the changeset viewer.