curvy 1.0.0 → 1.0.2

package/dist/curve/cubic2d.d.mts

@@ -1,7 +1,7 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.mjs';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Interval } from '../interval.mjs';
-import { C as CubicPolynomial } from '../shared/curvy.Bxxwly94.mjs';
+import { C as CubicPolynomial } from '../shared/curvy.DwHW5Fqt.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
@@ -23,5 +23,9 @@ declare const length: {
     (c: CubicCurve2d, i: Interval): number;
     (i: Interval): (c: CubicCurve2d) => number;
 };
+declare const curvature: {
+    (c: CubicCurve2d, t: number): number;
+    (t: number): (c: CubicCurve2d) => number;
+};
 
-export { type CubicCurve2d, fromPoints, fromPolynomials, isCubicCurve2d, length, solve };
+export { type CubicCurve2d, curvature, fromPoints, fromPolynomials, isCubicCurve2d, length, solve };

package/dist/curve/cubic2d.d.ts

@@ -1,7 +1,7 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.js';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Interval } from '../interval.js';
-import { C as CubicPolynomial } from '../shared/curvy.CO0Wclt0.js';
+import { C as CubicPolynomial } from '../shared/curvy.CkJRVnqj.js';
 import { Vector2 } from '../vector/vector2.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
@@ -23,5 +23,9 @@ declare const length: {
     (c: CubicCurve2d, i: Interval): number;
     (i: Interval): (c: CubicCurve2d) => number;
 };
+declare const curvature: {
+    (c: CubicCurve2d, t: number): number;
+    (t: number): (c: CubicCurve2d) => number;
+};
 
-export { type CubicCurve2d, fromPoints, fromPolynomials, isCubicCurve2d, length, solve };
+export { type CubicCurve2d, curvature, fromPoints, fromPolynomials, isCubicCurve2d, length, solve };

package/dist/curve/cubic2d.mjs

@@ -1,15 +1,15 @@
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
 import { size, scaleShift, biunit } from '../interval.mjs';
 import { G as GL32_W0, a as GL32_X0, b as GL32_W1, c as GL32_X1, d as GL32_W2, e as GL32_X2, f as GL32_W3, g as GL32_X3, h as GL32_W4, i as GL32_X4, j as GL32_W5, k as GL32_X5, l as GL32_W6, m as GL32_X6, n as GL32_W7, o as GL32_X7, p as GL32_W8, q as GL32_X8, r as GL32_W9, s as GL32_X9, t as GL32_W10, u as GL32_X10, v as GL32_W11, w as GL32_X11, x as GL32_W12, y as GL32_X12, z as GL32_W13, A as GL32_X13, B as GL32_W14, C as GL32_X14, D as GL32_W15, E as GL32_X15 } from '../shared/curvy.DmVFTi0N.mjs';
-import { make, solve as solve$2, derivative as derivative$1 } from '../polynomial/cubic.mjs';
+import { make, solve as solve$2, derivative as derivative$2 } from '../polynomial/cubic.mjs';
 import { make as make$2 } from '../polynomial/quadratic.mjs';
