ootk 2.0.1 → 4.0.0-0

package/lib/interpolator/LagrangeInterpolator.js

@@ -0,0 +1,156 @@
+/**
+ * @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 {
+    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/lib/interpolator/LagrangeInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"LagrangeInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/LagrangeInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAE,WAAW,EAAE,KAAK,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AACnE,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD,MAAM,OAAO,oBAAqB,SAAQ,iBAAiB;IAOzD,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;;;;;;OAMG;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,CAAC,EAC7C,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,EAC7C,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAC9C,CAAC;QACF,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAC3B,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,EAC7C,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,EAC7C,oBAAoB,CAAC,SAAS,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAC9C,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,CAAC,CAAC,EAAE,IAAI,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAC9F,CAAC;CACF"}

package/lib/interpolator/StateInterpolator.d.ts

@@ -0,0 +1,29 @@
+/**
+ * @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.
+     * @returns The interpolated state at the given epoch, or null if interpolation is not possible.
+     */
+    interpolate(epoch: EpochUTC): J2000 | null;
+    get sizeBytes(): number;
+}

package/lib/interpolator/StateInterpolator.js

@@ -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/>.
+ */
+/* 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.
+     * @returns The interpolated state at the given epoch, or null if interpolation is not possible.
+     */
+    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/lib/interpolator/StateInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"StateInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/StateInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAGH,sDAAsD;AACtD,2CAA2C;AAE3C,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAE9C,+CAA+C;AAC/C,MAAM,OAAgB,iBAAkB,SAAQ,YAAY;IAC1D;;;;;OAKG;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/lib/interpolator/VerletBlendInterpolator.d.ts

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

package/lib/interpolator/VerletBlendInterpolator.js

