In This Package:

MatrixManip.h

Go to the documentation of this file.

// $Id: MatrixManip.h,v 1.2 2007/07/10 07:25:53 cattanem Exp $
#ifndef LHCBMATH_MATRIXMANIP_H 
#define LHCBMATH_MATRIXMANIP_H 1

// Include files

#include "Math/SMatrix.h"

namespace Gaudi 
{
  
  namespace Math {

    using namespace ROOT::Math;

    template <typename M>  
    SMatrix<typename M::value_type, M::kRows, M::kRows, 
            MatRepSym<typename M::value_type, M::kRows> > Symmetrize(const M& rhs) {

      SMatrix<typename M::value_type, M::kRows, M::kRows, 
        MatRepSym<typename M::value_type, M::kRows> > result;
    
      for ( unsigned int i = 0; i < M::kRows; ++i ) {
        for ( unsigned int j = 0; j < M::kRows; ++j ) {
          if (i<=j) result(i,j) = rhs(i,j);
        }
      }
      return result;
    }

  }
}

#endif // LHCB_MATRIXMANIP_H

---