ootk 3.0.0 → 4.0.0-1

package/dist/interpolator/CubicSplineInterpolator.js

@@ -0,0 +1,80 @@
+/**
+ * @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 {
+    _splines;
+    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/dist/interpolator/CubicSplineInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CubicSplineInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/CubicSplineInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAA6C,MAAM,WAAW,CAAC;AACpG,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD;;;;;;;GAOG;AACH,MAAM,OAAO,uBAAwB,SAAQ,iBAAiB;IACxC;IAApB,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;aACpB,YAAY,CAAC,KAAK,CAAC;aACnB,WAAW,CAAC,KAAK,CAA0D,CAAC;QAE/E,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/dist/interpolator/FieldInterpolator.d.ts

@@ -0,0 +1,28 @@
+/**
+ * @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?: Float64Array;
+    /**
+     * Interpolate field values at the provided [epoch].
+     * @param final Final state vector.
+     * @param EpochUTC Epoch of the final state vector.
+     * @param epoch Epoch to interpolate field values at.
+     * @throws [Error] if the interpolator is not initialized.
+     */
+    interpolate(final: unknown, EpochUTC: unknown, epoch: unknown): void;
+}

package/dist/interpolator/FieldInterpolator.js

@@ -0,0 +1,34 @@
+/**
+ * @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 {
+    Float64List;
+    /**
+     * Interpolate field values at the provided [epoch].
+     * @param final Final state vector.
+     * @param EpochUTC Epoch of the final state vector.
+     * @param epoch Epoch to interpolate field values at.
+     * @throws [Error] if the interpolator is not initialized.
+     */
+    interpolate(final, EpochUTC, epoch) {
+        throw new Error('Method not implemented.');
+    }
+}
+//# sourceMappingURL=FieldInterpolator.js.map

package/dist/interpolator/FieldInterpolator.js.map

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

package/dist/interpolator/Interpolator.d.ts

