diff -r b2914aaeec1a -r e34d6a30b96d src/geometry.cpp
--- a/src/geometry.cpp	Mon Apr 10 14:25:19 2023 +0300
+++ b/src/geometry.cpp	Mon Apr 10 14:46:36 2023 +0300
@@ -256,43 +256,3 @@
 	}
 	return (sum < 0) ? winding_e::anticlockwise : winding_e::clockwise;
 }
-
-/**
- * @brief computes the point on a Bezier curve
- * @param curve
- * @param t scalar between 0 and 1, with t=0 being P0 and t=1 being P3
- * @return point on curve
- */
-glm::vec3 pointOnCurve(const BezierCurve &curve, float t)
-{
-	// clamp t as rounding errors might make it slightly out of bounds
-	t = std::clamp(t, 0.0f, 1.0f);
-	const float t_2 = t * t;
-	const float t_3 = t * t * t;
-	const float coeffs[3] = {
-		-1*t_3  +3*t_2  -3*t  +1,
-		+3*t_3  -6*t_2  +3*t,
-		-3*t_3  +3*t_2,
-	};
-	return coeffs[0] * curve[0] + coeffs[1] * curve[1] + coeffs[2] * curve[2] + t_3 * curve[3];
-}
-
-/**
- * @brief computes the derivative of a point on a Bezier curve
- * @param curve
- * @param t scalar between 0 and 1, with t=0 being P0 and t=1 being P3
- * @return point on curve
- */
-glm::vec3 derivativeOnCurve(const BezierCurve &curve, float t)
-{
-	// clamp t as rounding errors might make it slightly out of bounds
-	t = std::clamp(t, 0.0f, 1.0f);
-	const float t_2 = t * t;
-	const float coeffs[4] = {
-		-3*t_2  + 6*t  -3,
-		+9*t_2  -12*t  +3,
-		-9*t_2  + 6*t,
-		+3*t_2
-	};
-	return coeffs[0] * curve[0] + coeffs[1] * curve[1] + coeffs[2] * curve[2] + coeffs[3] * curve[3];
-}