From d183b813c0c3494972bff74da9c0052c2930203f Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Thu, 12 Sep 2024 09:42:49 +0200
Subject: [PATCH] Revert to original square radius implementation in POLAR
---
src/libnptm/clu_subs.cpp | 8 ++++++--
1 file changed, 6 insertions(+), 2 deletions(-)
diff --git a/src/libnptm/clu_subs.cpp b/src/libnptm/clu_subs.cpp
index 42d136b2..78bfc4af 100644
--- a/src/libnptm/clu_subs.cpp
+++ b/src/libnptm/clu_subs.cpp
@@ -1473,21 +1473,25 @@ void polar(
bool onx = (y == 0.0);
bool ony = (x == 0.0);
bool onz = (onx && ony);
+ double rhos = 0.0;
double rho = 0.0;
if (!onz) {
if (!onx) {
if (!ony) {
- rho = sqrt(x * x + y * y);
+ rhos = x * x + y * y;
+ rho = sqrt(rhos);
cph = x / rho;
sph = y / rho;
// goes to 25
} else { // label 20
+ rhos = y * y;
rho = (y > 0.0) ? y : -y;
cph = 0.0;
sph = (y > 0.0) ? 1.0 : -1.0;
// goes to 25
}
} else { // label 15
+ rhos = x * x;
rho = (x > 0.0) ? x : -x;
cph = (x > 0.0) ? 1.0 : -1.0;
sph = 0.0;
@@ -1513,7 +1517,7 @@ void polar(
}
} else { // label 35
if (!onz) {
- r = sqrt(rho * rho + z * z);
+ r = sqrt(rhos + z * z);
cth = z / r;
sth = rho / r;
// returns
--
GitLab