@lakuna/umath 1.0.0 → 1.1.0

package/dist/index.d.ts

@@ -14,6 +14,8 @@ export { type Matrix2Like, default as Matrix2 } from "@lakuna/umath/Matrix2";
 export { type Matrix3Like, default as Matrix3 } from "@lakuna/umath/Matrix3";
 export { type Matrix4Like, default as Matrix4 } from "@lakuna/umath/Matrix4";
 export { type QuaternionLike, default as Quaternion } from "@lakuna/umath/Quaternion";
+export { default as SlowMatrix } from "@lakuna/umath/SlowMatrix";
+export { default as SlowSquareMatrix } from "@lakuna/umath/SlowSquareMatrix";
 export type { default as SquareMatrix } from "./linear/SquareMatrix.js";
 export type { VectorLike, default as Vector } from "./linear/Vector.js";
 export { type Vector2Like, default as Vector2 } from "@lakuna/umath/Vector2";
@@ -25,5 +27,7 @@ export type { IntegerRepresentation } from "./types/IntegerRepresentation.js";
 export { default as BigNumber } from "./utility/BigNumber.js";
 export { default as epsilon } from "./utility/epsilon.js";
 export { default as MagnitudeError } from "./utility/MagnitudeError.js";
+export { default as MatrixSizeError } from "./utility/MatrixSizeError.js";
+export { default as PartialMatrixError } from "./utility/PartialMatrixError.js";
 export { default as SingularMatrixError } from "./utility/SingularMatrixError.js";
 //# sourceMappingURL=index.d.ts.map

package/dist/index.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,qBAAqB,EAAE,MAAM,uCAAuC,CAAC;AACzF,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,oCAAoC,CAAC;AACnF,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,KAAK,kBAAkB,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAClG,YAAY,EAAE,UAAU,EAAE,OAAO,IAAI,MAAM,EAAE,MAAM,oBAAoB,CAAC;AACxE,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,cAAc,EAAE,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtF,YAAY,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,0BAA0B,CAAC;AACxE,YAAY,EAAE,UAAU,EAAE,OAAO,IAAI,MAAM,EAAE,MAAM,oBAAoB,CAAC;AACxE,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,YAAY,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,sBAAsB,CAAC;AACjE,YAAY,EAAE,OAAO,IAAI,WAAW,EAAE,MAAM,wBAAwB,CAAC;AACrE,YAAY,EAAE,qBAAqB,EAAE,MAAM,kCAAkC,CAAC;AAC9E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,IAAI,mBAAmB,EAAE,MAAM,kCAAkC,CAAC"}
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,qBAAqB,EAAE,MAAM,uCAAuC,CAAC;AACzF,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,oCAAoC,CAAC;AACnF,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,KAAK,kBAAkB,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAClG,YAAY,EAAE,UAAU,EAAE,OAAO,IAAI,MAAM,EAAE,MAAM,oBAAoB,CAAC;AACxE,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,cAAc,EAAE,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtF,OAAO,EAAE,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAC7E,YAAY,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,0BAA0B,CAAC;AACxE,YAAY,EAAE,UAAU,EAAE,OAAO,IAAI,MAAM,EAAE,MAAM,oBAAoB,CAAC;AACxE,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAE,KAAK,WAAW,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,YAAY,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,sBAAsB,CAAC;AACjE,YAAY,EAAE,OAAO,IAAI,WAAW,EAAE,MAAM,wBAAwB,CAAC;AACrE,YAAY,EAAE,qBAAqB,EAAE,MAAM,kCAAkC,CAAC;AAC9E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,IAAI,eAAe,EAAE,MAAM,8BAA8B,CAAC;AAC1E,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,iCAAiC,CAAC;AAChF,OAAO,EAAE,OAAO,IAAI,mBAAmB,EAAE,MAAM,kCAAkC,CAAC"}

package/dist/index.js

@@ -13,11 +13,15 @@ export { default as Matrix2 } from "@lakuna/umath/Matrix2";
 export { default as Matrix3 } from "@lakuna/umath/Matrix3";
 export { default as Matrix4 } from "@lakuna/umath/Matrix4";
 export { default as Quaternion } from "@lakuna/umath/Quaternion";
+export { default as SlowMatrix } from "@lakuna/umath/SlowMatrix";
+export { default as SlowSquareMatrix } from "@lakuna/umath/SlowSquareMatrix";
 export { default as Vector2 } from "@lakuna/umath/Vector2";
 export { default as Vector3 } from "@lakuna/umath/Vector3";
 export { default as Vector4 } from "@lakuna/umath/Vector4";
 export { default as BigNumber } from "./utility/BigNumber.js";
 export { default as epsilon } from "./utility/epsilon.js";
 export { default as MagnitudeError } from "./utility/MagnitudeError.js";
+export { default as MatrixSizeError } from "./utility/MatrixSizeError.js";
+export { default as PartialMatrixError } from "./utility/PartialMatrixError.js";
 export { default as SingularMatrixError } from "./utility/SingularMatrixError.js";
 //# sourceMappingURL=index.js.map

package/dist/index.js.map

