From 0f9739c0e1e5ca451401ad40e8a242d0bf457460 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Fran=C3=A7ois=20LAURENT?= <francois.laurent@pasteur.fr>
Date: Wed, 12 Oct 2022 20:56:32 +0200
Subject: [PATCH] bug fix for time intervals with infinite number of decimals

---
 src/taggingbackends/features/skeleton.py | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/taggingbackends/features/skeleton.py b/src/taggingbackends/features/skeleton.py
index 959094a..e4e9539 100644
--- a/src/taggingbackends/features/skeleton.py
+++ b/src/taggingbackends/features/skeleton.py
@@ -48,10 +48,10 @@ def interpolate(times, spines, anchor, window_length,
     grid = range(istart, istop+1)
     series = []
     for i in grid:
-        t = round((anchor + i * frame_interval) * 1e4) * 1e-4
+        t = np.trunc((anchor + i * frame_interval) * 1e4) * 1e-4
         inext = np.flatnonzero(t <= times)[0]
         tnext, xnext = times[inext], spines[inext]
-        if tnext == t:
+        if tnext - t < 1e-4:
             x = xnext
         else:
             assert 0 < inext
-- 
GitLab