universal-physics-tensor 0.7.3 → 0.15.0

package/dist/numerical/formula-mathts.js

@@ -0,0 +1,98 @@
+/**
+ * MathTS-backed scalar-formula parser (Path A — see
+ * docs/planning/MathTS-Formula-Integration-Design-Note.md).
+ *
+ * Implements the same {@link FormulaParser} contract as the self-contained
+ * Path B parser (`formula.ts`), backed by `@danielsimonjr/mathts-functions`'s
+ * assembled mathjs-style engine (`parse(expr) → Node`, `node.evaluate(scope)`).
+ * The optional peer is loaded dynamically (it is NOT present at tsc time —
+ * the ambient declaration in `mathts-functions.ambient.d.ts` covers it), so
+ * the package still builds and runs without it; the registry
+ * (`formula-registry.ts`) falls back to Path B when it is absent.
+ *
+ * SCALAR-ONLY guard: the CLI/inference contract returns a `number`. If a
+ * formula evaluates to a non-number (matrix, complex, unit, function), this
+ * throws a {@link FormulaError} rather than leaking MathTS types through the
+ * seam — keeping the two parsers interchangeable.
+ *
+ * @module numerical/formula-mathts
+ */
+import { FormulaError } from './formula.js';
+/**
+ * Build a {@link FormulaParser} bound to an already-loaded mathts-functions
+ * module. Pure (no I/O) — the dynamic import lives in the registry.
+ */
+function createMathtsFormulaParser(mod) {
+    const builtinCache = new Map();
+    const isBuiltin = (name) => {
+        const cached = builtinCache.get(name);
+        if (cached !== undefined)
+            return cached;
+        let builtin;
+        try {
+            mod.parse(name).evaluate({});
+            builtin = true; // pi / e / sin / … resolve with an empty scope
+        }
+        catch {
+            builtin = false; // a free variable does not
+        }
+        builtinCache.set(name, builtin);
+        return builtin;
+    };
+    return {
+        parse(expr) {
+            if (!expr || !expr.trim())
+                throw new FormulaError('empty formula');
+            let node;
+            try {
+                node = mod.parse(expr);
+            }
+            catch (err) {
+                throw new FormulaError(`parse error: ${err instanceof Error ? err.message : String(err)}`);
+            }
+            // Free variables = symbol names − function callees − built-in
+            // constants/functions (MathTS's own namespace decides "built-in").
+            const callees = new Set(node
+                .filter((n) => n.isFunctionNode === true)
+                .map((n) => n.fn?.name)
+                .filter((n) => typeof n === 'string'));
+            const variables = [
+                ...new Set(node
+                    .filter((n) => n.isSymbolNode === true)
+                    .map((n) => n.name)
+                    .filter((n) => typeof n === 'string')
+                    .filter((n) => !callees.has(n) && !isBuiltin(n))),
+            ].sort();
+            return {
+                source: expr,
+                variables,
+                evaluate(scope) {
+                    let result;
+                    try {
+                        result = node.evaluate(scope);
+                    }
+                    catch (err) {
+                        throw new FormulaError(err instanceof Error ? err.message : String(err));
+                    }
+                    if (typeof result !== 'number' || !Number.isFinite(result)) {
+                        throw new FormulaError(`formula did not evaluate to a finite number (got ${typeof result})`);
+                    }
+                    return result;
+                },
+            };
+        },
+    };
+}
+/**
+ * Dynamically load the optional peer and build the MathTS-backed parser.
+ * Throws if the peer is absent or fails to assemble — the registry catches
+ * this and falls back to the self-contained Path B parser.
+ */
+export async function loadMathtsFormulaParser() {
+    const mod = (await import('@danielsimonjr/mathts-functions'));
+    if (typeof mod.parse !== 'function') {
+        throw new FormulaError('mathts-functions: no parse() export');
+    }
+    return createMathtsFormulaParser(mod);
+}
+//# sourceMappingURL=formula-mathts.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"formula-mathts.js","sourceRoot":"","sources":["../../src/numerical/formula-mathts.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;GAkBG;AAGH,OAAO,EAAE,YAAY,EAAE,MAAM,cAAc,CAAC;AAkB5C;;;GAGG;AACH,SAAS,yBAAyB,CAChC,GAA0B;IAE1B,MAAM,YAAY,GAAG,IAAI,GAAG,EAAmB,CAAC;IAChD,MAAM,SAAS,GAAG,CAAC,IAAY,EAAW,EAAE;QAC1C,MAAM,MAAM,GAAG,YAAY,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACtC,IAAI,MAAM,KAAK,SAAS;YAAE,OAAO,MAAM,CAAC;QACxC,IAAI,OAAgB,CAAC;QACrB,IAAI,CAAC;YACH,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;YAC7B,OAAO,GAAG,IAAI,CAAC,CAAC,+CAA+C;QACjE,CAAC;QAAC,MAAM,CAAC;YACP,OAAO,GAAG,KAAK,CAAC,CAAC,2BAA2B;QAC9C,CAAC;QACD,YAAY,CAAC,GAAG,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;QAChC,OAAO,OAAO,CAAC;IACjB,CAAC,CAAC;IAEF,OAAO;QACL,KAAK,CAAC,IAAY;YAChB,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE;gBAAE,MAAM,IAAI,YAAY,CAAC,eAAe,CAAC,CAAC;YACnE,IAAI,IAAc,CAAC;YACnB,IAAI,CAAC;gBACH,IAAI,GAAG,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;YACzB,CAAC;YAAC,OAAO,GAAG,EAAE,CAAC;gBACb,MAAM,IAAI,YAAY,CACpB,gBAAgB,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,EAAE,CACnE,CAAC;YACJ,CAAC;YAED,8DAA8D;YAC9D,mEAAmE;YACnE,MAAM,OAAO,GAAG,IAAI,GAAG,CACrB,IAAI;iBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,cAAc,KAAK,IAAI,CAAC;iBACxC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,IAAI,CAAC;iBACtB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CACrD,CAAC;YACF,MAAM,SAAS,GAAG;gBAChB,GAAG,IAAI,GAAG,CACR,IAAI;qBACD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,KAAK,IAAI,CAAC;qBACtC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;qBAClB,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC;qBACjD,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CACnD;aACF,CAAC,IAAI,EAAE,CAAC;YAET,OAAO;gBACL,MAAM,EAAE,IAAI;gBACZ,SAAS;gBACT,QAAQ,CAAC,KAA6B;oBACpC,IAAI,MAAe,CAAC;oBACpB,IAAI,CAAC;wBACH,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;oBAChC,CAAC;oBAAC,OAAO,GAAG,EAAE,CAAC;wBACb,MAAM,IAAI,YAAY,CACpB,GAAG,YAAY,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CACjD,CAAC;oBACJ,CAAC;oBACD,IAAI,OAAO,MAAM,KAAK,QAAQ,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC;wBAC3D,MAAM,IAAI,YAAY,CACpB,oDAAoD,OAAO,MAAM,GAAG,CACrE,CAAC;oBACJ,CAAC;oBACD,OAAO,MAAM,CAAC;gBAChB,CAAC;aACF,CAAC;QACJ,CAAC;KACF,CAAC;AACJ,CAAC;AAED;;;;GAIG;AACH,MAAM,CAAC,KAAK,UAAU,uBAAuB;IAC3C,MAAM,GAAG,GAAG,CAAC,MAAM,MAAM,CACvB,iCAAiC,CAClC,CAAqC,CAAC;IACvC,IAAI,OAAO,GAAG,CAAC,KAAK,KAAK,UAAU,EAAE,CAAC;QACpC,MAAM,IAAI,YAAY,CAAC,qCAAqC,CAAC,CAAC;IAChE,CAAC;IACD,OAAO,yBAAyB,CAAC,GAAG,CAAC,CAAC;AACxC,CAAC"}

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

