patched edge-case in LinearInterpolation

Author: packquickly

diffrax/global_interpolation.py

@@ -125,10 +125,11 @@ def _index(_ys):
         prev_t = self.ts[index]
         next_t = self.ts[index + 1]
         diff_t = next_t - prev_t
-
-        return (
-            prev_ys**ω + (next_ys**ω - prev_ys**ω) * (fractional_part / diff_t)
-        ).ω
+        return jnp.where(
+            diff_t >= jnp.finfo(diff_t.dtype).eps,
+            (prev_ys**ω + (next_ys**ω - prev_ys**ω) * (fractional_part / diff_t)).ω,
+            prev_ys
+        )
 
     @eqx.filter_jit
     def derivative(self, t: Scalar, left: bool = True) -> PyTree: