universal-physics-tensor 0.5.1 → 0.7.0

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

@@ -0,0 +1 @@
+{"version":3,"file":"kretschmann.d.ts","sourceRoot":"","sources":["../../src/numerical/kretschmann.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AACH,wBAAgB,kBAAkB,CAChC,YAAY,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,EAC5B,aAAa,EAAE,MAAM,EAAE,EAAE,GACxB,MAAM,CA0BR"}

package/dist/numerical/kretschmann.js

@@ -0,0 +1,88 @@
+/**
+ * Kretschmann scalar numerical contraction (v0.6.0 Phase 3, Task 3.6).
+ *
+ * Computes K = R_{ρσμν} R^{ρσμν} from the fully-lowered Riemann tensor and
+ * the inverse metric by raising all four indices on the second Riemann factor.
+ *
+ * **PERFORMANCE WARNING**: O(4⁸) = 65536 multiplications per call. Acceptable
+ * for diagnostic / sample-point use, but DO NOT call inside tight integration
+ * loops (10⁶ steps × 65k ops = ~6.5e10 ops). Future optimization (v0.7.0
+ * candidate): exploit Riemann symmetries to reduce to ~20 independent
+ * components × symmetric weighting; estimated speedup ~12×. Tracked in
+ * v0.6.0 brainstorm carry-forward.
+ *
+ * @module numerical/kretschmann
+ */
+/**
+ * Kretschmann invariant K = R_{ρσμν} R^{ρσμν}.
+ *
+ * Algorithm:
+ *   1. Raise all four indices on the second Riemann factor:
+ *      R^{ρσμν} = g^{ρα} g^{σβ} g^{μγ} g^{νδ} R_{αβγδ}
+ *   2. Contract:
+ *      K = Σ_{ρσμν} R_{ρσμν} · R^{ρσμν}
+ *
+ * The outer sum has 4⁴ = 256 terms; each inner raise requires 4⁴ = 256
+ * metric products, giving the O(4⁸) = 65536 total multiplications.
+ *
+ * @param riemannLower - R_{ρσμν}: all-lower Riemann tensor, shape [4][4][4][4].
+ *   Obtain by lowering R^ρ_{σμν} via g_{ρα}: R_{αβγδ} = g_{αρ} R^ρ_{βγδ}.
+ * @param metricInverse - g^{αβ}: inverse metric, shape [4][4].
+ * @returns K, a scalar. Dimension [L⁻⁴] (Riemann is [L⁻²] per F8/I3 convention).
+ *
+ * @example
+ * ```typescript
+ * import { computeKretschmann, G_SI, C_SI } from 'universal-physics-tensor';
+ * import {
+ *   schwarzschildGFn,
+ *   schwarzschildGInverseFn,
+ *   schwarzschildRs,
+ * } from '../tests/fixtures/schwarzschild.js';
+ * import { riemannLowerAt } from '../src/numerical/curvature-lowering-helpers.js';
+ * import { Float64ReferenceEngine } from '../src/numerical/float64-engine.js';
+ *
+ * const M = 1.989e30;
+ * const r_s = schwarzschildRs(M);
+ * const r = 5 * r_s;
+ * const x: [number, number, number, number] = [0, r, Math.PI / 2, 0];
+ *
+ * const engine = new Float64ReferenceEngine();
+ * const rLower = riemannLowerAt(x, schwarzschildGFn(M), schwarzschildGInverseFn(M), 4, engine);
+ * const gInv = schwarzschildGInverseFn(M)(x);
+ * const K = computeKretschmann(rLower, gInv);
+ *
+ * // Schwarzschild closed-form: K = 48 G² M² / (c⁴ r⁶)
+ * const K_analytic = 48 * G_SI**2 * M**2 / (C_SI**4 * r**6);
+ * // |K - K_analytic| / K_analytic < 1e-4
+ * ```
+ *
+ * @public
+ */
+export function computeKretschmann(riemannLower, metricInverse) {
+    // Raise: R^{ρσμν} = g^{ρα} g^{σβ} g^{μγ} g^{νδ} R_{αβγδ}
+    // Then K = R_{ρσμν} R^{ρσμν} (sum over all 4⁴ = 256 components)
+    let K = 0;
+    for (let rho = 0; rho < 4; rho++) {
+        for (let sigma = 0; sigma < 4; sigma++) {
+            for (let mu = 0; mu < 4; mu++) {
+                for (let nu = 0; nu < 4; nu++) {
+                    let raised = 0;
+                    for (let a = 0; a < 4; a++) {
+                        for (let b = 0; b < 4; b++) {
+                            for (let c = 0; c < 4; c++) {
+                                for (let d = 0; d < 4; d++) {
+                                    raised += metricInverse[rho][a] * metricInverse[sigma][b]
+                                        * metricInverse[mu][c] * metricInverse[nu][d]
+                                        * riemannLower[a][b][c][d];
+                                }
+                            }
+                        }
+                    }
+                    K += riemannLower[rho][sigma][mu][nu] * raised;
+                }
+            }
+        }
+    }
+    return K;
+}
+//# sourceMappingURL=kretschmann.js.map

package/dist/numerical/kretschmann.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"kretschmann.js","sourceRoot":"","sources":["../../src/numerical/kretschmann.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4CG;AACH,MAAM,UAAU,kBAAkB,CAChC,YAA4B,EAC5B,aAAyB;IAEzB,yDAAyD;IACzD,gEAAgE;IAChE,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC;YACvC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;oBAC9B,IAAI,MAAM,GAAG,CAAC,CAAC;oBACf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;wBAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;4BAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gCAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;oCAC3B,MAAM,IAAI,aAAa,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;0CAC/C,aAAa,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,aAAa,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;0CAC3C,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gCACrC,CAAC;4BACH,CAAC;wBACH,CAAC;oBACH,CAAC;oBACD,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC;gBACjD,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC"}

package/dist/numerical/lowering-utils.d.ts

@@ -0,0 +1,47 @@
+/**
+ * Shared private utilities for the `numerical/lowering*.ts` modules.
+ *
+ * Carved out of `lowering.ts` in the v0.7 follow-up to its v0.6.1
+ * sprint-target miss (LOC 903 vs ≤890). Centralizes the four helpers
+ * that both `lowering.ts` and the v0.7 `derivative-lowering.ts`
+ * need, avoiding both a forward-import cycle and code duplication.
+ *
+ * No public surface — every export carries an `@internal` JSDoc tag.
+ *
+ * @module numerical/lowering-utils
+ */
+import type { NumericalInputs, NestedArray } from './types.js';
+import type { MetricTensorNode } from '../dimensional/metric-validators.js';
+/**
+ * Module-private type predicate for the `metric-tensor` AST kind.
+ * Replaces the duplicated `(n as { kind?: unknown }).kind ===
+ * 'metric-tensor'` cast pattern. TypeScript narrows the operand to
+ * `MetricTensorNode` after the predicate, removing follow-up `as
+ * MetricTensorNode` casts.
+ *
+ * @internal
+ */
+export declare function isMetricTensorNode(n: unknown): n is MetricTensorNode;
+/**
+ * Coordinate-space dimension supplied via `NumericalInputs.dimension`,
+ * defaulting to 4 (spacetime). Read at most case-arms in lowering.
+ *
+ * @internal
+ */
+export declare function dimensionOf(inputs: NumericalInputs): number;
+/**
+ * Look up a named tensor's concrete value from
+ * `NumericalInputs.tensors`, or throw with a clear message.
+ *
+ * @internal
+ */
+export declare function requireValue(name: string, inputs: NumericalInputs): NestedArray;
+/**
+ * Flatten a `NestedArray` to a plain `number[]` and check the result
+ * matches `expectedSize`. Wraps the canonical `flattenNA()` with a
+ * size-mismatch guard so consumers don't silently swallow bad input.
+ *
+ * @internal
+ */
+export declare function flattenNestedArray(data: NestedArray, expectedSize: number): number[];
+//# sourceMappingURL=lowering-utils.d.ts.map

package/dist/numerical/lowering-utils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"lowering-utils.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering-utils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,eAAe,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAC/D,OAAO,KAAK,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AAI5E;;;;;;;;GAQG;AACH,wBAAgB,kBAAkB,CAAC,CAAC,EAAE,OAAO,GAAG,CAAC,IAAI,gBAAgB,CAGpE;AAED;;;;;GAKG;AACH,wBAAgB,WAAW,CAAC,MAAM,EAAE,eAAe,GAAG,MAAM,CAE3D;AAED;;;;;GAKG;AACH,wBAAgB,YAAY,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,eAAe,GAAG,WAAW,CAM/E;AAED;;;;;;GAMG;AACH,wBAAgB,kBAAkB,CAAC,IAAI,EAAE,WAAW,EAAE,YAAY,EAAE,MAAM,GAAG,MAAM,EAAE,CAQpF"}

