ootk 3.0.0 → 4.0.0-0

package/lib/propagator/RungeKutta4Propagator.d.ts

@@ -0,0 +1,46 @@
+/**
+ * @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 { ForceModel } from '../force/ForceModel';
+import { Thrust } from '../force/Thrust';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+export declare class RungeKutta4Propagator extends Propagator {
+    private _initState;
+    private _forceModel;
+    private _stepSize;
+    private _cacheState;
+    private _checkpoints;
+    /**
+     * Create a new [RungeKutta4Propagator] object from an initial state vector and
+     * along with an optional [ForceModel] and [stepSize] in seconds.
+     */
+    constructor(_initState: J2000, _forceModel?: ForceModel, _stepSize?: number, _cacheState?: J2000, _checkpoints?: J2000[]);
+    setStepSize(seconds: number): void;
+    setForceModel(forceModel: ForceModel): void;
+    ephemerisManeuver(start: EpochUTC, finish: EpochUTC, maneuvers: Thrust[], interval?: number): VerletBlendInterpolator;
+    maneuver(maneuver: Thrust, interval?: number): J2000[];
+    private _kFn;
+    private _integrate;
+    propagate(epoch: EpochUTC): J2000;
+    reset(): void;
+    get state(): J2000;
+    checkpoint(): number;
+    clearCheckpoints(): void;
+    restore(index: number): void;
+}

package/lib/propagator/RungeKutta4Propagator.js

@@ -0,0 +1,132 @@
+/**
+ * @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 { J2000, Vector } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+// / Runge-Kutta 4 fixed numerical propagator.
+export class RungeKutta4Propagator extends Propagator {
+    /**
+     * Create a new [RungeKutta4Propagator] object from an initial state vector and
+     * along with an optional [ForceModel] and [stepSize] in seconds.
+     */
+    constructor(_initState, _forceModel = new ForceModel().setGravity(), _stepSize = 15.0, _cacheState = _initState, _checkpoints = []) {
+        super();
+        this._initState = _initState;
+        this._forceModel = _forceModel;
+        this._stepSize = _stepSize;
+        this._cacheState = _cacheState;
+        this._checkpoints = _checkpoints;
+        this._stepSize = Math.abs(_stepSize);
+    }
+    // / Set the integrator step size to the provided number of [seconds].
+    setStepSize(seconds) {
+        this._stepSize = Math.abs(seconds);
+    }
+    // / Set numerical integration force model.
+    setForceModel(forceModel) {
+        this._forceModel = forceModel;
+    }
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        const tMvr = maneuvers.slice(0).filter((mvr) => mvr.start >= start || mvr.stop <= finish);
+        const ephemeris = [];
+        if (tMvr[0].start > start) {
+            ephemeris.push(this.propagate(start));
+        }
+        for (const mvr of tMvr) {
+            while (this._cacheState.epoch < mvr.start) {
+                const step = Math.min(mvr.start.difference(this._cacheState.epoch), interval);
+                this.propagate(this._cacheState.epoch.roll(step));
+                if (this._cacheState.epoch.posix !== mvr.start.posix) {
+                    ephemeris.push(this._cacheState);
+                }
+            }
+            ephemeris.push(...this.maneuver(mvr, interval));
+        }
+        while (this._cacheState.epoch.posix < finish.posix) {
+            const step = Math.min(finish.difference(this._cacheState.epoch), interval);
+            this.propagate(this._cacheState.epoch.roll(step));
+            ephemeris.push(this._cacheState);
+        }
+        return new VerletBlendInterpolator(ephemeris);
+    }
+    maneuver(maneuver, interval = 60.0) {
+        if (maneuver.isImpulsive) {
+            this._cacheState = maneuver.apply(this.propagate(maneuver.center));
+            return [this._cacheState];
+        }
+        let tState = this.propagate(maneuver.start);
+        this._forceModel.loadManeuver(maneuver);
+        const ephemeris = [tState];
+        while (tState.epoch < maneuver.stop) {
+            const step = Math.min(maneuver.stop.difference(tState.epoch), interval);
+            tState = this.propagate(tState.epoch.roll(step));
+            ephemeris.push(tState);
+        }
+        this._forceModel.clearManeuver();
+        return ephemeris;
+    }
+    _kFn(state, hArg, kArg) {
+        const epoch = state.epoch.roll(hArg);
+        const posvel = state.position.join(state.velocity);
+        const result = posvel.add(kArg);
+        const sample = new J2000(epoch, result.toVector3D(0), result.toVector3D(3));
+        return this._forceModel.derivative(sample);
+    }
+    _integrate(state, step) {
+        const k1 = this._kFn(state, 0, Vector.zero(6)).scale(step);
+        const k2 = this._kFn(state, 0.5 * step, k1.scale(0.5)).scale(step);
+        const k3 = this._kFn(state, 0.5 * step, k2.scale(0.5)).scale(step);
+        const k4 = this._kFn(state, step, k3).scale(step);
+        const v1 = k1;
+        const v2 = v1.add(k2.scale(2));
+        const v3 = v2.add(k3.scale(2));
+        const v4 = v3.add(k4);
+        const tNext = state.epoch.roll(step);
+        const posvel = state.position.join(state.velocity);
+        const result = posvel.add(v4.scale(1 / 6));
+        return new J2000(tNext, result.toVector3D(0), result.toVector3D(3));
+    }
+    propagate(epoch) {
+        let delta = epoch.difference(this._cacheState.epoch);
+        while (delta !== 0) {
+            const direction = delta >= 0 ? 1 : -1;
+            const dt = Math.min(Math.abs(delta), this._stepSize) * direction;
+            this._cacheState = this._integrate(this._cacheState, dt);
+            delta = epoch.difference(this._cacheState.epoch);
+        }
+        return this._cacheState;
+    }
+    reset() {
+        this._cacheState = this._initState;
+    }
+    get state() {
+        return this._cacheState;
+    }
+    checkpoint() {
+        this._checkpoints.push(this._cacheState);
+        return this._checkpoints.length - 1;
+    }
+    clearCheckpoints() {
+        this._checkpoints = [];
+    }
+    restore(index) {
+        this._cacheState = this._checkpoints[index];
+    }
+}
+//# sourceMappingURL=RungeKutta4Propagator.js.map

