universal-physics-tensor 9.0.0 → 11.0.0

package/dist/numerical/finite-difference.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"finite-difference.js","sourceRoot":"","sources":["../../src/numerical/finite-difference.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH;;;;;;GAMG;AACH,MAAM,UAAU,cAAc,CAC5B,CAA4C,EAC5C,KAAuC,EACvC,GAAW,EACX,IAAY;IAEZ,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAE,CAAC;IACtB,MAAM,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC;IACpB,MAAM,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC;IACtB,OAAO,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,CAAC;AACnF,CAAC","sourcesContent":["/**\n * One central-difference rule.\n *\n * Three readers differentiated a closed form by central differences, each\n * with its own step rule and its own denominator (`2h` at one, the\n * representable step at another). The derivative is here, once: the caller\n * passes the step its rule chose, and the denominator is the step the\n * floating-point point actually took, `dx = (x + h) − (x − h)`, so a\n * perturbation that rounds at an astrophysical `x` is divided by what was\n * applied, not by what was asked.\n *\n * @module numerical/finite-difference\n */\n\n/**\n * `∂f/∂x_key` at `point` by the central difference with half-step `step`:\n * `(f(x + h) − f(x − h)) / ((x + h) − (x − h))`. The other keys of `point` are\n * held. The step rule is the caller's; a zero or non-finite `step` makes the\n * denominator 0 or NaN, which the caller refuses before calling.\n * @internal\n */\nexport function centralPartial(\n  f: (point: Record<string, number>) => number,\n  point: Readonly<Record<string, number>>,\n  key: string,\n  step: number,\n): number {\n  const x = point[key]!;\n  const up = x + step;\n  const down = x - step;\n  return (f({ ...point, [key]: up }) - f({ ...point, [key]: down })) / (up - down);\n}\n"]}

package/dist/numerical/formula-contract.d.ts

@@ -1,17 +1,31 @@
 /**
- * Scalar-formula contract shared by the MathTS parser.
+ * Scalar-formula contract: the parser types, the refusal of `euler`, and THE
+ * ONE TABLE of documented scalar functions.
  *
- * Path B, the recursive-descent parser that used to live in `formula.ts`, is
- * gone. These types and the documented function table stay: the MathTS
- * adapter still returns a `number`, still refuses `euler`, and still supplies
- * `ln` when MathTS does not define it.
+ * `SCALAR_FUNCTIONS` carries, for each documented function, its arity, its
+ * numeric body, and its dimension rule (how the arguments' dimensions map to
+ * the result's). Three readers consume it and none keeps a table of its
+ * own: the MathTS formula parser supplies the bodies MathTS does not define
+ * (`ln`) and names the documented list in its diagnostics; the lowering pass
+ * evaluates a `transcendental`, `abs` or `^` node through it; the dimension
+ * checker and the binding reader apply the dimension rule. Path B, the
+ * recursive-descent parser that used to live in `formula.ts`, is gone.
  *
  * @module numerical/formula-contract
  */
+import type { TranscendentalFn } from '../dimensional/ast-types.js';
 /** A parse or evaluation failure (bad syntax, unknown symbol, arity). */
 export declare class FormulaError extends Error {
     constructor(message: string);
 }
+/**
+ * A formula that parsed and evaluated, whose value is not a finite real number (`ln(-1)`,
+ * `1/0`, `sqrt(-1)`). A bad value at the point, not a malformed formula: the CLI exits 1 on
+ * it where a `FormulaError` is exit 2 (Tom's review, second round).
+ */
+export declare class FormulaValueError extends FormulaError {
+    constructor(message: string);
+}
 /** A parsed formula: its free variables and a safe evaluator. @internal */
 export interface CompiledFormula {
     /** The original source string. */
@@ -35,8 +49,50 @@ export declare const EULER_NUMBER_ERROR = "euler is not Euler's number. Write it
 export declare class EulerNumberError extends FormulaError {
     constructor();
 }
+/**
+ * How a documented function maps its arguments' dimensions to its result's.
+ *
+ * - `transcendental`: the argument must be dimensionless and so is the result.
+ *   `node` is the dimensional grammar's `transcendental` node when the grammar
+ *   has one for this function; without it the checker keeps a dimensionless stub.
+ * - `same`: the result has the argument's dimension (`abs`).
+ * - `root`: the result is the argument's dimension to `power` (`sqrt`, `cbrt`).
+ * - `power`: `pow(base, exp)`; the exponent is a dimensionless constant and
+ *   the result is the base's dimension to that constant.
+ * - `ratio`: `atan2(y, x)`; both arguments share a dimension and the result is
+ *   dimensionless.
+ * @internal
+ */
+export type ScalarFunctionDimension = {
+    readonly kind: 'transcendental';
+    readonly node?: TranscendentalFn;
+} | {
+    readonly kind: 'same';
+} | {
+    readonly kind: 'root';
+    readonly power: number;
+} | {
+    readonly kind: 'power';
+} | {
+    readonly kind: 'ratio';
+};
+/** One documented scalar function: arity, numeric body, dimension rule. @internal */
+export interface ScalarFunction {
+    readonly arity: number;
+    readonly dimension: ScalarFunctionDimension;
+    /** The numeric body. Arity is checked by {@link callBuiltinFunction}, not here. */
+    readonly apply: (args: readonly number[]) => number;
+}
+/**
+ * The documented scalar functions, in the order help lists them. `log` is the
+ * natural logarithm (physics convention) and lowers to the grammar's `ln` node.
+ * @internal
+ */
+export declare const SCALAR_FUNCTIONS: Readonly<Record<string, ScalarFunction>>;
 /** Names of the built-in functions this parser documents. @internal */
 export declare const BUILTIN_FUNCTION_NAMES: readonly string[];
+/** `${name} expects ${arity} argument(s)`, the one arity sentence. @internal */
+export declare function arityMessage(name: string, arity: number): string;
 /**
  * Function names from other conventions, each with the documented function that computes the same
  * value (audit I4). An unknown name still fails: the equivalent is only named in the message, so a
@@ -56,9 +112,9 @@ export declare const BUILTIN_FUNCTION_LIST: string;
  */
 export declare function unknownFunctionMessage(name: string): string;
 /**
- * Call a documented function by name, with this parser's arity checks. Supplies
+ * Call a documented function by name, with the table's arity check. Supplies
  * a function MathTS does not define (`ln`).
  * @internal
  */
-export declare function callBuiltinFunction(name: string, args: number[]): number;
+export declare function callBuiltinFunction(name: string, args: readonly number[]): number;
 //# sourceMappingURL=formula-contract.d.ts.map

package/dist/numerical/formula-contract.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"formula-contract.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-contract.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,yEAAyE;AACzE,qBAAa,YAAa,SAAQ,KAAK;IACrC,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED,2EAA2E;AAC3E,MAAM,WAAW,eAAe;IAC9B,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,uEAAuE;IACvE,QAAQ,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,CAAC;IACtC,wEAAwE;IACxE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,CAAC;CACjD;AAED,0DAA0D;AAC1D,MAAM,WAAW,aAAa;IAC5B,KAAK,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAAC;CACtC;AAED;;;;GAIG;AACH,eAAO,MAAM,kBAAkB,yHACyF,CAAC;AAEzH,iEAAiE;AACjE,qBAAa,gBAAiB,SAAQ,YAAY;IAChD,cAGC;CACF;AAmCD,uEAAuE;AACvE,eAAO,MAAM,sBAAsB,EAAE,SAAS,MAAM,EAA2B,CAAC;AAEhF;;;;;;GAMG;AACH,eAAO,MAAM,oBAAoB,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE;IAAE,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAA;CAAE,CAAC,CAY7G,CAAC;AAEF,2GAA2G;AAC3G,eAAO,MAAM,qBAAqB,EAAE,MAExB,CAAC;AAEb;;;GAGG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,MAAM,GAAG,MAAM,CAW3D;AAED;;;;GAIG;AACH,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,GAAG,MAAM,CAIxE"}
+{"version":3,"file":"formula-contract.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-contract.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,6BAA6B,CAAC;AAEpE,yEAAyE;AACzE,qBAAa,YAAa,SAAQ,KAAK;IACrC,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED;;;;GAIG;AACH,qBAAa,iBAAkB,SAAQ,YAAY;IACjD,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED,2EAA2E;AAC3E,MAAM,WAAW,eAAe;IAC9B,kCAAkC;IAClC,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,uEAAuE;IACvE,QAAQ,CAAC,SAAS,EAAE,SAAS,MAAM,EAAE,CAAC;IACtC,wEAAwE;IACxE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,GAAG,MAAM,CAAC;CACjD;AAED,0DAA0D;AAC1D,MAAM,WAAW,aAAa;IAC5B,KAAK,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAAC;CACtC;AAED;;;;GAIG;AACH,eAAO,MAAM,kBAAkB,yHACyF,CAAC;AAEzH,iEAAiE;AACjE,qBAAa,gBAAiB,SAAQ,YAAY;IAChD,cAGC;CACF;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,MAAM,uBAAuB,GAC/B;IAAE,QAAQ,CAAC,IAAI,EAAE,gBAAgB,CAAC;IAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,gBAAgB,CAAA;CAAE,GACrE;IAAE,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAA;CAAE,GACzB;IAAE,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAA;CAAE,GACjD;IAAE,QAAQ,CAAC,IAAI,EAAE,OAAO,CAAA;CAAE,GAC1B;IAAE,QAAQ,CAAC,IAAI,EAAE,OAAO,CAAA;CAAE,CAAC;AAE/B,qFAAqF;AACrF,MAAM,WAAW,cAAc;IAC7B,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,QAAQ,CAAC,SAAS,EAAE,uBAAuB,CAAC;IAC5C,mFAAmF;IACnF,QAAQ,CAAC,KAAK,EAAE,CAAC,IAAI,EAAE,SAAS,MAAM,EAAE,KAAK,MAAM,CAAC;CACrD;AAQD;;;;GAIG;AACH,eAAO,MAAM,gBAAgB,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,cAAc,CAAC,CAoBrE,CAAC;AAEF,uEAAuE;AACvE,eAAO,MAAM,sBAAsB,EAAE,SAAS,MAAM,EAAkC,CAAC;AAEvF,gFAAgF;AAChF,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAEhE;AAED;;;;;;GAMG;AACH,eAAO,MAAM,oBAAoB,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE;IAAE,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IAAC,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAA;CAAE,CAAC,CAY7G,CAAC;AAEF,2GAA2G;AAC3G,eAAO,MAAM,qBAAqB,EAAE,MAExB,CAAC;AAEb;;;GAGG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,MAAM,GAAG,MAAM,CAW3D;AAED;;;;GAIG;AACH,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM,CAKjF"}

package/dist/numerical/formula-contract.js

@@ -1,10 +1,15 @@
 /**
- * Scalar-formula contract shared by the MathTS parser.
+ * Scalar-formula contract: the parser types, the refusal of `euler`, and THE
+ * ONE TABLE of documented scalar functions.
  *
- * Path B, the recursive-descent parser that used to live in `formula.ts`, is
- * gone. These types and the documented function table stay: the MathTS
- * adapter still returns a `number`, still refuses `euler`, and still supplies
- * `ln` when MathTS does not define it.
+ * `SCALAR_FUNCTIONS` carries, for each documented function, its arity, its
+ * numeric body, and its dimension rule (how the arguments' dimensions map to
+ * the result's). Three readers consume it and none keeps a table of its
+ * own: the MathTS formula parser supplies the bodies MathTS does not define
+ * (`ln`) and names the documented list in its diagnostics; the lowering pass
+ * evaluates a `transcendental`, `abs` or `^` node through it; the dimension
+ * checker and the binding reader apply the dimension rule. Path B, the
+ * recursive-descent parser that used to live in `formula.ts`, is gone.
  *
  * @module numerical/formula-contract
  */
@@ -15,6 +20,17 @@ export class FormulaError extends Error {
         this.name = 'FormulaError';
     }
 }
