universal-physics-tensor 10.0.0 → 11.0.0

package/dist/cli/commands/metric.js.map

@@ -1 +1 @@
-{"version":3,"file":"metric.js","sourceRoot":"","sources":["../../../src/cli/commands/metric.ts"],"names":[],"mappings":"AAUA,OAAO,EAAE,eAAe,EAAiC,MAAM,eAAe,CAAC;AAC/E,OAAO,EAAE,WAAW,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,cAAc,CAAC;AACxC,OAAO,EAAE,QAAQ,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AACpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,gBAAgB,CAAC;AAGlD,MAAM,KAAK,GAAe;IACxB,SAAS;IACT;QACE,IAAI,EAAE,YAAY;QAClB,UAAU,EAAE,MAAM;QAClB,WAAW,EAAE,oFAAoF;KAClG;CACF,CAAC;AAEF,MAAM,KAAK,GAAG,CAAC,WAAW,EAAE,eAAe,EAAE,MAAM,EAAE,MAAM,CAAU,CAAC;AAEtE,MAAM,IAAI,GAAG;;;;;;;;;;;;;;;;;;;wDAmB2C,CAAC;AAEzD,SAAS,QAAQ,CAAC,CAAqB;IACrC,OAAQ,KAA2B,CAAC,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;AACxD,CAAC;AAED;;;;;;GAMG;AACH,SAAS,gBAAgB,CAAC,GAAsB,EAAE,CAAU;IAC1D,IAAI,CAAC,YAAY,GAAG,CAAC,eAAe,IAAI,CAAC,YAAY,GAAG,CAAC,SAAS;QAAE,OAAO,IAAI,QAAQ,CAAC,eAAe,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;IACpH,IAAI,CAAC,YAAY,KAAK,IAAI,CAAC,CAAC,WAAW,KAAK,KAAK;QAAE,OAAO,IAAI,QAAQ,CAAC,eAAe,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;IACnG,OAAO,CAAC,CAAC;AACX,CAAC;AAED,KAAK,UAAU,GAAG,CAAC,GAAe;IAChC,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;IACjC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,2BAA2B,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,4BAA4B,CACxE,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACxC,sFAAsF;IACtF,gBAAgB,CAAC,QAAQ,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;IAC/C,IAAI,MAAM,CAAC;IACX,IAAI,CAAC;QACH,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IAC5C,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,MAAM,gBAAgB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC;IACD,IAAI,QAmBS,CAAC;IACd,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,IAAI,CAAC;YACH,IAAI,IAAI,KAAK,eAAe,EAAE,CAAC;gBAC7B,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,KAAK,GAAG,GAAG,CAAC,0BAA0B,CAAC;oBAC3C,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;oBACjC,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;iBAClC,CAAC,CAAC;gBACH,QAAQ,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,GAAG,KAAK,EAAE,CAAC;YAC5C,CAAC;iBAAM,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;gBAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC;gBACzC,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC;gBACnC,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,KAAK,IAAI,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;gBACrD,IAAI,KAAK,KAAK,SAAS,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;oBAAE,MAAM,IAAI,QAAQ,CAAC,qCAAqC,CAAC,CAAC;gBACnG,MAAM,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;gBACzC,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;gBACzB,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;oBAAE,MAAM,IAAI,QAAQ,CAAC,iCAAiC,CAAC,CAAC;gBACpF,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;gBACxD,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;oBACZ,wFAAwF;oBACxF,uFAAuF;oBACvF,sFAAsF;oBACtF,uFAAuF;oBACvF,uFAAuF;oBACvF,IAAI,CAAC,UAAU,EAAE,CAAC;wBAChB,MAAM,IAAI,QAAQ,CAChB,oKAAoK,CACrK,CAAC;oBACJ,CAAC;oBACD,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;oBAC9B,QAAQ,GAAG;wBACT,IAAI,EAAE,UAAmB;wBACzB,IAAI,EAAE,mHAAmH;wBACzH,GAAG,GAAG,CAAC,0BAA0B,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;qBAC5G,CAAC;gBACJ,CAAC;qBAAM,CAAC;oBACR,MAAM,MAAM,GAAG,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC;oBACxE,QAAQ;wBACN,UAAU;4BACR,CAAC,CAAC,EAAE,IAAI,EAAE,UAAmB,EAAE,GAAG,GAAG,CAAC,sBAAsB,CAAC,MAAM,CAAC,EAAE;4BACtE,CAAC,CAAC;gCACE,IAAI,EAAE,UAAmB;gCACzB,GAAG,GAAG,CAAC,YAAY,CAAC;oCAClB,GAAG,MAAM;oCACT,KAAK;oCACL,GAAG,GAAG,CAAC,qBAAqB,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC;oCACjE,GAAG,EAAE,CAAC;iCACP,CAAC;6BACH,CAAC;gBACR,CAAC;YACH,CAAC;iBAAM,CAAC;gBACN,MAAM,IAAI,UAAU,CAAC,sDAAsD,CAAC,CAAC;YAC/E,CAAC;QACH,CAAC;QAAC,OAAO,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,YAAY,UAAU,IAAI,CAAC,YAAY,QAAQ;gBAAE,MAAM,CAAC,CAAC;YAC9D,MAAM,gBAAgB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,QAAQ,CAAC,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,GAAG,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC;QACnH,OAAO,CAAC,CAAC;IACX,CAAC;IACD,GAAG,CAAC,KAAK,MAAM,CAAC,MAAM,eAAe,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC;IACzD,GAAG,CAAC,KAAK,MAAM,CAAC,aAAa,EAAE,CAAC,CAAC;IACjC,GAAG,CAAC,iBAAiB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAC3C,GAAG,CAAC,gBAAgB,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACjG,GAAG,CAAC,WAAW,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACvF,KAAK,MAAM,IAAI,IAAI,MAAM,CAAC,KAAK;QAAE,GAAG,CAAC,KAAK,IAAI,EAAE,CAAC,CAAC;IAClD,GAAG,CAAC,0BAA0B,CAAC,CAAC;IAChC,IAAI,MAAM,CAAC,WAAW,CAAC,MAAM,KAAK,CAAC;QAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACvD,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,WAAW;QAAE,GAAG,CAAC,OAAO,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;IACvE,GAAG,CAAC,oBAAoB,CAAC,CAAC;IAC1B,IAAI,MAAM,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACjD,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,KAAK;QAAE,GAAG,CAAC,OAAO,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;IACjE,GAAG,CAAC,sBAAsB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAChD,GAAG,CAAC,qBAAqB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAC/C,GAAG,CAAC,iBAAiB,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAClG,IAAI,QAAQ,EAAE,CAAC;QACb,MAAM,IAAI,GACR,QAAQ,CAAC,EAAE,KAAK,SAAS;YACvB,CAAC,CAAC,EAAE;YACJ,CAAC,CAAC,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,SAAS,QAAQ,CAAC,KAAK,IAAI,QAAQ,CAAC,OAAO,EAAE,CAAC;QAC1J,MAAM,KAAK,GACT,QAAQ,CAAC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,QAAQ,CAAC,MAAM,IAAI,QAAQ,CAAC,QAAQ,EAAE,CAAC;QACpF,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,KAAK,UAAU,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,gBAAgB,CAAC;QACjF,GAAG,CACD,eAAe,KAAK,OAAO,QAAQ,CAAC,EAAE,SAAS,QAAQ,CAAC,IAAI,GAAG,KAAK,OAAO,QAAQ,CAAC,UAAU,SAAS,QAAQ,CAAC,KAAK,SAAS,IAAI,EAAE,CACrI,CAAC;QACF,IAAI,QAAQ,CAAC,IAAI,KAAK,SAAS;YAAE,GAAG,CAAC,OAAO,QAAQ,CAAC,IAAI,EAAE,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,CAAC,MAAM,OAAO,GAAY;IAC9B,IAAI,EAAE,QAAQ;IACd,OAAO,EAAE,CAAC,WAAW,CAAC;IACtB,KAAK,EAAE,KAAK;IACZ,IAAI,EAAE,WAAW,CAAC,IAAI,EAAE,KAAK,CAAC;IAC9B,OAAO,EAAE,uEAAuE;IAChF,OAAO,EAAE,mDAAmD;IAC5D,KAAK,EAAE,UAAU;IACjB,GAAG;CACJ,CAAC;AAEF,eAAe,CAAC,OAAO,CAAC,CAAC","sourcesContent":["/**\n * `upt metric` — Christoffel symbols and curvature scalars for a named metric.\n *\n * `--geodesic` integrates a short Schwarzschild circular orbit, or a Kerr\n * geodesic. θ = π/2 is equatorial and circular. Any other θ is an inclined\n * spherical orbit with that polar turning point.\n *\n * @module cli/commands/metric\n */\nimport type { FlagSpec } from '../args.js';\nimport { registerCommand, type Command, type CommandCtx } from '../command.js';\nimport { commandHelp, JSON_FLAG } from '../flag-help.js';\nimport { emitJson } from '../output.js';\nimport { CliError, UsageError } from '../errors.js';\nimport { splitAssignments } from '../bindings.js';\nimport type { MetricId } from '../../cli-api.js';\n\nconst FLAGS: FlagSpec[] = [\n  JSON_FLAG,\n  {\n    name: '--geodesic',\n    valueStyle: 'none',\n    description: 'Integrate a short Schwarzschild circular orbit, or a Kerr geodesic at the given θ.',\n  },\n];\n\nconst NAMES = ['minkowski', 'schwarzschild', 'flrw', 'kerr'] as const;\n\nconst HELP = `upt metric <minkowski|schwarzschild|flrw|kerr> [key=value ...] [--geodesic] [--json]\n        Christoffel symbols, the Ricci tensor, the Ricci scalar and the\n        Kretschmann scalar of one exact metric. Alias: upt curvature.\n        Signature of the line element is (-,+,+,+), the same mostly-plus\n        signature as the canonical Einstein-equation metric node.\n        Schwarzschild Kretschmann is checked against 48 G² M² / (c⁴ r⁶).\n        FLRW prints the Friedmann equation including the curvature term\n        −k c²/a² and Λ. Kerr Kretschmann uses the Boyer–Lindquist closed form.\n        --geodesic integrates a short Schwarzschild circular orbit, or a Kerr\n        geodesic (equatorial and circular at θ = π/2; inclined, with that\n        polar turning point, at any other θ). Kerr's a is a length in metres;\n        the report prints a/M = a/(GM/c²). A Ricci scalar that has a closed\n        form is that form; a vacuum finite-difference residual is not listed\n        as Ricci. A parameter is a number, a unit (M=1Msun, r=1km) or a\n        constant expression (theta=pi/2). A bare number is already in the\n        parameter's SI unit; a dimensioned quantity must match it.\n        \\`upt help metric\\` describes every flag.\n        e.g.  upt metric schwarzschild M=1Msun r=1e8 theta=pi/2\n              upt metric flrw k=1 t=2\n              upt curvature kerr a=1000 r=1e8 theta=1.2`;\n\nfunction isMetric(s: string | undefined): s is MetricId {\n  return (NAMES as readonly string[]).includes(s ?? '');\n}\n\n/**\n * The metric library refuses a bad value (a non-positive mass, a point inside\n * the horizon, an unknown parameter, a start the integrator cannot take, a\n * parameter in the wrong unit) with `MetricMassError`, a `UnitError`, or a\n * plain `Error`. Each is a bad value, exit 1. Any other class (a TypeError, a\n * RangeError) is a defect and is not converted.\n */\nfunction metricValueError(api: CommandCtx['api'], e: unknown): unknown {\n  if (e instanceof api.MetricMassError || e instanceof api.UnitError) return new CliError(`upt metric: ${e.message}`);\n  if (e instanceof Error && e.constructor === Error) return new CliError(`upt metric: ${e.message}`);\n  return e;\n}\n\nasync function run(ctx: CommandCtx): Promise<number> {\n  const { args, api, out } = ctx;\n  const name = args.positionals[0];\n  if (!isMetric(name)) {\n    throw new UsageError(\n      `upt metric needs one of ${NAMES.join(', ')}. See \\`upt help metric\\`.`,\n    );\n  }\n  const pairs = args.positionals.slice(1);\n  // A token with no `=` is usage (exit 2); a parameter given twice is refused (exit 1).\n  splitAssignments('metric', pairs, 'key=value');\n  let report;\n  try {\n    report = api.curvatureReport(name, pairs);\n  } catch (e) {\n    throw metricValueError(api, e);\n  }\n  let geodesic:\n    | {\n        readonly kind: 'circular' | 'inclined';\n        readonly r0: number;\n        readonly rEnd: number;\n        readonly theta0?: number;\n        readonly thetaEnd?: number;\n        readonly phiAdvance: number;\n        readonly steps: number;\n        readonly E0?: number;\n        readonly EEnd?: number;\n        readonly L0?: number;\n        readonly LEnd?: number;\n        readonly Q0?: number;\n        readonly QEnd?: number;\n        readonly norm0?: number;\n        readonly normEnd?: number;\n        readonly note?: string;\n      }\n    | undefined;\n  if (args.flags.has('geodesic')) {\n    try {\n      if (name === 'schwarzschild') {\n        const M = report.parameters.M;\n        const r = report.parameters.r;\n        const orbit = api.schwarzschildCircularOrbit({\n          ...(M === undefined ? {} : { M }),\n          ...(r === undefined ? {} : { r }),\n        });\n        geodesic = { kind: 'circular', ...orbit };\n      } else if (name === 'kerr') {\n        const Mgeom = report.parameters.M_geom_m;\n        const a = report.parameters.a ?? 0;\n        const r = report.parameters.r;\n        const theta = report.parameters.theta ?? Math.PI / 2;\n        if (Mgeom === undefined || !(Mgeom > 0)) throw new CliError('Kerr geodesic needs a positive mass');\n        const rOverM = (r ?? 10 * Mgeom) / Mgeom;\n        const aOverM = a / Mgeom;\n        if (!(Math.abs(aOverM) <= 1)) throw new CliError('Kerr geodesic wants |a| ≤ GM/c²');\n        const equatorial = Math.abs(theta - Math.PI / 2) < 1e-6;\n        if (a === 0) {\n          // At a = 0 the Kerr metric is Schwarzschild. The Boyer–Lindquist integrator refuses the\n          // circular start there (\"radial potential is negative at the start\"): the two terms of\n          // its radial potential cancel to below its tolerance at a = 0, and only at a = 0. The\n          // equatorial circular orbit is the Schwarzschild one and is integrated as that, so the\n          // default parameters (a = 0) run. An inclined start has no Schwarzschild routine here.\n          if (!equatorial) {\n            throw new CliError(\n              'upt metric: at a = 0 an inclined Kerr orbit is a tilted Schwarzschild circular orbit, and this integrator needs a ≠ 0 for an inclined start; give a, or theta=pi/2',\n            );\n          }\n          const M = report.parameters.M;\n          geodesic = {\n            kind: 'circular' as const,\n            note: 'a = 0: the Kerr metric is Schwarzschild; the equatorial circular orbit is integrated by the Schwarzschild routine',\n            ...api.schwarzschildCircularOrbit({ ...(M === undefined ? {} : { M }), ...(r === undefined ? {} : { r }) }),\n          };\n        } else {\n        const shared = { M: Mgeom, aOverM, rOverM, fraction: 0.005, steps: 40 };\n        geodesic =\n          equatorial\n            ? { kind: 'circular' as const, ...api.kerrEquatorialCircular(shared) }\n            : {\n                kind: 'inclined' as const,\n                ...api.kerrGeodesic({\n                  ...shared,\n                  theta,\n                  ...api.kerrTurningPointOrbit({ M: Mgeom, aOverM, rOverM, theta }),\n                  mu2: 1,\n                }),\n              };\n        }\n      } else {\n        throw new UsageError('upt metric: --geodesic is for schwarzschild or kerr.');\n      }\n    } catch (e) {\n      if (e instanceof UsageError || e instanceof CliError) throw e;\n      throw metricValueError(api, e);\n    }\n  }\n  if (args.flags.has('json')) {\n    emitJson({ command: 'metric', result: { ...report, ...(geodesic === undefined ? {} : { geodesic }) } }, ctx.write);\n    return 0;\n  }\n  out(`\\n${report.metric}  signature ${report.signature}`);\n  out(`  ${report.signatureNote}`);\n  out(`  coordinates ${report.coordinates}`);\n  out('  parameters: ' + Object.entries(report.parameters).map(([k, v]) => `${k}=${v}`).join(' '));\n  out('  point: ' + Object.entries(report.point).map(([k, v]) => `${k}=${v}`).join(' '));\n  for (const note of report.notes) out(`  ${note}`);\n  out('  Christoffel (nonzero):');\n  if (report.christoffel.length === 0) out('    (none)');\n  for (const c of report.christoffel) out(`    ${c.index} = ${c.value}`);\n  out('  Ricci (nonzero):');\n  if (report.ricci.length === 0) out('    (none)');\n  for (const c of report.ricci) out(`    ${c.index} = ${c.value}`);\n  out(`  Ricci scalar R = ${report.ricciScalar}`);\n  out(`  Kretschmann K = ${report.kretschmann}`);\n  out('  closed form: ' + Object.entries(report.closedForm).map(([k, v]) => `${k}=${v}`).join(' '));\n  if (geodesic) {\n    const cons =\n      geodesic.E0 === undefined\n        ? ''\n        : ` E ${geodesic.E0}→${geodesic.EEnd} L ${geodesic.L0}→${geodesic.LEnd} Q ${geodesic.Q0}→${geodesic.QEnd} norm ${geodesic.norm0}→${geodesic.normEnd}`;\n    const polar =\n      geodesic.theta0 === undefined ? '' : ` θ ${geodesic.theta0}→${geodesic.thetaEnd}`;\n    const label = geodesic.kind === 'inclined' ? 'inclined orbit' : 'circular orbit';\n    out(\n      `  geodesic: ${label} r0=${geodesic.r0} rEnd=${geodesic.rEnd}${polar} Δφ=${geodesic.phiAdvance} over ${geodesic.steps} steps${cons}`,\n    );\n    if (geodesic.note !== undefined) out(`    ${geodesic.note}`);\n  }\n  return 0;\n}\n\nexport const command: Command = {\n  name: 'metric',\n  aliases: ['curvature'],\n  flags: FLAGS,\n  help: commandHelp(HELP, FLAGS),\n  summary: 'Print Christoffel symbols and curvature scalars for one exact metric.',\n  example: 'upt metric schwarzschild M=1Msun r=1e8 theta=pi/2',\n  group: 'evaluate',\n  run,\n};\n\nregisterCommand(command);\n"]}
+{"version":3,"file":"metric.js","sourceRoot":"","sources":["../../../src/cli/commands/metric.ts"],"names":[],"mappings":"AAUA,OAAO,EAAE,eAAe,EAAiC,MAAM,eAAe,CAAC;AAC/E,OAAO,EAAE,WAAW,EAAE,SAAS,EAAE,MAAM,iBAAiB,CAAC;AACzD,OAAO,EAAE,QAAQ,EAAE,MAAM,cAAc,CAAC;AACxC,OAAO,EAAE,QAAQ,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AACpD,OAAO,EAAE,gBAAgB,EAAE,MAAM,gBAAgB,CAAC;AAGlD,MAAM,KAAK,GAAe;IACxB,SAAS;IACT;QACE,IAAI,EAAE,YAAY;QAClB,UAAU,EAAE,MAAM;QAClB,WAAW,EAAE,oFAAoF;KAClG;CACF,CAAC;AAEF,MAAM,KAAK,GAAG,CAAC,WAAW,EAAE,eAAe,EAAE,MAAM,EAAE,MAAM,CAAU,CAAC;AAEtE,MAAM,IAAI,GAAG;;;;;;;;;;;;;;;;;;;wDAmB2C,CAAC;AAEzD,SAAS,QAAQ,CAAC,CAAqB;IACrC,OAAQ,KAA2B,CAAC,QAAQ,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;AACxD,CAAC;AAED;;;;;;;GAOG;AACH,SAAS,gBAAgB,CAAC,GAAsB,EAAE,CAAU;IAC1D,IAAI,CAAC,YAAY,GAAG,CAAC,eAAe,IAAI,CAAC,YAAY,GAAG,CAAC,gBAAgB,IAAI,CAAC,YAAY,GAAG,CAAC,SAAS,EAAE,CAAC;QACxG,OAAO,IAAI,QAAQ,CAAC,eAAe,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;IAClD,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,KAAK,UAAU,GAAG,CAAC,GAAe;IAChC,MAAM,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC;IAC/B,MAAM,IAAI,GAAG,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;IACjC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,MAAM,IAAI,UAAU,CAClB,2BAA2B,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,4BAA4B,CACxE,CAAC;IACJ,CAAC;IACD,MAAM,KAAK,GAAG,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACxC,sFAAsF;IACtF,gBAAgB,CAAC,QAAQ,EAAE,KAAK,EAAE,WAAW,CAAC,CAAC;IAC/C,IAAI,MAAM,CAAC;IACX,IAAI,CAAC;QACH,MAAM,GAAG,GAAG,CAAC,eAAe,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;IAC5C,CAAC;IAAC,OAAO,CAAC,EAAE,CAAC;QACX,MAAM,gBAAgB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;IACjC,CAAC;IACD,IAAI,QAmBS,CAAC;IACd,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,CAAC;QAC/B,IAAI,CAAC;YACH,IAAI,IAAI,KAAK,eAAe,EAAE,CAAC;gBAC7B,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,KAAK,GAAG,GAAG,CAAC,0BAA0B,CAAC;oBAC3C,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;oBACjC,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC;iBAClC,CAAC,CAAC;gBACH,QAAQ,GAAG,EAAE,IAAI,EAAE,UAAU,EAAE,GAAG,KAAK,EAAE,CAAC;YAC5C,CAAC;iBAAM,IAAI,IAAI,KAAK,MAAM,EAAE,CAAC;gBAC3B,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,QAAQ,CAAC;gBACzC,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC;gBACnC,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;gBAC9B,MAAM,KAAK,GAAG,MAAM,CAAC,UAAU,CAAC,KAAK,IAAI,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;gBACrD,oFAAoF;gBACpF,4EAA4E;gBAC5E,IAAI,KAAK,KAAK,SAAS,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;oBAAE,MAAM,IAAI,QAAQ,CAAC,iDAAiD,CAAC,CAAC;gBAC/G,MAAM,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC;gBACzC,MAAM,MAAM,GAAG,CAAC,GAAG,KAAK,CAAC;gBACzB,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;oBAAE,MAAM,IAAI,QAAQ,CAAC,6CAA6C,CAAC,CAAC;gBAChG,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;gBACxD,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;oBACZ,wFAAwF;oBACxF,uFAAuF;oBACvF,sFAAsF;oBACtF,uFAAuF;oBACvF,uFAAuF;oBACvF,IAAI,CAAC,UAAU,EAAE,CAAC;wBAChB,MAAM,IAAI,QAAQ,CAChB,oKAAoK,CACrK,CAAC;oBACJ,CAAC;oBACD,MAAM,CAAC,GAAG,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC;oBAC9B,QAAQ,GAAG;wBACT,IAAI,EAAE,UAAmB;wBACzB,IAAI,EAAE,mHAAmH;wBACzH,GAAG,GAAG,CAAC,0BAA0B,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;qBAC5G,CAAC;gBACJ,CAAC;qBAAM,CAAC;oBACR,MAAM,MAAM,GAAG,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,KAAK,EAAE,KAAK,EAAE,EAAE,EAAE,CAAC;oBACxE,QAAQ;wBACN,UAAU;4BACR,CAAC,CAAC,EAAE,IAAI,EAAE,UAAmB,EAAE,GAAG,GAAG,CAAC,sBAAsB,CAAC,MAAM,CAAC,EAAE;4BACtE,CAAC,CAAC;gCACE,IAAI,EAAE,UAAmB;gCACzB,GAAG,GAAG,CAAC,YAAY,CAAC;oCAClB,GAAG,MAAM;oCACT,KAAK;oCACL,GAAG,GAAG,CAAC,qBAAqB,CAAC,EAAE,CAAC,EAAE,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,KAAK,EAAE,CAAC;oCACjE,GAAG,EAAE,CAAC;iCACP,CAAC;6BACH,CAAC;gBACR,CAAC;YACH,CAAC;iBAAM,CAAC;gBACN,MAAM,IAAI,UAAU,CAAC,sDAAsD,CAAC,CAAC;YAC/E,CAAC;QACH,CAAC;QAAC,OAAO,CAAC,EAAE,CAAC;YACX,IAAI,CAAC,YAAY,UAAU,IAAI,CAAC,YAAY,QAAQ;gBAAE,MAAM,CAAC,CAAC;YAC9D,MAAM,gBAAgB,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACjC,CAAC;IACH,CAAC;IACD,IAAI,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,MAAM,CAAC,EAAE,CAAC;QAC3B,QAAQ,CAAC,EAAE,OAAO,EAAE,QAAQ,EAAE,MAAM,EAAE,EAAE,GAAG,MAAM,EAAE,GAAG,CAAC,QAAQ,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,EAAE,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC;QACnH,OAAO,CAAC,CAAC;IACX,CAAC;IACD,GAAG,CAAC,KAAK,MAAM,CAAC,MAAM,eAAe,MAAM,CAAC,SAAS,EAAE,CAAC,CAAC;IACzD,GAAG,CAAC,KAAK,MAAM,CAAC,aAAa,EAAE,CAAC,CAAC;IACjC,GAAG,CAAC,iBAAiB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAC3C,GAAG,CAAC,gBAAgB,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACjG,GAAG,CAAC,WAAW,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IACvF,KAAK,MAAM,IAAI,IAAI,MAAM,CAAC,KAAK;QAAE,GAAG,CAAC,KAAK,IAAI,EAAE,CAAC,CAAC;IAClD,GAAG,CAAC,0BAA0B,CAAC,CAAC;IAChC,IAAI,MAAM,CAAC,WAAW,CAAC,MAAM,KAAK,CAAC;QAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACvD,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,WAAW;QAAE,GAAG,CAAC,OAAO,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;IACvE,GAAG,CAAC,oBAAoB,CAAC,CAAC;IAC1B,IAAI,MAAM,CAAC,KAAK,CAAC,MAAM,KAAK,CAAC;QAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACjD,KAAK,MAAM,CAAC,IAAI,MAAM,CAAC,KAAK;QAAE,GAAG,CAAC,OAAO,CAAC,CAAC,KAAK,MAAM,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;IACjE,GAAG,CAAC,sBAAsB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAChD,GAAG,CAAC,qBAAqB,MAAM,CAAC,WAAW,EAAE,CAAC,CAAC;IAC/C,GAAG,CAAC,iBAAiB,GAAG,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;IAClG,IAAI,QAAQ,EAAE,CAAC;QACb,MAAM,IAAI,GACR,QAAQ,CAAC,EAAE,KAAK,SAAS;YACvB,CAAC,CAAC,EAAE;YACJ,CAAC,CAAC,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,MAAM,QAAQ,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,SAAS,QAAQ,CAAC,KAAK,IAAI,QAAQ,CAAC,OAAO,EAAE,CAAC;QAC1J,MAAM,KAAK,GACT,QAAQ,CAAC,MAAM,KAAK,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,QAAQ,CAAC,MAAM,IAAI,QAAQ,CAAC,QAAQ,EAAE,CAAC;QACpF,MAAM,KAAK,GAAG,QAAQ,CAAC,IAAI,KAAK,UAAU,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,CAAC,gBAAgB,CAAC;QACjF,GAAG,CACD,eAAe,KAAK,OAAO,QAAQ,CAAC,EAAE,SAAS,QAAQ,CAAC,IAAI,GAAG,KAAK,OAAO,QAAQ,CAAC,UAAU,SAAS,QAAQ,CAAC,KAAK,SAAS,IAAI,EAAE,CACrI,CAAC;QACF,IAAI,QAAQ,CAAC,IAAI,KAAK,SAAS;YAAE,GAAG,CAAC,OAAO,QAAQ,CAAC,IAAI,EAAE,CAAC,CAAC;IAC/D,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED,MAAM,CAAC,MAAM,OAAO,GAAY;IAC9B,IAAI,EAAE,QAAQ;IACd,OAAO,EAAE,CAAC,WAAW,CAAC;IACtB,WAAW,EAAE,EAAE,GAAG,EAAE,QAAQ,EAAE;IAC9B,KAAK,EAAE,KAAK;IACZ,IAAI,EAAE,WAAW,CAAC,IAAI,EAAE,KAAK,CAAC;IAC9B,OAAO,EAAE,uEAAuE;IAChF,OAAO,EAAE,mDAAmD;IAC5D,KAAK,EAAE,UAAU;IACjB,GAAG;CACJ,CAAC;AAEF,eAAe,CAAC,OAAO,CAAC,CAAC","sourcesContent":["/**\n * `upt metric` — Christoffel symbols and curvature scalars for a named metric.\n *\n * `--geodesic` integrates a short Schwarzschild circular orbit, or a Kerr\n * geodesic. θ = π/2 is equatorial and circular. Any other θ is an inclined\n * spherical orbit with that polar turning point.\n *\n * @module cli/commands/metric\n */\nimport type { FlagSpec } from '../args.js';\nimport { registerCommand, type Command, type CommandCtx } from '../command.js';\nimport { commandHelp, JSON_FLAG } from '../flag-help.js';\nimport { emitJson } from '../output.js';\nimport { CliError, UsageError } from '../errors.js';\nimport { splitAssignments } from '../bindings.js';\nimport type { MetricId } from '../../cli-api.js';\n\nconst FLAGS: FlagSpec[] = [\n  JSON_FLAG,\n  {\n    name: '--geodesic',\n    valueStyle: 'none',\n    description: 'Integrate a short Schwarzschild circular orbit, or a Kerr geodesic at the given θ.',\n  },\n];\n\nconst NAMES = ['minkowski', 'schwarzschild', 'flrw', 'kerr'] as const;\n\nconst HELP = `upt metric <minkowski|schwarzschild|flrw|kerr> [key=value ...] [--geodesic] [--json]\n        Christoffel symbols, the Ricci tensor, the Ricci scalar and the\n        Kretschmann scalar of one exact metric. Alias: upt curvature.\n        Signature of the line element is (-,+,+,+), the same mostly-plus\n        signature as the canonical Einstein-equation metric node.\n        Schwarzschild Kretschmann is checked against 48 G² M² / (c⁴ r⁶).\n        FLRW prints the Friedmann equation including the curvature term\n        −k c²/a² and Λ. Kerr Kretschmann uses the Boyer–Lindquist closed form.\n        --geodesic integrates a short Schwarzschild circular orbit, or a Kerr\n        geodesic (equatorial and circular at θ = π/2; inclined, with that\n        polar turning point, at any other θ). Kerr's a is a length in metres;\n        the report prints a/M = a/(GM/c²). A Ricci scalar that has a closed\n        form is that form; a vacuum finite-difference residual is not listed\n        as Ricci. A parameter is a number, a unit (M=1Msun, r=1km) or a\n        constant expression (theta=pi/2). A bare number is already in the\n        parameter's SI unit; a dimensioned quantity must match it.\n        \\`upt help metric\\` describes every flag.\n        e.g.  upt metric schwarzschild M=1Msun r=1e8 theta=pi/2\n              upt metric flrw k=1 t=2\n              upt curvature kerr a=1000 r=1e8 theta=1.2`;\n\nfunction isMetric(s: string | undefined): s is MetricId {\n  return (NAMES as readonly string[]).includes(s ?? '');\n}\n\n/**\n * The metric library refuses a bad value (a non-positive mass, a point inside\n * the horizon, an unknown or repeated parameter, a start the integrator cannot\n * take, a parameter in the wrong unit) with `MetricMassError`,\n * `MetricValueError` or a `UnitError`. Each is a bad value, exit 1. Any other\n * class (a plain `Error`, a TypeError, a RangeError) is a defect and is not\n * converted.\n */\nfunction metricValueError(api: CommandCtx['api'], e: unknown): unknown {\n  if (e instanceof api.MetricMassError || e instanceof api.MetricValueError || e instanceof api.UnitError) {\n    return new CliError(`upt metric: ${e.message}`);\n  }\n  return e;\n}\n\nasync function run(ctx: CommandCtx): Promise<number> {\n  const { args, api, out } = ctx;\n  const name = args.positionals[0];\n  if (!isMetric(name)) {\n    throw new UsageError(\n      `upt metric needs one of ${NAMES.join(', ')}. See \\`upt help metric\\`.`,\n    );\n  }\n  const pairs = args.positionals.slice(1);\n  // A token with no `=` is usage (exit 2); a parameter given twice is refused (exit 1).\n  splitAssignments('metric', pairs, 'key=value');\n  let report;\n  try {\n    report = api.curvatureReport(name, pairs);\n  } catch (e) {\n    throw metricValueError(api, e);\n  }\n  let geodesic:\n    | {\n        readonly kind: 'circular' | 'inclined';\n        readonly r0: number;\n        readonly rEnd: number;\n        readonly theta0?: number;\n        readonly thetaEnd?: number;\n        readonly phiAdvance: number;\n        readonly steps: number;\n        readonly E0?: number;\n        readonly EEnd?: number;\n        readonly L0?: number;\n        readonly LEnd?: number;\n        readonly Q0?: number;\n        readonly QEnd?: number;\n        readonly norm0?: number;\n        readonly normEnd?: number;\n        readonly note?: string;\n      }\n    | undefined;\n  if (args.flags.has('geodesic')) {\n    try {\n      if (name === 'schwarzschild') {\n        const M = report.parameters.M;\n        const r = report.parameters.r;\n        const orbit = api.schwarzschildCircularOrbit({\n          ...(M === undefined ? {} : { M }),\n          ...(r === undefined ? {} : { r }),\n        });\n        geodesic = { kind: 'circular', ...orbit };\n      } else if (name === 'kerr') {\n        const Mgeom = report.parameters.M_geom_m;\n        const a = report.parameters.a ?? 0;\n        const r = report.parameters.r;\n        const theta = report.parameters.theta ?? Math.PI / 2;\n        // Reachable: `curvatureReport` refuses M ≤ 0, but G=0 (or a c that overflows) makes\n        // M_geom = G M / c² zero with M > 0, and a supplied r passes its own check.\n        if (Mgeom === undefined || !(Mgeom > 0)) throw new CliError('upt metric: Kerr geodesic needs a positive mass');\n        const rOverM = (r ?? 10 * Mgeom) / Mgeom;\n        const aOverM = a / Mgeom;\n        if (!(Math.abs(aOverM) <= 1)) throw new CliError('upt metric: Kerr geodesic wants |a| ≤ GM/c²');\n        const equatorial = Math.abs(theta - Math.PI / 2) < 1e-6;\n        if (a === 0) {\n          // At a = 0 the Kerr metric is Schwarzschild. The Boyer–Lindquist integrator refuses the\n          // circular start there (\"radial potential is negative at the start\"): the two terms of\n          // its radial potential cancel to below its tolerance at a = 0, and only at a = 0. The\n          // equatorial circular orbit is the Schwarzschild one and is integrated as that, so the\n          // default parameters (a = 0) run. An inclined start has no Schwarzschild routine here.\n          if (!equatorial) {\n            throw new CliError(\n              'upt metric: at a = 0 an inclined Kerr orbit is a tilted Schwarzschild circular orbit, and this integrator needs a ≠ 0 for an inclined start; give a, or theta=pi/2',\n            );\n          }\n          const M = report.parameters.M;\n          geodesic = {\n            kind: 'circular' as const,\n            note: 'a = 0: the Kerr metric is Schwarzschild; the equatorial circular orbit is integrated by the Schwarzschild routine',\n            ...api.schwarzschildCircularOrbit({ ...(M === undefined ? {} : { M }), ...(r === undefined ? {} : { r }) }),\n          };\n        } else {\n        const shared = { M: Mgeom, aOverM, rOverM, fraction: 0.005, steps: 40 };\n        geodesic =\n          equatorial\n            ? { kind: 'circular' as const, ...api.kerrEquatorialCircular(shared) }\n            : {\n                kind: 'inclined' as const,\n                ...api.kerrGeodesic({\n                  ...shared,\n                  theta,\n                  ...api.kerrTurningPointOrbit({ M: Mgeom, aOverM, rOverM, theta }),\n                  mu2: 1,\n                }),\n              };\n        }\n      } else {\n        throw new UsageError('upt metric: --geodesic is for schwarzschild or kerr.');\n      }\n    } catch (e) {\n      if (e instanceof UsageError || e instanceof CliError) throw e;\n      throw metricValueError(api, e);\n    }\n  }\n  if (args.flags.has('json')) {\n    emitJson({ command: 'metric', result: { ...report, ...(geodesic === undefined ? {} : { geodesic }) } }, ctx.write);\n    return 0;\n  }\n  out(`\\n${report.metric}  signature ${report.signature}`);\n  out(`  ${report.signatureNote}`);\n  out(`  coordinates ${report.coordinates}`);\n  out('  parameters: ' + Object.entries(report.parameters).map(([k, v]) => `${k}=${v}`).join(' '));\n  out('  point: ' + Object.entries(report.point).map(([k, v]) => `${k}=${v}`).join(' '));\n  for (const note of report.notes) out(`  ${note}`);\n  out('  Christoffel (nonzero):');\n  if (report.christoffel.length === 0) out('    (none)');\n  for (const c of report.christoffel) out(`    ${c.index} = ${c.value}`);\n  out('  Ricci (nonzero):');\n  if (report.ricci.length === 0) out('    (none)');\n  for (const c of report.ricci) out(`    ${c.index} = ${c.value}`);\n  out(`  Ricci scalar R = ${report.ricciScalar}`);\n  out(`  Kretschmann K = ${report.kretschmann}`);\n  out('  closed form: ' + Object.entries(report.closedForm).map(([k, v]) => `${k}=${v}`).join(' '));\n  if (geodesic) {\n    const cons =\n      geodesic.E0 === undefined\n        ? ''\n        : ` E ${geodesic.E0}→${geodesic.EEnd} L ${geodesic.L0}→${geodesic.LEnd} Q ${geodesic.Q0}→${geodesic.QEnd} norm ${geodesic.norm0}→${geodesic.normEnd}`;\n    const polar =\n      geodesic.theta0 === undefined ? '' : ` θ ${geodesic.theta0}→${geodesic.thetaEnd}`;\n    const label = geodesic.kind === 'inclined' ? 'inclined orbit' : 'circular orbit';\n    out(\n      `  geodesic: ${label} r0=${geodesic.r0} rEnd=${geodesic.rEnd}${polar} Δφ=${geodesic.phiAdvance} over ${geodesic.steps} steps${cons}`,\n    );\n    if (geodesic.note !== undefined) out(`    ${geodesic.note}`);\n  }\n  return 0;\n}\n\nexport const command: Command = {\n  name: 'metric',\n  aliases: ['curvature'],\n  positionals: { max: Infinity },\n  flags: FLAGS,\n  help: commandHelp(HELP, FLAGS),\n  summary: 'Print Christoffel symbols and curvature scalars for one exact metric.',\n  example: 'upt metric schwarzschild M=1Msun r=1e8 theta=pi/2',\n  group: 'evaluate',\n  run,\n};\n\nregisterCommand(command);\n"]}

