Changeset 7af4ad9 in git
- Timestamp:
- Jan 7, 2001, 4:21:46 PM (23 years ago)
- Branches:
- (u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
- Children:
- 1adaa7cccf071b7264e7d562e31ff2fbcf28bda0
- Parents:
- 90611c96f29bf20282ba8567dc61e6aa0533155e
- Location:
- Singular/LIB
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/gaussman.lib
r90611c r7af4ad9 1 1 /////////////////////////////////////////////////////////////////////////////// 2 version="$Id: gaussman.lib,v 1.2 0 2000-12-24 15:39:11 mschulzeExp $";2 version="$Id: gaussman.lib,v 1.21 2001-01-07 15:21:09 greuel Exp $"; 3 3 category="Singularities"; 4 4 … … 147 147 RETURN: 148 148 @format 149 if mode=0:150 matrix M: exp(-2*pi*i*M) is a monodromy matrix of f151 if mode=1:152 ideal e: exp(-2*pi*i*e) is the spectrum of the monodromy of f149 if mode=0: 150 matrix M: exp(-2*pi*i*M) is a monodromy matrix of f 151 if mode=1: 152 ideal e: exp(-2*pi*i*e) is the spectrum of the monodromy of f 153 153 @end format 154 154 SEE ALSO: monodromy.lib, jordan.lib … … 371 371 372 372 proc vfiltration(poly f,list #) 373 "USAGE: vfiltration(f[,mode]); poly f, int mode[ =1]374 ASSUME: basering has local ordering, f hasisolated singularity at 0375 RETURN: 376 @format 377 list l:378 if mode=0 or mode=1:379 ideal l[1]:spectral numbers in increasing order380 intvec l[2]:381 int l[2][i]:multiplicity of spectral number l[1][i]382 if mode=1:383 list l[3]:384 module l[3][i]: vector space basis of l[1][i]-th graded part385 of the V-filtration on H''/H' in terms of l[4]386 ideal l[4]:monomial vector space basis of H''/H'387 ideal l[5]: standard basis of theJacobian ideal388 @end format389 NOTE: H' and H'' denote theBrieskorn lattices373 "USAGE: vfiltration(f[,mode]); poly f, int mode[default=1] 374 ASSUME: local ordering, f isolated singularity at 0 375 RETURN: list l: 376 @format 377 list l: 378 if mode=0 or mode=1: 379 l[1]: ideal, spectral numbers in increasing order 380 l[2]: intvec 381 l[2][i]: int, multiplicity of spectral number l[1][i] 382 if mode=1 : 383 l[3]: list 384 l[3][i]: module, vector space basis of l[1][i]-th graded 385 part of the V-filtration on H''/H' in terms of l[4] 386 l[4]: ideal, monomial vector space basis of H''/H' 387 l[5]: ideal, standard basis of Jacobian ideal 388 @end format 389 NOTE: H' and H'' denote Brieskorn lattices 390 390 SEE ALSO: spectrum.lib 391 391 KEYWORDS: singularities; Gauss-Manin connection; spectrum; … … 749 749 750 750 proc vfiltjacalg(list l) 751 "USAGE: vfiltjacalg(vfiltration(f)); poly f752 ASSUME: basering has local ordering, f hasisolated singularity at 0753 RETURN: 754 @format 755 list l:756 ideal l[1]: spectral numbers of the V-filtration757 on theJacobian algebra in increasing order758 intvec l[2]:759 int l[2][i]:multiplicity of spectral number l[1][i]760 list l[3]:761 module l[3][i]: vector space basis of thel[1][i]-th graded part762 of the V-filtration on the Jacobian algebra763 in termsof l[4]764 ideal l[4]:monomial vector space basis of the Jacobian algebra765 ideal l[5]: standard basis of theJacobian ideal751 "USAGE: vfiltjacalg(vfiltration(f)); 752 ASSUME: local ordering, f isolated singularity at 0 753 RETURN: 754 @format 755 list l: 756 l[1]: ideal, spectral numbers of the V-filtration on the 757 Jacobian algebra in increasing order 758 l[2]: intvec 759 l[2][i]: int, multiplicity of spectral number l[1][i] 760 l[3]: list 761 l[3][i]: module, vector space basis of l[1][i]-th graded part 762 of the V-filtration on the Jacobian algebra in terms 763 of l[4] 764 l[4]: ideal, monomial vector space basis of the Jacobian algebra 765 l[5]: ideal, standard basis of Jacobian ideal 766 766 @end format 767 767 EXAMPLE: example vfiltjacalg; shows an example -
Singular/LIB/matrix.lib
r90611c r7af4ad9 1 1 // GMG/BM, last modified: 8.10.98 2 2 /////////////////////////////////////////////////////////////////////////////// 3 version="$Id: matrix.lib,v 1. 19 2001-01-06 00:41:54greuel Exp $";3 version="$Id: matrix.lib,v 1.20 2001-01-07 15:21:46 greuel Exp $"; 4 4 category="Linear Algebra"; 5 5 info=" … … 503 503 504 504 proc gauss_col (matrix A, list #) 505 "USAGE: gauss_col(A[, e]); A matrix, eany type505 "USAGE: gauss_col(A[,any]); A a matrix, any any type 506 506 RETURN: - a matrix B, if called with one argument; B is the complete column- 507 507 reduced upper-triangular normal form of A if A is constant, … … 509 509 no division by polynomials). 510 510 @* - a list L of two matrices, if called with two arguments; 511 L sati dfies L[1] = A * L[2] with L[1] the column-reduced form of A511 L satisfies L[1] = A * L[2] with L[1] the column-reduced form of A 512 512 and L[2] the transformation matrix. 513 513 NOTE: * The procedure just applies interred to A with ordering (C,dp). … … 518 518 @* * Parameters are allowed. Hence, if the entries of A are parameters, 519 519 B is the column-reduced form of A over the rational function field. 520 SEE ALSO: colred 520 521 EXAMPLE: example gauss_col; shows an example 521 522 " … … 526 527 module M = A; 527 528 intvec v = option(get); 529 //------------------------ change ordering if necessary ---------------------- 528 530 if( ordstr(R) != "C,dp(nvars(R))" ) 529 531 { … … 532 534 module M = A; 533 535 } 536 //------------------------------ start computation --------------------------- 534 537 option(redSB); 535 538 M = simplify(interred(M),1); 536 /* 537 if (#[1] == 0) 538 { 539 M = simplify(std(M),1); 540 } 541 if (#[1] == 1) 542 { 543 M = simplify(interred(M),1); 544 } 545 */ 546 if (size(#) != 0) 547 { 548 module N = lift(A,M); 549 } 539 if(size(#) != 0) 540 { 541 module N = lift(A,M); 542 } 543 //--------------- reset ring and options and return -------------------------- 550 544 if ( u==1 ) 551 545 { … … 596 590 597 591 proc gauss_row (matrix A, list #) 598 "USAGE: gauss_row(A ); A matrix, e any type592 "USAGE: gauss_row(A [,e]); A matrix, e any type 599 593 RETURN: - a matrix B, if called with one argument; B is the complete row- 600 594 reduced lower-triangular normal form of A if A is constant, … … 611 605 @* * Parameters are allowed. Hence, if the entries of A are parameters, 612 606 B is the row-reduced form of A over the rational function field. 607 SEE ALSO: rowred 613 608 EXAMPLE: example gauss_row; shows an example 614 609 " … … 633 628 0, 0, 0, 0, 0, 0, 0, 2a;""; 634 629 print(gauss_row(m));""; 635 636 630 ring S=0,x,dp; 637 631 matrix A[4][5] = 3, 1,1,-1,2, … … 766 760 the row-reduced form of A and L[2] the transformation matrix 767 761 (if rowred is called with two arguments). 768 NOTE: * This procedure is designed for teaching purposes only.762 NOTE: * This procedure is designed for teaching purposes mainly. 769 763 @* * The straight forward Gaussian algorithm is implemented in the 770 library (no standard basis computation) , hence it is not very771 efficient. The transformation matrix is obtained by concatenating772 a unit matrix to A, proc gauss_redshould be faster.764 library (no standard basis computation). 765 The transformation matrix is obtained by concatenating a unit 766 matrix to A. proc gauss_row should be faster. 773 767 @* * It should only be used with exact coefficient field (there is no 774 768 pivoting) over the polynomial ring (ordering lp or dp). 775 769 @* * Parameters are allowed. Hence, if the entries of A are parameters 776 770 the computation takes place over the field of rational functions. 771 SEE ALSO: gauss_row 777 772 EXAMPLE: example rowred; shows an example 778 773 " … … 850 845 the column-reduced form of A and L[2] the transformation matrix 851 846 (if colred is called with two arguments). 852 NOTE: * This procedure is designed for teaching purposes only.853 @* * It applies rowred to the transposed matrix , hence it is not very854 efficient (cf. rowred);proc gauss_col should be faster.847 NOTE: * This procedure is designed for teaching purposes mainly. 848 @* * It applies rowred to the transposed matrix. 849 proc gauss_col should be faster. 855 850 @* * It should only be used with exact coefficient field (there is no 856 851 pivoting) over the polynomial ring (ordering lp or dp). 857 852 @* * Parameters are allowed. Hence, if the entries of A are parameters 858 853 the computation takes place over the field of rational functions. 854 SEE ALSO: gauss_col 859 855 EXAMPLE: example colred; shows an example 860 856 "
Note: See TracChangeset
for help on using the changeset viewer.