Run rustfmt.

Update glam requirement from >=0.10, <0.23 to >=0.10, <0.25

CHANGELOG.md

@@ -1,7 +1,7 @@
 # Changelog
 
-<!-- vim-markdown-toc GFM -->
+* [4.1.2](#412)
-
+* [4.1.1](#411)
 * [4.1](#41)
 * [4.0.3](#403)
 * [4.0.2](#402)
@@ -41,7 +41,25 @@
 * [0.1.1](#011)
 * [0.1](#01)
 
-<!-- vim-markdown-toc -->
+# 4.2.0
+
+> Feb 1, 2023
+
+- Add support for `glam-0.22`.
+- Add support for `nalgebra-0.32`.
+- Add deprecation lints for `impl-*` feature gates. Those shouldn’t be used anymore and the `*` variant should be 
+  preferred. For instance, if you used `impl-cgmath`, you should just use the `cgmath` feature gate now.
+
+# 4.1.1
+
+> Jul 27, 2022
+
+- Internal enhancement of sampling splines by looking for control points. That brings the lookup from _O(N)_ to
+  _O(log(N))_. That is super embarassing because it should have been the default from the very first commit. Sorry
+  about that.
+- Fix hermite cubic interpolation.
+- Add support for `glam-0.21`.
+- Add support for `nalgebra-0.31`.
 
 # 4.1
 

Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "splines"
-version = "4.1.0"
+version = "4.2.0"
 license = "BSD-3-Clause"
 authors = ["Dimitri Sabadie <dimitri.sabadie@gmail.com>"]
 description = "Spline interpolation made easy"
@@ -13,12 +13,6 @@ readme = "README.md"
 
 edition = "2021"
 
-[badges]
-travis-ci = { repository = "phaazon/splines", branch = "master" }
-is-it-maintained-issue-resolution = { repository = "phaazon/splines" }
-is-it-maintained-open-issues = { repository = "phaazon/splines" }
-maintenance = { status = "actively-developed" }
-
 [features]
 default = ["std"]
 impl-cgmath = ["cgmath"]
@@ -29,8 +23,8 @@ std = []
 
 [dependencies]
 cgmath = { version = ">=0.17, <0.19", optional = true }
-glam = { version = ">=0.10, <0.21", optional = true }
+glam = { version = ">=0.10, <0.25", optional = true }
-nalgebra = { version = ">=0.21, <0.31", optional = true }
+nalgebra = { version = ">=0.21, <0.33", optional = true }
 serde =  { version = "1", features = ["derive"], optional = true }
 
 [dev-dependencies]
@@ -38,7 +32,7 @@ float-cmp = ">=0.6, < 0.10"
 serde_json = "1"
 
 [package.metadata.docs.rs]
-all-features = true
+features = ["std", "cgmath", "glam", "nalgebra", "serde"]
 
 [[example]]
 name = "hello-world"

src/interpolate.rs

@@ -127,15 +127,16 @@ macro_rules! impl_Interpolate {
         let three_t = t * 3.;
         let t2 = t * t;
         let t3 = t2 * t;
-        let two_t3 = t3 * two_t;
+        let two_t3 = t2 * two_t;
-        let three_t2 = t2 * three_t;
+        let two_t2 = t * two_t;
+        let three_t2 = t * three_t;
 
         // tangents
-        let m0 = (b.1 - x.1) / (b.0 - x.0);
+        let m0 = (b.1 - x.1) / (b.0 - x.0) * (b.0 - a.0);
-        let m1 = (y.1 - a.1) / (y.0 - a.0);
+        let m1 = (y.1 - a.1) / (y.0 - a.0) * (b.0 - a.0);
 
         a.1 * (two_t3 - three_t2 + 1.)
-          + m0 * (t3 - t2 * two_t + t)
+          + m0 * (t3 - two_t2 + t)
           + b.1 * (three_t2 - two_t3)
           + m1 * (t3 - t2)
       }
@@ -192,15 +193,16 @@ macro_rules! impl_InterpolateT {
         let three_t = t * 3.;
         let t2 = t * t;
         let t3 = t2 * t;
-        let two_t3 = t3 * two_t;
+        let two_t3 = t2 * two_t;
-        let three_t2 = t2 * three_t;
+        let two_t2 = t * two_t;
+        let three_t2 = t * three_t;
 
         // tangents
-        let m0 = (b.1 - x.1) / (Self::from(b.0 - x.0));
+        let m0 = (b.1 - x.1) / (Self::from(b.0 - x.0)) * (Self::from(b.0 - a.0));
-        let m1 = (y.1 - a.1) / (Self::from(y.0 - a.0));
+        let m1 = (y.1 - a.1) / (Self::from(y.0 - a.0)) * (Self::from(b.0 - a.0));
 
         a.1 * (two_t3 - three_t2 + 1.)
-          + m0 * (t3 - t2 * two_t + t)
+          + m0 * (t3 - two_t2 + t)
           + b.1 * (three_t2 - two_t3)
           + m1 * (t3 - t2)
       }

src/lib.rs

@@ -108,6 +108,17 @@
 #![cfg_attr(not(feature = "std"), no_std)]
 #![cfg_attr(not(feature = "std"), feature(alloc))]
 #![cfg_attr(not(feature = "std"), feature(core_intrinsics))]
+#![cfg_attr(
+  any(
+    feature = "impl-cgmath",
+    feature = "impl-glam",
+    feature = "impl-nalgebra"
+  ),
+  deprecated(
+    since = "4.2.0",
+    note = "you are using an impl-* feature gate; please switch to * (e.g. impl-cgmath becomes cgmath)"
+  )
+)]
 
 #[cfg(not(feature = "std"))]
 extern crate alloc;

src/spline.rs

@@ -56,6 +56,13 @@ impl<T, V> Spline<T, V> {
     spline
   }
 
+  /// Clear the spline by removing all keys. Keeps the underlying allocated storage, so adding
+  /// new keys should be faster than creating a new [`Spline`]
+  #[inline]
+  pub fn clear(&mut self) {
+    self.0.clear()
+  }
+
   /// Create a new spline by consuming an `Iterater<Item = Key<T>>`. They keys don’t have to be
   /// sorted.
   ///
@@ -318,37 +325,20 @@ pub struct KeyMut<'a, T, V> {
 }
 
 // Find the lower control point corresponding to a given time.
+// It has the property to have a timestamp smaller or equal to t
 fn search_lower_cp<T, V>(cps: &[Key<T, V>], t: T) -> Option<usize>
 where
   T: PartialOrd,
 {
-  let mut i = 0;
   let len = cps.len();
-
   if len < 2 {
     return None;
   }
-
+  match cps.binary_search_by(|key| key.t.partial_cmp(&t).unwrap()) {
-  loop {
+    Err(i) if i >= len => None,
-    let cp = &cps[i];
+    Err(i) if i == 0 => None,
-    let cp1 = &cps[i + 1];
+    Err(i) => Some(i - 1),
-
+    Ok(i) if i == len - 1 => None,
-    if t >= cp1.t {
+    Ok(i) => Some(i),
-      if i >= len - 2 {
-        return None;
-      }
-
-      i += 1;
-    } else if t < cp.t {
-      if i == 0 {
-        return None;
-      }
-
-      i -= 1;
-    } else {
-      break; // found
-    }
   }
-
-  Some(i)
 }