document-compute.js 1.5.8 → 1.5.9

package/dist/index.cjs

@@ -11,20 +11,17 @@ function toRational(value) {
 	};
 }
 function gcd(a, b) {
-	let x = a < 0n ? -a : a;
-	let y = b;
+	let x = b;
+	let y = (a % b + b) % b;
 	while (y !== 0n) [x, y] = [y, x % y];
 	return x;
 }
 function reduce(n, d) {
 	if (d === 0n) throw new RangeError("rational.ts: denominator must not be zero");
-	const sign = d < 0n ? -1n : 1n;
-	const num = n * sign;
-	const den = d * sign;
-	const g = gcd(num, den);
+	const g = gcd(n, d);
 	return {
-		n: num / g,
-		d: den / g
+		n: n / g,
+		d: d / g
 	};
 }
 function toExactRational(value) {
@@ -267,9 +264,7 @@ function negateInterval(a) {
 	return interval(-a.max, -a.min, a.dimension);
 }
 function absInterval(a) {
-	if (a.min >= 0) return a;
-	if (a.max <= 0) return negateInterval(a);
-	return interval(0, Math.max(-a.min, a.max), a.dimension);
+	return interval(Math.max(0, a.min, -a.max), Math.max(-a.min, a.max), a.dimension);
 }
 //#endregion
 //#region src/compute/evaluate.ts
@@ -283,7 +278,7 @@ function isInterval(value) {
 function toInterval(value) {
 	return isInterval(value) ? value : pointInterval(value.magnitude, value.dimension);
 }
-function asQuantity$1(value, context) {
+function asQuantity(value, context) {
 	if (isInterval(value)) throw new UnsupportedExpressionError(context, "this position requires a plain Quantity, not an Interval");
 	return value;
 }
@@ -335,6 +330,9 @@ function evaluate(expression, bindings, context = EMPTY_SYMBOL_TABLE$2) {
 		case "unparsed": throw new UnsupportedExpressionError("evaluate", `this node is source LaTeX ("${expression.latex}") document-schema.js's lowering could not represent structurally, so there is nothing to evaluate`);
 	}
 }
+function evaluateQuantity(expression, bindings, context = EMPTY_SYMBOL_TABLE$2) {
+	return asQuantity(evaluate(expression, bindings, context), "evaluateQuantity");
+}
 function evaluateQty(node, context) {
 	const unit = context.units.find((entry) => entry.id === node.unit);
 	if (unit === void 0) throw new UnknownUnitError(node.unit);
@@ -371,13 +369,13 @@ function evaluateApp(node, bindings, context) {
 	}
 	if (node.operator === "math:pow") {
 		const [base, exponent] = expectTwoArgs(args, "'math:pow'");
-		return powQuantity(asQuantity$1(base, "evaluate"), asQuantity$1(exponent, "evaluate"));
+		return powQuantity(asQuantity(base, "evaluate"), asQuantity(exponent, "evaluate"));
 	}
 	throw new UnsupportedExpressionError("evaluate", `unknown operator '${node.operator}'`);
 }
 function evaluateBinder(node, bindings, context) {
-	const lower = asQuantity$1(evaluate(node.lower, bindings, context), `evaluate:${node.kind}`);
-	const upper = asQuantity$1(evaluate(node.upper, bindings, context), `evaluate:${node.kind}`);
+	const lower = asQuantity(evaluate(node.lower, bindings, context), `evaluate:${node.kind}`);
+	const upper = asQuantity(evaluate(node.upper, bindings, context), `evaluate:${node.kind}`);
 	if (!isDimensionless(lower.dimension) || !isDimensionless(upper.dimension)) throw new IncompatibleDimensionsError(`math:${node.kind}`, lower.dimension, upper.dimension, "binder bounds must be dimensionless");
 	if (!Number.isInteger(lower.magnitude) || !Number.isInteger(upper.magnitude)) throw new UnsupportedExpressionError(`evaluate:${node.kind}`, "binder bounds must evaluate to integers");
 	let accumulator = node.kind === "sum" ? quantity(0, {}) : quantity(1, {});
@@ -386,7 +384,7 @@ function evaluateBinder(node, bindings, context) {
 			...bindings,
 			[node.binder]: quantity(i, {})
 		};
-		const bodyValue = asQuantity$1(evaluate(node.body, bodyBindings, context), `evaluate:${node.kind}`);
+		const bodyValue = asQuantity(evaluate(node.body, bodyBindings, context), `evaluate:${node.kind}`);
 		accumulator = node.kind === "sum" ? addQuantities(accumulator, bodyValue) : multiplyQuantities(accumulator, bodyValue);
 	}
 	return accumulator;
@@ -430,7 +428,7 @@ function bisection(f, bracket, tolerance, maxIterations) {
 		const mid = (low + high) / 2;
 		const fMid = f(mid);
 		if (Math.abs(fMid) < tolerance) return mid;
-		if (fMid > 0 === fLow > 0) {
+		if (Math.sign(fMid) === Math.sign(fLow)) {
 			low = mid;
 			fLow = fMid;
 		} else high = mid;
@@ -445,9 +443,7 @@ function newton(f, initialGuess, tolerance, maxIterations, h) {
 		if (Math.abs(fx) < tolerance) return x;
 		const derivative = (f(x + h) - f(x - h)) / (2 * h);
 		if (!Number.isFinite(derivative) || Math.abs(derivative) < 1e-14) throw new NonConvergentSolveError("newton", i, `the numeric derivative vanished or diverged near x=${x}`);
-		const next = x - fx / derivative;
-		if (!Number.isFinite(next)) throw new NonConvergentSolveError("newton", i, `the iteration diverged to a non-finite value near x=${x}`);
-		x = next;
+		x = x - fx / derivative;
 	}
 	throw new NonConvergentSolveError("newton", maxIterations, `residual still exceeds tolerance ${tolerance} after ${maxIterations} iterations`);
 }
@@ -471,19 +467,13 @@ function asEquality(expression) {
 function containsSymbol(expression) {
 	switch (expression.kind) {
 		case "sym": return true;
-		case "num":
-		case "qty":
-		case "unparsed": return false;
 		case "app": return expression.args.some(containsSymbol);
 		case "sum":
 		case "prod": return containsSymbol(expression.lower) || containsSymbol(expression.upper) || containsSymbol(expression.body);
 		case "matrix": return expression.rows.some((row) => row.some(containsSymbol));
+		default: return false;
 	}
 }
-function asQuantity(value) {
-	if (isInterval(value)) throw new UnsupportedExpressionError("runWorkedExampleSequence", "this harness compares point-valued Quantity answers only; a symbol resolving to a range (Interval) has no stated-answer comparison defined yet");
-	return value;
-}
 function gapFromError(error) {
 	if (error instanceof UnboundSymbolError) return "unbound-symbol";
 	if (error instanceof UnknownUnitError) return "unknown-unit";
@@ -491,7 +481,6 @@ function gapFromError(error) {
 	if (error instanceof DivisionByZeroError) return "division-by-zero";
 	if (error instanceof UnsupportedExpressionError) return "unsupported-construct";
 	if (error instanceof NumericDomainError) return "numeric-domain";
-	if (error instanceof NonConvergentSolveError) return "non-convergent-solve";
 	return "other-evaluation-error";
 }
 function errorMessage(error) {
@@ -533,7 +522,7 @@ function runWorkedExampleSequence(formulas, symbolTable = EMPTY_SYMBOL_TABLE, op
 		}
 		let closedValue;
 		try {
-			closedValue = asQuantity(evaluate(rhs, EMPTY_BINDINGS, symbolTable));
+			closedValue = evaluateQuantity(rhs, EMPTY_BINDINGS, symbolTable);
 		} catch (error) {
 			outcomes.push({
 				outcome: "gap",
@@ -548,7 +537,7 @@ function runWorkedExampleSequence(formulas, symbolTable = EMPTY_SYMBOL_TABLE, op
 			pending = void 0;
 			let actual;
 			try {
-				actual = asQuantity(evaluate(definitionRhs, bindings, symbolTable));
+				actual = evaluateQuantity(definitionRhs, bindings, symbolTable);
 			} catch (error) {
 				outcomes.push({
 					outcome: "gap",
@@ -680,6 +669,7 @@ exports.divideIntervals = divideIntervals;
 exports.divideQuantities = divideQuantities;
 exports.divideRational = divideRational;
 exports.evaluate = evaluate;
+exports.evaluateQuantity = evaluateQuantity;
 exports.formatCorpusReport = formatCorpusReport;
 exports.interval = interval;
 exports.isDimensionless = isDimensionless;

package/dist/index.d.cts

@@ -82,6 +82,7 @@ declare function absInterval(a: Interval): Interval;
 type EvaluationResult = Quantity | Interval;
 declare function isInterval(value: EvaluationResult): value is Interval;
 declare function evaluate(expression: MathExpression, bindings: FormulaBindings, context?: SymbolTable): EvaluationResult;
+declare function evaluateQuantity(expression: MathExpression, bindings: FormulaBindings, context?: SymbolTable): Quantity;
 //#endregion
 //#region src/compute/solve.d.ts
 type SolveMethod = "bisection" | "newton";
@@ -104,7 +105,7 @@ interface SolveForOptions {
 declare function solveFor(expression: MathExpression, targetValue: number, unknownSymbol: string, bindings: FormulaBindings, options?: SolveForOptions, context?: SymbolTable): number;
 //#endregion
 //#region src/harness/worked-example.d.ts
-type WorkedExampleGap = "unbound-symbol" | "unknown-unit" | "incompatible-dimensions" | "division-by-zero" | "unsupported-construct" | "numeric-domain" | "non-convergent-solve" | "other-evaluation-error";
+type WorkedExampleGap = "unbound-symbol" | "unknown-unit" | "incompatible-dimensions" | "division-by-zero" | "unsupported-construct" | "numeric-domain" | "other-evaluation-error";
 interface WorkedExampleMatch {
   readonly outcome: "match";
   readonly targetSymbol: string;
@@ -165,4 +166,4 @@ interface CorpusReport {
 declare function runCorpus(documents: readonly CorpusDocument[], options?: WorkedExampleOptions): CorpusReport;
 declare function formatCorpusReport(report: CorpusReport): string;
 //#endregion
-export { CorpusDocument, CorpusDocumentReport, CorpusReport, DivisionByZeroError, EvaluationResult, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, Rational, SolveForOptions, SolveMethod, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, WorkedExampleGap, WorkedExampleGapResult, WorkedExampleMatch, WorkedExampleMismatch, WorkedExampleOptions, WorkedExampleOutcome, WorkedExampleReport, WorkedExampleUnresolved, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };
+export { CorpusDocument, CorpusDocumentReport, CorpusReport, DivisionByZeroError, EvaluationResult, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, Rational, SolveForOptions, SolveMethod, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, WorkedExampleGap, WorkedExampleGapResult, WorkedExampleMatch, WorkedExampleMismatch, WorkedExampleOptions, WorkedExampleOutcome, WorkedExampleReport, WorkedExampleUnresolved, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, evaluateQuantity, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };

package/dist/index.d.ts

@@ -82,6 +82,7 @@ declare function absInterval(a: Interval): Interval;
 type EvaluationResult = Quantity | Interval;
 declare function isInterval(value: EvaluationResult): value is Interval;
 declare function evaluate(expression: MathExpression, bindings: FormulaBindings, context?: SymbolTable): EvaluationResult;
+declare function evaluateQuantity(expression: MathExpression, bindings: FormulaBindings, context?: SymbolTable): Quantity;
 //#endregion
 //#region src/compute/solve.d.ts
 type SolveMethod = "bisection" | "newton";
@@ -104,7 +105,7 @@ interface SolveForOptions {
 declare function solveFor(expression: MathExpression, targetValue: number, unknownSymbol: string, bindings: FormulaBindings, options?: SolveForOptions, context?: SymbolTable): number;
 //#endregion
 //#region src/harness/worked-example.d.ts
-type WorkedExampleGap = "unbound-symbol" | "unknown-unit" | "incompatible-dimensions" | "division-by-zero" | "unsupported-construct" | "numeric-domain" | "non-convergent-solve" | "other-evaluation-error";
+type WorkedExampleGap = "unbound-symbol" | "unknown-unit" | "incompatible-dimensions" | "division-by-zero" | "unsupported-construct" | "numeric-domain" | "other-evaluation-error";
 interface WorkedExampleMatch {
   readonly outcome: "match";
   readonly targetSymbol: string;
@@ -165,4 +166,4 @@ interface CorpusReport {
 declare function runCorpus(documents: readonly CorpusDocument[], options?: WorkedExampleOptions): CorpusReport;
 declare function formatCorpusReport(report: CorpusReport): string;
 //#endregion
-export { CorpusDocument, CorpusDocumentReport, CorpusReport, DivisionByZeroError, EvaluationResult, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, Rational, SolveForOptions, SolveMethod, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, WorkedExampleGap, WorkedExampleGapResult, WorkedExampleMatch, WorkedExampleMismatch, WorkedExampleOptions, WorkedExampleOutcome, WorkedExampleReport, WorkedExampleUnresolved, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };
+export { CorpusDocument, CorpusDocumentReport, CorpusReport, DivisionByZeroError, EvaluationResult, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, Rational, SolveForOptions, SolveMethod, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, WorkedExampleGap, WorkedExampleGapResult, WorkedExampleMatch, WorkedExampleMismatch, WorkedExampleOptions, WorkedExampleOutcome, WorkedExampleReport, WorkedExampleUnresolved, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, evaluateQuantity, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };

package/dist/index.js

@@ -10,20 +10,17 @@ function toRational(value) {
 	};
 }
 function gcd(a, b) {
-	let x = a < 0n ? -a : a;
-	let y = b;
+	let x = b;
+	let y = (a % b + b) % b;
 	while (y !== 0n) [x, y] = [y, x % y];
 	return x;
 }
 function reduce(n, d) {
 	if (d === 0n) throw new RangeError("rational.ts: denominator must not be zero");
-	const sign = d < 0n ? -1n : 1n;
-	const num = n * sign;
-	const den = d * sign;
-	const g = gcd(num, den);
+	const g = gcd(n, d);
 	return {
-		n: num / g,
-		d: den / g
+		n: n / g,
+		d: d / g
 	};
 }
 function toExactRational(value) {
@@ -266,9 +263,7 @@ function negateInterval(a) {
 	return interval(-a.max, -a.min, a.dimension);
 }
 function absInterval(a) {
-	if (a.min >= 0) return a;
-	if (a.max <= 0) return negateInterval(a);
-	return interval(0, Math.max(-a.min, a.max), a.dimension);
+	return interval(Math.max(0, a.min, -a.max), Math.max(-a.min, a.max), a.dimension);
 }
 //#endregion
 //#region src/compute/evaluate.ts
@@ -282,7 +277,7 @@ function isInterval(value) {
 function toInterval(value) {
 	return isInterval(value) ? value : pointInterval(value.magnitude, value.dimension);
 }
-function asQuantity$1(value, context) {
+function asQuantity(value, context) {
 	if (isInterval(value)) throw new UnsupportedExpressionError(context, "this position requires a plain Quantity, not an Interval");
 	return value;
 }
@@ -334,6 +329,9 @@ function evaluate(expression, bindings, context = EMPTY_SYMBOL_TABLE$2) {
 		case "unparsed": throw new UnsupportedExpressionError("evaluate", `this node is source LaTeX ("${expression.latex}") document-schema.js's lowering could not represent structurally, so there is nothing to evaluate`);
 	}
 }
+function evaluateQuantity(expression, bindings, context = EMPTY_SYMBOL_TABLE$2) {
+	return asQuantity(evaluate(expression, bindings, context), "evaluateQuantity");
+}
 function evaluateQty(node, context) {
 	const unit = context.units.find((entry) => entry.id === node.unit);
 	if (unit === void 0) throw new UnknownUnitError(node.unit);
@@ -370,13 +368,13 @@ function evaluateApp(node, bindings, context) {
 	}
 	if (node.operator === "math:pow") {
 		const [base, exponent] = expectTwoArgs(args, "'math:pow'");
-		return powQuantity(asQuantity$1(base, "evaluate"), asQuantity$1(exponent, "evaluate"));
+		return powQuantity(asQuantity(base, "evaluate"), asQuantity(exponent, "evaluate"));
 	}
 	throw new UnsupportedExpressionError("evaluate", `unknown operator '${node.operator}'`);
 }
 function evaluateBinder(node, bindings, context) {
-	const lower = asQuantity$1(evaluate(node.lower, bindings, context), `evaluate:${node.kind}`);
-	const upper = asQuantity$1(evaluate(node.upper, bindings, context), `evaluate:${node.kind}`);
+	const lower = asQuantity(evaluate(node.lower, bindings, context), `evaluate:${node.kind}`);
+	const upper = asQuantity(evaluate(node.upper, bindings, context), `evaluate:${node.kind}`);
 	if (!isDimensionless(lower.dimension) || !isDimensionless(upper.dimension)) throw new IncompatibleDimensionsError(`math:${node.kind}`, lower.dimension, upper.dimension, "binder bounds must be dimensionless");
 	if (!Number.isInteger(lower.magnitude) || !Number.isInteger(upper.magnitude)) throw new UnsupportedExpressionError(`evaluate:${node.kind}`, "binder bounds must evaluate to integers");
 	let accumulator = node.kind === "sum" ? quantity(0, {}) : quantity(1, {});
@@ -385,7 +383,7 @@ function evaluateBinder(node, bindings, context) {
 			...bindings,
 			[node.binder]: quantity(i, {})
 		};
-		const bodyValue = asQuantity$1(evaluate(node.body, bodyBindings, context), `evaluate:${node.kind}`);
+		const bodyValue = asQuantity(evaluate(node.body, bodyBindings, context), `evaluate:${node.kind}`);
 		accumulator = node.kind === "sum" ? addQuantities(accumulator, bodyValue) : multiplyQuantities(accumulator, bodyValue);
 	}
 	return accumulator;
@@ -429,7 +427,7 @@ function bisection(f, bracket, tolerance, maxIterations) {
 		const mid = (low + high) / 2;
 		const fMid = f(mid);
 		if (Math.abs(fMid) < tolerance) return mid;
-		if (fMid > 0 === fLow > 0) {
+		if (Math.sign(fMid) === Math.sign(fLow)) {
 			low = mid;
 			fLow = fMid;
 		} else high = mid;
@@ -444,9 +442,7 @@ function newton(f, initialGuess, tolerance, maxIterations, h) {
 		if (Math.abs(fx) < tolerance) return x;
 		const derivative = (f(x + h) - f(x - h)) / (2 * h);
 		if (!Number.isFinite(derivative) || Math.abs(derivative) < 1e-14) throw new NonConvergentSolveError("newton", i, `the numeric derivative vanished or diverged near x=${x}`);
-		const next = x - fx / derivative;
-		if (!Number.isFinite(next)) throw new NonConvergentSolveError("newton", i, `the iteration diverged to a non-finite value near x=${x}`);
-		x = next;
+		x = x - fx / derivative;
 	}
 	throw new NonConvergentSolveError("newton", maxIterations, `residual still exceeds tolerance ${tolerance} after ${maxIterations} iterations`);
 }
@@ -470,19 +466,13 @@ function asEquality(expression) {
 function containsSymbol(expression) {
 	switch (expression.kind) {
 		case "sym": return true;
-		case "num":
-		case "qty":
-		case "unparsed": return false;
 		case "app": return expression.args.some(containsSymbol);
 		case "sum":
 		case "prod": return containsSymbol(expression.lower) || containsSymbol(expression.upper) || containsSymbol(expression.body);
 		case "matrix": return expression.rows.some((row) => row.some(containsSymbol));
+		default: return false;
 	}
 }
-function asQuantity(value) {
-	if (isInterval(value)) throw new UnsupportedExpressionError("runWorkedExampleSequence", "this harness compares point-valued Quantity answers only; a symbol resolving to a range (Interval) has no stated-answer comparison defined yet");
-	return value;
-}
 function gapFromError(error) {
 	if (error instanceof UnboundSymbolError) return "unbound-symbol";
 	if (error instanceof UnknownUnitError) return "unknown-unit";
@@ -490,7 +480,6 @@ function gapFromError(error) {
 	if (error instanceof DivisionByZeroError) return "division-by-zero";
 	if (error instanceof UnsupportedExpressionError) return "unsupported-construct";
 	if (error instanceof NumericDomainError) return "numeric-domain";
-	if (error instanceof NonConvergentSolveError) return "non-convergent-solve";
 	return "other-evaluation-error";
 }
 function errorMessage(error) {
@@ -532,7 +521,7 @@ function runWorkedExampleSequence(formulas, symbolTable = EMPTY_SYMBOL_TABLE, op
 		}
 		let closedValue;
 		try {
-			closedValue = asQuantity(evaluate(rhs, EMPTY_BINDINGS, symbolTable));
+			closedValue = evaluateQuantity(rhs, EMPTY_BINDINGS, symbolTable);
 		} catch (error) {
 			outcomes.push({
 				outcome: "gap",
@@ -547,7 +536,7 @@ function runWorkedExampleSequence(formulas, symbolTable = EMPTY_SYMBOL_TABLE, op
 			pending = void 0;
 			let actual;
 			try {
-				actual = asQuantity(evaluate(definitionRhs, bindings, symbolTable));
+				actual = evaluateQuantity(definitionRhs, bindings, symbolTable);
 			} catch (error) {
 				outcomes.push({
 					outcome: "gap",
@@ -657,4 +646,4 @@ function formatCorpusReport(report) {
 	return lines.join("\n");
 }
 //#endregion
-export { DivisionByZeroError, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };
+export { DivisionByZeroError, IncompatibleDimensionsError, NonConvergentSolveError, NumericDomainError, UnboundSymbolError, UnknownUnitError, UnsupportedExpressionError, absInterval, absQuantity, addIntervals, addQuantities, addRational, collectFormulas, cosQuantity, dimensionExponent, dimensionToString, dimensionsEqual, divideDimensions, divideIntervals, divideQuantities, divideRational, evaluate, evaluateQuantity, formatCorpusReport, interval, isDimensionless, isInterval, multiplyDimensions, multiplyIntervals, multiplyQuantities, multiplyRational, negateInterval, negateQuantity, pointInterval, powQuantity, quantity, rationalToNumber, runCorpus, runWorkedExampleSequence, scaleDimension, sinQuantity, solveFor, sqrtQuantity, subtractIntervals, subtractQuantities, subtractRational, tanQuantity, toExactRational, toRational };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "document-compute.js",
-  "version": "1.5.8",
+  "version": "1.5.9",
   "description": "Units-typed, tree-walking evaluator for document-schema.js's MathExpression -- exact-rational unit conversion, interval arithmetic, and bisection/Newton numeric solve-for, the compute package for the documents.js family.",
   "type": "module",
   "repository": {
@@ -79,7 +79,7 @@
     "@stryker-mutator/typescript-checker": "10.0.0",
     "@stryker-mutator/vitest-runner": "10.0.0",
     "@types/node": "26.2.0",
-    "documents.js": "7.20.5",
+    "documents.js": "7.20.6",
     "eslint": "10.8.1",
     "husky": "9.1.7",
     "markdown-codec": "6.7.1",