EvtMatrix< T > Class Template Reference

#include <EvtMatrix.hh>

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 ()

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 20 of file EvtMatrix.hh.

Constructor & Destructor Documentation

EvtMatrix()

template<class T>

EvtMatrix< T >::EvtMatrix ( )

inline

Definition at line 26 of file EvtMatrix.hh.

: _range( 0 ){};

Referenced by det(), inverse(), min(), and operator*.

~EvtMatrix()

template<class T>

EvtMatrix< T >::~EvtMatrix

(

)

Definition at line 64 of file EvtMatrix.hh.

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

Member Function Documentation

det()

template<class T>

T EvtMatrix< T >::det

(

)

Definition at line 82 of file EvtMatrix.hh.

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

Referenced by det(), and inverse().

dump()

template<class T>

std::string EvtMatrix< T >::dump

(

)

Definition at line 69 of file EvtMatrix.hh.

                                                {
  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();
}

inverse()

template<class T>

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

(

)

Definition at line 119 of file EvtMatrix.hh.

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

min()

template<class T>

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

Definition at line 102 of file EvtMatrix.hh.

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

operator()()

template<class T>

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

inline

Definition at line 30 of file EvtMatrix.hh.

{ return _mat[row][col]; }

operator[]()

template<class T>

T * EvtMatrix< T >::operator[] ( int row )

inline

Definition at line 31 of file EvtMatrix.hh.

{ return _mat[row]; }

setRange()

template<class T>

void EvtMatrix< T >::setRange ( int range )

inline

Definition at line 41 of file EvtMatrix.hh.

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

Referenced by inverse(), min(), and operator*().

Friends And Related Symbol Documentation

operator*

template<class T>

template<class M>

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

friend

---

The documentation for this class was generated from the following file:

• EvtMatrix.hh