@@ -0,0 +1,115 @@
+/**
+ * @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 { copySign, Earth, 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 class VerletBlendInterpolator extends StateInterpolator {
+    constructor(ephemeris) {
+        super();
+        this.ephemeris = ephemeris;
+    }
+    get sizeBytes() {
+        return (64 * 7 * this.ephemeris.length) / 8;
+    }
+    window() {
+        return new EpochWindow(this.ephemeris[0].epoch, this.ephemeris[this.ephemeris.length - 1].epoch);
+    }
+    static _getClosest(target, s1, s2) {
+        return target - s1.epoch.posix >= s2.epoch.posix - target ? s2 : s1;
+    }
+    _matchState(epoch) {
+        const target = epoch.posix;
+        if (target <= this.ephemeris[0].epoch.posix) {
+            return this.ephemeris[0];
+        }
+        if (target >= this.ephemeris[this.ephemeris.length - 1].epoch.posix) {
+            return this.ephemeris[this.ephemeris.length - 1];
+        }
+        let i = 0;
+        let j = this.ephemeris.length;
+        let mid = 0;
+        while (i < j) {
+            mid = (i + j) >> 1;
+            if (this.ephemeris[mid].epoch.posix === target) {
+                return this.ephemeris[mid];
+            }
+            if (target < this.ephemeris[mid].epoch.posix) {
+                if (mid > 0 && target > this.ephemeris[mid - 1].epoch.posix) {
+                    return VerletBlendInterpolator._getClosest(target, this.ephemeris[mid - 1], this.ephemeris[mid]);
+                }
+                j = mid;
+            }
+            else {
+                if (mid < this.ephemeris.length - 1 && target < this.ephemeris[mid + 1].epoch.posix) {
+                    return VerletBlendInterpolator._getClosest(target, this.ephemeris[mid], this.ephemeris[mid + 1]);
+                }
+                i = mid + 1;
+            }
+        }
+        return this.ephemeris[mid];
+    }
+    static _gravity(position) {
+        const r = position.magnitude();
+        return position.scale(-Earth.mu / (r * r * r));
+    }
+    static _integrate(state, step) {
+        const x0 = state.position;
+        const a0 = VerletBlendInterpolator._gravity(x0);
+        const v0 = state.velocity;
+        const x1 = x0.add(v0.scale(step)).add(a0.scale(0.5 * step * step));
+        const a1 = VerletBlendInterpolator._gravity(x1);
+        const v1 = v0.add(a0.add(a1).scale(0.5 * step));
+        return new J2000(state.epoch.roll(step), x1, v1);
+    }
+    interpolate(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        let state = this._matchState(epoch);
+        while (state.epoch.posix !== epoch.posix) {
+            const delta = epoch.posix - state.epoch.posix;
+            const stepMag = Math.min(5.0, Math.abs(delta));
+            const stepSize = copySign(stepMag, delta);
+            state = VerletBlendInterpolator._integrate(state, stepSize);
+        }
+        return state;
+    }
+    getCachedState(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        return this._matchState(epoch);
+    }
+    toCubicSpline() {
+        return CubicSplineInterpolator.fromEphemeris(this.ephemeris);
+    }
+    toLagrange(order = 10) {
+        return LagrangeInterpolator.fromEphemeris(this.ephemeris, order);
+    }
+}
+//# sourceMappingURL=VerletBlendInterpolator.js.map

package/lib/interpolator/VerletBlendInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"VerletBlendInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/VerletBlendInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAY,WAAW,EAAE,KAAK,EAAY,MAAM,WAAW,CAAC;AACpF,OAAO,EAAE,uBAAuB,EAAE,MAAM,2BAA2B,CAAC;AACpE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD;;;;;;;;GAQG;AACH,MAAM,OAAO,uBAAwB,SAAQ,iBAAiB;IAC5D,YAAmB,SAAkB;QACnC,KAAK,EAAE,CAAC;QADS,cAAS,GAAT,SAAS,CAAS;IAErC,CAAC;IAED,IAAI,SAAS;QACX,OAAO,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC9C,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,WAAW,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IACnG,CAAC;IAEO,MAAM,CAAC,WAAW,CAAC,MAAc,EAAE,EAAS,EAAE,EAAS;QAC7D,OAAO,MAAM,GAAG,EAAE,CAAC,KAAK,CAAC,KAAK,IAAI,EAAE,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;IACtE,CAAC;IAEO,WAAW,CAAC,KAAe;QACjC,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC;QAE3B,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;YAC3C,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;SAC1B;QACD,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;YACnE,OAAO,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;SAClD;QAED,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC;QAC9B,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,SAAS,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,KAAK,MAAM,EAAE;gBAC9C,OAAO,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;aAC5B;YACD,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;gBAC5C,IAAI,GAAG,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;oBAC3D,OAAO,uBAAuB,CAAC,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC;iBAClG;gBACD,CAAC,GAAG,GAAG,CAAC;aACT;iBAAM;gBACL,IAAI,GAAG,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;oBACnF,OAAO,uBAAuB,CAAC,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;iBAClG;gBACD,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;aACb;SACF;QAED,OAAO,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;IAC7B,CAAC;IAEO,MAAM,CAAC,QAAQ,CAAC,QAAkB;QACxC,MAAM,CAAC,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;QAE/B,OAAO,QAAQ,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC;IAEO,MAAM,CAAC,UAAU,CAAC,KAAY,EAAE,IAAY;QAClD,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,MAAM,EAAE,GAAG,uBAAuB,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAChD,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;QACnE,MAAM,EAAE,GAAG,uBAAuB,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAChD,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC;QAEhD,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACnD,CAAC;IAED,WAAW,CAAC,KAAe;QACzB,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QACD,IAAI,KAAK,GAAG,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;QAEpC,OAAO,KAAK,CAAC,KAAK,CAAC,KAAK,KAAK,KAAK,CAAC,KAAK,EAAE;YACxC,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC;YAC9C,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;YAC/C,MAAM,QAAQ,GAAG,QAAQ,CAAC,OAAO,EAAE,KAAK,CAAC,CAAC;YAE1C,KAAK,GAAG,uBAAuB,CAAC,UAAU,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;SAC7D;QAED,OAAO,KAAK,CAAC;IACf,CAAC;IAED,cAAc,CAAC,KAAe;QAC5B,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QAED,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;IACjC,CAAC;IAED,aAAa;QACX,OAAO,uBAAuB,CAAC,aAAa,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;IAED,UAAU,CAAC,KAAK,GAAG,EAAE;QACnB,OAAO,oBAAoB,CAAC,aAAa,CAAC,IAAI,CAAC,SAAS,EAAE,KAAK,CAAC,CAAC;IACnE,CAAC;CACF"}

package/lib/maneuver/TwoBurnOrbitTransfer.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 } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+export declare class TwoBurnOrbitTransfer {
+    vInit: number;
+    vFinal: number;
+    vTransA: number;
+    vTransB: number;
+    tTrans: number;
+    constructor(vInit: number, vFinal: number, vTransA: number, vTransB: number, tTrans: number);
+    /**
+     * Calculates the parameters for a Hohmann transfer orbit between two circular orbits.
+     *
+     * @param rInit The initial radius of the orbit. (km)
+     * @param rFinal The final radius of the orbit. (km)
+     * @returns An instance of TwoBurnOrbitTransfer containing the calculated parameters.
+     */
+    static hohmannTransfer(rInit: number, rFinal: number): TwoBurnOrbitTransfer;
+    get deltaV(): number;
+    /**
+     * Calculates the two maneuver thrusts required for a two-burn orbit transfer.
+     * @param epoch The epoch of the maneuver.
+     * @param durationRate The duration rate of the maneuver (s/m/s).
+     * @returns An array containing the two thrust objects representing the maneuvers.
+     */
+    toManeuvers(epoch: EpochUTC, durationRate?: number): [Thrust, Thrust];
+}

