universal-physics-tensor 0.4.5

package/dist/numerical/mathts-engine.js

@@ -0,0 +1,164 @@
+/**
+ * MathTSEngine — a TensorEngine implementation backed by
+ * @danielsimonjr/mathts-tensor's rank-N Tensor. The second v0.3.5 engine;
+ * becomes UPT's default in v0.4.0 (v0.3.5-Design.md §12).
+ *
+ * Thin adapter: it translates the TensorEngine contract onto the MathTS
+ * Tensor's methods. Both engines pass the identical engine-conformance
+ * suite, which is what guarantees behavioural parity across the two repos.
+ *
+ * @module numerical/mathts-engine
+ */
+import { Tensor } from '@danielsimonjr/mathts-tensor';
+import { NumericalBackendError } from './errors.js';
+import { EngineCapabilityError } from './tensor-engine.js';
+/** EngineTensor handle wrapping a MathTS Tensor. */
+class MathTSEngineTensor {
+    inner;
+    constructor(inner) {
+        this.inner = inner;
+    }
+    get shape() { return this.inner.shape; }
+}
+function unwrap(t, op) {
+    if (!(t instanceof MathTSEngineTensor)) {
+        throw new NumericalBackendError(`MathTSEngine.${op}: operand is not a MathTSEngineTensor`);
+    }
+    return t.inner;
+}
+/**
+ * `TensorEngine` backed by `@danielsimonjr/mathts-tensor`'s rank-N Tensor —
+ * the second v0.3.5 engine; becomes UPT's default in v0.4.0.
+ *
+ * @public — reachable only via the
+ * `universal-physics-tensor/numerical/mathts-engine` exports subpath; requires
+ * the `@danielsimonjr/mathts-tensor` optional dependency. Intentionally NOT
+ * re-exported from the root barrel.
+ */
+export class MathTSEngine {
+    name = 'MathTSEngine';
+    fromNested(data, shape) {
+        return new MathTSEngineTensor(Tensor.fromNested(data, shape));
+    }
+    toNested(t) {
+        return unwrap(t, 'toNested').toNested();
+    }
+    einsum(spec, ...operands) {
+        const inner = operands.map((o, i) => unwrap(o, `einsum (operand ${i})`));
+        return new MathTSEngineTensor(Tensor.einsum(spec, ...inner));
+    }
+    matMul(a, b) {
+        return new MathTSEngineTensor(unwrap(a, 'matMul').matMul(unwrap(b, 'matMul')));
+    }
+    transpose(t, perm) {
+        return new MathTSEngineTensor(unwrap(t, 'transpose').transpose(perm));
+    }
+    reshape(t, shape) {
+        return new MathTSEngineTensor(unwrap(t, 'reshape').reshape(shape));
+    }
+    add(a, b) {
+        return new MathTSEngineTensor(unwrap(a, 'add').add(unwrap(b, 'add')));
+    }
+    sub(a, b) {
+        return new MathTSEngineTensor(unwrap(a, 'sub').sub(unwrap(b, 'sub')));
+    }
+    mul(a, b) {
+        return new MathTSEngineTensor(unwrap(a, 'mul').mul(unwrap(b, 'mul')));
+    }
+    scale(t, k) {
+        return new MathTSEngineTensor(unwrap(t, 'scale').scale(k));
+    }
+    identity(n) {
+        return new MathTSEngineTensor(Tensor.identity(n));
+    }
+    normInf(t) {
+        return unwrap(t, 'normInf').normInf();
+    }
+    // ---------------------------------------------------------------------------
+    // Internal helpers for AD boundary conversions
+    // ---------------------------------------------------------------------------
+    /** Unwrap an EngineTensor to its underlying MathTS Tensor (for the AD boundary). */
+    toMathTSTensor(t) {
+        return unwrap(t, 'toMathTSTensor');
+    }
+    /** Wrap a MathTS Tensor (returned by autograd) back as an EngineTensor. */
+    fromMathTSTensor(t) {
+        return new MathTSEngineTensor(t);
+    }
+    // ---------------------------------------------------------------------------
+    // Forward-mode AD (Jacobian-vector product)
+    // ---------------------------------------------------------------------------
+    /**
+     * Forward-mode automatic differentiation via a lazy-imported
+     * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if the
+     * optional dependency is absent.
+     *
+     * S1 fix: `fn` is passed UNCHANGED to `autograd.forwardGrad`. The autograd
+     * package wraps `x` as a DualTensor internally; MathTSEngine's arithmetic
+     * methods (mul/add/sub/scale) must dispatch DualTensor inputs through to
+     * mathts-autograd's dual arithmetic — there must be no fn-wrapping at the
+     * boundary, which would strip DualTensor instrumentation and silence the
+     * AD trace.
+     */
+    async forwardGrad(fn, x) {
+        // eslint-disable-next-line @typescript-eslint/no-explicit-any
+        let autograd;
+        // @ts-ignore — @danielsimonjr/mathts-autograd is an optional dependency;
+        // tsc cannot resolve it when not installed. The try/catch handles absence
+        // at runtime; EngineCapabilityError is thrown when the package is absent.
+        try {
+            autograd = await import('@danielsimonjr/mathts-autograd');
+        }
+        catch {
+            throw new EngineCapabilityError('MathTSEngine', 'forwardGrad');
+        }
+        // S1 fix: pass fn UNCHANGED. autograd.forwardGrad wraps x as a DualTensor;
+        // MathTSEngine's arithmetic methods (mul/add/sub/scale) MUST dispatch
+        // DualTensor inputs to mathts-autograd's dual arithmetic (same dispatch
+        // story Float64ReferenceEngine has, via `'tangent' in arg`). Wrapping fn
+        // at the boundary strips the instrumentation — the v0 sketch's bug.
+        const xMathts = this.toMathTSTensor(x);
+        const { value, jacobian } = await autograd.forwardGrad(fn, xMathts);
+        return {
+            value: this.fromMathTSTensor(value),
+            jacobian: this.fromMathTSTensor(jacobian),
+        };
+    }
+    // ---------------------------------------------------------------------------
+    // Reverse-mode AD (vector-Jacobian product)
+    // ---------------------------------------------------------------------------
+    /**
+     * Reverse-mode automatic differentiation via a lazy-imported
+     * `@danielsimonjr/mathts-autograd`. Throws `EngineCapabilityError` if the
+     * optional dependency is absent.
+     *
+     * S1 fix: `fn` is passed UNCHANGED (see forwardGrad note above). The
+     * autograd package wraps `x` as a TapedTensor; MathTSEngine's op methods
+     * must dispatch TapedTensor inputs to mathts-autograd's tape arithmetic via
+     * `'tape' in arg` branching.
+     */
+    async reverseGrad(fn, x, cotangent) {
+        // eslint-disable-next-line @typescript-eslint/no-explicit-any
+        let autograd;
+        // @ts-ignore — @danielsimonjr/mathts-autograd is an optional dependency;
+        // tsc cannot resolve it when not installed. The try/catch handles absence
+        // at runtime; EngineCapabilityError is thrown when the package is absent.
+        try {
+            autograd = await import('@danielsimonjr/mathts-autograd');
+        }
+        catch {
+            throw new EngineCapabilityError('MathTSEngine', 'reverseGrad');
+        }
+        // S1 fix: pass fn UNCHANGED (see forwardGrad note above). autograd.reverseGrad
+        // wraps x as a TapedTensor; MathTSEngine's mul/add/sub/scale dispatch
+        // TapedTensor inputs to mathts-autograd's tape arithmetic via `'tape' in arg`.
+        const xMathts = this.toMathTSTensor(x);
+        const ctMathts = cotangent ? this.toMathTSTensor(cotangent) : undefined;
+        const { value, gradient } = await autograd.reverseGrad(fn, xMathts, ctMathts);
+        return {
+            value: this.fromMathTSTensor(value),
+            gradient: this.fromMathTSTensor(gradient),
+        };
+    }
+}
+//# sourceMappingURL=mathts-engine.js.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"mathts-engine.js","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,OAAO,EAAE,MAAM,EAAE,MAAM,8BAA8B,CAAC;AAGtD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,qBAAqB,EAAE,MAAM,oBAAoB,CAAC;AAE3D,oDAAoD;AACpD,MAAM,kBAAkB;IACD;IAArB,YAAqB,KAAa;QAAb,UAAK,GAAL,KAAK,CAAQ;IAAG,CAAC;IACtC,IAAI,KAAK,KAA4B,OAAO,IAAI,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC,CAAC;CAChE;AAED,SAAS,MAAM,CAAC,CAAe,EAAE,EAAU;IACzC,IAAI,CAAC,CAAC,CAAC,YAAY,kBAAkB,CAAC,EAAE,CAAC;QACvC,MAAM,IAAI,qBAAqB,CAAC,gBAAgB,EAAE,uCAAuC,CAAC,CAAC;IAC7F,CAAC;IACD,OAAO,CAAC,CAAC,KAAK,CAAC;AACjB,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,OAAO,YAAY;IACd,IAAI,GAAG,cAAc,CAAC;IAE/B,UAAU,CAAC,IAAiB,EAAE,KAA4B;QACxD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,CAAC;IAChE,CAAC;IACD,QAAQ,CAAC,CAAe;QACtB,OAAO,MAAM,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,QAAQ,EAAE,CAAC;IAC1C,CAAC;IAED,MAAM,CAAC,IAAgB,EAAE,GAAG,QAAwB;QAClD,MAAM,KAAK,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,EAAE,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC;QACzE,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC;IAC/D,CAAC;IACD,MAAM,CAAC,CAAe,EAAE,CAAe;QACrC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC;IACjF,CAAC;IACD,SAAS,CAAC,CAAe,EAAE,IAA4B;QACrD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,OAAO,CAAC,CAAe,EAAE,KAA4B;QACnD,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;IACrE,CAAC;IAED,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,GAAG,CAAC,CAAe,EAAE,CAAe;QAClC,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IACD,KAAK,CAAC,CAAe,EAAE,CAAS;QAC9B,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7D,CAAC;IAED,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,kBAAkB,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACpD,CAAC;IACD,OAAO,CAAC,CAAe;QACrB,OAAO,MAAM,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,OAAO,EAAE,CAAC;IACxC,CAAC;IAED,8EAA8E;IAC9E,+CAA+C;IAC/C,8EAA8E;IAE9E,oFAAoF;IAC5E,cAAc,CAAC,CAAe;QACpC,OAAO,MAAM,CAAC,CAAC,EAAE,gBAAgB,CAAC,CAAC;IACrC,CAAC;IAED,2EAA2E;IACnE,gBAAgB,CAAC,CAAS;QAChC,OAAO,IAAI,kBAAkB,CAAC,CAAC,CAAC,CAAC;IACnC,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;;;OAWG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe;QAEf,8DAA8D;QAC9D,IAAI,QAAa,CAAC;QAClB,yEAAyE;QACzE,0EAA0E;QAC1E,0EAA0E;QAC1E,IAAI,CAAC;YAAC,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAAC,CAAC;QAAC,CAAC;QAClE,MAAM,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,CAAC,CAAC;QAAC,CAAC;QAEzE,2EAA2E;QAC3E,sEAAsE;QACtE,wEAAwE;QACxE,yEAAyE;QACzE,oEAAoE;QACpE,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAS,EAAE,OAAc,CAAC,CAAC;QAClF,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;IAED,8EAA8E;IAC9E,4CAA4C;IAC5C,8EAA8E;IAE9E;;;;;;;;;OASG;IACH,KAAK,CAAC,WAAW,CACf,EAAqC,EACrC,CAAe,EACf,SAAwB;QAExB,8DAA8D;QAC9D,IAAI,QAAa,CAAC;QAClB,yEAAyE;QACzE,0EAA0E;QAC1E,0EAA0E;QAC1E,IAAI,CAAC;YAAC,QAAQ,GAAG,MAAM,MAAM,CAAC,gCAAgC,CAAC,CAAC;QAAC,CAAC;QAClE,MAAM,CAAC;YAAC,MAAM,IAAI,qBAAqB,CAAC,cAAc,EAAE,aAAa,CAAC,CAAC;QAAC,CAAC;QAEzE,+EAA+E;QAC/E,sEAAsE;QACtE,+EAA+E;QAC/E,MAAM,OAAO,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC;QACvC,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,cAAc,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,SAAS,CAAC;QACxE,MAAM,EAAE,KAAK,EAAE,QAAQ,EAAE,GAAG,MAAM,QAAQ,CAAC,WAAW,CAAC,EAAS,EAAE,OAAc,EAAE,QAAe,CAAC,CAAC;QACnG,OAAO;YACL,KAAK,EAAE,IAAI,CAAC,gBAAgB,CAAC,KAAK,CAAC;YACnC,QAAQ,EAAE,IAAI,CAAC,gBAAgB,CAAC,QAAQ,CAAC;SAC1C,CAAC;IACJ,CAAC;CACF"}