@@ -0,0 +1,30 @@
+/**
+ * Formula-parser registry (Path A selector — mirrors `engine-registry.ts`).
+ *
+ * `getFormulaParser()` returns the MathTS-backed parser when the optional
+ * peer is installed AND actually assembles and evaluates (smoke-tested),
+ * otherwise the self-contained Path B parser. The choice is detected once
+ * and cached. A broken or absent MathTS never breaks the caller — it falls
+ * back silently to Path B.
+ *
+ * The CLI and any inference code depend only on the {@link FormulaParser}
+ * interface, so swapping Path B ↔ Path A is transparent.
+ *
+ * @module numerical/formula-registry
+ */
+import type { FormulaParser } from './formula.js';
+import type { FormulaDimensionChecker } from './formula-dimension.js';
+type FormulaParserKind = 'mathts' | 'builtin';
+/** Resolve the active formula parser (cached). @internal */
+export declare function getFormulaParser(): Promise<FormulaParser>;
+/** Which parser is active — `mathts` (Path A) or `builtin` (Path B). @internal */
+export declare function getFormulaParserKind(): Promise<FormulaParserKind>;
+/**
+ * The dimensional checker for user formulas (MathTS Phase 2). Always
+ * available: it uses the MathTS AST when the peer is installed and
+ * smoke-tests clean, else the self-contained Path B AST — both transpile to
+ * the same `ExprNode`, so the dimensional verdict is identical. @internal
+ */
+export declare function getFormulaDimensionChecker(): Promise<FormulaDimensionChecker>;
+export {};
+//# sourceMappingURL=formula-registry.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"formula-registry.d.ts","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAEH,OAAO,KAAK,EAAE,aAAa,EAAE,MAAM,cAAc,CAAC;AAGlD,OAAO,KAAK,EAAE,uBAAuB,EAAE,MAAM,wBAAwB,CAAC;AAGtE,KAAK,iBAAiB,GAAG,QAAQ,GAAG,SAAS,CAAC;AA8C9C,4DAA4D;AAC5D,wBAAsB,gBAAgB,IAAI,OAAO,CAAC,aAAa,CAAC,CAG/D;AAED,kFAAkF;AAClF,wBAAsB,oBAAoB,IAAI,OAAO,CAAC,iBAAiB,CAAC,CAGvE;AAMD;;;;;GAKG;AACH,wBAAsB,0BAA0B,IAAI,OAAO,CAAC,uBAAuB,CAAC,CAYnF"}

package/dist/numerical/formula-registry.js

