curvy 2.0.0-alpha.13 → 2.0.0-alpha.15

package/dist/matrix/index.mjs

@@ -1,5 +1,6 @@
 import { t as __exportAll } from "../shared/chunk-25Hzd3z8.mjs";
-import { A as isMatrix2x2$1, B as transpose$2, C as vectorProductLeft$3, D as equals$2, E as determinant$2, F as setColumn$2, H as vectorProductRight$2, I as setRow$2, L as solveSystem$2, M as reverseColumns$2, N as reverseRows$2, O as fromColumns$2, P as rowVector$2, R as toColumns$2, S as transpose$3, T as columnVector$2, V as vectorProductLeft$2, _ as setColumn$3, b as toColumns$3, c as equals$3, d as isMatrix3x3$1, f as make$3, g as rowVector$3, h as reverseRows$3, j as make$2, k as fromRows$2, l as fromColumns$3, m as reverseColumns$3, o as columnVector$3, p as minor$1, r as matrix4x4_exports, s as determinant$3, u as fromRows$3, v as setRow$3, w as vectorProductRight$3, x as toRows$3, y as solveSystem$3, z as toRows$2 } from "../shared/matrix4x4-BWmf4lZM.mjs";
+import { A as fromRows$2, B as toRows$2, C as transpose$3, D as determinant$2, E as columnVector$2, F as rowVector$2, H as vectorProductLeft$2, I as setColumn$2, L as setRow$2, M as make$2, N as reverseColumns$2, O as equals$2, P as reverseRows$2, R as solveSystem$2, S as toRows$3, T as vectorProductRight$3, U as vectorProductRight$2, V as transpose$2, _ as rowVector$3, a as fromColumns$3, b as solveSystem$3, c as inverse$1, d as isMatrix3x3$1, f as make$3, g as reverseRows$3, h as reverseColumns$3, i as equals$3, j as isMatrix2x2$1, k as fromColumns$2, l as inverseUnsafe$1, m as multiply$1, n as columnVector$3, o as fromRows$3, p as minor$1, r as determinant$3, s as identity$1, t as asInvertible$1, u as isInvertible$1, v as setColumn$3, w as vectorProductLeft$3, x as toColumns$3, y as setRow$3, z as toColumns$2 } from "../shared/matrix3x3.internal-Cod-r7H3.mjs";
+import { r as matrix4x4_exports } from "../shared/matrix4x4-BrqgYZCD.mjs";
 //#region src/matrix/matrix2x2.ts
 var matrix2x2_exports = /* @__PURE__ */ __exportAll({
 	columnVector: () => columnVector$1,
@@ -131,14 +132,20 @@ const reverseColumns$1 = reverseColumns$2;
 //#endregion
 //#region src/matrix/matrix3x3.ts
 var matrix3x3_exports = /* @__PURE__ */ __exportAll({
+	asInvertible: () => asInvertible,
 	columnVector: () => columnVector,
 	determinant: () => determinant,
 	equals: () => equals,
 	fromColumns: () => fromColumns,
 	fromRows: () => fromRows,
+	identity: () => identity,
+	inverse: () => inverse,
+	inverseUnsafe: () => inverseUnsafe,
+	isInvertible: () => isInvertible,
 	isMatrix3x3: () => isMatrix3x3,
 	make: () => make,
 	minor: () => minor,
+	multiply: () => multiply,
 	reverseColumns: () => reverseColumns,
 	reverseRows: () => reverseRows,
 	rowVector: () => rowVector,
@@ -261,6 +268,45 @@ const reverseRows = reverseRows$3;
 * @since 1.0.0
 */
 const reverseColumns = reverseColumns$3;
+/**
+* The 3x3 identity matrix.
+*
+* @since 2.0.0
+*/
+const identity = identity$1;
+const multiply = multiply$1;
+const inverse = inverse$1;
+/**
+* Returns the inverse of a `Matrix3x3`. Throws when the matrix is singular
+* (determinant approximately zero). Useful when you've already validated
+* invertibility out-of-band and want the unwrapped matrix directly; otherwise
+* prefer {@link inverse}.
+*
+* @param m - The matrix to invert.
+* @returns The inverse matrix.
+* @throws If `m` is singular.
+* @since 2.0.0
+*/
+const inverseUnsafe = inverseUnsafe$1;
+/**
+* Type-narrowing predicate: refines a `Matrix3x3<T>` to
+* `Matrix3x3<T & Invertible>` when the matrix's determinant is non-zero
+* (within the default absolute tolerance).
+*
+* @param m - The matrix to check.
+* @returns `true` when `m` is invertible.
+* @since 2.0.0
+*/
+const isInvertible = isInvertible$1;
+/**
+* Asserts that the matrix is invertible, throwing on a singular matrix.
+*
+* @param m - The matrix to assert against.
+* @returns The same matrix, typed with the `Invertible` brand.
+* @throws When `m` is singular.
+* @since 2.0.0
+*/
+const asInvertible = asInvertible$1;
 //#endregion
 export { matrix2x2_exports as Matrix2x2, matrix3x3_exports as Matrix3x3, matrix4x4_exports as Matrix4x4 };
 

package/dist/matrix/index.mjs.map

@@ -1 +1 @@
-{"version":3,"file":"index.mjs","names":["internal.isMatrix2x2","equals","internal.equals","make","internal.make","fromRows","internal.fromRows","fromColumns","internal.fromColumns","setRow","internal.setRow","setColumn","internal.setColumn","determinant","internal.determinant","vectorProductLeft","internal.vectorProductLeft","vectorProductRight","internal.vectorProductRight","solveSystem","internal.solveSystem","toRows","internal.toRows","toColumns","internal.toColumns","rowVector","internal.rowVector","columnVector","internal.columnVector","transpose","internal.transpose","reverseRows","internal.reverseRows","reverseColumns","internal.reverseColumns","internal.isMatrix3x3","internal.equals","internal.make","internal.fromRows","internal.fromColumns","internal.setRow","internal.setColumn","internal.determinant","internal.minor","internal.vectorProductLeft","internal.vectorProductRight","internal.solveSystem","internal.toRows","internal.toColumns","internal.rowVector","internal.columnVector","internal.transpose","internal.reverseRows","internal.reverseColumns"],"sources":["../../src/matrix/matrix2x2.ts","../../src/matrix/matrix3x3.ts"],"sourcesContent":["import type { TwoDimensionalComponent, TwoDimensionalIndex } from '../dimensions.ts'\nimport type { Pipeable } from '../utils.ts'\nimport type { Vector2 } from '../vector/vector2.ts'\nimport type { Matrix2x2TypeId } from './matrix2x2.internal.ts'\nimport * as internal from './matrix2x2.internal.ts'\n\n/**\n * A type representing the coordinate of a 2D matrix.\n *\n * @since 1.0.0\n */\nexport type Matrix2x2Coordinate = TwoDimensionalIndex | TwoDimensionalComponent\n\n/**\n * A 2D matrix with 2 rows and 2 columns.\n *\n * All fields are readonly and immutable, and all operations create new instances.\n *\n * @since 1.0.0\n */\nexport interface Matrix2x2 extends Pipeable {\n  readonly [Matrix2x2TypeId]: Matrix2x2TypeId\n\n  /**\n   * The value in the first row and first column.\n   */\n  readonly m00: number\n  /**\n   * The value in the first row and second column.\n   */\n  readonly m01: number\n  /**\n   * The value in the second row and first column.\n   */\n  readonly m10: number\n  /**\n   * The value in the second row and second column.\n   */\n  readonly m11: number\n}\n\n/**\n * Checks if a value is a `Matrix2x2`.\n *\n * @param m - The value to check.\n * @returns `true` if the value is a `Matrix2x2`, `false` otherwise.\n * @since 1.0.0\n */\nexport const isMatrix2x2: (m: unknown) => m is Matrix2x2 = internal.isMatrix2x2\n\nexport const equals: {\n  /**\n   * Checks if two `Matrix2x2` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param a - The first matrix.\n   * @param b - The second matrix.\n   * @returns `true` when each pair of components is within tolerance.\n   * @since 1.1.0\n   */\n  (a: Matrix2x2, b: Matrix2x2): boolean\n  /**\n   * Checks if two `Matrix2x2` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param b - The second matrix.\n   * @returns A function that takes the first matrix and returns the comparison result.\n   * @since 1.1.0\n   */\n  (b: Matrix2x2): (a: Matrix2x2) => boolean\n} = internal.equals\n\n/**\n * Creates a new `Matrix2x2` instance.\n *\n * @param m00 - The value in the first row and first column. Defaults to 0.\n * @param m01 - The value in the first row and second column. Defaults to 0.\n * @param m10 - The value in the second row and first column. Defaults to 0.\n * @param m11 - The value in the second row and second column. Defaults to 0.\n */\nexport const make: (m00?: number, m01?: number, m10?: number, m11?: number) => Matrix2x2 =\n  internal.make\n\n/**\n * Creates a new `Matrix2x2` instance from two vectors representing rows.\n *\n * @param v0 - The first row vector.\n * @param v1 - The second row vector.\n * @returns A new `Matrix2x2` instance.\n * @since 1.0.0\n */\nexport const fromRows: (v0: Vector2, v1: Vector2) => Matrix2x2 = internal.fromRows\n\n/**\n * Creates a new `Matrix2x2` instance from two vectors representing columns.\n *\n * @param v0 - The first column vector.\n * @param v1 - The second column vector.\n * @returns A new `Matrix2x2` instance.\n * @since 1.0.0\n */\nexport const fromColumns: (v0: Vector2, v1: Vector2) => Matrix2x2 = internal.fromColumns\n\nexport const setRow: {\n  /**\n   * Sets the specified row of the matrix to the given vector.\n   *\n   * @param m - The matrix to modify.\n   * @param row - The row index (0 or 1).\n   * @param v - The vector to set the row to.\n   * @returns A new `Matrix2x2` instance with the specified row set to the given vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, row: Matrix2x2Coordinate, v: Vector2): Matrix2x2\n  /**\n   * Sets the specified row of the matrix to the given vector.\n   *\n   * @param row - The row index (0 or 1).\n   * @param v - The vector to set the row to.\n   * @returns A function that takes a matrix and returns a new `Matrix2x2` instance with the specified row set to the given vector.\n   * @since 1.0.0\n   */\n  (row: Matrix2x2Coordinate, v: Vector2): (m: Matrix2x2) => Matrix2x2\n} = internal.setRow\n\nexport const setColumn: {\n  /**\n   * Sets the specified column of the matrix to the given vector.\n   *\n   * @param m - The matrix to modify.\n   * @param column - The column index (0 or 1).\n   * @param v - The vector to set the column to.\n   * @returns A new `Matrix2x2` instance with the specified column set to the given vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, column: Matrix2x2Coordinate, v: Vector2): Matrix2x2\n  /**\n   * Sets the specified column of the matrix to the given vector.\n   *\n   * @param column - The column index (0 or 1).\n   * @param v - The vector to set the column to.\n   * @returns A function that takes a matrix and returns a new `Matrix2x2` instance with the specified column set to the given vector.\n   * @since 1.0.0\n   */\n  (column: Matrix2x2Coordinate, v: Vector2): (m: Matrix2x2) => Matrix2x2\n} = internal.setColumn\n\n/**\n * Calculates the determinant of a 2D matrix.\n *\n * @param m - The matrix to calculate the determinant for.\n * @returns The determinant of the matrix.\n * @since 1.0.0\n */\nexport const determinant: (m: Matrix2x2) => number = internal.determinant\n\nexport const vectorProductLeft: {\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 2 * 6 = 17\n   * // y = 3 * 5 + 4 * 6 = 39\n   * const result = Matrix2x2.vectorProductLeft(matrix, vector)\n   * ```\n   * @param m - The matrix to multiply.\n   * @param v - The vector to multiply.\n   * @returns The resulting vector after the left vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 2 * 6 = 17\n   * // y = 3 * 5 + 4 * 6 = 39\n   * const result = Matrix2x2.vectorProductLeft(vector)(matrix)\n   * ```\n   * @param v - The vector to multiply.\n   * @returns A function that takes a matrix and returns the resulting vector after the left vector product.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.vectorProductLeft\n\nexport const vectorProductRight: {\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 3 * 6 = 33\n   * // y = 2 * 5 + 4 * 6 = 62\n   * const result = Matrix2x2.vectorProductRight(matrix, vector)\n   * ```\n   * @param m - The matrix to multiply.\n   * @param v - The vector to multiply.\n   * @returns The resulting vector after the right vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 3 * 6 = 33\n   * // y = 2 * 5 + 4 * 6 = 62\n   * const result = Matrix2x2.vectorProductRight(vector)(matrix)\n   * ```\n   * @param v - The vector to multiply.\n   * @returns A function that takes a matrix and returns the resulting vector after the right vector product.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.vectorProductRight\n\nexport const solveSystem: {\n  /**\n   * Solves a system of linear equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Solves the system of equations:\n   * // 1x + 2y = 5\n   * // 3x + 4y = 6\n   * const result = Matrix2x2.solveSystem(matrix, vector)\n   * ```\n   * @param m - The matrix representing the coefficients of the system.\n   * @param v - The vector representing the constants of the system.\n   * @returns The solution vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Solves a system of linear equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Solves the system of equations:\n   * // 1x + 2y = 5\n   * // 3x + 4y = 6\n   * const result = Matrix2x2.solveSystem(vector)(matrix)\n   * ```\n   * @param v - The vector representing the constants of the system.\n   * @returns A function that takes a matrix and returns the solution vector.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.solveSystem\n\n/**\n * Returns the rows of a matrix as an array of vectors.\n *\n * @param m - The matrix to convert.\n * @returns An array of vectors representing the rows of the matrix.\n * @since 1.0.0\n */\nexport const toRows: (m: Matrix2x2) => [Vector2, Vector2] = internal.toRows\n\n/**\n * Returns the columns of a matrix as an array of vectors.\n *\n * @param m - The matrix to convert.\n * @returns An array of vectors representing the columns of the matrix.\n * @since 1.0.0\n */\nexport const toColumns: (m: Matrix2x2) => [Vector2, Vector2] = internal.toColumns\n\n/**\n * Returns the specified row of a matrix as a vector.\n *\n * @param m - The matrix to get the row from.\n * @param row - The row index (0 or 1).\n * @returns A vector representing the specified row of the matrix.\n * @since 1.0.0\n */\nexport const rowVector: (m: Matrix2x2, row: Matrix2x2Coordinate) => Vector2 = internal.rowVector\n\n/**\n * Returns the specified column of a matrix as a vector.\n *\n * @param m - The matrix to get the column from.\n * @param column - The column index (0 or 1).\n * @returns A vector representing the specified column of the matrix.\n * @since 1.0.0\n */\nexport const columnVector: (m: Matrix2x2, column: Matrix2x2Coordinate) => Vector2 =\n  internal.columnVector\n\n/**\n * Transposes a matrix.\n *\n * @param m - The matrix to transpose.\n * @returns A new `Matrix2x2` instance representing the transposed matrix.\n * @since 1.0.0\n */\nexport const transpose: (m: Matrix2x2) => Matrix2x2 = internal.transpose\n\n/**\n * Reverses the rows of a matrix.\n *\n * @param m - The matrix to reverse the rows of.\n * @returns A new `Matrix2x2` instance with the rows reversed.\n * @since 1.0.0\n */\nexport const reverseRows: (m: Matrix2x2) => Matrix2x2 = internal.reverseRows\n\n/**\n * Reverses the columns of a matrix.\n *\n * @param m - The matrix to reverse the columns of.\n * @returns A new `Matrix2x2` instance with the columns reversed.\n * @since 1.0.0\n */\nexport const reverseColumns: (m: Matrix2x2) => Matrix2x2 = internal.reverseColumns\n","import type { ThreeDimensionalComponent, ThreeDimensionalIndex } from '../dimensions.ts'\nimport type { Pipeable } from '../utils.ts'\nimport type { Vector3 } from '../vector/vector3.ts'\nimport type { Matrix2x2 } from './matrix2x2.ts'\nimport type { Matrix3x3TypeId } from './matrix3x3.internal.ts'\nimport * as internal from './matrix3x3.internal.ts'\n\n/**\n * A coordinate in a 3x3 matrix.\n *\n * @since 1.0.0\n */\nexport type Matrix3x3Coordinate = ThreeDimensionalIndex | ThreeDimensionalComponent\n\n/**\n * A 3x3 matrix.\n *\n * All fields are readonly and immutable, and all operations create new instances.\n *\n * @since 1.0.0\n */\nexport interface Matrix3x3 extends Pipeable {\n  readonly [Matrix3x3TypeId]: Matrix3x3TypeId\n  /**\n   * The value of the first row and first column.\n   */\n  readonly m00: number\n  /**\n   * The value of the first row and second column.\n   */\n  readonly m01: number\n  /**\n   * The value of the first row and third column.\n   */\n  readonly m02: number\n  /**\n   * The value of the second row and first column.\n   */\n  readonly m10: number\n  /**\n   * The value of the second row and second column.\n   */\n  readonly m11: number\n  /**\n   * The value of the second row and third column.\n   */\n  readonly m12: number\n  /**\n   * The value of the third row and first column.\n   */\n  readonly m20: number\n  /**\n   * The value of the third row and second column.\n   */\n  readonly m21: number\n  /**\n   * The value of the third row and third column.\n   */\n  readonly m22: number\n}\n\n/**\n * Checks if a value is a `Matrix3x3`.\n *\n * @param m - The value to check.\n * @returns `true` if the value is a `Matrix3x3`, `false` otherwise.\n * @since 1.0.0\n */\nexport const isMatrix3x3: (m: unknown) => m is Matrix3x3 = internal.isMatrix3x3\n\nexport const equals: {\n  /**\n   * Checks if two `Matrix3x3` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param a - The first matrix.\n   * @param b - The second matrix.\n   * @returns `true` when each pair of components is within tolerance.\n   * @since 1.1.0\n   */\n  (a: Matrix3x3, b: Matrix3x3): boolean\n  /**\n   * Checks if two `Matrix3x3` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param b - The second matrix.\n   * @returns A function that takes the first matrix and returns the comparison result.\n   * @since 1.1.0\n   */\n  (b: Matrix3x3): (a: Matrix3x3) => boolean\n} = internal.equals\n\n/**\n * Creates a new `Matrix3x3` instance.\n *\n * @param m00 - The value of the first row and first column.\n * @param m01 - The value of the first row and second column.\n * @param m02 - The value of the first row and third column.\n * @param m10 - The value of the second row and first column.\n * @param m11 - The value of the second row and second column.\n * @param m12 - The value of the second row and third column.\n * @param m20 - The value of the third row and first column.\n * @param m21 - The value of the third row and second column.\n * @param m22 - The value of the third row and third column.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const make: (\n  m00?: number,\n  m01?: number,\n  m02?: number,\n  m10?: number,\n  m11?: number,\n  m12?: number,\n  m20?: number,\n  m21?: number,\n  m22?: number,\n) => Matrix3x3 = internal.make\n\n/**\n * Creates a new `Matrix3x3` instance from three row vectors.\n *\n * @param v0 - The first row vector.\n * @param v1 - The second row vector.\n * @param v2 - The third row vector.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const fromRows: (v0: Vector3, v1: Vector3, v2: Vector3) => Matrix3x3 = internal.fromRows\n\n/**\n * Creates a new `Matrix3x3` instance from three column vectors.\n *\n * @param v0 - The first column vector.\n * @param v1 - The second column vector.\n * @param v2 - The third column vector.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const fromColumns: (v0: Vector3, v1: Vector3, v2: Vector3) => Matrix3x3 =\n  internal.fromColumns\n\nexport const setRow: {\n  /**\n   * Sets a row of the matrix to a new vector.\n   *\n   * @param m - The matrix to set the row of.\n   * @param row - The row index to set.\n   * @param v - The new vector for the row.\n   * @returns A new matrix with the specified row set to the new vector.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, row: Matrix3x3Coordinate, v: Vector3): Matrix3x3\n  /**\n   * Sets a row of the matrix to a new vector.\n   *\n   * @param row - The row index to set.\n   * @param v - The new vector for the row.\n   * @returns A function that takes a matrix and returns a new matrix with the specified row set to the new vector.\n   * @since 1.0.0\n   */\n  (row: Matrix3x3Coordinate, v: Vector3): (m: Matrix3x3) => Matrix3x3\n} = internal.setRow\n\nexport const setColumn: {\n  /**\n   * Sets a column of the matrix to a new vector.\n   *\n   * @param m - The matrix to set the column of.\n   * @param column - The column index to set.\n   * @param v - The new vector for the column.\n   * @returns A new matrix with the specified column set to the new vector.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, column: Matrix3x3Coordinate, v: Vector3): Matrix3x3\n  /**\n   * Sets a column of the matrix to a new vector.\n   *\n   * @param column - The column index to set.\n   * @param v - The new vector for the column.\n   * @returns A function that takes a matrix and returns a new matrix with the specified column set to the new vector.\n   * @since 1.0.0\n   */\n  (column: Matrix3x3Coordinate, v: Vector3): (m: Matrix3x3) => Matrix3x3\n} = internal.setColumn\n\n/**\n * Calculates the determinant of a 3x3 matrix.\n *\n * @param m - The matrix to calculate the determinant of.\n * @returns The determinant of the matrix.\n * @since 1.0.0\n */\nexport const determinant: (m: Matrix3x3) => number = internal.determinant\n\n/**\n * Calculates the minor of a 3x3 matrix.\n *\n * @param m - The matrix to calculate the minor of.\n * @param row - The row index of the element to calculate the minor for.\n * @param column - The column index of the element to calculate the minor for.\n * @returns The minor of the matrix at the specified row and column.\n * @since 1.0.0\n */\nexport const minor: (\n  m: Matrix3x3,\n  row: Matrix3x3Coordinate,\n  column: Matrix3x3Coordinate,\n) => Matrix2x2 = internal.minor\n\nexport const vectorProductLeft: {\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 2 * 11 + 3 * 12\n   * // y = 4 * 10 + 5 * 11 + 6 * 12\n   * // z = 7 * 10 + 8 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductLeft(m, v)\n   * ```\n   * @param m - The matrix to calculate the left vector product with.\n   * @param v - The vector to calculate the left vector product with.\n   * @returns The resulting vector from the left vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 2 * 11 + 3 * 12\n   * // y = 4 * 10 + 5 * 11 + 6 * 12\n   * // z = 7 * 10 + 8 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductLeft(v)(m)\n   * ```\n   * @param v - The vector to calculate the left vector product with.\n   * @returns A function that takes a matrix and returns the resulting vector from the left vector product.\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.vectorProductLeft\n\nexport const vectorProductRight: {\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 4 * 11 + 7 * 12\n   * // y = 2 * 10 + 5 * 11 + 8 * 12\n   * // z = 3 * 10 + 6 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductRight(m, v)\n   * ```\n   * @param m - The matrix to calculate the right vector product with.\n   * @param v - The vector to calculate the right vector product with.\n   * @returns The resulting vector from the right vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 4 * 11 + 7 * 12\n   * // y = 2 * 10 + 5 * 11 + 8 * 12\n   * // z = 3 * 10 + 6 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductRight(v)(m)\n   * ```\n   * @param v - The vector to calculate the right vector product with.\n   * @returns A function that takes a matrix and returns the resulting vector from the right vector product.\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.vectorProductRight\n\nexport const solveSystem: {\n  /**\n   * Solves a system of equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = m00 * x + m01 * y + m02 * z\n   * // y = m10 * x + m11 * y + m12 * z\n   * // z = m20 * x + m21 * y + m22 * z\n   * const result = Matrix3x3.solveSystem(m, v)\n   * ```\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Solves a system of equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = m00 * x + m01 * y + m02 * z\n   * // y = m10 * x + m11 * y + m12 * z\n   * // z = m20 * x + m21 * y + m22 * z\n   * const result = Matrix3x3.solveSystem(v)(m)\n   * ```\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.solveSystem\n\n/**\n * Returns the rows of a `Matrix3x3` as an array of vectors.\n *\n * @param m - The matrix to get the rows from.\n * @returns An array of vectors representing the rows of the matrix.\n * @since 1.0.0\n */\nexport const toRows: (m: Matrix3x3) => [Vector3, Vector3, Vector3] = internal.toRows as (\n  m: Matrix3x3,\n) => [Vector3, Vector3, Vector3]\n\n/**\n * Returns the columns of a `Matrix3x3` as an array of vectors.\n *\n * @param m - The matrix to get the columns from.\n * @returns An array of vectors representing the columns of the matrix.\n * @since 1.0.0\n */\nexport const toColumns: (m: Matrix3x3) => [Vector3, Vector3, Vector3] = internal.toColumns as (\n  m: Matrix3x3,\n) => [Vector3, Vector3, Vector3]\n\nexport const rowVector: {\n  /**\n   * Gets a row vector from the matrix.\n   *\n   * @param m - The matrix to get the row vector from.\n   * @param row - The row index to get the vector from.\n   * @returns The row vector at the specified index.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, row: Matrix3x3Coordinate): Vector3\n  /**\n   * Gets a row vector from the matrix.\n   *\n   * @param row - The row index to get the vector from.\n   * @returns A function that takes a matrix and returns the row vector at the specified index.\n   * @since 1.0.0\n   */\n  (row: Matrix3x3Coordinate): (m: Matrix3x3) => Vector3\n} = internal.rowVector\n\nexport const columnVector: {\n  /**\n   * Gets a column vector from the matrix.\n   *\n   * @param m - The matrix to get the column vector from.\n   * @param column - The column index to get the vector from.\n   * @returns The column vector at the specified index.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, column: Matrix3x3Coordinate): Vector3\n  /**\n   * Gets a column vector from the matrix.\n   *\n   * @param column - The column index to get the vector from.\n   * @returns A function that takes a matrix and returns the column vector at the specified index.\n   * @since 1.0.0\n   */\n  (column: Matrix3x3Coordinate): (m: Matrix3x3) => Vector3\n} = internal.columnVector\n\n/**\n * Transposes a `Matrix3x3`.\n *\n * @param m - The matrix to transpose.\n * @returns A new `Matrix3x3` instance that is the transpose of the input matrix.\n * @since 1.0.0\n */\nexport const transpose: (m: Matrix3x3) => Matrix3x3 = internal.transpose\n\n/**\n * Reverses the rows of a `Matrix3x3`.\n *\n * @param m - The matrix to reverse the rows of.\n * @returns A new `Matrix3x3` instance with the rows reversed.\n * @since 1.0.0\n */\nexport const reverseRows: (m: Matrix3x3) => Matrix3x3 = internal.reverseRows\n\n/**\n * Reverses the columns of a `Matrix3x3`.\n *\n * @param m - The matrix to reverse the columns of.\n * @returns A new `Matrix3x3` instance with the columns reversed.\n * @since 1.0.0\n */\nexport const reverseColumns: (m: Matrix3x3) => Matrix3x3 = internal.reverseColumns\n"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAgDA,MAAa,cAA8CA;AAE3D,MAAaC,WAoBTC;;;;;;;;;AAUJ,MAAaC,SACXC;;;;;;;;;AAUF,MAAaC,aAAoDC;;;;;;;;;AAUjE,MAAaC,gBAAuDC;AAEpE,MAAaC,WAoBTC;AAEJ,MAAaC,cAoBTC;;;;;;;;AASJ,MAAaC,gBAAwCC;AAErD,MAAaC,sBA4CTC;AAEJ,MAAaC,uBA4CTC;AAEJ,MAAaC,gBA4CTC;;;;;;;;AASJ,MAAaC,WAA+CC;;;;;;;;AAS5D,MAAaC,cAAkDC;;;;;;;;;AAU/D,MAAaC,cAAiEC;;;;;;;;;AAU9E,MAAaC,iBACXC;;;;;;;;AASF,MAAaC,cAAyCC;;;;;;;;AAStD,MAAaC,gBAA2CC;;;;;;;;AASxD,MAAaC,mBAA8CC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AClS3D,MAAa,cAA8CC;AAE3D,MAAa,SAoBTC;;;;;;;;;;;;;;;;AAiBJ,MAAa,OAUIC;;;;;;;;;;AAWjB,MAAa,WAAiEC;;;;;;;;;;AAW9E,MAAa,cACXC;AAEF,MAAa,SAoBTC;AAEJ,MAAa,YAoBTC;;;;;;;;AASJ,MAAa,cAAwCC;;;;;;;;;;AAWrD,MAAa,QAIIC;AAEjB,MAAa,oBA8CTC;AAEJ,MAAa,qBA8CTC;AAEJ,MAAa,cAyCTC;;;;;;;;AASJ,MAAa,SAAwDC;;;;;;;;AAWrE,MAAa,YAA2DC;AAIxE,MAAa,YAkBTC;AAEJ,MAAa,eAkBTC;;;;;;;;AASJ,MAAa,YAAyCC;;;;;;;;AAStD,MAAa,cAA2CC;;;;;;;;AASxD,MAAa,iBAA8CC"}
+{"version":3,"file":"index.mjs","names":["internal.isMatrix2x2","equals","internal.equals","make","internal.make","fromRows","internal.fromRows","fromColumns","internal.fromColumns","setRow","internal.setRow","setColumn","internal.setColumn","determinant","internal.determinant","vectorProductLeft","internal.vectorProductLeft","vectorProductRight","internal.vectorProductRight","solveSystem","internal.solveSystem","toRows","internal.toRows","toColumns","internal.toColumns","rowVector","internal.rowVector","columnVector","internal.columnVector","transpose","internal.transpose","reverseRows","internal.reverseRows","reverseColumns","internal.reverseColumns","internal.isMatrix3x3","internal.equals","internal.make","internal.fromRows","internal.fromColumns","internal.setRow","internal.setColumn","internal.determinant","internal.minor","internal.vectorProductLeft","internal.vectorProductRight","internal.solveSystem","internal.toRows","internal.toColumns","internal.rowVector","internal.columnVector","internal.transpose","internal.reverseRows","internal.reverseColumns","internal.identity","internal.multiply","internal.inverse","internal.inverseUnsafe","internal.isInvertible","internal.asInvertible"],"sources":["../../src/matrix/matrix2x2.ts","../../src/matrix/matrix3x3.ts"],"sourcesContent":["import type { TwoDimensionalComponent, TwoDimensionalIndex } from '../dimensions.ts'\nimport type { Pipeable } from '../utils.ts'\nimport type { Vector2 } from '../vector/vector2.ts'\nimport type { Matrix2x2TypeId } from './matrix2x2.internal.ts'\nimport * as internal from './matrix2x2.internal.ts'\n\n/**\n * A type representing the coordinate of a 2D matrix.\n *\n * @since 1.0.0\n */\nexport type Matrix2x2Coordinate = TwoDimensionalIndex | TwoDimensionalComponent\n\n/**\n * A 2D matrix with 2 rows and 2 columns.\n *\n * All fields are readonly and immutable, and all operations create new instances.\n *\n * @since 1.0.0\n */\nexport interface Matrix2x2 extends Pipeable {\n  readonly [Matrix2x2TypeId]: Matrix2x2TypeId\n\n  /**\n   * The value in the first row and first column.\n   */\n  readonly m00: number\n  /**\n   * The value in the first row and second column.\n   */\n  readonly m01: number\n  /**\n   * The value in the second row and first column.\n   */\n  readonly m10: number\n  /**\n   * The value in the second row and second column.\n   */\n  readonly m11: number\n}\n\n/**\n * Checks if a value is a `Matrix2x2`.\n *\n * @param m - The value to check.\n * @returns `true` if the value is a `Matrix2x2`, `false` otherwise.\n * @since 1.0.0\n */\nexport const isMatrix2x2: (m: unknown) => m is Matrix2x2 = internal.isMatrix2x2\n\nexport const equals: {\n  /**\n   * Checks if two `Matrix2x2` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param a - The first matrix.\n   * @param b - The second matrix.\n   * @returns `true` when each pair of components is within tolerance.\n   * @since 1.1.0\n   */\n  (a: Matrix2x2, b: Matrix2x2): boolean\n  /**\n   * Checks if two `Matrix2x2` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param b - The second matrix.\n   * @returns A function that takes the first matrix and returns the comparison result.\n   * @since 1.1.0\n   */\n  (b: Matrix2x2): (a: Matrix2x2) => boolean\n} = internal.equals\n\n/**\n * Creates a new `Matrix2x2` instance.\n *\n * @param m00 - The value in the first row and first column. Defaults to 0.\n * @param m01 - The value in the first row and second column. Defaults to 0.\n * @param m10 - The value in the second row and first column. Defaults to 0.\n * @param m11 - The value in the second row and second column. Defaults to 0.\n */\nexport const make: (m00?: number, m01?: number, m10?: number, m11?: number) => Matrix2x2 =\n  internal.make\n\n/**\n * Creates a new `Matrix2x2` instance from two vectors representing rows.\n *\n * @param v0 - The first row vector.\n * @param v1 - The second row vector.\n * @returns A new `Matrix2x2` instance.\n * @since 1.0.0\n */\nexport const fromRows: (v0: Vector2, v1: Vector2) => Matrix2x2 = internal.fromRows\n\n/**\n * Creates a new `Matrix2x2` instance from two vectors representing columns.\n *\n * @param v0 - The first column vector.\n * @param v1 - The second column vector.\n * @returns A new `Matrix2x2` instance.\n * @since 1.0.0\n */\nexport const fromColumns: (v0: Vector2, v1: Vector2) => Matrix2x2 = internal.fromColumns\n\nexport const setRow: {\n  /**\n   * Sets the specified row of the matrix to the given vector.\n   *\n   * @param m - The matrix to modify.\n   * @param row - The row index (0 or 1).\n   * @param v - The vector to set the row to.\n   * @returns A new `Matrix2x2` instance with the specified row set to the given vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, row: Matrix2x2Coordinate, v: Vector2): Matrix2x2\n  /**\n   * Sets the specified row of the matrix to the given vector.\n   *\n   * @param row - The row index (0 or 1).\n   * @param v - The vector to set the row to.\n   * @returns A function that takes a matrix and returns a new `Matrix2x2` instance with the specified row set to the given vector.\n   * @since 1.0.0\n   */\n  (row: Matrix2x2Coordinate, v: Vector2): (m: Matrix2x2) => Matrix2x2\n} = internal.setRow\n\nexport const setColumn: {\n  /**\n   * Sets the specified column of the matrix to the given vector.\n   *\n   * @param m - The matrix to modify.\n   * @param column - The column index (0 or 1).\n   * @param v - The vector to set the column to.\n   * @returns A new `Matrix2x2` instance with the specified column set to the given vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, column: Matrix2x2Coordinate, v: Vector2): Matrix2x2\n  /**\n   * Sets the specified column of the matrix to the given vector.\n   *\n   * @param column - The column index (0 or 1).\n   * @param v - The vector to set the column to.\n   * @returns A function that takes a matrix and returns a new `Matrix2x2` instance with the specified column set to the given vector.\n   * @since 1.0.0\n   */\n  (column: Matrix2x2Coordinate, v: Vector2): (m: Matrix2x2) => Matrix2x2\n} = internal.setColumn\n\n/**\n * Calculates the determinant of a 2D matrix.\n *\n * @param m - The matrix to calculate the determinant for.\n * @returns The determinant of the matrix.\n * @since 1.0.0\n */\nexport const determinant: (m: Matrix2x2) => number = internal.determinant\n\nexport const vectorProductLeft: {\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 2 * 6 = 17\n   * // y = 3 * 5 + 4 * 6 = 39\n   * const result = Matrix2x2.vectorProductLeft(matrix, vector)\n   * ```\n   * @param m - The matrix to multiply.\n   * @param v - The vector to multiply.\n   * @returns The resulting vector after the left vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 2 * 6 = 17\n   * // y = 3 * 5 + 4 * 6 = 39\n   * const result = Matrix2x2.vectorProductLeft(vector)(matrix)\n   * ```\n   * @param v - The vector to multiply.\n   * @returns A function that takes a matrix and returns the resulting vector after the left vector product.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.vectorProductLeft\n\nexport const vectorProductRight: {\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 3 * 6 = 33\n   * // y = 2 * 5 + 4 * 6 = 62\n   * const result = Matrix2x2.vectorProductRight(matrix, vector)\n   * ```\n   * @param m - The matrix to multiply.\n   * @param v - The vector to multiply.\n   * @returns The resulting vector after the right vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Performs:\n   * // x = 1 * 5 + 3 * 6 = 33\n   * // y = 2 * 5 + 4 * 6 = 62\n   * const result = Matrix2x2.vectorProductRight(vector)(matrix)\n   * ```\n   * @param v - The vector to multiply.\n   * @returns A function that takes a matrix and returns the resulting vector after the right vector product.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.vectorProductRight\n\nexport const solveSystem: {\n  /**\n   * Solves a system of linear equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Solves the system of equations:\n   * // 1x + 2y = 5\n   * // 3x + 4y = 6\n   * const result = Matrix2x2.solveSystem(matrix, vector)\n   * ```\n   * @param m - The matrix representing the coefficients of the system.\n   * @param v - The vector representing the constants of the system.\n   * @returns The solution vector.\n   * @since 1.0.0\n   */\n  (m: Matrix2x2, v: Vector2): Vector2\n  /**\n   * Solves a system of linear equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix2x2 from 'curvy/matrix2x2'\n   * import * as Vector2 from 'curvy/vector2'\n   *\n   * const matrix = Matrix2x2.make(1, 2, 3, 4)\n   * const vector = Vector2.make(5, 6)\n   *\n   * // Solves the system of equations:\n   * // 1x + 2y = 5\n   * // 3x + 4y = 6\n   * const result = Matrix2x2.solveSystem(vector)(matrix)\n   * ```\n   * @param v - The vector representing the constants of the system.\n   * @returns A function that takes a matrix and returns the solution vector.\n   * @since 1.0.0\n   */\n  (v: Vector2): (m: Matrix2x2) => Vector2\n} = internal.solveSystem\n\n/**\n * Returns the rows of a matrix as an array of vectors.\n *\n * @param m - The matrix to convert.\n * @returns An array of vectors representing the rows of the matrix.\n * @since 1.0.0\n */\nexport const toRows: (m: Matrix2x2) => [Vector2, Vector2] = internal.toRows\n\n/**\n * Returns the columns of a matrix as an array of vectors.\n *\n * @param m - The matrix to convert.\n * @returns An array of vectors representing the columns of the matrix.\n * @since 1.0.0\n */\nexport const toColumns: (m: Matrix2x2) => [Vector2, Vector2] = internal.toColumns\n\n/**\n * Returns the specified row of a matrix as a vector.\n *\n * @param m - The matrix to get the row from.\n * @param row - The row index (0 or 1).\n * @returns A vector representing the specified row of the matrix.\n * @since 1.0.0\n */\nexport const rowVector: (m: Matrix2x2, row: Matrix2x2Coordinate) => Vector2 = internal.rowVector\n\n/**\n * Returns the specified column of a matrix as a vector.\n *\n * @param m - The matrix to get the column from.\n * @param column - The column index (0 or 1).\n * @returns A vector representing the specified column of the matrix.\n * @since 1.0.0\n */\nexport const columnVector: (m: Matrix2x2, column: Matrix2x2Coordinate) => Vector2 =\n  internal.columnVector\n\n/**\n * Transposes a matrix.\n *\n * @param m - The matrix to transpose.\n * @returns A new `Matrix2x2` instance representing the transposed matrix.\n * @since 1.0.0\n */\nexport const transpose: (m: Matrix2x2) => Matrix2x2 = internal.transpose\n\n/**\n * Reverses the rows of a matrix.\n *\n * @param m - The matrix to reverse the rows of.\n * @returns A new `Matrix2x2` instance with the rows reversed.\n * @since 1.0.0\n */\nexport const reverseRows: (m: Matrix2x2) => Matrix2x2 = internal.reverseRows\n\n/**\n * Reverses the columns of a matrix.\n *\n * @param m - The matrix to reverse the columns of.\n * @returns A new `Matrix2x2` instance with the columns reversed.\n * @since 1.0.0\n */\nexport const reverseColumns: (m: Matrix2x2) => Matrix2x2 = internal.reverseColumns\n","import type { ThreeDimensionalComponent, ThreeDimensionalIndex } from '../dimensions.ts'\nimport type * as Solution from '../solution/solution.ts'\nimport type { Pipeable } from '../utils.ts'\nimport type { Vector3 } from '../vector/vector3.ts'\nimport type { Matrix2x2 } from './matrix2x2.ts'\nimport type { Matrix3x3TypeId } from './matrix3x3.internal.ts'\nimport * as internal from './matrix3x3.internal.ts'\nimport type { Invertible, MatrixTraits } from './traits.ts'\n\nexport type { Invertible } from './traits.ts'\n\n/**\n * A coordinate in a 3x3 matrix.\n *\n * @since 1.0.0\n */\nexport type Matrix3x3Coordinate = ThreeDimensionalIndex | ThreeDimensionalComponent\n\n/**\n * A 3x3 matrix.\n *\n * All fields are readonly and immutable, and all operations create new instances.\n *\n * The `Traits` type parameter is a phantom marker that accumulates trait\n * brands as the matrix is refined via `isInvertible` / `asInvertible`. The\n * runtime value carries no trait data — `Traits` only flows through the type\n * system to enable tighter return types on operations like {@link inverse}.\n *\n * @since 1.0.0\n */\nexport interface Matrix3x3<out Traits = unknown> extends Pipeable {\n  readonly [Matrix3x3TypeId]: Matrix3x3TypeId\n  readonly [MatrixTraits]: Traits\n  /**\n   * The value of the first row and first column.\n   */\n  readonly m00: number\n  /**\n   * The value of the first row and second column.\n   */\n  readonly m01: number\n  /**\n   * The value of the first row and third column.\n   */\n  readonly m02: number\n  /**\n   * The value of the second row and first column.\n   */\n  readonly m10: number\n  /**\n   * The value of the second row and second column.\n   */\n  readonly m11: number\n  /**\n   * The value of the second row and third column.\n   */\n  readonly m12: number\n  /**\n   * The value of the third row and first column.\n   */\n  readonly m20: number\n  /**\n   * The value of the third row and second column.\n   */\n  readonly m21: number\n  /**\n   * The value of the third row and third column.\n   */\n  readonly m22: number\n}\n\n/**\n * Checks if a value is a `Matrix3x3`.\n *\n * @param m - The value to check.\n * @returns `true` if the value is a `Matrix3x3`, `false` otherwise.\n * @since 1.0.0\n */\nexport const isMatrix3x3: (m: unknown) => m is Matrix3x3 = internal.isMatrix3x3\n\nexport const equals: {\n  /**\n   * Checks if two `Matrix3x3` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param a - The first matrix.\n   * @param b - The second matrix.\n   * @returns `true` when each pair of components is within tolerance.\n   * @since 1.1.0\n   */\n  (a: Matrix3x3, b: Matrix3x3): boolean\n  /**\n   * Checks if two `Matrix3x3` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param b - The second matrix.\n   * @returns A function that takes the first matrix and returns the comparison result.\n   * @since 1.1.0\n   */\n  (b: Matrix3x3): (a: Matrix3x3) => boolean\n} = internal.equals\n\n/**\n * Creates a new `Matrix3x3` instance.\n *\n * @param m00 - The value of the first row and first column.\n * @param m01 - The value of the first row and second column.\n * @param m02 - The value of the first row and third column.\n * @param m10 - The value of the second row and first column.\n * @param m11 - The value of the second row and second column.\n * @param m12 - The value of the second row and third column.\n * @param m20 - The value of the third row and first column.\n * @param m21 - The value of the third row and second column.\n * @param m22 - The value of the third row and third column.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const make: (\n  m00?: number,\n  m01?: number,\n  m02?: number,\n  m10?: number,\n  m11?: number,\n  m12?: number,\n  m20?: number,\n  m21?: number,\n  m22?: number,\n) => Matrix3x3 = internal.make\n\n/**\n * Creates a new `Matrix3x3` instance from three row vectors.\n *\n * @param v0 - The first row vector.\n * @param v1 - The second row vector.\n * @param v2 - The third row vector.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const fromRows: (v0: Vector3, v1: Vector3, v2: Vector3) => Matrix3x3 = internal.fromRows\n\n/**\n * Creates a new `Matrix3x3` instance from three column vectors.\n *\n * @param v0 - The first column vector.\n * @param v1 - The second column vector.\n * @param v2 - The third column vector.\n * @returns A new `Matrix3x3` instance.\n * @since 1.0.0\n */\nexport const fromColumns: (v0: Vector3, v1: Vector3, v2: Vector3) => Matrix3x3 =\n  internal.fromColumns\n\nexport const setRow: {\n  /**\n   * Sets a row of the matrix to a new vector.\n   *\n   * @param m - The matrix to set the row of.\n   * @param row - The row index to set.\n   * @param v - The new vector for the row.\n   * @returns A new matrix with the specified row set to the new vector.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, row: Matrix3x3Coordinate, v: Vector3): Matrix3x3\n  /**\n   * Sets a row of the matrix to a new vector.\n   *\n   * @param row - The row index to set.\n   * @param v - The new vector for the row.\n   * @returns A function that takes a matrix and returns a new matrix with the specified row set to the new vector.\n   * @since 1.0.0\n   */\n  (row: Matrix3x3Coordinate, v: Vector3): (m: Matrix3x3) => Matrix3x3\n} = internal.setRow\n\nexport const setColumn: {\n  /**\n   * Sets a column of the matrix to a new vector.\n   *\n   * @param m - The matrix to set the column of.\n   * @param column - The column index to set.\n   * @param v - The new vector for the column.\n   * @returns A new matrix with the specified column set to the new vector.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, column: Matrix3x3Coordinate, v: Vector3): Matrix3x3\n  /**\n   * Sets a column of the matrix to a new vector.\n   *\n   * @param column - The column index to set.\n   * @param v - The new vector for the column.\n   * @returns A function that takes a matrix and returns a new matrix with the specified column set to the new vector.\n   * @since 1.0.0\n   */\n  (column: Matrix3x3Coordinate, v: Vector3): (m: Matrix3x3) => Matrix3x3\n} = internal.setColumn\n\n/**\n * Calculates the determinant of a 3x3 matrix.\n *\n * @param m - The matrix to calculate the determinant of.\n * @returns The determinant of the matrix.\n * @since 1.0.0\n */\nexport const determinant: (m: Matrix3x3) => number = internal.determinant\n\n/**\n * Calculates the minor of a 3x3 matrix.\n *\n * @param m - The matrix to calculate the minor of.\n * @param row - The row index of the element to calculate the minor for.\n * @param column - The column index of the element to calculate the minor for.\n * @returns The minor of the matrix at the specified row and column.\n * @since 1.0.0\n */\nexport const minor: (\n  m: Matrix3x3,\n  row: Matrix3x3Coordinate,\n  column: Matrix3x3Coordinate,\n) => Matrix2x2 = internal.minor\n\nexport const vectorProductLeft: {\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 2 * 11 + 3 * 12\n   * // y = 4 * 10 + 5 * 11 + 6 * 12\n   * // z = 7 * 10 + 8 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductLeft(m, v)\n   * ```\n   * @param m - The matrix to calculate the left vector product with.\n   * @param v - The vector to calculate the left vector product with.\n   * @returns The resulting vector from the left vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Calculates the left vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 2 * 11 + 3 * 12\n   * // y = 4 * 10 + 5 * 11 + 6 * 12\n   * // z = 7 * 10 + 8 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductLeft(v)(m)\n   * ```\n   * @param v - The vector to calculate the left vector product with.\n   * @returns A function that takes a matrix and returns the resulting vector from the left vector product.\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.vectorProductLeft\n\nexport const vectorProductRight: {\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 4 * 11 + 7 * 12\n   * // y = 2 * 10 + 5 * 11 + 8 * 12\n   * // z = 3 * 10 + 6 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductRight(m, v)\n   * ```\n   * @param m - The matrix to calculate the right vector product with.\n   * @param v - The vector to calculate the right vector product with.\n   * @returns The resulting vector from the right vector product.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Calculates the right vector product of a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = 1 * 10 + 4 * 11 + 7 * 12\n   * // y = 2 * 10 + 5 * 11 + 8 * 12\n   * // z = 3 * 10 + 6 * 11 + 9 * 12\n   * const result = Matrix3x3.vectorProductRight(v)(m)\n   * ```\n   * @param v - The vector to calculate the right vector product with.\n   * @returns A function that takes a matrix and returns the resulting vector from the right vector product.\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.vectorProductRight\n\nexport const solveSystem: {\n  /**\n   * Solves a system of equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = m00 * x + m01 * y + m02 * z\n   * // y = m10 * x + m11 * y + m12 * z\n   * // z = m20 * x + m21 * y + m22 * z\n   * const result = Matrix3x3.solveSystem(m, v)\n   * ```\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, v: Vector3): Vector3\n  /**\n   * Solves a system of equations represented by a matrix and a vector.\n   *\n   * @example\n   * ```ts\n   * import * as Matrix3x3 from 'curvy/matrix3x3'\n   * import * as Vector3 from 'curvy/vector3'\n   *\n   * const m = Matrix3x3.make(1, 2, 3, 4, 5, 6, 7, 8, 9)\n   * const v = Vector3.vector3(10, 11, 12)\n   *\n   * // Performs:\n   * // x = m00 * x + m01 * y + m02 * z\n   * // y = m10 * x + m11 * y + m12 * z\n   * // z = m20 * x + m21 * y + m22 * z\n   * const result = Matrix3x3.solveSystem(v)(m)\n   * ```\n   * @since 1.0.0\n   */\n  (v: Vector3): (m: Matrix3x3) => Vector3\n} = internal.solveSystem\n\n/**\n * Returns the rows of a `Matrix3x3` as an array of vectors.\n *\n * @param m - The matrix to get the rows from.\n * @returns An array of vectors representing the rows of the matrix.\n * @since 1.0.0\n */\nexport const toRows: (m: Matrix3x3) => [Vector3, Vector3, Vector3] = internal.toRows as (\n  m: Matrix3x3,\n) => [Vector3, Vector3, Vector3]\n\n/**\n * Returns the columns of a `Matrix3x3` as an array of vectors.\n *\n * @param m - The matrix to get the columns from.\n * @returns An array of vectors representing the columns of the matrix.\n * @since 1.0.0\n */\nexport const toColumns: (m: Matrix3x3) => [Vector3, Vector3, Vector3] = internal.toColumns as (\n  m: Matrix3x3,\n) => [Vector3, Vector3, Vector3]\n\nexport const rowVector: {\n  /**\n   * Gets a row vector from the matrix.\n   *\n   * @param m - The matrix to get the row vector from.\n   * @param row - The row index to get the vector from.\n   * @returns The row vector at the specified index.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, row: Matrix3x3Coordinate): Vector3\n  /**\n   * Gets a row vector from the matrix.\n   *\n   * @param row - The row index to get the vector from.\n   * @returns A function that takes a matrix and returns the row vector at the specified index.\n   * @since 1.0.0\n   */\n  (row: Matrix3x3Coordinate): (m: Matrix3x3) => Vector3\n} = internal.rowVector\n\nexport const columnVector: {\n  /**\n   * Gets a column vector from the matrix.\n   *\n   * @param m - The matrix to get the column vector from.\n   * @param column - The column index to get the vector from.\n   * @returns The column vector at the specified index.\n   * @since 1.0.0\n   */\n  (m: Matrix3x3, column: Matrix3x3Coordinate): Vector3\n  /**\n   * Gets a column vector from the matrix.\n   *\n   * @param column - The column index to get the vector from.\n   * @returns A function that takes a matrix and returns the column vector at the specified index.\n   * @since 1.0.0\n   */\n  (column: Matrix3x3Coordinate): (m: Matrix3x3) => Vector3\n} = internal.columnVector\n\n/**\n * Transposes a `Matrix3x3`.\n *\n * @param m - The matrix to transpose.\n * @returns A new `Matrix3x3` instance that is the transpose of the input matrix.\n * @since 1.0.0\n */\nexport const transpose: (m: Matrix3x3) => Matrix3x3 = internal.transpose\n\n/**\n * Reverses the rows of a `Matrix3x3`.\n *\n * @param m - The matrix to reverse the rows of.\n * @returns A new `Matrix3x3` instance with the rows reversed.\n * @since 1.0.0\n */\nexport const reverseRows: (m: Matrix3x3) => Matrix3x3 = internal.reverseRows\n\n/**\n * Reverses the columns of a `Matrix3x3`.\n *\n * @param m - The matrix to reverse the columns of.\n * @returns A new `Matrix3x3` instance with the columns reversed.\n * @since 1.0.0\n */\nexport const reverseColumns: (m: Matrix3x3) => Matrix3x3 = internal.reverseColumns\n\n/**\n * The 3x3 identity matrix.\n *\n * @since 2.0.0\n */\nexport const identity: Matrix3x3 = internal.identity\n\nexport const multiply: {\n  /**\n   * Multiplies two `Matrix3x3` instances. Result is `a · b`.\n   *\n   * @param a - The left-hand matrix.\n   * @param b - The right-hand matrix.\n   * @returns The product matrix.\n   * @since 2.0.0\n   */\n  (a: Matrix3x3, b: Matrix3x3): Matrix3x3\n  /**\n   * Multiplies two `Matrix3x3` instances. Result is `a · b`.\n   *\n   * @param b - The right-hand matrix.\n   * @returns A function that takes the left-hand matrix and returns the product.\n   * @since 2.0.0\n   */\n  (b: Matrix3x3): (a: Matrix3x3) => Matrix3x3\n} = internal.multiply\n\n// The implementation always runs the same det/threshold branch. For\n// `Invertible`-branded matrices the empty branch is provably unreachable,\n// but TS can't see that — the cast pairs the type-level promise with the\n// runtime brand check.\nexport const inverse: {\n  /**\n   * Returns the inverse of an `Invertible`-branded `Matrix3x3`. Because the\n   * brand guarantees non-singularity, the result is always a single matrix —\n   * wrapped in `Solution.One`.\n   *\n   * @param m - The matrix to invert.\n   * @since 2.0.0\n   */\n  <T extends Invertible>(m: Matrix3x3<T>): Solution.One<Matrix3x3<T>>\n  /**\n   * Returns the inverse of a `Matrix3x3`. Returns `Solution.none` if the\n   * matrix is singular (determinant approximately zero). Use {@link asInvertible}\n   * or {@link isInvertible} to refine the type and tighten the return to\n   * `Solution.One`, or use {@link inverseUnsafe} if you want the throwing form.\n   *\n   * @param m - The matrix to invert.\n   * @since 2.0.0\n   */\n  <T>(m: Matrix3x3<T>): Solution.AtMostOne<Matrix3x3<T>>\n} = internal.inverse as never\n\n/**\n * Returns the inverse of a `Matrix3x3`. Throws when the matrix is singular\n * (determinant approximately zero). Useful when you've already validated\n * invertibility out-of-band and want the unwrapped matrix directly; otherwise\n * prefer {@link inverse}.\n *\n * @param m - The matrix to invert.\n * @returns The inverse matrix.\n * @throws If `m` is singular.\n * @since 2.0.0\n */\nexport const inverseUnsafe: (m: Matrix3x3) => Matrix3x3 = internal.inverseUnsafe\n\n/**\n * Type-narrowing predicate: refines a `Matrix3x3<T>` to\n * `Matrix3x3<T & Invertible>` when the matrix's determinant is non-zero\n * (within the default absolute tolerance).\n *\n * @param m - The matrix to check.\n * @returns `true` when `m` is invertible.\n * @since 2.0.0\n */\nexport const isInvertible: <T>(m: Matrix3x3<T>) => m is Matrix3x3<T & Invertible> =\n  internal.isInvertible\n\n/**\n * Asserts that the matrix is invertible, throwing on a singular matrix.\n *\n * @param m - The matrix to assert against.\n * @returns The same matrix, typed with the `Invertible` brand.\n * @throws When `m` is singular.\n * @since 2.0.0\n */\nexport const asInvertible: <T>(m: Matrix3x3<T>) => Matrix3x3<T & Invertible> = internal.asInvertible\n"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;AAgDA,MAAa,cAA8CA;AAE3D,MAAaC,WAoBTC;;;;;;;;;AAUJ,MAAaC,SACXC;;;;;;;;;AAUF,MAAaC,aAAoDC;;;;;;;;;AAUjE,MAAaC,gBAAuDC;AAEpE,MAAaC,WAoBTC;AAEJ,MAAaC,cAoBTC;;;;;;;;AASJ,MAAaC,gBAAwCC;AAErD,MAAaC,sBA4CTC;AAEJ,MAAaC,uBA4CTC;AAEJ,MAAaC,gBA4CTC;;;;;;;;AASJ,MAAaC,WAA+CC;;;;;;;;AAS5D,MAAaC,cAAkDC;;;;;;;;;AAU/D,MAAaC,cAAiEC;;;;;;;;;AAU9E,MAAaC,iBACXC;;;;;;;;AASF,MAAaC,cAAyCC;;;;;;;;AAStD,MAAaC,gBAA2CC;;;;;;;;AASxD,MAAaC,mBAA8CC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;ACxR3D,MAAa,cAA8CC;AAE3D,MAAa,SAoBTC;;;;;;;;;;;;;;;;AAiBJ,MAAa,OAUIC;;;;;;;;;;AAWjB,MAAa,WAAiEC;;;;;;;;;;AAW9E,MAAa,cACXC;AAEF,MAAa,SAoBTC;AAEJ,MAAa,YAoBTC;;;;;;;;AASJ,MAAa,cAAwCC;;;;;;;;;;AAWrD,MAAa,QAIIC;AAEjB,MAAa,oBA8CTC;AAEJ,MAAa,qBA8CTC;AAEJ,MAAa,cAyCTC;;;;;;;;AASJ,MAAa,SAAwDC;;;;;;;;AAWrE,MAAa,YAA2DC;AAIxE,MAAa,YAkBTC;AAEJ,MAAa,eAkBTC;;;;;;;;AASJ,MAAa,YAAyCC;;;;;;;;AAStD,MAAa,cAA2CC;;;;;;;;AASxD,MAAa,iBAA8CC;;;;;;AAO3D,MAAa,WAAsBC;AAEnC,MAAa,WAkBTC;AAMJ,MAAa,UAoBTC;;;;;;;;;;;;AAaJ,MAAa,gBAA6CC;;;;;;;;;;AAW1D,MAAa,eACXC;;;;;;;;;AAUF,MAAa,eAAkEC"}

package/dist/matrix/matrix3x3.d.ts

@@ -1,8 +1,11 @@
 import type { ThreeDimensionalComponent, ThreeDimensionalIndex } from '../dimensions.ts';
+import type * as Solution from '../solution/solution.ts';
 import type { Pipeable } from '../utils.ts';
 import type { Vector3 } from '../vector/vector3.ts';
 import type { Matrix2x2 } from './matrix2x2.ts';
 import type { Matrix3x3TypeId } from './matrix3x3.internal.ts';
+import type { Invertible, MatrixTraits } from './traits.ts';
+export type { Invertible } from './traits.ts';
 /**
  * A coordinate in a 3x3 matrix.
  *
@@ -14,10 +17,16 @@ export type Matrix3x3Coordinate = ThreeDimensionalIndex | ThreeDimensionalCompon
  *
  * All fields are readonly and immutable, and all operations create new instances.
  *
+ * The `Traits` type parameter is a phantom marker that accumulates trait
+ * brands as the matrix is refined via `isInvertible` / `asInvertible`. The
+ * runtime value carries no trait data — `Traits` only flows through the type
+ * system to enable tighter return types on operations like {@link inverse}.
+ *
  * @since 1.0.0
  */
-export interface Matrix3x3 extends Pipeable {
+export interface Matrix3x3<out Traits = unknown> extends Pipeable {
     readonly [Matrix3x3TypeId]: Matrix3x3TypeId;
+    readonly [MatrixTraits]: Traits;
     /**
      * The value of the first row and first column.
      */
@@ -394,3 +403,80 @@ export declare const reverseRows: (m: Matrix3x3) => Matrix3x3;
  * @since 1.0.0
  */
 export declare const reverseColumns: (m: Matrix3x3) => Matrix3x3;
+/**
+ * The 3x3 identity matrix.
+ *
+ * @since 2.0.0
+ */
+export declare const identity: Matrix3x3;
+export declare const multiply: {
+    /**
+     * Multiplies two `Matrix3x3` instances. Result is `a · b`.
+     *
+     * @param a - The left-hand matrix.
+     * @param b - The right-hand matrix.
+     * @returns The product matrix.
+     * @since 2.0.0
+     */
+    (a: Matrix3x3, b: Matrix3x3): Matrix3x3;
+    /**
+     * Multiplies two `Matrix3x3` instances. Result is `a · b`.
+     *
+     * @param b - The right-hand matrix.
+     * @returns A function that takes the left-hand matrix and returns the product.
+     * @since 2.0.0
+     */
+    (b: Matrix3x3): (a: Matrix3x3) => Matrix3x3;
+};
+export declare const inverse: {
+    /**
+     * Returns the inverse of an `Invertible`-branded `Matrix3x3`. Because the
+     * brand guarantees non-singularity, the result is always a single matrix —
+     * wrapped in `Solution.One`.
+     *
+     * @param m - The matrix to invert.
+     * @since 2.0.0
+     */
+    <T extends Invertible>(m: Matrix3x3<T>): Solution.One<Matrix3x3<T>>;
+    /**
+     * Returns the inverse of a `Matrix3x3`. Returns `Solution.none` if the
+     * matrix is singular (determinant approximately zero). Use {@link asInvertible}
+     * or {@link isInvertible} to refine the type and tighten the return to
+     * `Solution.One`, or use {@link inverseUnsafe} if you want the throwing form.
+     *
+     * @param m - The matrix to invert.
+     * @since 2.0.0
+     */
+    <T>(m: Matrix3x3<T>): Solution.AtMostOne<Matrix3x3<T>>;
+};
+/**
+ * Returns the inverse of a `Matrix3x3`. Throws when the matrix is singular
+ * (determinant approximately zero). Useful when you've already validated
+ * invertibility out-of-band and want the unwrapped matrix directly; otherwise
+ * prefer {@link inverse}.
+ *
+ * @param m - The matrix to invert.
+ * @returns The inverse matrix.
+ * @throws If `m` is singular.
+ * @since 2.0.0
+ */
+export declare const inverseUnsafe: (m: Matrix3x3) => Matrix3x3;
+/**
+ * Type-narrowing predicate: refines a `Matrix3x3<T>` to
+ * `Matrix3x3<T & Invertible>` when the matrix's determinant is non-zero
+ * (within the default absolute tolerance).
+ *
+ * @param m - The matrix to check.
+ * @returns `true` when `m` is invertible.
+ * @since 2.0.0
+ */
+export declare const isInvertible: <T>(m: Matrix3x3<T>) => m is Matrix3x3<T & Invertible>;
+/**
+ * Asserts that the matrix is invertible, throwing on a singular matrix.
+ *
+ * @param m - The matrix to assert against.
+ * @returns The same matrix, typed with the `Invertible` brand.
+ * @throws When `m` is singular.
+ * @since 2.0.0
+ */
+export declare const asInvertible: <T>(m: Matrix3x3<T>) => Matrix3x3<T & Invertible>;

package/dist/matrix/matrix4x4.d.ts

@@ -1,8 +1,11 @@
 import type { FourDimensionalComponent, FourDimensionalIndex } from '../dimensions.ts';
+import type * as Solution from '../solution/solution.ts';
 import type { Pipeable } from '../utils.ts';
 import type { Vector4 } from '../vector/vector4.ts';
 import type { Matrix3x3 } from './matrix3x3.ts';
 import type { Matrix4x4TypeId } from './matrix4x4.internal.ts';
+import type { Invertible, MatrixTraits } from './traits.ts';
+export type { Invertible } from './traits.ts';
 /**
  * A coordinate in a 4x4 matrix.
  *
@@ -14,10 +17,16 @@ export type Matrix4x4Coordinate = FourDimensionalIndex | FourDimensionalComponen
  *
  * All fields are readonly and immutable, and all operations create new instances.
  *
+ * The `Traits` type parameter is a phantom marker that accumulates trait
+ * brands as the matrix is refined via `isInvertible` / `asInvertible`. The
+ * runtime value carries no trait data — `Traits` only flows through the type
+ * system to enable tighter return types on operations like {@link inverse}.
+ *
  * @since 1.0.0
  */
-export interface Matrix4x4 extends Pipeable {
+export interface Matrix4x4<out Traits = unknown> extends Pipeable {
     readonly [Matrix4x4TypeId]: Matrix4x4TypeId;
+    readonly [MatrixTraits]: Traits;
     /**
      * The value of the first row and first column.
      */
@@ -330,16 +339,58 @@ export declare const multiply: {
      */
     (b: Matrix4x4): (a: Matrix4x4) => Matrix4x4;
 };
+export declare const inverse: {
+    /**
+     * Returns the inverse of an `Invertible`-branded `Matrix4x4`. Because the
+     * brand guarantees non-singularity, the result is always a single matrix —
+     * wrapped in `Solution.One`.
+     *
+     * @param m - The matrix to invert.
+     * @since 2.0.0
+     */
+    <T extends Invertible>(m: Matrix4x4<T>): Solution.One<Matrix4x4<T>>;
+    /**
+     * Returns the inverse of a `Matrix4x4`. Returns `Solution.none` if the
+     * matrix is singular (determinant approximately zero). Use {@link asInvertible}
+     * or {@link isInvertible} to refine the type and tighten the return to
+     * `Solution.One`, or use {@link inverseUnsafe} if you want the throwing form.
+     *
+     * @param m - The matrix to invert.
+     * @since 2.0.0
+     */
+    <T>(m: Matrix4x4<T>): Solution.AtMostOne<Matrix4x4<T>>;
+};
 /**
- * Returns the inverse of a `Matrix4x4`.
- *
- * Throws when the matrix is singular (determinant is zero).
+ * Returns the inverse of a `Matrix4x4`. Throws when the matrix is singular
+ * (determinant approximately zero). Useful when you've already validated
+ * invertibility out-of-band and want the unwrapped matrix directly; otherwise
+ * prefer {@link inverse}.
  *
  * @param m - The matrix to invert.
  * @returns The inverse matrix.
- * @since 1.1.0
+ * @throws If `m` is singular.
+ * @since 2.0.0
+ */
+export declare const inverseUnsafe: (m: Matrix4x4) => Matrix4x4;
+/**
+ * Type-narrowing predicate: refines a `Matrix4x4<T>` to
+ * `Matrix4x4<T & Invertible>` when the matrix's determinant is non-zero
+ * (within the default absolute tolerance).
+ *
+ * @param m - The matrix to check.
+ * @returns `true` when `m` is invertible.
+ * @since 2.0.0
+ */
+export declare const isInvertible: <T>(m: Matrix4x4<T>) => m is Matrix4x4<T & Invertible>;
+/**
+ * Asserts that the matrix is invertible, throwing on a singular matrix.
+ *
+ * @param m - The matrix to assert against.
+ * @returns The same matrix, typed with the `Invertible` brand.
+ * @throws When `m` is singular.
+ * @since 2.0.0
  */
-export declare const inverse: (m: Matrix4x4) => Matrix4x4;
+export declare const asInvertible: <T>(m: Matrix4x4<T>) => Matrix4x4<T & Invertible>;
 export declare const solveSystem: {
     /**
      * Solves a system of linear equations represented by a `Matrix4x4` and a `Vector4`.

package/dist/matrix/traits.d.ts

@@ -0,0 +1,34 @@
+/**
+ * Trait brands for matrices.
+ *
+ * Brands are phantom types — they exist only at the type level and carry no
+ * runtime data. A matrix value at runtime is identical regardless of the
+ * traits attached to its type. Refiners (`isInvertible`) and assertions
+ * (`asInvertible`) on each matrix module attach these brands to the matrix's
+ * `Traits` parameter at the type level.
+ *
+ * @since 2.0.0
+ */
+declare const InvertibleId: unique symbol;
+/**
+ * Brand for invertible matrices — matrices whose determinant is non-zero
+ * (within the default absolute tolerance). An `Invertible`-branded matrix
+ * narrows {@link inverse} from `Solution.AtMostOne` to `Solution.One`.
+ *
+ * @since 2.0.0
+ */
+export type Invertible = {
+    readonly [InvertibleId]: 'invertible';
+};
+/**
+ * Phantom property key used by matrix interfaces to hold their `Traits` type
+ * parameter. The runtime object never has a value at this key — the `declare`
+ * modifier on the impl class makes the property type-visible only. The symbol
+ * itself is exported as a real value so that bundlers can resolve the
+ * computed key reference even though the property is never assigned.
+ *
+ * @since 2.0.0
+ */
+export declare const MatrixTraits: unique symbol;
+export type MatrixTraits = typeof MatrixTraits;
+export {};

package/dist/monotonicity/index.mjs

@@ -1,2 +1,2 @@
-import { r as monotonicity_exports } from "../shared/monotonicity-BdebFOtv.mjs";
+import { r as monotonicity_exports } from "../shared/monotonicity-DNxBoC7F.mjs";
 export { monotonicity_exports as Monotonicity };

package/dist/monotonicity/monotonicity.d.ts

@@ -78,3 +78,28 @@ export declare const isStrict: (m: Monotonicity) => m is 1 | 2;
  * @since 2.0.0
  */
 export declare const fromComparison: (s0: number, s1: number) => GuaranteedMonotonicity;
+/**
+ * Classifies monotonicity from the extreme values of a function's derivative
+ * over an interval, tolerating rounding noise in the derivative.
+ *
+ * `min` and `max` must be the true range extremes of the derivative on the
+ * interval (for a polynomial derivative: its values at the interval endpoints
+ * and at any interior stationary point). Each bound sets its sign-coverage
+ * bit only when it exceeds the tolerance, so a derivative that merely grazes
+ * zero — or dips past it by an ulp-scale amount that construction arithmetic
+ * cannot distinguish from zero — does not register as a crossing.
+ *
+ * The default tolerance is {@link RELATIVE_TOLERANCE} scaled by the larger
+ * bound magnitude, making the classification invariant under scaling of the
+ * input data. Derivative values are used rather than derivative *root
+ * locations* because values are well-conditioned where roots are not: a root
+ * computed from rounded coefficients can land on either side of an interval
+ * boundary, while the value it corresponds to stays pinned near zero.
+ *
+ * @param min - The smallest derivative value attained on the interval.
+ * @param max - The largest derivative value attained on the interval.
+ * @param tolerance - Absolute threshold below which a derivative value is treated as zero. Defaults to `RELATIVE_TOLERANCE * max(|min|, |max|)`.
+ * @returns The sign-coverage classification.
+ * @since 2.0.0
+ */
+export declare const fromDerivativeRange: (min: number, max: number, tolerance?: number) => Monotonicity;

package/dist/number.d.ts

@@ -4,9 +4,32 @@
  * @since 1.1.0
  */
 export declare const EPSILON = 1e-10;
+/**
+ * Default relative tolerance for sign and degeneracy decisions about
+ * *computed* quantities — derivative values, determinants, discriminants.
+ *
+ * Scale-free: multiply by the natural magnitude of the quantity under test
+ * (e.g. the largest derivative value over the interval) to get an absolute
+ * threshold. Values inside the resulting band are indistinguishable from
+ * accumulated rounding noise and must not decide a sign question.
+ *
+ * The value sits ~4,500× above IEEE double rounding noise (2⁻⁵² ≈ 2.2e-16
+ * relative per operation) — generous headroom for noise accumulated across a
+ * construction pipeline — and far below any geometrically meaningful feature.
+ *
+ * @since 2.0.0
+ */
+export declare const RELATIVE_TOLERANCE = 1e-12;
 /**
  * Approximate equality between two numbers, within an absolute tolerance.
  *
+ * This is the raw absolute comparison — the right tool when the domain is
+ * normalized by construction (curve parameter space, where `[0, 1]` makes an
+ * absolute band meaningful) or when a caller-chosen absolute band is the
+ * actual question. For identity questions about computed values of arbitrary
+ * magnitude ("is this the same point / value?"), use {@link coincident},
+ * which adds a magnitude-relative term.
+ *
  * @param a - The first value.
  * @param b - The second value.
  * @param eps - The maximum allowed absolute difference. Defaults to {@link EPSILON}.
@@ -14,6 +37,36 @@ export declare const EPSILON = 1e-10;
  * @since 1.1.0
  */
 export declare function epsEquals(a: number, b: number, eps?: number): boolean;
+/**
+ * Library-policy approximate equality: `true` when
+ * `|a - b| <= absolute + relative * max(|a|, |b|)`.
+ *
+ * The relative term covers rounding noise accumulated by arithmetic, which
+ * grows with the magnitude of the values; the absolute term covers the
+ * neighborhood of zero, where relative comparison degenerates (no nonzero
+ * value is relatively close to `0`). The band is strictly wider than the
+ * plain absolute comparison, never narrower.
+ *
+ * Defaults: `absolute = ` {@link EPSILON} — a modeling floor that assumes
+ * coordinates of order ~1; pass an explicit value when your data has a
+ * different natural resolution — and `relative = ` {@link RELATIVE_TOLERANCE}.
+ *
+ * Not transitive: `coincident(a, b)` and `coincident(b, c)` do not imply
+ * `coincident(a, c)` — chains of pairwise-coincident values can drift beyond
+ * the band. Pairwise semantics are intentional (e.g. path continuity checks
+ * each junction locally).
+ *
+ * The tolerance parameters are positional, mirroring {@link epsEquals}: the
+ * common override is `absolute` alone, for a domain with a known resolution.
+ *
+ * @param a - The first value.
+ * @param b - The second value.
+ * @param absolute - Absolute floor of the band. Defaults to {@link EPSILON}.
+ * @param relative - Magnitude-proportional term of the band. Defaults to {@link RELATIVE_TOLERANCE}.
+ * @returns `true` when the values are equal to within the combined band.
+ * @since 2.0.0
+ */
+export declare function coincident(a: number, b: number, absolute?: number, relative?: number): boolean;
 /**
  * Snaps near-zero values to exactly zero for stable sign tests.
  *

package/dist/number.mjs

@@ -8,8 +8,31 @@ const PRECISION = 8;
 */
 const EPSILON = 1e-10;
 /**
+* Default relative tolerance for sign and degeneracy decisions about
+* *computed* quantities — derivative values, determinants, discriminants.
+*
+* Scale-free: multiply by the natural magnitude of the quantity under test
+* (e.g. the largest derivative value over the interval) to get an absolute
+* threshold. Values inside the resulting band are indistinguishable from
+* accumulated rounding noise and must not decide a sign question.
+*
+* The value sits ~4,500× above IEEE double rounding noise (2⁻⁵² ≈ 2.2e-16
+* relative per operation) — generous headroom for noise accumulated across a
+* construction pipeline — and far below any geometrically meaningful feature.
+*
+* @since 2.0.0
+*/
+const RELATIVE_TOLERANCE = 1e-12;
+/**
 * Approximate equality between two numbers, within an absolute tolerance.
 *
+* This is the raw absolute comparison — the right tool when the domain is
+* normalized by construction (curve parameter space, where `[0, 1]` makes an
+* absolute band meaningful) or when a caller-chosen absolute band is the
+* actual question. For identity questions about computed values of arbitrary
+* magnitude ("is this the same point / value?"), use {@link coincident},
+* which adds a magnitude-relative term.
+*
 * @param a - The first value.
 * @param b - The second value.
 * @param eps - The maximum allowed absolute difference. Defaults to {@link EPSILON}.
@@ -20,6 +43,38 @@ function epsEquals(a, b, eps = EPSILON) {
 	return Math.abs(a - b) <= eps;
 }
 /**
+* Library-policy approximate equality: `true` when
+* `|a - b| <= absolute + relative * max(|a|, |b|)`.
+*
+* The relative term covers rounding noise accumulated by arithmetic, which
+* grows with the magnitude of the values; the absolute term covers the
+* neighborhood of zero, where relative comparison degenerates (no nonzero
+* value is relatively close to `0`). The band is strictly wider than the
+* plain absolute comparison, never narrower.
+*
+* Defaults: `absolute = ` {@link EPSILON} — a modeling floor that assumes
+* coordinates of order ~1; pass an explicit value when your data has a
+* different natural resolution — and `relative = ` {@link RELATIVE_TOLERANCE}.
+*
+* Not transitive: `coincident(a, b)` and `coincident(b, c)` do not imply
+* `coincident(a, c)` — chains of pairwise-coincident values can drift beyond
+* the band. Pairwise semantics are intentional (e.g. path continuity checks
+* each junction locally).
+*
+* The tolerance parameters are positional, mirroring {@link epsEquals}: the
+* common override is `absolute` alone, for a domain with a known resolution.
+*
+* @param a - The first value.
+* @param b - The second value.
+* @param absolute - Absolute floor of the band. Defaults to {@link EPSILON}.
+* @param relative - Magnitude-proportional term of the band. Defaults to {@link RELATIVE_TOLERANCE}.
+* @returns `true` when the values are equal to within the combined band.
+* @since 2.0.0
+*/
+function coincident(a, b, absolute = EPSILON, relative = RELATIVE_TOLERANCE) {
+	return Math.abs(a - b) <= absolute + relative * Math.max(Math.abs(a), Math.abs(b));
+}
+/**
 * Snaps near-zero values to exactly zero for stable sign tests.
 *
 * @param value - The value to test.
@@ -135,6 +190,6 @@ function mod(n, m) {
 	return (n % m + m) % m;
 }
 //#endregion
-export { EPSILON, clamp, clampToZero, clip, epsEquals, lerp, minMax, mod, normalize, remap, round, roundDown, roundUp };
+export { EPSILON, RELATIVE_TOLERANCE, clamp, clampToZero, clip, coincident, epsEquals, lerp, minMax, mod, normalize, remap, round, roundDown, roundUp };
 
 //# sourceMappingURL=number.mjs.map