invertPosDefMatrix.h File Reference

#include "Math/SMatrix.h"
#include "Math/CholeskyDecomp.h"
#include <type_traits>

Go to the source code of this file.

Macros
#define	SMATRIX_USE_CONSTEXPR

Functions
template<typename PDM2 >
void	fastInvertPDM2 (PDM2 &mm)
template<typename T , unsigned int N>
bool	invertPosDefMatrix (ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > &m)
template<typename T , unsigned int N>
bool	invertPosDefMatrix (ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &mIn, ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > &mOut)
template<>
bool	invertPosDefMatrix< double, 1 > (ROOT::Math::SMatrix< double, 1, 1, ROOT::Math::MatRepSym< double, 1 > > &m)
template<>
bool	invertPosDefMatrix< double, 2 > (ROOT::Math::SMatrix< double, 2, 2, ROOT::Math::MatRepSym< double, 2 > > &m)
template<>
bool	invertPosDefMatrix< float, 1 > (ROOT::Math::SMatrix< float, 1, 1, ROOT::Math::MatRepSym< float, 1 > > &m)
template<>
bool	invertPosDefMatrix< float, 2 > (ROOT::Math::SMatrix< float, 2, 2, ROOT::Math::MatRepSym< float, 2 > > &m)

Macro Definition Documentation

SMATRIX_USE_CONSTEXPR

#define SMATRIX_USE_CONSTEXPR

Definition at line 4 of file invertPosDefMatrix.h.

Function Documentation

fastInvertPDM2()

template<typename PDM2 >

void fastInvertPDM2

(

PDM2 &

mm

)

Definition at line 20 of file invertPosDefMatrix.h.

                              {
  auto m = mm.Array();
 
  constexpr typename std::remove_reference<decltype(m[0])>::type one = 1.;
  auto c0 = one / m[0];
  auto c1 = m[1] * m[1] * c0;
  auto c2 = one / (m[2] - c1);
 
  auto li21 = c1 * c0 * c2;
  m[0] = li21 + c0;
  m[1] = -m[1] * c0 * c2;
  m[2] = c2;
}

References fftjetpileupestimator_calo_uncalib_cfi::c0, alignmentValidation::c1, math::cholesky::li21, visualization-live-secondInstance_cfg::m, and SiPixelPI::one.

Referenced by invertPosDefMatrix< double, 2 >(), and invertPosDefMatrix< float, 2 >().

invertPosDefMatrix() [1/2]

template<typename T , unsigned int N>

bool invertPosDefMatrix ( ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > &  m )

inline

Definition at line 10 of file invertPosDefMatrix.h.

                                                                                      {
  ROOT::Math::CholeskyDecomp<T, N> decomp(m);
  if (!decomp) {
    return m.Invert();
  } else
    decomp.Invert(m);
  return true;
}

References visualization-live-secondInstance_cfg::m.

Referenced by SiTrackerMultiRecHitUpdator::calcParameters(), PerigeeTrajectoryError::calculateWeightMatrix(), KalmanVertexUpdator< 5 >::chi2Increment(), TrajectoryStateCombiner::combine(), SiTrackerMultiRecHitUpdator::ComputeWeight(), MRHChi2MeasurementEstimator::estimate(), KalmanVertexUpdator< 5 >::positionUpdate(), KalmanVertexTrackUpdator< 5 >::trackRefit(), KinematicConstrainedVertexUpdatorT< nTrk, nConstraint >::update(), LocalTrajectoryError::weightMatrix(), and PosteriorWeightsCalculator::weights().

invertPosDefMatrix() [2/2]

template<typename T , unsigned int N>

bool invertPosDefMatrix ( ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > const &  mIn, ROOT::Math::SMatrix< T, N, N, ROOT::Math::MatRepSym< T, N > > &  mOut  )

inline

Definition at line 57 of file invertPosDefMatrix.h.

                                                                                         {
  ROOT::Math::CholeskyDecomp<T, N> decomp(mIn);
  if (!decomp) {
    mOut = mIn;
    return mOut.Invert();
  } else
    decomp.Invert(mOut);
  return true;
}

invertPosDefMatrix< double, 1 >()

template<>

bool invertPosDefMatrix< double, 1 > ( ROOT::Math::SMatrix< double, 1, 1, ROOT::Math::MatRepSym< double, 1 > > &  m )

inline

Definition at line 35 of file invertPosDefMatrix.h.

                                                                                                           {
  m(0, 0) = 1. / m(0, 0);
  return true;
}

References visualization-live-secondInstance_cfg::m.

invertPosDefMatrix< double, 2 >()

template<>

bool invertPosDefMatrix< double, 2 > ( ROOT::Math::SMatrix< double, 2, 2, ROOT::Math::MatRepSym< double, 2 > > &  m )

inline

Definition at line 46 of file invertPosDefMatrix.h.

                                                                                                           {
  fastInvertPDM2(m);
  return true;
}

References fastInvertPDM2(), and visualization-live-secondInstance_cfg::m.

invertPosDefMatrix< float, 1 >()

template<>

bool invertPosDefMatrix< float, 1 > ( ROOT::Math::SMatrix< float, 1, 1, ROOT::Math::MatRepSym< float, 1 > > &  m )

inline

Definition at line 40 of file invertPosDefMatrix.h.

                                                                                                        {
  m(0, 0) = 1.f / m(0, 0);
  return true;
}

References visualization-live-secondInstance_cfg::m.

invertPosDefMatrix< float, 2 >()

template<>

bool invertPosDefMatrix< float, 2 > ( ROOT::Math::SMatrix< float, 2, 2, ROOT::Math::MatRepSym< float, 2 > > &  m )

inline

Definition at line 51 of file invertPosDefMatrix.h.

                                                                                                        {
  fastInvertPDM2(m);
  return true;
}

References fastInvertPDM2(), and visualization-live-secondInstance_cfg::m.