@@ -1 +1 @@
-{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,qBAAqB,EAAE,MAAM,uCAAuC,CAAC;AACzF,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,oCAAoC,CAAC;AACnF,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAA2B,OAAO,IAAI,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAElG,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAuB,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AAGtF,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAI7E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,IAAI,mBAAmB,EAAE,MAAM,kCAAkC,CAAC"}
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../src/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,qBAAqB,EAAE,MAAM,uCAAuC,CAAC;AACzF,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,yBAAyB,CAAC;AAC7D,OAAO,EAAE,OAAO,IAAI,YAAY,EAAE,MAAM,8BAA8B,CAAC;AACvE,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,oCAAoC,CAAC;AACnF,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AAC/E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,2BAA2B,CAAC;AACjE,OAAO,EAA2B,OAAO,IAAI,cAAc,EAAE,MAAM,8BAA8B,CAAC;AAElG,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAuB,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACtF,OAAO,EAAE,OAAO,IAAI,UAAU,EAAE,MAAM,0BAA0B,CAAC;AACjE,OAAO,EAAE,OAAO,IAAI,gBAAgB,EAAE,MAAM,gCAAgC,CAAC;AAG7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAC7E,OAAO,EAAoB,OAAO,IAAI,OAAO,EAAE,MAAM,uBAAuB,CAAC;AAI7E,OAAO,EAAE,OAAO,IAAI,SAAS,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,OAAO,IAAI,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC1D,OAAO,EAAE,OAAO,IAAI,cAAc,EAAE,MAAM,6BAA6B,CAAC;AACxE,OAAO,EAAE,OAAO,IAAI,eAAe,EAAE,MAAM,8BAA8B,CAAC;AAC1E,OAAO,EAAE,OAAO,IAAI,kBAAkB,EAAE,MAAM,iCAAiC,CAAC;AAChF,OAAO,EAAE,OAAO,IAAI,mBAAmB,EAAE,MAAM,kCAAkC,CAAC"}

package/dist/linear/Matrix.d.ts

@@ -1,14 +1,13 @@
 export type MatrixLike = Matrix | ArrayLike<number>;
 export default interface Matrix extends ArrayLike<number> {
-    width: number;
     height: number;
-    equals(matrix: MatrixLike): boolean;
-    exactEquals(matrix: MatrixLike): boolean;
+    width: number;
+    get frob(): number;
     add(matrix: MatrixLike): Matrix;
-    adjoint(): Matrix;
     clone(): Matrix;
     copy(matrix: MatrixLike): this;
-    get frob(): number;
+    equals(matrix: MatrixLike): boolean;
+    exactEquals(matrix: MatrixLike): boolean;
     multiply(matrix: MatrixLike): Matrix;
     multiplyScalar(scalar: number): Matrix;
     multiplyScalarAndAdd(matrix: MatrixLike, scalar: number): Matrix;

package/dist/linear/Matrix.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Matrix.d.ts","sourceRoot":"","sources":["../../src/linear/Matrix.ts"],"names":[],"mappings":"AAIA,MAAM,MAAM,UAAU,GAAG,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;AAMpD,MAAM,CAAC,OAAO,WAAW,MAAO,SAAQ,SAAS,CAAC,MAAM,CAAC;IAExD,KAAK,EAAE,MAAM,CAAC;IAGd,MAAM,EAAE,MAAM,CAAC;IAOf,MAAM,CAAC,MAAM,EAAE,UAAU,GAAG,OAAO,CAAC;IAOpC,WAAW,CAAC,MAAM,EAAE,UAAU,GAAG,OAAO,CAAC;IAQzC,GAAG,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAOhC,OAAO,IAAI,MAAM,CAAC;IAMlB,KAAK,IAAI,MAAM,CAAC;IAOhB,IAAI,CAAC,MAAM,EAAE,UAAU,GAAG,IAAI,CAAC;IAM/B,IAAI,IAAI,IAAI,MAAM,CAAC;IAQnB,QAAQ,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAQrC,cAAc,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC;IAUvC,oBAAoB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC;IAQjE,QAAQ,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAOrC,SAAS,IAAI,MAAM,CAAC;CACpB"}
+{"version":3,"file":"Matrix.d.ts","sourceRoot":"","sources":["../../src/linear/Matrix.ts"],"names":[],"mappings":"AAIA,MAAM,MAAM,UAAU,GAAG,MAAM,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;AAMpD,MAAM,CAAC,OAAO,WAAW,MAAO,SAAQ,SAAS,CAAC,MAAM,CAAC;IAExD,MAAM,EAAE,MAAM,CAAC;IAGf,KAAK,EAAE,MAAM,CAAC;IAMd,IAAI,IAAI,IAAI,MAAM,CAAC;IAQnB,GAAG,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAMhC,KAAK,IAAI,MAAM,CAAC;IAOhB,IAAI,CAAC,MAAM,EAAE,UAAU,GAAG,IAAI,CAAC;IAO/B,MAAM,CAAC,MAAM,EAAE,UAAU,GAAG,OAAO,CAAC;IAOpC,WAAW,CAAC,MAAM,EAAE,UAAU,GAAG,OAAO,CAAC;IAQzC,QAAQ,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAQrC,cAAc,CAAC,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC;IAUvC,oBAAoB,CAAC,MAAM,EAAE,UAAU,EAAE,MAAM,EAAE,MAAM,GAAG,MAAM,CAAC;IAQjE,QAAQ,CAAC,MAAM,EAAE,UAAU,GAAG,MAAM,CAAC;IAOrC,SAAS,IAAI,MAAM,CAAC;CACpB"}

package/dist/linear/SlowMatrix.d.ts

@@ -0,0 +1,18 @@
+import { type Matrix } from "@lakuna/umath";
+export default class SlowMatrix extends Float32Array implements Matrix {
+    constructor(...cols: Array<Array<number>>);
+    readonly height: number;
+    readonly width: number;
+    get frob(): number;
+    add(matrix: Matrix): SlowMatrix;
+    clone(): SlowMatrix;
+    copy(matrix: Matrix): this;
+    equals(matrix: Matrix): boolean;
+    exactEquals(matrix: Matrix): boolean;
+    multiply(matrix: Matrix): SlowMatrix;
+    multiplyScalar(scalar: number): SlowMatrix;
+    multiplyScalarAndAdd(matrix: Matrix, scalar: number): SlowMatrix;
+    subtract(matrix: Matrix): SlowMatrix;
+    transpose(): SlowMatrix;
+}
+//# sourceMappingURL=SlowMatrix.d.ts.map

package/dist/linear/SlowMatrix.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SlowMatrix.d.ts","sourceRoot":"","sources":["../../src/linear/SlowMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,MAAM,EAAgD,MAAM,eAAe,CAAC;AAM1F,MAAM,CAAC,OAAO,OAAO,UAAW,SAAQ,YAAa,YAAW,MAAM;gBAKlD,GAAG,IAAI,EAAE,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IAkBhD,SAAgB,MAAM,EAAE,MAAM,CAAC;IAG/B,SAAgB,KAAK,EAAE,MAAM,CAAC;IAM9B,IAAW,IAAI,IAAI,MAAM,CAExB;IAQM,GAAG,CAAC,MAAM,EAAE,MAAM,GAAG,UAAU;IAmB/B,KAAK,IAAI,UAAU;IAgBnB,IAAI,CAAC,MAAM,EAAE,MAAM,GAAG,IAAI;IAmB1B,MAAM,CAAC,MAAM,EAAE,MAAM,GAAG,OAAO;IAmB/B,WAAW,CAAC,MAAM,EAAE,MAAM,GAAG,OAAO;IAoBpC,QAAQ,CAAC,MAAM,EAAE,MAAM,GAAG,UAAU;IA6BpC,cAAc,CAAC,MAAM,EAAE,MAAM,GAAG,UAAU;IAiB1C,oBAAoB,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,GAAG,UAAU;IAUhE,QAAQ,CAAC,MAAM,EAAE,MAAM,GAAG,UAAU;IAoBpC,SAAS,IAAI,UAAU;CAU9B"}

package/dist/linear/SlowMatrix.js

@@ -0,0 +1,126 @@
+import { PartialMatrixError, epsilon, MatrixSizeError } from "@lakuna/umath";
+export default class SlowMatrix extends Float32Array {
+    constructor(...cols) {
+        super(cols.flat());
+        this.width = cols.length;
+        this.height = this.width && cols[0].length;
+        if (this.width < 1 || this.height < 1) {
+            throw new MatrixSizeError();
+        }
+        for (let i = 0; i < this.width; i++) {
+            if ((cols[i]?.length ?? 0) != this.height) {
+                throw new PartialMatrixError();
+            }
+        }
+    }
+    height;
+    width;
+    get frob() {
+        return Math.hypot(...this);
+    }
+    add(matrix) {
+        if (this.width != matrix.width || this.height != matrix.height) {
+            throw new MatrixSizeError();
+        }
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[x] ??= [])[y] = this[x * this.height + y] + matrix[x * matrix.height + y];
+            }
+        }
+        return new SlowMatrix(...cols);
+    }
+    clone() {
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[x] ??= [])[y] = this[x * this.height + y];
+            }
+        }
+        return new SlowMatrix(...cols);
+    }
+    copy(matrix) {
+        if (this.width != matrix.width || this.height != matrix.height) {
+            throw new MatrixSizeError();
+        }
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                this[x * this.width + y] = matrix[x * this.height + y];
+            }
+        }
+        return this;
+    }
+    equals(matrix) {
+        if (this.width != matrix.width || this.height != matrix.height) {
+            return false;
+        }
+        for (let i = 0; i < this.length; i++) {
+            if (Math.abs(this[i] - matrix[i]) > epsilon) {
+                return false;
+            }
+        }
+        return true;
+    }
+    exactEquals(matrix) {
+        if (this.width != matrix.width || this.height != matrix.height) {
+            return false;
+        }
+        for (let i = 0; i < this.length; i++) {
+            if (this[i] != matrix[i]) {
+                return false;
+            }
+        }
+        return true;
+    }
+    multiply(matrix) {
+        if (this.width != matrix.height) {
+            throw new MatrixSizeError();
+        }
+        const n = this.height;
+        const m = this.width;
+        const p = matrix.width;
+        const out = [];
+        for (let i = 0; i < n; i++) {
+            for (let j = 0; j < p; j++) {
+                let sum = 0;
+                for (let k = 0; k < m; k++) {
+                    sum += this[k * this.height + i] * matrix[j * matrix.height + k];
+                }
+                (out[j] ??= [])[i] = sum;
+            }
+        }
+        return new SlowMatrix(...out);
+    }
+    multiplyScalar(scalar) {
+        const out = this.clone();
+        for (let i = 0; i < out.length; i++) {
+            out[i] *= scalar;
+        }
+        return out;
+    }
+    multiplyScalarAndAdd(matrix, scalar) {
+        return this.add(matrix.multiplyScalar(scalar));
+    }
+    subtract(matrix) {
+        if (this.width != matrix.width || this.height != matrix.height) {
+            throw new MatrixSizeError();
+        }
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[x] ??= [])[y] = this[x * this.height + y] - matrix[x * matrix.height + y];
+            }
+        }
+        return new SlowMatrix(...cols);
+    }
+    transpose() {
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[y] ??= [])[x] = this[x * this.height + y];
+            }
+        }
+        return new SlowMatrix(...cols);
+    }
+}
+//# sourceMappingURL=SlowMatrix.js.map

