Home Online Manual
Top
Back: displayHNE
Forward: displayInvariants
FastBack: gmspoly_lib
FastForward: kskernel_lib
Up: hnoether_lib
Top: Singular Manual
Contents: Table of Contents
Index: Index
About: About this document

D.6.9.6 invariants

Procedure from library hnoether.lib (see hnoether_lib).

Usage:
invariants(INPUT); INPUT list or poly

Assume:
INPUT is the output of develop(f), or of extdevelop(develop(f),n), or one entry of the list of HN data computed by hnexpansion(f[,"ess"]).

Return:
list INV of the following format:
 
    INV[1]:  intvec    (characteristic exponents)
    INV[2]:  intvec    (generators of the semigroup)
    INV[3]:  intvec    (Puiseux pairs, 1st components)
    INV[4]:  intvec    (Puiseux pairs, 2nd components)
    INV[5]:  int       (degree of the conductor)
    INV[6]:  intvec    (sequence of multiplicities)
If INPUT contains no valid HN expansion, the empty list is returned.

Assume:
INPUT is a bivariate polynomial f, or the output of hnexpansion(f), or the list of HN data computed by hnexpansion(f [,"ess"]).

Return:
list INV, such that INV[i] coincides with the output of invariants(develop(f[i])), where f[i] is the i-th branch of f, and the last entry of INV contains further invariants of f in the format:
 
    INV[last][1] : intmat    (contact matrix of the branches)
    INV[last][2] : intmat    (intersection multiplicities of the branches)
    INV[last][3] : int       (delta invariant of f)

Note:
In case the Hamburger-Noether expansion of the curve f is needed for other purposes as well it is better to calculate this first with the aid of hnexpansion and use it as input instead of the polynomial itself.

Example:
 
LIB "hnoether.lib";
ring exring=0,(x,y),dp;
list Hne=develop(y4+2x3y2+x6+x5y);
list INV=invariants(Hne);
INV[1];                   // the characteristic exponents
==> 4,6,7
INV[2];                   // the generators of the semigroup of values
==> 4,6,13
INV[3],INV[4];            // the Puiseux pairs in packed form
==> 3,7 2,2
INV[5] div 2;             // the delta-invariant
==> 8
INV[6];                   // the sequence of multiplicities
==> 4,2,2,1,1
// To display the invariants more 'nicely':
displayInvariants(Hne);
==>  characteristic exponents  : 4,6,7
==>  generators of semigroup   : 4,6,13
==>  Puiseux pairs             : (3,2)(7,2)
==>  degree of the conductor   : 16
==>  delta invariant           : 8
==>  sequence of multiplicities: 4,2,2,1,1
/////////////////////////////
INV=invariants((x2-y3)*(x3-y5));
INV[1][1];                // the characteristic exponents of the first branch
==> 2,3
INV[2][6];                // the sequence of multiplicities of the second branch
==> 3,2,1,1
print(INV[size(INV)][1]);         // the contact matrix of the branches
==>      0     3
==>      3     0
print(INV[size(INV)][2]);         // the intersection numbers of the branches
==>      0     9
==>      9     0
INV[size(INV)][3];                // the delta invariant of the curve
==> 14
See also: develop; displayInvariants; hnexpansion; intersection; multsequence.