@danielsimonjr/mathts-compat 0.4.29 → 0.4.31

package/dist/index.js

@@ -1,621 +1,376 @@
-// src/shims.ts
-import {
-  Complex,
-  Fraction,
-  BigNumber,
-  I,
-  COMPLEX_ZERO,
-  isComplex,
-  isFraction,
-  isBigNumber,
-  isNumber,
-  LN2 as CORE_LN2,
-  LN10 as CORE_LN10,
-  LOG2E as CORE_LOG2E,
-  LOG10E as CORE_LOG10E,
-  SQRT2 as CORE_SQRT2,
-  SQRT1_2 as CORE_SQRT1_2
-} from "@danielsimonjr/mathts-core";
-import {
-  add as _add,
-  subtract as _subtract,
-  multiply as _multiply,
-  divide as _divide,
-  pow as _pow,
-  sqrt as _sqrt,
-  abs as _abs,
-  exp as _exp,
-  log as _log,
-  sin as _sin,
-  cos as _cos,
-  tan as _tan,
-  sum as _sum,
-  mean as _mean,
-  min as _min,
-  max as _max,
-  gcd as _gcd,
-  lcm as _lcm,
-  round as _round,
-  floor as _floor,
-  ceil as _ceil
-} from "@danielsimonjr/mathts-functions";
-import {
-  DenseMatrix,
-  SparseMatrix,
-  add as _matAdd,
-  subtract as _matSubtract,
-  multiply as _matMultiply
-} from "@danielsimonjr/mathts-matrix";
-function complex(re2, im2) {
-  if (re2 === void 0) {
-    return COMPLEX_ZERO;
-  }
-  if (re2 instanceof Complex) {
-    return re2;
-  }
-  if (typeof re2 === "string") {
-    return Complex.parse(re2);
-  }
-  return new Complex(re2, im2 ?? 0);
+import { BIGNUMBER_E, BIGNUMBER_ONE, BIGNUMBER_PI, BIGNUMBER_ZERO, BigNumber, BigNumber as BigNumber$1, COMPLEX_ONE, COMPLEX_ZERO, COMPLEX_ZERO as COMPLEX_ZERO$1, Complex, Complex as Complex$1, FRACTION_ONE, FRACTION_ZERO, Fraction, Fraction as Fraction$1, I, I as I$1, LN10 as LN10$1, LN2 as LN2$1, LOG10E as LOG10E$1, LOG2E as LOG2E$1, SQRT1_2 as SQRT1_2$1, SQRT2 as SQRT2$1, isBigNumber, isComplex, isFraction, isNumber } from "@danielsimonjr/mathts-core";
+import * as mathFunctions from "@danielsimonjr/mathts-functions";
+import { abs as abs$1, add as add$1, ceil as ceil$1, cos as cos$1, divide as divide$1, exp as exp$1, floor as floor$1, gcd as gcd$1, lcm as lcm$1, log as log$1, max as max$1, mean as mean$1, min as min$1, multiply as multiply$1, pow as pow$1, round as round$1, sin as sin$1, sqrt as sqrt$1, subtract as subtract$1, sum as sum$1, tan as tan$1 } from "@danielsimonjr/mathts-functions";
+import { DenseMatrix, DenseMatrix as DenseMatrix$1, SparseMatrix, SparseMatrix as SparseMatrix$1, add as add$2, multiply as multiply$2, subtract as subtract$2 } from "@danielsimonjr/mathts-matrix";
+import { computePool } from "@danielsimonjr/mathts-parallel";
+//#region src/shims.ts
+function complex(re, im) {
+	if (re === void 0) return COMPLEX_ZERO$1;
+	if (re instanceof Complex$1) return re;
+	if (typeof re === "string") return Complex$1.parse(re);
+	return new Complex$1(re, im ?? 0);
 }
 function fraction(numerator, denominator) {
-  if (numerator === void 0) {
-    return new Fraction(0, 1);
-  }
-  if (numerator instanceof Fraction) {
-    return numerator;
-  }
-  if (typeof numerator === "string") {
-    return Fraction.parse(numerator);
-  }
-  if (denominator !== void 0) {
-    return new Fraction(numerator, denominator);
-  }
-  return Fraction.fromNumber(numerator);
+	if (numerator === void 0) return new Fraction$1(0, 1);
+	if (numerator instanceof Fraction$1) return numerator;
+	if (typeof numerator === "string") return Fraction$1.parse(numerator);
+	if (denominator !== void 0) return new Fraction$1(numerator, denominator);
+	return Fraction$1.fromNumber(numerator);
 }
 function bignumber(value) {
-  if (value === void 0) {
-    return BigNumber.fromNumber(0);
-  }
-  if (value instanceof BigNumber) {
-    return value;
-  }
-  if (typeof value === "string") {
-    return BigNumber.parse(value);
-  }
-  return BigNumber.fromNumber(value);
+	if (value === void 0) return BigNumber$1.fromNumber(0);
+	if (value instanceof BigNumber$1) return value;
+	if (typeof value === "string") return BigNumber$1.parse(value);
+	return BigNumber$1.fromNumber(value);
 }
 function matrix(data, format) {
-  if (data === void 0) {
-    return DenseMatrix.zeros(0, 0);
-  }
-  if (data instanceof DenseMatrix || data instanceof SparseMatrix) {
-    return data;
-  }
-  const dense = DenseMatrix.fromArray(data);
-  if (format === "sparse") {
-    return SparseMatrix.fromDense(dense);
-  }
-  return dense;
+	if (data === void 0) return DenseMatrix$1.zeros(0, 0);
+	if (data instanceof DenseMatrix$1 || data instanceof SparseMatrix$1) return data;
+	const dense = DenseMatrix$1.fromArray(data);
+	if (format === "sparse") return SparseMatrix$1.fromDense(dense);
+	return dense;
 }
 function sparse(data) {
-  if (data === void 0) {
-    return SparseMatrix.zeros(0, 0);
-  }
-  return SparseMatrix.fromDense(DenseMatrix.fromArray(data));
+	if (data === void 0) return SparseMatrix$1.zeros(0, 0);
+	return SparseMatrix$1.fromDense(DenseMatrix$1.fromArray(data));
 }
-var isMatrixLike = (x) => x instanceof DenseMatrix || Array.isArray(x) && x.length > 0 && Array.isArray(x[0]);
-var toDenseMatrix = (x) => x instanceof DenseMatrix ? x : DenseMatrix.fromArray(x);
-var unwrap = (r) => r instanceof DenseMatrix ? r.toArray() : r;
-var applyArithmetic = (scalarFn, matrixFn, a, b) => {
-  if (isMatrixLike(a) || isMatrixLike(b)) {
-    const aa = isMatrixLike(a) ? toDenseMatrix(a) : a;
-    const bb = isMatrixLike(b) ? toDenseMatrix(b) : b;
-    return unwrap(matrixFn(aa, bb));
-  }
-  return scalarFn(a, b);
+const isMatrixLike = (x) => x instanceof DenseMatrix$1 || Array.isArray(x) && x.length > 0 && Array.isArray(x[0]);
+const toDenseMatrix = (x) => x instanceof DenseMatrix$1 ? x : DenseMatrix$1.fromArray(x);
+const unwrap = (r) => r instanceof DenseMatrix$1 ? r.toArray() : r;
+const applyArithmetic = (scalarFn, matrixFn, a, b) => {
+	if (isMatrixLike(a) || isMatrixLike(b)) {
+		const aa = isMatrixLike(a) ? toDenseMatrix(a) : a;
+		const bb = isMatrixLike(b) ? toDenseMatrix(b) : b;
+		return unwrap(matrixFn(aa, bb));
+	}
+	return scalarFn(a, b);
 };
-var makeArithmetic = (scalarFn, matrixFn) => (...args) => {
-  if (args.length === 0) {
-    throw new TypeError("Too few arguments (expected at least 1)");
-  }
-  if (args.length === 1) {
-    return args[0];
-  }
-  return args.reduce((acc, next) => applyArithmetic(scalarFn, matrixFn, acc, next));
+const makeArithmetic = (scalarFn, matrixFn) => (...args) => {
+	if (args.length === 0) throw new TypeError("Too few arguments (expected at least 1)");
+	if (args.length === 1) return args[0];
+	return args.reduce((acc, next) => applyArithmetic(scalarFn, matrixFn, acc, next));
 };
-var add = makeArithmetic(
-  _add,
-  _matAdd
-);
-var subtract = makeArithmetic(
-  _subtract,
-  _matSubtract
-);
-var multiply = makeArithmetic(
-  _multiply,
-  _matMultiply
-);
-var divide = _divide;
-var pow = _pow;
-var sqrt = _sqrt;
-var abs = _abs;
-var exp = _exp;
-var log = _log;
-var sin = _sin;
-var cos = _cos;
-var tan = _tan;
+const add = makeArithmetic(add$1, add$2);
+const subtract = makeArithmetic(subtract$1, subtract$2);
+const multiply = makeArithmetic(multiply$1, multiply$2);
+const divide = divide$1;
+const pow = pow$1;
+const sqrt = sqrt$1;
+const abs = abs$1;
+const exp = exp$1;
+const log = log$1;
+const sin = sin$1;
+const cos = cos$1;
+const tan = tan$1;
 function asin(x) {
-  return Math.asin(x);
+	return Math.asin(x);
 }
 function acos(x) {
-  return Math.acos(x);
+	return Math.acos(x);
 }
 function atan(x) {
-  return Math.atan(x);
+	return Math.atan(x);
 }
 function atan2(y, x) {
-  return Math.atan2(y, x);
+	return Math.atan2(y, x);
 }
-var sum = _sum;
-var mean = _mean;
+const sum = sum$1;
+const mean = mean$1;
 function toNumericArray(data) {
-  if (data instanceof DenseMatrix) return data.toArray().flat(Infinity);
-  if (Array.isArray(data)) return data.flat(Infinity);
-  return [data];
+	if (data instanceof DenseMatrix$1) return data.toArray().flat(Infinity);
+	if (Array.isArray(data)) return data.flat(Infinity);
+	return [data];
 }
 function variance(data, normalization = "unbiased") {
-  const arr = toNumericArray(data);
-  const n = arr.length;
-  if (n === 0) return NaN;
-  const mu = arr.reduce((s, x) => s + x, 0) / n;
-  const ss = arr.reduce((s, x) => s + (x - mu) * (x - mu), 0);
-  const denom = normalization === "uncorrected" ? n : normalization === "biased" ? n + 1 : n > 1 ? n - 1 : n;
-  return ss / denom;
+	const arr = toNumericArray(data);
+	const n = arr.length;
+	if (n === 0) return NaN;
+	const mu = arr.reduce((s, x) => s + x, 0) / n;
+	return arr.reduce((s, x) => s + (x - mu) * (x - mu), 0) / (normalization === "uncorrected" ? n : normalization === "biased" ? n + 1 : n > 1 ? n - 1 : n);
 }
 function std(data, normalization = "unbiased") {
-  return Math.sqrt(variance(data, normalization));
+	return Math.sqrt(variance(data, normalization));
 }
-var min = _min;
-var max = _max;
-var gcd = _gcd;
-var lcm = _lcm;
-var round = _round;
-var floor = _floor;
-var ceil = _ceil;
+const min = min$1;
+const max = max$1;
+const gcd = gcd$1;
+const lcm = lcm$1;
+const round = round$1;
+const floor = floor$1;
+const ceil = ceil$1;
 function conj(c) {
-  if (!isComplex(c)) {
-    throw new TypeError("Expected Complex number");
-  }
-  return c.conjugate();
+	if (!isComplex(c)) throw new TypeError("Expected Complex number");
+	return c.conjugate();
 }
 function re(c) {
-  if (typeof c === "number") return c;
-  if (isComplex(c)) return c.re;
-  throw new TypeError("Expected Complex or number");
+	if (typeof c === "number") return c;
+	if (isComplex(c)) return c.re;
+	throw new TypeError("Expected Complex or number");
 }
 function im(c) {
-  if (typeof c === "number") return 0;
-  if (isComplex(c)) return c.im;
-  throw new TypeError("Expected Complex or number");
+	if (typeof c === "number") return 0;
+	if (isComplex(c)) return c.im;
+	throw new TypeError("Expected Complex or number");
 }
 function arg(c) {
-  if (!isComplex(c)) {
-    throw new TypeError("Expected Complex number");
-  }
-  return c.arg();
+	if (!isComplex(c)) throw new TypeError("Expected Complex number");
+	return c.arg();
 }
 function transpose(m) {
-  if (Array.isArray(m)) {
-    return DenseMatrix.fromArray(m).transpose().toArray();
-  }
-  return m.transpose();
+	if (Array.isArray(m)) return DenseMatrix$1.fromArray(m).transpose().toArray();
+	return m.transpose();
 }
 function det(m) {
-  const matrix2 = Array.isArray(m) ? DenseMatrix.fromArray(m) : m;
-  if (!(matrix2 instanceof DenseMatrix)) {
-    throw new TypeError("Expected DenseMatrix or 2D array");
-  }
-  if (matrix2.rows !== matrix2.cols) {
-    throw new Error("Determinant requires a square matrix");
-  }
-  const n = matrix2.rows;
-  if (n === 0) return 1;
-  if (n === 1) return matrix2.get(0, 0);
-  if (n === 2) {
-    return matrix2.get(0, 0) * matrix2.get(1, 1) - matrix2.get(0, 1) * matrix2.get(1, 0);
-  }
-  if (n === 3) {
-    const a = matrix2.get(0, 0), b = matrix2.get(0, 1), c = matrix2.get(0, 2);
-    const d = matrix2.get(1, 0), e2 = matrix2.get(1, 1), f = matrix2.get(1, 2);
-    const g = matrix2.get(2, 0), h = matrix2.get(2, 1), i2 = matrix2.get(2, 2);
-    return a * (e2 * i2 - f * h) - b * (d * i2 - f * g) + c * (d * h - e2 * g);
-  }
-  const lu = [];
-  for (let i2 = 0; i2 < n; i2++) {
-    lu[i2] = [];
-    for (let j = 0; j < n; j++) {
-      lu[i2][j] = matrix2.get(i2, j);
-    }
-  }
-  let swaps = 0;
-  for (let k = 0; k < n - 1; k++) {
-    let maxVal = Math.abs(lu[k][k]);
-    let maxRow = k;
-    for (let i2 = k + 1; i2 < n; i2++) {
-      if (Math.abs(lu[i2][k]) > maxVal) {
-        maxVal = Math.abs(lu[i2][k]);
-        maxRow = i2;
-      }
-    }
-    if (maxVal < 1e-12) {
-      return 0;
-    }
-    if (maxRow !== k) {
-      [lu[k], lu[maxRow]] = [lu[maxRow], lu[k]];
-      swaps++;
-    }
-    for (let i2 = k + 1; i2 < n; i2++) {
-      lu[i2][k] /= lu[k][k];
-      for (let j = k + 1; j < n; j++) {
-        lu[i2][j] -= lu[i2][k] * lu[k][j];
-      }
-    }
-  }
-  let result = 1;
-  for (let i2 = 0; i2 < n; i2++) {
-    result *= lu[i2][i2];
-  }
-  return swaps % 2 === 0 ? result : -result;
+	const matrix = Array.isArray(m) ? DenseMatrix$1.fromArray(m) : m;
+	if (!(matrix instanceof DenseMatrix$1)) throw new TypeError("Expected DenseMatrix or 2D array");
+	if (matrix.rows !== matrix.cols) throw new Error("Determinant requires a square matrix");
+	const n = matrix.rows;
+	if (n === 0) return 1;
+	if (n === 1) return matrix.get(0, 0);
+	if (n === 2) return matrix.get(0, 0) * matrix.get(1, 1) - matrix.get(0, 1) * matrix.get(1, 0);
+	if (n === 3) {
+		const a = matrix.get(0, 0), b = matrix.get(0, 1), c = matrix.get(0, 2);
+		const d = matrix.get(1, 0), e = matrix.get(1, 1), f = matrix.get(1, 2);
+		const g = matrix.get(2, 0), h = matrix.get(2, 1), i = matrix.get(2, 2);
+		return a * (e * i - f * h) - b * (d * i - f * g) + c * (d * h - e * g);
+	}
+	const lu = [];
+	for (let i = 0; i < n; i++) {
+		lu[i] = [];
+		for (let j = 0; j < n; j++) lu[i][j] = matrix.get(i, j);
+	}
+	let swaps = 0;
+	for (let k = 0; k < n - 1; k++) {
+		let maxVal = Math.abs(lu[k][k]);
+		let maxRow = k;
+		for (let i = k + 1; i < n; i++) if (Math.abs(lu[i][k]) > maxVal) {
+			maxVal = Math.abs(lu[i][k]);
+			maxRow = i;
+		}
+		if (maxVal < 1e-12) return 0;
+		if (maxRow !== k) {
+			[lu[k], lu[maxRow]] = [lu[maxRow], lu[k]];
+			swaps++;
+		}
+		for (let i = k + 1; i < n; i++) {
+			lu[i][k] /= lu[k][k];
+			for (let j = k + 1; j < n; j++) lu[i][j] -= lu[i][k] * lu[k][j];
+		}
+	}
+	let result = 1;
+	for (let i = 0; i < n; i++) result *= lu[i][i];
+	return swaps % 2 === 0 ? result : -result;
 }
 function identity(n) {
-  return DenseMatrix.identity(n);
+	return DenseMatrix$1.identity(n);
 }
 function zeros(rows, cols) {
-  if (cols === void 0) {
-    return new Array(rows).fill(0);
-  }
-  return DenseMatrix.zeros(rows, cols);
+	if (cols === void 0) return new Array(rows).fill(0);
+	return DenseMatrix$1.zeros(rows, cols);
 }
 function ones(rows, cols) {
-  if (cols === void 0) {
-    return new Array(rows).fill(1);
-  }
-  return DenseMatrix.ones(rows, cols);
+	if (cols === void 0) return new Array(rows).fill(1);
+	return DenseMatrix$1.ones(rows, cols);
 }
 function size(m) {
-  if (Array.isArray(m)) {
-    return [m.length, m[0]?.length ?? 0];
-  }
-  return [m.rows, m.cols];
+	if (Array.isArray(m)) return [m.length, m[0]?.length ?? 0];
+	return [m.rows, m.cols];
 }
 function isComplex_(x) {
-  return isComplex(x);
+	return isComplex(x);
 }
 function isFraction_(x) {
-  return isFraction(x);
+	return isFraction(x);
 }
 function isBigNumber_(x) {
-  return isBigNumber(x);
+	return isBigNumber(x);
 }
 function isNumber_(x) {
-  return isNumber(x);
+	return isNumber(x);
 }
 function isMatrix(x) {
-  return x instanceof DenseMatrix || x instanceof SparseMatrix;
+	return x instanceof DenseMatrix$1 || x instanceof SparseMatrix$1;
 }
-var i = I;
-var pi = Math.PI;
-var e = Math.E;
-var phi = (1 + Math.sqrt(5)) / 2;
-var tau = 2 * Math.PI;
-var LN2 = CORE_LN2;
-var LN10 = CORE_LN10;
-var LOG2E = CORE_LOG2E;
-var LOG10E = CORE_LOG10E;
-var SQRT2 = CORE_SQRT2;
-var SQRT1_2 = CORE_SQRT1_2;
-var Infinity_ = Infinity;
-var NaN_ = NaN;
-var shims = {
-  // Type creation
-  complex,
-  fraction,
-  bignumber,
-  matrix,
-  sparse,
-  // Arithmetic
-  add,
-  subtract,
-  multiply,
-  divide,
-  pow,
-  sqrt,
-  abs,
-  exp,
-  log,
-  // Trigonometry
-  sin,
-  cos,
-  tan,
-  asin,
-  acos,
-  atan,
-  atan2,
-  // Statistics
-  sum,
-  mean,
-  std,
-  variance,
-  min,
-  max,
-  // Number theory
-  gcd,
-  lcm,
-  // Rounding
-  round,
-  floor,
-  ceil,
-  // Complex-specific
-  conj,
-  re,
-  im,
-  arg,
-  // Matrix-specific
-  transpose,
-  det,
-  identity,
-  zeros,
-  ones,
-  size,
-  // Type checking
-  isComplex: isComplex_,
-  isFraction: isFraction_,
-  isBigNumber: isBigNumber_,
-  isNumber: isNumber_,
-  isMatrix,
-  // Constants
-  i,
-  pi,
-  e,
-  phi,
-  tau,
-  LN2,
-  LN10,
-  LOG2E,
-  LOG10E,
-  SQRT2,
-  SQRT1_2,
-  Infinity: Infinity_,
-  NaN: NaN_
+const i = I$1;
+const pi = Math.PI;
+const e = Math.E;
+const phi = (1 + Math.sqrt(5)) / 2;
+const tau = 2 * Math.PI;
+const LN2 = LN2$1;
+const LN10 = LN10$1;
+const LOG2E = LOG2E$1;
+const LOG10E = LOG10E$1;
+const SQRT2 = SQRT2$1;
+const SQRT1_2 = SQRT1_2$1;
+const Infinity_ = Infinity;
+const NaN_ = NaN;
+const shims = {
+	complex,
+	fraction,
+	bignumber,
+	matrix,
+	sparse,
+	add,
+	subtract,
+	multiply,
+	divide,
+	pow,
+	sqrt,
+	abs,
+	exp,
+	log,
+	sin,
+	cos,
+	tan,
+	asin,
+	acos,
+	atan,
+	atan2,
+	sum,
+	mean,
+	std,
+	variance,
+	min,
+	max,
+	gcd,
+	lcm,
+	round,
+	floor,
+	ceil,
+	conj,
+	re,
+	im,
+	arg,
+	transpose,
+	det,
+	identity,
+	zeros,
+	ones,
+	size,
+	isComplex: isComplex_,
+	isFraction: isFraction_,
+	isBigNumber: isBigNumber_,
+	isNumber: isNumber_,
+	isMatrix,
+	i,
+	pi,
+	e,
+	phi,
+	tau,
+	LN2,
+	LN10,
+	LOG2E,
+	LOG10E,
+	SQRT2,
+	SQRT1_2,
+	Infinity: Infinity_,
+	NaN: NaN_
 };
-
-// src/chain.ts
+//#endregion
+//#region src/chain.ts
 function createChain(math) {
-  function chain(value) {
-    const base = {
-      done: () => value,
-      valueOf: () => value,
-      toString: () => String(value)
-    };
-    return new Proxy(base, {
-      get(target, prop) {
-        if (prop in target || typeof prop === "symbol") {
-          return target[prop];
-        }
-        const fn = math[prop];
-        if (typeof fn === "function") {
-          return (...args) => chain(fn(value, ...args));
-        }
-        return void 0;
-      }
-    });
-  }
-  return chain;
+	function chain(value) {
+		return new Proxy({
+			done: () => value,
+			valueOf: () => value,
+			toString: () => String(value)
+		}, { get(target, prop) {
+			if (prop in target || typeof prop === "symbol") return target[prop];
+			const fn = math[prop];
+			if (typeof fn === "function") return (...args) => chain(fn(value, ...args));
+		} });
+	}
+	return chain;
 }
-
-// src/index.ts
-import * as mathFunctions from "@danielsimonjr/mathts-functions";
-import {
-  Complex as Complex2,
-  Fraction as Fraction2,
-  BigNumber as BigNumber2,
-  I as I2,
-  COMPLEX_ZERO as COMPLEX_ZERO2,
-  COMPLEX_ONE,
-  FRACTION_ZERO,
-  FRACTION_ONE,
-  BIGNUMBER_ZERO,
-  BIGNUMBER_ONE,
-  BIGNUMBER_PI,
-  BIGNUMBER_E
-} from "@danielsimonjr/mathts-core";
-import { DenseMatrix as DenseMatrix2, SparseMatrix as SparseMatrix2 } from "@danielsimonjr/mathts-matrix";
-import { computePool } from "@danielsimonjr/mathts-parallel";
-var defaultConfig = {
-  precision: 64,
-  matrix: "Matrix",
-  number: "number",
-  epsilon: 1e-12,
-  randomSeed: null
+//#endregion
+//#region src/index.ts
+const defaultConfig = {
+	precision: 64,
+	matrix: "Matrix",
+	number: "number",
+	epsilon: 1e-12,
+	randomSeed: null
 };
-function create(factories = all, config2 = {}) {
-  const currentConfig = { ...defaultConfig, ...config2 };
-  if (Object.keys(config2).length > 0) {
-    mathFunctions.config(config2);
-  }
-  const instance = {
-    // Full @danielsimonjr/mathts-functions namespace — det, integrate, eigs,
-    // simplify, and the rest. The compat-specific shims below overlay it.
-    ...factories,
-    // Configuration — forwards to the shared functions runtime config so it
-    // actually DRIVES behavior (e.g. `config({ matrix: 'Array' })` changes the
-    // return type of matrix-producing ops). `currentConfig` mirrors it for the
-    // read path and seeds the functions config with this instance's overrides.
-    config: (newConfig) => {
-      const merged = mathFunctions.config(newConfig);
-      Object.assign(currentConfig, merged);
-      return { ...currentConfig };
-    },
-    // Type creation
-    complex: shims.complex,
-    fraction: shims.fraction,
-    bignumber: shims.bignumber,
-    matrix: shims.matrix,
-    sparse: shims.sparse,
-    // Arithmetic
-    add: shims.add,
-    subtract: shims.subtract,
-    multiply: shims.multiply,
-    divide: shims.divide,
-    pow: shims.pow,
-    sqrt: shims.sqrt,
-    abs: shims.abs,
-    exp: shims.exp,
-    log: shims.log,
-    // Trigonometry
-    sin: shims.sin,
-    cos: shims.cos,
-    tan: shims.tan,
-    asin: shims.asin,
-    acos: shims.acos,
-    atan: shims.atan,
-    atan2: shims.atan2,
-    // Statistics
-    sum: shims.sum,
-    mean: shims.mean,
-    std: shims.std,
-    variance: shims.variance,
-    min: shims.min,
-    max: shims.max,
-    // Number theory
-    gcd: shims.gcd,
-    lcm: shims.lcm,
-    // Rounding
-    round: shims.round,
-    floor: shims.floor,
-    ceil: shims.ceil,
-    // Complex-specific
-    conj: shims.conj,
-    re: shims.re,
-    im: shims.im,
-    arg: shims.arg,
-    // Matrix-specific
-    transpose: shims.transpose,
-    det: shims.det,
-    identity: shims.identity,
-    zeros: shims.zeros,
-    ones: shims.ones,
-    size: shims.size,
-    // Type checking
-    isComplex: shims.isComplex,
-    isFraction: shims.isFraction,
-    isBigNumber: shims.isBigNumber,
-    isNumber: shims.isNumber,
-    isMatrix: shims.isMatrix,
-    // Constants
-    i: shims.i,
-    pi: shims.pi,
-    e: shims.e,
-    phi: shims.phi,
-    tau: shims.tau,
-    LN2: shims.LN2,
-    LN10: shims.LN10,
-    LOG2E: shims.LOG2E,
-    LOG10E: shims.LOG10E,
-    SQRT2: shims.SQRT2,
-    SQRT1_2: shims.SQRT1_2,
-    Infinity: shims.Infinity,
-    NaN: shims.NaN,
-    // Placeholder; rebound below once `instance` exists so chain methods can
-    // resolve against the full surfaced namespace.
-    chain: (() => {
-      throw new Error("chain not initialized");
-    })
-  };
-  instance.chain = createChain(instance);
-  return instance;
+function create(factories = all, config = {}) {
+	const currentConfig = {
+		...defaultConfig,
+		...config
+	};
+	if (Object.keys(config).length > 0) mathFunctions.config(config);
+	const instance = {
+		...factories,
+		config: (newConfig) => {
+			const merged = mathFunctions.config(newConfig);
+			Object.assign(currentConfig, merged);
+			return { ...currentConfig };
+		},
+		complex: shims.complex,
+		fraction: shims.fraction,
+		bignumber: shims.bignumber,
+		matrix: shims.matrix,
+		sparse: shims.sparse,
+		add: shims.add,
+		subtract: shims.subtract,
+		multiply: shims.multiply,
+		divide: shims.divide,
+		pow: shims.pow,
+		sqrt: shims.sqrt,
+		abs: shims.abs,
+		exp: shims.exp,
+		log: shims.log,
+		sin: shims.sin,
+		cos: shims.cos,
+		tan: shims.tan,
+		asin: shims.asin,
+		acos: shims.acos,
+		atan: shims.atan,
+		atan2: shims.atan2,
+		sum: shims.sum,
+		mean: shims.mean,
+		std: shims.std,
+		variance: shims.variance,
+		min: shims.min,
+		max: shims.max,
+		gcd: shims.gcd,
+		lcm: shims.lcm,
+		round: shims.round,
+		floor: shims.floor,
+		ceil: shims.ceil,
+		conj: shims.conj,
+		re: shims.re,
+		im: shims.im,
+		arg: shims.arg,
+		transpose: shims.transpose,
+		det: shims.det,
+		identity: shims.identity,
+		zeros: shims.zeros,
+		ones: shims.ones,
+		size: shims.size,
+		isComplex: shims.isComplex,
+		isFraction: shims.isFraction,
+		isBigNumber: shims.isBigNumber,
+		isNumber: shims.isNumber,
+		isMatrix: shims.isMatrix,
+		i: shims.i,
+		pi: shims.pi,
+		e: shims.e,
+		phi: shims.phi,
+		tau: shims.tau,
+		LN2: shims.LN2,
+		LN10: shims.LN10,
+		LOG2E: shims.LOG2E,
+		LOG10E: shims.LOG10E,
+		SQRT2: shims.SQRT2,
+		SQRT1_2: shims.SQRT1_2,
+		Infinity: shims.Infinity,
+		NaN: shims.NaN,
+		chain: (() => {
+			throw new Error("chain not initialized");
+		})
+	};
+	instance.chain = createChain(instance);
+	return instance;
 }
-var all = {
-  ...mathFunctions
-};
-export {
-  BIGNUMBER_E,
-  BIGNUMBER_ONE,
-  BIGNUMBER_PI,
-  BIGNUMBER_ZERO,
-  BigNumber2 as BigNumber,
-  COMPLEX_ONE,
-  COMPLEX_ZERO2 as COMPLEX_ZERO,
-  Complex2 as Complex,
-  DenseMatrix2 as DenseMatrix,
-  FRACTION_ONE,
-  FRACTION_ZERO,
-  Fraction2 as Fraction,
-  I2 as I,
-  Infinity_,
-  LN10,
-  LN2,
-  LOG10E,
-  LOG2E,
-  NaN_,
-  SQRT1_2,
-  SQRT2,
-  SparseMatrix2 as SparseMatrix,
-  abs,
-  acos,
-  add,
-  all,
-  arg,
-  asin,
-  atan,
-  atan2,
-  bignumber,
-  ceil,
-  complex,
-  computePool,
-  conj,
-  cos,
-  create,
-  det,
-  divide,
-  e,
-  exp,
-  floor,
-  fraction,
-  gcd,
-  i,
-  identity,
-  im,
-  isBigNumber_,
-  isComplex_,
-  isFraction_,
-  isMatrix,
-  isNumber_,
-  lcm,
-  log,
-  matrix,
-  max,
-  mean,
-  min,
-  multiply,
-  ones,
-  phi,
-  pi,
-  pow,
-  re,
-  round,
-  shims,
-  sin,
-  size,
-  sparse,
-  sqrt,
-  std,
-  subtract,
-  sum,
-  tan,
-  tau,
-  transpose,
-  variance,
-  zeros
-};
+const all = { ...mathFunctions };
+//#endregion
+export { BIGNUMBER_E, BIGNUMBER_ONE, BIGNUMBER_PI, BIGNUMBER_ZERO, BigNumber, COMPLEX_ONE, COMPLEX_ZERO, Complex, DenseMatrix, FRACTION_ONE, FRACTION_ZERO, Fraction, I, Infinity_, LN10, LN2, LOG10E, LOG2E, NaN_, SQRT1_2, SQRT2, SparseMatrix, abs, acos, add, all, arg, asin, atan, atan2, bignumber, ceil, complex, computePool, conj, cos, create, det, divide, e, exp, floor, fraction, gcd, i, identity, im, isBigNumber_, isComplex_, isFraction_, isMatrix, isNumber_, lcm, log, matrix, max, mean, min, multiply, ones, phi, pi, pow, re, round, shims, sin, size, sparse, sqrt, std, subtract, sum, tan, tau, transpose, variance, zeros };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-compat",
-  "version": "0.4.29",
+  "version": "0.4.31",
   "description": "mathjs compatibility layer for MathTS",
   "type": "module",
   "main": "./dist/index.js",
@@ -16,22 +16,33 @@
     "README.md"
   ],
   "scripts": {
-    "build": "tsup src/index.ts --format esm --clean && tsc -p tsconfig.dts.json",
-    "dev": "tsup src/index.ts --format esm --watch --onSuccess \"tsc -p tsconfig.dts.json\"",
+    "build": "tsdown src/index.ts --format esm && tsc -p tsconfig.dts.json",
+    "dev": "tsdown src/index.ts --format esm --watch --on-success \"tsc -p tsconfig.dts.json\"",
     "typecheck": "tsc --noEmit && tsc -p tsconfig.test.json",
     "test": "bun test --isolate --timeout 30000",
     "test:coverage": "bun test --isolate --timeout 30000 --coverage",
-    "build:prod": "tsup src/index.ts --format esm --clean --minify --treeshake && tsc -p tsconfig.dts.json"
+    "build:prod": "tsdown src/index.ts --format esm --minify && tsc -p tsconfig.dts.json"
+  },
+  "tsdown": {
+    "fixedExtension": false,
+    "dts": false,
+    "outputOptions": {
+      "comments": {
+        "legal": true,
+        "annotation": true,
+        "jsdoc": false
+      }
+    }
   },
   "dependencies": {
-    "@danielsimonjr/mathts-core": "^0.17.0",
-    "@danielsimonjr/mathts-functions": "^0.68.1",
-    "@danielsimonjr/mathts-matrix": "^0.7.7",
-    "@danielsimonjr/mathts-parallel": "^0.6.9"
+    "@danielsimonjr/mathts-core": "^0.18.0",
+    "@danielsimonjr/mathts-functions": "^0.69.0",
+    "@danielsimonjr/mathts-matrix": "^0.7.8",
+    "@danielsimonjr/mathts-parallel": "^0.6.10"
   },
   "devDependencies": {
     "@types/node": "^26.4.0",
-    "tsup": "^8.0.0",
+    "tsdown": "^0.23.0",
     "typescript": "^7.0.2",
     "vitest": "^5.0.0"
   },