@@ -0,0 +1,88 @@
+/**
+ * Formula-parser registry (Path A selector — mirrors `engine-registry.ts`).
+ *
+ * `getFormulaParser()` returns the MathTS-backed parser when the optional
+ * peer is installed AND actually assembles and evaluates (smoke-tested),
+ * otherwise the self-contained Path B parser. The choice is detected once
+ * and cached. A broken or absent MathTS never breaks the caller — it falls
+ * back silently to Path B.
+ *
+ * The CLI and any inference code depend only on the {@link FormulaParser}
+ * interface, so swapping Path B ↔ Path A is transparent.
+ *
+ * @module numerical/formula-registry
+ */
+import { defaultFormulaParser } from './formula.js';
+import { loadMathtsFormulaParser } from './formula-mathts.js';
+import { loadFormulaDimensionChecker, builtinFormulaDimensionChecker } from './formula-dimension.js';
+let cached;
+/** A formula MathTS must reproduce to be accepted (the §1-gate smoke test). */
+const SMOKE_EXPR = 'a*b^2 + 1';
+const SMOKE_SCOPE = { a: 3, b: 4 };
+const SMOKE_EXPECTED = 49; // 3*16 + 1
+/** Run `fn` with MathTS's import-time WASM-fallback chatter suppressed. */
+async function quietly(fn) {
+    const origWarn = console.warn;
+    const origError = console.error;
+    const origWrite = process.stderr.write.bind(process.stderr);
+    console.warn = () => { };
+    console.error = () => { };
+    // mathts prints some notices via process.stderr.write directly.
+    process.stderr.write = () => true;
+    try {
+        return await fn();
+    }
+    finally {
+        console.warn = origWarn;
+        console.error = origError;
+        process.stderr.write = origWrite;
+    }
+}
+async function detect() {
+    try {
+        const parser = await quietly(loadMathtsFormulaParser);
+        // Smoke test: it must actually assemble AND evaluate correctly.
+        const v = parser.parse(SMOKE_EXPR).evaluate(SMOKE_SCOPE);
+        if (Math.abs(v - SMOKE_EXPECTED) < 1e-9) {
+            return { parser, kind: 'mathts' };
+        }
+    }
+    catch {
+        /* peer absent or failed to assemble — fall through */
+    }
+    return { parser: defaultFormulaParser, kind: 'builtin' };
+}
+/** Resolve the active formula parser (cached). @internal */
+export async function getFormulaParser() {
+    cached ??= detect();
+    return (await cached).parser;
+}
+/** Which parser is active — `mathts` (Path A) or `builtin` (Path B). @internal */
+export async function getFormulaParserKind() {
+    cached ??= detect();
+    return (await cached).kind;
+}
+let cachedChecker;
+const LENGTH_DIM = { L: 1, M: 0, T: 0, I: 0, Theta: 0, N: 0, J: 0 };
+/**
+ * The dimensional checker for user formulas (MathTS Phase 2). Always
+ * available: it uses the MathTS AST when the peer is installed and
+ * smoke-tests clean, else the self-contained Path B AST — both transpile to
+ * the same `ExprNode`, so the dimensional verdict is identical. @internal
+ */
+export async function getFormulaDimensionChecker() {
+    cachedChecker ??= (async () => {
+        try {
+            const checker = await quietly(loadFormulaDimensionChecker);
+            // smoke test: a dimensionless ratio is homogeneous and dimensionless.
+            if (checker.check('a/a', { a: LENGTH_DIM }).ok)
+                return checker;
+        }
+        catch {
+            /* peer absent or failed — fall back to the built-in checker */
+        }
+        return builtinFormulaDimensionChecker();
+    })();
+    return cachedChecker;
+}
+//# sourceMappingURL=formula-registry.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"formula-registry.js","sourceRoot":"","sources":["../../src/numerical/formula-registry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AAGH,OAAO,EAAE,oBAAoB,EAAE,MAAM,cAAc,CAAC;AACpD,OAAO,EAAE,uBAAuB,EAAE,MAAM,qBAAqB,CAAC;AAE9D,OAAO,EAAE,2BAA2B,EAAE,8BAA8B,EAAE,MAAM,wBAAwB,CAAC;AASrG,IAAI,MAAqC,CAAC;AAE1C,+EAA+E;AAC/E,MAAM,UAAU,GAAG,WAAW,CAAC;AAC/B,MAAM,WAAW,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AACnC,MAAM,cAAc,GAAG,EAAE,CAAC,CAAC,WAAW;AAEtC,2EAA2E;AAC3E,KAAK,UAAU,OAAO,CAAI,EAAoB;IAC5C,MAAM,QAAQ,GAAG,OAAO,CAAC,IAAI,CAAC;IAC9B,MAAM,SAAS,GAAG,OAAO,CAAC,KAAK,CAAC;IAChC,MAAM,SAAS,GAAG,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC;IAC5D,OAAO,CAAC,IAAI,GAAG,GAAG,EAAE,GAAE,CAAC,CAAC;IACxB,OAAO,CAAC,KAAK,GAAG,GAAG,EAAE,GAAE,CAAC,CAAC;IACzB,gEAAgE;IAC/D,OAAO,CAAC,MAA6B,CAAC,KAAK,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC;IAC1D,IAAI,CAAC;QACH,OAAO,MAAM,EAAE,EAAE,CAAC;IACpB,CAAC;YAAS,CAAC;QACT,OAAO,CAAC,IAAI,GAAG,QAAQ,CAAC;QACxB,OAAO,CAAC,KAAK,GAAG,SAAS,CAAC;QACzB,OAAO,CAAC,MAA6B,CAAC,KAAK,GAAG,SAAS,CAAC;IAC3D,CAAC;AACH,CAAC;AAED,KAAK,UAAU,MAAM;IACnB,IAAI,CAAC;QACH,MAAM,MAAM,GAAG,MAAM,OAAO,CAAC,uBAAuB,CAAC,CAAC;QACtD,gEAAgE;QAChE,MAAM,CAAC,GAAG,MAAM,CAAC,KAAK,CAAC,UAAU,CAAC,CAAC,QAAQ,CAAC,WAAW,CAAC,CAAC;QACzD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,cAAc,CAAC,GAAG,IAAI,EAAE,CAAC;YACxC,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,CAAC;QACpC,CAAC;IACH,CAAC;IAAC,MAAM,CAAC;QACP,sDAAsD;IACxD,CAAC;IACD,OAAO,EAAE,MAAM,EAAE,oBAAoB,EAAE,IAAI,EAAE,SAAS,EAAE,CAAC;AAC3D,CAAC;AAED,4DAA4D;AAC5D,MAAM,CAAC,KAAK,UAAU,gBAAgB;IACpC,MAAM,KAAK,MAAM,EAAE,CAAC;IACpB,OAAO,CAAC,MAAM,MAAM,CAAC,CAAC,MAAM,CAAC;AAC/B,CAAC;AAED,kFAAkF;AAClF,MAAM,CAAC,KAAK,UAAU,oBAAoB;IACxC,MAAM,KAAK,MAAM,EAAE,CAAC;IACpB,OAAO,CAAC,MAAM,MAAM,CAAC,CAAC,IAAI,CAAC;AAC7B,CAAC;AAED,IAAI,aAA2D,CAAC;AAEhE,MAAM,UAAU,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;AAEpE;;;;;GAKG;AACH,MAAM,CAAC,KAAK,UAAU,0BAA0B;IAC9C,aAAa,KAAK,CAAC,KAAK,IAAI,EAAE;QAC5B,IAAI,CAAC;YACH,MAAM,OAAO,GAAG,MAAM,OAAO,CAAC,2BAA2B,CAAC,CAAC;YAC3D,sEAAsE;YACtE,IAAI,OAAO,CAAC,KAAK,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,CAAC,EAAE;gBAAE,OAAO,OAAO,CAAC;QACjE,CAAC;QAAC,MAAM,CAAC;YACP,+DAA+D;QACjE,CAAC;QACD,OAAO,8BAA8B,EAAE,CAAC;IAC1C,CAAC,CAAC,EAAE,CAAC;IACL,OAAO,aAAa,CAAC;AACvB,CAAC"}

