From d346f2ff37545b0a419cc85bdacc7e4e7ebab67b Mon Sep 17 00:00:00 2001 From: Adrian Croucher Date: Thu, 4 Apr 2024 16:40:19 +1300 Subject: [PATCH] Modify tolerance checking in column_track() --- mulgrids.py | 10 ++++------ 1 file changed, 4 insertions(+), 6 deletions(-) diff --git a/mulgrids.py b/mulgrids.py index fe128e4..03bdc17 100755 --- a/mulgrids.py +++ b/mulgrids.py @@ -2003,12 +2003,10 @@ def column_track(self, line): line. """ - dl = max(norm(line[1] - line[0]), 1) tol = 1e-3 - - def track_dist(p): - """Non-dimensionalised distance of point along track""" - return norm(p - line[0]) / dl + def track_dist(p): return norm(p - line[0]) + dl = track_dist(line[1]) + tol_dl = tol * dl track, dist = [], [] start_col, end_col = None, None @@ -2034,7 +2032,7 @@ def track_dist(p): elif col == end_col: pts = [pts[0], line[-1]] din, dout = track_dist(pts[0]), track_dist(pts[-1]) - if abs(dout - din) > tol: + if abs(dout - din) > tol_dl: track.append((col, pts[0], pts[-1])) dist.append(din)