ootk 3.0.0 → 4.0.0-0

package/lib/interpolator/ChebyshevInterpolator.d.ts

@@ -0,0 +1,39 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow, J2000 } from 'ootk-core';
+import { ChebyshevCoefficients } from './ChebyshevCoefficients';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Compressed Chebyshev ephemeris interpolator.
+ *
+ * The ChebyshevInterpolator sacrifices state accuracy in order to gain
+ * speed and memory requirements, so you can store many ephemerides in RAM
+ * and interpolate through them quickly.
+ *
+ * Using more coefficients per revolution during lossy compression results
+ * in increased accuracy and decreased performance.
+ */
+export declare class ChebyshevInterpolator extends StateInterpolator {
+    private _coefficients;
+    constructor(coefficients: ChebyshevCoefficients[]);
+    private _calcSizeBytes;
+    get sizeBytes(): number;
+    interpolate(epoch: EpochUTC): J2000 | null;
+    window(): EpochWindow;
+    private _matchCoefficients;
+}

package/lib/interpolator/ChebyshevInterpolator.js

@@ -0,0 +1,71 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow, J2000 } from 'ootk-core';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Compressed Chebyshev ephemeris interpolator.
+ *
+ * The ChebyshevInterpolator sacrifices state accuracy in order to gain
+ * speed and memory requirements, so you can store many ephemerides in RAM
+ * and interpolate through them quickly.
+ *
+ * Using more coefficients per revolution during lossy compression results
+ * in increased accuracy and decreased performance.
+ */
+export class ChebyshevInterpolator extends StateInterpolator {
+    constructor(coefficients) {
+        super();
+        this._coefficients = coefficients;
+    }
+    _calcSizeBytes() {
+        let output = 0;
+        for (const coeffs of this._coefficients) {
+            output += coeffs.sizeBytes;
+        }
+        return output;
+    }
+    get sizeBytes() {
+        return this._calcSizeBytes();
+    }
+    interpolate(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        const coeffs = this._matchCoefficients(epoch.posix);
+        const { position: pos, velocity: vel } = coeffs.interpolate(epoch.posix);
+        return new J2000(epoch, pos, vel);
+    }
+    window() {
+        return new EpochWindow(new EpochUTC(this._coefficients[0].a), new EpochUTC(this._coefficients[this._coefficients.length - 1].b));
+    }
+    _matchCoefficients(posix) {
+        let left = 0;
+        let right = this._coefficients.length;
+        while (left < right) {
+            const middle = (left + right) >> 1;
+            if (this._coefficients[middle].b < posix) {
+                left = middle + 1;
+            }
+            else {
+                right = middle;
+            }
+        }
+        return this._coefficients[left];
+    }
+}
+//# sourceMappingURL=ChebyshevInterpolator.js.map

package/lib/interpolator/ChebyshevInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ChebyshevInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/ChebyshevInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,WAAW,CAAC;AAEzD,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD;;;;;;;;;GASG;AACH,MAAM,OAAO,qBAAsB,SAAQ,iBAAiB;IAG1D,YAAY,YAAqC;QAC/C,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;IACpC,CAAC;IAEO,cAAc;QACpB,IAAI,MAAM,GAAG,CAAC,CAAC;QAEf,KAAK,MAAM,MAAM,IAAI,IAAI,CAAC,aAAa,EAAE;YACvC,MAAM,IAAI,MAAM,CAAC,SAAS,CAAC;SAC5B;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,IAAI,SAAS;QACX,OAAO,IAAI,CAAC,cAAc,EAAE,CAAC;IAC/B,CAAC;IAED,WAAW,CAAC,KAAe;QACzB,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,kBAAkB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACpD,MAAM,EAAE,QAAQ,EAAE,GAAG,EAAE,QAAQ,EAAE,GAAG,EAAE,GAAG,MAAM,CAAC,WAAW,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QAEzE,OAAO,IAAI,KAAK,CAAC,KAAK,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;IACpC,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,WAAW,CACpB,IAAI,QAAQ,CAAC,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EACrC,IAAI,QAAQ,CAAC,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAClE,CAAC;IACJ,CAAC;IAEO,kBAAkB,CAAC,KAAa;QACtC,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,KAAK,GAAG,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC;QAEtC,OAAO,IAAI,GAAG,KAAK,EAAE;YACnB,MAAM,MAAM,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YAEnC,IAAI,IAAI,CAAC,aAAa,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE;gBACxC,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC;aACnB;iBAAM;gBACL,KAAK,GAAG,MAAM,CAAC;aAChB;SACF;QAED,OAAO,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC;IAClC,CAAC;CACF"}

