@essrt/physics 20.0.1 → 21.0.0

package/src/layer0/index.js

@@ -43,11 +43,11 @@
  * @property {number} moonEccentricityFactor     sqrt(1 - e^2)
  * @property {number} earthMassKg
  * @property {number} earthRadiusM
- * @property {number} holisticYearJ2000
  * @property {number} lodNowH13Seconds
  * @property {number} siderealYearJ2000Seconds
  * @property {number} solarMassLossFracPerYear
  * @property {number} precessionSolarShareJ2000  f_S, the solar fraction of the J2000 precession torque
+ * @property {number} precessionRateJ2000ArcsecPerYr  p₀, the composed precession rate's J2000 anchor (1,296,000/axial0)
  */
 
 /**
@@ -55,12 +55,8 @@
  * @property {(year: number) => number} tMa
  * @property {(year: number) => number} moonDistanceMetres
  * @property {(year: number) => number|null} lodSeconds
- * @property {(year: number) => number|null} holisticH               the UNIT: 13 × the composed lunisolar precession period
  * @property {(year: number) => number} siderealYearSeconds
- * @property {(year: number) => number} tropicalYearSeconds
- * @property {(year: number) => number|null} anomalisticYearSeconds
- * @property {(year: number) => number|null} eraClockH               the frozen era clock's phase convention H₀·LOD/LOD₀ (device tier)
- * @property {(year: number) => number} eraClockTropicalYearSeconds  T_sid·(1 − 13/H_era) — the frozen comb family's base
+ * @property {(year: number) => number} tropicalYearSeconds      sidereal less one equinox turn per T_p(t) (the composed precession period)
  */
 
 /**
@@ -100,35 +96,10 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
          - p.moonMassKg * Math.sqrt(p.gmEarthMoonM3S2 * a) * p.moonEccentricityFactor);
   };
 
-  /**
-   * The frozen era clock's phase convention: H_era(t) = H₀·LOD(t)/LOD₀ —
-   * pure spin scaling. This is the counter the FROZEN devices (the cardinal
-   * era clock, the year-length comb family) were fitted against; it ships
-   * with their coefficients as a named device constant (plan 06 Phase 3,
-   * D8), never as "H(t)".
-   * @param {number} t age in Ma
-   * @returns {number|null} years
-   */
-  const eraClockHCore = (t) => {
-    const lod = lodSecondsCore(t);
-    return lod === null ? null : p.holisticYearJ2000 * lod / p.lodNowH13Seconds;
-  };
-
-  /**
-   * The unit H(t) ≡ 13 × the composed lunisolar precession period:
-   * H_era(t) / [f_S + (1 − f_S)·(a₀/a_M(t))³] — the spin scaling AND the
-   * lunar torque growing on the recession history (plan 06 D6/Phase 3; the
-   * formula's home is earth/precession-composed, this is its Layer-0 twin,
-   * held bit-identical to deltat/deep-time.cjs by the layer0 gate).
-   * @param {number} t age in Ma
-   * @returns {number|null} years
-   */
-  const holisticHCore = (t) => {
-    const hEra = eraClockHCore(t);
-    if (hEra === null) return null;
-    const lf = Math.pow(p.moonDistanceNowM / moonDistanceMetresCore(t), 3);
-    return hEra / (p.precessionSolarShareJ2000 + (1 - p.precessionSolarShareJ2000) * lf);
-  };
+  // (Plan 09 step 7b: the unit H(t) = H₀·(LOD/LOD₀)/[f_S + (1 − f_S)(a₀/a_M)³]
+  // — `holisticH`, 13.011 composed precession periods, a fit constant's
+  // multiple of T_p(t) — and its spin-scaling factor (step 7a's `eraClockH`)
+  // are deleted. The physics is the composed precession rate below.)
 
   /**
    * Sidereal year. Driver 2 alone — solar mass loss via Kepler, `dT/T = -2 dM/M`
@@ -149,57 +120,40 @@ export const createEpochPrimitives = ({ params: p, alphaAtAgeMa }) => {
   };
 
   /**
-   * Tropical year — sidereal less one calendar turn of the unit, 13/H(t): the
-   * unit's CALENDAR convention (plan 06 S5 — a device identity, 0.086 % from
-   * one turn per the published precession period; D2 decided: it stays a
-   * named device, never re-based).
-   * Falls back to the J2000 H past the asymptote, matching the shipped chain.
+   * Tropical year of date — the sidereal year less the equinox's own motion
+   * over one year: one turn per the composed lunisolar precession period
+   * T_p(t), i.e. 1/T_trop = 1/T_sid + ψ̇(t)/(1,296,000″ · Julian-year s) with
+   * ψ̇(t) = [ω₀/ω(t)]⁻¹… spelled as (LOD₀/LOD(t))·p₀·[f_S + (1 − f_S)(a₀/a_M)³].
+   * Plan 09 step 5d: it was the unit's calendar convention sid·(1 − 13/H(t))
+   * (plan 06 D2 kept it as a named device), 0.086 % slow on the equinox term
+   * — 1.0 s of year at J2000, integrated by the deep JD↔year calendar. The
+   * SAME operations as deltat/deep-time.cjs tropicalYearSecondsFromLod (the
+   * layer0 gate holds the twins bit-identical). Past the tidal chain's
+   * domain the J2000 rate stands in.
+   * (The former anomalistic identity T_sid·(H − 13)/(H − 16) left with it:
+   * the model publishes the anomalistic year of date inside the one-family
+   * window only, from the series' own apsidal rate.)
    * @param {number} t age in Ma
    * @returns {number} seconds
    */
   const tropicalYearSecondsCore = (t) => {
     const sid = siderealYearSecondsCore(t);
-    const H = holisticHCore(t);
-    return H === null ? sid * (1 - 13 / p.holisticYearJ2000) : sid * (1 - 13 / H);
-  };
-
-  /**
-   * Anomalistic year.
-   *
-   * NOT `T_sid * (1 + 3/H)`. That is the FIRST-ORDER form and the plan quotes
-   * it, but the shipped chain composes the lattice steps exactly: tropical is
-   * `(H-13)/H` of sidereal, then anomalistic is `H/(H-16)` of tropical, giving
-   * `T_sid * (H-13)/(H-16)`. The two differ by `48/H^2 * T_sid` ~ 13.5 ms —
-   * comparable to 6d's own 0.0178 s anomalistic RMSE, so adopting the
-   * first-order form would move a shipped value for no reason. Exact
-   * composition it is; the identity 13 + 3 = 16 still holds to first order.
-   *
-   * @param {number} t age in Ma
-   * @returns {number|null} seconds
-   */
-  const anomalisticYearSecondsCore = (t) => {
-    const H = holisticHCore(t);
-    if (H === null) return null;
-    const sid = siderealYearSecondsCore(t);
-    return sid * (H - 13) / H * H / (H - 16);
+    const lod = lodSecondsCore(t);
+    const a = moonDistanceMetresCore(t);
+    const lf = (a === null || !(a > 0)) ? null : Math.pow(p.moonDistanceNowM / a, 3);
+    const term = lf === null ? null : p.precessionSolarShareJ2000 + (1 - p.precessionSolarShareJ2000) * lf;
+    const p0 = p.precessionRateJ2000ArcsecPerYr;
+    const rate = (lod === null || !(lod > 0) || term === null) ? p0 : (p.lodNowH13Seconds / lod) * p0 * term;
+    return 1 / (1 / sid + rate / (1296000 * 365.25 * 86400));
   };
 