package/dist/numerical/metric-inverse.d.ts

@@ -0,0 +1,31 @@
+/**
+ * InverseMetricInconsistencyWarning — numerical path. Builds g⁻¹ and g as
+ * EngineTensors, computes ‖matMul(g⁻¹, g) − I‖_∞, and emits a
+ * 'warning'-severity Violation when it exceeds tolerance. v0.3.5-Design.md
+ * §7. Resolves the v0.3.0 §13 Q2 deferral. `scanForMetricPair` finds an
+ * identifiable lower/upper metric pair in an AST so evaluateNumerical can
+ * run the check automatically.
+ *
+ * @module numerical/metric-inverse
+ */
+import type { ExprNode, Violation } from '../dimensional/validator.js';
+import type { MetricTensorNode } from '../dimensional/metric-validators.js';
+import type { TensorEngine } from './tensor-engine.js';
+import type { NumericalInputs } from './types.js';
+/** Compute ‖g⁻¹g − I‖_∞ and, if it exceeds tolerance, return a warning.
+ *  @public */
+export declare function evaluateMetricInverse(gInverse: MetricTensorNode, g: MetricTensorNode, inputs: NumericalInputs, tolerance?: number, options?: {
+    engine?: TensorEngine;
+}): Promise<{
+    residualNorm: number;
+    warning?: Violation;
+}>;
+/** Walk an ExprNode tree; return the first all-lower / all-upper
+ *  metric-tensor pair found, or null. Cheap — evaluateNumerical already
+ *  walks the tree to lower it.
+ *  @internal — cross-module use only; not part of the consumer surface. */
+export declare function scanForMetricPair(node: ExprNode): {
+    gLower: MetricTensorNode;
+    gUpper: MetricTensorNode;
+} | null;
+//# sourceMappingURL=metric-inverse.d.ts.map

