/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Generator/BesEvtGen/BesEvtGen-00-03-58/src/EvtGen/EvtGenBase/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

template<class T>

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

Definition at line 169 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, EvtMatrix< T >::_range, genRecEmupikp::line, and EvtMatrix< T >::setRange().

{
  // 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;
}

---