-  /** The frozen comb family's tropical base: T_sid·(1 − 13/H_era). @param {number} t */
-  const eraClockTropicalYearSecondsCore = (t) => {
-    const sid = siderealYearSecondsCore(t);
-    const H = eraClockHCore(t);
-    return H === null ? sid * (1 - 13 / p.holisticYearJ2000) : sid * (1 - 13 / H);
-  };
+  // (plan 09 step 7a: the frozen comb family's tropical base T_sid·(1 − 13/H_era) is deleted)
 
   return Object.freeze({
     tMa,
     moonDistanceMetres: (year) => moonDistanceMetresCore(tMa(year)),
     lodSeconds: (year) => lodSecondsCore(tMa(year)),
-    holisticH: (year) => holisticHCore(tMa(year)),
     siderealYearSeconds: (year) => siderealYearSecondsCore(tMa(year)),
     tropicalYearSeconds: (year) => tropicalYearSecondsCore(tMa(year)),
-    anomalisticYearSeconds: (year) => anomalisticYearSecondsCore(tMa(year)),
-    eraClockH: (year) => eraClockHCore(tMa(year)),
-    eraClockTropicalYearSeconds: (year) => eraClockTropicalYearSecondsCore(tMa(year)),
   });
 };

package/src/layer1/index.js

