universal-physics-tensor 0.7.3 → 0.14.0

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"}

package/dist/numerical/geometrized.js

@@ -0,0 +1,66 @@
+/**
+ * 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 { UPTError } from '../dimensional/errors.js';
+import { C_SI, G_SI } from '../core/constants.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 class NonGeometrizableDimensionError extends UPTError {
+    constructor(message) {
+        super(message);
+        this.name = 'NonGeometrizableDimensionError';
+        Object.setPrototypeOf(this, NonGeometrizableDimensionError.prototype);
+    }
+}
+/**
+ * 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 function geometrizedFactor(dim) {
+    if (dim.I !== 0 || dim.Theta !== 0 || dim.N !== 0 || dim.J !== 0) {
+        throw new NonGeometrizableDimensionError(`geometrizedFactor: only c and G are in scope; this dimension has a ` +
+            `nonzero non-mechanizable exponent (I=${dim.I}, Θ=${dim.Theta}, ` +
+            `N=${dim.N}, J=${dim.J}). Geometrization handles length / mass / time.`);
+    }
+    return Math.pow(G_SI, dim.M) * Math.pow(C_SI, dim.T - 2 * dim.M);
+}
+/** SI → geometrized (G = c = 1). @public */
+export function toGeometrized(valueSI, dim) {
+    return valueSI * geometrizedFactor(dim);
+}
+/** Geometrized (G = c = 1) → SI. @public */
+export function fromGeometrized(valueGeom, dim) {
+    return valueGeom / geometrizedFactor(dim);
+}
+//# sourceMappingURL=geometrized.js.map

package/dist/numerical/geometrized.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"geometrized.js","sourceRoot":"","sources":["../../src/numerical/geometrized.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,0BAA0B,CAAC;AACpD,OAAO,EAAE,IAAI,EAAE,IAAI,EAAE,MAAM,sBAAsB,CAAC;AAElD;;;;;GAKG;AACH,MAAM,OAAO,8BAA+B,SAAQ,QAAQ;IAC1D,YAAY,OAAe;QACzB,KAAK,CAAC,OAAO,CAAC,CAAC;QACf,IAAI,CAAC,IAAI,GAAG,gCAAgC,CAAC;QAC7C,MAAM,CAAC,cAAc,CAAC,IAAI,EAAE,8BAA8B,CAAC,SAAS,CAAC,CAAC;IACxE,CAAC;CACF;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,iBAAiB,CAAC,GAAc;IAC9C,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC;QACjE,MAAM,IAAI,8BAA8B,CACtC,qEAAqE;YACnE,wCAAwC,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,KAAK,IAAI;YACjE,KAAK,GAAG,CAAC,CAAC,OAAO,GAAG,CAAC,CAAC,iDAAiD,CAC1E,CAAC;IACJ,CAAC;IACD,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;AACnE,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,aAAa,CAAC,OAAe,EAAE,GAAc;IAC3D,OAAO,OAAO,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,4CAA4C;AAC5C,MAAM,UAAU,eAAe,CAAC,SAAiB,EAAE,GAAc;IAC/D,OAAO,SAAS,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC;AAC5C,CAAC"}

package/dist/numerical/gl4-integrator.d.ts