package/lib/interpolator/CubicSpline.d.ts

@@ -0,0 +1,36 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Vector3D } from 'ootk-core';
+export declare class CubicSpline {
+    t0: number;
+    p0: Vector3D;
+    m0: Vector3D;
+    t1: number;
+    p1: Vector3D;
+    m1: Vector3D;
+    constructor(t0: number, p0: Vector3D, m0: Vector3D, t1: number, p1: Vector3D, m1: Vector3D);
+    private _position;
+    private _velocity;
+    /**
+     * Interpolates the position and velocity at a given time.
+     * (km) and velocity (km/s) vectors at the provided time.
+     * @param t The time value to interpolate at _(POSIX seconds)_.
+     * @returns An array containing the interpolated position and velocity as Vector3D objects.
+     */
+    interpolate(t: number): Vector3D[];
+}

package/lib/interpolator/CubicSpline.js

@@ -0,0 +1,64 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+// / Container for cubic spline data.
+export class CubicSpline {
+    // / Create a new [CubicSpline] object.
+    constructor(t0, p0, m0, t1, p1, m1) {
+        this.t0 = t0;
+        this.p0 = p0;
+        this.m0 = m0;
+        this.t1 = t1;
+        this.p1 = p1;
+        this.m1 = m1;
+        // Nothing to do here.
+    }
+    // / Interpolate position at the provided time [t] _(POSIX seconds)_.
+    _position(t) {
+        const t2 = t * t;
+        const t3 = t2 * t;
+        const r0 = this.p0.scale(2 * t3 - 3 * t2 + 1);
+        const v0 = this.m0.scale((t3 - 2 * t2 + t) * (this.t1 - this.t0));
+        const r1 = this.p1.scale(-2 * t3 + 3 * t2);
+        const v1 = this.m1.scale((t3 - t2) * (this.t1 - this.t0));
+        return r0.add(v0).add(r1).add(v1);
+    }
+    // / Interpolate velocity at the provided time [t] _(POSIX seconds)_.
+    _velocity(t) {
+        const t2 = t * t;
+        const r0 = this.p0.scale(6 * t2 - 6 * t);
+        const v0 = this.m0.scale((3 * t2 - 4 * t + 1) * (this.t1 - this.t0));
+        const r1 = this.p1.scale(-6 * t2 + 6 * t);
+        const v1 = this.m1.scale((3 * t2 - 2 * t) * (this.t1 - this.t0));
+        return r0
+            .add(v0)
+            .add(r1)
+            .add(v1)
+            .scale(1 / (this.t1 - this.t0));
+    }
+    /**
+     * Interpolates the position and velocity at a given time.
+     * (km) and velocity (km/s) vectors at the provided time.
+     * @param t The time value to interpolate at _(POSIX seconds)_.
+     * @returns An array containing the interpolated position and velocity as Vector3D objects.
+     */
+    interpolate(t) {
+        const n = (t - this.t0) / (this.t1 - this.t0);
+        return [this._position(n), this._velocity(n)];
+    }
+}
+//# sourceMappingURL=CubicSpline.js.map

