@essrt/physics 4.7.0 → 4.8.0

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "4.7.0",
+  "version": "4.8.0",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core \u2014 pure, injectable-constants, no I/O, no globals, no DOM. The complete model incl. fitted coefficients; every published version immutably ships one recorded model identity.",
   "essrt": {
     "modelVersion": "v13.0",

package/src/chain-cycles/index.cjs

@@ -73,6 +73,29 @@ function createChainCycleIntegrator({
   /** @type {Map<PeriodSecondsAtAge, Float64Array | null>} */
   const tables = new Map();
 
+  // Exact-argument memo for the year-length integrand of the out-of-range
+  // Simpson fallback. Within one animation frame every integrator node
+  // (the six-plus Moon chains + the seven planets) integrates the SAME
+  // (anchor, currentYear) pair, so the Simpson sample lattice — and with
+  // it every tropicalYearSecondsAtAge(tMa) — repeats across chains, up to
+  // 13× per sample. Exact-key reuse returns the bit-identical value: pure
+  // dedupe, not a resolution change (the sampler-purity trap class is
+  // about grids serving NEARBY arguments; this serves only the same one).
+  // Bounded by wholesale clear — one frame's lattice is ≤ ~1025 entries.
+  /** @type {Map<number, number | null>} */
+  const yearSecondsMemo = new Map();
+  const YEAR_SECONDS_MEMO_CAP = 4096;
+  /** @param {number} tMa @returns {number | null} */
+  function tropicalYearSecondsMemo(tMa) {
+    let v = yearSecondsMemo.get(tMa);
+    if (v === undefined) {
+      v = tropicalYearSecondsAtAge(tMa);
+      if (yearSecondsMemo.size >= YEAR_SECONDS_MEMO_CAP) yearSecondsMemo.clear();
+      yearSecondsMemo.set(tMa, v);
+    }
+    return v;
+  }
+
   /** Lazy cumulative table for one chain; null when the chain is undefined
    *  anywhere in range (→ Simpson fallback).
    *  @param {PeriodSecondsAtAge} periodFnSeconds
@@ -167,7 +190,7 @@ function createChainCycleIntegrator({
       const tMa = (ageAnchorYear - y) / 1e6;
       const tPeriodS = periodFnSeconds(tMa);
       if (tPeriodS === null) return null;
-      const tYrS = tropicalYearSecondsAtAge(tMa);
+      const tYrS = tropicalYearSecondsMemo(tMa);
       if (tYrS === null) return null;
       const integrand = tYrS / tPeriodS;  // cycles per SI year
       const w = (i === 0 || i === n) ? 1 : (i % 2 === 1 ? 4 : 2);

package/src/earth/year-lengths.cjs

@@ -97,7 +97,7 @@ function createYearLengths({ sampleAt, massLossSiderealSecondsAtYearFn, apsidalS
     axialPrecessionYearsAtYear: (y) => { const s = sid(y), t = trop(y); return s / (s - t); },
     /** Perihelion (apsidal-vs-equinox) precession period of date, years. @param {number} y */
     perihelionPrecessionYearsAtYear: (y) => { const a = anom(y), t = trop(y); return a / (a - t); },
-    /** Inclination precession period of date, years (hypersensitive: ±1 s ↔ ≈400 yr). @param {number} y */
+    /** Apsidal precession period of date vs the fixed stars, years (hypersensitive: ±1 s ↔ ≈400 yr; the identifier keeps the historical "inclination" name). @param {number} y */
     inclinationPrecessionYearsAtYear: (y) => { const a = anom(y), s = sid(y); return a / (a - s); },
     /** The per-cardinal layer (D6 rules: year lengths λ̇-corrected; offsets/spreads raw — μs-class/second-order there). */
     cardinal: Object.freeze({

package/src/model.js

@@ -1239,6 +1239,20 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     }),
     eclipse: Object.freeze({
       sunLonDegAtJD: /** @param {number} jd @returns {number} */ (jd) => eclipseFinders.sunLonDegAt(jd),
+      // The COMPLETED tier Sun (K8b follow-up, owner-approved injection):
+      // the bare finder Sun minus the DERIVED planetary-completion table —
+      // 70 framework-carrier terms plus the 6.44″ Earth-around-EMB wobble
+      // (the "lunar equation"; eclipse/sun-planetary-completion.cjs, zero
+      // fitted constants). The finders themselves deliberately stay bare
+      // (their certified canon statistics were produced on the bare form;
+      // besselian.cjs applies the same subtraction internally). This
+      // surface exists for the SCENE tier's E5 wheel-Sun bridge, so the
+      // rendered Sun carries the model's own derived apparent-class terms.
+      // Clock: the completion argument rides TT centuries from the same
+      // framework ΔT the besselian composition uses.
+      sunLonCompletedDegAtJD: /** @param {number} jd @returns {number} */ (jd) =>
+        eclipseFinders.sunLonDegAt(jd)
+        - sunPlanetaryCompletionDeg((jdTTFromUT(jd) - j2000JD) / julianCenturyDays),
       findLunarInRange: /** @param {number} jdStart @param {number} jdEnd */ (jdStart, jdEnd) => eclipseFinders.findLunarEclipsesInRange(jdStart, jdEnd),
       findSolarInRange: /** @param {number} jdStart @param {number} jdEnd */ (jdStart, jdEnd) => eclipseFinders.findSolarEclipsesInRange(jdStart, jdEnd),
       deltaTSecondsAtJD: frameworkDeltaTSecondsAtJD,