ootk 2.0.1 → 4.0.0-0

package/lib/orbit_determination/ModifiedGoodingIOD.js

@@ -0,0 +1,82 @@
+/**
+ * @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, J2000, RadecTopocentric, Vector3D } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator';
+import { DownhillSimplex } from './../optimize/DownhillSimplex';
+import { GoodingIOD } from './GoodingIOD';
+/**
+ * Gooding angles-only initial orbit determination.
+ *
+ * Used for orbit determination from multiple optical observations.
+ */
+export class ModifiedGoodingIOD {
+    constructor(observations, mu = Earth.mu) {
+        this._observations = observations;
+        this._mu = mu;
+    }
+    _createInitial(r0, rN, nRev, direction) {
+        const iod = new GoodingIOD(this._observations[0], this._observations[Math.floor(this._observations.length / 2)], this._observations[this._observations.length - 1], this._mu);
+        return iod.solve(r0, rN, nRev, direction);
+    }
+    _createErrorFunction(aprioriEpoch) {
+        const forceModel = new ForceModel().setGravity(this._mu);
+        const scoreFn = (x) => {
+            const position = new Vector3D(x[0], x[1], x[2]);
+            const velocity = new Vector3D(x[3], x[4], x[5]);
+            const state = new J2000(aprioriEpoch, position, velocity);
+            const propagator = new RungeKutta89Propagator(state, forceModel);
+            let total = 0;
+            for (const oC of this._observations) {
+                const sC = propagator.propagate(oC.epoch);
+                const pC = oC.site;
+                const expected = oC.observation.lineOfSight();
+                const actual = RadecTopocentric.fromStateVectors(sC, pC).lineOfSight();
+                const error = expected.angle(actual);
+                total += error;
+            }
+            return total;
+        };
+        return scoreFn;
+    }
+    solve(r0, rN, { nRev = 0, direction = true, posSearch = 10.0, velSearch = 0.1, tolerance = 1e-6, printIter = false, }) {
+        if (this._observations.length < 3) {
+            throw new Error('At least 3 observations required for Gooding IOD.');
+        }
+        const init = this._createInitial(r0, rN, nRev, direction);
+        const guess = Float64Array.from([...init.position.toArray(), ...init.velocity.toArray()]);
+        const solveFn = this._createErrorFunction(init.epoch);
+        const simplex = [
+            Float64Array.from(guess),
+            Float64Array.from([guess[0] + posSearch, guess[1], guess[2], guess[3], guess[4], guess[5]]),
+            Float64Array.from([guess[0], guess[1] + posSearch, guess[2], guess[3], guess[4], guess[5]]),
+            Float64Array.from([guess[0], guess[1], guess[2] + posSearch, guess[3], guess[4], guess[5]]),
+            Float64Array.from([guess[0], guess[1], guess[2], guess[3] + velSearch, guess[4], guess[5]]),
+            Float64Array.from([guess[0], guess[1], guess[2], guess[3], guess[4] + velSearch, guess[5]]),
+            Float64Array.from([guess[0], guess[1], guess[2], guess[3], guess[4], guess[5] + velSearch]),
+        ];
+        const result = DownhillSimplex.solveSimplex(solveFn, simplex, {
+            adaptive: true,
+            xTolerance: tolerance,
+            fTolerance: tolerance,
+            printIter,
+        });
+        return new J2000(init.epoch, new Vector3D(result[0], result[1], result[2]), new Vector3D(result[3], result[4], result[5]));
+    }
+}
+//# sourceMappingURL=ModifiedGoodingIOD.js.map

