EvtMatrix.hh File Reference

#include <vector>
#include <sstream>

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 169 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;
}