@danielsimonjr/mathts-functions 0.3.0 → 0.5.0

package/dist/index.js

@@ -468,7 +468,7 @@ __export(typed_exports, {
 });
 
 // src/typed/arithmetic.ts
-import { mathTyped, Complex, Fraction, BigNumber } from "@danielsimonjr/mathts-core";
+import { mathTyped, Complex, Fraction, BigNumber, Dual } from "@danielsimonjr/mathts-core";
 import { DenseMatrix, backendManager } from "@danielsimonjr/mathts-matrix";
 import { computePool } from "@danielsimonjr/mathts-parallel";
 
@@ -510,6 +510,10 @@ var add = mathTyped("add", {
   "BigNumber, BigNumber": (a, b) => a.add(b),
   // Unit arithmetic (mathjs parity): same-dimension units add (e.g. 5 cm + 3 mm).
   "Unit, Unit": (a, b) => a.add(b),
+  // Dual numbers (forward-mode AD).
+  "Dual, Dual": (a, b) => a.add(b),
+  "Dual, number": (a, b) => a.add(Dual.constant(b)),
+  "number, Dual": (a, b) => Dual.constant(a).add(b),
   // Mixed type operations (auto-coercion via conversions)
   "number, Complex": (a, b) => Complex.fromNumber(a).add(b),
   "Complex, number": (a, b) => a.add(Complex.fromNumber(b)),
@@ -544,6 +548,10 @@ var subtract = mathTyped("subtract", {
   "BigNumber, BigNumber": (a, b) => a.subtract(b),
   // Unit subtraction (same dimension): 5 cm − 3 mm.
   "Unit, Unit": (a, b) => a.sub(b),
+  // Dual numbers (forward-mode AD).
+  "Dual, Dual": (a, b) => a.sub(b),
+  "Dual, number": (a, b) => a.sub(Dual.constant(b)),
+  "number, Dual": (a, b) => Dual.constant(a).sub(b),
   "number, Complex": (a, b) => Complex.fromNumber(a).subtract(b),
   "Complex, number": (a, b) => a.subtract(Complex.fromNumber(b)),
   "number, Fraction": (a, b) => Fraction.fromNumber(a).subtract(b),
@@ -566,6 +574,10 @@ var multiply = mathTyped("multiply", {
   "Unit, Unit": (a, b) => a.mul(b),
   "Unit, number": (a, b) => a.mul(b),
   "number, Unit": (a, b) => b.mul(a),
+  // Dual numbers (forward-mode AD).
+  "Dual, Dual": (a, b) => a.mul(b),
+  "Dual, number": (a, b) => a.mul(Dual.constant(b)),
+  "number, Dual": (a, b) => Dual.constant(a).mul(b),
   // Matrix × matrix (2-D arrays) routed through the native DenseMatrix +
   // BackendManager — i.e. WASM/GPU acceleration for large matrices (GC7).
   // Previously `multiply(2DArray, 2DArray)` threw; this is additive. Non-2-D
@@ -621,6 +633,10 @@ var divide = mathTyped("divide", {
   // Unit division / scaling (mathjs parity). Unit ÷ Unit may be dimensionless.
   "Unit, Unit": (a, b) => a.div(b),
   "Unit, number": (a, b) => a.div(b),
+  // Dual numbers (forward-mode AD).
+  "Dual, Dual": (a, b) => a.div(b),
+  "Dual, number": (a, b) => a.div(Dual.constant(b)),
+  "number, Dual": (a, b) => Dual.constant(a).div(b),
   "number, Complex": (a, b) => Complex.fromNumber(a).divide(b),
   "Complex, number": (a, b) => a.divide(Complex.fromNumber(b)),
   "number, Fraction": (a, b) => Fraction.fromNumber(a).divide(b),
@@ -660,6 +676,7 @@ var abs = mathTyped("abs", {
   Complex: (a) => a.abs(),
   Fraction: (a) => a.abs(),
   BigNumber: (a) => a.abs(),
+  Dual: (a) => a.abs(),
   // Parallel array abs
   Float64Array: async (a) => {
     const wasm = elementwiseUnaryDispatch("abs", a);
@@ -696,7 +713,10 @@ var pow = mathTyped("pow", {
   "Complex, Complex": (a, b) => a.pow(b),
   "Fraction, number": (a, b) => a.pow(Math.floor(b)),
   "BigNumber, number": (a, b) => a.pow(b),
-  "BigNumber, BigNumber": (a, b) => a.pow(b.valueOf())
+  "BigNumber, BigNumber": (a, b) => a.pow(b.valueOf()),
+  // Dual numbers (forward-mode AD): constant or variable exponent.
+  "Dual, number": (a, b) => a.powConst(b),
+  "Dual, Dual": (a, b) => a.pow(b)
 });
 var sqrt = mathTyped("sqrt", {
   number: (a) => {
@@ -707,6 +727,7 @@ var sqrt = mathTyped("sqrt", {
   },
   Complex: (a) => a.sqrt(),
   BigNumber: (a) => a.sqrt(),
+  Dual: (a) => a.sqrt(),
   // Parallel array sqrt
   Float64Array: async (a) => {
     const result = await computePool.sqrt(a);
@@ -719,6 +740,7 @@ var square = mathTyped("square", {
   Complex: (a) => a.multiply(a),
   Fraction: (a) => a.multiply(a),
   BigNumber: (a) => a.multiply(a),
+  Dual: (a) => a.square(),
   // Parallel array square
   Float64Array: async (a) => {
     const result = await computePool.square(a);
@@ -731,6 +753,7 @@ var cube = mathTyped("cube", {
   Complex: (a) => a.multiply(a).multiply(a),
   Fraction: (a) => a.multiply(a).multiply(a),
   BigNumber: (a) => a.multiply(a).multiply(a),
+  Dual: (a) => a.powConst(3),
   // Parallel array cube
   Float64Array: async (a) => {
     if (computePool.shouldParallelize(a.length)) {
@@ -746,6 +769,7 @@ var cbrt = mathTyped("cbrt", {
   number: (a) => Math.cbrt(a),
   Complex: (a) => a.pow(1 / 3),
   BigNumber: (a) => a.cbrt(),
+  Dual: (a) => a.powConst(1 / 3),
   // Parallel array cbrt
   Float64Array: async (a) => {
     if (computePool.shouldParallelize(a.length)) {
@@ -766,6 +790,7 @@ var exp = mathTyped("exp", {
   number: (a) => Math.exp(a),
   Complex: (a) => a.exp(),
   BigNumber: (a) => a.exp(),
+  Dual: (a) => a.exp(),
   // Parallel array exp
   Float64Array: async (a) => {
     const wasm = elementwiseUnaryDispatch("exp", a);
@@ -778,6 +803,7 @@ var log = mathTyped("log", {
   number: (a) => Math.log(a),
   Complex: (a) => a.log(),
   BigNumber: (a) => a.ln(),
+  Dual: (a) => a.log(),
   "number, number": (a, base) => Math.log(a) / Math.log(base),
   // Parallel array log
   Float64Array: async (a) => {
@@ -1349,6 +1375,7 @@ var sin = mathTyped2("sin", {
   number: (a) => Math.sin(a),
   Complex: (a) => a.sin(),
   BigNumber: (a) => a.sin(),
+  Dual: (a) => a.sin(),
   // Parallel Float64Array sin
   Float64Array: async (a) => {
     const wasm = elementwiseUnaryDispatch("sin", a);
@@ -1361,6 +1388,7 @@ var cos = mathTyped2("cos", {
   number: (a) => Math.cos(a),
   Complex: (a) => a.cos(),
   BigNumber: (a) => a.cos(),
+  Dual: (a) => a.cos(),
   // Parallel Float64Array cos
   Float64Array: async (a) => {
     const wasm = elementwiseUnaryDispatch("cos", a);
@@ -1373,6 +1401,7 @@ var tan = mathTyped2("tan", {
   number: (a) => Math.tan(a),
   Complex: (a) => a.tan(),
   BigNumber: (a) => a.tan(),
+  Dual: (a) => a.tan(),
   // Parallel Float64Array tan
   Float64Array: async (a) => {
     const wasm = elementwiseUnaryDispatch("tan", a);
@@ -14910,7 +14939,18 @@ import {
   createEvaluate,
   compileExpression as _compileExpression
 } from "@danielsimonjr/mathts-expression";
-import { Complex as Complex8, Fraction as Fraction4 } from "@danielsimonjr/mathts-core";
+import {
+  Complex as Complex8,
+  Fraction as Fraction4,
+  I,
+  PHI,
+  SQRT2,
+  SQRT1_2,
+  LN2 as LN22,
+  LN10,
+  LOG2E,
+  LOG10E
+} from "@danielsimonjr/mathts-core";
 
 // src/factories/scope.ts
 import { mathTyped as mathTyped17, Complex as Complex7, Fraction as Fraction3, BigNumber as BigNumber8 } from "@danielsimonjr/mathts-core";
@@ -16761,13 +16801,13 @@ function clone2(array) {
 // src/utils/switch.ts
 function _switch(mat) {
   const rows = mat;
-  const I = rows.length;
+  const I2 = rows.length;
   const J = rows[0].length;
   let i, j;
   const ret = [];
   for (j = 0; j < J; j++) {
     const tmp = [];
-    for (i = 0; i < I; i++) {
+    for (i = 0; i < I2; i++) {
       tmp.push(rows[i][j]);
     }
     ret.push(tmp);
@@ -19717,17 +19757,17 @@ var createDot = /* @__PURE__ */ factory(
       let i = 0;
       let j = 0;
       while (i < xindex.length && j < yindex.length) {
-        const I = xindex[i];
+        const I2 = xindex[i];
         const J = yindex[j];
-        if (I < J) {
+        if (I2 < J) {
           i++;
           continue;
         }
-        if (I > J) {
+        if (I2 > J) {
           j++;
           continue;
         }
-        if (I === J) {
+        if (I2 === J) {
           c = add2(c, mul(xvalues[i], yvalues[j]));
           i++;
           j++;
@@ -20260,7 +20300,7 @@ var dependencies85 = ["typed", "config", "matrix", "BigNumber", "DenseMatrix", "
 var createIdentity = /* @__PURE__ */ factory(
   name85,
   dependencies85,
-  ({ typed: typed3, config, matrix: matrix2, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4, SparseMatrix }) => {
+  ({ typed: typed3, config, matrix: matrix2, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5, SparseMatrix }) => {
     return typed3(name85, {
       "": function() {
         return config.matrix === "Matrix" ? matrix2([]) : [];
@@ -20323,7 +20363,7 @@ var createIdentity = /* @__PURE__ */ factory(
           return SparseMatrix.diagonal(size2, one, 0, defaultValue);
         }
         if (format6 === "dense") {
-          return DenseMatrix4.diagonal(size2, one, 0, defaultValue);
+          return DenseMatrix5.diagonal(size2, one, 0, defaultValue);
         }
         throw new TypeError(`Unknown matrix type "${format6}"`);
       }
@@ -20431,7 +20471,7 @@ var dependencies88 = ["typed", "matrix", "DenseMatrix", "SparseMatrix"];
 var createDiag = /* @__PURE__ */ factory(
   name88,
   dependencies88,
-  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, SparseMatrix }) => {
+  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, SparseMatrix }) => {
     return typed3(name88, {
       // FIXME: simplify this huge amount of signatures as soon as typed-function supports optional arguments
       Array: function(x) {
@@ -20490,7 +20530,7 @@ var createDiag = /* @__PURE__ */ factory(
       if (format6 && format6 !== "sparse" && format6 !== "dense") {
         throw new TypeError(`Unknown matrix type ${format6}"`);
       }
-      const m = format6 === "sparse" ? SparseMatrix.diagonal(ms, x, k) : DenseMatrix4.diagonal(ms, x, k);
+      const m = format6 === "sparse" ? SparseMatrix.diagonal(ms, x, k) : DenseMatrix5.diagonal(ms, x, k);
       return format6 !== null ? m : m.valueOf();
     }
     function _getDiagonal(x, k, format6, s, kSub, kSuper) {
@@ -21052,6 +21092,79 @@ var createDet = /* @__PURE__ */ factory(
   }
 );
 
+// src/matrix/native-accel.ts
+import { DenseMatrix as DenseMatrix4, lu } from "@danielsimonjr/mathts-matrix";
+var NATIVE_MATRIX_THRESHOLD = 8;
+function isLargeNumericSquare(a) {
+  if (!Array.isArray(a) || a.length < NATIVE_MATRIX_THRESHOLD) return false;
+  const n = a.length;
+  for (const row2 of a) {
+    if (!Array.isArray(row2) || row2.length !== n) return false;
+    for (let j = 0; j < n; j++) if (typeof row2[j] !== "number") return false;
+  }
+  return true;
+}
+function permutationParity(P2) {
+  let sign3 = 1;
+  const seen = new Array(P2.length).fill(false);
+  for (let i = 0; i < P2.length; i++) {
+    if (seen[i]) continue;
+    let j = i;
+    let len = 0;
+    while (!seen[j]) {
+      seen[j] = true;
+      j = P2[j];
+      len++;
+    }
+    if (len % 2 === 0) sign3 = -sign3;
+  }
+  return sign3;
+}
+var SINGULAR = /singular|zero pivot/i;
+function acceleratedDet(a, factoryDet2) {
+  if (!isLargeNumericSquare(a)) return factoryDet2(a);
+  try {
+    const { U, P: P2 } = lu(DenseMatrix4.fromArray(a));
+    const Ud = U.toArray();
+    let d = permutationParity(P2);
+    for (let i = 0; i < a.length; i++) d *= Ud[i][i];
+    return d;
+  } catch (e) {
+    if (e instanceof Error && SINGULAR.test(e.message)) return 0;
+    return factoryDet2(a);
+  }
+}
+function acceleratedInv(a, factoryInv2) {
+  if (!isLargeNumericSquare(a)) return factoryInv2(a);
+  try {
+    const { L, U, P: P2 } = lu(DenseMatrix4.fromArray(a));
+    const Ld = L.toArray();
+    const Ud = U.toArray();
+    const n = a.length;
+    const inv2 = Array.from({ length: n }, () => new Array(n).fill(0));
+    for (let j = 0; j < n; j++) {
+      const b = new Array(n).fill(0);
+      for (let i = 0; i < n; i++) if (P2[i] === j) b[i] = 1;
+      const y = new Array(n).fill(0);
+      for (let i = 0; i < n; i++) {
+        let s = b[i];
+        for (let k = 0; k < i; k++) s -= Ld[i][k] * y[k];
+        y[i] = s / Ld[i][i];
+      }
+      const x = new Array(n).fill(0);
+      for (let i = n - 1; i >= 0; i--) {
+        let s = y[i];
+        for (let k = i + 1; k < n; k++) s -= Ud[i][k] * x[k];
+        x[i] = s / Ud[i][i];
+      }
+      for (let i = 0; i < n; i++) inv2[i][j] = x[i];
+    }
+    return inv2;
+  } catch {
+    return factoryInv2(a);
+  }
+}
+
 // src/matrix/reshape.ts
 var name96 = "reshape";
 var dependencies96 = ["typed", "isInteger", "matrix"];
@@ -21958,7 +22071,7 @@ var createCsSymperm = /* @__PURE__ */ factory(
 );
 
 // src/algebra/solver/utils/solveValidation.ts
-function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
+function createSolveValidation({ DenseMatrix: DenseMatrix5 }) {
   return function solveValidation(m, b, copy) {
     const mSize = m.size();
     if (mSize.length !== 2) {
@@ -21981,7 +22094,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
         for (let i = 0; i < rows; i++) {
           data[i] = [bdata[i]];
         }
-        return new DenseMatrix4({
+        return new DenseMatrix5({
           data,
           size: [rows, 1],
           datatype: bm._datatype
@@ -21997,7 +22110,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
             for (let i = 0; i < rows; i++) {
               data[i] = [bdata[i][0]];
             }
-            return new DenseMatrix4({
+            return new DenseMatrix5({
               data,
               size: [rows, 1],
               datatype: bm._datatype
@@ -22016,7 +22129,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
             const i = index[k];
             data[i][0] = values[k];
           }
-          return new DenseMatrix4({
+          return new DenseMatrix5({
             data,
             size: [rows, 1],
             datatype: bm._datatype
@@ -22036,7 +22149,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
         for (let i = 0; i < rows; i++) {
           data[i] = [b[i]];
         }
-        return new DenseMatrix4({
+        return new DenseMatrix5({
           data,
           size: [rows, 1]
         });
@@ -22048,7 +22161,7 @@ function createSolveValidation({ DenseMatrix: DenseMatrix4 }) {
         for (let i = 0; i < rows; i++) {
           data[i] = [b[i][0]];
         }
-        return new DenseMatrix4({
+        return new DenseMatrix5({
           data,
           size: [rows, 1]
         });
@@ -22118,10 +22231,10 @@ var createLsolve = /* @__PURE__ */ factory(
     multiplyScalar: multiplyScalar2,
     subtractScalar: subtractScalar2,
     equalScalar: equalScalar3,
-    DenseMatrix: DenseMatrix4
+    DenseMatrix: DenseMatrix5
   }) => {
     const solveValidation = createSolveValidation({
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     return typed3(name101, {
       "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22159,7 +22272,7 @@ var createLsolve = /* @__PURE__ */ factory(
             for (let i = 0; i < rows; i++) {
               x2[i] = [resultAlloc.array[i]];
             }
-            return new DenseMatrix4({
+            return new DenseMatrix5({
               data: x2,
               size: [rows, 1]
             });
@@ -22192,7 +22305,7 @@ var createLsolve = /* @__PURE__ */ factory(
         }
         x[j] = [xj];
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: x,
         size: [rows, 1]
       });
@@ -22236,7 +22349,7 @@ var createLsolve = /* @__PURE__ */ factory(
           x[j] = [0];
         }
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: x,
         size: [rows, 1]
       });
@@ -22265,10 +22378,10 @@ var createLsolveAll = /* @__PURE__ */ factory(
     multiplyScalar: multiplyScalar2,
     subtractScalar: subtractScalar2,
     equalScalar: equalScalar3,
-    DenseMatrix: DenseMatrix4
+    DenseMatrix: DenseMatrix5
   }) => {
     const solveValidation = createSolveValidation({
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     return typed3(name102, {
       "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22318,7 +22431,7 @@ var createLsolveAll = /* @__PURE__ */ factory(
         }
       }
       return B.map(
-        (x) => new DenseMatrix4({
+        (x) => new DenseMatrix5({
           data: x.map((e) => [e]),
           size: [rows, 1]
         })
@@ -22377,7 +22490,7 @@ var createLsolveAll = /* @__PURE__ */ factory(
         }
       }
       return B.map(
-        (x) => new DenseMatrix4({
+        (x) => new DenseMatrix5({
           data: x.map((e) => [e]),
           size: [rows, 1]
         })
@@ -22443,10 +22556,10 @@ var createUsolve = /* @__PURE__ */ factory(
     multiplyScalar: multiplyScalar2,
     subtractScalar: subtractScalar2,
     equalScalar: equalScalar3,
-    DenseMatrix: DenseMatrix4
+    DenseMatrix: DenseMatrix5
   }) => {
     const solveValidation = createSolveValidation({
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     return typed3(name103, {
       "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22484,7 +22597,7 @@ var createUsolve = /* @__PURE__ */ factory(
             for (let i = 0; i < rows; i++) {
               x2[i] = [resultAlloc.array[i]];
             }
-            return new DenseMatrix4({
+            return new DenseMatrix5({
               data: x2,
               size: [rows, 1]
             });
@@ -22517,7 +22630,7 @@ var createUsolve = /* @__PURE__ */ factory(
         }
         x[j] = [xj];
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: x,
         size: [rows, 1]
       });
@@ -22561,7 +22674,7 @@ var createUsolve = /* @__PURE__ */ factory(
           x[j] = [0];
         }
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: x,
         size: [rows, 1]
       });
@@ -22590,10 +22703,10 @@ var createUsolveAll = /* @__PURE__ */ factory(
     multiplyScalar: multiplyScalar2,
     subtractScalar: subtractScalar2,
     equalScalar: equalScalar3,
-    DenseMatrix: DenseMatrix4
+    DenseMatrix: DenseMatrix5
   }) => {
     const solveValidation = createSolveValidation({
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     return typed3(name104, {
       "SparseMatrix, Array | Matrix": function(m, b) {
@@ -22643,7 +22756,7 @@ var createUsolveAll = /* @__PURE__ */ factory(
         }
       }
       return B.map(
-        (x) => new DenseMatrix4({
+        (x) => new DenseMatrix5({
           data: x.map((e) => [e]),
           size: [rows, 1]
         })
@@ -22702,7 +22815,7 @@ var createUsolveAll = /* @__PURE__ */ factory(
         }
       }
       return B.map(
-        (x) => new DenseMatrix4({
+        (x) => new DenseMatrix5({
           data: x.map((e) => [e]),
           size: [rows, 1]
         })
@@ -22920,7 +23033,7 @@ var dependencies108 = ["typed", "DenseMatrix"];
 var createMatAlgo12xSfs = /* @__PURE__ */ factory(
   name108,
   dependencies108,
-  ({ typed: typed3, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, DenseMatrix: DenseMatrix5 }) => {
     return function matAlgo12xSfs(s, b, callback, inverse) {
       const avalues = s._values;
       const aindex = s._index;
@@ -22960,7 +23073,7 @@ var createMatAlgo12xSfs = /* @__PURE__ */ factory(
           }
         }
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: cdata,
         size: [rows, columns],
         datatype: dt
@@ -23033,9 +23146,9 @@ var dependencies110 = [
 var createRound = /* @__PURE__ */ factory(
   name110,
   dependencies110,
-  ({ typed: typed3, config, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, config, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5 }) => {
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     function toExponent(epsilon) {
       return Math.abs(splitNumber(epsilon).exponent);
@@ -23619,13 +23732,13 @@ function _mapSlices(mat, dim, callback) {
   }
 }
 function _switch2(mat) {
-  const I = mat.length;
+  const I2 = mat.length;
   const J = mat[0].length;
   let i, j;
   const ret = [];
   for (j = 0; j < J; j++) {
     const tmp = [];
-    for (i = 0; i < I; i++) {
+    for (i = 0; i < I2; i++) {
       tmp.push(mat[i][j]);
     }
     ret.push(tmp);
@@ -24444,13 +24557,13 @@ var createEqual = /* @__PURE__ */ factory(
     typed: typed3,
     matrix: matrix2,
     equalScalar: equalScalar3,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: _concat,
     SparseMatrix
   }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({ typed: typed3, matrix: matrix2 });
     return typed3(
       name124,
@@ -24840,7 +24953,7 @@ var createDotDivide = /* @__PURE__ */ factory(
     matrix: matrix2,
     equalScalar: equalScalar3,
     divideScalar: divideScalar2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3,
     SparseMatrix
   }) => {
@@ -24848,7 +24961,7 @@ var createDotDivide = /* @__PURE__ */ factory(
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -25017,9 +25130,9 @@ var createFloorNumber = /* @__PURE__ */ factory(
 var createFloor = /* @__PURE__ */ factory(
   name133,
   dependencies133,
-  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const floorNumber = createFloorNumber({ typed: typed3, config, round: round2 });
     function _bigFloor(x) {
@@ -25230,7 +25343,7 @@ var dependencies135 = [
 var createMod = /* @__PURE__ */ factory(
   name135,
   dependencies135,
-  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const floor2 = createFloor({
       typed: typed3,
       config,
@@ -25238,13 +25351,13 @@ var createMod = /* @__PURE__ */ factory(
       matrix: matrix2,
       equalScalar: equalScalar3,
       zeros: zeros3,
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -25457,7 +25570,7 @@ var dependencies138 = ["typed", "DenseMatrix"];
 var createMatAlgo10xSids = /* @__PURE__ */ factory(
   name138,
   dependencies138,
-  ({ typed: typed3, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, DenseMatrix: DenseMatrix5 }) => {
     return function matAlgo10xSids(s, b, callback, inverse) {
       const avalues = s._values;
       const aindex = s._index;
@@ -25497,7 +25610,7 @@ var createMatAlgo10xSids = /* @__PURE__ */ factory(
           }
         }
       }
-      return new DenseMatrix4({
+      return new DenseMatrix5({
         data: cdata,
         size: [rows, columns],
         datatype: dt
@@ -25536,7 +25649,7 @@ var createGcd = /* @__PURE__ */ factory(
     equalScalar: equalScalar3,
     zeros: zeros3,
     BigNumber: BigNumber9,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3
   }) => {
     const mod2 = createMod({
@@ -25546,12 +25659,12 @@ var createGcd = /* @__PURE__ */ factory(
       matrix: matrix2,
       equalScalar: equalScalar3,
       zeros: zeros3,
-      DenseMatrix: DenseMatrix4,
+      DenseMatrix: DenseMatrix5,
       concat: concat3
     });
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo04xSidSid = createMatAlgo04xSidSid({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -26118,9 +26231,9 @@ var createCeilNumber = /* @__PURE__ */ factory(
 var createCeil = /* @__PURE__ */ factory(
   name145,
   dependencies145,
-  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, config, round: round2, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const ceilNumber = createCeilNumber({ typed: typed3, config, round: round2 });
     function _bigCeil(x) {
@@ -26254,10 +26367,10 @@ var dependencies147 = ["typed", "matrix", "equalScalar", "DenseMatrix", "concat"
 var createBitOr = /* @__PURE__ */ factory(
   name147,
   dependencies147,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo04xSidSid = createMatAlgo04xSidSid({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -26285,10 +26398,10 @@ var dependencies148 = ["typed", "matrix", "DenseMatrix", "concat", "SparseMatrix
 var createBitXor = /* @__PURE__ */ factory(
   name148,
   dependencies148,
-  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
+  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -26429,11 +26542,11 @@ var dependencies150 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
 var createLeftShift = /* @__PURE__ */ factory(
   name150,
   dependencies150,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26500,11 +26613,11 @@ var dependencies151 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
 var createRightArithShift = /* @__PURE__ */ factory(
   name151,
   dependencies151,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26571,11 +26684,11 @@ var dependencies152 = ["typed", "matrix", "equalScalar", "zeros", "DenseMatrix",
 var createRightLogShift = /* @__PURE__ */ factory(
   name152,
   dependencies152,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo08xS0Sid = createMatAlgo08xS0Sid({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
@@ -26641,10 +26754,10 @@ var dependencies153 = ["typed", "matrix", "equalScalar", "DenseMatrix", "concat"
 var createOr = /* @__PURE__ */ factory(
   name153,
   dependencies153,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -26680,10 +26793,10 @@ var dependencies154 = ["typed", "matrix", "DenseMatrix", "concat", "SparseMatrix
 var createXor = /* @__PURE__ */ factory(
   name154,
   dependencies154,
-  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
+  ({ typed: typed3, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -26870,12 +26983,12 @@ var createCompare = /* @__PURE__ */ factory(
     matrix: matrix2,
     BigNumber: BigNumber9,
     Fraction: Fraction5,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3
   }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27000,10 +27113,10 @@ var dependencies159 = [
 var createLarger = /* @__PURE__ */ factory(
   name159,
   dependencies159,
-  ({ typed: typed3, config, bignumber: bignumber2, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
+  ({ typed: typed3, config, bignumber: bignumber2, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27058,10 +27171,10 @@ var dependencies160 = ["typed", "config", "matrix", "DenseMatrix", "concat", "Sp
 var createLargerEq = /* @__PURE__ */ factory(
   name160,
   dependencies160,
-  ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
+  ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27124,13 +27237,13 @@ var createSmaller = /* @__PURE__ */ factory(
     config,
     bignumber: bignumber2,
     matrix: matrix2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3,
     SparseMatrix
   }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27185,10 +27298,10 @@ var dependencies162 = ["typed", "config", "matrix", "DenseMatrix", "concat", "Sp
 var createSmallerEq = /* @__PURE__ */ factory(
   name162,
   dependencies162,
-  ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix4, concat: concat3, SparseMatrix }) => {
+  ({ typed: typed3, config, matrix: matrix2, DenseMatrix: DenseMatrix5, concat: concat3, SparseMatrix }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27249,13 +27362,13 @@ var createUnequal = /* @__PURE__ */ factory(
     config: _config,
     equalScalar: equalScalar3,
     matrix: matrix2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3,
     SparseMatrix
   }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27306,12 +27419,12 @@ var dependencies164 = ["typed", "matrix", "equalScalar", "BigNumber", "DenseMatr
 var createAtan2 = /* @__PURE__ */ factory(
   name164,
   dependencies164,
-  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix4, concat: concat3 }) => {
+  ({ typed: typed3, matrix: matrix2, equalScalar: equalScalar3, BigNumber: BigNumber9, DenseMatrix: DenseMatrix5, concat: concat3 }) => {
     const matAlgo02xDS0 = createMatAlgo02xDS0({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo09xS0Sf = createMatAlgo09xS0Sf({ typed: typed3, equalScalar: equalScalar3 });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27380,14 +27493,14 @@ var createDotPow = /* @__PURE__ */ factory(
     equalScalar: equalScalar3,
     matrix: matrix2,
     pow: pow2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3,
     SparseMatrix
   }) => {
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo07xSSf = createMatAlgo07xSSf({ typed: typed3, SparseMatrix });
     const matAlgo11xS0s = createMatAlgo11xS0s({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -27444,8 +27557,8 @@ var createFixNumber = /* @__PURE__ */ factory(
 var createFix = /* @__PURE__ */ factory(
   name167,
   dependencies167,
-  ({ typed: typed3, Complex: Complex10, matrix: matrix2, ceil: ceil2, floor: floor2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix4 }) => {
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+  ({ typed: typed3, Complex: Complex10, matrix: matrix2, ceil: ceil2, floor: floor2, equalScalar: equalScalar3, zeros: zeros3, DenseMatrix: DenseMatrix5 }) => {
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matAlgo14xDs = createMatAlgo14xDs({ typed: typed3 });
     const fixNumber = createFixNumber({ typed: typed3, ceil: ceil2, floor: floor2 });
     return typed3("fix", {
@@ -30448,7 +30561,7 @@ var createFft = /* @__PURE__ */ factory(
     divideScalar: divideScalar2,
     exp: exp2,
     tau: tau2,
-    i: I,
+    i: I2,
     dotDivide: dotDivide2,
     conj: conj3,
     pow: pow2,
@@ -30496,7 +30609,7 @@ var createFft = /* @__PURE__ */ factory(
     }
     function _czt(arr) {
       const n = arr.length;
-      const w = exp2(divideScalar2(multiplyScalar2(-1, multiplyScalar2(I, tau2)), n));
+      const w = exp2(divideScalar2(multiplyScalar2(-1, multiplyScalar2(I2, tau2)), n));
       const chirp = [];
       for (let i = 1 - n; i < n; i++) {
         chirp.push(pow2(w, divideScalar2(pow2(i, 2), 2)));
@@ -30563,7 +30676,7 @@ var createFft = /* @__PURE__ */ factory(
           const p = ret[k];
           const q = multiplyScalar2(
             ret[k + length / 2],
-            exp2(multiplyScalar2(multiplyScalar2(tau2, I), divideScalar2(-k, length)))
+            exp2(multiplyScalar2(multiplyScalar2(tau2, I2), divideScalar2(-k, length)))
           );
           ret[k] = addScalar2(p, q);
           ret[k + length / 2] = addScalar2(p, multiplyScalar2(-1, q));
@@ -30709,15 +30822,15 @@ var createSubtract = /* @__PURE__ */ factory(
     equalScalar: equalScalar3,
     subtractScalar: subtractScalar2,
     unaryMinus: _unaryMinus,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     concat: concat3,
     nodeOperations: nodeOperations2
   }) => {
     const matAlgo01xDSid = createMatAlgo01xDSid({ typed: typed3 });
     const matAlgo03xDSf = createMatAlgo03xDSf({ typed: typed3 });
     const matAlgo05xSfSf = createMatAlgo05xSfSf({ typed: typed3, equalScalar: equalScalar3 });
-    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix4 });
-    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix4 });
+    const matAlgo10xSids = createMatAlgo10xSids({ typed: typed3, DenseMatrix: DenseMatrix5 });
+    const matAlgo12xSfs = createMatAlgo12xSfs({ typed: typed3, DenseMatrix: DenseMatrix5 });
     const matrixAlgorithmSuite = createMatrixAlgorithmSuite({
       typed: typed3,
       matrix: matrix2,
@@ -34457,8 +34570,8 @@ var dependencies207 = ["smaller", "DenseMatrix"];
 var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
   name207,
   dependencies207,
-  ({ smaller: smaller2, DenseMatrix: DenseMatrix4 }) => {
-    class ImmutableDenseMatrix extends DenseMatrix4 {
+  ({ smaller: smaller2, DenseMatrix: DenseMatrix5 }) => {
+    class ImmutableDenseMatrix extends DenseMatrix5 {
       /**
        * Type identifier
        */
@@ -34498,7 +34611,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
           throw new Error("Invalid datatype: " + datatype);
         }
         if (isMatrix(data) || isArray(data)) {
-          const matrix2 = new DenseMatrix4(data, datatype);
+          const matrix2 = new DenseMatrix5(data, datatype);
           this._data = matrix2._data;
           this._size = matrix2._size;
           this._datatype = matrix2._datatype;
@@ -34537,7 +34650,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
       subset(index, _replacement, _defaultValue) {
         switch (arguments.length) {
           case 1: {
-            const m = DenseMatrix4.prototype.subset.call(this, index);
+            const m = DenseMatrix5.prototype.subset.call(this, index);
             if (isMatrix(m)) {
               const mm = m;
               return new ImmutableDenseMatrix({
@@ -34649,7 +34762,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
         if (this._min === null) {
           let m = null;
           const smallerFn = smaller2;
-          DenseMatrix4.prototype.forEach.call(this, function(v) {
+          DenseMatrix5.prototype.forEach.call(this, function(v) {
             if (m === null || smallerFn(v, m)) {
               m = v;
             }
@@ -34666,7 +34779,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
         if (this._max === null) {
           let m = null;
           const smallerFn = smaller2;
-          DenseMatrix4.prototype.forEach.call(this, function(v) {
+          DenseMatrix5.prototype.forEach.call(this, function(v) {
             if (m === null || smallerFn(m, v)) {
               m = v;
             }
@@ -34676,7 +34789,7 @@ var createImmutableDenseMatrixClass = /* @__PURE__ */ factory(
         return this._max ?? void 0;
       }
     }
-    Object.setPrototypeOf(ImmutableDenseMatrix.prototype, DenseMatrix4.prototype);
+    Object.setPrototypeOf(ImmutableDenseMatrix.prototype, DenseMatrix5.prototype);
     const proto = ImmutableDenseMatrix.prototype;
     proto.constructor = ImmutableDenseMatrix;
     proto.type = "ImmutableDenseMatrix";
@@ -35056,7 +35169,7 @@ var createLup = /* @__PURE__ */ factory(
     larger: larger2,
     equalScalar: equalScalar3,
     unaryMinus: unaryMinus2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     SparseMatrix,
     Spa
   }) => {
@@ -35115,11 +35228,11 @@ var createLup = /* @__PURE__ */ factory(
             for (let i2 = 0; i2 < rows; i2++) {
               p2[i2] = permAlloc.array[i2];
             }
-            const l2 = new DenseMatrix4({
+            const l2 = new DenseMatrix5({
               data: ldata2,
               size: [rows, n]
             });
-            const u2 = new DenseMatrix4({
+            const u2 = new DenseMatrix5({
               data: udata2,
               size: [n, columns]
             });
@@ -35176,7 +35289,7 @@ var createLup = /* @__PURE__ */ factory(
         }
         if (j !== pi) {
           p[j] = [p[pi], p[pi] = p[j]][0];
-          DenseMatrix4._swapRows(j, pi, data);
+          DenseMatrix5._swapRows(j, pi, data);
         }
         if (j < rows) {
           for (i = j + 1; i < rows; i++) {
@@ -35221,11 +35334,11 @@ var createLup = /* @__PURE__ */ factory(
           }
         }
       }
-      const l = new DenseMatrix4({
+      const l = new DenseMatrix5({
         data: ldata,
         size: lsize
       });
-      const u = new DenseMatrix4({
+      const u = new DenseMatrix5({
         data: udata,
         size: usize
       });
@@ -36763,7 +36876,7 @@ var dependencies227 = ["typed", "size", "subset", "compareNatural", "Index", "De
 var createSetCartesian = /* @__PURE__ */ factory(
   name227,
   dependencies227,
-  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
     return typed3(name227, {
       "Array | Matrix, Array | Matrix": function(a1, a2) {
         let result = [];
@@ -36780,7 +36893,7 @@ var createSetCartesian = /* @__PURE__ */ factory(
         if (Array.isArray(a1) && Array.isArray(a2)) {
           return result;
         }
-        return new DenseMatrix4(result);
+        return new DenseMatrix5(result);
       }
     });
   }
@@ -36792,7 +36905,7 @@ var dependencies228 = ["typed", "size", "subset", "compareNatural", "Index", "De
 var createSetDifference = /* @__PURE__ */ factory(
   name228,
   dependencies228,
-  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
     return typed3(name228, {
       "Array | Matrix, Array | Matrix": function(a1, a2) {
         let result;
@@ -36821,7 +36934,7 @@ var createSetDifference = /* @__PURE__ */ factory(
         if (Array.isArray(a1) && Array.isArray(a2)) {
           return generalize(result);
         }
-        return new DenseMatrix4(generalize(result));
+        return new DenseMatrix5(generalize(result));
       }
     });
   }
@@ -36833,7 +36946,7 @@ var dependencies229 = ["typed", "size", "subset", "compareNatural", "Index", "De
 var createSetDistinct = /* @__PURE__ */ factory(
   name229,
   dependencies229,
-  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
     return typed3(name229, {
       "Array | Matrix": function(a) {
         let result;
@@ -36852,7 +36965,7 @@ var createSetDistinct = /* @__PURE__ */ factory(
         if (Array.isArray(a)) {
           return result;
         }
-        return new DenseMatrix4(result);
+        return new DenseMatrix5(result);
       }
     });
   }
@@ -36864,7 +36977,7 @@ var dependencies230 = ["typed", "size", "subset", "compareNatural", "Index", "De
 var createSetIntersect = /* @__PURE__ */ factory(
   name230,
   dependencies230,
-  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, size: size2, subset: subset2, compareNatural: compareNatural3, Index: Index2, DenseMatrix: DenseMatrix5 }) => {
     return typed3(name230, {
       "Array | Matrix, Array | Matrix": function(a1, a2) {
         let result;
@@ -36886,7 +36999,7 @@ var createSetIntersect = /* @__PURE__ */ factory(
         if (Array.isArray(a1) && Array.isArray(a2)) {
           return generalize(result);
         }
-        return new DenseMatrix4(generalize(result));
+        return new DenseMatrix5(generalize(result));
       }
     });
   }
@@ -38552,15 +38665,15 @@ function createComplexEigs({
     for (let i = 0; i < N; i++) {
       M[i] = Array(N).fill(0);
     }
-    let I = 0;
+    let I2 = 0;
     for (const sub of arr) {
       const n = sub.length;
       for (let i = 0; i < n; i++) {
         for (let j = 0; j < n; j++) {
-          M[I + i][I + j] = sub[i][j];
+          M[I2 + i][I2 + j] = sub[i][j];
         }
       }
-      I += n;
+      I2 += n;
     }
     return M;
   }
@@ -39309,9 +39422,9 @@ var dependencies240 = ["typed", "matrix", "lup", "slu", "usolve", "lsolve", "Den
 var createLusolve = /* @__PURE__ */ factory(
   name240,
   dependencies240,
-  ({ typed: typed3, matrix: matrix2, lup: lup2, slu: slu2, usolve: usolve2, lsolve: lsolve2, DenseMatrix: DenseMatrix4 }) => {
+  ({ typed: typed3, matrix: matrix2, lup: lup2, slu: slu2, usolve: usolve2, lsolve: lsolve2, DenseMatrix: DenseMatrix5 }) => {
     const solveValidation = createSolveValidation({
-      DenseMatrix: DenseMatrix4
+      DenseMatrix: DenseMatrix5
     });
     return typed3(name240, {
       "Array, Array | Matrix": function(a, b) {
@@ -39394,7 +39507,7 @@ var createLusolve = /* @__PURE__ */ factory(
                 for (let i = 0; i < rows; i++) {
                   x[i] = [resultAlloc.array[i]];
                 }
-                return new DenseMatrix4({
+                return new DenseMatrix5({
                   data: x,
                   size: [rows, 1]
                 });
@@ -41856,7 +41969,7 @@ var createRotationMatrix = /* @__PURE__ */ factory(
     norm: norm3,
     BigNumber: BigNumber9,
     matrix: matrix2,
-    DenseMatrix: DenseMatrix4,
+    DenseMatrix: DenseMatrix5,
     SparseMatrix,
     cos: cos2,
     sin: sin2
@@ -41920,7 +42033,7 @@ var createRotationMatrix = /* @__PURE__ */ factory(
           return new SparseMatrix(data);
         }
         if (format6 === "dense") {
-          return new DenseMatrix4(data);
+          return new DenseMatrix5(data);
         }
         throw new TypeError(`Unknown matrix type "${format6}"`);
       }
@@ -42687,8 +42800,9 @@ factoryScope.add = factoryScope.addScalar;
 var trace2 = createTrace(factoryScope);
 factoryScope.trace = trace2;
 factoryScope.multiply = factoryScope.multiplyScalar;
-var det = createDet(factoryScope);
-factoryScope.det = det;
+var factoryDet = createDet(factoryScope);
+factoryScope.det = factoryDet;
+var det = ((m) => acceleratedDet(m, factoryDet));
 factoryScope.isInteger = (x) => typeof x === "number" && Number.isInteger(x);
 var reshape2 = createReshape(factoryScope);
 factoryScope.reshape = reshape2;
@@ -42704,8 +42818,9 @@ var resize2 = createResize(factoryScope);
 factoryScope.resize = resize2;
 var subset = createSubset(factoryScope);
 factoryScope.subset = subset;
-var inv = createInv(factoryScope);
-factoryScope.inv = inv;
+var factoryInv = createInv(factoryScope);
+factoryScope.inv = factoryInv;
+var inv = ((m) => acceleratedInv(m, factoryInv));
 var factory_cbrt = createCbrt(factoryScope);
 factoryScope.cbrt = factory_cbrt;
 var nthRoots = createNthRoots(factoryScope);
@@ -43145,6 +43260,14 @@ var mathScope = {
   pi: Math.PI,
   e: Math.E,
   tau: 2 * Math.PI,
+  phi: PHI,
+  i: I,
+  SQRT2,
+  SQRT1_2,
+  LN2: LN22,
+  LN10,
+  LOG2E,
+  LOG10E,
   Infinity: Infinity,
   NaN: NaN,
   null: null,
@@ -44950,6 +45073,24 @@ function help(search) {
   }
   return new Help(doc);
 }
+
+// src/grad-forward.ts
+import { Dual as Dual3 } from "@danielsimonjr/mathts-core";
+function derivativeAt(fn, x0) {
+  return fn(Dual3.variable(x0)).deriv;
+}
+function valueAndDerivativeAt(fn, x0) {
+  const r = fn(Dual3.variable(x0));
+  return { value: r.value, deriv: r.deriv };
+}
+function gradientAt(fn, x) {
+  const grad = new Array(x.length);
+  for (let k = 0; k < x.length; k++) {
+    const duals = x.map((xi, i) => i === k ? Dual3.variable(xi) : Dual3.constant(xi));
+    grad[k] = fn(duals).deriv;
+  }
+  return grad;
+}
 export {
   ARRAY_WORKER_THRESHOLD,
   CAS_BATCH_THRESHOLD,
@@ -45109,6 +45250,7 @@ export {
   degree,
   delaunayTriangulation,
   derivative,
+  derivativeAt,
   det,
   deuteronMass,
   diag,
@@ -45279,6 +45421,7 @@ export {
   gpuMatmul,
   gpuScale,
   gpuTranspose,
+  gradientAt,
   gradientSymbolic,
   graphDistance,
   gravitationConstant,
@@ -45678,6 +45821,7 @@ export {
   usolveAll,
   vacuumImpedance,
   validateClosureSource,
+  valueAndDerivativeAt,
   variables,
   variance,
   voronoiDiagram,