Changeset 8942a5 in git for Singular/LIB/matrix.lib
- Timestamp:
- Dec 22, 2000, 3:33:13 PM (23 years ago)
- Branches:
- (u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
- Children:
- 8bb77bfee92151ac9ed599c9d789653dc8a3570c
- Parents:
- 803c5a17d620ffa41d14efaab791dcf17ad1f7a4
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/matrix.lib
r803c5a1 r8942a5 1 // $Id: matrix.lib,v 1.14 2000-12-19 14:37:26 anne Exp $2 1 // GMG/BM, last modified: 8.10.98 3 2 /////////////////////////////////////////////////////////////////////////////// 4 5 version="$Id: matrix.lib,v 1.14 2000-12-19 14:37:26 anne Exp $"; 3 version="$Id: matrix.lib,v 1.15 2000-12-22 14:15:05 greuel Exp $"; 6 4 category="Linear Algebra"; 7 5 info=" 8 LIBRARY: matrix.lib PROCEDURES FOR MATRIX OPERATIONS6 LIBRARY: matrix.lib Elementary Matrix Operations 9 7 10 8 PROCEDURES: … … 463 461 { 464 462 int i,j; 465 matrix C=B; 466 for( i=2; i<=nrows(A); i=i+1 ) { C=dsum(C,B); } 467 matrix D[nrows(C)][ncols(A)*nrows(B)]; 468 for( j=1; j<=nrows(B); j=j+1 ) 469 { 470 for( i=1; i<=nrows(A); i=i+1 ) 471 { 472 D[(i-1)*nrows(B)+j,(j-1)*ncols(A)+1..j*ncols(A)]=A[i,1..ncols(A)]; 473 } 474 } 475 return(concat(C,D)); 463 matrix C,D; 464 for( i=1; i<=nrows(A); i++ ) 465 { 466 C = A[i,1]*B; 467 for( j=2; j<=ncols(A); j++ ) 468 { 469 C = concat(C,A[i,j]*B); 470 } 471 D = concat(D,transpose(C)); 472 } 473 D = transpose(D); 474 return(submat(D,2..nrows(D),1..ncols(D))); 476 475 } 477 476 example
Note: See TracChangeset
for help on using the changeset viewer.