package/lib/maneuver/TwoBurnOrbitTransfer.js

@@ -0,0 +1,62 @@
+/**
+ * @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 { Earth } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+// / Container for a two-burn orbit transfer.
+export class TwoBurnOrbitTransfer {
+    // / Create a new [TwoBurnOrbitTransfer] object.
+    constructor(vInit, vFinal, vTransA, vTransB, tTrans) {
+        this.vInit = vInit;
+        this.vFinal = vFinal;
+        this.vTransA = vTransA;
+        this.vTransB = vTransB;
+        this.tTrans = tTrans;
+        // Nothing to do here.
+    }
+    /**
+     * Calculates the parameters for a Hohmann transfer orbit between two circular orbits.
+     *
+     * @param rInit The initial radius of the orbit. (km)
+     * @param rFinal The final radius of the orbit. (km)
+     * @returns An instance of TwoBurnOrbitTransfer containing the calculated parameters.
+     */
+    static hohmannTransfer(rInit, rFinal) {
+        const vInit = Math.sqrt(Earth.mu / rInit);
+        const vFinal = Math.sqrt(Earth.mu / rFinal);
+        const vTransA = vInit * (Math.sqrt((2.0 * rFinal) / (rInit + rFinal)) - 1.0);
+        const vTransB = vFinal * (1.0 - Math.sqrt((2 * rInit) / (rInit + rFinal)));
+        const tTrans = Math.PI * Math.sqrt((rInit + rFinal) ** 3 / (8.0 * Earth.mu));
+        return new TwoBurnOrbitTransfer(vInit, vFinal, vTransA, vTransB, tTrans);
+    }
+    // / Return the total delta-velocity magnitude for both maneuvers _(km/s)_.
+    get deltaV() {
+        return Math.abs(this.vTransA) + Math.abs(this.vTransB);
+    }
+    /**
+     * Calculates the two maneuver thrusts required for a two-burn orbit transfer.
+     * @param epoch The epoch of the maneuver.
+     * @param durationRate The duration rate of the maneuver (s/m/s).
+     * @returns An array containing the two thrust objects representing the maneuvers.
+     */
+    toManeuvers(epoch, durationRate = 0.0) {
+        const mA = new Thrust(epoch, 0.0, this.vTransA * 1000.0, 0.0, durationRate);
+        const mB = new Thrust(epoch.roll(this.tTrans), 0.0, this.vTransB * 1000.0, 0.0, durationRate);
+        return [mA, mB];
+    }
+}
+//# sourceMappingURL=TwoBurnOrbitTransfer.js.map