package/dist/numerical/formula.d.ts

@@ -0,0 +1,75 @@
+/**
+ * Self-contained scalar-formula parser/evaluator (Path B).
+ *
+ * A small, dependency-free recursive-descent parser for closed-form
+ * scalar physics expressions — `hbar*c^3/(8*pi*G*M*k_B)` and the like —
+ * so non-TypeScript users can supply their OWN equation to the CLI
+ * without UPT shipping a code evaluator. It is SAFE: there is no `eval`,
+ * no function constructor, no property access — only arithmetic over
+ * numbers, a fixed whitelist of math functions, and the variables the
+ * caller explicitly supplies. An unknown symbol is an error, never an
+ * implicit global.
+ *
+ * It sits behind the {@link FormulaParser} interface so a MathTS-backed
+ * parser (Path A) can be dropped in later by implementing the same
+ * contract — the inference/CLI code depends only on the interface.
+ *
+ * Grammar (standard precedence; `^` is right-associative; unary minus
+ * binds looser than `^`, so `-2^2 = -4` and `2^-2` parses):
+ *   expr  → term (('+'|'-') term)*
+ *   term  → unary (('*'|'/') unary)*
+ *   unary → ('+'|'-') unary | power
+ *   power → atom ('^' unary)?
+ *   atom  → NUMBER | NAME | NAME '(' args ')' | '(' expr ')'
+ *
+ * @module numerical/formula
+ */
+/** A parse or evaluation failure (bad syntax, unknown symbol, arity). */
+export declare class FormulaError extends Error {
+    constructor(message: string);
+}
+/** A parsed formula: its free variables and a safe evaluator. @internal */
+export interface CompiledFormula {
+    /** The original source string. */
+    readonly source: string;
+    /** Free variable names (excludes built-in constants and functions). */
+    readonly variables: readonly string[];
+    /** Evaluate against a `name → value` scope. Throws on a missing var. */
+    evaluate(scope: Record<string, number>): number;
+}
+/** The swap point for Path A (a MathTS-backed parser implements this). */
+export interface FormulaParser {
+    parse(expr: string): CompiledFormula;
+}
+type Node = {
+    kind: 'num';
+    value: number;
+} | {
+    kind: 'sym';
+    name: string;
+} | {
+    kind: 'unary';
+    op: '+' | '-';
+    arg: Node;
+} | {
+    kind: 'bin';
+    op: '+' | '-' | '*' | '/' | '^';
+    left: Node;
+    right: Node;
+} | {
+    kind: 'call';
+    fn: string;
+    args: Node[];
+};
+/** The self-contained Path B parser. @internal */
+export declare const defaultFormulaParser: FormulaParser;
+/** Parse a scalar formula with the default (self-contained) parser. @internal */
+export declare function parseFormula(expr: string): CompiledFormula;
+/** The Path B formula AST node. @internal */
+export type FormulaAstNode = Node;
+/** Parse a formula to its Path B AST. @internal */
+export declare const parseFormulaToAst: (expr: string) => FormulaAstNode;
+/** Evaluate a Path B AST node against a scope. @internal */
+export declare const evalFormulaAst: (node: FormulaAstNode, scope: Record<string, number>) => number;
+export {};
+//# sourceMappingURL=formula.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"formula.d.ts","sourceRoot":"","sources":["../../src/numerical/formula.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAEH,yEAAyE;AACzE,qBAAa,YAAa,SAAQ,KAAK;gBACzB,OAAO,EAAE,MAAM;CAI5B;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,0EAA0E;AAC1E,MAAM,WAAW,aAAa;IAC5B,KAAK,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAAC;CACtC;AA4CD,KAAK,IAAI,GACL;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,OAAO,CAAC;IAAC,EAAE,EAAE,GAAG,GAAG,GAAG,CAAC;IAAC,GAAG,EAAE,IAAI,CAAA;CAAE,GAC3C;IAAE,IAAI,EAAE,KAAK,CAAC;IAAC,EAAE,EAAE,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;IAAC,IAAI,EAAE,IAAI,CAAC;IAAC,KAAK,EAAE,IAAI,CAAA;CAAE,GACzE;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,EAAE,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,IAAI,EAAE,CAAA;CAAE,CAAC;AA4L/C,kDAAkD;AAClD,eAAO,MAAM,oBAAoB,EAAE,aAalC,CAAC;AAEF,iFAAiF;AACjF,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,GAAG,eAAe,CAE1D;AAQD,6CAA6C;AAC7C,MAAM,MAAM,cAAc,GAAG,IAAI,CAAC;AAClC,mDAAmD;AACnD,eAAO,MAAM,iBAAiB,EAAE,CAAC,IAAI,EAAE,MAAM,KAAK,cAA2B,CAAC;AAC9E,4DAA4D;AAC5D,eAAO,MAAM,cAAc,EAAE,CAAC,IAAI,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,KAAK,MAAiB,CAAC"}

