From 4e3be44203ea6b6de2cc27c2c579d30edf1cf643 Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Wed, 20 May 2026 17:32:40 +0200 Subject: [PATCH 1/4] Update besselnumba.py --- timflow/bessel/besselnumba.py | 1 + 1 file changed, 1 insertion(+) diff --git a/timflow/bessel/besselnumba.py b/timflow/bessel/besselnumba.py index 2d9c13d..3d648b9 100644 --- a/timflow/bessel/besselnumba.py +++ b/timflow/bessel/besselnumba.py @@ -22,6 +22,7 @@ fac = 1.0 +# max number of terms fixed to 20 nrange = np.arange(21, dtype=np.float64) a = np.zeros(21, dtype=np.float64) From ee64d67b44459153d3c1d8594a807f468e1e15a2 Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Wed, 20 May 2026 17:42:45 +0200 Subject: [PATCH 2/4] fix param_3d return values --- timflow/steady/inhomogeneity.py | 10 ++++++++-- 1 file changed, 8 insertions(+), 2 deletions(-) diff --git a/timflow/steady/inhomogeneity.py b/timflow/steady/inhomogeneity.py index 7744d4c..e20d93d 100644 --- a/timflow/steady/inhomogeneity.py +++ b/timflow/steady/inhomogeneity.py @@ -687,7 +687,9 @@ def __init__( layers = [0] if z is None: z = [1, 0] - (kaq, c, npor, ltype) = param_3d(kaq, z, kzoverkh, npor, topboundary, topres) + kaq, kzoverkh, c, npor, ltype = param_3d( + kaq, z, kzoverkh, npor, topboundary, topres + ) if topboundary == "semi": z = np.hstack((z[0] + topthick, z)) super().__init__( @@ -703,6 +705,7 @@ def __init__( ndeg=ndeg, layers=layers, ) + self.kzoverkh = kzoverkh # add kzoverkh to inhomogeneity object class LeakyBuildingPit(BuildingPit): @@ -1060,7 +1063,9 @@ def __init__( layers = [0] if z is None: z = [1, 0] - (kaq, c, npor, ltype) = param_3d(kaq, z, kzoverkh, npor, topboundary, topres) + kaq, kzoverkh, c, npor, ltype = param_3d( + kaq, z, kzoverkh, npor, topboundary, topres + ) if topboundary == "semi": z = np.hstack((z[0] + topthick, z)) super().__init__( @@ -1077,6 +1082,7 @@ def __init__( layers=layers, res=res, ) + self.kzoverkh = kzoverkh # add kzoverkh to inhomogeneity object class AreaSinkInhom(Element): From 16747f19d422156d276af376907c511b9fca55a1 Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Wed, 20 May 2026 18:20:04 +0200 Subject: [PATCH 3/4] long line to print to readable string --- timflow/steady/inhomogeneity.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/timflow/steady/inhomogeneity.py b/timflow/steady/inhomogeneity.py index e20d93d..2f1dbbb 100644 --- a/timflow/steady/inhomogeneity.py +++ b/timflow/steady/inhomogeneity.py @@ -424,7 +424,7 @@ def __repr__(self): "BuildingPit: layers " + str(list(self.layers)) + ", " - + str(list(self.x, self.y)) + + str(list(zip(self.x.tolist(), self.y.tolist()))) ) def isinside(self, x, y): @@ -815,7 +815,7 @@ def __repr__(self): "LeakyBuildingPit: layers " + str(list(self.layers)) + ", " - + str(list(self.x, self.y)) + + str(list(zip(self.x.tolist(), self.y.tolist()))) ) def create_elements(self): From 84fd0a27ec8f484fcc32800d84726089dd8fb6a8 Mon Sep 17 00:00:00 2001 From: Mark Bakker Date: Wed, 20 May 2026 18:28:12 +0200 Subject: [PATCH 4/4] How dumb is that? You must supply a strict=False part. --- timflow/steady/inhomogeneity.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/timflow/steady/inhomogeneity.py b/timflow/steady/inhomogeneity.py index 2f1dbbb..a9021cc 100644 --- a/timflow/steady/inhomogeneity.py +++ b/timflow/steady/inhomogeneity.py @@ -424,7 +424,7 @@ def __repr__(self): "BuildingPit: layers " + str(list(self.layers)) + ", " - + str(list(zip(self.x.tolist(), self.y.tolist()))) + + str(list(zip(self.x.tolist(), self.y.tolist(), strict=False))) ) def isinside(self, x, y): @@ -815,7 +815,7 @@ def __repr__(self): "LeakyBuildingPit: layers " + str(list(self.layers)) + ", " - + str(list(zip(self.x.tolist(), self.y.tolist()))) + + str(list(zip(self.x.tolist(), self.y.tolist(), strict=False))) ) def create_elements(self):