package/lib/orbit_determination/ModifiedGoodingIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ModifiedGoodingIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/ModifiedGoodingIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,KAAK,EAAY,KAAK,EAAE,gBAAgB,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAC/E,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AACjD,OAAO,EAAE,sBAAsB,EAAE,MAAM,sCAAsC,CAAC;AAE9E,OAAO,EAAE,eAAe,EAAE,MAAM,+BAA+B,CAAC;AAEhE,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAW1C;;;;GAIG;AACH,MAAM,OAAO,kBAAkB;IAI7B,YAAY,YAAkC,EAAE,KAAa,KAAK,CAAC,EAAE;QACnE,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;IAChB,CAAC;IAEO,cAAc,CAAC,EAAU,EAAE,EAAU,EAAE,IAAY,EAAE,SAAkB;QAC7E,MAAM,GAAG,GAAG,IAAI,UAAU,CACxB,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,EACrB,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,EAC7D,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,aAAa,CAAC,MAAM,GAAG,CAAC,CAAC,EACjD,IAAI,CAAC,GAAG,CACT,CAAC;QAEF,OAAO,GAAG,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,SAAS,CAAC,CAAC;IAC5C,CAAC;IAEO,oBAAoB,CAAC,YAAsB;QACjD,MAAM,UAAU,GAAG,IAAI,UAAU,EAAE,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACzD,MAAM,OAAO,GAAG,CAAC,CAAe,EAAU,EAAE;YAC1C,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,MAAM,QAAQ,GAAG,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,YAAY,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;YAC1D,MAAM,UAAU,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;YACjE,IAAI,KAAK,GAAG,CAAC,CAAC;YAEd,KAAK,MAAM,EAAE,IAAI,IAAI,CAAC,aAAa,EAAE;gBACnC,MAAM,EAAE,GAAG,UAAU,CAAC,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;gBAC1C,MAAM,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC;gBACnB,MAAM,QAAQ,GAAG,EAAE,CAAC,WAAW,CAAC,WAAW,EAAE,CAAC;gBAC9C,MAAM,MAAM,GAAG,gBAAgB,CAAC,gBAAgB,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,WAAW,EAAE,CAAC;gBACvE,MAAM,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;gBAErC,KAAK,IAAI,KAAK,CAAC;aAChB;YAED,OAAO,KAAK,CAAC;QACf,CAAC,CAAC;QAEF,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,KAAK,CACH,EAAU,EACV,EAAU,EACV,EACE,IAAI,GAAG,CAAC,EACR,SAAS,GAAG,IAAI,EAChB,SAAS,GAAG,IAAI,EAChB,SAAS,GAAG,GAAG,EACf,SAAS,GAAG,IAAI,EAChB,SAAS,GAAG,KAAK,GACJ;QAEf,IAAI,IAAI,CAAC,aAAa,CAAC,MAAM,GAAG,CAAC,EAAE;YACjC,MAAM,IAAI,KAAK,CAAC,mDAAmD,CAAC,CAAC;SACtE;QACD,MAAM,IAAI,GAAG,IAAI,CAAC,cAAc,CAAC,EAAE,EAAE,EAAE,EAAE,IAAI,EAAE,SAAS,CAAC,CAAC;QAC1D,MAAM,KAAK,GAAG,YAAY,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,EAAE,GAAG,IAAI,CAAC,QAAQ,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC;QAC1F,MAAM,OAAO,GAAG,IAAI,CAAC,oBAAoB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACtD,MAAM,OAAO,GAAG;YACd,YAAY,CAAC,IAAI,CAAC,KAAK,CAAC;YACxB,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3F,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3F,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3F,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3F,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;YAC3F,YAAY,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,CAAC;SAC5F,CAAC;QACF,MAAM,MAAM,GAAG,eAAe,CAAC,YAAY,CAAC,OAAO,EAAE,OAAO,EAAE;YAC5D,QAAQ,EAAE,IAAI;YACd,UAAU,EAAE,SAAS;YACrB,UAAU,EAAE,SAAS;YACrB,SAAS;SACV,CAAC,CAAC;QAEH,OAAO,IAAI,KAAK,CACd,IAAI,CAAC,KAAK,EACV,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,EAC7C,IAAI,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAC9C,CAAC;IACJ,CAAC;CACF"}

package/lib/propagator/DormandPrince54Propagator.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 DormandPrince54Propagator 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/DormandPrince54Propagator.js

@@ -0,0 +1,68 @@
+/**
+ * @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 class-methods-use-this */
+import { RungeKuttaAdaptive } from './RungeKuttaAdaptive';
+export class DormandPrince54Propagator extends RungeKuttaAdaptive {
+    constructor() {
+        super(...arguments);
+        this._a = Float64Array.from([0.0, 1.0 / 5.0, 3.0 / 10.0, 4.0 / 5.0, 8.0 / 9.0, 1.0, 1.0]);
+        this._b = [
+            new Float64Array(),
+            Float64Array.from([1.0 / 5.0]),
+            Float64Array.from([3.0 / 40.0, 9.0 / 40.0]),
+            Float64Array.from([44.0 / 45.0, -56.0 / 15.0, 32.0 / 9.0]),
+            Float64Array.from([19372.0 / 6561.0, -25360.0 / 2187.0, 64448.0 / 6561.0, -212.0 / 729.0]),
+            Float64Array.from([9017.0 / 3168.0, -355.0 / 33.0, 46732.0 / 5247.0, 49.0 / 176.0, -5103.0 / 18656.0]),
+            Float64Array.from([35.0 / 384.0, 0.0, 500.0 / 1113.0, 125.0 / 192.0, -2187.0 / 6784.0, 11.0 / 84.0]),
+        ];
+        this._ch = Float64Array.from([
+            35.0 / 384.0,
+            0.0,
+            500.0 / 1113.0,
+            125.0 / 192.0,
+            -2187.0 / 6784.0,
+            11.0 / 84.0,
+            0.0,
+        ]);
+        this._c = Float64Array.from([
+            5179.0 / 57600.0,
+            0.0,
+            7571.0 / 16695.0,
+            393.0 / 640.0,
+            -92097.0 / 339200.0,
+            187.0 / 2100.0,
+            1.0 / 40.0,
+        ]);
+    }
+    get a() {
+        return this._a;
+    }
+    get b() {
+        return this._b;
+    }
+    get ch() {
+        return this._ch;
+    }
+    get c() {
+        return this._c;
+    }
+    get order() {
+        return 5;
+    }
+}
+//# sourceMappingURL=DormandPrince54Propagator.js.map