package/lib/propagator/RungeKutta4Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKutta4Propagator.js","sourceRoot":"","sources":["../../src/propagator/RungeKutta4Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAY,KAAK,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AACpD,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AAEjD,OAAO,EAAE,uBAAuB,EAAE,MAAM,yCAAyC,CAAC;AAClF,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C,8CAA8C;AAC9C,MAAM,OAAO,qBAAsB,SAAQ,UAAU;IACnD;;;OAGG;IACH,YACU,UAAiB,EACjB,cAA0B,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,EACvD,YAAoB,IAAI,EACxB,cAAqB,UAAU,EAC/B,eAAwB,EAAE;QAElC,KAAK,EAAE,CAAC;QANA,eAAU,GAAV,UAAU,CAAO;QACjB,gBAAW,GAAX,WAAW,CAA4C;QACvD,cAAS,GAAT,SAAS,CAAe;QACxB,gBAAW,GAAX,WAAW,CAAoB;QAC/B,iBAAY,GAAZ,YAAY,CAAc;QAGlC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;IACvC,CAAC;IAED,sEAAsE;IACtE,WAAW,CAAC,OAAe;QACzB,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC;IACrC,CAAC;IAED,2CAA2C;IAC3C,aAAa,CAAC,UAAsB;QAClC,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;IAChC,CAAC;IAED,iBAAiB,CAAC,KAAe,EAAE,MAAgB,EAAE,SAAmB,EAAE,QAAQ,GAAG,IAAI;QACvF,MAAM,IAAI,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,IAAI,KAAK,IAAI,GAAG,CAAC,IAAI,IAAI,MAAM,CAAC,CAAC;QAC1F,MAAM,SAAS,GAAY,EAAE,CAAC;QAE9B,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,EAAE;YACzB,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;SACvC;QACD,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;YACtB,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,GAAG,GAAG,CAAC,KAAK,EAAE;gBACzC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAE9E,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;gBAClD,IAAI,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,KAAK,KAAK,GAAG,CAAC,KAAK,CAAC,KAAK,EAAE;oBACpD,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;iBAClC;aACF;YACD,SAAS,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC,CAAC;SACjD;QACD,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,EAAE;YAClD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;YAE3E,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;YAClD,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;SAClC;QAED,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IAED,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,IAAI;QACxC,IAAI,QAAQ,CAAC,WAAW,EAAE;YACxB,IAAI,CAAC,WAAW,GAAG,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;YAEnE,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;SAC3B;QACD,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QAE5C,IAAI,CAAC,WAAW,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACxC,MAAM,SAAS,GAAY,CAAC,MAAM,CAAC,CAAC;QAEpC,OAAO,MAAM,CAAC,KAAK,GAAG,QAAQ,CAAC,IAAI,EAAE;YACnC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,UAAU,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;YAExE,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;YACjD,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;SACxB;QACD,IAAI,CAAC,WAAW,CAAC,aAAa,EAAE,CAAC;QAEjC,OAAO,SAAS,CAAC;IACnB,CAAC;IAEO,IAAI,CAAC,KAAY,EAAE,IAAY,EAAE,IAAY;QACnD,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACrC,MAAM,MAAM,GAAG,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACnD,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAChC,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAE5E,OAAO,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IAEO,UAAU,CAAC,KAAY,EAAE,IAAY;QAC3C,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC3D,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,GAAG,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACnE,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,GAAG,IAAI,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACnE,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAClD,MAAM,EAAE,GAAG,EAAE,CAAC;QACd,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QACtB,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QACrC,MAAM,MAAM,GAAG,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QACnD,MAAM,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAE3C,OAAO,IAAI,KAAK,CAAC,KAAK,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;IACtE,CAAC;IAED,SAAS,CAAC,KAAe;QACvB,IAAI,KAAK,GAAG,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;QAErD,OAAO,KAAK,KAAK,CAAC,EAAE;YAClB,MAAM,SAAS,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,SAAS,CAAC;YAEjE,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC;YACzD,KAAK,GAAG,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;SAClD;QAED,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,KAAK;QACH,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC;IACrC,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,UAAU;QACR,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;QAEzC,OAAO,IAAI,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;IACtC,CAAC;IAED,gBAAgB;QACd,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC;IACzB,CAAC;IAED,OAAO,CAAC,KAAa;QACnB,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC;IAC9C,CAAC;CACF"}