package/dist/linear/SlowMatrix.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SlowMatrix.js","sourceRoot":"","sources":["../../src/linear/SlowMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAe,kBAAkB,EAAE,OAAO,EAAE,eAAe,EAAE,MAAM,eAAe,CAAC;AAM1F,MAAM,CAAC,OAAO,OAAO,UAAW,SAAQ,YAAY;IAKnD,YAAmB,GAAG,IAA0B;QAC/C,KAAK,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;QACnB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC;QACzB,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,KAAK,IAAK,IAAI,CAAC,CAAC,CAAmB,CAAC,MAAM,CAAC;QAE9D,IAAI,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;YACtC,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAGD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,IAAI,CAAC,CAAC,IAAI,IAAI,CAAC,MAAM,EAAE;gBAC1C,MAAM,IAAI,kBAAkB,EAAE,CAAC;aAC/B;SACD;IACF,CAAC;IAGe,MAAM,CAAS;IAGf,KAAK,CAAS;IAM9B,IAAW,IAAI;QACd,OAAO,IAAI,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,CAAC;IAC5B,CAAC;IAQM,GAAG,CAAC,MAAc;QACxB,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,EAAE;YAC/D,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAED,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAY,GAAI,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAY,CAAC;aACxG;SACD;QAED,OAAO,IAAI,UAAU,CAAC,GAAG,IAAI,CAAC,CAAC;IAChC,CAAC;IAMM,KAAK;QACX,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAW,CAAC;aAC1D;SACD;QAED,OAAO,IAAI,UAAU,CAAC,GAAG,IAAI,CAAC,CAAC;IAChC,CAAC;IAOM,IAAI,CAAC,MAAc;QACzB,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,EAAE;YAC/D,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACpC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,GAAG,CAAC,CAAY,GAAG,MAAM,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAW,CAAC;aAC7E;SACD;QAED,OAAO,IAAI,CAAC;IACb,CAAC;IAOM,MAAM,CAAC,MAAc;QAC3B,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,EAAE;YAC/D,OAAO,KAAK,CAAC;SACb;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACrC,IAAI,IAAI,CAAC,GAAG,CAAE,IAAI,CAAC,CAAC,CAAY,GAAI,MAAM,CAAC,CAAC,CAAY,CAAC,GAAG,OAAO,EAAE;gBACpE,OAAO,KAAK,CAAC;aACb;SACD;QAED,OAAO,IAAI,CAAC;IACb,CAAC;IAOM,WAAW,CAAC,MAAc;QAChC,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,EAAE;YAC/D,OAAO,KAAK,CAAC;SACb;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACrC,IAAI,IAAI,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,EAAE;gBACzB,OAAO,KAAK,CAAC;aACb;SACD;QAED,OAAO,IAAI,CAAC;IACb,CAAC;IAQM,QAAQ,CAAC,MAAc;QAC7B,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,MAAM,EAAE;YAChC,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAED,MAAM,CAAC,GAAW,IAAI,CAAC,MAAM,CAAC;QAC9B,MAAM,CAAC,GAAW,IAAI,CAAC,KAAK,CAAC;QAC7B,MAAM,CAAC,GAAW,MAAM,CAAC,KAAK,CAAC;QAE/B,MAAM,GAAG,GAAyB,EAAE,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC3B,IAAI,GAAG,GAAG,CAAC,CAAC;gBACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;oBAC3B,GAAG,IAAK,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAY,GAAI,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAY,CAAA;iBACxF;gBACD,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;aACzB;SACD;QAED,OAAO,IAAI,UAAU,CAAC,GAAG,GAAG,CAAC,CAAC;IAC/B,CAAC;IAQM,cAAc,CAAC,MAAc;QACnC,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACpC,GAAG,CAAC,CAAC,CAAC,IAAI,MAAM,CAAC;SACjB;QAED,OAAO,GAAG,CAAC;IACZ,CAAC;IAUM,oBAAoB,CAAC,MAAc,EAAE,MAAc;QACzD,OAAO,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC,CAAC;IAChD,CAAC;IAQM,QAAQ,CAAC,MAAc;QAC7B,IAAI,IAAI,CAAC,KAAK,IAAI,MAAM,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,IAAI,MAAM,CAAC,MAAM,EAAE;YAC/D,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAED,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAY,GAAI,MAAM,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,CAAY,CAAC;aACxG;SACD;QAED,OAAO,IAAI,UAAU,CAAC,GAAG,IAAI,CAAC,CAAC;IAChC,CAAC;IAOM,SAAS;QACf,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAY,CAAC;aAC5D;SACD;QAED,OAAO,IAAI,UAAU,CAAC,GAAG,IAAI,CAAC,CAAC;IAChC,CAAC;CACD"}

package/dist/linear/SlowSquareMatrix.d.ts