package/dist/numerical/lowering-utils.js

@@ -0,0 +1,64 @@
+/**
+ * Shared private utilities for the `numerical/lowering*.ts` modules.
+ *
+ * Carved out of `lowering.ts` in the v0.7 follow-up to its v0.6.1
+ * sprint-target miss (LOC 903 vs ≤890). Centralizes the four helpers
+ * that both `lowering.ts` and the v0.7 `derivative-lowering.ts`
+ * need, avoiding both a forward-import cycle and code duplication.
+ *
+ * No public surface — every export carries an `@internal` JSDoc tag.
+ *
+ * @module numerical/lowering-utils
+ */
+import { NumericalBackendError } from './errors.js';
+import { flattenNA } from './connection-lowering-helpers.js';
+/**
+ * Module-private type predicate for the `metric-tensor` AST kind.
+ * Replaces the duplicated `(n as { kind?: unknown }).kind ===
+ * 'metric-tensor'` cast pattern. TypeScript narrows the operand to
+ * `MetricTensorNode` after the predicate, removing follow-up `as
+ * MetricTensorNode` casts.
+ *
+ * @internal
+ */
+export function isMetricTensorNode(n) {
+    return typeof n === 'object' && n !== null
+        && n.kind === 'metric-tensor';
+}
+/**
+ * Coordinate-space dimension supplied via `NumericalInputs.dimension`,
+ * defaulting to 4 (spacetime). Read at most case-arms in lowering.
+ *
+ * @internal
+ */
+export function dimensionOf(inputs) {
+    return inputs.dimension ?? 4;
+}
+/**
+ * Look up a named tensor's concrete value from
+ * `NumericalInputs.tensors`, or throw with a clear message.
+ *
+ * @internal
+ */
+export function requireValue(name, inputs) {
+    const v = inputs.tensors.get(name);
+    if (v === undefined) {
+        throw new NumericalBackendError(`lowering: no value supplied for "${name}" in inputs.tensors`);
+    }
+    return v;
+}
+/**
+ * Flatten a `NestedArray` to a plain `number[]` and check the result
+ * matches `expectedSize`. Wraps the canonical `flattenNA()` with a
+ * size-mismatch guard so consumers don't silently swallow bad input.
+ *
+ * @internal
+ */
+export function flattenNestedArray(data, expectedSize) {
+    const out = flattenNA(data);
+    if (out.length !== expectedSize) {
+        throw new NumericalBackendError(`lowering: flattenNestedArray: got ${out.length} elements, expected ${expectedSize}`);
+    }
+    return out;
+}
+//# sourceMappingURL=lowering-utils.js.map