package/lib/propagator/RungeKutta89Propagator.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 { RungeKuttaAdaptive } from './RungeKuttaAdaptive';
+export declare class RungeKutta89Propagator extends RungeKuttaAdaptive {
+    private a_;
+    private b_;
+    private ch_;
+    private c_;
+    get a(): Float64Array;
+    get b(): Float64Array[];
+    get ch(): Float64Array;
+    get c(): Float64Array;
+    get order(): number;
+}

package/lib/propagator/RungeKutta89Propagator.js

@@ -0,0 +1,130 @@
+/**
+ * @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-loss-of-precision */
+/* eslint-disable class-methods-use-this */
+import { RungeKuttaAdaptive } from './RungeKuttaAdaptive';
+// / Runge-Kutta 8(9) adaptive numerical propagator.
+export class RungeKutta89Propagator extends RungeKuttaAdaptive {
+    constructor() {
+        super(...arguments);
+        this.a_ = Float64Array.from([
+            0.0, 0.3462e-1, 0.9702435063878044594828361677100617517633e-1, 0.1455365259581706689224254251565092627645, 0.561,
+            0.2290079115904850126662751771814700052182, 0.5449920884095149873337248228185299947818, 0.645, 0.48375, 0.6757e-1,
+            0.25, 0.6590650618730998549405331618649220295334, 0.8206, 0.9012, 1.0, 1.0,
+        ]);
+        this.b_ = [
+            new Float64Array(0),
+            Float64Array.from([0.3462e-1]),
+            Float64Array.from([-0.389335438857287327017042687229284478532e-1, 0.1359578945245091786499878854939346230295]),
+            Float64Array.from([0.3638413148954266723060635628912731569111e-1, 0.0, 0.1091523944686280016918190688673819470733]),
+            Float64Array.from([
+                2.025763914393969636805657604282571047511, 0.0, -7.638023836496292020387602153091964592952,
+                6.173259922102322383581944548809393545442,
+            ]),
+            Float64Array.from([
+                0.5112275589406060872792270881648288397197e-1, 0.0, 0.0, 0.1770823794555021537929910813839068684087,
+                0.80277624092225014536138698108025283759e-3,
+            ]),
+            Float64Array.from([
+                0.1316006357975216279279871693164256985334, 0.0, 0.0, -0.2957276252669636417685183174672273730699,
+                0.878137803564295237421124704053886667082e-1, 0.62130529752252747743214350056394300261,
+            ]),
+            Float64Array.from([
+                0.7166666666666666666666666666666666666667e-1, 0.0, 0.0, 0.0, 0.0, 0.3305533578915319409260346730051472207728,
+                0.2427799754418013924072986603281861125606,
+            ]),
+            Float64Array.from([
+                0.71806640625e-1, 0.0, 0.0, 0.0, 0.0, 0.3294380283228177160744825466257672816401,
+                0.1165190029271822839255174533742327183599, -0.34013671875e-1,
+            ]),
+            Float64Array.from([
+                0.4836757646340646986611287718844085773549e-1, 0.0, 0.0, 0.0, 0.0, 0.3928989925676163974333190042057047002852e-1,
+                0.1054740945890344608263649267140088017604, -0.2143865284648312665982642293830533996214e-1,
+                -0.1041229174627194437759832813847147895623,
+            ]),
+            Float64Array.from([
+                -0.2664561487201478635337289243849737340534e-1, 0.0, 0.0, 0.0, 0.0, 0.3333333333333333333333333333333333333333e-1,
+                -0.1631072244872467239162704487554706387141, 0.3396081684127761199487954930015522928244e-1,
+                0.1572319413814626097110769806810024118077, 0.215226747803187955230353477879477037696,
+            ]),
+            Float64Array.from([
+                0.3689009248708622334786359863227633989718e-1, 0.0, 0.0, 0.0, 0.0, -0.1465181576725542928653609891758501156785,
+                0.2242577768172024345345469822625833796001, 0.2294405717066072637090897902753790803034e-1,
+                -0.35850052905728761357394424889330334334e-2, 0.8669223316444385506869203619044453906053e-1,
+                0.4383840651968337846196219974168630120572,
+            ]),
+            Float64Array.from([
+                -0.4866012215113340846662212357570395295088, 0.0, 0.0, 0.0, 0.0, -6.304602650282852990657772792012007122988,
+                -0.281245618289472564778284183790118418111, -2.679019236219849057687906597489223155566,
+                0.518815663924157511565311164615012522024, 1.365353187603341710683633635235238678626,
+                5.885091088503946585721274891680604830712, 2.802808786272062889819965117517532194812,
+            ]),
+            Float64Array.from([
+                0.4185367457753471441471025246471931649633, 0.0, 0.0, 0.0, 0.0, 6.724547581906459363100870806514855026676,
+                -0.425444280164611790606983409697113064616, 3.343279153001265577811816947557982637749,
+                0.617081663117537759528421117507709784737, -0.929966123939932833937749523988800852013,
+                -6.099948804751010722472962837945508844846, -3.002206187889399044804158084895173690015,
+                0.2553202529443445472336424602988558373637,
+            ]),
+            Float64Array.from([
+                -0.779374086122884664644623040843840506343, 0.0, 0.0, 0.0, 0.0, -13.93734253810777678786523664804936051203,
+                1.252048853379357320949735183924200895136, -14.69150040801686878191527989293072091588,
+                -0.494705058533141685655191992136962873577, 2.242974909146236657906984549543692874755,
+                13.36789380382864375813864978592679139881, 14.39665048665068644512236935340272139005,
+                -0.7975813331776800379127866056663258667437, 0.4409353709534277758753793068298041158235,
+            ]),
+            Float64Array.from([
+                2.058051337466886442151242368989994043993, 0.0, 0.0, 0.0, 0.0, 22.35793772796803295519317565842520212899,
+                0.90949810997556332745009198137971890783, 35.89110098240264104710550686568482456493,
+                -3.442515027624453437985000403608480262211, -4.865481358036368826566013387928704014496,
+                -18.90980381354342625688427480879773032857, -34.26354448030451782929251177395134170515,
+                1.264756521695642578827783499806516664686, 0.0, 0.0,
+            ]),
+        ];
+        this.ch_ = Float64Array.from([
+            0.1996996514886773085518508418098868756464e-1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+            2.191499304949330054530747099310837524864, 0.8857071848208438030833722031786358862953e-1,
+            0.1140560234865965622484956605091432032674, 0.2533163805345107065564577734569651977347,
+            -2.056564386240941011158999594595981300493, 0.340809679901311993516009489422454381283, 0.0, 0.0,
+            0.4834231373823958314376726739772871714902e-1,
+        ]);
+        this.c_ = Float64Array.from([
+            0.1461197685842315252051541915018784713459e-1, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0,
+            -0.391521186233133908941022826728824203081, 0.2310932500289506415909675644868993669908,
+            0.1274766769992852382560589467488989175618, 0.2246434176204157731566981937082069688984,
+            0.5684352689748512932705226972873692126743, 0.5825871557215827200814768021863420902155e-1,
+            0.1364317403482215641609022744494239843327, 0.3057013983082797397721005067920369646664e-1, 0.0,
+        ]);
+    }
+    get a() {
+        return this.a_;
+    }
+    get b() {
+        return this.b_;
+    }
+    get ch() {
+        return this.ch_;
+    }
+    get c() {
+        return this.c_;
+    }
+    get order() {
+        return 8;
+    }
+}
+//# sourceMappingURL=RungeKutta89Propagator.js.map