@@ -0,0 +1,15 @@
+import { type SquareMatrix } from "@lakuna/umath";
+import SlowMatrix from "@lakuna/umath/SlowMatrix";
+export default class SlowSquareMatrix extends SlowMatrix implements SquareMatrix {
+    constructor(...cols: Array<Array<number>>);
+    get determinant(): number;
+    adjoint(): SlowSquareMatrix;
+    clone(): SlowSquareMatrix;
+    cofactor(): SlowSquareMatrix;
+    identity(): this;
+    invert(): SlowSquareMatrix;
+    minor(row: number, col: number): number;
+    submatrix(row: number, col: number): SlowSquareMatrix;
+    transpose(): SlowSquareMatrix;
+}
+//# sourceMappingURL=SlowSquareMatrix.d.ts.map

package/dist/linear/SlowSquareMatrix.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SlowSquareMatrix.d.ts","sourceRoot":"","sources":["../../src/linear/SlowSquareMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,KAAK,YAAY,EAAwC,MAAM,eAAe,CAAC;AACxF,OAAO,UAAU,MAAM,0BAA0B,CAAC;AAMlD,MAAM,CAAC,OAAO,OAAO,gBAAiB,SAAQ,UAAW,YAAW,YAAY;gBAM5D,GAAG,IAAI,EAAE,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IAehD,IAAW,WAAW,IAAI,MAAM,CAe/B;IAOM,OAAO,IAAI,gBAAgB;IAQlB,KAAK,IAAI,gBAAgB;IAgBlC,QAAQ,IAAI,gBAAgB;IAgB5B,QAAQ,IAAI,IAAI;IAehB,MAAM,IAAI,gBAAgB;IAoD1B,KAAK,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,MAAM;IAUvC,SAAS,CAAC,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,GAAG,gBAAgB;IAyB5C,SAAS,IAAI,gBAAgB;CAU7C"}

package/dist/linear/SlowSquareMatrix.js

@@ -0,0 +1,118 @@
+import { MatrixSizeError, SingularMatrixError } from "@lakuna/umath";
+import SlowMatrix from "@lakuna/umath/SlowMatrix";
+export default class SlowSquareMatrix extends SlowMatrix {
+    constructor(...cols) {
+        super(...cols);
+        for (let i = 0; i < this.width; i++) {
+            if (cols[i].length != this.height) {
+                throw new MatrixSizeError();
+            }
+        }
+    }
+    get determinant() {
+        if (this.length < 1) {
+            throw new MatrixSizeError();
+        }
+        if (this.length == 1) {
+            return this[0] ?? 0;
+        }
+        let out = 0;
+        for (let i = 0; i < this.width; i++) {
+            out += this[i * this.height] * (i % 2 ? -1 : 1) * this.minor(0, i);
+        }
+        return out;
+    }
+    adjoint() {
+        return this.cofactor().transpose();
+    }
+    clone() {
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[x] ??= [])[y] = this[x * this.height + y];
+            }
+        }
+        return new SlowSquareMatrix(...cols);
+    }
+    cofactor() {
+        const out = [];
+        for (let i = 0; i < this.width; i++) {
+            for (let j = 0; j < this.height; j++) {
+                (out[j] ??= [])[i] = ((i + j) % 2 ? -1 : 1) * this.minor(i, j);
+            }
+        }
+        return new SlowSquareMatrix(...out);
+    }
+    identity() {
+        for (let i = 0; i < this.width; i++) {
+            for (let j = 0; j < this.height; j++) {
+                this[i * this.height + j] = i == j ? 1 : 0;
+            }
+        }
+        return this;
+    }
+    invert() {
+        const dim = this.width;
+        const clone = this.clone();
+        const identity = this.clone().identity();
+        for (let i = 0; i < dim; i++) {
+            let diagonal = clone[i * dim + i];
+            if (!diagonal) {
+                for (let ii = i + 1; ii < dim; ii++) {
+                    if (clone[ii * dim + i]) {
+                        for (let j = 0; j < dim; j++) {
+                            for (const matrix of [clone, identity]) {
+                                [matrix[i * dim + j], matrix[ii * dim + j]] = [matrix[ii * dim + j], matrix[i * dim + j]];
+                            }
+                        }
+                        break;
+                    }
+                }
+                diagonal = clone[i * dim + i];
+                if (!diagonal) {
+                    throw new SingularMatrixError();
+                }
+            }
+            for (let ii = 0; ii < dim; ii++) {
+                if (ii == i) {
+                    continue;
+                }
+                const temp = clone[ii * dim + i];
+                for (let j = 0; j < dim; j++) {
+                    for (const matrix of [clone, identity]) {
+                        matrix[ii * dim + j] -= temp * matrix[i * dim + j];
+                    }
+                }
+            }
+        }
+        return identity;
+    }
+    minor(row, col) {
+        return this.submatrix(row, col).determinant;
+    }
+    submatrix(row, col) {
+        const cols = [];
+        for (let i = 0; i < this.width; i++) {
+            if (i == col) {
+                continue;
+            }
+            for (let j = 0; j < this.height; j++) {
+                if (j == row) {
+                    continue;
+                }
+                (cols[i < col ? i : i - 1] ??= [])[j < row ? j : j - 1] = this[i * this.height + j];
+            }
+        }
+        return new SlowSquareMatrix(...cols);
+    }
+    transpose() {
+        const cols = [];
+        for (let x = 0; x < this.width; x++) {
+            for (let y = 0; y < this.height; y++) {
+                (cols[y] ??= [])[x] = this[x * this.height + y];
+            }
+        }
+        return new SlowSquareMatrix(...cols);
+    }
+}
+//# sourceMappingURL=SlowSquareMatrix.js.map

