universal-physics-tensor 0.4.5

package/dist/numerical/lowering.js

@@ -0,0 +1,443 @@
+/**
+ * AST → EngineTensor lowering. Walks a validated ExprNode tree and emits
+ * TensorEngine calls. The tensor-product case turns computeContraction()'s
+ * contractionPairs into an EinsumSpec (v0.3.5-Design.md §5).
+ *
+ * v0.3.5 scope: flat tensor-products of contractable operands (the three
+ * index-carrying nodes + tensor-partial-derivative) + scalar operands.
+ * Nested products, integral/derivative numerical eval are out of scope and
+ * throw clearly.
+ *
+ * @module numerical/lowering
+ */
+import { validate } from '../dimensional/validator.js';
+import { computeContraction, validateTensorSymbol } from '../dimensional/tensor.js';
+import { pderivGrid, pderivNumericalFn, pderivSymbolic } from './pderiv.js';
+import { validateMetricTensor, validateKroneckerDelta, validatePartialDerivative, } from '../dimensional/metric-validators.js';
+import { NumericalBackendError } from './errors.js';
+import { zeroTensor, zeroTensorLike, flatToNested, flattenNA, tensorAdd, tensorAddScaled, computeChristoffelTensor, contractChristoffelWithOperand, getMetricDerivFlat, } from './connection-lowering-helpers.js';
+function isContractable(node) {
+    return node.kind === 'tensor-symbol'
+        || node.kind === 'metric-tensor'
+        || node.kind === 'kronecker-delta'
+        || node.kind === 'tensor-partial-derivative';
+}
+/** The effective index list of a contractable operand, in the SAME axis
+ *  order as its lowered EngineTensor. For the three index-carrying nodes
+ *  that is `node.indices`. For ∂_μ(of) the lowered tensor has shape
+ *  [...ofShape, N] (Task 10), so the effective indices are `of`'s indices
+ *  (v0.3.5: `of` is a tensor-symbol or metric-tensor) followed by the wrtIndex. */
+function operandIndices(node) {
+    if (node.kind === 'tensor-partial-derivative') {
+        const of = node.of;
+        if (of.kind !== 'tensor-symbol' && of.kind !== 'metric-tensor') {
+            throw new NumericalBackendError(`lowering: a tensor-partial-derivative operand requires a tensor-symbol or `
+                + `metric-tensor 'of' in v0.3.5/v0.4.0 — got '${of.kind}'`);
+        }
+        return [...of.indices, node.wrtIndex];
+    }
+    return node.indices;
+}
+function dimensionOf(inputs) {
+    return inputs.dimension ?? 4;
+}
+/** Look up a named tensor's concrete value, or throw. */
+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 expected size.
+ *  Delegates to the canonical flattenNA() from connection-lowering-helpers. */
+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;
+}
+/**
+ * Build the EinsumSpec for a flat tensor-product.
+ *
+ * CRITICAL (finding #1 of the v0.3.5 adversarial review): this function does
+ * NOT decide which indices contract. `computeContraction()` — the v0.2.0
+ * symbolic-layer authority, which is variance-aware and applies the
+ * implicit-identity-metric rule — 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 function buildEinsumSpec(operands, contractionPairs, freeIndices) {
+    // Map every label to its (operand, axis) sites — via operandIndices() so a
+    // tensor-partial-derivative operand contributes its [...of.indices, wrtIndex]
+    // effective axes.
+    const sites = new Map();
+    operands.forEach((op, opIndex) => {
+        operandIndices(op).forEach((idx, axis) => {
+            const list = sites.get(idx.label) ?? [];
+            list.push([opIndex, axis]);
+            sites.set(idx.label, list);
+        });
+    });
+    // Contractions: computeContraction told us exactly which labels contract.
+    const contractions = [];
+    for (const { label } of contractionPairs) {
+        const locs = sites.get(label);
+        if (!locs || locs.length !== 2) {
+            throw new NumericalBackendError(`buildEinsumSpec: contracted label "${label}" must occur at exactly 2 `
+                + `operand sites, found ${locs?.length ?? 0}`);
+        }
+        contractions.push({ pair: [locs[0], locs[1]] });
+    }
+    // Free axes in computeContraction's freeIndices order — the symbolic
+    // layer's canonical ordering, so the numerical result's axis order
+    // matches what the rest of UPT expects.
+    const free = [];
+    for (const label of freeIndices.keys()) {
+        const locs = sites.get(label);
+        if (!locs || locs.length !== 1) {
+            throw new NumericalBackendError(`buildEinsumSpec: free label "${label}" must occur at exactly 1 `
+                + `operand site, found ${locs?.length ?? 0}`);
+        }
+        free.push({ operand: locs[0][0], axis: locs[0][1] });
+    }
+    return { contractions, free };
+}
+/** Lower a contractable operand to a concrete EngineTensor. */
+function lowerContractable(node, inputs, engine) {
+    const N = dimensionOf(inputs);
+    if (node.kind === 'kronecker-delta')
+        return engine.identity(N);
+    if (node.kind === 'tensor-partial-derivative') {
+        // Delegate to lowerNode, which handles all three numericalForm paths.
+        return lowerNode(node, inputs, engine);
+    }
+    // tensor-symbol / metric-tensor: shape is N per index.
+    const shape = node.indices.map(() => N);
+    return engine.fromNested(requireValue(node.name, inputs), shape);
+}
+/** Lower a validated ExprNode to an EngineTensor.
+ *  @internal — cross-module/test use only; not part of the consumer surface. */
+export function lowerNode(node, inputs, engine) {
+    switch (node.kind) {
+        case 'symbol':
+            return engine.fromNested(requireValue(node.name, inputs), []);
+        case 'tensor-symbol':
+        case 'metric-tensor':
+        case 'kronecker-delta':
+            return lowerContractable(node, inputs, engine);
+        case 'op': {
+            if (node.op === '+' || node.op === '-') {
+                if (node.args.length === 0)
+                    return engine.fromNested(0, []);
+                let acc = lowerNode(node.args[0], inputs, engine);
+                for (let i = 1; i < node.args.length; i++) {
+                    const next = lowerNode(node.args[i], inputs, engine);
+                    acc = node.op === '+' ? engine.add(acc, next) : engine.sub(acc, next);
+                }
+                return acc;
+            }
+            // '*' / '/' / '^' are scalar-only. Guard arity so an unvalidated AST
+            // surfaces a clean NumericalBackendError instead of a raw TypeError
+            // (zero-operand '/') or a silent NaN (wrong-arity '^'). The '*' case
+            // with zero args is already fine — reduce(..., 1) returns 1.
+            if (node.op === '/' && node.args.length === 0) {
+                throw new NumericalBackendError("lowering: op '/' requires at least one operand");
+            }
+            if (node.op === '^' && node.args.length !== 2) {
+                throw new NumericalBackendError(`lowering: op '^' requires exactly 2 operands (base, exponent), got ${node.args.length}`);
+            }
+            // '*' / '/' / '^' are scalar-only (the validator rejects tensor
+            // operands). Lower each to rank-0, do the arithmetic in JS, lift back.
+            const scalars = node.args.map((a) => {
+                const t = lowerNode(a, inputs, engine);
+                if (t.shape.length !== 0) {
+                    throw new NumericalBackendError(`lowering: op '${node.op}' got a rank-${t.shape.length} operand — scalar ops require rank-0`);
+                }
+                return engine.toNested(t);
+            });
+            let value;
+            if (node.op === '*')
+                value = scalars.reduce((a, b) => a * b, 1);
+            else if (node.op === '/')
+                value = scalars.reduce((a, b) => a / b);
+            else
+                value = Math.pow(scalars[0], scalars[1]); // '^'
+            return engine.fromNested(value, []);
+        }
+        case 'tensor-product': {
+            for (const arg of node.args) {
+                if (arg.kind === 'tensor-product') {
+                    throw new NumericalBackendError('lowering: nested tensor-product numerical evaluation is not supported in v0.3.5 — '
+                        + 'flatten the contraction into a single product');
+                }
+            }
+            // Contractable operands (tensor-symbol / metric-tensor / kronecker-delta
+            // / tensor-partial-derivative) participate in the einsum; everything
+            // else must be a rank-0 scalar factor.
+            const operands = node.args.filter(isContractable);
+            const scalarArgs = node.args.filter((a) => !isContractable(a));
+            // computeContraction is the single authority on WHICH indices contract
+            // (variance-aware, implicit-metric rule — finding #1). buildEinsumSpec
+            // only maps the labels it classified to (operand, axis) sites. The
+            // recursive validateContractionChild resolves tensor-partial-derivative
+            // operands via validatePartialDerivative.
+            function validateContractionChild(child) {
+                if (child.kind === 'tensor-symbol')
+                    return validateTensorSymbol(child);
+                if (child.kind === 'metric-tensor')
+                    return validateMetricTensor(child);
+                if (child.kind === 'kronecker-delta')
+                    return validateKroneckerDelta(child);
+                if (child.kind === 'tensor-partial-derivative') {
+                    const r = validatePartialDerivative(child, (g) => validateContractionChild(g));
+                    return { dim: r.dim, freeIndices: r.freeIndices };
+                }
+                throw new NumericalBackendError(`lowering: unexpected operand '${child.kind}' in tensor-product einsum`);
+            }
+            const { contractionPairs, freeIndices } = computeContraction(operands, validateContractionChild);
+            const spec = buildEinsumSpec(operands, contractionPairs, freeIndices);
+            const operandTensors = operands.map((n) => lowerContractable(n, inputs, engine));
+            let result = engine.einsum(spec, ...operandTensors);
+            // Scalar operands multiply the whole contraction.
+            for (const s of scalarArgs) {
+                const st = lowerNode(s, inputs, engine);
+                if (st.shape.length !== 0) {
+                    throw new NumericalBackendError('lowering: non-scalar non-contractable operand in tensor-product');
+                }
+                result = engine.scale(result, engine.toNested(st));
+            }
+            return result;
+        }
+        case 'tensor-partial-derivative': {
+            // v0.3.5/v0.4.0 scope: `of` is a tensor-symbol or metric-tensor.
+            // ∂_μ(of) adds the wrtIndex as a trailing axis — the result shape is
+            // [...ofShape, N], NOT ofShape. (For BE-37, `of` = the scalar S is
+            // rank-0, so ∂_μ S is the rank-1 wave covector k_μ, shape [N].)
+            const of = node.of;
+            // v0.4.0 extension: metric-tensor pderiv dispatch.
+            if (of.kind === 'metric-tensor') {
+                const mNode = of;
+                const strategy = mNode.derivativeStrategy ?? 'computed';
+                const N = dimensionOf(inputs);
+                const coordLabel = node.wrtIndex.label;
+                const ofShape = mNode.indices.map(() => N);
+                const resultShape = [...ofShape, N];
+                if (strategy === 'zero') {
+                    // ∂g = 0 everywhere (constant/flat metric).
+                    return zeroTensor(resultShape, engine);
+                }
+                if (strategy === 'supplied') {
+                    // Look up the N slices ∂_mu g for mu=0..N-1 and stack them as the
+                    // trailing axis. Key format: `${metricName}/${coordLabel}_${mu}`.
+                    // Result shape: [...ofShape, N] = [N, N, N] for a rank-2 metric.
+                    // Build as a flat array then convert to nested for engine.fromNested().
+                    const size = resultShape.reduce((a, b) => a * b, 1);
+                    const flat = new Array(size).fill(0);
+                    // ofShape = [N, N], resultShape = [N, N, N]
+                    // flat[i*N*N + j*N + mu] = (∂_mu g)[i][j]
+                    for (let mu = 0; mu < N; mu++) {
+                        const key = `${mNode.name}/${coordLabel}_${mu}`;
+                        const slice = inputs.metricDerivatives?.get(key);
+                        if (slice === undefined) {
+                            throw new NumericalBackendError(`lowering: metric-tensor pderiv with strategy='supplied': ` +
+                                `no metricDerivatives entry for "${key}"`);
+                        }
+                        // Flatten the slice (shape [N,N]) and write into flat at stride N (last axis)
+                        const flatSlice = flattenNestedArray(slice, N * N);
+                        for (let ij = 0; ij < N * N; ij++) {
+                            flat[ij * N + mu] = flatSlice[ij];
+                        }
+                    }
+                    return engine.fromNested(flatToNested(flat, resultShape), resultShape);
+                }
+                // strategy === 'computed': treat as zero for constant-tensor metrics
+                // (the raw metric supplied via inputs.tensors has no coordinate dependence).
+                return zeroTensor(resultShape, engine);
+            }
+            if (of.kind !== 'tensor-symbol') {
+                throw new NumericalBackendError(`lowering: tensor-partial-derivative numerical eval requires a tensor-symbol `
+                    + `or metric-tensor 'of' operand in v0.3.5/v0.4.0 — got '${of.kind}'`);
+            }
+            const sym = of;
+            const form = sym.numericalForm ?? 'symbolic';
+            const coordLabel = node.wrtIndex.label;
+            const N = dimensionOf(inputs);
+            const ofShape = sym.indices.map(() => N);
+            const resultShape = [...ofShape, N];
+            if (form === 'symbolic') {
+                // pderivSymbolic returns the caller-supplied full ∂_μ(of) tensor,
+                // which must already be of shape [...ofShape, N].
+                const d = pderivSymbolic(sym.name, coordLabel, inputs.derivatives ?? new Map());
+                return engine.fromNested(d, resultShape);
+            }
+            if (form === 'numerical-fn') {
+                // v0.3.5: 'numerical-fn' lowering is scoped to a rank-0 `of` (scalar
+                // field). ∂_μ ranges over all N coordinate axes — stack the N
+                // single-axis derivatives into the rank-1 result. Higher-rank fields
+                // under 'numerical-fn' are a v0.4.0 concern.
+                if (ofShape.length !== 0) {
+                    throw new NumericalBackendError(`lowering: 'numerical-fn' pderiv lowering supports a rank-0 'of' in v0.3.5; `
+                        + `"${sym.name}" is rank ${ofShape.length}`);
+                }
+                const fn = inputs.fields?.get(sym.name);
+                if (!fn) {
+                    throw new NumericalBackendError(`lowering: 'numerical-fn' tensor-symbol "${sym.name}" has no field fn in inputs.fields`);
+                }
+                const coordValues = inputs.coords ? [...inputs.coords.values()] : [];
+                const components = [];
+                for (let axis = 0; axis < N; axis++) {
+                    components.push(pderivNumericalFn(fn, coordValues, axis));
+                }
+                return engine.fromNested(components, [N]);
+            }
+            // form === 'grid': the GridField is the field sampled over space, and
+            // pderivGrid returns the derivative field sampled on that same grid —
+            // result shape is grid.shape. This is a distinct semantic from
+            // symbolic/numerical-fn (a sampled derivative field, not a single
+            // tensor), kept as the v0.5.0 BSSN forward-compat path. v0.3.5 has no
+            // release test driving 'grid' through lowering; pderivGrid itself is
+            // unit-tested in Task 10's pderiv.test.ts.
+            const grid = inputs.grids?.get(sym.name);
+            if (!grid) {
+                throw new NumericalBackendError(`lowering: 'grid' tensor-symbol "${sym.name}" has no GridField in inputs.grids`);
+            }
+            const gridAxis = inputs.coords ? [...inputs.coords.keys()].indexOf(coordLabel) : 0;
+            const flat = pderivGrid(grid, gridAxis < 0 ? 0 : gridAxis);
+            return engine.fromNested(flat.length === 1 ? flat[0] : flat, grid.shape);
+        }
+        case 'covariant-derivative': {
+            // S2(a) fix: `of.freeIndices` does NOT exist on the raw ExprNode.
+            // Re-validate the `of` subtree to obtain its free-index structure.
+            const covNode = node;
+            const ofExpr = covNode.of;
+            const ofValidation = validate(ofExpr);
+            // Build ordered list of free indices: [{label, variance}].
+            // Iterate of.indices (NOT validation.freeIndices Map) — declaration order IS
+            // the axis layout. For tensor-symbol and metric-tensor, of.indices is the
+            // canonical axis order; the freeIndices Map insertion order is also
+            // of.indices order (validateTensorSymbol iterates node.indices), but relying
+            // on that is an undocumented invariant. Iterating of.indices directly makes
+            // the axis ordering guarantee explicit and safe for future of-kinds.
+            const ofFreeIndices = [];
+            const ofIndices = ofExpr.indices;
+            if (ofIndices) {
+                // tensor-symbol / metric-tensor: iterate the declared indices in order.
+                // Each index is either free (present in ofValidation.freeIndices) or
+                // contracted (absent — skip). In practice, of.indices for a simple
+                // tensor-symbol or metric-tensor has no contracted indices, but we
+                // guard with the Map lookup for safety.
+                for (const idx of ofIndices) {
+                    const counts = ofValidation.freeIndices.get(idx.label);
+                    if (counts === undefined)
+                        continue; // contracted — not a free axis
+                    ofFreeIndices.push({
+                        label: idx.label,
+                        variance: idx.variance,
+                        pos: ofFreeIndices.length,
+                    });
+                }
+            }
+            else {
+                // Fallback for future of-kinds without .indices (e.g. tensor-product).
+                // Map iteration order is insertion order — a best-effort axis ordering.
+                let axisPos = 0;
+                for (const [label, counts] of ofValidation.freeIndices) {
+                    for (let i = 0; i < counts.upper; i++) {
+                        ofFreeIndices.push({ label, variance: 'upper', pos: axisPos++ });
+                    }
+                    for (let i = 0; i < counts.lower; i++) {
+                        ofFreeIndices.push({ label, variance: 'lower', pos: axisPos++ });
+                    }
+                }
+            }
+            // Lower the operand tensor.
+            const ofTensor = lowerNode(ofExpr, inputs, engine);
+            const N = dimensionOf(inputs);
+            const strategy = covNode.gLower.derivativeStrategy ?? 'computed';
+            // S2(b): strategy='zero' → flat space, Γ=0, ∇_μ T = ∂_μ T.
+            // For constant tensors (like a flat metric), ∂_μ T = 0, so result is all zeros.
+            // We return a zero tensor of shape [...ofShape, N] (wrt axis appended last).
+            if (strategy === 'zero') {
+                const outShape = [...ofTensor.shape, N];
+                return zeroTensor(outShape, engine);
+            }
+            // v0.4.0 CRITICAL FIX (Finding #1): 'computed' on a raw-tensor metric means
+            // constant metric → Γ = 0 → covariant-derivative = partial derivative only.
+            // v0.5.0 will replace this with coordinate-grid finite-difference.
+            // This early return MUST come before the Christoffel construction below,
+            // because getMetricDerivFlat only accepts 'zero' | 'supplied'; if 'computed'
+            // fell through, the type cast was silently wrong at runtime and the function
+            // threw NumericalBackendError for every coordinate.
+            //
+            // Compute partial first so we can return it directly.
+            // (partial is also needed by the 'supplied' path below, so we compute it
+            //  unconditionally and use the early return only for 'computed'.)
+            let partial;
+            if (ofExpr.kind === 'tensor-symbol' || ofExpr.kind === 'metric-tensor') {
+                const pdNode = {
+                    kind: 'tensor-partial-derivative',
+                    of: ofExpr,
+                    wrt: covNode.wrt,
+                    wrtIndex: covNode.wrtIndex,
+                };
+                partial = lowerNode(pdNode, inputs, engine);
+            }
+            else {
+                // Scalar or other: partial derivative is zero.
+                const outShape = [...ofTensor.shape, N];
+                partial = zeroTensor(outShape, engine);
+            }
+            // v0.4.0 spec: 'computed' on a raw-tensor metric = constant metric → Γ = 0.
+            // The covariant-derivative reduces to the ordinary partial derivative.
+            // v0.5.0 will add coordinate-grid finite-difference here.
+            if (strategy === 'computed') {
+                return partial;
+            }
+            // S2(c) + S2(d): Build Christoffel Γ^α_{μν} from metric data and apply
+            // the sign rule to all free indices of `of`:
+            //   ∇_μ T^α_β = ∂_μ T^α_β + Γ^α_{μλ} T^λ_β − Γ^λ_{μβ} T^α_λ
+            //
+            // If `of` has no free indices (scalar), no correction needed.
+            if (ofFreeIndices.length === 0) {
+                return partial;
+            }
+            // Get metric and inverse metric data.
+            const gLowerNode = covNode.gLower;
+            const gInverseNode = covNode.gInverse;
+            const gInverseData = flattenNestedArray(requireValue(gInverseNode.name, inputs), N * N);
+            // coordLabel for metricDerivatives keys (wrtIndex.label of the covariant derivative)
+            const coordLabel = covNode.wrtIndex.label;
+            // getMetricDeriv(mu): returns flat [N*N] of ∂_{mu} g
+            const getMetricDeriv = (mu) => getMetricDerivFlat(gLowerNode.name, coordLabel, mu, strategy, N, inputs.metricDerivatives);
+            // Compute Γ^α_{μν} from metric data.
+            const GammaTensor = computeChristoffelTensor(gInverseData, getMetricDeriv, N, engine);
+            const GammaFlat = flattenNestedArray(engine.toNested(GammaTensor), N * N * N);
+            // Apply the Christoffel correction for each free index of `of`.
+            // S2(d) sign rule: upper → +Γ, lower → −Γ.
+            const ofFlat = flattenNestedArray(engine.toNested(ofTensor), ofTensor.shape.reduce((a, b) => a * b, 1));
+            const ofShapeArr = [...ofTensor.shape];
+            let correction = zeroTensorLike(partial, engine);
+            for (const freeIdx of ofFreeIndices) {
+                const term = contractChristoffelWithOperand(GammaFlat, ofFlat, ofShapeArr, freeIdx.pos, freeIdx.variance, N, engine);
+                const sign = freeIdx.variance === 'upper' ? 1 : -1;
+                correction = tensorAddScaled(correction, term, sign, engine);
+            }
+            return tensorAdd(partial, correction, engine);
+        }
+        case 'integral':
+        case 'derivative':
+            throw new NumericalBackendError(`lowering: '${node.kind}' is not numerically evaluated in v0.3.5 — `
+                + 'use tensor-partial-derivative for differentiation');
+        default: {
+            const _exhaustive = node;
+            void _exhaustive;
+            throw new NumericalBackendError(`lowering: unknown ExprNode.kind ${JSON.stringify(node.kind)}`);
+        }
+    }
+}
+//# sourceMappingURL=lowering.js.map

