Changeset 8942a5 in git for Singular/LIB/matrix.lib


Ignore:
Timestamp:
Dec 22, 2000, 3:33:13 PM (23 years ago)
Author:
Gert-Martin Greuel <greuel@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
8bb77bfee92151ac9ed599c9d789653dc8a3570c
Parents:
803c5a17d620ffa41d14efaab791dcf17ad1f7a4
Message:
* GMG: Kosmetik


git-svn-id: file:///usr/local/Singular/svn/trunk@4982 2c84dea3-7e68-4137-9b89-c4e89433aadc
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 $
    21// GMG/BM, last modified: 8.10.98
    32///////////////////////////////////////////////////////////////////////////////
    4 
    5 version="$Id: matrix.lib,v 1.14 2000-12-19 14:37:26 anne Exp $";
     3version="$Id: matrix.lib,v 1.15 2000-12-22 14:15:05 greuel Exp $";
    64category="Linear Algebra";
    75info="
    8 LIBRARY:  matrix.lib    PROCEDURES FOR MATRIX OPERATIONS
     6LIBRARY:  matrix.lib    Elementary Matrix Operations
    97
    108PROCEDURES:
     
    463461{
    464462   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)));
    476475}
    477476example
Note: See TracChangeset for help on using the changeset viewer.