package/dist/linear/SlowSquareMatrix.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SlowSquareMatrix.js","sourceRoot":"","sources":["../../src/linear/SlowSquareMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,EAAqB,eAAe,EAAE,mBAAmB,EAAE,MAAM,eAAe,CAAC;AACxF,OAAO,UAAU,MAAM,0BAA0B,CAAC;AAMlD,MAAM,CAAC,OAAO,OAAO,gBAAiB,SAAQ,UAAU;IAMvD,YAAmB,GAAG,IAA0B;QAC/C,KAAK,CAAC,GAAG,IAAI,CAAC,CAAC;QAGf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,IAAK,IAAI,CAAC,CAAC,CAAmB,CAAC,MAAM,IAAI,IAAI,CAAC,MAAM,EAAE;gBACrD,MAAM,IAAI,eAAe,EAAE,CAAC;aAC5B;SACD;IACF,CAAC;IAMD,IAAW,WAAW;QACrB,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;YACpB,MAAM,IAAI,eAAe,EAAE,CAAC;SAC5B;QAED,IAAI,IAAI,CAAC,MAAM,IAAI,CAAC,EAAE;YACrB,OAAO,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;SACpB;QAED,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,GAAG,IAAK,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,CAAY,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;SAC/E;QAED,OAAO,GAAG,CAAC;IACZ,CAAC;IAOM,OAAO;QACb,OAAO,IAAI,CAAC,QAAQ,EAAE,CAAC,SAAS,EAAE,CAAC;IACpC,CAAC;IAMe,KAAK;QACpB,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAW,CAAC;aAC1D;SACD;QAED,OAAO,IAAI,gBAAgB,CAAC,GAAG,IAAI,CAAC,CAAC;IACtC,CAAC;IAOM,QAAQ;QACd,MAAM,GAAG,GAAyB,EAAE,CAAC;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;aAC/D;SACD;QAED,OAAO,IAAI,gBAAgB,CAAC,GAAG,GAAG,CAAC,CAAC;IACrC,CAAC;IAOM,QAAQ;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aAC3C;SACD;QAED,OAAO,IAAI,CAAC;IACb,CAAC;IAOM,MAAM;QACZ,MAAM,GAAG,GAAW,IAAI,CAAC,KAAK,CAAC;QAE/B,MAAM,KAAK,GAAqB,IAAI,CAAC,KAAK,EAAE,CAAC;QAC7C,MAAM,QAAQ,GAAqB,IAAI,CAAC,KAAK,EAAE,CAAC,QAAQ,EAAE,CAAC;QAE3D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE;YAC7B,IAAI,QAAQ,GAAW,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAW,CAAC;YACpD,IAAI,CAAC,QAAQ,EAAE;gBACd,KAAK,IAAI,EAAE,GAAW,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE;oBAC5C,IAAI,KAAK,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE;wBACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE;4BAC7B,KAAK,MAAM,MAAM,IAAI,CAAC,KAAK,EAAE,QAAQ,CAAC,EAAE;gCACvC,CAAC,MAAM,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,CAAW,EAAE,MAAM,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAW,CAAC,CAAC;6BAC9G;yBACD;wBAED,MAAM;qBACN;iBACD;gBAED,QAAQ,GAAG,KAAK,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAW,CAAC;gBACxC,IAAI,CAAC,QAAQ,EAAE;oBACd,MAAM,IAAI,mBAAmB,EAAE,CAAC;iBAChC;aACD;YAED,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE;gBAChC,IAAI,EAAE,IAAI,CAAC,EAAE;oBACZ,SAAS;iBACT;gBAED,MAAM,IAAI,GAAW,KAAK,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,CAAW,CAAC;gBAEnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE;oBAC7B,KAAK,MAAM,MAAM,IAAI,CAAC,KAAK,EAAE,QAAQ,CAAC,EAAE;wBACvC,MAAM,CAAC,EAAE,GAAG,GAAG,GAAG,CAAC,CAAC,IAAI,IAAI,GAAI,MAAM,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAY,CAAC;qBAC/D;iBACD;aACD;SACD;QAED,OAAO,QAAQ,CAAC;IACjB,CAAC;IASM,KAAK,CAAC,GAAW,EAAE,GAAW;QACpC,OAAO,IAAI,CAAC,SAAS,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,WAAW,CAAC;IAC7C,CAAC;IAQM,SAAS,CAAC,GAAW,EAAE,GAAW;QACxC,MAAM,IAAI,GAAyB,EAAE,CAAC;QAEtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,IAAI,CAAC,IAAI,GAAG,EAAE;gBACb,SAAS;aACT;YAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,IAAI,CAAC,IAAI,GAAG,EAAE;oBACb,SAAS;iBACT;gBAED,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAW,CAAC;aAC9F;SACD;QAED,OAAO,IAAI,gBAAgB,CAAC,GAAG,IAAI,CAAC,CAAC;IACtC,CAAC;IAOe,SAAS;QACxB,MAAM,IAAI,GAAyB,EAAE,CAAC;QACtC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE;YACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACrC,CAAC,IAAI,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC,GAAI,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,MAAM,GAAG,CAAC,CAAY,CAAC;aAC5D;SACD;QAED,OAAO,IAAI,gBAAgB,CAAC,GAAG,IAAI,CAAC,CAAC;IACtC,CAAC;CACD"}

package/dist/linear/SquareMatrix.d.ts

@@ -1,6 +1,7 @@
 import type { Matrix } from "@lakuna/umath";
 export default interface SquareMatrix extends Matrix {
     get determinant(): number;
+    adjoint(): Matrix;
     identity(): this;
     invert(): SquareMatrix;
 }

