In This Package:

MatrixInversion.h

Go to the documentation of this file.

#ifndef LHCBMATH_MATRIXINVERSION_H
#define LHCBMATH_MATRIXINVERSION_H

// ============================================================================
// ROOT matrix types
// ============================================================================
#include "GaudiKernel/SymmetricMatrixTypes.h"

// ============================================================================
// matrix inversion implementation
// ============================================================================
#include "SymMatrixInverter.h"
#include "SymPosDefMatrixInverter.h"

// ============================================================================
namespace Gaudi
{
  namespace Math 
  {
    template<class T> inline bool invertPosDefSymMatrix(T& matrix)
    {
      enum { N = T::kRows };
      typedef typename T::value_type F;
      typedef typename ROOT::Math::MatRepSym<F, N> R;
      typedef typename ROOT::Math::SMatrix<F, N, N, R> M;
      return Gaudi::Math::SymPosDefMatrixInverter::inverter<M, F, N>()(matrix);
    }

    template<class T> inline bool invertSymMatrix(T& matrix)
    {
      enum { N = T::kRows };
      typedef typename T::value_type F;
      typedef typename ROOT::Math::MatRepSym<F, N> R;
      typedef typename ROOT::Math::SMatrix<F, N, N, R> M;
      return Gaudi::Math::SymMatrixInverter::inverter<M, F, N>()(matrix);
    }
  }
}

#endif // LHCBMATH_MATRIXINVERSION_H

// vim:sw=2:tw=78:ft=cpp

---