universal-physics-tensor 10.0.0 → 12.0.0

package/dist/atlas/benchmark/stats.js

@@ -13,8 +13,10 @@
  * @module atlas/benchmark/stats
  * @internal
  */
-/** The two-sided 95% normal quantile. @internal */
-export const Z95 = 1.959963984540054;
+import { binomialTest, normalQuantile } from '@danielsimonjr/mathts-functions';
+import { numberFrom } from '../../core/mathts-read.js';
+/** The two-sided 95% normal quantile, MathTS `normalQuantile(0.975)` (1.959963984540054). @internal */
+export const Z95 = normalQuantile(0.975);
 /**
  * Wilson score interval for a binomial proportion.
  *
@@ -39,13 +41,6 @@ export function wilsonInterval(successes, n, z = Z95) {
         upper: successes === n ? 1 : Math.min(1, centre + half),
     };
 }
-/** `ln C(n, k)`, summed; exact enough for the n here. */
-function logChoose(n, k) {
-    let s = 0;
-    for (let i = 1; i <= k; i++)
-        s += Math.log(n - k + i) - Math.log(i);
-    return s;
-}
 /**
  * McNemar's paired test. The EXACT binomial p-value is the one to report at
  * these sample sizes; the χ² is given because readers expect it.
@@ -54,14 +49,13 @@ function logChoose(n, k) {
  */
 export function mcnemar(table) {
     const n = table.b + table.c;
-    if (n === 0)
+    // No discordant pair, or an equal split: every outcome is at most as likely as the observed one, so p = 1 exactly.
+    if (table.b === table.c)
         return { chiSquared: 0, exactP: 1 };
     const chiSquared = Math.max(0, Math.abs(table.b - table.c) - 1) ** 2 / n;
-    const k = Math.min(table.b, table.c);
-    let tail = 0;
-    for (let i = 0; i <= k; i++)
-        tail += Math.exp(logChoose(n, i) - n * Math.LN2);
-    return { chiSquared, exactP: Math.min(1, 2 * tail) };
+    // MathTS `binomialTest` sums every outcome no more likely than the observed one; at p = ½ that is twice the smaller tail.
+    const exactP = numberFrom('binomialTest', binomialTest(Math.min(table.b, table.c), n).pValue);
+    return { chiSquared, exactP };
 }
 /**
  * Newcombe's method 10: a 95% interval for the PAIRED difference in proportions

package/dist/atlas/benchmark/stats.js.map

@@ -1 +1 @@
-{"version":3,"file":"stats.js","sourceRoot":"","sources":["../../../src/atlas/benchmark/stats.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,mDAAmD;AACnD,MAAM,CAAC,MAAM,GAAG,GAAG,iBAAiB,CAAC;AAQrC;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,SAAiB,EAAE,CAAS,EAAE,CAAC,GAAG,GAAG;IAClE,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,EAAE,CAAC,CAAC;IAC9G,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,SAAS,CAAC,IAAI,SAAS,GAAG,CAAC,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;QACnE,MAAM,IAAI,KAAK,CAAC,+DAA+D,SAAS,EAAE,CAAC,CAAC;IAC9F,CAAC;IACD,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,KAAK,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACzB,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3E,8EAA8E;IAC9E,wEAAwE;IACxE,OAAO;QACL,KAAK,EAAE,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC;QACvD,KAAK,EAAE,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC;KACxD,CAAC;AACJ,CAAC;AAkBD,yDAAyD;AACzD,SAAS,SAAS,CAAC,CAAS,EAAE,CAAS;IACrC,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;QAAE,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpE,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,OAAO,CAAC,KAAkB;IACxC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC5B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IACjD,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACzE,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IACrC,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE;QAAE,IAAI,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;IAC9E,OAAO,EAAE,UAAU,EAAE,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;AACvD,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,wBAAwB,CAAC,KAAkB,EAAE,CAAC,GAAG,GAAG;IAClE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,KAAK,CAAC;IAC7B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxB,IAAI,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;IACtE,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/D,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;QACd,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC5B,oEAAoE;QACpE,GAAG,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC;IACjE,CAAC;IACD,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IACrB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,GAAG,EAAE,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,GAAG,EAAE,CAAC;IAC1B,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;IACzF,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;IACzF,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;AACzE,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,WAAW,CAAC,MAAsC;IAChE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,IAAI,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,KAAK,CAAC,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IAC3H,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC,CAAC,CAAE,CAAC;YACzB,CAAC,IAAI,CAAC,CAAC;YACP,IAAI,CAAC,CAAC,CAAE,IAAI,CAAC,CAAC;YACd,IAAI,CAAC,CAAC,CAAE,IAAI,CAAC,CAAC;YACd,IAAI,CAAC,KAAK,CAAC;gBAAE,KAAK,IAAI,CAAC,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,IAAI,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;IACxD,MAAM,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC;IACrB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC;IAClE,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAChC,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;AAC9B,CAAC;AAaD;;;;;;GAMG;AACH,MAAM,UAAU,WAAW,CAAC,SAAiB,EAAE,eAAe,GAAG,GAAG;IAClE,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,eAAe,GAAG,SAAS,CAAC,CAAC;IAC1D,MAAM,QAAQ,GAAG,cAAc,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IACtD,MAAM,SAAS,GAAG,CAAC,QAAQ,CAAC,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACxD,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC;IAC9D,OAAO;QACL,SAAS;QACT,eAAe;QACf,QAAQ;QACR,SAAS;QACT,SAAS,EACP,QAAQ,SAAS,uBAAuB,GAAG,CAAC,eAAe,CAAC,qCAAqC;YACjG,OAAO,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,OAAO,GAAG,CAAC,SAAS,CAAC,uBAAuB;YAC9F,wCAAwC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,mCAAmC;YAC7F,sEAAsE;KACzE,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Atlas Phase 5, S5.4 — the benchmark's statistics and power report.\n *\n * Design note: `docs/planning/Atlas-Phase-5-Design.md` §8.\n *\n * Every function here is tested against a value computed independently, not\n * against its own output: the Wilson intervals against the plan's textbook\n * values, McNemar against a hand-computed 2×2 table, and κ against a hand-computed\n * confusion matrix. The paired-difference interval (Newcombe's method 10) has no\n * textbook value in the repository; it is tested by properties that pin it to\n * Wilson and to the unpaired square-and-add interval where those must agree.\n *\n * @module atlas/benchmark/stats\n * @internal\n */\n\n/** The two-sided 95% normal quantile. @internal */\nexport const Z95 = 1.959963984540054;\n\n/** A closed interval. @internal */\nexport interface Interval {\n  readonly lower: number;\n  readonly upper: number;\n}\n\n/**\n * Wilson score interval for a binomial proportion.\n *\n * @throws Error when `n` is not a positive integer or `successes` is outside `[0, n]`.\n * @internal\n */\nexport function wilsonInterval(successes: number, n: number, z = Z95): Interval {\n  if (!Number.isInteger(n) || n <= 0) throw new Error(`wilsonInterval: n must be a positive integer, got ${n}`);\n  if (!Number.isInteger(successes) || successes < 0 || successes > n) {\n    throw new Error(`wilsonInterval: successes must be an integer in [0, n], got ${successes}`);\n  }\n  const p = successes / n;\n  const z2 = z * z;\n  const denom = 1 + z2 / n;\n  const centre = (p + z2 / (2 * n)) / denom;\n  const half = (z / denom) * Math.sqrt((p * (1 - p)) / n + z2 / (4 * n * n));\n  // At x = 0 the lower bound, and at x = n the upper bound, are EXACTLY 0 and 1\n  // analytically; the formula reaches them only to within float rounding.\n  return {\n    lower: successes === 0 ? 0 : Math.max(0, centre - half),\n    upper: successes === n ? 1 : Math.min(1, centre + half),\n  };\n}\n\n/** A paired 2×2 table: a = both right, b = only A right, c = only B right, d = both wrong. @internal */\nexport interface PairedTable {\n  readonly a: number;\n  readonly b: number;\n  readonly c: number;\n  readonly d: number;\n}\n\n/** McNemar's test on the discordant cells. @internal */\nexport interface McNemarResult {\n  /** Continuity-corrected χ² = max(0, |b − c| − 1)² / (b + c); 0 when b + c = 0. The correction never goes negative (Edwards 1948). */\n  readonly chiSquared: number;\n  /** Exact two-sided p-value from Binomial(b + c, ½); 1 when b + c = 0. */\n  readonly exactP: number;\n}\n\n/** `ln C(n, k)`, summed; exact enough for the n here. */\nfunction logChoose(n: number, k: number): number {\n  let s = 0;\n  for (let i = 1; i <= k; i++) s += Math.log(n - k + i) - Math.log(i);\n  return s;\n}\n\n/**\n * McNemar's paired test. The EXACT binomial p-value is the one to report at\n * these sample sizes; the χ² is given because readers expect it.\n *\n * @internal\n */\nexport function mcnemar(table: PairedTable): McNemarResult {\n  const n = table.b + table.c;\n  if (n === 0) return { chiSquared: 0, exactP: 1 };\n  const chiSquared = Math.max(0, Math.abs(table.b - table.c) - 1) ** 2 / n;\n  const k = Math.min(table.b, table.c);\n  let tail = 0;\n  for (let i = 0; i <= k; i++) tail += Math.exp(logChoose(n, i) - n * Math.LN2);\n  return { chiSquared, exactP: Math.min(1, 2 * tail) };\n}\n\n/**\n * Newcombe's method 10: a 95% interval for the PAIRED difference in proportions\n * p_A − p_B = (b − c)/n, built from the two marginal Wilson intervals and a\n * correlation correction φ. Pre-registered criterion: the interval for the\n * atlas-minus-baseline difference must exclude zero.\n *\n * Newcombe, Statistics in Medicine 17 (1998) 2635.\n *\n * @internal\n */\nexport function pairedDifferenceInterval(table: PairedTable, z = Z95): Interval & { readonly diff: number } {\n  const { a, b, c, d } = table;\n  const n = a + b + c + d;\n  if (n === 0) throw new Error('pairedDifferenceInterval: empty table');\n  const p1 = (a + b) / n;\n  const p2 = (a + c) / n;\n  const w1 = wilsonInterval(a + b, n, z);\n  const w2 = wilsonInterval(a + c, n, z);\n  const denom = Math.sqrt((a + b) * (c + d) * (a + c) * (b + d));\n  let phi = 0;\n  if (denom > 0) {\n    const cross = a * d - b * c;\n    // Newcombe's continuity adjustment of φ, applied only when ad > bc.\n    phi = (cross > 0 ? Math.max(cross - n / 2, 0) : cross) / denom;\n  }\n  const diff = p1 - p2;\n  const dl1 = p1 - w1.lower;\n  const du1 = w1.upper - p1;\n  const dl2 = p2 - w2.lower;\n  const du2 = w2.upper - p2;\n  const lower = diff - Math.sqrt(Math.max(0, dl1 * dl1 - 2 * phi * dl1 * du2 + du2 * du2));\n  const upper = diff + Math.sqrt(Math.max(0, du1 * du1 - 2 * phi * du1 * dl2 + dl2 * dl2));\n  return { diff, lower: Math.max(-1, lower), upper: Math.min(1, upper) };\n}\n\n/**\n * Cohen's κ for two raters over the same items, from a square confusion matrix\n * (`matrix[i][j]` = items rater 1 put in category i and rater 2 in category j).\n *\n * @returns NaN when chance agreement is 1 (κ is undefined there), never 1.\n * @internal\n */\nexport function cohensKappa(matrix: readonly (readonly number[])[]): number {\n  const k = matrix.length;\n  if (k === 0 || matrix.some((row) => row.length !== k)) throw new Error('cohensKappa: matrix must be square and non-empty');\n  let n = 0;\n  let agree = 0;\n  const rows = new Array<number>(k).fill(0);\n  const cols = new Array<number>(k).fill(0);\n  for (let i = 0; i < k; i++) {\n    for (let j = 0; j < k; j++) {\n      const v = matrix[i]![j]!;\n      n += v;\n      rows[i]! += v;\n      cols[j]! += v;\n      if (i === j) agree += v;\n    }\n  }\n  if (n === 0) throw new Error('cohensKappa: no ratings');\n  const po = agree / n;\n  let pe = 0;\n  for (let i = 0; i < k; i++) pe += (rows[i]! / n) * (cols[i]! / n);\n  if (pe === 1) return Number.NaN;\n  return (po - pe) / (1 - pe);\n}\n\n/** What a frozen set of a given size can and cannot resolve. @internal */\nexport interface PowerReport {\n  readonly nPerClass: number;\n  /** The assumed per-class accuracy the interval is computed at. */\n  readonly assumedAccuracy: number;\n  readonly interval: Interval;\n  /** Half the interval width: differences smaller than about twice this are not resolvable. */\n  readonly halfWidth: number;\n  readonly statement: string;\n}\n\n/**\n * The interval a frozen set of `nPerClass` items per class can afford at an\n * assumed accuracy. This is the \"report the interval it can afford rather than\n * the claim it cannot\" line of ROADMAP Phase 5.\n *\n * @internal\n */\nexport function powerReport(nPerClass: number, assumedAccuracy = 0.8): PowerReport {\n  const successes = Math.round(assumedAccuracy * nPerClass);\n  const interval = wilsonInterval(successes, nPerClass);\n  const halfWidth = (interval.upper - interval.lower) / 2;\n  const pct = (x: number): string => `${(100 * x).toFixed(1)}%`;\n  return {\n    nPerClass,\n    assumedAccuracy,\n    interval,\n    halfWidth,\n    statement:\n      `With ${nPerClass} items per class at ${pct(assumedAccuracy)} accuracy, the 95% Wilson interval ` +\n      `is [${pct(interval.lower)}, ${pct(interval.upper)}] (±${pct(halfWidth)}). Two methods whose ` +\n      `accuracies differ by less than about ${pct(2 * halfWidth)} cannot be separated on this set ` +\n      'by unpaired intervals; the paired McNemar design is what narrows it.',\n  };\n}\n"]}
+{"version":3,"file":"stats.js","sourceRoot":"","sources":["../../../src/atlas/benchmark/stats.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,YAAY,EAAE,cAAc,EAAE,MAAM,iCAAiC,CAAC;AAC/E,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AAEvD,uGAAuG;AACvG,MAAM,CAAC,MAAM,GAAG,GAAW,cAAc,CAAC,KAAK,CAAC,CAAC;AAQjD;;;;;GAKG;AACH,MAAM,UAAU,cAAc,CAAC,SAAiB,EAAE,CAAS,EAAE,CAAC,GAAG,GAAG;IAClE,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,qDAAqD,CAAC,EAAE,CAAC,CAAC;IAC9G,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,SAAS,CAAC,IAAI,SAAS,GAAG,CAAC,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;QACnE,MAAM,IAAI,KAAK,CAAC,+DAA+D,SAAS,EAAE,CAAC,CAAC;IAC9F,CAAC;IACD,MAAM,CAAC,GAAG,SAAS,GAAG,CAAC,CAAC;IACxB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;IACjB,MAAM,KAAK,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;IACzB,MAAM,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IAC1C,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3E,8EAA8E;IAC9E,wEAAwE;IACxE,OAAO;QACL,KAAK,EAAE,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC;QACvD,KAAK,EAAE,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC;KACxD,CAAC;AACJ,CAAC;AAkBD;;;;;GAKG;AACH,MAAM,UAAU,OAAO,CAAC,KAAkB;IACxC,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAC5B,mHAAmH;IACnH,IAAI,KAAK,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC;QAAE,OAAO,EAAE,UAAU,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IAC7D,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IACzE,0HAA0H;IAC1H,MAAM,MAAM,GAAG,UAAU,CAAC,cAAc,EAAE,YAAY,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC;IAC9F,OAAO,EAAE,UAAU,EAAE,MAAM,EAAE,CAAC;AAChC,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,wBAAwB,CAAC,KAAkB,EAAE,CAAC,GAAG,GAAG;IAClE,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,KAAK,CAAC;IAC7B,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACxB,IAAI,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,uCAAuC,CAAC,CAAC;IACtE,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IACvB,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,EAAE,GAAG,cAAc,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACvC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/D,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,KAAK,GAAG,CAAC,EAAE,CAAC;QACd,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAC5B,oEAAoE;QACpE,GAAG,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,KAAK,CAAC;IACjE,CAAC;IACD,MAAM,IAAI,GAAG,EAAE,GAAG,EAAE,CAAC;IACrB,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,GAAG,EAAE,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC;IAC1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,GAAG,EAAE,CAAC;IAC1B,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;IACzF,MAAM,KAAK,GAAG,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC;IACzF,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC;AACzE,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,WAAW,CAAC,MAAsC;IAChE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;IACxB,IAAI,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,KAAK,CAAC,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;IAC3H,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAE,CAAC,CAAC,CAAE,CAAC;YACzB,CAAC,IAAI,CAAC,CAAC;YACP,IAAI,CAAC,CAAC,CAAE,IAAI,CAAC,CAAC;YACd,IAAI,CAAC,CAAC,CAAE,IAAI,CAAC,CAAC;YACd,IAAI,CAAC,KAAK,CAAC;gBAAE,KAAK,IAAI,CAAC,CAAC;QAC1B,CAAC;IACH,CAAC;IACD,IAAI,CAAC,KAAK,CAAC;QAAE,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;IACxD,MAAM,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC;IACrB,IAAI,EAAE,GAAG,CAAC,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;QAAE,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAE,GAAG,CAAC,CAAC,CAAC;IAClE,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IAChC,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;AAC9B,CAAC;AAaD;;;;;;GAMG;AACH,MAAM,UAAU,WAAW,CAAC,SAAiB,EAAE,eAAe,GAAG,GAAG;IAClE,MAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,eAAe,GAAG,SAAS,CAAC,CAAC;IAC1D,MAAM,QAAQ,GAAG,cAAc,CAAC,SAAS,EAAE,SAAS,CAAC,CAAC;IACtD,MAAM,SAAS,GAAG,CAAC,QAAQ,CAAC,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;IACxD,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,CAAC;IAC9D,OAAO;QACL,SAAS;QACT,eAAe;QACf,QAAQ;QACR,SAAS;QACT,SAAS,EACP,QAAQ,SAAS,uBAAuB,GAAG,CAAC,eAAe,CAAC,qCAAqC;YACjG,OAAO,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,QAAQ,CAAC,KAAK,CAAC,OAAO,GAAG,CAAC,SAAS,CAAC,uBAAuB;YAC9F,wCAAwC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,mCAAmC;YAC7F,sEAAsE;KACzE,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Atlas Phase 5, S5.4 — the benchmark's statistics and power report.\n *\n * Design note: `docs/planning/Atlas-Phase-5-Design.md` §8.\n *\n * Every function here is tested against a value computed independently, not\n * against its own output: the Wilson intervals against the plan's textbook\n * values, McNemar against a hand-computed 2×2 table, and κ against a hand-computed\n * confusion matrix. The paired-difference interval (Newcombe's method 10) has no\n * textbook value in the repository; it is tested by properties that pin it to\n * Wilson and to the unpaired square-and-add interval where those must agree.\n *\n * @module atlas/benchmark/stats\n * @internal\n */\n\nimport { binomialTest, normalQuantile } from '@danielsimonjr/mathts-functions';\nimport { numberFrom } from '../../core/mathts-read.js';\n\n/** The two-sided 95% normal quantile, MathTS `normalQuantile(0.975)` (1.959963984540054). @internal */\nexport const Z95: number = normalQuantile(0.975);\n\n/** A closed interval. @internal */\nexport interface Interval {\n  readonly lower: number;\n  readonly upper: number;\n}\n\n/**\n * Wilson score interval for a binomial proportion.\n *\n * @throws Error when `n` is not a positive integer or `successes` is outside `[0, n]`.\n * @internal\n */\nexport function wilsonInterval(successes: number, n: number, z = Z95): Interval {\n  if (!Number.isInteger(n) || n <= 0) throw new Error(`wilsonInterval: n must be a positive integer, got ${n}`);\n  if (!Number.isInteger(successes) || successes < 0 || successes > n) {\n    throw new Error(`wilsonInterval: successes must be an integer in [0, n], got ${successes}`);\n  }\n  const p = successes / n;\n  const z2 = z * z;\n  const denom = 1 + z2 / n;\n  const centre = (p + z2 / (2 * n)) / denom;\n  const half = (z / denom) * Math.sqrt((p * (1 - p)) / n + z2 / (4 * n * n));\n  // At x = 0 the lower bound, and at x = n the upper bound, are EXACTLY 0 and 1\n  // analytically; the formula reaches them only to within float rounding.\n  return {\n    lower: successes === 0 ? 0 : Math.max(0, centre - half),\n    upper: successes === n ? 1 : Math.min(1, centre + half),\n  };\n}\n\n/** A paired 2×2 table: a = both right, b = only A right, c = only B right, d = both wrong. @internal */\nexport interface PairedTable {\n  readonly a: number;\n  readonly b: number;\n  readonly c: number;\n  readonly d: number;\n}\n\n/** McNemar's test on the discordant cells. @internal */\nexport interface McNemarResult {\n  /** Continuity-corrected χ² = max(0, |b − c| − 1)² / (b + c); 0 when b + c = 0. The correction never goes negative (Edwards 1948). */\n  readonly chiSquared: number;\n  /** Exact two-sided p-value from Binomial(b + c, ½); 1 when b + c = 0. */\n  readonly exactP: number;\n}\n\n/**\n * McNemar's paired test. The EXACT binomial p-value is the one to report at\n * these sample sizes; the χ² is given because readers expect it.\n *\n * @internal\n */\nexport function mcnemar(table: PairedTable): McNemarResult {\n  const n = table.b + table.c;\n  // No discordant pair, or an equal split: every outcome is at most as likely as the observed one, so p = 1 exactly.\n  if (table.b === table.c) return { chiSquared: 0, exactP: 1 };\n  const chiSquared = Math.max(0, Math.abs(table.b - table.c) - 1) ** 2 / n;\n  // MathTS `binomialTest` sums every outcome no more likely than the observed one; at p = ½ that is twice the smaller tail.\n  const exactP = numberFrom('binomialTest', binomialTest(Math.min(table.b, table.c), n).pValue);\n  return { chiSquared, exactP };\n}\n\n/**\n * Newcombe's method 10: a 95% interval for the PAIRED difference in proportions\n * p_A − p_B = (b − c)/n, built from the two marginal Wilson intervals and a\n * correlation correction φ. Pre-registered criterion: the interval for the\n * atlas-minus-baseline difference must exclude zero.\n *\n * Newcombe, Statistics in Medicine 17 (1998) 2635.\n *\n * @internal\n */\nexport function pairedDifferenceInterval(table: PairedTable, z = Z95): Interval & { readonly diff: number } {\n  const { a, b, c, d } = table;\n  const n = a + b + c + d;\n  if (n === 0) throw new Error('pairedDifferenceInterval: empty table');\n  const p1 = (a + b) / n;\n  const p2 = (a + c) / n;\n  const w1 = wilsonInterval(a + b, n, z);\n  const w2 = wilsonInterval(a + c, n, z);\n  const denom = Math.sqrt((a + b) * (c + d) * (a + c) * (b + d));\n  let phi = 0;\n  if (denom > 0) {\n    const cross = a * d - b * c;\n    // Newcombe's continuity adjustment of φ, applied only when ad > bc.\n    phi = (cross > 0 ? Math.max(cross - n / 2, 0) : cross) / denom;\n  }\n  const diff = p1 - p2;\n  const dl1 = p1 - w1.lower;\n  const du1 = w1.upper - p1;\n  const dl2 = p2 - w2.lower;\n  const du2 = w2.upper - p2;\n  const lower = diff - Math.sqrt(Math.max(0, dl1 * dl1 - 2 * phi * dl1 * du2 + du2 * du2));\n  const upper = diff + Math.sqrt(Math.max(0, du1 * du1 - 2 * phi * du1 * dl2 + dl2 * dl2));\n  return { diff, lower: Math.max(-1, lower), upper: Math.min(1, upper) };\n}\n\n/**\n * Cohen's κ for two raters over the same items, from a square confusion matrix\n * (`matrix[i][j]` = items rater 1 put in category i and rater 2 in category j).\n *\n * @returns NaN when chance agreement is 1 (κ is undefined there), never 1.\n * @internal\n */\nexport function cohensKappa(matrix: readonly (readonly number[])[]): number {\n  const k = matrix.length;\n  if (k === 0 || matrix.some((row) => row.length !== k)) throw new Error('cohensKappa: matrix must be square and non-empty');\n  let n = 0;\n  let agree = 0;\n  const rows = new Array<number>(k).fill(0);\n  const cols = new Array<number>(k).fill(0);\n  for (let i = 0; i < k; i++) {\n    for (let j = 0; j < k; j++) {\n      const v = matrix[i]![j]!;\n      n += v;\n      rows[i]! += v;\n      cols[j]! += v;\n      if (i === j) agree += v;\n    }\n  }\n  if (n === 0) throw new Error('cohensKappa: no ratings');\n  const po = agree / n;\n  let pe = 0;\n  for (let i = 0; i < k; i++) pe += (rows[i]! / n) * (cols[i]! / n);\n  if (pe === 1) return Number.NaN;\n  return (po - pe) / (1 - pe);\n}\n\n/** What a frozen set of a given size can and cannot resolve. @internal */\nexport interface PowerReport {\n  readonly nPerClass: number;\n  /** The assumed per-class accuracy the interval is computed at. */\n  readonly assumedAccuracy: number;\n  readonly interval: Interval;\n  /** Half the interval width: differences smaller than about twice this are not resolvable. */\n  readonly halfWidth: number;\n  readonly statement: string;\n}\n\n/**\n * The interval a frozen set of `nPerClass` items per class can afford at an\n * assumed accuracy. This is the \"report the interval it can afford rather than\n * the claim it cannot\" line of ROADMAP Phase 5.\n *\n * @internal\n */\nexport function powerReport(nPerClass: number, assumedAccuracy = 0.8): PowerReport {\n  const successes = Math.round(assumedAccuracy * nPerClass);\n  const interval = wilsonInterval(successes, nPerClass);\n  const halfWidth = (interval.upper - interval.lower) / 2;\n  const pct = (x: number): string => `${(100 * x).toFixed(1)}%`;\n  return {\n    nPerClass,\n    assumedAccuracy,\n    interval,\n    halfWidth,\n    statement:\n      `With ${nPerClass} items per class at ${pct(assumedAccuracy)} accuracy, the 95% Wilson interval ` +\n      `is [${pct(interval.lower)}, ${pct(interval.upper)}] (±${pct(halfWidth)}). Two methods whose ` +\n      `accuracies differ by less than about ${pct(2 * halfWidth)} cannot be separated on this set ` +\n      'by unpaired intervals; the paired McNemar design is what narrows it.',\n  };\n}\n"]}

package/dist/atlas/catalog-formal-ref.d.ts

@@ -1,13 +1,17 @@
 /**
  * Catalog formal references.
  *
- * The catalog row stores the manifest key. This module asks PhysJS for that
- * key. It does not parse an id out of text, and it does not store an evidence
- * tag.
+ * The catalog row stores the manifest key, and beside it two projections of
+ * the manifest: the theorem that key names and the record's proof status.
+ * This module builds the reference from the key, derives the status from
+ * that reference, and reports a row whose stored projections disagree with
+ * the derivation. It does not parse an id out of text, and it does not store
+ * an evidence tag.
  *
  * @module atlas/catalog-formal-ref
  */
 import type { FormalRef } from '../relations/types.js';
+import type { CatalogEntry, ProofStatus } from '../bridges/catalog-types.js';
 /**
  * The reviewed PhysJS reference for a catalog id, or `undefined` when the
  * row records no manifest key.
@@ -15,4 +19,23 @@ import type { FormalRef } from '../relations/types.js';
  * @internal
  */
 export declare function catalogFormalRef(id: number): FormalRef | undefined;
+/**
+ * The proof status a record earns, derived from its reference: `pending` with
+ * no formalKey, `full` for a reviewed kind-`bridge` reference, `weak` for any
+ * other reference. `reference` is the record's {@link catalogFormalRef}.
+ *
+ * @internal
+ */
+export declare function deriveProofStatus(reference: FormalRef | undefined): ProofStatus;
+/**
+ * Records whose stored `theorem` or `proofStatus` disagree with the manifest:
+ * a theorem that is not the one the manifest entry of `formalKey` names, a
+ * theorem on a record with no `formalKey` (or none on a record with one), and
+ * a status other than {@link deriveProofStatus}. Empty means every record's
+ * projections are the manifest's. `referenceOf` defaults to the manifest's
+ * reference for the key; a test passes a changed one.
+ *
+ * @internal
+ */
+export declare function catalogProofProblems(entries?: readonly CatalogEntry[], referenceOf?: (key: string) => FormalRef): string[];
 //# sourceMappingURL=catalog-formal-ref.d.ts.map

package/dist/atlas/catalog-formal-ref.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"catalog-formal-ref.d.ts","sourceRoot":"","sources":["../../src/atlas/catalog-formal-ref.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,uBAAuB,CAAC;AAOvD;;;;;GAKG;AACH,wBAAgB,gBAAgB,CAAC,EAAE,EAAE,MAAM,GAAG,SAAS,GAAG,SAAS,CAMlE"}
+{"version":3,"file":"catalog-formal-ref.d.ts","sourceRoot":"","sources":["../../src/atlas/catalog-formal-ref.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,uBAAuB,CAAC;AACvD,OAAO,KAAK,EAAE,YAAY,EAAE,WAAW,EAAE,MAAM,6BAA6B,CAAC;AAO7E;;;;;GAKG;AACH,wBAAgB,gBAAgB,CAAC,EAAE,EAAE,MAAM,GAAG,SAAS,GAAG,SAAS,CAUlE;AAED;;;;;;GAMG;AACH,wBAAgB,iBAAiB,CAAC,SAAS,EAAE,SAAS,GAAG,SAAS,GAAG,WAAW,CAG/E;AAED;;;;;;;;;GASG;AACH,wBAAgB,oBAAoB,CAClC,OAAO,GAAE,SAAS,YAAY,EAAqB,EACnD,WAAW,GAAE,CAAC,GAAG,EAAE,MAAM,KAAK,SAA2B,GACxD,MAAM,EAAE,CAeV"}

package/dist/atlas/catalog-formal-ref.js

@@ -1,13 +1,16 @@
 /**
  * Catalog formal references.
  *
- * The catalog row stores the manifest key. This module asks PhysJS for that
- * key. It does not parse an id out of text, and it does not store an evidence
- * tag.
+ * The catalog row stores the manifest key, and beside it two projections of
+ * the manifest: the theorem that key names and the record's proof status.
+ * This module builds the reference from the key, derives the status from
+ * that reference, and reports a row whose stored projections disagree with
+ * the derivation. It does not parse an id out of text, and it does not store
+ * an evidence tag.
  *
  * @module atlas/catalog-formal-ref
  */
-import { catalogEntry } from '../bridges/catalog-load.js';
+import { catalogEntries, catalogEntry } from '../bridges/catalog-load.js';
 import { physjsFormalRef } from './physjs-ref.js';
 /** Memo of the PhysJS reference already built for a catalog id. */
 const byId = new Map();
@@ -20,9 +23,53 @@ const byId = new Map();
 export function catalogFormalRef(id) {
     if (byId.has(id))
         return byId.get(id);
-    const key = catalogEntry(id)?.formalKey;
+    const entry = catalogEntry(id);
+    const key = entry?.formalKey;
     const ref = key === undefined ? undefined : physjsFormalRef(key);
+    if (ref !== undefined && entry?.theorem !== ref.statement) {
+        throw new Error(`catalog id ${id} names theorem '${entry?.theorem}', the manifest entry '${key}' names '${ref.statement}'`);
+    }
     byId.set(id, ref);
     return ref;
 }
+/**
+ * The proof status a record earns, derived from its reference: `pending` with
+ * no formalKey, `full` for a reviewed kind-`bridge` reference, `weak` for any
+ * other reference. `reference` is the record's {@link catalogFormalRef}.
+ *
+ * @internal
+ */
+export function deriveProofStatus(reference) {
+    if (reference === undefined)
+        return 'pending';
+    return reference.kind === 'bridge' && reference.fidelity !== 'unreviewed' ? 'full' : 'weak';
+}
+/**
+ * Records whose stored `theorem` or `proofStatus` disagree with the manifest:
+ * a theorem that is not the one the manifest entry of `formalKey` names, a
+ * theorem on a record with no `formalKey` (or none on a record with one), and
+ * a status other than {@link deriveProofStatus}. Empty means every record's
+ * projections are the manifest's. `referenceOf` defaults to the manifest's
+ * reference for the key; a test passes a changed one.
+ *
+ * @internal
+ */
+export function catalogProofProblems(entries = catalogEntries(), referenceOf = physjsFormalRef) {
+    const problems = [];
+    for (const entry of entries) {
+        const reference = entry.formalKey === undefined ? undefined : referenceOf(entry.formalKey);
+        if (reference === undefined) {
+            if (entry.theorem !== undefined)
+                problems.push(`catalog id ${entry.id} names theorem '${entry.theorem}' and has no formalKey`);
+        }
+        else if (entry.theorem !== reference.statement) {
+            problems.push(`catalog id ${entry.id} theorem is '${entry.theorem ?? '(none)'}', the manifest entry '${entry.formalKey}' names '${reference.statement}'`);
+        }
+        const derived = deriveProofStatus(reference);
+        if (entry.proofStatus !== derived) {
+            problems.push(`catalog id ${entry.id} proofStatus is '${entry.proofStatus}', derived '${derived}'`);
+        }
+    }
+    return problems;
+}
 //# sourceMappingURL=catalog-formal-ref.js.map

package/dist/atlas/catalog-formal-ref.js.map

@@ -1 +1 @@
-{"version":3,"file":"catalog-formal-ref.js","sourceRoot":"","sources":["../../src/atlas/catalog-formal-ref.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AAGH,OAAO,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAC1D,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAElD,mEAAmE;AACnE,MAAM,IAAI,GAAG,IAAI,GAAG,EAAiC,CAAC;AAEtD;;;;;GAKG;AACH,MAAM,UAAU,gBAAgB,CAAC,EAAU;IACzC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAAE,OAAO,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACtC,MAAM,GAAG,GAAG,YAAY,CAAC,EAAE,CAAC,EAAE,SAAS,CAAC;IACxC,MAAM,GAAG,GAAG,GAAG,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;IACjE,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAClB,OAAO,GAAG,CAAC;AACb,CAAC","sourcesContent":["/**\n * Catalog formal references.\n *\n * The catalog row stores the manifest key. This module asks PhysJS for that\n * key. It does not parse an id out of text, and it does not store an evidence\n * tag.\n *\n * @module atlas/catalog-formal-ref\n */\n\nimport type { FormalRef } from '../relations/types.js';\nimport { catalogEntry } from '../bridges/catalog-load.js';\nimport { physjsFormalRef } from './physjs-ref.js';\n\n/** Memo of the PhysJS reference already built for a catalog id. */\nconst byId = new Map<number, FormalRef | undefined>();\n\n/**\n * The reviewed PhysJS reference for a catalog id, or `undefined` when the\n * row records no manifest key.\n *\n * @internal\n */\nexport function catalogFormalRef(id: number): FormalRef | undefined {\n  if (byId.has(id)) return byId.get(id);\n  const key = catalogEntry(id)?.formalKey;\n  const ref = key === undefined ? undefined : physjsFormalRef(key);\n  byId.set(id, ref);\n  return ref;\n}\n"]}
+{"version":3,"file":"catalog-formal-ref.js","sourceRoot":"","sources":["../../src/atlas/catalog-formal-ref.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,EAAE,cAAc,EAAE,YAAY,EAAE,MAAM,4BAA4B,CAAC;AAC1E,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAElD,mEAAmE;AACnE,MAAM,IAAI,GAAG,IAAI,GAAG,EAAiC,CAAC;AAEtD;;;;;GAKG;AACH,MAAM,UAAU,gBAAgB,CAAC,EAAU;IACzC,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;QAAE,OAAO,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACtC,MAAM,KAAK,GAAG,YAAY,CAAC,EAAE,CAAC,CAAC;IAC/B,MAAM,GAAG,GAAG,KAAK,EAAE,SAAS,CAAC;IAC7B,MAAM,GAAG,GAAG,GAAG,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,eAAe,CAAC,GAAG,CAAC,CAAC;IACjE,IAAI,GAAG,KAAK,SAAS,IAAI,KAAK,EAAE,OAAO,KAAK,GAAG,CAAC,SAAS,EAAE,CAAC;QAC1D,MAAM,IAAI,KAAK,CAAC,cAAc,EAAE,mBAAmB,KAAK,EAAE,OAAO,0BAA0B,GAAG,YAAY,GAAG,CAAC,SAAS,GAAG,CAAC,CAAC;IAC9H,CAAC;IACD,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;IAClB,OAAO,GAAG,CAAC;AACb,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,SAAgC;IAChE,IAAI,SAAS,KAAK,SAAS;QAAE,OAAO,SAAS,CAAC;IAC9C,OAAO,SAAS,CAAC,IAAI,KAAK,QAAQ,IAAI,SAAS,CAAC,QAAQ,KAAK,YAAY,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC;AAC9F,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,oBAAoB,CAClC,OAAO,GAA4B,cAAc,EAAE,EACnD,WAAW,GAA+B,eAAe;IAEzD,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,KAAK,MAAM,KAAK,IAAI,OAAO,EAAE,CAAC;QAC5B,MAAM,SAAS,GAAG,KAAK,CAAC,SAAS,KAAK,SAAS,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,WAAW,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;QAC3F,IAAI,SAAS,KAAK,SAAS,EAAE,CAAC;YAC5B,IAAI,KAAK,CAAC,OAAO,KAAK,SAAS;gBAAE,QAAQ,CAAC,IAAI,CAAC,cAAc,KAAK,CAAC,EAAE,mBAAmB,KAAK,CAAC,OAAO,wBAAwB,CAAC,CAAC;QACjI,CAAC;aAAM,IAAI,KAAK,CAAC,OAAO,KAAK,SAAS,CAAC,SAAS,EAAE,CAAC;YACjD,QAAQ,CAAC,IAAI,CAAC,cAAc,KAAK,CAAC,EAAE,gBAAgB,KAAK,CAAC,OAAO,IAAI,QAAQ,0BAA0B,KAAK,CAAC,SAAS,YAAY,SAAS,CAAC,SAAS,GAAG,CAAC,CAAC;QAC5J,CAAC;QACD,MAAM,OAAO,GAAG,iBAAiB,CAAC,SAAS,CAAC,CAAC;QAC7C,IAAI,KAAK,CAAC,WAAW,KAAK,OAAO,EAAE,CAAC;YAClC,QAAQ,CAAC,IAAI,CAAC,cAAc,KAAK,CAAC,EAAE,oBAAoB,KAAK,CAAC,WAAW,eAAe,OAAO,GAAG,CAAC,CAAC;QACtG,CAAC;IACH,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC","sourcesContent":["/**\n * Catalog formal references.\n *\n * The catalog row stores the manifest key, and beside it two projections of\n * the manifest: the theorem that key names and the record's proof status.\n * This module builds the reference from the key, derives the status from\n * that reference, and reports a row whose stored projections disagree with\n * the derivation. It does not parse an id out of text, and it does not store\n * an evidence tag.\n *\n * @module atlas/catalog-formal-ref\n */\n\nimport type { FormalRef } from '../relations/types.js';\nimport type { CatalogEntry, ProofStatus } from '../bridges/catalog-types.js';\nimport { catalogEntries, catalogEntry } from '../bridges/catalog-load.js';\nimport { physjsFormalRef } from './physjs-ref.js';\n\n/** Memo of the PhysJS reference already built for a catalog id. */\nconst byId = new Map<number, FormalRef | undefined>();\n\n/**\n * The reviewed PhysJS reference for a catalog id, or `undefined` when the\n * row records no manifest key.\n *\n * @internal\n */\nexport function catalogFormalRef(id: number): FormalRef | undefined {\n  if (byId.has(id)) return byId.get(id);\n  const entry = catalogEntry(id);\n  const key = entry?.formalKey;\n  const ref = key === undefined ? undefined : physjsFormalRef(key);\n  if (ref !== undefined && entry?.theorem !== ref.statement) {\n    throw new Error(`catalog id ${id} names theorem '${entry?.theorem}', the manifest entry '${key}' names '${ref.statement}'`);\n  }\n  byId.set(id, ref);\n  return ref;\n}\n\n/**\n * The proof status a record earns, derived from its reference: `pending` with\n * no formalKey, `full` for a reviewed kind-`bridge` reference, `weak` for any\n * other reference. `reference` is the record's {@link catalogFormalRef}.\n *\n * @internal\n */\nexport function deriveProofStatus(reference: FormalRef | undefined): ProofStatus {\n  if (reference === undefined) return 'pending';\n  return reference.kind === 'bridge' && reference.fidelity !== 'unreviewed' ? 'full' : 'weak';\n}\n\n/**\n * Records whose stored `theorem` or `proofStatus` disagree with the manifest:\n * a theorem that is not the one the manifest entry of `formalKey` names, a\n * theorem on a record with no `formalKey` (or none on a record with one), and\n * a status other than {@link deriveProofStatus}. Empty means every record's\n * projections are the manifest's. `referenceOf` defaults to the manifest's\n * reference for the key; a test passes a changed one.\n *\n * @internal\n */\nexport function catalogProofProblems(\n  entries: readonly CatalogEntry[] = catalogEntries(),\n  referenceOf: (key: string) => FormalRef = physjsFormalRef,\n): string[] {\n  const problems: string[] = [];\n  for (const entry of entries) {\n    const reference = entry.formalKey === undefined ? undefined : referenceOf(entry.formalKey);\n    if (reference === undefined) {\n      if (entry.theorem !== undefined) problems.push(`catalog id ${entry.id} names theorem '${entry.theorem}' and has no formalKey`);\n    } else if (entry.theorem !== reference.statement) {\n      problems.push(`catalog id ${entry.id} theorem is '${entry.theorem ?? '(none)'}', the manifest entry '${entry.formalKey}' names '${reference.statement}'`);\n    }\n    const derived = deriveProofStatus(reference);\n    if (entry.proofStatus !== derived) {\n      problems.push(`catalog id ${entry.id} proofStatus is '${entry.proofStatus}', derived '${derived}'`);\n    }\n  }\n  return problems;\n}\n"]}

package/dist/atlas/chain-pipeline.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"chain-pipeline.d.ts","sourceRoot":"","sources":["../../src/atlas/chain-pipeline.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAC5D,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAIzD,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,wBAAwB,CAAC;AAKzD,OAAO,EAGL,KAAK,WAAW,EACjB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAsB,KAAK,mBAAmB,EAAE,MAAM,gCAAgC,CAAC;AAK9F,2FAA2F;AAC3F,MAAM,WAAW,uBAAuB;IACtC,QAAQ,CAAC,IAAI,EAAE,cAAc,CAAC;IAC9B,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;CACrC;AAED,6FAA6F;AAC7F,MAAM,WAAW,sBAAsB;IACrC,QAAQ,CAAC,IAAI,EAAE,aAAa,CAAC;IAC7B,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,QAAQ,CAAC,cAAc,EAAE,MAAM,CAAC;IAChC,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;CACrC;AAED,mEAAmE;AACnE,MAAM,WAAW,eAAe;IAC9B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;IACpC,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;CACvB;AAED;;;;;;;GAOG;AACH,MAAM,WAAW,uBAAuB;IACtC,QAAQ,CAAC,IAAI,EAAE,6BAA6B,CAAC;IAC7C,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC5C,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;CAC1B;AAED,0FAA0F;AAC1F,MAAM,MAAM,mBAAmB,GAC3B,uBAAuB,GACvB,sBAAsB,GACtB,eAAe,GACf,mBAAmB,GACnB,uBAAuB,CAAC;AAY5B;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,IAAI,EAAE,QAAQ,EAAE,UAAU,EAAE,MAAM,GAAG;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,SAAS,CAAA;CAAE,EAAE,CASlG;AAyCD;;;;;;;;;GASG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,WAAW,EACnB,YAAY,GAAE,SAAS,MAAM,EAAO,GACnC,mBAAmB,CAyBrB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gBAAgB,CAAC,KAAK,EAAE,SAAS,UAAU,EAAE,GAAG,SAAS,mBAAmB,EAAE,CAqD7F"}
+{"version":3,"file":"chain-pipeline.d.ts","sourceRoot":"","sources":["../../src/atlas/chain-pipeline.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAGH,OAAO,KAAK,EAAE,QAAQ,EAAE,MAAM,6BAA6B,CAAC;AAC5D,OAAO,KAAK,EAAE,SAAS,EAAE,MAAM,yBAAyB,CAAC;AAKzD,OAAO,KAAK,EAAE,UAAU,EAAE,MAAM,wBAAwB,CAAC;AAKzD,OAAO,EAGL,KAAK,WAAW,EACjB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAsB,KAAK,mBAAmB,EAAE,MAAM,gCAAgC,CAAC;AAK9F,2FAA2F;AAC3F,MAAM,WAAW,uBAAuB;IACtC,QAAQ,CAAC,IAAI,EAAE,cAAc,CAAC;IAC9B,QAAQ,CAAC,SAAS,EAAE,MAAM,CAAC;IAC3B,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;CACrC;AAED,6FAA6F;AAC7F,MAAM,WAAW,sBAAsB;IACrC,QAAQ,CAAC,IAAI,EAAE,aAAa,CAAC;IAC7B,QAAQ,CAAC,WAAW,EAAE,MAAM,CAAC;IAC7B,QAAQ,CAAC,cAAc,EAAE,MAAM,CAAC;IAChC,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;CACrC;AAED,mEAAmE;AACnE,MAAM,WAAW,eAAe;IAC9B,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;IACtB,QAAQ,CAAC,EAAE,EAAE,MAAM,CAAC;IACpB,QAAQ,CAAC,OAAO,EAAE,SAAS,MAAM,EAAE,CAAC;IACpC,QAAQ,CAAC,IAAI,EAAE,MAAM,CAAC;CACvB;AAED;;;;;;;GAOG;AACH,MAAM,WAAW,uBAAuB;IACtC,QAAQ,CAAC,IAAI,EAAE,6BAA6B,CAAC;IAC7C,QAAQ,CAAC,OAAO,EAAE,SAAS,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC5C,QAAQ,CAAC,OAAO,EAAE,MAAM,CAAC;CAC1B;AAED,0FAA0F;AAC1F,MAAM,MAAM,mBAAmB,GAC3B,uBAAuB,GACvB,sBAAsB,GACtB,eAAe,GACf,mBAAmB,GACnB,uBAAuB,CAAC;AAQ5B;;;;;;;GAOG;AACH,wBAAgB,WAAW,CAAC,IAAI,EAAE,QAAQ,EAAE,UAAU,EAAE,MAAM,GAAG;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,SAAS,CAAA;CAAE,EAAE,CASlG;AAyCD;;;;;;;;;GASG;AACH,wBAAgB,iBAAiB,CAC/B,MAAM,EAAE,WAAW,EACnB,YAAY,GAAE,SAAS,MAAM,EAAO,GACnC,mBAAmB,CAyBrB;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,gBAAgB,CAAC,KAAK,EAAE,SAAS,UAAU,EAAE,GAAG,SAAS,mBAAmB,EAAE,CAqD7F"}

package/dist/atlas/chain-pipeline.js

@@ -13,9 +13,11 @@
  *
  * @module atlas/chain-pipeline
  */
+import { byKey } from '../core/compare.js';
 import { DIMENSIONLESS } from '../dimensional/types.js';
 import { equals } from '../dimensional/algebra.js';
-import { CONSTANTS } from '../dimensional/symbolic-constants.js';
+import { CONSTANTS, isRegisteredConstant } from '../dimensional/symbolic-constants.js';
+import { isDecimalLiteral } from '../dimensional/ast-leaves.js';
 import { enumerateCompositionsWithRefusals } from '../composition/enumerate.js';
 import { buckinghamFilter } from '../composition/buckingham-filter.js';
 import { matchChain } from '../composition/chain-match.js';
@@ -25,14 +27,11 @@ import { joinRegimeMismatch } from '../composition/chain-regime.js';
 import { bridgeSeedKeys, physjsTheorem } from './physjs-ref.js';
 import { emitProofTarget } from './proof-target.js';
 import { scalarSymbolsFromMathTs } from '../composition/mathts-scalar-symbols.js';
-function isNumericName(name) {
-    return /^\d+(\.\d+)?$/.test(name);
-}
+/** A registered constant whose dimension is none, or a `pi` multiple (`2pi`). Own keys only. */
 function isDimensionlessConstant(name) {
     if (/^\d*pi$/.test(name))
         return true;
-    const registered = CONSTANTS[name];
-    return registered !== undefined && equals(registered.dim, DIMENSIONLESS);
+    return isRegisteredConstant(name) && equals(CONSTANTS[name].dim, DIMENSIONLESS);
 }
 /**
  * Dimensioned leaves of the composed formula, excluding the target name.
@@ -47,11 +46,11 @@ export function governingOf(expr, targetName) {
     for (const leaf of scalarSymbolsFromMathTs(expr)) {
         if (leaf.name === targetName)
             continue;
-        if (isNumericName(leaf.name) || isDimensionlessConstant(leaf.name))
+        if (isDecimalLiteral(leaf.name) || isDimensionlessConstant(leaf.name))
             continue;
         vars.push(leaf);
     }
-    vars.sort((a, b) => (a.name < b.name ? -1 : a.name > b.name ? 1 : 0));
+    vars.sort(byKey((v) => v.name));
     return vars;
 }
 function theoremsFor(edgeIds) {

package/dist/atlas/chain-pipeline.js.map

@@ -1 +1 @@
-{"version":3,"file":"chain-pipeline.js","sourceRoot":"","sources":["../../src/atlas/chain-pipeline.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAIH,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AACxD,OAAO,EAAE,MAAM,EAAE,MAAM,2BAA2B,CAAC;AACnD,OAAO,EAAE,SAAS,EAAE,MAAM,sCAAsC,CAAC;AAEjE,OAAO,EAAE,iCAAiC,EAAE,MAAM,6BAA6B,CAAC;AAChF,OAAO,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AACvE,OAAO,EAAE,UAAU,EAAE,MAAM,+BAA+B,CAAC;AAC3D,OAAO,EAAE,mBAAmB,EAAuB,MAAM,mCAAmC,CAAC;AAC7F,OAAO,EACL,aAAa,EACb,iBAAiB,GAElB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAE,kBAAkB,EAA4B,MAAM,gCAAgC,CAAC;AAC9F,OAAO,EAAE,cAAc,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChE,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,uBAAuB,EAAE,MAAM,yCAAyC,CAAC;AA+ClF,SAAS,aAAa,CAAC,IAAY;IACjC,OAAO,eAAe,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;AACpC,CAAC;AAED,SAAS,uBAAuB,CAAC,IAAY;IAC3C,IAAI,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,IAAI,CAAC;IACtC,MAAM,UAAU,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;IACnC,OAAO,UAAU,KAAK,SAAS,IAAI,MAAM,CAAC,UAAU,CAAC,GAAG,EAAE,aAAa,CAAC,CAAC;AAC3E,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,WAAW,CAAC,IAAc,EAAE,UAAkB;IAC5D,MAAM,IAAI,GAAuC,EAAE,CAAC;IACpD,KAAK,MAAM,IAAI,IAAI,uBAAuB,CAAC,IAAI,CAAC,EAAE,CAAC;QACjD,IAAI,IAAI,CAAC,IAAI,KAAK,UAAU;YAAE,SAAS;QACvC,IAAI,aAAa,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,uBAAuB,CAAC,IAAI,CAAC,IAAI,CAAC;YAAE,SAAS;QAC7E,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAClB,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACtE,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,WAAW,CAAC,OAA0B;IAC7C,OAAO,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE;QACxB,MAAM,IAAI,GAAG,aAAa,CAAC,EAAE,CAAC,CAAC;QAC/B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAM,IAAI,KAAK,CAAC,0DAA0D,EAAE,GAAG,CAAC,CAAC;QACnF,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC,CAAC,CAAC;AACL,CAAC;AAED,SAAS,cAAc,CAAC,MAAmB;IACzC,MAAM,IAAI,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACnC,MAAM,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,SAAS,CAAC;IAC5C,MAAM,OAAO,GAAG,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC;IACzD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,cAAc,EAAE,CAAC;QAClD,OAAO;YACL,IAAI;YACJ,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,SAAS,EAAE,MAAM,CAAC,cAAc,CAAC,SAAS;YAC1C,GAAG,OAAO;SACX,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;QACjD,OAAO;YACL,IAAI;YACJ,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,WAAW,EAAE,MAAM,CAAC,cAAc,CAAC,WAAW;YAC9C,cAAc,EAAE,MAAM,CAAC,cAAc,CAAC,cAAc;YACpD,GAAG,OAAO;SACX,CAAC;IACJ,CAAC;IACD,OAAO;QACL,IAAI;QACJ,OAAO,EAAE,MAAM,CAAC,OAAO;QACvB,EAAE,EAAE,MAAM,CAAC,cAAc,CAAC,EAAE;QAC5B,GAAG,OAAO;KACX,CAAC;AACJ,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAmB,EACnB,YAAY,GAAsB,EAAE;IAEpC,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,cAAc,EAAE,CAAC;QAClD,OAAO;YACL,IAAI,EAAE,cAAc;YACpB,SAAS,EAAE,MAAM,CAAC,cAAc,CAAC,SAAS;YAC1C,OAAO,EAAE,MAAM,CAAC,OAAO;SACxB,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;QACjD,OAAO;YACL,IAAI,EAAE,aAAa;YACnB,WAAW,EAAE,MAAM,CAAC,cAAc,CAAC,WAAW;YAC9C,cAAc,EAAE,MAAM,CAAC,cAAc,CAAC,cAAc;YACpD,OAAO,EAAE,MAAM,CAAC,OAAO;SACxB,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,QAAQ,KAAK,SAAS;QAAE,OAAO,MAAM,CAAC,QAAQ,CAAC;IAC1D,MAAM,SAAS,GAAG,cAAc,CAAC,MAAM,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,SAAS,CAAC,EAAE,IAAI,SAAS,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;IAC/D,OAAO;QACL,IAAI,EAAE,MAAM;QACZ,EAAE;QACF,OAAO,EAAE,MAAM,CAAC,OAAO;QACvB,IAAI,EAAE,eAAe,CAAC,EAAE,GAAG,SAAS,EAAE,EAAE,EAAE,EAAE,YAAY,CAAC;KAC1D,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,gBAAgB,CAAC,KAA4B;IAC3D,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,cAAc,EAAE,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,iCAAiC,CAAC,KAAK,EAAE,EAAE,OAAO,EAAE,CAAC,CAAC;IACpE,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC;IAC5B,MAAM,OAAO,GAAkB,EAAE,CAAC;IAClC,MAAM,UAAU,GAAkB,EAAE,CAAC;IACrC,KAAK,MAAM,MAAM,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC;QACzC,MAAM,QAAQ,GAAG,gBAAgB,CAC/B,EAAE,IAAI,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,EAClE,WAAW,CAAC,MAAM,CAAC,IAAI,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CACpD,CAAC;QACF,IAAI,QAAQ,KAAK,SAAS;YAAE,SAAS;QACrC,MAAM,cAAc,GAAG,UAAU,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;QACpF,MAAM,OAAO,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;QACpD,MAAM,MAAM,GAAgB;YAC1B,OAAO;YACP,cAAc;YACd,OAAO,EAAE,QAAQ,CAAC,OAAO;YACzB,QAAQ,EAAE,SAAS;YACnB,uEAAuE;YACvE,uEAAuE;YACvE,sEAAsE;YACtE,8DAA8D;YAC9D,mBAAmB,EAAE,SAAS;SAC/B,CAAC;QACF,IAAI,cAAc,CAAC,IAAI,KAAK,cAAc,IAAI,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;YACpF,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YACrB,SAAS;QACX,CAAC;QACD,MAAM,QAAQ,GAAG,kBAAkB,CAAC,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QACjE,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;YAC3B,UAAU,CAAC,IAAI,CAAC,EAAE,GAAG,MAAM,EAAE,QAAQ,EAAE,CAAC,CAAC;YACzC,SAAS;QACX,CAAC;QACD,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACvB,CAAC;IACD,MAAM,iBAAiB,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CACtD,mBAAmB,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC,OAAO,CAAC,CAC1C,CAAC;IACF,MAAM,aAAa,GAA8B,KAAK,CAAC,gBAAgB;SACpE,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QACb,IAAI,EAAE,6BAAsC;QAC5C,OAAO,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,CAAC,MAAM,CAAC,EAAE,CAAqB;QAC1D,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;SACF,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;IAC7D,OAAO;QACL,GAAG,iBAAiB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAC3C,iBAAiB,CAAC,MAAM,EAAE,WAAW,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CACvD;QACD,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC/D,GAAG,aAAa;KACjB,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Internal bridge-discovery pipeline.\n *\n * Seeds, then seed-filtered symbolic enumeration, the Buckingham filter,\n * the structural classifier, the regime join gate, the chain order,\n * and the proof-target stub.\n * A confirmation reports the catalog id. A restatement reports the\n * pre-declared canonical equation. A regime mismatch is recorded and\n * is not a stub. Any other survivor is a statement skeleton. Nothing\n * is written into the catalog.\n *\n * No command calls this. `discovery.ts` and the probe are not inputs.\n *\n * @module atlas/chain-pipeline\n */\n\nimport type { ExprNode } from '../dimensional/validator.js';\nimport type { Dimension } from '../dimensional/types.js';\nimport { DIMENSIONLESS } from '../dimensional/types.js';\nimport { equals } from '../dimensional/algebra.js';\nimport { CONSTANTS } from '../dimensional/symbolic-constants.js';\nimport type { BridgeEdge } from '../composition/edge.js';\nimport { enumerateCompositionsWithRefusals } from '../composition/enumerate.js';\nimport { buckinghamFilter } from '../composition/buckingham-filter.js';\nimport { matchChain } from '../composition/chain-match.js';\nimport { compareChainEdgeIds, type ChainCandidate } from '../composition/chain-candidate.js';\nimport {\n  chainOrderKey,\n  orderChainRecords,\n  type ChainRecord,\n} from '../composition/chain-result.js';\nimport { joinRegimeMismatch, type ChainRegimeMismatch } from '../composition/chain-regime.js';\nimport { bridgeSeedKeys, physjsTheorem } from './physjs-ref.js';\nimport { emitProofTarget } from './proof-target.js';\nimport { scalarSymbolsFromMathTs } from '../composition/mathts-scalar-symbols.js';\n\n/** A chain whose normal form matches a catalog equation. Nothing was written. @internal */\nexport interface ChainConfirmationRecord {\n  readonly kind: 'confirmation';\n  readonly catalogId: number;\n  readonly edgeIds: readonly string[];\n}\n\n/** A chain the registry pre-declared with `restatesBridge`. Not a new equation. @internal */\nexport interface ChainRestatementRecord {\n  readonly kind: 'restatement';\n  readonly canonicalId: string;\n  readonly restatesBridge: string;\n  readonly edgeIds: readonly string[];\n}\n\n/** A provisional chain, as a Lean statement skeleton. @internal */\nexport interface ChainStubRecord {\n  readonly kind: 'stub';\n  readonly id: string;\n  readonly edgeIds: readonly string[];\n  readonly text: string;\n}\n\n/**\n * A pair the composition table refused.\n *\n * The message is `UndefinedCompositionError`. A dimension mismatch is not\n * this record.\n *\n * @internal\n */\nexport interface ChainCompositionRefusal {\n  readonly kind: 'rejected: composition table';\n  readonly edgeIds: readonly [string, string];\n  readonly message: string;\n}\n\n/** One ordered survivor, a regime rejection, or a composition-table refusal. @internal */\nexport type ChainPipelineResult =\n  | ChainConfirmationRecord\n  | ChainRestatementRecord\n  | ChainStubRecord\n  | ChainRegimeMismatch\n  | ChainCompositionRefusal;\n\nfunction isNumericName(name: string): boolean {\n  return /^\\d+(\\.\\d+)?$/.test(name);\n}\n\nfunction isDimensionlessConstant(name: string): boolean {\n  if (/^\\d*pi$/.test(name)) return true;\n  const registered = CONSTANTS[name];\n  return registered !== undefined && equals(registered.dim, DIMENSIONLESS);\n}\n\n/**\n * Dimensioned leaves of the composed formula, excluding the target name.\n * The names are the MathTS symbol filter. A dimensionless constant is not\n * a governing symbol. The elementary charge is, because its dimension is\n * charge.\n *\n * @internal\n */\nexport function governingOf(expr: ExprNode, targetName: string): { name: string; dim: Dimension }[] {\n  const vars: { name: string; dim: Dimension }[] = [];\n  for (const leaf of scalarSymbolsFromMathTs(expr)) {\n    if (leaf.name === targetName) continue;\n    if (isNumericName(leaf.name) || isDimensionlessConstant(leaf.name)) continue;\n    vars.push(leaf);\n  }\n  vars.sort((a, b) => (a.name < b.name ? -1 : a.name > b.name ? 1 : 0));\n  return vars;\n}\n\nfunction theoremsFor(edgeIds: readonly string[]): string[] {\n  return edgeIds.map((id) => {\n    const name = physjsTheorem(id);\n    if (name === undefined) {\n      throw new Error(`chain pipeline: manifest copy has no theorem for seed '${id}'`);\n    }\n    return name;\n  });\n}\n\nfunction proofCandidate(record: ChainRecord): ChainCandidate {\n  const kind = chainOrderKey(record);\n  const theorem = record.theorem ?? undefined;\n  const carried = theorem === undefined ? {} : { theorem };\n  if (record.classification.kind === 'confirmation') {\n    return {\n      kind,\n      edgeIds: record.edgeIds,\n      catalogId: record.classification.catalogId,\n      ...carried,\n    };\n  }\n  if (record.classification.kind === 'restatement') {\n    return {\n      kind,\n      edgeIds: record.edgeIds,\n      canonicalId: record.classification.canonicalId,\n      restatesBridge: record.classification.restatesBridge,\n      ...carried,\n    };\n  }\n  return {\n    kind,\n    edgeIds: record.edgeIds,\n    id: record.classification.id,\n    ...carried,\n  };\n}\n\n/**\n * The orchestrator view of one {@link ChainRecord}.\n *\n * A confirmation and a restatement ignore `mismatch`. A provisional\n * record with a mismatch is that rejection. Any other provisional record\n * is a stub whose id and theorem come from the same record the order key\n * reads.\n *\n * @internal\n */\nexport function renderChainRecord(\n  record: ChainRecord,\n  seedTheorems: readonly string[] = [],\n): ChainPipelineResult {\n  if (record.classification.kind === 'confirmation') {\n    return {\n      kind: 'confirmation',\n      catalogId: record.classification.catalogId,\n      edgeIds: record.edgeIds,\n    };\n  }\n  if (record.classification.kind === 'restatement') {\n    return {\n      kind: 'restatement',\n      canonicalId: record.classification.canonicalId,\n      restatesBridge: record.classification.restatesBridge,\n      edgeIds: record.edgeIds,\n    };\n  }\n  if (record.mismatch !== undefined) return record.mismatch;\n  const candidate = proofCandidate(record);\n  const id = candidate.id ?? `chain-${record.edgeIds.join('-')}`;\n  return {\n    kind: 'stub',\n    id,\n    edgeIds: record.edgeIds,\n    text: emitProofTarget({ ...candidate, id }, seedTheorems),\n  };\n}\n\n/**\n * Run the pipeline on `edges`.\n *\n * A manifest key whose derived kind is not `bridge` is not a premise.\n * A chain the Buckingham filter drops is absent. A confirmation or a\n * restatement is not sent to the regime gate. Any other mismatch is\n * `rejected: regime mismatch` and is not a stub. A composition-table\n * refusal is `rejected: composition table` and is not a dimension\n * failure. The catalog array is not modified.\n *\n * @internal\n */\nexport function runChainPipeline(edges: readonly BridgeEdge[]): readonly ChainPipelineResult[] {\n  const seedIds = new Set(bridgeSeedKeys());\n  const found = enumerateCompositionsWithRefusals(edges, { seedIds });\n  const report = found.report;\n  const records: ChainRecord[] = [];\n  const rejections: ChainRecord[] = [];\n  for (const target of report.proofTargets) {\n    const filtered = buckinghamFilter(\n      { name: target.second.target.name, dim: target.second.target.dim },\n      governingOf(target.expr, target.second.target.name),\n    );\n    if (filtered === undefined) continue;\n    const classification = matchChain(target.expr, [target.first.id, target.second.id]);\n    const edgeIds = [target.first.id, target.second.id];\n    const record: ChainRecord = {\n      edgeIds,\n      classification,\n      theorem: filtered.theorem,\n      mismatch: undefined,\n      // A quantity edge stores a quantity name, not a category object id, so\n      // no step can state a category claim for the pair. Unset, and the pair\n      // is not dropped. The sentence that a function computed this value is\n      // the record from before it was found to return nothing else.\n      categoryComposition: undefined,\n    };\n    if (classification.kind === 'confirmation' || classification.kind === 'restatement') {\n      records.push(record);\n      continue;\n    }\n    const mismatch = joinRegimeMismatch(target.first, target.second);\n    if (mismatch !== undefined) {\n      rejections.push({ ...record, mismatch });\n      continue;\n    }\n    records.push(record);\n  }\n  const orderedRejections = [...rejections].sort((a, b) =>\n    compareChainEdgeIds(a.edgeIds, b.edgeIds),\n  );\n  const tableRefusals: ChainCompositionRefusal[] = found.relationRefusals\n    .map((row) => ({\n      kind: 'rejected: composition table' as const,\n      edgeIds: [row.first.id, row.second.id] as [string, string],\n      message: row.message,\n    }))\n    .sort((a, b) => compareChainEdgeIds(a.edgeIds, b.edgeIds));\n  return [\n    ...orderChainRecords(records).map((record) =>\n      renderChainRecord(record, theoremsFor(record.edgeIds)),\n    ),\n    ...orderedRejections.map((record) => renderChainRecord(record)),\n    ...tableRefusals,\n  ];\n}\n"]}
+{"version":3,"file":"chain-pipeline.js","sourceRoot":"","sources":["../../src/atlas/chain-pipeline.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;GAcG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,oBAAoB,CAAC;AAG3C,OAAO,EAAE,aAAa,EAAE,MAAM,yBAAyB,CAAC;AACxD,OAAO,EAAE,MAAM,EAAE,MAAM,2BAA2B,CAAC;AACnD,OAAO,EAAE,SAAS,EAAE,oBAAoB,EAAE,MAAM,sCAAsC,CAAC;AACvF,OAAO,EAAE,gBAAgB,EAAE,MAAM,8BAA8B,CAAC;AAEhE,OAAO,EAAE,iCAAiC,EAAE,MAAM,6BAA6B,CAAC;AAChF,OAAO,EAAE,gBAAgB,EAAE,MAAM,qCAAqC,CAAC;AACvE,OAAO,EAAE,UAAU,EAAE,MAAM,+BAA+B,CAAC;AAC3D,OAAO,EAAE,mBAAmB,EAAuB,MAAM,mCAAmC,CAAC;AAC7F,OAAO,EACL,aAAa,EACb,iBAAiB,GAElB,MAAM,gCAAgC,CAAC;AACxC,OAAO,EAAE,kBAAkB,EAA4B,MAAM,gCAAgC,CAAC;AAC9F,OAAO,EAAE,cAAc,EAAE,aAAa,EAAE,MAAM,iBAAiB,CAAC;AAChE,OAAO,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AACpD,OAAO,EAAE,uBAAuB,EAAE,MAAM,yCAAyC,CAAC;AA+ClF,gGAAgG;AAChG,SAAS,uBAAuB,CAAC,IAAY;IAC3C,IAAI,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,IAAI,CAAC;IACtC,OAAO,oBAAoB,CAAC,IAAI,CAAC,IAAI,MAAM,CAAC,SAAS,CAAC,IAAI,CAAE,CAAC,GAAG,EAAE,aAAa,CAAC,CAAC;AACnF,CAAC;AAED;;;;;;;GAOG;AACH,MAAM,UAAU,WAAW,CAAC,IAAc,EAAE,UAAkB;IAC5D,MAAM,IAAI,GAAuC,EAAE,CAAC;IACpD,KAAK,MAAM,IAAI,IAAI,uBAAuB,CAAC,IAAI,CAAC,EAAE,CAAC;QACjD,IAAI,IAAI,CAAC,IAAI,KAAK,UAAU;YAAE,SAAS;QACvC,IAAI,gBAAgB,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,uBAAuB,CAAC,IAAI,CAAC,IAAI,CAAC;YAAE,SAAS;QAChF,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IAClB,CAAC;IACD,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC;IAChC,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,WAAW,CAAC,OAA0B;IAC7C,OAAO,OAAO,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE;QACxB,MAAM,IAAI,GAAG,aAAa,CAAC,EAAE,CAAC,CAAC;QAC/B,IAAI,IAAI,KAAK,SAAS,EAAE,CAAC;YACvB,MAAM,IAAI,KAAK,CAAC,0DAA0D,EAAE,GAAG,CAAC,CAAC;QACnF,CAAC;QACD,OAAO,IAAI,CAAC;IACd,CAAC,CAAC,CAAC;AACL,CAAC;AAED,SAAS,cAAc,CAAC,MAAmB;IACzC,MAAM,IAAI,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC;IACnC,MAAM,OAAO,GAAG,MAAM,CAAC,OAAO,IAAI,SAAS,CAAC;IAC5C,MAAM,OAAO,GAAG,OAAO,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC;IACzD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,cAAc,EAAE,CAAC;QAClD,OAAO;YACL,IAAI;YACJ,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,SAAS,EAAE,MAAM,CAAC,cAAc,CAAC,SAAS;YAC1C,GAAG,OAAO;SACX,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;QACjD,OAAO;YACL,IAAI;YACJ,OAAO,EAAE,MAAM,CAAC,OAAO;YACvB,WAAW,EAAE,MAAM,CAAC,cAAc,CAAC,WAAW;YAC9C,cAAc,EAAE,MAAM,CAAC,cAAc,CAAC,cAAc;YACpD,GAAG,OAAO;SACX,CAAC;IACJ,CAAC;IACD,OAAO;QACL,IAAI;QACJ,OAAO,EAAE,MAAM,CAAC,OAAO;QACvB,EAAE,EAAE,MAAM,CAAC,cAAc,CAAC,EAAE;QAC5B,GAAG,OAAO;KACX,CAAC;AACJ,CAAC;AAED;;;;;;;;;GASG;AACH,MAAM,UAAU,iBAAiB,CAC/B,MAAmB,EACnB,YAAY,GAAsB,EAAE;IAEpC,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,cAAc,EAAE,CAAC;QAClD,OAAO;YACL,IAAI,EAAE,cAAc;YACpB,SAAS,EAAE,MAAM,CAAC,cAAc,CAAC,SAAS;YAC1C,OAAO,EAAE,MAAM,CAAC,OAAO;SACxB,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;QACjD,OAAO;YACL,IAAI,EAAE,aAAa;YACnB,WAAW,EAAE,MAAM,CAAC,cAAc,CAAC,WAAW;YAC9C,cAAc,EAAE,MAAM,CAAC,cAAc,CAAC,cAAc;YACpD,OAAO,EAAE,MAAM,CAAC,OAAO;SACxB,CAAC;IACJ,CAAC;IACD,IAAI,MAAM,CAAC,QAAQ,KAAK,SAAS;QAAE,OAAO,MAAM,CAAC,QAAQ,CAAC;IAC1D,MAAM,SAAS,GAAG,cAAc,CAAC,MAAM,CAAC,CAAC;IACzC,MAAM,EAAE,GAAG,SAAS,CAAC,EAAE,IAAI,SAAS,MAAM,CAAC,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;IAC/D,OAAO;QACL,IAAI,EAAE,MAAM;QACZ,EAAE;QACF,OAAO,EAAE,MAAM,CAAC,OAAO;QACvB,IAAI,EAAE,eAAe,CAAC,EAAE,GAAG,SAAS,EAAE,EAAE,EAAE,EAAE,YAAY,CAAC;KAC1D,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;GAWG;AACH,MAAM,UAAU,gBAAgB,CAAC,KAA4B;IAC3D,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,cAAc,EAAE,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,iCAAiC,CAAC,KAAK,EAAE,EAAE,OAAO,EAAE,CAAC,CAAC;IACpE,MAAM,MAAM,GAAG,KAAK,CAAC,MAAM,CAAC;IAC5B,MAAM,OAAO,GAAkB,EAAE,CAAC;IAClC,MAAM,UAAU,GAAkB,EAAE,CAAC;IACrC,KAAK,MAAM,MAAM,IAAI,MAAM,CAAC,YAAY,EAAE,CAAC;QACzC,MAAM,QAAQ,GAAG,gBAAgB,CAC/B,EAAE,IAAI,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,GAAG,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,EAAE,EAClE,WAAW,CAAC,MAAM,CAAC,IAAI,EAAE,MAAM,CAAC,MAAM,CAAC,MAAM,CAAC,IAAI,CAAC,CACpD,CAAC;QACF,IAAI,QAAQ,KAAK,SAAS;YAAE,SAAS;QACrC,MAAM,cAAc,GAAG,UAAU,CAAC,MAAM,CAAC,IAAI,EAAE,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;QACpF,MAAM,OAAO,GAAG,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,EAAE,MAAM,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC;QACpD,MAAM,MAAM,GAAgB;YAC1B,OAAO;YACP,cAAc;YACd,OAAO,EAAE,QAAQ,CAAC,OAAO;YACzB,QAAQ,EAAE,SAAS;YACnB,uEAAuE;YACvE,uEAAuE;YACvE,sEAAsE;YACtE,8DAA8D;YAC9D,mBAAmB,EAAE,SAAS;SAC/B,CAAC;QACF,IAAI,cAAc,CAAC,IAAI,KAAK,cAAc,IAAI,cAAc,CAAC,IAAI,KAAK,aAAa,EAAE,CAAC;YACpF,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;YACrB,SAAS;QACX,CAAC;QACD,MAAM,QAAQ,GAAG,kBAAkB,CAAC,MAAM,CAAC,KAAK,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC;QACjE,IAAI,QAAQ,KAAK,SAAS,EAAE,CAAC;YAC3B,UAAU,CAAC,IAAI,CAAC,EAAE,GAAG,MAAM,EAAE,QAAQ,EAAE,CAAC,CAAC;YACzC,SAAS;QACX,CAAC;QACD,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACvB,CAAC;IACD,MAAM,iBAAiB,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CACtD,mBAAmB,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC,OAAO,CAAC,CAC1C,CAAC;IACF,MAAM,aAAa,GAA8B,KAAK,CAAC,gBAAgB;SACpE,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QACb,IAAI,EAAE,6BAAsC;QAC5C,OAAO,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,GAAG,CAAC,MAAM,CAAC,EAAE,CAAqB;QAC1D,OAAO,EAAE,GAAG,CAAC,OAAO;KACrB,CAAC,CAAC;SACF,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,mBAAmB,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC;IAC7D,OAAO;QACL,GAAG,iBAAiB,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAC3C,iBAAiB,CAAC,MAAM,EAAE,WAAW,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC,CACvD;QACD,GAAG,iBAAiB,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,EAAE,CAAC,iBAAiB,CAAC,MAAM,CAAC,CAAC;QAC/D,GAAG,aAAa;KACjB,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Internal bridge-discovery pipeline.\n *\n * Seeds, then seed-filtered symbolic enumeration, the Buckingham filter,\n * the structural classifier, the regime join gate, the chain order,\n * and the proof-target stub.\n * A confirmation reports the catalog id. A restatement reports the\n * pre-declared canonical equation. A regime mismatch is recorded and\n * is not a stub. Any other survivor is a statement skeleton. Nothing\n * is written into the catalog.\n *\n * No command calls this. `discovery.ts` and the probe are not inputs.\n *\n * @module atlas/chain-pipeline\n */\n\nimport { byKey } from '../core/compare.js';\nimport type { ExprNode } from '../dimensional/validator.js';\nimport type { Dimension } from '../dimensional/types.js';\nimport { DIMENSIONLESS } from '../dimensional/types.js';\nimport { equals } from '../dimensional/algebra.js';\nimport { CONSTANTS, isRegisteredConstant } from '../dimensional/symbolic-constants.js';\nimport { isDecimalLiteral } from '../dimensional/ast-leaves.js';\nimport type { BridgeEdge } from '../composition/edge.js';\nimport { enumerateCompositionsWithRefusals } from '../composition/enumerate.js';\nimport { buckinghamFilter } from '../composition/buckingham-filter.js';\nimport { matchChain } from '../composition/chain-match.js';\nimport { compareChainEdgeIds, type ChainCandidate } from '../composition/chain-candidate.js';\nimport {\n  chainOrderKey,\n  orderChainRecords,\n  type ChainRecord,\n} from '../composition/chain-result.js';\nimport { joinRegimeMismatch, type ChainRegimeMismatch } from '../composition/chain-regime.js';\nimport { bridgeSeedKeys, physjsTheorem } from './physjs-ref.js';\nimport { emitProofTarget } from './proof-target.js';\nimport { scalarSymbolsFromMathTs } from '../composition/mathts-scalar-symbols.js';\n\n/** A chain whose normal form matches a catalog equation. Nothing was written. @internal */\nexport interface ChainConfirmationRecord {\n  readonly kind: 'confirmation';\n  readonly catalogId: number;\n  readonly edgeIds: readonly string[];\n}\n\n/** A chain the registry pre-declared with `restatesBridge`. Not a new equation. @internal */\nexport interface ChainRestatementRecord {\n  readonly kind: 'restatement';\n  readonly canonicalId: string;\n  readonly restatesBridge: string;\n  readonly edgeIds: readonly string[];\n}\n\n/** A provisional chain, as a Lean statement skeleton. @internal */\nexport interface ChainStubRecord {\n  readonly kind: 'stub';\n  readonly id: string;\n  readonly edgeIds: readonly string[];\n  readonly text: string;\n}\n\n/**\n * A pair the composition table refused.\n *\n * The message is `UndefinedCompositionError`. A dimension mismatch is not\n * this record.\n *\n * @internal\n */\nexport interface ChainCompositionRefusal {\n  readonly kind: 'rejected: composition table';\n  readonly edgeIds: readonly [string, string];\n  readonly message: string;\n}\n\n/** One ordered survivor, a regime rejection, or a composition-table refusal. @internal */\nexport type ChainPipelineResult =\n  | ChainConfirmationRecord\n  | ChainRestatementRecord\n  | ChainStubRecord\n  | ChainRegimeMismatch\n  | ChainCompositionRefusal;\n\n/** A registered constant whose dimension is none, or a `pi` multiple (`2pi`). Own keys only. */\nfunction isDimensionlessConstant(name: string): boolean {\n  if (/^\\d*pi$/.test(name)) return true;\n  return isRegisteredConstant(name) && equals(CONSTANTS[name]!.dim, DIMENSIONLESS);\n}\n\n/**\n * Dimensioned leaves of the composed formula, excluding the target name.\n * The names are the MathTS symbol filter. A dimensionless constant is not\n * a governing symbol. The elementary charge is, because its dimension is\n * charge.\n *\n * @internal\n */\nexport function governingOf(expr: ExprNode, targetName: string): { name: string; dim: Dimension }[] {\n  const vars: { name: string; dim: Dimension }[] = [];\n  for (const leaf of scalarSymbolsFromMathTs(expr)) {\n    if (leaf.name === targetName) continue;\n    if (isDecimalLiteral(leaf.name) || isDimensionlessConstant(leaf.name)) continue;\n    vars.push(leaf);\n  }\n  vars.sort(byKey((v) => v.name));\n  return vars;\n}\n\nfunction theoremsFor(edgeIds: readonly string[]): string[] {\n  return edgeIds.map((id) => {\n    const name = physjsTheorem(id);\n    if (name === undefined) {\n      throw new Error(`chain pipeline: manifest copy has no theorem for seed '${id}'`);\n    }\n    return name;\n  });\n}\n\nfunction proofCandidate(record: ChainRecord): ChainCandidate {\n  const kind = chainOrderKey(record);\n  const theorem = record.theorem ?? undefined;\n  const carried = theorem === undefined ? {} : { theorem };\n  if (record.classification.kind === 'confirmation') {\n    return {\n      kind,\n      edgeIds: record.edgeIds,\n      catalogId: record.classification.catalogId,\n      ...carried,\n    };\n  }\n  if (record.classification.kind === 'restatement') {\n    return {\n      kind,\n      edgeIds: record.edgeIds,\n      canonicalId: record.classification.canonicalId,\n      restatesBridge: record.classification.restatesBridge,\n      ...carried,\n    };\n  }\n  return {\n    kind,\n    edgeIds: record.edgeIds,\n    id: record.classification.id,\n    ...carried,\n  };\n}\n\n/**\n * The orchestrator view of one {@link ChainRecord}.\n *\n * A confirmation and a restatement ignore `mismatch`. A provisional\n * record with a mismatch is that rejection. Any other provisional record\n * is a stub whose id and theorem come from the same record the order key\n * reads.\n *\n * @internal\n */\nexport function renderChainRecord(\n  record: ChainRecord,\n  seedTheorems: readonly string[] = [],\n): ChainPipelineResult {\n  if (record.classification.kind === 'confirmation') {\n    return {\n      kind: 'confirmation',\n      catalogId: record.classification.catalogId,\n      edgeIds: record.edgeIds,\n    };\n  }\n  if (record.classification.kind === 'restatement') {\n    return {\n      kind: 'restatement',\n      canonicalId: record.classification.canonicalId,\n      restatesBridge: record.classification.restatesBridge,\n      edgeIds: record.edgeIds,\n    };\n  }\n  if (record.mismatch !== undefined) return record.mismatch;\n  const candidate = proofCandidate(record);\n  const id = candidate.id ?? `chain-${record.edgeIds.join('-')}`;\n  return {\n    kind: 'stub',\n    id,\n    edgeIds: record.edgeIds,\n    text: emitProofTarget({ ...candidate, id }, seedTheorems),\n  };\n}\n\n/**\n * Run the pipeline on `edges`.\n *\n * A manifest key whose derived kind is not `bridge` is not a premise.\n * A chain the Buckingham filter drops is absent. A confirmation or a\n * restatement is not sent to the regime gate. Any other mismatch is\n * `rejected: regime mismatch` and is not a stub. A composition-table\n * refusal is `rejected: composition table` and is not a dimension\n * failure. The catalog array is not modified.\n *\n * @internal\n */\nexport function runChainPipeline(edges: readonly BridgeEdge[]): readonly ChainPipelineResult[] {\n  const seedIds = new Set(bridgeSeedKeys());\n  const found = enumerateCompositionsWithRefusals(edges, { seedIds });\n  const report = found.report;\n  const records: ChainRecord[] = [];\n  const rejections: ChainRecord[] = [];\n  for (const target of report.proofTargets) {\n    const filtered = buckinghamFilter(\n      { name: target.second.target.name, dim: target.second.target.dim },\n      governingOf(target.expr, target.second.target.name),\n    );\n    if (filtered === undefined) continue;\n    const classification = matchChain(target.expr, [target.first.id, target.second.id]);\n    const edgeIds = [target.first.id, target.second.id];\n    const record: ChainRecord = {\n      edgeIds,\n      classification,\n      theorem: filtered.theorem,\n      mismatch: undefined,\n      // A quantity edge stores a quantity name, not a category object id, so\n      // no step can state a category claim for the pair. Unset, and the pair\n      // is not dropped. The sentence that a function computed this value is\n      // the record from before it was found to return nothing else.\n      categoryComposition: undefined,\n    };\n    if (classification.kind === 'confirmation' || classification.kind === 'restatement') {\n      records.push(record);\n      continue;\n    }\n    const mismatch = joinRegimeMismatch(target.first, target.second);\n    if (mismatch !== undefined) {\n      rejections.push({ ...record, mismatch });\n      continue;\n    }\n    records.push(record);\n  }\n  const orderedRejections = [...rejections].sort((a, b) =>\n    compareChainEdgeIds(a.edgeIds, b.edgeIds),\n  );\n  const tableRefusals: ChainCompositionRefusal[] = found.relationRefusals\n    .map((row) => ({\n      kind: 'rejected: composition table' as const,\n      edgeIds: [row.first.id, row.second.id] as [string, string],\n      message: row.message,\n    }))\n    .sort((a, b) => compareChainEdgeIds(a.edgeIds, b.edgeIds));\n  return [\n    ...orderChainRecords(records).map((record) =>\n      renderChainRecord(record, theoremsFor(record.edgeIds)),\n    ),\n    ...orderedRejections.map((record) => renderChainRecord(record)),\n    ...tableRefusals,\n  ];\n}\n"]}

package/dist/atlas/diffusion/bridges-closure.d.ts

@@ -20,14 +20,10 @@
  * @module atlas/diffusion/bridges-closure
  */
 import type { AtlasBridge } from '../types.js';
-/** τ_p/t ceiling of the Langevin → Fick regime: WD4's fine resolution, t = 100 τ_p. @internal */
-export declare const LANGEVIN_MAX_TAU_RATIO = 0.01;
 /** ε = τDq² ceiling of the telegraph → Fick domain. @internal */
 export declare const TELEGRAPH_FICK_MAX_EPS = 0.05;
 /** ε = τDq² floor of the telegraph → wave domain. @internal */
 export declare const TELEGRAPH_WAVE_MIN_EPS = 25;
-/** Fourier-number floor of the steady-state restriction. @internal */
-export declare const STEADY_MIN_FOURIER = 1;
 /** Bridge: Langevin → Fick, the Einstein coarse-graining D = k_BT/γ. @internal */
 export declare const BRIDGE_LANGEVIN_DIFFUSION: AtlasBridge;
 /**

package/dist/atlas/diffusion/bridges-closure.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-closure.d.ts","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges-closure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;GAoBG;AAOH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAS/C,iGAAiG;AACjG,eAAO,MAAM,sBAAsB,OAAO,CAAC;AAE3C,iEAAiE;AACjE,eAAO,MAAM,sBAAsB,OAAO,CAAC;AAE3C,+DAA+D;AAC/D,eAAO,MAAM,sBAAsB,KAAK,CAAC;AAEzC,sEAAsE;AACtE,eAAO,MAAM,kBAAkB,IAAI,CAAC;AAwBpC,kFAAkF;AAClF,eAAO,MAAM,yBAAyB,EAAE,WAyDvC,CAAC;AAEF;;;;;;;GAOG;AACH,eAAO,MAAM,sBAAsB,EAAE,WA4DpC,CAAC;AAEF;;;;;;;;;GASG;AACH,eAAO,MAAM,0BAA0B,EAAE,WA0ExC,CAAC;AAEF;;;;GAIG;AACH,eAAO,MAAM,qBAAqB,EAAE,WAkDnC,CAAC;AAEF,sEAAsE;AACtE,eAAO,MAAM,mBAAmB,EAAE,WAqCjC,CAAC;AAEF,sDAAsD;AACtD,eAAO,MAAM,yBAAyB,EAAE,SAAS,WAAW,EAM3D,CAAC"}
+{"version":3,"file":"bridges-closure.d.ts","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges-closure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;GAoBG;AAOH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAY/C,iEAAiE;AACjE,eAAO,MAAM,sBAAsB,OAAO,CAAC;AAE3C,+DAA+D;AAC/D,eAAO,MAAM,sBAAsB,KAAK,CAAC;AAuBzC,kFAAkF;AAClF,eAAO,MAAM,yBAAyB,EAAE,WAyDvC,CAAC;AAEF;;;;;;;GAOG;AACH,eAAO,MAAM,sBAAsB,EAAE,WA4DpC,CAAC;AAEF;;;;;;;;;GASG;AACH,eAAO,MAAM,0BAA0B,EAAE,WA0ExC,CAAC;AAEF;;;;GAIG;AACH,eAAO,MAAM,qBAAqB,EAAE,WAkDnC,CAAC;AAEF,sEAAsE;AACtE,eAAO,MAAM,mBAAmB,EAAE,WAqCjC,CAAC;AAEF,sDAAsD;AACtD,eAAO,MAAM,yBAAyB,EAAE,SAAS,WAAW,EAM3D,CAAC"}

package/dist/atlas/diffusion/bridges-closure.js

@@ -22,29 +22,25 @@
 import { ENERGY, LENGTH, MASS, TIME, VELOCITY } from '../../dimensional/types.js';
 import { makeApproximation } from '../oscillators/bridges-limits.js';
 import { DAMPING } from '../oscillators/dimensions.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { physjsFormalRef } from '../physjs-ref.js';
 import { WAVENUMBER } from '../waves/models.js';
 import { DENSITY, DIFFUSIVITY, SPECIFIC_HEAT, THERMAL_CONDUCTIVITY } from './dimensions.js';
 import { DIFFUSION_FAMILY_NAME, VISCOSITY } from './models.js';
-import { telegraphSlowRateRatio, telegraphWaveFrequencyRatio } from './numerics.js';
+import { telegraphSlowRateExcess, telegraphWaveFrequencyRatio } from './numerics.js';
 const TEST = 'tests/atlas/closure.test.ts';
 const WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';
 /** τ_p/t ceiling of the Langevin → Fick regime: WD4's fine resolution, t = 100 τ_p. @internal */
-export const LANGEVIN_MAX_TAU_RATIO = 0.01;
+const LANGEVIN_MAX_TAU_RATIO = 0.01;
 /** ε = τDq² ceiling of the telegraph → Fick domain. @internal */
 export const TELEGRAPH_FICK_MAX_EPS = 0.05;
 /** ε = τDq² floor of the telegraph → wave domain. @internal */
 export const TELEGRAPH_WAVE_MIN_EPS = 25;
 /** Fourier-number floor of the steady-state restriction. @internal */
-export const STEADY_MIN_FOURIER = 1;
+const STEADY_MIN_FOURIER = 1;
 const regime = (parameters, inequalities) => regimeWithInputs(parameters, [], inequalities);
 /** A regime whose coordinates also include dimensionless model inputs (such as Re). */
-const regimeWithInputs = (parameters, dimensionlessInputs, inequalities) => ({
-    family: DIFFUSION_FAMILY_NAME,
-    inequalities,
-    groupDefinitions: deriveRegimeGroups(DIFFUSION_FAMILY_NAME, parameters, dimensionlessInputs),
-});
+const regimeWithInputs = (parameters, dimensionlessInputs, inequalities) => makeRegime(DIFFUSION_FAMILY_NAME, parameters, dimensionlessInputs, inequalities);
 /** Reynolds-number ceiling of the Stokes–Einstein derivation: a chosen threshold for Re ≪ 1. @internal */
 const STOKES_MAX_RE = 0.1;
 /** Ceiling of τ_p/t = m/(γt) for the overdamped limit: a chosen threshold for ≪ 1. @internal */
@@ -193,8 +189,8 @@ export const BRIDGE_TELEGRAPH_DIFFUSION = {
     sideConditions: ['ε = τDq² ≤ 0.05', 'after the initial layer: t ≫ τ'],
     bound: makeApproximation({
         K: 1,
-        delta: telegraphSlowRateRatio(TELEGRAPH_FICK_MAX_EPS, 1, 1) - 1,
-        deltaAt: (p) => telegraphSlowRateRatio(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN) - 1,
+        delta: telegraphSlowRateExcess(TELEGRAPH_FICK_MAX_EPS, 1, 1),
+        deltaAt: (p) => telegraphSlowRateExcess(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN),
         deltaAtBasis: 'closed-form',
         norm: 'relative error of the slow-mode decay rate of a Fourier mode, normalized by the value of the reduced model',
         domain: 'ε = τDq² ≤ 0.05',
@@ -204,7 +200,7 @@ export const BRIDGE_TELEGRAPH_DIFFUSION = {
             const { tau, D, q } = p;
             if (tau === undefined || D === undefined || q === undefined)
                 return false;
-            const delta = telegraphSlowRateRatio(tau, D, q) - 1;
+            const delta = telegraphSlowRateExcess(tau, D, q);
             if (!Number.isFinite(delta))
                 return false;
             if (t <= 5 * tau)

package/dist/atlas/diffusion/bridges-closure.js.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-closure.js","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges-closure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;GAoBG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AAClF,OAAO,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC;AACrE,OAAO,EAAE,OAAO,EAAE,MAAM,8BAA8B,CAAC;AACvD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClD,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AAEnD,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,WAAW,EAAE,aAAa,EAAE,oBAAoB,EAAE,MAAM,iBAAiB,CAAC;AAC5F,OAAO,EAAE,qBAAqB,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAC/D,OAAO,EAAE,sBAAsB,EAAE,2BAA2B,EAAE,MAAM,eAAe,CAAC;AAEpF,MAAM,IAAI,GAAG,6BAA6B,CAAC;AAC3C,MAAM,oBAAoB,GAAG,qCAAqC,CAAC;AAEnE,iGAAiG;AACjG,MAAM,CAAC,MAAM,sBAAsB,GAAG,IAAI,CAAC;AAE3C,iEAAiE;AACjE,MAAM,CAAC,MAAM,sBAAsB,GAAG,IAAI,CAAC;AAE3C,+DAA+D;AAC/D,MAAM,CAAC,MAAM,sBAAsB,GAAG,EAAE,CAAC;AAEzC,sEAAsE;AACtE,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,CAAC;AAEpC,MAAM,MAAM,GAAG,CACb,UAA4F,EAC5F,YAAmD,EAC5B,EAAE,CAAC,gBAAgB,CAAC,UAAU,EAAE,EAAE,EAAE,YAAY,CAAC,CAAC;AAE3E,uFAAuF;AACvF,MAAM,gBAAgB,GAAG,CACvB,UAA4F,EAC5F,mBAAsC,EACtC,YAAmD,EAC5B,EAAE,CAAC,CAAC;IAC3B,MAAM,EAAE,qBAAqB;IAC7B,YAAY;IACZ,gBAAgB,EAAE,kBAAkB,CAAC,qBAAqB,EAAE,UAAU,EAAE,mBAAmB,CAAC;CAC7F,CAAC,CAAC;AAEH,0GAA0G;AAC1G,MAAM,aAAa,GAAG,GAAG,CAAC;AAE1B,gGAAgG;AAChG,MAAM,oBAAoB,GAAG,IAAI,CAAC;AAElC,kFAAkF;AAClF,MAAM,CAAC,MAAM,yBAAyB,GAAgB;IACpD,EAAE,EAAE,uBAAuB;IAC3B,QAAQ,EAAE,iBAAiB;IAC3B,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,0EAA0E;IAC1F,SAAS,EAAE,CAAC,wCAAwC,EAAE,2CAA2C,CAAC;IAClG,eAAe,EAAE;QACf,kCAAkC;QAClC,yDAAyD;QACzD,0DAA0D;KAC3D;IACD,cAAc,EAAE,CAAC,yCAAyC,EAAE,qDAAqD,CAAC;IAClH,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE;QAC/B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;QAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;KACzB,EACD;QACE;YACE,KAAK,EAAE,qBAAqB;YAC5B,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,sBAAsB;YAC7B,KAAK,EAAE,iCAAiC;SACzC;KACF,CACF;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,wFAAwF;gBACxF,mFAAmF;YACrF,OAAO,EAAE,MAAM;SAChB;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,uFAAuF;SACnG;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,oCAAoC;SAChD;KACF;IACD,SAAS,EAAE;QACT,uFAAuF;QACvF,+DAA+D;KAChE;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,sBAAsB,GAAgB;IACjD,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,YAAY;IACtB,QAAQ,EAAE,CAAC,gBAAgB,EAAE,mBAAmB,CAAC;IACjD,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,wDAAwD;IACxE,SAAS,EAAE,CAAC,iDAAiD,CAAC;IAC9D,eAAe,EAAE,CAAC,kCAAkC,EAAE,+BAA+B,CAAC;IACtF,cAAc,EAAE;QACd,+FAA+F;QAC/F,kBAAkB;QAClB,kHAAkH;KACnH;IACD,+EAA+E;IAC/E,wEAAwE;IACxE,gFAAgF;IAChF,8EAA8E;IAC9E,MAAM,EAAE,gBAAgB,CACtB;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;QAC3B,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,SAAS,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;KAChC,EACD,CAAC,IAAI,CAAC,EACN;QACE,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,aAAa,EAAE,KAAK,EAAE,gCAAgC,EAAE;QACxF;YACE,KAAK,EAAE,qBAAqB;YAC5B,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,oBAAoB;YAC3B,KAAK,EAAE,wCAAwC;SAChD;KACF,CACF;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EACP,2GAA2G;gBAC3G,6DAA6D;SAChE;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,oFAAoF;KACrF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,0BAA0B,GAAgB;IACrD,EAAE,EAAE,wBAAwB;IAC5B,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,iBAAiB,CAAC;IAC7B,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,2EAA2E;IAC3F,SAAS,EAAE,CAAC,kDAAkD,EAAE,wBAAwB,CAAC;IACzF,eAAe,EAAE;QACf,gCAAgC;QAChC,wCAAwC;QACxC,4CAA4C;KAC7C;IACD,cAAc,EAAE,CAAC,iBAAiB,EAAE,gCAAgC,CAAC;IACrE,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,sBAAsB,CAAC,sBAAsB,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC;QAC/D,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,CACb,sBAAsB,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC;QAChG,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,4GAA4G;QAClH,MAAM,EAAE,iBAAiB;QACzB,OAAO,EACL,2FAA2F;YAC3F,oCAAoC;QACtC,YAAY,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACrB,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;YACxB,IAAI,GAAG,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;gBAAE,OAAO,KAAK,CAAC;YAC1E,MAAM,KAAK,GAAG,sBAAsB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YACpD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YAC1C,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG;gBAAE,OAAO,KAAK,CAAC;YAC/B,OAAO,KAAK,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACtD,CAAC;QACD,cAAc,EAAE,iBAAiB;QACjC,cAAc,EAAE,UAAU;QAC1B,yEAAyE;QACzE,yBAAyB;QACzB,UAAU,EAAE,CAAC,+DAA+D,CAAC;KAC9E,CAAC;IACF,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;KAC/B,EACD,CAAC,EAAE,KAAK,EAAE,eAAe,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,sBAAsB,EAAE,KAAK,EAAE,iBAAiB,EAAE,CAAC,CAChG;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,4FAA4F;gBAC5F,+BAA+B;YACjC,OAAO,EAAE,MAAM;SAChB;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,gFAAgF;SAC5F;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,4BAA4B;SACxC;KACF;IACD,SAAS,EAAE;QACT,iGAAiG;QACjG,0HAA0H;KAC3H;IACD,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,wBAAwB,CAAC;CACrD,CAAC;AAEF;;;;GAIG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAgB;IAChD,EAAE,EAAE,mBAAmB;IACvB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,iBAAiB,CAAC;IAC7B,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,iDAAiD;IACjE,SAAS,EAAE,CAAC,yBAAyB,EAAE,qCAAqC,CAAC;IAC7E,eAAe,EAAE,CAAC,+CAA+C,EAAE,uCAAuC,CAAC;IAC3G,cAAc,EAAE,CAAC,eAAe,EAAE,sDAAsD,CAAC;IACzF,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,CAAC,GAAG,2BAA2B,CAAC,sBAAsB,EAAE,CAAC,EAAE,CAAC,CAAC;QACpE,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,CACb,CAAC,GAAG,2BAA2B,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC;QACrG,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,6GAA6G;QACnH,MAAM,EAAE,eAAe;QACvB,OAAO,EACL,8FAA8F;QAChG,YAAY,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACrB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAClB,IAAI,GAAG,KAAK,SAAS,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;gBAAE,OAAO,KAAK,CAAC;YAClD,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAClD,CAAC;QACD,cAAc,EAAE,eAAe;QAC/B,cAAc,EAAE,SAAS;QACzB,6EAA6E;QAC7E,UAAU,EAAE,CAAC,8DAA8D,CAAC;KAC7E,CAAC;IACF,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;QAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;KAC7B,EACD,CAAC,EAAE,KAAK,EAAE,eAAe,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,sBAAsB,EAAE,KAAK,EAAE,eAAe,EAAE,CAAC,CAC9F;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,uEAAuE;SACnF;KACF;IACD,SAAS,EAAE,CAAC,0HAA0H,CAAC;IACvI,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,mBAAmB,CAAC;CAChD,CAAC;AAEF,sEAAsE;AACtE,MAAM,CAAC,MAAM,mBAAmB,GAAgB;IAC9C,EAAE,EAAE,iBAAiB;IACrB,QAAQ,EAAE,aAAa;IACvB,QAAQ,EAAE,CAAC,YAAY,CAAC;IACxB,UAAU,EAAE,kBAAkB;IAC9B,cAAc,EAAE,8DAA8D;IAC9E,SAAS,EAAE,CAAC,sBAAsB,EAAE,qCAAqC,CAAC;IAC1E,eAAe,EAAE,CAAC,eAAe,EAAE,+CAA+C,CAAC;IACnF,cAAc,EAAE,CAAC,wBAAwB,EAAE,+CAA+C,CAAC;IAC3F,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,oBAAoB,EAAE;QAC5C,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,OAAO,EAAE;QAC7B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,aAAa,EAAE;QAClC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;QAC5B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;KACzB,EACD;QACE;YACE,KAAK,EAAE,qCAAqC;YAC5C,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,kBAAkB;YACzB,KAAK,EAAE,gBAAgB;SACxB;KACF,CACF;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,qFAAqF;SACjG;KACF;IACD,SAAS,EAAE,CAAC,+EAA+E,CAAC;IAC5F,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,sDAAsD;AACtD,MAAM,CAAC,MAAM,yBAAyB,GAA2B;IAC/D,yBAAyB;IACzB,sBAAsB;IACtB,0BAA0B;IAC1B,qBAAqB;IACrB,mBAAmB;CACpB,CAAC","sourcesContent":["/**\n * Atlas Phase 4 — the five diffusion-family bridges added to close Sprint 4's\n * \"≥ 20 bridges\" exit criterion.\n *\n * Design note: `docs/planning/Atlas-Phase-4-Design.md` §4 (Sprint 4 closure).\n *\n * | Bridge | Relation | Witnesses |\n * |---|---|---|\n * | Langevin → Fick | `coarse-graining` | WD4 (convergence), WD4b (ballistic regime) |\n * | Langevin + Stokes drag → Fick | `derivation`, **hyperedge** | WD5 (canonical + water), WD5s (CAS) |\n * | telegraph → Fick | `approximation` (singular, τ → 0) | WD6 (convergence), WD6b (oscillatory modes) |\n * | telegraph → 1-D wave | `approximation` (τDq² → ∞, t ≪ τ) | WD7 (convergence) |\n * | heat → Laplace | `restriction` (steady state) | WD8 (convergence in time) |\n *\n * Every regime inequality keys a `PiGroup.formula` derived from the dimension\n * matrix. Every `delta` is the EXACT error at its domain edge; every error\n * function used for it is monotone on the domain, so the edge value is the\n * supremum.\n *\n * @module atlas/diffusion/bridges-closure\n */\n\nimport { ENERGY, LENGTH, MASS, TIME, VELOCITY } from '../../dimensional/types.js';\nimport { makeApproximation } from '../oscillators/bridges-limits.js';\nimport { DAMPING } from '../oscillators/dimensions.js';\nimport { deriveRegimeGroups } from '../regime.js';\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport type { AtlasBridge } from '../types.js';\nimport { WAVENUMBER } from '../waves/models.js';\nimport { DENSITY, DIFFUSIVITY, SPECIFIC_HEAT, THERMAL_CONDUCTIVITY } from './dimensions.js';\nimport { DIFFUSION_FAMILY_NAME, VISCOSITY } from './models.js';\nimport { telegraphSlowRateRatio, telegraphWaveFrequencyRatio } from './numerics.js';\n\nconst TEST = 'tests/atlas/closure.test.ts';\nconst WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';\n\n/** τ_p/t ceiling of the Langevin → Fick regime: WD4's fine resolution, t = 100 τ_p. @internal */\nexport const LANGEVIN_MAX_TAU_RATIO = 0.01;\n\n/** ε = τDq² ceiling of the telegraph → Fick domain. @internal */\nexport const TELEGRAPH_FICK_MAX_EPS = 0.05;\n\n/** ε = τDq² floor of the telegraph → wave domain. @internal */\nexport const TELEGRAPH_WAVE_MIN_EPS = 25;\n\n/** Fourier-number floor of the steady-state restriction. @internal */\nexport const STEADY_MIN_FOURIER = 1;\n\nconst regime = (\n  parameters: readonly { name: string; dim: import('../../dimensional/types.js').Dimension }[],\n  inequalities: AtlasBridge['regime']['inequalities'],\n): AtlasBridge['regime'] => regimeWithInputs(parameters, [], inequalities);\n\n/** A regime whose coordinates also include dimensionless model inputs (such as Re). */\nconst regimeWithInputs = (\n  parameters: readonly { name: string; dim: import('../../dimensional/types.js').Dimension }[],\n  dimensionlessInputs: readonly string[],\n  inequalities: AtlasBridge['regime']['inequalities'],\n): AtlasBridge['regime'] => ({\n  family: DIFFUSION_FAMILY_NAME,\n  inequalities,\n  groupDefinitions: deriveRegimeGroups(DIFFUSION_FAMILY_NAME, parameters, dimensionlessInputs),\n});\n\n/** Reynolds-number ceiling of the Stokes–Einstein derivation: a chosen threshold for Re ≪ 1. @internal */\nconst STOKES_MAX_RE = 0.1;\n\n/** Ceiling of τ_p/t = m/(γt) for the overdamped limit: a chosen threshold for ≪ 1. @internal */\nconst STOKES_MAX_TAU_RATIO = 0.01;\n\n/** Bridge: Langevin → Fick, the Einstein coarse-graining D = k_BT/γ. @internal */\nexport const BRIDGE_LANGEVIN_DIFFUSION: AtlasBridge = {\n  id: 'ab-langevin-diffusion',\n  relation: 'coarse-graining',\n  premises: ['model-langevin'],\n  conclusion: 'model-fick',\n  transformation: 'average over times t ≫ τ_p = m/γ: ⟨x²⟩ → 2Dt with D = k_B T/γ (Einstein)',\n  preserves: ['the long-time mean-square displacement', 'the equilibrium velocity variance k_B T/m'],\n  doesNotPreserve: [\n    'the velocity as a state variable',\n    'the ballistic regime t ≲ τ_p, where ⟨x²⟩ ≈ (k_B T/m) t²',\n    'the momentum relaxation time τ_p as an independent scale',\n  ],\n  sideConditions: ['overdamped observation times: τ_p/t ≪ 1', 'white thermal noise, fluctuation–dissipation 2γk_BT'],\n  regime: regime(\n    [\n      { name: 'm', dim: MASS },\n      { name: 'gamma', dim: DAMPING },\n      { name: 'kT', dim: ENERGY },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 't', dim: TIME },\n    ],\n    [\n      {\n        group: 'm · gamma^-1 · t^-1',\n        op: '<=',\n        bound: LANGEVIN_MAX_TAU_RATIO,\n        alias: 'τ_p/t ≪ 1 (machine form ≤ 0.01)',\n      },\n    ],\n  ),\n  counterexamples: [\n    {\n      description:\n        'At t = 0.1 τ_p the Langevin mean-square displacement is 0.048 of 2Dt: the particle is ' +\n        'still ballistic, ⟨x²⟩ ≈ (k_BT/m)t², and the diffusion model does not describe it.',\n      witness: 'WD4b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'WD4',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '⟨x²⟩/(2Dt) within 0.02 of 1 at t = 100 τ_p, from the Langevin moment equations by RK4',\n    },\n    {\n      id: 'WD4b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '⟨x²⟩/(2Dt) = 0.0484 at t = 0.1 τ_p',\n    },\n  ],\n  citations: [\n    'Uhlenbeck & Ornstein, Phys. Rev. 36 (1930) 823 - On the theory of the Brownian motion',\n    'Einstein, Ann. Phys. 17 (1905) 549 - the relation D = k_B T/γ',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge: Langevin + Stokes drag → Fick with the Stokes–Einstein coefficient.\n *\n * A HYPEREDGE: the Langevin model supplies D = k_BT/γ with γ free; Stokes drag\n * supplies γ = 6πηa for a sphere. Neither alone gives D in terms of the fluid.\n *\n * @internal\n */\nexport const BRIDGE_STOKES_EINSTEIN: AtlasBridge = {\n  id: 'ab-stokes-einstein',\n  relation: 'derivation',\n  premises: ['model-langevin', 'model-stokes-drag'],\n  conclusion: 'model-fick',\n  transformation: 'substitute γ = 6πηa into D = k_B T/γ: D = k_B T/(6πηa)',\n  preserves: ['the long-time diffusion coefficient of a sphere'],\n  doesNotPreserve: ['the particle mass (it drops out)', 'the shape beyond the radius a'],\n  sideConditions: [\n    'creeping flow around the sphere, Re ≪ 1 (machine form Re ≤ 0.1: a chosen threshold for \"≪ 1\")',\n    'no-slip boundary',\n    'overdamped times t ≫ m/γ (machine form m/(γt) ≤ 0.01: a chosen threshold for \"≪ 1\", as in ab-langevin-diffusion)',\n  ],\n  // The prose conditions were not machine-checked, so `upt regime` reported this\n  // bridge VACUOUS and \"valid\" at Re = 1000 (persona finding L8). Re is a\n  // dimensionless input of the flow; m/(γt) is listed first so its π-group key is\n  // the one ab-langevin-diffusion uses, and one --at value checks both bridges.\n  regime: regimeWithInputs(\n    [\n      { name: 'm', dim: MASS },\n      { name: 'gamma', dim: DAMPING },\n      { name: 't', dim: TIME },\n      { name: 'kT', dim: ENERGY },\n      { name: 'eta', dim: VISCOSITY },\n      { name: 'a', dim: LENGTH },\n      { name: 'D', dim: DIFFUSIVITY },\n    ],\n    ['Re'],\n    [\n      { group: 'Re', op: '<=', bound: STOKES_MAX_RE, alias: 'Re ≪ 1 (machine form Re ≤ 0.1)' },\n      {\n        group: 'm · gamma^-1 · t^-1',\n        op: '<=',\n        bound: STOKES_MAX_TAU_RATIO,\n        alias: 'τ_p/t ≪ 1 (machine form m/(γt) ≤ 0.01)',\n      },\n    ],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD5',\n      kind: 'numeric',\n      test: TEST,\n      tolerance:\n        'agrees with CE-stokes-einstein to 1e-12 relative; a sphere of radius a = 0.5 µm (1 µm diameter) in water ' +\n        'at T = 293.15 K, η = 1.0016e-3 Pa·s gives D = 4.29e-13 m²/s',\n    },\n    {\n      id: 'WD5s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: k_BT/γ under γ ↦ 6πηa minus k_BT/(6πηa) simplifies to literal 0',\n    },\n  ],\n  citations: [\n    'Einstein, Ann. Phys. 17 (1905) 549 - Brownian motion of a sphere, D = k_B T/(6πηa)',\n    'Stokes, Trans. Camb. Phil. Soc. 9 (1851) 8 - the drag on a sphere in creeping flow',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge: telegraph → Fick, the SINGULAR limit τ → 0.\n *\n * Singular because the highest time derivative is lost: the telegraph equation\n * needs two initial conditions and diffusion takes one, so a fast mode\n * `e^{−t/τ}` forms an initial layer. The bound therefore holds only AFTER that\n * layer, and the machine horizon has a LOWER edge as well as an upper one.\n *\n * @internal\n */\nexport const BRIDGE_TELEGRAPH_DIFFUSION: AtlasBridge = {\n  id: 'ab-telegraph-diffusion',\n  relation: 'approximation',\n  premises: ['model-telegraph'],\n  conclusion: 'model-fick',\n  transformation: 'drop τ u_tt: the slow mode decays at Dq²(1 + ε + 2ε² + …) → Dq², ε = τDq²',\n  preserves: ['the slow decay rate of each Fourier mode to O(ε)', 'the total amount ∫u dx'],\n  doesNotPreserve: [\n    'the finite signal speed √(D/τ)',\n    'the second initial condition u_t(x, 0)',\n    'oscillatory modes, which exist for ε > 1/4',\n  ],\n  sideConditions: ['ε = τDq² ≤ 0.05', 'after the initial layer: t ≫ τ'],\n  bound: makeApproximation({\n    K: 1,\n    delta: telegraphSlowRateRatio(TELEGRAPH_FICK_MAX_EPS, 1, 1) - 1,\n    deltaAt: (p) =>\n      telegraphSlowRateRatio(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN) - 1,\n    deltaAtBasis: 'closed-form',\n    norm: 'relative error of the slow-mode decay rate of a Fourier mode, normalized by the value of the reduced model',\n    domain: 'ε = τDq² ≤ 0.05',\n    horizon:\n      'τ ≪ t ≪ 1/(δ D q²): after the initial layer and before the decay-rate error accumulates; ' +\n      'machine form 5τ < t < 0.1/(δ D q²)',\n    horizonHolds: (t, p) => {\n      const { tau, D, q } = p;\n      if (tau === undefined || D === undefined || q === undefined) return false;\n      const delta = telegraphSlowRateRatio(tau, D, q) - 1;\n      if (!Number.isFinite(delta)) return false;\n      if (t <= 5 * tau) return false;\n      return delta === 0 || t < 0.1 / (delta * D * q * q);\n    },\n    parameterRange: 'ε = τDq² ≤ 0.05',\n    limitCharacter: 'singular',\n    // A slow-mode decay-rate error. The horizon says that error accumulates,\n    // so time is not listed.\n    uniformity: ['slow-mode decay rate of one Fourier mode, for ε = τDq² ≤ 0.05'],\n  }),\n  regime: regime(\n    [\n      { name: 'tau', dim: TIME },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 'q', dim: WAVENUMBER },\n    ],\n    [{ group: 'tau · D · q^2', op: '<=', bound: TELEGRAPH_FICK_MAX_EPS, alias: 'ε = τDq² ≤ 0.05' }],\n  ),\n  counterexamples: [\n    {\n      description:\n        'At ε = τDq² = 1 the telegraph mode is OSCILLATORY (complex rate, NaN on the slow branch): ' +\n        'no diffusion mode oscillates.',\n      witness: 'WD6b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'WD6',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'slow rate / Dq² within 0.03 of 1 at τ = 0.025 (D = q = 1); error halves with τ',\n    },\n    {\n      id: 'WD6b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'no real slow rate at ε = 1',\n    },\n  ],\n  citations: [\n    'Cattaneo, Atti Sem. Mat. Fis. Univ. Modena 3 (1948) 83 - heat conduction with a relaxation time',\n    'Goldstein, Q. J. Mech. Appl. Math. 4 (1951) 129 - On diffusion by discontinuous movements, and on the telegraph equation',\n  ],\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-telegraph-diffusion'),\n};\n\n/**\n * Bridge: telegraph → 1-D wave, for strongly underdamped modes over short times.\n *\n * @internal\n */\nexport const BRIDGE_TELEGRAPH_WAVE: AtlasBridge = {\n  id: 'ab-telegraph-wave',\n  relation: 'approximation',\n  premises: ['model-telegraph'],\n  conclusion: 'model-wave-1d',\n  transformation: 'drop u_t: τ u_tt = D u_xx, a wave with c² = D/τ',\n  preserves: ['the signal speed √(D/τ)', 'the oscillation frequency to O(1/ε)'],\n  doesNotPreserve: ['damping: telegraph modes decay as e^{−t/(2τ)}', 'relaxation to diffusion at long times'],\n  sideConditions: ['ε = τDq² ≥ 25', 'short times t ≪ τ (the bound is not uniform in time)'],\n  bound: makeApproximation({\n    K: 1,\n    delta: 1 - telegraphWaveFrequencyRatio(TELEGRAPH_WAVE_MIN_EPS, 1, 1),\n    deltaAt: (p) =>\n      1 - telegraphWaveFrequencyRatio(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN),\n    deltaAtBasis: 'closed-form',\n    norm: 'relative error of the oscillation frequency of a Fourier mode, normalized by the value of the reduced model',\n    domain: 'ε = τDq² ≥ 25',\n    horizon:\n      't ≪ τ: before damping e^{−t/(2τ)} removes 10% of the amplitude; machine form t < 2τ ln(10/9)',\n    horizonHolds: (t, p) => {\n      const { tau } = p;\n      if (tau === undefined || !(tau > 0)) return false;\n      return t >= 0 && t < 2 * tau * Math.log(10 / 9);\n    },\n    parameterRange: 'ε = τDq² ≥ 25',\n    limitCharacter: 'regular',\n    // An oscillation-frequency error. The record says it is not uniform in time.\n    uniformity: ['oscillation frequency of one Fourier mode, for ε = τDq² ≥ 25'],\n  }),\n  regime: regime(\n    [\n      { name: 'tau', dim: TIME },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 'q', dim: WAVENUMBER },\n      { name: 'c', dim: VELOCITY },\n    ],\n    [{ group: 'tau · D · q^2', op: '>=', bound: TELEGRAPH_WAVE_MIN_EPS, alias: 'ε = τDq² ≥ 25' }],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD7',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'frequency ratio within 3e-3 of 1 at ε = 50; error halves as ε doubles',\n    },\n  ],\n  citations: ['Goldstein, Q. J. Mech. Appl. Math. 4 (1951) 129 - On diffusion by discontinuous movements, and on the telegraph equation'],\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-telegraph-wave'),\n};\n\n/** Bridge: heat → Laplace, the steady-state restriction. @internal */\nexport const BRIDGE_HEAT_LAPLACE: AtlasBridge = {\n  id: 'ab-heat-laplace',\n  relation: 'restriction',\n  premises: ['model-heat'],\n  conclusion: 'model-laplace-1d',\n  transformation: '∂T/∂t = 0: κ T_xx = 0, T linear between the end temperatures',\n  preserves: ['the end temperatures', 'the steady heat flux κ(T_b − T_a)/ℓ'],\n  doesNotPreserve: ['the transient', 'ρ and c_p, which drop out of the steady state'],\n  sideConditions: ['fixed end temperatures', 'Fourier number αt/ℓ² ≥ 1 (transients decayed)'],\n  regime: regime(\n    [\n      { name: 'kappa', dim: THERMAL_CONDUCTIVITY },\n      { name: 'rho', dim: DENSITY },\n      { name: 'cp', dim: SPECIFIC_HEAT },\n      { name: 'ell', dim: LENGTH },\n      { name: 't', dim: TIME },\n    ],\n    [\n      {\n        group: 'kappa · rho^-1 · cp^-1 · ell^-2 · t',\n        op: '>=',\n        bound: STEADY_MIN_FOURIER,\n        alias: 'Fo = αt/ℓ² ≥ 1',\n      },\n    ],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD8',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'max deviation from the linear profile below 0.025 at αt/ℓ² = 0.4; shrinks with time',\n    },\n  ],\n  citations: ['Carslaw & Jaeger, Conduction of Heat in Solids - steady linear flow in a slab'],\n  reviewStatus: 'proposed',\n};\n\n/** The diffusion-family closure bridges. @internal */\nexport const DIFFUSION_CLOSURE_BRIDGES: readonly AtlasBridge[] = [\n  BRIDGE_LANGEVIN_DIFFUSION,\n  BRIDGE_STOKES_EINSTEIN,\n  BRIDGE_TELEGRAPH_DIFFUSION,\n  BRIDGE_TELEGRAPH_WAVE,\n  BRIDGE_HEAT_LAPLACE,\n];\n"]}
+{"version":3,"file":"bridges-closure.js","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges-closure.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;GAoBG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,IAAI,EAAE,IAAI,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AAClF,OAAO,EAAE,iBAAiB,EAAE,MAAM,kCAAkC,CAAC;AACrE,OAAO,EAAE,OAAO,EAAE,MAAM,8BAA8B,CAAC;AACvD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AAEnD,OAAO,EAAE,UAAU,EAAE,MAAM,oBAAoB,CAAC;AAChD,OAAO,EAAE,OAAO,EAAE,WAAW,EAAE,aAAa,EAAE,oBAAoB,EAAE,MAAM,iBAAiB,CAAC;AAC5F,OAAO,EAAE,qBAAqB,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAC/D,OAAO,EAAE,uBAAuB,EAAE,2BAA2B,EAAE,MAAM,eAAe,CAAC;AAErF,MAAM,IAAI,GAAG,6BAA6B,CAAC;AAC3C,MAAM,oBAAoB,GAAG,qCAAqC,CAAC;AAEnE,iGAAiG;AACjG,MAAM,sBAAsB,GAAG,IAAI,CAAC;AAEpC,iEAAiE;AACjE,MAAM,CAAC,MAAM,sBAAsB,GAAG,IAAI,CAAC;AAE3C,+DAA+D;AAC/D,MAAM,CAAC,MAAM,sBAAsB,GAAG,EAAE,CAAC;AAEzC,sEAAsE;AACtE,MAAM,kBAAkB,GAAG,CAAC,CAAC;AAE7B,MAAM,MAAM,GAAG,CACb,UAA4F,EAC5F,YAAmD,EAC5B,EAAE,CAAC,gBAAgB,CAAC,UAAU,EAAE,EAAE,EAAE,YAAY,CAAC,CAAC;AAE3E,uFAAuF;AACvF,MAAM,gBAAgB,GAAG,CACvB,UAA4F,EAC5F,mBAAsC,EACtC,YAAmD,EAC5B,EAAE,CAAC,UAAU,CAAC,qBAAqB,EAAE,UAAU,EAAE,mBAAmB,EAAE,YAAY,CAAC,CAAC;AAE7G,0GAA0G;AAC1G,MAAM,aAAa,GAAG,GAAG,CAAC;AAE1B,gGAAgG;AAChG,MAAM,oBAAoB,GAAG,IAAI,CAAC;AAElC,kFAAkF;AAClF,MAAM,CAAC,MAAM,yBAAyB,GAAgB;IACpD,EAAE,EAAE,uBAAuB;IAC3B,QAAQ,EAAE,iBAAiB;IAC3B,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,0EAA0E;IAC1F,SAAS,EAAE,CAAC,wCAAwC,EAAE,2CAA2C,CAAC;IAClG,eAAe,EAAE;QACf,kCAAkC;QAClC,yDAAyD;QACzD,0DAA0D;KAC3D;IACD,cAAc,EAAE,CAAC,yCAAyC,EAAE,qDAAqD,CAAC;IAClH,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE;QAC/B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;QAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;KACzB,EACD;QACE;YACE,KAAK,EAAE,qBAAqB;YAC5B,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,sBAAsB;YAC7B,KAAK,EAAE,iCAAiC;SACzC;KACF,CACF;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,wFAAwF;gBACxF,mFAAmF;YACrF,OAAO,EAAE,MAAM;SAChB;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,uFAAuF;SACnG;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,oCAAoC;SAChD;KACF;IACD,SAAS,EAAE;QACT,uFAAuF;QACvF,+DAA+D;KAChE;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;GAOG;AACH,MAAM,CAAC,MAAM,sBAAsB,GAAgB;IACjD,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,YAAY;IACtB,QAAQ,EAAE,CAAC,gBAAgB,EAAE,mBAAmB,CAAC;IACjD,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,wDAAwD;IACxE,SAAS,EAAE,CAAC,iDAAiD,CAAC;IAC9D,eAAe,EAAE,CAAC,kCAAkC,EAAE,+BAA+B,CAAC;IACtF,cAAc,EAAE;QACd,+FAA+F;QAC/F,kBAAkB;QAClB,kHAAkH;KACnH;IACD,+EAA+E;IAC/E,wEAAwE;IACxE,gFAAgF;IAChF,8EAA8E;IAC9E,MAAM,EAAE,gBAAgB,CACtB;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,OAAO,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,MAAM,EAAE;QAC3B,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,SAAS,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;KAChC,EACD,CAAC,IAAI,CAAC,EACN;QACE,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,aAAa,EAAE,KAAK,EAAE,gCAAgC,EAAE;QACxF;YACE,KAAK,EAAE,qBAAqB;YAC5B,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,oBAAoB;YAC3B,KAAK,EAAE,wCAAwC;SAChD;KACF,CACF;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EACP,2GAA2G;gBAC3G,6DAA6D;SAChE;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,oFAAoF;KACrF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,0BAA0B,GAAgB;IACrD,EAAE,EAAE,wBAAwB;IAC5B,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,iBAAiB,CAAC;IAC7B,UAAU,EAAE,YAAY;IACxB,cAAc,EAAE,2EAA2E;IAC3F,SAAS,EAAE,CAAC,kDAAkD,EAAE,wBAAwB,CAAC;IACzF,eAAe,EAAE;QACf,gCAAgC;QAChC,wCAAwC;QACxC,4CAA4C;KAC7C;IACD,cAAc,EAAE,CAAC,iBAAiB,EAAE,gCAAgC,CAAC;IACrE,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,uBAAuB,CAAC,sBAAsB,EAAE,CAAC,EAAE,CAAC,CAAC;QAC5D,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,CACb,uBAAuB,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC;QAC7F,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,4GAA4G;QAClH,MAAM,EAAE,iBAAiB;QACzB,OAAO,EACL,2FAA2F;YAC3F,oCAAoC;QACtC,YAAY,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACrB,MAAM,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;YACxB,IAAI,GAAG,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS,IAAI,CAAC,KAAK,SAAS;gBAAE,OAAO,KAAK,CAAC;YAC1E,MAAM,KAAK,GAAG,uBAAuB,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YACjD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YAC1C,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG;gBAAE,OAAO,KAAK,CAAC;YAC/B,OAAO,KAAK,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACtD,CAAC;QACD,cAAc,EAAE,iBAAiB;QACjC,cAAc,EAAE,UAAU;QAC1B,yEAAyE;QACzE,yBAAyB;QACzB,UAAU,EAAE,CAAC,+DAA+D,CAAC;KAC9E,CAAC;IACF,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;KAC/B,EACD,CAAC,EAAE,KAAK,EAAE,eAAe,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,sBAAsB,EAAE,KAAK,EAAE,iBAAiB,EAAE,CAAC,CAChG;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,4FAA4F;gBAC5F,+BAA+B;YACjC,OAAO,EAAE,MAAM;SAChB;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,gFAAgF;SAC5F;QACD;YACE,EAAE,EAAE,MAAM;YACV,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,4BAA4B;SACxC;KACF;IACD,SAAS,EAAE;QACT,iGAAiG;QACjG,0HAA0H;KAC3H;IACD,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,wBAAwB,CAAC;CACrD,CAAC;AAEF;;;;GAIG;AACH,MAAM,CAAC,MAAM,qBAAqB,GAAgB;IAChD,EAAE,EAAE,mBAAmB;IACvB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,iBAAiB,CAAC;IAC7B,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,iDAAiD;IACjE,SAAS,EAAE,CAAC,yBAAyB,EAAE,qCAAqC,CAAC;IAC7E,eAAe,EAAE,CAAC,+CAA+C,EAAE,uCAAuC,CAAC;IAC3G,cAAc,EAAE,CAAC,eAAe,EAAE,sDAAsD,CAAC;IACzF,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,KAAK,EAAE,CAAC,GAAG,2BAA2B,CAAC,sBAAsB,EAAE,CAAC,EAAE,CAAC,CAAC;QACpE,OAAO,EAAE,CAAC,CAAC,EAAE,EAAE,CACb,CAAC,GAAG,2BAA2B,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,GAAG,CAAC;QACrG,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,6GAA6G;QACnH,MAAM,EAAE,eAAe;QACvB,OAAO,EACL,8FAA8F;QAChG,YAAY,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;YACrB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YAClB,IAAI,GAAG,KAAK,SAAS,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;gBAAE,OAAO,KAAK,CAAC;YAClD,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAClD,CAAC;QACD,cAAc,EAAE,eAAe;QAC/B,cAAc,EAAE,SAAS;QACzB,6EAA6E;QAC7E,UAAU,EAAE,CAAC,8DAA8D,CAAC;KAC7E,CAAC;IACF,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,IAAI,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,WAAW,EAAE;QAC/B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,UAAU,EAAE;QAC9B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;KAC7B,EACD,CAAC,EAAE,KAAK,EAAE,eAAe,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,sBAAsB,EAAE,KAAK,EAAE,eAAe,EAAE,CAAC,CAC9F;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,uEAAuE;SACnF;KACF;IACD,SAAS,EAAE,CAAC,0HAA0H,CAAC;IACvI,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,mBAAmB,CAAC;CAChD,CAAC;AAEF,sEAAsE;AACtE,MAAM,CAAC,MAAM,mBAAmB,GAAgB;IAC9C,EAAE,EAAE,iBAAiB;IACrB,QAAQ,EAAE,aAAa;IACvB,QAAQ,EAAE,CAAC,YAAY,CAAC;IACxB,UAAU,EAAE,kBAAkB;IAC9B,cAAc,EAAE,8DAA8D;IAC9E,SAAS,EAAE,CAAC,sBAAsB,EAAE,qCAAqC,CAAC;IAC1E,eAAe,EAAE,CAAC,eAAe,EAAE,+CAA+C,CAAC;IACnF,cAAc,EAAE,CAAC,wBAAwB,EAAE,+CAA+C,CAAC;IAC3F,MAAM,EAAE,MAAM,CACZ;QACE,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,oBAAoB,EAAE;QAC5C,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,OAAO,EAAE;QAC7B,EAAE,IAAI,EAAE,IAAI,EAAE,GAAG,EAAE,aAAa,EAAE;QAClC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;QAC5B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;KACzB,EACD;QACE;YACE,KAAK,EAAE,qCAAqC;YAC5C,EAAE,EAAE,IAAI;YACR,KAAK,EAAE,kBAAkB;YACzB,KAAK,EAAE,gBAAgB;SACxB;KACF,CACF;IACD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,qFAAqF;SACjG;KACF;IACD,SAAS,EAAE,CAAC,+EAA+E,CAAC;IAC5F,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,sDAAsD;AACtD,MAAM,CAAC,MAAM,yBAAyB,GAA2B;IAC/D,yBAAyB;IACzB,sBAAsB;IACtB,0BAA0B;IAC1B,qBAAqB;IACrB,mBAAmB;CACpB,CAAC","sourcesContent":["/**\n * Atlas Phase 4 — the five diffusion-family bridges added to close Sprint 4's\n * \"≥ 20 bridges\" exit criterion.\n *\n * Design note: `docs/planning/Atlas-Phase-4-Design.md` §4 (Sprint 4 closure).\n *\n * | Bridge | Relation | Witnesses |\n * |---|---|---|\n * | Langevin → Fick | `coarse-graining` | WD4 (convergence), WD4b (ballistic regime) |\n * | Langevin + Stokes drag → Fick | `derivation`, **hyperedge** | WD5 (canonical + water), WD5s (CAS) |\n * | telegraph → Fick | `approximation` (singular, τ → 0) | WD6 (convergence), WD6b (oscillatory modes) |\n * | telegraph → 1-D wave | `approximation` (τDq² → ∞, t ≪ τ) | WD7 (convergence) |\n * | heat → Laplace | `restriction` (steady state) | WD8 (convergence in time) |\n *\n * Every regime inequality keys a `PiGroup.formula` derived from the dimension\n * matrix. Every `delta` is the EXACT error at its domain edge; every error\n * function used for it is monotone on the domain, so the edge value is the\n * supremum.\n *\n * @module atlas/diffusion/bridges-closure\n */\n\nimport { ENERGY, LENGTH, MASS, TIME, VELOCITY } from '../../dimensional/types.js';\nimport { makeApproximation } from '../oscillators/bridges-limits.js';\nimport { DAMPING } from '../oscillators/dimensions.js';\nimport { makeRegime } from '../regime.js';\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport type { AtlasBridge } from '../types.js';\nimport { WAVENUMBER } from '../waves/models.js';\nimport { DENSITY, DIFFUSIVITY, SPECIFIC_HEAT, THERMAL_CONDUCTIVITY } from './dimensions.js';\nimport { DIFFUSION_FAMILY_NAME, VISCOSITY } from './models.js';\nimport { telegraphSlowRateExcess, telegraphWaveFrequencyRatio } from './numerics.js';\n\nconst TEST = 'tests/atlas/closure.test.ts';\nconst WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';\n\n/** τ_p/t ceiling of the Langevin → Fick regime: WD4's fine resolution, t = 100 τ_p. @internal */\nconst LANGEVIN_MAX_TAU_RATIO = 0.01;\n\n/** ε = τDq² ceiling of the telegraph → Fick domain. @internal */\nexport const TELEGRAPH_FICK_MAX_EPS = 0.05;\n\n/** ε = τDq² floor of the telegraph → wave domain. @internal */\nexport const TELEGRAPH_WAVE_MIN_EPS = 25;\n\n/** Fourier-number floor of the steady-state restriction. @internal */\nconst STEADY_MIN_FOURIER = 1;\n\nconst regime = (\n  parameters: readonly { name: string; dim: import('../../dimensional/types.js').Dimension }[],\n  inequalities: AtlasBridge['regime']['inequalities'],\n): AtlasBridge['regime'] => regimeWithInputs(parameters, [], inequalities);\n\n/** A regime whose coordinates also include dimensionless model inputs (such as Re). */\nconst regimeWithInputs = (\n  parameters: readonly { name: string; dim: import('../../dimensional/types.js').Dimension }[],\n  dimensionlessInputs: readonly string[],\n  inequalities: AtlasBridge['regime']['inequalities'],\n): AtlasBridge['regime'] => makeRegime(DIFFUSION_FAMILY_NAME, parameters, dimensionlessInputs, inequalities);\n\n/** Reynolds-number ceiling of the Stokes–Einstein derivation: a chosen threshold for Re ≪ 1. @internal */\nconst STOKES_MAX_RE = 0.1;\n\n/** Ceiling of τ_p/t = m/(γt) for the overdamped limit: a chosen threshold for ≪ 1. @internal */\nconst STOKES_MAX_TAU_RATIO = 0.01;\n\n/** Bridge: Langevin → Fick, the Einstein coarse-graining D = k_BT/γ. @internal */\nexport const BRIDGE_LANGEVIN_DIFFUSION: AtlasBridge = {\n  id: 'ab-langevin-diffusion',\n  relation: 'coarse-graining',\n  premises: ['model-langevin'],\n  conclusion: 'model-fick',\n  transformation: 'average over times t ≫ τ_p = m/γ: ⟨x²⟩ → 2Dt with D = k_B T/γ (Einstein)',\n  preserves: ['the long-time mean-square displacement', 'the equilibrium velocity variance k_B T/m'],\n  doesNotPreserve: [\n    'the velocity as a state variable',\n    'the ballistic regime t ≲ τ_p, where ⟨x²⟩ ≈ (k_B T/m) t²',\n    'the momentum relaxation time τ_p as an independent scale',\n  ],\n  sideConditions: ['overdamped observation times: τ_p/t ≪ 1', 'white thermal noise, fluctuation–dissipation 2γk_BT'],\n  regime: regime(\n    [\n      { name: 'm', dim: MASS },\n      { name: 'gamma', dim: DAMPING },\n      { name: 'kT', dim: ENERGY },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 't', dim: TIME },\n    ],\n    [\n      {\n        group: 'm · gamma^-1 · t^-1',\n        op: '<=',\n        bound: LANGEVIN_MAX_TAU_RATIO,\n        alias: 'τ_p/t ≪ 1 (machine form ≤ 0.01)',\n      },\n    ],\n  ),\n  counterexamples: [\n    {\n      description:\n        'At t = 0.1 τ_p the Langevin mean-square displacement is 0.048 of 2Dt: the particle is ' +\n        'still ballistic, ⟨x²⟩ ≈ (k_BT/m)t², and the diffusion model does not describe it.',\n      witness: 'WD4b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'WD4',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '⟨x²⟩/(2Dt) within 0.02 of 1 at t = 100 τ_p, from the Langevin moment equations by RK4',\n    },\n    {\n      id: 'WD4b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '⟨x²⟩/(2Dt) = 0.0484 at t = 0.1 τ_p',\n    },\n  ],\n  citations: [\n    'Uhlenbeck & Ornstein, Phys. Rev. 36 (1930) 823 - On the theory of the Brownian motion',\n    'Einstein, Ann. Phys. 17 (1905) 549 - the relation D = k_B T/γ',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge: Langevin + Stokes drag → Fick with the Stokes–Einstein coefficient.\n *\n * A HYPEREDGE: the Langevin model supplies D = k_BT/γ with γ free; Stokes drag\n * supplies γ = 6πηa for a sphere. Neither alone gives D in terms of the fluid.\n *\n * @internal\n */\nexport const BRIDGE_STOKES_EINSTEIN: AtlasBridge = {\n  id: 'ab-stokes-einstein',\n  relation: 'derivation',\n  premises: ['model-langevin', 'model-stokes-drag'],\n  conclusion: 'model-fick',\n  transformation: 'substitute γ = 6πηa into D = k_B T/γ: D = k_B T/(6πηa)',\n  preserves: ['the long-time diffusion coefficient of a sphere'],\n  doesNotPreserve: ['the particle mass (it drops out)', 'the shape beyond the radius a'],\n  sideConditions: [\n    'creeping flow around the sphere, Re ≪ 1 (machine form Re ≤ 0.1: a chosen threshold for \"≪ 1\")',\n    'no-slip boundary',\n    'overdamped times t ≫ m/γ (machine form m/(γt) ≤ 0.01: a chosen threshold for \"≪ 1\", as in ab-langevin-diffusion)',\n  ],\n  // The prose conditions were not machine-checked, so `upt regime` reported this\n  // bridge VACUOUS and \"valid\" at Re = 1000 (persona finding L8). Re is a\n  // dimensionless input of the flow; m/(γt) is listed first so its π-group key is\n  // the one ab-langevin-diffusion uses, and one --at value checks both bridges.\n  regime: regimeWithInputs(\n    [\n      { name: 'm', dim: MASS },\n      { name: 'gamma', dim: DAMPING },\n      { name: 't', dim: TIME },\n      { name: 'kT', dim: ENERGY },\n      { name: 'eta', dim: VISCOSITY },\n      { name: 'a', dim: LENGTH },\n      { name: 'D', dim: DIFFUSIVITY },\n    ],\n    ['Re'],\n    [\n      { group: 'Re', op: '<=', bound: STOKES_MAX_RE, alias: 'Re ≪ 1 (machine form Re ≤ 0.1)' },\n      {\n        group: 'm · gamma^-1 · t^-1',\n        op: '<=',\n        bound: STOKES_MAX_TAU_RATIO,\n        alias: 'τ_p/t ≪ 1 (machine form m/(γt) ≤ 0.01)',\n      },\n    ],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD5',\n      kind: 'numeric',\n      test: TEST,\n      tolerance:\n        'agrees with CE-stokes-einstein to 1e-12 relative; a sphere of radius a = 0.5 µm (1 µm diameter) in water ' +\n        'at T = 293.15 K, η = 1.0016e-3 Pa·s gives D = 4.29e-13 m²/s',\n    },\n    {\n      id: 'WD5s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: k_BT/γ under γ ↦ 6πηa minus k_BT/(6πηa) simplifies to literal 0',\n    },\n  ],\n  citations: [\n    'Einstein, Ann. Phys. 17 (1905) 549 - Brownian motion of a sphere, D = k_B T/(6πηa)',\n    'Stokes, Trans. Camb. Phil. Soc. 9 (1851) 8 - the drag on a sphere in creeping flow',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge: telegraph → Fick, the SINGULAR limit τ → 0.\n *\n * Singular because the highest time derivative is lost: the telegraph equation\n * needs two initial conditions and diffusion takes one, so a fast mode\n * `e^{−t/τ}` forms an initial layer. The bound therefore holds only AFTER that\n * layer, and the machine horizon has a LOWER edge as well as an upper one.\n *\n * @internal\n */\nexport const BRIDGE_TELEGRAPH_DIFFUSION: AtlasBridge = {\n  id: 'ab-telegraph-diffusion',\n  relation: 'approximation',\n  premises: ['model-telegraph'],\n  conclusion: 'model-fick',\n  transformation: 'drop τ u_tt: the slow mode decays at Dq²(1 + ε + 2ε² + …) → Dq², ε = τDq²',\n  preserves: ['the slow decay rate of each Fourier mode to O(ε)', 'the total amount ∫u dx'],\n  doesNotPreserve: [\n    'the finite signal speed √(D/τ)',\n    'the second initial condition u_t(x, 0)',\n    'oscillatory modes, which exist for ε > 1/4',\n  ],\n  sideConditions: ['ε = τDq² ≤ 0.05', 'after the initial layer: t ≫ τ'],\n  bound: makeApproximation({\n    K: 1,\n    delta: telegraphSlowRateExcess(TELEGRAPH_FICK_MAX_EPS, 1, 1),\n    deltaAt: (p) =>\n      telegraphSlowRateExcess(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN),\n    deltaAtBasis: 'closed-form',\n    norm: 'relative error of the slow-mode decay rate of a Fourier mode, normalized by the value of the reduced model',\n    domain: 'ε = τDq² ≤ 0.05',\n    horizon:\n      'τ ≪ t ≪ 1/(δ D q²): after the initial layer and before the decay-rate error accumulates; ' +\n      'machine form 5τ < t < 0.1/(δ D q²)',\n    horizonHolds: (t, p) => {\n      const { tau, D, q } = p;\n      if (tau === undefined || D === undefined || q === undefined) return false;\n      const delta = telegraphSlowRateExcess(tau, D, q);\n      if (!Number.isFinite(delta)) return false;\n      if (t <= 5 * tau) return false;\n      return delta === 0 || t < 0.1 / (delta * D * q * q);\n    },\n    parameterRange: 'ε = τDq² ≤ 0.05',\n    limitCharacter: 'singular',\n    // A slow-mode decay-rate error. The horizon says that error accumulates,\n    // so time is not listed.\n    uniformity: ['slow-mode decay rate of one Fourier mode, for ε = τDq² ≤ 0.05'],\n  }),\n  regime: regime(\n    [\n      { name: 'tau', dim: TIME },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 'q', dim: WAVENUMBER },\n    ],\n    [{ group: 'tau · D · q^2', op: '<=', bound: TELEGRAPH_FICK_MAX_EPS, alias: 'ε = τDq² ≤ 0.05' }],\n  ),\n  counterexamples: [\n    {\n      description:\n        'At ε = τDq² = 1 the telegraph mode is OSCILLATORY (complex rate, NaN on the slow branch): ' +\n        'no diffusion mode oscillates.',\n      witness: 'WD6b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'WD6',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'slow rate / Dq² within 0.03 of 1 at τ = 0.025 (D = q = 1); error halves with τ',\n    },\n    {\n      id: 'WD6b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'no real slow rate at ε = 1',\n    },\n  ],\n  citations: [\n    'Cattaneo, Atti Sem. Mat. Fis. Univ. Modena 3 (1948) 83 - heat conduction with a relaxation time',\n    'Goldstein, Q. J. Mech. Appl. Math. 4 (1951) 129 - On diffusion by discontinuous movements, and on the telegraph equation',\n  ],\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-telegraph-diffusion'),\n};\n\n/**\n * Bridge: telegraph → 1-D wave, for strongly underdamped modes over short times.\n *\n * @internal\n */\nexport const BRIDGE_TELEGRAPH_WAVE: AtlasBridge = {\n  id: 'ab-telegraph-wave',\n  relation: 'approximation',\n  premises: ['model-telegraph'],\n  conclusion: 'model-wave-1d',\n  transformation: 'drop u_t: τ u_tt = D u_xx, a wave with c² = D/τ',\n  preserves: ['the signal speed √(D/τ)', 'the oscillation frequency to O(1/ε)'],\n  doesNotPreserve: ['damping: telegraph modes decay as e^{−t/(2τ)}', 'relaxation to diffusion at long times'],\n  sideConditions: ['ε = τDq² ≥ 25', 'short times t ≪ τ (the bound is not uniform in time)'],\n  bound: makeApproximation({\n    K: 1,\n    delta: 1 - telegraphWaveFrequencyRatio(TELEGRAPH_WAVE_MIN_EPS, 1, 1),\n    deltaAt: (p) =>\n      1 - telegraphWaveFrequencyRatio(p['tau'] ?? Number.NaN, p['D'] ?? Number.NaN, p['q'] ?? Number.NaN),\n    deltaAtBasis: 'closed-form',\n    norm: 'relative error of the oscillation frequency of a Fourier mode, normalized by the value of the reduced model',\n    domain: 'ε = τDq² ≥ 25',\n    horizon:\n      't ≪ τ: before damping e^{−t/(2τ)} removes 10% of the amplitude; machine form t < 2τ ln(10/9)',\n    horizonHolds: (t, p) => {\n      const { tau } = p;\n      if (tau === undefined || !(tau > 0)) return false;\n      return t >= 0 && t < 2 * tau * Math.log(10 / 9);\n    },\n    parameterRange: 'ε = τDq² ≥ 25',\n    limitCharacter: 'regular',\n    // An oscillation-frequency error. The record says it is not uniform in time.\n    uniformity: ['oscillation frequency of one Fourier mode, for ε = τDq² ≥ 25'],\n  }),\n  regime: regime(\n    [\n      { name: 'tau', dim: TIME },\n      { name: 'D', dim: DIFFUSIVITY },\n      { name: 'q', dim: WAVENUMBER },\n      { name: 'c', dim: VELOCITY },\n    ],\n    [{ group: 'tau · D · q^2', op: '>=', bound: TELEGRAPH_WAVE_MIN_EPS, alias: 'ε = τDq² ≥ 25' }],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD7',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'frequency ratio within 3e-3 of 1 at ε = 50; error halves as ε doubles',\n    },\n  ],\n  citations: ['Goldstein, Q. J. Mech. Appl. Math. 4 (1951) 129 - On diffusion by discontinuous movements, and on the telegraph equation'],\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-telegraph-wave'),\n};\n\n/** Bridge: heat → Laplace, the steady-state restriction. @internal */\nexport const BRIDGE_HEAT_LAPLACE: AtlasBridge = {\n  id: 'ab-heat-laplace',\n  relation: 'restriction',\n  premises: ['model-heat'],\n  conclusion: 'model-laplace-1d',\n  transformation: '∂T/∂t = 0: κ T_xx = 0, T linear between the end temperatures',\n  preserves: ['the end temperatures', 'the steady heat flux κ(T_b − T_a)/ℓ'],\n  doesNotPreserve: ['the transient', 'ρ and c_p, which drop out of the steady state'],\n  sideConditions: ['fixed end temperatures', 'Fourier number αt/ℓ² ≥ 1 (transients decayed)'],\n  regime: regime(\n    [\n      { name: 'kappa', dim: THERMAL_CONDUCTIVITY },\n      { name: 'rho', dim: DENSITY },\n      { name: 'cp', dim: SPECIFIC_HEAT },\n      { name: 'ell', dim: LENGTH },\n      { name: 't', dim: TIME },\n    ],\n    [\n      {\n        group: 'kappa · rho^-1 · cp^-1 · ell^-2 · t',\n        op: '>=',\n        bound: STEADY_MIN_FOURIER,\n        alias: 'Fo = αt/ℓ² ≥ 1',\n      },\n    ],\n  ),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'WD8',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'max deviation from the linear profile below 0.025 at αt/ℓ² = 0.4; shrinks with time',\n    },\n  ],\n  citations: ['Carslaw & Jaeger, Conduction of Heat in Solids - steady linear flow in a slab'],\n  reviewStatus: 'proposed',\n};\n\n/** The diffusion-family closure bridges. @internal */\nexport const DIFFUSION_CLOSURE_BRIDGES: readonly AtlasBridge[] = [\n  BRIDGE_LANGEVIN_DIFFUSION,\n  BRIDGE_STOKES_EINSTEIN,\n  BRIDGE_TELEGRAPH_DIFFUSION,\n  BRIDGE_TELEGRAPH_WAVE,\n  BRIDGE_HEAT_LAPLACE,\n];\n"]}

package/dist/atlas/diffusion/bridges.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges.d.ts","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AAIH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAO/C,qGAAqG;AACrG,eAAO,MAAM,sBAAsB,OAAO,CAAC;AAE3C,kFAAkF;AAClF,eAAO,MAAM,qBAAqB,EAAE,WAuEnC,CAAC;AAEF,gFAAgF;AAChF,eAAO,MAAM,qBAAqB,EAAE,WA0DnC,CAAC;AAEF,6EAA6E;AAC7E,eAAO,MAAM,4BAA4B,EAAE,WAkE1C,CAAC;AAEF,sEAAsE;AACtE,eAAO,MAAM,iBAAiB,EAAE,SAAS,WAAW,EAInD,CAAC"}
+{"version":3,"file":"bridges.d.ts","sourceRoot":"","sources":["../../../src/atlas/diffusion/bridges.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;GAuBG;AAIH,OAAO,KAAK,EAAE,WAAW,EAAE,MAAM,aAAa,CAAC;AAO/C,qGAAqG;AACrG,eAAO,MAAM,sBAAsB,OAAO,CAAC;AAE3C,kFAAkF;AAClF,eAAO,MAAM,qBAAqB,EAAE,WAoEnC,CAAC;AAEF,gFAAgF;AAChF,eAAO,MAAM,qBAAqB,EAAE,WAuDnC,CAAC;AAEF,6EAA6E;AAC7E,eAAO,MAAM,4BAA4B,EAAE,WA+D1C,CAAC;AAEF,sEAAsE;AACtE,eAAO,MAAM,iBAAiB,EAAE,SAAS,WAAW,EAInD,CAAC"}