@@ -0,0 +1,22 @@
+/**
+ * @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/dist/interpolator/Interpolator.js

@@ -0,0 +1,47 @@
+/**
+ * @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/dist/interpolator/Interpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Interpolator.js","sourceRoot":"","sources":["../../src/interpolator/Interpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAW,MAAM,WAAW,CAAC;AAE3D,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,CAAY,CAAC,CAAC;YACjE,MAAM,EAAE,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAY,CAAC,CAAC;YAEjE,OAAO,IAAI,WAAW,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;SAChC;QAED,OAAO,IAAI,CAAC;IACd,CAAC;CACF"}

package/dist/interpolator/LagrangeInterpolator.d.ts

@@ -0,0 +1,40 @@
+/**
+ * @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;
+}

package/dist/interpolator/LagrangeInterpolator.js

@@ -0,0 +1,159 @@
+/**
+ * @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, Vector3D } from 'ootk-core';
+import { StateInterpolator } from './StateInterpolator';
+export class LagrangeInterpolator extends StateInterpolator {
+    t_;
+    x_;
+    y_;
+    z_;
+    order;
+    constructor(t, x, y, z, order = 10) {
+        super();
+        this.t_ = t;
+        this.x_ = x;
+        this.y_ = y;
+        this.z_ = z;
+        this.order = order;
+    }
+    /**
+     * 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, order = 10) {
+        const k = ephemeris.length;
+        const t = new Float64Array(k);
+        const x = new Float64Array(k);
+        const y = new Float64Array(k);
+        const z = new Float64Array(k);
+        for (let i = 0; i < k; i++) {
+            const state = ephemeris[i];
+            t[i] = state.epoch.posix;
+            x[i] = state.position.x;
+            y[i] = state.position.y;
+            z[i] = state.position.z;
+        }
+        return new LagrangeInterpolator(t, x, y, z, order);
+    }
+    get sizeBytes() {
+        return (64 * 4 * this.t_.length) / 8;
+    }
+    interpolate(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        const posix = epoch.posix;
+        const subDex = this.slice_(posix);
+        const start = subDex.left;
+        const stop = subDex.right;
+        const ts = this.t_.subarray(start, stop);
+        const xs = this.x_.subarray(start, stop);
+        const ys = this.y_.subarray(start, stop);
+        const zs = this.z_.subarray(start, stop);
+        const position = new Vector3D(LagrangeInterpolator.position_(ts, xs, posix), LagrangeInterpolator.position_(ts, ys, posix), LagrangeInterpolator.position_(ts, zs, posix));
+        const velocity = new Vector3D(LagrangeInterpolator.velocity_(ts, xs, posix), LagrangeInterpolator.velocity_(ts, ys, posix), LagrangeInterpolator.velocity_(ts, zs, posix));
+        return new J2000(epoch, position, velocity);
+    }
+    static position_(xs, ys, x) {
+        const k = xs.length - 1;
+        let result = 0.0;
+        for (let j = 0; j < k; j++) {
+            let product = ys[j];
+            for (let m = 0; m < k; m++) {
+                if (j === m) {
+                    continue;
+                }
+                product *= (x - xs[m]) / (xs[j] - xs[m]);
+            }
+            result += product;
+        }
+        return result;
+    }
+    static velocity_(xs, ys, x) {
+        const k = xs.length;
+        let result = 0.0;
+        for (let j = 0; j < k; j++) {
+            let total = 0.0;
+            for (let i = 0; i < k; i++) {
+                if (i === j) {
+                    continue;
+                }
+                let product = 1 / (xs[j] - xs[i]);
+                for (let m = 0; m < k; m++) {
+                    if (m === i || m === j) {
+                        continue;
+                    }
+                    product *= (x - xs[m]) / (xs[j] - xs[m]);
+                }
+                total += product;
+            }
+            result += ys[j] * total;
+        }
+        return result;
+    }
+    static _getClosest(target, t1, d1, t2, d2) {
+        return target - t1 >= t2 - target ? d2 : d1;
+    }
+    slice_(posix) {
+        const n = this.t_.length;
+        if (posix <= this.t_[0]) {
+            return { left: 0, right: this.order };
+        }
+        if (posix >= this.t_[n - 1]) {
+            return { left: n - this.order, right: n };
+        }
+        let i = 0;
+        let j = this.t_.length;
+        let mid = 0;
+        while (i < j) {
+            mid = (i + j) >> 1;
+            if (this.t_[mid] === posix) {
+                break;
+            }
+            if (posix < this.t_[mid]) {
+                if (mid > 0 && posix > this.t_[mid - 1]) {
+                    mid = LagrangeInterpolator._getClosest(posix, this.t_[mid - 1], mid - 1, this.t_[mid], mid);
+                    break;
+                }
+                j = mid;
+            }
+            else {
+                if (mid < this.t_.length - 1 && posix < this.t_[mid + 1]) {
+                    mid = LagrangeInterpolator._getClosest(posix, this.t_[mid], mid, this.t_[mid + 1], mid + 1);
+                    break;
+                }
+                i = mid + 1;
+            }
+        }
+        const offset = Math.floor(this.order / 2);
+        const left = mid - offset;
+        const right = mid + offset - (this.order % 2 === 1 ? 1 : 0);
+        if (left < 0) {
+            return { left: 0, right: this.order };
+        }
+        if (right > n) {
+            return { left: n - this.order, right: n };
+        }
+        return { left, right };
+    }
+    window() {
+        return new EpochWindow(new EpochUTC(this.t_[0]), new EpochUTC(this.t_[this.t_.length - 1]));
+    }
+}
+//# sourceMappingURL=LagrangeInterpolator.js.map

package/dist/interpolator/LagrangeInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"LagrangeInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/LagrangeInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAA4C,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC7G,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD,MAAM,OAAO,oBAAqB,SAAQ,iBAAiB;IACjD,EAAE,CAAe;IACjB,EAAE,CAAe;IACjB,EAAE,CAAe;IACjB,EAAE,CAAe;IACjB,KAAK,CAAS;IAEtB,YAAY,CAAe,EAAE,CAAe,EAAE,CAAe,EAAE,CAAe,EAAE,KAAK,GAAG,EAAE;QACxF,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACZ,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;IACrB,CAAC;IAED;;;;;OAKG;IACH,MAAM,CAAC,aAAa,CAAC,SAAkB,EAAE,KAAK,GAAG,EAAE;QACjD,MAAM,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC;QAC3B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC9B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAE9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAE3B,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC;YACzB,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YACxB,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YACxB,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;SACzB;QAED,OAAO,IAAI,oBAAoB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,KAAK,CAAC,CAAC;IACrD,CAAC;IAED,IAAI,SAAS;QACX,OAAO,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IACvC,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,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC;QAClC,MAAM,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC;QAC1B,MAAM,IAAI,GAAG,MAAM,CAAC,KAAK,CAAC;QAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;QACzC,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;QACzC,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAC3B,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAe,EAC3D,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAe,EAC3D,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAe,CAC5D,CAAC;QACF,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAC3B,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAwB,EACpE,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAwB,EACpE,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAwB,CACrE,CAAC;QAEF,OAAO,IAAI,KAAK,CAAC,KAAK,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAC9C,CAAC;IAEO,MAAM,CAAC,SAAS,CAAC,EAAgB,EAAE,EAAgB,EAAE,CAAS;QACpE,MAAM,CAAC,GAAG,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC;QACxB,IAAI,MAAM,GAAG,GAAG,CAAC;QAEjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,IAAI,OAAO,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;YAEpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,IAAI,CAAC,KAAK,CAAC,EAAE;oBACX,SAAS;iBACV;gBACD,OAAO,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aAC1C;YACD,MAAM,IAAI,OAAO,CAAC;SACnB;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAEO,MAAM,CAAC,SAAS,CAAC,EAAgB,EAAE,EAAgB,EAAE,CAAS;QACpE,MAAM,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC;QACpB,IAAI,MAAM,GAAG,GAAG,CAAC;QAEjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,IAAI,KAAK,GAAG,GAAG,CAAC;YAEhB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,IAAI,CAAC,KAAK,CAAC,EAAE;oBACX,SAAS;iBACV;gBACD,IAAI,OAAO,GAAG,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAElC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;oBAC1B,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE;wBACtB,SAAS;qBACV;oBACD,OAAO,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;iBAC1C;gBACD,KAAK,IAAI,OAAO,CAAC;aAClB;YACD,MAAM,IAAI,EAAE,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;SACzB;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAEO,MAAM,CAAC,WAAW,CAAC,MAAc,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU,EAAE,EAAU;QACvF,OAAO,MAAM,GAAG,EAAE,IAAI,EAAE,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;IAC9C,CAAC;IAEO,MAAM,CAAC,KAAa;QAC1B,MAAM,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC;QAEzB,IAAI,KAAK,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE;YACvB,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC;SACvC;QACD,IAAI,KAAK,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE;YAC3B,OAAO,EAAE,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;SAC3C;QAED,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC;QACvB,IAAI,GAAG,GAAG,CAAC,CAAC;QAEZ,OAAO,CAAC,GAAG,CAAC,EAAE;YACZ,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;YACnB,IAAI,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,KAAK,KAAK,EAAE;gBAC1B,MAAM;aACP;YACD,IAAI,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE;gBACxB,IAAI,GAAG,GAAG,CAAC,IAAI,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE;oBACvC,GAAG,GAAG,oBAAoB,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;oBAC5F,MAAM;iBACP;gBACD,CAAC,GAAG,GAAG,CAAC;aACT;iBAAM;gBACL,IAAI,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,IAAI,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE;oBACxD,GAAG,GAAG,oBAAoB,CAAC,WAAW,CAAC,KAAK,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,CAAC;oBAC5F,MAAM;iBACP;gBACD,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;aACb;SACF;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC1C,MAAM,IAAI,GAAG,GAAG,GAAG,MAAM,CAAC;QAC1B,MAAM,KAAK,GAAG,GAAG,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAE5D,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,OAAO,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC;SACvC;QACD,IAAI,KAAK,GAAG,CAAC,EAAE;YACb,OAAO,EAAE,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC;SAC3C;QAED,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,CAAC;IACzB,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,WAAW,CAAC,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAY,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAY,CAAC,CAAC,CAAC;IACpH,CAAC;CACF"}

package/dist/interpolator/StateInterpolator.d.ts

@@ -0,0 +1,27 @@
+/**
+ * @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, J2000 } from 'ootk-core';
+import { Interpolator } from './Interpolator';
+export declare abstract class StateInterpolator extends Interpolator {
+    /**
+     * Interpolates the state at the given epoch.
+     * @param epoch The epoch in UTC format.
+     * @throws If the interpolator has not been initialized.
+     */
+    interpolate(epoch: EpochUTC): J2000 | null;
+    get sizeBytes(): number;
+}