+/**
+ * A formula that parsed and evaluated, whose value is not a finite real number (`ln(-1)`,
+ * `1/0`, `sqrt(-1)`). A bad value at the point, not a malformed formula: the CLI exits 1 on
+ * it where a `FormulaError` is exit 2 (Tom's review, second round).
+ */
+export class FormulaValueError extends FormulaError {
+    constructor(message) {
+        super(message);
+        this.name = 'FormulaValueError';
+    }
+}
 /**
  * The name `euler` is not Euler's number. Write `exp(x)`, for example `exp(1)`.
  * A bare `e` is the elementary charge. `E` is energy.
@@ -28,42 +44,43 @@ export class EulerNumberError extends FormulaError {
         this.name = 'EulerNumberError';
     }
 }
-const arity1 = (f) => (a) => {
-    if (a.length !== 1)
-        throw new FormulaError('expected 1 argument');
-    return f(a[0]);
-};
-const FUNCTIONS = {
-    sqrt: arity1(Math.sqrt),
-    cbrt: arity1(Math.cbrt),
-    exp: arity1(Math.exp),
-    ln: arity1(Math.log),
-    log: arity1(Math.log), // natural log (physics convention)
-    log10: arity1(Math.log10),
-    log2: arity1(Math.log2),
-    abs: arity1(Math.abs),
-    sin: arity1(Math.sin),
-    cos: arity1(Math.cos),
-    tan: arity1(Math.tan),
-    asin: arity1(Math.asin),
-    acos: arity1(Math.acos),
-    atan: arity1(Math.atan),
-    sinh: arity1(Math.sinh),
-    cosh: arity1(Math.cosh),
-    tanh: arity1(Math.tanh),
-    pow: (a) => {
-        if (a.length !== 2)
-            throw new FormulaError('pow expects 2 arguments');
-        return Math.pow(a[0], a[1]);
-    },
-    atan2: (a) => {
-        if (a.length !== 2)
-            throw new FormulaError('atan2 expects 2 arguments');
-        return Math.atan2(a[0], a[1]);
-    },
+const transcendental = (node, f) => ({
+    arity: 1,
+    dimension: node === undefined ? { kind: 'transcendental' } : { kind: 'transcendental', node },
+    apply: (a) => f(a[0]),
+});
+/**
+ * The documented scalar functions, in the order help lists them. `log` is the
+ * natural logarithm (physics convention) and lowers to the grammar's `ln` node.
+ * @internal
+ */
+export const SCALAR_FUNCTIONS = {
+    sqrt: { arity: 1, dimension: { kind: 'root', power: 0.5 }, apply: (a) => Math.sqrt(a[0]) },
+    cbrt: { arity: 1, dimension: { kind: 'root', power: 1 / 3 }, apply: (a) => Math.cbrt(a[0]) },
+    exp: transcendental('exp', Math.exp),
+    ln: transcendental('ln', Math.log),
+    log: transcendental('ln', Math.log),
+    log10: transcendental('log10', Math.log10),
+    log2: transcendental('log2', Math.log2),
+    abs: { arity: 1, dimension: { kind: 'same' }, apply: (a) => Math.abs(a[0]) },
+    sin: transcendental('sin', Math.sin),
+    cos: transcendental('cos', Math.cos),
+    tan: transcendental('tan', Math.tan),
+    asin: transcendental(undefined, Math.asin),
+    acos: transcendental(undefined, Math.acos),
+    atan: transcendental(undefined, Math.atan),
+    sinh: transcendental('sinh', Math.sinh),
+    cosh: transcendental('cosh', Math.cosh),
+    tanh: transcendental('tanh', Math.tanh),
+    pow: { arity: 2, dimension: { kind: 'power' }, apply: (a) => Math.pow(a[0], a[1]) },
+    atan2: { arity: 2, dimension: { kind: 'ratio' }, apply: (a) => Math.atan2(a[0], a[1]) },
 };
 /** Names of the built-in functions this parser documents. @internal */
-export const BUILTIN_FUNCTION_NAMES = Object.keys(FUNCTIONS);
+export const BUILTIN_FUNCTION_NAMES = Object.keys(SCALAR_FUNCTIONS);
+/** `${name} expects ${arity} argument(s)`, the one arity sentence. @internal */
+export function arityMessage(name, arity) {
+    return `${name} expects ${arity} argument${arity === 1 ? '' : 's'}`;
+}
 /**
  * Function names from other conventions, each with the documented function that computes the same
  * value (audit I4). An unknown name still fails: the equivalent is only named in the message, so a
@@ -101,14 +118,16 @@ export function unknownFunctionMessage(name) {
     return `unknown function '${name}'.${hint} Documented functions: ${BUILTIN_FUNCTION_LIST}.`;
 }
 /**
- * Call a documented function by name, with this parser's arity checks. Supplies
+ * Call a documented function by name, with the table's arity check. Supplies
  * a function MathTS does not define (`ln`).
  * @internal
  */
 export function callBuiltinFunction(name, args) {
-    const fn = FUNCTIONS[name];
+    const fn = SCALAR_FUNCTIONS[name];
     if (fn === undefined)
         throw new FormulaError(unknownFunctionMessage(name));
-    return fn(args);
+    if (args.length !== fn.arity)
+        throw new FormulaError(arityMessage(name, fn.arity));
+    return fn.apply(args);
 }
 //# sourceMappingURL=formula-contract.js.map

package/dist/numerical/formula-contract.js.map