package/dist/numerical/lowering-utils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lowering-utils.js","sourceRoot":"","sources":["../../src/numerical/lowering-utils.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EAAE,SAAS,EAAE,MAAM,kCAAkC,CAAC;AAE7D;;;;;;;;GAQG;AACH,MAAM,UAAU,kBAAkB,CAAC,CAAU;IAC3C,OAAO,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,KAAK,IAAI;WACpC,CAAwB,CAAC,IAAI,KAAK,eAAe,CAAC;AAC1D,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,WAAW,CAAC,MAAuB;IACjD,OAAO,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,YAAY,CAAC,IAAY,EAAE,MAAuB;IAChE,MAAM,CAAC,GAAG,MAAM,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;IACnC,IAAI,CAAC,KAAK,SAAS,EAAE,CAAC;QACpB,MAAM,IAAI,qBAAqB,CAAC,oCAAoC,IAAI,qBAAqB,CAAC,CAAC;IACjG,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,kBAAkB,CAAC,IAAiB,EAAE,YAAoB;IACxE,MAAM,GAAG,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;IAC5B,IAAI,GAAG,CAAC,MAAM,KAAK,YAAY,EAAE,CAAC;QAChC,MAAM,IAAI,qBAAqB,CAC7B,qCAAqC,GAAG,CAAC,MAAM,uBAAuB,YAAY,EAAE,CACrF,CAAC;IACJ,CAAC;IACD,OAAO,GAAG,CAAC;AACb,CAAC"}