package/lib/propagator/DormandPrince54Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"DormandPrince54Propagator.js","sourceRoot":"","sources":["../../src/propagator/DormandPrince54Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,2CAA2C;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D,MAAM,OAAO,yBAA0B,SAAQ,kBAAkB;IAAjE;;QACU,OAAE,GAAiB,YAAY,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC;QAEnG,OAAE,GAAmB;YAC3B,IAAI,YAAY,EAAE;YAClB,YAAY,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;YAC9B,YAAY,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,EAAE,GAAG,GAAG,IAAI,CAAC,CAAC;YAC3C,YAAY,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,IAAI,EAAE,CAAC,IAAI,GAAG,IAAI,EAAE,IAAI,GAAG,GAAG,CAAC,CAAC;YAC1D,YAAY,CAAC,IAAI,CAAC,CAAC,OAAO,GAAG,MAAM,EAAE,CAAC,OAAO,GAAG,MAAM,EAAE,OAAO,GAAG,MAAM,EAAE,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;YAC1F,YAAY,CAAC,IAAI,CAAC,CAAC,MAAM,GAAG,MAAM,EAAE,CAAC,KAAK,GAAG,IAAI,EAAE,OAAO,GAAG,MAAM,EAAE,IAAI,GAAG,KAAK,EAAE,CAAC,MAAM,GAAG,OAAO,CAAC,CAAC;YACtG,YAAY,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,KAAK,EAAE,GAAG,EAAE,KAAK,GAAG,MAAM,EAAE,KAAK,GAAG,KAAK,EAAE,CAAC,MAAM,GAAG,MAAM,EAAE,IAAI,GAAG,IAAI,CAAC,CAAC;SACrG,CAAC;QAEM,QAAG,GAAiB,YAAY,CAAC,IAAI,CAAC;YAC5C,IAAI,GAAG,KAAK;YACZ,GAAG;YACH,KAAK,GAAG,MAAM;YACd,KAAK,GAAG,KAAK;YACb,CAAC,MAAM,GAAG,MAAM;YAChB,IAAI,GAAG,IAAI;YACX,GAAG;SACJ,CAAC,CAAC;QAEK,OAAE,GAAiB,YAAY,CAAC,IAAI,CAAC;YAC3C,MAAM,GAAG,OAAO;YAChB,GAAG;YACH,MAAM,GAAG,OAAO;YAChB,KAAK,GAAG,KAAK;YACb,CAAC,OAAO,GAAG,QAAQ;YACnB,KAAK,GAAG,MAAM;YACd,GAAG,GAAG,IAAI;SACX,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/KeplerPropagator.d.ts

@@ -0,0 +1,36 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { ClassicalElements, EpochUTC, J2000 } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+export declare class KeplerPropagator extends Propagator {
+    private _initElements;
+    private _elements;
+    private _cacheState;
+    private _checkpoints;
+    constructor(initElements: ClassicalElements);
+    get state(): J2000;
+    propagate(epoch: EpochUTC): J2000;
+    reset(): void;
+    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/KeplerPropagator.js

@@ -0,0 +1,81 @@
+/**
+ * @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 } from 'ootk-core';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { Propagator } from './Propagator';
+// / Kepler analytical two-body propagator.
+export class KeplerPropagator extends Propagator {
+    constructor(initElements) {
+        super();
+        this._initElements = initElements;
+        this._elements = initElements;
+        this._cacheState = J2000.fromClassicalElements(initElements);
+        this._checkpoints = [];
+    }
+    get state() {
+        return this._cacheState;
+    }
+    propagate(epoch) {
+        this._cacheState = J2000.fromClassicalElements(this._elements.propagate(epoch));
+        return this._cacheState;
+    }
+    reset() {
+        this._elements = this._initElements;
+        this._cacheState = J2000.fromClassicalElements(this._elements);
+    }
+    // eslint-disable-next-line @typescript-eslint/no-unused-vars
+    maneuver(maneuver, interval = 60) {
+        this._cacheState = maneuver.apply(this.propagate(maneuver.center));
+        this._elements = this._cacheState.toClassicalElements();
+        return [this._cacheState];
+    }
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        const tMvr = maneuvers.slice(0).filter((mvr) => mvr.center >= start || mvr.center <= finish);
+        const ephemeris = [];
+        if (tMvr[0].start > start) {
+            ephemeris.push(this.propagate(start));
+        }
+        for (const mvr of tMvr) {
+            while (this._cacheState.epoch < mvr.center) {
+                const step = Math.min(mvr.center.difference(this._cacheState.epoch), interval);
+                this.propagate(this._cacheState.epoch.roll(step));
+                if (this._cacheState.epoch.posix !== mvr.center.posix) {
+                    ephemeris.push(this._cacheState);
+                }
+            }
+            ephemeris.push(...this.maneuver(mvr, interval));
+        }
+        while (this._cacheState.epoch < finish) {
+            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(this._cacheState);
+        return this._checkpoints.length - 1;
+    }
+    clearCheckpoints() {
+        this._checkpoints = [];
+    }
+    restore(index) {
+        this._cacheState = this._checkpoints[index];
+    }
+}
+//# sourceMappingURL=KeplerPropagator.js.map

package/lib/propagator/KeplerPropagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"KeplerPropagator.js","sourceRoot":"","sources":["../../src/propagator/KeplerPropagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAA+B,KAAK,EAAE,MAAM,WAAW,CAAC;AAE/D,OAAO,EAAE,uBAAuB,EAAE,MAAM,yCAAyC,CAAC;AAClF,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C,2CAA2C;AAC3C,MAAM,OAAO,gBAAiB,SAAQ,UAAU;IAM9C,YAAY,YAA+B;QACzC,KAAK,EAAE,CAAC;QACR,IAAI,CAAC,aAAa,GAAG,YAAY,CAAC;QAClC,IAAI,CAAC,SAAS,GAAG,YAAY,CAAC;QAC9B,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,qBAAqB,CAAC,YAAY,CAAC,CAAC;QAC7D,IAAI,CAAC,YAAY,GAAG,EAAE,CAAC;IACzB,CAAC;IAED,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,SAAS,CAAC,KAAe;QACvB,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,qBAAqB,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;QAEhF,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,KAAK;QACH,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,aAAa,CAAC;QACpC,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,qBAAqB,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IACjE,CAAC;IAED,6DAA6D;IAC7D,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,EAAE;QACtC,IAAI,CAAC,WAAW,GAAG,QAAQ,CAAC,KAAK,CAAC,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC;QACnE,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,WAAW,CAAC,mBAAmB,EAAE,CAAC;QAExD,OAAO,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;IAC5B,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,MAAM,IAAI,KAAK,IAAI,GAAG,CAAC,MAAM,IAAI,MAAM,CAAC,CAAC;QAC7F,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,MAAM,EAAE;gBAC1C,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,UAAU,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,EAAE,QAAQ,CAAC,CAAC;gBAE/E,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,MAAM,CAAC,KAAK,EAAE;oBACrD,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,GAAG,MAAM,EAAE;YACtC,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,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/Propagator.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 } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+import { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+export declare abstract class Propagator {
+    abstract propagate(epoch: EpochUTC): J2000;
+    ephemeris(start: EpochUTC, stop: EpochUTC, interval?: number): VerletBlendInterpolator;
+    abstract reset(): void;
+    abstract checkpoint(): number;
+    abstract restore(index: number): void;
+    abstract clearCheckpoints(): void;
+    abstract get state(): J2000;
+    maneuver(maneuver: Thrust, interval?: number): J2000[];
+    ephemerisManeuver(start: EpochUTC, finish: EpochUTC, maneuvers: Thrust[], interval?: number): VerletBlendInterpolator;
+    ascendingNodeEpoch(start: EpochUTC): EpochUTC;
+    descendingNodeEpoch(start: EpochUTC): EpochUTC;
+    apogeeEpoch(start: EpochUTC): EpochUTC;
+    perigeeEpoch(start: EpochUTC): EpochUTC;
+}

package/lib/propagator/Propagator.js

@@ -0,0 +1,163 @@
+/**
+ * @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 { VerletBlendInterpolator } from '../interpolator/VerletBlendInterpolator';
+import { GoldenSection } from './../optimize/GoldenSection';
+// Propagator base class.
+export class Propagator {
+    /*
+     * Generate a [VerletBlendInterpolator] containing ephemeris over the
+     * [start] and [stop] propagation period, with an optional
+     * ephemeris [interval].
+     */
+    ephemeris(start, stop, interval = 60.0) {
+        const output = [this.propagate(start)];
+        let tempEpoch = start;
+        while (tempEpoch <= stop) {
+            tempEpoch = tempEpoch.roll(interval);
+            output.push(this.propagate(tempEpoch));
+        }
+        return new VerletBlendInterpolator(output);
+    }
+    /*
+     * Generate a list of [J2000] states integrating over a maneuver.
+     *
+     * If the maneuver is impulsive, the list will only contain a single state.
+     */
+    maneuver(maneuver, interval = 60.0) {
+        const output = [];
+        const tempEpoch = maneuver.start;
+        const stop = maneuver.stop;
+        output.push(this.propagate(tempEpoch));
+        let current = tempEpoch;
+        while (current <= stop) {
+            current = current.roll(interval);
+            output.push(this.propagate(current));
+        }
+        return output;
+    }
+    /*
+     * Generate a [VerletBlendInterpolator] containing maneuver ephemeris over
+     * the [start] and [finish] interval, with an optional ephemeris [interval].
+     */
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        const output = [];
+        const tempEpoch = start;
+        output.push(this.propagate(tempEpoch));
+        const stop = finish;
+        let current = tempEpoch;
+        while (current <= stop) {
+            current = current.roll(interval);
+            output.push(this.propagate(current));
+        }
+        return new VerletBlendInterpolator(output);
+    }
+    // Return the epoch of the ascending node after the [start] epoch.
+    ascendingNodeEpoch(start) {
+        const period = this.state.period();
+        const step = period / 8;
+        let current = start;
+        const stop = current.roll(period);
+        this.propagate(current);
+        let previous = this.state.position.z;
+        while (current <= stop) {
+            current = current.roll(step);
+            this.propagate(current);
+            if (Math.sign(this.state.position.z) === Math.sign(-previous) && this.state.velocity.z > 0) {
+                break;
+            }
+            previous = this.state.position.z;
+        }
+        const t = GoldenSection.search((x) => {
+            this.propagate(new EpochUTC(x));
+            return Math.abs(this.state.position.z);
+        }, current.posix - step, current.posix, { tolerance: 1e-3 });
+        return new EpochUTC(t);
+    }
+    // Return the epoch of the descending node after the [start] epoch.
+    descendingNodeEpoch(start) {
+        const period = this.state.period();
+        const step = period / 8;
+        let current = start;
+        const stop = current.roll(period);
+        this.propagate(current);
+        let previous = this.state.position.z;
+        while (current <= stop) {
+            current = current.roll(step);
+            this.propagate(current);
+            if (Math.sign(this.state.position.z) === Math.sign(-previous) && this.state.velocity.z < 0) {
+                break;
+            }
+            previous = this.state.position.z;
+        }
+        const t = GoldenSection.search((x) => {
+            this.propagate(new EpochUTC(x));
+            return Math.abs(this.state.position.z);
+        }, current.posix - step, current.posix, { tolerance: 1e-3 });
+        return new EpochUTC(t);
+    }
+    // Return the epoch of apogee after the [start] epoch.
+    apogeeEpoch(start) {
+        const slice = 8;
+        const period = this.state.period();
+        const step = period / slice;
+        let current = start;
+        this.propagate(current);
+        let tCache = current;
+        let rCache = this.state.position.magnitude();
+        for (let i = 0; i < slice; i++) {
+            current = current.roll(step);
+            const t = new EpochUTC(GoldenSection.search((x) => {
+                this.propagate(new EpochUTC(x));
+                return this.state.position.magnitude();
+            }, current.posix - step, current.posix, { tolerance: 1e-3, solveMax: true }));
+            this.propagate(t);
+            const r = this.state.position.magnitude();
+            if (r > rCache) {
+                tCache = t;
+                rCache = r;
+            }
+        }
+        return tCache;
+    }
+    // Return the epoch of perigee after the [start] epoch.
+    perigeeEpoch(start) {
+        const slice = 8;
+        const period = this.state.period();
+        const step = period / slice;
+        let current = start;
+        this.propagate(current);
+        let tCache = current;
+        let rCache = this.state.position.magnitude();
+        for (let i = 0; i < slice; i++) {
+            current = current.roll(step);
+            const t = new EpochUTC(GoldenSection.search((x) => {
+                this.propagate(new EpochUTC(x));
+                return this.state.position.magnitude();
+            }, current.posix - step, current.posix, { tolerance: 1e-3, solveMax: false }));
+            this.propagate(t);
+            const r = this.state.position.magnitude();
+            if (r < rCache) {
+                tCache = t;
+                rCache = r;
+            }
+        }
+        return tCache;
+    }
+}
+//# sourceMappingURL=Propagator.js.map