@@ -1 +1 @@
-{"version":3,"file":"formula-contract.js","sourceRoot":"","sources":["../../src/numerical/formula-contract.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AAEH,yEAAyE;AACzE,MAAM,OAAO,YAAa,SAAQ,KAAK;IACrC,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,cAAc,CAAC;IAC7B,CAAC;CACF;AAiBD;;;;GAIG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAC7B,sHAAsH,CAAC;AAEzH,iEAAiE;AACjE,MAAM,OAAO,gBAAiB,SAAQ,YAAY;IAChD;QACE,KAAK,CAAC,kBAAkB,CAAC,CAAC;QAC1B,IAAI,CAAC,IAAI,GAAG,kBAAkB,CAAC;IACjC,CAAC;CACF;AAGD,MAAM,MAAM,GAAG,CAAC,CAAwB,EAAM,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE;IACrD,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC;QAAE,MAAM,IAAI,YAAY,CAAC,qBAAqB,CAAC,CAAC;IAClE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;AAClB,CAAC,CAAC;AACF,MAAM,SAAS,GAAiC;IAC9C,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACrB,EAAE,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACpB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,mCAAmC;IAC1D,KAAK,EAAE,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC;IACzB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACrB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACrB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACrB,GAAG,EAAE,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC;IACrB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,IAAI,EAAE,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC;IACvB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;QACT,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,YAAY,CAAC,yBAAyB,CAAC,CAAC;QACtE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;IAChC,CAAC;IACD,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE;QACX,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,YAAY,CAAC,2BAA2B,CAAC,CAAC;QACxE,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,CAAC;IAClC,CAAC;CACF,CAAC;AAEF,uEAAuE;AACvE,MAAM,CAAC,MAAM,sBAAsB,GAAsB,MAAM,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;AAEhF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAiF;IAChH,EAAE,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,uBAAuB,EAAE;IACtD,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,sBAAsB,EAAE;IACpD,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,sBAAsB,EAAE;IACpD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,kBAAkB,EAAE;IACpD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,oBAAoB,EAAE;IACtD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,qBAAqB,EAAE;IACvD,KAAK,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,qBAAqB,EAAE;IACtD,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,aAAa,EAAE;IAC1C,IAAI,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,oBAAoB,EAAE;IACnD,KAAK,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,SAAS,EAAE;IACzC,OAAO,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,kCAAkC,EAAE;CACvE,CAAC;AAEF,2GAA2G;AAC3G,MAAM,CAAC,MAAM,qBAAqB,GAAW,sBAAsB,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAC5E,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,qBAAqB,CAAC,CAAC,CAAC,CAAC,CACxC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAEb;;;GAGG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAAY;IACjD,MAAM,KAAK,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,CAAC,WAAW,EAAE,CAAC;IACjC,MAAM,GAAG,GACP,KAAK,EAAE,GAAG,IAAI,CAAC,KAAK,KAAK,IAAI,IAAI,sBAAsB,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,oBAAoB,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;IACtH,MAAM,OAAO,GAAG,KAAK,EAAE,OAAO,IAAI,oBAAoB,CAAC,KAAK,CAAC,EAAE,OAAO,CAAC;IACvE,MAAM,IAAI,GACR,GAAG,KAAK,SAAS;QACf,CAAC,CAAC,EAAE;QACJ,CAAC,CAAC,QAAQ,OAAO,IAAI,gBAAgB,GAAG,EAAE,QAAQ,GAAG,0CAA0C,CAAC;IACpG,OAAO,qBAAqB,IAAI,KAAK,IAAI,0BAA0B,qBAAqB,GAAG,CAAC;AAC9F,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,IAAY,EAAE,IAAc;IAC9D,MAAM,EAAE,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;IAC3B,IAAI,EAAE,KAAK,SAAS;QAAE,MAAM,IAAI,YAAY,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC,CAAC;IAC3E,OAAO,EAAE,CAAC,IAAI,CAAC,CAAC;AAClB,CAAC"}
+{"version":3,"file":"formula-contract.js","sourceRoot":"","sources":["../../src/numerical/formula-contract.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAIH,yEAAyE;AACzE,MAAM,OAAO,YAAa,SAAQ,KAAK;IACrC,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,cAAc,CAAC;IAC7B,CAAC;CACF;AAED;;;;GAIG;AACH,MAAM,OAAO,iBAAkB,SAAQ,YAAY;IACjD,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,mBAAmB,CAAC;IAClC,CAAC;CACF;AAiBD;;;;GAIG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAC7B,sHAAsH,CAAC;AAEzH,iEAAiE;AACjE,MAAM,OAAO,gBAAiB,SAAQ,YAAY;IAChD;QACE,KAAK,CAAC,kBAAkB,CAAC,CAAC;QAC1B,IAAI,CAAC,IAAI,GAAG,kBAAkB,CAAC;IACjC,CAAC;CACF;AA+BD,MAAM,cAAc,GAAG,CAAC,IAAkC,EAAE,CAAwB,EAAkB,EAAE,CAAC,CAAC;IACxG,KAAK,EAAE,CAAC;IACR,SAAS,EAAE,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,gBAAgB,EAAE,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,gBAAgB,EAAE,IAAI,EAAE;IAC7F,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC;CACvB,CAAC,CAAC;AAEH;;;;GAIG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAA6C;IACxE,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,EAAE;IAC3F,IAAI,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,EAAE;IAC7F,GAAG,EAAE,cAAc,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC;IACpC,EAAE,EAAE,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC;IAClC,GAAG,EAAE,cAAc,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC;IACnC,KAAK,EAAE,cAAc,CAAC,OAAO,EAAE,IAAI,CAAC,KAAK,CAAC;IAC1C,IAAI,EAAE,cAAc,CAAC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC;IACvC,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,CAAC,EAAE;IAC7E,GAAG,EAAE,cAAc,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC;IACpC,GAAG,EAAE,cAAc,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC;IACpC,GAAG,EAAE,cAAc,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC;IACpC,IAAI,EAAE,cAAc,CAAC,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC;IAC1C,IAAI,EAAE,cAAc,CAAC,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC;IAC1C,IAAI,EAAE,cAAc,CAAC,SAAS,EAAE,IAAI,CAAC,IAAI,CAAC;IAC1C,IAAI,EAAE,cAAc,CAAC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC;IACvC,IAAI,EAAE,cAAc,CAAC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC;IACvC,IAAI,EAAE,cAAc,CAAC,MAAM,EAAE,IAAI,CAAC,IAAI,CAAC;IACvC,GAAG,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,EAAE;IACrF,KAAK,EAAE,EAAE,KAAK,EAAE,CAAC,EAAE,SAAS,EAAE,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAE,EAAE,CAAC,CAAC,CAAC,CAAE,CAAC,EAAE;CAC1F,CAAC;AAEF,uEAAuE;AACvE,MAAM,CAAC,MAAM,sBAAsB,GAAsB,MAAM,CAAC,IAAI,CAAC,gBAAgB,CAAC,CAAC;AAEvF,gFAAgF;AAChF,MAAM,UAAU,YAAY,CAAC,IAAY,EAAE,KAAa;IACtD,OAAO,GAAG,IAAI,YAAY,KAAK,YAAY,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;AACtE,CAAC;AAED;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,oBAAoB,GAAiF;IAChH,EAAE,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,uBAAuB,EAAE;IACtD,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,sBAAsB,EAAE;IACpD,EAAE,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,sBAAsB,EAAE;IACpD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,kBAAkB,EAAE;IACpD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,oBAAoB,EAAE;IACtD,MAAM,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,qBAAqB,EAAE;IACvD,KAAK,EAAE,EAAE,GAAG,EAAE,MAAM,EAAE,OAAO,EAAE,qBAAqB,EAAE;IACtD,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,aAAa,EAAE;IAC1C,IAAI,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,oBAAoB,EAAE;IACnD,KAAK,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,OAAO,EAAE,SAAS,EAAE;IACzC,OAAO,EAAE,EAAE,GAAG,EAAE,OAAO,EAAE,OAAO,EAAE,kCAAkC,EAAE;CACvE,CAAC;AAEF,2GAA2G;AAC3G,MAAM,CAAC,MAAM,qBAAqB,GAAW,sBAAsB,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAC5E,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,qBAAqB,CAAC,CAAC,CAAC,CAAC,CACxC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AAEb;;;GAGG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAAY;IACjD,MAAM,KAAK,GAAG,oBAAoB,CAAC,IAAI,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,CAAC,WAAW,EAAE,CAAC;IACjC,MAAM,GAAG,GACP,KAAK,EAAE,GAAG,IAAI,CAAC,KAAK,KAAK,IAAI,IAAI,sBAAsB,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,oBAAoB,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC;IACtH,MAAM,OAAO,GAAG,KAAK,EAAE,OAAO,IAAI,oBAAoB,CAAC,KAAK,CAAC,EAAE,OAAO,CAAC;IACvE,MAAM,IAAI,GACR,GAAG,KAAK,SAAS;QACf,CAAC,CAAC,EAAE;QACJ,CAAC,CAAC,QAAQ,OAAO,IAAI,gBAAgB,GAAG,EAAE,QAAQ,GAAG,0CAA0C,CAAC;IACpG,OAAO,qBAAqB,IAAI,KAAK,IAAI,0BAA0B,qBAAqB,GAAG,CAAC;AAC9F,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,mBAAmB,CAAC,IAAY,EAAE,IAAuB;IACvE,MAAM,EAAE,GAAG,gBAAgB,CAAC,IAAI,CAAC,CAAC;IAClC,IAAI,EAAE,KAAK,SAAS;QAAE,MAAM,IAAI,YAAY,CAAC,sBAAsB,CAAC,IAAI,CAAC,CAAC,CAAC;IAC3E,IAAI,IAAI,CAAC,MAAM,KAAK,EAAE,CAAC,KAAK;QAAE,MAAM,IAAI,YAAY,CAAC,YAAY,CAAC,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC;IACnF,OAAO,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;AACxB,CAAC","sourcesContent":["/**\n * Scalar-formula contract: the parser types, the refusal of `euler`, and THE\n * ONE TABLE of documented scalar functions.\n *\n * `SCALAR_FUNCTIONS` carries, for each documented function, its arity, its\n * numeric body, and its dimension rule (how the arguments' dimensions map to\n * the result's). Three readers consume it and none keeps a table of its\n * own: the MathTS formula parser supplies the bodies MathTS does not define\n * (`ln`) and names the documented list in its diagnostics; the lowering pass\n * evaluates a `transcendental`, `abs` or `^` node through it; the dimension\n * checker and the binding reader apply the dimension rule. Path B, the\n * recursive-descent parser that used to live in `formula.ts`, is gone.\n *\n * @module numerical/formula-contract\n */\n\nimport type { TranscendentalFn } from '../dimensional/ast-types.js';\n\n/** A parse or evaluation failure (bad syntax, unknown symbol, arity). */\nexport class FormulaError extends Error {\n  constructor(message: string) {\n    super(message);\n    this.name = 'FormulaError';\n  }\n}\n\n/**\n * A formula that parsed and evaluated, whose value is not a finite real number (`ln(-1)`,\n * `1/0`, `sqrt(-1)`). A bad value at the point, not a malformed formula: the CLI exits 1 on\n * it where a `FormulaError` is exit 2 (Tom's review, second round).\n */\nexport class FormulaValueError extends FormulaError {\n  constructor(message: string) {\n    super(message);\n    this.name = 'FormulaValueError';\n  }\n}\n\n/** A parsed formula: its free variables and a safe evaluator. @internal */\nexport interface CompiledFormula {\n  /** The original source string. */\n  readonly source: string;\n  /** Free variable names (excludes built-in constants and functions). */\n  readonly variables: readonly string[];\n  /** Evaluate against a `name → value` scope. Throws on a missing var. */\n  evaluate(scope: Record<string, number>): number;\n}\n\n/** The MathTS-backed parser implements this. @internal */\nexport interface FormulaParser {\n  parse(expr: string): CompiledFormula;\n}\n\n/**\n * The name `euler` is not Euler's number. Write `exp(x)`, for example `exp(1)`.\n * A bare `e` is the elementary charge. `E` is energy.\n * @internal\n */\nexport const EULER_NUMBER_ERROR =\n  \"euler is not Euler's number. Write it as exp(x), for example exp(1). A bare e is the elementary charge. E is energy.\";\n\n/** The name `euler` was written for Euler's number. @internal */\nexport class EulerNumberError extends FormulaError {\n  constructor() {\n    super(EULER_NUMBER_ERROR);\n    this.name = 'EulerNumberError';\n  }\n}\n\n/**\n * How a documented function maps its arguments' dimensions to its result's.\n *\n * - `transcendental`: the argument must be dimensionless and so is the result.\n *   `node` is the dimensional grammar's `transcendental` node when the grammar\n *   has one for this function; without it the checker keeps a dimensionless stub.\n * - `same`: the result has the argument's dimension (`abs`).\n * - `root`: the result is the argument's dimension to `power` (`sqrt`, `cbrt`).\n * - `power`: `pow(base, exp)`; the exponent is a dimensionless constant and\n *   the result is the base's dimension to that constant.\n * - `ratio`: `atan2(y, x)`; both arguments share a dimension and the result is\n *   dimensionless.\n * @internal\n */\nexport type ScalarFunctionDimension =\n  | { readonly kind: 'transcendental'; readonly node?: TranscendentalFn }\n  | { readonly kind: 'same' }\n  | { readonly kind: 'root'; readonly power: number }\n  | { readonly kind: 'power' }\n  | { readonly kind: 'ratio' };\n\n/** One documented scalar function: arity, numeric body, dimension rule. @internal */\nexport interface ScalarFunction {\n  readonly arity: number;\n  readonly dimension: ScalarFunctionDimension;\n  /** The numeric body. Arity is checked by {@link callBuiltinFunction}, not here. */\n  readonly apply: (args: readonly number[]) => number;\n}\n\nconst transcendental = (node: TranscendentalFn | undefined, f: (x: number) => number): ScalarFunction => ({\n  arity: 1,\n  dimension: node === undefined ? { kind: 'transcendental' } : { kind: 'transcendental', node },\n  apply: (a) => f(a[0]!),\n});\n\n/**\n * The documented scalar functions, in the order help lists them. `log` is the\n * natural logarithm (physics convention) and lowers to the grammar's `ln` node.\n * @internal\n */\nexport const SCALAR_FUNCTIONS: Readonly<Record<string, ScalarFunction>> = {\n  sqrt: { arity: 1, dimension: { kind: 'root', power: 0.5 }, apply: (a) => Math.sqrt(a[0]!) },\n  cbrt: { arity: 1, dimension: { kind: 'root', power: 1 / 3 }, apply: (a) => Math.cbrt(a[0]!) },\n  exp: transcendental('exp', Math.exp),\n  ln: transcendental('ln', Math.log),\n  log: transcendental('ln', Math.log),\n  log10: transcendental('log10', Math.log10),\n  log2: transcendental('log2', Math.log2),\n  abs: { arity: 1, dimension: { kind: 'same' }, apply: (a) => Math.abs(a[0]!) },\n  sin: transcendental('sin', Math.sin),\n  cos: transcendental('cos', Math.cos),\n  tan: transcendental('tan', Math.tan),\n  asin: transcendental(undefined, Math.asin),\n  acos: transcendental(undefined, Math.acos),\n  atan: transcendental(undefined, Math.atan),\n  sinh: transcendental('sinh', Math.sinh),\n  cosh: transcendental('cosh', Math.cosh),\n  tanh: transcendental('tanh', Math.tanh),\n  pow: { arity: 2, dimension: { kind: 'power' }, apply: (a) => Math.pow(a[0]!, a[1]!) },\n  atan2: { arity: 2, dimension: { kind: 'ratio' }, apply: (a) => Math.atan2(a[0]!, a[1]!) },\n};\n\n/** Names of the built-in functions this parser documents. @internal */\nexport const BUILTIN_FUNCTION_NAMES: readonly string[] = Object.keys(SCALAR_FUNCTIONS);\n\n/** `${name} expects ${arity} argument(s)`, the one arity sentence. @internal */\nexport function arityMessage(name: string, arity: number): string {\n  return `${name} expects ${arity} argument${arity === 1 ? '' : 's'}`;\n}\n\n/**\n * Function names from other conventions, each with the documented function that computes the same\n * value (audit I4). An unknown name still fails: the equivalent is only named in the message, so a\n * formula never means something its author did not write. Every `use` is in\n * {@link BUILTIN_FUNCTION_NAMES}.\n * @internal\n */\nexport const FUNCTION_EQUIVALENTS: Readonly<Record<string, { readonly use: string; readonly meaning: string }>> = {\n  lg: { use: 'log10', meaning: 'the base-10 logarithm' },\n  lb: { use: 'log2', meaning: 'the base-2 logarithm' },\n  ld: { use: 'log2', meaning: 'the base-2 logarithm' },\n  arcsin: { use: 'asin', meaning: 'the inverse sine' },\n  arccos: { use: 'acos', meaning: 'the inverse cosine' },\n  arctan: { use: 'atan', meaning: 'the inverse tangent' },\n  arctg: { use: 'atan', meaning: 'the inverse tangent' },\n  tg: { use: 'tan', meaning: 'the tangent' },\n  fabs: { use: 'abs', meaning: 'the absolute value' },\n  power: { use: 'pow', meaning: 'a power' },\n  arctan2: { use: 'atan2', meaning: 'the two-argument inverse tangent' },\n};\n\n/** The documented functions as help and diagnostics list them, with the base of `log` stated. @internal */\nexport const BUILTIN_FUNCTION_LIST: string = BUILTIN_FUNCTION_NAMES.map((n) =>\n  n === 'log' ? 'log (natural, = ln)' : n,\n).join(', ');\n\n/**\n * The diagnostic for a function this parser does not document: the documented equivalent when the\n * name is a known alias or differs only in case, then the documented list. @internal\n */\nexport function unknownFunctionMessage(name: string): string {\n  const alias = FUNCTION_EQUIVALENTS[name];\n  const lower = name.toLowerCase();\n  const use =\n    alias?.use ?? (lower !== name && BUILTIN_FUNCTION_NAMES.includes(lower) ? lower : FUNCTION_EQUIVALENTS[lower]?.use);\n  const meaning = alias?.meaning ?? FUNCTION_EQUIVALENTS[lower]?.meaning;\n  const hint =\n    use === undefined\n      ? ''\n      : ` For ${meaning ?? `the function ${use}`} use ${use}(…), which the formula parser evaluates.`;\n  return `unknown function '${name}'.${hint} Documented functions: ${BUILTIN_FUNCTION_LIST}.`;\n}\n\n/**\n * Call a documented function by name, with the table's arity check. Supplies\n * a function MathTS does not define (`ln`).\n * @internal\n */\nexport function callBuiltinFunction(name: string, args: readonly number[]): number {\n  const fn = SCALAR_FUNCTIONS[name];\n  if (fn === undefined) throw new FormulaError(unknownFunctionMessage(name));\n  if (args.length !== fn.arity) throw new FormulaError(arityMessage(name, fn.arity));\n  return fn.apply(args);\n}\n"]}

package/dist/numerical/formula-dimension.d.ts

@@ -32,13 +32,6 @@ export declare class UnsupportedSyntaxError extends FormulaDimensionError {
  * @internal
  */
 export declare function formulaSymbolDimension(name: string): Dimension | undefined;
-/**
- * What to say when a bare `e` was the elementary charge inside a sum that is
- * not homogeneous. Declaring or binding `e` keeps a different quantity.
- * Euler's number is `exp(x)`.
- * @internal
- */
-export declare const ELEMENTARY_CHARGE_MIX_MESSAGE = "e is the elementary charge and is not dimensionless here. Declare it or bind a value (e=<number>) if you mean a different quantity, or write Euler's number as exp(x). The catalog writes the eccentricity factor as one_minus_e_sq.";
 /** Normalized parse node — the shape the MathTS AST adapts to. @internal */
 type PNode = {
     kind: 'num';
@@ -64,6 +57,14 @@ type PNode = {
 };
 /** The normalized node {@link parseFormulaPNode} returns. @internal */
 export type FormulaPNode = PNode;
+/**
+ * The number a numeric-constant exponent folds to (`2`, `-1`, `1/2`,
+ * `(2+1)/3`), or undefined when the node names a symbol or a call. The
+ * binding reader uses the same folding for its exponents, so the two
+ * readers agree on what an exponent may be.
+ * @internal
+ */
+export declare function constantExponent(node: PNode): number | undefined;
 /** The result of a formula dimensional check. */
 interface FormulaDimensionResult {
     /** True iff the formula is dimensionally homogeneous and well-formed. */

package/dist/numerical/formula-dimension.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"formula-dimension.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-dimension.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAGH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAGzD,OAAO,KAAK,EAAE,QAAQ,EAAoB,MAAM,6BAA6B,CAAC;AAK9E;;cAEc;AACd,qBAAa,qBAAsB,SAAQ,KAAK;IAC9C,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED,6GAA6G;AAC7G,qBAAa,sBAAuB,SAAQ,qBAAqB;IACnD,QAAQ,CAAC,QAAQ,EAAE,MAAM;IAArC,YAAqB,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAGpD;CACF;AA4BD;;;;;GAKG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,MAAM,GAAG,SAAS,GAAG,SAAS,CAM1E;AAuDD;;;;;GAKG;AACH,eAAO,MAAM,6BAA6B,yOAC+L,CAAC;AA2C1O,4EAA4E;AAC5E,KAAK,KAAK,GACN;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,GAC9B;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,GAC7B;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,GAAG,EAAE,KAAK,CAAA;CAAE,GAC3B;IAAE,IAAI,EAAE,IAAI,CAAC;IAAC,EAAE,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,IAAI,EAAE,KAAK,EAAE,CAAA;CAAE,GACxD;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,IAAI,EAAE,KAAK,CAAC;IAAC,GAAG,EAAE,KAAK,CAAA;CAAE,GACxC;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,KAAK,EAAE,CAAA;CAAE,CAAC;AAEhD,uEAAuE;AACvE,MAAM,MAAM,YAAY,GAAG,KAAK,CAAC;AAmGjC,iDAAiD;AACjD,UAAU,sBAAsB;IAC9B,yEAAyE;IACzE,QAAQ,CAAC,EAAE,EAAE,OAAO,CAAC;IACrB,wCAAwC;IACxC,QAAQ,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC;IACzB,wCAAwC;IACxC,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;IACxB,0FAA0F;IAC1F,QAAQ,CAAC,qBAAqB,CAAC,EAAE,IAAI,CAAC;IACtC,gGAAgG;IAChG,QAAQ,CAAC,gBAAgB,CAAC,EAAE,MAAM,CAAC;CACpC;AAED,oFAAoF;AACpF,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,IAAI,EAAE,QAAQ,CAAC;IACxB,QAAQ,CAAC,SAAS,EAAE,SAAS,CAAC;CAC/B;AAED,kFAAkF;AAClF,MAAM,WAAW,uBAAuB;IACtC,KAAK,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GAAG,sBAAsB,CAAC;IACvF;+EAC2E;IAC3E,KAAK,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GAAG,aAAa,CAAC;CAC/E;AAmDD;;;;GAIG;AACH,wBAAgB,iBAAiB,CAAC,IAAI,EAAE,MAAM,GAAG,KAAK,CAIrD;AA0BD,8EAA8E;AAC9E,wBAAgB,8BAA8B,IAAI,uBAAuB,CAExE"}
+{"version":3,"file":"formula-dimension.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-dimension.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAGH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAGzD,OAAO,KAAK,EAAE,QAAQ,EAAoB,MAAM,6BAA6B,CAAC;AAK9E;;cAEc;AACd,qBAAa,qBAAsB,SAAQ,KAAK;IAC9C,YAAY,OAAO,EAAE,MAAM,EAG1B;CACF;AAED,6GAA6G;AAC7G,qBAAa,sBAAuB,SAAQ,qBAAqB;IACnD,QAAQ,CAAC,QAAQ,EAAE,MAAM;IAArC,YAAqB,QAAQ,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAGpD;CACF;AA4BD;;;;;GAKG;AACH,wBAAgB,sBAAsB,CAAC,IAAI,EAAE,MAAM,GAAG,SAAS,GAAG,SAAS,CAM1E;AAwHD,4EAA4E;AAC5E,KAAK,KAAK,GACN;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,GAC9B;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,GAC7B;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,GAAG,EAAE,KAAK,CAAA;CAAE,GAC3B;IAAE,IAAI,EAAE,IAAI,CAAC;IAAC,EAAE,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,IAAI,EAAE,KAAK,EAAE,CAAA;CAAE,GACxD;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,IAAI,EAAE,KAAK,CAAC;IAAC,GAAG,EAAE,KAAK,CAAA;CAAE,GACxC;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,KAAK,EAAE,CAAA;CAAE,CAAC;AAEhD,uEAAuE;AACvE,MAAM,MAAM,YAAY,GAAG,KAAK,CAAC;AAkCjC;;;;;;GAMG;AACH,wBAAgB,gBAAgB,CAAC,IAAI,EAAE,KAAK,GAAG,MAAM,GAAG,SAAS,CAOhE;AA6DD,iDAAiD;AACjD,UAAU,sBAAsB;IAC9B,yEAAyE;IACzE,QAAQ,CAAC,EAAE,EAAE,OAAO,CAAC;IACrB,wCAAwC;IACxC,QAAQ,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC;IACzB,wCAAwC;IACxC,QAAQ,CAAC,KAAK,CAAC,EAAE,MAAM,CAAC;IACxB,0FAA0F;IAC1F,QAAQ,CAAC,qBAAqB,CAAC,EAAE,IAAI,CAAC;IACtC,gGAAgG;IAChG,QAAQ,CAAC,gBAAgB,CAAC,EAAE,MAAM,CAAC;CACpC;AAED,oFAAoF;AACpF,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,IAAI,EAAE,QAAQ,CAAC;IACxB,QAAQ,CAAC,SAAS,EAAE,SAAS,CAAC;CAC/B;AAED,kFAAkF;AAClF,MAAM,WAAW,uBAAuB;IACtC,KAAK,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GAAG,sBAAsB,CAAC;IACvF;+EAC2E;IAC3E,KAAK,CAAC,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,SAAS,CAAC,CAAC,GAAG,aAAa,CAAC;CAC/E;AAsDD;;;;GAIG;AACH,wBAAgB,iBAAiB,CAAC,IAAI,EAAE,MAAM,GAAG,KAAK,CAIrD;AA0BD,8EAA8E;AAC9E,wBAAgB,8BAA8B,IAAI,uBAAuB,CAExE"}

package/dist/numerical/formula-dimension.js

@@ -17,7 +17,7 @@ import { CHARGE, DIMENSIONLESS, ENERGY } from '../dimensional/types.js';
 import { equals, format, multiply } from '../dimensional/algebra.js';
 import { validate } from '../dimensional/validator.js';
 import { sym } from '../dimensional/ast-builders.js';
-import { EULER_NUMBER_ERROR, EulerNumberError, FormulaError } from './formula-contract.js';
+import { EULER_NUMBER_ERROR, EulerNumberError, FormulaError, SCALAR_FUNCTIONS, unknownFunctionMessage } from './formula-contract.js';
 /** A formula cannot be dimensionally analyzed (undeclared symbol, variable
  *  exponent, transcendental of a dimensional argument, unsupported node).
  *  @public */
@@ -58,7 +58,7 @@ class ElementaryChargeMixError extends FormulaDimensionError {
         this.name = 'ElementaryChargeMixError';
     }
 }
-/** Dimensionless math constants both parsers recognize. Bare `e` is not here:
+/** Dimensionless math constants the parser recognizes. Bare `e` is not here:
  *  ISO 80000 names that symbol the elementary charge. Euler's number is `exp(x)`. */
 const MATH_CONSTANTS = new Set(['pi', 'tau', 'phi', 'Infinity', 'NaN']);
 /**
@@ -78,29 +78,10 @@ export function formulaSymbolDimension(name) {
         return DIMENSIONLESS;
     return undefined;
 }
-/** Parser function names → the dimensional grammar's `transcendental` node fn
- *  (dimensionless → dimensionless). `log` is natural log (mathjs convention). */
-const TRANSCENDENTAL_FN = {
-    exp: 'exp', ln: 'ln', log: 'ln', log10: 'log10', log2: 'log2',
-    sin: 'sin', cos: 'cos', tan: 'tan', sinh: 'sinh', cosh: 'cosh', tanh: 'tanh',
-};
-/** Transcendentals WITHOUT a grammar node — kept as dimensionless stubs (still
- *  dimensionless→dimensionless, just not a faithful node). */
-const TRANSCENDENTAL_STUB = new Set(['asin', 'acos', 'atan', 'sec', 'csc', 'cot']);
 const op = (o, args) => ({ kind: 'op', op: o, args });
 const powExpr = (base, exp) => op('^', [base, sym(String(exp), DIMENSIONLESS)]);
 const transcendental = (fn, arg) => ({ kind: 'transcendental', fn, arg });
 const absNode = (arg) => ({ kind: 'abs', arg });
-/** Dispatch a parsed function name to its `ExprNode`. Shared by both transpilers. */
-function transpileFunction(fn, argExpr) {
-    if (fn === 'abs')
-        return absNode(argExpr);
-    if (fn in TRANSCENDENTAL_FN)
-        return transcendental(TRANSCENDENTAL_FN[fn], argExpr);
-    if (TRANSCENDENTAL_STUB.has(fn))
-        return transcendentalStub(fn, argExpr);
-    return null;
-}
 /** Resolve a symbol to a dimensioned `ExprNode` (declared dim, or a
  *  dimensionless math constant, else an error). */
 function resolveSymbol(name, dims) {
@@ -132,9 +113,8 @@ function dimensionIsCharge(dim) {
  * What to say when a bare `e` was the elementary charge inside a sum that is
  * not homogeneous. Declaring or binding `e` keeps a different quantity.
  * Euler's number is `exp(x)`.
- * @internal
  */
-export const ELEMENTARY_CHARGE_MIX_MESSAGE = 'e is the elementary charge and is not dimensionless here. Declare it or bind a value (e=<number>) if you mean a different quantity, or write Euler\'s number as exp(x). The catalog writes the eccentricity factor as one_minus_e_sq.';
+const ELEMENTARY_CHARGE_MIX_MESSAGE = 'e is the elementary charge and is not dimensionless here. Declare it or bind a value (e=<number>) if you mean a different quantity, or write Euler\'s number as exp(x). The catalog writes the eccentricity factor as one_minus_e_sq.';
 /** Inferred dimension of an `ExprNode`, or throw if not homogeneous. */
 function dimensionOf(node) {
     const r = validate(node);
@@ -152,6 +132,39 @@ function transcendentalStub(fn, argExprNode) {
     }
     return sym(`${fn}(...)`, DIMENSIONLESS);
 }
+/**
+ * A call, dispatched by the one function table's dimension rule. The arity is
+ * the table's: a call with another number of arguments is a dimension error,
+ * not a crash on a missing argument or a silently dropped extra one.
+ */
+function callToExpr(fn, args, dims) {
+    const spec = SCALAR_FUNCTIONS[fn];
+    if (spec === undefined)
+        throw new FormulaDimensionError(unknownFunctionMessage(fn));
+    if (args.length !== spec.arity) {
+        throw new FormulaDimensionError(`${fn}() expects ${spec.arity} argument${spec.arity === 1 ? '' : 's'}, got ${args.length}`);
+    }
+    const rule = spec.dimension;
+    switch (rule.kind) {
+        case 'root':
+            return powExpr(normToExpr(args[0], dims), rule.power);
+        case 'power':
+            return powOf(normToExpr(args[0], dims), args[1], dims);
+        case 'same':
+            return absNode(normToExpr(args[0], dims));
+        case 'transcendental': {
+            const arg = normToExpr(args[0], dims);
+            return rule.node === undefined ? transcendentalStub(fn, arg) : transcendental(rule.node, arg);
+        }
+        case 'ratio': {
+            const [y, x] = [dimensionOf(normToExpr(args[0], dims)), dimensionOf(normToExpr(args[1], dims))];
+            if (!equals(y, x)) {
+                throw new FormulaDimensionError(`${fn}() arguments must have the same dimension, got ${format(y)} and ${format(x)}`);
+            }
+            return sym(`${fn}(...)`, DIMENSIONLESS);
+        }
+    }
+}
 /** Adapt a MathTS AST node to a `PNode` (pure shape map; no dimensions). */
 function mathtsToPNode(node) {
     switch (node.type) {
@@ -185,6 +198,23 @@ function pnodeConstant(node) {
         throw new FormulaDimensionError('exponent must be a finite numeric constant');
     return v;
 }
+/**
+ * The number a numeric-constant exponent folds to (`2`, `-1`, `1/2`,
+ * `(2+1)/3`), or undefined when the node names a symbol or a call. The
+ * binding reader uses the same folding for its exponents, so the two
+ * readers agree on what an exponent may be.
+ * @internal
+ */
+export function constantExponent(node) {
+    try {
+        return pnodeConstant(node);
+    }
+    catch (error) {
+        if (error instanceof FormulaDimensionError)
+            return undefined;
+        throw error;
+    }
+}
 function evalConstPNode(node) {
     switch (node.kind) {
         case 'num':
@@ -238,19 +268,8 @@ function normToExpr(node, dims) {
             return op(node.op, node.args.map((a) => normToExpr(a, dims)));
         case 'pow':
             return powOf(normToExpr(node.base, dims), node.exp, dims);
-        case 'call': {
-            const fn = node.fn;
-            if (fn === 'sqrt')
-                return powExpr(normToExpr(node.args[0], dims), 0.5);
-            if (fn === 'cbrt')
-                return powExpr(normToExpr(node.args[0], dims), 1 / 3);
-            if (fn === 'pow')
-                return powOf(normToExpr(node.args[0], dims), node.args[1], dims);
-            const fnNode = transpileFunction(fn, normToExpr(node.args[0], dims));
-            if (fnNode)
-                return fnNode;
-            throw new FormulaDimensionError(`unsupported function '${fn}'`);
-        }
+        case 'call':
+            return callToExpr(node.fn, node.args, dims);
     }
 }
 /** Build a checker from an `(expr, dims) → ExprNode` transpile (parse +
@@ -295,6 +314,9 @@ function createFormulaDimensionChecker(toExprNode) {
 }
 function parsedMathNode(expr) {
     try {
+        // MathTS's Node type is not imported; `MathNode` is the structural subset
+        // this adapter reads (`type`, `value`, `name`, `op`, `args`, `content`,
+        // `fn`), every field of which exists on a MathTS node of that `type`.
         return parseMathTs(expr);
     }
     catch (e) {

package/dist/numerical/formula-dimension.js.map

@@ -1 +1 @@
-{"version":3,"file":"formula-dimension.js","sourceRoot":"","sources":["../../src/numerical/formula-dimension.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,KAAK,IAAI,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAEvE,OAAO,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,MAAM,yBAAyB,CAAC;AACxE,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAErE,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AACvD,OAAO,EAAE,GAAG,EAAE,MAAM,gCAAgC,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,gBAAgB,EAAE,YAAY,EAAE,MAAM,uBAAuB,CAAC;AAE3F;;cAEc;AACd,MAAM,OAAO,qBAAsB,SAAQ,KAAK;IAC9C,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;IACtC,CAAC;CACF;AAED,6GAA6G;AAC7G,MAAM,OAAO,sBAAuB,SAAQ,qBAAqB;IAC1C,QAAQ;IAA7B,YAAqB,QAAgB,EAAE,MAAc;QACnD,KAAK,CAAC,MAAM,CAAC,CAAC;wBADK,QAAQ;QAE3B,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;IACvC,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,qBAAsB,SAAQ,qBAAqB;IAClC,MAAM;IAA3B,YAAqB,MAAc;QACjC,KAAK,CAAC,sBAAsB,MAAM,sFAAsF,CAAC,CAAC;sBADvG,MAAM;QAEzB,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;IACtC,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,wBAAyB,SAAQ,qBAAqB;IAC1D,YAAY,IAAY;QACtB,KAAK,CAAC,GAAG,6BAA6B,IAAI,IAAI,EAAE,CAAC,CAAC;QAClD,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;IACzC,CAAC;CACF;AAED;qFACqF;AACrF,MAAM,cAAc,GAAG,IAAI,GAAG,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,KAAK,CAAC,CAAC,CAAC;AAExE;;;;;GAKG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAAY;IACjD,IAAI,IAAI,KAAK,GAAG;QAAE,OAAO,MAAM,CAAC;IAChC,IAAI,IAAI,KAAK,GAAG;QAAE,OAAO,MAAM,CAAC;IAChC,IAAI,IAAI,KAAK,OAAO;QAAE,OAAO,SAAS,CAAC;IACvC,IAAI,cAAc,CAAC,GAAG,CAAC,IAAI,CAAC;QAAE,OAAO,aAAa,CAAC;IACnD,OAAO,SAAS,CAAC;AACnB,CAAC;AAED;iFACiF;AACjF,MAAM,iBAAiB,GAA+C;IACpE,GAAG,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,IAAI,EAAE,MAAM;IAC7D,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM;CAC7E,CAAC;AAEF;8DAC8D;AAC9D,MAAM,mBAAmB,GAAG,IAAI,GAAG,CAAC,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC,CAAC;AAEnF,MAAM,EAAE,GAAG,CAAC,CAA8B,EAAE,IAAgB,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC;AACzG,MAAM,OAAO,GAAG,CAAC,IAAc,EAAE,GAAW,EAAY,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;AAC5G,MAAM,cAAc,GAAG,CAAC,EAAoB,EAAE,GAAa,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,gBAAgB,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;AAChH,MAAM,OAAO,GAAG,CAAC,GAAa,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC;AAEpE,qFAAqF;AACrF,SAAS,iBAAiB,CAAC,EAAU,EAAE,OAAiB;IACtD,IAAI,EAAE,KAAK,KAAK;QAAE,OAAO,OAAO,CAAC,OAAO,CAAC,CAAC;IAC1C,IAAI,EAAE,IAAI,iBAAiB;QAAE,OAAO,cAAc,CAAC,iBAAiB,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;IACnF,IAAI,mBAAmB,CAAC,GAAG,CAAC,EAAE,CAAC;QAAE,OAAO,kBAAkB,CAAC,EAAE,EAAE,OAAO,CAAC,CAAC;IACxE,OAAO,IAAI,CAAC;AACd,CAAC;AAED;mDACmD;AACnD,SAAS,aAAa,CAAC,IAAY,EAAE,IAAyC;IAC5E,IAAI,IAAI,KAAK,OAAO,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC;QAAE,MAAM,IAAI,qBAAqB,CAAC,kBAAkB,CAAC,CAAC;IAC7F,IAAI,IAAI,IAAI,IAAI;QAAE,OAAO,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,sBAAsB,CAAC,IAAI,CAAC,CAAC;IAC3C,IAAI,KAAK,KAAK,SAAS;QAAE,OAAO,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IACjD,MAAM,IAAI,qBAAqB,CAAC,IAAI,CAAC,CAAC;AACxC,CAAC;AAED,MAAM,cAAc,GAAG,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;AAEhD,6EAA6E;AAC7E,SAAS,eAAe,CAAC,IAAY;IACnC,OAAO,oCAAoC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AACzD,CAAC;AAED;;;GAGG;AACH,SAAS,mBAAmB,CAAC,IAAyC;IACpE,OAAO,CAAC,CAAC,GAAG,IAAI,IAAI,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;AAClD,CAAC;AAED,SAAS,iBAAiB,CAAC,GAAc;IACvC,OAAO,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,cAAc,CAAC,CAAC;AAC5D,CAAC;AAED;;;;;GAKG;AACH,MAAM,CAAC,MAAM,6BAA6B,GACxC,uOAAuO,CAAC;AAE1O,wEAAwE;AACxE,SAAS,WAAW,CAAC,IAAc;IACjC,MAAM,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzB,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,iBAAiB,KAAK,IAAI,EAAE,CAAC;QAC1C,MAAM,IAAI,qBAAqB,CAC7B,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,IAAI,IAAI,4CAA4C,CACtE,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC,iBAAiB,CAAC;AAC7B,CAAC;AAED;0EAC0E;AAC1E,SAAS,kBAAkB,CAAC,EAAU,EAAE,WAAqB;IAC3D,MAAM,MAAM,GAAG,WAAW,CAAC,WAAW,CAAC,CAAC;IACxC,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,aAAa,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,qBAAqB,CAC7B,GAAG,EAAE,6CAA6C,MAAM,CAAC,MAAM,CAAC,EAAE,CACnE,CAAC;IACJ,CAAC;IACD,OAAO,GAAG,CAAC,GAAG,EAAE,OAAO,EAAE,aAAa,CAAC,CAAC;AAC1C,CAAC;AAgCD,4EAA4E;AAC5E,SAAS,aAAa,CAAC,IAAc;IACnC,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,cAAc;YACjB,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACpD,KAAK,YAAY;YACf,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,IAAK,EAAE,CAAC;QAC3C,KAAK,iBAAiB;YACpB,OAAO,aAAa,CAAC,IAAI,CAAC,OAAQ,CAAC,CAAC;QACtC,KAAK,cAAc,EAAE,CAAC;YACpB,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;YAClD,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC;gBAAE,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;YAC/E,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG;gBAAE,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;YACzE,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,EAAE,CAAC;gBAC7E,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,EAAE,CAAC;YAC3C,CAAC;YACD,MAAM,IAAI,sBAAsB,CAAC,IAAI,CAAC,IAAI,EAAE,yBAAyB,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACnF,CAAC;QACD,KAAK,cAAc;YACjB,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,IAAI,IAAI,CAAC,IAAI,IAAI,EAAE,EAAE,IAAI,EAAE,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,aAAa,CAAC,EAAE,CAAC;QAC5G;YACE,MAAM,IAAI,sBAAsB,CAAC,IAAI,CAAC,IAAI,EAAE,qBAAqB,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC;IACnF,CAAC;AACH,CAAC;AAED,yEAAyE;AACzE,SAAS,aAAa,CAAC,IAAW;IAChC,MAAM,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,CAAC;IAC/B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,qBAAqB,CAAC,4CAA4C,CAAC,CAAC;IACvG,OAAO,CAAC,CAAC;AACX,CAAC;AACD,SAAS,cAAc,CAAC,IAAW;IACjC,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,KAAK,CAAC;QACpB,KAAK,KAAK;YACR,OAAO,CAAC,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACnC,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QACvE,KAAK,IAAI,EAAE,CAAC;YACV,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,CAAC;YACxC,QAAQ,IAAI,CAAC,EAAE,EAAE,CAAC;gBAChB,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC9C,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACpE,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC9C,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC7C,CAAC;QACH,CAAC;QACD,6DAA6D;QAC7D;YACE,MAAM,IAAI,qBAAqB,CAAC,qCAAqC,CAAC,CAAC;IAC3E,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,SAAS,KAAK,CAAC,IAAc,EAAE,QAAe,EAAE,IAAyC;IACvF,IAAI,QAA4B,CAAC;IACjC,IAAI,CAAC;QACH,QAAQ,GAAG,aAAa,CAAC,QAAQ,CAAC,CAAC;IACrC,CAAC;IAAC,OAAO,KAAK,EAAE,CAAC;QACf,IAAI,CAAC,CAAC,KAAK,YAAY,qBAAqB,CAAC;YAAE,MAAM,KAAK,CAAC;IAC7D,CAAC;IACD,IAAI,QAAQ,KAAK,SAAS;QAAE,OAAO,OAAO,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;IAC3D,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,UAAU,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;AACrD,CAAC;AAED,0EAA0E;AAC1E,SAAS,UAAU,CAAC,IAAW,EAAE,IAAyC;IACxE,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,aAAa,CAAC,CAAC;QAChD,KAAK,KAAK;YACR,OAAO,aAAa,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACxC,KAAK,KAAK;YACR,OAAO,UAAU,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,4BAA4B;QACjE,KAAK,IAAI;YACP,OAAO,EAAE,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAChE,KAAK,KAAK;YACR,OAAO,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QAC5D,KAAK,MAAM,EAAE,CAAC;YACZ,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;YACnB,IAAI,EAAE,KAAK,MAAM;gBAAE,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;YACvE,IAAI,EAAE,KAAK,MAAM;gBAAE,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACzE,IAAI,EAAE,KAAK,KAAK;gBAAE,OAAO,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC;YACnF,MAAM,MAAM,GAAG,iBAAiB,CAAC,EAAE,EAAE,UAAU,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;YACrE,IAAI,MAAM;gBAAE,OAAO,MAAM,CAAC;YAC1B,MAAM,IAAI,qBAAqB,CAAC,yBAAyB,EAAE,GAAG,CAAC,CAAC;QAClE,CAAC;IACH,CAAC;AACH,CAAC;AAgCD;4EAC4E;AAC5E,SAAS,6BAA6B,CACpC,UAAiF;IAEjF,MAAM,KAAK,GAAG,CAAC,IAAY,EAAE,IAAyC,EAAiB,EAAE;QACvF,IAAI,QAAkB,CAAC;QACvB,IAAI,CAAC;YACH,QAAQ,GAAG,UAAU,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACpC,CAAC;QAAC,OAAO,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,YAAY,qBAAqB;gBAAE,MAAM,CAAC,CAAC;YAChD,MAAM,IAAI,qBAAqB,CAAC,gBAAgB,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QAChG,CAAC;QACD,MAAM,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;QAC7B,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,iBAAiB,KAAK,IAAI,EAAE,CAAC;YAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,IAAI,IAAI,+BAA+B,CAAC;YACtE,MAAM,WAAW,GACf,eAAe,CAAC,IAAI,CAAC;gBACrB,mBAAmB,CAAC,IAAI,CAAC;gBACzB,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,iBAAiB,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,iBAAiB,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;YACzF,MAAM,WAAW,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,qBAAqB,CAAC,IAAI,CAAC,CAAC;QAC3F,CAAC;QACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,SAAS,EAAE,CAAC,CAAC,iBAAiB,EAAE,CAAC;IAC5D,CAAC,CAAC;IACF,OAAO;QACL,KAAK;QACL,KAAK,CAAC,IAAI,EAAE,IAAI;YACd,IAAI,CAAC;gBACH,OAAO,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,SAAS,EAAE,CAAC;YACxD,CAAC;YAAC,OAAO,CAAC,EAAE,CAAC;gBACX,OAAO;oBACL,EAAE,EAAE,KAAK;oBACT,KAAK,EAAE,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;oBACjD,GAAG,CAAC,CAAC,YAAY,wBAAwB,CAAC,CAAC,CAAC,EAAE,qBAAqB,EAAE,IAAa,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;oBAC1F,GAAG,CAAC,CAAC,YAAY,qBAAqB,CAAC,CAAC,CAAC,EAAE,gBAAgB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;iBAC9E,CAAC;YACJ,CAAC;QACH,CAAC;KACF,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CAAC,IAAY;IAClC,IAAI,CAAC;QACH,OAAO,WAAW,CAAC,IAAI,CAAwB,CAAC;IAClD,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,YAAY,CAAC,gBAAgB,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACvF,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,iBAAiB,CAAC,IAAY;IAC5C,MAAM,IAAI,GAAG,aAAa,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC;IACjD,WAAW,CAAC,IAAI,CAAC,CAAC;IAClB,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,WAAW,CAAC,IAAW;IAC9B,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO;gBAAE,MAAM,IAAI,gBAAgB,EAAE,CAAC;YACxD,OAAO;QACT,KAAK,KAAK;YACR,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACtB,OAAO;QACT,KAAK,IAAI;YACP,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI;gBAAE,WAAW,CAAC,CAAC,CAAC,CAAC;YAC1C,OAAO;QACT,KAAK,KAAK;YACR,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YACvB,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACtB,OAAO;QACT,KAAK,MAAM;YACT,IAAI,IAAI,CAAC,EAAE,KAAK,OAAO;gBAAE,MAAM,IAAI,gBAAgB,EAAE,CAAC;YACtD,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI;gBAAE,WAAW,CAAC,CAAC,CAAC,CAAC;YAC1C,OAAO;QACT,KAAK,KAAK;YACR,OAAO;IACX,CAAC;AACH,CAAC;AAED,8EAA8E;AAC9E,MAAM,UAAU,8BAA8B;IAC5C,OAAO,6BAA6B,CAAC,CAAC,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC,UAAU,CAAC,aAAa,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;AAC9G,CAAC"}
+{"version":3,"file":"formula-dimension.js","sourceRoot":"","sources":["../../src/numerical/formula-dimension.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,EAAE,KAAK,IAAI,WAAW,EAAE,MAAM,iCAAiC,CAAC;AAEvE,OAAO,EAAE,MAAM,EAAE,aAAa,EAAE,MAAM,EAAE,MAAM,yBAAyB,CAAC;AACxE,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,2BAA2B,CAAC;AAErE,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AACvD,OAAO,EAAE,GAAG,EAAE,MAAM,gCAAgC,CAAC;AACrD,OAAO,EAAE,kBAAkB,EAAE,gBAAgB,EAAE,YAAY,EAAE,gBAAgB,EAAE,sBAAsB,EAAE,MAAM,uBAAuB,CAAC;AAErI;;cAEc;AACd,MAAM,OAAO,qBAAsB,SAAQ,KAAK;IAC9C,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;IACtC,CAAC;CACF;AAED,6GAA6G;AAC7G,MAAM,OAAO,sBAAuB,SAAQ,qBAAqB;IAC1C,QAAQ;IAA7B,YAAqB,QAAgB,EAAE,MAAc;QACnD,KAAK,CAAC,MAAM,CAAC,CAAC;wBADK,QAAQ;QAE3B,IAAI,CAAC,IAAI,GAAG,wBAAwB,CAAC;IACvC,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,qBAAsB,SAAQ,qBAAqB;IAClC,MAAM;IAA3B,YAAqB,MAAc;QACjC,KAAK,CAAC,sBAAsB,MAAM,sFAAsF,CAAC,CAAC;sBADvG,MAAM;QAEzB,IAAI,CAAC,IAAI,GAAG,uBAAuB,CAAC;IACtC,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,wBAAyB,SAAQ,qBAAqB;IAC1D,YAAY,IAAY;QACtB,KAAK,CAAC,GAAG,6BAA6B,IAAI,IAAI,EAAE,CAAC,CAAC;QAClD,IAAI,CAAC,IAAI,GAAG,0BAA0B,CAAC;IACzC,CAAC;CACF;AAED;qFACqF;AACrF,MAAM,cAAc,GAAG,IAAI,GAAG,CAAC,CAAC,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,UAAU,EAAE,KAAK,CAAC,CAAC,CAAC;AAExE;;;;;GAKG;AACH,MAAM,UAAU,sBAAsB,CAAC,IAAY;IACjD,IAAI,IAAI,KAAK,GAAG;QAAE,OAAO,MAAM,CAAC;IAChC,IAAI,IAAI,KAAK,GAAG;QAAE,OAAO,MAAM,CAAC;IAChC,IAAI,IAAI,KAAK,OAAO;QAAE,OAAO,SAAS,CAAC;IACvC,IAAI,cAAc,CAAC,GAAG,CAAC,IAAI,CAAC;QAAE,OAAO,aAAa,CAAC;IACnD,OAAO,SAAS,CAAC;AACnB,CAAC;AAED,MAAM,EAAE,GAAG,CAAC,CAA8B,EAAE,IAAgB,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC;AACzG,MAAM,OAAO,GAAG,CAAC,IAAc,EAAE,GAAW,EAAY,EAAE,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC;AAC5G,MAAM,cAAc,GAAG,CAAC,EAAoB,EAAE,GAAa,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,gBAAgB,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;AAChH,MAAM,OAAO,GAAG,CAAC,GAAa,EAAY,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,CAAC,CAAC;AAEpE;mDACmD;AACnD,SAAS,aAAa,CAAC,IAAY,EAAE,IAAyC;IAC5E,IAAI,IAAI,KAAK,OAAO,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,CAAC;QAAE,MAAM,IAAI,qBAAqB,CAAC,kBAAkB,CAAC,CAAC;IAC7F,IAAI,IAAI,IAAI,IAAI;QAAE,OAAO,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;IAC/C,MAAM,KAAK,GAAG,sBAAsB,CAAC,IAAI,CAAC,CAAC;IAC3C,IAAI,KAAK,KAAK,SAAS;QAAE,OAAO,GAAG,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IACjD,MAAM,IAAI,qBAAqB,CAAC,IAAI,CAAC,CAAC;AACxC,CAAC;AAED,MAAM,cAAc,GAAG,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;AAEhD,6EAA6E;AAC7E,SAAS,eAAe,CAAC,IAAY;IACnC,OAAO,oCAAoC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AACzD,CAAC;AAED;;;GAGG;AACH,SAAS,mBAAmB,CAAC,IAAyC;IACpE,OAAO,CAAC,CAAC,GAAG,IAAI,IAAI,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC;AAClD,CAAC;AAED,SAAS,iBAAiB,CAAC,GAAc;IACvC,OAAO,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,cAAc,CAAC,CAAC;AAC5D,CAAC;AAED;;;;GAIG;AACH,MAAM,6BAA6B,GACjC,uOAAuO,CAAC;AAE1O,wEAAwE;AACxE,SAAS,WAAW,CAAC,IAAc;IACjC,MAAM,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,CAAC;IACzB,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,iBAAiB,KAAK,IAAI,EAAE,CAAC;QAC1C,MAAM,IAAI,qBAAqB,CAC7B,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,IAAI,IAAI,4CAA4C,CACtE,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC,iBAAiB,CAAC;AAC7B,CAAC;AAED;0EAC0E;AAC1E,SAAS,kBAAkB,CAAC,EAAU,EAAE,WAAqB;IAC3D,MAAM,MAAM,GAAG,WAAW,CAAC,WAAW,CAAC,CAAC;IACxC,IAAI,CAAC,MAAM,CAAC,MAAM,EAAE,aAAa,CAAC,EAAE,CAAC;QACnC,MAAM,IAAI,qBAAqB,CAC7B,GAAG,EAAE,6CAA6C,MAAM,CAAC,MAAM,CAAC,EAAE,CACnE,CAAC;IACJ,CAAC;IACD,OAAO,GAAG,CAAC,GAAG,EAAE,OAAO,EAAE,aAAa,CAAC,CAAC;AAC1C,CAAC;AAED;;;;GAIG;AACH,SAAS,UAAU,CAAC,EAAU,EAAE,IAAsB,EAAE,IAAyC;IAC/F,MAAM,IAAI,GAAG,gBAAgB,CAAC,EAAE,CAAC,CAAC;IAClC,IAAI,IAAI,KAAK,SAAS;QAAE,MAAM,IAAI,qBAAqB,CAAC,sBAAsB,CAAC,EAAE,CAAC,CAAC,CAAC;IACpF,IAAI,IAAI,CAAC,MAAM,KAAK,IAAI,CAAC,KAAK,EAAE,CAAC;QAC/B,MAAM,IAAI,qBAAqB,CAC7B,GAAG,EAAE,cAAc,IAAI,CAAC,KAAK,YAAY,IAAI,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,SAAS,IAAI,CAAC,MAAM,EAAE,CAC3F,CAAC;IACJ,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,CAAC,SAAS,CAAC;IAC5B,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,MAAM;YACT,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QACzD,KAAK,OAAO;YACV,OAAO,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,CAAC;QAC3D,KAAK,MAAM;YACT,OAAO,OAAO,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,CAAC,CAAC;QAC7C,KAAK,gBAAgB,EAAE,CAAC;YACtB,MAAM,GAAG,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,CAAC;YACvC,OAAO,IAAI,CAAC,IAAI,KAAK,SAAS,CAAC,CAAC,CAAC,kBAAkB,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,cAAc,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;QAChG,CAAC;QACD,KAAK,OAAO,EAAE,CAAC;YACb,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,WAAW,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,CAAC,EAAE,WAAW,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAE,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;YAClG,IAAI,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;gBAClB,MAAM,IAAI,qBAAqB,CAAC,GAAG,EAAE,kDAAkD,MAAM,CAAC,CAAC,CAAC,QAAQ,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACvH,CAAC;YACD,OAAO,GAAG,CAAC,GAAG,EAAE,OAAO,EAAE,aAAa,CAAC,CAAC;QAC1C,CAAC;IACH,CAAC;AACH,CAAC;AAgCD,4EAA4E;AAC5E,SAAS,aAAa,CAAC,IAAc;IACnC,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,cAAc;YACjB,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;QACpD,KAAK,YAAY;YACf,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,IAAK,EAAE,CAAC;QAC3C,KAAK,iBAAiB;YACpB,OAAO,aAAa,CAAC,IAAI,CAAC,OAAQ,CAAC,CAAC;QACtC,KAAK,cAAc,EAAE,CAAC;YACpB,MAAM,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,aAAa,CAAC,CAAC;YAClD,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC;gBAAE,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;YAC/E,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG;gBAAE,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;YACzE,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,EAAE,CAAC;gBAC7E,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,EAAE,CAAC;YAC3C,CAAC;YACD,MAAM,IAAI,sBAAsB,CAAC,IAAI,CAAC,IAAI,EAAE,yBAAyB,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACnF,CAAC;QACD,KAAK,cAAc;YACjB,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,EAAE,IAAI,IAAI,IAAI,CAAC,IAAI,IAAI,EAAE,EAAE,IAAI,EAAE,CAAC,IAAI,CAAC,IAAI,IAAI,EAAE,CAAC,CAAC,GAAG,CAAC,aAAa,CAAC,EAAE,CAAC;QAC5G;YACE,MAAM,IAAI,sBAAsB,CAAC,IAAI,CAAC,IAAI,EAAE,qBAAqB,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC;IACnF,CAAC;AACH,CAAC;AAED,yEAAyE;AACzE,SAAS,aAAa,CAAC,IAAW;IAChC,MAAM,CAAC,GAAG,cAAc,CAAC,IAAI,CAAC,CAAC;IAC/B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAAE,MAAM,IAAI,qBAAqB,CAAC,4CAA4C,CAAC,CAAC;IACvG,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,gBAAgB,CAAC,IAAW;IAC1C,IAAI,CAAC;QACH,OAAO,aAAa,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IAAC,OAAO,KAAK,EAAE,CAAC;QACf,IAAI,KAAK,YAAY,qBAAqB;YAAE,OAAO,SAAS,CAAC;QAC7D,MAAM,KAAK,CAAC;IACd,CAAC;AACH,CAAC;AAED,SAAS,cAAc,CAAC,IAAW;IACjC,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,KAAK,CAAC;QACpB,KAAK,KAAK;YACR,OAAO,CAAC,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACnC,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,cAAc,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;QACvE,KAAK,IAAI,EAAE,CAAC;YACV,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,cAAc,CAAC,CAAC;YACxC,QAAQ,IAAI,CAAC,EAAE,EAAE,CAAC;gBAChB,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC9C,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;gBACpE,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;gBAC9C,KAAK,GAAG,EAAE,OAAO,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC7C,CAAC;QACH,CAAC;QACD,6DAA6D;QAC7D;YACE,MAAM,IAAI,qBAAqB,CAAC,qCAAqC,CAAC,CAAC;IAC3E,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,SAAS,KAAK,CAAC,IAAc,EAAE,QAAe,EAAE,IAAyC;IACvF,IAAI,QAA4B,CAAC;IACjC,IAAI,CAAC;QACH,QAAQ,GAAG,aAAa,CAAC,QAAQ,CAAC,CAAC;IACrC,CAAC;IAAC,OAAO,KAAK,EAAE,CAAC;QACf,IAAI,CAAC,CAAC,KAAK,YAAY,qBAAqB,CAAC;YAAE,MAAM,KAAK,CAAC;IAC7D,CAAC;IACD,IAAI,QAAQ,KAAK,SAAS;QAAE,OAAO,OAAO,CAAC,IAAI,EAAE,QAAQ,CAAC,CAAC;IAC3D,OAAO,EAAE,CAAC,GAAG,EAAE,CAAC,IAAI,EAAE,UAAU,CAAC,QAAQ,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;AACrD,CAAC;AAED,0EAA0E;AAC1E,SAAS,UAAU,CAAC,IAAW,EAAE,IAAyC;IACxE,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,aAAa,CAAC,CAAC;QAChD,KAAK,KAAK;YACR,OAAO,aAAa,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACxC,KAAK,KAAK;YACR,OAAO,UAAU,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,4BAA4B;QACjE,KAAK,IAAI;YACP,OAAO,EAAE,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;QAChE,KAAK,KAAK;YACR,OAAO,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QAC5D,KAAK,MAAM;YACT,OAAO,UAAU,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;IAChD,CAAC;AACH,CAAC;AAgCD;4EAC4E;AAC5E,SAAS,6BAA6B,CACpC,UAAiF;IAEjF,MAAM,KAAK,GAAG,CAAC,IAAY,EAAE,IAAyC,EAAiB,EAAE;QACvF,IAAI,QAAkB,CAAC;QACvB,IAAI,CAAC;YACH,QAAQ,GAAG,UAAU,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC;QACpC,CAAC;QAAC,OAAO,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,YAAY,qBAAqB;gBAAE,MAAM,CAAC,CAAC;YAChD,MAAM,IAAI,qBAAqB,CAAC,gBAAgB,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QAChG,CAAC;QACD,MAAM,CAAC,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;QAC7B,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,iBAAiB,KAAK,IAAI,EAAE,CAAC;YAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,IAAI,IAAI,+BAA+B,CAAC;YACtE,MAAM,WAAW,GACf,eAAe,CAAC,IAAI,CAAC;gBACrB,mBAAmB,CAAC,IAAI,CAAC;gBACzB,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,iBAAiB,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,iBAAiB,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;YACzF,MAAM,WAAW,CAAC,CAAC,CAAC,IAAI,wBAAwB,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,qBAAqB,CAAC,IAAI,CAAC,CAAC;QAC3F,CAAC;QACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,SAAS,EAAE,CAAC,CAAC,iBAAiB,EAAE,CAAC;IAC5D,CAAC,CAAC;IACF,OAAO;QACL,KAAK;QACL,KAAK,CAAC,IAAI,EAAE,IAAI;YACd,IAAI,CAAC;gBACH,OAAO,EAAE,EAAE,EAAE,IAAI,EAAE,GAAG,EAAE,KAAK,CAAC,IAAI,EAAE,IAAI,CAAC,CAAC,SAAS,EAAE,CAAC;YACxD,CAAC;YAAC,OAAO,CAAC,EAAE,CAAC;gBACX,OAAO;oBACL,EAAE,EAAE,KAAK;oBACT,KAAK,EAAE,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;oBACjD,GAAG,CAAC,CAAC,YAAY,wBAAwB,CAAC,CAAC,CAAC,EAAE,qBAAqB,EAAE,IAAa,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;oBAC1F,GAAG,CAAC,CAAC,YAAY,qBAAqB,CAAC,CAAC,CAAC,EAAE,gBAAgB,EAAE,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;iBAC9E,CAAC;YACJ,CAAC;QACH,CAAC;KACF,CAAC;AACJ,CAAC;AAED,SAAS,cAAc,CAAC,IAAY;IAClC,IAAI,CAAC;QACH,0EAA0E;QAC1E,wEAAwE;QACxE,sEAAsE;QACtE,OAAO,WAAW,CAAC,IAAI,CAAwB,CAAC;IAClD,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,MAAM,IAAI,YAAY,CAAC,gBAAgB,CAAC,YAAY,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACvF,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,iBAAiB,CAAC,IAAY;IAC5C,MAAM,IAAI,GAAG,aAAa,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC;IACjD,WAAW,CAAC,IAAI,CAAC,CAAC;IAClB,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,WAAW,CAAC,IAAW;IAC9B,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,IAAI,IAAI,CAAC,IAAI,KAAK,OAAO;gBAAE,MAAM,IAAI,gBAAgB,EAAE,CAAC;YACxD,OAAO;QACT,KAAK,KAAK;YACR,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACtB,OAAO;QACT,KAAK,IAAI;YACP,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI;gBAAE,WAAW,CAAC,CAAC,CAAC,CAAC;YAC1C,OAAO;QACT,KAAK,KAAK;YACR,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YACvB,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YACtB,OAAO;QACT,KAAK,MAAM;YACT,IAAI,IAAI,CAAC,EAAE,KAAK,OAAO;gBAAE,MAAM,IAAI,gBAAgB,EAAE,CAAC;YACtD,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI;gBAAE,WAAW,CAAC,CAAC,CAAC,CAAC;YAC1C,OAAO;QACT,KAAK,KAAK;YACR,OAAO;IACX,CAAC;AACH,CAAC;AAED,8EAA8E;AAC9E,MAAM,UAAU,8BAA8B;IAC5C,OAAO,6BAA6B,CAAC,CAAC,IAAI,EAAE,IAAI,EAAE,EAAE,CAAC,UAAU,CAAC,aAAa,CAAC,cAAc,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;AAC9G,CAAC","sourcesContent":["/**\n * Formula dimensional check (MathTS Phase 2 — see\n * docs/planning/Formula-Dimensional-Check-Design-Note.md).\n *\n * Transpiles a parsed formula AST into UPT's own dimensional `ExprNode` and\n * runs `validate()`, so the CLI can report whether a user's formula is\n * dimensionally HOMOGENEOUS and what dimension it has — unifying string→AST\n * with AST→dimension.\n *\n * The MathTS AST is the only front-end. It adapts to a normalized `PNode`,\n * and the dimensional transpilation lives on that node.\n *\n * @module numerical/formula-dimension\n */\n\nimport { parse as parseMathTs } from '@danielsimonjr/mathts-functions';\nimport type { Dimension } from '../dimensional/types.js';\nimport { CHARGE, DIMENSIONLESS, ENERGY } from '../dimensional/types.js';\nimport { equals, format, multiply } from '../dimensional/algebra.js';\nimport type { ExprNode, TranscendentalFn } from '../dimensional/validator.js';\nimport { validate } from '../dimensional/validator.js';\nimport { sym } from '../dimensional/ast-builders.js';\nimport { EULER_NUMBER_ERROR, EulerNumberError, FormulaError, SCALAR_FUNCTIONS, unknownFunctionMessage } from './formula-contract.js';\n\n/** A formula cannot be dimensionally analyzed (undeclared symbol, variable\n *  exponent, transcendental of a dimensional argument, unsupported node).\n *  @public */\nexport class FormulaDimensionError extends Error {\n  constructor(message: string) {\n    super(message);\n    this.name = 'FormulaDimensionError';\n  }\n}\n\n/** The parse tree holds a node the normalized form has no shape for (`a.b`, `[1, 2]`, `x > 2`). @internal */\nexport class UnsupportedSyntaxError extends FormulaDimensionError {\n  constructor(readonly nodeType: string, detail: string) {\n    super(detail);\n    this.name = 'UnsupportedSyntaxError';\n  }\n}\n\n/**\n * A formula names a symbol with no declared dimension and no known one. A\n * caller reads it through the check result's `undeclaredSymbol`. @internal\n */\nclass UndeclaredSymbolError extends FormulaDimensionError {\n  constructor(readonly symbol: string) {\n    super(`undeclared symbol '${symbol}' — declare its dimension (a separate name:dimension argument, for example x:length)`);\n    this.name = 'UndeclaredSymbolError';\n  }\n}\n\n/**\n * A sum is not homogeneous because a bare `e`, read as the elementary charge,\n * meets a dimensionless term (`1 - e^2`). @internal\n */\nclass ElementaryChargeMixError extends FormulaDimensionError {\n  constructor(note: string) {\n    super(`${ELEMENTARY_CHARGE_MIX_MESSAGE} ${note}`);\n    this.name = 'ElementaryChargeMixError';\n  }\n}\n\n/** Dimensionless math constants the parser recognizes. Bare `e` is not here:\n *  ISO 80000 names that symbol the elementary charge. Euler's number is `exp(x)`. */\nconst MATH_CONSTANTS = new Set(['pi', 'tau', 'phi', 'Infinity', 'NaN']);\n\n/**\n * Dimension of a symbol the physics parser knows when the caller did not\n * declare one. A caller-supplied dimension still wins. `e` is the elementary\n * charge. `E` is energy. The name `euler` is not a constant.\n * @internal\n */\nexport function formulaSymbolDimension(name: string): Dimension | undefined {\n  if (name === 'e') return CHARGE;\n  if (name === 'E') return ENERGY;\n  if (name === 'euler') return undefined;\n  if (MATH_CONSTANTS.has(name)) return DIMENSIONLESS;\n  return undefined;\n}\n\nconst op = (o: '+' | '-' | '*' | '/' | '^', args: ExprNode[]): ExprNode => ({ kind: 'op', op: o, args });\nconst powExpr = (base: ExprNode, exp: number): ExprNode => op('^', [base, sym(String(exp), DIMENSIONLESS)]);\nconst transcendental = (fn: TranscendentalFn, arg: ExprNode): ExprNode => ({ kind: 'transcendental', fn, arg });\nconst absNode = (arg: ExprNode): ExprNode => ({ kind: 'abs', arg });\n\n/** Resolve a symbol to a dimensioned `ExprNode` (declared dim, or a\n *  dimensionless math constant, else an error). */\nfunction resolveSymbol(name: string, dims: Readonly<Record<string, Dimension>>): ExprNode {\n  if (name === 'euler' && !(name in dims)) throw new FormulaDimensionError(EULER_NUMBER_ERROR);\n  if (name in dims) return sym(name, dims[name]);\n  const known = formulaSymbolDimension(name);\n  if (known !== undefined) return sym(name, known);\n  throw new UndeclaredSymbolError(name);\n}\n\nconst CHARGE_SQUARED = multiply(CHARGE, CHARGE);\n\n/** A bare `e`, not `exp`, `e_charge`, `eps0`, or the exponent in `1e-19`. */\nfunction formulaHasBareE(expr: string): boolean {\n  return /(?<![A-Za-z0-9_])e(?![A-Za-z0-9_])/.test(expr);\n}\n\n/**\n * `e` is the elementary charge when the caller did not give it another\n * dimension. A declared dimensionless `e` is not this case.\n */\nfunction eIsElementaryCharge(dims: Readonly<Record<string, Dimension>>): boolean {\n  return !('e' in dims) || equals(dims.e, CHARGE);\n}\n\nfunction dimensionIsCharge(dim: Dimension): boolean {\n  return equals(dim, CHARGE) || equals(dim, CHARGE_SQUARED);\n}\n\n/**\n * What to say when a bare `e` was the elementary charge inside a sum that is\n * not homogeneous. Declaring or binding `e` keeps a different quantity.\n * Euler's number is `exp(x)`.\n */\nconst ELEMENTARY_CHARGE_MIX_MESSAGE =\n  'e is the elementary charge and is not dimensionless here. Declare it or bind a value (e=<number>) if you mean a different quantity, or write Euler\\'s number as exp(x). The catalog writes the eccentricity factor as one_minus_e_sq.';\n\n/** Inferred dimension of an `ExprNode`, or throw if not homogeneous. */\nfunction dimensionOf(node: ExprNode): Dimension {\n  const r = validate(node);\n  if (!r.ok || r.inferredDimension === null) {\n    throw new FormulaDimensionError(\n      r.violations[0]?.note ?? 'expression is not dimensionally consistent',\n    );\n  }\n  return r.inferredDimension;\n}\n\n/** Transcendental typed-stub: assert the argument is dimensionless, collapse\n *  to a dimensionless symbol (keeps the `ExprNode` grammar unchanged). */\nfunction transcendentalStub(fn: string, argExprNode: ExprNode): ExprNode {\n  const argDim = dimensionOf(argExprNode);\n  if (!equals(argDim, DIMENSIONLESS)) {\n    throw new FormulaDimensionError(\n      `${fn}() requires a dimensionless argument, got ${format(argDim)}`,\n    );\n  }\n  return sym(`${fn}(...)`, DIMENSIONLESS);\n}\n\n/**\n * A call, dispatched by the one function table's dimension rule. The arity is\n * the table's: a call with another number of arguments is a dimension error,\n * not a crash on a missing argument or a silently dropped extra one.\n */\nfunction callToExpr(fn: string, args: readonly PNode[], dims: Readonly<Record<string, Dimension>>): ExprNode {\n  const spec = SCALAR_FUNCTIONS[fn];\n  if (spec === undefined) throw new FormulaDimensionError(unknownFunctionMessage(fn));\n  if (args.length !== spec.arity) {\n    throw new FormulaDimensionError(\n      `${fn}() expects ${spec.arity} argument${spec.arity === 1 ? '' : 's'}, got ${args.length}`,\n    );\n  }\n  const rule = spec.dimension;\n  switch (rule.kind) {\n    case 'root':\n      return powExpr(normToExpr(args[0]!, dims), rule.power);\n    case 'power':\n      return powOf(normToExpr(args[0]!, dims), args[1]!, dims);\n    case 'same':\n      return absNode(normToExpr(args[0]!, dims));\n    case 'transcendental': {\n      const arg = normToExpr(args[0]!, dims);\n      return rule.node === undefined ? transcendentalStub(fn, arg) : transcendental(rule.node, arg);\n    }\n    case 'ratio': {\n      const [y, x] = [dimensionOf(normToExpr(args[0]!, dims)), dimensionOf(normToExpr(args[1]!, dims))];\n      if (!equals(y, x)) {\n        throw new FormulaDimensionError(`${fn}() arguments must have the same dimension, got ${format(y)} and ${format(x)}`);\n      }\n      return sym(`${fn}(...)`, DIMENSIONLESS);\n    }\n  }\n}\n\n// --- one transpiler over a normalized parse node (Phase 2) ----------------\n//\n// The MathTS AST, the one front-end, is adapted to a tiny normalized `PNode`;\n// the dimensional transpilation (`normToExpr`) and the constant-exponent\n// folding (`pnodeConstant`) live in ONE place. `ExprNode` is the single\n// semantic IR; the parse-tree is transient.\n\n/** Structural shape of a MathTS AST node (the bits the adapter reads). */\ninterface MathNode {\n  readonly type: string;\n  readonly value?: number;\n  readonly name?: string;\n  readonly op?: string;\n  readonly args?: MathNode[];\n  readonly content?: MathNode;\n  readonly fn?: { readonly name?: string };\n}\n\n/** Normalized parse node — the shape the MathTS AST adapts to. @internal */\ntype PNode =\n  | { kind: 'num'; value: number }\n  | { kind: 'sym'; name: string }\n  | { kind: 'neg'; arg: PNode } // unary minus (dimension-neutral)\n  | { kind: 'op'; op: '+' | '-' | '*' | '/'; args: PNode[] }\n  | { kind: 'pow'; base: PNode; exp: PNode } // exp must fold to a constant\n  | { kind: 'call'; fn: string; args: PNode[] };\n\n/** The normalized node {@link parseFormulaPNode} returns. @internal */\nexport type FormulaPNode = PNode;\n\n/** Adapt a MathTS AST node to a `PNode` (pure shape map; no dimensions). */\nfunction mathtsToPNode(node: MathNode): PNode {\n  switch (node.type) {\n    case 'ConstantNode':\n      return { kind: 'num', value: Number(node.value) };\n    case 'SymbolNode':\n      return { kind: 'sym', name: node.name! };\n    case 'ParenthesisNode':\n      return mathtsToPNode(node.content!);\n    case 'OperatorNode': {\n      const args = (node.args ?? []).map(mathtsToPNode);\n      if (node.op === '-' && args.length === 1) return { kind: 'neg', arg: args[0] };\n      if (node.op === '^') return { kind: 'pow', base: args[0], exp: args[1] };\n      if (node.op === '+' || node.op === '-' || node.op === '*' || node.op === '/') {\n        return { kind: 'op', op: node.op, args };\n      }\n      throw new UnsupportedSyntaxError(node.type, `unsupported operator '${node.op}'`);\n    }\n    case 'FunctionNode':\n      return { kind: 'call', fn: node.fn?.name ?? node.name ?? '', args: (node.args ?? []).map(mathtsToPNode) };\n    default:\n      throw new UnsupportedSyntaxError(node.type, `unsupported node '${node.type}'`);\n  }\n}\n\n/** Fold a (presumed constant) `PNode` to a finite number, else throw. */\nfunction pnodeConstant(node: PNode): number {\n  const v = evalConstPNode(node);\n  if (!Number.isFinite(v)) throw new FormulaDimensionError('exponent must be a finite numeric constant');\n  return v;\n}\n\n/**\n * The number a numeric-constant exponent folds to (`2`, `-1`, `1/2`,\n * `(2+1)/3`), or undefined when the node names a symbol or a call. The\n * binding reader uses the same folding for its exponents, so the two\n * readers agree on what an exponent may be.\n * @internal\n */\nexport function constantExponent(node: PNode): number | undefined {\n  try {\n    return pnodeConstant(node);\n  } catch (error) {\n    if (error instanceof FormulaDimensionError) return undefined;\n    throw error;\n  }\n}\n\nfunction evalConstPNode(node: PNode): number {\n  switch (node.kind) {\n    case 'num':\n      return node.value;\n    case 'neg':\n      return -evalConstPNode(node.arg);\n    case 'pow':\n      return Math.pow(evalConstPNode(node.base), evalConstPNode(node.exp));\n    case 'op': {\n      const v = node.args.map(evalConstPNode);\n      switch (node.op) {\n        case '+': return v.reduce((a, b) => a + b, 0);\n        case '-': return v.length === 1 ? -v[0] : v.reduce((a, b) => a - b);\n        case '*': return v.reduce((a, b) => a * b, 1);\n        case '/': return v.reduce((a, b) => a / b);\n      }\n    }\n    // falls through: a 'sym' or 'call' is not a numeric constant\n    default:\n      throw new FormulaDimensionError('exponent must be a numeric constant');\n  }\n}\n\n/**\n * `base ^ exponent`. A numeric-constant exponent is the dimensional power. An\n * input-dependent exponent is built as a `^` node and the validator decides:\n * it is legal only on a dimensionless base, with a dimensionless exponent.\n */\nfunction powOf(base: ExprNode, exponent: PNode, dims: Readonly<Record<string, Dimension>>): ExprNode {\n  let constant: number | undefined;\n  try {\n    constant = pnodeConstant(exponent);\n  } catch (error) {\n    if (!(error instanceof FormulaDimensionError)) throw error;\n  }\n  if (constant !== undefined) return powExpr(base, constant);\n  return op('^', [base, normToExpr(exponent, dims)]);\n}\n\n/** The ONE transpiler: normalized parse node → dimensional `ExprNode`. */\nfunction normToExpr(node: PNode, dims: Readonly<Record<string, Dimension>>): ExprNode {\n  switch (node.kind) {\n    case 'num':\n      return sym(String(node.value), DIMENSIONLESS);\n    case 'sym':\n      return resolveSymbol(node.name, dims);\n    case 'neg':\n      return normToExpr(node.arg, dims); // sign is dimension-neutral\n    case 'op':\n      return op(node.op, node.args.map((a) => normToExpr(a, dims)));\n    case 'pow':\n      return powOf(normToExpr(node.base, dims), node.exp, dims);\n    case 'call':\n      return callToExpr(node.fn, node.args, dims);\n  }\n}\n\n// --- the shared checker ---------------------------------------------------\n\n/** The result of a formula dimensional check. */\ninterface FormulaDimensionResult {\n  /** True iff the formula is dimensionally homogeneous and well-formed. */\n  readonly ok: boolean;\n  /** The inferred dimension when `ok`. */\n  readonly dim?: Dimension;\n  /** Human-readable reason when `!ok`. */\n  readonly error?: string;\n  /** Set when `!ok` because a bare `e` (the elementary charge) met a dimensionless term. */\n  readonly elementaryChargeMixed?: true;\n  /** Set when `!ok` because the formula names this symbol with no dimension declared or known. */\n  readonly undeclaredSymbol?: string;\n}\n\n/** A parsed physics expression: its dimensional `ExprNode` + inferred dimension. */\nexport interface ParsedPhysics {\n  readonly expr: ExprNode;\n  readonly dimension: Dimension;\n}\n\n/** A bound checker: non-throwing `check`, and throwing `parse` (the ExprNode). */\nexport interface FormulaDimensionChecker {\n  check(expr: string, dims: Readonly<Record<string, Dimension>>): FormulaDimensionResult;\n  /** Parse + transpile + validate → `{ expr, dimension }`.\n   *  @throws {FormulaDimensionError} on a parse error or non-homogeneity. */\n  parse(expr: string, dims: Readonly<Record<string, Dimension>>): ParsedPhysics;\n}\n\n/** Build a checker from an `(expr, dims) → ExprNode` transpile (parse +\n *  transpile, which may throw a parse error or `FormulaDimensionError`). */\nfunction createFormulaDimensionChecker(\n  toExprNode: (expr: string, dims: Readonly<Record<string, Dimension>>) => ExprNode,\n): FormulaDimensionChecker {\n  const parse = (expr: string, dims: Readonly<Record<string, Dimension>>): ParsedPhysics => {\n    let exprNode: ExprNode;\n    try {\n      exprNode = toExprNode(expr, dims);\n    } catch (e) {\n      if (e instanceof FormulaDimensionError) throw e;\n      throw new FormulaDimensionError(`parse error: ${e instanceof Error ? e.message : String(e)}`);\n    }\n    const r = validate(exprNode);\n    if (!r.ok || r.inferredDimension === null) {\n      const note = r.violations[0]?.note ?? 'not dimensionally homogeneous';\n      const chargeMixed =\n        formulaHasBareE(expr) &&\n        eIsElementaryCharge(dims) &&\n        r.violations.some((v) => dimensionIsCharge(v.expected) || dimensionIsCharge(v.actual));\n      throw chargeMixed ? new ElementaryChargeMixError(note) : new FormulaDimensionError(note);\n    }\n    return { expr: exprNode, dimension: r.inferredDimension };\n  };\n  return {\n    parse,\n    check(expr, dims) {\n      try {\n        return { ok: true, dim: parse(expr, dims).dimension };\n      } catch (e) {\n        return {\n          ok: false,\n          error: e instanceof Error ? e.message : String(e),\n          ...(e instanceof ElementaryChargeMixError ? { elementaryChargeMixed: true as const } : {}),\n          ...(e instanceof UndeclaredSymbolError ? { undeclaredSymbol: e.symbol } : {}),\n        };\n      }\n    },\n  };\n}\n\nfunction parsedMathNode(expr: string): MathNode {\n  try {\n    // MathTS's Node type is not imported; `MathNode` is the structural subset\n    // this adapter reads (`type`, `value`, `name`, `op`, `args`, `content`,\n    // `fn`), every field of which exists on a MathTS node of that `type`.\n    return parseMathTs(expr) as unknown as MathNode;\n  } catch (e) {\n    throw new FormulaError(`parse error: ${e instanceof Error ? e.message : String(e)}`);\n  }\n}\n\n/**\n * Parse a formula to the normalized node the binding reader evaluates.\n * `euler` is refused here so a binding never becomes Math.E.\n * @internal\n */\nexport function parseFormulaPNode(expr: string): PNode {\n  const node = mathtsToPNode(parsedMathNode(expr));\n  refuseEuler(node);\n  return node;\n}\n\nfunction refuseEuler(node: PNode): void {\n  switch (node.kind) {\n    case 'sym':\n      if (node.name === 'euler') throw new EulerNumberError();\n      return;\n    case 'neg':\n      refuseEuler(node.arg);\n      return;\n    case 'op':\n      for (const a of node.args) refuseEuler(a);\n      return;\n    case 'pow':\n      refuseEuler(node.base);\n      refuseEuler(node.exp);\n      return;\n    case 'call':\n      if (node.fn === 'euler') throw new EulerNumberError();\n      for (const a of node.args) refuseEuler(a);\n      return;\n    case 'num':\n      return;\n  }\n}\n\n/** The dimensional checker over the MathTS AST. The historical name stays. */\nexport function builtinFormulaDimensionChecker(): FormulaDimensionChecker {\n  return createFormulaDimensionChecker((expr, dims) => normToExpr(mathtsToPNode(parsedMathNode(expr)), dims));\n}\n"]}

package/dist/numerical/formula-mathts.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"formula-mathts.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-mathts.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAGH,OAAO,KAAK,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AAuK5E,2CAA2C;AAC3C,eAAO,MAAM,mBAAmB,EAAE,aAAuD,CAAC;AAE1F,wCAAwC;AACxC,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAE1D"}
+{"version":3,"file":"formula-mathts.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-mathts.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAGH,OAAO,KAAK,EAAE,eAAe,EAAE,aAAa,EAAE,MAAM,uBAAuB,CAAC;AA0K5E,2CAA2C;AAC3C,eAAO,MAAM,mBAAmB,EAAE,aAAuD,CAAC;AAE1F,wCAAwC;AACxC,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAE1D"}

package/dist/numerical/formula-mathts.js

@@ -18,7 +18,7 @@
  * @module numerical/formula-mathts
  */
 import { compileExpr, parse as parseMathTs } from '@danielsimonjr/mathts-functions';
-import { BUILTIN_FUNCTION_NAMES, callBuiltinFunction, EulerNumberError, FormulaError, unknownFunctionMessage } from './formula-contract.js';
+import { BUILTIN_FUNCTION_NAMES, callBuiltinFunction, EulerNumberError, FormulaError, unknownFunctionMessage, FormulaValueError } from './formula-contract.js';
 /** Bare `e` is the SI magnitude, not the coulomb Unit. */
 const PHYSICS_SCALAR = { physics: true, charge: 'scalar' };
 /**
@@ -136,7 +136,7 @@ function createMathtsFormulaParser(mod) {
                     if (typeof result !== 'number' || !Number.isFinite(result)) {
                         // typeof Infinity and NaN is "number", which reads as a type error.
                         // A Complex is an object; "got object" hides that sqrt(-1) left the reals.
-                        throw new FormulaError(`formula did not evaluate to a finite number (got ${describeNonFinite(result)})`);
+                        throw new FormulaValueError(`formula did not evaluate to a finite number (got ${describeNonFinite(result)})`);
                     }
                     return result;
                 },
@@ -158,6 +158,9 @@ function describeNonFinite(result) {
     return typeof result;
 }
 const mathTsModule = {
+    // MathTS's Node type is not imported; `MathNode` is the structural subset
+    // this parser reads (`evaluate`, `filter`, `isSymbolNode`, `isFunctionNode`,
+    // `name`, `fn`), every member of which exists on a MathTS node.
     parse: (expr) => parseMathTs(expr),
 };
 /** The MathTS formula parser. @internal */