package/dist/numerical/formula.js

@@ -0,0 +1,275 @@
+/**
+ * Self-contained scalar-formula parser/evaluator (Path B).
+ *
+ * A small, dependency-free recursive-descent parser for closed-form
+ * scalar physics expressions — `hbar*c^3/(8*pi*G*M*k_B)` and the like —
+ * so non-TypeScript users can supply their OWN equation to the CLI
+ * without UPT shipping a code evaluator. It is SAFE: there is no `eval`,
+ * no function constructor, no property access — only arithmetic over
+ * numbers, a fixed whitelist of math functions, and the variables the
+ * caller explicitly supplies. An unknown symbol is an error, never an
+ * implicit global.
+ *
+ * It sits behind the {@link FormulaParser} interface so a MathTS-backed
+ * parser (Path A) can be dropped in later by implementing the same
+ * contract — the inference/CLI code depends only on the interface.
+ *
+ * Grammar (standard precedence; `^` is right-associative; unary minus
+ * binds looser than `^`, so `-2^2 = -4` and `2^-2` parses):
+ *   expr  → term (('+'|'-') term)*
+ *   term  → unary (('*'|'/') unary)*
+ *   unary → ('+'|'-') unary | power
+ *   power → atom ('^' unary)?
+ *   atom  → NUMBER | NAME | NAME '(' args ')' | '(' expr ')'
+ *
+ * @module numerical/formula
+ */
+/** A parse or evaluation failure (bad syntax, unknown symbol, arity). */
+export class FormulaError extends Error {
+    constructor(message) {
+        super(message);
+        this.name = 'FormulaError';
+    }
+}
+// --- built-ins ------------------------------------------------------------
+const CONSTANTS = {
+    pi: Math.PI,
+    tau: 2 * Math.PI,
+};
+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 NUM_RE = /^(?:\d+\.?\d*|\.\d+)(?:[eE][+-]?\d+)?/;
+const NAME_RE = /^[A-Za-z_][A-Za-z0-9_]*/;
+function tokenize(src) {
+    const toks = [];
+    let s = src;
+    while (s.length) {
+        const ch = s[0];
+        if (ch === ' ' || ch === '\t' || ch === '\n') {
+            s = s.slice(1);
+            continue;
+        }
+        if ('+-*/^(),'.includes(ch)) {
+            toks.push({ t: 'op', v: ch });
+            s = s.slice(1);
+            continue;
+        }
+        const num = NUM_RE.exec(s);
+        if (num) {
+            toks.push({ t: 'num', v: Number(num[0]) });
+            s = s.slice(num[0].length);
+            continue;
+        }
+        const name = NAME_RE.exec(s);
+        if (name) {
+            toks.push({ t: 'name', v: name[0] });
+            s = s.slice(name[0].length);
+            continue;
+        }
+        throw new FormulaError(`unexpected character '${ch}' in formula`);
+    }
+    return toks;
+}
+// --- parser (recursive descent) ------------------------------------------
+function parseToAst(src) {
+    const toks = tokenize(src);
+    let i = 0;
+    const peek = () => toks[i];
+    const eat = (v) => {
+        const tk = toks[i];
+        if (!tk)
+            throw new FormulaError('unexpected end of formula');
+        if (v !== undefined && !(tk.t === 'op' && tk.v === v)) {
+            throw new FormulaError(`expected '${v}'`);
+        }
+        i++;
+        return tk;
+    };
+    const isOp = (v) => {
+        const tk = peek();
+        return !!tk && tk.t === 'op' && tk.v === v;
+    };
+    function expr() {
+        let node = term();
+        while (isOp('+') || isOp('-')) {
+            const op = eat().v;
+            node = { kind: 'bin', op, left: node, right: term() };
+        }
+        return node;
+    }
+    function term() {
+        let node = unary();
+        while (isOp('*') || isOp('/')) {
+            const op = eat().v;
+            node = { kind: 'bin', op, left: node, right: unary() };
+        }
+        return node;
+    }
+    function unary() {
+        if (isOp('+') || isOp('-')) {
+            const op = eat().v;
+            return { kind: 'unary', op, arg: unary() };
+        }
+        return power();
+    }
+    function power() {
+        const base = atom();
+        if (isOp('^')) {
+            eat('^');
+            return { kind: 'bin', op: '^', left: base, right: unary() };
+        }
+        return base;
+    }
+    function atom() {
+        const tk = peek();
+        if (!tk)
+            throw new FormulaError('unexpected end of formula');
+        if (tk.t === 'num') {
+            eat();
+            return { kind: 'num', value: tk.v };
+        }
+        if (tk.t === 'name') {
+            eat();
+            if (isOp('(')) {
+                eat('(');
+                const args = [];
+                if (!isOp(')')) {
+                    args.push(expr());
+                    while (isOp(',')) {
+                        eat(',');
+                        args.push(expr());
+                    }
+                }
+                eat(')');
+                return { kind: 'call', fn: tk.v, args };
+            }
+            return { kind: 'sym', name: tk.v };
+        }
+        if (tk.t === 'op' && tk.v === '(') {
+            eat('(');
+            const node = expr();
+            eat(')');
+            return node;
+        }
+        throw new FormulaError(`unexpected token '${tk.v}'`);
+    }
+    const node = expr();
+    if (i !== toks.length) {
+        throw new FormulaError(`unexpected trailing input in formula`);
+    }
+    return node;
+}
+function collectSymbols(node, out) {
+    switch (node.kind) {
+        case 'num':
+            return;
+        case 'sym':
+            if (!(node.name in CONSTANTS))
+                out.add(node.name);
+            return;
+        case 'unary':
+            collectSymbols(node.arg, out);
+            return;
+        case 'bin':
+            collectSymbols(node.left, out);
+            collectSymbols(node.right, out);
+            return;
+        case 'call':
+            for (const a of node.args)
+                collectSymbols(a, out);
+            return;
+    }
+}
+function evalNode(node, scope) {
+    switch (node.kind) {
+        case 'num':
+            return node.value;
+        case 'sym': {
+            if (node.name in CONSTANTS)
+                return CONSTANTS[node.name];
+            const v = scope[node.name];
+            if (v === undefined)
+                throw new FormulaError(`unknown variable '${node.name}'`);
+            return v;
+        }
+        case 'unary': {
+            const a = evalNode(node.arg, scope);
+            return node.op === '-' ? -a : a;
+        }
+        case 'bin': {
+            const l = evalNode(node.left, scope);
+            const r = evalNode(node.right, scope);
+            switch (node.op) {
+                case '+': return l + r;
+                case '-': return l - r;
+                case '*': return l * r;
+                case '/': return l / r;
+                case '^': return Math.pow(l, r);
+            }
+        }
+        // eslint-disable-next-line no-fallthrough
+        case 'call': {
+            const fn = FUNCTIONS[node.fn];
+            if (!fn)
+                throw new FormulaError(`unknown function '${node.fn}'`);
+            return fn(node.args.map((a) => evalNode(a, scope)));
+        }
+    }
+}
+/** The self-contained Path B parser. @internal */
+export const defaultFormulaParser = {
+    parse(expr) {
+        if (!expr || !expr.trim())
+            throw new FormulaError('empty formula');
+        const ast = parseToAst(expr);
+        const syms = new Set();
+        collectSymbols(ast, syms);
+        const variables = [...syms].sort();
+        return {
+            source: expr,
+            variables,
+            evaluate: (scope) => evalNode(ast, scope),
+        };
+    },
+};
+/** Parse a scalar formula with the default (self-contained) parser. @internal */
+export function parseFormula(expr) {
+    return defaultFormulaParser.parse(expr);
+}
+/** Parse a formula to its Path B AST. @internal */
+export const parseFormulaToAst = parseToAst;
+/** Evaluate a Path B AST node against a scope. @internal */
+export const evalFormulaAst = evalNode;
+//# sourceMappingURL=formula.js.map