package/dist/linear/SquareMatrix.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"SquareMatrix.d.ts","sourceRoot":"","sources":["../../src/linear/SquareMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,eAAe,CAAC;AAM5C,MAAM,CAAC,OAAO,WAAW,YAAa,SAAQ,MAAM;IAKnD,IAAI,WAAW,IAAI,MAAM,CAAC;IAO1B,QAAQ,IAAI,IAAI,CAAC;IAOjB,MAAM,IAAI,YAAY,CAAC;CACvB"}
+{"version":3,"file":"SquareMatrix.d.ts","sourceRoot":"","sources":["../../src/linear/SquareMatrix.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,MAAM,EAAE,MAAM,eAAe,CAAC;AAM5C,MAAM,CAAC,OAAO,WAAW,YAAa,SAAQ,MAAM;IAKnD,IAAI,WAAW,IAAI,MAAM,CAAC;IAO1B,OAAO,IAAI,MAAM,CAAC;IAOlB,QAAQ,IAAI,IAAI,CAAC;IAOjB,MAAM,IAAI,YAAY,CAAC;CACvB"}

package/dist/utility/MagnitudeError.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"MagnitudeError.d.ts","sourceRoot":"","sources":["../../src/utility/MagnitudeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,KAAK;gBAK7B,OAAO,SAA6B;CAGvD"}
+{"version":3,"file":"MagnitudeError.d.ts","sourceRoot":"","sources":["../../src/utility/MagnitudeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,KAAK;gBAK7B,OAAO,SAA6B;CAIvD"}

package/dist/utility/MagnitudeError.js

@@ -1,6 +1,7 @@
 export default class MagnitudeError extends Error {
     constructor(message = "The vector is too small.") {
         super(message);
+        this.name = "MagnitudeError";
     }
 }
 //# sourceMappingURL=MagnitudeError.js.map

package/dist/utility/MagnitudeError.js.map

@@ -1 +1 @@
-{"version":3,"file":"MagnitudeError.js","sourceRoot":"","sources":["../../src/utility/MagnitudeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,KAAK;IAKhD,YAAmB,OAAO,GAAG,0BAA0B;QACtD,KAAK,CAAC,OAAO,CAAC,CAAC;IAChB,CAAC;CACD"}
+{"version":3,"file":"MagnitudeError.js","sourceRoot":"","sources":["../../src/utility/MagnitudeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,cAAe,SAAQ,KAAK;IAKhD,YAAmB,OAAO,GAAG,0BAA0B;QACtD,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gBAAgB,CAAC;IAC9B,CAAC;CACD"}

package/dist/utility/MatrixSizeError.d.ts

@@ -0,0 +1,4 @@
+export default class PartialMatrixError extends Error {
+    constructor(message?: string);
+}
+//# sourceMappingURL=MatrixSizeError.d.ts.map

package/dist/utility/MatrixSizeError.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"MatrixSizeError.d.ts","sourceRoot":"","sources":["../../src/utility/MatrixSizeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,kBAAmB,SAAQ,KAAK;gBAKjC,OAAO,SAA+B;CAIzD"}

package/dist/utility/MatrixSizeError.js

@@ -0,0 +1,7 @@
+export default class PartialMatrixError extends Error {
+    constructor(message = "Invalid matrix dimensions.") {
+        super(message);
+        this.name = "MatrixSizeError";
+    }
+}
+//# sourceMappingURL=MatrixSizeError.js.map

package/dist/utility/MatrixSizeError.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"MatrixSizeError.js","sourceRoot":"","sources":["../../src/utility/MatrixSizeError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,kBAAmB,SAAQ,KAAK;IAKpD,YAAmB,OAAO,GAAG,4BAA4B;QACxD,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,iBAAiB,CAAC;IAC/B,CAAC;CACD"}

package/dist/utility/PartialMatrixError.d.ts

@@ -0,0 +1,4 @@
+export default class PartialMatrixError extends Error {
+    constructor(message?: string);
+}
+//# sourceMappingURL=PartialMatrixError.d.ts.map

package/dist/utility/PartialMatrixError.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"PartialMatrixError.d.ts","sourceRoot":"","sources":["../../src/utility/PartialMatrixError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,kBAAmB,SAAQ,KAAK;gBAKjC,OAAO,SAAmC;CAI7D"}

package/dist/utility/PartialMatrixError.js

@@ -0,0 +1,7 @@
+export default class PartialMatrixError extends Error {
+    constructor(message = "The matrix is not rectangular.") {
+        super(message);
+        this.name = "PartialMatrixError";
+    }
+}
+//# sourceMappingURL=PartialMatrixError.js.map

package/dist/utility/PartialMatrixError.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"PartialMatrixError.js","sourceRoot":"","sources":["../../src/utility/PartialMatrixError.ts"],"names":[],"mappings":"AACA,MAAM,CAAC,OAAO,OAAO,kBAAmB,SAAQ,KAAK;IAKpD,YAAmB,OAAO,GAAG,gCAAgC;QAC5D,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,oBAAoB,CAAC;IAClC,CAAC;CACD"}

package/package.json

@@ -1,6 +1,6 @@
 {
 	"name": "@lakuna/umath",
-	"version": "1.0.0",
+	"version": "1.1.0",
 	"description": "A lightweight math library.",
 	"keywords": [
 		"math"
@@ -53,6 +53,8 @@
 		"./Matrix3": "./dist/linear/Matrix3.js",
 		"./Matrix4": "./dist/linear/Matrix4.js",
 		"./Quaternion": "./dist/linear/Quaternion.js",
+		"./SlowMatrix": "./dist/linear/SlowMatrix.js",
+		"./SlowSquareMatrix": "./dist/linear/SlowSquareMatrix.js",
 		"./Vector2": "./dist/linear/Vector2.js",
 		"./Vector3": "./dist/linear/Vector3.js",
 		"./Vector4": "./dist/linear/Vector4.js"

package/src/index.ts

@@ -14,6 +14,8 @@ export { type Matrix2Like, default as Matrix2 } from "@lakuna/umath/Matrix2";
 export { type Matrix3Like, default as Matrix3 } from "@lakuna/umath/Matrix3";
 export { type Matrix4Like, default as Matrix4 } from "@lakuna/umath/Matrix4";
 export { type QuaternionLike, default as Quaternion } from "@lakuna/umath/Quaternion";
+export { default as SlowMatrix } from "@lakuna/umath/SlowMatrix";
+export { default as SlowSquareMatrix } from "@lakuna/umath/SlowSquareMatrix";
 export type { default as SquareMatrix } from "./linear/SquareMatrix.js";
 export type { VectorLike, default as Vector } from "./linear/Vector.js";
 export { type Vector2Like, default as Vector2 } from "@lakuna/umath/Vector2";
@@ -25,4 +27,6 @@ export type { IntegerRepresentation } from "./types/IntegerRepresentation.js";
 export { default as BigNumber } from "./utility/BigNumber.js";
 export { default as epsilon } from "./utility/epsilon.js";
 export { default as MagnitudeError } from "./utility/MagnitudeError.js";
+export { default as MatrixSizeError } from "./utility/MatrixSizeError.js";
+export { default as PartialMatrixError } from "./utility/PartialMatrixError.js";
 export { default as SingularMatrixError } from "./utility/SingularMatrixError.js";

package/src/linear/Matrix.ts

@@ -9,26 +9,18 @@ export type MatrixLike = Matrix | ArrayLike<number>;
  * @see [Matrix](https://en.wikipedia.org/wiki/Matrix_(mathematics))
  */
 export default interface Matrix extends ArrayLike<number> {
-	/** The number of columns in this matrix. */
-	width: number;
-
 	/** The number of rows in this matrix. */
 	height: number;
 
-	/**
-	 * Determines whether this matrix is roughly equivalent to another.
-	 * @param matrix The other matrix.
-	 * @returns Whether the matrices are equivalent.
-	 */
-	equals(matrix: MatrixLike): boolean;
+	/** The number of columns in this matrix. */
+	width: number;
 
 	/**
-	 * Determines whether this matrix is exactly equivalent to another.
-	 * @param matrix The other matrix.
-	 * @returns Whether the matrices are equivalent.
+	 * The Frobenius norm of this matrix.
+	 * @see [Matrix norm](https://en.wikipedia.org/wiki/Matrix_norm)
 	 */
-	exactEquals(matrix: MatrixLike): boolean;
-
+	get frob(): number;
+	
 	/**
 	 * Adds two matrices of the same size.
 	 * @param matrix The other matrix.
@@ -37,13 +29,6 @@ export default interface Matrix extends ArrayLike<number> {
 	 */
 	add(matrix: MatrixLike): Matrix;
 
-	/**
-	 * Calculates the adjugate of this matrix.
-	 * @returns The adjugate of this matrix.
-	 * @see [Adjugate matrix](https://en.wikipedia.org/wiki/Adjugate_matrix)
-	 */
-	adjoint(): Matrix;
-
 	/**
 	 * Creates a copy of this matrix.
 	 * @returns A copy of this matrix.
@@ -58,10 +43,18 @@ export default interface Matrix extends ArrayLike<number> {
 	copy(matrix: MatrixLike): this;
 
 	/**
-	 * The Frobenius norm of this matrix.
-	 * @see [Matrix norm](https://en.wikipedia.org/wiki/Matrix_norm)
+	 * Determines whether this matrix is roughly equivalent to another.
+	 * @param matrix The other matrix.
+	 * @returns Whether the matrices are equivalent.
 	 */
-	get frob(): number;
+	equals(matrix: MatrixLike): boolean;
+
+	/**
+	 * Determines whether this matrix is exactly equivalent to another.
+	 * @param matrix The other matrix.
+	 * @returns Whether the matrices are equivalent.
+	 */
+	exactEquals(matrix: MatrixLike): boolean;
 
 	/**
 	 * Multiplies this matrix by another.

package/src/linear/SlowMatrix.ts

@@ -0,0 +1,228 @@
+import { type Matrix, PartialMatrixError, epsilon, MatrixSizeError } from "@lakuna/umath";
+
+/**
+ * A variable-size matrix.
+ * @see [Matrix](https://en.wikipedia.org/wiki/Matrix_(mathematics))
+ */
+export default class SlowMatrix extends Float32Array implements Matrix {
+	/**
+	 * Creates a variable-size matrix from the given columns.
+	 * @param cols The columns in the matrix.
+	 */
+	public constructor(...cols: Array<Array<number>>) {
+		super(cols.flat());
+		this.width = cols.length;
+		this.height = this.width && (cols[0] as Array<number>).length;
+
+		if (this.width < 1 || this.height < 1) {
+			throw new MatrixSizeError();
+		}
+
+		// Ensure that every column is the same height.
+		for (let i = 0; i < this.width; i++) {
+			if ((cols[i]?.length ?? 0) != this.height) {
+				throw new PartialMatrixError();
+			}
+		}
+	}
+
+	/** The number of rows in this matrix. */
+	public readonly height: number;
+
+	/** The number of columns in this matrix. */
+	public readonly width: number;
+
+	/**
+	 * The Frobenius norm of this matrix.
+	 * @see [Matrix norm](https://en.wikipedia.org/wiki/Matrix_norm)
+	 */
+	public get frob(): number {
+		return Math.hypot(...this);
+	}
+
+	/**
+	 * Adds two matrices of the same size.
+	 * @param matrix The other matrix.
+	 * @returns The sum of the matrices.
+	 * @see [Matrix addition](https://en.wikipedia.org/wiki/Matrix_addition)
+	 */
+	public add(matrix: Matrix): SlowMatrix {
+		if (this.width != matrix.width || this.height != matrix.height) {
+			throw new MatrixSizeError();
+		}
+
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[x] ??= [])[y] = (this[x * this.height + y] as number) + (matrix[x * matrix.height + y] as number);
+			}
+		}
+
+		return new SlowMatrix(...cols);
+	}
+
+	/**
+	 * Creates a copy of this matrix.
+	 * @returns A copy of this matrix.
+	 */
+	public clone(): SlowMatrix {
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[x] ??= [])[y] = this[x * this.height + y] as number;
+			}
+		}
+
+		return new SlowMatrix(...cols);
+	}
+
+	/**
+	 * Copies the values of another matrix into this one.
+	 * @param matrix The matrix to copy.
+	 * @returns This matrix.
+	 */
+	public copy(matrix: Matrix): this {
+		if (this.width != matrix.width || this.height != matrix.height) {
+			throw new MatrixSizeError();
+		}
+
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(this[x * this.width + y] as number) = matrix[x * this.height + y] as number;
+			}
+		}
+
+		return this;
+	}
+
+	/**
+	 * Determines whether this matrix is roughly equivalent to another.
+	 * @param matrix The other matrix.
+	 * @returns Whether the matrices are equivalent.
+	 */
+	public equals(matrix: Matrix): boolean {
+		if (this.width != matrix.width || this.height != matrix.height) {
+			return false;
+		}
+
+		for (let i = 0; i < this.length; i++) {
+			if (Math.abs((this[i] as number) - (matrix[i] as number)) > epsilon) {
+				return false;
+			}
+		}
+
+		return true;
+	}
+
+	/**
+	 * Determines whether this matrix is exactly equivalent to another.
+	 * @param matrix The other matrix.
+	 * @returns Whether the matrices are equivalent.
+	 */
+	public exactEquals(matrix: Matrix): boolean {
+		if (this.width != matrix.width || this.height != matrix.height) {
+			return false;
+		}
+
+		for (let i = 0; i < this.length; i++) {
+			if (this[i] != matrix[i]) {
+				return false;
+			}
+		}
+
+		return true;
+	}
+
+	/**
+	 * Multiplies this matrix by another.
+	 * @param matrix The other matrix.
+	 * @returns The product of the matrices.
+	 * @see [Matrix multiplication](https://en.wikipedia.org/wiki/Matrix_multiplication)
+	 */
+	public multiply(matrix: Matrix): SlowMatrix {
+		if (this.width != matrix.height) {
+			throw new MatrixSizeError();
+		}
+
+		const n: number = this.height;
+		const m: number = this.width;
+		const p: number = matrix.width;
+
+		const out: Array<Array<number>> = [];
+		for (let i = 0; i < n; i++) {
+			for (let j = 0; j < p; j++) {
+				let sum = 0;
+				for (let k = 0; k < m; k++) {
+					sum += (this[k * this.height + i] as number) * (matrix[j * matrix.height + k] as number)
+				}
+				(out[j] ??= [])[i] = sum;
+			}
+		}
+
+		return new SlowMatrix(...out);
+	}
+
+	/**
+	 * Multiplies this matrix by a scalar value.
+	 * @param scalar The scalar value.
+	 * @returns The product of the matrix and the scalar value.
+	 * @see [Matrix multiplication](https://en.wikipedia.org/wiki/Matrix_multiplication)
+	 */
+	public multiplyScalar(scalar: number): SlowMatrix {
+		const out = this.clone();
+		for (let i = 0; i < out.length; i++) {
+			out[i] *= scalar;
+		}
+
+		return out;
+	}
+
+	/**
+	 * Adds this matrix to another after multiplying the other by a scalar.
+	 * @param matrix The other matrix.
+	 * @param scalar The scalar.
+	 * @returns The sum.
+	 * @see [Matrix addition](https://en.wikipedia.org/wiki/Matrix_addition)
+	 * @see [Matrix multiplication](https://en.wikipedia.org/wiki/Matrix_multiplication)
+	 */
+	public multiplyScalarAndAdd(matrix: Matrix, scalar: number): SlowMatrix {
+		return this.add(matrix.multiplyScalar(scalar));
+	}
+
+	/**
+	 * Subtracts another matrix from this one.
+	 * @param matrix The other matrix.
+	 * @returns The difference between the matrices.
+	 * @see [Matrix addition](https://en.wikipedia.org/wiki/Matrix_addition)
+	 */
+	public subtract(matrix: Matrix): SlowMatrix {
+		if (this.width != matrix.width || this.height != matrix.height) {
+			throw new MatrixSizeError();
+		}
+
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[x] ??= [])[y] = (this[x * this.height + y] as number) - (matrix[x * matrix.height + y] as number);
+			}
+		}
+
+		return new SlowMatrix(...cols);
+	}
+
+	/**
+	 * Transposes this matrix.
+	 * @returns The transpose of this matrix.
+	 * @see [Transpose](https://en.wikipedia.org/wiki/Transpose)
+	 */
+	public transpose(): SlowMatrix {
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[y] ??= [])[x] = (this[x * this.height + y] as number);
+			}
+		}
+
+		return new SlowMatrix(...cols);
+	}
+}

