Changeset 4ad2bb in git


Ignore:
Timestamp:
Jun 22, 1999, 5:39:03 PM (25 years ago)
Author:
Mathias Schulze <mschulze@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
97a7b44e9344c736a265288ae4f2fa5ab105da59
Parents:
1eb7af12ba2eda02be4f55b12705297ff94d3b13
Message:
*** empty log message ***


git-svn-id: file:///usr/local/Singular/svn/trunk@3161 2c84dea3-7e68-4137-9b89-c4e89433aadc
Location:
Singular/LIB
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • Singular/LIB/jordan.lib

    r1eb7af r4ad2bb  
    11///////////////////////////////////////////////////////////////////////////////
    22
    3 version="$Id: jordan.lib,v 1.12 1999-06-21 08:45:41 mschulze Exp $";
     3version="$Id: jordan.lib,v 1.13 1999-06-22 15:39:03 mschulze Exp $";
    44info="
    55LIBRARY: jordan.lib  PROCEDURES TO COMPUTE THE JORDAN NORMAL FORM
    … …  
    7070         By default, opt=0.
    7171NOTE:    A non constant polynomial matrix M is replaced by its constant part.
     72DISPLAY: The procedure displays comments if printlevel>=1.
    7273EXAMPLE: example jordan; shows an example.
    7374"
    … …  
    386387RETURN:  The procedure returns the Jordan matrix J with eigenvalues jd[1] and
    387388         size jd[2][i][j] of j-th Jordan block with eigenvalue jd[1][i].
     389DISPLAY: The procedure displays comments if printlevel>=1.
    388390EXAMPLE: example jordanmatrix; shows an example.
    389391"
    … …  
    477479RETURN:  The procedure returns the Jordan normal form of M.
    478480NOTE:    A non constant polynomial matrix M is replaced by its constant part.
     481DISPLAY: The procedure displays more comments for higher printlevel.
    479482EXAMPLE: example jordanform; shows an example.
    480483"

  • Singular/LIB/mondromy.lib

    r1eb7af r4ad2bb  
    11///////////////////////////////////////////////////////////////////////////////
    22
    3 version="$Id: mondromy.lib,v 1.2 1999-06-22 14:50:49 mschulze Exp $";
     3version="$Id: mondromy.lib,v 1.3 1999-06-22 15:38:54 mschulze Exp $";
    44info="
    55LIBRARY: mondromy.lib  PROCEDURES TO COMPUTE THE MONODROMY OF A SINGULARITY
    … …  
    149149RETURN:  The procedure returns the series inverse of u up to order n
    150150         or a zero polynomial if u is no series unit.
    151 DISPLAY: printlevel>=1; shows comments.
     151DISPLAY: The procedure displays comments if printlevel>=1.
    152152EXAMPLE: example invunit; shows an example.
    153153"
    … …  
    164164    poly v=lift(fetch(br,u),1)[1,1];
    165165    dbprint(printlevel-voice+2,"//...inverse computed ["+string(timer-t)+
    166       " secs]");
     166      " secs, "+string(memory(1))+" bytes]");
    167167
    168168    setring br;
    … …  
    190190         If U is a square matrix and the determinant of U not zero,
    191191         then the second entry is the adjoint matrix of U.
    192 DISPLAY: printlevel>=1; shows comments.
     192DISPLAY: The procedure displays comments if printlevel>=1.
    193193EXAMPLE: example detadj; shows an example.
    194194"
    … …  
    200200    poly detU=det(U);
    201201    dbprint(printlevel-voice+2,"//...determinant computed ["+string(timer-t)+
    202       " secs]");
     202      " secs, "+string(memory(1))+" bytes]");
    203203
    204204    if(detU==0)
    … …  
    218218      matrix adjU=lift(U,detU*freemodule(nrows(U)));
    219219      dbprint(printlevel-voice+2,"//...adjoint matrix computed ["
    220         +string(timer-t)+" secs]");
     220        +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
    221221
    222222      setring br;
    … …  
    250250         int, vector, poly such that kappa is minimal with f^kappa in jacob(f),
    251251         u is a unit, and u*f^kappa=(matrix(jacob(f))*xi)[1,1].
    252 DISPLAY: printlevel>=1; shows comments.
     252DISPLAY: The procedure displays comments if printlevel>=1.
    253253EXAMPLE: example jacoblift; shows an example.
    254254"
    … …  
    267267  }
    268268  dbprint(printlevel-voice+2,"//kappa="+string(kappa));
    269   dbprint(printlevel-voice+2,"//...kappa computed ["+string(timer-t)+" secs]");
     269  dbprint(printlevel-voice+2,"//...kappa computed ["+string(timer-t)+" secs, "
     270    +string(memory(1))+" bytes]");
    270271
    271272  dbprint(printlevel-voice+2,"//computing xi...");
    272273  t=timer;
    273274  vector xi=lift(jf,fkappa)[1];
    274   dbprint(printlevel-voice+2,"//...xi computed ["+string(timer-t)+" secs]");
     275  dbprint(printlevel-voice+2,"//...xi computed ["+string(timer-t)+" secs, "
     276    +string(memory(1))+" bytes]");
    275277
    276278  dbprint(printlevel-voice+2,"//computing u...");
    277279  t=timer;
    278280  poly u=(matrix(jf)*xi)[1,1]/fkappa;
    279   dbprint(printlevel-voice+2,"//...u computed ["+string(timer-t)+" secs]");
     281  dbprint(printlevel-voice+2,"//...u computed ["+string(timer-t)+" secs, "
     282    +string(memory(1))+" bytes]");
    280283
    281284  return(list(kappa,xi,u));
    … …  
    362365  {
    363366    dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t)
    364       +" secs]");
     367      +" secs, "+string(memory(1))+" bytes]");
    365368
    366369    deltaP1=getdeltaP1(f,K,N,deltaN);
    … …  
    383386  }
    384387  dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t)
    385     +" secs]");
     388    +" secs, "+string(memory(1))+" bytes]");
    386389
    387390  return(K,N,P1,P2,Pe,V1,V2,Ve);
    … …  
    423426  Vnablae=nablaK(f,kappa,xi,u,N,prevN,Vnablae,e);
    424427  dbprint(printlevel-voice+2,"//...nabla(e) computed ["+string(timer-t)
    425     +" secs]");
     428    +" secs, "+string(memory(1))+" bytes]");
    426429
    427430  dbprint(printlevel-voice+2,
    … …  
    434437  matrix C=lift(W,module(Vnablae[1..size(Vnablae)]));
    435438  dbprint(printlevel-voice+2,"//...nabla(e) lifted ["+string(timer-t)
    436     +" secs]");
     439    +" secs, "+string(memory(1))+" bytes]");
    437440
    438441  dbprint(printlevel-voice+2,"//computing e-lift of nabla(e)...");
    … …  
    451454  }
    452455  dbprint(printlevel-voice+2,"//...e-lift of nabla(e) computed ["
    453     +string(timer-t)+" secs]");
     456    +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
    454457
    455458  return(M,N,Vnablae);
    … …  
    493496    U=syz(M0e)+U;
    494497  }
    495   dbprint(printlevel-voice+2,"//...U computed ["+string(timer-t)+" secs]");
     498  dbprint(printlevel-voice+2,"//...U computed ["+string(timer-t)+" secs, "
     499    +string(memory(1))+" bytes]");
    496500
    497501  dbprint(printlevel-voice+2,"//transforming M to U...");
    … …  
    500504  daU[2]=(1/number(daU[1]))*daU[2];
    501505  M=daU[2]*M*U;
    502   dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs]");
     506  dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, "
     507    +string(memory(1))+" bytes]");
    503508
    504509  dbprint(printlevel-voice+2,
    … …  
    536541  dbprint(printlevel-voice+2,
    537542    "//...integer differences of eigenvalues of M0 computed ["+string(timer-t)
    538     +" secs]");
     543    +" secs, "+string(memory(1))+" bytes]");
    539544
    540545  dbprint(printlevel-voice+2,"//transforming M...");
    … …  
    561566    }
    562567  }
    563   dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs]");
     568  dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, "
     569    +string(memory(1))+" bytes]");
    564570
    565571  return(M);
    … …  
    602608  int t=timer;
    603609  e=pcvcv2p(quotV(V1+V2,N),0,N);
    604   dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs]");
     610  dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
     611    +string(memory(1))+" bytes]");
    605612
    606613  dbprint(printlevel-voice+2,"//e=");
    … …  
    631638    {
    632639      dbprint(printlevel-voice+2,"//...compared with previous lattice ["
    633         +string(timer-t)+" secs]");
     640        +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
    634641
    635642      dbprint(printlevel-voice+2,"//increasing K and N...");
    … …  
    649656        module(mdivp(M*U,var(1)^(kappa-1))),(kappa-1)*(mu-1)));
    650657      dbprint(printlevel-voice+2,"//...lattice enlarged ["+string(timer-t)
    651         +" secs]");
     658        +" secs, "+string(memory(1))+" bytes]");
    652659
    653660      dbprint(printlevel-voice+2,"//comparing with previous lattice...");
    … …  
    655662    }
    656663    dbprint(printlevel-voice+2,"//...compared with previous lattice ["
    657       +string(timer-t)+" secs]");
     664      +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
    658665    dbprint(printlevel-voice+2,"//...t*nabla-stable lattice computed");
    659666
    … …  
    666673      dbprint(printlevel-voice+2,
    667674        "//...C{f}-basis of t*nabla-stable lattice computed ["+string(timer-t)
    668         +" secs]");
     675        +" secs, "+string(memory(1))+" bytes]");
    669676    }
    670677
    … …  
    699706      leadcoef(daU[1])*var(1)^(kappa+pcvmindeg(daU[1])-1));
    700707    dbprint(printlevel-voice+2,"//...M/t^kappa transformed to simple pole ["
    701       +string(timer-t)+" secs]");
     708      +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
    702709  }
    703710
    … …  
    711718    int delta=mid(eM0);
    712719    dbprint(printlevel-voice+2,"//...delta computed ["+string(timer-t)
    713       +" secs]");
     720      +" secs, "+string(memory(1))+" bytes]");
    714721
    715722    dbprint(printlevel-voice+2,"//delta="+string(delta));
    … …  
    743750        dbprint(printlevel-voice+2,
    744751          "//...M/t^kappa transformed to simple pole ["+string(timer-t)
    745           +" secs]");
     752          +" secs, "+string(memory(1))+" bytes]");
    746753      }
    747754
    … …  
    776783  int t=timer;
    777784  ideal e=kbase(std(jacob(f)));
    778   dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs]");
     785  dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
     786    +string(memory(1))+" bytes]");
    779787
    780788  dbprint(printlevel-voice+2,
    … …  
    812820         Ann. Inst. Fourier, Grenoble 23,1 (1973), pp. 157-195) to transform
    813821         it to a simple pole.
    814 DISPLAY: printlevel>=1; shows comments.
     822DISPLAY: The procedure displays more comments for higher printlevel.
    815823EXAMPLE: example monodromy; shows an example.
    816824"
    … …  
    856864    int t=timer;
    857865    int mu=milnor(f);
    858     dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs]");
     866    dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, "
     867      +string(memory(1))+" bytes]");
    859868
    860869    dbprint(printlevel-voice+2,"//mu="+string(mu));
    … …  
    913922         Brieskorn lattice H''=Omega^(n+1)/df^dOmega^(n-1).
    914923THEORY:  H'' is a free C{f}-module of rank milnor(f).
    915 DISPLAY: printlevel>=1; shows comments.
     924DISPLAY: The procedure displays more comments for higher printlevel.
    916925EXAMPLE: example H''basis; shows an example.
    917926"
    … …  
    943952    int t=timer;
    944953    int mu=milnor(f);
    945     dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs]");
     954    dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, "
     955      +string(memory(1))+" bytes]");
    946956
    947957    dbprint(printlevel-voice+2,"//mu="+string(mu));
    … …  
    9941004      t=timer;
    9951005      e=pcvcv2p(quotV(V1+V2,N),0,N);
    996       dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs]");
     1006      dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
     1007        +string(memory(1))+" bytes]");
    9971008
    9981009      dbprint(printlevel-voice+2,"//e=");
Note: See TracChangeset for help on using the changeset viewer.