ootk 3.0.0 → 4.0.0-1

package/dist/propagator/RungeKuttaAdaptive.d.ts

@@ -0,0 +1,54 @@
+/**
+ * @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 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].
+     * @param initState_ Initial state vector.
+     * @param forceModel_ Numerical integration force model.
+     * @param tolerance_ Minimum allowable local error 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?: Seconds): J2000[];
+    ephemerisManeuver(start: EpochUTC, finish: EpochUTC, maneuvers: Thrust[], interval?: Seconds): VerletBlendInterpolator;
+    checkpoint(): number;
+    clearCheckpoints(): void;
+    restore(index: number): void;
+}

package/dist/propagator/RungeKuttaAdaptive.js

@@ -0,0 +1,159 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { 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 {
+    initState_;
+    forceModel_;
+    tolerance_;
+    /**
+     * Create a new [RungeKuttaAdaptive] object from an initial state vector
+     * along with an optional [ForceModel] and [tolerance].
+     * @param initState_ Initial state vector.
+     * @param forceModel_ Numerical integration force model.
+     * @param tolerance_ Minimum allowable local error tolerance.
+     */
+    constructor(initState_, forceModel_ = new ForceModel().setGravity(), tolerance_ = 1e-9) {
+        super();
+        this.initState_ = initState_;
+        this.forceModel_ = forceModel_;
+        this.tolerance_ = tolerance_;
+        this._cacheState = this.initState_;
+        this.tolerance_ = Math.max(RungeKuttaAdaptive._minTolerance, Math.abs(tolerance_));
+    }
+    // / Initial state vector.
+    _cacheState;
+    _checkpoints = [];
+    // / Integration step size _(seconds)_.
+    _stepSize = 60.0;
+    // / Minimum allowable local error tolerance.
+    static _minTolerance = 1e-15;
+    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;
+    }
+}
+//# sourceMappingURL=RungeKuttaAdaptive.js.map

package/dist/propagator/RungeKuttaAdaptive.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RungeKuttaAdaptive.js","sourceRoot":"","sources":["../../src/propagator/RungeKuttaAdaptive.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,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;IAS/C;IACA;IACA;IAVV;;;;;;OAMG;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;QAGjC,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;IAED,0BAA0B;IAClB,WAAW,CAAQ;IAEV,YAAY,GAAmB,EAAE,CAAC;IAEnD,uCAAuC;IAC/B,SAAS,GAAG,IAAI,CAAC;IAEzB,6CAA6C;IACrC,MAAM,CAAU,aAAa,GAAW,KAAK,CAAC;IAiBtD,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,CACV,KAAe,EACf,EAA0C,EAC1C,IAAa,EACb,IAAwB,EACxB,IAAa;QAEb,MAAM,CAAC,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,IAAe,CAAC,CAAC;QAC7C,MAAM,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,KAAK,CACtB,CAAC,EACD,KAAK,CAAC,UAAU,CAAC,CAAC,CAAyB,EAC3C,KAAK,CAAC,UAAU,CAAC,CAAC,CAAkC,CACrD,CAAC;QAEF,OAAO,IAAI,CAAC,WAAW,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC;IACzD,CAAC;IAEO,UAAU,CAAC,KAAY,EAAE,IAAa;QAC5C,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,CAAuB,CAAC;QAEpE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACtC,IAAI,IAAI,GAAG,MAAM,CAAC,OAA6B,CAAC;YAEhD,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,CAAuB,CAAC;iBACjE;aACF;YACD,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAY,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;SACpE;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,CAAuB,CAAC;YAC1D,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAuB,CAAC;SAC1D;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,CACjB,IAAI,KAAK,CACP,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EACtB,EAAE,CAAC,UAAU,CAAC,CAAC,CAAyB,EACxC,EAAE,CAAC,UAAU,CAAC,CAAC,CAAkC,CAClD,EACD,KAAK,EACL,IAAI,CACL,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;YAC5E,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,WAAW,IAAe;QACnD,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;IAED,iBAAiB,CACf,KAAe,EACf,MAAgB,EAChB,SAAmB,EACnB,WAAW,IAAe;QAE1B,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,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,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"}

package/dist/propagator/Sgp4Propagator.d.ts

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

package/dist/propagator/Sgp4Propagator.js

