EvtMatrix.hh File Reference

#include <sstream>
#include <vector>

Go to the source code of this file.

Classes
class	EvtMatrix< T >

Functions
template<class T>
EvtMatrix< T > *	operator* (const EvtMatrix< T > &left, const EvtMatrix< T > &right)

Function Documentation

operator*()

template<class T>

EvtMatrix< T > * operator*	(	const EvtMatrix< T > &	left,
		const EvtMatrix< T > &	right )

Definition at line 147 of file EvtMatrix.hh.

                                                                               {
  // Chech that the matrices have the correct range.
  if ( left._range != right._range )
  {
    std::cerr << "These matrices cannot be multiplied." << std::endl;
    return new EvtMatrix<T>();
  }
 
  EvtMatrix<T>* mat = new EvtMatrix<T>();
  mat->setRange( left._range );
 
  // Initialize the elements of the matrix.
  for ( int row = 0; row < left._range; row++ )
    for ( int col = 0; col < right._range; col++ ) ( *mat )[row][col] = 0;
 
  for ( int row = 0; row < left._range; row++ )
    for ( int col = 0; col < right._range; col++ )
      for ( int line = 0; line < right._range; line++ )
        ( *mat )[row][col] += left._mat[row][line] * right._mat[line][col];
 
  return mat;
}