package/dist/numerical/lowering.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"lowering.js","sourceRoot":"","sources":["../../src/numerical/lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAGH,OAAO,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAGvD,OAAO,EAAE,kBAAkB,EAAE,oBAAoB,EAAE,MAAM,0BAA0B,CAAC;AACpF,OAAO,EAAE,UAAU,EAAE,iBAAiB,EAAE,cAAc,EAAE,MAAM,aAAa,CAAC;AAC5E,OAAO,EACL,oBAAoB,EACpB,sBAAsB,EACtB,yBAAyB,GAC1B,MAAM,qCAAqC,CAAC;AAO7C,OAAO,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AACpD,OAAO,EACL,UAAU,EACV,cAAc,EACd,YAAY,EACZ,SAAS,EACT,SAAS,EACT,eAAe,EACf,wBAAwB,EACxB,8BAA8B,EAC9B,kBAAkB,GACnB,MAAM,kCAAkC,CAAC;AAS1C,SAAS,cAAc,CAAC,IAAc;IACpC,OAAO,IAAI,CAAC,IAAI,KAAK,eAAe;WAC/B,IAAI,CAAC,IAAI,KAAK,eAAe;WAC7B,IAAI,CAAC,IAAI,KAAK,iBAAiB;WAC/B,IAAI,CAAC,IAAI,KAAK,2BAA2B,CAAC;AACjD,CAAC;AAED;;;;mFAImF;AACnF,SAAS,cAAc,CAAC,IAAsB;IAC5C,IAAI,IAAI,CAAC,IAAI,KAAK,2BAA2B,EAAE,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,EAAc,CAAC;QAC/B,IAAI,EAAE,CAAC,IAAI,KAAK,eAAe,IAAI,EAAE,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;YAC/D,MAAM,IAAI,qBAAqB,CAC7B,4EAA4E;kBAC1E,8CAA8C,EAAE,CAAC,IAAI,GAAG,CAC3D,CAAC;QACJ,CAAC;QACD,OAAO,CAAC,GAAG,EAAE,CAAC,OAAO,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;IACxC,CAAC;IACD,OAAO,IAAI,CAAC,OAAO,CAAC;AACtB,CAAC;AAED,SAAS,WAAW,CAAC,MAAuB;IAC1C,OAAO,MAAM,CAAC,SAAS,IAAI,CAAC,CAAC;AAC/B,CAAC;AAED,yDAAyD;AACzD,SAAS,YAAY,CAAC,IAAY,EAAE,MAAuB;IACzD,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;+EAC+E;AAC/E,SAAS,kBAAkB,CAAC,IAAiB,EAAE,YAAoB;IACjE,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;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,eAAe,CAC7B,QAAyC,EACzC,gBAAkD,EAClD,WAAkE;IAElE,2EAA2E;IAC3E,8EAA8E;IAC9E,kBAAkB;IAClB,MAAM,KAAK,GAAG,IAAI,GAAG,EAAmC,CAAC;IACzD,QAAQ,CAAC,OAAO,CAAC,CAAC,EAAE,EAAE,OAAO,EAAE,EAAE;QAC/B,cAAc,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,GAAG,EAAE,IAAI,EAAE,EAAE;YACvC,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,EAAE,CAAC;YACxC,IAAI,CAAC,IAAI,CAAC,CAAC,OAAO,EAAE,IAAI,CAAC,CAAC,CAAC;YAC3B,KAAK,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;QAC7B,CAAC,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IAEH,0EAA0E;IAC1E,MAAM,YAAY,GAAwB,EAAE,CAAC;IAC7C,KAAK,MAAM,EAAE,KAAK,EAAE,IAAI,gBAAgB,EAAE,CAAC;QACzC,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC9B,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC/B,MAAM,IAAI,qBAAqB,CAC7B,sCAAsC,KAAK,4BAA4B;kBACrE,wBAAwB,IAAI,EAAE,MAAM,IAAI,CAAC,EAAE,CAC9C,CAAC;QACJ,CAAC;QACD,YAAY,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IAClD,CAAC;IAED,qEAAqE;IACrE,mEAAmE;IACnE,wCAAwC;IACxC,MAAM,IAAI,GAA6C,EAAE,CAAC;IAC1D,KAAK,MAAM,KAAK,IAAI,WAAW,CAAC,IAAI,EAAE,EAAE,CAAC;QACvC,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC9B,IAAI,CAAC,IAAI,IAAI,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;YAC/B,MAAM,IAAI,qBAAqB,CAC7B,gCAAgC,KAAK,4BAA4B;kBAC/D,uBAAuB,IAAI,EAAE,MAAM,IAAI,CAAC,EAAE,CAC7C,CAAC;QACJ,CAAC;QACD,IAAI,CAAC,IAAI,CAAC,EAAE,OAAO,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;IACvD,CAAC;IAED,OAAO,EAAE,YAAY,EAAE,IAAI,EAAE,CAAC;AAChC,CAAC;AAED,+DAA+D;AAC/D,SAAS,iBAAiB,CACxB,IAAsB,EAAE,MAAuB,EAAE,MAAoB;IAErE,MAAM,CAAC,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;IAC9B,IAAI,IAAI,CAAC,IAAI,KAAK,iBAAiB;QAAE,OAAO,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IAC/D,IAAI,IAAI,CAAC,IAAI,KAAK,2BAA2B,EAAE,CAAC;QAC9C,sEAAsE;QACtE,OAAO,SAAS,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACzC,CAAC;IACD,uDAAuD;IACvD,MAAM,KAAK,GAAG,IAAI,CAAC,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;IACxC,OAAO,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,IAAI,CAAC,IAAI,EAAE,MAAM,CAAC,EAAE,KAAK,CAAC,CAAC;AACnE,CAAC;AAED;gFACgF;AAChF,MAAM,UAAU,SAAS,CACvB,IAAc,EACd,MAAuB,EACvB,MAAoB;IAEpB,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,QAAQ;YACX,OAAO,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,IAAI,CAAC,IAAI,EAAE,MAAM,CAAC,EAAE,EAAE,CAAC,CAAC;QAEhE,KAAK,eAAe,CAAC;QACrB,KAAK,eAAe,CAAC;QACrB,KAAK,iBAAiB;YACpB,OAAO,iBAAiB,CAAC,IAAI,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;QAEjD,KAAK,IAAI,CAAC,CAAC,CAAC;YACV,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,EAAE,CAAC;gBACvC,IAAI,IAAI,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC;oBAAE,OAAO,MAAM,CAAC,UAAU,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;gBAC5D,IAAI,GAAG,GAAG,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;gBAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC1C,MAAM,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;oBACrD,GAAG,GAAG,IAAI,CAAC,EAAE,KAAK,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;gBACxE,CAAC;gBACD,OAAO,GAAG,CAAC;YACb,CAAC;YACD,qEAAqE;YACrE,oEAAoE;YACpE,qEAAqE;YACrE,6DAA6D;YAC7D,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAC9C,MAAM,IAAI,qBAAqB,CAAC,gDAAgD,CAAC,CAAC;YACpF,CAAC;YACD,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAC9C,MAAM,IAAI,qBAAqB,CAC7B,sEAAsE,IAAI,CAAC,IAAI,CAAC,MAAM,EAAE,CACzF,CAAC;YACJ,CAAC;YACD,gEAAgE;YAChE,uEAAuE;YACvE,MAAM,OAAO,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;gBAClC,MAAM,CAAC,GAAG,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;gBACvC,IAAI,CAAC,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBACzB,MAAM,IAAI,qBAAqB,CAC7B,iBAAiB,IAAI,CAAC,EAAE,gBAAgB,CAAC,CAAC,KAAK,CAAC,MAAM,sCAAsC,CAC7F,CAAC;gBACJ,CAAC;gBACD,OAAO,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAW,CAAC;YACtC,CAAC,CAAC,CAAC;YACH,IAAI,KAAa,CAAC;YAClB,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG;gBAAE,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;iBAC3D,IAAI,IAAI,CAAC,EAAE,KAAK,GAAG;gBAAE,KAAK,GAAG,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;;gBAC7D,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM;YACrD,OAAO,MAAM,CAAC,UAAU,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACtC,CAAC;QAED,KAAK,gBAAgB,CAAC,CAAC,CAAC;YACtB,KAAK,MAAM,GAAG,IAAI,IAAI,CAAC,IAAI,EAAE,CAAC;gBAC5B,IAAI,GAAG,CAAC,IAAI,KAAK,gBAAgB,EAAE,CAAC;oBAClC,MAAM,IAAI,qBAAqB,CAC7B,oFAAoF;0BAClF,+CAA+C,CAClD,CAAC;gBACJ,CAAC;YACH,CAAC;YACD,yEAAyE;YACzE,qEAAqE;YACrE,uCAAuC;YACvC,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,cAAc,CAAC,CAAC;YAClD,MAAM,UAAU,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/D,uEAAuE;YACvE,uEAAuE;YACvE,mEAAmE;YACnE,wEAAwE;YACxE,0CAA0C;YAC1C,SAAS,wBAAwB,CAAC,KAAe;gBAI/C,IAAI,KAAK,CAAC,IAAI,KAAK,eAAe;oBAAE,OAAO,oBAAoB,CAAC,KAAK,CAAC,CAAC;gBACvE,IAAI,KAAK,CAAC,IAAI,KAAK,eAAe;oBAAE,OAAO,oBAAoB,CAAC,KAAK,CAAC,CAAC;gBACvE,IAAI,KAAK,CAAC,IAAI,KAAK,iBAAiB;oBAAE,OAAO,sBAAsB,CAAC,KAAK,CAAC,CAAC;gBAC3E,IAAI,KAAK,CAAC,IAAI,KAAK,2BAA2B,EAAE,CAAC;oBAC/C,MAAM,CAAC,GAAG,yBAAyB,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAC/C,wBAAwB,CAAC,CAAa,CAAC,CACxC,CAAC;oBACF,OAAO,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,WAAW,EAAE,CAAC,CAAC,WAAW,EAAE,CAAC;gBACpD,CAAC;gBACD,MAAM,IAAI,qBAAqB,CAC7B,iCAAiC,KAAK,CAAC,IAAI,4BAA4B,CACxE,CAAC;YACJ,CAAC;YACD,MAAM,EAAE,gBAAgB,EAAE,WAAW,EAAE,GAAG,kBAAkB,CAC1D,QAAQ,EAAE,wBAAwB,CACnC,CAAC;YACF,MAAM,IAAI,GAAG,eAAe,CAAC,QAAQ,EAAE,gBAAgB,EAAE,WAAW,CAAC,CAAC;YACtE,MAAM,cAAc,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,CAAC;YACjF,IAAI,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,cAAc,CAAC,CAAC;YACpD,kDAAkD;YAClD,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;gBACxC,IAAI,EAAE,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBAC1B,MAAM,IAAI,qBAAqB,CAC7B,iEAAiE,CAClE,CAAC;gBACJ,CAAC;gBACD,MAAM,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAW,CAAC,CAAC;YAC/D,CAAC;YACD,OAAO,MAAM,CAAC;QAChB,CAAC;QAED,KAAK,2BAA2B,CAAC,CAAC,CAAC;YACjC,iEAAiE;YACjE,qEAAqE;YACrE,mEAAmE;YACnE,gEAAgE;YAChE,MAAM,EAAE,GAAG,IAAI,CAAC,EAAc,CAAC;YAE/B,mDAAmD;YACnD,IAAI,EAAE,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;gBAChC,MAAM,KAAK,GAAG,EAAsB,CAAC;gBACrC,MAAM,QAAQ,GAAG,KAAK,CAAC,kBAAkB,IAAI,UAAU,CAAC;gBACxD,MAAM,CAAC,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;gBAC9B,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC;gBACvC,MAAM,OAAO,GAAG,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAC3C,MAAM,WAAW,GAAG,CAAC,GAAG,OAAO,EAAE,CAAC,CAAC,CAAC;gBAEpC,IAAI,QAAQ,KAAK,MAAM,EAAE,CAAC;oBACxB,4CAA4C;oBAC5C,OAAO,UAAU,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;gBACzC,CAAC;gBACD,IAAI,QAAQ,KAAK,UAAU,EAAE,CAAC;oBAC5B,kEAAkE;oBAClE,kEAAkE;oBAClE,iEAAiE;oBACjE,wEAAwE;oBACxE,MAAM,IAAI,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;oBACpD,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC7C,4CAA4C;oBAC5C,0CAA0C;oBAC1C,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;wBAC9B,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,IAAI,IAAI,UAAU,IAAI,EAAE,EAAE,CAAC;wBAChD,MAAM,KAAK,GAAG,MAAM,CAAC,iBAAiB,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;wBACjD,IAAI,KAAK,KAAK,SAAS,EAAE,CAAC;4BACxB,MAAM,IAAI,qBAAqB,CAC7B,2DAA2D;gCAC3D,mCAAmC,GAAG,GAAG,CAC1C,CAAC;wBACJ,CAAC;wBACD,8EAA8E;wBAC9E,MAAM,SAAS,GAAG,kBAAkB,CAAC,KAAK,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;wBACnD,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;4BAClC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,GAAG,SAAS,CAAC,EAAE,CAAC,CAAC;wBACpC,CAAC;oBACH,CAAC;oBACD,OAAO,MAAM,CAAC,UAAU,CAAC,YAAY,CAAC,IAAI,EAAE,WAAW,CAAC,EAAE,WAAW,CAAC,CAAC;gBACzE,CAAC;gBAED,qEAAqE;gBACrE,6EAA6E;gBAC7E,OAAO,UAAU,CAAC,WAAW,EAAE,MAAM,CAAC,CAAC;YACzC,CAAC;YAED,IAAI,EAAE,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;gBAChC,MAAM,IAAI,qBAAqB,CAC7B,8EAA8E;sBAC5E,yDAAyD,EAAE,CAAC,IAAI,GAAG,CACtE,CAAC;YACJ,CAAC;YACD,MAAM,GAAG,GAAG,EAAsB,CAAC;YACnC,MAAM,IAAI,GAAG,GAAG,CAAC,aAAa,IAAI,UAAU,CAAC;YAC7C,MAAM,UAAU,GAAG,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC;YACvC,MAAM,CAAC,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;YAC9B,MAAM,OAAO,GAAG,GAAG,CAAC,OAAO,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YACzC,MAAM,WAAW,GAAG,CAAC,GAAG,OAAO,EAAE,CAAC,CAAC,CAAC;YAEpC,IAAI,IAAI,KAAK,UAAU,EAAE,CAAC;gBACxB,kEAAkE;gBAClE,kDAAkD;gBAClD,MAAM,CAAC,GAAG,cAAc,CAAC,GAAG,CAAC,IAAI,EAAE,UAAU,EAAE,MAAM,CAAC,WAAW,IAAI,IAAI,GAAG,EAAE,CAAC,CAAC;gBAChF,OAAO,MAAM,CAAC,UAAU,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC;YAC3C,CAAC;YAED,IAAI,IAAI,KAAK,cAAc,EAAE,CAAC;gBAC5B,qEAAqE;gBACrE,8DAA8D;gBAC9D,qEAAqE;gBACrE,6CAA6C;gBAC7C,IAAI,OAAO,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;oBACzB,MAAM,IAAI,qBAAqB,CAC7B,6EAA6E;0BAC3E,IAAI,GAAG,CAAC,IAAI,aAAa,OAAO,CAAC,MAAM,EAAE,CAC5C,CAAC;gBACJ,CAAC;gBACD,MAAM,EAAE,GAAG,MAAM,CAAC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;gBACxC,IAAI,CAAC,EAAE,EAAE,CAAC;oBACR,MAAM,IAAI,qBAAqB,CAC7B,2CAA2C,GAAG,CAAC,IAAI,oCAAoC,CACxF,CAAC;gBACJ,CAAC;gBACD,MAAM,WAAW,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACrE,MAAM,UAAU,GAAa,EAAE,CAAC;gBAChC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;oBACpC,UAAU,CAAC,IAAI,CAAC,iBAAiB,CAAC,EAAE,EAAE,WAAW,EAAE,IAAI,CAAW,CAAC,CAAC;gBACtE,CAAC;gBACD,OAAO,MAAM,CAAC,UAAU,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5C,CAAC;YAED,sEAAsE;YACtE,sEAAsE;YACtE,+DAA+D;YAC/D,kEAAkE;YAClE,sEAAsE;YACtE,qEAAqE;YACrE,2CAA2C;YAC3C,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;YACzC,IAAI,CAAC,IAAI,EAAE,CAAC;gBACV,MAAM,IAAI,qBAAqB,CAC7B,mCAAmC,GAAG,CAAC,IAAI,oCAAoC,CAChF,CAAC;YACJ,CAAC;YACD,MAAM,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnF,MAAM,IAAI,GAAG,UAAU,CAAC,IAAI,EAAE,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;YAC3D,OAAO,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC;QAC3E,CAAC;QAED,KAAK,sBAAsB,CAAC,CAAC,CAAC;YAC5B,kEAAkE;YAClE,mEAAmE;YACnE,MAAM,OAAO,GAAG,IAA+B,CAAC;YAChD,MAAM,MAAM,GAAG,OAAO,CAAC,EAAc,CAAC;YACtC,MAAM,YAAY,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC;YACtC,2DAA2D;YAC3D,6EAA6E;YAC7E,0EAA0E;YAC1E,oEAAoE;YACpE,6EAA6E;YAC7E,4EAA4E;YAC5E,qEAAqE;YACrE,MAAM,aAAa,GAAuE,EAAE,CAAC;YAC7F,MAAM,SAAS,GAAI,MAA2E,CAAC,OAAO,CAAC;YACvG,IAAI,SAAS,EAAE,CAAC;gBACd,wEAAwE;gBACxE,qEAAqE;gBACrE,mEAAmE;gBACnE,mEAAmE;gBACnE,wCAAwC;gBACxC,KAAK,MAAM,GAAG,IAAI,SAAS,EAAE,CAAC;oBAC5B,MAAM,MAAM,GAAG,YAAY,CAAC,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;oBACvD,IAAI,MAAM,KAAK,SAAS;wBAAE,SAAS,CAAC,+BAA+B;oBACnE,aAAa,CAAC,IAAI,CAAC;wBACjB,KAAK,EAAE,GAAG,CAAC,KAAK;wBAChB,QAAQ,EAAE,GAAG,CAAC,QAA6B;wBAC3C,GAAG,EAAE,aAAa,CAAC,MAAM;qBAC1B,CAAC,CAAC;gBACL,CAAC;YACH,CAAC;iBAAM,CAAC;gBACN,uEAAuE;gBACvE,wEAAwE;gBACxE,IAAI,OAAO,GAAG,CAAC,CAAC;gBAChB,KAAK,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,IAAI,YAAY,CAAC,WAAW,EAAE,CAAC;oBACvD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;wBACtC,aAAa,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,EAAE,CAAC,CAAC;oBACnE,CAAC;oBACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;wBACtC,aAAa,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE,EAAE,CAAC,CAAC;oBACnE,CAAC;gBACH,CAAC;YACH,CAAC;YAED,4BAA4B;YAC5B,MAAM,QAAQ,GAAG,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;YACnD,MAAM,CAAC,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC;YAE9B,MAAM,QAAQ,GAAI,OAAO,CAAC,MAA2B,CAAC,kBAAkB,IAAI,UAAU,CAAC;YAEvF,2DAA2D;YAC3D,gFAAgF;YAChF,6EAA6E;YAC7E,IAAI,QAAQ,KAAK,MAAM,EAAE,CAAC;gBACxB,MAAM,QAAQ,GAAG,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBACxC,OAAO,UAAU,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;YACtC,CAAC;YAED,4EAA4E;YAC5E,4EAA4E;YAC5E,mEAAmE;YACnE,yEAAyE;YACzE,6EAA6E;YAC7E,6EAA6E;YAC7E,oDAAoD;YACpD,EAAE;YACF,sDAAsD;YACtD,yEAAyE;YACzE,kEAAkE;YAClE,IAAI,OAAqB,CAAC;YAC1B,IAAI,MAAM,CAAC,IAAI,KAAK,eAAe,IAAI,MAAM,CAAC,IAAI,KAAK,eAAe,EAAE,CAAC;gBACvE,MAAM,MAAM,GAAa;oBACvB,IAAI,EAAE,2BAA2B;oBACjC,EAAE,EAAE,MAAM;oBACV,GAAG,EAAE,OAAO,CAAC,GAAe;oBAC5B,QAAQ,EAAE,OAAO,CAAC,QAAQ;iBAC3B,CAAC;gBACF,OAAO,GAAG,SAAS,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;YAC9C,CAAC;iBAAM,CAAC;gBACN,+CAA+C;gBAC/C,MAAM,QAAQ,GAAG,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBACxC,OAAO,GAAG,UAAU,CAAC,QAAQ,EAAE,MAAM,CAAC,CAAC;YACzC,CAAC;YAED,4EAA4E;YAC5E,uEAAuE;YACvE,0DAA0D;YAC1D,IAAI,QAAQ,KAAK,UAAU,EAAE,CAAC;gBAC5B,OAAO,OAAO,CAAC;YACjB,CAAC;YAED,uEAAuE;YACvE,6CAA6C;YAC7C,4DAA4D;YAC5D,EAAE;YACF,8DAA8D;YAC9D,IAAI,aAAa,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAC/B,OAAO,OAAO,CAAC;YACjB,CAAC;YAED,sCAAsC;YACtC,MAAM,UAAU,GAAG,OAAO,CAAC,MAA0B,CAAC;YACtD,MAAM,YAAY,GAAG,OAAO,CAAC,QAA4B,CAAC;YAC1D,MAAM,YAAY,GAAG,kBAAkB,CAAC,YAAY,CAAC,YAAY,CAAC,IAAI,EAAE,MAAM,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YAExF,qFAAqF;YACrF,MAAM,UAAU,GAAG,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAC;YAE1C,qDAAqD;YACrD,MAAM,cAAc,GAAG,CAAC,EAAU,EAAY,EAAE,CAC9C,kBAAkB,CAChB,UAAU,CAAC,IAAI,EACf,UAAU,EACV,EAAE,EACF,QAA+B,EAC/B,CAAC,EACD,MAAM,CAAC,iBAAiB,CACzB,CAAC;YAEJ,qCAAqC;YACrC,MAAM,WAAW,GAAG,wBAAwB,CAAC,YAAY,EAAE,cAAc,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;YACtF,MAAM,SAAS,GAAG,kBAAkB,CAAC,MAAM,CAAC,QAAQ,CAAC,WAAW,CAAgB,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAE7F,gEAAgE;YAChE,2CAA2C;YAC3C,MAAM,MAAM,GAAG,kBAAkB,CAAC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAgB,EACxE,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAC7C,MAAM,UAAU,GAAG,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;YAEvC,IAAI,UAAU,GAAG,cAAc,CAAC,OAAO,EAAE,MAAM,CAAC,CAAC;YACjD,KAAK,MAAM,OAAO,IAAI,aAAa,EAAE,CAAC;gBACpC,MAAM,IAAI,GAAG,8BAA8B,CACzC,SAAS,EACT,MAAM,EACN,UAAU,EACV,OAAO,CAAC,GAAG,EACX,OAAO,CAAC,QAAQ,EAChB,CAAC,EACD,MAAM,CACP,CAAC;gBACF,MAAM,IAAI,GAAG,OAAO,CAAC,QAAQ,KAAK,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnD,UAAU,GAAG,eAAe,CAAC,UAAU,EAAE,IAAI,EAAE,IAAc,EAAE,MAAM,CAAC,CAAC;YACzE,CAAC;YAED,OAAO,SAAS,CAAC,OAAO,EAAE,UAAU,EAAE,MAAM,CAAC,CAAC;QAChD,CAAC;QAED,KAAK,UAAU,CAAC;QAChB,KAAK,YAAY;YACf,MAAM,IAAI,qBAAqB,CAC7B,cAAc,IAAI,CAAC,IAAI,6CAA6C;kBAClE,mDAAmD,CACtD,CAAC;QAEJ,OAAO,CAAC,CAAC,CAAC;YACR,MAAM,WAAW,GAAU,IAAI,CAAC;YAChC,KAAK,WAAW,CAAC;YACjB,MAAM,IAAI,qBAAqB,CAC7B,mCAAmC,IAAI,CAAC,SAAS,CAAE,IAA2B,CAAC,IAAI,CAAC,EAAE,CACvF,CAAC;QACJ,CAAC;IACH,CAAC;AACH,CAAC"}