-import { make as make$1, magnitude } from '../vector/vector2.mjs';
-import { length as length$2, fromPolynomials as fromPolynomials$2, solve as solve$3 } from './quadratic2d.mjs';
-import '../shared/curvy.CLzGuqek.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import { make as make$1, magnitude, cross } from '../vector/vector2.mjs';
+import { solve as solve$4 } from './linear2d.mjs';
+import { length as length$2, fromPolynomials as fromPolynomials$2, solve as solve$3, derivative as derivative$1 } from './quadratic2d.mjs';
+import '../shared/curvy.ZjG-q57I.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 import '../polynomial/linear.mjs';
-import './linear2d.mjs';
 
 const CubicCurve2dTypeId = Symbol("curvy/curve/cubic2d");
 class CubicCurve2dImpl extends Pipeable {
@@ -39,8 +39,8 @@ const solve$1 = dual(
   (c, t) => make$1(solve$2(c.x, t), solve$2(c.y, t))
 );
 const derivative = (c) => fromPolynomials$2(
-  derivative$1(c.x),
-  derivative$1(c.y)
+  derivative$2(c.x),
+  derivative$2(c.y)
 );
 const length$1 = dual(2, (c, i) => {
   if (size(i) === 0) {
@@ -75,11 +75,22 @@ const length$1 = dual(2, (c, i) => {
     )) * scale
   );
 });
+const curvature$1 = dual(2, (c, t) => {
+  const d = derivative(c);
+  const v = solve$3(d, t);
+  const a = solve$4(derivative$1(d), t);
+  const vMag = magnitude(v);
+  if (vMag === 0) {
+    return 0;
+  }
+  return round(Math.abs(cross(v, a)) / vMag ** 3);
+});
 
 const fromPolynomials = fromPolynomials$1;
 const fromPoints = fromPoints$1;
 const isCubicCurve2d = isCubicCurve2d$1;
 const solve = solve$1;
 const length = length$1;
+const curvature = curvature$1;
 
-export { fromPoints, fromPolynomials, isCubicCurve2d, length, solve };
+export { curvature, fromPoints, fromPolynomials, isCubicCurve2d, length, solve };

package/dist/curve/linear2d.d.mts

@@ -1,7 +1,7 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.mjs';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Interval } from '../interval.mjs';
-import { L as LinearPolynomial } from '../shared/curvy.Bxxwly94.mjs';
+import { L as LinearPolynomial } from '../shared/curvy.DwHW5Fqt.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
@@ -15,6 +15,7 @@ interface LinearCurve2d extends Pipeable, TwoDimensional<LinearPolynomial> {
 declare const fromPolynomials: (x: LinearPolynomial, y: LinearPolynomial) => LinearCurve2d;
 declare const fromPoints: (p0: Vector2, p1: Vector2) => LinearCurve2d;
 declare const isLinearCurve2d: (c: unknown) => c is LinearCurve2d;
+declare const derivative: (c: LinearCurve2d) => Vector2;
 declare const solve: {
     (t: number): (c: LinearCurve2d) => Vector2;
     (c: LinearCurve2d, t: number): Vector2;
@@ -24,4 +25,4 @@ declare const length: {
     (i: Interval): (c: LinearCurve2d) => number;
 };
 
-export { type LinearCurve2d, fromPoints, fromPolynomials, isLinearCurve2d, length, solve };
+export { type LinearCurve2d, derivative, fromPoints, fromPolynomials, isLinearCurve2d, length, solve };

package/dist/curve/linear2d.d.ts

@@ -1,7 +1,7 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.js';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Interval } from '../interval.js';
-import { L as LinearPolynomial } from '../shared/curvy.CO0Wclt0.js';
+import { L as LinearPolynomial } from '../shared/curvy.CkJRVnqj.js';
 import { Vector2 } from '../vector/vector2.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
@@ -15,6 +15,7 @@ interface LinearCurve2d extends Pipeable, TwoDimensional<LinearPolynomial> {
 declare const fromPolynomials: (x: LinearPolynomial, y: LinearPolynomial) => LinearCurve2d;
 declare const fromPoints: (p0: Vector2, p1: Vector2) => LinearCurve2d;
 declare const isLinearCurve2d: (c: unknown) => c is LinearCurve2d;
+declare const derivative: (c: LinearCurve2d) => Vector2;
 declare const solve: {
     (t: number): (c: LinearCurve2d) => Vector2;
     (c: LinearCurve2d, t: number): Vector2;
@@ -24,4 +25,4 @@ declare const length: {
     (i: Interval): (c: LinearCurve2d) => number;
 };
 
-export { type LinearCurve2d, fromPoints, fromPolynomials, isLinearCurve2d, length, solve };
+export { type LinearCurve2d, derivative, fromPoints, fromPolynomials, isLinearCurve2d, length, solve };

package/dist/curve/linear2d.mjs

@@ -1,9 +1,9 @@
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
 import { make, solve as solve$2 } from '../polynomial/linear.mjs';
 import { make as make$1 } from '../vector/vector2.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
 import '../interval.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 
 const LinearCurve2dTypeId = Symbol("curvy/curve/linear2d");
 class LinearCurve2dImpl extends Pipeable {
@@ -28,6 +28,7 @@ const fromPoints$1 = (p0, p1) => new LinearCurve2dImpl(
   make(p0.y, p1.y)
 );
 const isLinearCurve2d$1 = (c) => typeof c === "object" && c !== null && LinearCurve2dTypeId in c;
+const derivative$1 = (c) => make$1(c.x.c1, c.y.c1);
 const solve$1 = dual(
   2,
   (c, t) => make$1(solve$2(c.x, t), solve$2(c.y, t))
@@ -40,7 +41,8 @@ const length$1 = dual(
 const fromPolynomials = fromPolynomials$1;
 const fromPoints = fromPoints$1;
 const isLinearCurve2d = isLinearCurve2d$1;
+const derivative = derivative$1;
 const solve = solve$1;
 const length = length$1;
 
-export { fromPoints, fromPolynomials, isLinearCurve2d, length, solve };
+export { derivative, fromPoints, fromPolynomials, isLinearCurve2d, length, solve };

package/dist/curve/quadratic2d.d.mts

@@ -1,8 +1,9 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.mjs';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Interval } from '../interval.mjs';
-import { Q as QuadraticPolynomial } from '../shared/curvy.Bxxwly94.mjs';
+import { Q as QuadraticPolynomial } from '../shared/curvy.DwHW5Fqt.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
+import { LinearCurve2d } from './linear2d.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
 
@@ -23,5 +24,10 @@ declare const length: {
     (c: QuadraticCurve2d, i: Interval): number;
     (i: Interval): (c: QuadraticCurve2d) => number;
 };
+declare const derivative: (c: QuadraticCurve2d) => LinearCurve2d;
+declare const curvature: {
+    (c: QuadraticCurve2d, t: number): number;
+    (t: number): (c: QuadraticCurve2d) => number;
+};
 
-export { type QuadraticCurve2d, fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };
+export { type QuadraticCurve2d, curvature, derivative, fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };

package/dist/curve/quadratic2d.d.ts

@@ -1,8 +1,9 @@
-import { e as TwoDimensional } from '../shared/curvy.CT38g_II.js';
+import { T as TwoDimensional } from '../shared/curvy.Bf5TIk31.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Interval } from '../interval.js';
-import { Q as QuadraticPolynomial } from '../shared/curvy.CO0Wclt0.js';
+import { Q as QuadraticPolynomial } from '../shared/curvy.CkJRVnqj.js';
 import { Vector2 } from '../vector/vector2.js';
+import { LinearCurve2d } from './linear2d.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
 
@@ -23,5 +24,10 @@ declare const length: {
     (c: QuadraticCurve2d, i: Interval): number;
     (i: Interval): (c: QuadraticCurve2d) => number;
 };
+declare const derivative: (c: QuadraticCurve2d) => LinearCurve2d;
+declare const curvature: {
+    (c: QuadraticCurve2d, t: number): number;
+    (t: number): (c: QuadraticCurve2d) => number;
+};
 
-export { type QuadraticCurve2d, fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };
+export { type QuadraticCurve2d, curvature, derivative, fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };

package/dist/curve/quadratic2d.mjs

@@ -1,13 +1,13 @@
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
 import { size, scaleShift, biunit } from '../interval.mjs';
 import { G as GL32_W0, a as GL32_X0, b as GL32_W1, c as GL32_X1, d as GL32_W2, e as GL32_X2, f as GL32_W3, g as GL32_X3, h as GL32_W4, i as GL32_X4, j as GL32_W5, k as GL32_X5, l as GL32_W6, m as GL32_X6, n as GL32_W7, o as GL32_X7, p as GL32_W8, q as GL32_X8, r as GL32_W9, s as GL32_X9, t as GL32_W10, u as GL32_X10, v as GL32_W11, w as GL32_X11, x as GL32_W12, y as GL32_X12, z as GL32_W13, A as GL32_X13, B as GL32_W14, C as GL32_X14, D as GL32_W15, E as GL32_X15 } from '../shared/curvy.DmVFTi0N.mjs';
 import { make as make$2 } from '../polynomial/linear.mjs';
-import { make, solve as solve$2, derivative as derivative$1 } from '../polynomial/quadratic.mjs';
-import { make as make$1, magnitude } from '../vector/vector2.mjs';
-import { length as length$2, fromPolynomials as fromPolynomials$2, solve as solve$3 } from './linear2d.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
-import '../shared/curvy.CLzGuqek.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import { make, solve as solve$2, derivative as derivative$2 } from '../polynomial/quadratic.mjs';
+import { make as make$1, magnitude, cross } from '../vector/vector2.mjs';
+import { fromPolynomials as fromPolynomials$2, length as length$2, solve as solve$3, derivative as derivative$3 } from './linear2d.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
+import '../shared/curvy.ZjG-q57I.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 
 const QuadraticCurve2dTypeId = Symbol(
   "curvy/curve/quadratic2d"
@@ -41,9 +41,9 @@ const solve$1 = dual(
     solve$2(c.y, t)
   )
 );
-const derivative = (c) => fromPolynomials$2(
-  derivative$1(c.x),
-  derivative$1(c.y)
+const derivative$1 = (c) => fromPolynomials$2(
+  derivative$2(c.x),
+  derivative$2(c.y)
 );
 const length$1 = dual(2, (c, i) => {
   if (size(i) === 0) {
@@ -58,17 +58,29 @@ const length$1 = dual(2, (c, i) => {
       i
     );
   }
-  const d = derivative(c);
+  const d = derivative$1(c);
   const { scale, shift } = scaleShift(biunit, i);
   return round(
     (GL32_W0 * magnitude(solve$3(d, shift + scale * -GL32_X0)) + GL32_W0 * magnitude(solve$3(d, shift + scale * GL32_X0)) + GL32_W1 * magnitude(solve$3(d, shift + scale * -GL32_X1)) + GL32_W1 * magnitude(solve$3(d, shift + scale * GL32_X1)) + GL32_W2 * magnitude(solve$3(d, shift + scale * -GL32_X2)) + GL32_W2 * magnitude(solve$3(d, shift + scale * GL32_X2)) + GL32_W3 * magnitude(solve$3(d, shift + scale * -GL32_X3)) + GL32_W3 * magnitude(solve$3(d, shift + scale * GL32_X3)) + GL32_W4 * magnitude(solve$3(d, shift + scale * -GL32_X4)) + GL32_W4 * magnitude(solve$3(d, shift + scale * GL32_X4)) + GL32_W5 * magnitude(solve$3(d, shift + scale * -GL32_X5)) + GL32_W5 * magnitude(solve$3(d, shift + scale * GL32_X5)) + GL32_W6 * magnitude(solve$3(d, shift + scale * -GL32_X6)) + GL32_W6 * magnitude(solve$3(d, shift + scale * GL32_X6)) + GL32_W7 * magnitude(solve$3(d, shift + scale * -GL32_X7)) + GL32_W7 * magnitude(solve$3(d, shift + scale * GL32_X7)) + GL32_W8 * magnitude(solve$3(d, shift + scale * -GL32_X8)) + GL32_W8 * magnitude(solve$3(d, shift + scale * GL32_X8)) + GL32_W9 * magnitude(solve$3(d, shift + scale * -GL32_X9)) + GL32_W9 * magnitude(solve$3(d, shift + scale * GL32_X9)) + GL32_W10 * magnitude(solve$3(d, shift + scale * -GL32_X10)) + GL32_W10 * magnitude(solve$3(d, shift + scale * GL32_X10)) + GL32_W11 * magnitude(solve$3(d, shift + scale * -GL32_X11)) + GL32_W11 * magnitude(solve$3(d, shift + scale * GL32_X11)) + GL32_W12 * magnitude(solve$3(d, shift + scale * -GL32_X12)) + GL32_W12 * magnitude(solve$3(d, shift + scale * GL32_X12)) + GL32_W13 * magnitude(solve$3(d, shift + scale * -GL32_X13)) + GL32_W13 * magnitude(solve$3(d, shift + scale * GL32_X13)) + GL32_W14 * magnitude(solve$3(d, shift + scale * -GL32_X14)) + GL32_W14 * magnitude(solve$3(d, shift + scale * GL32_X14)) + GL32_W15 * magnitude(solve$3(d, shift + scale * -GL32_X15)) + GL32_W15 * magnitude(solve$3(d, shift + scale * GL32_X15))) * scale
   );
 });
+const curvature$1 = dual(2, (c, t) => {
+  const d = derivative$1(c);
+  const v = solve$3(d, t);
+  const a = derivative$3(d);
+  const vMag = magnitude(v);
+  if (vMag === 0) {
+    return 0;
+  }
+  return round(Math.abs(cross(v, a)) / vMag ** 3);
+});
 
 const fromPolynomials = fromPolynomials$1;
 const fromPoints = fromPoints$1;
 const isQuadraticCurve2d = isQuadraticCurve2d$1;
 const solve = solve$1;
 const length = length$1;
+const derivative = derivative$1;
+const curvature = curvature$1;
 
-export { fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };
+export { curvature, derivative, fromPoints, fromPolynomials, isQuadraticCurve2d, length, solve };

package/dist/interval.mjs

@@ -1,4 +1,4 @@
-import { d as dual, P as Pipeable, r as round } from './shared/curvy.ChkqgKuv.mjs';
+import { d as dual, P as Pipeable, r as round } from './shared/curvy.DLOeble8.mjs';
 
 const TypeId = Symbol.for("curvy/interval");
 class IntervalImpl extends Pipeable {

package/dist/matrix/matrix2x2.d.mts

@@ -1,4 +1,4 @@
-import { c as TwoDimensionalIndex, d as TwoDimensionalComponent } from '../shared/curvy.CT38g_II.mjs';
+import { d as TwoDimensionalIndex, e as TwoDimensionalComponent } from '../shared/curvy.Bf5TIk31.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
 

package/dist/matrix/matrix2x2.d.ts

@@ -1,4 +1,4 @@
-import { c as TwoDimensionalIndex, d as TwoDimensionalComponent } from '../shared/curvy.CT38g_II.js';
+import { d as TwoDimensionalIndex, e as TwoDimensionalComponent } from '../shared/curvy.Bf5TIk31.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Vector2 } from '../vector/vector2.js';
 

package/dist/matrix/matrix2x2.mjs

@@ -1,6 +1,6 @@
-import { i as isMatrix2x2$1, m as make, f as fromRows$1, a as fromColumns$1, s as setRow$1, b as setColumn$1, d as determinant$1, v as vectorProductLeft$1, c as vectorProductRight$1, e as solveSystem$1, t as toRows$1, g as toColumns$1, r as rowVector$1, h as columnVector$1, j as transpose$1, k as reverseRows$1, l as reverseColumns$1 } from '../shared/curvy.C-sUk7ux.mjs';
-import '../shared/curvy.ChkqgKuv.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import { i as isMatrix2x2$1, m as make, f as fromRows$1, a as fromColumns$1, s as setRow$1, b as setColumn$1, d as determinant$1, v as vectorProductLeft$1, c as vectorProductRight$1, e as solveSystem$1, t as toRows$1, g as toColumns$1, r as rowVector$1, h as columnVector$1, j as transpose$1, k as reverseRows$1, l as reverseColumns$1 } from '../shared/curvy.BrrC9QVB.mjs';
+import '../shared/curvy.DLOeble8.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 
 const isMatrix2x2 = isMatrix2x2$1;
 const matrix2x2 = make;

package/dist/matrix/matrix3x3.d.mts

@@ -1,4 +1,4 @@
-import { T as ThreeDimensionalIndex, b as ThreeDimensionalComponent } from '../shared/curvy.CT38g_II.mjs';
+import { b as ThreeDimensionalIndex, c as ThreeDimensionalComponent } from '../shared/curvy.Bf5TIk31.mjs';
 import { Vector3 } from '../vector/vector3.mjs';
 import { Matrix2x2 } from './matrix2x2.mjs';
 import '../shared/curvy.BN-oFBsc.mjs';

package/dist/matrix/matrix3x3.d.ts

@@ -1,4 +1,4 @@
-import { T as ThreeDimensionalIndex, b as ThreeDimensionalComponent } from '../shared/curvy.CT38g_II.js';
+import { b as ThreeDimensionalIndex, c as ThreeDimensionalComponent } from '../shared/curvy.Bf5TIk31.js';
 import { Vector3 } from '../vector/vector3.js';
 import { Matrix2x2 } from './matrix2x2.js';
 import '../shared/curvy.BN-oFBsc.js';

package/dist/matrix/matrix3x3.mjs

@@ -1,8 +1,8 @@
-import { i as isMatrix3x3$1, m as make, f as fromRows$1, a as fromColumns$1, s as setRow$1, b as setColumn$1, d as determinant$1, c as minor$1, v as vectorProductLeft$1, e as vectorProductRight$1, g as solveSystem$1, t as toRows$1, h as toColumns$1, r as rowVector$1, j as columnVector$1, k as transpose$1, l as reverseRows$1, n as reverseColumns$1 } from '../shared/curvy.BAW_cKLu.mjs';
-import '../shared/curvy.C-sUk7ux.mjs';
-import '../shared/curvy.ChkqgKuv.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
-import '../shared/curvy.D7bhnwZy.mjs';
+import { i as isMatrix3x3$1, m as make, f as fromRows$1, a as fromColumns$1, s as setRow$1, b as setColumn$1, d as determinant$1, c as minor$1, v as vectorProductLeft$1, e as vectorProductRight$1, g as solveSystem$1, t as toRows$1, h as toColumns$1, r as rowVector$1, j as columnVector$1, k as transpose$1, l as reverseRows$1, n as reverseColumns$1 } from '../shared/curvy.DTco5vku.mjs';
+import '../shared/curvy.BrrC9QVB.mjs';
+import '../shared/curvy.DLOeble8.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
+import '../shared/curvy.DV4S13Cq.mjs';
 
 const isMatrix3x3 = isMatrix3x3$1;
 const matrix3x3 = make;

package/dist/matrix/matrix4x4.d.mts

@@ -1,4 +1,4 @@
-import { F as FourDimensionalIndex, a as FourDimensionalComponent } from '../shared/curvy.CT38g_II.mjs';
+import { F as FourDimensionalIndex, a as FourDimensionalComponent } from '../shared/curvy.Bf5TIk31.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Vector4 } from '../vector/vector4.mjs';
 import { Matrix3x3 } from './matrix3x3.mjs';

package/dist/matrix/matrix4x4.d.ts

@@ -1,4 +1,4 @@
-import { F as FourDimensionalIndex, a as FourDimensionalComponent } from '../shared/curvy.CT38g_II.js';
+import { F as FourDimensionalIndex, a as FourDimensionalComponent } from '../shared/curvy.Bf5TIk31.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Vector4 } from '../vector/vector4.js';
 import { Matrix3x3 } from './matrix3x3.js';

package/dist/matrix/matrix4x4.mjs

@@ -1,9 +1,9 @@
-import { n as toFourDimensionalIndex } from '../shared/curvy.C-sUk7ux.mjs';
-import { d as dual, r as round, i as invariant, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
-import { c as components, m as make$3, d as dot } from '../shared/curvy.Cty3LrOR.mjs';
-import { d as determinant$2, m as make$2 } from '../shared/curvy.BAW_cKLu.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
-import '../shared/curvy.D7bhnwZy.mjs';
+import { n as toFourDimensionalIndex } from '../shared/curvy.BrrC9QVB.mjs';
+import { d as dual, r as round, i as invariant, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
+import { c as components, m as make$3, d as dot } from '../shared/curvy.B1zfBfvf.mjs';
+import { d as determinant$2, m as make$2 } from '../shared/curvy.DTco5vku.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
+import '../shared/curvy.DV4S13Cq.mjs';
 
 const Matrix4x4TypeId = Symbol.for("curvy/matrix4x4");
 class Matrix4x4Impl extends Pipeable {

package/dist/path/cubic2d.d.mts

@@ -1,9 +1,9 @@
 import { CubicCurve2d } from '../curve/cubic2d.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
-import '../shared/curvy.CT38g_II.mjs';
+import '../shared/curvy.Bf5TIk31.mjs';
 import '../interval.mjs';
-import '../shared/curvy.Bxxwly94.mjs';
+import '../shared/curvy.DwHW5Fqt.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
 

package/dist/path/cubic2d.d.ts

@@ -1,9 +1,9 @@
 import { CubicCurve2d } from '../curve/cubic2d.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Vector2 } from '../vector/vector2.js';
-import '../shared/curvy.CT38g_II.js';
+import '../shared/curvy.Bf5TIk31.js';
 import '../interval.js';
-import '../shared/curvy.CO0Wclt0.js';
+import '../shared/curvy.CkJRVnqj.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
 

package/dist/path/cubic2d.mjs

@@ -1,22 +1,23 @@
 import { length as length$2, solve as solve$2 } from '../curve/cubic2d.mjs';
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable, i as invariant } from '../shared/curvy.DLOeble8.mjs';
 import { unit } from '../interval.mjs';
 import '../shared/curvy.DmVFTi0N.mjs';
 import '../polynomial/cubic.mjs';
-import '../shared/curvy.CLzGuqek.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
+import '../shared/curvy.ZjG-q57I.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
 import '../polynomial/quadratic.mjs';
 import '../vector/vector2.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
-import '../curve/quadratic2d.mjs';
-import '../polynomial/linear.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 import '../curve/linear2d.mjs';
+import '../polynomial/linear.mjs';
+import '../curve/quadratic2d.mjs';
 
 const CubicPath2dTypeId = Symbol("curvy/path/cubic2d");
 class CubicPath2dImpl extends Pipeable {
   [CubicPath2dTypeId] = CubicPath2dTypeId;
   curves;
   constructor(curves) {
+    invariant(curves.length > 0, "cubic path requires at least one curve");
     super();
     this.curves = curves;
   }

package/dist/path/linear2d.d.mts

@@ -1,9 +1,9 @@
 import { LinearCurve2d } from '../curve/linear2d.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
-import '../shared/curvy.CT38g_II.mjs';
+import '../shared/curvy.Bf5TIk31.mjs';
 import '../interval.mjs';
-import '../shared/curvy.Bxxwly94.mjs';
+import '../shared/curvy.DwHW5Fqt.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
 

package/dist/path/linear2d.d.ts

@@ -1,9 +1,9 @@
 import { LinearCurve2d } from '../curve/linear2d.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Vector2 } from '../vector/vector2.js';
-import '../shared/curvy.CT38g_II.js';
+import '../shared/curvy.Bf5TIk31.js';
 import '../interval.js';
-import '../shared/curvy.CO0Wclt0.js';
+import '../shared/curvy.CkJRVnqj.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
 

package/dist/path/linear2d.mjs

@@ -1,10 +1,10 @@
 import { length as length$2, solve as solve$2 } from '../curve/linear2d.mjs';
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
 import { unit } from '../interval.mjs';
 import '../polynomial/linear.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
 import '../vector/vector2.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 
 const LinearPath2dTypeId = Symbol("curvy/path/linear2d");
 class LinearPath2dImpl extends Pipeable {

package/dist/path/quadratic2d.d.mts

@@ -1,11 +1,12 @@
 import { QuadraticCurve2d } from '../curve/quadratic2d.mjs';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.mjs';
 import { Vector2 } from '../vector/vector2.mjs';
-import '../shared/curvy.CT38g_II.mjs';
+import '../shared/curvy.Bf5TIk31.mjs';
 import '../interval.mjs';
-import '../shared/curvy.Bxxwly94.mjs';
+import '../shared/curvy.DwHW5Fqt.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';
+import '../curve/linear2d.mjs';
 
 declare const QuadraticPath2dTypeId: unique symbol;
 type QuadraticPath2dTypeId = typeof QuadraticPath2dTypeId;

package/dist/path/quadratic2d.d.ts

@@ -1,11 +1,12 @@
 import { QuadraticCurve2d } from '../curve/quadratic2d.js';
 import { P as Pipeable } from '../shared/curvy.BN-oFBsc.js';
 import { Vector2 } from '../vector/vector2.js';
-import '../shared/curvy.CT38g_II.js';
+import '../shared/curvy.Bf5TIk31.js';
 import '../interval.js';
-import '../shared/curvy.CO0Wclt0.js';
+import '../shared/curvy.CkJRVnqj.js';
 import '../vector/vector3.js';
 import '../vector/vector4.js';
+import '../curve/linear2d.js';
 
 declare const QuadraticPath2dTypeId: unique symbol;
 type QuadraticPath2dTypeId = typeof QuadraticPath2dTypeId;

package/dist/path/quadratic2d.mjs

@@ -1,13 +1,13 @@
 import { length as length$2, solve as solve$2 } from '../curve/quadratic2d.mjs';
-import { d as dual, r as round, P as Pipeable } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round, P as Pipeable } from '../shared/curvy.DLOeble8.mjs';
 import { unit } from '../interval.mjs';
 import '../shared/curvy.DmVFTi0N.mjs';
 import '../polynomial/linear.mjs';
-import '../shared/curvy.DgRsieUW.mjs';
+import '../shared/curvy.O57gIOmy.mjs';
 import '../polynomial/quadratic.mjs';
-import '../shared/curvy.CLzGuqek.mjs';
+import '../shared/curvy.ZjG-q57I.mjs';
 import '../vector/vector2.mjs';
-import '../shared/curvy.DOLqYbKE.mjs';
+import '../shared/curvy.BPOevbwG.mjs';
 import '../curve/linear2d.mjs';
 
 const QuadraticPath2dTypeId = Symbol(

package/dist/polynomial/cubic.d.mts

@@ -1,7 +1,7 @@
 import '../shared/curvy.BN-oFBsc.mjs';
 import '../interval.mjs';
 import '../vector/vector4.mjs';
-export { C as CubicPolynomial, d as derivative, g as domain, e as extrema, f as fromVector, i as isCubicPolynomial, l as length, m as make, c as monotonicity, h as range, r as roots, s as solve, a as solveInverse, b as toInverseSolver, t as toSolver } from '../shared/curvy.Bxxwly94.mjs';
-import '../shared/curvy.CT38g_II.mjs';
+export { C as CubicPolynomial, j as curvature, d as derivative, g as domain, e as extrema, f as fromVector, i as isCubicPolynomial, l as length, m as make, c as monotonicity, h as range, r as roots, s as solve, a as solveInverse, b as toInverseSolver, t as toSolver } from '../shared/curvy.DwHW5Fqt.mjs';
+import '../shared/curvy.Bf5TIk31.mjs';
 import '../vector/vector2.mjs';
 import '../vector/vector3.mjs';

package/dist/polynomial/cubic.d.ts

@@ -1,7 +1,7 @@
 import '../shared/curvy.BN-oFBsc.js';
 import '../interval.js';
 import '../vector/vector4.js';
-export { C as CubicPolynomial, d as derivative, g as domain, e as extrema, f as fromVector, i as isCubicPolynomial, l as length, m as make, c as monotonicity, h as range, r as roots, s as solve, a as solveInverse, b as toInverseSolver, t as toSolver } from '../shared/curvy.CO0Wclt0.js';
-import '../shared/curvy.CT38g_II.js';
+export { C as CubicPolynomial, j as curvature, d as derivative, g as domain, e as extrema, f as fromVector, i as isCubicPolynomial, l as length, m as make, c as monotonicity, h as range, r as roots, s as solve, a as solveInverse, b as toInverseSolver, t as toSolver } from '../shared/curvy.CkJRVnqj.js';
+import '../shared/curvy.Bf5TIk31.js';
 import '../vector/vector2.js';
 import '../vector/vector3.js';

package/dist/polynomial/cubic.mjs

@@ -1,8 +1,8 @@
-import { d as dual, r as round } from '../shared/curvy.ChkqgKuv.mjs';
+import { d as dual, r as round } from '../shared/curvy.DLOeble8.mjs';
 import { size, filter, fromMinMax, scaleShift, biunit } from '../interval.mjs';
 import { G as GL32_W0, a as GL32_X0, b as GL32_W1, c as GL32_X1, d as GL32_W2, e as GL32_X2, f as GL32_W3, g as GL32_X3, h as GL32_W4, i as GL32_X4, j as GL32_W5, k as GL32_X5, l as GL32_W6, m as GL32_X6, n as GL32_W7, o as GL32_X7, p as GL32_W8, q as GL32_X8, r as GL32_W9, s as GL32_X9, t as GL32_W10, u as GL32_X10, v as GL32_W11, w as GL32_X11, x as GL32_W12, y as GL32_X12, z as GL32_W13, A as GL32_X13, B as GL32_W14, C as GL32_X14, D as GL32_W15, E as GL32_X15 } from '../shared/curvy.DmVFTi0N.mjs';
-import { C as CubicPolynomialImpl, a as CubicPolynomialTypeId, s as solveInverse$2, m as make$2, r as roots$2, b as monotonicity$3, g as guaranteedMonotonicityFromComparison, l as length$3, c as solve$2 } from '../shared/curvy.CLzGuqek.mjs';
-import { m as monotonicity$2, a as make$3, l as length$2 } from '../shared/curvy.DgRsieUW.mjs';
+import { C as CubicPolynomialImpl, a as CubicPolynomialTypeId, s as solveInverse$2, m as make$2, r as roots$2, b as monotonicity$3, g as guaranteedMonotonicityFromComparison, l as length$3, c as solve$2, d as derivative$2 } from '../shared/curvy.ZjG-q57I.mjs';
+import { m as monotonicity$2, a as make$3, l as length$2, s as solve$3 } from '../shared/curvy.O57gIOmy.mjs';
 
 const make$1 = (c0 = 0, c1 = 0, c2 = 0, c3 = 0) => new CubicPolynomialImpl(c0, c1, c2, c3);
 const isCubicPolynomial$1 = (v) => typeof v === "object" && v !== null && CubicPolynomialTypeId in v;
@@ -120,6 +120,13 @@ const length$1 = dual(2, (p, domain2) => {
     (GL32_W0 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X0) ** 2) + GL32_W0 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X0) ** 2) + GL32_W1 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X1) ** 2) + GL32_W1 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X1) ** 2) + GL32_W2 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X2) ** 2) + GL32_W2 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X2) ** 2) + GL32_W3 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X3) ** 2) + GL32_W3 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X3) ** 2) + GL32_W4 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X4) ** 2) + GL32_W4 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X4) ** 2) + GL32_W5 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X5) ** 2) + GL32_W5 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X5) ** 2) + GL32_W6 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X6) ** 2) + GL32_W6 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X6) ** 2) + GL32_W7 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X7) ** 2) + GL32_W7 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X7) ** 2) + GL32_W8 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X8) ** 2) + GL32_W8 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X8) ** 2) + GL32_W9 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X9) ** 2) + GL32_W9 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X9) ** 2) + GL32_W10 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X10) ** 2) + GL32_W10 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X10) ** 2) + GL32_W11 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X11) ** 2) + GL32_W11 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X11) ** 2) + GL32_W12 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X12) ** 2) + GL32_W12 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X12) ** 2) + GL32_W13 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X13) ** 2) + GL32_W13 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X13) ** 2) + GL32_W14 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X14) ** 2) + GL32_W14 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X14) ** 2) + GL32_W15 * Math.sqrt(1 + solve$2(d, shift + scale * -GL32_X15) ** 2) + GL32_W15 * Math.sqrt(1 + solve$2(d, shift + scale * GL32_X15) ** 2)) * scale
   );
 });
