Convergence line elements large order

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

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:

timml/version.py

@@ -1 +1 @@
-__version__ = "6.4.1"
+__version__ = "6.5.0.dev0"