@@ -71,14 +71,16 @@ export interface GL4Options {
     readonly steps: number;
     /** Final proper time τ_max (initial τ = 0). */
     readonly tauMax: number;
-    /** Inverse-metric closure: `gInverseFn(x)[μ][ν] = g^{μν}(x)`. */
-    readonly gInverseFn: (x: readonly number[]) => readonly (readonly number[])[];
+    /** Inverse-metric closure (v0.9.0 O-1 flat layout):
+     *  `gInverseFn(x)[μ*dim + ν] = g^{μν}(x)`, row-major Float64Array(dim²). */
+    readonly gInverseFn: (x: readonly number[]) => Float64Array;
     /**
      * Partial derivatives of the inverse metric.
-     * Index order: `dgInverseFn(x)[lambda][mu][nu] = ∂_lambda g^{mu nu}` at coords x.
+     * Index order (v0.9.0 O-1 flat layout):
+     * `dgInverseFn(x)[lambda*dim² + mu*dim + nu] = ∂_lambda g^{mu nu}` at coords x.
      * (I2: axis semantics pinned here to prevent silent transposition bugs.)
      */
-    readonly dgInverseFn: (x: readonly number[]) => readonly (readonly (readonly number[])[])[];
+    readonly dgInverseFn: (x: readonly number[]) => Float64Array;
     /** Picard fixed-point tolerance (default chosen in Task 2). */
     readonly picardTol?: number;
     /** Picard fixed-point iteration cap (default chosen in Task 2). */
@@ -141,7 +143,7 @@ interface StageSolveResult {
  * The `dgInverseFn` index order is `dg[λ][μ][ν] = ∂_λ g^{μν}` (Task 0 I2
  * pin, also recorded on `GL4Options.dgInverseFn`). When we evaluate
  * `dp_μ = −½ (∂_μ g^{νρ}) P_ν P_ρ` we therefore read
- * `dgInvAtXj[mu][nu][rho]` — `mu` is the differentiation axis (λ in the
+ * `dgInvAtXj[mu*dim²+nu*dim+rho]` — `mu` is the differentiation axis (λ in the
  * pinned order) and `(nu, rho)` are the upper metric indices.
  *
  * Throws `GL4ConvergenceError` with message matching
@@ -149,7 +151,7 @@ interface StageSolveResult {
  *
  * @internal
  */
-export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => readonly (readonly number[])[], dgInverseFn: (x: readonly number[]) => readonly (readonly (readonly number[])[])[], opts: {
+export declare function solveGL4Stage(state: GL4State, h: number, gInverseFn: (x: readonly number[]) => Float64Array, dgInverseFn: (x: readonly number[]) => Float64Array, opts: {
     picardTol: number;
     picardMaxIter: number;
 }): StageSolveResult;

package/dist/numerical/gl4-integrator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA6BA;;;;;;;;GAQG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAA4C,CAAC;AAEzF;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAGjF,CAAC;AAEF;;;;;;GAMG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAc,CAAC;AAE3D;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;GAMG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB,iEAAiE;IACjE,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC;IAC9E;;;;OAIG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,CAAC,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;IAC5F,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAED;;;;;;;GAOG;AACH,UAAU,gBAAgB;IACxB,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AACH,wBAAgB,aAAa,CAC3B,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,EACpE,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,SAAS,CAAC,SAAS,CAAC,SAAS,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,EAClF,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,gBAAgB,CAgFlB;AA8ED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAwDG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CA4DxB"}
+{"version":3,"file":"gl4-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AA6BA;;;;;;;;GAQG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAA4C,CAAC;AAEzF;;;;;;;;;;;GAWG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAGjF,CAAC;AAEF;;;;;;GAMG;AACH,eAAO,MAAM,KAAK,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAc,CAAC;AAE3D;;;;;;;;;;;GAWG;AACH,MAAM,WAAW,QAAQ;IACvB,+BAA+B;IAC/B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,yCAAyC;IACzC,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAC/B;AAED;;;;;;GAMG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,CAAC,GAAG,EAAE,MAAM,CAAC;IACrB,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;IAC9B,QAAQ,CAAC,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,CAAC;CAChC;AAED;;;;GAIG;AACH,MAAM,WAAW,UAAU;IACzB;qCACiC;IACjC,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;IACvB,+CAA+C;IAC/C,QAAQ,CAAC,MAAM,EAAE,MAAM,CAAC;IACxB;gFAC4E;IAC5E,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC5D;;;;;OAKG;IACH,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,CAAC;IAC7D,+DAA+D;IAC/D,QAAQ,CAAC,SAAS,CAAC,EAAE,MAAM,CAAC;IAC5B,mEAAmE;IACnE,QAAQ,CAAC,aAAa,CAAC,EAAE,MAAM,CAAC;IAChC,+FAA+F;IAC/F,QAAQ,CAAC,IAAI,CAAC,EAAE,MAAM,CAAC;IACvB;;qCAEiC;IACjC,QAAQ,CAAC,eAAe,CAAC,EAAE,MAAM,CAAC;IAClC;;;;;;;;;;;;;;OAcG;IACH,QAAQ,CAAC,MAAM,CAAC,EAAE,CAAC,KAAK,EAAE;QAAE,IAAI,EAAE,MAAM,CAAC;QAAC,UAAU,EAAE,MAAM,CAAC;QAAC,QAAQ,EAAE,MAAM,CAAA;KAAE,KAAK,IAAI,CAAC;CAC3F;AAED;;;;;;;GAOG;AACH,UAAU,gBAAgB;IACxB,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,MAAM,EAAE,SAAS,CAAC,SAAS,MAAM,EAAE,EAAE,SAAS,MAAM,EAAE,CAAC,CAAC;IACjE,QAAQ,CAAC,UAAU,EAAE,MAAM,CAAC;CAC7B;AAED;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AACH,wBAAgB,aAAa,CAC3B,KAAK,EAAE,QAAQ,EACf,CAAC,EAAE,MAAM,EACT,UAAU,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EAClD,WAAW,EAAE,CAAC,CAAC,EAAE,SAAS,MAAM,EAAE,KAAK,YAAY,EACnD,IAAI,EAAE;IAAE,SAAS,EAAE,MAAM,CAAC;IAAC,aAAa,EAAE,MAAM,CAAA;CAAE,GACjD,gBAAgB,CAgFlB;AA8ED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAwDG;AACH,wBAAgB,oBAAoB,CAClC,YAAY,EAAE,QAAQ,EACtB,OAAO,EAAE,UAAU,GAClB,SAAS,WAAW,EAAE,CA4DxB"}