From f3ca711bd589beaca2a1d363c2bdd813fef5cfb6 Mon Sep 17 00:00:00 2001
From: Giovanni La Mura <giovanni.lamura@inaf.it>
Date: Thu, 12 Sep 2024 09:03:42 +0200
Subject: [PATCH] Fix f1 initialization bug in CBF and RBF

---
 src/libnptm/sph_subs.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/libnptm/sph_subs.cpp b/src/libnptm/sph_subs.cpp
index 4a639480..046c280d 100644
--- a/src/libnptm/sph_subs.cpp
+++ b/src/libnptm/sph_subs.cpp
@@ -123,7 +123,7 @@ void cbf(int n, dcomplex z, int &nm, dcomplex *csj) {
   if (m < n) nm = m;
   else m = msta2(a0, n, 15);
   dcomplex cf0 = 0.0;
-  dcomplex cf1 = 1.0e-100;
+  dcomplex cf1 = 1.0e+00 - 100.0;
   dcomplex cf, cs;
   for (int k = m; k >= 0; k--) {
     cf = (2.0 * k + 3.0) * cf1 / z - cf0;
@@ -668,7 +668,7 @@ void rbf(int n, double x, int &nm, double sj[]) {
   if (m < n) nm = m;
   else m = msta2(a0, n, 15);
   double f0 = 0.0;
-  double f1 = 1.0e-100;
+  double f1 = 1.0e+00 - 100.0;
   double f;
   for (int k = m; k >= 0; k--) {
     f = (2.0 * k +3.0) * f1 / x - f0;
-- 
GitLab