// Matrix.cpp: implementation of the Matrix class.
//
//////////////////////////////////////////////////////////////////////
#include <iostream.h>
#include <iomanip.h>
#include "Matrix.h"
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
Matrix::Matrix()
{
m = 0;
n = 0;
elem = NULL;
}
Matrix::Matrix( int rSize, int cSize, float init )
{
m = rSize;
n = cSize;
elem = new float[ m * n ];
for ( int i = 0; i < m * n; i++ )
{
elem[i] = init;
}
}
Matrix::Matrix( Matrix & A )
{
m = A.m;
n = A.n;
elem = new float[ m * n ];
for ( int i = 0; i < m * n; i++ )
{
elem[i] = A.elem[i];
}
}
Matrix::~Matrix()
{
if ( !IsEmpty() )
{
delete [] elem;
}
}
bool Matrix::IsEmpty()
{
if ( elem == NULL )
{
return true;
}
else
{
return false;
}
}
float & Matrix::operator ()( int mid,int nid )
{
if ( IsEmpty() || mid >= m || nid >= n)
{
elemEmpty = 0;
return elemEmpty;
}
else
{
return elem[ mid * m + nid ];
}
}
ostream & operator << ( ostream & out, Matrix & A )
{
int width = out.width();
for ( int i = 0; i < A.m; i++ )
{
for ( int j = 0; j < A.n; j++ )
{
out << setw( width ) << A( i, j );
}
out << endl;
}
return out;
}
Matrix operator + ( Matrix & A, Matrix & B )
{
Matrix result( A.m, A.n );
for ( int i = 0; i < result.m; i++ )
{
for ( int j = 0; j < result.n; j++ )
{
result( i, j ) = A( i, j ) + B( i, j );
}
}
return result;
}
Matrix operator - ( Matrix & A, Matrix & B )
{
Matrix result( A.m, A.n );
for ( int i = 0; i < result.m; i++ )
{
for ( int j = 0; j < result.n; j++ )
{
result( i, j ) = A( i, j ) - B( i, j );
}
}
return result;
}
Matrix operator * ( Matrix & A, Matrix & B )
{
Matrix result( A.m, B.n );
int i, j, k;
int maxN = A.n > B.m ? A.n : B.m;
for ( i = 0; i < result.m; i++ )
{
for ( j = 0; j < result.n; j++ )
{
for( k = 0; k < maxN; k++ )
{
result( i, j ) += A( i, k ) * B( k, j );
}
}
}
return result;
}
Matrix & Matrix::operator = ( Matrix & A )
{
m = A.m;
n = A.n;
if ( !IsEmpty() )
{
delete [] elem;
}
elem = new float[ m * n ];
for ( int i = 0; i < m * n; i++ )
{
elem[i] = A.elem[i];
}
return *this;
}
Matrix Matrix::Transpose()
{
Matrix result( n, m );
for( int i = 0; i < n; i++ )
{
for ( int j = 0; j < m; j++ )
{
result( i, j ) = (*this)( j, i );
}
}
return result;
}
这个是Matrix.cpp中的实现细节