package/dist/numerical/formula.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"formula.js","sourceRoot":"","sources":["../../src/numerical/formula.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;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,6EAA6E;AAE7E,MAAM,SAAS,GAAqC;IAClD,EAAE,EAAE,IAAI,CAAC,EAAE;IACX,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE;CACjB,CAAC;AAGF,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,CAAC,CAAC,CAAC;AACjB,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,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9B,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,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAChC,CAAC;CACF,CAAC;AAkBF,MAAM,MAAM,GAAG,uCAAuC,CAAC;AACvD,MAAM,OAAO,GAAG,yBAAyB,CAAC;AAE1C,SAAS,QAAQ,CAAC,GAAW;IAC3B,MAAM,IAAI,GAAU,EAAE,CAAC;IACvB,IAAI,CAAC,GAAG,GAAG,CAAC;IACZ,OAAO,CAAC,CAAC,MAAM,EAAE,CAAC;QAChB,MAAM,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAChB,IAAI,EAAE,KAAK,GAAG,IAAI,EAAE,KAAK,IAAI,IAAI,EAAE,KAAK,IAAI,EAAE,CAAC;YAC7C,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACf,SAAS;QACX,CAAC;QACD,IAAI,UAAU,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;YAC5B,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;YAC9B,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACf,SAAS;QACX,CAAC;QACD,MAAM,GAAG,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,GAAG,EAAE,CAAC;YACR,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC3C,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YAC3B,SAAS;QACX,CAAC;QACD,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;QAC7B,IAAI,IAAI,EAAE,CAAC;YACT,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACrC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;YAC5B,SAAS;QACX,CAAC;QACD,MAAM,IAAI,YAAY,CAAC,yBAAyB,EAAE,cAAc,CAAC,CAAC;IACpE,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,4EAA4E;AAE5E,SAAS,UAAU,CAAC,GAAW;IAC7B,MAAM,IAAI,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;IAC3B,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,MAAM,IAAI,GAAG,GAAoB,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,GAAG,GAAG,CAAC,CAAU,EAAO,EAAE;QAC9B,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACnB,IAAI,CAAC,EAAE;YAAE,MAAM,IAAI,YAAY,CAAC,2BAA2B,CAAC,CAAC;QAC7D,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,KAAK,IAAI,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;YACtD,MAAM,IAAI,YAAY,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;QAC5C,CAAC;QACD,CAAC,EAAE,CAAC;QACJ,OAAO,EAAE,CAAC;IACZ,CAAC,CAAC;IACF,MAAM,IAAI,GAAG,CAAC,CAAS,EAAW,EAAE;QAClC,MAAM,EAAE,GAAG,IAAI,EAAE,CAAC;QAClB,OAAO,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,CAAC,KAAK,IAAI,IAAI,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;IAC7C,CAAC,CAAC;IAEF,SAAS,IAAI;QACX,IAAI,IAAI,GAAG,IAAI,EAAE,CAAC;QAClB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAI,GAAG,EAAuB,CAAC,CAAC,CAAC;YACzC,IAAI,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC;QACxD,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,SAAS,IAAI;QACX,IAAI,IAAI,GAAG,KAAK,EAAE,CAAC;QACnB,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAC9B,MAAM,EAAE,GAAI,GAAG,EAAuB,CAAC,CAAC,CAAC;YACzC,IAAI,GAAG,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QACzD,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,SAAS,KAAK;QACZ,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YAC3B,MAAM,EAAE,GAAI,GAAG,EAAuB,CAAC,CAAC,CAAC;YACzC,OAAO,EAAE,IAAI,EAAE,OAAO,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC;QAC7C,CAAC;QACD,OAAO,KAAK,EAAE,CAAC;IACjB,CAAC;IACD,SAAS,KAAK;QACZ,MAAM,IAAI,GAAG,IAAI,EAAE,CAAC;QACpB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACd,GAAG,CAAC,GAAG,CAAC,CAAC;YACT,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,EAAE,EAAE,GAAG,EAAE,IAAI,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC;QAC9D,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC;IACD,SAAS,IAAI;QACX,MAAM,EAAE,GAAG,IAAI,EAAE,CAAC;QAClB,IAAI,CAAC,EAAE;YAAE,MAAM,IAAI,YAAY,CAAC,2BAA2B,CAAC,CAAC;QAC7D,IAAI,EAAE,CAAC,CAAC,KAAK,KAAK,EAAE,CAAC;YACnB,GAAG,EAAE,CAAC;YACN,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC;QACtC,CAAC;QACD,IAAI,EAAE,CAAC,CAAC,KAAK,MAAM,EAAE,CAAC;YACpB,GAAG,EAAE,CAAC;YACN,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;gBACd,GAAG,CAAC,GAAG,CAAC,CAAC;gBACT,MAAM,IAAI,GAAW,EAAE,CAAC;gBACxB,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;oBACf,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;oBAClB,OAAO,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;wBACjB,GAAG,CAAC,GAAG,CAAC,CAAC;wBACT,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE,CAAC,CAAC;oBACpB,CAAC;gBACH,CAAC;gBACD,GAAG,CAAC,GAAG,CAAC,CAAC;gBACT,OAAO,EAAE,IAAI,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC;YAC1C,CAAC;YACD,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC;QACrC,CAAC;QACD,IAAI,EAAE,CAAC,CAAC,KAAK,IAAI,IAAI,EAAE,CAAC,CAAC,KAAK,GAAG,EAAE,CAAC;YAClC,GAAG,CAAC,GAAG,CAAC,CAAC;YACT,MAAM,IAAI,GAAG,IAAI,EAAE,CAAC;YACpB,GAAG,CAAC,GAAG,CAAC,CAAC;YACT,OAAO,IAAI,CAAC;QACd,CAAC;QACD,MAAM,IAAI,YAAY,CAAC,qBAAsB,EAAqB,CAAC,CAAC,GAAG,CAAC,CAAC;IAC3E,CAAC;IAED,MAAM,IAAI,GAAG,IAAI,EAAE,CAAC;IACpB,IAAI,CAAC,KAAK,IAAI,CAAC,MAAM,EAAE,CAAC;QACtB,MAAM,IAAI,YAAY,CAAC,sCAAsC,CAAC,CAAC;IACjE,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,cAAc,CAAC,IAAU,EAAE,GAAgB;IAClD,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO;QACT,KAAK,KAAK;YACR,IAAI,CAAC,CAAC,IAAI,CAAC,IAAI,IAAI,SAAS,CAAC;gBAAE,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAClD,OAAO;QACT,KAAK,OAAO;YACV,cAAc,CAAC,IAAI,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YAC9B,OAAO;QACT,KAAK,KAAK;YACR,cAAc,CAAC,IAAI,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC;YAC/B,cAAc,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;YAChC,OAAO;QACT,KAAK,MAAM;YACT,KAAK,MAAM,CAAC,IAAI,IAAI,CAAC,IAAI;gBAAE,cAAc,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;YAClD,OAAO;IACX,CAAC;AACH,CAAC;AAED,SAAS,QAAQ,CAAC,IAAU,EAAE,KAA6B;IACzD,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,KAAK;YACR,OAAO,IAAI,CAAC,KAAK,CAAC;QACpB,KAAK,KAAK,CAAC,CAAC,CAAC;YACX,IAAI,IAAI,CAAC,IAAI,IAAI,SAAS;gBAAE,OAAO,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YACxD,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC3B,IAAI,CAAC,KAAK,SAAS;gBAAE,MAAM,IAAI,YAAY,CAAC,qBAAqB,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC;YAC/E,OAAO,CAAC,CAAC;QACX,CAAC;QACD,KAAK,OAAO,CAAC,CAAC,CAAC;YACb,MAAM,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;YACpC,OAAO,IAAI,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,CAAC;QACD,KAAK,KAAK,CAAC,CAAC,CAAC;YACX,MAAM,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;YACrC,MAAM,CAAC,GAAG,QAAQ,CAAC,IAAI,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;YACtC,QAAQ,IAAI,CAAC,EAAE,EAAE,CAAC;gBAChB,KAAK,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;gBACvB,KAAK,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;gBACvB,KAAK,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;gBACvB,KAAK,GAAG,CAAC,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;gBACvB,KAAK,GAAG,CAAC,CAAC,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAClC,CAAC;QACH,CAAC;QACD,0CAA0C;QAC1C,KAAK,MAAM,CAAC,CAAC,CAAC;YACZ,MAAM,EAAE,GAAG,SAAS,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;YAC9B,IAAI,CAAC,EAAE;gBAAE,MAAM,IAAI,YAAY,CAAC,qBAAqB,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;YACjE,OAAO,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,QAAQ,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;QACtD,CAAC;IACH,CAAC;AACH,CAAC;AAED,kDAAkD;AAClD,MAAM,CAAC,MAAM,oBAAoB,GAAkB;IACjD,KAAK,CAAC,IAAY;QAChB,IAAI,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,EAAE;YAAE,MAAM,IAAI,YAAY,CAAC,eAAe,CAAC,CAAC;QACnE,MAAM,GAAG,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAG,IAAI,GAAG,EAAU,CAAC;QAC/B,cAAc,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QAC1B,MAAM,SAAS,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,IAAI,EAAE,CAAC;QACnC,OAAO;YACL,MAAM,EAAE,IAAI;YACZ,SAAS;YACT,QAAQ,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,QAAQ,CAAC,GAAG,EAAE,KAAK,CAAC;SAC1C,CAAC;IACJ,CAAC;CACF,CAAC;AAEF,iFAAiF;AACjF,MAAM,UAAU,YAAY,CAAC,IAAY;IACvC,OAAO,oBAAoB,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;AAC1C,CAAC;AAUD,mDAAmD;AACnD,MAAM,CAAC,MAAM,iBAAiB,GAAqC,UAAU,CAAC;AAC9E,4DAA4D;AAC5D,MAAM,CAAC,MAAM,cAAc,GAAoE,QAAQ,CAAC"}

