ootk 3.0.0 → 4.0.0-1

package/dist/propagator/Propagator.js

@@ -0,0 +1,162 @@
+/**
+ * @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 / 60;
+        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 / 60;
+        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 / 60;
+        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 / 60;
+        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/dist/propagator/Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Propagator.js","sourceRoot":"","sources":["../../src/propagator/Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,QAAQ,EAAkB,MAAM,WAAW,CAAC;AAErD,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,WAAW,IAAe;QACnE,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,WAAW,IAAe;QACnD,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,CACf,KAAe,EACf,MAAgB,EAChB,SAAmB,EACnB,WAAW,IAAe;QAE1B,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,GAAG,EAAa,CAAC;QACjD,MAAM,IAAI,GAAG,MAAM,GAAG,CAAY,CAAC;QACnC,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,CAAC,EAAE,EAAE;YACJ,IAAI,CAAC,SAAS,CAAC,IAAI,QAAQ,CAAC,CAAY,CAAC,CAAC,CAAC;YAE3C,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAY,CAAC;QACpD,CAAC,EACD,OAAO,CAAC,KAAK,GAAG,IAAI,EACpB,OAAO,CAAC,KAAK,EACb,EAAE,SAAS,EAAE,IAAI,EAAE,CACT,CAAC;QAEb,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC;IACzB,CAAC;IAED,mEAAmE;IACnE,mBAAmB,CAAC,KAAe;QACjC,MAAM,MAAM,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,EAAa,CAAC;QACjD,MAAM,IAAI,GAAG,MAAM,GAAG,CAAY,CAAC;QACnC,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,CAAY,CAAC,CAAC,CAAC;YAE3C,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,CACT,CAAC;QAEb,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,GAAG,EAAa,CAAC;QACjD,MAAM,IAAI,GAAG,MAAM,GAAG,KAAgB,CAAC;QACvC,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,CAAY,CAAC,CAAC,CAAC;gBAE3C,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,CACzB,CACb,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,GAAG,EAAa,CAAC;QACjD,MAAM,IAAI,GAAG,MAAM,GAAG,KAAgB,CAAC;QACvC,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,CAAY,CAAC,CAAC,CAAC;gBAE3C,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,CAC1B,CACb,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/dist/propagator/RkCheckpoint.d.ts

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

package/dist/propagator/RkCheckpoint.js

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

package/dist/propagator/RkCheckpoint.js.map

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

package/dist/propagator/RkResult.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 RkResult {
+    state: J2000;
+    error: number;
+    newStep: number;
+    constructor(state: J2000, error: number, newStep: number);
+}

package/dist/propagator/RkResult.js

@@ -0,0 +1,30 @@
+/**
+ * @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 {
+    state;
+    error;
+    newStep;
+    // / 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/dist/propagator/RkResult.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RkResult.js","sourceRoot":"","sources":["../../src/propagator/RkResult.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAIH,8CAA8C;AAC9C,MAAM,OAAO,QAAQ;IAEA;IAAqB;IAAsB;IAD9D,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"}

package/dist/propagator/RungeKutta4Propagator.d.ts

@@ -0,0 +1,50 @@
+/**
+ * @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, Seconds } 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.
+     * @param initState_ Initial state vector.
+     * @param forceModel_ Numerical integration force model.
+     * @param stepSize_ Integration step size _(seconds)_.
+     * @param cacheState_ Cached state vector.
+     * @param checkpoints_ Cached state vector checkpoints.
+     */
+    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?: Seconds): 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/dist/propagator/RungeKutta4Propagator.js

@@ -0,0 +1,141 @@
+/**
+ * @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 {
+    initState_;
+    forceModel_;
+    stepSize_;
+    cacheState_;
+    checkpoints_;
+    /**
+     * Create a new [RungeKutta4Propagator] object from an initial state vector and
+     * along with an optional [ForceModel] and [stepSize] in seconds.
+     * @param initState_ Initial state vector.
+     * @param forceModel_ Numerical integration force model.
+     * @param stepSize_ Integration step size _(seconds)_.
+     * @param cacheState_ Cached state vector.
+     * @param checkpoints_ Cached state vector checkpoints.
+     */
+    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/dist/propagator/RungeKutta4Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKutta4Propagator.js","sourceRoot":"","sources":["../../src/propagator/RungeKutta4Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAY,KAAK,EAA4C,MAAM,EAAY,MAAM,WAAW,CAAC;AACxG,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;IAWzC;IACA;IACA;IACA;IACA;IAdV;;;;;;;;OAQG;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,CACf,KAAe,EACf,MAAgB,EAChB,SAAmB,EACnB,WAAW,IAAe;QAE1B,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,CAAY,CAAC;gBAEzF,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,CAAY,CAAC;YAEtF,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,CAAY,CAAC;YAEnF,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,IAAa,EAAE,IAAY;QACpD,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,CACtB,KAAK,EACL,MAAM,CAAC,UAAU,CAAC,CAAC,CAAyB,EAC5C,MAAM,CAAC,UAAU,CAAC,CAAC,CAAkC,CACtD,CAAC;QAEF,OAAO,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC;IAC7C,CAAC;IAEO,UAAU,CAAC,KAAY,EAAE,IAAa;QAC5C,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,CAAY,EAAE,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QACtE,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,GAAG,IAAe,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9E,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,GAAG,GAAG,IAAe,EAAE,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;QAC9E,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,CACd,KAAK,EACL,MAAM,CAAC,UAAU,CAAC,CAAC,CAAyB,EAC5C,MAAM,CAAC,UAAU,CAAC,CAAC,CAAkC,CACtD,CAAC;IACJ,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,SAAoB,CAAC;YAE5E,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/dist/propagator/RungeKutta89Propagator.d.ts

