Changeset cb1e43b in git
- Timestamp:
- Jan 3, 2001, 2:28:48 AM (23 years ago)
- Branches:
- (u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')
- Children:
- f875fae3f8c98402b705e50855bccd0ce2c50d90
- Parents:
- 9cc5d37df563ce9febe8a449f39cf9986d705e1e
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/matrix.lib
r9cc5d37 rcb1e43b 1 1 // GMG/BM, last modified: 8.10.98 2 2 /////////////////////////////////////////////////////////////////////////////// 3 version="$Id: matrix.lib,v 1.1 7 2001-01-02 13:10:59greuel Exp $";3 version="$Id: matrix.lib,v 1.18 2001-01-03 01:28:48 greuel Exp $"; 4 4 category="Linear Algebra"; 5 5 info=" … … 505 505 "USAGE: gauss_col(A[,e]); A matrix, e any type 506 506 RETURN: - a matrix B, if called with one argument; B is the complete column- 507 reduced upper-triangular normal form of A if A is constant, resp.508 as far as this is possible if A is a polynomial matrix509 (no division by polynomials).507 reduced upper-triangular normal form of A if A is constant, 508 (resp. as far as this is possible if A is a polynomial matrix; 509 no division by polynomials). 510 510 @* - a list L of two matrices, if called with two arguments; 511 511 L satidfies L[1] = A * L[2] with L[1] the column-reduced form of A … … 534 534 option(redSB); 535 535 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 /* 536 546 if (size(#) != 0) 537 547 { … … 588 598 "USAGE: gauss_row(A); A matrix, e any type 589 599 RETURN: - a matrix B, if called with one argument; B is the complete row- 590 reduced lower-triangular normal form of A if A is constant, resp.591 as far as this is possible if A is a polynomial matrix592 (no division by polynomials).600 reduced lower-triangular normal form of A if A is constant, 601 (resp. as far as this is possible if A is a polynomial matrix; 602 no division by polynomials). 593 603 @* - a list L of two matrices, if called with two arguments; 594 604 L satidfies L[1] = L[2] * A with L[1] the row-reduced form of A … … 749 759 proc rowred (matrix A,list #) 750 760 "USAGE: rowred(A[,e]); A matrix, e any type 751 RETURN: - a matrix B, being the row reduced form of A, as far as this 752 is possible over the polynomial ring (no division by polynomials) 753 (if rowred is called with one argument). 761 RETURN: - a matrix B, being the row reduced form of A, if rowred is called 762 with one argument. 763 (as far as this is possible over the polynomial ring; no division 764 by polynomials) 754 765 @* - a list L of two matrices, such that L[1] = L[2] * A with L[1] 755 766 the row-reduced form of A and L[2] the transformation matrix … … 759 770 library (no standard basis computation), hence it is not very 760 771 efficient. The transformation matrix is obtained by concatenating 761 a unit matrix to A .772 a unit matrix to A, proc gauss_red should be faster. 762 773 @* * It should only be used with exact coefficient field (there is no 763 774 pivoting) over the polynomial ring (ordering lp or dp). … … 832 843 proc colred (matrix A,list #) 833 844 "USAGE: colred(A[,e]); A matrix, e any type 834 RETURN: - a matrix B, being the column reduced form of A, as far as this 835 is possible over the polynomial ring (no division by polynomials) 836 (if colred is called with one argument). 845 RETURN: - a matrix B, being the column reduced form of A, if colred is 846 called with one argument. 847 (as far as this is possible over the polynomial ring; 848 no division by polynomials) 837 849 @* - a list L of two matrices, such that L[1] = A * L[2] with L[1] 838 850 the column-reduced form of A and L[2] the transformation matrix … … 840 852 NOTE: * This procedure is designed for teaching purposes only. 841 853 @* * It applies rowred to the transposed matrix, hence it is not very 842 efficient (cf. rowred) .854 efficient (cf. rowred); proc gauss_col should be faster. 843 855 @* * It should only be used with exact coefficient field (there is no 844 856 pivoting) over the polynomial ring (ordering lp or dp).
Note: See TracChangeset
for help on using the changeset viewer.