package/src/linear/SlowSquareMatrix.ts

@@ -0,0 +1,203 @@
+import { type SquareMatrix, MatrixSizeError, SingularMatrixError } from "@lakuna/umath";
+import SlowMatrix from "@lakuna/umath/SlowMatrix";
+
+/**
+ * A variable-size matrix with the same number of rows and columns.
+ * @see [Square matrix](https://en.wikipedia.org/wiki/Square_matrix)
+ */
+export default class SlowSquareMatrix extends SlowMatrix implements SquareMatrix {
+	/**
+	 * Creates a variable-size matrix with the same number of rows and columns from the given columns.
+	 * @param cols The columns in the matrix.
+	 * @see [Square matrix](https://en.wikipedia.org/wiki/Square_matrix)
+	 */
+	public constructor(...cols: Array<Array<number>>) {
+		super(...cols);
+
+		// Ensure that every column is the same height and that the height is equal to the width.
+		for (let i = 0; i < this.width; i++) {
+			if ((cols[i] as Array<number>).length != this.height) {
+				throw new MatrixSizeError();
+			}
+		}
+	}
+
+	/**
+	 * The determinant of this matrix.
+	 * @see [Determinant](https://en.wikipedia.org/wiki/Determinant)
+	 */
+	public get determinant(): number {
+		if (this.length < 1) {
+			throw new MatrixSizeError();
+		}
+
+		if (this.length == 1) {
+			return this[0] ?? 0;
+		}
+
+		let out = 0;
+		for (let i = 0; i < this.width; i++) {
+			out += (this[i * this.height] as number) * (i % 2 ? -1 : 1) * this.minor(0, i);
+		}
+		
+		return out;
+	}
+
+	/**
+	 * Calculates the adjugate of this matrix.
+	 * @returns The adjugate of this matrix.
+	 * @see [Adjugate matrix](https://en.wikipedia.org/wiki/Adjugate_matrix)
+	 */
+	public adjoint(): SlowSquareMatrix {
+		return this.cofactor().transpose();
+	}
+
+	/**
+	 * Creates a copy of this matrix.
+	 * @returns A copy of this matrix.
+	 */
+	public override clone(): SlowSquareMatrix {
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[x] ??= [])[y] = this[x * this.height + y] as number;
+			}
+		}
+
+		return new SlowSquareMatrix(...cols);
+	}
+
+	/**
+	 * Calculates the cofactor of this matrix.
+	 * @returns The cofactor.
+	 * @see [Minor](https://en.wikipedia.org/wiki/Minor_(linear_algebra))
+	 */
+	public cofactor(): SlowSquareMatrix {
+		const out: Array<Array<number>> = [];
+		for (let i = 0; i < this.width; i++) { // col i
+			for (let j = 0; j < this.height; j++) { // row j
+				(out[j] ??= [])[i] = ((i + j) % 2 ? -1 : 1) * this.minor(i, j);
+			}
+		}
+
+		return new SlowSquareMatrix(...out);
+	}
+
+	/**
+	 * Resets this matrix to identity.
+	 * @returns This matrix.
+	 * @see [Identity matrix](https://en.wikipedia.org/wiki/Identity_matrix)
+	 */
+	public identity(): this {
+		for (let i = 0; i < this.width; i++) {
+			for (let j = 0; j < this.height; j++) {
+				this[i * this.height + j] = i == j ? 1 : 0;
+			}
+		}
+		
+		return this;
+	}
+
+	/**
+	 * Inverts this matrix.
+	 * @returns The inverted matrix.
+	 * @see [Invertible matrix](https://en.wikipedia.org/wiki/Invertible_matrix)
+	 */
+	public invert(): SlowSquareMatrix {
+		const dim: number = this.width;
+
+		const clone: SlowSquareMatrix = this.clone();
+		const identity: SlowSquareMatrix = this.clone().identity();
+
+		for (let i = 0; i < dim; i++) {
+			let diagonal: number = clone[i * dim + i] as number;
+			if (!diagonal) {
+				for (let ii: number = i + 1; ii < dim; ii++) {
+					if (clone[ii * dim + i]) {
+						for (let j = 0; j < dim; j++) {
+							for (const matrix of [clone, identity]) {
+								[matrix[i * dim + j], matrix[ii * dim + j]] = [matrix[ii * dim + j] as number, matrix[i * dim + j] as number];
+							}
+						}
+
+						break;
+					}
+				}
+
+				diagonal = clone[i * dim + i] as number;
+				if (!diagonal) {
+					throw new SingularMatrixError();
+				}
+			}
+
+			for (let ii = 0; ii < dim; ii++) {
+				if (ii == i) {
+					continue;
+				}
+
+				const temp: number = clone[ii * dim + i] as number;
+
+				for (let j = 0; j < dim; j++) {
+					for (const matrix of [clone, identity]) {
+						matrix[ii * dim + j] -= temp * (matrix[i * dim + j] as number);
+					}
+				}
+			}
+		}
+
+		return identity;
+	}
+
+	/**
+	 * Calculates the minor of this matrix which results from removing the given row and column.
+	 * @param row The row to remove.
+	 * @param col The column to remove.
+	 * @returns The minor.
+	 * @see [Minor](https://en.wikipedia.org/wiki/Minor_(linear_algebra))
+	 */
+	public minor(row: number, col: number): number {
+		return this.submatrix(row, col).determinant;
+	}
+
+	/**
+	 * Creates a submatrix by removing the given row and column from this matrix.
+	 * @param row The row to remove.
+	 * @param col The column to remove.
+	 * @returns The submatrix.
+	 */
+	public submatrix(row: number, col: number): SlowSquareMatrix {
+		const cols: Array<Array<number>> = [];
+
+		for (let i = 0; i < this.width; i++) {
+			if (i == col) {
+				continue;
+			}
+
+			for (let j = 0; j < this.height; j++) {
+				if (j == row) {
+					continue;
+				}
+
+				(cols[i < col ? i : i - 1] ??= [])[j < row ? j : j - 1] = this[i * this.height + j] as number;
+			}
+		}
+
+		return new SlowSquareMatrix(...cols);
+	}
+
+	/**
+	 * Transposes this matrix.
+	 * @returns The transpose of this matrix.
+	 * @see [Transpose](https://en.wikipedia.org/wiki/Transpose)
+	 */
+	public override transpose(): SlowSquareMatrix {
+		const cols: Array<Array<number>> = [];
+		for (let x = 0; x < this.width; x++) {
+			for (let y = 0; y < this.height; y++) {
+				(cols[y] ??= [])[x] = (this[x * this.height + y] as number);
+			}
+		}
+
+		return new SlowSquareMatrix(...cols);
+	}
+}

