From f507c0718d919a669613c6ab31d8e1c238f08988 Mon Sep 17 00:00:00 2001
From: LunaTheFoxgirl <cli-luna@protonmail.com>
Date: Sun, 23 Oct 2022 18:12:35 +0200
Subject: [PATCH] Cubic interpolation

---
 inmath/interpolate.d | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/inmath/interpolate.d b/inmath/interpolate.d
index 9385a77..a76243c 100644
--- a/inmath/interpolate.d
+++ b/inmath/interpolate.d
@@ -151,4 +151,14 @@ T hermite(T)(T x, T tx, T y, T ty, float t) {
     float h3 = t^^3 - 2 * t^^2 + t;
     float h4 = t^^3 - t^^2;
     return h1 * x + h3 * tx + h2 * y + h4 * ty;
+}
+
+/// Cubic interpolation
+T cubic(T)(T p0, T p1, T p2, T p3, float t) {
+    T a = -0.5 * p0 + 1.5 * p1 - 1.5 * p2 + 0.5 * p3;
+    T b = p0 - 2.5 * p1 + 2 * p2 - 0.5 * p3;
+    T c = -0.5 * p0 + 0.5 * p2;
+    T d = p1;
+    
+    return a * (t ^^ 3) + b * (t ^^ 2) + c * t + d;
 }
\ No newline at end of file