package/dist/numerical/mathts-engine.d.ts

@@ -0,0 +1,55 @@
+import type { EngineTensor, TensorEngine, EinsumSpec, ForwardGradResult, ReverseGradResult } from './tensor-engine.js';
+import type { NestedArray } from './types.js';
+/**
+ * `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 declare class MathTSEngine implements TensorEngine {
+    readonly name = "MathTSEngine";
+    fromNested(data: NestedArray, shape: ReadonlyArray<number>): EngineTensor;
+    toNested(t: EngineTensor): NestedArray;
+    einsum(spec: EinsumSpec, ...operands: EngineTensor[]): EngineTensor;
+    matMul(a: EngineTensor, b: EngineTensor): EngineTensor;
+    transpose(t: EngineTensor, perm?: ReadonlyArray<number>): EngineTensor;
+    reshape(t: EngineTensor, shape: ReadonlyArray<number>): EngineTensor;
+    add(a: EngineTensor, b: EngineTensor): EngineTensor;
+    sub(a: EngineTensor, b: EngineTensor): EngineTensor;
+    mul(a: EngineTensor, b: EngineTensor): EngineTensor;
+    scale(t: EngineTensor, k: number): EngineTensor;
+    identity(n: number): EngineTensor;
+    normInf(t: EngineTensor): number;
+    /** Unwrap an EngineTensor to its underlying MathTS Tensor (for the AD boundary). */
+    private toMathTSTensor;
+    /** Wrap a MathTS Tensor (returned by autograd) back as an EngineTensor. */
+    private fromMathTSTensor;
+    /**
+     * 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.
+     */
+    forwardGrad(fn: (x: EngineTensor) => EngineTensor, x: EngineTensor): Promise<ForwardGradResult>;
+    /**
+     * 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.
+     */
+    reverseGrad(fn: (x: EngineTensor) => EngineTensor, x: EngineTensor, cotangent?: EngineTensor): Promise<ReverseGradResult>;
+}
+//# sourceMappingURL=mathts-engine.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"mathts-engine.d.ts","sourceRoot":"","sources":["../../src/numerical/mathts-engine.ts"],"names":[],"mappings":"AAYA,OAAO,KAAK,EAAE,YAAY,EAAE,YAAY,EAAE,UAAU,EAAE,iBAAiB,EAAE,iBAAiB,EAAE,MAAM,oBAAoB,CAAC;AACvH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,YAAY,CAAC;AAiB9C;;;;;;;;GAQG;AACH,qBAAa,YAAa,YAAW,YAAY;IAC/C,QAAQ,CAAC,IAAI,kBAAkB;IAE/B,UAAU,CAAC,IAAI,EAAE,WAAW,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY;IAGzE,QAAQ,CAAC,CAAC,EAAE,YAAY,GAAG,WAAW;IAItC,MAAM,CAAC,IAAI,EAAE,UAAU,EAAE,GAAG,QAAQ,EAAE,YAAY,EAAE,GAAG,YAAY;IAInE,MAAM,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY;IAGtD,SAAS,CAAC,CAAC,EAAE,YAAY,EAAE,IAAI,CAAC,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY;IAGtE,OAAO,CAAC,CAAC,EAAE,YAAY,EAAE,KAAK,EAAE,aAAa,CAAC,MAAM,CAAC,GAAG,YAAY;IAIpE,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY;IAGnD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY;IAGnD,GAAG,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,YAAY;IAGnD,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,MAAM,GAAG,YAAY;IAI/C,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,YAAY;IAGjC,OAAO,CAAC,CAAC,EAAE,YAAY,GAAG,MAAM;IAQhC,oFAAoF;IACpF,OAAO,CAAC,cAAc;IAItB,2EAA2E;IAC3E,OAAO,CAAC,gBAAgB;IAQxB;;;;;;;;;;;OAWG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,GACd,OAAO,CAAC,iBAAiB,CAAC;IA0B7B;;;;;;;;;OASG;IACG,WAAW,CACf,EAAE,EAAE,CAAC,CAAC,EAAE,YAAY,KAAK,YAAY,EACrC,CAAC,EAAE,YAAY,EACf,SAAS,CAAC,EAAE,YAAY,GACvB,OAAO,CAAC,iBAAiB,CAAC;CAoB9B"}