@@ -0,0 +1,56 @@
+/**
+ * @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 { Propagator } from './Propagator';
+// / SGP4 propagator.
+export class Sgp4Propagator extends Propagator {
+    tle;
+    // / Create a new [Sgp4Propagator] object from a [TLE].
+    constructor(tle) {
+        super();
+        this.tle = tle;
+        this._cacheState = tle.state.toJ2000();
+    }
+    _cacheState;
+    _checkpoints = [];
+    get state() {
+        return this._cacheState;
+    }
+    ephemerisManeuver(start, finish, maneuvers, interval = 60.0) {
+        throw new Error('Maneuvers cannot be modelled with SGP4.');
+    }
+    maneuver(maneuver, interval = 60.0) {
+        throw new Error('Maneuvers cannot be modelled with SGP4.');
+    }
+    propagate(epoch) {
+        this._cacheState = this.tle.propagate(epoch).toJ2000();
+        return this._cacheState;
+    }
+    reset() {
+        this._cacheState = this.tle.state.toJ2000();
+    }
+    checkpoint() {
+        this._checkpoints.push(this._cacheState);
+        return this._checkpoints.length - 1;
+    }
+    clearCheckpoints() {
+        this._checkpoints = [];
+    }
+    restore(index) {
+        this._cacheState = this._checkpoints[index];
+    }
+}
+//# sourceMappingURL=Sgp4Propagator.js.map

package/dist/propagator/Sgp4Propagator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Sgp4Propagator.js","sourceRoot":"","sources":["../../src/propagator/Sgp4Propagator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAOH,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAE1C,qBAAqB;AACrB,MAAM,OAAO,cAAe,SAAQ,UAAU;IAExB;IADpB,uDAAuD;IACvD,YAAoB,GAAQ;QAC1B,KAAK,EAAE,CAAC;QADU,QAAG,GAAH,GAAG,CAAK;QAE1B,IAAI,CAAC,WAAW,GAAG,GAAG,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IACzC,CAAC;IAEO,WAAW,CAAQ;IACnB,YAAY,GAAY,EAAE,CAAC;IAEnC,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,iBAAiB,CAAC,KAAe,EAAE,MAAgB,EAAE,SAAmB,EAAE,QAAQ,GAAG,IAAI;QACvF,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,QAAQ,CAAC,QAAgB,EAAE,QAAQ,GAAG,IAAI;QACxC,MAAM,IAAI,KAAK,CAAC,yCAAyC,CAAC,CAAC;IAC7D,CAAC;IAED,SAAS,CAAC,KAAe;QACvB,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,OAAO,EAAE,CAAC;QAEvD,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,KAAK;QACH,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC;IAC9C,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/smoothing/ExponentialSmoothing.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 { EpochUTC } from 'ootk-core';
+export declare class ExponentialSmoothing {
+    private constructor();
+    static smooth(xs: number[], alpha: number): number[];
+    static smoothDouble(xs: number[], alpha: number, beta: number): number[];
+    static smoothTime(epochs: EpochUTC[], xs: number[], timeConstant: number): number[];
+}

package/dist/smoothing/ExponentialSmoothing.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/>.
+ */
+// / Exponential data smoothing methods.
+export class ExponentialSmoothing {
+    constructor() {
+        // disable constructor
+    }
+    /*
+     * Perform exponential smoothing on data set [xs] using the provided
+     * data smoothing factor [alpha] _(0.0 <= alpha <= 1.0)_.
+     */
+    static smooth(xs, alpha) {
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+            }
+            else {
+                ss.push(alpha * xs[i] + (1 - alpha) * ss[i - 1]);
+            }
+        }
+        return ss;
+    }
+    /*
+     * Perform exponential smoothing on data set [xs] using the provided data
+     * smoothing factor [alpha] _(0.0 <= alpha <= 1.0)_ and trend smoothing
+     * factor [beta] _(0.0 <= beta <= 1.0)_.
+     */
+    static smoothDouble(xs, alpha, beta) {
+        const bs = [];
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+                bs.push(xs[1] - xs[0]);
+            }
+            else {
+                ss.push(alpha * xs[i] + (1 - alpha) * (ss[i - 1] + bs[i - 1]));
+                bs.push(beta * (ss[i] - ss[i - 1]) + (1 - beta) * bs[i - 1]);
+            }
+        }
+        return ss;
+    }
+    /*
+     * Perform exponential smoothing on time series data set [xs] correlated
+     * with the [epochs] array using the provided [timeConstant]
+     * value _(seconds)_.
+     *
+     * Note: the [timeConstant] is the amount of time for the smoothed response
+     * of a unit step function to reach ~63.2% of the original signal.
+     */
+    static smoothTime(epochs, xs, timeConstant) {
+        const ts = epochs.map((e) => e.posix);
+        const ss = [];
+        for (let i = 0; i < xs.length; i++) {
+            if (i === 0) {
+                ss.push(xs[0]);
+            }
+            else {
+                const a = 1 - Math.exp(-(ts[i] - ts[i - 1]) / timeConstant);
+                ss.push(a * xs[i] + (1 - a) * ss[i - 1]);
+            }
+        }
+        return ss;
+    }
+}
+//# sourceMappingURL=ExponentialSmoothing.js.map