package/lib/propagator/RungeKutta89Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKutta89Propagator.js","sourceRoot":"","sources":["../../src/propagator/RungeKutta89Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,4DAA4D;AAC5D,2CAA2C;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D,oDAAoD;AACpD,MAAM,OAAO,sBAAuB,SAAQ,kBAAkB;IAA9D;;QACU,OAAE,GAAG,YAAY,CAAC,IAAI,CAAC;YAC7B,GAAG,EAAE,SAAS,EAAE,6CAA6C,EAAE,0CAA0C,EAAE,KAAK;YAChH,0CAA0C,EAAE,0CAA0C,EAAE,KAAK,EAAE,OAAO,EAAE,SAAS;YACjH,IAAI,EAAE,0CAA0C,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,GAAG;SAC3E,CAAC,CAAC;QACK,OAAE,GAAG;YACX,IAAI,YAAY,CAAC,CAAC,CAAC;YACnB,YAAY,CAAC,IAAI,CAAC,CAAC,SAAS,CAAC,CAAC;YAC9B,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,4CAA4C,EAAE,0CAA0C,CAAC,CAAC;YAC9G,YAAY,CAAC,IAAI,CAAC,CAAC,6CAA6C,EAAE,GAAG,EAAE,0CAA0C,CAAC,CAAC;YACnH,YAAY,CAAC,IAAI,CAAC;gBAChB,yCAAyC,EAAE,GAAG,EAAE,CAAC,yCAAyC;gBAC1F,yCAAyC;aAC1C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;gBACnG,2CAA2C;aAC5C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,0CAA0C;gBACjG,4CAA4C,EAAE,wCAAwC;aACvF,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;gBAC7G,0CAA0C;aAC3C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,gBAAgB,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;gBAChF,0CAA0C,EAAE,CAAC,gBAAgB;aAC9D,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,6CAA6C;gBAChH,0CAA0C,EAAE,CAAC,6CAA6C;gBAC1F,CAAC,0CAA0C;aAC5C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,CAAC,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,6CAA6C;gBACjH,CAAC,0CAA0C,EAAE,6CAA6C;gBAC1F,0CAA0C,EAAE,yCAAyC;aACtF,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,0CAA0C;gBAC9G,0CAA0C,EAAE,6CAA6C;gBACzF,CAAC,2CAA2C,EAAE,6CAA6C;gBAC3F,0CAA0C;aAC3C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,CAAC,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,yCAAyC;gBAC3G,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;gBACtF,yCAAyC,EAAE,yCAAyC;gBACpF,yCAAyC,EAAE,yCAAyC;aACrF,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,yCAAyC;gBACzG,CAAC,yCAAyC,EAAE,yCAAyC;gBACrF,yCAAyC,EAAE,CAAC,yCAAyC;gBACrF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;gBACtF,0CAA0C;aAC3C,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,CAAC,yCAAyC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,yCAAyC;gBAC1G,yCAAyC,EAAE,CAAC,yCAAyC;gBACrF,CAAC,yCAAyC,EAAE,yCAAyC;gBACrF,yCAAyC,EAAE,yCAAyC;gBACpF,CAAC,0CAA0C,EAAE,0CAA0C;aACxF,CAAC;YACF,YAAY,CAAC,IAAI,CAAC;gBAChB,yCAAyC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,yCAAyC;gBACxG,wCAAwC,EAAE,yCAAyC;gBACnF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;gBACtF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;gBACtF,yCAAyC,EAAE,GAAG,EAAE,GAAG;aACpD,CAAC;SACH,CAAC;QACM,QAAG,GAAG,YAAY,CAAC,IAAI,CAAC;YAC9B,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG;YAC3E,yCAAyC,EAAE,6CAA6C;YACxF,0CAA0C,EAAE,0CAA0C;YACtF,CAAC,yCAAyC,EAAE,yCAAyC,EAAE,GAAG,EAAE,GAAG;YAC/F,6CAA6C;SAC9C,CAAC,CAAC;QACK,OAAE,GAAG,YAAY,CAAC,IAAI,CAAC;YAC7B,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG;YAC3E,CAAC,yCAAyC,EAAE,0CAA0C;YACtF,0CAA0C,EAAE,0CAA0C;YACtF,0CAA0C,EAAE,6CAA6C;YACzF,0CAA0C,EAAE,6CAA6C,EAAE,GAAG;SAC/F,CAAC,CAAC;IAqBL,CAAC;IAnBC,IAAI,CAAC;QACH,OAAO,IAAI,CAAC,EAAE,CAAC;IACjB,CAAC;IAED,IAAI,CAAC;QACH,OAAO,IAAI,CAAC,EAAE,CAAC;IACjB,CAAC;IAED,IAAI,EAAE;QACJ,OAAO,IAAI,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,IAAI,CAAC;QACH,OAAO,IAAI,CAAC,EAAE,CAAC;IACjB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,CAAC,CAAC;IACX,CAAC;CACF"}

package/lib/propagator/RungeKuttaAdaptive.d.ts

@@ -0,0 +1,52 @@
+/**
+ * @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 { ForceModel } from '../force/ForceModel';
+import { Thrust } from '../force/Thrust';
+import { VerletBlendInterpolator } from './../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+export declare abstract class RungeKuttaAdaptive extends Propagator {
+    private _initState;
+    private _forceModel;
+    private _tolerance;
+    /**
+     * Create a new [RungeKuttaAdaptive] object from an initial state vector
+     * along with an optional [ForceModel] and [tolerance].
+     */
+    constructor(_initState: J2000, _forceModel?: ForceModel, _tolerance?: number);
+    private _cacheState;
+    private readonly _checkpoints;
+    private _stepSize;
+    private static readonly _minTolerance;
+    protected abstract get a(): Float64Array;
+    protected abstract get b(): Float64Array[];
+    protected abstract get ch(): Float64Array;
+    protected abstract get c(): Float64Array;
+    protected abstract get order(): number;
+    get state(): J2000;
+    reset(): void;
+    setForceModel(forceModel: ForceModel): void;
+    private _kfn;
+    private _integrate;
+    propagate(epoch: EpochUTC): J2000;
+    maneuver(maneuver: Thrust, interval?: number): J2000[];
+    ephemerisManeuver(start: EpochUTC, finish: EpochUTC, maneuvers: Thrust[], interval?: number): VerletBlendInterpolator;
+    checkpoint(): number;
+    clearCheckpoints(): void;
+    restore(index: number): void;
+}

