Ignore:
Timestamp:
Jan 16, 2019, 6:43:39 PM (5 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'spielwiese', '8e0ad00ce244dfd0756200662572aef8402f13d5')
Children:
ecf1dbe0e95a9a6bc8b3ea517b7f9643198a5c5f
Parents:
3fffaf51c086050aa295bf603497ddd25fcdc14d
Message:
format (arnoldclassify.lib,SingularityDBM.lib)
File:
1 edited

Legend:

Unmodified
Added
Removed

  • Singular/LIB/SingularityDBM.lib

    r3fffaf5 rab758e  
    11/////////////////////////////////////////////////////////////////////////////
    2 version="version singularityDBM.lib 4.1.0.0 Sep_2017 "; //
     2version="version singularityDBM.lib 4.1.0.0 Sep_2017 "; // $Id$
    33category="Singularities";
    44info="
    … …  
    1818static proc makedbm_init()
    1919{
    20 //* Generates file containing a data base for singularities up to corank 2 
     20//* Generates file containing a data base for singularities up to corank 2
    2121//* listed by Arnol'd. This file is needed for arnoldclassify.lib.
    2222
    … …  
    5858"USAGE:   create_singularity_dbm();
    5959COMPUTE: Generates two files,  Singularitylist.dir and Singularitylist.pag
    60          containing a data base for singularities up to corank 2 listed 
     60         containing a data base for singularities up to corank 2 listed
    6161                                 by Arnol'd.
    6262RETURN:  Nothing
    … …  
    6868  string s;
    6969
    70 //* A[k]       
    71         s = "singseries f;
    72         f.Series = \"A[k]\"; 
     70//* A[k]
     71        s = "singseries f;
     72        f.Series = \"A[k]\";
    7373        f.Modality = \"0\";
    7474        f.Corank = \"1\";
    7575        f.MilnorNumber = \"k\";
    76         f.MilnorCode = \"k\"; 
     76        f.MilnorCode = \"k\";
    7777        f.NormalForm = \"x^(k+1)\";
    7878        f.SpecialForm = \"x^(k+1)\";
    7979        f.Restrictions = \"(k>1)\";";
    8080  write(l, "A[k]", s);
    81        
    82 //* D[k]       
    83         s = "singseries f;
    84         f.Series = \"D[k]\"; 
     81
     82//* D[k]
     83        s = "singseries f;
     84        f.Series = \"D[k]\";
    8585        f.Modality = \"0\";
    8686        f.Corank = \"2\";
    8787        f.MilnorNumber = \"k\";
    88         f.MilnorCode = \"1,1,k-3\"; 
     88        f.MilnorCode = \"1,1,k-3\";
    8989        f.NormalForm = \"x^2*y+y^(k-1)\";
    9090        f.SpecialForm = \"x^2*y+y^(k-1)\";
    … …  
    9494//* J[k,0]
    9595        s = "singseries f;
    96         f.Series = \"J[k,0]\"; 
     96        f.Series = \"J[k,0]\";
    9797        f.Modality = \"0\";
    9898        f.Corank = \"2\";
    9999        f.MilnorNumber = \"6*k-2\";
    100         f.MilnorCode = \"1,2*k+j,2*k-2*j-3\"; 
     100        f.MilnorCode = \"1,2*k+j,2*k-2*j-3\";
    101101        f.NormalForm = \"x^3 + b(y)*x^2*y^k+c(y)*x*y^(2*k+1)+y^(3*k)\";
    102102        f.SpecialForm = \"x^3 + x^2*y^k+y^(3*k)\";
    … …  
    107107//* J[k,r]
    108108        s = "singseries f;
    109         f.Series = \"J[k,r]\"; 
     109        f.Series = \"J[k,r]\";
    110110        f.Modality = \"0\";
    111111        f.Corank = \"2\";
    112112        f.MilnorNumber = \"6*k-2+r\";
    113         f.MilnorCode = \"1,2*k-1,2*k+r-1\"; 
     113        f.MilnorCode = \"1,2*k-1,2*k+r-1\";
    114114        f.NormalForm = \"x^3 + x^2*y^k+a(y)*y^(3*k+r)\";
    115115        f.SpecialForm = \"x^3 + x^2*y^k+y^(3*k+r)\";
    … …  
    119119//* E[6k]
    120120        s = "singseries f;
    121         f.Series = \"E[6k]\"; 
     121        f.Series = \"E[6k]\";
    122122        f.Modality = \"0\";
    123123        f.Corank = \"2\";
    124124        f.MilnorNumber = \"6*k\";
    125         f.MilnorCode = \"1,2*k+j,2*k-2j-1\"; 
     125        f.MilnorCode = \"1,2*k+j,2*k-2j-1\";
    126126        f.NormalForm = \"x^3 + a(y)*x*y^(2*k+1)+y^(3*k+1)\";
    127127        f.SpecialForm = \"x^3+y^(3*k+1)\";
    … …  
    129129  write(l, "E[6k]", s);
    130130
    131 //* E[6k+1]     
    132         s = "singseries f;
    133         f.Series = \"E[6k+1]\"; 
     131//* E[6k+1]
     132        s = "singseries f;
     133        f.Series = \"E[6k+1]\";
    134134        f.Modality = \"0\";
    135135        f.Corank = \"2\";
    136136        f.MilnorNumber = \"6*k+1\";
    137         f.MilnorCode = \"1,2*k,2*k\"; 
     137        f.MilnorCode = \"1,2*k,2*k\";
    138138        f.NormalForm = \"x^3 + x*y^(2*k+1)+a(y)*y^(3*k+2)\";
    139139        f.SpecialForm = \"x^3 + x*y^(2*k+1)\";
    … …  
    141141  write(l, "E[6k+1]", s);
    142142
    143 //* E[6k+2]     
    144         s = "singseries f;
    145         f.Series = \"E[6k+2]\"; 
     143//* E[6k+2]
     144        s = "singseries f;
     145        f.Series = \"E[6k+2]\";
    146146        f.Modality = \"0\";
    147147        f.Corank = \"2\";
    148148        f.MilnorNumber = \"6*k+2\";
    149         f.MilnorCode = \"1,2*k+j+1,2*k-2j-1\"; 
     149        f.MilnorCode = \"1,2*k+j+1,2*k-2j-1\";
    150150        f.NormalForm = \"x^3 + a(y)*x*y^(2*k+2)+y^(3*k+2)\";
    151151        f.SpecialForm = \"x^3 +y^(3*k+2)\";
    … …  
    155155//* X[k,0]
    156156        s = "singseries f;
    157         f.Series = \"X[k,0]\"; 
     157        f.Series = \"X[k,0]\";
    158158        f.Modality = \"3*k-2\";
    159159        f.Corank = \"2\";
    160160        f.MilnorNumber = \"12*k-3\";
    161         f.MilnorCode = \"1,1,2*k-1+j,2k-1-2*j+t,2*k-1+j-2t\"; 
     161        f.MilnorCode = \"1,1,2*k-1+j,2k-1-2*j+t,2*k-1+j-2t\";
    162162        f.NormalForm = \"x^4 + b(y)*x^3*y^k + a(y)*x^2*y^(2*k) + x*y^(3*k)\";
    163163        f.SpecialForm = \"x^4 + x^3*y^k + x*y^(3*k)\";
    164164        f.Restrictions = \"(jet(a,0)*jet(b,0)!=9)&&(k>1)&&(4*(jet(a,0)^3+jet(b,0)^3)
    165          - jet(a,0)^2*jet(b,0)^2-18* jet(a,0)*jet(b,0) + 27 !=0)&&(deg(a)<=(k-2)) 
     165         - jet(a,0)^2*jet(b,0)^2-18* jet(a,0)*jet(b,0) + 27 !=0)&&(deg(a)<=(k-2))
    166166         &&(deg(b)<=(2*k-2))\";";
    167167  write(l, "X[k,0]", s);
    168        
     168
    169169        //* X[1,0]
    170170                s = "singseries f;
    171                 f.Series = \"X[1,0]\"; 
     171                f.Series = \"X[1,0]\";
    172172                f.Modality = \"1\";
    173173                f.Corank = \"2\";
    174174                f.MilnorNumber = \"9\";
    175                 f.MilnorCode = \"1,1,1+j,1-2*j+t,1+j-2t\"; 
     175                f.MilnorCode = \"1,1,1+j,1-2*j+t,1+j-2t\";
    176176                f.NormalForm = \"x^4 + a(y)*x^2*y^2 + y^4\";
    177177                f.SpecialForm = \"x^4 + x^2*y^2 + y^4\";
    178178                f.Restrictions = \"(deg(a)==0)&&(jet(a,0)^2!=4)\";";
    179179          write(l, "X[1,0]", s);
    180        
    181         //* X[k,r]     
    182                 s = "singseries f;
    183                 f.Series = \"X[k,r]\"; 
     180
     181        //* X[k,r]
     182                s = "singseries f;
     183                f.Series = \"X[k,r]\";
    184184                f.Modality = \"3*k-2\";
    185185                f.Corank = \"2\";
    186186                f.MilnorNumber = \"12*k-3+r\";
    187                 f.MilnorCode = \"1,1,2*k-1+j,2k-1-2*j,2*k-1+j+r\"; 
     187                f.MilnorCode = \"1,1,2*k-1+j,2k-1-2*j,2*k-1+j+r\";
    188188                f.NormalForm = \"x4+a(y)*x3*y^(k)+x^2*y^(2*k)+b(y)*y^(4*k+r)\";
    189189                f.SpecialForm = \"x4+x3*y^(k)+x^2*y^(2*k)+y^(4*k+r)\";
    … …  
    191191                (jet(b,0)!=0)&&(deg(b)<=(2*k-2))\";";
    192192          write(l, "X[k,r]", s);
    193                
     193
    194194        //* X[1,r]
    195195                s = "singseries f;
    196                 f.Series = \"X[1,r]\"; 
     196                f.Series = \"X[1,r]\";
    197197                f.Modality = \"1\";
    198198                f.Corank = \"2\";
    199199                f.MilnorNumber = \"9+r\";
    200                 f.MilnorCode = \"1,1,1+j,1-2*j,1+j+r\"; 
     200                f.MilnorCode = \"1,1,1+j,1-2*j,1+j+r\";
    201201                f.NormalForm = \"x4+x^2*y^2+a(y)*y^(4+r)\";
    202202                f.SpecialForm = \"x4+x^2*y^2+y^(4+r)\";
    203203                f.Restrictions = \"(deg(a)==0)&&(jet(a,0)!=0)\";";
    204204                write(l, "X[1,r]", s);
    205                
     205
    206206        //* Y[k,r,s]
    207207                s = "singseries f;
    208                 f.Series = \"Y[k,r,s]\"; 
     208                f.Series = \"Y[k,r,s]\";
    209209                f.Modality = \"3*k-2\";
    210210                f.Corank = \"2\";
    211211                f.MilnorNumber = \"12*k-3+r+s\";
    212                 f.MilnorCode = \"1,1,2*k-1,2*k-1+j,2*k-1-2*j+r+s\"; 
     212                f.MilnorCode = \"1,1,2*k-1,2*k-1+j,2*k-1-2*j+r+s\";
    213213                f.NormalForm = \"((x + a(y)*y^k)^2 + b(y)*y^(2*k+s))*(x2 + y^(2*k+r))\";
    214214                f.SpecialForm = \"((x + y^k)^2 + y^(2*k+s))*(x2 + y^(2*k+r))\";
    … …  
    216216                &&(1<=s)&&(s<=7)\";";
    217217                write(l, "Y[k,r,s]", s);
    218        
     218
    219219        //* Y[1,r,s]
    220220                s = "singseries f;
    221                 f.Series = \"Y[1,r,s]\"; 
     221                f.Series = \"Y[1,r,s]\";
    222222                f.Modality = \"1\";
    223223                f.Corank = \"2\";
    224224                f.MilnorNumber = \"9+r+s\";
    225                 f.MilnorCode = \"1,1,1,1+j,1-2*j+r+s\"; 
     225                f.MilnorCode = \"1,1,1,1+j,1-2*j+r+s\";
    226226                f.NormalForm = \" x^(4+r)+ a(y)*x2*y2 + y^(4+s)\";
    227227                f.SpecialForm = \" x^(4+r)+ x2*y2 + y^(4+s)\";
    228228                f.Restrictions = \"(deg(a)==0)&&(jet(a,0)!=0)&&(1<=s)&&(s<=7)\";";
    229229                write(l, "Y[1,r,s]", s);
    230                
     230
    231231        //* Z[k,r]
    232232                s = "singseries f;
    233                 f.Series = \"Z[k,r]\"; 
     233                f.Series = \"Z[k,r]\";
    234234                f.Modality = \"3*k+r-2\";
    235235                f.Corank = \"2\";
    236236                f.MilnorNumber = \"12*k-3+6*r\";
    237                 f.MilnorCode = \"1,1,2*k-1,2*k-1+j,2*k-1+6*r-2*j\"; 
    238                 f.NormalForm = \"(x + a(y)*y^k)*(x^3 + d(y)*x2*y^(k+1) + 
     237                f.MilnorCode = \"1,1,2*k-1,2*k-1+j,2*k-1+6*r-2*j\";
     238                f.NormalForm = \"(x + a(y)*y^k)*(x^3 + d(y)*x2*y^(k+1) +
    239239                c(y)*x*y^(2*k+2*r+1) + y^(3*k+3*r))\";
    240                 f.SpecialForm = \"(x + y^k)*(x^3 + 2*y^(k+1) + x*y^(2*k+2*r+1) + 
     240                f.SpecialForm = \"(x + y^k)*(x^3 + 2*y^(k+1) + x*y^(2*k+2*r+1) +
    241241                y^(3*k+3*r))\";
    242242                f.Restrictions = \"(k>1)&&(r>=0)&&(4*d^3+27!=0)&&(deg(d)==0)&&
    243243                (deg(c)<=(2*k+r-3))&&(deg(a)<=(k-2))\";";
    244244                write(l, "Z[k,r]", s);
    245                
     245
    246246        //* Z[1,r]
    247247                s = "singseries f;
    248                 f.Series = \"Z[1,r]\"; 
     248                f.Series = \"Z[1,r]\";
    249249                f.Modality = \"1+r\";
    250250                f.Corank = \"2\";
    251251                f.MilnorNumber = \"9+6*r\";
    252                 f.MilnorCode = \"1,1,1,1+j,1+6*r-2*j\"; 
    253                 f.NormalForm = \"y*(x^3 + d(y)*x^2*y^(2) + c(y)*x*y^(2+2*r+1) + 
     252                f.MilnorCode = \"1,1,1,1+j,1+6*r-2*j\";
     253                f.NormalForm = \"y*(x^3 + d(y)*x^2*y^(2) + c(y)*x*y^(2+2*r+1) +
    254254                y^(3+3*r))\";
    255                 f.SpecialForm = \"y*(x^3 + x^2*y^(2) + x*y^(2+2*r+1) + 
     255                f.SpecialForm = \"y*(x^3 + x^2*y^(2) + x*y^(2+2*r+1) +
    256256                y^(3+3*r))\";
    257257                f.Restrictions = \"(r>=0)&&(4*d^3+27!=0)&&(deg(d)==0)
    … …  
    261261        //* Z[k,r,s]
    262262                s = "singseries f;
    263                 f.Series = \"Z[k,r,s]\"; 
     263                f.Series = \"Z[k,r,s]\";
    264264                f.Modality = \"3*k+r-2\";
    265265                f.Corank = \"2\";
    266266                f.MilnorNumber = \"12*k+6*r+s-3\";
    267                 f.MilnorCode = \"1,1,2*k-1,2*k-1+2*r,2*k-1+2*r-s\"; 
     267                f.MilnorCode = \"1,1,2*k-1,2*k-1+2*r,2*k-1+2*r-s\";
    268268                f.NormalForm = \"(x^2 + a(y)*x*y^k + b(y)*y^(2*k+r))*
    269269                (x^2 + y^(2*k+2*r+s))\";
    … …  
    272272                (jet(b,0)!=0)&&(deg(b)<=(2*k+r-2))\";";
    273273                write(l, "Z[k,r,s]", s);
    274                
     274
    275275        //* Z[1,r,s]
    276276                s = "singseries f;
    277                 f.Series = \"Z[1,r,s]\"; 
     277                f.Series = \"Z[1,r,s]\";
    278278                f.Modality = \"1+r\";
    279279                f.Corank = \"2\";
    280280                f.MilnorNumber = \"9+6*r+s\";
    281                 f.MilnorCode = \"1,1,1,1+2*r,1+2*r-s\"; 
     281                f.MilnorCode = \"1,1,1,1+2*r,1+2*r-s\";
    282282                f.NormalForm = \"y*(x^3 + x^2*y^(r+1) + b(y)*y^(3*r+s+3))\";
    283283                f.SpecialForm = \"y*(x^3 + x^2*y^(r+1) + y^(3*r+s+3))\";
    284284                f.Restrictions = \"(r>=0)&&(jet(b,0)!=0)&&(deg(b)<=(2*k+r-2))\";";
    285285                write(l, "Z[1,r,s]", s);
    286                
    287         //* Z[k,12k+6r-1]       
    288                 s = "singseries f;
    289                 f.Series = \"Z[k,12k+6r-1]\"; 
     286
     287        //* Z[k,12k+6r-1]
     288                s = "singseries f;
     289                f.Series = \"Z[k,12k+6r-1]\";
    290290                f.Modality = \"3*k+r-2\";
    291291                f.Corank = \"2\";
    292292                f.MilnorNumber = \"12*k+6r-1\";
    293                 f.MilnorCode = \"1,1,2k-1,2k-1+j,2k+1+6*r-2*j\"; 
     293                f.MilnorCode = \"1,1,2k-1,2k-1+j,2k+1+6*r-2*j\";
    294294                f.NormalForm = \"(x + a(y)*y^k)*(x^3 + b(y)*x*y^(2*k+2*r+1) +
    295295                 y^(3*k+3*r+1))\";
    … …  
    298298                (jet(b,0)!=0)&&(deg(b)<=(2*k+r-2))\";";
    299299                write(l, "Z[k,12k+6r-1]", s);
    300                
    301                 //* Z[1,6r+11] 
     300
     301                //* Z[1,6r+11]
    302302                        s = "singseries f;
    303                         f.Series = \"Z[1,6r+11]\"; 
     303                        f.Series = \"Z[1,6r+11]\";
    304304                        f.Modality = \"1+r\";
    305305                        f.Corank = \"2\";
    306306                        f.MilnorNumber = \"6r+11\";
    307                         f.MilnorCode = \"1,1,1,1+j,3+6*r-2*j\"; 
     307                        f.MilnorCode = \"1,1,1,1+j,3+6*r-2*j\";
    308308                        f.NormalForm = \"y*(x^3 + b(y)*x*y^(2+2*r+1) + y^(3+3*r+1))\";
    309309                        f.SpecialForm = \"y*(x^3 + x*y^(2+2*r+1) + y^(3+3*r+1))\";
    310310                        f.Restrictions = \"(r>=0)&&(deg(b)<=(r))\";";
    311311                        write(l, "Z[1,6r+11]", s);
    312                
    313         //* Z[k,12k+6r+1]       
    314                 s = "singseries f;
    315                 f.Series = \"Z[k,12k+6r+1]\"; 
     312
     313        //* Z[k,12k+6r+1]
     314                s = "singseries f;
     315                f.Series = \"Z[k,12k+6r+1]\";
    316316                f.Modality = \"3*k+r-2\";
    317317                f.Corank = \"2\";
    318318                f.MilnorNumber = \"12*k+6r+1\";
    319                 f.MilnorCode = \"1,1,2k-1,2k-1+j,2k+3+6*r-2*j\"; 
     319                f.MilnorCode = \"1,1,2k-1,2k-1+j,2k+3+6*r-2*j\";
    320320                f.NormalForm = \"(x + a(y)*y^k)*(x^3 + b(y)*x*y^(2*k+2*r+2) +
    321321                 y^(3*k+3*r+2))\";
    … …  
    324324                (jet(b,0)!=0)&&(deg(b)<=(2*k+r-2))\";";
    325325                write(l, "Z[k,12k+6r+1]", s);
    326                
     326
    327327        //* Z[1,6r+13]
    328328                s = "singseries f;
    329                 f.Series = \"Z[1,6r+13]\"; 
     329                f.Series = \"Z[1,6r+13]\";
    330330                f.Modality = \"1+r\";
    331331                f.Corank = \"2\";
    332332                f.MilnorNumber = \"6r+13\";
    333                 f.MilnorCode = \"1,1,1,1+j,5+6*r-2*j\"; 
     333                f.MilnorCode = \"1,1,1,1+j,5+6*r-2*j\";
    334334                f.NormalForm = \"y*(x^3 + b(y)*x*y^(2*r+4) + y^(3*r+5))\";
    335335                f.SpecialForm = \"y*(x^3 + x*y^(2*r+4) + y^(3*r+5))\";
    336336                f.Restrictions = \" (r>=0)&&(deg(b)<=(r))\";";
    337337                write(l, "Z[1,6r+13]", s);
    338                
    339         //* Z[k,12k+6r] 
    340                 s = "singseries f;
    341                 f.Series = \"Z[k,12k+6r]\"; 
     338
     339        //* Z[k,12k+6r]
     340                s = "singseries f;
     341                f.Series = \"Z[k,12k+6r]\";
    342342                f.Modality = \"3*k+r-2\";
    343343                f.Corank = \"2\";
    344344                f.MilnorNumber = \"12*k+6r\";
    345                 f.MilnorCode = \"1,1,2k-1,2k-1+2*r,2k+2*r\"; 
     345                f.MilnorCode = \"1,1,2k-1,2k-1+2*r,2k+2*r\";
    346346                f.NormalForm = \"(x + a(y)*y^k)*(x^3 + x*y^(2*k+2*r+1) +
    347347                b(y)* y^(3*k+3*r+2))\";
    … …  
    350350                (jet(b,0)!=0)&&(deg(b)<=(2*k+r-2))\";";
    351351                write(l, "Z[k,12k+6r]", s);
    352                
    353                
     352
     353
    354354        //* Z[1,6r+12]
    355355                s = "singseries f;
    356                 f.Series = \"Z[1,6r+12]\"; 
     356                f.Series = \"Z[1,6r+12]\";
    357357                f.Modality = \"1+r\";
    358358                f.Corank = \"2\";
    359359                f.MilnorNumber = \"6*r+12\";
    360                 f.MilnorCode = \"1,1,1,1+2*r,2+2*r\"; 
     360                f.MilnorCode = \"1,1,1,1+2*r,2+2*r\";
    361361                f.NormalForm = \"y*(x^3 + x*y^(2*r+3) +b(y)* y^(3*r+5))\";
    362362                f.SpecialForm = \"y*(x^3 + x*y^(2*r+3) +y^(3*r+5))\";
    363363                f.Restrictions = \"(r>=0)&&(deg(b)<=(r))\";";
    364364                write(l, "Z[1,6r+12]", s);
    365                
    366                
     365
     366
    367367        //* W[k,r]
    368368                s = "singseries f;
    369                 f.Series = \"W[k,r]\"; 
     369                f.Series = \"W[k,r]\";
    370370                f.Modality = \"3*k-1\";
    371371                f.Corank = \"2\";
    372372                f.MilnorNumber = \"12*k+3+r\";
    373                 f.MilnorCode = \"1,1,2k,2k,2k+r\"; 
     373                f.MilnorCode = \"1,1,2k,2k,2k+r\";
    374374                f.NormalForm = \"x4+a(y)*x^3*y^(k+1)+x^2*y^(2*k+1)+b(y)*y^(4*k+2+r) \";
    375375                f.SpecialForm = \"x4+x^2*y^(2*k+1)+y^(4*k+2+r) \";
    … …  
    380380        //* W[k,0]
    381381                s = "singseries f;
    382                 f.Series = \"W[k,0]\"; 
     382                f.Series = \"W[k,0]\";
    383383                f.Modality = \"3*k-1\";
    384384                f.Corank = \"2\";
    385385                f.MilnorNumber = \"12*k+3\";
    386                 f.MilnorCode = \"1,1,2k+j,2k-2-2*j+t,2k+6+j+2*t\"; 
     386                f.MilnorCode = \"1,1,2k+j,2k-2-2*j+t,2k+6+j+2*t\";
    387387                f.NormalForm = \"x4+b(y)*x2*y^(2*k+1)+a(y)*x*y^(3*k+2)+y^(4*k+2) \";
    388388                f.SpecialForm = \"x4+x2*y^(2*k+1)+y^(4*k+2) \";
    … …  
    393393        //* W[12k]
    394394                s = "singseries f;
    395                 f.Series = \"W[12k]\"; 
     395                f.Series = \"W[12k]\";
    396396                f.Modality = \"3*k-2\";
    397397                f.Corank = \"2\";
    398398                f.MilnorNumber = \"12*k\";
    399                 f.MilnorCode = \"1,1,2k+j,2k-3-2*j+t,2k+3+j-2*t\"; 
     399                f.MilnorCode = \"1,1,2k+j,2k-3-2*j+t,2k+3+j-2*t\";
    400400                f.NormalForm = \"x4+a(y)*x*y^(3*k+1)+c(y)*x^2*y^(2*k+1)+y^(4*k+1)\";
    401401                f.SpecialForm = \"x4+x^2*y^(2*k+1)+y^(4*k+1)\";
    … …  
    403403                (deg(c)<=(2*k-2))\";";
    404404                write(l, "W[12k]", s);
    405                
     405
    406406        //* W[12k+1]
    407407                s = "singseries f;
    408                 f.Series = \"W[12k+1]\"; 
     408                f.Series = \"W[12k+1]\";
    409409                f.Modality = \"3*k-2\";
    410410                f.Corank = \"2\";
    411411                f.MilnorNumber = \"12*k+1\";
    412                 f.MilnorCode = \"1,1,2k+j,2k-1-2*j,2k+j\"; 
     412                f.MilnorCode = \"1,1,2k+j,2k-1-2*j,2k+j\";
    413413                f.NormalForm = \"x4+x*y^(3*k+1)+a(y)*x^2*y^(2*k+1)+c(y)*y^(4*k+2) \";
    414414                f.SpecialForm = \"x4+x*y^(3*k+1)+y^(4*k+2) \";
    … …  
    419419        //* W[12k+5]
    420420                s = "singseries f;
    421                 f.Series = \"W[12k+5]\"; 
     421                f.Series = \"W[12k+5]\";
    422422                f.Modality = \"3*k-1\";
    423423                f.Corank = \"2\";
    424424                f.MilnorNumber = \"12*k+5\";
    425                 f.MilnorCode = \"1,1,2k+j,2k+1-2*j,2k+j\"; 
     425                f.MilnorCode = \"1,1,2k+j,2k+1-2*j,2k+j\";
    426426                f.NormalForm = \"x4+x*y^(3*k+2)+a(y)*x^2*y^(2*k+2)+b(y)*y^(4*k+3) \";
    427427                f.SpecialForm = \"x4+x*y^(3*k+2)+y^(4*k+3) \";
    … …  
    429429                (deg(b)<=(2*k-1))\";";
    430430                write(l, "W[12k+5]", s);
    431                
     431
    432432        //* W[12k+6]
    433433                s = "singseries f;
    434                 f.Series = \"W[12k+6]\"; 
     434                f.Series = \"W[12k+6]\";
    435435                f.Modality = \"3*k-1\";
    436436                f.Corank = \"2\";
    437437                f.MilnorNumber = \"12*k+6\";
    438                 f.MilnorCode = \"1,1,2k+j,2k-3-2*j+t,2k+9+j-2*t\"; 
     438                f.MilnorCode = \"1,1,2k+j,2k-3-2*j+t,2k+9+j-2*t\";
    439439                f.NormalForm = \"x4+a(y)*x*y^(3*k+3)+b(y)*x^2*y^(2*k+2)+y^(4*k+3) \";
    440440                f.SpecialForm = \"x4+x^2*y^(2*k+2)+y^(4*k+3) \";
    … …  
    442442                (deg(b)<=(2*k-1))\";";
    443443                write(l, "W[12k+6]", s);
    444                
     444
    445445        //* W#[k,2r]
    446446                s = "singseries f;
    447                 f.Series = \"W#[k,2r]\"; 
     447                f.Series = \"W#[k,2r]\";
    448448                f.Modality = \"3*k-1\";
    449449                f.Corank = \"2\";
    450450                f.MilnorNumber = \"12*k+3+2*r\";
    451                 f.MilnorCode = \"1,1,2k,2k+r,2k\"; 
     451                f.MilnorCode = \"1,1,2k,2k+r,2k\";
    452452                f.NormalForm = \"(x2+y^(2*k+1))^2+b(y)*x^2*y^(2*k+1+r)+
    453453                a(y)*x*y^(3*k+2+r) \";
    … …  
    455455                f.Restrictions = \"(k>=1)&&(r>0)&&(k>1||a==0)&&(deg(a)<=(k-2))&&
    456456                (jet(b,0)!=0)&&(deg(b)<=(2*k-1))\";";
    457                 write(l, "W#[k,2r]", s);       
     457                write(l, "W#[k,2r]", s);
    458458
    459459        //* W#[k,2r-1]
    460460                s = "singseries f;
    461                 f.Series = \"W#[k,2r-1]\"; 
     461                f.Series = \"W#[k,2r-1]\";
    462462                f.Modality = \"3*k-1\";
    463463                f.Corank = \"2\";
    464464                f.MilnorNumber = \"12*k+2+2*r\";
    465                 f.MilnorCode = \"1,1,2k,2k-3+j,2*k+5+2*r-2*j\"; 
     465                f.MilnorCode = \"1,1,2k,2k-3+j,2*k+5+2*r-2*j\";
    466466                f.NormalForm = \"(x2+y^(2*k+1))^2+b(y)*x*y^(3*k+1+r)+
    467467                a(y)*y^(4*k+2+r)\";
    … …  
    469469                f.Restrictions = \"(k>=1)&&(r>0)&&(k>1||a==0)&&(deg(a)<=(k-2))
    470470                &&(jet(b,0)!=0)&&(deg(b)<=(2*k-1))\";";
    471                 write(l, "W#[k,2r-1]", s);             
    472                
     471                write(l, "W#[k,2r-1]", s);
     472
    473473 write(l,"VERSION", "1.0");
    474474  close(l);
    … …  
    484484  else { DatabasePath = "Singularitylist"; }
    485485  Database="DBM: ",DatabasePath;
    486        
     486
    487487  link dbmLink=Database;
    488488  Tp = read(dbmLink, typ);
Note: See TracChangeset for help on using the changeset viewer.