@essrt/physics 4.7.0 → 4.9.0

package/types/planets/secular-series.d.cts

@@ -0,0 +1,73 @@
+/**
+ * @param {{
+ *   series: {t0Yr: number,
+ *            bodies: Object<string, {stepYr: number,
+ *              zetaQ: ReadonlyArray<number>, zetaP: ReadonlyArray<number>,
+ *              zQ: ReadonlyArray<number>, zP: ReadonlyArray<number>}>,
+ *            verdict?: {planetHandover?: {rows?: Object<string, {pastKyr: number, futureKyr: number}>}}},
+ *   anchorElements: Object<string, {e: number, lonPeriEclipticDeg: number,
+ *     inclEclipticDeg: number, ascNodeEclipticDeg: number}>,
+ *   invariablePlane?: {inclEclipticDeg: number, ascNodeEclipticDeg: number},
+ *   planetZModes: Object<string, ReadonlyArray<{omegaRadPerYr:number,re:number,im:number}>>,
+ *   planetZetaModes: Object<string, ReadonlyArray<{omegaRadPerYr:number,re:number,im:number}>>,
+ * }} deps — the parsed series artifact, the chain's J2000 anchors + banked
+ *   invariable plane, and the deep mode tables (the beyond-span tail).
+ */
+export function createSecularSeriesOverride({ series, anchorElements, invariablePlane, planetZModes, planetZetaModes }: {
+    series: {
+        t0Yr: number;
+        bodies: {
+            [x: string]: {
+                stepYr: number;
+                zetaQ: ReadonlyArray<number>;
+                zetaP: ReadonlyArray<number>;
+                zQ: ReadonlyArray<number>;
+                zP: ReadonlyArray<number>;
+            };
+        };
+        verdict?: {
+            planetHandover?: {
+                rows?: {
+                    [x: string]: {
+                        pastKyr: number;
+                        futureKyr: number;
+                    };
+                };
+            };
+        };
+    };
+    anchorElements: {
+        [x: string]: {
+            e: number;
+            lonPeriEclipticDeg: number;
+            inclEclipticDeg: number;
+            ascNodeEclipticDeg: number;
+        };
+    };
+    invariablePlane?: {
+        inclEclipticDeg: number;
+        ascNodeEclipticDeg: number;
+    };
+    planetZModes: {
+        [x: string]: ReadonlyArray<{
+            omegaRadPerYr: number;
+            re: number;
+            im: number;
+        }>;
+    };
+    planetZetaModes: {
+        [x: string]: ReadonlyArray<{
+            omegaRadPerYr: number;
+            re: number;
+            im: number;
+        }>;
+    };
+}): {
+    applyToElements: (nm: string, year: number, el: any) => any;
+};
+/** Anchored mode-sum: z(t) = Σ (re + i·im)·e^{iωt}, t years from J2000. */
+export function modeSum(modes: ReadonlyArray<{
+    omegaRadPerYr: number;
+    re: number;
+    im: number;
+}>, t: number): number[];

package/types/sun/longitude-correction.d.cts

@@ -0,0 +1,41 @@
+export type SunLongitudeCorrectionDeps = {
+    /**
+     * - H, J2000-fixed (the fitted axis)
+     */
+    hYears: number;
+    /**
+     * - J2000-fixed
+     */
+    balancedYear: number;
+    j2000JD: number;
+    /**
+     * - the fitted mean offset
+     */
+    meanDeg: number;
+    /**
+     * - [divisor, sinC, cosC]
+     */
+    harmonics: Array<[number, number, number]>;
+    /**
+     * - lunar nodal divisor (per 8H convention)
+     */
+    nNodalJ2000: number;
+    /**
+     * - lunar apsidal divisor
+     */
+    nApsidalJ2000: number;
+};
+/**
+ * @typedef {Object} SunLongitudeCorrectionDeps
+ * @property {number} hYears - H, J2000-fixed (the fitted axis)
+ * @property {number} balancedYear - J2000-fixed
+ * @property {number} j2000JD
+ * @property {number} meanDeg - the fitted mean offset
+ * @property {Array<[number, number, number]>} harmonics - [divisor, sinC, cosC]
+ * @property {number} nNodalJ2000 - lunar nodal divisor (per 8H convention)
+ * @property {number} nApsidalJ2000 - lunar apsidal divisor
+ */
+/** @param {SunLongitudeCorrectionDeps} deps */
+export function createSunLongitudeCorrection(deps: SunLongitudeCorrectionDeps): {
+    correctionDegAt: (jd: number) => number;
+};