Skip to content
Snippets Groups Projects
Commit 75b55554 authored by Wuttke, Joachim's avatar Wuttke, Joachim
Browse files

ditto

parent 1c6c103c
No related branches found
No related tags found
1 merge request!145new class SliceStack
...@@ -15,9 +15,15 @@ ...@@ -15,9 +15,15 @@
#include "Sample/Specular/SpecularMagneticNCStrategy.h" #include "Sample/Specular/SpecularMagneticNCStrategy.h"
namespace { namespace {
complex_t checkForUnderflow(complex_t val);
complex_t checkForUnderflow(complex_t val)
{
return std::abs(val.imag()) < 1e-80 && val.real() < 0 ? complex_t(val.real(), 1e-40) : val;
} }
} // namespace
std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> SpecularMagneticNCStrategy::computeRoughnessMatrices( std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> SpecularMagneticNCStrategy::computeRoughnessMatrices(
const MatrixRTCoefficients& coeff_i, const MatrixRTCoefficients& coeff_i1, double sigma) const const MatrixRTCoefficients& coeff_i, const MatrixRTCoefficients& coeff_i1, double sigma) const
{ {
...@@ -82,10 +88,3 @@ SpecularMagneticNCStrategy::computeBackwardsSubmatrices(const MatrixRTCoefficien ...@@ -82,10 +88,3 @@ SpecularMagneticNCStrategy::computeBackwardsSubmatrices(const MatrixRTCoefficien
return {mp, mm}; return {mp, mm};
} }
namespace {
complex_t checkForUnderflow(complex_t val)
{
return std::abs(val.imag()) < 1e-80 && val.real() < 0 ? complex_t(val.real(), 1e-40) : val;
}
} // namespace
...@@ -17,9 +17,12 @@ ...@@ -17,9 +17,12 @@
#include "Base/Math/Functions.h" #include "Base/Math/Functions.h"
namespace { namespace {
const double pi2_15 = std::pow(M_PI_2, 1.5); const double pi2_15 = std::pow(M_PI_2, 1.5);
} // namespace } // namespace
Eigen::Matrix2cd Eigen::Matrix2cd
SpecularMagneticTanhStrategy::computeRoughnessMatrix(const MatrixRTCoefficients& coeff, SpecularMagneticTanhStrategy::computeRoughnessMatrix(const MatrixRTCoefficients& coeff,
double sigma, bool inverse) const double sigma, bool inverse) const
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment