EvtMatrix< T > Class Template Reference

#include <EvtMatrix.hh>

List of all members.

Public Member Functions
	EvtMatrix ()
	~EvtMatrix ()
void	setRange (int range)
T &	operator() (int row, int col)
T *	operator[] (int row)
T	det ()
EvtMatrix *	min (int row, int col)
EvtMatrix *	inverse ()
std::string	dump ()
Private Attributes
T **	_mat
int	_range
Friends
template<class M>
EvtMatrix< M > *	operator * (const EvtMatrix< M > &left, const EvtMatrix< M > &right)

---

Detailed Description

template<class T>
class EvtMatrix< T >

Definition at line 22 of file EvtMatrix.hh.

---

Constructor & Destructor Documentation

template<class T>

EvtMatrix< T >::EvtMatrix

(

)

[inline]

Definition at line 28 of file EvtMatrix.hh.

: _range( 0 ) {};

template<class T>

EvtMatrix< T >::~EvtMatrix

(

)

Definition at line 73 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, and EvtMatrix< T >::_range.

{
  for( int row = 0; row < _range; row++ )
    delete[] _mat[ row ];
  delete[] _mat;
}

---

Member Function Documentation

template<class T>

T EvtMatrix< T >::det

(

)

Definition at line 97 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, EvtMatrix< T >::_range, EvtMatrix< T >::det(), EvtMatrix< T >::min(), and T.

Referenced by EvtMatrix< T >::det(), and EvtMatrix< T >::inverse().

{
  if ( _range == 1 )
    return _mat[ 0 ][ 0 ];

  // There's no need to define the range 2 determinant manually, but it may
  //    speed up the calculation.
  if ( _range == 2 )
    return _mat[ 0 ][ 0 ] * _mat[ 1 ][ 1 ] - _mat[ 0 ][ 1 ] * _mat[ 1 ][ 0 ];

  T sum = 0.;

  for ( int col = 0; col < _range; col++ )
    {
      EvtMatrix< T >* minor = min( 0, col );
      sum += std::pow( -1., col ) * _mat[ 0 ][ col ] * minor->det();
      delete minor;
    }

  return sum;
}

template<class T>

std::string EvtMatrix< T >::dump

(

)

Definition at line 81 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, and EvtMatrix< T >::_range.

{
  std::ostringstream str;

  for ( int row = 0; row < _range; row++ )
    {
      str << "|";
      for ( int col = 0; col < _range; col++ )
        str << "\t" << _mat[ row ][ col ];
      str << "\t|" << std::endl;
    }

  return str.str();
}

template<class T>

EvtMatrix< T > * EvtMatrix< T >::inverse

(

)

Definition at line 140 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, EvtMatrix< T >::_range, EvtMatrix< T >::det(), EvtMatrix< T >::min(), EvtMatrix< T >::setRange(), and T.

Referenced by EvtDalitzReso::Fvector().

{
  EvtMatrix< T >* inv = new EvtMatrix< T >();
  inv->setRange( _range );

  if ( det() == 0 )
    {
      std::cerr << "This matrix has a null determinant and cannot be inverted. Returning zero matrix." << std::endl;
      for ( int row = 0; row < _range; row++ )
        for ( int col = 0; col < _range; col++ )
          (*inv)( row, col ) = 0.;
      return inv;
    }

  T determinant = det();

  for ( int row = 0; row < _range; row++ )
    for ( int col = 0; col < _range; col++ )
      {
        EvtMatrix< T >* minor = min( row, col );
        inv->_mat[col][row] = std::pow( -1., row + col ) * minor->det() / determinant;
        delete minor;
      }

  return inv;
}

template<class T>

EvtMatrix< T > * EvtMatrix< T >::min	(	int	row,
		int	col
	)

Definition at line 121 of file EvtMatrix.hh.

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

Referenced by EvtMatrix< T >::det(), and EvtMatrix< T >::inverse().

{
  EvtMatrix< T >* minor = new EvtMatrix< T >();
  minor->setRange( _range - 1 );

  int minIndex = 0;

  for ( int r = 0; r < _range; r++ )
    for ( int c = 0; c < _range; c++ )
      if ( ( r != row ) && ( c != col ) )
        {
          (*minor)( minIndex / ( _range - 1 ), minIndex % ( _range - 1 ) ) = _mat[ r ][ c ];
          minIndex++;
        }

  return minor;
}

template<class T>

T& EvtMatrix< T >::operator()	(	int	row,
		int	col
	)			[inline]

Definition at line 32 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat.

{ return _mat[ row ][ col ]; }

template<class T>

T* EvtMatrix< T >::operator[]

(

int

row

)

[inline]

Definition at line 33 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat.

{ return _mat[ row ];        }

template<class T>

void EvtMatrix< T >::setRange

(

int

range

)

[inline]

Definition at line 45 of file EvtMatrix.hh.

References EvtMatrix< T >::_mat, EvtMatrix< T >::_range, and T.

Referenced by EvtDalitzReso::Fvector(), EvtMatrix< T >::inverse(), EvtMatrix< T >::min(), and operator *().

{
  // If the range is changed, delete any previous matrix stored
  //    and allocate elements with the newly specified range.
  if ( _range != range )
    {
      if ( _range )
        {
          for ( int row = 0; row < _range; row++ )
            delete[] _mat[ row ];
          delete[] _mat;
        }

      _mat = new T*[ range ];
      for ( int row = 0; row < range; row++ )
        _mat[ row ] = new T[ range ];

      // Set the new range.
      _range = range;
    }

  // Since user is willing to change the range, reset the matrix elements.
  for ( int row = 0; row < _range; row++ )
    for ( int col = 0; col < _range; col++ )
      _mat[ row ][ col ] = 0.;
}

---

Friends And Related Function Documentation

template<class T>

template<class M>

EvtMatrix< M >* operator *	(	const EvtMatrix< M > &	left,
		const EvtMatrix< M > &	right
	)			[friend]

---

Member Data Documentation

template<class T>

T** EvtMatrix< T >::_mat [private]

Definition at line 25 of file EvtMatrix.hh.

Referenced by EvtMatrix< T >::det(), EvtMatrix< T >::dump(), EvtMatrix< T >::inverse(), EvtMatrix< T >::min(), operator *(), EvtMatrix< T >::operator()(), EvtMatrix< T >::operator[](), EvtMatrix< T >::setRange(), and EvtMatrix< T >::~EvtMatrix().

template<class T>

int EvtMatrix< T >::_range [private]

Definition at line 26 of file EvtMatrix.hh.

Referenced by EvtMatrix< T >::det(), EvtMatrix< T >::dump(), EvtMatrix< T >::inverse(), EvtMatrix< T >::min(), operator *(), EvtMatrix< T >::setRange(), and EvtMatrix< T >::~EvtMatrix().

---