From aa3e78342c98070b82dd5766d70956d6c052c5e6 Mon Sep 17 00:00:00 2001
From: Randolf Beerwerth <r.beerwerth@fz-juelich.de>
Date: Tue, 22 Sep 2020 12:47:37 +0200
Subject: [PATCH] Remove unnecessary cutoff
---
Core/Multilayer/SpecularMagneticNewStrategy.cpp | 3 ++-
Core/RT/MatrixRTCoefficients_v3.cpp | 3 +--
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/Core/Multilayer/SpecularMagneticNewStrategy.cpp b/Core/Multilayer/SpecularMagneticNewStrategy.cpp
index 6dac8e0bf6e..79dc7f6225d 100644
--- a/Core/Multilayer/SpecularMagneticNewStrategy.cpp
+++ b/Core/Multilayer/SpecularMagneticNewStrategy.cpp
@@ -82,7 +82,8 @@ SpecularMagneticNewStrategy::computeTR(const std::vector<Slice>& slices,
result[0].m_T = Eigen::Matrix2cd::Identity();
result[0].m_R = Eigen::Matrix2cd::Zero();
return result;
- }else if(kzs[0] == 0.0){
+
+ }else if( kzs[0] == 0. ){
result[0].m_T = Eigen::Matrix2cd::Identity();
result[0].m_R = -Eigen::Matrix2cd::Identity();
for (size_t i = 1; i < N; ++i) {
diff --git a/Core/RT/MatrixRTCoefficients_v3.cpp b/Core/RT/MatrixRTCoefficients_v3.cpp
index b72cfead875..98919c7a34b 100644
--- a/Core/RT/MatrixRTCoefficients_v3.cpp
+++ b/Core/RT/MatrixRTCoefficients_v3.cpp
@@ -144,8 +144,7 @@ Eigen::Matrix2cd MatrixRTCoefficients_v3::computeInverseP() const
const complex_t alpha = m_lambda(1) + m_lambda(0);
const complex_t beta = m_lambda(1) - m_lambda(0);
- if (std::abs(alpha * alpha - beta * beta) < eps)
-// throw std::runtime_error("Singular p_m");
+ if (std::abs(alpha * alpha - beta * beta) == 0.)
return Eigen::Matrix2cd::Zero();
Eigen::Matrix2cd result = pMatrixHelper(-1.);
--
GitLab