package/src/linear/SquareMatrix.ts

@@ -11,6 +11,13 @@ export default interface SquareMatrix extends Matrix {
 	 */
 	get determinant(): number;
 
+	/**
+	 * Calculates the adjugate of this matrix.
+	 * @returns The adjugate of this matrix.
+	 * @see [Adjugate matrix](https://en.wikipedia.org/wiki/Adjugate_matrix)
+	 */
+	adjoint(): Matrix;
+
 	/**
 	 * Resets this matrix to identity.
 	 * @returns This matrix.

package/src/utility/MagnitudeError.ts

@@ -6,5 +6,6 @@ export default class MagnitudeError extends Error {
 	 */
 	public constructor(message = "The vector is too small.") {
 		super(message);
+		this.name = "MagnitudeError";
 	}
 }

package/src/utility/MatrixSizeError.ts

@@ -0,0 +1,11 @@
+/** An error resulting from trying to use a matrix that is the wrong size. */
+export default class PartialMatrixError extends Error {
+	/**
+	 * Creates an error resulting from trying to use a matrix that is the wrong size.
+	 * @param message The message of the error.
+	 */
+	public constructor(message = "Invalid matrix dimensions.") {
+		super(message);
+		this.name = "MatrixSizeError";
+	}
+}

package/src/utility/PartialMatrixError.ts

@@ -0,0 +1,11 @@
+/** An error resulting from creating a matrix that isn't a perfect rectangle. */
+export default class PartialMatrixError extends Error {
+	/**
+	 * Creates an error resulting from creating a matrix that isn't a perfect rectangle.
+	 * @param message The message of the error.
+	 */
+	public constructor(message = "The matrix is not rectangular.") {
+		super(message);
+		this.name = "PartialMatrixError";
+	}
+}