solve-engine 1.0.1 → 1.1.0

package/dist/{chunk-O3ANBHSA.js → chunk-NNQ2TYDF.js}

@@ -1,5 +1,123 @@
-import { boolValue, symbolicValue, numberValue, matrixValue, errorValue, checkedArray, faultedOperand } from './chunk-GQCOSXMG.js';
-import { ErrorFactory } from './chunk-3D7V24DG.js';
+import { boolValue, symbolicValue, numberValue, matrixValue, errorValue, checkedArray, faultedOperand } from './chunk-267JPOTF.js';
+import { ErrorFactory } from './chunk-ERCOHGXD.js';
+
+// src/decimal/Decimal.ts
+var DEFAULT_DIVISION_SCALE = 20;
+var POW10 = [1n];
+function pow10(n) {
+  if (n < 0) {
+    throw ErrorFactory.internal("DECIMAL_NEGATIVE_POWER", `pow10 needs a non-negative exponent, got ${n}`, { exponent: n });
+  }
+  for (let i = POW10.length; i <= n; i++) POW10[i] = POW10[i - 1] * 10n;
+  return POW10[n];
+}
+function decimalFromInteger(n) {
+  return { coef: typeof n === "bigint" ? n : BigInt(n), scale: 0 };
+}
+function decimalFromNumberIfExact(n) {
+  if (!Number.isFinite(n) || !Number.isInteger(n)) return null;
+  return decimalFromInteger(BigInt(n));
+}
+function decimalFromLiteral(text) {
+  let s = text;
+  let negative = false;
+  if (s.startsWith("-")) {
+    negative = true;
+    s = s.slice(1);
+  } else if (s.startsWith("+")) {
+    s = s.slice(1);
+  }
+  const dot = s.indexOf(".");
+  const intPart = dot < 0 ? s : s.slice(0, dot);
+  const fracPart = dot < 0 ? "" : s.slice(dot + 1);
+  const digits = intPart + fracPart || "0";
+  if (!/^[0-9]+$/.test(digits)) {
+    throw ErrorFactory.parsing("INVALID_DECIMAL_LITERAL", `Not a decimal literal: "${text}"`, { raw: text });
+  }
+  let coef = BigInt(digits);
+  if (negative) coef = -coef;
+  return { coef, scale: fracPart.length };
+}
+function align(a, b) {
+  if (a.scale === b.scale) return { ca: a.coef, cb: b.coef, scale: a.scale };
+  if (a.scale > b.scale) return { ca: a.coef, cb: b.coef * pow10(a.scale - b.scale), scale: a.scale };
+  return { ca: a.coef * pow10(b.scale - a.scale), cb: b.coef, scale: b.scale };
+}
+function decimalAdd(a, b) {
+  const { ca, cb, scale } = align(a, b);
+  return { coef: ca + cb, scale };
+}
+function decimalSubtract(a, b) {
+  const { ca, cb, scale } = align(a, b);
+  return { coef: ca - cb, scale };
+}
+function decimalMultiply(a, b) {
+  return { coef: a.coef * b.coef, scale: a.scale + b.scale };
+}
+function decimalNegate(a) {
+  return { coef: -a.coef, scale: a.scale };
+}
+function decimalIsZero(a) {
+  return a.coef === 0n;
+}
+function roundCoefToScale(coef, scale, targetScale) {
+  if (targetScale >= scale) return coef * pow10(targetScale - scale);
+  const divisor = pow10(scale - targetScale);
+  const negative = coef < 0n;
+  const magnitude = negative ? -coef : coef;
+  let q = magnitude / divisor;
+  const r = magnitude % divisor;
+  if (r * 2n >= divisor) q += 1n;
+  return negative ? -q : q;
+}
+function decimalRound(a, targetScale) {
+  return { coef: roundCoefToScale(a.coef, a.scale, targetScale), scale: targetScale };
+}
+function decimalDivide(a, b, maxFractionDigits = DEFAULT_DIVISION_SCALE) {
+  if (decimalIsZero(b)) {
+    throw ErrorFactory.internal("DECIMAL_DIVISION_BY_ZERO", "decimalDivide was called with a zero divisor", {});
+  }
+  const s = maxFractionDigits;
+  const numerator = a.coef * pow10(s + b.scale);
+  const denominator = b.coef * pow10(a.scale);
+  const negative = numerator < 0n !== denominator < 0n;
+  const num = numerator < 0n ? -numerator : numerator;
+  const den = denominator < 0n ? -denominator : denominator;
+  let q = num / den;
+  const r = num % den;
+  if (r * 2n >= den) q += 1n;
+  return { coef: negative ? -q : q, scale: s };
+}
+function decimalCompare(a, b) {
+  const { ca, cb } = align(a, b);
+  if (ca < cb) return -1;
+  if (ca > cb) return 1;
+  return 0;
+}
+function decimalToString(a) {
+  const negative = a.coef < 0n;
+  let digits = (negative ? -a.coef : a.coef).toString();
+  if (a.scale === 0) return negative ? `-${digits}` : digits;
+  if (digits.length <= a.scale) digits = digits.padStart(a.scale + 1, "0");
+  const cut = digits.length - a.scale;
+  const whole = digits.slice(0, cut);
+  const frac = digits.slice(cut);
+  return `${negative ? "-" : ""}${whole}.${frac}`;
+}
+function decimalToNumber(a) {
+  return Number(decimalToString(a));
+}
+function decimalToFixed(a, fractionDigits) {
+  const rounded = decimalRound(a, fractionDigits);
+  const negative = rounded.coef < 0n && rounded.coef !== 0n;
+  let digits = (rounded.coef < 0n ? -rounded.coef : rounded.coef).toString();
+  if (fractionDigits === 0) return negative ? `-${digits}` : digits;
+  if (digits.length <= fractionDigits) digits = digits.padStart(fractionDigits + 1, "0");
+  const cut = digits.length - fractionDigits;
+  const whole = digits.slice(0, cut);
+  const frac = digits.slice(cut);
+  return `${negative ? "-" : ""}${whole}.${frac}`;
+}
 
 // src/symbolic/Rational.ts
 var RATIONAL_ZERO = { n: 0n, d: 1n };
