From 26b6d687437db5e0287196eafd5ee8618df4013a Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Tue, 16 Apr 2024 19:18:13 +0200 Subject: [PATCH] Set Rconv based on order If order > 5, Rconv =5. Otherwise it is 7 --- timml/besselaesnumba/besselaesnumba.py | 20 ++++++++++++++++---- 1 file changed, 16 insertions(+), 4 deletions(-) diff --git a/timml/besselaesnumba/besselaesnumba.py b/timml/besselaesnumba/besselaesnumba.py index 46e9877d..4c1cef62 100644 --- 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