@@ -21,8 +21,6 @@
  * @typedef {Object} DerivedViews
  * @property {(year: number) => number|null} tropicalYearDays
  * @property {(year: number) => number|null} siderealYearDays
- * @property {(year: number) => number|null} anomalisticYearDays
- * @property {(year: number) => number|null} siderealYearDaysViaLattice
  */
 
 /**
@@ -44,45 +42,16 @@ export const createDerivedViews = ({ primitives: L0 }) => {
 
   const tropicalYearDays = perDay(L0.tropicalYearSeconds);
   const siderealYearDays = perDay(L0.siderealYearSeconds);
-  const anomalisticYearDays = perDay(L0.anomalisticYearSeconds);
 
-  /**
-   * The sidereal year in days, reached the long way — from the tropical year
-   * through the H/13 lattice identity instead of from sidereal seconds.
-   *
-   * ALGEBRAICALLY THE SAME NUMBER. `T_trop = T_sid·(H−13)/H`, so
-   * `(T_trop/LOD)·H/(H−13)` cancels to `T_sid/LOD` — but the operation ORDER
-   * differs, so it is not bit-identical everywhere. Measured across the fixture
-   * grid: identical at six of seven epochs, 0.7 ULP apart at +0.3 Ma. The
-   * browser shows exactly the same one-epoch split, which is how we know the
-   * port is faithful rather than merely close.
-   *
-   * It exists because `src/script.js` carried BOTH as separate mutable globals
-   * — `meansiderealyearlengthinDays` and `meansiderealyearlengthinDays_kinematic`
-   * — holding the same value at every epoch once `recomputeEpochAnchors` had
-   * run, diverging only at module load, where one is the IAU anchor
-   * (365.256363004) and the other lattice-derived: the 118 ms of R16.
-   *
-   * So the "kinematic" variant was never a different quantity, only a
-   * different arrival route with a different starting constant. The year/day
-   * collapse retired the mutable global: `meansiderealyearlengthinDays_kinematic`
-   * is now a J2000 CONSTANT in the browser (this route at f(2000)), and the
-   * epoch value is read from `meansiderealyearlengthinDays`. This view stays
-   * so the equality is measured, not asserted.
-   *
-   * @param {number} year
-   * @returns {number|null}
-   */
-  const siderealYearDaysViaLattice = (year) => {
-    const trop = tropicalYearDays(year);
-    const H = L0.holisticH(year);
-    return (trop === null || H === null) ? null : trop * H / (H - 13);
-  };
+  // (Plan 09 step 5d: `anomalisticYearDays` — the H/16 device identity — and
+  // `siderealYearDaysViaLattice` — the tropical year inverted through the
+  // H/(H − 13) identity, kept to measure the browser's former two sidereal
+  // globals against each other — left with the identities: the tropical year
+  // of date is one equinox turn per T_p(t), and the anomalistic year of date
+  // is published from the series' apsidal rate inside the one-family window.)
 
   return Object.freeze({
     tropicalYearDays,
     siderealYearDays,
-    anomalisticYearDays,
-    siderealYearDaysViaLattice,
   });
 };