package/lib/maneuver/TwoBurnOrbitTransfer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"TwoBurnOrbitTransfer.js","sourceRoot":"","sources":["../../src/maneuver/TwoBurnOrbitTransfer.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAY,MAAM,WAAW,CAAC;AAC5C,OAAO,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAEzC,6CAA6C;AAC7C,MAAM,OAAO,oBAAoB;IAC/B,gDAAgD;IAChD,YACS,KAAa,EACb,MAAc,EACd,OAAe,EACf,OAAe,EACf,MAAc;QAJd,UAAK,GAAL,KAAK,CAAQ;QACb,WAAM,GAAN,MAAM,CAAQ;QACd,YAAO,GAAP,OAAO,CAAQ;QACf,YAAO,GAAP,OAAO,CAAQ;QACf,WAAM,GAAN,MAAM,CAAQ;QAErB,sBAAsB;IACxB,CAAC;IAED;;;;;;OAMG;IACH,MAAM,CAAC,eAAe,CAAC,KAAa,EAAE,MAAc;QAClD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC;QAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QAC5C,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC7E,MAAM,OAAO,GAAG,MAAM,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3E,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;QAE7E,OAAO,IAAI,oBAAoB,CAAC,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;IAC3E,CAAC;IAED,2EAA2E;IAC3E,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IACzD,CAAC;IAED;;;;;OAKG;IACH,WAAW,CAAC,KAAe,EAAE,YAAY,GAAG,GAAG;QAC7C,MAAM,EAAE,GAAG,IAAI,MAAM,CAAC,KAAK,EAAE,GAAG,EAAE,IAAI,CAAC,OAAO,GAAG,MAAM,EAAE,GAAG,EAAE,YAAY,CAAC,CAAC;QAC5E,MAAM,EAAE,GAAG,IAAI,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,EAAE,GAAG,EAAE,IAAI,CAAC,OAAO,GAAG,MAAM,EAAE,GAAG,EAAE,YAAY,CAAC,CAAC;QAE9F,OAAO,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAClB,CAAC;CACF"}

package/lib/maneuver/Waypoint.d.ts

@@ -0,0 +1,66 @@
+/**
+ * @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, Vector3D } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+import { StateInterpolator } from '../interpolator/StateInterpolator';
+import { ForceModel } from './../force/ForceModel';
+export declare class Waypoint {
+    epoch: EpochUTC;
+    relativePosition: Vector3D;
+    constructor(epoch: EpochUTC, relativePosition: Vector3D);
+    /**
+     * Return the perturbed error in a [maneuver] when compared against the
+     * target [waypoint] given an initial [state], [forceModel],
+     * [target] interpolator, and speculative relative maneuver
+     * [components] _(m/s)_.
+     */
+    static _error(waypoint: Waypoint, maneuver: Thrust, state: J2000, forceModel: ForceModel, target: StateInterpolator, components: Float64Array): number;
+    /**
+     * Generate a score function for refining perturbed [waypoint] maneuvers.
+     *
+     * The score function takes an array of speculative radial, intrack, and
+     * crosstrack components _(m/s)_ and returns the propagated error from the
+     * desired waypoint target.
+     */
+    static _refineManeuverScore(waypoint: Waypoint, maneuver: Thrust, state: J2000, forceModel: ForceModel, target: StateInterpolator): (components: Float64Array) => number;
+    /**
+     * Convert an array of [waypoints] into a maneuver sequence given an
+     * [interceptor] state, [pivot] epoch for the first burn to arrive at the
+     * first waypoint, and [target] ephemeris interpolator.
+     *
+     * Optional arguments are as follows:
+     *  - `preManeuvers`: maneuvers to execute before the pivot burn
+     * - `postManeuvers`: maneuvers to execute after the last pivot burn
+     * - `durationRate`: thruster duration rate _(s/m/s)_
+     * - `forceModel`: interceptor force model, defaults to two-body
+     * - `refine`: refine maneuvers to account for perturbations if `true`
+     * - `maxIter`: maximum refinement iterations per maneuver
+     * - `printIter`: print debug information on each refinement iteration
+     */
+    static toManeuvers(interceptor: J2000, pivot: EpochUTC, waypoints: Waypoint[], target: StateInterpolator, preManeuvers: Thrust[] | null, postManeuvers: Thrust[] | null, { durationRate, forceModel, refine, maxIter, printIter, }?: {
+        durationRate?: number;
+        forceModel?: ForceModel;
+        refine?: boolean;
+        maxIter?: number;
+        printIter?: boolean;
+    }): Thrust[];
+    static _refineManeuvers(waypoints: Waypoint[], maneuvers: Thrust[], interceptor: J2000, forceModel: ForceModel, target: StateInterpolator, { maxIter, printIter, }?: {
+        maxIter?: number;
+        printIter?: boolean;
+    }): Thrust[];
+}