package/lib/propagator/Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Propagator.js","sourceRoot":"","sources":["../../src/propagator/Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,QAAQ,EAAS,MAAM,WAAW,CAAC;AAE5C,OAAO,EAAE,uBAAuB,EAAE,MAAM,yCAAyC,CAAC;AAClF,OAAO,EAAE,aAAa,EAAE,MAAM,6BAA6B,CAAC;AAE5D,yBAAyB;AACzB,MAAM,OAAgB,UAAU;IAI9B;;;;OAIG;IACH,SAAS,CAAC,KAAe,EAAE,IAAc,EAAE,QAAQ,GAAG,IAAI;QACxD,MAAM,MAAM,GAAY,CAAC,IAAI,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,CAAC;QAChD,IAAI,SAAS,GAAG,KAAK,CAAC;QAEtB,OAAO,SAAS,IAAI,IAAI,EAAE;YACxB,SAAS,GAAG,SAAS,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;YACrC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC,CAAC;SACxC;QAED,OAAO,IAAI,uBAAuB,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IAyBD;;;;OAIG;IACH,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,IAAI;QACxC,MAAM,MAAM,GAAY,EAAE,CAAC;QAC3B,MAAM,SAAS,GAAG,QAAQ,CAAC,KAAK,CAAC;QACjC,MAAM,IAAI,GAAG,QAAQ,CAAC,IAAI,CAAC;QAE3B,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC,CAAC;QACvC,IAAI,OAAO,GAAG,SAAS,CAAC;QAExB,OAAO,OAAO,IAAI,IAAI,EAAE;YACtB,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;YACjC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;SACtC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;OAGG;IACH,iBAAiB,CAAC,KAAe,EAAE,MAAgB,EAAE,SAAmB,EAAE,QAAQ,GAAG,IAAI;QACvF,MAAM,MAAM,GAAY,EAAE,CAAC;QAC3B,MAAM,SAAS,GAAG,KAAK,CAAC;QAExB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,MAAM,CAAC;QACpB,IAAI,OAAO,GAAG,SAAS,CAAC;QAExB,OAAO,OAAO,IAAI,IAAI,EAAE;YACtB,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;YACjC,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;SACtC;QAED,OAAO,IAAI,uBAAuB,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IAED,kEAAkE;IAClE,kBAAkB,CAAC,KAAe;QAChC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACnC,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC;QACxB,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAElC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;QACxB,IAAI,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAErC,OAAO,OAAO,IAAI,IAAI,EAAE;YACtB,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC7B,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;YACxB,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE;gBAC1F,MAAM;aACP;YACD,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;SAClC;QACD,MAAM,CAAC,GAAG,aAAa,CAAC,MAAM,CAC5B,CAAC,CAAS,EAAE,EAAE;YACZ,IAAI,CAAC,SAAS,CAAC,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;YAEhC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QACzC,CAAC,EACD,OAAO,CAAC,KAAK,GAAG,IAAI,EACpB,OAAO,CAAC,KAAK,EACb,EAAE,SAAS,EAAE,IAAI,EAAE,CACpB,CAAC;QAEF,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC;IAED,mEAAmE;IACnE,mBAAmB,CAAC,KAAe;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACnC,MAAM,IAAI,GAAG,MAAM,GAAG,CAAC,CAAC;QACxB,IAAI,OAAO,GAAG,KAAK,CAAC;QACpB,MAAM,IAAI,GAAG,OAAO,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAElC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;QACxB,IAAI,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;QAErC,OAAO,OAAO,IAAI,IAAI,EAAE;YACtB,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC7B,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;YACxB,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,IAAI,CAAC,CAAC,QAAQ,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,EAAE;gBAC1F,MAAM;aACP;YACD,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;SAClC;QACD,MAAM,CAAC,GAAG,aAAa,CAAC,MAAM,CAC5B,CAAC,CAAS,EAAE,EAAE;YACZ,IAAI,CAAC,SAAS,CAAC,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;YAEhC,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QACzC,CAAC,EACD,OAAO,CAAC,KAAK,GAAG,IAAI,EACpB,OAAO,CAAC,KAAK,EACb,EAAE,SAAS,EAAE,IAAI,EAAE,CACpB,CAAC;QAEF,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC;IAED,sDAAsD;IACtD,WAAW,CAAC,KAAe;QACzB,MAAM,KAAK,GAAG,CAAC,CAAC;QAChB,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACnC,MAAM,IAAI,GAAG,MAAM,GAAG,KAAK,CAAC;QAC5B,IAAI,OAAO,GAAG,KAAK,CAAC;QAEpB,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;QACxB,IAAI,MAAM,GAAG,OAAO,CAAC;QACrB,IAAI,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAE7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE;YAC9B,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC7B,MAAM,CAAC,GAAG,IAAI,QAAQ,CACpB,aAAa,CAAC,MAAM,CAClB,CAAC,CAAS,EAAE,EAAE;gBACZ,IAAI,CAAC,SAAS,CAAC,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;gBAEhC,OAAO,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;YACzC,CAAC,EACD,OAAO,CAAC,KAAK,GAAG,IAAI,EACpB,OAAO,CAAC,KAAK,EACb,EAAE,SAAS,EAAE,IAAI,EAAE,QAAQ,EAAE,IAAI,EAAE,CACpC,CACF,CAAC;YAEF,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;YAClB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;YAE1C,IAAI,CAAC,GAAG,MAAM,EAAE;gBACd,MAAM,GAAG,CAAC,CAAC;gBACX,MAAM,GAAG,CAAC,CAAC;aACZ;SACF;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,uDAAuD;IACvD,YAAY,CAAC,KAAe;QAC1B,MAAM,KAAK,GAAG,CAAC,CAAC;QAChB,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC;QACnC,MAAM,IAAI,GAAG,MAAM,GAAG,KAAK,CAAC;QAC5B,IAAI,OAAO,GAAG,KAAK,CAAC;QAEpB,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC;QACxB,IAAI,MAAM,GAAG,OAAO,CAAC;QACrB,IAAI,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAE7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE;YAC9B,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;YAC7B,MAAM,CAAC,GAAG,IAAI,QAAQ,CACpB,aAAa,CAAC,MAAM,CAClB,CAAC,CAAS,EAAE,EAAE;gBACZ,IAAI,CAAC,SAAS,CAAC,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;gBAEhC,OAAO,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;YACzC,CAAC,EACD,OAAO,CAAC,KAAK,GAAG,IAAI,EACpB,OAAO,CAAC,KAAK,EACb,EAAE,SAAS,EAAE,IAAI,EAAE,QAAQ,EAAE,KAAK,EAAE,CACrC,CACF,CAAC;YAEF,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;YAClB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;YAE1C,IAAI,CAAC,GAAG,MAAM,EAAE;gBACd,MAAM,GAAG,CAAC,CAAC;gBACX,MAAM,GAAG,CAAC,CAAC;aACZ;SACF;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/lib/propagator/RkCheckpoint.d.ts

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