package/dist/numerical/lowering.d.ts

@@ -11,37 +11,9 @@
  * @module numerical/lowering
  */
 import type { ExprNode } from '../dimensional/validator.js';
-import type { TensorIndex } from '../dimensional/tensor.js';
-import type { EngineTensor, TensorEngine, EinsumSpec } from './tensor-engine.js';
+import type { EngineTensor, TensorEngine } from './tensor-engine.js';
 import type { NumericalInputs } from './types.js';
-/** Operand kinds a flat tensor-product can lower in v0.3.5: the three
- *  index-carrying nodes plus tensor-partial-derivative (whose effective
- *  indices are `of`'s indices followed by its wrtIndex). */
-type ContractableNode = Extract<ExprNode, {
-    indices: ReadonlyArray<TensorIndex>;
-}> | Extract<ExprNode, {
-    kind: 'tensor-partial-derivative';
-}>;
-/**
- * Build the EinsumSpec for a flat tensor-product.
- *
- * NOTE: this function does NOT decide which indices contract. That authority
- * belongs to `computeContraction()` (v0.2.0 symbolic-layer) — variance-aware,
- * implicit-identity-metric rule — which already classified every label as
- * contracted or free. buildEinsumSpec only maps those already-classified
- * labels to their (operand, axis) sites. There is exactly one
- * contraction-decision implementation in the codebase.
- *
- * @internal — cross-module/test use only; not part of the consumer surface.
- */
-export declare function buildEinsumSpec(operands: ReadonlyArray<ContractableNode>, contractionPairs: ReadonlyArray<{
-    label: string;
-}>, freeIndices: ReadonlyMap<string, {
-    upper: number;
-    lower: number;
-}>): EinsumSpec;
 /** Lower a validated ExprNode to an EngineTensor.
  *  @internal — cross-module/test use only; not part of the consumer surface. */
 export declare function lowerNode(node: ExprNode, inputs: NumericalInputs, engine: TensorEngine): EngineTensor;
-export {};
 //# sourceMappingURL=lowering.d.ts.map

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

@@ -1 +1 @@
-{"version":3,"file":"lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAE5D,OAAO,KAAK,EAAE,WAAW,EAAoB,MAAM,0BAA0B,CAAC;AAY9E,OAAO,KAAK,EACV,YAAY,EAAE,YAAY,EAAE,UAAU,EACvC,MAAM,oBAAoB,CAAC;AAC5B,OAAO,KAAK,EAAE,eAAe,EAAe,MAAM,YAAY,CAAC;AAkC/D;;4DAE4D;AAC5D,KAAK,gBAAgB,GACjB,OAAO,CAAC,QAAQ,EAAE;IAAE,OAAO,EAAE,aAAa,CAAC,WAAW,CAAC,CAAA;CAAE,CAAC,GAC1D,OAAO,CAAC,QAAQ,EAAE;IAAE,IAAI,EAAE,2BAA2B,CAAA;CAAE,CAAC,CAAC;AAqD7D;;;;;;;;;;;GAWG;AACH,wBAAgB,eAAe,CAC7B,QAAQ,EAAE,aAAa,CAAC,gBAAgB,CAAC,EACzC,gBAAgB,EAAE,aAAa,CAAC;IAAE,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC,EAClD,WAAW,EAAE,WAAW,CAAC,MAAM,EAAE;IAAE,KAAK,EAAE,MAAM,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,CAAC,GACjE,UAAU,CA0CZ;AAiBD;gFACgF;AAChF,wBAAgB,SAAS,CACvB,IAAI,EAAE,QAAQ,EACd,MAAM,EAAE,eAAe,EACvB,MAAM,EAAE,YAAY,GACnB,YAAY,CAulBd"}
+{"version":3,"file":"lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAiB5D,OAAO,KAAK,EACV,YAAY,EAAE,YAAY,EAC3B,MAAM,oBAAoB,CAAC;AAC5B,OAAO,KAAK,EAAE,eAAe,EAAe,MAAM,YAAY,CAAC;AAwW/D;gFACgF;AAChF,wBAAgB,SAAS,CACvB,IAAI,EAAE,QAAQ,EACd,MAAM,EAAE,eAAe,EACvB,MAAM,EAAE,YAAY,GACnB,YAAY,CAqMd"}