@@ -0,0 +1,28 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { 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/dist/propagator/RungeKutta89Propagator.js

@@ -0,0 +1,126 @@
+/**
+ * @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 {
+    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,
+    ]);
+    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,
+        ]),
+    ];
+    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,
+    ]);
+    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/dist/propagator/RungeKutta89Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKutta89Propagator.js","sourceRoot":"","sources":["../../src/propagator/RungeKutta89Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,4DAA4D;AAC5D,2CAA2C;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,sBAAsB,CAAC;AAE1D,oDAAoD;AACpD,MAAM,OAAO,sBAAuB,SAAQ,kBAAkB;IACpD,EAAE,GAAG,YAAY,CAAC,IAAI,CAAC;QAC7B,GAAG,EAAE,SAAS,EAAE,6CAA6C,EAAE,0CAA0C,EAAE,KAAK;QAChH,0CAA0C,EAAE,0CAA0C,EAAE,KAAK,EAAE,OAAO,EAAE,SAAS;QACjH,IAAI,EAAE,0CAA0C,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,EAAE,GAAG;KAC3E,CAAC,CAAC;IACK,EAAE,GAAG;QACX,IAAI,YAAY,CAAC,CAAC,CAAC;QACnB,YAAY,CAAC,IAAI,CAAC,CAAC,SAAS,CAAC,CAAC;QAC9B,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,4CAA4C,EAAE,0CAA0C,CAAC,CAAC;QAC9G,YAAY,CAAC,IAAI,CAAC,CAAC,6CAA6C,EAAE,GAAG,EAAE,0CAA0C,CAAC,CAAC;QACnH,YAAY,CAAC,IAAI,CAAC;YAChB,yCAAyC,EAAE,GAAG,EAAE,CAAC,yCAAyC;YAC1F,yCAAyC;SAC1C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;YACnG,2CAA2C;SAC5C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,0CAA0C;YACjG,4CAA4C,EAAE,wCAAwC;SACvF,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;YAC7G,0CAA0C;SAC3C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,gBAAgB,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,0CAA0C;YAChF,0CAA0C,EAAE,CAAC,gBAAgB;SAC9D,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,6CAA6C;YAChH,0CAA0C,EAAE,CAAC,6CAA6C;YAC1F,CAAC,0CAA0C;SAC5C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,CAAC,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,6CAA6C;YACjH,CAAC,0CAA0C,EAAE,6CAA6C;YAC1F,0CAA0C,EAAE,yCAAyC;SACtF,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,0CAA0C;YAC9G,0CAA0C,EAAE,6CAA6C;YACzF,CAAC,2CAA2C,EAAE,6CAA6C;YAC3F,0CAA0C;SAC3C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,CAAC,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,yCAAyC;YAC3G,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;YACtF,yCAAyC,EAAE,yCAAyC;YACpF,yCAAyC,EAAE,yCAAyC;SACrF,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,0CAA0C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,yCAAyC;YACzG,CAAC,yCAAyC,EAAE,yCAAyC;YACrF,yCAAyC,EAAE,CAAC,yCAAyC;YACrF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;YACtF,0CAA0C;SAC3C,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,CAAC,yCAAyC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,yCAAyC;YAC1G,yCAAyC,EAAE,CAAC,yCAAyC;YACrF,CAAC,yCAAyC,EAAE,yCAAyC;YACrF,yCAAyC,EAAE,yCAAyC;YACpF,CAAC,0CAA0C,EAAE,0CAA0C;SACxF,CAAC;QACF,YAAY,CAAC,IAAI,CAAC;YAChB,yCAAyC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,yCAAyC;YACxG,wCAAwC,EAAE,yCAAyC;YACnF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;YACtF,CAAC,yCAAyC,EAAE,CAAC,yCAAyC;YACtF,yCAAyC,EAAE,GAAG,EAAE,GAAG;SACpD,CAAC;KACH,CAAC;IACM,GAAG,GAAG,YAAY,CAAC,IAAI,CAAC;QAC9B,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG;QAC3E,yCAAyC,EAAE,6CAA6C;QACxF,0CAA0C,EAAE,0CAA0C;QACtF,CAAC,yCAAyC,EAAE,yCAAyC,EAAE,GAAG,EAAE,GAAG;QAC/F,6CAA6C;KAC9C,CAAC,CAAC;IACK,EAAE,GAAG,YAAY,CAAC,IAAI,CAAC;QAC7B,6CAA6C,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG;QAC3E,CAAC,yCAAyC,EAAE,0CAA0C;QACtF,0CAA0C,EAAE,0CAA0C;QACtF,0CAA0C,EAAE,6CAA6C;QACzF,0CAA0C,EAAE,6CAA6C,EAAE,GAAG;KAC/F,CAAC,CAAC;IAEH,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"}