From f7ad16a17b1e4d3f649e53943b3c09af7199bb39 Mon Sep 17 00:00:00 2001
From: "Joachim Wuttke (h)" <j.wuttke@fz-juelich.de>
Date: Tue, 22 Nov 2022 10:46:37 +0100
Subject: [PATCH] reuse sqrt(1-z^2)
---
Sample/HardParticle/HemiEllipsoid.cpp | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/Sample/HardParticle/HemiEllipsoid.cpp b/Sample/HardParticle/HemiEllipsoid.cpp
index 33549ed7d91..1714be15a79 100644
--- a/Sample/HardParticle/HemiEllipsoid.cpp
+++ b/Sample/HardParticle/HemiEllipsoid.cpp
@@ -59,7 +59,8 @@ complex_t HemiEllipsoid::formfactor_at_bottom(C3 q) const
complex_t qxRz = q.x() * Rz;
complex_t qyWz = q.y() * Wz;
- complex_t gamma = std::sqrt(qxRz * qxRz + qyWz * qyWz);
+ complex_t gamma = std::sqrt((q.x() * R)*(q.x() * R)
+ + (q.y() * W) * (q.y() * W)) * std::sqrt(1.0 - z * z);
complex_t J1_gamma_div_gamma = Math::Bessel::J1c(gamma);
return Rz * Wz * J1_gamma_div_gamma * exp_I(q.z() * H * z);
--
GitLab