Set Rconv based on order

If order > 5, Rconv =5. Otherwise it is 7
mbakker7 · Apr 16, 2024 · 26b6d68 · 26b6d68
1 parent d534aa0
commit 26b6d68
Showing 1 changed file with 16 additions and 4 deletions.
diff --git a/timml/besselaesnumba/besselaesnumba.py b/timml/besselaesnumba/besselaesnumba.py
@@ -121,7 +121,10 @@ def potbeslsho(x, y, z1, z2, labda, order, ilap, naq):
     lstype = 1
 
     # Radius of convergence
-    Rconv = 5.0
+    if order > 5:
+        Rconv = 5.0
+    else:
+        Rconv = 7.0
 
     # if (ilap==1) :
     #    istart = 1
@@ -183,7 +186,10 @@ def disbeslsho(x, y, z1, z2, labda, order, ilap, naq):
 
     rv = np.zeros((2, naq))
     # Radius of convergence
-    Rconv = 5.0
+    if order > 5:
+        Rconv = 5.0
+    else:
+        Rconv = 7.0
 
     # lstype = 1 means line-sink
     lstype = 1
@@ -267,7 +273,10 @@ def potbesldho(x, y, z1, z2, labda, order, ilap, naq):
     rv = np.zeros(naq)
 
     # Radius of convergence
-    Rconv = 5.0
+    if order > 5:
+        Rconv = 5.0
+    else:
+        Rconv = 7.0
 
     # lstype=2 means line-doublet
     lstype = 2
@@ -334,7 +343,10 @@ def disbesldho(x, y, z1, z2, labda, order, ilap, naq):
     #                   and mod.Helmholtz potentials in remaining spots
     rv = np.zeros((2, naq))
     # Radius of convergence
-    Rconv = 5.0
+    if order > 5:
+        Rconv = 5.0
+    else:
+        Rconv = 7.0
 
     # lstype=2 means line-doublet
     lstype = 2