solve-engine 1.1.1 → 1.2.0

package/dist/{chunk-HIQ5HSZL.cjs → chunk-A2DZ5BEJ.cjs}

@@ -1,18 +1,24 @@
 'use strict';
 
-var chunk2TZKENDH_cjs = require('./chunk-2TZKENDH.cjs');
-var chunkB4HBEFTB_cjs = require('./chunk-B4HBEFTB.cjs');
+var chunkIT74O4BE_cjs = require('./chunk-IT74O4BE.cjs');
+var chunkMIGXCBXT_cjs = require('./chunk-MIGXCBXT.cjs');
 
 // src/decimal/Decimal.ts
 var DEFAULT_DIVISION_SCALE = 20;
 var POW10 = [1n];
 function pow10(n) {
   if (n < 0) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.internal("DECIMAL_NEGATIVE_POWER", `pow10 needs a non-negative exponent, got ${n}`, { exponent: n });
+    throw chunkMIGXCBXT_cjs.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 makeDecimal(coef, scale) {
+  if (!Number.isInteger(scale) || scale < 0) {
+    throw chunkMIGXCBXT_cjs.ErrorFactory.internal("DECIMAL_INVALID_SCALE", `A decimal scale must be a non-negative integer, got ${scale}`, { scale });
+  }
+  return { coef, scale };
+}
 function decimalFromInteger(n) {
   return { coef: typeof n === "bigint" ? n : BigInt(n), scale: 0 };
 }
@@ -34,7 +40,7 @@ function decimalFromLiteral(text) {
   const fracPart = dot < 0 ? "" : s.slice(dot + 1);
   const digits = intPart + fracPart || "0";
   if (!/^[0-9]+$/.test(digits)) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.parsing("INVALID_DECIMAL_LITERAL", `Not a decimal literal: "${text}"`, { raw: text });
+    throw chunkMIGXCBXT_cjs.ErrorFactory.parsing("INVALID_DECIMAL_LITERAL", `Not a decimal literal: "${text}"`, { raw: text });
   }
   let coef = BigInt(digits);
   if (negative) coef = -coef;
@@ -77,7 +83,7 @@ function decimalRound(a, targetScale) {
 }
 function decimalDivide(a, b, maxFractionDigits = DEFAULT_DIVISION_SCALE) {
   if (decimalIsZero(b)) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.internal("DECIMAL_DIVISION_BY_ZERO", "decimalDivide was called with a zero divisor", {});
+    throw chunkMIGXCBXT_cjs.ErrorFactory.internal("DECIMAL_DIVISION_BY_ZERO", "decimalDivide was called with a zero divisor", {});
   }
   const s = maxFractionDigits;
   const numerator = a.coef * pow10(s + b.scale);
@@ -151,7 +157,7 @@ function gcd(a, b) {
 function guardMagnitude(n, d) {
   const magnitude = n < 0n ? -n : n;
   if (magnitude >= MAX_MAGNITUDE || d >= MAX_MAGNITUDE) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution(
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution(
       "SYMBOLIC_RATIONAL_OVERFLOW",
       `This calculation produced a fraction too large to represent exactly (over ${RATIONAL_MAX_BITS} bits).`,
       { limitBits: RATIONAL_MAX_BITS }
@@ -160,7 +166,7 @@ function guardMagnitude(n, d) {
 }
 function normalize(n, d) {
   if (d === 0n) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a symbolic expression.");
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a symbolic expression.");
   }
   if (n === 0n) return RATIONAL_ZERO;
   let numerator = d < 0n ? -n : n;
@@ -177,7 +183,7 @@ function rational(numerator, denominator = 1n) {
 var DECIMAL_FORM = /^(-?)(\d+)(?:\.(\d+))?(?:[eE]([+-]?\d+))?$/;
 function rationalFromNumber(value) {
   if (!Number.isFinite(value)) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution(
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution(
       "SYMBOLIC_NONFINITE_OPERAND",
       `"${String(value)}" has no exact value, so it cannot appear in a symbolic expression.`,
       { value: String(value) }
@@ -186,7 +192,7 @@ function rationalFromNumber(value) {
   if (Number.isInteger(value)) return normalize(BigInt(value), 1n);
   const parts = DECIMAL_FORM.exec(String(value));
   if (parts === null) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.internal(
+    throw chunkMIGXCBXT_cjs.ErrorFactory.internal(
       "INTERNAL_RATIONAL_PARSE",
       `Internal error: could not read "${String(value)}" as a decimal number.`,
       { value: String(value) }
@@ -218,7 +224,7 @@ function rationalMul(a, b) {
 }
 function rationalDiv(a, b) {
   if (b.n === 0n) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a symbolic expression.");
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a symbolic expression.");
   }
   return normalize(a.n * b.d, a.d * b.n);
 }
@@ -229,7 +235,7 @@ function rationalPow(base, exponent) {
   if (exponent === 0n) return RATIONAL_ONE;
   if (exponent < 0n) {
     if (base.n === 0n) {
-      throw chunkB4HBEFTB_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Zero raised to a negative power is undefined.");
+      throw chunkMIGXCBXT_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Zero raised to a negative power is undefined.");
     }
     const positive = -exponent;
     return normalize(bigintPow(base.d, positive), bigintPow(base.n, positive));
@@ -309,7 +315,7 @@ function complexMul(a, b) {
 }
 function complexDiv(a, b) {
   if (isComplexZero(b)) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a complex expression.");
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Division by zero in a complex expression.");
   }
   const denominator = rationalAdd(rationalMul(b.re, b.re), rationalMul(b.im, b.im));
   return {
@@ -330,7 +336,7 @@ function complexPow(base, exponent) {
   if (exponent === 0n) return COMPLEX_ONE;
   if (exponent < 0n) {
     if (isComplexZero(base)) {
-      throw chunkB4HBEFTB_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Zero raised to a negative power is undefined.");
+      throw chunkMIGXCBXT_cjs.ErrorFactory.execution("SYMBOLIC_DIVISION_BY_ZERO", "Zero raised to a negative power is undefined.");
     }
     return complexDiv(COMPLEX_ONE, complexPow(base, -exponent));
   }
@@ -1043,7 +1049,7 @@ function foldCall(name, args) {
 }
 function simplifySymbolic(node) {
   if (nodeCount(node, SYMBOLIC_MAX_NODES) > SYMBOLIC_MAX_NODES) {
-    throw chunkB4HBEFTB_cjs.ErrorFactory.execution(
+    throw chunkMIGXCBXT_cjs.ErrorFactory.execution(
       "SYMBOLIC_NODE_LIMIT_EXCEEDED",
       `This symbolic expression is too large to simplify (over ${SYMBOLIC_MAX_NODES} terms).`,
       { limit: SYMBOLIC_MAX_NODES }
@@ -1395,9 +1401,9 @@ function symbolicToEntry(node) {
   return simplified.kind === "const" ? rationalToNumber(simplified.value) : simplified;
 }
 function matrixEntryToValue(cell) {
-  if (typeof cell === "boolean") return chunk2TZKENDH_cjs.boolValue(cell);
-  if (typeof cell === "object" && cell !== null) return chunk2TZKENDH_cjs.symbolicValue(cell);
-  return chunk2TZKENDH_cjs.numberValue(cell);
+  if (typeof cell === "boolean") return chunkIT74O4BE_cjs.boolValue(cell);
+  if (typeof cell === "object" && cell !== null) return chunkIT74O4BE_cjs.symbolicValue(cell);
+  return chunkIT74O4BE_cjs.numberValue(cell);
 }
 function rowMajorToColumnMajor(rows, cols, rowMajorData) {
   const data = new Array(rows * cols);
@@ -1426,21 +1432,21 @@ function matrixMultiply(l, r) {
     if (useSymbolic) {
       const scalarNode = entryToSymbolic(scalarOperand.data[0]);
       const result2 = matrixOperand.data.map((v) => symbolicToEntry(simplifySymbolic({ kind: "mul", left: entryToSymbolic(v), right: scalarNode })));
-      return chunk2TZKENDH_cjs.matrixValue(matrixOperand.rows, matrixOperand.cols, result2);
+      return chunkIT74O4BE_cjs.matrixValue(matrixOperand.rows, matrixOperand.cols, result2);
     }
     const scalar = scalarOperand.data[0];
     const result = matrixOperand.data.map((v) => v * scalar);
-    return chunk2TZKENDH_cjs.matrixValue(matrixOperand.rows, matrixOperand.cols, result);
+    return chunkIT74O4BE_cjs.matrixValue(matrixOperand.rows, matrixOperand.cols, result);
   }
   if (l.cols !== r.rows) {
-    return chunk2TZKENDH_cjs.errorValue(
+    return chunkIT74O4BE_cjs.errorValue(
       "DIMENSION_MISMATCH",
       `Cannot multiply a ${l.rows}x${l.cols} matrix by a ${r.rows}x${r.cols} matrix \u2014 inner dimensions must match (${l.cols} !== ${r.rows}).`
     );
   }
   const resultRows = l.rows;
   const resultCols = r.cols;
-  const data = chunk2TZKENDH_cjs.checkedArray(resultRows * resultCols, "matrix cells");
+  const data = chunkIT74O4BE_cjs.checkedArray(resultRows * resultCols, "matrix cells");
   for (let row = 0; row < resultRows; row++) {
     for (let col = 0; col < resultCols; col++) {
       if (useSymbolic) {
@@ -1459,14 +1465,14 @@ function matrixMultiply(l, r) {
       }
     }
   }
-  return chunk2TZKENDH_cjs.matrixValue(resultRows, resultCols, data);
+  return chunkIT74O4BE_cjs.matrixValue(resultRows, resultCols, data);
 }
 function matrixPower(m, exponent) {
   if (!isSquare(m)) {
-    return chunk2TZKENDH_cjs.errorValue("MATRIX_POWER_REQUIRES_SQUARE_MATRIX", `^: only a square matrix can be raised to a power (got ${m.rows}x${m.cols}).`);
+    return chunkIT74O4BE_cjs.errorValue("MATRIX_POWER_REQUIRES_SQUARE_MATRIX", `^: only a square matrix can be raised to a power (got ${m.rows}x${m.cols}).`);
   }
   if (!Number.isInteger(exponent) || exponent < 0) {
-    return chunk2TZKENDH_cjs.errorValue(
+    return chunkIT74O4BE_cjs.errorValue(
       "MATRIX_POWER_REQUIRES_WHOLE_EXPONENT",
       `^: a matrix can only be raised to a whole number of 0 or more (got ${exponent}); "^-1" means the inverse and "^T" the transpose.`
     );
@@ -1490,18 +1496,18 @@ function matrixPower(m, exponent) {
       squared = doubled.value;
     }
   }
-  return chunk2TZKENDH_cjs.matrixValue(accumulated.rows, accumulated.cols, accumulated.data);
+  return chunkIT74O4BE_cjs.matrixValue(accumulated.rows, accumulated.cols, accumulated.data);
 }
 function matrixCompare(l, r, cmp) {
   if (!sameShape(l, r)) {
-    return chunk2TZKENDH_cjs.errorValue(
+    return chunkIT74O4BE_cjs.errorValue(
       "DIMENSION_MISMATCH",
       `Cannot compare matrices of different shapes: ${l.rows}x${l.cols} and ${r.rows}x${r.cols}`
     );
   }
   const data = new Array(l.data.length);
   for (let i = 0; i < l.data.length; i++) data[i] = cmp(l.data[i], r.data[i]);
-  return chunk2TZKENDH_cjs.matrixValue(l.rows, l.cols, data);
+  return chunkIT74O4BE_cjs.matrixValue(l.rows, l.cols, data);
 }
 function transpose(m) {
   const data = new Array(m.rows * m.cols);
@@ -1510,7 +1516,7 @@ function transpose(m) {
       data[r + c * m.cols] = matAt(m, c, r);
     }
   }
-  return chunk2TZKENDH_cjs.matrixValue(m.cols, m.rows, data);
+  return chunkIT74O4BE_cjs.matrixValue(m.cols, m.rows, data);
 }
 function toWorkingRows(m) {
   const rows = [];
@@ -1566,7 +1572,7 @@ function integerDeterminant(rows) {
 }
 function numericDeterminant(m) {
   const integers = toIntegerRows(m);
-  if (integers !== null) return chunk2TZKENDH_cjs.numberValue(Number(integerDeterminant(integers)));
+  if (integers !== null) return chunkIT74O4BE_cjs.numberValue(Number(integerDeterminant(integers)));
   const n = m.rows;
   const a = toWorkingRows(m);
   let det = 1;
@@ -1579,7 +1585,7 @@ function numericDeterminant(m) {
         pivotRow = r;
       }
     }
-    if (maxAbs === 0) return chunk2TZKENDH_cjs.numberValue(0);
+    if (maxAbs === 0) return chunkIT74O4BE_cjs.numberValue(0);
     if (pivotRow !== col) {
       const tmp = a[col];
       a[col] = a[pivotRow];
@@ -1593,7 +1599,7 @@ function numericDeterminant(m) {
       for (let c = col; c < n; c++) a[r][c] -= factor * a[col][c];
     }
   }
-  return chunk2TZKENDH_cjs.numberValue(det);
+  return chunkIT74O4BE_cjs.numberValue(det);
 }
 function symbolicDeterminant(m) {
   const n = m.rows;
@@ -1606,7 +1612,7 @@ function symbolicDeterminant(m) {
   let det = constNode(1);
   for (let col = 0; col < n; col++) {
     const pivot = simplifySymbolic(a[col][col]);
-    if (pivot.kind === "const" && isRationalZero(pivot.value)) return chunk2TZKENDH_cjs.numberValue(0);
+    if (pivot.kind === "const" && isRationalZero(pivot.value)) return chunkIT74O4BE_cjs.numberValue(0);
     det = simplifySymbolic({ kind: "mul", left: det, right: pivot });
     for (let r = col + 1; r < n; r++) {
       const factor = simplifySymbolic({ kind: "div", left: a[r][col], right: pivot });
@@ -1615,11 +1621,11 @@ function symbolicDeterminant(m) {
       }
     }
   }
-  return det.kind === "const" ? chunk2TZKENDH_cjs.numberValue(rationalToNumber(det.value)) : chunk2TZKENDH_cjs.symbolicValue(det);
+  return det.kind === "const" ? chunkIT74O4BE_cjs.numberValue(rationalToNumber(det.value)) : chunkIT74O4BE_cjs.symbolicValue(det);
 }
 function determinant(m) {
   if (!isSquare(m)) {
-    return chunk2TZKENDH_cjs.errorValue("DETERMINANT_REQUIRES_SQUARE_MATRIX", `det: matrix must be square (got ${m.rows}x${m.cols}).`);
+    return chunkIT74O4BE_cjs.errorValue("DETERMINANT_REQUIRES_SQUARE_MATRIX", `det: matrix must be square (got ${m.rows}x${m.cols}).`);
   }
   return m.hasSymbolic ? symbolicDeterminant(m) : numericDeterminant(m);
 }
@@ -1639,7 +1645,7 @@ function numericInverse(m) {
       }
     }
     if (maxAbs < 1e-10) {
-      return chunk2TZKENDH_cjs.errorValue("SINGULAR_MATRIX", `inv: a ${n}x${n} matrix with a zero (or near-zero) pivot column is not invertible.`);
+      return chunkIT74O4BE_cjs.errorValue("SINGULAR_MATRIX", `inv: a ${n}x${n} matrix with a zero (or near-zero) pivot column is not invertible.`);
     }
     if (pivotRow !== col) {
       const tmp = a[col];
@@ -1659,13 +1665,13 @@ function numericInverse(m) {
   for (let r = 0; r < n; r++) {
     for (let c = 0; c < n; c++) data[r + c * n] = a[r][n + c];
   }
-  return chunk2TZKENDH_cjs.matrixValue(n, n, data);
+  return chunkIT74O4BE_cjs.matrixValue(n, n, data);
 }
 var SYMBOLIC_INVERSE_DIMENSION_LIMIT = 9;
 function symbolicInverse(m) {
   const n = m.rows;
   if (n >= SYMBOLIC_INVERSE_DIMENSION_LIMIT) {
-    return chunk2TZKENDH_cjs.errorValue(
+    return chunkIT74O4BE_cjs.errorValue(
       "SYMBOLIC_INVERSE_DIMENSION_LIMIT",
       `Symbolic matrix inversion is limited to matrices smaller than ${SYMBOLIC_INVERSE_DIMENSION_LIMIT}x${SYMBOLIC_INVERSE_DIMENSION_LIMIT} (got ${n}x${n}).`
     );
@@ -1680,7 +1686,7 @@ function symbolicInverse(m) {
   for (let col = 0; col < n; col++) {
     const pivot = simplifySymbolic(a[col][col]);
     if (pivot.kind === "const" && isRationalZero(pivot.value)) {
-      return chunk2TZKENDH_cjs.errorValue(
+      return chunkIT74O4BE_cjs.errorValue(
         "SYMBOLIC_SINGULAR_OR_UNSUPPORTED_PIVOT",
         `Symbolic inverse: the pivot at row/col ${col} is exactly zero \u2014 this matrix's structure isn't invertible via this engine's diagonal-first symbolic elimination (no row-swapping); try reordering rows manually.`
       );
@@ -1701,22 +1707,22 @@ function symbolicInverse(m) {
   for (let r = 0; r < n; r++) {
     for (let c = 0; c < n; c++) data[r + c * n] = symbolicToEntry(a[r][n + c]);
   }
-  return chunk2TZKENDH_cjs.matrixValue(n, n, data);
+  return chunkIT74O4BE_cjs.matrixValue(n, n, data);
 }
 function inverse(m) {
   if (!isSquare(m)) {
-    return chunk2TZKENDH_cjs.errorValue("INVERSE_REQUIRES_SQUARE_MATRIX", `inv: matrix must be square (got ${m.rows}x${m.cols}).`);
+    return chunkIT74O4BE_cjs.errorValue("INVERSE_REQUIRES_SQUARE_MATRIX", `inv: matrix must be square (got ${m.rows}x${m.cols}).`);
   }
   return m.hasSymbolic ? symbolicInverse(m) : numericInverse(m);
 }
 function tooLarge(length, maxElements) {
-  return chunk2TZKENDH_cjs.errorValue(
+  return chunkIT74O4BE_cjs.errorValue(
     "COLLECTION_TOO_LARGE",
     `This collection has ${length} elements, past the limit of ${maxElements} (see the engine's vm.maxCollectionSize setting).`
   );
 }
 function collectionToValues(v, maxElements = Number.POSITIVE_INFINITY) {
-  const faulted = chunk2TZKENDH_cjs.faultedOperand(v);
+  const faulted = chunkIT74O4BE_cjs.faultedOperand(v);
   if (faulted) return faulted;
   if (v.type === 7 /* Matrix */) {
     const m = v.value;
@@ -1728,10 +1734,10 @@ function collectionToValues(v, maxElements = Number.POSITIVE_INFINITY) {
     const length = r.max - r.min + 1;
     if (length > maxElements) return tooLarge(length, maxElements);
     const out = new Array(length);
-    for (let i = 0; i < length; i++) out[i] = chunk2TZKENDH_cjs.numberValue(r.min + i);
+    for (let i = 0; i < length; i++) out[i] = chunkIT74O4BE_cjs.numberValue(r.min + i);
     return out;
   }
-  return chunk2TZKENDH_cjs.errorValue(
+  return chunkIT74O4BE_cjs.errorValue(
     "MAP_REDUCE_REQUIRES_COLLECTION",
     `map/reduce requires a Matrix or Range collection (e.g. "[1,2,3]" or "0:3").`
   );
@@ -1884,6 +1890,7 @@ exports.inverse = inverse;
 exports.isRationalOne = isRationalOne;
 exports.isRationalZero = isRationalZero;
 exports.isZeroPolynomial = isZeroPolynomial;
+exports.makeDecimal = makeDecimal;
 exports.matAt = matAt;
 exports.matIndex = matIndex;
 exports.matrixCompare = matrixCompare;
@@ -1919,5 +1926,5 @@ exports.symbolicToEntry = symbolicToEntry;
 exports.toPolynomial = toPolynomial;
 exports.transpose = transpose;
 exports.varNode = varNode;
-//# sourceMappingURL=chunk-HIQ5HSZL.cjs.map
-//# sourceMappingURL=chunk-HIQ5HSZL.cjs.map
+//# sourceMappingURL=chunk-A2DZ5BEJ.cjs.map
+//# sourceMappingURL=chunk-A2DZ5BEJ.cjs.map