package/dist/smoothing/ExponentialSmoothing.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ExponentialSmoothing.js","sourceRoot":"","sources":["../../src/smoothing/ExponentialSmoothing.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAIH,wCAAwC;AACxC,MAAM,OAAO,oBAAoB;IAC/B;QACE,sBAAsB;IACxB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,MAAM,CAAC,EAAY,EAAE,KAAa;QACvC,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aAChB;iBAAM;gBACL,EAAE,CAAC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAClD;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,YAAY,CAAC,EAAY,EAAE,KAAa,EAAE,IAAY;QAC3D,MAAM,EAAE,GAAa,EAAE,CAAC;QACxB,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACf,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aACxB;iBAAM;gBACL,EAAE,CAAC,IAAI,CAAC,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/D,EAAE,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAC9D;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;IAED;;;;;;;OAOG;IACH,MAAM,CAAC,UAAU,CAAC,MAAkB,EAAE,EAAY,EAAE,YAAoB;QACtE,MAAM,EAAE,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;QACtC,MAAM,EAAE,GAAa,EAAE,CAAC;QAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YAClC,IAAI,CAAC,KAAK,CAAC,EAAE;gBACX,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aAChB;iBAAM;gBACL,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,YAAY,CAAC,CAAC;gBAE5D,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aAC1C;SACF;QAED,OAAO,EAAE,CAAC;IACZ,CAAC;CACF"}

package/dist/transforms/index.d.ts

@@ -0,0 +1 @@
+export * from './transforms';

package/dist/transforms/index.js

@@ -0,0 +1,2 @@
+export * from './transforms';
+//# sourceMappingURL=index.js.map

package/dist/transforms/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/transforms/index.ts"],"names":[],"mappings":"AAAA,cAAc,cAAc,CAAC"}

package/dist/transforms/transforms.d.ts

@@ -1,94 +1,58 @@
+import { Degrees, EnuVec3, Radians, RaeVec3 } from 'ootk-core';
+import { RfSensor } from 'src/objects/RfSensor';
+import { RfVec3, RuvVec3 } from 'src/types/types';
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
- *
- * @file The Transforms module contains a collection of conversions not contained
- * in the original SGP4 library such as ECI to ECF and ECF to RAE. This was based
- * on some of the functions in satellite.js.
- *
- * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 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/>.
+ * Converts Azimuth and Elevation to U and V.
+ * Azimuth is the angle off of boresight in the horizontal plane.
+ * Elevation is the angle off of boresight in the vertical plane.
+ * Cone half angle is the angle of the cone of the radar max field of view.
+ * @param az - Azimuth in radians
+ * @param el - Elevation in radians
+ * @param coneHalfAngle - Cone half angle in radians
+ * @returns U and V in radians
  */
