Changeset 97a7b44 in git
- Timestamp:
- Jun 22, 1999, 5:47:35 PM (24 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- 3e9ae9c78c81633291fa286789448020d1984470
- Parents:
- 4ad2bb592a6ffe7f195320a5d4bec4f9ef0c57dd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/mondromy.lib
r4ad2bb r97a7b44 1 1 /////////////////////////////////////////////////////////////////////////////// 2 2 3 version="$Id: mondromy.lib,v 1. 3 1999-06-22 15:38:54mschulze Exp $";3 version="$Id: mondromy.lib,v 1.4 1999-06-22 15:47:35 mschulze Exp $"; 4 4 info=" 5 5 LIBRARY: mondromy.lib PROCEDURES TO COMPUTE THE MONODROMY OF A SINGULARITY … … 164 164 poly v=lift(fetch(br,u),1)[1,1]; 165 165 dbprint(printlevel-voice+2,"//...inverse computed ["+string(timer-t)+ 166 " secs, "+string( memory(1))+" bytes]");166 " secs, "+string((memory(1)+1023)/1024)+" K]"); 167 167 168 168 setring br; … … 200 200 poly detU=det(U); 201 201 dbprint(printlevel-voice+2,"//...determinant computed ["+string(timer-t)+ 202 " secs, "+string( memory(1))+" bytes]");202 " secs, "+string((memory(1)+1023)/1024)+" K]"); 203 203 204 204 if(detU==0) … … 218 218 matrix adjU=lift(U,detU*freemodule(nrows(U))); 219 219 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]"); 221 221 222 222 setring br; … … 268 268 dbprint(printlevel-voice+2,"//kappa="+string(kappa)); 269 269 dbprint(printlevel-voice+2,"//...kappa computed ["+string(timer-t)+" secs, " 270 +string( memory(1))+" bytes]");270 +string((memory(1)+1023)/1024)+" K]"); 271 271 272 272 dbprint(printlevel-voice+2,"//computing xi..."); … … 274 274 vector xi=lift(jf,fkappa)[1]; 275 275 dbprint(printlevel-voice+2,"//...xi computed ["+string(timer-t)+" secs, " 276 +string( memory(1))+" bytes]");276 +string((memory(1)+1023)/1024)+" K]"); 277 277 278 278 dbprint(printlevel-voice+2,"//computing u..."); … … 280 280 poly u=(matrix(jf)*xi)[1,1]/fkappa; 281 281 dbprint(printlevel-voice+2,"//...u computed ["+string(timer-t)+" secs, " 282 +string( memory(1))+" bytes]");282 +string((memory(1)+1023)/1024)+" K]"); 283 283 284 284 return(list(kappa,xi,u)); … … 365 365 { 366 366 dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t) 367 +" secs, "+string( memory(1))+" bytes]");367 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 368 368 369 369 deltaP1=getdeltaP1(f,K,N,deltaN); … … 386 386 } 387 387 dbprint(printlevel-voice+2,"//...codimension computed ["+string(timer-t) 388 +" secs, "+string( memory(1))+" bytes]");388 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 389 389 390 390 return(K,N,P1,P2,Pe,V1,V2,Ve); … … 426 426 Vnablae=nablaK(f,kappa,xi,u,N,prevN,Vnablae,e); 427 427 dbprint(printlevel-voice+2,"//...nabla(e) computed ["+string(timer-t) 428 +" secs, "+string( memory(1))+" bytes]");428 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 429 429 430 430 dbprint(printlevel-voice+2, … … 437 437 matrix C=lift(W,module(Vnablae[1..size(Vnablae)])); 438 438 dbprint(printlevel-voice+2,"//...nabla(e) lifted ["+string(timer-t) 439 +" secs, "+string( memory(1))+" bytes]");439 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 440 440 441 441 dbprint(printlevel-voice+2,"//computing e-lift of nabla(e)..."); … … 454 454 } 455 455 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]"); 457 457 458 458 return(M,N,Vnablae); … … 497 497 } 498 498 dbprint(printlevel-voice+2,"//...U computed ["+string(timer-t)+" secs, " 499 +string( memory(1))+" bytes]");499 +string((memory(1)+1023)/1024)+" K]"); 500 500 501 501 dbprint(printlevel-voice+2,"//transforming M to U..."); … … 505 505 M=daU[2]*M*U; 506 506 dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, " 507 +string( memory(1))+" bytes]");507 +string((memory(1)+1023)/1024)+" K]"); 508 508 509 509 dbprint(printlevel-voice+2, … … 541 541 dbprint(printlevel-voice+2, 542 542 "//...integer differences of eigenvalues of M0 computed ["+string(timer-t) 543 +" secs, "+string( memory(1))+" bytes]");543 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 544 544 545 545 dbprint(printlevel-voice+2,"//transforming M..."); … … 567 567 } 568 568 dbprint(printlevel-voice+2,"//...M transformed ["+string(timer-t)+" secs, " 569 +string( memory(1))+" bytes]");569 +string((memory(1)+1023)/1024)+" K]"); 570 570 571 571 return(M); … … 609 609 e=pcvcv2p(quotV(V1+V2,N),0,N); 610 610 dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, " 611 +string( memory(1))+" bytes]");611 +string((memory(1)+1023)/1024)+" K]"); 612 612 613 613 dbprint(printlevel-voice+2,"//e="); … … 638 638 { 639 639 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]"); 641 641 642 642 dbprint(printlevel-voice+2,"//increasing K and N..."); … … 656 656 module(mdivp(M*U,var(1)^(kappa-1))),(kappa-1)*(mu-1))); 657 657 dbprint(printlevel-voice+2,"//...lattice enlarged ["+string(timer-t) 658 +" secs, "+string( memory(1))+" bytes]");658 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 659 659 660 660 dbprint(printlevel-voice+2,"//comparing with previous lattice..."); … … 662 662 } 663 663 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]"); 665 665 dbprint(printlevel-voice+2,"//...t*nabla-stable lattice computed"); 666 666 … … 673 673 dbprint(printlevel-voice+2, 674 674 "//...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]"); 676 676 } 677 677 … … 706 706 leadcoef(daU[1])*var(1)^(kappa+pcvmindeg(daU[1])-1)); 707 707 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]"); 709 709 } 710 710 … … 718 718 int delta=mid(eM0); 719 719 dbprint(printlevel-voice+2,"//...delta computed ["+string(timer-t) 720 +" secs, "+string( memory(1))+" bytes]");720 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 721 721 722 722 dbprint(printlevel-voice+2,"//delta="+string(delta)); … … 750 750 dbprint(printlevel-voice+2, 751 751 "//...M/t^kappa transformed to simple pole ["+string(timer-t) 752 +" secs, "+string( memory(1))+" bytes]");752 +" secs, "+string((memory(1)+1023)/1024)+" K]"); 753 753 } 754 754 … … 784 784 ideal e=kbase(std(jacob(f))); 785 785 dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, " 786 +string( memory(1))+" bytes]");786 +string((memory(1)+1023)/1024)+" K]"); 787 787 788 788 dbprint(printlevel-voice+2, … … 865 865 int mu=milnor(f); 866 866 dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, " 867 +string( memory(1))+" bytes]");867 +string((memory(1)+1023)/1024)+" K]"); 868 868 869 869 dbprint(printlevel-voice+2,"//mu="+string(mu)); … … 953 953 int mu=milnor(f); 954 954 dbprint(printlevel-voice+2,"//...mu computed ["+string(timer-t)+" secs, " 955 +string( memory(1))+" bytes]");955 +string((memory(1)+1023)/1024)+" K]"); 956 956 957 957 dbprint(printlevel-voice+2,"//mu="+string(mu)); … … 1005 1005 e=pcvcv2p(quotV(V1+V2,N),0,N); 1006 1006 dbprint(printlevel-voice+2,"//...e computed ["+string(timer-t)+" secs, " 1007 +string( memory(1))+" bytes]");1007 +string((memory(1)+1023)/1024)+" K]"); 1008 1008 1009 1009 dbprint(printlevel-voice+2,"//e=");
Note: See TracChangeset
for help on using the changeset viewer.