IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 59e82949 authored by Maximilian Reininghaus's avatar Maximilian Reininghaus :vulcan: Committed by ralfulrich
Browse files

removed utl::sgn usage

parent 1eca04b6
No related branches found
No related tags found
No related merge requests found
......@@ -14,7 +14,6 @@
#include <Eigen/Dense>
#include <stdexcept>
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/utility/sgn.hpp>
namespace corsika {
......@@ -38,14 +37,12 @@ namespace corsika {
auto CoordinateSystem::RotateToZ(Vector<TDim> vVec) const
{
auto const a = vVec.normalized().GetComponents(*this).eVector;
auto const a1 = a(0), a2 = a(1);
auto const s = corsika::sgn(a(2));
auto const c = 1 / (1 + s * a(2));
auto const a1 = a(0), a2 = a(1), a3 = a(2);
Eigen::Matrix3d A, B;
if (s > 0) {
if (a3 > 0) {
auto const c = 1 / (1 + a3);
A << 1, 0, -a1, // comment to prevent clang-format
0, 1, -a2, // .
a1, a2, 1; // .
......@@ -54,6 +51,7 @@ namespace corsika {
0, 0, -(a1 * a1 + a2 * a2) * c; // .
} else {
auto const c = 1 / (1 - a3);
A << 1, 0, a1, // .
0, -1, -a2, // .
a1, a2, -1; // .
......
......@@ -18,7 +18,6 @@
#include <corsika/framework/geometry/FourVector.hpp>
#include <corsika/framework/geometry/Vector.hpp>
#include <corsika/framework/logging/Logging.h>
#include <corsika/framework/utility/sgn.hpp>
// using namespace corsika::units::si;
......
......@@ -10,7 +10,6 @@ n/*
#include <corsika/framework/core/PhysicalUnits.hpp>
#include <corsika/framework/geometry/QuantityVector.hpp>
#include <corsika/framework/utility/sgn.hpp>
#include <Eigen/Dense>
#include <stdexcept>
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment