From c3e17e03c88d084d0b8e95ffc190df62b49cb5ef Mon Sep 17 00:00:00 2001
From: Mikhail Svechnikov <m.svechnikov@fz-juelich.de>
Date: Fri, 9 Aug 2024 14:24:52 +0200
Subject: [PATCH] use magxy
---
Sample/Interface/AutoCorrModels.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/Sample/Interface/AutoCorrModels.cpp b/Sample/Interface/AutoCorrModels.cpp
index 799a64fe97f..63ee626ea3a 100644
--- a/Sample/Interface/AutoCorrModels.cpp
+++ b/Sample/Interface/AutoCorrModels.cpp
@@ -67,7 +67,7 @@ double BasicAutoCorr::spectralFunction(const R3& k) const
ASSERT(m_validated);
double H = m_hurst_parameter;
double clength2 = m_lateral_corr_length * m_lateral_corr_length;
- double Qpar2 = k.x() * k.x() + k.y() * k.y();
+ double Qpar2 = k.magxy2();
return 4.0 * pi * H * m_sigma * m_sigma * clength2 * std::pow(1 + Qpar2 * clength2, -1 - H);
}
@@ -77,7 +77,7 @@ double BasicAutoCorr::corrFunction(const R3& k) const
ASSERT(m_validated);
double H = m_hurst_parameter;
double clength = m_lateral_corr_length;
- double R = sqrt(k.x() * k.x() + k.y() * k.y());
+ double R = k.magxy();
return m_sigma * m_sigma * std::pow(2., 1 - H) / tgamma(H) * std::pow(R / clength, H)
* gsl_sf_bessel_Knu(H, R / clength);
}
--
GitLab