From 82f85f58a94baa02e32860732ddc43b6c689551f Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Tue, 19 Feb 2019 15:53:23 +0100
Subject: [PATCH] typo
---
Processes/TrackingLine/TrackingLine.cc | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/Processes/TrackingLine/TrackingLine.cc b/Processes/TrackingLine/TrackingLine.cc
index e19d9077a..b540dfcc9 100644
--- a/Processes/TrackingLine/TrackingLine.cc
+++ b/Processes/TrackingLine/TrackingLine.cc
@@ -32,11 +32,12 @@ namespace corsika::process::tracking_line {
geometry::Sphere const& sphere) {
auto const delta = line.GetR0() - sphere.GetCenter();
auto const v = line.GetV0();
+ auto const vSqNorm = v.squaredNorm();
auto const R = sphere.GetRadius();
auto const vDotDelta = v.dot(delta);
auto const discriminant =
- vDotDelta * vDotDelta - v.squaredNorm() * (delta.squaredNorm() - R * R);
+ vDotDelta * vDotDelta - vSqNorm * (delta.squaredNorm() - R * R);
//~ std::cout << "discriminant: " << discriminant << std::endl;
//~ std::cout << "alpha: " << alpha << std::endl;
@@ -44,7 +45,7 @@ namespace corsika::process::tracking_line {
if (discriminant.magnitude() > 0) {
auto const sqDisc = sqrt(discriminant);
- auto const invDenom = 1 / v0.squaredNorm();
+ auto const invDenom = 1 / vSqNorm;
return std::make_pair((vDotDelta - sqDisc) * invDenom),
(vDotDelta + sqDisc) * invDenom));
} else {
--
GitLab