package/lib/propagator/RungeKuttaAdaptive.js

@@ -0,0 +1,152 @@
+/**
+ * @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 { J2000, Vector } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { VerletBlendInterpolator } from './../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+import { RkCheckpoint } from './RkCheckpoint';
+import { RkResult } from './RkResult';
+// / Adaptive Runge-Kutta propagator base class.
+export class RungeKuttaAdaptive extends Propagator {
+    /**
+     * Create a new [RungeKuttaAdaptive] object from an initial state vector
+     * along with an optional [ForceModel] and [tolerance].
+     */
+    constructor(_initState, _forceModel = new ForceModel().setGravity(), _tolerance = 1e-9) {
+        super();
+        this._initState = _initState;
+        this._forceModel = _forceModel;
+        this._tolerance = _tolerance;
+        this._checkpoints = [];
+        // / Integration step size _(seconds)_.
+        this._stepSize = 60.0;
+        this._cacheState = this._initState;
+        this._tolerance = Math.max(RungeKuttaAdaptive._minTolerance, Math.abs(_tolerance));
+    }
+    get state() {
+        return this._cacheState;
+    }
+    reset() {
+        this._cacheState = this._initState;
+        this._stepSize = 60.0;
+    }
+    // / Set numerical integration force model.
+    setForceModel(forceModel) {
+        this._forceModel = forceModel;
+    }
+    _kfn(epoch, rv, hArg, kArg, step) {
+        const t = epoch.roll(hArg * step);
+        const rvNew = rv.add(kArg);
+        const sample = new J2000(t, rvNew.toVector3D(0), rvNew.toVector3D(3));
+        return this._forceModel.derivative(sample).scale(step);
+    }
+    _integrate(state, step) {
+        const k = new Array(this.a.length).fill(Vector.origin3);
+        const y = state.position.join(state.velocity);
+        for (let i = 0; i < this.a.length; i++) {
+            let kArg = Vector.origin6;
+            if (i !== 0) {
+                for (let j = 0; j < i; j++) {
+                    kArg = kArg.add(k[j].scale(this.b[i][j]));
+                }
+            }
+            k[i] = this._kfn(state.epoch, y, this.a[i], kArg, step);
+        }
+        let y1 = y;
+        let y2 = y;
+        for (let i = 0; i < k.length; i++) {
+            y1 = y1.add(k[i].scale(this.ch[i]));
+            y2 = y2.add(k[i].scale(this.c[i]));
+        }
+        const teVal = y1.distance(y2);
+        let hNew = 0.9 * step * (this._tolerance / teVal) ** (1.0 / this.order);
+        const hOld = Math.abs(step);
+        hNew = Math.max(0.2 * hOld, Math.min(5.0 * hOld, hNew));
+        hNew = Math.max(1e-5, Math.min(1000.0, hNew));
+        return new RkResult(new J2000(state.epoch.roll(step), y1.toVector3D(0), y1.toVector3D(3)), teVal, hNew);
+    }
+    propagate(epoch) {
+        let delta = epoch.difference(this._cacheState.epoch);
+        while (delta !== 0) {
+            const direction = delta >= 0 ? 1 : -1;
+            const dt = Math.min(Math.abs(delta), this._stepSize) * direction;
+            const result = this._integrate(this._cacheState, dt);
+            this._stepSize = result.newStep;
+            if (result.error > this._tolerance) {
+                continue;
+            }
+            this._cacheState = result.state;
+            delta = epoch.difference(this._cacheState.epoch);
+        }
+        return this._cacheState;
+    }
+    maneuver(maneuver, interval = 60.0) {
+        if (maneuver.isImpulsive) {
+            this._cacheState = maneuver.apply(this.propagate(maneuver.center));
+            return [this._cacheState];
+        }
+        let tState = this.propagate(maneuver.start);
+        this._forceModel.loadManeuver(maneuver);
+        const ephemeris = [tState];
+        while (tState.epoch < maneuver.stop) {
+            const step = Math.min(maneuver.stop.difference(tState.epoch), interval);
+            tState = this.propagate(tState.epoch.roll(step));
+            ephemeris.push(tState);
+        }
+        this._forceModel.clearManeuver();
+        return ephemeris;
+    }
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        const tMvr = maneuvers.filter((mvr) => mvr.start >= start || mvr.stop <= finish);
+        const ephemeris = [];
+        if (tMvr[0].start > start) {
+            ephemeris.push(this.propagate(start));
+        }
+        for (const mvr of tMvr) {
+            while (this._cacheState.epoch < mvr.start) {
+                const step = Math.min(mvr.start.difference(this._cacheState.epoch), interval);
+                this.propagate(this._cacheState.epoch.roll(step));
+                if (this._cacheState.epoch.posix !== mvr.start.posix) {
+                    ephemeris.push(this._cacheState);
+                }
+            }
+            ephemeris.push(...this.maneuver(mvr, interval));
+        }
+        while (this._cacheState.epoch.posix < finish.posix) {
+            const step = Math.min(finish.difference(this._cacheState.epoch), interval);
+            this.propagate(this._cacheState.epoch.roll(step));
+            ephemeris.push(this._cacheState);
+        }
+        return new VerletBlendInterpolator(ephemeris);
+    }
+    checkpoint() {
+        this._checkpoints.push(new RkCheckpoint(this._cacheState, this._stepSize));
+        return this._checkpoints.length - 1;
+    }
+    clearCheckpoints() {
+        this._checkpoints.length = 0;
+    }
+    restore(index) {
+        const checkpoint = this._checkpoints[index];
+        this._cacheState = checkpoint.cacheState;
+        this._stepSize = checkpoint.stepSize;
+    }
+}
+// / Minimum allowable local error tolerance.
+RungeKuttaAdaptive._minTolerance = 1e-15;
+//# sourceMappingURL=RungeKuttaAdaptive.js.map

