@essrt/physics 9.1.0 → 9.1.1

package/README.md

@@ -1,7 +1,7 @@
 # @essrt/physics
 
 The physics core of **ESSRT** — the Expanding Solar System Resonance Theory —
-as implemented by the [Holistic Universe Model](https://3d.holisticuniverse.com)
+as implemented by the [ESSRT simulator](https://3d.holisticuniverse.com)
 ([source, AGPL-3.0](https://github.com/dvansonsbeek/3d)). Pure computation:
 no I/O, no globals, no DOM, constants injected. The package ships the
 **complete model**, including the fitted coefficients, so an installed copy

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@essrt/physics",
-  "version": "9.1.0",
+  "version": "9.1.1",
   "description": "ESSRT (Expanding Solar System Resonance Theory) physics core — 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": "v14.0",

package/src/deltat/deep-time.cjs

@@ -125,7 +125,15 @@ function createDeepTimeLod(deps) {
    * the unit H(t) below.
    * @param {number} t_Ma @returns {number|null} */
   function eraClockHAtAge(t_Ma) {
-    const LOD_s = lodSecondsAtAge(t_Ma);
+    return eraClockHFromLod(lodSecondsAtAge(t_Ma));
+  }
+  // The `…FromLod` twins take an ALREADY-COMPUTED LOD(t): the Simpson ΔT
+  // integrand below evaluated the LOD four times and the tropical year twice
+  // per node through these helpers (measured 8.8 → 4.8 ms per cardinal panel
+  // year with the model's ΔT memo). Same operations, same inputs — every
+  // value bit-identical (11,958-value sweep, 0 differences).
+  /** @param {number|null} LOD_s @returns {number|null} */
+  function eraClockHFromLod(LOD_s) {
     if (LOD_s === null) return null;
     return K.holisticYearJ2000 * LOD_s / K.lodNowH13Seconds;
   }
@@ -156,7 +164,11 @@ function createDeepTimeLod(deps) {
    * ≡ H_era wherever (a₀/a_M)³ ≈ 1. Identifier kept (plan 06 P4).
    * @param {number} t_Ma @returns {number|null} */
   function hAtAge(t_Ma) {
-    const hEra = eraClockHAtAge(t_Ma);
+    return hFromLod(t_Ma, lodSecondsAtAge(t_Ma));
+  }
+  /** @param {number} t_Ma @param {number|null} LOD_s @returns {number|null} */
+  function hFromLod(t_Ma, LOD_s) {
+    const hEra = eraClockHFromLod(LOD_s);
     const term = composed.torqueTermAtAge(t_Ma);
     return hEra === null || term === null ? null : hEra / term;
   }
@@ -202,8 +214,12 @@ function createDeepTimeLod(deps) {
    * ±2 Myr, measured), and re-basing would move every comb coefficient
    * for a rounding-level gain. This tier stays a named device. */
   function tropicalYearSecondsAtAge(t_Ma) {
+    return tropicalYearSecondsFromLod(t_Ma, lodSecondsAtAge(t_Ma));
+  }
+  /** @param {number} t_Ma @param {number|null} LOD_s @returns {number} */
+  function tropicalYearSecondsFromLod(t_Ma, LOD_s) {
     const sidSec = siderealYearSecondsAtAge(t_Ma);
-    const Ht = hAtAge(t_Ma);
+    const Ht = hFromLod(t_Ma, LOD_s);
     if (Ht === null) return sidSec * (1 - 13 / K.holisticYearJ2000);
     return sidSec * (1 - 13 / Ht);
   }
@@ -281,15 +297,16 @@ function createDeepTimeLod(deps) {
       const tau = i * h;
       const lodMean = lodSecondsAtAge(tau);
       if (lodMean === null) return NaN;
-      const yearS = tropicalYearSecondsAtAge(tau);
+      const yearS = tropicalYearSecondsFromLod(tau, lodMean);
       // Ecliptic "missing motion" — the solar day is measured against the
       // Sun on the ECLIPTIC, whose plane turns on the invariable plane at the
       // nodal rate s₃ (period ≈ 68,751 yr, an orbital quantity), not in the
       // inclination frame. Adds ~3.4 ms at J2000; the fitted cycle stack
       // closes Layer-4 LOD_real onto the USNO anchor. Non-null: the helper
       // can only be null when lodSecondsAtAge(tau) is null, which already
-      // returned NaN above — the checker can't see the chain.
-      const lodH5Raw = lodMean + /** @type {number} */ (eclipticLodCorrectionSecondsAtAge(tau));
+      // returned NaN above. Spelled inline from the node's own lodMean and
+      // yearS — the same operations as eclipticLodCorrectionSecondsAtAge(tau).
+      const lodH5Raw = lodMean + lodMean / (deps.nodalPeriodYearsFn() * (yearS / 86400));
       const integrand = (86400 - lodH5Raw) * yearS * 1e6 / 86400;
       const w = (i === 0 || i === n) ? 1 : (i % 2 === 1 ? 4 : 2);
       sum += w * integrand;

package/src/model.js

@@ -411,9 +411,9 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
   // own ±10-Myr mode table, anchored form — feeds the ENTIRE lunar chain
   // (modulation, cycle counts, arguments eccAt/channelIntegral, E-factor).
   // The Sun/clock machinery (eclipse Sun EoC, besselian Sun distance,
-  // cardinal braid) stays on the H/3 channel above — a certification
-  // split, not a physics one (the two agree within 4.2e-5 in-era). The
-  // laws hook deliberately does NOT reach this channel.
+  // cardinal braid) reads the one-source series below (plan 06 layer B —
+  // the two agree within 4.2e-5 in-era). The laws hook deliberately does
+  // NOT reach this channel.
   const deepEcc = createDeepEccChannel(DEEP_MODES_ARTIFACT);
   /** @param {number} year @returns {number} */
   // Plan 06 layer B: the PUBLISHED Earth eccentricity is the one-source series
@@ -673,8 +673,17 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
 
   // ── ΔT (TT − UT1): raw Simpson + sequential stack adds ────────────────────
   /** @param {number} tMa @returns {number} */
+  // Exact-argument memo (the alphaMemo pattern above): the cardinal Newton
+  // solve asks for ΔT at the SAME jd twice per step (jdTTFromUT in the
+  // apparent Sun and again under the planetary completion) — measured 8.8 →
+  // 5.6 ms per panel year (8 solves) at modern epochs, values bit-identical.
+  /** @type {Map<number, number>} */
+  const dtMemo = new Map();
+  /** @param {number} tMa @returns {number} */
   const meanDeltaTSecondsAtAge = (tMa) => {
     if (tMa === 0) return 0;
+    const hit = dtMemo.get(tMa);
+    if (hit !== undefined) return hit;
     let result = deepLod.deltaTRawSecondsAtAge(tMa);
     const yearY = 2000 - tMa * 1e6;
     result += dtCycles.cycleDeltaTSecondsAt('bond', yearY);
@@ -682,6 +691,8 @@ export function assembleModel(C, F, laws = {}, secularSeriesArtifact = /** @type
     result += dtCycles.cycleDeltaTSecondsAt('jose5', yearY);
     result += dtCycles.cycleDeltaTSecondsAt('jose4', yearY);
     result += dtCycles.swingDeltaTSecondsAt(yearY);
+    if (dtMemo.size >= 8192) dtMemo.clear();
+    dtMemo.set(tMa, result);
     return result;
   };
   /** @param {number} year @returns {number} */

package/src/planets/orientation.cjs

@@ -14,7 +14,7 @@
  *
  * Two ascending-node conventions exist and BOTH are quantities (the S-P4
  * lesson): the year-2000-anchored linear node here
- * (ascendingNodeInvPlaneLinearAt — the node-integrator/dashboard mirror)
+ * (ascendingNodeInvPlaneLinearAt — the node-integrator mirror)
  * and the balanced-year −8H/N node inside the scene tilt
  * (eclipticInclinationFromBalanced). Do not merge them.
  */
@@ -54,7 +54,7 @@ function invPlaneInclinationAt(body, yearsSinceBalanced, env) {
 
 /**
  * Ascending node on the invariable plane — the LINEAR year-2000-anchored
- * convention (node-integrator/dashboard; see the header).
+ * convention (node-integrator; see the header).
  * @param {{ ascendingNodeInvPlane?: number, ascendingNodePeriod?: number,
  *   perihelionEclipticYears: number }} body
  * @param {number} year @returns {number} degrees 0–360 */

package/types/planets/orientation.d.cts

@@ -25,7 +25,7 @@ export function invPlaneInclinationAt(body: {
 }): number;
 /**
  * Ascending node on the invariable plane — the LINEAR year-2000-anchored
- * convention (node-integrator/dashboard; see the header).
+ * convention (node-integrator; see the header).
  * @param {{ ascendingNodeInvPlane?: number, ascendingNodePeriod?: number,
  *   perihelionEclipticYears: number }} body
  * @param {number} year @returns {number} degrees 0–360 */