Changeset 97a7b44 in git


Ignore:
Timestamp:
Jun 22, 1999, 5:47:35 PM (25 years ago)
Author:
Mathias Schulze <mschulze@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b52fc4b2495505785981d640dcf7eb3e456778ef')
Children:
3e9ae9c78c81633291fa286789448020d1984470
Parents:
4ad2bb592a6ffe7f195320a5d4bec4f9ef0c57dd
Message:
*** empty log message ***


git-svn-id: file:///usr/local/Singular/svn/trunk@3162 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/mondromy.lib

    r4ad2bb r97a7b44  
    11///////////////////////////////////////////////////////////////////////////////
    22
    3 version="$Id: mondromy.lib,v 1.3 1999-06-22 15:38:54 mschulze Exp $";
     3version="$Id: mondromy.lib,v 1.4 1999-06-22 15:47:35 mschulze Exp $";
    44info="
    55LIBRARY: mondromy.lib  PROCEDURES TO COMPUTE THE MONODROMY OF A SINGULARITY
     
    164164    poly v=lift(fetch(br,u),1)[1,1];
    165165    dbprint(printlevel-voice+2,"//...inverse computed ["+string(timer-t)+
    166       " secs, "+string(memory(1))+" bytes]");
     166      " secs, "+string((memory(1)+1023)/1024)+" K]");
    167167
    168168    setring br;
     
    200200    poly detU=det(U);
    201201    dbprint(printlevel-voice+2,"//...determinant computed ["+string(timer-t)+
    202       " secs, "+string(memory(1))+" bytes]");
     202      " secs, "+string((memory(1)+1023)/1024)+" K]");
    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, "+string(memory(1))+" bytes]");
     220        +string(timer-t)+" secs, "+string((memory(1)+1023)/1024)+" K]");
    221221
    222222      setring br;
     
    268268  dbprint(printlevel-voice+2,"//kappa="+string(kappa));
    269269  dbprint(printlevel-voice+2,"//...kappa computed ["+string(timer-t)+" secs, "
    270     +string(memory(1))+" bytes]");
     270    +string((memory(1)+1023)/1024)+" K]");
    271271
    272272  dbprint(printlevel-voice+2,"//computing xi...");
     
    274274  vector xi=lift(jf,fkappa)[1];
    275275  dbprint(printlevel-voice+2,"//...xi computed ["+string(timer-t)+" secs, "
    276     +string(memory(1))+" bytes]");
     276    +string((memory(1)+1023)/1024)+" K]");
    277277
    278278  dbprint(printlevel-voice+2,"//computing u...");
     
    280280  poly u=(matrix(jf)*xi)[1,1]/fkappa;
    281281  dbprint(printlevel-voice+2,"//...u computed ["+string(timer-t)+" secs, "
    282     +string(memory(1))+" bytes]");
     282    +string((memory(1)+1023)/1024)+" K]");
    283283
    284284  return(list(kappa,xi,u));
     
    365365  {
    366366    dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t)
    367       +" secs, "+string(memory(1))+" bytes]");
     367      +" secs, "+string((memory(1)+1023)/1024)+" K]");
    368368
    369369    deltaP1=getdeltaP1(f,K,N,deltaN);
     
    386386  }
    387387  dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t)
    388     +" secs, "+string(memory(1))+" bytes]");
     388    +" secs, "+string((memory(1)+1023)/1024)+" K]");
    389389
    390390  return(K,N,P1,P2,Pe,V1,V2,Ve);
     
    426426  Vnablae=nablaK(f,kappa,xi,u,N,prevN,Vnablae,e);
    427427  dbprint(printlevel-voice+2,"//...nabla(e) computed ["+string(timer-t)
    428     +" secs, "+string(memory(1))+" bytes]");
     428    +" secs, "+string((memory(1)+1023)/1024)+" K]");
    429429
    430430  dbprint(printlevel-voice+2,
     
    437437  matrix C=lift(W,module(Vnablae[1..size(Vnablae)]));
    438438  dbprint(printlevel-voice+2,"//...nabla(e) lifted ["+string(timer-t)
    439     +" secs, "+string(memory(1))+" bytes]");
     439    +" secs, "+string((memory(1)+1023)/1024)+" K]");
    440440
    441441  dbprint(printlevel-voice+2,"//computing e-lift of nabla(e)...");
     
    454454  }
    455455  dbprint(printlevel-voice+2,"//...e-lift of nabla(e) computed ["
    456     +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
     456    +string(timer-t)+" secs, "+string((memory(1)+1023)/1024)+" K]");
    457457
    458458  return(M,N,Vnablae);
     
    497497  }
    498498  dbprint(printlevel-voice+2,"//...U computed ["+string(timer-t)+" secs, "
    499     +string(memory(1))+" bytes]");
     499    +string((memory(1)+1023)/1024)+" K]");
    500500
    501501  dbprint(printlevel-voice+2,"//transforming M to U...");
     
    505505  M=daU[2]*M*U;
    506506  dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, "
    507     +string(memory(1))+" bytes]");
     507    +string((memory(1)+1023)/1024)+" K]");
    508508
    509509  dbprint(printlevel-voice+2,
     
    541541  dbprint(printlevel-voice+2,
    542542    "//...integer differences of eigenvalues of M0 computed ["+string(timer-t)
    543     +" secs, "+string(memory(1))+" bytes]");
     543    +" secs, "+string((memory(1)+1023)/1024)+" K]");
    544544
    545545  dbprint(printlevel-voice+2,"//transforming M...");
     
    567567  }
    568568  dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, "
    569     +string(memory(1))+" bytes]");
     569    +string((memory(1)+1023)/1024)+" K]");
    570570
    571571  return(M);
     
    609609  e=pcvcv2p(quotV(V1+V2,N),0,N);
    610610  dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
    611     +string(memory(1))+" bytes]");
     611    +string((memory(1)+1023)/1024)+" K]");
    612612
    613613  dbprint(printlevel-voice+2,"//e=");
     
    638638    {
    639639      dbprint(printlevel-voice+2,"//...compared with previous lattice ["
    640         +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
     640        +string(timer-t)+" secs, "+string((memory(1)+1023)/1024)+" K]");
    641641
    642642      dbprint(printlevel-voice+2,"//increasing K and N...");
     
    656656        module(mdivp(M*U,var(1)^(kappa-1))),(kappa-1)*(mu-1)));
    657657      dbprint(printlevel-voice+2,"//...lattice enlarged ["+string(timer-t)
    658         +" secs, "+string(memory(1))+" bytes]");
     658        +" secs, "+string((memory(1)+1023)/1024)+" K]");
    659659
    660660      dbprint(printlevel-voice+2,"//comparing with previous lattice...");
     
    662662    }
    663663    dbprint(printlevel-voice+2,"//...compared with previous lattice ["
    664       +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
     664      +string(timer-t)+" secs, "+string((memory(1)+1023)/1024)+" K]");
    665665    dbprint(printlevel-voice+2,"//...t*nabla-stable lattice computed");
    666666
     
    673673      dbprint(printlevel-voice+2,
    674674        "//...C{f}-basis of t*nabla-stable lattice computed ["+string(timer-t)
    675         +" secs, "+string(memory(1))+" bytes]");
     675        +" secs, "+string((memory(1)+1023)/1024)+" K]");
    676676    }
    677677
     
    706706      leadcoef(daU[1])*var(1)^(kappa+pcvmindeg(daU[1])-1));
    707707    dbprint(printlevel-voice+2,"//...M/t^kappa transformed to simple pole ["
    708       +string(timer-t)+" secs, "+string(memory(1))+" bytes]");
     708      +string(timer-t)+" secs, "+string((memory(1)+1023)/1024)+" K]");
    709709  }
    710710
     
    718718    int delta=mid(eM0);
    719719    dbprint(printlevel-voice+2,"//...delta computed ["+string(timer-t)
    720       +" secs, "+string(memory(1))+" bytes]");
     720      +" secs, "+string((memory(1)+1023)/1024)+" K]");
    721721
    722722    dbprint(printlevel-voice+2,"//delta="+string(delta));
     
    750750        dbprint(printlevel-voice+2,
    751751          "//...M/t^kappa transformed to simple pole ["+string(timer-t)
    752           +" secs, "+string(memory(1))+" bytes]");
     752          +" secs, "+string((memory(1)+1023)/1024)+" K]");
    753753      }
    754754
     
    784784  ideal e=kbase(std(jacob(f)));
    785785  dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
    786     +string(memory(1))+" bytes]");
     786    +string((memory(1)+1023)/1024)+" K]");
    787787
    788788  dbprint(printlevel-voice+2,
     
    865865    int mu=milnor(f);
    866866    dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, "
    867       +string(memory(1))+" bytes]");
     867      +string((memory(1)+1023)/1024)+" K]");
    868868
    869869    dbprint(printlevel-voice+2,"//mu="+string(mu));
     
    953953    int mu=milnor(f);
    954954    dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, "
    955       +string(memory(1))+" bytes]");
     955      +string((memory(1)+1023)/1024)+" K]");
    956956
    957957    dbprint(printlevel-voice+2,"//mu="+string(mu));
     
    10051005      e=pcvcv2p(quotV(V1+V2,N),0,N);
    10061006      dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, "
    1007         +string(memory(1))+" bytes]");
     1007        +string((memory(1)+1023)/1024)+" K]");
    10081008
    10091009      dbprint(printlevel-voice+2,"//e=");
Note: See TracChangeset for help on using the changeset viewer.