package/lib/interpolator/CubicSpline.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CubicSpline.js","sourceRoot":"","sources":["../../src/interpolator/CubicSpline.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAIH,qCAAqC;AACrC,MAAM,OAAO,WAAW;IACtB,uCAAuC;IACvC,YACS,EAAU,EACV,EAAY,EACZ,EAAY,EACZ,EAAU,EACV,EAAY,EACZ,EAAY;QALZ,OAAE,GAAF,EAAE,CAAQ;QACV,OAAE,GAAF,EAAE,CAAU;QACZ,OAAE,GAAF,EAAE,CAAU;QACZ,OAAE,GAAF,EAAE,CAAQ;QACV,OAAE,GAAF,EAAE,CAAU;QACZ,OAAE,GAAF,EAAE,CAAU;QAEnB,sBAAsB;IACxB,CAAC;IAED,qEAAqE;IAC7D,SAAS,CAAC,CAAS;QACzB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAClB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAClE,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;QAC3C,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1D,OAAO,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;IACpC,CAAC;IAED,qEAAqE;IAC7D,SAAS,CAAC,CAAS;QACzB,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QACrE,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjE,OAAO,EAAE;aACN,GAAG,CAAC,EAAE,CAAC;aACP,GAAG,CAAC,EAAE,CAAC;aACP,GAAG,CAAC,EAAE,CAAC;aACP,KAAK,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC;IACpC,CAAC;IAED;;;;;OAKG;IACH,WAAW,CAAC,CAAS;QACnB,MAAM,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAE9C,OAAO,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC;IAChD,CAAC;CACF"}

package/lib/interpolator/CubicSplineInterpolator.d.ts

@@ -0,0 +1,37 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow, J2000 } from 'ootk-core';
+import { CubicSpline } from './CubicSpline';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Cubic spline ephemeris interpolator.
+ *
+ * The [CubicSplineInterpolator] is a very fast and accurate interpolator
+ * at the expense of memory due to the cached spline pairs used in the
+ * interpolation operation. Accuracy is significantly impacted when using
+ * sparse ephemerides.
+ */
+export declare class CubicSplineInterpolator extends StateInterpolator {
+    private _splines;
+    constructor(_splines: CubicSpline[]);
+    static fromEphemeris(ephemeris: J2000[]): CubicSplineInterpolator;
+    get sizeBytes(): number;
+    private _matchSpline;
+    interpolate(epoch: EpochUTC): J2000 | null;
+    window(): EpochWindow;
+}

package/lib/interpolator/CubicSplineInterpolator.js

@@ -0,0 +1,78 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow, J2000 } from 'ootk-core';
+import { CubicSpline } from './CubicSpline';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Cubic spline ephemeris interpolator.
+ *
+ * The [CubicSplineInterpolator] is a very fast and accurate interpolator
+ * at the expense of memory due to the cached spline pairs used in the
+ * interpolation operation. Accuracy is significantly impacted when using
+ * sparse ephemerides.
+ */
+export class CubicSplineInterpolator extends StateInterpolator {
+    constructor(_splines) {
+        super();
+        this._splines = _splines;
+    }
+    static fromEphemeris(ephemeris) {
+        const splines = [];
+        for (let i = 0; i < ephemeris.length - 1; i++) {
+            const e0 = ephemeris[i];
+            const t0 = e0.epoch.posix;
+            const p0 = e0.position;
+            const m0 = e0.velocity;
+            const e1 = ephemeris[i + 1];
+            const t1 = e1.epoch.posix;
+            const p1 = e1.position;
+            const m1 = e1.velocity;
+            splines.push(new CubicSpline(t0, p0, m0, t1, p1, m1));
+        }
+        return new CubicSplineInterpolator(splines);
+    }
+    get sizeBytes() {
+        return (64 * 14 * this._splines.length) / 8;
+    }
+    _matchSpline(posix) {
+        let left = 0;
+        let right = this._splines.length;
+        while (left < right) {
+            const middle = (left + right) >> 1;
+            if (this._splines[middle].t1 < posix) {
+                left = middle + 1;
+            }
+            else {
+                right = middle;
+            }
+        }
+        return this._splines[left];
+    }
+    interpolate(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        const posix = epoch.posix;
+        const splineVecs = this._matchSpline(posix).interpolate(posix);
+        return new J2000(epoch, splineVecs[0], splineVecs[1]);
+    }
+    window() {
+        return new EpochWindow(new EpochUTC(this._splines[0].t0), new EpochUTC(this._splines[this._splines.length - 1].t1));
+    }
+}
+//# sourceMappingURL=CubicSplineInterpolator.js.map

