From 9a6e2760ec8ab6bdf01019371d8895d9989d770a Mon Sep 17 00:00:00 2001 From: Randolf Beerwerth <r.beerwerth@fz-juelich.de> Date: Mon, 30 Nov 2020 17:42:21 +0100 Subject: [PATCH] Rename MatrixRTCoefficients_v3 -> MatrixRTCoefficients --- ...cients_v3.cpp => MatrixRTCoefficients.cpp} | 48 +++++++++---------- ...efficients_v3.h => MatrixRTCoefficients.h} | 20 ++++---- .../Specular/SpecularMagneticNCStrategy.cpp | 8 ++-- Sample/Specular/SpecularMagneticNCStrategy.h | 8 ++-- Sample/Specular/SpecularMagneticStrategy.cpp | 8 ++-- Sample/Specular/SpecularMagneticStrategy.h | 10 ++-- .../Specular/SpecularMagneticTanhStrategy.cpp | 4 +- .../Specular/SpecularMagneticTanhStrategy.h | 6 +-- ...3Test.cpp => MatrixRTCoefficientsTest.cpp} | 24 +++++----- 9 files changed, 68 insertions(+), 68 deletions(-) rename Sample/RT/{MatrixRTCoefficients_v3.cpp => MatrixRTCoefficients.cpp} (70%) rename Sample/RT/{MatrixRTCoefficients_v3.h => MatrixRTCoefficients.h} (81%) rename Tests/UnitTests/Core/Fresnel/{MatrixRTCoefficients_v3Test.cpp => MatrixRTCoefficientsTest.cpp} (66%) diff --git a/Sample/RT/MatrixRTCoefficients_v3.cpp b/Sample/RT/MatrixRTCoefficients.cpp similarity index 70% rename from Sample/RT/MatrixRTCoefficients_v3.cpp rename to Sample/RT/MatrixRTCoefficients.cpp index 0e3ccbf5bbd..cced7f1f1de 100644 --- a/Sample/RT/MatrixRTCoefficients_v3.cpp +++ b/Sample/RT/MatrixRTCoefficients.cpp @@ -2,8 +2,8 @@ // // BornAgain: simulate and fit scattering at grazing incidence // -//! @file Sample/RT/MatrixRTCoefficients_v3.cpp -//! @brief Implements class MatrixRTCoefficients_v3. +//! @file Sample/RT/MatrixRTCoefficients.cpp +//! @brief Implements class MatrixRTCoefficients. //! //! @homepage http://www.bornagainproject.org //! @license GNU General Public License v3 or higher (see COPYING) @@ -12,7 +12,7 @@ // // ************************************************************************************************ -#include "Sample/RT/MatrixRTCoefficients_v3.h" +#include "Sample/RT/MatrixRTCoefficients.h" #include "Base/Utils/Assert.h" namespace { @@ -20,7 +20,7 @@ complex_t GetImExponential(complex_t exponent); const auto eps = std::numeric_limits<double>::epsilon() * 10.; } // namespace -MatrixRTCoefficients_v3::MatrixRTCoefficients_v3(double kz_sign, Eigen::Vector2cd eigenvalues, +MatrixRTCoefficients::MatrixRTCoefficients(double kz_sign, Eigen::Vector2cd eigenvalues, kvector_t b, double magnetic_SLD) : m_kz_sign(kz_sign) , m_lambda(std::move(eigenvalues)) @@ -32,15 +32,15 @@ MatrixRTCoefficients_v3::MatrixRTCoefficients_v3(double kz_sign, Eigen::Vector2c m_R << -1, 0, 0, -1; } -MatrixRTCoefficients_v3::MatrixRTCoefficients_v3(const MatrixRTCoefficients_v3& other) = default; +MatrixRTCoefficients::MatrixRTCoefficients(const MatrixRTCoefficients& other) = default; -MatrixRTCoefficients_v3::~MatrixRTCoefficients_v3() = default; +MatrixRTCoefficients::~MatrixRTCoefficients() = default; -MatrixRTCoefficients_v3* MatrixRTCoefficients_v3::clone() const { - return new MatrixRTCoefficients_v3(*this); +MatrixRTCoefficients* MatrixRTCoefficients::clone() const { + return new MatrixRTCoefficients(*this); } -Eigen::Matrix2cd MatrixRTCoefficients_v3::TransformationMatrix(Eigen::Vector2d selection) const { +Eigen::Matrix2cd MatrixRTCoefficients::TransformationMatrix(Eigen::Vector2d selection) const { const Eigen::Matrix2cd exp2 = Eigen::DiagonalMatrix<complex_t, 2>(selection); if (std::abs(m_b.mag() - 1.) < eps) { @@ -55,51 +55,51 @@ Eigen::Matrix2cd MatrixRTCoefficients_v3::TransformationMatrix(Eigen::Vector2d s throw std::runtime_error("Broken magnetic field vector"); } -Eigen::Matrix2cd MatrixRTCoefficients_v3::T1Matrix() const { +Eigen::Matrix2cd MatrixRTCoefficients::T1Matrix() const { return TransformationMatrix({0., 1.}); } -Eigen::Matrix2cd MatrixRTCoefficients_v3::T2Matrix() const { +Eigen::Matrix2cd MatrixRTCoefficients::T2Matrix() const { return TransformationMatrix({1., 0.}); } -Eigen::Vector2cd MatrixRTCoefficients_v3::T1plus() const { +Eigen::Vector2cd MatrixRTCoefficients::T1plus() const { return T1Matrix() * m_T.col(0); } -Eigen::Vector2cd MatrixRTCoefficients_v3::R1plus() const { +Eigen::Vector2cd MatrixRTCoefficients::R1plus() const { return T1Matrix() * m_R.col(0); } -Eigen::Vector2cd MatrixRTCoefficients_v3::T2plus() const { +Eigen::Vector2cd MatrixRTCoefficients::T2plus() const { return T2Matrix() * m_T.col(0); } -Eigen::Vector2cd MatrixRTCoefficients_v3::R2plus() const { +Eigen::Vector2cd MatrixRTCoefficients::R2plus() const { return T2Matrix() * m_R.col(0); } -Eigen::Vector2cd MatrixRTCoefficients_v3::T1min() const { +Eigen::Vector2cd MatrixRTCoefficients::T1min() const { return T1Matrix() * m_T.col(1); } -Eigen::Vector2cd MatrixRTCoefficients_v3::R1min() const { +Eigen::Vector2cd MatrixRTCoefficients::R1min() const { return T1Matrix() * m_R.col(1); } -Eigen::Vector2cd MatrixRTCoefficients_v3::T2min() const { +Eigen::Vector2cd MatrixRTCoefficients::T2min() const { return T2Matrix() * m_T.col(1); } -Eigen::Vector2cd MatrixRTCoefficients_v3::R2min() const { +Eigen::Vector2cd MatrixRTCoefficients::R2min() const { return T2Matrix() * m_R.col(1); } -Eigen::Vector2cd MatrixRTCoefficients_v3::getKz() const { +Eigen::Vector2cd MatrixRTCoefficients::getKz() const { return m_kz_sign * m_lambda; } -Eigen::Matrix2cd MatrixRTCoefficients_v3::pMatrixHelper(double sign) const { +Eigen::Matrix2cd MatrixRTCoefficients::pMatrixHelper(double sign) const { const complex_t alpha = m_lambda(1) + m_lambda(0); const complex_t beta = m_lambda(1) - m_lambda(0); @@ -113,14 +113,14 @@ Eigen::Matrix2cd MatrixRTCoefficients_v3::pMatrixHelper(double sign) const { return result; } -Eigen::Matrix2cd MatrixRTCoefficients_v3::computeP() const { +Eigen::Matrix2cd MatrixRTCoefficients::computeP() const { Eigen::Matrix2cd result = pMatrixHelper(1.); result *= 0.5; return result; } -Eigen::Matrix2cd MatrixRTCoefficients_v3::computeInverseP() const { +Eigen::Matrix2cd MatrixRTCoefficients::computeInverseP() const { const complex_t alpha = m_lambda(1) + m_lambda(0); const complex_t beta = m_lambda(1) - m_lambda(0); @@ -133,7 +133,7 @@ Eigen::Matrix2cd MatrixRTCoefficients_v3::computeInverseP() const { return result; } -Eigen::Matrix2cd MatrixRTCoefficients_v3::computeDeltaMatrix(double thickness) { +Eigen::Matrix2cd MatrixRTCoefficients::computeDeltaMatrix(double thickness) { Eigen::Matrix2cd result; const complex_t alpha = 0.5 * thickness * (m_lambda(1) + m_lambda(0)); diff --git a/Sample/RT/MatrixRTCoefficients_v3.h b/Sample/RT/MatrixRTCoefficients.h similarity index 81% rename from Sample/RT/MatrixRTCoefficients_v3.h rename to Sample/RT/MatrixRTCoefficients.h index b7f104e35ae..a7c3858ebb9 100644 --- a/Sample/RT/MatrixRTCoefficients_v3.h +++ b/Sample/RT/MatrixRTCoefficients.h @@ -2,8 +2,8 @@ // // BornAgain: simulate and fit scattering at grazing incidence // -//! @file Sample/RT/MatrixRTCoefficients_v3.h -//! @brief Defines class MatrixRTCoefficients_v3. +//! @file Sample/RT/MatrixRTCoefficients.h +//! @brief Defines class MatrixRTCoefficients. //! //! @homepage http://www.bornagainproject.org //! @license GNU General Public License v3 or higher (see COPYING) @@ -12,8 +12,8 @@ // // ************************************************************************************************ -#ifndef BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_V3_H -#define BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_V3_H +#ifndef BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_H +#define BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_H #include "Base/Vector/Vectors3D.h" #include "Sample/RT/ILayerRTCoefficients.h" @@ -23,19 +23,19 @@ //! of magnetic interactions between the scattered particle and the layer. //! @ingroup algorithms_internal -class MatrixRTCoefficients_v3 : public ILayerRTCoefficients { +class MatrixRTCoefficients : public ILayerRTCoefficients { public: friend class SpecularMagneticStrategy; friend class SpecularMagneticNCStrategy; friend class SpecularMagneticNCTestingStrategy; friend class SpecularMagneticTanhStrategy; - MatrixRTCoefficients_v3(double kz_sign, Eigen::Vector2cd eigenvalues, kvector_t b, + MatrixRTCoefficients(double kz_sign, Eigen::Vector2cd eigenvalues, kvector_t b, double magnetic_SLD); - MatrixRTCoefficients_v3(const MatrixRTCoefficients_v3& other); - ~MatrixRTCoefficients_v3() override; + MatrixRTCoefficients(const MatrixRTCoefficients& other); + ~MatrixRTCoefficients() override; - MatrixRTCoefficients_v3* clone() const override; + MatrixRTCoefficients* clone() const override; //! The following functions return the transmitted and reflected amplitudes //! for different incoming beam polarizations and eigenmodes @@ -75,4 +75,4 @@ private: Eigen::Matrix2cd pMatrixHelper(double sign) const; }; -#endif // BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_V3_H +#endif // BORNAGAIN_SAMPLE_RT_MATRIXRTCOEFFICIENTS_H diff --git a/Sample/Specular/SpecularMagneticNCStrategy.cpp b/Sample/Specular/SpecularMagneticNCStrategy.cpp index f7495c9f9be..85746b1327d 100644 --- a/Sample/Specular/SpecularMagneticNCStrategy.cpp +++ b/Sample/Specular/SpecularMagneticNCStrategy.cpp @@ -19,8 +19,8 @@ complex_t checkForUnderflow(complex_t val); } std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> -SpecularMagneticNCStrategy::computeRoughnessMatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, +SpecularMagneticNCStrategy::computeRoughnessMatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const { complex_t beta_i = coeff_i.m_lambda(1) - coeff_i.m_lambda(0); complex_t beta_i1 = coeff_i1.m_lambda(1) - coeff_i1.m_lambda(0); @@ -65,8 +65,8 @@ SpecularMagneticNCStrategy::computeRoughnessMatrices(const MatrixRTCoefficients_ } std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> -SpecularMagneticNCStrategy::computeBackwardsSubmatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, +SpecularMagneticNCStrategy::computeBackwardsSubmatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const { Eigen::Matrix2cd roughness_sum{Eigen::Matrix2cd::Identity()}; Eigen::Matrix2cd roughness_diff{Eigen::Matrix2cd::Identity()}; diff --git a/Sample/Specular/SpecularMagneticNCStrategy.h b/Sample/Specular/SpecularMagneticNCStrategy.h index dd04616317b..573c5b9e7a3 100644 --- a/Sample/Specular/SpecularMagneticNCStrategy.h +++ b/Sample/Specular/SpecularMagneticNCStrategy.h @@ -30,12 +30,12 @@ class SpecularMagneticNCStrategy : public SpecularMagneticStrategy { private: std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> - computeRoughnessMatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, double sigma) const; + computeRoughnessMatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const; virtual std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> - computeBackwardsSubmatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, double sigma) const; + computeBackwardsSubmatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const; }; #endif // BORNAGAIN_SAMPLE_SPECULAR_SPECULARMAGNETICNCSTRATEGY_H diff --git a/Sample/Specular/SpecularMagneticStrategy.cpp b/Sample/Specular/SpecularMagneticStrategy.cpp index db5049df637..e940ff0093f 100644 --- a/Sample/Specular/SpecularMagneticStrategy.cpp +++ b/Sample/Specular/SpecularMagneticStrategy.cpp @@ -44,12 +44,12 @@ SpecularMagneticStrategy::Execute(const std::vector<Slice>& slices, ISpecularStrategy::coeffs_t result; for (auto& coeff : computeTR(slices, kz)) - result.push_back(std::make_unique<MatrixRTCoefficients_v3>(coeff)); + result.push_back(std::make_unique<MatrixRTCoefficients>(coeff)); return result; } -std::vector<MatrixRTCoefficients_v3> +std::vector<MatrixRTCoefficients> SpecularMagneticStrategy::computeTR(const std::vector<Slice>& slices, const std::vector<complex_t>& kzs) const { const size_t N = slices.size(); @@ -60,7 +60,7 @@ SpecularMagneticStrategy::computeTR(const std::vector<Slice>& slices, if (slices.empty()) return {}; - std::vector<MatrixRTCoefficients_v3> result; + std::vector<MatrixRTCoefficients> result; result.reserve(N); const double kz_sign = kzs.front().real() >= 0.0 ? 1.0 : -1.0; // save sign to restore it later @@ -94,7 +94,7 @@ SpecularMagneticStrategy::computeTR(const std::vector<Slice>& slices, return result; } -void SpecularMagneticStrategy::calculateUpwards(std::vector<MatrixRTCoefficients_v3>& coeff, +void SpecularMagneticStrategy::calculateUpwards(std::vector<MatrixRTCoefficients>& coeff, const std::vector<Slice>& slices) const { const auto N = slices.size(); std::vector<Eigen::Matrix2cd> SMatrices(N - 1); diff --git a/Sample/Specular/SpecularMagneticStrategy.h b/Sample/Specular/SpecularMagneticStrategy.h index 573e4f6e434..d9fcc4a9c84 100644 --- a/Sample/Specular/SpecularMagneticStrategy.h +++ b/Sample/Specular/SpecularMagneticStrategy.h @@ -15,7 +15,7 @@ #ifndef BORNAGAIN_SAMPLE_SPECULAR_SPECULARMAGNETICSTRATEGY_H #define BORNAGAIN_SAMPLE_SPECULAR_SPECULARMAGNETICSTRATEGY_H -#include "Sample/RT/MatrixRTCoefficients_v3.h" +#include "Sample/RT/MatrixRTCoefficients.h" #include "Sample/Specular/ISpecularStrategy.h" #include <memory> #include <vector> @@ -42,14 +42,14 @@ public: const std::vector<complex_t>& kz) const; private: - std::vector<MatrixRTCoefficients_v3> computeTR(const std::vector<Slice>& slices, + std::vector<MatrixRTCoefficients> computeTR(const std::vector<Slice>& slices, const std::vector<complex_t>& kzs) const; virtual std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> - computeBackwardsSubmatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, double sigma) const = 0; + computeBackwardsSubmatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const = 0; - void calculateUpwards(std::vector<MatrixRTCoefficients_v3>& coeff, + void calculateUpwards(std::vector<MatrixRTCoefficients>& coeff, const std::vector<Slice>& slices) const; }; diff --git a/Sample/Specular/SpecularMagneticTanhStrategy.cpp b/Sample/Specular/SpecularMagneticTanhStrategy.cpp index b964ba6b435..6ddc26a6b29 100644 --- a/Sample/Specular/SpecularMagneticTanhStrategy.cpp +++ b/Sample/Specular/SpecularMagneticTanhStrategy.cpp @@ -21,7 +21,7 @@ const double pi2_15 = std::pow(M_PI_2, 1.5); } // namespace Eigen::Matrix2cd -SpecularMagneticTanhStrategy::computeRoughnessMatrix(const MatrixRTCoefficients_v3& coeff, +SpecularMagneticTanhStrategy::computeRoughnessMatrix(const MatrixRTCoefficients& coeff, double sigma, bool inverse) const { if (sigma < 10 * std::numeric_limits<double>::epsilon()) return Eigen::Matrix2cd{Eigen::Matrix2cd::Identity()}; @@ -61,7 +61,7 @@ SpecularMagneticTanhStrategy::computeRoughnessMatrix(const MatrixRTCoefficients_ std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> SpecularMagneticTanhStrategy::computeBackwardsSubmatrices( - const MatrixRTCoefficients_v3& coeff_i, const MatrixRTCoefficients_v3& coeff_i1, + const MatrixRTCoefficients& coeff_i, const MatrixRTCoefficients& coeff_i1, double sigma) const { Eigen::Matrix2cd R{Eigen::Matrix2cd::Identity()}; Eigen::Matrix2cd RInv{Eigen::Matrix2cd::Identity()}; diff --git a/Sample/Specular/SpecularMagneticTanhStrategy.h b/Sample/Specular/SpecularMagneticTanhStrategy.h index 6b13a7d930e..cb10b10636b 100644 --- a/Sample/Specular/SpecularMagneticTanhStrategy.h +++ b/Sample/Specular/SpecularMagneticTanhStrategy.h @@ -28,10 +28,10 @@ class SpecularMagneticTanhStrategy : public SpecularMagneticStrategy { private: virtual std::pair<Eigen::Matrix2cd, Eigen::Matrix2cd> - computeBackwardsSubmatrices(const MatrixRTCoefficients_v3& coeff_i, - const MatrixRTCoefficients_v3& coeff_i1, double sigma) const; + computeBackwardsSubmatrices(const MatrixRTCoefficients& coeff_i, + const MatrixRTCoefficients& coeff_i1, double sigma) const; - Eigen::Matrix2cd computeRoughnessMatrix(const MatrixRTCoefficients_v3& coeff, double sigma, + Eigen::Matrix2cd computeRoughnessMatrix(const MatrixRTCoefficients& coeff, double sigma, bool inverse = false) const; }; diff --git a/Tests/UnitTests/Core/Fresnel/MatrixRTCoefficients_v3Test.cpp b/Tests/UnitTests/Core/Fresnel/MatrixRTCoefficientsTest.cpp similarity index 66% rename from Tests/UnitTests/Core/Fresnel/MatrixRTCoefficients_v3Test.cpp rename to Tests/UnitTests/Core/Fresnel/MatrixRTCoefficientsTest.cpp index e237eb7ba99..ba2843d4118 100644 --- a/Tests/UnitTests/Core/Fresnel/MatrixRTCoefficients_v3Test.cpp +++ b/Tests/UnitTests/Core/Fresnel/MatrixRTCoefficientsTest.cpp @@ -1,60 +1,60 @@ -#include "Sample/RT/MatrixRTCoefficients_v3.h" +#include "Sample/RT/MatrixRTCoefficients.h" #include "Tests/GTestWrapper/google_test.h" -class MatrixRTCoefficients_v3Test : public ::testing::Test { +class MatrixRTCoefficientsTest : public ::testing::Test { protected: - MatrixRTCoefficients_v3 mrtcDefault{-1., {0., 0.}, kvector_t{0.0, 0.0, 0.0}, 0.}; + MatrixRTCoefficients mrtcDefault{-1., {0., 0.}, kvector_t{0.0, 0.0, 0.0}, 0.}; }; -TEST_F(MatrixRTCoefficients_v3Test, T1plus) { +TEST_F(MatrixRTCoefficientsTest, T1plus) { Eigen::Vector2cd vector = mrtcDefault.T1plus(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, T1min) { +TEST_F(MatrixRTCoefficientsTest, T1min) { Eigen::Vector2cd vector = mrtcDefault.T1min(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(complex_t(1.0, 0.0), vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, T2plus) { +TEST_F(MatrixRTCoefficientsTest, T2plus) { Eigen::Vector2cd vector = mrtcDefault.T2plus(); EXPECT_EQ(complex_t(1.0, 0.0), vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, T2min) { +TEST_F(MatrixRTCoefficientsTest, T2min) { Eigen::Vector2cd vector = mrtcDefault.T2min(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, R1plus) { +TEST_F(MatrixRTCoefficientsTest, R1plus) { Eigen::Vector2cd vector = mrtcDefault.R1plus(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, R1min) { +TEST_F(MatrixRTCoefficientsTest, R1min) { Eigen::Vector2cd vector = mrtcDefault.R1min(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(complex_t(-1.0, 0.0), vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, R2plus) { +TEST_F(MatrixRTCoefficientsTest, R2plus) { Eigen::Vector2cd vector = mrtcDefault.R2plus(); EXPECT_EQ(complex_t(-1.0, 0.0), vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, R2min) { +TEST_F(MatrixRTCoefficientsTest, R2min) { Eigen::Vector2cd vector = mrtcDefault.R2min(); EXPECT_EQ(0.0, vector(0)); EXPECT_EQ(0.0, vector(1)); } -TEST_F(MatrixRTCoefficients_v3Test, getKz) { +TEST_F(MatrixRTCoefficientsTest, getKz) { Eigen::Vector2cd vector = mrtcDefault.getKz(); EXPECT_EQ(complex_t(0.0, 0.0), vector(0)); EXPECT_EQ(complex_t(0.0, 0.0), vector(1)); -- GitLab