package/dist/numerical/metric-inverse.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"metric-inverse.d.ts","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AACH,OAAO,KAAK,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM,6BAA6B,CAAC;AACvE,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AAE5E,OAAO,KAAK,EAAE,YAAY,EAAE,MAAM,oBAAoB,CAAC;AACvD,OAAO,KAAK,EAAE,eAAe,EAAE,MAAM,YAAY,CAAC;AAMlD;cACc;AACd,wBAAsB,qBAAqB,CACzC,QAAQ,EAAE,gBAAgB,EAC1B,CAAC,EAAE,gBAAgB,EACnB,MAAM,EAAE,eAAe,EACvB,SAAS,GAAE,MAA0B,EACrC,OAAO,CAAC,EAAE;IAAE,MAAM,CAAC,EAAE,YAAY,CAAA;CAAE,GAClC,OAAO,CAAC;IAAE,YAAY,EAAE,MAAM,CAAC;IAAC,OAAO,CAAC,EAAE,SAAS,CAAA;CAAE,CAAC,CAkCxD;AAED;;;2EAG2E;AAC3E,wBAAgB,iBAAiB,CAC/B,IAAI,EAAE,QAAQ,GACb;IAAE,MAAM,EAAE,gBAAgB,CAAC;IAAC,MAAM,EAAE,gBAAgB,CAAA;CAAE,GAAG,IAAI,CAiB/D"}