package/lib/propagator/RkCheckpoint.js

@@ -0,0 +1,26 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+// / Runge-Kutta adaptive state checkpoint.
+export class RkCheckpoint {
+    constructor(cacheState, stepSize) {
+        this.cacheState = cacheState;
+        this.stepSize = stepSize;
+        // Nothing to do here.
+    }
+}
+//# sourceMappingURL=RkCheckpoint.js.map

package/lib/propagator/RkCheckpoint.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RkCheckpoint.js","sourceRoot":"","sources":["../../src/propagator/RkCheckpoint.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAIH,2CAA2C;AAC3C,MAAM,OAAO,YAAY;IACvB,YAAmB,UAAiB,EAAS,QAAgB;QAA1C,eAAU,GAAV,UAAU,CAAO;QAAS,aAAQ,GAAR,QAAQ,CAAQ;QAC3D,sBAAsB;IACxB,CAAC;CACF"}

package/lib/propagator/RkResult.d.ts

@@ -0,0 +1,24 @@
+/**
+ * @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 } from 'ootk-core';
+export declare class RkResult {
+    state: J2000;
+    error: number;
+    newStep: number;
+    constructor(state: J2000, error: number, newStep: number);
+}

package/lib/propagator/RkResult.js

@@ -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/>.
+ */
+// / Result of adaptive numerical integration.
+export class RkResult {
+    // / Create a new [RkResult] object.
+    constructor(state, error, newStep) {
+        this.state = state;
+        this.error = error;
+        this.newStep = newStep;
+        // Nothing to do here.
+    }
+}
+//# sourceMappingURL=RkResult.js.map

package/lib/propagator/RkResult.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RkResult.js","sourceRoot":"","sources":["../../src/propagator/RkResult.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAIH,8CAA8C;AAC9C,MAAM,OAAO,QAAQ;IACnB,oCAAoC;IACpC,YAAmB,KAAY,EAAS,KAAa,EAAS,OAAe;QAA1D,UAAK,GAAL,KAAK,CAAO;QAAS,UAAK,GAAL,KAAK,CAAQ;QAAS,YAAO,GAAP,OAAO,CAAQ;QAC3E,sBAAsB;IACxB,CAAC;CACF"}