package/lib/interpolator/CubicSplineInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CubicSplineInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/CubicSplineInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,WAAW,CAAC;AACzD,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD;;;;;;;GAOG;AACH,MAAM,OAAO,uBAAwB,SAAQ,iBAAiB;IAC5D,YAAoB,QAAuB;QACzC,KAAK,EAAE,CAAC;QADU,aAAQ,GAAR,QAAQ,CAAe;IAE3C,CAAC;IAED,MAAM,CAAC,aAAa,CAAC,SAAkB;QACrC,MAAM,OAAO,GAAkB,EAAE,CAAC;QAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC7C,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC;YAC1B,MAAM,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC;YACvB,MAAM,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC;YACvB,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC5B,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC;YAC1B,MAAM,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC;YACvB,MAAM,EAAE,GAAG,EAAE,CAAC,QAAQ,CAAC;YAEvB,OAAO,CAAC,IAAI,CAAC,IAAI,WAAW,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;SACvD;QAED,OAAO,IAAI,uBAAuB,CAAC,OAAO,CAAC,CAAC;IAC9C,CAAC;IAED,IAAI,SAAS;QACX,OAAO,CAAC,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC9C,CAAC;IAEO,YAAY,CAAC,KAAa;QAChC,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC;QAEjC,OAAO,IAAI,GAAG,KAAK,EAAE;YACnB,MAAM,MAAM,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;YAEnC,IAAI,IAAI,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,EAAE,GAAG,KAAK,EAAE;gBACpC,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC;aACnB;iBAAM;gBACL,KAAK,GAAG,MAAM,CAAC;aAChB;SACF;QAED,OAAO,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC;IAC7B,CAAC;IAED,WAAW,CAAC,KAAe;QACzB,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QACD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC;QAC1B,MAAM,UAAU,GAAG,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;QAE/D,OAAO,IAAI,KAAK,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;IACxD,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,WAAW,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IACtH,CAAC;CACF"}

package/lib/interpolator/FieldInterpolator.d.ts

@@ -0,0 +1,26 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Interpolator } from './Interpolator';
+export declare abstract class FieldInterpolator extends Interpolator {
+    Float64List?: any;
+    /**
+     * Interpolate field values at the provided [epoch].
+     * @returns Returns `null` if the epoch is outside the cached field range.
+     */
+    interpolate(final: any, EpochUTC: any, epoch: any): void;
+}

package/lib/interpolator/FieldInterpolator.js

@@ -0,0 +1,31 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+/* eslint-disable @typescript-eslint/no-unused-vars */
+/* eslint-disable class-methods-use-this */
+import { Interpolator } from './Interpolator';
+// / Base class for arbitrary field interpolators.
+export class FieldInterpolator extends Interpolator {
+    /**
+     * Interpolate field values at the provided [epoch].
+     * @returns Returns `null` if the epoch is outside the cached field range.
+     */
+    interpolate(final, EpochUTC, epoch) {
+        throw new Error('Method not implemented.');
+    }
+}
+//# sourceMappingURL=FieldInterpolator.js.map

package/lib/interpolator/FieldInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FieldInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/FieldInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,sDAAsD;AACtD,2CAA2C;AAE3C,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,kDAAkD;AAClD,MAAM,OAAgB,iBAAkB,SAAQ,YAAY;IAG1D;;;OAGG;IACH,WAAW,CAAC,KAAK,EAAE,QAAQ,EAAE,KAAK;QAChC,MAAM,IAAI,KAAK,CAAC,yBAAyB,CAAC,CAAC;IAC7C,CAAC;CACF"}