package/dist/numerical/metric-inverse.js

@@ -0,0 +1,68 @@
+import { DIMENSIONLESS } from '../dimensional/types.js';
+import { getActiveEngine } from './engine-registry.js';
+import { NumericalBackendError } from './errors.js';
+const DEFAULT_TOLERANCE = 1e-10;
+/** Compute ‖g⁻¹g − I‖_∞ and, if it exceeds tolerance, return a warning.
+ *  @public */
+export async function evaluateMetricInverse(gInverse, g, inputs, tolerance = DEFAULT_TOLERANCE, options) {
+    const engine = options?.engine ?? await getActiveEngine();
+    const N = inputs.dimension ?? 4;
+    const gInvData = inputs.tensors.get(gInverse.name);
+    const gData = inputs.tensors.get(g.name);
+    if (gInvData === undefined) {
+        throw new NumericalBackendError(`evaluateMetricInverse: missing tensor components for "${gInverse.name}" in inputs.tensors`);
+    }
+    if (gData === undefined) {
+        throw new NumericalBackendError(`evaluateMetricInverse: missing tensor components for "${g.name}" in inputs.tensors`);
+    }
+    const gInvT = engine.fromNested(gInvData, [N, N]);
+    const gT = engine.fromNested(gData, [N, N]);
+    const residual = engine.sub(engine.matMul(gInvT, gT), engine.identity(N));
+    const residualNorm = engine.normInf(residual);
+    if (residualNorm > tolerance) {
+        const warning = {
+            location: `${gInverse.name}·${g.name}`,
+            expected: DIMENSIONLESS,
+            actual: DIMENSIONLESS,
+            note: `InverseMetricInconsistencyWarning: ‖${gInverse.name}·${g.name} − I‖_∞ = `
+                + `${residualNorm.toExponential(3)} exceeds tolerance ${tolerance.toExponential(1)} — `
+                + `the supplied metrics are not numerical inverses`,
+            severity: 'warning',
+        };
+        return { residualNorm, warning };
+    }
+    return { residualNorm };
+}
+/** Walk an ExprNode tree; return the first all-lower / all-upper
+ *  metric-tensor pair found, or null. Cheap — evaluateNumerical already
+ *  walks the tree to lower it.
+ *  @internal — cross-module use only; not part of the consumer surface. */
+export function scanForMetricPair(node) {
+    const metrics = [];
+    const walk = (n) => {
+        if (n.kind === 'metric-tensor')
+            metrics.push(n);
+        else if (n.kind === 'op')
+            n.args.forEach(walk);
+        else if (n.kind === 'tensor-product')
+            n.args.forEach(walk);
+        else if (n.kind === 'integral') {
+            walk(n.over);
+            walk(n.integrand);
+        }
+        else if (n.kind === 'derivative') {
+            walk(n.of);
+            walk(n.wrt);
+        }
+        else if (n.kind === 'tensor-partial-derivative') {
+            walk(n.of);
+            walk(n.wrt);
+        }
+        // tensor-symbol / kronecker-delta / symbol: leaves, no metric inside
+    };
+    walk(node);
+    const lower = metrics.find((m) => m.indices.every((i) => i.variance === 'lower'));
+    const upper = metrics.find((m) => m.indices.every((i) => i.variance === 'upper'));
+    return lower && upper ? { gLower: lower, gUpper: upper } : null;
+}
+//# sourceMappingURL=metric-inverse.js.map

