RungeKuttaUtils.cxx File Reference

#include "TrkExUtils/RungeKuttaUtils.h"
#include "TrkPatternParameters/PatternTrackParameters.h"
#include "TrkSurfaces/ConeSurface.h"
#include "TrkSurfaces/CylinderSurface.h"
#include "TrkSurfaces/DiscSurface.h"
#include "TrkSurfaces/DistanceSolution.h"
#include "TrkSurfaces/PerigeeSurface.h"
#include "TrkSurfaces/PlaneSurface.h"
#include "TrkSurfaces/StraightLineSurface.h"
#include "CxxUtils/vec.h"
#include "CxxUtils/inline_hints.h"

Go to the source code of this file.

Functions
ATH_FLATTEN	AmgSymMatrix (5) Trk

Function Documentation

AmgSymMatrix()

ATH_FLATTEN AmgSymMatrix

(

5

)

Definition at line 1391 of file RungeKuttaUtils.cxx.

{
  AmgSymMatrix(5) nM;
  AmgSymMatrix(5)& m = nM;
 
  Eigen::Map<const AmgVector(5)> JacMap0(&J[0], 5, 1);
  const AmgVector(5) a1 = M * JacMap0; //(5x5 * 5x1)
  m(0, 0) = a1.dot(JacMap0);           // dot product
 
  Eigen::Map<const AmgVector(5)> JacMap5(&J[5]);
  const AmgVector(5) a2 = M * JacMap5;
  m(1, 0) = a2.dot(JacMap0);
  m(1, 1) = a2.dot(JacMap5);
  m(0, 1) = m(1, 0);
 
  Eigen::Map<const AmgVector(5)> JacMap10(&J[10], 5, 1);
  const AmgVector(5) a3 = M * JacMap10;
  m(2, 0) = a3.dot(JacMap0);
  m(2, 1) = a3.dot(JacMap5);
  m(2, 2) = a3.dot(JacMap10);
  m(0, 2) = m(2, 0);
  m(1, 2) = m(2, 1);
 
  Eigen::Map<const AmgVector(5)> JacMap15(&J[15], 5, 1);
  const AmgVector(5) a4 = M * JacMap15;
  m(3, 0) = a4.dot(JacMap0);
  m(3, 1) = a4.dot(JacMap5);
  m(3, 2) = a4.dot(JacMap10);
  m(3, 3) = a4.dot(JacMap15);
  m(0, 3) = m(3, 0);
  m(1, 3) = m(3, 1);
  m(2, 3) = m(3, 2);
 
  const AmgVector(5) a5 = M.row(4) * J[20];
  m(4, 0) = a5.dot(JacMap0);
  m(4, 1) = a5.dot(JacMap5);
  m(4, 2) = a5.dot(JacMap10);
  m(4, 3) = a5.dot(JacMap15);
  m(4, 4) = a5[4] * J[20];
  m(0, 4) = m(4, 0);
  m(1, 4) = m(4, 1);
  m(2, 4) = m(4, 2);
  m(3, 4) = m(4, 3);
 
  return nM;
}