package/lib/interpolator/Interpolator.d.ts

@@ -0,0 +1,23 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow } from 'ootk-core';
+export declare abstract class Interpolator {
+    abstract window(): EpochWindow;
+    inWindow(epoch: EpochUTC): boolean;
+    overlap(interpolator: Interpolator): EpochWindow | null;
+}

package/lib/interpolator/Interpolator.js

@@ -0,0 +1,48 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow } from 'ootk-core';
+// / Interpolator base class.
+export class Interpolator {
+    /*
+     * Return `true` if the provided [epoch] is within this interpolator's
+     * cached value range.
+     */
+    inWindow(epoch) {
+        const start = this.window().start;
+        const stop = this.window().end;
+        return start <= epoch && epoch <= stop;
+    }
+    /*
+     * Calculate the start/stop epoch between this and another [Interpolator].
+     *
+     * Returns `null` if there is no overlap between interpolators.
+     */
+    overlap(interpolator) {
+        const x1 = this.window().start;
+        const x2 = this.window().end;
+        const y1 = interpolator.window().start;
+        const y2 = interpolator.window().end;
+        if (x1 <= y2 && y1 <= x2) {
+            const e1 = new EpochUTC(Math.max(x1.posix, y1.posix));
+            const e2 = new EpochUTC(Math.min(x2.posix, y2.posix));
+            return new EpochWindow(e1, e2);
+        }
+        return null;
+    }
+}
+//# sourceMappingURL=Interpolator.js.map

package/lib/interpolator/Interpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Interpolator.js","sourceRoot":"","sources":["../../src/interpolator/Interpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,MAAM,WAAW,CAAC;AAElD,6BAA6B;AAC7B,MAAM,OAAgB,YAAY;IAIhC;;;OAGG;IACH,QAAQ,CAAC,KAAe;QACtB,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,KAAK,CAAC;QAClC,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC;QAE/B,OAAO,KAAK,IAAI,KAAK,IAAI,KAAK,IAAI,IAAI,CAAC;IACzC,CAAC;IAED;;;;OAIG;IACH,OAAO,CAAC,YAA0B;QAChC,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,KAAK,CAAC;QAC/B,MAAM,EAAE,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC;QAC7B,MAAM,EAAE,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC,KAAK,CAAC;QACvC,MAAM,EAAE,GAAG,YAAY,CAAC,MAAM,EAAE,CAAC,GAAG,CAAC;QAErC,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE;YACxB,MAAM,EAAE,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC;YACtD,MAAM,EAAE,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC;YAEtD,OAAO,IAAI,WAAW,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;SAChC;QAED,OAAO,IAAI,CAAC;IACd,CAAC;CACF"}

package/lib/interpolator/LagrangeInterpolator.d.ts

@@ -0,0 +1,42 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { EpochUTC, EpochWindow, J2000 } from 'ootk-core';
+import { StateInterpolator } from './StateInterpolator';
+export declare class LagrangeInterpolator extends StateInterpolator {
+    private _t;
+    private _x;
+    private _y;
+    private _z;
+    private order;
+    constructor(t: Float64Array, x: Float64Array, y: Float64Array, z: Float64Array, order?: number);
+    /**
+     * Creates a LagrangeInterpolator from an array of J2000 ephemeris data.
+     *
+     * @param ephemeris - The array of J2000 ephemeris data.
+     * @param order - The order of the LagrangeInterpolator. Default is 10.
+     * @returns A new LagrangeInterpolator instance.
+     */
+    static fromEphemeris(ephemeris: J2000[], order?: number): LagrangeInterpolator;
+    get sizeBytes(): number;
+    interpolate(epoch: EpochUTC): J2000 | null;
+    private static _position;
+    private static _velocity;
+    private static _getClosest;
+    private _slice;
+    window(): EpochWindow;
+}