-import { Degrees, EcefVec3, EcfVec3, EciVec3, EnuVec3, Kilometers, LlaVec3, RadarSensor, Radians, RaeVec3, RfVec3, RuvVec3, SezVec3 } from '../types/types';
-declare class Transforms {
-    static getDegLat(radians: Radians): Degrees;
-    static getDegLon(radians: Radians): Degrees;
-    static getRadLat(degrees: Degrees): Radians;
-    static getRadLon(degrees: Degrees): Radians;
-    static rad2deg(radians: Radians): Degrees;
-    static deg2rad(degrees: Degrees): Radians;
-    static ecf2eci(ecf: EcfVec3, gmst: number): EciVec3;
-    static ecf2rae(lla: LlaVec3, ecf: EcfVec3): RaeVec3;
-    /**
-     * EciToGeodetic converts eci coordinates to lla coordinates
-     * @param {vec3} eci takes xyz coordinates
-     * @param {number} gmst takes a number in gmst time
-     * @returns {array} array containing lla coordinates
-     */
-    static eci2lla(eci: EciVec3, gmst: number): LlaVec3;
-    static eci2ecf(eci: EciVec3, gmst: number): EcfVec3;
-    static lla2ecf(lla: LlaVec3): EcfVec3;
-    static lla2sez(lla: LlaVec3, ecf: EcfVec3): SezVec3;
-    static rae2sez(rae: RaeVec3): SezVec3;
-    static rae2ecf(rae: RaeVec3, lla: LlaVec3): EcfVec3;
-    /**
-     * @param {Object} sez Containing SEZ coordinates
-     * @param {Number} sez.s Positive horizontal vector S due south.
-     * @param {Number} sez.e Positive horizontal vector E due east.
-     * @param {Number} sez.z Vector Z normal to the surface of the earth (up).
-     * @returns {RaeVec3} Rng, Az, El array
-     */
-    static sez2rae(sez: SezVec3): RaeVec3;
-    /**
-     * Calculates Geodetic Lat Lon Alt to ECEF coordinates.
-     */
-    static lla2ecef(lla: LlaVec3): EcefVec3<Kilometers>;
-    static ecf2enu(ecf: EcefVec3<Kilometers>, lla: LlaVec3): EnuVec3<Kilometers>;
-    static rae2enu(rae: RaeVec3): EnuVec3<Kilometers>;
-    static enu2rf({ x, y, z }: EnuVec3<Kilometers>, az: Radians, el: Radians): RfVec3;
-    /**
-     * Converts Azimuth and Elevation to U and V.
-     * Azimuth is the angle off of boresight in the horizontal plane.
-     * Elevation is the angle off of boresight in the vertical plane.
-     * Cone half angle is the angle of the cone of the radar max field of view.
-     */
-    static azel2uv(az: Radians, el: Radians, coneHalfAngle: Radians): {
-        u: Radians;
-        v: Radians;
-    };
-    /**
-     * Converts U and V to Azimuth and Elevation off of boresight.
-     */
-    static uv2azel(u: Radians, v: Radians, coneHalfAngle: Radians): {
-        az: Radians;
-        el: Radians;
-    };
-    /**
-     * Converts Range Az El to Range U V.
-     */
-    static rae2ruv(rae: RaeVec3, sensor: RadarSensor, maxSensorAz: Degrees): RuvVec3;
-    /**
-     * Determine azimuth and elevation off of boresight based on sensor orientation and RAE.
-     *
-     * @param {rae} rae Range, Azimuth, Elevation
-     * @param {RadarSensor} sensor Radar sensor object
-     * @returns {az, el} Azimuth and Elevation off of boresight
-     */
-    static rae2raeOffBoresight(rae: RaeVec3, sensor: RadarSensor, maxSensorAz: Degrees): {
-        az: Radians;
-        el: Radians;
-    };
-}
-export { Transforms };
+export declare function azel2uv(az: Radians, el: Radians, coneHalfAngle: Radians): {
+    u: number;
+    v: number;
+};
+/**
+ * Determine azimuth and elevation off of boresight based on sensor orientation and RAE.
+ * @param rae Range, Azimuth, Elevation
+ * @param sensor Radar sensor object
+ * @param maxSensorAz Maximum sensor azimuth
+ * @returns Azimuth and Elevation off of boresight
+ */
+export declare function rae2raeOffBoresight(rae: RaeVec3, sensor: RfSensor, maxSensorAz: Degrees): {
+    az: Radians;
+    el: Radians;
+};
+/**
+ * Converts Range Az El to Range U V.
+ * @param rae Range, Azimuth, Elevation
+ * @param sensor Radar sensor object
+ * @param maxSensorAz Maximum sensor azimuth
+ * @returns Range, U, V
+ */
+export declare function rae2ruv(rae: RaeVec3, sensor: RfSensor, maxSensorAz: Degrees): RuvVec3;
+/**
+ * Converts U and V to Azimuth and Elevation off of boresight.
+ * @param u The U coordinate.
+ * @param v The V coordinate.
+ * @param coneHalfAngle The cone half angle of the radar.
+ * @returns Azimuth and Elevation off of boresight.
+ */
+export declare function uv2azel(u: number, v: number, coneHalfAngle: Radians): {
+    az: Radians;
+    el: Radians;
+};
+/**
+ * Converts coordinates from East-North-Up (ENU) to Right-Front-Up (RF) coordinate system.
+ * @param enu - The ENU coordinates to be converted.
+ * @param enu.x - The east coordinate.
+ * @param enu.y - The north coordinate.
+ * @param enu.z - The up coordinate.
+ * @param az - The azimuth angle in radians.
+ * @param el - The elevation angle in radians.
+ * @returns The converted RF coordinates.
+ */
+export declare function enu2rf<D extends number, A extends number = Radians>({ x, y, z }: EnuVec3<D>, az: A, el: A): RfVec3<D>;