+const curvature$1 = dual(2, (p, x) => {
+  const d = derivative$1(p);
+  const dd = derivative$2(d);
+  return round(
+    Math.abs(solve$3(dd, x)) / (1 + solve$2(d, x) ** 2) ** 1.5
+  );
+});
 
 const isCubicPolynomial = isCubicPolynomial$1;
 const make = make$1;
@@ -135,5 +142,6 @@ const monotonicity = monotonicity$1;
 const domain = domain$1;
 const range = range$1;
 const length = length$1;
+const curvature = curvature$1;
 
-export { derivative, domain, extrema, fromVector, isCubicPolynomial, length, make, monotonicity, range, roots, solve, solveInverse, toInverseSolver, toSolver };
+export { curvature, derivative, domain, extrema, fromVector, isCubicPolynomial, length, make, monotonicity, range, roots, solve, solveInverse, toInverseSolver, toSolver };

package/dist/polynomial/linear.d.mts

@@ -1,7 +1,7 @@
 import '../shared/curvy.BN-oFBsc.mjs';
 import '../interval.mjs';
 import '../vector/vector2.mjs';
-export { L as LinearPolynomial, O as antiderivative, P as domain, G as fromVector, E as isLinearPolynomial, S as length, F as make, N as monotonicity, R as range, M as root, H as solve, J as solveInverse, K as toInverseSolver, I as toSolver } from '../shared/curvy.Bxxwly94.mjs';
-import '../shared/curvy.CT38g_II.mjs';
+export { L as LinearPolynomial, U as antiderivative, T as derivative, V as domain, J as fromPointSlope, K as fromPoints, I as fromVector, G as isLinearPolynomial, X as length, H as make, S as monotonicity, W as range, R as root, M as solve, O as solveInverse, P as toInverseSolver, N as toSolver } from '../shared/curvy.DwHW5Fqt.mjs';
+import '../shared/curvy.Bf5TIk31.mjs';
 import '../vector/vector3.mjs';
 import '../vector/vector4.mjs';