package/dist/cli/commands/path.d.ts

@@ -21,6 +21,16 @@ type AtlasBridge = import('../../cli-api.js').AtlasBridge;
 type Bridges = readonly AtlasBridge[];
 type Translation = import('../../cli-api.js').ObservableTranslation;
 type Carriage = import('../../cli-api.js').ObservableCarriage;
+/**
+ * A translation declared on a bridge that states no bound. A translation restates the bound's
+ * error in another observable, so it presupposes the bound; a bridge without one declaring a
+ * translation is a defect of the atlas record, not a user error, and this is not a `CliError`.
+ */
+export declare class TranslationWithoutBoundError extends Error {
+    constructor(bridgeId: string, observable: string);
+}
+/** The first bridge's bound, which the translation `tr` presupposes. @internal */
+export declare function translatedBound(bridges: Bridges, tr: Translation): NonNullable<AtlasBridge['bound']>;
 /**
  * Whether the path is adequate for a requested tolerance at one point: every
  * regime holds, every horizon holds at the given t, and the closed-form point

package/dist/cli/commands/path.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"path.d.ts","sourceRoot":"","sources":["../../../src/cli/commands/path.ts"],"names":[],"mappings":"AA2BA,OAAO,EAAmB,KAAK,OAAO,EAAE,KAAK,UAAU,EAAE,MAAM,eAAe,CAAC;AA0I/E,UAAU,aAAa;IACrB,QAAQ,EAAE,MAAM,CAAC;IACjB,OAAO,EAAE,MAAM,CAAC;IAChB,KAAK,EAAE,OAAO,GAAG,IAAI,CAAC;IACtB,gGAAgG;IAChG,UAAU,CAAC,EAAE;QAAE,SAAS,EAAE,MAAM,CAAC;QAAC,OAAO,EAAE,MAAM,CAAA;KAAE,EAAE,CAAC;IACtD,8FAA8F;IAC9F,WAAW,CAAC,EAAE,cAAc,CAAC;CAC9B;AAID,qFAAqF;AACrF,wBAAgB,UAAU,CACxB,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,IAAI,EAAE,MAAM,GACX;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,EAAE,CAAC;IAAC,OAAO,EAAE,QAAQ,GAAG,KAAK,CAAA;CAAE,CA8B/D;AAID,KAAK,WAAW,GAAG,OAAO,kBAAkB,EAAE,WAAW,CAAC;AAE1D,KAAK,OAAO,GAAG,SAAS,WAAW,EAAE,CAAC;AACtC,KAAK,WAAW,GAAG,OAAO,kBAAkB,EAAE,qBAAqB,CAAC;AACpE,KAAK,QAAQ,GAAG,OAAO,kBAAkB,EAAE,kBAAkB,CAAC;AAiH9D;;;;;;GAMG;AACH,wBAAgB,cAAc,CAC5B,SAAS,EAAE,MAAM,EACjB,CAAC,EAAE,UAAU,EACb,MAAM,EAAE,OAAO,GACd;IAAE,OAAO,EAAE,UAAU,GAAG,YAAY,GAAG,cAAc,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAazE;AAED,iFAAiF;AACjF,MAAM,WAAW,gBAAgB;IAC/B,UAAU,EAAE,MAAM,GAAG,IAAI,CAAC;IAC1B,KAAK,EAAE,MAAM,CAAC;CACf;AAED,0KAA0K;AAC1K,wBAAgB,cAAc,CAAC,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,MAAM,GAAG,SAAS,GAAG,gBAAgB,GAAG,IAAI,CAgBvG;AAKD,yFAAyF;AACzF,UAAU,mBAAmB;IAC3B,OAAO,EAAE,UAAU,GAAG,YAAY,GAAG,cAAc,CAAC;IACpD,MAAM,EAAE,MAAM,CAAC;IACf,WAAW,EAAE,WAAW,GAAG,IAAI,CAAC;IAChC,8EAA8E;IAC9E,SAAS,EAAE,QAAQ,EAAE,CAAC;IACtB,2EAA2E;IAC3E,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,0EAA0E;IAC1E,cAAc,EAAE,MAAM,GAAG,IAAI,CAAC;IAC9B,2GAA2G;IAC3G,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,2EAA2E;IAC3E,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,0EAA0E;IAC1E,aAAa,EAAE,MAAM,GAAG,IAAI,CAAC;IAC7B,6EAA6E;IAC7E,kBAAkB,EAAE,OAAO,GAAG,IAAI,CAAC;CACpC;AAED;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAC7B,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,OAAO,EAAE,gBAAgB,GAAG;IAAE,UAAU,EAAE,MAAM,CAAA;CAAE,EAClD,OAAO,EAAE,OAAO,EAChB,CAAC,EAAE,UAAU,EACb,EAAE,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GACnC,mBAAmB,CAkFrB;AAkCD,KAAK,YAAY,GACb;IAAE,WAAW,EAAE,MAAM,CAAA;CAAE,GACvB;IACE,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IACf,MAAM,EAAE,MAAM,CAAC;IACf,OAAO,CAAC,EAAE;QAAE,EAAE,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,aAAa,EAAE,OAAO,CAAA;KAAE,CAAC;CACjG,CAAC;AAqBN,kHAAkH;AAClH,wBAAgB,YAAY,CAC1B,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,OAAO,EAAE,gBAAgB,GAAG;IAAE,UAAU,EAAE,MAAM,CAAA;CAAE,EAClD,OAAO,EAAE,OAAO,EAChB,CAAC,EAAE,UAAU,EACb,EAAE,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACpC,WAAW,EAAE,OAAO,GACnB,mBAAmB,GAAG;IAAE,YAAY,EAAE,YAAY,GAAG,IAAI,CAAA;CAAE,CAa7D;AAwHD,KAAK,UAAU,GAAG;IAChB,cAAc,EAAE,OAAO,GAAG,SAAS,CAAC;IACpC,QAAQ,EAAE,SAAS,aAAa,EAAE,CAAC;IACnC,OAAO,EAAE,OAAO,CAAC;IACjB,UAAU,EAAE;QAAE,CAAC,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAE,GAAG,IAAI,CAAC;IAChD,gBAAgB,EAAE,MAAM,GAAG,IAAI,CAAC;IAChC,gBAAgB,EAAE,SAAS;QAAE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;QAAC,QAAQ,CAAC,QAAQ,EAAE,SAAS,MAAM,EAAE,CAAA;KAAE,EAAE,CAAC;CAC/F,CAAC;AA4qBF,eAAO,MAAM,OAAO,EAAE,OASrB,CAAC"}
+{"version":3,"file":"path.d.ts","sourceRoot":"","sources":["../../../src/cli/commands/path.ts"],"names":[],"mappings":"AA2BA,OAAO,EAAmB,KAAK,OAAO,EAAE,KAAK,UAAU,EAAE,MAAM,eAAe,CAAC;AA2I/E,UAAU,aAAa;IACrB,QAAQ,EAAE,MAAM,CAAC;IACjB,OAAO,EAAE,MAAM,CAAC;IAChB,KAAK,EAAE,OAAO,GAAG,IAAI,CAAC;IACtB,gGAAgG;IAChG,UAAU,CAAC,EAAE;QAAE,SAAS,EAAE,MAAM,CAAC;QAAC,OAAO,EAAE,MAAM,CAAA;KAAE,EAAE,CAAC;IACtD,8FAA8F;IAC9F,WAAW,CAAC,EAAE,cAAc,CAAC;CAC9B;AAID,qFAAqF;AACrF,wBAAgB,UAAU,CACxB,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,IAAI,EAAE,MAAM,GACX;IAAE,IAAI,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,EAAE,CAAC;IAAC,OAAO,EAAE,QAAQ,GAAG,KAAK,CAAA;CAAE,CA+B/D;AAID,KAAK,WAAW,GAAG,OAAO,kBAAkB,EAAE,WAAW,CAAC;AAE1D,KAAK,OAAO,GAAG,SAAS,WAAW,EAAE,CAAC;AACtC,KAAK,WAAW,GAAG,OAAO,kBAAkB,EAAE,qBAAqB,CAAC;AACpE,KAAK,QAAQ,GAAG,OAAO,kBAAkB,EAAE,kBAAkB,CAAC;AAE9D;;;;GAIG;AACH,qBAAa,4BAA6B,SAAQ,KAAK;IACrD,YAAY,QAAQ,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,EAM/C;CACF;AAED,kFAAkF;AAClF,wBAAgB,eAAe,CAAC,OAAO,EAAE,OAAO,EAAE,EAAE,EAAE,WAAW,GAAG,WAAW,CAAC,WAAW,CAAC,OAAO,CAAC,CAAC,CAIpG;AAiHD;;;;;;GAMG;AACH,wBAAgB,cAAc,CAC5B,SAAS,EAAE,MAAM,EACjB,CAAC,EAAE,UAAU,EACb,MAAM,EAAE,OAAO,GACd;IAAE,OAAO,EAAE,UAAU,GAAG,YAAY,GAAG,cAAc,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAazE;AAED,iFAAiF;AACjF,MAAM,WAAW,gBAAgB;IAC/B,UAAU,EAAE,MAAM,GAAG,IAAI,CAAC;IAC1B,KAAK,EAAE,MAAM,CAAC;CACf;AAED,0KAA0K;AAC1K,wBAAgB,cAAc,CAAC,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EAAE,GAAG,EAAE,MAAM,GAAG,SAAS,GAAG,gBAAgB,GAAG,IAAI,CAgBvG;AAKD,yFAAyF;AACzF,UAAU,mBAAmB;IAC3B,OAAO,EAAE,UAAU,GAAG,YAAY,GAAG,cAAc,CAAC;IACpD,MAAM,EAAE,MAAM,CAAC;IACf,WAAW,EAAE,WAAW,GAAG,IAAI,CAAC;IAChC,8EAA8E;IAC9E,SAAS,EAAE,QAAQ,EAAE,CAAC;IACtB,2EAA2E;IAC3E,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,0EAA0E;IAC1E,cAAc,EAAE,MAAM,GAAG,IAAI,CAAC;IAC9B,2GAA2G;IAC3G,iBAAiB,EAAE,MAAM,GAAG,IAAI,CAAC;IACjC,2EAA2E;IAC3E,OAAO,EAAE,MAAM,GAAG,IAAI,CAAC;IACvB,0EAA0E;IAC1E,aAAa,EAAE,MAAM,GAAG,IAAI,CAAC;IAC7B,6EAA6E;IAC7E,kBAAkB,EAAE,OAAO,GAAG,IAAI,CAAC;CACpC;AAED;;;;;;;;;GASG;AACH,wBAAgB,eAAe,CAC7B,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,OAAO,EAAE,gBAAgB,GAAG;IAAE,UAAU,EAAE,MAAM,CAAA;CAAE,EAClD,OAAO,EAAE,OAAO,EAChB,CAAC,EAAE,UAAU,EACb,EAAE,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,GACnC,mBAAmB,CAkFrB;AAkCD,KAAK,YAAY,GACb;IAAE,WAAW,EAAE,MAAM,CAAA;CAAE,GACvB;IACE,EAAE,EAAE,MAAM,CAAC;IACX,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IACf,MAAM,EAAE,MAAM,CAAC;IACf,OAAO,CAAC,EAAE;QAAE,EAAE,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,MAAM,EAAE,MAAM,CAAC;QAAC,aAAa,EAAE,OAAO,CAAA;KAAE,CAAC;CACjG,CAAC;AAqBN,kHAAkH;AAClH,wBAAgB,YAAY,CAC1B,GAAG,EAAE,UAAU,CAAC,KAAK,CAAC,EACtB,OAAO,EAAE,gBAAgB,GAAG;IAAE,UAAU,EAAE,MAAM,CAAA;CAAE,EAClD,OAAO,EAAE,OAAO,EAChB,CAAC,EAAE,UAAU,EACb,EAAE,EAAE,QAAQ,CAAC,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC,EACpC,WAAW,EAAE,OAAO,GACnB,mBAAmB,GAAG;IAAE,YAAY,EAAE,YAAY,GAAG,IAAI,CAAA;CAAE,CAa7D;AAwHD,KAAK,UAAU,GAAG;IAChB,cAAc,EAAE,OAAO,GAAG,SAAS,CAAC;IACpC,QAAQ,EAAE,SAAS,aAAa,EAAE,CAAC;IACnC,OAAO,EAAE,OAAO,CAAC;IACjB,UAAU,EAAE;QAAE,CAAC,EAAE,MAAM,CAAC;QAAC,KAAK,EAAE,MAAM,CAAA;KAAE,GAAG,IAAI,CAAC;IAChD,gBAAgB,EAAE,MAAM,GAAG,IAAI,CAAC;IAChC,gBAAgB,EAAE,SAAS;QAAE,QAAQ,CAAC,KAAK,EAAE,MAAM,CAAC;QAAC,QAAQ,CAAC,QAAQ,EAAE,SAAS,MAAM,EAAE,CAAA;KAAE,EAAE,CAAC;CAC/F,CAAC;AA4qBF,eAAO,MAAM,OAAO,EAAE,OAUrB,CAAC"}

package/dist/cli/commands/path.js

@@ -2,6 +2,7 @@ import { registerCommand } from '../command.js';
 import { commandHelp, JSON_FLAG } from '../flag-help.js';
 import { CliError, EXIT_CHECK_FAILED, UsageError } from '../errors.js';
 import { emitJson } from '../output.js';
+import { integerFlag } from '../numeric-flag.js';
 import { parseAt, resolveAtPoint, showInequality } from './regime.js';
 import { explainsRefusal, missingForComposite, routeClaim, selectRoute, transportReport } from './_atlas-route.js';
 const FLAGS = [
@@ -136,7 +137,8 @@ export function parseSweep(api, spec) {
     };
     const lo = endpoint(loS);
     const hi = endpoint(hiS);
-    const n = Number(nS);
+    // The grammar first (`0x5` is usage, exit 2); the 2..MAX_SAMPLES range after (exit 1).
+    const n = integerFlag('path', '--sweep sample count', nS);
     const spacing = sp === 'log' ? 'log' : 'linear';
     if (!Number.isFinite(lo) || !Number.isFinite(hi) || !(lo < hi)) {
         throw new CliError(`upt path: --sweep '${spec}' needs finite lo < hi`);
@@ -149,6 +151,25 @@ export function parseSweep(api, spec) {
     const values = Array.from({ length: n }, (_, i) => spacing === 'log' ? lo * (hi / lo) ** (i / (n - 1)) : lo + ((hi - lo) * i) / (n - 1));
     return { name, values, spacing };
 }
+/**
+ * A translation declared on a bridge that states no bound. A translation restates the bound's
+ * error in another observable, so it presupposes the bound; a bridge without one declaring a
+ * translation is a defect of the atlas record, not a user error, and this is not a `CliError`.
+ */
+export class TranslationWithoutBoundError extends Error {
+    constructor(bridgeId, observable) {
+        super(`upt path: '${bridgeId}' declares a translation to ${observable} but states no bound; ` +
+            'a translation presupposes the bound it translates (a defect of the record, not a usage error)');
+        this.name = 'TranslationWithoutBoundError';
+    }
+}
+/** The first bridge's bound, which the translation `tr` presupposes. @internal */
+export function translatedBound(bridges, tr) {
+    const first = bridges[0];
+    if (first.bound === undefined)
+        throw new TranslationWithoutBoundError(first.id, tr.observable);
+    return first.bound;
+}
 function modelFamilyOf(api, id) {
     if (id === null)
         return undefined;
@@ -448,7 +469,7 @@ export function judgeAtPoint(api, request, bridges, e, at, withControl) {
 }
 function thresholdOrigin(request, tr, bridges) {
     return (`t* is derived from the requested tolerance ${request.value} ${tr.unit} through the translation; ` +
-        `${bridges[0].id}'s own horizon (${bridges[0].bound.horizon}) is the record's fixed threshold, judged separately`);
+        `${bridges[0].id}'s own horizon (${translatedBound(bridges, tr).horizon}) is the record's fixed threshold, judged separately`);
 }
 function observableReport(request, j, bridges, evidence) {
     const tr = j.translation;
@@ -479,7 +500,7 @@ function observableReport(request, j, bridges, evidence) {
                 domainHorizon: j.domainHorizon,
                 horizonConvention: tr.boundary,
                 thresholdOrigin: thresholdOrigin(request, tr, bridges),
-                bridgeHorizon: { horizon: bridges[0].bound.horizon, holds: j.bridgeHorizonHolds },
+                bridgeHorizon: { horizon: translatedBound(bridges, tr).horizon, holds: j.bridgeHorizonHolds },
                 evidence,
                 pointWitness: j.pointWitness,
             }),
@@ -1161,6 +1182,7 @@ async function run(ctx) {
 export const command = {
     name: 'path',
     aliases: [],
+    positionals: { max: Infinity },
     flags: FLAGS,
     help: commandHelp(HELP, FLAGS),
     summary: 'Show the bridge chain between two models and whether a bound is claimed there.',