package/dist/numerical/geometrized.d.ts

@@ -0,0 +1,51 @@
+/**
+ * Geometrized-units boundary adapters (G-9 increment 1).
+ *
+ * Converts a scalar between SI and geometrized (G = c = 1) units, driven
+ * MECHANICALLY by the `Dimension` exponent vector — the dimension functor
+ * earning its keep. For an SI quantity of dimension L^l M^m T^t the single
+ * conversion factor is `G^m · c^(t − 2m)` (each kg → G/c² metres, each second
+ * → c metres); the geometrized quantity is a pure length L^(l+m+t).
+ *
+ * This is the self-contained FOUNDATION of the units-normalization layer
+ * (design: docs/planning/v0.10.0-Units-Normalization-Design-Note.md, Adam-
+ * vetted r2; plan: docs/planning/v0.13-G9-Adapters-Plan.md, Eve-vetted). Only c
+ * and G are in scope: a nonzero electromagnetic / thermal / molar / luminous
+ * exponent (I/Θ/N/J) throws.
+ *
+ * PUBLIC boundary API (G-9 increment 2): a consumer running the GR pipeline in
+ * geometrized (G = c = 1) units converts scalars at the boundary with these
+ * adapters. The consumer-wide DEFAULT-pipeline migration (increment 3) was
+ * DECLINED after adversarial review — it is measured to give no precision win
+ * (and far worse far-field; see docs/planning/v0.14-G9-Increment3-Disposition.md).
+ * The SI pipeline stays the default; this geometrized layer is the optional
+ * boundary API.
+ *
+ * @module numerical/geometrized
+ */
+import type { Dimension } from '../dimensional/types.js';
+import { UPTError } from '../dimensional/errors.js';
+/**
+ * A dimension carries a non-c/G base (I/Θ/N/J ≠ 0) and cannot be geometrized
+ * by c and G alone (k_B / ε₀ extensions are out of scope for this increment).
+ *
+ * @public
+ */
+export declare class NonGeometrizableDimensionError extends UPTError {
+    constructor(message: string);
+}
+/**
+ * The single SI→geometrized conversion factor for `dim`: `G^M · c^(T − 2M)`
+ * (only the mass and time exponents enter). Shared by `toGeometrized` /
+ * `fromGeometrized` so the round-trip is a clean inverse. A nonzero I/Θ/N/J
+ * exponent throws `NonGeometrizableDimensionError` BEFORE any factor is
+ * computed.
+ *
+ * @public
+ */
+export declare function geometrizedFactor(dim: Dimension): number;
+/** SI → geometrized (G = c = 1). @public */
+export declare function toGeometrized(valueSI: number, dim: Dimension): number;
+/** Geometrized (G = c = 1) → SI. @public */
+export declare function fromGeometrized(valueGeom: number, dim: Dimension): number;
+//# sourceMappingURL=geometrized.d.ts.map

package/dist/numerical/geometrized.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"geometrized.d.ts","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AAGpD;;;;;GAKG;AACH,qBAAa,8BAA+B,SAAQ,QAAQ;gBAC9C,OAAO,EAAE,MAAM;CAK5B;AAED;;;;;;;;GAQG;AACH,wBAAgB,iBAAiB,CAAC,GAAG,EAAE,SAAS,GAAG,MAAM,CASxD;AAED,4CAA4C;AAC5C,wBAAgB,aAAa,CAAC,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAErE;AAED,4CAA4C;AAC5C,wBAAgB,eAAe,CAAC,SAAS,EAAE,MAAM,EAAE,GAAG,EAAE,SAAS,GAAG,MAAM,CAEzE"}