package/dist/numerical/metric-inverse.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"metric-inverse.js","sourceRoot":"","sources":["../../src/numerical/metric-inverse.ts"],"names":[],"mappings":"AAYA,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AAGxD,OAAO,EAAE,eAAe,EAAE,MAAM,sBAAsB,CAAC;AACvD,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEpD,MAAM,iBAAiB,GAAG,KAAK,CAAC;AAEhC;cACc;AACd,MAAM,CAAC,KAAK,UAAU,qBAAqB,CACzC,QAA0B,EAC1B,CAAmB,EACnB,MAAuB,EACvB,YAAoB,iBAAiB,EACrC,OAAmC;IAEnC,MAAM,MAAM,GAAG,OAAO,EAAE,MAAM,IAAI,MAAM,eAAe,EAAE,CAAC;IAC1D,MAAM,CAAC,GAAG,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;IAEhC,MAAM,QAAQ,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IACnD,MAAM,KAAK,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACzC,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;QAC3B,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,QAAQ,CAAC,IAAI,qBAAqB,CAC5F,CAAC;IACJ,CAAC;IACD,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;QACxB,MAAM,IAAI,qBAAqB,CAC7B,yDAAyD,CAAC,CAAC,IAAI,qBAAqB,CACrF,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAClD,MAAM,EAAE,GAAG,MAAM,CAAC,UAAU,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IAC5C,MAAM,QAAQ,GAAG,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1E,MAAM,YAAY,GAAG,MAAM,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC;IAE9C,IAAI,YAAY,GAAG,SAAS,EAAE,CAAC;QAC7B,MAAM,OAAO,GAAc;YACzB,QAAQ,EAAE,GAAG,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,EAAE;YACtC,QAAQ,EAAE,aAAa;YACvB,MAAM,EAAE,aAAa;YACrB,IAAI,EAAE,uCAAuC,QAAQ,CAAC,IAAI,IAAI,CAAC,CAAC,IAAI,YAAY;kBAC5E,GAAG,YAAY,CAAC,aAAa,CAAC,CAAC,CAAC,sBAAsB,SAAS,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK;kBACrF,iDAAiD;YACrD,QAAQ,EAAE,SAAS;SACpB,CAAC;QACF,OAAO,EAAE,YAAY,EAAE,OAAO,EAAE,CAAC;IACnC,CAAC;IACD,OAAO,EAAE,YAAY,EAAE,CAAC;AAC1B,CAAC;AAED;;;2EAG2E;AAC3E,MAAM,UAAU,iBAAiB,CAC/B,IAAc;IAEd,MAAM,OAAO,GAAuB,EAAE,CAAC;IACvC,MAAM,IAAI,GAAG,CAAC,CAAW,EAAQ,EAAE;QACjC,IAAI,CAAC,CAAC,IAAI,KAAK,eAAe;YAAE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;aAC3C,IAAI,CAAC,CAAC,IAAI,KAAK,IAAI;YAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;aAC1C,IAAI,CAAC,CAAC,IAAI,KAAK,gBAAgB;YAAE,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;aACtD,IAAI,CAAC,CAAC,IAAI,KAAK,UAAU,EAAE,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC;QAAC,CAAC;aAC/D,IAAI,CAAC,CAAC,IAAI,KAAK,YAAY,EAAE,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAAC,CAAC;aACzD,IAAI,CAAC,CAAC,IAAI,KAAK,2BAA2B,EAAE,CAAC;YAChD,IAAI,CAAC,CAAC,CAAC,EAAc,CAAC,CAAC;YAAC,IAAI,CAAC,CAAC,CAAC,GAAe,CAAC,CAAC;QAClD,CAAC;QACD,qEAAqE;IACvE,CAAC,CAAC;IACF,IAAI,CAAC,IAAI,CAAC,CAAC;IACX,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,MAAM,KAAK,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC;IAClF,OAAO,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC;AAClE,CAAC"}

package/dist/numerical/null-ray-integrator.d.ts

@@ -0,0 +1,13 @@
+/** A first-order ODE system: dy/dλ = f(λ, y). `y` and the return are
+ *  state vectors of equal length.
+ *  @internal — cross-module/test use only; not part of the consumer surface. */
+export type ODESystem = (lambda: number, y: ReadonlyArray<number>) => number[];
+/**
+ * Integrate `system` from affine parameter `lambda0` to `lambda1` in
+ * `steps` fixed RK4 steps, starting from state `y0`. Returns the final
+ * state vector. Classical 4th-order Runge-Kutta — global error O(h⁴).
+ *
+ * @internal — cross-module/test use only; not part of the consumer surface.
+ */
+export declare function integrateRK4(system: ODESystem, y0: ReadonlyArray<number>, lambda0: number, lambda1: number, steps: number): number[];
+//# sourceMappingURL=null-ray-integrator.d.ts.map

package/dist/numerical/null-ray-integrator.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"null-ray-integrator.d.ts","sourceRoot":"","sources":["../../src/numerical/null-ray-integrator.ts"],"names":[],"mappings":"AAUA;;gFAEgF;AAChF,MAAM,MAAM,SAAS,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,MAAM,EAAE,CAAC;AAgB/E;;;;;;GAMG;AACH,wBAAgB,YAAY,CAC1B,MAAM,EAAE,SAAS,EACjB,EAAE,EAAE,aAAa,CAAC,MAAM,CAAC,EACzB,OAAO,EAAE,MAAM,EACf,OAAO,EAAE,MAAM,EACf,KAAK,EAAE,MAAM,GACZ,MAAM,EAAE,CAyBV"}

package/dist/numerical/null-ray-integrator.js

@@ -0,0 +1,53 @@
+/**
+ * Fixed-step classical RK4 integrator for affine-parameterized null
+ * geodesics — the numerical core of the BE-37 Shapiro-delay evaluator
+ * (v0.3.5-Design.md §8). Self-contained: operates on plain number[] state
+ * vectors, no TensorEngine dependency.
+ *
+ * @module numerical/null-ray-integrator
+ */
+import { NumericalBackendError } from './errors.js';
+function addScaled(a, b, k) {
+    const out = new Array(a.length);
+    for (let i = 0; i < a.length; i++)
+        out[i] = a[i] + k * b[i];
+    return out;
+}
+function checkLength(k, expected) {
+    if (k.length !== expected) {
+        throw new NumericalBackendError(`integrateRK4: ODE system returned a state vector of length ${k.length}, expected ${expected}`);
+    }
+}
+/**
+ * Integrate `system` from affine parameter `lambda0` to `lambda1` in
+ * `steps` fixed RK4 steps, starting from state `y0`. Returns the final
+ * state vector. Classical 4th-order Runge-Kutta — global error O(h⁴).
+ *
+ * @internal — cross-module/test use only; not part of the consumer surface.
+ */
+export function integrateRK4(system, y0, lambda0, lambda1, steps) {
+    if (!Number.isInteger(steps) || steps <= 0) {
+        throw new NumericalBackendError(`integrateRK4: step count must be a positive integer, got ${steps}`);
+    }
+    const h = (lambda1 - lambda0) / steps;
+    let lambda = lambda0;
+    let y = [...y0];
+    for (let n = 0; n < steps; n++) {
+        const k1 = system(lambda, y);
+        checkLength(k1, y.length);
+        const k2 = system(lambda + h / 2, addScaled(y, k1, h / 2));
+        checkLength(k2, y.length);
+        const k3 = system(lambda + h / 2, addScaled(y, k2, h / 2));
+        checkLength(k3, y.length);
+        const k4 = system(lambda + h, addScaled(y, k3, h));
+        checkLength(k4, y.length);
+        const next = new Array(y.length);
+        for (let i = 0; i < y.length; i++) {
+            next[i] = y[i] + (h / 6) * (k1[i] + 2 * k2[i] + 2 * k3[i] + k4[i]);
+        }
+        y = next;
+        lambda += h;
+    }
+    return y;
+}
+//# sourceMappingURL=null-ray-integrator.js.map