package/dist/interpolator/StateInterpolator.js

@@ -0,0 +1,35 @@
+/**
+ * @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 state vector interpolators.
+export class StateInterpolator extends Interpolator {
+    /**
+     * Interpolates the state at the given epoch.
+     * @param epoch The epoch in UTC format.
+     * @throws If the interpolator has not been initialized.
+     */
+    interpolate(epoch) {
+        throw new Error('Not implemented.');
+    }
+    // / Return the size _(bytes)_ of this interpolator's cached data.
+    get sizeBytes() {
+        throw new Error('Not implemented.');
+    }
+}
+//# sourceMappingURL=StateInterpolator.js.map

package/dist/interpolator/StateInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"StateInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/StateInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAGH,sDAAsD;AACtD,2CAA2C;AAE3C,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,+CAA+C;AAC/C,MAAM,OAAgB,iBAAkB,SAAQ,YAAY;IAC1D;;;;OAIG;IACH,WAAW,CAAC,KAAe;QACzB,MAAM,IAAI,KAAK,CAAC,kBAAkB,CAAC,CAAC;IACtC,CAAC;IAED,kEAAkE;IAClE,IAAI,SAAS;QACX,MAAM,IAAI,KAAK,CAAC,kBAAkB,CAAC,CAAC;IACtC,CAAC;CACF"}

package/dist/interpolator/VerletBlendInterpolator.d.ts

@@ -0,0 +1,43 @@
+/**
+ * @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 { CubicSplineInterpolator } from './CubicSplineInterpolator';
+import { LagrangeInterpolator } from './LagrangeInterpolator';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Two-body Velocity Verlet Blend interpolator.
+ *
+ * The [VerletBlendInterpolator] retains the original ephemerides, so the
+ * original _"truth"_ states can be retrieved if needed without imparting any
+ * additional error, so this can be used to build other interpolator types.
+ * The implementation is simple and very tolerant when working with sparse
+ * ephemerides.
+ */
+export declare class VerletBlendInterpolator extends StateInterpolator {
+    ephemeris: J2000[];
+    constructor(ephemeris: J2000[]);
+    get sizeBytes(): number;
+    window(): EpochWindow;
+    private static getClosest_;
+    private matchState_;
+    private static _gravity;
+    private static integrate_;
+    interpolate(epoch: EpochUTC): J2000 | null;
+    getCachedState(epoch: EpochUTC): J2000 | null;
+    toCubicSpline(): CubicSplineInterpolator;
+    toLagrange(order?: number): LagrangeInterpolator;
+}