package/lib/propagator/RungeKuttaAdaptive.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKuttaAdaptive.js","sourceRoot":"","sources":["../../src/propagator/RungeKuttaAdaptive.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAY,KAAK,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AACpD,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AAEjD,OAAO,EAAE,uBAAuB,EAAE,MAAM,2CAA2C,CAAC;AACpF,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,YAAY,EAAE,MAAM,gBAAgB,CAAC;AAC9C,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAEtC,gDAAgD;AAChD,MAAM,OAAgB,kBAAmB,SAAQ,UAAU;IACzD;;;OAGG;IACH,YACU,UAAiB,EACjB,cAA0B,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,EACvD,aAAqB,IAAI;QAEjC,KAAK,EAAE,CAAC;QAJA,eAAU,GAAV,UAAU,CAAO;QACjB,gBAAW,GAAX,WAAW,CAA4C;QACvD,eAAU,GAAV,UAAU,CAAe;QAUlB,iBAAY,GAAmB,EAAE,CAAC;QAEnD,uCAAuC;QAC/B,cAAS,GAAG,IAAI,CAAC;QAVvB,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC;QACnC,IAAI,CAAC,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,kBAAkB,CAAC,aAAa,EAAE,IAAI,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC;IACrF,CAAC;IA4BD,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,KAAK;QACH,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,UAAU,CAAC;QACnC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC;IACxB,CAAC;IAED,2CAA2C;IAC3C,aAAa,CAAC,UAAsB;QAClC,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC;IAChC,CAAC;IAEO,IAAI,CAAC,KAAe,EAAE,EAAU,EAAE,IAAY,EAAE,IAAY,EAAE,IAAY;QAChF,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;QAClC,MAAM,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,CAAC,EAAE,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC;QAEtE,OAAO,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACzD,CAAC;IAEO,UAAU,CAAC,KAAY,EAAE,IAAY;QAC3C,MAAM,CAAC,GAAa,IAAI,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC,CAAC;QAClE,MAAM,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QAE9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACtC,IAAI,IAAI,GAAG,MAAM,CAAC,OAAyB,CAAC;YAE5C,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;oBAC1B,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;iBAC3C;aACF;YACD,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;SACzD;QACD,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,EAAE,GAAG,CAAC,CAAC;QAEX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACjC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SACpC;QACD,MAAM,KAAK,GAAG,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAC9B,IAAI,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,CAAC,IAAI,CAAC,UAAU,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;QACxE,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAE5B,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,EAAE,IAAI,CAAC,CAAC,CAAC;QACxD,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC;QAE9C,OAAO,IAAI,QAAQ,CAAC,IAAI,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC;IAC1G,CAAC;IAED,SAAS,CAAC,KAAe;QACvB,IAAI,KAAK,GAAG,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;QAErD,OAAO,KAAK,KAAK,CAAC,EAAE;YAClB,MAAM,SAAS,GAAG,KAAK,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,SAAS,CAAC;YACjE,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,EAAE,EAAE,CAAC,CAAC;YAErD,IAAI,CAAC,SAAS,GAAG,MAAM,CAAC,OAAO,CAAC;YAChC,IAAI,MAAM,CAAC,KAAK,GAAG,IAAI,CAAC,UAAU,EAAE;gBAClC,SAAS;aACV;YACD,IAAI,CAAC,WAAW,GAAG,MAAM,CAAC,KAAK,CAAC;YAChC,KAAK,GAAG,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;SAClD;QAED,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,IAAI;QACxC,IAAI,QAAQ,CAAC,WAAW,EAAE;YACxB,IAAI,CAAC,WAAW,GAAG,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;YAEnE,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;SAC3B;QACD,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QAE5C,IAAI,CAAC,WAAW,CAAC,YAAY,CAAC,QAAQ,CAAC,CAAC;QACxC,MAAM,SAAS,GAAY,CAAC,MAAM,CAAC,CAAC;QAEpC,OAAO,MAAM,CAAC,KAAK,GAAG,QAAQ,CAAC,IAAI,EAAE;YACnC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,UAAU,CAAC,MAAM,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;YAExE,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;YACjD,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;SACxB;QACD,IAAI,CAAC,WAAW,CAAC,aAAa,EAAE,CAAC;QAEjC,OAAO,SAAS,CAAC;IACnB,CAAC;IAED,iBAAiB,CAAC,KAAe,EAAE,MAAgB,EAAE,SAAmB,EAAE,QAAQ,GAAG,IAAI;QACvF,MAAM,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,KAAK,IAAI,KAAK,IAAI,GAAG,CAAC,IAAI,IAAI,MAAM,CAAC,CAAC;QACjF,MAAM,SAAS,GAAY,EAAE,CAAC;QAE9B,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,EAAE;YACzB,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;SACvC;QACD,KAAK,MAAM,GAAG,IAAI,IAAI,EAAE;YACtB,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,GAAG,GAAG,CAAC,KAAK,EAAE;gBACzC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAE9E,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;gBAClD,IAAI,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,KAAK,KAAK,GAAG,CAAC,KAAK,CAAC,KAAK,EAAE;oBACpD,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;iBAClC;aACF;YACD,SAAS,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,QAAQ,CAAC,CAAC,CAAC;SACjD;QACD,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,KAAK,EAAE;YAClD,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;YAE3E,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;YAClD,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;SAClC;QAED,OAAO,IAAI,uBAAuB,CAAC,SAAS,CAAC,CAAC;IAChD,CAAC;IAED,UAAU;QACR,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,IAAI,YAAY,CAAC,IAAI,CAAC,WAAW,EAAE,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC;QAE3E,OAAO,IAAI,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;IACtC,CAAC;IAED,gBAAgB;QACd,IAAI,CAAC,YAAY,CAAC,MAAM,GAAG,CAAC,CAAC;IAC/B,CAAC;IAED,OAAO,CAAC,KAAa;QACnB,MAAM,UAAU,GAAG,IAAI,CAAC,YAAY,CAAC,KAAK,CAAC,CAAC;QAE5C,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,UAAU,CAAC;QACzC,IAAI,CAAC,SAAS,GAAG,UAAU,CAAC,QAAQ,CAAC;IACvC,CAAC;;AA3JD,6CAA6C;AACrB,gCAAa,GAAW,KAAK,CAAC"}

package/lib/propagator/Sgp4Propagator.d.ts

@@ -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/>.
+ */
+import { EpochUTC, J2000, TLE } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+export declare class Sgp4Propagator extends Propagator {
+    private tle;
+    constructor(tle: TLE);
+    private _cacheState;
+    private _checkpoints;
+    get state(): J2000;
+    ephemerisManeuver(start: EpochUTC, finish: EpochUTC, maneuvers: Thrust[], interval?: number): VerletBlendInterpolator;
+    maneuver(maneuver: Thrust, interval?: number): J2000[];
+    propagate(epoch: EpochUTC): J2000;
+    reset(): void;
+    checkpoint(): number;
+    clearCheckpoints(): void;
+    restore(index: number): void;
+}