package/dist/numerical/null-ray-integrator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"null-ray-integrator.js","sourceRoot":"","sources":["../../src/numerical/null-ray-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AACH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAOpD,SAAS,SAAS,CAAC,CAAwB,EAAE,CAAwB,EAAE,CAAS;IAC9E,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,MAAM,CAAC,CAAC;IACxC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,OAAO,GAAG,CAAC;AACb,CAAC;AAED,SAAS,WAAW,CAAC,CAAwB,EAAE,QAAgB;IAC7D,IAAI,CAAC,CAAC,MAAM,KAAK,QAAQ,EAAE,CAAC;QAC1B,MAAM,IAAI,qBAAqB,CAC7B,8DAA8D,CAAC,CAAC,MAAM,cAAc,QAAQ,EAAE,CAC/F,CAAC;IACJ,CAAC;AACH,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,YAAY,CAC1B,MAAiB,EACjB,EAAyB,EACzB,OAAe,EACf,OAAe,EACf,KAAa;IAEb,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QAC3C,MAAM,IAAI,qBAAqB,CAAC,4DAA4D,KAAK,EAAE,CAAC,CAAC;IACvG,CAAC;IACD,MAAM,CAAC,GAAG,CAAC,OAAO,GAAG,OAAO,CAAC,GAAG,KAAK,CAAC;IACtC,IAAI,MAAM,GAAG,OAAO,CAAC;IACrB,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;IAEhB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;QAC7B,WAAW,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC3D,WAAW,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,EAAE,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC3D,WAAW,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;QAC1B,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,SAAS,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC;QACnD,WAAW,CAAC,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC;QAC1B,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,MAAM,CAAC,CAAC;QACzC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,CAAC;QACD,CAAC,GAAG,IAAI,CAAC;QACT,MAAM,IAAI,CAAC,CAAC;IACd,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/dist/numerical/pderiv.d.ts

@@ -0,0 +1,43 @@
+/**
+ * Numerical partial derivative — two-way dispatch (v0.3.5-Design.md §6).
+ *   'grid'         — centered finite-difference over a GridField.
+ *   'numerical-fn' — centered finite-difference over a caller-supplied fn.
+ *   'symbolic'     — no CAS in v0.3.5; consumes an explicit pre-supplied
+ *                    derivative from NumericalInputs.derivatives.
+ *
+ * @module numerical/pderiv
+ */
+import type { GridField } from './grid-field.js';
+import type { NestedArray } from './types.js';
+/**
+ * Centered finite-difference of a GridField along `axis`. Returns a flat
+ * number[] of the same total size as the grid. O(h²) interior stencil;
+ * edges follow `grid.boundary`:
+ *   - 'clamp'    — one-sided (forward/backward) difference at the edge.
+ *   - 'periodic' — wraps to the opposite edge.
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export declare function pderivGrid(grid: GridField, axis: number): number[];
+/**
+ * Centered finite-difference of a caller-supplied scalar field with respect
+ * to coordinate `axis`. Step h = 1e-6 · max(|x|, 1) (v0.3.5-Design.md §13 Q3).
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export declare function pderivNumericalFn(fn: (coords: ReadonlyArray<number>) => NestedArray, coords: ReadonlyArray<number>, axis: number): NestedArray;
+/**
+ * 'symbolic' numericalForm path: v0.3.5 has no CAS, so the caller supplies
+ * the derivative explicitly, keyed `${symbolName}/${coordLabel}`.
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export declare function pderivSymbolic(symbolName: string, coordLabel: string, derivatives: ReadonlyMap<string, NestedArray>): NestedArray;
+/**
+ * 'supplied' derivativeStrategy path: v0.4.0 metric derivative looked up
+ * from explicit caller-supplied components. Keyed `${metricName}/${coordLabel}`.
+ * Mirror of `pderivSymbolic` for the metric-tensor lowering path.
+ * Throws NumericalBackendError with a clear message when absent.
+ */
+export declare function metricDerivSupplied(metricName: string, coordLabel: string, metricDerivatives: ReadonlyMap<string, NestedArray>): NestedArray;
+//# sourceMappingURL=pderiv.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"pderiv.d.ts","sourceRoot":"","sources":["../../src/numerical/pderiv.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AACH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AACjD,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAe9C;;;;;;;;GAQG;AACH,wBAAgB,UAAU,CAAC,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,MAAM,GAAG,MAAM,EAAE,CAiClE;AAED;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAC/B,EAAE,EAAE,CAAC,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,KAAK,WAAW,EAClD,MAAM,EAAE,aAAa,CAAC,MAAM,CAAC,EAC7B,IAAI,EAAE,MAAM,GACX,WAAW,CAeb;AAED;;;;;GAKG;AACH,wBAAgB,cAAc,CAC5B,UAAU,EAAE,MAAM,EAClB,UAAU,EAAE,MAAM,EAClB,WAAW,EAAE,WAAW,CAAC,MAAM,EAAE,WAAW,CAAC,GAC5C,WAAW,CAWb;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CACjC,UAAU,EAAE,MAAM,EAClB,UAAU,EAAE,MAAM,EAClB,iBAAiB,EAAE,WAAW,CAAC,MAAM,EAAE,WAAW,CAAC,GAClD,WAAW,CAYb"}

package/dist/numerical/pderiv.js

@@ -0,0 +1,121 @@
+import { NumericalBackendError } from './errors.js';
+/** Flatten a (possibly scalar) NestedArray to a number[] in row-major order. */
+function flattenToNumbers(data) {
+    if (typeof data === 'number')
+        return [data];
+    const out = [];
+    const walk = (n) => {
+        if (typeof n === 'number')
+            out.push(n);
+        else
+            for (const c of n)
+                walk(c);
+    };
+    walk(data);
+    return out;
+}
+/**
+ * Centered finite-difference of a GridField along `axis`. Returns a flat
+ * number[] of the same total size as the grid. O(h²) interior stencil;
+ * edges follow `grid.boundary`:
+ *   - 'clamp'    — one-sided (forward/backward) difference at the edge.
+ *   - 'periodic' — wraps to the opposite edge.
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export function pderivGrid(grid, axis) {
+    if (axis < 0 || axis >= grid.shape.length) {
+        throw new NumericalBackendError(`pderivGrid: axis ${axis} out of range for shape [${grid.shape}]`);
+    }
+    const h = grid.spacing[axis];
+    if (!(h > 0)) {
+        throw new NumericalBackendError(`pderivGrid: non-positive spacing ${h} on axis ${axis}`);
+    }
+    const flat = flattenToNumbers(grid.data);
+    const shape = grid.shape;
+    const strides = new Array(shape.length);
+    let acc = 1;
+    for (let k = shape.length - 1; k >= 0; k--) {
+        strides[k] = acc;
+        acc *= shape[k];
+    }
+    const out = new Array(flat.length).fill(0);
+    const n = shape[axis];
+    const step = strides[axis];
+    // Walk every multi-index; compute the axis position from the flat index.
+    for (let f = 0; f < flat.length; f++) {
+        const pos = Math.floor(f / step) % n;
+        if (pos > 0 && pos < n - 1) {
+            out[f] = (flat[f + step] - flat[f - step]) / (2 * h);
+        }
+        else if (grid.boundary === 'periodic') {
+            const next = pos === n - 1 ? f - (n - 1) * step : f + step;
+            const prev = pos === 0 ? f + (n - 1) * step : f - step;
+            out[f] = (flat[next] - flat[prev]) / (2 * h);
+        }
+        else { // 'clamp' — one-sided difference at the edge
+            out[f] = pos === 0
+                ? (flat[f + step] - flat[f]) / h
+                : (flat[f] - flat[f - step]) / h;
+        }
+    }
+    return out;
+}
+/**
+ * Centered finite-difference of a caller-supplied scalar field with respect
+ * to coordinate `axis`. Step h = 1e-6 · max(|x|, 1) (v0.3.5-Design.md §13 Q3).
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export function pderivNumericalFn(fn, coords, axis) {
+    if (axis < 0 || axis >= coords.length) {
+        throw new NumericalBackendError(`pderivNumericalFn: axis ${axis} out of range`);
+    }
+    const x = coords[axis];
+    const h = 1e-6 * Math.max(Math.abs(x), 1);
+    const plus = [...coords];
+    plus[axis] = x + h;
+    const minus = [...coords];
+    minus[axis] = x - h;
+    const fp = flattenToNumbers(fn(plus));
+    const fm = flattenToNumbers(fn(minus));
+    if (fp.length !== fm.length) {
+        throw new NumericalBackendError('pderivNumericalFn: field returned inconsistent shapes');
+    }
+    const d = fp.map((v, i) => (v - fm[i]) / (2 * h));
+    return d.length === 1 ? d[0] : d;
+}
+/**
+ * 'symbolic' numericalForm path: v0.3.5 has no CAS, so the caller supplies
+ * the derivative explicitly, keyed `${symbolName}/${coordLabel}`.
+ *
+ * @internal — consumed by the lowering pass; not part of the consumer surface.
+ */
+export function pderivSymbolic(symbolName, coordLabel, derivatives) {
+    const key = `${symbolName}/${coordLabel}`;
+    const d = derivatives.get(key);
+    if (d === undefined) {
+        throw new NumericalBackendError(`pderivSymbolic: no explicit derivative supplied for "${key}" — `
+            + `a 'symbolic' tensor-symbol under a partial-derivative requires `
+            + `inputs.derivatives to contain its pre-computed components (v0.3.5 has no CAS)`);
+    }
+    return d;
+}
+/**
+ * 'supplied' derivativeStrategy path: v0.4.0 metric derivative looked up
+ * from explicit caller-supplied components. Keyed `${metricName}/${coordLabel}`.
+ * Mirror of `pderivSymbolic` for the metric-tensor lowering path.
+ * Throws NumericalBackendError with a clear message when absent.
+ */
+export function metricDerivSupplied(metricName, coordLabel, metricDerivatives) {
+    const key = `${metricName}/${coordLabel}`;
+    const d = metricDerivatives.get(key);
+    if (d === undefined) {
+        throw new NumericalBackendError(`metricDerivSupplied: no metric derivative supplied for "${key}" — `
+            + `a metric-tensor with derivativeStrategy='supplied' under a `
+            + `christoffel/covariant-derivative requires inputs.metricDerivatives `
+            + `to contain its pre-computed ∂g components`);
+    }
+    return d;
+}
+//# sourceMappingURL=pderiv.js.map

package/dist/numerical/pderiv.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"pderiv.js","sourceRoot":"","sources":["../../src/numerical/pderiv.ts"],"names":[],"mappings":"AAWA,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEpD,gFAAgF;AAChF,SAAS,gBAAgB,CAAC,IAAiB;IACzC,IAAI,OAAO,IAAI,KAAK,QAAQ;QAAE,OAAO,CAAC,IAAI,CAAC,CAAC;IAC5C,MAAM,GAAG,GAAa,EAAE,CAAC;IACzB,MAAM,IAAI,GAAG,CAAC,CAAc,EAAQ,EAAE;QACpC,IAAI,OAAO,CAAC,KAAK,QAAQ;YAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;;YAClC,KAAK,MAAM,CAAC,IAAI,CAAC;gBAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC;IACF,IAAI,CAAC,IAAI,CAAC,CAAC;IACX,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,UAAU,CAAC,IAAe,EAAE,IAAY;IACtD,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,IAAI,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QAC1C,MAAM,IAAI,qBAAqB,CAAC,oBAAoB,IAAI,4BAA4B,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC;IACrG,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC;IAC7B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACb,MAAM,IAAI,qBAAqB,CAAC,oCAAoC,CAAC,YAAY,IAAI,EAAE,CAAC,CAAC;IAC3F,CAAC;IACD,MAAM,IAAI,GAAG,gBAAgB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC;IACzB,MAAM,OAAO,GAAa,IAAI,KAAK,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IAClD,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,KAAK,IAAI,CAAC,GAAG,KAAK,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAAC,OAAO,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAAC,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC;IAClF,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;IACtB,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,CAAC;IAE3B,yEAAyE;IACzE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACrC,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;QACrC,IAAI,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC;YAC3B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QACvD,CAAC;aAAM,IAAI,IAAI,CAAC,QAAQ,KAAK,UAAU,EAAE,CAAC;YACxC,MAAM,IAAI,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YAC3D,MAAM,IAAI,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YACvD,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/C,CAAC;aAAM,CAAC,CAAC,6CAA6C;YACpD,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,KAAK,CAAC;gBAChB,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;gBAChC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACrC,CAAC;IACH,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,iBAAiB,CAC/B,EAAkD,EAClD,MAA6B,EAC7B,IAAY;IAEZ,IAAI,IAAI,GAAG,CAAC,IAAI,IAAI,IAAI,MAAM,CAAC,MAAM,EAAE,CAAC;QACtC,MAAM,IAAI,qBAAqB,CAAC,2BAA2B,IAAI,eAAe,CAAC,CAAC;IAClF,CAAC;IACD,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC;IACvB,MAAM,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;IAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC7C,MAAM,KAAK,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;IAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC/C,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;IACtC,MAAM,EAAE,GAAG,gBAAgB,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC;IACvC,IAAI,EAAE,CAAC,MAAM,KAAK,EAAE,CAAC,MAAM,EAAE,CAAC;QAC5B,MAAM,IAAI,qBAAqB,CAAC,uDAAuD,CAAC,CAAC;IAC3F,CAAC;IACD,MAAM,CAAC,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAClD,OAAO,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACnC,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAC5B,UAAkB,EAClB,UAAkB,EAClB,WAA6C;IAE7C,MAAM,GAAG,GAAG,GAAG,UAAU,IAAI,UAAU,EAAE,CAAC;IAC1C,MAAM,CAAC,GAAG,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC/B,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;QACpB,MAAM,IAAI,qBAAqB,CAC7B,wDAAwD,GAAG,MAAM;cAC/D,iEAAiE;cACjE,+EAA+E,CAClF,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,mBAAmB,CACjC,UAAkB,EAClB,UAAkB,EAClB,iBAAmD;IAEnD,MAAM,GAAG,GAAG,GAAG,UAAU,IAAI,UAAU,EAAE,CAAC;IAC1C,MAAM,CAAC,GAAG,iBAAiB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACrC,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;QACpB,MAAM,IAAI,qBAAqB,CAC7B,2DAA2D,GAAG,MAAM;cAClE,6DAA6D;cAC7D,qEAAqE;cACrE,2CAA2C,CAC9C,CAAC;IACJ,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}