universal-physics-tensor 10.0.0 → 11.0.0

package/dist/atlas/link-prediction.js.map

@@ -1 +1 @@
-{"version":3,"file":"link-prediction.js","sourceRoot":"","sources":["../../src/atlas/link-prediction.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;GA2BG;AAwCH,MAAM,MAAM,GAAG,CAAC,IAAY,EAAe,EAAE,CAC3C,IAAI,GAAG,CAAC,IAAI,CAAC,WAAW,EAAE,CAAC,KAAK,CAAC,iBAAiB,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;AAEnF,MAAM,OAAO,GAAG,CAAC,CAAsB,EAAE,CAAsB,EAAU,EAAE;IACzE,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC;QAAE,OAAO,CAAC,CAAC;IAC3C,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,KAAK,MAAM,CAAC,IAAI,CAAC;QAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAAE,KAAK,EAAE,CAAC;IACzC,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;AAC3C,CAAC,CAAC;AAEF,mFAAmF;AACnF,SAAS,IAAI,CAAC,UAA6B,EAAE,KAAwC;IACnF,OAAO,CAAC,GAAG,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACnC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;YAAE,IAAI,EAAE,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;gBAAE,OAAO,EAAE,CAAC,CAAC,CAAE,GAAG,EAAE,CAAC,CAAC,CAAE,CAAC;QAChF,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpC,CAAC,CAAC,CAAC;AACL,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,iBAAiB,CAAC,QAAgC,EAAE,CAAC,GAAG,EAAE;IACxE,MAAM,MAAM,GAAG,IAAI,GAAG,EAAsB,CAAC;IAC7C,KAAK,MAAM,CAAC,IAAI,QAAQ;QAAE,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,MAAM;YAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACxE,MAAM,OAAO,GAAG,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;IAEnD,MAAM,OAAO,GAAG,CAAC,CAA2B,EAAoC,EAAE,CAChF,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,CAAU,CAAC,CAAC;IAEpD,MAAM,MAAM,GAA0B,EAAE,CAAC;IACzC,IAAI,cAAc,GAAG,CAAC,CAAC;IACvB,KAAK,MAAM,IAAI,IAAI,OAAO,EAAE,CAAC;QAC3B,yCAAyC;QACzC,MAAM,GAAG,GAAG,IAAI,GAAG,EAAuB,CAAC;QAC3C,KAAK,MAAM,EAAE,IAAI,MAAM,CAAC,IAAI,EAAE;YAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,GAAG,EAAE,CAAC,CAAC;QACvD,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;YACxB,IAAI,CAAC,KAAK,IAAI;gBAAE,SAAS;YACzB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC;gBAChC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;gBACnB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACrB,CAAC;QACH,CAAC;QACD,KAAK,MAAM,CAAC,OAAO,EAAE,UAAU,CAAC,IAAI,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;YAClD,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YACnC,IAAI,MAAM,KAAK,SAAS;gBAAE,SAAS;YACnC,MAAM,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,GAAG,EAAU,CAAC;YACzD,IAAI,UAAU,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,CAAC;gBAC/B,cAAc,EAAE,CAAC;gBACjB,SAAS;YACX,CAAC;YACD,MAAM,UAAU,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,OAAO,IAAI,CAAC,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5F,MAAM,MAAM,GAAG,CAAC,EAAU,EAAU,EAAE;gBACpC,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,MAAM,CAAC,IAAI,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE;oBAAE,IAAI,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC;wBAAE,CAAC,EAAE,CAAC;gBAC9D,OAAO,CAAC,CAAC;YACX,CAAC,CAAC;YACF,MAAM,UAAU,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,MAAM,KAAK,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC5G,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC,QAAQ,IAAI,MAAM,CAAC,UAAU,EAAE,CAAC,CAAC;YAChE,MAAM,IAAI,GAAG,IAAI,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE;gBACnC,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC;gBAC1B,OAAO,CAAC,OAAO,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,QAAQ,IAAI,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC;YACnE,CAAC,CAAC,CAAC;YACH,MAAM,EAAE,GAAG,CAAC,IAAc,EAAU,EAAE;gBACpC,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC;gBACnC,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAClD,CAAC,CAAC;YACF,MAAM,CAAC,IAAI,CAAC;gBACV,QAAQ,EAAE,IAAI,CAAC,EAAE;gBACjB,OAAO;gBACP,UAAU;gBACV,cAAc,EAAE,EAAE,CAAC,UAAU,CAAC;gBAC9B,QAAQ,EAAE,EAAE,CAAC,IAAI,CAAC;gBAClB,UAAU,EAAE,UAAU,CAAC,MAAM;aAC9B,CAAC,CAAC;QACL,CAAC;IACH,CAAC;IACD,MAAM,IAAI,GAAG,CAAC,EAAqB,EAAU,EAAE,CAAC,CAAC,EAAE,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC,CAAC;IAC3H,OAAO;QACL,CAAC;QACD,MAAM;QACN,cAAc;QACd,cAAc,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,cAAc,IAAI,CAAC,CAAC,CAAC,MAAM;QAClE,QAAQ,EAAE,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,IAAI,CAAC,CAAC,CAAC,MAAM;QACtD,YAAY,EAAE,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC;QAC/E,aAAa,EAAE,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,cAAc,CAAC,CAAC;QAC5D,OAAO,EAAE,IAAI,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;KACjD,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Atlas Phase 6, S6.3 — link prediction over the typed model graph: ONE result,\n * stated as a hypothesis about the representation.\n *\n * Design note: `docs/planning/Atlas-Phase-6-Design.md` §1.\n *\n * The graph: every model in every registered family is a node; every\n * premise→conclusion pair of every bridge is an undirected edge. The test is\n * leave-one-bridge-out: remove a bridge's edges, then ask each predictor to rank\n * the conclusion among every model not already adjacent to the premise. A hit is\n * the held-out conclusion inside the top k.\n *\n * - **structural predictor** — common neighbours in the remaining graph, ties\n *   broken by sharing the premise's family, then by id;\n * - **text baseline** — word-token overlap between the models' dynamics and\n *   state-space descriptions, ties broken by id.\n *\n * ⚠ **What this result cannot carry.** The graph is small (24 models, 22 pairs)\n * and was curated by the same agent that runs the test, in-distribution. It is a\n * statement about how this representation behaves on its own content — not a\n * discovery claim, not a hub or gap claim, and not evidence about physics the\n * atlas does not contain. ROADMAP Phase 6 asks for it as a hypothesis test, and\n * that is its whole scope. Product A (`src/composition/discovery.ts`) is not\n * touched: a test pins that nothing under `src/atlas/` imports it.\n *\n * @module atlas/link-prediction\n * @internal\n */\n\nimport type { AtlasFamily } from './family.js';\nimport type { AtlasModel } from './model.js';\n\n/** One held-out pair and whether each predictor recovered it. @internal */\nexport interface LinkPredictionTrial {\n  readonly bridgeId: string;\n  readonly premise: string;\n  readonly conclusion: string;\n  /** 1-based rank of the held-out conclusion; Infinity when it was not a candidate. */\n  readonly structuralRank: number;\n  readonly textRank: number;\n  /** How many candidates the conclusion competed against. */\n  readonly candidates: number;\n}\n\n/** The whole leave-one-bridge-out run. @internal */\nexport interface LinkPredictionResult {\n  readonly k: number;\n  readonly trials: readonly LinkPredictionTrial[];\n  /**\n   * Held-out pairs whose endpoints stay adjacent through ANOTHER bridge. Such a\n   * pair is not a missing link — the conclusion cannot even be a candidate — so\n   * it is excluded from the trials and counted here instead of silently.\n   */\n  readonly stillConnected: number;\n  readonly structuralHits: number;\n  readonly textHits: number;\n  /**\n   * The hit rate a RANDOM ranking would get: the mean over trials of\n   * min(k, candidates)/candidates. On a small graph this is high, and a\n   * predictor's recall means nothing until it is read against it.\n   */\n  readonly chanceRecall: number;\n  /** Mean reciprocal rank of the held-out conclusion. */\n  readonly structuralMrr: number;\n  readonly textMrr: number;\n}\n\nconst tokens = (text: string): Set<string> =>\n  new Set(text.toLowerCase().split(/[^\\p{L}\\p{N}]+/u).filter((t) => t.length > 0));\n\nconst jaccard = (a: ReadonlySet<string>, b: ReadonlySet<string>): number => {\n  if (a.size === 0 && b.size === 0) return 0;\n  let inter = 0;\n  for (const x of a) if (b.has(x)) inter++;\n  return inter / (a.size + b.size - inter);\n};\n\n/** Rank candidates by descending score list (lexicographic), then ascending id. */\nfunction rank(candidates: readonly string[], score: (id: string) => readonly number[]): string[] {\n  return [...candidates].sort((x, y) => {\n    const sx = score(x);\n    const sy = score(y);\n    for (let i = 0; i < sx.length; i++) if (sx[i] !== sy[i]) return sy[i]! - sx[i]!;\n    return x < y ? -1 : x > y ? 1 : 0;\n  });\n}\n\n/**\n * Run leave-one-bridge-out link prediction over the families' model graph.\n *\n * @param families - the atlas families (normally `ATLAS_FAMILIES`).\n * @param k - the depth a hit must fall within (10, as pre-registered for recall).\n * @internal\n */\nexport function runLinkPrediction(families: readonly AtlasFamily[], k = 10): LinkPredictionResult {\n  const models = new Map<string, AtlasModel>();\n  for (const f of families) for (const m of f.models) models.set(m.id, m);\n  const bridges = families.flatMap((f) => f.bridges);\n\n  const edgesOf = (b: (typeof bridges)[number]): Array<readonly [string, string]> =>\n    b.premises.map((p) => [p, b.conclusion] as const);\n\n  const trials: LinkPredictionTrial[] = [];\n  let stillConnected = 0;\n  for (const held of bridges) {\n    // The graph WITHOUT the held-out bridge.\n    const adj = new Map<string, Set<string>>();\n    for (const id of models.keys()) adj.set(id, new Set());\n    for (const b of bridges) {\n      if (b === held) continue;\n      for (const [u, v] of edgesOf(b)) {\n        adj.get(u)?.add(v);\n        adj.get(v)?.add(u);\n      }\n    }\n    for (const [premise, conclusion] of edgesOf(held)) {\n      const pModel = models.get(premise);\n      if (pModel === undefined) continue;\n      const neighbours = adj.get(premise) ?? new Set<string>();\n      if (neighbours.has(conclusion)) {\n        stillConnected++;\n        continue;\n      }\n      const candidates = [...models.keys()].filter((id) => id !== premise && !neighbours.has(id));\n      const common = (id: string): number => {\n        let n = 0;\n        for (const x of adj.get(id) ?? []) if (neighbours.has(x)) n++;\n        return n;\n      };\n      const structural = rank(candidates, (id) => [common(id), models.get(id)!.family === pModel.family ? 1 : 0]);\n      const pText = tokens(`${pModel.dynamics} ${pModel.stateSpace}`);\n      const text = rank(candidates, (id) => {\n        const m = models.get(id)!;\n        return [jaccard(pText, tokens(`${m.dynamics} ${m.stateSpace}`))];\n      });\n      const at = (list: string[]): number => {\n        const i = list.indexOf(conclusion);\n        return i < 0 ? Number.POSITIVE_INFINITY : i + 1;\n      };\n      trials.push({\n        bridgeId: held.id,\n        premise,\n        conclusion,\n        structuralRank: at(structural),\n        textRank: at(text),\n        candidates: candidates.length,\n      });\n    }\n  }\n  const mean = (xs: readonly number[]): number => (xs.length === 0 ? Number.NaN : xs.reduce((a, b) => a + b, 0) / xs.length);\n  return {\n    k,\n    trials,\n    stillConnected,\n    structuralHits: trials.filter((t) => t.structuralRank <= k).length,\n    textHits: trials.filter((t) => t.textRank <= k).length,\n    chanceRecall: mean(trials.map((t) => Math.min(k, t.candidates) / t.candidates)),\n    structuralMrr: mean(trials.map((t) => 1 / t.structuralRank)),\n    textMrr: mean(trials.map((t) => 1 / t.textRank)),\n  };\n}\n"]}
+{"version":3,"file":"link-prediction.js","sourceRoot":"","sources":["../../src/atlas/link-prediction.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAuCG;AAEH,OAAO,EAAE,WAAW,EAAE,GAAG,EAAE,eAAe,EAAE,MAAM,iCAAiC,CAAC;AAIpF,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAC5D,OAAO,EAAE,OAAO,EAAE,MAAM,sBAAsB,CAAC;AAC/C,OAAO,EAAE,MAAM,EAAE,MAAM,4BAA4B,CAAC;AAkDpD,uGAAuG;AACvG,MAAM,CAAC,MAAM,gBAAgB,GAC3B,kGAA2G,CAAC;AA+B9G,qDAAqD;AACrD,MAAM,CAAC,MAAM,cAAc,GAAG,wDAAiE,CAAC;AAyChG,mGAAmG;AACnG,MAAM,aAAa,GAAG,QAAQ,CAAC;AAC/B,MAAM,SAAS,GAAG,KAAK,CAAC;AAaxB,uFAAuF;AACvF,SAAS,QAAQ,CAAC,QAAgC;IAKhD,MAAM,MAAM,GAAG,IAAI,GAAG,EAAsB,CAAC;IAC7C,KAAK,MAAM,CAAC,IAAI,QAAQ;QAAE,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,MAAM;YAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;IACxE,MAAM,OAAO,GAAG,QAAQ,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;IACnD,MAAM,OAAO,GAAG,CAAC,CAAc,EAAoC,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,CAAU,CAAC,CAAC;IAExH,MAAM,OAAO,GAAc,EAAE,CAAC;IAC9B,IAAI,cAAc,GAAG,CAAC,CAAC;IACvB,KAAK,MAAM,IAAI,IAAI,OAAO,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,IAAI,GAAG,EAAuB,CAAC;QAC3C,KAAK,MAAM,EAAE,IAAI,MAAM,CAAC,IAAI,EAAE;YAAE,GAAG,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,GAAG,EAAE,CAAC,CAAC;QACvD,KAAK,MAAM,CAAC,IAAI,OAAO,EAAE,CAAC;YACxB,IAAI,CAAC,KAAK,IAAI;gBAAE,SAAS;YACzB,KAAK,MAAM,CAAC,CAAC,EAAE,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC;gBAChC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;gBACnB,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;YACrB,CAAC;QACH,CAAC;QACD,KAAK,MAAM,CAAC,OAAO,EAAE,UAAU,CAAC,IAAI,OAAO,CAAC,IAAI,CAAC,EAAE,CAAC;YAClD,MAAM,YAAY,GAAG,MAAM,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;YACzC,IAAI,YAAY,KAAK,SAAS;gBAAE,SAAS;YACzC,MAAM,UAAU,GAAG,GAAG,CAAC,GAAG,CAAC,OAAO,CAAC,IAAI,IAAI,GAAG,EAAU,CAAC;YACzD,IAAI,UAAU,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,CAAC;gBAC/B,cAAc,EAAE,CAAC;gBACjB,SAAS;YACX,CAAC;YACD,MAAM,UAAU,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,MAAM,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,KAAK,OAAO,IAAI,CAAC,UAAU,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5F,OAAO,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,OAAO,EAAE,UAAU,EAAE,YAAY,EAAE,GAAG,EAAE,UAAU,EAAE,CAAC,CAAC;QAC7E,CAAC;IACH,CAAC;IACD,OAAO,EAAE,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,CAAC;AAC7C,CAAC;AAED,mEAAmE;AACnE,MAAM,MAAM,GAAG,CAAC,IAAuB,EAAE,EAAU,EAAU,EAAE;IAC7D,MAAM,CAAC,GAAG,IAAI,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC;IAC3B,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AAClD,CAAC,CAAC;AAEF,qGAAqG;AACrG,SAAS,aAAa,CAAC,MAAyB;IAC9C,MAAM,EAAE,GAAG,WAAW,CAAC,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,MAAM,EAAE,EAAE,CAAC,MAAM,CAAC,MAAM,CAAC,EAAE;QAC9D,UAAU,EAAE,IAAI;QAChB,SAAS,EAAE,SAAS;QACpB,IAAI,EAAE,aAAa;KACpB,CAAC,CAAC;IACH,OAAO,EAAE,QAAQ,EAAE,EAAE,CAAC,QAAQ,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,EAAE,CAAC,KAAK,EAAE,UAAU,EAAE,EAAE,CAAC,UAAU,EAAE,CAAC;AAChG,CAAC;AAED,gGAAgG;AAChG,MAAM,SAAS,GAAG,CAAC,EAAqB,EAAU,EAAE,CAAC,CAAC,EAAE,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;AAEjG;;;;;;;GAOG;AACH,MAAM,UAAU,iBAAiB,CAAC,QAAgC,EAAE,CAAC,GAAG,EAAE;IACxE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC/D,MAAM,MAAM,GAA0B,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;QACtD,MAAM,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAC,IAAI,IAAI,GAAG,EAAU,CAAC;QAC7D,MAAM,MAAM,GAAG,CAAC,EAAU,EAAU,EAAE;YACpC,IAAI,CAAC,GAAG,CAAC,CAAC;YACV,KAAK,MAAM,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE;gBAAE,IAAI,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC;oBAAE,CAAC,EAAE,CAAC;YAChE,OAAO,CAAC,CAAC;QACX,CAAC,CAAC;QACF,MAAM,UAAU,GAAG,MAAM,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,MAAM,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,MAAM,KAAK,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACxH,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,YAAY,CAAC,QAAQ,IAAI,CAAC,CAAC,YAAY,CAAC,UAAU,EAAE,CAAC,CAAC;QAChF,MAAM,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE;YACvC,MAAM,CAAC,GAAG,MAAM,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC;YAC1B,OAAO,CAAC,OAAO,CAAC,KAAK,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,QAAQ,IAAI,CAAC,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC;QACnE,CAAC,CAAC,CAAC;QACH,OAAO;YACL,QAAQ,EAAE,CAAC,CAAC,IAAI,CAAC,EAAE;YACnB,OAAO,EAAE,CAAC,CAAC,OAAO;YAClB,UAAU,EAAE,CAAC,CAAC,UAAU;YACxB,cAAc,EAAE,MAAM,CAAC,UAAU,EAAE,CAAC,CAAC,UAAU,CAAC;YAChD,QAAQ,EAAE,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC,UAAU,CAAC;YACpC,UAAU,EAAE,CAAC,CAAC,UAAU,CAAC,MAAM;SAChC,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,MAAM,aAAa,GAAa,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,cAAc,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnF,MAAM,OAAO,GAAa,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACvE,MAAM,YAAY,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,cAAc,CAAC,CAAC;IAC7D,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;IAEjD,mEAAmE;IACnE,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,IAAI,CAAC,GAAG,CAAC,CAAC;IACV,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE;QACtB,IAAI,aAAa,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC;YAAE,CAAC,EAAE,CAAC;aACnC,IAAI,aAAa,CAAC,CAAC,CAAC;YAAE,CAAC,EAAE,CAAC;aAC1B,IAAI,OAAO,CAAC,CAAC,CAAC;YAAE,CAAC,EAAE,CAAC;;YACpB,CAAC,EAAE,CAAC;IACX,CAAC,CAAC,CAAC;IAEH,MAAM,QAAQ,GAAG,eAAe,CAAC,YAAY,EAAE,MAAM,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC,CAAC,EAAE;QAC5F,SAAS,EAAE,SAAS;QACpB,IAAI,EAAE,aAAa;KACpB,CAAC,CAAC;IAEH,OAAO;QACL,CAAC;QACD,MAAM;QACN,cAAc;QACd,cAAc,EAAE,aAAa,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;QAC5D,QAAQ,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;QAChD,YAAY,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC;QACpF,aAAa,EAAE,SAAS,CAAC,YAAY,CAAC;QACtC,OAAO,EAAE,SAAS,CAAC,MAAM,CAAC;QAC1B,SAAS,EAAE;YACT,IAAI,EAAE,aAAa;YACnB,SAAS,EAAE,SAAS;YACpB,gBAAgB,EAAE,aAAa,CAAC,aAAa,CAAC;YAC9C,UAAU,EAAE,aAAa,CAAC,OAAO,CAAC;YAClC,aAAa,EAAE,aAAa,CAAC,YAAY,CAAC;YAC1C,OAAO,EAAE,aAAa,CAAC,MAAM,CAAC;SAC/B;QACD,WAAW,EAAE;YACX,SAAS,EAAE,QAAQ,CAAC,SAAS;YAC7B,MAAM,EAAE,QAAQ,CAAC,MAAM;YACvB,IAAI,EAAE,aAAa;YACnB,SAAS,EAAE,SAAS;YACpB,IAAI,EAAE,gBAAgB;SACvB;QACD,aAAa,EAAE,OAAO,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,MAAM;KAC9C,CAAC;AACJ,CAAC;AAED,wGAAwG;AACxG,SAAS,YAAY,CAAC,KAAc,EAAE,CAAS;IAC7C,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAChD,MAAM,IAAI,SAAS,CAAC,wBAAwB,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;IACD,KAAK,MAAM,GAAG,IAAI,KAAK,EAAE,CAAC;QACxB,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;YAC9G,MAAM,IAAI,SAAS,CAAC,wBAAwB,CAAC,IAAI,CAAC,2BAA2B,CAAC,CAAC;QACjF,CAAC;IACH,CAAC;IACD,OAAO,KAAmB,CAAC;AAC7B,CAAC;AAED;;;;;GAKG;AACH,SAAS,SAAS,CAAC,GAA6C,EAAE,GAAsB;IACtF,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC;IACrB,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,CAAC,GAAe,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACpF,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,KAAK,MAAM,CAAC,CAAC,EAAE,UAAU,CAAC,IAAI,GAAG,EAAE,CAAC;QAClC,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvB,IAAI,CAAC,KAAK,SAAS;YAAE,SAAS;QAC9B,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,UAAU,CAAC,IAAI,CAAC,CAAC;QACjD,KAAK,MAAM,CAAC,IAAI,UAAU,EAAE,CAAC;YAC3B,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YACvB,IAAI,CAAC,KAAK,SAAS;gBAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACnC,CAAC;IACH,CAAC;IACD,MAAM,IAAI,GAAG,SAAS,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC;IACzD,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;IAC7E,MAAM,OAAO,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACxC,MAAM,CAAC,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5E,OAAO,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC;AACrB,CAAC;AAED;;;;;;;;GAQG;AACH,MAAM,UAAU,qBAAqB,CAAC,QAAgC,EAAE,CAAC,GAAG,EAAE;IAC5E,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,GAAG,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC/D,MAAM,GAAG,GAAG,CAAC,GAAG,MAAM,CAAC,IAAI,EAAE,CAAC,CAAC;IAC/B,MAAM,KAAK,GAAG,IAAI,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACnD,MAAM,MAAM,GAAgB,OAAO,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE;QAC5C,MAAM,IAAI,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,OAAO,CAAE,CAAC;QAChC,MAAM,KAAK,GAAG,CAAC,EAAU,EAAU,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAE,CAAC,CAAC;QAChE,MAAM,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QACzD,OAAO;YACL,QAAQ,EAAE,CAAC,CAAC,IAAI,CAAC,EAAE;YACnB,OAAO,EAAE,CAAC,CAAC,OAAO;YAClB,UAAU,EAAE,CAAC,CAAC,UAAU;YACxB,QAAQ,EAAE,MAAM,CAAC,MAAM,EAAE,CAAC,CAAC,UAAU,CAAC;YACtC,SAAS,EAAE,CAAC,CAAC,UAAU,CAAC,QAAQ,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;YACxE,UAAU,EAAE,CAAC,CAAC,UAAU,CAAC,MAAM;YAC/B,IAAI,EAAE,IAAI,CAAC,IAAI;YACf,OAAO,EAAE,IAAI;SACd,CAAC;IACJ,CAAC,CAAC,CAAC;IACH,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;IACxC,MAAM,GAAG,GAAa,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACnE,MAAM,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,QAAQ,CAAC,CAAC;IAC7C,OAAO;QACL,OAAO,EAAE,IAAI;QACb,KAAK,EAAE,cAAc;QACrB,CAAC;QACD,SAAS,EAAE,KAAK,CAAC,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,EAAE,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC;QACnG,MAAM;QACN,cAAc;QACd,QAAQ,EAAE,GAAG,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;QAC5C,OAAO,EAAE,SAAS,CAAC,EAAE,CAAC;QACtB,YAAY,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC;QACpF,SAAS,EAAE,EAAE,IAAI,EAAE,aAAa,EAAE,SAAS,EAAE,SAAS,EAAE,MAAM,EAAE,aAAa,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,aAAa,CAAC,EAAE,CAAC,EAAE;KAC7G,CAAC;AACJ,CAAC","sourcesContent":["/**\n * Atlas Phase 6, S6.3 — link prediction over the typed model graph: ONE result,\n * stated as a hypothesis about the representation.\n *\n * Design note: `docs/planning/Atlas-Phase-6-Design.md` §1.\n *\n * The graph: every model in every registered family is a node; every\n * premise→conclusion pair of every bridge is an undirected edge. The test is\n * leave-one-bridge-out: remove a bridge's edges, then ask each predictor to rank\n * the conclusion among every model not already adjacent to the premise. A hit is\n * the held-out conclusion inside the top k.\n *\n * - **structural predictor** — common neighbours in the remaining graph, ties\n *   broken by sharing the premise's family, then by id;\n * - **text baseline** — word-token overlap between the models' dynamics and\n *   state-space descriptions (`text-overlap.ts`), ties broken by id.\n *\n * **The recorded result stands.** The negative result (typed neighbourhood below text\n * overlap on twenty trials; NOTES.md, \"Link prediction is NEGATIVE\") is reproduced by\n * {@link runLinkPrediction} exactly as recorded and is not re-run for a better answer. What\n * was added beside it, never inside it: a seeded percentile-bootstrap interval on every\n * recall and MRR (MathTS `bootstrapCI`), a seeded permutation test on the two per-trial\n * reciprocal-rank vectors (MathTS `permutationTest`, stated for what it is), and McNemar's\n * exact p on the paired hits (`benchmark/stats.ts`). The Katz-index predictor\n * ({@link runKatzLinkPrediction}) is a SEPARATE result: every record it returns carries\n * `postHoc: true` and the label, and it is never a field of the recorded comparison.\n *\n * ⚠ **What this result cannot carry.** The graph is small (24 models, 22 pairs)\n * and was curated by the same agent that runs the test, in-distribution. It is a\n * statement about how this representation behaves on its own content — not a\n * discovery claim, not a hub or gap claim, and not evidence about physics the\n * atlas does not contain. ROADMAP Phase 6 asks for it as a hypothesis test, and\n * that is its whole scope. Product A (`src/composition/discovery.ts`) is not\n * touched: a test pins that nothing under `src/atlas/` imports it, and the Katz\n * predictor is not wired into discovery, whose candidates are in different\n * components by construction (every path score there is 0).\n *\n * @module atlas/link-prediction\n * @internal\n */\n\nimport { bootstrapCI, inv, permutationTest } from '@danielsimonjr/mathts-functions';\nimport type { AtlasFamily } from './family.js';\nimport type { AtlasModel } from './model.js';\nimport type { AtlasBridge } from './types.js';\nimport { jaccard, rankBy, tokens } from './text-overlap.js';\nimport { mcnemar } from './benchmark/stats.js';\nimport { meanOf } from '../numerical/statistics.js';\n\n/** One held-out pair and whether each predictor recovered it. @internal */\nexport interface LinkPredictionTrial {\n  readonly bridgeId: string;\n  readonly premise: string;\n  readonly conclusion: string;\n  /** 1-based rank of the held-out conclusion; Infinity when it was not a candidate. */\n  readonly structuralRank: number;\n  readonly textRank: number;\n  /** How many candidates the conclusion competed against. */\n  readonly candidates: number;\n}\n\n/** A percentile-bootstrap interval with the point estimate it brackets (MathTS `bootstrapCI`). @internal */\nexport interface BootstrapInterval {\n  readonly estimate: number;\n  readonly lower: number;\n  readonly upper: number;\n  /** The confidence level, 0.95. */\n  readonly confidence: number;\n}\n\n/** The seeded intervals of the recorded comparison. @internal */\nexport interface LinkPredictionIntervals {\n  /** The seed every resample used; stated so the interval is reproducible. */\n  readonly seed: number;\n  readonly resamples: number;\n  readonly structuralRecall: BootstrapInterval;\n  readonly textRecall: BootstrapInterval;\n  readonly structuralMrr: BootstrapInterval;\n  readonly textMrr: BootstrapInterval;\n}\n\n/** The seeded permutation test on the two per-trial reciprocal-rank vectors. @internal */\nexport interface ReciprocalRankPermutation {\n  /** Mean structural reciprocal rank minus mean text reciprocal rank (the observed MRR difference). */\n  readonly statistic: number;\n  /** Two-tailed p from MathTS `permutationTest`. */\n  readonly pValue: number;\n  readonly seed: number;\n  readonly resamples: number;\n  /**\n   * What the test is. MathTS has no paired permutation: `permutationTest` pools the two vectors\n   * and re-splits them, so the trial pairing is not kept and the null is looser than a sign-flip\n   * test's. The paired test on the hits is `mcnemarExactP`.\n   */\n  readonly kind: typeof PERMUTATION_KIND;\n}\n\n/** The statement of the permutation test's kind, so no reader takes it for a paired test. @internal */\nexport const PERMUTATION_KIND =\n  'two-sample permutation of the reciprocal-rank vectors (pooled re-split; the pairing is not kept)' as const;\n\n/** The whole leave-one-bridge-out run. @internal */\nexport interface LinkPredictionResult {\n  readonly k: number;\n  readonly trials: readonly LinkPredictionTrial[];\n  /**\n   * Held-out pairs whose endpoints stay adjacent through ANOTHER bridge. Such a\n   * pair is not a missing link — the conclusion cannot even be a candidate — so\n   * it is excluded from the trials and counted here instead of silently.\n   */\n  readonly stillConnected: number;\n  readonly structuralHits: number;\n  readonly textHits: number;\n  /**\n   * The hit rate a RANDOM ranking would get: the mean over trials of\n   * min(k, candidates)/candidates. On a small graph this is high, and a\n   * predictor's recall means nothing until it is read against it.\n   */\n  readonly chanceRecall: number;\n  /** Mean reciprocal rank of the held-out conclusion. */\n  readonly structuralMrr: number;\n  readonly textMrr: number;\n  /** Seeded bootstrap intervals on the four numbers above. */\n  readonly intervals: LinkPredictionIntervals;\n  /** The seeded permutation test on the per-trial reciprocal ranks. */\n  readonly permutation: ReciprocalRankPermutation;\n  /** McNemar's exact p on the paired hit table (both, structural only, text only, neither). */\n  readonly mcnemarExactP: number;\n}\n\n/** The label every Katz record carries. @internal */\nexport const POST_HOC_LABEL = 'post hoc: added after the negative result was recorded' as const;\n\n/** One held-out pair under the Katz-index predictor. @internal */\nexport interface KatzTrial {\n  readonly bridgeId: string;\n  readonly premise: string;\n  readonly conclusion: string;\n  /** 1-based rank of the held-out conclusion by Katz index; Infinity when it was not a candidate. */\n  readonly katzRank: number;\n  /** The Katz index K[premise][conclusion] = ((I − βA)⁻¹ − I)[premise][conclusion]; 0 when no path. */\n  readonly katzScore: number;\n  readonly candidates: number;\n  /** This trial's attenuation β = 1/(2·max degree) of its own held-out graph. */\n  readonly beta: number;\n  readonly postHoc: true;\n}\n\n/** The Katz-index run: a separate, post-hoc result, never merged into the recorded comparison. @internal */\nexport interface KatzLinkPredictionResult {\n  readonly postHoc: true;\n  readonly label: typeof POST_HOC_LABEL;\n  readonly k: number;\n  /**\n   * `[min, max]` of the trials' attenuations. Each held-out graph has its own maximum degree, so\n   * β = 1/(2·max degree) is per trial (`KatzTrial.beta`), each below its graph's spectral bound\n   * 1/λ_max since λ_max ≤ max degree. The result once reported one β, the last trial's.\n   */\n  readonly betaRange: readonly [number, number];\n  readonly trials: readonly KatzTrial[];\n  readonly stillConnected: number;\n  readonly katzHits: number;\n  readonly katzMrr: number;\n  readonly chanceRecall: number;\n  readonly intervals: {\n    readonly seed: number;\n    readonly resamples: number;\n    readonly recall: BootstrapInterval;\n    readonly mrr: BootstrapInterval;\n  };\n}\n\n/** The one seed every resample in this module uses, and the resample count. Stated, not hidden. */\nconst RESAMPLE_SEED = 20260924;\nconst RESAMPLES = 10000;\n\n/** One leave-one-bridge-out query: the graph without the bridge, and the candidates. */\ninterface HoldOut {\n  readonly held: AtlasBridge;\n  readonly premise: string;\n  readonly conclusion: string;\n  readonly premiseModel: AtlasModel;\n  /** The undirected graph WITHOUT the held-out bridge. */\n  readonly adj: ReadonlyMap<string, ReadonlySet<string>>;\n  readonly candidates: readonly string[];\n}\n\n/** Every held-out query of the families, and the pairs excluded as still connected. */\nfunction holdOuts(families: readonly AtlasFamily[]): {\n  readonly models: ReadonlyMap<string, AtlasModel>;\n  readonly queries: readonly HoldOut[];\n  readonly stillConnected: number;\n} {\n  const models = new Map<string, AtlasModel>();\n  for (const f of families) for (const m of f.models) models.set(m.id, m);\n  const bridges = families.flatMap((f) => f.bridges);\n  const edgesOf = (b: AtlasBridge): Array<readonly [string, string]> => b.premises.map((p) => [p, b.conclusion] as const);\n\n  const queries: HoldOut[] = [];\n  let stillConnected = 0;\n  for (const held of bridges) {\n    const adj = new Map<string, Set<string>>();\n    for (const id of models.keys()) adj.set(id, new Set());\n    for (const b of bridges) {\n      if (b === held) continue;\n      for (const [u, v] of edgesOf(b)) {\n        adj.get(u)?.add(v);\n        adj.get(v)?.add(u);\n      }\n    }\n    for (const [premise, conclusion] of edgesOf(held)) {\n      const premiseModel = models.get(premise);\n      if (premiseModel === undefined) continue;\n      const neighbours = adj.get(premise) ?? new Set<string>();\n      if (neighbours.has(conclusion)) {\n        stillConnected++;\n        continue;\n      }\n      const candidates = [...models.keys()].filter((id) => id !== premise && !neighbours.has(id));\n      queries.push({ held, premise, conclusion, premiseModel, adj, candidates });\n    }\n  }\n  return { models, queries, stillConnected };\n}\n\n/** 1-based position of `id` in a ranking; Infinity when absent. */\nconst rankOf = (list: readonly string[], id: string): number => {\n  const i = list.indexOf(id);\n  return i < 0 ? Number.POSITIVE_INFINITY : i + 1;\n};\n\n/** Seeded percentile bootstrap of the mean of per-trial values (a 0/1 hit, or a reciprocal rank). */\nfunction bootstrapMean(values: readonly number[]): BootstrapInterval {\n  const ci = bootstrapCI([...values], (sample) => meanOf(sample), {\n    confidence: 0.95,\n    resamples: RESAMPLES,\n    seed: RESAMPLE_SEED,\n  });\n  return { estimate: ci.estimate, lower: ci.lower, upper: ci.upper, confidence: ci.confidence };\n}\n\n/** The mean of a trial statistic; NaN with no trials, as the recorded result always read it. */\nconst meanOrNaN = (xs: readonly number[]): number => (xs.length === 0 ? Number.NaN : meanOf(xs));\n\n/**\n * Run leave-one-bridge-out link prediction over the families' model graph. The recorded\n * comparison (structural against text), with its seeded intervals and paired statistics.\n *\n * @param families - the atlas families (normally `ATLAS_FAMILIES`).\n * @param k - the depth a hit must fall within (10, as pre-registered for recall).\n * @internal\n */\nexport function runLinkPrediction(families: readonly AtlasFamily[], k = 10): LinkPredictionResult {\n  const { models, queries, stillConnected } = holdOuts(families);\n  const trials: LinkPredictionTrial[] = queries.map((q) => {\n    const neighbours = q.adj.get(q.premise) ?? new Set<string>();\n    const common = (id: string): number => {\n      let n = 0;\n      for (const x of q.adj.get(id) ?? []) if (neighbours.has(x)) n++;\n      return n;\n    };\n    const structural = rankBy(q.candidates, (id) => [common(id), models.get(id)!.family === q.premiseModel.family ? 1 : 0]);\n    const pText = tokens(`${q.premiseModel.dynamics} ${q.premiseModel.stateSpace}`);\n    const text = rankBy(q.candidates, (id) => {\n      const m = models.get(id)!;\n      return [jaccard(pText, tokens(`${m.dynamics} ${m.stateSpace}`))];\n    });\n    return {\n      bridgeId: q.held.id,\n      premise: q.premise,\n      conclusion: q.conclusion,\n      structuralRank: rankOf(structural, q.conclusion),\n      textRank: rankOf(text, q.conclusion),\n      candidates: q.candidates.length,\n    };\n  });\n\n  const structuralHit: number[] = trials.map((t) => (t.structuralRank <= k ? 1 : 0));\n  const textHit: number[] = trials.map((t) => (t.textRank <= k ? 1 : 0));\n  const structuralRR = trials.map((t) => 1 / t.structuralRank);\n  const textRR = trials.map((t) => 1 / t.textRank);\n\n  // The paired hit table: both, structural only, text only, neither.\n  let a = 0;\n  let b = 0;\n  let c = 0;\n  let d = 0;\n  trials.forEach((_, i) => {\n    if (structuralHit[i] && textHit[i]) a++;\n    else if (structuralHit[i]) b++;\n    else if (textHit[i]) c++;\n    else d++;\n  });\n\n  const permuted = permutationTest(structuralRR, textRR, (x, y) => meanOrNaN(x) - meanOrNaN(y), {\n    resamples: RESAMPLES,\n    seed: RESAMPLE_SEED,\n  });\n\n  return {\n    k,\n    trials,\n    stillConnected,\n    structuralHits: structuralHit.reduce((acc, x) => acc + x, 0),\n    textHits: textHit.reduce((acc, x) => acc + x, 0),\n    chanceRecall: meanOrNaN(trials.map((t) => Math.min(k, t.candidates) / t.candidates)),\n    structuralMrr: meanOrNaN(structuralRR),\n    textMrr: meanOrNaN(textRR),\n    intervals: {\n      seed: RESAMPLE_SEED,\n      resamples: RESAMPLES,\n      structuralRecall: bootstrapMean(structuralHit),\n      textRecall: bootstrapMean(textHit),\n      structuralMrr: bootstrapMean(structuralRR),\n      textMrr: bootstrapMean(textRR),\n    },\n    permutation: {\n      statistic: permuted.statistic,\n      pValue: permuted.pValue,\n      seed: RESAMPLE_SEED,\n      resamples: RESAMPLES,\n      kind: PERMUTATION_KIND,\n    },\n    mcnemarExactP: mcnemar({ a, b, c, d }).exactP,\n  };\n}\n\n/** A square number matrix of order n, read back from MathTS `inv`, which types its result `unknown`. */\nfunction squareMatrix(value: unknown, n: number): number[][] {\n  if (!Array.isArray(value) || value.length !== n) {\n    throw new TypeError(`inv did not return a ${n}×${n} matrix`);\n  }\n  for (const row of value) {\n    if (!Array.isArray(row) || row.length !== n || !row.every((x) => typeof x === 'number' && Number.isFinite(x))) {\n      throw new TypeError(`inv did not return a ${n}×${n} matrix of finite numbers`);\n    }\n  }\n  return value as number[][];\n}\n\n/**\n * The Katz index of the graph without the held-out bridge: K = (I − βA)⁻¹ − I through MathTS\n * `inv`, with β = 1/(2·max degree). The largest eigenvalue of a 0/1 adjacency is at most the\n * maximum degree, so β is below the spectral bound 1/λ_max and the series Σ βˡ Aˡ converges.\n * A graph with no edge has K = 0 everywhere, and β is stated as 0.5 (there is nothing to attenuate).\n */\nfunction katzIndex(adj: ReadonlyMap<string, ReadonlySet<string>>, ids: readonly string[]): { K: number[][]; beta: number } {\n  const n = ids.length;\n  const index = new Map(ids.map((id, i) => [id, i]));\n  const A: number[][] = Array.from({ length: n }, () => new Array<number>(n).fill(0));\n  let maxDegree = 0;\n  for (const [u, neighbours] of adj) {\n    const i = index.get(u);\n    if (i === undefined) continue;\n    maxDegree = Math.max(maxDegree, neighbours.size);\n    for (const v of neighbours) {\n      const j = index.get(v);\n      if (j !== undefined) A[i][j] = 1;\n    }\n  }\n  const beta = maxDegree === 0 ? 0.5 : 1 / (2 * maxDegree);\n  const M = A.map((row, i) => row.map((x, j) => (i === j ? 1 : 0) - beta * x));\n  const inverse = squareMatrix(inv(M), n);\n  const K = inverse.map((row, i) => row.map((x, j) => (i === j ? x - 1 : x)));\n  return { K, beta };\n}\n\n/**\n * The Katz-index predictor, SEPARATE from the recorded comparison and labelled post hoc: it was\n * added after the negative result was recorded, so its numbers are exploratory and are never\n * merged into that result. Ranks each candidate by K[premise][candidate] (paths of every length,\n * attenuated by β per step), ties by id. Not wired into discovery (`src/composition/`): Product A\n * candidates are in different components by construction, so every path score there is 0.\n *\n * @internal\n */\nexport function runKatzLinkPrediction(families: readonly AtlasFamily[], k = 10): KatzLinkPredictionResult {\n  const { models, queries, stillConnected } = holdOuts(families);\n  const ids = [...models.keys()];\n  const index = new Map(ids.map((id, i) => [id, i]));\n  const trials: KatzTrial[] = queries.map((q) => {\n    const katz = katzIndex(q.adj, ids);\n    const p = index.get(q.premise)!;\n    const score = (id: string): number => katz.K[p][index.get(id)!];\n    const ranked = rankBy(q.candidates, (id) => [score(id)]);\n    return {\n      bridgeId: q.held.id,\n      premise: q.premise,\n      conclusion: q.conclusion,\n      katzRank: rankOf(ranked, q.conclusion),\n      katzScore: q.candidates.includes(q.conclusion) ? score(q.conclusion) : 0,\n      candidates: q.candidates.length,\n      beta: katz.beta,\n      postHoc: true,\n    };\n  });\n  const betas = trials.map((t) => t.beta);\n  const hit: number[] = trials.map((t) => (t.katzRank <= k ? 1 : 0));\n  const rr = trials.map((t) => 1 / t.katzRank);\n  return {\n    postHoc: true,\n    label: POST_HOC_LABEL,\n    k,\n    betaRange: betas.length === 0 ? [Number.NaN, Number.NaN] : [Math.min(...betas), Math.max(...betas)],\n    trials,\n    stillConnected,\n    katzHits: hit.reduce((acc, x) => acc + x, 0),\n    katzMrr: meanOrNaN(rr),\n    chanceRecall: meanOrNaN(trials.map((t) => Math.min(k, t.candidates) / t.candidates)),\n    intervals: { seed: RESAMPLE_SEED, resamples: RESAMPLES, recall: bootstrapMean(hit), mrr: bootstrapMean(rr) },\n  };\n}\n"]}

package/dist/atlas/oscillators/bridges-coarse.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-coarse.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-coarse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,KAAK,EAAE,WAAW,EAAE,gBAAgB,EAAE,MAAM,aAAa,CAAC;AAMjE;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAExF;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE1F;AAED,+DAA+D;AAC/D,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE1E;AAED;;;;;;GAMG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEhE;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,4BAA4B,CAAC,EAAE,EAAE,MAAM,GAAG,MAAM,CAE/D;AAED,oEAAoE;AACpE,eAAO,MAAM,iBAAiB,EAAE,WAkE/B,CAAC;AAEF;;;;;;;;;;;;;;GAcG;AACH,eAAO,MAAM,mBAAmB,EAAE,gBAAwD,CAAC"}
+{"version":3,"file":"bridges-coarse.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-coarse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAIH,OAAO,KAAK,EAAE,WAAW,EAAE,gBAAgB,EAAE,MAAM,aAAa,CAAC;AAMjE;;;;;GAKG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAExF;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE1F;AAED,+DAA+D;AAC/D,wBAAgB,cAAc,CAAC,CAAC,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAE1E;AAED;;;;;;GAMG;AACH,wBAAgB,eAAe,CAAC,KAAK,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,MAAM,CAEhE;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,4BAA4B,CAAC,EAAE,EAAE,MAAM,GAAG,MAAM,CAE/D;AAED,oEAAoE;AACpE,eAAO,MAAM,iBAAiB,EAAE,WA+D/B,CAAC;AAEF;;;;;;;;;;;;;;GAcG;AACH,eAAO,MAAM,mBAAmB,EAAE,gBAAwD,CAAC"}

package/dist/atlas/oscillators/bridges-coarse.js

@@ -11,7 +11,7 @@
  * @module atlas/oscillators/bridges-coarse
  */
 import { LENGTH, MASS, VELOCITY } from '../../dimensional/types.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { relationContractOf } from './bridges-exact.js';
 import { SPRING_CONSTANT } from './dimensions.js';
 const FAMILY = 'oscillators';
@@ -82,16 +82,12 @@ export const BRIDGE_CHAIN_WAVE = {
         'the discrete lattice spacing a as an independent scale',
     ],
     sideConditions: ['long-wavelength, qa ≪ 1'],
-    regime: {
-        family: FAMILY,
-        inequalities: [{ group: 'qa', op: '<', bound: 1, alias: 'qa ≪ 1' }],
-        groupDefinitions: deriveRegimeGroups(FAMILY, [
-            { name: 'm', dim: MASS },
-            { name: 'kappa', dim: SPRING_CONSTANT },
-            { name: 'a', dim: LENGTH },
-            { name: 'c', dim: VELOCITY },
-        ], ['qa']),
-    },
+    regime: makeRegime(FAMILY, [
+        { name: 'm', dim: MASS },
+        { name: 'kappa', dim: SPRING_CONSTANT },
+        { name: 'a', dim: LENGTH },
+        { name: 'c', dim: VELOCITY },
+    ], ['qa'], [{ group: 'qa', op: '<', bound: 1, alias: 'qa ≪ 1' }]),
     counterexamples: [
         {
             description: 'A mode at the band edge qa = π has ω = 2√(κ/m) on the lattice, while the ' +

package/dist/atlas/oscillators/bridges-coarse.js.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-coarse.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-coarse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AACpE,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAElD,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAElD,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B;;;;;GAKG;AACH,MAAM,UAAU,iBAAiB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAa,EAAE,CAAS;IAC9E,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACpE,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,mBAAmB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAa,EAAE,CAAS;IAChF,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtD,CAAC;AAED,+DAA+D;AAC/D,MAAM,UAAU,cAAc,CAAC,CAAS,EAAE,KAAa,EAAE,CAAS;IAChE,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;AAClC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,eAAe,CAAC,KAAa,EAAE,CAAS;IACtD,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;AAClC,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,4BAA4B,CAAC,EAAU;IACrD,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;AACxB,CAAC;AAED,oEAAoE;AACpE,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,EAAE,EAAE,eAAe;IACnB,QAAQ,EAAE,iBAAiB;IAC3B,QAAQ,EAAE,CAAC,aAAa,CAAC;IACzB,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,4CAA4C;IAC5D,SAAS,EAAE;QACT,oCAAoC;QACpC,yBAAyB;QACzB,WAAW;KACZ;IACD,eAAe,EAAE;QACf,oBAAoB;QACpB,0CAA0C;QAC1C,wDAAwD;KACzD;IACD,cAAc,EAAE,CAAC,yBAAyB,CAAC;IAC3C,MAAM,EAAE;QACN,MAAM,EAAE,MAAM;QACd,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC;QACnE,gBAAgB,EAAE,kBAAkB,CAClC,MAAM,EACN;YACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;YACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,eAAe,EAAE;YACvC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;YAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;SAC7B,EACD,CAAC,IAAI,CAAC,CACP;KACF;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,2EAA2E;gBAC3E,2EAA2E;gBAC3E,mBAAmB;YACrB,OAAO,EAAE,IAAI;SACd;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,wCAAwC;YAC9C,SAAS,EAAE,kEAAkE;SAC9E;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,wCAAwC;YAC9C,SAAS,EACP,kGAAkG;gBAClG,+GAA+G;SAClH;KACF;IACD,gFAAgF;IAChF,uFAAuF;IACvF,mFAAmF;IACnF,wFAAwF;IACxF,qBAAqB;IACrB,SAAS,EAAE;QACT,8IAA8I;QAC9I,0IAA0I;KAC3I;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;;;;;;;;GAcG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAqB,kBAAkB,CAAC,iBAAiB,CAAC,CAAC","sourcesContent":["/**\n * Bridge 5 of the oscillator pilot: the COARSE-GRAINING bridge from the\n * discrete mass–spring chain to the 1-D wave equation.\n *\n * This is the only Phase 0 bridge whose relation is `coarse-graining`: the\n * reduction map `u_n(t) ↦ u(x = na, t)` is not invertible, and the\n * short-wavelength half of the lattice Brillouin zone (`q > π/a`) has no image\n * in the continuum model. The information loss is therefore recorded in\n * `doesNotPreserve` and witnessed numerically by `W9`.\n *\n * @module atlas/oscillators/bridges-coarse\n */\n\nimport { LENGTH, MASS, VELOCITY } from '../../dimensional/types.js';\nimport { deriveRegimeGroups } from '../regime.js';\nimport type { AtlasBridge, RelationContract } from '../types.js';\nimport { relationContractOf } from './bridges-exact.js';\nimport { SPRING_CONSTANT } from './dimensions.js';\n\nconst FAMILY = 'oscillators';\n\n/**\n * The lattice dispersion relation of `m u_n″ = κ(u_{n+1} − 2u_n + u_{n−1})`:\n * `ω(q) = 2 √(κ/m) |sin(qa/2)|`.\n *\n * @internal\n */\nexport function latticeDispersion(q: number, a: number, kappa: number, m: number): number {\n  return 2 * Math.sqrt(kappa / m) * Math.abs(Math.sin((q * a) / 2));\n}\n\n/**\n * The continuum dispersion relation of `u_tt = c² u_xx` with `c² = κ a² / m`:\n * `ω = c q`.\n *\n * @internal\n */\nexport function continuumDispersion(q: number, a: number, kappa: number, m: number): number {\n  return Math.sqrt((kappa * a * a) / m) * Math.abs(q);\n}\n\n/** The long-wavelength wave speed `c = a √(κ/m)`. @internal */\nexport function chainWaveSpeed(a: number, kappa: number, m: number): number {\n  return a * Math.sqrt(kappa / m);\n}\n\n/**\n * The band edge of the lattice, at `qa = π`: `ω_max = 2 √(κ/m)`.\n * The continuum relation `ω = c q` is unbounded and has no counterpart —\n * that absence IS the information the coarse-graining discards.\n *\n * @internal\n */\nexport function latticeBandEdge(kappa: number, m: number): number {\n  return 2 * Math.sqrt(kappa / m);\n}\n\n/**\n * Leading-order relative dispersion error `1 − ω_lattice/ω_continuum`,\n * approximated as `(qa)²/24`.\n *\n * Convention note (fixed for the whole sprint): with `x = qa/2`,\n * `1 − sin x / x = x²/6 − x⁴/120 + …`, so the true relative error is\n * `(qa)²/24 − (qa)⁴/1920`. The quartic term makes `(qa)²/24` an\n * OVERESTIMATE; the signed deviation `(true − approx)/approx` is therefore\n * NEGATIVE and equals `−(qa)²/80` to leading order. The W9 assertion is on\n * the MAGNITUDE, so the sign convention does not change the test — it is\n * stated here so the next reader does not rederive it.\n *\n * @internal\n */\nexport function dispersionErrorApproximation(qa: number): number {\n  return (qa * qa) / 24;\n}\n\n/** Bridge 5: chain → wave equation by coarse-graining. @internal */\nexport const BRIDGE_CHAIN_WAVE: AtlasBridge = {\n  id: 'ab-chain-wave',\n  relation: 'coarse-graining',\n  premises: ['model-chain'],\n  conclusion: 'model-wave-1d',\n  transformation: 'u_n(t) ↦ u(x = n a, t), with c² = κ a² / m',\n  preserves: [\n    'long-wavelength dispersion ω ≈ c q',\n    'wave speed c = a √(κ/m)',\n    'linearity',\n  ],\n  doesNotPreserve: [\n    'modes with q > π/a',\n    'the band edge ω_max = 2 √(κ/m) at qa = π',\n    'the discrete lattice spacing a as an independent scale',\n  ],\n  sideConditions: ['long-wavelength, qa ≪ 1'],\n  regime: {\n    family: FAMILY,\n    inequalities: [{ group: 'qa', op: '<', bound: 1, alias: 'qa ≪ 1' }],\n    groupDefinitions: deriveRegimeGroups(\n      FAMILY,\n      [\n        { name: 'm', dim: MASS },\n        { name: 'kappa', dim: SPRING_CONSTANT },\n        { name: 'a', dim: LENGTH },\n        { name: 'c', dim: VELOCITY },\n      ],\n      ['qa'],\n    ),\n  },\n  counterexamples: [\n    {\n      description:\n        'A mode at the band edge qa = π has ω = 2√(κ/m) on the lattice, while the ' +\n        'continuum relation ω = c q is unbounded: the coarse-grained model has no ' +\n        'band edge at all.',\n      witness: 'W9',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W9',\n      kind: 'numeric',\n      test: 'tests/atlas/oscillators-coarse.test.ts',\n      tolerance: 'relative dispersion error matches (qa)²/24 within 0.5% of itself',\n    },\n    {\n      id: 'W9b',\n      kind: 'numeric',\n      test: 'tests/atlas/oscillators-coarse.test.ts',\n      tolerance:\n        'ring of 64 masses integrated: ω matches the lattice dispersion within 1e-9 and the coarse error ' +\n        '(qa)²/24 within 1%; superposition within 1e-10 (linearity); a 10% wrong κ and a cubic on-site force each fail',\n    },\n  ],\n  // The monatomic-chain dispersion omega(q) = 2 sqrt(kappa/m) |sin(qa/2)| and its\n  // long-wavelength continuum limit are standard solid-state results. Work and topic are\n  // cited rather than a precise section number, which would be fabrication in a more\n  // confident costume. This array was EMPTY, which is honest but leaves a textbook result\n  // looking unsourced.\n  citations: [\n    'Kittel, Introduction to Solid State Physics - Phonons I: the monatomic linear chain, its dispersion relation and the qa << 1 continuum limit',\n    'Ashcroft & Mermin, Solid State Physics - the one-dimensional monatomic Bravais lattice and the first Brillouin zone band edge at qa = pi',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge 5 re-registered through `RelationContract`.\n *\n * `coarse-graining` leaves `bound` OPTIONAL, and this record carries none: the\n * `qa²/24` dispersion residual is a per-mode series term, not a Lipschitz\n * constant with a uniform offset in a stated norm. The union permits the\n * absence rather than forcing an invented `K` and `delta`.\n *\n * No `deltaAt` either, for the same reason and by the same rule: `deltaAt` is\n * the machine form OF a `delta`, so a record with no `delta` has nothing for\n * it to be the machine form of. `admitApproximation` requires `deltaAt` only\n * of an `approximation`, so this record passes admission unchanged.\n *\n * @internal\n */\nexport const CONTRACT_CHAIN_WAVE: RelationContract = relationContractOf(BRIDGE_CHAIN_WAVE);\n"]}
+{"version":3,"file":"bridges-coarse.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-coarse.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AAEH,OAAO,EAAE,MAAM,EAAE,IAAI,EAAE,QAAQ,EAAE,MAAM,4BAA4B,CAAC;AACpE,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAElD,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B;;;;;GAKG;AACH,MAAM,UAAU,iBAAiB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAa,EAAE,CAAS;IAC9E,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACpE,CAAC;AAED;;;;;GAKG;AACH,MAAM,UAAU,mBAAmB,CAAC,CAAS,EAAE,CAAS,EAAE,KAAa,EAAE,CAAS;IAChF,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtD,CAAC;AAED,+DAA+D;AAC/D,MAAM,UAAU,cAAc,CAAC,CAAS,EAAE,KAAa,EAAE,CAAS;IAChE,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;AAClC,CAAC;AAED;;;;;;GAMG;AACH,MAAM,UAAU,eAAe,CAAC,KAAa,EAAE,CAAS;IACtD,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;AAClC,CAAC;AAED;;;;;;;;;;;;;GAaG;AACH,MAAM,UAAU,4BAA4B,CAAC,EAAU;IACrD,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;AACxB,CAAC;AAED,oEAAoE;AACpE,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,EAAE,EAAE,eAAe;IACnB,QAAQ,EAAE,iBAAiB;IAC3B,QAAQ,EAAE,CAAC,aAAa,CAAC;IACzB,UAAU,EAAE,eAAe;IAC3B,cAAc,EAAE,4CAA4C;IAC5D,SAAS,EAAE;QACT,oCAAoC;QACpC,yBAAyB;QACzB,WAAW;KACZ;IACD,eAAe,EAAE;QACf,oBAAoB;QACpB,0CAA0C;QAC1C,wDAAwD;KACzD;IACD,cAAc,EAAE,CAAC,yBAAyB,CAAC;IAC3C,MAAM,EAAE,UAAU,CAChB,MAAM,EACN;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,EAAE,eAAe,EAAE;QACvC,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE;QAC1B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,QAAQ,EAAE;KAC7B,EACD,CAAC,IAAI,CAAC,EACN,CAAC,EAAE,KAAK,EAAE,IAAI,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC,CACtD;IACD,eAAe,EAAE;QACf;YACE,WAAW,EACT,2EAA2E;gBAC3E,2EAA2E;gBAC3E,mBAAmB;YACrB,OAAO,EAAE,IAAI;SACd;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,wCAAwC;YAC9C,SAAS,EAAE,kEAAkE;SAC9E;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,wCAAwC;YAC9C,SAAS,EACP,kGAAkG;gBAClG,+GAA+G;SAClH;KACF;IACD,gFAAgF;IAChF,uFAAuF;IACvF,mFAAmF;IACnF,wFAAwF;IACxF,qBAAqB;IACrB,SAAS,EAAE;QACT,8IAA8I;QAC9I,0IAA0I;KAC3I;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF;;;;;;;;;;;;;;GAcG;AACH,MAAM,CAAC,MAAM,mBAAmB,GAAqB,kBAAkB,CAAC,iBAAiB,CAAC,CAAC","sourcesContent":["/**\n * Bridge 5 of the oscillator pilot: the COARSE-GRAINING bridge from the\n * discrete mass–spring chain to the 1-D wave equation.\n *\n * This is the only Phase 0 bridge whose relation is `coarse-graining`: the\n * reduction map `u_n(t) ↦ u(x = na, t)` is not invertible, and the\n * short-wavelength half of the lattice Brillouin zone (`q > π/a`) has no image\n * in the continuum model. The information loss is therefore recorded in\n * `doesNotPreserve` and witnessed numerically by `W9`.\n *\n * @module atlas/oscillators/bridges-coarse\n */\n\nimport { LENGTH, MASS, VELOCITY } from '../../dimensional/types.js';\nimport { makeRegime } from '../regime.js';\nimport type { AtlasBridge, RelationContract } from '../types.js';\nimport { relationContractOf } from './bridges-exact.js';\nimport { SPRING_CONSTANT } from './dimensions.js';\n\nconst FAMILY = 'oscillators';\n\n/**\n * The lattice dispersion relation of `m u_n″ = κ(u_{n+1} − 2u_n + u_{n−1})`:\n * `ω(q) = 2 √(κ/m) |sin(qa/2)|`.\n *\n * @internal\n */\nexport function latticeDispersion(q: number, a: number, kappa: number, m: number): number {\n  return 2 * Math.sqrt(kappa / m) * Math.abs(Math.sin((q * a) / 2));\n}\n\n/**\n * The continuum dispersion relation of `u_tt = c² u_xx` with `c² = κ a² / m`:\n * `ω = c q`.\n *\n * @internal\n */\nexport function continuumDispersion(q: number, a: number, kappa: number, m: number): number {\n  return Math.sqrt((kappa * a * a) / m) * Math.abs(q);\n}\n\n/** The long-wavelength wave speed `c = a √(κ/m)`. @internal */\nexport function chainWaveSpeed(a: number, kappa: number, m: number): number {\n  return a * Math.sqrt(kappa / m);\n}\n\n/**\n * The band edge of the lattice, at `qa = π`: `ω_max = 2 √(κ/m)`.\n * The continuum relation `ω = c q` is unbounded and has no counterpart —\n * that absence IS the information the coarse-graining discards.\n *\n * @internal\n */\nexport function latticeBandEdge(kappa: number, m: number): number {\n  return 2 * Math.sqrt(kappa / m);\n}\n\n/**\n * Leading-order relative dispersion error `1 − ω_lattice/ω_continuum`,\n * approximated as `(qa)²/24`.\n *\n * Convention note (fixed for the whole sprint): with `x = qa/2`,\n * `1 − sin x / x = x²/6 − x⁴/120 + …`, so the true relative error is\n * `(qa)²/24 − (qa)⁴/1920`. The quartic term makes `(qa)²/24` an\n * OVERESTIMATE; the signed deviation `(true − approx)/approx` is therefore\n * NEGATIVE and equals `−(qa)²/80` to leading order. The W9 assertion is on\n * the MAGNITUDE, so the sign convention does not change the test — it is\n * stated here so the next reader does not rederive it.\n *\n * @internal\n */\nexport function dispersionErrorApproximation(qa: number): number {\n  return (qa * qa) / 24;\n}\n\n/** Bridge 5: chain → wave equation by coarse-graining. @internal */\nexport const BRIDGE_CHAIN_WAVE: AtlasBridge = {\n  id: 'ab-chain-wave',\n  relation: 'coarse-graining',\n  premises: ['model-chain'],\n  conclusion: 'model-wave-1d',\n  transformation: 'u_n(t) ↦ u(x = n a, t), with c² = κ a² / m',\n  preserves: [\n    'long-wavelength dispersion ω ≈ c q',\n    'wave speed c = a √(κ/m)',\n    'linearity',\n  ],\n  doesNotPreserve: [\n    'modes with q > π/a',\n    'the band edge ω_max = 2 √(κ/m) at qa = π',\n    'the discrete lattice spacing a as an independent scale',\n  ],\n  sideConditions: ['long-wavelength, qa ≪ 1'],\n  regime: makeRegime(\n    FAMILY,\n    [\n      { name: 'm', dim: MASS },\n      { name: 'kappa', dim: SPRING_CONSTANT },\n      { name: 'a', dim: LENGTH },\n      { name: 'c', dim: VELOCITY },\n    ],\n    ['qa'],\n    [{ group: 'qa', op: '<', bound: 1, alias: 'qa ≪ 1' }],\n  ),\n  counterexamples: [\n    {\n      description:\n        'A mode at the band edge qa = π has ω = 2√(κ/m) on the lattice, while the ' +\n        'continuum relation ω = c q is unbounded: the coarse-grained model has no ' +\n        'band edge at all.',\n      witness: 'W9',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W9',\n      kind: 'numeric',\n      test: 'tests/atlas/oscillators-coarse.test.ts',\n      tolerance: 'relative dispersion error matches (qa)²/24 within 0.5% of itself',\n    },\n    {\n      id: 'W9b',\n      kind: 'numeric',\n      test: 'tests/atlas/oscillators-coarse.test.ts',\n      tolerance:\n        'ring of 64 masses integrated: ω matches the lattice dispersion within 1e-9 and the coarse error ' +\n        '(qa)²/24 within 1%; superposition within 1e-10 (linearity); a 10% wrong κ and a cubic on-site force each fail',\n    },\n  ],\n  // The monatomic-chain dispersion omega(q) = 2 sqrt(kappa/m) |sin(qa/2)| and its\n  // long-wavelength continuum limit are standard solid-state results. Work and topic are\n  // cited rather than a precise section number, which would be fabrication in a more\n  // confident costume. This array was EMPTY, which is honest but leaves a textbook result\n  // looking unsourced.\n  citations: [\n    'Kittel, Introduction to Solid State Physics - Phonons I: the monatomic linear chain, its dispersion relation and the qa << 1 continuum limit',\n    'Ashcroft & Mermin, Solid State Physics - the one-dimensional monatomic Bravais lattice and the first Brillouin zone band edge at qa = pi',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/**\n * Bridge 5 re-registered through `RelationContract`.\n *\n * `coarse-graining` leaves `bound` OPTIONAL, and this record carries none: the\n * `qa²/24` dispersion residual is a per-mode series term, not a Lipschitz\n * constant with a uniform offset in a stated norm. The union permits the\n * absence rather than forcing an invented `K` and `delta`.\n *\n * No `deltaAt` either, for the same reason and by the same rule: `deltaAt` is\n * the machine form OF a `delta`, so a record with no `delta` has nothing for\n * it to be the machine form of. `admitApproximation` requires `deltaAt` only\n * of an `approximation`, so this record passes admission unchanged.\n *\n * @internal\n */\nexport const CONTRACT_CHAIN_WAVE: RelationContract = relationContractOf(BRIDGE_CHAIN_WAVE);\n"]}

package/dist/atlas/oscillators/bridges-exact.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-exact.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-exact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAKH,OAAO,KAAK,EAAE,WAAW,EAAU,gBAAgB,EAAE,MAAM,aAAa,CAAC;AA8CzE;;;;;;;;;;;;GAYG;AACH,eAAO,MAAM,gBAAgB,EAAE,WA6D9B,CAAC;AAEF;;;;;;;;;GASG;AACH,eAAO,MAAM,iBAAiB,EAAE,WAoC/B,CAAC;AAMF;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,WAAW,GAAG,gBAAgB,CAwBxE;AAED,qEAAqE;AACrE,eAAO,MAAM,kBAAkB,EAAE,gBAAuD,CAAC;AAEzF,qEAAqE;AACrE,eAAO,MAAM,mBAAmB,EAAE,gBAAwD,CAAC"}
+{"version":3,"file":"bridges-exact.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-exact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAKH,OAAO,KAAK,EAAE,WAAW,EAAU,gBAAgB,EAAE,MAAM,aAAa,CAAC;AA6CzE;;;;;;;;;;;;GAYG;AACH,eAAO,MAAM,gBAAgB,EAAE,WA6D9B,CAAC;AAEF;;;;;;;;;GASG;AACH,eAAO,MAAM,iBAAiB,EAAE,WAoC/B,CAAC;AAMF;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,kBAAkB,CAAC,MAAM,EAAE,WAAW,GAAG,gBAAgB,CAwBxE;AAED,qEAAqE;AACrE,eAAO,MAAM,kBAAkB,EAAE,gBAAuD,CAAC;AAEzF,qEAAqE;AACrE,eAAO,MAAM,mBAAmB,EAAE,gBAAwD,CAAC"}

package/dist/atlas/oscillators/bridges-exact.js

@@ -21,7 +21,7 @@
  * @module atlas/oscillators/bridges-exact
  */
 import { physjsFormalRef } from '../physjs-ref.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { getAtlasModel } from './models.js';
 import { SPRING_LC_RELATIVE_PERIOD_TRANSPORT } from './norm-transport.js';
 const FAMILY = 'oscillators';
@@ -40,7 +40,10 @@ const WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';
  *
  * Work and topic are cited rather than precise section numbers, because a section number
  * stated without checking the edition in hand is the same fabrication in a more confident
- * costume. `reviewStatus` stays 'proposed' until the independent physicist review lands.
+ * costume. `reviewStatus` is the record's declared review state and derives no evidence
+ * tag (`deriveEvidence` never reads it). It is 'proposed' here because no independent
+ * physicist review of these records exists; the atlas no longer gates on one (the ruling
+ * is recorded in `ACTIVE.md` and `NOTES.md`).
  */
 const CITATIONS = [
     'Olson, Dynamical Analogies - the force-voltage (impedance) analogy between a mass-spring–damper and a series RLC circuit',
@@ -57,11 +60,7 @@ function bridgeRegime(premiseId, conclusionId) {
         ...getAtlasModel(premiseId).parameters,
         ...getAtlasModel(conclusionId).parameters,
     ];
-    return {
-        family: FAMILY,
-        inequalities: [],
-        groupDefinitions: deriveRegimeGroups(FAMILY, parameters, []),
-    };
+    return makeRegime(FAMILY, parameters, [], []);
 }
 /**
  * Bridge 1 — the lossless spring ↔ LC dictionary.

package/dist/atlas/oscillators/bridges-exact.js.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-exact.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-exact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClD,OAAO,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAE5C,OAAO,EAAE,mCAAmC,EAAE,MAAM,qBAAqB,CAAC;AAE1E,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B,yDAAyD;AACzD,MAAM,IAAI,GAAG,uCAAuC,CAAC;AACrD,iFAAiF;AACjF,MAAM,oBAAoB,GAAG,qCAAqC,CAAC;AAEnE;;;;;;;;;;;;GAYG;AACH,MAAM,SAAS,GAAsB;IACnC,0HAA0H;IAC1H,+GAA+G;CAChH,CAAC;AAEF;;;;;GAKG;AACH,SAAS,YAAY,CAAC,SAAiB,EAAE,YAAoB;IAC3D,MAAM,UAAU,GAAG;QACjB,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC,UAAU;QACtC,GAAG,aAAa,CAAC,YAAY,CAAC,CAAC,UAAU;KAC1C,CAAC;IACF,OAAO;QACL,MAAM,EAAE,MAAM;QACd,YAAY,EAAE,EAAE;QAChB,gBAAgB,EAAE,kBAAkB,CAAC,MAAM,EAAE,UAAU,EAAE,EAAE,CAAC;KAC7D,CAAC;AACJ,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAgB;IAC3C,EAAE,EAAE,cAAc;IAClB,QAAQ,EAAE,mBAAmB;IAC7B,QAAQ,EAAE,CAAC,cAAc,CAAC;IAC1B,UAAU,EAAE,UAAU;IACtB,mFAAmF;IACnF,oFAAoF;IACpF,kEAAkE;IAClE,cAAc,EACZ,uEAAuE;QACvE,sEAAsE;IACxE,OAAO,EAAE,oBAAoB;IAC7B,SAAS,EAAE,CAAC,0DAA0D,EAAE,oBAAoB,EAAE,gBAAgB,CAAC;IAC/G,eAAe,EAAE,CAAC,yBAAyB,EAAE,OAAO,CAAC;IACrD,cAAc,EAAE,CAAC,gBAAgB,EAAE,UAAU,EAAE,UAAU,EAAE,gBAAgB,CAAC;IAC5E,MAAM,EAAE,YAAY,CAAC,cAAc,EAAE,UAAU,CAAC;IAChD,eAAe,EAAE;QACf;YACE,WAAW,EACT,6FAA6F;gBAC7F,6BAA6B;gBAC7B,yFAAyF;YAC3F,OAAO,EAAE,KAAK;SACf;KACF;IACD,2EAA2E;IAC3E,4EAA4E;IAC5E,4EAA4E;IAC5E,sCAAsC;IACtC,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,yDAAyD;SACrE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,qDAAqD;SACjE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,kFAAkF;SAC9F;QACD,EAAE,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,EAAE,SAAS,EAAE,4BAA4B,EAAE;QACnF;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,oEAAoE;SAChF;KACF;IACD,SAAS,EAAE,SAAS;IACpB,YAAY,EAAE,UAAU;IACxB,+FAA+F;IAC/F,cAAc,EAAE,CAAC,mCAAmC,CAAC;IACrD,SAAS,EAAE,eAAe,CAAC,cAAc,CAAC;CAC3C,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,EAAE,EAAE,eAAe;IACnB,QAAQ,EAAE,mBAAmB;IAC7B,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,WAAW;IACvB,cAAc,EACZ,sDAAsD;QACtD,0EAA0E;IAC5E,OAAO,EAAE,oBAAoB;IAC7B,SAAS,EAAE,CAAC,0DAA0D,EAAE,eAAe,EAAE,gBAAgB,CAAC;IAC1G,eAAe,EAAE,CAAC,yBAAyB,EAAE,OAAO,CAAC;IACrD,cAAc,EAAE;QACd,gFAAgF;QAChF,sBAAsB;QACtB,UAAU;QACV,gBAAgB;KACjB;IACD,MAAM,EAAE,YAAY,CAAC,qBAAqB,EAAE,WAAW,CAAC;IACxD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,sDAAsD;SAClE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,kFAAkF;SAC9F;KACF;IACD,SAAS,EAAE,SAAS;IACpB,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,eAAe,CAAC;CAC5C,CAAC;AAEF,gFAAgF;AAChF,2EAA2E;AAC3E,gFAAgF;AAEhF;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAmB;IACpD,MAAM,EAAE,QAAQ,EAAE,cAAc,EAAE,GAAG,MAAM,CAAC;IAC5C,QAAQ,QAAQ,EAAE,CAAC;QACjB,KAAK,mBAAmB,EAAE,CAAC;YACzB,IAAI,MAAM,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBACjC,MAAM,IAAI,SAAS,CAAC,GAAG,MAAM,CAAC,EAAE,wCAAwC,CAAC,CAAC;YAC5E,CAAC;YACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,OAAO,EAAE,MAAM,CAAC,OAAO,EAAE,CAAC;QACrE,CAAC;QACD,KAAK,eAAe,EAAE,CAAC;YACrB,IAAI,MAAM,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;gBAC/B,MAAM,IAAI,SAAS,CAAC,GAAG,MAAM,CAAC,EAAE,iCAAiC,CAAC,CAAC;YACrE,CAAC;YACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,CAAC;QACjE,CAAC;QACD,KAAK,iBAAiB;YACpB,wEAAwE;YACxE,yDAAyD;YACzD,OAAO,MAAM,CAAC,KAAK,KAAK,SAAS;gBAC/B,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE;gBACpC,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,CAAC;QAC9D;YACE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,CAAC;IAC9C,CAAC;AACH,CAAC;AAED,qEAAqE;AACrE,MAAM,CAAC,MAAM,kBAAkB,GAAqB,kBAAkB,CAAC,gBAAgB,CAAC,CAAC;AAEzF,qEAAqE;AACrE,MAAM,CAAC,MAAM,mBAAmB,GAAqB,kBAAkB,CAAC,iBAAiB,CAAC,CAAC","sourcesContent":["/**\n * The two Phase 0 `exact-equivalence` bridges of design note §3.\n *\n * - `ab-spring-lc` — the lossless mechanical/electrical dictionary, exact with\n *   no side condition beyond positivity and losslessness.\n * - `ab-damped-rlc` — exact only ON a side condition: the two damping ratios\n *   must agree, `b/(2√(mk)) = (R/2)√(C/L)`.\n *\n * **Regimes carry no inequalities.** `RegimeInequality` expresses `<`, `<=`,\n * `>`, `>=` only, so the bridge-2 side condition — an EQUALITY between two\n * derived ratios — cannot be written as one. It is carried in\n * `sideConditions` as prose instead, and the equality is what witness W2\n * executes. Writing an inequality that merely resembles it would be a\n * fabricated constraint, so none is written. `groupDefinitions` is derived\n * from the premise and conclusion parameters by `deriveRegimeGroups`, exactly\n * as the models do; nothing there is hand-authored.\n *\n * **Evidence tags are witness-gated.** A tag appears only because a witness in\n * `tests/atlas/oscillators-exact.test.ts` passes and shows it.\n *\n * @module atlas/oscillators/bridges-exact\n */\n\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport { deriveRegimeGroups } from '../regime.js';\nimport { getAtlasModel } from './models.js';\nimport type { AtlasBridge, Regime, RelationContract } from '../types.js';\nimport { SPRING_LC_RELATIVE_PERIOD_TRANSPORT } from './norm-transport.js';\n\nconst FAMILY = 'oscillators';\n\n/** The witness test file every witness below runs in. */\nconst TEST = 'tests/atlas/oscillators-exact.test.ts';\n/** Where the CAS witnesses W1s and W2s are named; the artifact pin runs them. */\nconst WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';\n\n/**\n * Literature for the electromechanical analogy, NOT a pointer at our own plan.\n *\n * The previous value here was `docs/planning/Atlas-Phase-0-Design.md §3, §6` - a citation\n * naming the document that ASSERTED the claim. That is circular, and it is worse than an\n * empty array because it LOOKS sourced and would satisfy any check that only asks whether\n * `citations[]` is non-empty. Same failure class as the evidence-tag theatre that\n * `tests/atlas/evidence-rule.test.ts` exists to catch.\n *\n * Work and topic are cited rather than precise section numbers, because a section number\n * stated without checking the edition in hand is the same fabrication in a more confident\n * costume. `reviewStatus` stays 'proposed' until the independent physicist review lands.\n */\nconst CITATIONS: readonly string[] = [\n  'Olson, Dynamical Analogies - the force-voltage (impedance) analogy between a mass-spring–damper and a series RLC circuit',\n  'Feynman, Lectures on Physics Vol. II, chapter on resonance - the LC oscillator and its mechanical counterpart',\n];\n\n/**\n * Regime coordinates spanning a bridge's premise and conclusion parameters.\n *\n * The two models' parameter names are disjoint in both bridges, so the union\n * is a legal `deriveRegimeGroups` input.\n */\nfunction bridgeRegime(premiseId: string, conclusionId: string): Regime {\n  const parameters = [\n    ...getAtlasModel(premiseId).parameters,\n    ...getAtlasModel(conclusionId).parameters,\n  ];\n  return {\n    family: FAMILY,\n    inequalities: [],\n    groupDefinitions: deriveRegimeGroups(FAMILY, parameters, []),\n  };\n}\n\n/**\n * Bridge 1 — the lossless spring ↔ LC dictionary.\n *\n * `m ↔ L`, `k ↔ 1/C`, `x ↔ q`. Under `u = x/x0` (or `q/q0`) and `τ = ω0 t`\n * both systems become `u'' + u = 0`, so the equivalence is exact and\n * invertible.\n *\n * Its counterexample is W2b, which lives on THIS bridge: it shows that adding\n * a resistor to the LC side breaks the equivalence, so bridge 1 is a statement\n * about the LOSSLESS pair only.\n *\n * @internal\n */\nexport const BRIDGE_SPRING_LC: AtlasBridge = {\n  id: 'ab-spring-lc',\n  relation: 'exact-equivalence',\n  premises: ['model-spring'],\n  conclusion: 'model-lc',\n  // Two nondimensionalizations, each with its own ω0; composed, the spring→LC map is\n  // q(t) = (q0/x0)·x(ω_LC t / ω_s). The dimensioned frequencies differ (W1a: ω_s = 2,\n  // ω_LC = 2√2); what is preserved is the frequency in units of ω0.\n  transformation:\n    'force–voltage analogy m ↔ L, k ↔ 1/C, x ↔ q (so ω0² = k/m ↔ 1/(LC)); ' +\n    'u = x/x0 or q/q0, τ = ω0 t; composed, q(t) = (q0/x0)·x(ω_LC t / ω_s)',\n  inverse: 'x = x0 u, t = τ/ω0',\n  preserves: ['the natural frequency in units of ω0 (ω = 1 in τ = ω0 t)', 'energy up to scale', 'phase portrait'],\n  doesNotPreserve: ['physical interpretation', 'units'],\n  sideConditions: ['m, k, L, C > 0', 'lossless', 'unforced', 'x0, q0 nonzero'],\n  regime: bridgeRegime('model-spring', 'model-lc'),\n  counterexamples: [\n    {\n      description:\n        'adding R to bridge 1 breaks it: the circuit L = 2, C = 0.125 of witness W2b with R = 4 has ' +\n        'ζ_RLC = (R/2)√(C/L) = 0.5, ' +\n        'which no lossless spring matches — the trajectories separate by more than 1e-2 at τ = π',\n      witness: 'W2b',\n    },\n  ],\n  // `symbolically-checked` is NOT stored here any more (Phase 4 S4.3): it is\n  // DERIVED from `data/atlas/witness-results.json`, where W1s records the CAS\n  // check of this bridge's dictionary. A stored copy would be a second source\n  // of truth that no run could falsify.\n  witnesses: [\n    {\n      id: 'W1',\n      kind: 'symbolic',\n      test: TEST,\n      tolerance: 'exact — rational exponent arithmetic, no floating point',\n    },\n    {\n      id: 'W1a',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '|u_spring − u_lc| < 1e-8; each within 1e-8 of cos τ',\n    },\n    {\n      id: 'W1b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'inverse-mapped trajectory solves m x″ + k x = 0 within 1e-8; a 1% wrong ω0 fails',\n    },\n    { id: 'W2b', kind: 'numeric', test: TEST, tolerance: 'separation > 1e-2 at τ = π' },\n    {\n      id: 'W1s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: k/m under m ↔ L, k ↔ 1/C minus 1/(LC) simplifies to literal 0',\n    },\n  ],\n  citations: CITATIONS,\n  reviewStatus: 'proposed',\n  // spring → LC in relative period error only; see ./norm-transport.ts for what is not declared.\n  normTransports: [SPRING_LC_RELATIVE_PERIOD_TRANSPORT],\n  formalRef: physjsFormalRef('ab-spring-lc'),\n};\n\n/**\n * Bridge 2 — damped spring ↔ RLC, exact ON a side condition.\n *\n * The undamped dictionary of bridge 1 does not extend for free: `b` and `R`\n * are independent parameters, and the nondimensional systems coincide only\n * when the damping ratios match. The side condition is that equality, and W2\n * DERIVES `R` from it rather than assuming a value.\n *\n * @internal\n */\nexport const BRIDGE_DAMPED_RLC: AtlasBridge = {\n  id: 'ab-damped-rlc',\n  relation: 'exact-equivalence',\n  premises: ['model-damped-spring'],\n  conclusion: 'model-rlc',\n  transformation:\n    'force–voltage analogy m ↔ L, k ↔ 1/C, b ↔ R, x ↔ q; ' +\n    'u = x/x0 or q/q0, τ = ω0 t; composed, q(t) = (q0/x0)·x(ω_RLC t / ω_mech)',\n  inverse: 'x = x0 u, t = τ/ω0',\n  preserves: ['the natural frequency in units of ω0 (ω = 1 in τ = ω0 t)', 'damping ratio', 'phase portrait'],\n  doesNotPreserve: ['physical interpretation', 'units'],\n  sideConditions: [\n    'b/√(mk) = R√(C/L), equivalently ζ_mech = b/(2√(mk)) equals ζ_RLC = (R/2)√(C/L)',\n    'm, b, k, L, R, C > 0',\n    'unforced',\n    'x0, q0 nonzero',\n  ],\n  regime: bridgeRegime('model-damped-spring', 'model-rlc'),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'W2',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'R derived exactly; |u_mech − u_rlc| < 1e-8 at four τ',\n    },\n    {\n      id: 'W2s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: b²/(4mk) under m ↔ L, k ↔ 1/C, b ↔ R minus R²C/(4L) simplifies to literal 0',\n    },\n  ],\n  citations: CITATIONS,\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-damped-rlc'),\n};\n\n// ─────────────────────────────────────────────────────────────────────────────\n// Phase 1 re-registration: the Sprint 0 bridges through `RelationContract`\n// ─────────────────────────────────────────────────────────────────────────────\n\n/**\n * Re-register an `AtlasBridge` through the Phase 1 `RelationContract` union.\n *\n * `AtlasBridge` states its two conditional requirements — `exact-equivalence`\n * needs `inverse`, `approximation` needs `bound` — only in doc comments, so a\n * record that breaks them type-checks. `RelationContract` states them in the\n * type system. This function is the one place the Sprint 0 records cross from\n * the weaker shape to the stronger one, and it THROWS on the two cases the\n * doc comments could only describe.\n *\n * It copies rather than restates: every field comes off the bridge, so a\n * contract cannot drift from the record it re-registers.\n *\n * @throws TypeError when the bridge omits the field its relation requires.\n * @internal\n */\nexport function relationContractOf(bridge: AtlasBridge): RelationContract {\n  const { relation, transformation } = bridge;\n  switch (relation) {\n    case 'exact-equivalence': {\n      if (bridge.inverse === undefined) {\n        throw new TypeError(`${bridge.id}: exact-equivalence without an inverse`);\n      }\n      return { type: relation, transformation, inverse: bridge.inverse };\n    }\n    case 'approximation': {\n      if (bridge.bound === undefined) {\n        throw new TypeError(`${bridge.id}: approximation without a bound`);\n      }\n      return { type: relation, transformation, bound: bridge.bound };\n    }\n    case 'coarse-graining':\n      // `bound` is OPTIONAL here: the chain/wave bridge carries none, and the\n      // union allows that rather than forcing an invented one.\n      return bridge.bound === undefined\n        ? { type: relation, transformation }\n        : { type: relation, transformation, bound: bridge.bound };\n    default:\n      return { type: relation, transformation };\n  }\n}\n\n/** Bridge 1 re-registered through the Phase 1 contract. @internal */\nexport const CONTRACT_SPRING_LC: RelationContract = relationContractOf(BRIDGE_SPRING_LC);\n\n/** Bridge 2 re-registered through the Phase 1 contract. @internal */\nexport const CONTRACT_DAMPED_RLC: RelationContract = relationContractOf(BRIDGE_DAMPED_RLC);\n"]}
+{"version":3,"file":"bridges-exact.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-exact.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,aAAa,EAAE,MAAM,aAAa,CAAC;AAE5C,OAAO,EAAE,mCAAmC,EAAE,MAAM,qBAAqB,CAAC;AAE1E,MAAM,MAAM,GAAG,aAAa,CAAC;AAE7B,yDAAyD;AACzD,MAAM,IAAI,GAAG,uCAAuC,CAAC;AACrD,iFAAiF;AACjF,MAAM,oBAAoB,GAAG,qCAAqC,CAAC;AAEnE;;;;;;;;;;;;;;;GAeG;AACH,MAAM,SAAS,GAAsB;IACnC,0HAA0H;IAC1H,+GAA+G;CAChH,CAAC;AAEF;;;;;GAKG;AACH,SAAS,YAAY,CAAC,SAAiB,EAAE,YAAoB;IAC3D,MAAM,UAAU,GAAG;QACjB,GAAG,aAAa,CAAC,SAAS,CAAC,CAAC,UAAU;QACtC,GAAG,aAAa,CAAC,YAAY,CAAC,CAAC,UAAU;KAC1C,CAAC;IACF,OAAO,UAAU,CAAC,MAAM,EAAE,UAAU,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;AAChD,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,CAAC,MAAM,gBAAgB,GAAgB;IAC3C,EAAE,EAAE,cAAc;IAClB,QAAQ,EAAE,mBAAmB;IAC7B,QAAQ,EAAE,CAAC,cAAc,CAAC;IAC1B,UAAU,EAAE,UAAU;IACtB,mFAAmF;IACnF,oFAAoF;IACpF,kEAAkE;IAClE,cAAc,EACZ,uEAAuE;QACvE,sEAAsE;IACxE,OAAO,EAAE,oBAAoB;IAC7B,SAAS,EAAE,CAAC,0DAA0D,EAAE,oBAAoB,EAAE,gBAAgB,CAAC;IAC/G,eAAe,EAAE,CAAC,yBAAyB,EAAE,OAAO,CAAC;IACrD,cAAc,EAAE,CAAC,gBAAgB,EAAE,UAAU,EAAE,UAAU,EAAE,gBAAgB,CAAC;IAC5E,MAAM,EAAE,YAAY,CAAC,cAAc,EAAE,UAAU,CAAC;IAChD,eAAe,EAAE;QACf;YACE,WAAW,EACT,6FAA6F;gBAC7F,6BAA6B;gBAC7B,yFAAyF;YAC3F,OAAO,EAAE,KAAK;SACf;KACF;IACD,2EAA2E;IAC3E,4EAA4E;IAC5E,4EAA4E;IAC5E,sCAAsC;IACtC,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,yDAAyD;SACrE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,qDAAqD;SACjE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,kFAAkF;SAC9F;QACD,EAAE,EAAE,EAAE,KAAK,EAAE,IAAI,EAAE,SAAS,EAAE,IAAI,EAAE,IAAI,EAAE,SAAS,EAAE,4BAA4B,EAAE;QACnF;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,oEAAoE;SAChF;KACF;IACD,SAAS,EAAE,SAAS;IACpB,YAAY,EAAE,UAAU;IACxB,+FAA+F;IAC/F,cAAc,EAAE,CAAC,mCAAmC,CAAC;IACrD,SAAS,EAAE,eAAe,CAAC,cAAc,CAAC;CAC3C,CAAC;AAEF;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,iBAAiB,GAAgB;IAC5C,EAAE,EAAE,eAAe;IACnB,QAAQ,EAAE,mBAAmB;IAC7B,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,WAAW;IACvB,cAAc,EACZ,sDAAsD;QACtD,0EAA0E;IAC5E,OAAO,EAAE,oBAAoB;IAC7B,SAAS,EAAE,CAAC,0DAA0D,EAAE,eAAe,EAAE,gBAAgB,CAAC;IAC1G,eAAe,EAAE,CAAC,yBAAyB,EAAE,OAAO,CAAC;IACrD,cAAc,EAAE;QACd,gFAAgF;QAChF,sBAAsB;QACtB,UAAU;QACV,gBAAgB;KACjB;IACD,MAAM,EAAE,YAAY,CAAC,qBAAqB,EAAE,WAAW,CAAC;IACxD,eAAe,EAAE,EAAE;IACnB,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,sDAAsD;SAClE;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,UAAU;YAChB,IAAI,EAAE,oBAAoB;YAC1B,SAAS,EAAE,kFAAkF;SAC9F;KACF;IACD,SAAS,EAAE,SAAS;IACpB,YAAY,EAAE,UAAU;IACxB,SAAS,EAAE,eAAe,CAAC,eAAe,CAAC;CAC5C,CAAC;AAEF,gFAAgF;AAChF,2EAA2E;AAC3E,gFAAgF;AAEhF;;;;;;;;;;;;;;;GAeG;AACH,MAAM,UAAU,kBAAkB,CAAC,MAAmB;IACpD,MAAM,EAAE,QAAQ,EAAE,cAAc,EAAE,GAAG,MAAM,CAAC;IAC5C,QAAQ,QAAQ,EAAE,CAAC;QACjB,KAAK,mBAAmB,EAAE,CAAC;YACzB,IAAI,MAAM,CAAC,OAAO,KAAK,SAAS,EAAE,CAAC;gBACjC,MAAM,IAAI,SAAS,CAAC,GAAG,MAAM,CAAC,EAAE,wCAAwC,CAAC,CAAC;YAC5E,CAAC;YACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,OAAO,EAAE,MAAM,CAAC,OAAO,EAAE,CAAC;QACrE,CAAC;QACD,KAAK,eAAe,EAAE,CAAC;YACrB,IAAI,MAAM,CAAC,KAAK,KAAK,SAAS,EAAE,CAAC;gBAC/B,MAAM,IAAI,SAAS,CAAC,GAAG,MAAM,CAAC,EAAE,iCAAiC,CAAC,CAAC;YACrE,CAAC;YACD,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,CAAC;QACjE,CAAC;QACD,KAAK,iBAAiB;YACpB,wEAAwE;YACxE,yDAAyD;YACzD,OAAO,MAAM,CAAC,KAAK,KAAK,SAAS;gBAC/B,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE;gBACpC,CAAC,CAAC,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,CAAC;QAC9D;YACE,OAAO,EAAE,IAAI,EAAE,QAAQ,EAAE,cAAc,EAAE,CAAC;IAC9C,CAAC;AACH,CAAC;AAED,qEAAqE;AACrE,MAAM,CAAC,MAAM,kBAAkB,GAAqB,kBAAkB,CAAC,gBAAgB,CAAC,CAAC;AAEzF,qEAAqE;AACrE,MAAM,CAAC,MAAM,mBAAmB,GAAqB,kBAAkB,CAAC,iBAAiB,CAAC,CAAC","sourcesContent":["/**\n * The two Phase 0 `exact-equivalence` bridges of design note §3.\n *\n * - `ab-spring-lc` — the lossless mechanical/electrical dictionary, exact with\n *   no side condition beyond positivity and losslessness.\n * - `ab-damped-rlc` — exact only ON a side condition: the two damping ratios\n *   must agree, `b/(2√(mk)) = (R/2)√(C/L)`.\n *\n * **Regimes carry no inequalities.** `RegimeInequality` expresses `<`, `<=`,\n * `>`, `>=` only, so the bridge-2 side condition — an EQUALITY between two\n * derived ratios — cannot be written as one. It is carried in\n * `sideConditions` as prose instead, and the equality is what witness W2\n * executes. Writing an inequality that merely resembles it would be a\n * fabricated constraint, so none is written. `groupDefinitions` is derived\n * from the premise and conclusion parameters by `deriveRegimeGroups`, exactly\n * as the models do; nothing there is hand-authored.\n *\n * **Evidence tags are witness-gated.** A tag appears only because a witness in\n * `tests/atlas/oscillators-exact.test.ts` passes and shows it.\n *\n * @module atlas/oscillators/bridges-exact\n */\n\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport { makeRegime } from '../regime.js';\nimport { getAtlasModel } from './models.js';\nimport type { AtlasBridge, Regime, RelationContract } from '../types.js';\nimport { SPRING_LC_RELATIVE_PERIOD_TRANSPORT } from './norm-transport.js';\n\nconst FAMILY = 'oscillators';\n\n/** The witness test file every witness below runs in. */\nconst TEST = 'tests/atlas/oscillators-exact.test.ts';\n/** Where the CAS witnesses W1s and W2s are named; the artifact pin runs them. */\nconst WITNESS_RESULTS_TEST = 'tests/atlas/witness-results.test.ts';\n\n/**\n * Literature for the electromechanical analogy, NOT a pointer at our own plan.\n *\n * The previous value here was `docs/planning/Atlas-Phase-0-Design.md §3, §6` - a citation\n * naming the document that ASSERTED the claim. That is circular, and it is worse than an\n * empty array because it LOOKS sourced and would satisfy any check that only asks whether\n * `citations[]` is non-empty. Same failure class as the evidence-tag theatre that\n * `tests/atlas/evidence-rule.test.ts` exists to catch.\n *\n * Work and topic are cited rather than precise section numbers, because a section number\n * stated without checking the edition in hand is the same fabrication in a more confident\n * costume. `reviewStatus` is the record's declared review state and derives no evidence\n * tag (`deriveEvidence` never reads it). It is 'proposed' here because no independent\n * physicist review of these records exists; the atlas no longer gates on one (the ruling\n * is recorded in `ACTIVE.md` and `NOTES.md`).\n */\nconst CITATIONS: readonly string[] = [\n  'Olson, Dynamical Analogies - the force-voltage (impedance) analogy between a mass-spring–damper and a series RLC circuit',\n  'Feynman, Lectures on Physics Vol. II, chapter on resonance - the LC oscillator and its mechanical counterpart',\n];\n\n/**\n * Regime coordinates spanning a bridge's premise and conclusion parameters.\n *\n * The two models' parameter names are disjoint in both bridges, so the union\n * is a legal `deriveRegimeGroups` input.\n */\nfunction bridgeRegime(premiseId: string, conclusionId: string): Regime {\n  const parameters = [\n    ...getAtlasModel(premiseId).parameters,\n    ...getAtlasModel(conclusionId).parameters,\n  ];\n  return makeRegime(FAMILY, parameters, [], []);\n}\n\n/**\n * Bridge 1 — the lossless spring ↔ LC dictionary.\n *\n * `m ↔ L`, `k ↔ 1/C`, `x ↔ q`. Under `u = x/x0` (or `q/q0`) and `τ = ω0 t`\n * both systems become `u'' + u = 0`, so the equivalence is exact and\n * invertible.\n *\n * Its counterexample is W2b, which lives on THIS bridge: it shows that adding\n * a resistor to the LC side breaks the equivalence, so bridge 1 is a statement\n * about the LOSSLESS pair only.\n *\n * @internal\n */\nexport const BRIDGE_SPRING_LC: AtlasBridge = {\n  id: 'ab-spring-lc',\n  relation: 'exact-equivalence',\n  premises: ['model-spring'],\n  conclusion: 'model-lc',\n  // Two nondimensionalizations, each with its own ω0; composed, the spring→LC map is\n  // q(t) = (q0/x0)·x(ω_LC t / ω_s). The dimensioned frequencies differ (W1a: ω_s = 2,\n  // ω_LC = 2√2); what is preserved is the frequency in units of ω0.\n  transformation:\n    'force–voltage analogy m ↔ L, k ↔ 1/C, x ↔ q (so ω0² = k/m ↔ 1/(LC)); ' +\n    'u = x/x0 or q/q0, τ = ω0 t; composed, q(t) = (q0/x0)·x(ω_LC t / ω_s)',\n  inverse: 'x = x0 u, t = τ/ω0',\n  preserves: ['the natural frequency in units of ω0 (ω = 1 in τ = ω0 t)', 'energy up to scale', 'phase portrait'],\n  doesNotPreserve: ['physical interpretation', 'units'],\n  sideConditions: ['m, k, L, C > 0', 'lossless', 'unforced', 'x0, q0 nonzero'],\n  regime: bridgeRegime('model-spring', 'model-lc'),\n  counterexamples: [\n    {\n      description:\n        'adding R to bridge 1 breaks it: the circuit L = 2, C = 0.125 of witness W2b with R = 4 has ' +\n        'ζ_RLC = (R/2)√(C/L) = 0.5, ' +\n        'which no lossless spring matches — the trajectories separate by more than 1e-2 at τ = π',\n      witness: 'W2b',\n    },\n  ],\n  // `symbolically-checked` is NOT stored here any more (Phase 4 S4.3): it is\n  // DERIVED from `data/atlas/witness-results.json`, where W1s records the CAS\n  // check of this bridge's dictionary. A stored copy would be a second source\n  // of truth that no run could falsify.\n  witnesses: [\n    {\n      id: 'W1',\n      kind: 'symbolic',\n      test: TEST,\n      tolerance: 'exact — rational exponent arithmetic, no floating point',\n    },\n    {\n      id: 'W1a',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: '|u_spring − u_lc| < 1e-8; each within 1e-8 of cos τ',\n    },\n    {\n      id: 'W1b',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'inverse-mapped trajectory solves m x″ + k x = 0 within 1e-8; a 1% wrong ω0 fails',\n    },\n    { id: 'W2b', kind: 'numeric', test: TEST, tolerance: 'separation > 1e-2 at τ = π' },\n    {\n      id: 'W1s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: k/m under m ↔ L, k ↔ 1/C minus 1/(LC) simplifies to literal 0',\n    },\n  ],\n  citations: CITATIONS,\n  reviewStatus: 'proposed',\n  // spring → LC in relative period error only; see ./norm-transport.ts for what is not declared.\n  normTransports: [SPRING_LC_RELATIVE_PERIOD_TRANSPORT],\n  formalRef: physjsFormalRef('ab-spring-lc'),\n};\n\n/**\n * Bridge 2 — damped spring ↔ RLC, exact ON a side condition.\n *\n * The undamped dictionary of bridge 1 does not extend for free: `b` and `R`\n * are independent parameters, and the nondimensional systems coincide only\n * when the damping ratios match. The side condition is that equality, and W2\n * DERIVES `R` from it rather than assuming a value.\n *\n * @internal\n */\nexport const BRIDGE_DAMPED_RLC: AtlasBridge = {\n  id: 'ab-damped-rlc',\n  relation: 'exact-equivalence',\n  premises: ['model-damped-spring'],\n  conclusion: 'model-rlc',\n  transformation:\n    'force–voltage analogy m ↔ L, k ↔ 1/C, b ↔ R, x ↔ q; ' +\n    'u = x/x0 or q/q0, τ = ω0 t; composed, q(t) = (q0/x0)·x(ω_RLC t / ω_mech)',\n  inverse: 'x = x0 u, t = τ/ω0',\n  preserves: ['the natural frequency in units of ω0 (ω = 1 in τ = ω0 t)', 'damping ratio', 'phase portrait'],\n  doesNotPreserve: ['physical interpretation', 'units'],\n  sideConditions: [\n    'b/√(mk) = R√(C/L), equivalently ζ_mech = b/(2√(mk)) equals ζ_RLC = (R/2)√(C/L)',\n    'm, b, k, L, R, C > 0',\n    'unforced',\n    'x0, q0 nonzero',\n  ],\n  regime: bridgeRegime('model-damped-spring', 'model-rlc'),\n  counterexamples: [],\n  witnesses: [\n    {\n      id: 'W2',\n      kind: 'numeric',\n      test: TEST,\n      tolerance: 'R derived exactly; |u_mech − u_rlc| < 1e-8 at four τ',\n    },\n    {\n      id: 'W2s',\n      kind: 'symbolic',\n      test: WITNESS_RESULTS_TEST,\n      tolerance: 'CAS: b²/(4mk) under m ↔ L, k ↔ 1/C, b ↔ R minus R²C/(4L) simplifies to literal 0',\n    },\n  ],\n  citations: CITATIONS,\n  reviewStatus: 'proposed',\n  formalRef: physjsFormalRef('ab-damped-rlc'),\n};\n\n// ─────────────────────────────────────────────────────────────────────────────\n// Phase 1 re-registration: the Sprint 0 bridges through `RelationContract`\n// ─────────────────────────────────────────────────────────────────────────────\n\n/**\n * Re-register an `AtlasBridge` through the Phase 1 `RelationContract` union.\n *\n * `AtlasBridge` states its two conditional requirements — `exact-equivalence`\n * needs `inverse`, `approximation` needs `bound` — only in doc comments, so a\n * record that breaks them type-checks. `RelationContract` states them in the\n * type system. This function is the one place the Sprint 0 records cross from\n * the weaker shape to the stronger one, and it THROWS on the two cases the\n * doc comments could only describe.\n *\n * It copies rather than restates: every field comes off the bridge, so a\n * contract cannot drift from the record it re-registers.\n *\n * @throws TypeError when the bridge omits the field its relation requires.\n * @internal\n */\nexport function relationContractOf(bridge: AtlasBridge): RelationContract {\n  const { relation, transformation } = bridge;\n  switch (relation) {\n    case 'exact-equivalence': {\n      if (bridge.inverse === undefined) {\n        throw new TypeError(`${bridge.id}: exact-equivalence without an inverse`);\n      }\n      return { type: relation, transformation, inverse: bridge.inverse };\n    }\n    case 'approximation': {\n      if (bridge.bound === undefined) {\n        throw new TypeError(`${bridge.id}: approximation without a bound`);\n      }\n      return { type: relation, transformation, bound: bridge.bound };\n    }\n    case 'coarse-graining':\n      // `bound` is OPTIONAL here: the chain/wave bridge carries none, and the\n      // union allows that rather than forcing an invented one.\n      return bridge.bound === undefined\n        ? { type: relation, transformation }\n        : { type: relation, transformation, bound: bridge.bound };\n    default:\n      return { type: relation, transformation };\n  }\n}\n\n/** Bridge 1 re-registered through the Phase 1 contract. @internal */\nexport const CONTRACT_SPRING_LC: RelationContract = relationContractOf(BRIDGE_SPRING_LC);\n\n/** Bridge 2 re-registered through the Phase 1 contract. @internal */\nexport const CONTRACT_DAMPED_RLC: RelationContract = relationContractOf(BRIDGE_DAMPED_RLC);\n"]}

package/dist/atlas/oscillators/bridges-limits.d.ts

@@ -35,11 +35,11 @@ import type { ApproximationBound, AtlasBridge, RelationContract } from '../types
  * @internal
  */
 export declare function makeApproximation(bound: ApproximationBound): ApproximationBound;
-/** Arithmetic–geometric mean, the standard evaluation route for `K`. @internal */
-export declare function agm(a0: number, b0: number): number;
 /**
  * The EXACT relative period error of the pendulum at amplitude `theta0`:
- * `(2/π) K(sin(θ0/2)) − 1 = 1/AGM(1, cos(θ0/2)) − 1`.
+ * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the
+ * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM
+ * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).
  *
  * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —
  * the bound equals the error it bounds — and it is strictly increasing in
@@ -102,19 +102,6 @@ export declare const AB_PENDULUM_LINEAR: AtlasBridge;
 export declare const AB_DAMPED_MASSLESS: AtlasBridge;
 /** The two approximation bridges of this module, in design-note order. @internal */
 export declare const LIMIT_BRIDGES: readonly AtlasBridge[];
-/**
- * The two approximation bridges re-registered through `RelationContract`.
- *
- * The union REQUIRES `bound` on an `approximation`, and `ApproximationBound`
- * requires its own machine `horizonHolds`, so these two lines are what turn
- * Sprint 0's doc-comment requirement into a compile-and-throw one. Both bounds
- * are the records' own, measured ones — nothing is restated here.
- *
- * @internal
- */
-export declare const CONTRACT_PENDULUM_LINEAR: RelationContract;
-/** @internal */
-export declare const CONTRACT_DAMPED_MASSLESS: RelationContract;
 /** The two limit contracts, in design-note order. @internal */
 export declare const LIMIT_CONTRACTS: readonly RelationContract[];
 //# sourceMappingURL=bridges-limits.d.ts.map

package/dist/atlas/oscillators/bridges-limits.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-limits.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAKH,OAAO,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,gBAAgB,EAAU,MAAM,aAAa,CAAC;AAS7F;;;;;;;;;;;;GAYG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,kBAAkB,GAAG,kBAAkB,CAO/E;AAED,kFAAkF;AAClF,wBAAgB,GAAG,CAAC,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,GAAG,MAAM,CAOlD;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,qBAAqB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAItF;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAKpF;AAuCD;;;;;;;;;;;;;;;;GAgBG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAgFhC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAmFhC,CAAC;AAEF,oFAAoF;AACpF,eAAO,MAAM,aAAa,EAAE,SAAS,WAAW,EAA6C,CAAC;AAE9F;;;;;;;;;GASG;AACH,eAAO,MAAM,wBAAwB,EAAE,gBACC,CAAC;AAEzC,gBAAgB;AAChB,eAAO,MAAM,wBAAwB,EAAE,gBACC,CAAC;AAEzC,+DAA+D;AAC/D,eAAO,MAAM,eAAe,EAAE,SAAS,gBAAgB,EAGtD,CAAC"}
+{"version":3,"file":"bridges-limits.d.ts","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAOH,OAAO,KAAK,EAAE,kBAAkB,EAAE,WAAW,EAAE,gBAAgB,EAAU,MAAM,aAAa,CAAC;AAS7F;;;;;;;;;;;;GAYG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,kBAAkB,GAAG,kBAAkB,CAO/E;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,qBAAqB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAMtF;AAED;;;;;;;;;;;;GAYG;AACH,wBAAgB,mBAAmB,CAAC,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,MAAM,CAKpF;AAiCD;;;;;;;;;;;;;;;;GAgBG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAgFhC,CAAC;AAEF;;;;;;;;;;GAUG;AACH,eAAO,MAAM,kBAAkB,EAAE,WAmFhC,CAAC;AAEF,oFAAoF;AACpF,eAAO,MAAM,aAAa,EAAE,SAAS,WAAW,EAA6C,CAAC;AAmB9F,+DAA+D;AAC/D,eAAO,MAAM,eAAe,EAAE,SAAS,gBAAgB,EAGtD,CAAC"}

package/dist/atlas/oscillators/bridges-limits.js

@@ -20,8 +20,10 @@
  *
  * @module atlas/oscillators/bridges-limits
  */
+import { ellipticK } from '@danielsimonjr/mathts-functions';
+import { numberFrom } from '../../core/mathts-read.js';
 import { physjsFormalRef } from '../physjs-ref.js';
-import { deriveRegimeGroups } from '../regime.js';
+import { makeRegime } from '../regime.js';
 import { MissingHorizonError } from '../types.js';
 import { relationContractOf } from './bridges-exact.js';
 import { DAMPING, SPRING_CONSTANT } from './dimensions.js';
@@ -48,18 +50,11 @@ export function makeApproximation(bound) {
     }
     return bound;
 }
-/** Arithmetic–geometric mean, the standard evaluation route for `K`. @internal */
-export function agm(a0, b0) {
-    let a = a0;
-    let b = b0;
-    for (let i = 0; i < 60 && a !== b; i++) {
-        [a, b] = [(a + b) / 2, Math.sqrt(a * b)];
-    }
-    return a;
-}
 /**
  * The EXACT relative period error of the pendulum at amplitude `theta0`:
- * `(2/π) K(sin(θ0/2)) − 1 = 1/AGM(1, cos(θ0/2)) − 1`.
+ * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the
+ * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM
+ * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).
  *
  * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —
  * the bound equals the error it bounds — and it is strictly increasing in
@@ -79,7 +74,9 @@ export function pendulumPeriodErrorAt(params) {
     const { theta0 } = params;
     if (!Number.isFinite(theta0) || Math.abs(theta0) >= Math.PI)
         return Infinity;
-    return 1 / agm(1, Math.cos(theta0 / 2)) - 1;
+    const m = Math.sin(theta0 / 2) ** 2;
+    const k = numberFrom('ellipticK', ellipticK(m));
+    return (2 / Math.PI) * k - 1;
 }
 /**
  * The position-offset bound of `ab-damped-massless` at a point:
@@ -103,14 +100,10 @@ export function dampedOffsetBoundAt(params) {
     return (2 * (1 + Math.abs(v0)) * m) / b;
 }
 /** Pendulum regime: `θ0 ≤ 0.5 rad`, stated on the dimensionless input. @internal */
-const PENDULUM_REGIME = {
-    family: FAMILY,
-    inequalities: [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }],
-    groupDefinitions: deriveRegimeGroups(FAMILY, [
-        { name: 'g', dim: ACCELERATION },
-        { name: 'ell', dim: LENGTH },
-    ], ['theta0']),
-};
+const PENDULUM_REGIME = makeRegime(FAMILY, [
+    { name: 'g', dim: ACCELERATION },
+    { name: 'ell', dim: LENGTH },
+], ['theta0'], [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }]);
 /**
  * Damped-spring regime: overdamped, `m k / b² < 1/4`.
  *
@@ -120,15 +113,11 @@ const PENDULUM_REGIME = {
  *
  * @internal
  */
-const DAMPED_REGIME = {
-    family: FAMILY,
-    inequalities: [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }],
-    groupDefinitions: deriveRegimeGroups(FAMILY, [
-        { name: 'm', dim: MASS },
-        { name: 'b', dim: DAMPING },
-        { name: 'k', dim: SPRING_CONSTANT },
-    ], []),
-};
+const DAMPED_REGIME = makeRegime(FAMILY, [
+    { name: 'm', dim: MASS },
+    { name: 'b', dim: DAMPING },
+    { name: 'k', dim: SPRING_CONSTANT },
+], [], [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }]);
 /**
  * `ab-pendulum-linear` — small-angle pendulum → harmonic spring.
  *
@@ -158,8 +147,8 @@ export const AB_PENDULUM_LINEAR = {
     bound: makeApproximation({
         K: 1,
         // Sup over `θ0 ≤ 0.5`: the error is strictly increasing in |θ0|, so the
-        // edge value IS the supremum. 0.0158525311014... (AGM), independently
-        // reproduced by quadrature of the complete elliptic integral.
+        // edge value IS the supremum. 0.0158525311014... (MathTS ellipticK, by the
+        // AGM), independently reproduced by quadrature of the complete elliptic integral.
         delta: pendulumPeriodErrorAt({ theta0: 0.5 }),
         deltaAt: pendulumPeriodErrorAt,
         deltaAtBasis: 'closed-form',
@@ -224,7 +213,7 @@ export const AB_PENDULUM_LINEAR = {
     // certifies the TRANSFORMATION (mass mℓ², spring constant mgℓ, ω0² = g/ℓ),
     // not `bound.delta`. Physlib's period results concern `periodFormula`, which
     // its own TODO has not yet tied to the motion. The sanity lemmas instantiate
-    // that imported statement. Keyed by `formal/physjs/manifest.json`.
+    // that imported statement. Keyed by `formal/lean/manifest.json`.
     formalRef: physjsFormalRef('ab-pendulum-linear'),
 };
 /**
@@ -333,9 +322,9 @@ export const LIMIT_BRIDGES = [AB_PENDULUM_LINEAR, AB_DAMPED_MASSLESS];
  *
  * @internal
  */
-export const CONTRACT_PENDULUM_LINEAR = relationContractOf(AB_PENDULUM_LINEAR);
+const CONTRACT_PENDULUM_LINEAR = relationContractOf(AB_PENDULUM_LINEAR);
 /** @internal */
-export const CONTRACT_DAMPED_MASSLESS = relationContractOf(AB_DAMPED_MASSLESS);
+const CONTRACT_DAMPED_MASSLESS = relationContractOf(AB_DAMPED_MASSLESS);
 /** The two limit contracts, in design-note order. @internal */
 export const LIMIT_CONTRACTS = [
     CONTRACT_PENDULUM_LINEAR,

package/dist/atlas/oscillators/bridges-limits.js.map

@@ -1 +1 @@
-{"version":3,"file":"bridges-limits.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,kBAAkB,EAAE,MAAM,cAAc,CAAC;AAClD,OAAO,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAElD,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAC3D,OAAO,EAAE,oBAAoB,EAAE,MAAM,YAAY,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAExE,MAAM,MAAM,GAAG,aAAa,CAAC;AAC7B,MAAM,SAAS,GAAG,wCAAwC,CAAC;AAE3D;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,iBAAiB,CAAC,KAAyB;IACzD,IAAI,KAAK,CAAC,OAAO,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QAChC,MAAM,IAAI,mBAAmB,CAC3B,+EAA+E,CAChF,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED,kFAAkF;AAClF,MAAM,UAAU,GAAG,CAAC,EAAU,EAAE,EAAU;IACxC,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACvC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3C,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,MAAM,UAAU,qBAAqB,CAAC,MAAwC;IAC5E,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;IAC1B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,EAAE;QAAE,OAAO,QAAQ,CAAC;IAC7E,OAAO,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AAC9C,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAwC;IAC1E,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,MAAM,CAAC;IACnC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/D,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,oFAAoF;AACpF,MAAM,eAAe,GAAW;IAC9B,MAAM,EAAE,MAAM;IACd,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,cAAc,EAAE,CAAC;IAChF,gBAAgB,EAAE,kBAAkB,CAClC,MAAM,EACN;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;QAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;KAC7B,EACD,CAAC,QAAQ,CAAC,CACX;CACF,CAAC;AAEF;;;;;;;;GAQG;AACH,MAAM,aAAa,GAAW;IAC5B,MAAM,EAAE,MAAM;IACd,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,cAAc,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC;IAC/E,gBAAgB,EAAE,kBAAkB,CAClC,MAAM,EACN;QACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;QACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;QAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;KACpC,EACD,EAAE,CACH;CACF,CAAC;AAEF;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,cAAc;IAC1B,cAAc,EAAE,uCAAuC;IACvD,SAAS,EAAE,CAAC,iCAAiC,EAAE,qBAAqB,EAAE,wBAAwB,CAAC;IAC/F,eAAe,EAAE,CAAC,wCAAwC,EAAE,gCAAgC,CAAC;IAC7F,cAAc,EAAE,CAAC,cAAc,EAAE,wDAAwD,CAAC;IAC1F,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,wEAAwE;QACxE,sEAAsE;QACtE,8DAA8D;QAC9D,KAAK,EAAE,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;QAC7C,OAAO,EAAE,qBAAqB;QAC9B,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,oBAAoB;QAC1B,MAAM,EAAE,cAAc;QACtB,OAAO,EAAE,8DAA8D;QACvE,sEAAsE;QACtE,yEAAyE;QACzE,yEAAyE;QACzE,wEAAwE;QACxE,uEAAuE;QACvE,uEAAuE;QACvE,kEAAkE;QAClE,+DAA+D;QAC/D,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;YAC9B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACnF,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,CAAC;QACpC,CAAC;QACD,cAAc,EAAE,cAAc;QAC9B,cAAc,EAAE,SAAS;QACzB,wEAAwE;QACxE,cAAc;QACd,UAAU,EAAE,CAAC,yCAAyC,CAAC;KACxD,CAAC;IACF,MAAM,EAAE,eAAe;IACvB,eAAe,EAAE;QACf;YACE,WAAW,EACT,oFAAoF;gBACpF,yDAAyD;YAC3D,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,4EAA4E;SACxF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,6CAA6C;SACzD;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,iDAAiD;KAClD;IACD,YAAY,EAAE,UAAU;IACxB,4EAA4E;IAC5E,yEAAyE;IACzE,2EAA2E;IAC3E,6EAA6E;IAC7E,6EAA6E;IAC7E,mEAAmE;IACnE,SAAS,EAAE,eAAe,CAAC,oBAAoB,CAAC;CACjD,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,mBAAmB;IAC/B,cAAc,EAAE,+BAA+B;IAC/C,SAAS,EAAE,CAAC,sCAAsC,EAAE,wCAAwC,CAAC;IAC7F,eAAe,EAAE;QACf,qCAAqC;QACrC,qEAAqE;QACrE,wBAAwB;KACzB;IACD,cAAc,EAAE;QACd,4BAA4B;QAC5B,4DAA4D;KAC7D;IACD,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,yEAAyE;QACzE,oEAAoE;QACpE,0EAA0E;QAC1E,0EAA0E;QAC1E,uEAAuE;QACvE,0EAA0E;QAC1E,mEAAmE;QACnE,kBAAkB;QAClB,wEAAwE;QACxE,4EAA4E;QAC5E,yEAAyE;QACzE,kEAAkE;QAClE,2EAA2E;QAC3E,6EAA6E;QAC7E,2DAA2D;QAC3D,yDAAyD;QACzD,yDAAyD;QACzD,KAAK,EAAE,mBAAmB,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;QACjE,OAAO,EAAE,mBAAmB;QAC5B,sFAAsF;QACtF,YAAY,EAAE,uBAAuB;QACrC,IAAI,EAAE,mCAAmC;QACzC,MAAM,EAAE,+DAA+D;QACvE,OAAO,EAAE,wCAAwC;QACjD,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YACxB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACxE,OAAO,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,CAAC;QACD,cAAc,EAAE,0DAA0D;QAC1E,cAAc,EAAE,UAAU;QAC1B,UAAU,EAAE;YACV,+DAA+D;YAC/D,0BAA0B;SAC3B;KACF,CAAC;IACF,MAAM,EAAE,aAAa;IACrB,eAAe,EAAE;QACf;YACE,WAAW,EACT,8EAA8E;gBAC9E,wEAAwE;YAC1E,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,wEAAwE;SACpF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EACP,2FAA2F;SAC9F;KACF;IACD,SAAS,EAAE;QACT,4FAA4F;QAC5F,mFAAmF;KACpF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,oFAAoF;AACpF,MAAM,CAAC,MAAM,aAAa,GAA2B,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAE9F;;;;;;;;;GASG;AACH,MAAM,CAAC,MAAM,wBAAwB,GACnC,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,gBAAgB;AAChB,MAAM,CAAC,MAAM,wBAAwB,GACnC,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,+DAA+D;AAC/D,MAAM,CAAC,MAAM,eAAe,GAAgC;IAC1D,wBAAwB;IACxB,wBAAwB;CACzB,CAAC","sourcesContent":["/**\n * The two Phase 0 APPROXIMATION bridges — design note §3, bridges 3 and 4.\n *\n * - `ab-pendulum-linear` — a REGULAR limit. `θ'' + ω0² sin θ = 0` reduces to\n *   `θ'' + ω0² θ = 0` as `θ0 → 0`. The order of the system does not drop, the\n *   error is `O(θ0²)` in the relative period, and the approximation is\n *   NON-UNIFORM in time: a fixed relative period error accumulates phase\n *   without bound, so the bound carries the horizon `t ≪ 16 T0/θ0²`.\n * - `ab-damped-massless` — a SINGULAR limit. `m x'' + b x' + k x = 0` reduces\n *   to `b x' + k x = 0` as `m → 0`. The order DROPS from two to one, one\n *   initial condition is lost, and the reduced model cannot satisfy `x'(0)`.\n *   The bound therefore holds only OUTSIDE the boundary layer, `t ≥ 5 m/b`.\n *\n * On the position offset of the singular bridge: `x_full − x_red` is an\n * `O((1 + |v0|) m)` deviation that PERSISTS THROUGHOUT THE OUTER REGION and\n * decays only on the slow time scale `b/k`. It is the matched-asymptotics\n * correction, not a boundary-layer transient, and waiting out the fast mode\n * does not remove it. (Design note §9 YELLOW (5) — \"non-vanishing\" was the\n * imprecise wording this replaces: the offset does go to zero as `t → ∞`.)\n *\n * @module atlas/oscillators/bridges-limits\n */\n\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport { deriveRegimeGroups } from '../regime.js';\nimport { MissingHorizonError } from '../types.js';\nimport type { ApproximationBound, AtlasBridge, RelationContract, Regime } from '../types.js';\nimport { relationContractOf } from './bridges-exact.js';\nimport { DAMPING, SPRING_CONSTANT } from './dimensions.js';\nimport { RELATIVE_PERIOD_NORM } from './norms.js';\nimport { ACCELERATION, LENGTH, MASS } from '../../dimensional/types.js';\n\nconst FAMILY = 'oscillators';\nconst TEST_FILE = 'tests/atlas/oscillators-limits.test.ts';\n\n/**\n * Build an `ApproximationBound`, rejecting one with no prose horizon.\n *\n * Every `approximation` bridge in this module is constructed through here, so\n * the mandatory-horizon rule of `ApproximationBound` has something that FAILS\n * rather than a comment asking for compliance. `horizonHolds` is mandatory by\n * the type; the prose `horizon` is the field a caller can satisfy with `''`,\n * so it is the one checked at run time. `uniformity: null` is accepted: not\n * yet analysed is a real state, and the refusal lives in `boundPath`.\n *\n * @throws MissingHorizonError if `horizon` is empty or whitespace.\n * @internal\n */\nexport function makeApproximation(bound: ApproximationBound): ApproximationBound {\n  if (bound.horizon.trim() === '') {\n    throw new MissingHorizonError(\n      'an approximation bound requires a non-empty prose horizon beside horizonHolds',\n    );\n  }\n  return bound;\n}\n\n/** Arithmetic–geometric mean, the standard evaluation route for `K`. @internal */\nexport function agm(a0: number, b0: number): number {\n  let a = a0;\n  let b = b0;\n  for (let i = 0; i < 60 && a !== b; i++) {\n    [a, b] = [(a + b) / 2, Math.sqrt(a * b)];\n  }\n  return a;\n}\n\n/**\n * The EXACT relative period error of the pendulum at amplitude `theta0`:\n * `(2/π) K(sin(θ0/2)) − 1 = 1/AGM(1, cos(θ0/2)) − 1`.\n *\n * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —\n * the bound equals the error it bounds — and it is strictly increasing in\n * `|θ0|`, so its value at the edge of a declared amplitude range is the\n * supremum over that range. That last property is what lets the record state\n * one scalar for the whole domain honestly.\n *\n * The series `θ0²/16` is the LEADING TERM of this quantity, not an upper bound\n * on it: at the domain edge `θ0 = 0.5` the series gives `0.0156250` while the\n * exact error is `0.0158525311`, so the series UNDERSTATES by 1.456%. A record\n * that declares the series value is violated at its own boundary.\n *\n * @returns `Infinity` outside `|θ0| < π`, where no finite period error exists.\n * @internal\n */\nexport function pendulumPeriodErrorAt(params: Readonly<Record<string, number>>): number {\n  const { theta0 } = params;\n  if (!Number.isFinite(theta0) || Math.abs(theta0) >= Math.PI) return Infinity;\n  return 1 / agm(1, Math.cos(theta0 / 2)) - 1;\n}\n\n/**\n * The position-offset bound of `ab-damped-massless` at a point:\n * `2 (1 + |v0|) m / b`, the formula this record's docstring has always\n * carried. It is a bound ONLY at the witness normalisation `b = k = x0 = 1`:\n * the leading offset is `(m/b)(v0 + k x0/b)`, and the formula has no `k` or\n * `x0` in it. At `m = 1e-3, b = 1, k = 100` (inside `m k/b² < 1/4`) the true\n * error is more than 20 times the formula.\n *\n * @returns `Infinity` unless `m`, `b`, `k`, `x0` and `v0` are all finite and\n *   `b = k = x0 = 1` — a missing parameter would otherwise silently return a\n *   bound for a normalisation the caller did not ask about.\n * @internal\n */\nexport function dampedOffsetBoundAt(params: Readonly<Record<string, number>>): number {\n  const { m, b, k, x0, v0 } = params;\n  if (![m, b, k, x0, v0].every(Number.isFinite)) return Infinity;\n  if (b !== 1 || k !== 1 || x0 !== 1) return Infinity;\n  return (2 * (1 + Math.abs(v0)) * m) / b;\n}\n\n/** Pendulum regime: `θ0 ≤ 0.5 rad`, stated on the dimensionless input. @internal */\nconst PENDULUM_REGIME: Regime = {\n  family: FAMILY,\n  inequalities: [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }],\n  groupDefinitions: deriveRegimeGroups(\n    FAMILY,\n    [\n      { name: 'g', dim: ACCELERATION },\n      { name: 'ell', dim: LENGTH },\n    ],\n    ['theta0'],\n  ),\n};\n\n/**\n * Damped-spring regime: overdamped, `m k / b² < 1/4`.\n *\n * The inequality is written on the π-group, not on `ζ = ½ (m k / b²)^(−½)`,\n * which is a derived quantity and not a group (`RegimeInequality` doc).\n * `m k / b² < 1/4` is exactly `ζ > 1`, recorded as the display alias.\n *\n * @internal\n */\nconst DAMPED_REGIME: Regime = {\n  family: FAMILY,\n  inequalities: [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }],\n  groupDefinitions: deriveRegimeGroups(\n    FAMILY,\n    [\n      { name: 'm', dim: MASS },\n      { name: 'b', dim: DAMPING },\n      { name: 'k', dim: SPRING_CONSTANT },\n    ],\n    [],\n  ),\n};\n\n/**\n * `ab-pendulum-linear` — small-angle pendulum → harmonic spring.\n *\n * `θ0²/16` is the leading term of the exact relative period error\n * `(2/π) K(sin(θ0/2)) − 1`. At `θ0 = 0.2` the residual is `5.744e-6`, which\n * the next series term `11 θ0⁴/3072 = 5.729e-6` accounts for to `1.5e-8`.\n *\n * `delta` is the EXACT error at the edge of the declared range, not that\n * series term. The record previously declared `0.5²/16 = 0.0156250`, which is\n * the series evaluated at `θ0 = 0.5` — and the exact error there is\n * `0.0158525311`, so the declared bound was VIOLATED AT ITS OWN BOUNDARY by\n * 1.456%. A truncated series is an approximation of the error, not a bound on\n * it, and the sign of its omitted tail decides which. Here the tail is\n * positive, so the series sits below what it was asked to cover.\n *\n * @internal\n */\nexport const AB_PENDULUM_LINEAR: AtlasBridge = {\n  id: 'ab-pendulum-linear',\n  relation: 'approximation',\n  premises: ['model-pendulum'],\n  conclusion: 'model-spring',\n  transformation: 'sin θ → θ, with x ↔ ℓ θ and ω0² = g/ℓ',\n  preserves: ['harmonic frequency ω0 to O(θ0²)', 'energy conservation', 'time-reversal symmetry'],\n  doesNotPreserve: ['the amplitude dependence of the period', 'phase over times t ≳ 16 T0/θ0²'],\n  sideConditions: ['θ0 ≤ 0.5 rad', 'the bound is a PERIOD error and is not uniform in time'],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over `θ0 ≤ 0.5`: the error is strictly increasing in |θ0|, so the\n    // edge value IS the supremum. 0.0158525311014... (AGM), independently\n    // reproduced by quadrature of the complete elliptic integral.\n    delta: pendulumPeriodErrorAt({ theta0: 0.5 }),\n    deltaAt: pendulumPeriodErrorAt,\n    deltaAtBasis: 'closed-form',\n    norm: RELATIVE_PERIOD_NORM,\n    domain: 'θ0 ≤ 0.5 rad',\n    horizon: 't ≪ 16 T0/θ0²; machine form t < 4 T0/θ0², the π/2-drift time',\n    // ⚠ The machine form is a QUARTER of the prose scale, and it must be.\n    // `16 T0/θ0²` is the time at which the accumulated phase drift reaches a\n    // FULL 2π (measured: 400 T0 at θ0 = 0.2, drift 6.2976 rad) — by then the\n    // approximation has lapped, so `t < 16 T0/θ0²` is not a horizon at all.\n    // The `≪` is what carries the whole claim, and `<` does not render it.\n    // `4 T0/θ0²` is the π/2-drift time, which is exactly what W7b measures\n    // independently (99.77 cycles at θ0 = 0.2). It is also what makes\n    // `horizonHolds(200 T0)` false and `horizonHolds(10 T0)` true.\n    horizonHolds: (t, params) => {\n      const { T0, theta0 } = params;\n      if (!Number.isFinite(T0) || !Number.isFinite(theta0) || theta0 === 0) return false;\n      return t < (4 * T0) / theta0 ** 2;\n    },\n    parameterRange: 'θ0 ≤ 0.5 rad',\n    limitCharacter: 'regular',\n    // A period error. The record says it is not uniform in time, so time is\n    // not listed.\n    uniformity: ['one period, for θ0 in the stated domain'],\n  }),\n  regime: PENDULUM_REGIME,\n  counterexamples: [\n    {\n      description:\n        'the approximation is non-uniform in time: at θ0 = 0.2 the phase drift reaches π/2 ' +\n        'after ~100 cycles, however small the per-cycle error is',\n      witness: 'W7b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W7',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'T/T0 − 1 ∈ [0.002505, 0.002507]; residual ∈ [5.70e-6, 5.76e-6] at θ0 = 0.2',\n    },\n    {\n      id: 'W7b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'cycles to π/2 drift ∈ [99, 101] at θ0 = 0.2',\n    },\n    {\n      id: 'W7c',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'RK4 zero-crossing lag within 0.05° of the elliptic prediction 89.98°',\n    },\n  ],\n  citations: [\n    'Landau & Lifshitz, Mechanics §11 (pendulum period as a complete elliptic integral)',\n    'Abramowitz & Stegun §17.6 (AGM evaluation of K)',\n  ],\n  reviewStatus: 'proposed',\n  // Milestone 1. The reference is the PhysJS theorem, which imports Physlib's\n  // `ClassicalMechanics.SimplePendulum.linearizedEquationOfMotion_iff`. It\n  // certifies the TRANSFORMATION (mass mℓ², spring constant mgℓ, ω0² = g/ℓ),\n  // not `bound.delta`. Physlib's period results concern `periodFormula`, which\n  // its own TODO has not yet tied to the motion. The sanity lemmas instantiate\n  // that imported statement. Keyed by `formal/physjs/manifest.json`.\n  formalRef: physjsFormalRef('ab-pendulum-linear'),\n};\n\n/**\n * `ab-damped-massless` — the SINGULAR `m → 0` limit of the damped spring.\n *\n * `deltaAt = 2(1 + |v0|) m/b` bounds the position offset outside the boundary\n * layer, and `delta` is its supremum over the declared range rather than its\n * value at one fixture. The lost initial condition is visible in the VELOCITY: inside the\n * layer the reduced model's `x'` is wrong by O(1/m·…), and it is only after\n * `t ≈ 5 m/b` that the fast mode has decayed.\n *\n * @internal\n */\nexport const AB_DAMPED_MASSLESS: AtlasBridge = {\n  id: 'ab-damped-massless',\n  relation: 'approximation',\n  premises: ['model-damped-spring'],\n  conclusion: 'model-first-order',\n  transformation: 'm → 0, dropping the m x″ term',\n  preserves: ['the slow relaxation rate k/b to O(m)', 'the sign and monotonicity of the decay'],\n  doesNotPreserve: [\n    'the order of the system (two → one)',\n    \"the initial condition x'(0), which the reduced model cannot satisfy\",\n    'the fast mode r ≈ −b/m',\n  ],\n  sideConditions: [\n    'overdamped: m k / b² < 1/4',\n    'the bound holds only outside the boundary layer, t ≥ 5 m/b',\n  ],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over the DECLARED range, not the fixture. `2(1+|v0|)m/b` increases\n    // in both m and |v0|, and the declared range is `m k/b² < 1/4` with\n    // `|v0| ≤ 5`; at the witness normalisation `b = k = 1` that is `m < 1/4`,\n    // so the supremum is `2·(1+5)·(1/4) = 3`, approached at the open edge and\n    // not attained. The record previously declared `2·(1+5)·1e-3 = 0.012`,\n    // the same formula frozen at the ONE fixture mass `m = 1e-3` — true error\n    // there is 5.96e-3, but at m = 0.24 it is 5.19e-1, forty times the\n    // declared bound.\n    // KNOWN LOOSENESS (0.47.0 persona re-test, L9 residual; documented, not\n    // tightened): for x0 = 1 the measured sup of |x − x_reduced| over the whole\n    // declared range (m < 1/4, |v0| ≤ 5, t ≥ 5m/b) is 0.524, at m ≈ 0.21 and\n    // v0 = +5, and it does not grow toward critical damping (0.514 at\n    // m = 0.249999). So δ = 3 holds but is 5.7× loose, and it exceeds the unit\n    // signal it bounds. A tighter δ needs a bound with its own basis, not a grid\n    // maximum; until then `upt path` withholds the point bound\n    // (deltaAtBasis below). Independent check: the persona's\n    // checks/recheck_047.py and a 115 × 2 grid over (m, ±5).\n    delta: dampedOffsetBoundAt({ m: 0.25, b: 1, k: 1, x0: 1, v0: 5 }),\n    deltaAt: dampedOffsetBoundAt,\n    // Witness W8b supports it numerically for m ∈ {1e-1, 1e-2, 1e-3}; no proof covers it.\n    deltaAtBasis: 'numerically-supported',\n    norm: 'sup |x − x_reduced| for t ≥ 5 m/b',\n    domain: 't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n    horizon: 't ≥ 5 m/b (outside the boundary layer)',\n    horizonHolds: (t, params) => {\n      const { m, b } = params;\n      if (!Number.isFinite(m) || !Number.isFinite(b) || b === 0) return false;\n      return t >= (5 * m) / b;\n    },\n    parameterRange: 'm k / b² < 1/4, |v0| ≤ 5; with b = k = 1 this is m < 1/4',\n    limitCharacter: 'singular',\n    uniformity: [\n      't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n      'm k / b² < 1/4, |v0| ≤ 5',\n    ],\n  }),\n  regime: DAMPED_REGIME,\n  counterexamples: [\n    {\n      description:\n        \"inside the boundary layer the reduced model's velocity is wrong by O(1): at \" +\n        \"t = 0.5 m/b with v0 = 5 the velocity error is 3.6, against x'_red = −1\",\n      witness: 'W8b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W8',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: '|r_slow + k/b| < 2m and |r_fast·m + b| < 2m for m ∈ {1e-1, 1e-2, 1e-3}',\n    },\n    {\n      id: 'W8b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance:\n        \"|x_full − x_red| < 2(1+v0)·m for t ≥ 5 m/b; |x'| error > 1 at 0.5 m/b and < 0.05 at 5 m/b\",\n    },\n  ],\n  citations: [\n    \"Kevorkian & Cole, Multiple Scale and Singular Perturbation Methods §2 (Tikhonov's theorem)\",\n    'Verhulst, Methods and Applications of Singular Perturbations §8 (boundary layers)',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/** The two approximation bridges of this module, in design-note order. @internal */\nexport const LIMIT_BRIDGES: readonly AtlasBridge[] = [AB_PENDULUM_LINEAR, AB_DAMPED_MASSLESS];\n\n/**\n * The two approximation bridges re-registered through `RelationContract`.\n *\n * The union REQUIRES `bound` on an `approximation`, and `ApproximationBound`\n * requires its own machine `horizonHolds`, so these two lines are what turn\n * Sprint 0's doc-comment requirement into a compile-and-throw one. Both bounds\n * are the records' own, measured ones — nothing is restated here.\n *\n * @internal\n */\nexport const CONTRACT_PENDULUM_LINEAR: RelationContract =\n  relationContractOf(AB_PENDULUM_LINEAR);\n\n/** @internal */\nexport const CONTRACT_DAMPED_MASSLESS: RelationContract =\n  relationContractOf(AB_DAMPED_MASSLESS);\n\n/** The two limit contracts, in design-note order. @internal */\nexport const LIMIT_CONTRACTS: readonly RelationContract[] = [\n  CONTRACT_PENDULUM_LINEAR,\n  CONTRACT_DAMPED_MASSLESS,\n];\n"]}
+{"version":3,"file":"bridges-limits.js","sourceRoot":"","sources":["../../../src/atlas/oscillators/bridges-limits.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;GAqBG;AAEH,OAAO,EAAE,SAAS,EAAE,MAAM,iCAAiC,CAAC;AAC5D,OAAO,EAAE,UAAU,EAAE,MAAM,2BAA2B,CAAC;AACvD,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,mBAAmB,EAAE,MAAM,aAAa,CAAC;AAElD,OAAO,EAAE,kBAAkB,EAAE,MAAM,oBAAoB,CAAC;AACxD,OAAO,EAAE,OAAO,EAAE,eAAe,EAAE,MAAM,iBAAiB,CAAC;AAC3D,OAAO,EAAE,oBAAoB,EAAE,MAAM,YAAY,CAAC;AAClD,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,4BAA4B,CAAC;AAExE,MAAM,MAAM,GAAG,aAAa,CAAC;AAC7B,MAAM,SAAS,GAAG,wCAAwC,CAAC;AAE3D;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,iBAAiB,CAAC,KAAyB;IACzD,IAAI,KAAK,CAAC,OAAO,CAAC,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;QAChC,MAAM,IAAI,mBAAmB,CAC3B,+EAA+E,CAChF,CAAC;IACJ,CAAC;IACD,OAAO,KAAK,CAAC;AACf,CAAC;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,MAAM,UAAU,qBAAqB,CAAC,MAAwC;IAC5E,MAAM,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;IAC1B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,IAAI,CAAC,EAAE;QAAE,OAAO,QAAQ,CAAC;IAC7E,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,UAAU,CAAC,WAAW,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAC/B,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,mBAAmB,CAAC,MAAwC;IAC1E,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,MAAM,CAAC;IACnC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,QAAQ,CAAC;QAAE,OAAO,QAAQ,CAAC;IAC/D,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,EAAE,KAAK,CAAC;QAAE,OAAO,QAAQ,CAAC;IACpD,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;AAC1C,CAAC;AAED,oFAAoF;AACpF,MAAM,eAAe,GAAW,UAAU,CACxC,MAAM,EACN;IACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,YAAY,EAAE;IAChC,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,MAAM,EAAE;CAC7B,EACD,CAAC,QAAQ,CAAC,EACV,CAAC,EAAE,KAAK,EAAE,QAAQ,EAAE,EAAE,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,KAAK,EAAE,cAAc,EAAE,CAAC,CACnE,CAAC;AAEF;;;;;;;;GAQG;AACH,MAAM,aAAa,GAAW,UAAU,CACtC,MAAM,EACN;IACE,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE;IACxB,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,OAAO,EAAE;IAC3B,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,eAAe,EAAE;CACpC,EACD,EAAE,EACF,CAAC,EAAE,KAAK,EAAE,cAAc,EAAE,EAAE,EAAE,GAAG,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,OAAO,EAAE,CAAC,CAClE,CAAC;AAEF;;;;;;;;;;;;;;;;GAgBG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,gBAAgB,CAAC;IAC5B,UAAU,EAAE,cAAc;IAC1B,cAAc,EAAE,uCAAuC;IACvD,SAAS,EAAE,CAAC,iCAAiC,EAAE,qBAAqB,EAAE,wBAAwB,CAAC;IAC/F,eAAe,EAAE,CAAC,wCAAwC,EAAE,gCAAgC,CAAC;IAC7F,cAAc,EAAE,CAAC,cAAc,EAAE,wDAAwD,CAAC;IAC1F,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,wEAAwE;QACxE,2EAA2E;QAC3E,kFAAkF;QAClF,KAAK,EAAE,qBAAqB,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,CAAC;QAC7C,OAAO,EAAE,qBAAqB;QAC9B,YAAY,EAAE,aAAa;QAC3B,IAAI,EAAE,oBAAoB;QAC1B,MAAM,EAAE,cAAc;QACtB,OAAO,EAAE,8DAA8D;QACvE,sEAAsE;QACtE,yEAAyE;QACzE,yEAAyE;QACzE,wEAAwE;QACxE,uEAAuE;QACvE,uEAAuE;QACvE,kEAAkE;QAClE,+DAA+D;QAC/D,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,MAAM,CAAC;YAC9B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACnF,OAAO,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,MAAM,IAAI,CAAC,CAAC;QACpC,CAAC;QACD,cAAc,EAAE,cAAc;QAC9B,cAAc,EAAE,SAAS;QACzB,wEAAwE;QACxE,cAAc;QACd,UAAU,EAAE,CAAC,yCAAyC,CAAC;KACxD,CAAC;IACF,MAAM,EAAE,eAAe;IACvB,eAAe,EAAE;QACf;YACE,WAAW,EACT,oFAAoF;gBACpF,yDAAyD;YAC3D,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,4EAA4E;SACxF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,6CAA6C;SACzD;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,sEAAsE;SAClF;KACF;IACD,SAAS,EAAE;QACT,oFAAoF;QACpF,iDAAiD;KAClD;IACD,YAAY,EAAE,UAAU;IACxB,4EAA4E;IAC5E,yEAAyE;IACzE,2EAA2E;IAC3E,6EAA6E;IAC7E,6EAA6E;IAC7E,iEAAiE;IACjE,SAAS,EAAE,eAAe,CAAC,oBAAoB,CAAC;CACjD,CAAC;AAEF;;;;;;;;;;GAUG;AACH,MAAM,CAAC,MAAM,kBAAkB,GAAgB;IAC7C,EAAE,EAAE,oBAAoB;IACxB,QAAQ,EAAE,eAAe;IACzB,QAAQ,EAAE,CAAC,qBAAqB,CAAC;IACjC,UAAU,EAAE,mBAAmB;IAC/B,cAAc,EAAE,+BAA+B;IAC/C,SAAS,EAAE,CAAC,sCAAsC,EAAE,wCAAwC,CAAC;IAC7F,eAAe,EAAE;QACf,qCAAqC;QACrC,qEAAqE;QACrE,wBAAwB;KACzB;IACD,cAAc,EAAE;QACd,4BAA4B;QAC5B,4DAA4D;KAC7D;IACD,KAAK,EAAE,iBAAiB,CAAC;QACvB,CAAC,EAAE,CAAC;QACJ,yEAAyE;QACzE,oEAAoE;QACpE,0EAA0E;QAC1E,0EAA0E;QAC1E,uEAAuE;QACvE,0EAA0E;QAC1E,mEAAmE;QACnE,kBAAkB;QAClB,wEAAwE;QACxE,4EAA4E;QAC5E,yEAAyE;QACzE,kEAAkE;QAClE,2EAA2E;QAC3E,6EAA6E;QAC7E,2DAA2D;QAC3D,yDAAyD;QACzD,yDAAyD;QACzD,KAAK,EAAE,mBAAmB,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,CAAC;QACjE,OAAO,EAAE,mBAAmB;QAC5B,sFAAsF;QACtF,YAAY,EAAE,uBAAuB;QACrC,IAAI,EAAE,mCAAmC;QACzC,MAAM,EAAE,+DAA+D;QACvE,OAAO,EAAE,wCAAwC;QACjD,YAAY,EAAE,CAAC,CAAC,EAAE,MAAM,EAAE,EAAE;YAC1B,MAAM,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,MAAM,CAAC;YACxB,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,KAAK,CAAC;YACxE,OAAO,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC1B,CAAC;QACD,cAAc,EAAE,0DAA0D;QAC1E,cAAc,EAAE,UAAU;QAC1B,UAAU,EAAE;YACV,+DAA+D;YAC/D,0BAA0B;SAC3B;KACF,CAAC;IACF,MAAM,EAAE,aAAa;IACrB,eAAe,EAAE;QACf;YACE,WAAW,EACT,8EAA8E;gBAC9E,wEAAwE;YAC1E,OAAO,EAAE,KAAK;SACf;KACF;IACD,SAAS,EAAE;QACT;YACE,EAAE,EAAE,IAAI;YACR,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EAAE,wEAAwE;SACpF;QACD;YACE,EAAE,EAAE,KAAK;YACT,IAAI,EAAE,SAAS;YACf,IAAI,EAAE,SAAS;YACf,SAAS,EACP,2FAA2F;SAC9F;KACF;IACD,SAAS,EAAE;QACT,4FAA4F;QAC5F,mFAAmF;KACpF;IACD,YAAY,EAAE,UAAU;CACzB,CAAC;AAEF,oFAAoF;AACpF,MAAM,CAAC,MAAM,aAAa,GAA2B,CAAC,kBAAkB,EAAE,kBAAkB,CAAC,CAAC;AAE9F;;;;;;;;;GASG;AACH,MAAM,wBAAwB,GAC5B,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,gBAAgB;AAChB,MAAM,wBAAwB,GAC5B,kBAAkB,CAAC,kBAAkB,CAAC,CAAC;AAEzC,+DAA+D;AAC/D,MAAM,CAAC,MAAM,eAAe,GAAgC;IAC1D,wBAAwB;IACxB,wBAAwB;CACzB,CAAC","sourcesContent":["/**\n * The two Phase 0 APPROXIMATION bridges — design note §3, bridges 3 and 4.\n *\n * - `ab-pendulum-linear` — a REGULAR limit. `θ'' + ω0² sin θ = 0` reduces to\n *   `θ'' + ω0² θ = 0` as `θ0 → 0`. The order of the system does not drop, the\n *   error is `O(θ0²)` in the relative period, and the approximation is\n *   NON-UNIFORM in time: a fixed relative period error accumulates phase\n *   without bound, so the bound carries the horizon `t ≪ 16 T0/θ0²`.\n * - `ab-damped-massless` — a SINGULAR limit. `m x'' + b x' + k x = 0` reduces\n *   to `b x' + k x = 0` as `m → 0`. The order DROPS from two to one, one\n *   initial condition is lost, and the reduced model cannot satisfy `x'(0)`.\n *   The bound therefore holds only OUTSIDE the boundary layer, `t ≥ 5 m/b`.\n *\n * On the position offset of the singular bridge: `x_full − x_red` is an\n * `O((1 + |v0|) m)` deviation that PERSISTS THROUGHOUT THE OUTER REGION and\n * decays only on the slow time scale `b/k`. It is the matched-asymptotics\n * correction, not a boundary-layer transient, and waiting out the fast mode\n * does not remove it. (Design note §9 YELLOW (5) — \"non-vanishing\" was the\n * imprecise wording this replaces: the offset does go to zero as `t → ∞`.)\n *\n * @module atlas/oscillators/bridges-limits\n */\n\nimport { ellipticK } from '@danielsimonjr/mathts-functions';\nimport { numberFrom } from '../../core/mathts-read.js';\nimport { physjsFormalRef } from '../physjs-ref.js';\nimport { makeRegime } from '../regime.js';\nimport { MissingHorizonError } from '../types.js';\nimport type { ApproximationBound, AtlasBridge, RelationContract, Regime } from '../types.js';\nimport { relationContractOf } from './bridges-exact.js';\nimport { DAMPING, SPRING_CONSTANT } from './dimensions.js';\nimport { RELATIVE_PERIOD_NORM } from './norms.js';\nimport { ACCELERATION, LENGTH, MASS } from '../../dimensional/types.js';\n\nconst FAMILY = 'oscillators';\nconst TEST_FILE = 'tests/atlas/oscillators-limits.test.ts';\n\n/**\n * Build an `ApproximationBound`, rejecting one with no prose horizon.\n *\n * Every `approximation` bridge in this module is constructed through here, so\n * the mandatory-horizon rule of `ApproximationBound` has something that FAILS\n * rather than a comment asking for compliance. `horizonHolds` is mandatory by\n * the type; the prose `horizon` is the field a caller can satisfy with `''`,\n * so it is the one checked at run time. `uniformity: null` is accepted: not\n * yet analysed is a real state, and the refusal lives in `boundPath`.\n *\n * @throws MissingHorizonError if `horizon` is empty or whitespace.\n * @internal\n */\nexport function makeApproximation(bound: ApproximationBound): ApproximationBound {\n  if (bound.horizon.trim() === '') {\n    throw new MissingHorizonError(\n      'an approximation bound requires a non-empty prose horizon beside horizonHolds',\n    );\n  }\n  return bound;\n}\n\n/**\n * The EXACT relative period error of the pendulum at amplitude `theta0`:\n * `(2/π) K(sin(θ0/2)) − 1`, with `K` MathTS `ellipticK`, which takes the\n * PARAMETER `m = k² = sin²(θ0/2)` and evaluates it by the AGM\n * (`K = π / (2·AGM(1, √(1−m)))`, so this is `1/AGM(1, cos(θ0/2)) − 1`).\n *\n * This is the machine form of `ab-pendulum-linear`'s `delta`. It is SHARP —\n * the bound equals the error it bounds — and it is strictly increasing in\n * `|θ0|`, so its value at the edge of a declared amplitude range is the\n * supremum over that range. That last property is what lets the record state\n * one scalar for the whole domain honestly.\n *\n * The series `θ0²/16` is the LEADING TERM of this quantity, not an upper bound\n * on it: at the domain edge `θ0 = 0.5` the series gives `0.0156250` while the\n * exact error is `0.0158525311`, so the series UNDERSTATES by 1.456%. A record\n * that declares the series value is violated at its own boundary.\n *\n * @returns `Infinity` outside `|θ0| < π`, where no finite period error exists.\n * @internal\n */\nexport function pendulumPeriodErrorAt(params: Readonly<Record<string, number>>): number {\n  const { theta0 } = params;\n  if (!Number.isFinite(theta0) || Math.abs(theta0) >= Math.PI) return Infinity;\n  const m = Math.sin(theta0 / 2) ** 2;\n  const k = numberFrom('ellipticK', ellipticK(m));\n  return (2 / Math.PI) * k - 1;\n}\n\n/**\n * The position-offset bound of `ab-damped-massless` at a point:\n * `2 (1 + |v0|) m / b`, the formula this record's docstring has always\n * carried. It is a bound ONLY at the witness normalisation `b = k = x0 = 1`:\n * the leading offset is `(m/b)(v0 + k x0/b)`, and the formula has no `k` or\n * `x0` in it. At `m = 1e-3, b = 1, k = 100` (inside `m k/b² < 1/4`) the true\n * error is more than 20 times the formula.\n *\n * @returns `Infinity` unless `m`, `b`, `k`, `x0` and `v0` are all finite and\n *   `b = k = x0 = 1` — a missing parameter would otherwise silently return a\n *   bound for a normalisation the caller did not ask about.\n * @internal\n */\nexport function dampedOffsetBoundAt(params: Readonly<Record<string, number>>): number {\n  const { m, b, k, x0, v0 } = params;\n  if (![m, b, k, x0, v0].every(Number.isFinite)) return Infinity;\n  if (b !== 1 || k !== 1 || x0 !== 1) return Infinity;\n  return (2 * (1 + Math.abs(v0)) * m) / b;\n}\n\n/** Pendulum regime: `θ0 ≤ 0.5 rad`, stated on the dimensionless input. @internal */\nconst PENDULUM_REGIME: Regime = makeRegime(\n  FAMILY,\n  [\n    { name: 'g', dim: ACCELERATION },\n    { name: 'ell', dim: LENGTH },\n  ],\n  ['theta0'],\n  [{ group: 'theta0', op: '<=', bound: 0.5, alias: 'θ0 ≤ 0.5 rad' }],\n);\n\n/**\n * Damped-spring regime: overdamped, `m k / b² < 1/4`.\n *\n * The inequality is written on the π-group, not on `ζ = ½ (m k / b²)^(−½)`,\n * which is a derived quantity and not a group (`RegimeInequality` doc).\n * `m k / b² < 1/4` is exactly `ζ > 1`, recorded as the display alias.\n *\n * @internal\n */\nconst DAMPED_REGIME: Regime = makeRegime(\n  FAMILY,\n  [\n    { name: 'm', dim: MASS },\n    { name: 'b', dim: DAMPING },\n    { name: 'k', dim: SPRING_CONSTANT },\n  ],\n  [],\n  [{ group: 'm · b^-2 · k', op: '<', bound: 0.25, alias: 'ζ > 1' }],\n);\n\n/**\n * `ab-pendulum-linear` — small-angle pendulum → harmonic spring.\n *\n * `θ0²/16` is the leading term of the exact relative period error\n * `(2/π) K(sin(θ0/2)) − 1`. At `θ0 = 0.2` the residual is `5.744e-6`, which\n * the next series term `11 θ0⁴/3072 = 5.729e-6` accounts for to `1.5e-8`.\n *\n * `delta` is the EXACT error at the edge of the declared range, not that\n * series term. The record previously declared `0.5²/16 = 0.0156250`, which is\n * the series evaluated at `θ0 = 0.5` — and the exact error there is\n * `0.0158525311`, so the declared bound was VIOLATED AT ITS OWN BOUNDARY by\n * 1.456%. A truncated series is an approximation of the error, not a bound on\n * it, and the sign of its omitted tail decides which. Here the tail is\n * positive, so the series sits below what it was asked to cover.\n *\n * @internal\n */\nexport const AB_PENDULUM_LINEAR: AtlasBridge = {\n  id: 'ab-pendulum-linear',\n  relation: 'approximation',\n  premises: ['model-pendulum'],\n  conclusion: 'model-spring',\n  transformation: 'sin θ → θ, with x ↔ ℓ θ and ω0² = g/ℓ',\n  preserves: ['harmonic frequency ω0 to O(θ0²)', 'energy conservation', 'time-reversal symmetry'],\n  doesNotPreserve: ['the amplitude dependence of the period', 'phase over times t ≳ 16 T0/θ0²'],\n  sideConditions: ['θ0 ≤ 0.5 rad', 'the bound is a PERIOD error and is not uniform in time'],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over `θ0 ≤ 0.5`: the error is strictly increasing in |θ0|, so the\n    // edge value IS the supremum. 0.0158525311014... (MathTS ellipticK, by the\n    // AGM), independently reproduced by quadrature of the complete elliptic integral.\n    delta: pendulumPeriodErrorAt({ theta0: 0.5 }),\n    deltaAt: pendulumPeriodErrorAt,\n    deltaAtBasis: 'closed-form',\n    norm: RELATIVE_PERIOD_NORM,\n    domain: 'θ0 ≤ 0.5 rad',\n    horizon: 't ≪ 16 T0/θ0²; machine form t < 4 T0/θ0², the π/2-drift time',\n    // ⚠ The machine form is a QUARTER of the prose scale, and it must be.\n    // `16 T0/θ0²` is the time at which the accumulated phase drift reaches a\n    // FULL 2π (measured: 400 T0 at θ0 = 0.2, drift 6.2976 rad) — by then the\n    // approximation has lapped, so `t < 16 T0/θ0²` is not a horizon at all.\n    // The `≪` is what carries the whole claim, and `<` does not render it.\n    // `4 T0/θ0²` is the π/2-drift time, which is exactly what W7b measures\n    // independently (99.77 cycles at θ0 = 0.2). It is also what makes\n    // `horizonHolds(200 T0)` false and `horizonHolds(10 T0)` true.\n    horizonHolds: (t, params) => {\n      const { T0, theta0 } = params;\n      if (!Number.isFinite(T0) || !Number.isFinite(theta0) || theta0 === 0) return false;\n      return t < (4 * T0) / theta0 ** 2;\n    },\n    parameterRange: 'θ0 ≤ 0.5 rad',\n    limitCharacter: 'regular',\n    // A period error. The record says it is not uniform in time, so time is\n    // not listed.\n    uniformity: ['one period, for θ0 in the stated domain'],\n  }),\n  regime: PENDULUM_REGIME,\n  counterexamples: [\n    {\n      description:\n        'the approximation is non-uniform in time: at θ0 = 0.2 the phase drift reaches π/2 ' +\n        'after ~100 cycles, however small the per-cycle error is',\n      witness: 'W7b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W7',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'T/T0 − 1 ∈ [0.002505, 0.002507]; residual ∈ [5.70e-6, 5.76e-6] at θ0 = 0.2',\n    },\n    {\n      id: 'W7b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'cycles to π/2 drift ∈ [99, 101] at θ0 = 0.2',\n    },\n    {\n      id: 'W7c',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: 'RK4 zero-crossing lag within 0.05° of the elliptic prediction 89.98°',\n    },\n  ],\n  citations: [\n    'Landau & Lifshitz, Mechanics §11 (pendulum period as a complete elliptic integral)',\n    'Abramowitz & Stegun §17.6 (AGM evaluation of K)',\n  ],\n  reviewStatus: 'proposed',\n  // Milestone 1. The reference is the PhysJS theorem, which imports Physlib's\n  // `ClassicalMechanics.SimplePendulum.linearizedEquationOfMotion_iff`. It\n  // certifies the TRANSFORMATION (mass mℓ², spring constant mgℓ, ω0² = g/ℓ),\n  // not `bound.delta`. Physlib's period results concern `periodFormula`, which\n  // its own TODO has not yet tied to the motion. The sanity lemmas instantiate\n  // that imported statement. Keyed by `formal/lean/manifest.json`.\n  formalRef: physjsFormalRef('ab-pendulum-linear'),\n};\n\n/**\n * `ab-damped-massless` — the SINGULAR `m → 0` limit of the damped spring.\n *\n * `deltaAt = 2(1 + |v0|) m/b` bounds the position offset outside the boundary\n * layer, and `delta` is its supremum over the declared range rather than its\n * value at one fixture. The lost initial condition is visible in the VELOCITY: inside the\n * layer the reduced model's `x'` is wrong by O(1/m·…), and it is only after\n * `t ≈ 5 m/b` that the fast mode has decayed.\n *\n * @internal\n */\nexport const AB_DAMPED_MASSLESS: AtlasBridge = {\n  id: 'ab-damped-massless',\n  relation: 'approximation',\n  premises: ['model-damped-spring'],\n  conclusion: 'model-first-order',\n  transformation: 'm → 0, dropping the m x″ term',\n  preserves: ['the slow relaxation rate k/b to O(m)', 'the sign and monotonicity of the decay'],\n  doesNotPreserve: [\n    'the order of the system (two → one)',\n    \"the initial condition x'(0), which the reduced model cannot satisfy\",\n    'the fast mode r ≈ −b/m',\n  ],\n  sideConditions: [\n    'overdamped: m k / b² < 1/4',\n    'the bound holds only outside the boundary layer, t ≥ 5 m/b',\n  ],\n  bound: makeApproximation({\n    K: 1,\n    // Sup over the DECLARED range, not the fixture. `2(1+|v0|)m/b` increases\n    // in both m and |v0|, and the declared range is `m k/b² < 1/4` with\n    // `|v0| ≤ 5`; at the witness normalisation `b = k = 1` that is `m < 1/4`,\n    // so the supremum is `2·(1+5)·(1/4) = 3`, approached at the open edge and\n    // not attained. The record previously declared `2·(1+5)·1e-3 = 0.012`,\n    // the same formula frozen at the ONE fixture mass `m = 1e-3` — true error\n    // there is 5.96e-3, but at m = 0.24 it is 5.19e-1, forty times the\n    // declared bound.\n    // KNOWN LOOSENESS (0.47.0 persona re-test, L9 residual; documented, not\n    // tightened): for x0 = 1 the measured sup of |x − x_reduced| over the whole\n    // declared range (m < 1/4, |v0| ≤ 5, t ≥ 5m/b) is 0.524, at m ≈ 0.21 and\n    // v0 = +5, and it does not grow toward critical damping (0.514 at\n    // m = 0.249999). So δ = 3 holds but is 5.7× loose, and it exceeds the unit\n    // signal it bounds. A tighter δ needs a bound with its own basis, not a grid\n    // maximum; until then `upt path` withholds the point bound\n    // (deltaAtBasis below). Independent check: the persona's\n    // checks/recheck_047.py and a 115 × 2 grid over (m, ±5).\n    delta: dampedOffsetBoundAt({ m: 0.25, b: 1, k: 1, x0: 1, v0: 5 }),\n    deltaAt: dampedOffsetBoundAt,\n    // Witness W8b supports it numerically for m ∈ {1e-1, 1e-2, 1e-3}; no proof covers it.\n    deltaAtBasis: 'numerically-supported',\n    norm: 'sup |x − x_reduced| for t ≥ 5 m/b',\n    domain: 't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n    horizon: 't ≥ 5 m/b (outside the boundary layer)',\n    horizonHolds: (t, params) => {\n      const { m, b } = params;\n      if (!Number.isFinite(m) || !Number.isFinite(b) || b === 0) return false;\n      return t >= (5 * m) / b;\n    },\n    parameterRange: 'm k / b² < 1/4, |v0| ≤ 5; with b = k = 1 this is m < 1/4',\n    limitCharacter: 'singular',\n    uniformity: [\n      't ≥ 5 m/b, overdamped, at the witness normalisation b = k = 1',\n      'm k / b² < 1/4, |v0| ≤ 5',\n    ],\n  }),\n  regime: DAMPED_REGIME,\n  counterexamples: [\n    {\n      description:\n        \"inside the boundary layer the reduced model's velocity is wrong by O(1): at \" +\n        \"t = 0.5 m/b with v0 = 5 the velocity error is 3.6, against x'_red = −1\",\n      witness: 'W8b',\n    },\n  ],\n  witnesses: [\n    {\n      id: 'W8',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance: '|r_slow + k/b| < 2m and |r_fast·m + b| < 2m for m ∈ {1e-1, 1e-2, 1e-3}',\n    },\n    {\n      id: 'W8b',\n      kind: 'numeric',\n      test: TEST_FILE,\n      tolerance:\n        \"|x_full − x_red| < 2(1+v0)·m for t ≥ 5 m/b; |x'| error > 1 at 0.5 m/b and < 0.05 at 5 m/b\",\n    },\n  ],\n  citations: [\n    \"Kevorkian & Cole, Multiple Scale and Singular Perturbation Methods §2 (Tikhonov's theorem)\",\n    'Verhulst, Methods and Applications of Singular Perturbations §8 (boundary layers)',\n  ],\n  reviewStatus: 'proposed',\n};\n\n/** The two approximation bridges of this module, in design-note order. @internal */\nexport const LIMIT_BRIDGES: readonly AtlasBridge[] = [AB_PENDULUM_LINEAR, AB_DAMPED_MASSLESS];\n\n/**\n * The two approximation bridges re-registered through `RelationContract`.\n *\n * The union REQUIRES `bound` on an `approximation`, and `ApproximationBound`\n * requires its own machine `horizonHolds`, so these two lines are what turn\n * Sprint 0's doc-comment requirement into a compile-and-throw one. Both bounds\n * are the records' own, measured ones — nothing is restated here.\n *\n * @internal\n */\nconst CONTRACT_PENDULUM_LINEAR: RelationContract =\n  relationContractOf(AB_PENDULUM_LINEAR);\n\n/** @internal */\nconst CONTRACT_DAMPED_MASSLESS: RelationContract =\n  relationContractOf(AB_DAMPED_MASSLESS);\n\n/** The two limit contracts, in design-note order. @internal */\nexport const LIMIT_CONTRACTS: readonly RelationContract[] = [\n  CONTRACT_PENDULUM_LINEAR,\n  CONTRACT_DAMPED_MASSLESS,\n];\n"]}