Skip to content

Commit 2dbf448

Browse files
committed
ITS3: exploit symmetric matrices for interpolation
Signed-off-by: Felix Schlepper <felix.schlepper@cern.ch>
1 parent 2cec687 commit 2dbf448

File tree

1 file changed

+5
-3
lines changed
  • Detectors/Upgrades/ITS3/alignment/include/ITS3Align

1 file changed

+5
-3
lines changed

Detectors/Upgrades/ITS3/alignment/include/ITS3Align/TrackFit.h

Lines changed: 5 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -76,10 +76,12 @@ o2::track::TrackParametrizationWithError<T> interpolateTrackParCov(
7676
};
7777
Mat55 cA = unpack(tA.getCov());
7878
Mat55 cB = unpack(tB.getCov());
79-
Mat55 wA = cA.inverse();
80-
Mat55 wB = cB.inverse();
79+
Eigen::LLT<Mat55> lltA(cA), lltB(cB);
80+
Mat55 wA = lltA.solve(Mat55::Identity());
81+
Mat55 wB = lltB.solve(Mat55::Identity());
8182
Mat55 wTot = wA + wB;
82-
Mat55 cTot = wTot.inverse();
83+
Eigen::LLT<Mat55> lltTot(wTot);
84+
Mat55 cTot = lltTot.solve(Mat55::Identity());
8385
Mat51 pA, pB;
8486
for (int i = 0; i < 5; ++i) {
8587
pA(i) = tA.getParam(i);

0 commit comments

Comments
 (0)