Merge pull request #113 from mbakker7/convergence_large_radius

Convergence large radius
mbakker7 · Sep 26, 2024 · fd0d1d3 · fd0d1d3
2 parents a7cb4c9 + 26b6d68
commit fd0d1d3
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 5 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 = 7.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 = 7.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 = 7.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 = 7.0
+    if order > 5:
+        Rconv = 5.0
+    else:
+        Rconv = 7.0
 
     # lstype=2 means line-doublet
     lstype = 2

diff --git a/timml/constant.py b/timml/constant.py
@@ -153,7 +153,7 @@ def equation(self):
                         mat[0:, ieq : ieq + e.nunknowns] += e.potinflayers(
                             self.xc[icp], self.yc[icp], self.layers
                         ).sum(0)
-                        ieq += e.nunknowns
+                    ieq += e.nunknowns # I decreased the tab here
                         # else:
                         #    mat[0, ieq:ieq+e. nunknowns] += -1
                 else: