@pawells/math-extended 1.0.1

package/build/matrices/core.d.ts

@@ -0,0 +1,168 @@
+import { IMatrix, IMatrix1, IMatrix2, IMatrix3, IMatrix4, TMatrixResult } from './types.js';
+export declare function MatrixCreate(): IMatrix1;
+export declare function MatrixCreate(size: 1): IMatrix1;
+export declare function MatrixCreate(size: 2): IMatrix2;
+export declare function MatrixCreate(size: 3): IMatrix3;
+export declare function MatrixCreate(size: 4): IMatrix4;
+export declare function MatrixCreate(size: number): IMatrix;
+export declare function MatrixCreate(rows: number, cols: number): IMatrix;
+/**
+ * Returns the dimensions of a matrix as [rows, columns].
+ * @param matrix - The matrix to measure
+ * @returns {[number, number]} Tuple [rows, columns]. Returns [0, 0] for empty matrices.
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixSize([[1, 2, 3], [4, 5, 6]]) // [2, 3]
+ */
+export declare function MatrixSize(matrix: IMatrix): [number, number];
+/**
+ * Returns the size of a square matrix (number of rows/columns).
+ * @param matrix - The square matrix to measure
+ * @returns {number} The size (n for an n×n matrix)
+ * @throws {Error} If the matrix is not square
+ * @example MatrixSizeSquare([[1, 2], [3, 4]]) // 2
+ */
+export declare function MatrixSizeSquare(matrix: IMatrix): number;
+/**
+ * Validates that the input is a well-formed matrix.
+ * @param matrix - The value to validate as a matrix
+ * @returns {boolean} True if the input is a valid matrix, false otherwise
+ * @example MatrixIsValid([[1, 2], [3, 4]]) // true
+ * @example MatrixIsValid([[1, 'invalid'], [3, 4]]) // false
+ */
+export declare function MatrixIsValid(matrix: unknown): boolean;
+/**
+ * Checks if a matrix is square (equal number of rows and columns).
+ * @param matrix - The matrix to check
+ * @returns {boolean} True if the matrix is square, false otherwise
+ * @example MatrixIsSquare([[1, 2], [3, 4]]) // true (2×2)
+ * @example MatrixIsSquare([[1, 2, 3], [4, 5, 6]]) // false (3×2)
+ */
+export declare function MatrixIsSquare(matrix: IMatrix): boolean;
+/**
+ * Checks if a matrix is a zero matrix (all elements are zero within tolerance).
+ * @param matrix - The matrix to check
+ * @param threshold - Tolerance for considering values as zero (default: 1e-14)
+ * @returns {boolean} True if all matrix elements are within threshold of zero
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsZero([[0, 0], [0, 0]]) // true
+ * @example MatrixIsZero([[1e-15, 0], [0, 0]]) // true (within default threshold)
+ * @example MatrixIsZero([[0.1, 0], [0, 0]]) // false
+ */
+export declare function MatrixIsZero(matrix: IMatrix, threshold?: number): boolean;
+/**
+ * Checks if a matrix is an identity matrix (1s on the main diagonal, 0s elsewhere).
+ * @param matrix - The matrix to check
+ * @param threshold - Tolerance for floating-point comparisons (default: 1e-14)
+ * @returns {boolean} True if the matrix is an identity matrix within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsIdentity([[1, 0], [0, 1]]) // true
+ * @example MatrixIsIdentity([[1, 0, 0], [0, 1, 0], [0, 0, 1]]) // true
+ * @example MatrixIsIdentity([[1, 1], [0, 1]]) // false
+ */
+export declare function MatrixIsIdentity(matrix: IMatrix, threshold?: number): boolean;
+/**
+ * Checks if a matrix is symmetric (A = Aᵀ, i.e. A[i][j] === A[j][i] for all i, j).
+ * @param matrix - The matrix to check (must be square)
+ * @param threshold - Tolerance for floating-point comparisons (default: 1e-14)
+ * @returns {boolean} True if the matrix is symmetric within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsSymmetric([[1, 2], [2, 1]]) // true
+ * @example MatrixIsSymmetric([[1, 2], [3, 4]]) // false
+ * @example MatrixIsSymmetric([[1, 2, 3], [2, 5, 4], [3, 4, 6]]) // true
+ */
+export declare function MatrixIsSymmetric(matrix: IMatrix, threshold?: number): boolean;
+/**
+ * Checks if a matrix is diagonal (all off-diagonal elements are zero within tolerance).
+ * A diagonal matrix has non-zero values only on its main diagonal.
+ * @param matrix - The matrix to check (must be square)
+ * @param threshold - Tolerance for considering values as zero (default: 1e-14)
+ * @returns {boolean} True if the matrix is diagonal within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsDiagonal([[3, 0], [0, 7]]) // true
+ * @example MatrixIsDiagonal([[1, 0, 0], [0, 5, 0], [0, 0, 2]]) // true
+ * @example MatrixIsDiagonal([[1, 2], [0, 1]]) // false
+ */
+export declare function MatrixIsDiagonal(matrix: IMatrix, threshold?: number): boolean;
+/**
+ * Creates an identity matrix of the specified size.
+ * @param size - The dimensions of the square identity matrix (must be non-negative integer)
+ * @returns {IMatrix} A square identity matrix of size n×n
+ * @throws {Error} If size is negative or not an integer
+ * @example MatrixIdentity(2) // [[1, 0], [0, 1]]
+ * @example MatrixIdentity(3) // [[1, 0, 0], [0, 1, 0], [0, 0, 1]]
+ */
+export declare function MatrixIdentity(size: 0): IMatrix;
+export declare function MatrixIdentity(size: 1): IMatrix1;
+export declare function MatrixIdentity(size: 2): IMatrix2;
+export declare function MatrixIdentity(size: 3): IMatrix3;
+export declare function MatrixIdentity(size: 4): IMatrix4;
+export declare function MatrixIdentity(size: number): IMatrix;
+/**
+ * Creates a deep copy of the given matrix.
+ * @param matrix - The matrix to clone (supports both IMatrix and IMatrixFlat formats)
+ * @returns {IMatrix | IMatrixFlat} A new matrix with identical values but independent memory allocation
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixClone([[1, 2], [3, 4]]) // [[1, 2], [3, 4]] (independent copy)
+ * @example MatrixClone({data: Float64Array[1,2,3,4], rows: 2, columns: 2}) // New IMatrixFlat with copied data
+ */
+export declare function MatrixClone<T extends IMatrix>(matrix: T): TMatrixResult<T>;
+/**
+ * Checks if two matrices are equal within a specified tolerance.
+ * @param a - First matrix to compare (supports both IMatrix and IMatrixFlat formats)
+ * @param b - Second matrix to compare (supports both IMatrix and IMatrixFlat formats)
+ * @param tolerance - Maximum allowed difference between corresponding elements (default: 1e-8)
+ * @returns {boolean} True if matrices are equal within tolerance, false otherwise
+ * @throws {Error} If either input is not a valid matrix or if matrix data is corrupted
+ * @example MatrixEquals([[1, 2]], [[1.0001, 2]], 0.001) // true
+ * @example MatrixEquals([[1, 2]], [[1, 3]]) // false
+ */
+export declare function MatrixEquals(a: IMatrix, b: IMatrix, tolerance?: number): boolean;
+/**
+ * Converts a matrix to a formatted string representation.
+ * @param matrix - The matrix to convert to string (supports both IMatrix and IMatrixFlat formats)
+ * @param precision - Number of decimal places for formatting (default: 2)
+ * @returns {string} A formatted string representation of the matrix
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixToString([[1.23, 2.7]]) // "[ 1.23, 2.70 ]"
+ * @example MatrixToString([[1, 2], [3, 4]], 0) // "[ 1, 2 ]\n[ 3, 4 ]"
+ */
+export declare function MatrixToString(matrix: IMatrix, precision?: number): string;
+/**
+ * Computes the rank of a matrix using Gaussian elimination.
+ * @param matrix - The input matrix (any dimensions, supports both IMatrix and IMatrixFlat formats)
+ * @param tolerance - Numerical tolerance for zero detection (default: 1e-10)
+ * @returns {number} The rank of the matrix (0 ≤ rank ≤ min(rows, columns))
+ * @throws {Error} If the matrix contains invalid values
+ * @example MatrixRank([[1, 2], [2, 4]]) // 1 (second row = 2 × first row)
+ * @example MatrixRank([[1, 0], [0, 1]]) // 2 (full rank)
+ */
+export declare function MatrixRank(matrix: IMatrix, tolerance?: number): number;
+/**
+ * Computes the trace of a matrix (sum of main diagonal elements).
+ * @param matrix - The input matrix (can be square or rectangular, supports both IMatrix and IMatrixFlat formats)
+ * @returns {number} The trace value (sum of diagonal elements)
+ * @throws {Error} If the matrix contains invalid values
+ * @example MatrixTrace([[1, 2, 3], [4, 5, 6], [7, 8, 9]]) // 15 (1 + 5 + 9)
+ * @example MatrixTrace([[1, 2], [3, 4], [5, 6]]) // 5 (1 + 4, rectangular matrix)
+ */
+export declare function MatrixTrace(matrix: IMatrix): number;
+/**
+ * Returns the transpose of a matrix (rows become columns and vice versa).
+ * @param matrix - The matrix to transpose (can be any m×n matrix, supports both IMatrix and IMatrixFlat formats)
+ * @returns {IMatrix | IMatrixFlat} The transposed matrix with dimensions n×m (same format as input)
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixTranspose([[1, 2, 3], [4, 5, 6]]) // [[1, 4], [2, 5], [3, 6]]
+ * @example MatrixTranspose([[1, 2], [3, 4]]) // [[1, 3], [2, 4]]
+ */
+export declare function MatrixTranspose<T extends IMatrix>(matrix: T): TMatrixResult<T>;
+/**
+ * Applies a transformation function to each element of the matrix.
+ * @param matrix - The input matrix to transform (supports both IMatrix and IMatrixFlat formats)
+ * @param fn - Transformation function: (value, row, col) => transformedValue
+ * @returns {IMatrix | IMatrixFlat} A new matrix with transformed values (same dimensions and format as input)
+ * @throws {Error} If the input matrix is invalid or transformation function throws
+ * @example MatrixMap([[1, 2], [3, 4]], (value) => value * value) // [[1, 4], [9, 16]]
+ * @example MatrixMap([[1, 2], [3, 4]], (value, row, col) => value + row + col) // [[1, 3], [4, 6]]
+ */
+export declare function MatrixMap<T extends IMatrix>(matrix: T, fn: (value: number, row: number, col: number) => number): TMatrixResult<T>;
+//# sourceMappingURL=core.d.ts.map

package/build/matrices/core.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"core.d.ts","sourceRoot":"","sources":["../../src/matrices/core.ts"],"names":[],"mappings":"AAEA,OAAO,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,QAAQ,EAAE,aAAa,EAAE,MAAM,YAAY,CAAC;AAI5F,wBAAgB,YAAY,IAAI,QAAQ,CAAC;AAEzC,wBAAgB,YAAY,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAEhD,wBAAgB,YAAY,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAEhD,wBAAgB,YAAY,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAEhD,wBAAgB,YAAY,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAEhD,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,OAAO,CAAC;AAEpD,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,GAAG,OAAO,CAAC;AA4DlE;;;;;;GAMG;AACH,wBAAgB,UAAU,CAAC,MAAM,EAAE,OAAO,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAS5D;AAED;;;;;;GAMG;AACH,wBAAgB,gBAAgB,CAAC,MAAM,EAAE,OAAO,GAAG,MAAM,CAKxD;AAED;;;;;;GAMG;AACH,wBAAgB,aAAa,CAAC,MAAM,EAAE,OAAO,GAAG,OAAO,CAOtD;AAED;;;;;;GAMG;AACH,wBAAgB,cAAc,CAAC,MAAM,EAAE,OAAO,GAAG,OAAO,CAOvD;AAED;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAc,GAAG,OAAO,CAoBhF;AAED;;;;;;;;;GASG;AACH,wBAAgB,gBAAgB,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAc,GAAG,OAAO,CAoBpF;AAED;;;;;;;;;GASG;AACH,wBAAgB,iBAAiB,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAc,GAAG,OAAO,CAwBrF;AAED;;;;;;;;;;GAUG;AACH,wBAAgB,gBAAgB,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAc,GAAG,OAAO,CAqBpF;AAED;;;;;;;GAOG;AAEH,wBAAgB,cAAc,CAAC,IAAI,EAAE,CAAC,GAAG,OAAO,CAAC;AAEjD,wBAAgB,cAAc,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAElD,wBAAgB,cAAc,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAElD,wBAAgB,cAAc,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAElD,wBAAgB,cAAc,CAAC,IAAI,EAAE,CAAC,GAAG,QAAQ,CAAC;AAElD,wBAAgB,cAAc,CAAC,IAAI,EAAE,MAAM,GAAG,OAAO,CAAC;AAgBtD;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,OAAO,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAI1E;AAED;;;;;;;;;GASG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,EAAE,SAAS,GAAE,MAAa,GAAG,OAAO,CA+BtF;AAED;;;;;;;;GAQG;AACH,wBAAgB,cAAc,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAU,GAAG,MAAM,CAW7E;AAED;;;;;;;;GAQG;AACH,wBAAgB,UAAU,CAAC,MAAM,EAAE,OAAO,EAAE,SAAS,GAAE,MAAc,GAAG,MAAM,CAiE7E;AAED;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,MAAM,EAAE,OAAO,GAAG,MAAM,CAmBnD;AAED;;;;;;;GAOG;AACH,wBAAgB,eAAe,CAAC,CAAC,SAAS,OAAO,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CA4B9E;AAED;;;;;;;;GAQG;AACH,wBAAgB,SAAS,CAAC,CAAC,SAAS,OAAO,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,KAAK,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CAqBjI"}

package/build/matrices/core.js

@@ -0,0 +1,457 @@
+import { AssertNumber } from '@pawells/typescript-common';
+import { AssertMatrix, AssertMatrixRow, AssertMatrixValue } from './asserts.js';
+/**
+ * Creates a matrix with specified dimensions, initialized with zeros.
+ * @param rows - Number of rows (non-negative integer). If only this parameter is provided, creates a square matrix.
+ * @param cols - Number of columns (non-negative integer). Optional for square matrices.
+ * @returns {IMatrix | IMatrix1 | IMatrix2 | IMatrix3 | IMatrix4} Zero-filled matrix with specified dimensions
+ * @throws {Error} If rows or cols are negative or not integers
+ * @example MatrixCreate() // [[0]] (1x1 matrix)
+ * @example MatrixCreate(2) // [[0, 0], [0, 0]] (2x2 matrix)
+ * @example MatrixCreate(2, 3) // [[0, 0, 0], [0, 0, 0]] (2x3 matrix)
+ */
+export function MatrixCreate(rows, cols) {
+    // Handle no parameters - default to 1x1
+    if (rows === undefined) {
+        return [[0]];
+    }
+    // Handle single parameter - create square matrix
+    if (cols === undefined) {
+        // Return specific types for common square matrices
+        switch (rows) {
+            case 1:
+                return [[0]];
+            case 2:
+                return [
+                    [0, 0],
+                    [0, 0],
+                ];
+            case 3:
+                return [
+                    [0, 0, 0],
+                    [0, 0, 0],
+                    [0, 0, 0],
+                ];
+            case 4:
+                return [
+                    [0, 0, 0, 0],
+                    [0, 0, 0, 0],
+                    [0, 0, 0, 0],
+                    [0, 0, 0, 0],
+                ];
+        }
+    }
+    // Validate parameters
+    const effectiveCols = cols ?? rows;
+    AssertNumber(rows, { integer: true, gte: 0 }, { message: 'Rows must be a non-negative integer' });
+    AssertNumber(effectiveCols, { integer: true, gte: 0 }, { message: 'Columns must be a non-negative integer' });
+    const result = [];
+    // Create each row filled with zeros
+    for (let i = 0; i < rows; i++) {
+        result.push(new Array(effectiveCols).fill(0));
+    }
+    return result;
+}
+/**
+ * Returns the dimensions of a matrix as [rows, columns].
+ * @param matrix - The matrix to measure
+ * @returns {[number, number]} Tuple [rows, columns]. Returns [0, 0] for empty matrices.
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixSize([[1, 2, 3], [4, 5, 6]]) // [2, 3]
+ */
+export function MatrixSize(matrix) {
+    // Basic validation without calling AssertMatrix to avoid circular dependency
+    if (!Array.isArray(matrix)) {
+        throw new Error('Input must be an array');
+    }
+    // Handle empty matrix or matrix with no columns
+    if (matrix.length === 0 || !Array.isArray(matrix[0]))
+        return [0, 0];
+    return [matrix.length, matrix[0]?.length ?? 0];
+}
+/**
+ * Returns the size of a square matrix (number of rows/columns).
+ * @param matrix - The square matrix to measure
+ * @returns {number} The size (n for an n×n matrix)
+ * @throws {Error} If the matrix is not square
+ * @example MatrixSizeSquare([[1, 2], [3, 4]]) // 2
+ */
+export function MatrixSizeSquare(matrix) {
+    AssertMatrix(matrix, { square: true });
+    const [rows] = MatrixSize(matrix);
+    return rows;
+}
+/**
+ * Validates that the input is a well-formed matrix.
+ * @param matrix - The value to validate as a matrix
+ * @returns {boolean} True if the input is a valid matrix, false otherwise
+ * @example MatrixIsValid([[1, 2], [3, 4]]) // true
+ * @example MatrixIsValid([[1, 'invalid'], [3, 4]]) // false
+ */
+export function MatrixIsValid(matrix) {
+    try {
+        AssertMatrix(matrix);
+        return true;
+    }
+    catch {
+        return false;
+    }
+}
+/**
+ * Checks if a matrix is square (equal number of rows and columns).
+ * @param matrix - The matrix to check
+ * @returns {boolean} True if the matrix is square, false otherwise
+ * @example MatrixIsSquare([[1, 2], [3, 4]]) // true (2×2)
+ * @example MatrixIsSquare([[1, 2, 3], [4, 5, 6]]) // false (3×2)
+ */
+export function MatrixIsSquare(matrix) {
+    try {
+        MatrixSizeSquare(matrix);
+        return true;
+    }
+    catch {
+        return false;
+    }
+}
+/**
+ * Checks if a matrix is a zero matrix (all elements are zero within tolerance).
+ * @param matrix - The matrix to check
+ * @param threshold - Tolerance for considering values as zero (default: 1e-14)
+ * @returns {boolean} True if all matrix elements are within threshold of zero
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsZero([[0, 0], [0, 0]]) // true
+ * @example MatrixIsZero([[1e-15, 0], [0, 0]]) // true (within default threshold)
+ * @example MatrixIsZero([[0.1, 0], [0, 0]]) // false
+ */
+export function MatrixIsZero(matrix, threshold = 1e-14) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    if (rows === 0 || cols === 0)
+        return true; // Empty matrix is considered zero
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        for (let col = 0; col < cols; col++) {
+            const value = matrixRow[col];
+            AssertMatrixValue(value, { rowIndex: row, columnIndex: col });
+            if (Math.abs(value) > threshold) {
+                return false; // Found a non-zero value
+            }
+        }
+    }
+    return true; // All values are within the threshold
+}
+/**
+ * Checks if a matrix is an identity matrix (1s on the main diagonal, 0s elsewhere).
+ * @param matrix - The matrix to check
+ * @param threshold - Tolerance for floating-point comparisons (default: 1e-14)
+ * @returns {boolean} True if the matrix is an identity matrix within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsIdentity([[1, 0], [0, 1]]) // true
+ * @example MatrixIsIdentity([[1, 0, 0], [0, 1, 0], [0, 0, 1]]) // true
+ * @example MatrixIsIdentity([[1, 1], [0, 1]]) // false
+ */
+export function MatrixIsIdentity(matrix, threshold = 1e-14) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    if (rows !== cols)
+        return false; // Identity must be square
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        for (let col = 0; col < cols; col++) {
+            const value = matrixRow[col];
+            AssertMatrixValue(value, { rowIndex: row, columnIndex: col });
+            const expected = row === col ? 1 : 0;
+            if (Math.abs(value - expected) > threshold)
+                return false;
+        }
+    }
+    return true;
+}
+/**
+ * Checks if a matrix is symmetric (A = Aᵀ, i.e. A[i][j] === A[j][i] for all i, j).
+ * @param matrix - The matrix to check (must be square)
+ * @param threshold - Tolerance for floating-point comparisons (default: 1e-14)
+ * @returns {boolean} True if the matrix is symmetric within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsSymmetric([[1, 2], [2, 1]]) // true
+ * @example MatrixIsSymmetric([[1, 2], [3, 4]]) // false
+ * @example MatrixIsSymmetric([[1, 2, 3], [2, 5, 4], [3, 4, 6]]) // true
+ */
+export function MatrixIsSymmetric(matrix, threshold = 1e-14) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    if (rows !== cols)
+        return false; // Symmetry requires a square matrix
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        for (let col = row + 1; col < cols; col++) {
+            const upper = matrixRow[col];
+            AssertMatrixValue(upper, { rowIndex: row, columnIndex: col });
+            const lowerRow = matrix[col];
+            AssertMatrixRow(lowerRow, { rowIndex: col });
+            const lower = lowerRow[row];
+            AssertMatrixValue(lower, { rowIndex: col, columnIndex: row });
+            if (Math.abs(upper - lower) > threshold)
+                return false;
+        }
+    }
+    return true;
+}
+/**
+ * Checks if a matrix is diagonal (all off-diagonal elements are zero within tolerance).
+ * A diagonal matrix has non-zero values only on its main diagonal.
+ * @param matrix - The matrix to check (must be square)
+ * @param threshold - Tolerance for considering values as zero (default: 1e-14)
+ * @returns {boolean} True if the matrix is diagonal within tolerance
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixIsDiagonal([[3, 0], [0, 7]]) // true
+ * @example MatrixIsDiagonal([[1, 0, 0], [0, 5, 0], [0, 0, 2]]) // true
+ * @example MatrixIsDiagonal([[1, 2], [0, 1]]) // false
+ */
+export function MatrixIsDiagonal(matrix, threshold = 1e-14) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    if (rows !== cols)
+        return false; // Diagonality requires a square matrix
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        for (let col = 0; col < cols; col++) {
+            if (row === col)
+                continue; // Skip diagonal elements
+            const value = matrixRow[col];
+            AssertMatrixValue(value, { rowIndex: row, columnIndex: col });
+            if (Math.abs(value) > threshold)
+                return false;
+        }
+    }
+    return true;
+}
+export function MatrixIdentity(size) {
+    AssertNumber(size, { integer: true, gte: 0 }, { message: 'Size must be a non-negative integer' });
+    const result = MatrixCreate(size, size);
+    // Set diagonal elements to 1
+    for (let i = 0; i < size; i++) {
+        const row = result[i];
+        AssertMatrixRow(row, { rowIndex: i });
+        row[i] = 1;
+    }
+    return result;
+}
+/**
+ * Creates a deep copy of the given matrix.
+ * @param matrix - The matrix to clone (supports both IMatrix and IMatrixFlat formats)
+ * @returns {IMatrix | IMatrixFlat} A new matrix with identical values but independent memory allocation
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixClone([[1, 2], [3, 4]]) // [[1, 2], [3, 4]] (independent copy)
+ * @example MatrixClone({data: Float64Array[1,2,3,4], rows: 2, columns: 2}) // New IMatrixFlat with copied data
+ */
+export function MatrixClone(matrix) {
+    AssertMatrix(matrix);
+    // Create a deep copy by mapping each row to a new array
+    return matrix.map((row) => [...row]);
+}
+/**
+ * Checks if two matrices are equal within a specified tolerance.
+ * @param a - First matrix to compare (supports both IMatrix and IMatrixFlat formats)
+ * @param b - Second matrix to compare (supports both IMatrix and IMatrixFlat formats)
+ * @param tolerance - Maximum allowed difference between corresponding elements (default: 1e-8)
+ * @returns {boolean} True if matrices are equal within tolerance, false otherwise
+ * @throws {Error} If either input is not a valid matrix or if matrix data is corrupted
+ * @example MatrixEquals([[1, 2]], [[1.0001, 2]], 0.001) // true
+ * @example MatrixEquals([[1, 2]], [[1, 3]]) // false
+ */
+export function MatrixEquals(a, b, tolerance = 1e-8) {
+    AssertMatrix(a);
+    AssertMatrix(b);
+    AssertNumber(tolerance, { gte: 0 }, { message: 'Tolerance must be a non-negative number' });
+    const [rowsA, colsA] = MatrixSize(a);
+    const [rowsB, colsB] = MatrixSize(b);
+    if (rowsA !== rowsB || colsA !== colsB) {
+        return false;
+    }
+    for (let row = 0; row < rowsA; row++) {
+        const rowA = a[row];
+        const rowB = b[row];
+        AssertMatrixRow(rowA, { rowIndex: row });
+        AssertMatrixRow(rowB, { rowIndex: row });
+        for (let col = 0; col < colsA; col++) {
+            const valA = rowA[col];
+            const valB = rowB[col];
+            AssertMatrixValue(valA, { rowIndex: row, columnIndex: col });
+            AssertMatrixValue(valB, { rowIndex: row, columnIndex: col });
+            if (Math.abs(valA - valB) > tolerance) {
+                return false;
+            }
+        }
+    }
+    return true;
+}
+/**
+ * Converts a matrix to a formatted string representation.
+ * @param matrix - The matrix to convert to string (supports both IMatrix and IMatrixFlat formats)
+ * @param precision - Number of decimal places for formatting (default: 2)
+ * @returns {string} A formatted string representation of the matrix
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixToString([[1.23, 2.7]]) // "[ 1.23, 2.70 ]"
+ * @example MatrixToString([[1, 2], [3, 4]], 0) // "[ 1, 2 ]\n[ 3, 4 ]"
+ */
+export function MatrixToString(matrix, precision = 2) {
+    AssertMatrix(matrix);
+    AssertNumber(precision, { integer: true, gte: 0 }, { message: 'Precision must be a non-negative integer' });
+    return matrix
+        .map((row) => '[ ' + row
+        .map((val) => {
+        return typeof val === 'number' ? val.toFixed(precision) : String(val);
+    })
+        .join(', ') + ' ]')
+        .join('\n');
+}
+/**
+ * Computes the rank of a matrix using Gaussian elimination.
+ * @param matrix - The input matrix (any dimensions, supports both IMatrix and IMatrixFlat formats)
+ * @param tolerance - Numerical tolerance for zero detection (default: 1e-10)
+ * @returns {number} The rank of the matrix (0 ≤ rank ≤ min(rows, columns))
+ * @throws {Error} If the matrix contains invalid values
+ * @example MatrixRank([[1, 2], [2, 4]]) // 1 (second row = 2 × first row)
+ * @example MatrixRank([[1, 0], [0, 1]]) // 2 (full rank)
+ */
+export function MatrixRank(matrix, tolerance = 1e-10) {
+    AssertMatrix(matrix);
+    AssertNumber(tolerance, { gte: 0 }, { message: 'Tolerance must be a non-negative number' });
+    const [rows, cols] = MatrixSize(matrix);
+    if (rows === 0 || cols === 0)
+        return 0;
+    // Make a deep copy to avoid mutating the input
+    const mat = matrix.map((row) => [...row]);
+    let rank = 0;
+    const rowUsed = new Array(rows).fill(false);
+    for (let col = 0; col < cols; col++) {
+        let pivotRow = -1;
+        for (let row = 0; row < rows; row++) {
+            const matrixRow = mat[row];
+            AssertMatrixRow(matrixRow, { rowIndex: row });
+            const value = matrixRow[col];
+            AssertMatrixValue(value, { rowIndex: row, columnIndex: col });
+            if (!rowUsed[row] && Math.abs(value) > tolerance) {
+                pivotRow = row;
+                break;
+            }
+        }
+        if (pivotRow === -1)
+            continue;
+        rowUsed[pivotRow] = true;
+        rank++;
+        // Eliminate this column in all other rows
+        for (let row = 0; row < rows; row++) {
+            if (row !== pivotRow) {
+                const currentRow = mat[row];
+                const pivotRowData = mat[pivotRow];
+                AssertMatrixRow(currentRow, { rowIndex: row });
+                AssertMatrixRow(pivotRowData, { rowIndex: pivotRow });
+                const pivotValue = pivotRowData[col];
+                const currentValue = currentRow[col];
+                AssertMatrixValue(pivotValue, { rowIndex: pivotRow, columnIndex: col });
+                AssertMatrixValue(currentValue, { rowIndex: row, columnIndex: col });
+                const factor = currentValue / pivotValue;
+                for (let k = col; k < cols; k++) {
+                    const pivotK = pivotRowData[k];
+                    const currentK = currentRow[k];
+                    AssertMatrixValue(pivotK, { rowIndex: pivotRow, columnIndex: k });
+                    AssertMatrixValue(currentK, { rowIndex: row, columnIndex: k });
+                    currentRow[k] = currentK - (factor * pivotK);
+                    const updatedValue = currentRow[k];
+                    AssertMatrixValue(updatedValue, { rowIndex: row, columnIndex: k });
+                    if (Math.abs(updatedValue) < tolerance) {
+                        currentRow[k] = 0;
+                    }
+                }
+            }
+        }
+    }
+    return rank;
+}
+/**
+ * Computes the trace of a matrix (sum of main diagonal elements).
+ * @param matrix - The input matrix (can be square or rectangular, supports both IMatrix and IMatrixFlat formats)
+ * @returns {number} The trace value (sum of diagonal elements)
+ * @throws {Error} If the matrix contains invalid values
+ * @example MatrixTrace([[1, 2, 3], [4, 5, 6], [7, 8, 9]]) // 15 (1 + 5 + 9)
+ * @example MatrixTrace([[1, 2], [3, 4], [5, 6]]) // 5 (1 + 4, rectangular matrix)
+ */
+export function MatrixTrace(matrix) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    let trace = 0;
+    // The trace is the sum of the diagonal elements (i.e., elements where row === col)
+    const n = Math.min(rows, cols);
+    for (let i = 0; i < n; i++) {
+        const row = matrix[i];
+        AssertMatrixRow(row, { rowIndex: i });
+        const val = row[i];
+        AssertMatrixValue(val, { rowIndex: i, columnIndex: i });
+        trace += val;
+    }
+    return trace;
+}
+/**
+ * Returns the transpose of a matrix (rows become columns and vice versa).
+ * @param matrix - The matrix to transpose (can be any m×n matrix, supports both IMatrix and IMatrixFlat formats)
+ * @returns {IMatrix | IMatrixFlat} The transposed matrix with dimensions n×m (same format as input)
+ * @throws {Error} If the input is not a valid matrix
+ * @example MatrixTranspose([[1, 2, 3], [4, 5, 6]]) // [[1, 4], [2, 5], [3, 6]]
+ * @example MatrixTranspose([[1, 2], [3, 4]]) // [[1, 3], [2, 4]]
+ */
+export function MatrixTranspose(matrix) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    // Note: result dimensions are swapped (columns become rows)
+    const result = [];
+    // Initialize result matrix with swapped dimensions
+    for (let i = 0; i < cols; i++) {
+        result.push(new Array(rows).fill(0));
+    }
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        for (let col = 0; col < cols; col++) {
+            const val = matrixRow[col];
+            AssertMatrixValue(val, { rowIndex: row, columnIndex: col });
+            // Swap row and column indices in the result
+            const resultCol = result[col];
+            AssertMatrixRow(resultCol, { rowIndex: col });
+            resultCol[row] = val;
+        }
+    }
+    return result;
+}
+/**
+ * Applies a transformation function to each element of the matrix.
+ * @param matrix - The input matrix to transform (supports both IMatrix and IMatrixFlat formats)
+ * @param fn - Transformation function: (value, row, col) => transformedValue
+ * @returns {IMatrix | IMatrixFlat} A new matrix with transformed values (same dimensions and format as input)
+ * @throws {Error} If the input matrix is invalid or transformation function throws
+ * @example MatrixMap([[1, 2], [3, 4]], (value) => value * value) // [[1, 4], [9, 16]]
+ * @example MatrixMap([[1, 2], [3, 4]], (value, row, col) => value + row + col) // [[1, 3], [4, 6]]
+ */
+export function MatrixMap(matrix, fn) {
+    AssertMatrix(matrix);
+    const [rows, cols] = MatrixSize(matrix);
+    const result = [];
+    for (let row = 0; row < rows; row++) {
+        const matrixRow = matrix[row];
+        AssertMatrixRow(matrixRow, { rowIndex: row });
+        const newRow = [];
+        for (let col = 0; col < cols; col++) {
+            const val = matrixRow[col];
+            AssertMatrixValue(val, { rowIndex: row, columnIndex: col });
+            newRow.push(fn(val, row, col));
+        }
+        result.push(newRow);
+    }
+    return result;
+}
+//# sourceMappingURL=core.js.map

package/build/matrices/core.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"core.js","sourceRoot":"","sources":["../../src/matrices/core.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAC1D,OAAO,EAAE,YAAY,EAAE,eAAe,EAAE,iBAAiB,EAAE,MAAM,cAAc,CAAC;AAmBhF;;;;;;;;;GASG;AACH,MAAM,UAAU,YAAY,CAAC,IAAa,EAAE,IAAa;IACxD,wCAAwC;IACxC,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;QACxB,OAAO,CAAC,CAAC,CAAC,CAAC,CAAa,CAAC;IAC1B,CAAC;IAED,iDAAiD;IACjD,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;QACxB,mDAAmD;QACnD,QAAQ,IAAI,EAAE,CAAC;YACd,KAAK,CAAC;gBACL,OAAO,CAAC,CAAC,CAAC,CAAC,CAAa,CAAC;YAC1B,KAAK,CAAC;gBACL,OAAO;oBACN,CAAC,CAAC,EAAE,CAAC,CAAC;oBACN,CAAC,CAAC,EAAE,CAAC,CAAC;iBACM,CAAC;YACf,KAAK,CAAC;gBACL,OAAO;oBACN,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBACT,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;iBACG,CAAC;YACf,KAAK,CAAC;gBACL,OAAO;oBACN,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;oBACZ,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;iBACA,CAAC;QAChB,CAAC;IACF,CAAC;IAED,sBAAsB;IACtB,MAAM,aAAa,GAAG,IAAI,IAAI,IAAI,CAAC;IACnC,YAAY,CAAC,IAAI,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,qCAAqC,EAAE,CAAC,CAAC;IAClG,YAAY,CAAC,aAAa,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,wCAAwC,EAAE,CAAC,CAAC;IAE9G,MAAM,MAAM,GAAe,EAAE,CAAC;IAE9B,oCAAoC;IACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,aAAa,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAC/C,CAAC;IAED,OAAO,MAAM,CAAC;AACf,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,UAAU,CAAC,MAAe;IACzC,6EAA6E;IAC7E,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,EAAE,CAAC;QAC5B,MAAM,IAAI,KAAK,CAAC,wBAAwB,CAAC,CAAC;IAC3C,CAAC;IAED,gDAAgD;IAChD,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACpE,OAAO,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,IAAI,CAAC,CAAC,CAAC;AAChD,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAe;IAC/C,YAAY,CAAC,MAAM,EAAE,EAAE,MAAM,EAAE,IAAI,EAAE,CAAC,CAAC;IAEvC,MAAM,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IAClC,OAAO,IAAI,CAAC;AACb,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,aAAa,CAAC,MAAe;IAC5C,IAAI,CAAC;QACJ,YAAY,CAAC,MAAM,CAAC,CAAC;QACrB,OAAO,IAAI,CAAC;IACb,CAAC;IAAC,MAAM,CAAC;QACR,OAAO,KAAK,CAAC;IACd,CAAC;AACF,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAAC,MAAe;IAC7C,IAAI,CAAC;QACJ,gBAAgB,CAAC,MAAM,CAAC,CAAC;QACzB,OAAO,IAAI,CAAC;IACb,CAAC;IAAC,MAAM,CAAC;QACR,OAAO,KAAK,CAAC;IACd,CAAC;AACF,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,YAAY,CAAC,MAAe,EAAE,YAAoB,KAAK;IACtE,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,IAAI,CAAC,CAAC,kCAAkC;IAC7E,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC7B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9D,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,SAAS,EAAE,CAAC;gBACjC,OAAO,KAAK,CAAC,CAAC,yBAAyB;YACxC,CAAC;QACF,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC,CAAC,sCAAsC;AACpD,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAe,EAAE,YAAoB,KAAK;IAC1E,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,IAAI,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC,CAAC,0BAA0B;IAE3D,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC7B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9D,MAAM,QAAQ,GAAG,GAAG,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,QAAQ,CAAC,GAAG,SAAS;gBAAE,OAAO,KAAK,CAAC;QAC1D,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC;AACb,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,iBAAiB,CAAC,MAAe,EAAE,YAAoB,KAAK;IAC3E,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,IAAI,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC,CAAC,oCAAoC;IAErE,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,KAAK,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YAC3C,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC7B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9D,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;YAC7B,eAAe,CAAC,QAAQ,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;YAC7C,MAAM,KAAK,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC5B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9D,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,SAAS;gBAAE,OAAO,KAAK,CAAC;QACvD,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC;AACb,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,gBAAgB,CAAC,MAAe,EAAE,YAAoB,KAAK;IAC1E,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,IAAI,KAAK,IAAI;QAAE,OAAO,KAAK,CAAC,CAAC,uCAAuC;IAExE,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,IAAI,GAAG,KAAK,GAAG;gBAAE,SAAS,CAAC,yBAAyB;YAEpD,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC7B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9D,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,SAAS;gBAAE,OAAO,KAAK,CAAC;QAC/C,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC;AACb,CAAC;AAsBD,MAAM,UAAU,cAAc,CAAC,IAAY;IAC1C,YAAY,CAAC,IAAI,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,qCAAqC,EAAE,CAAC,CAAC;IAElG,MAAM,MAAM,GAAG,YAAY,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IAExC,6BAA6B;IAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACtB,eAAe,CAAC,GAAG,EAAE,EAAE,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC;QACtC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IACZ,CAAC;IAED,OAAO,MAAM,CAAC;AACf,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,WAAW,CAAoB,MAAS;IACvD,YAAY,CAAC,MAAM,CAAC,CAAC;IACrB,wDAAwD;IACxD,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAqB,CAAC;AAC1D,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,YAAY,CAAC,CAAU,EAAE,CAAU,EAAE,YAAoB,IAAI;IAC5E,YAAY,CAAC,CAAC,CAAC,CAAC;IAChB,YAAY,CAAC,CAAC,CAAC,CAAC;IAChB,YAAY,CAAC,SAAS,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,yCAAyC,EAAE,CAAC,CAAC;IAE5F,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IACrC,MAAM,CAAC,KAAK,EAAE,KAAK,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAErC,IAAI,KAAK,KAAK,KAAK,IAAI,KAAK,KAAK,KAAK,EAAE,CAAC;QACxC,OAAO,KAAK,CAAC;IACd,CAAC;IAED,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC;QACtC,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACpB,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QACpB,eAAe,CAAC,IAAI,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QACzC,eAAe,CAAC,IAAI,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAEzC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,KAAK,EAAE,GAAG,EAAE,EAAE,CAAC;YACtC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YACvB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YACvB,iBAAiB,CAAC,IAAI,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAC7D,iBAAiB,CAAC,IAAI,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE7D,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,SAAS,EAAE,CAAC;gBACvC,OAAO,KAAK,CAAC;YACd,CAAC;QACF,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC;AACb,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,cAAc,CAAC,MAAe,EAAE,YAAoB,CAAC;IACpE,YAAY,CAAC,MAAM,CAAC,CAAC;IACrB,YAAY,CAAC,SAAS,EAAE,EAAE,OAAO,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,0CAA0C,EAAE,CAAC,CAAC;IAE5G,OAAO,MAAM;SACX,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,GAAG,GAAG;SACtB,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE;QACZ,OAAO,OAAO,GAAG,KAAK,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACvE,CAAC,CAAC;SACD,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC;SACnB,IAAI,CAAC,IAAI,CAAC,CAAC;AACd,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,UAAU,CAAC,MAAe,EAAE,YAAoB,KAAK;IACpE,YAAY,CAAC,MAAM,CAAC,CAAC;IACrB,YAAY,CAAC,SAAS,EAAE,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,EAAE,OAAO,EAAE,yCAAyC,EAAE,CAAC,CAAC;IAE5F,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,IAAI,KAAK,CAAC,IAAI,IAAI,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAEvC,+CAA+C;IAC/C,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC1C,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,MAAM,OAAO,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;IAE5C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,IAAI,QAAQ,GAAG,CAAC,CAAC,CAAC;QAElB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,MAAM,SAAS,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;YAC3B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;YAE9C,MAAM,KAAK,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC7B,iBAAiB,CAAC,KAAK,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAC9D,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,SAAS,EAAE,CAAC;gBAClD,QAAQ,GAAG,GAAG,CAAC;gBACf,MAAM;YACP,CAAC;QACF,CAAC;QACD,IAAI,QAAQ,KAAK,CAAC,CAAC;YAAE,SAAS;QAE9B,OAAO,CAAC,QAAQ,CAAC,GAAG,IAAI,CAAC;QACzB,IAAI,EAAE,CAAC;QAEP,0CAA0C;QAC1C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,IAAI,GAAG,KAAK,QAAQ,EAAE,CAAC;gBACtB,MAAM,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBAC5B,MAAM,YAAY,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC;gBACnC,eAAe,CAAC,UAAU,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;gBAC/C,eAAe,CAAC,YAAY,EAAE,EAAE,QAAQ,EAAE,QAAQ,EAAE,CAAC,CAAC;gBAEtD,MAAM,UAAU,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC;gBACrC,MAAM,YAAY,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;gBACrC,iBAAiB,CAAC,UAAU,EAAE,EAAE,QAAQ,EAAE,QAAQ,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;gBACxE,iBAAiB,CAAC,YAAY,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;gBAErE,MAAM,MAAM,GAAG,YAAY,GAAG,UAAU,CAAC;gBAEzC,KAAK,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;oBACjC,MAAM,MAAM,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC;oBAC/B,MAAM,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;oBAC/B,iBAAiB,CAAC,MAAM,EAAE,EAAE,QAAQ,EAAE,QAAQ,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC;oBAClE,iBAAiB,CAAC,QAAQ,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC;oBAE/D,UAAU,CAAC,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;oBAE7C,MAAM,YAAY,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;oBACnC,iBAAiB,CAAC,YAAY,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC;oBACnE,IAAI,IAAI,CAAC,GAAG,CAAC,YAAY,CAAC,GAAG,SAAS,EAAE,CAAC;wBACxC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBACnB,CAAC;gBACF,CAAC;YACF,CAAC;QACF,CAAC;IACF,CAAC;IAED,OAAO,IAAI,CAAC;AACb,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,WAAW,CAAC,MAAe;IAC1C,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,mFAAmF;IACnF,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IAE/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5B,MAAM,GAAG,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QACtB,eAAe,CAAC,GAAG,EAAE,EAAE,QAAQ,EAAE,CAAC,EAAE,CAAC,CAAC;QAEtC,MAAM,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACnB,iBAAiB,CAAC,GAAG,EAAE,EAAE,QAAQ,EAAE,CAAC,EAAE,WAAW,EAAE,CAAC,EAAE,CAAC,CAAC;QACxD,KAAK,IAAI,GAAG,CAAC;IACd,CAAC;IAED,OAAO,KAAK,CAAC;AACd,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,eAAe,CAAoB,MAAS;IAC3D,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,4DAA4D;IAC5D,MAAM,MAAM,GAAe,EAAE,CAAC;IAE9B,mDAAmD;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACtC,CAAC;IAED,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAI,MAAkB,CAAC,GAAG,CAAC,CAAC;QAC3C,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC3B,iBAAiB,CAAC,GAAG,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAE5D,4CAA4C;YAC5C,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;YAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;YAC9C,SAAS,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QACtB,CAAC;IACF,CAAC;IAED,OAAO,MAA0B,CAAC;AACnC,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,SAAS,CAAoB,MAAS,EAAE,EAAuD;IAC9G,YAAY,CAAC,MAAM,CAAC,CAAC;IAErB,MAAM,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC;IACxC,MAAM,MAAM,GAAe,EAAE,CAAC;IAE9B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;QACrC,MAAM,SAAS,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC;QAC9B,eAAe,CAAC,SAAS,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC;QAE9C,MAAM,MAAM,GAAa,EAAE,CAAC;QAE5B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,EAAE,EAAE,CAAC;YACrC,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC;YAC3B,iBAAiB,CAAC,GAAG,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,WAAW,EAAE,GAAG,EAAE,CAAC,CAAC;YAC5D,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAChC,CAAC;QACD,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACrB,CAAC;IAED,OAAO,MAA0B,CAAC;AACnC,CAAC"}

package/build/matrices/core.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=core.spec.d.ts.map

package/build/matrices/core.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"core.spec.d.ts","sourceRoot":"","sources":["../../src/matrices/core.spec.ts"],"names":[],"mappings":""}