@@ -1725,6 +1843,6 @@ var CURRENCY_DISPLAY = {
   CZK: { symbol: "K\u010D", position: "suffix", spaced: true }
 };
 
-export { COMPLEX_I, CURRENCY_DISPLAY, RATIONAL_ONE, RATIONAL_ZERO, callNode, cancelSymbolic, coefficientGcd, collectionToValues, columnMajorToRowMajor, compareMonomials, complex, complexNode, constNode, determinant, exactIntegerCbrt, exactIntegerSqrt, exactRationalSqrt2 as exactRationalSqrt, expandSymbolic, formatSymbolic, freeVariables, fromPolynomial, inBounds, inverse, isRationalOne, isRationalZero, isZeroPolynomial, matAt, matIndex, matrixCompare, matrixEntryToValue, matrixMultiply, matrixPower, monomialFromKey, monomialToKey, polyCoefficients, polyDegree, polyUnivariateVar, polynomialDivide, powNode, rationalAdd, rationalCompare, rationalDiv, rationalFromNumber, rationalMul, rationalNeg, rationalSub, rationalToNumber, resolveCurrencyAlias, rowMajorToColumnMajor, sameShape, signedSurdNode, simplifySymbolic, squareFreePart, substitute, surdNode, symbolicKey, symbolicToEntry, toPolynomial, transpose, varNode };
-//# sourceMappingURL=chunk-O3ANBHSA.js.map
-//# sourceMappingURL=chunk-O3ANBHSA.js.map
+export { COMPLEX_I, CURRENCY_DISPLAY, RATIONAL_ONE, RATIONAL_ZERO, callNode, cancelSymbolic, coefficientGcd, collectionToValues, columnMajorToRowMajor, compareMonomials, complex, complexNode, constNode, decimalAdd, decimalCompare, decimalDivide, decimalFromInteger, decimalFromLiteral, decimalFromNumberIfExact, decimalIsZero, decimalMultiply, decimalNegate, decimalRound, decimalSubtract, decimalToFixed, decimalToNumber, determinant, exactIntegerCbrt, exactIntegerSqrt, exactRationalSqrt2 as exactRationalSqrt, expandSymbolic, formatSymbolic, freeVariables, fromPolynomial, inBounds, inverse, isRationalOne, isRationalZero, isZeroPolynomial, matAt, matIndex, matrixCompare, matrixEntryToValue, matrixMultiply, matrixPower, monomialFromKey, monomialToKey, polyCoefficients, polyDegree, polyUnivariateVar, polynomialDivide, powNode, rational, rationalAdd, rationalCompare, rationalDiv, rationalFromNumber, rationalMul, rationalNeg, rationalSub, rationalToNumber, resolveCurrencyAlias, rowMajorToColumnMajor, sameShape, signedSurdNode, simplifySymbolic, squareFreePart, substitute, surdNode, symbolicKey, symbolicToEntry, toPolynomial, transpose, varNode };
+//# sourceMappingURL=chunk-NNQ2TYDF.js.map
+//# sourceMappingURL=chunk-NNQ2TYDF.js.map