ootk 3.0.0 → 4.0.0-0

package/lib/covariance/CovarianceSample.d.ts

@@ -0,0 +1,44 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Epoch, EpochUTC, J2000 } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { Thrust } from '../force/Thrust';
+import { StateCovariance } from './StateCovariance';
+export declare class CovarianceSample {
+    private _origin;
+    private _samples;
+    private _pts;
+    /**
+     * Create a new [CovarianceSample] object from an inertial state, covariance
+     * and optional force models for the origin state and samples.
+     *
+     * Two-body physics will be used if a force model is not provided.
+     * @param state
+     * @param covariance
+     * @param originForceModel
+     * @param sampleForceModel
+     */
+    constructor(state: J2000, covariance: StateCovariance, originForceModel?: ForceModel, sampleForceModel?: ForceModel);
+    get epoch(): Epoch;
+    get state(): J2000;
+    private _rebuildCovariance;
+    propagate(epoch: EpochUTC): void;
+    maneuver(maneuver: Thrust): void;
+    desampleJ2000(): StateCovariance;
+    desampleRIC(): StateCovariance;
+}

package/lib/covariance/CovarianceSample.js

@@ -0,0 +1,146 @@
+/**
+ * @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 { J2000, Matrix, RelativeState, RIC, Vector, Vector3D } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator';
+import { CovarianceFrame, StateCovariance } from './StateCovariance';
+// / Sigma point covariance sample.
+export class CovarianceSample {
+    /**
+     * Create a new [CovarianceSample] object from an inertial state, covariance
+     * and optional force models for the origin state and samples.
+     *
+     * Two-body physics will be used if a force model is not provided.
+     * @param state
+     * @param covariance
+     * @param originForceModel
+     * @param sampleForceModel
+     */
+    constructor(state, covariance, originForceModel, sampleForceModel) {
+        this._samples = [];
+        this._pts = Matrix.zero(6, 12);
+        originForceModel ??= new ForceModel().setGravity();
+        sampleForceModel ??= new ForceModel().setGravity();
+        this._origin = new RungeKutta89Propagator(state, originForceModel);
+        const s = covariance.matrix.cholesky();
+        const sqrt6 = Math.sqrt(6.0);
+        for (let i = 0; i < 6; i++) {
+            for (let j = 0; j < 6; j++) {
+                s[i][j] *= sqrt6;
+            }
+        }
+        // 6 x 12 matrix
+        const sigmapts = Matrix.zero(6, 12);
+        for (let i = 0; i < 6; i++) {
+            const jj = (i - 1) * 2 + 2;
+            for (let j = 0; j < 3; j++) {
+                sigmapts[j][jj] = s[j][i];
+                sigmapts[j + 3][jj] = s[j + 3][i];
+                sigmapts[j][jj + 1] = -s[j][i];
+                sigmapts[j + 3][jj + 1] = -s[j + 3][i];
+            }
+        }
+        for (let i = 0; i < 12; i++) {
+            const sampleR = new Vector3D(sigmapts[0][i], sigmapts[1][i], sigmapts[2][i]);
+            const sampleV = new Vector3D(sigmapts[3][i], sigmapts[4][i], sigmapts[5][i]);
+            if (covariance.frame === CovarianceFrame.ECI) {
+                const sample = new J2000(state.epoch, state.position.add(sampleR), state.velocity.add(sampleV));
+                this._samples.push(new RungeKutta89Propagator(sample, sampleForceModel));
+            }
+            else if (covariance.frame === CovarianceFrame.RIC) {
+                const sample = new RIC(sampleR, sampleV).toJ2000(state);
+                this._samples.push(new RungeKutta89Propagator(sample, sampleForceModel));
+            }
+        }
+    }
+    // / Current covariance sample epoch.
+    get epoch() {
+        return this._origin.state.epoch;
+    }
+    // / Current covariance sample origin state.
+    get state() {
+        return this._origin.state;
+    }
+    // / Rebuild covariance from sigma points.
+    _rebuildCovariance(pts) {
+        const c = 1.0 / 12.0;
+        const yu = new Vector([0, 0, 0, 0, 0, 0]);
+        const y = Matrix.zero(6, 12);
+        for (let i = 0; i < 12; i++) {
+            for (let j = 0; j < 6; j++) {
+                yu[j] += pts[j][i];
+            }
+        }
+        for (let j = 0; j < 6; j++) {
+            yu[j] *= c;
+        }
+        for (let i = 0; i < 12; i++) {
+            for (let j = 0; j < 6; j++) {
+                y[j][i] = pts[j][i] - yu[j];
+            }
+        }
+        const yt = y.transpose();
+        const tmp = y.multiply(yt);
+        return tmp.scale(c);
+    }
+    // / Propagate covariance to a new epoch.
+    propagate(epoch) {
+        this._origin.propagate(epoch);
+        for (const sample of this._samples) {
+            sample.propagate(epoch);
+        }
+    }
+    // / Apply a maneuver to this covariance.
+    maneuver(maneuver) {
+        this._origin.maneuver(maneuver);
+        for (const sample of this._samples) {
+            sample.maneuver(maneuver);
+        }
+    }
+    // / Desample covariance in J2000 frame.
+    desampleJ2000() {
+        for (let i = 0; i < 12; i++) {
+            const state = this._samples[i].state;
+            this._pts[0][i] = state.position[0];
+            this._pts[1][i] = state.position[1];
+            this._pts[2][i] = state.position[2];
+            this._pts[3][i] = state.velocity[0];
+            this._pts[4][i] = state.velocity[1];
+            this._pts[5][i] = state.velocity[2];
+        }
+        const matrix = this._rebuildCovariance(this._pts);
+        return new StateCovariance(matrix, CovarianceFrame.ECI);
+    }
+    // / Desample covariance in RIC frame.
+    desampleRIC() {
+        const rot = RelativeState.createMatrix(this._origin.state.position, this._origin.state.velocity);
+        for (let i = 0; i < 12; i++) {
+            const state = RIC.fromJ2000Matrix(this._samples[i].state, this._origin.state, rot);
+            this._pts[0][i] = state.position[0];
+            this._pts[1][i] = state.position[1];
+            this._pts[2][i] = state.position[2];
+            this._pts[3][i] = state.velocity[0];
+            this._pts[4][i] = state.velocity[1];
+            this._pts[5][i] = state.velocity[2];
+        }
+        const matrix = this._rebuildCovariance(this._pts);
+        return new StateCovariance(matrix, CovarianceFrame.RIC);
+    }
+}
+//# sourceMappingURL=CovarianceSample.js.map

package/lib/covariance/CovarianceSample.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CovarianceSample.js","sourceRoot":"","sources":["../../src/covariance/CovarianceSample.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,2CAA2C;AAC3C,OAAO,EAAmB,KAAK,EAAE,MAAM,EAAE,aAAa,EAAE,GAAG,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AACjG,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AAEjD,OAAO,EAAE,sBAAsB,EAAE,MAAM,sCAAsC,CAAC;AAC9E,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,mBAAmB,CAAC;AAErE,mCAAmC;AACnC,MAAM,OAAO,gBAAgB;IAK3B;;;;;;;;;OASG;IACH,YAAY,KAAY,EAAE,UAA2B,EAAE,gBAA6B,EAAE,gBAA6B;QAb3G,aAAQ,GAA6B,EAAE,CAAC;QACxC,SAAI,GAAW,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAaxC,gBAAgB,KAAK,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;QACnD,gBAAgB,KAAK,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;QACnD,IAAI,CAAC,OAAO,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,gBAAgB,CAAC,CAAC;QACnE,MAAM,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,QAAQ,EAAE,CAAC;QACvC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAE7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC;aAClB;SACF;QAED,gBAAgB;QAChB,MAAM,QAAQ,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAEpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAE3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAElC,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC/B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACxC;SACF;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;YAC3B,MAAM,OAAO,GAAG,IAAI,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAC7E,MAAM,OAAO,GAAG,IAAI,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAE7E,IAAI,UAAU,CAAC,KAAK,KAAK,eAAe,CAAC,GAAG,EAAE;gBAC5C,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC;gBAEhG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,sBAAsB,CAAC,MAAM,EAAE,gBAAgB,CAAC,CAAC,CAAC;aAC1E;iBAAM,IAAI,UAAU,CAAC,KAAK,KAAK,eAAe,CAAC,GAAG,EAAE;gBACnD,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,OAAO,EAAE,OAAO,CAAC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC;gBAExD,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,sBAAsB,CAAC,MAAM,EAAE,gBAAgB,CAAC,CAAC,CAAC;aAC1E;SACF;IACH,CAAC;IAED,qCAAqC;IACrC,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,CAAC;IAClC,CAAC;IAED,4CAA4C;IAC5C,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC;IAC5B,CAAC;IAED,0CAA0C;IAClC,kBAAkB,CAAC,GAAW;QACpC,MAAM,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC;QACrB,MAAM,EAAE,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC1C,MAAM,CAAC,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAE7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,EAAE,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACpB;SACF;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC;SACZ;QAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;aAC7B;SACF;QAED,MAAM,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QACzB,MAAM,GAAG,GAAG,CAAC,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;QAE3B,OAAO,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IACtB,CAAC;IAED,yCAAyC;IACzC,SAAS,CAAC,KAAe;QACvB,IAAI,CAAC,OAAO,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;QAC9B,KAAK,MAAM,MAAM,IAAI,IAAI,CAAC,QAAQ,EAAE;YAClC,MAAM,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;SACzB;IACH,CAAC;IAED,yCAAyC;IACzC,QAAQ,CAAC,QAAgB;QACvB,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;QAChC,KAAK,MAAM,MAAM,IAAI,IAAI,CAAC,QAAQ,EAAE;YAClC,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;SAC3B;IACH,CAAC;IAED,wCAAwC;IACxC,aAAa;QACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;YAC3B,MAAM,KAAK,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;YAErC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;SACrC;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAElD,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;IAC1D,CAAC;IAED,sCAAsC;IACtC,WAAW;QACT,MAAM,GAAG,GAAG,aAAa,CAAC,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC;QAEjG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE;YAC3B,MAAM,KAAK,GAAG,GAAG,CAAC,eAAe,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,IAAI,CAAC,OAAO,CAAC,KAAK,EAAE,GAAG,CAAC,CAAC;YAEnF,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;SACrC;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAElD,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;IAC1D,CAAC;CACF"}

package/lib/covariance/StateCovariance.d.ts

@@ -0,0 +1,41 @@
+/**
+ * @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 { Matrix, Vector } from 'ootk-core';
+/** Covariance Frame */
+export declare enum CovarianceFrame {
+    /** Earth-centered inertial */
+    ECI = "eci",
+    /** Radial-Intrack-Crosstrack */
+    RIC = "ric"
+}
+export declare class StateCovariance {
+    matrix: Matrix;
+    frame: CovarianceFrame;
+    /**
+     * Create a new [StateCovariance] object given its covariance [matrix] and
+     * [CovarianceFrame].
+     */
+    constructor(matrix: Matrix, frame: CovarianceFrame);
+    static fromSigmas(sigmas: number[], frame: CovarianceFrame): StateCovariance;
+    /**
+     * Calculates the standard deviations (sigmas) of each element in the covariance matrix.
+     *
+     * @returns A vector containing the standard deviations of each element in the covariance matrix.
+     */
+    sigmas(): Vector;
+}

package/lib/covariance/StateCovariance.js

@@ -0,0 +1,62 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Matrix, Vector } from 'ootk-core';
+/** Covariance Frame */
+export var CovarianceFrame;
+(function (CovarianceFrame) {
+    /** Earth-centered inertial */
+    CovarianceFrame["ECI"] = "eci";
+    /** Radial-Intrack-Crosstrack */
+    CovarianceFrame["RIC"] = "ric";
+})(CovarianceFrame || (CovarianceFrame = {}));
+// / State covariance.
+export class StateCovariance {
+    /**
+     * Create a new [StateCovariance] object given its covariance [matrix] and
+     * [CovarianceFrame].
+     */
+    constructor(matrix, frame) {
+        this.matrix = matrix;
+        this.frame = frame;
+        // Nothing to do here.
+    }
+    // / Create a new [StateCovariance] object from 1-sigma values.
+    static fromSigmas(sigmas, frame) {
+        const n = sigmas.length;
+        const output = Matrix.zero(n, n);
+        for (let i = 0; i < n; i++) {
+            output[i][i] = Math.max(sigmas[i] * sigmas[i], 1e-32);
+        }
+        return new StateCovariance(output, frame);
+    }
+    /**
+     * Calculates the standard deviations (sigmas) of each element in the covariance matrix.
+     *
+     * @returns A vector containing the standard deviations of each element in the covariance matrix.
+     */
+    sigmas() {
+        const c = this.matrix.columns;
+        const result = new Float64Array(c);
+        for (let i = 0; i < c; i++) {
+            const variance = this.matrix[i][i];
+            result[i] = Math.sqrt(variance);
+        }
+        return new Vector(result);
+    }
+}
+//# sourceMappingURL=StateCovariance.js.map

package/lib/covariance/StateCovariance.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"StateCovariance.js","sourceRoot":"","sources":["../../src/covariance/StateCovariance.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAE3C,uBAAuB;AACvB,MAAM,CAAN,IAAY,eAKX;AALD,WAAY,eAAe;IACzB,8BAA8B;IAC9B,8BAAW,CAAA;IACX,gCAAgC;IAChC,8BAAW,CAAA;AACb,CAAC,EALW,eAAe,KAAf,eAAe,QAK1B;AAED,sBAAsB;AACtB,MAAM,OAAO,eAAe;IAC1B;;;OAGG;IACH,YAAmB,MAAc,EAAS,KAAsB;QAA7C,WAAM,GAAN,MAAM,CAAQ;QAAS,UAAK,GAAL,KAAK,CAAiB;QAC9D,sBAAsB;IACxB,CAAC;IAED,+DAA+D;IAC/D,MAAM,CAAC,UAAU,CAAC,MAAgB,EAAE,KAAsB;QACxD,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;SACvD;QAED,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;IAC5C,CAAC;IAED;;;;OAIG;IACH,MAAM;QACJ,MAAM,CAAC,GAAG,IAAI,CAAC,MAAM,CAAC,OAAO,CAAC;QAC9B,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAEnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,QAAQ,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEnC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SACjC;QAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;CACF"}

package/lib/force/AtmosphericDrag.d.ts

@@ -0,0 +1,32 @@
+/**
+ * @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, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+/**
+ * Harris-Priester atmospheric drag force model.
+ * Atmospheric density model assumes mean solar flux.
+ */
+export declare class AtmosphericDrag implements Force {
+    mass: number;
+    area: number;
+    dragCoeff: number;
+    cosine: number;
+    constructor(mass: number, area: number, dragCoeff: number, cosine: number);
+    private static _getHPDensity;
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/AtmosphericDrag.js

@@ -0,0 +1,66 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { DataHandler, DEG2RAD, Earth, ITRF, J2000, Sun, Vector3D } from 'ootk-core';
+/**
+ * Harris-Priester atmospheric drag force model.
+ * Atmospheric density model assumes mean solar flux.
+ */
+export class AtmosphericDrag {
+    constructor(mass, area, dragCoeff, cosine) {
+        this.mass = mass;
+        this.area = area;
+        this.dragCoeff = dragCoeff;
+        this.cosine = cosine;
+    }
+    static _getHPDensity(state, n) {
+        const hpa = DataHandler.getInstance().getHpAtmosphere(state.getHeight());
+        if (hpa === null) {
+            return 0.0;
+        }
+        const sunPos = Sun.positionApparent(state.epoch);
+        const sunVec = new J2000(state.epoch, sunPos, Vector3D.origin).toITRF().position.normalize();
+        const bulVec = sunVec.rotZ(-30.0 * DEG2RAD);
+        const cosPsi = bulVec.normalize().dot(state.position.normalize());
+        const c2Psi2 = 0.5 * (1.0 + cosPsi);
+        const cPsi2 = Math.sqrt(c2Psi2);
+        const cosPow = cPsi2 > 1e-12 ? c2Psi2 * cPsi2 ** (n - 2) : 0.0;
+        const altitude = hpa.height;
+        const [h0, min0, max0] = hpa.hp0;
+        const [h1, min1, max1] = hpa.hp1;
+        const dH = (h0 - altitude) / (h0 - h1);
+        const rhoMin = min0 * (min1 / min0) ** dH;
+        if (cosPow === 0) {
+            return rhoMin;
+        }
+        const rhoMax = max0 * (max1 / max0) ** dH;
+        return rhoMin + (rhoMax - rhoMin) * cosPow;
+    }
+    acceleration(state) {
+        const itrfState = state.toITRF();
+        const density = AtmosphericDrag._getHPDensity(itrfState, this.cosine);
+        if (density === 0) {
+            return Vector3D.origin;
+        }
+        const rotation = new ITRF(state.epoch, Earth.rotation, Vector3D.origin).toJ2000().position;
+        const vRel = state.velocity.subtract(rotation.cross(state.position)).scale(1000.0);
+        const vm = vRel.magnitude();
+        const fScale = -0.5 * density * ((this.dragCoeff * this.area) / this.mass) * vm;
+        return vRel.scale(fScale / 1000.0);
+    }
+}
+//# sourceMappingURL=AtmosphericDrag.js.map

package/lib/force/AtmosphericDrag.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"AtmosphericDrag.js","sourceRoot":"","sources":["../../src/force/AtmosphericDrag.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,KAAK,EAAE,IAAI,EAAE,KAAK,EAAE,GAAG,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAGpF;;;GAGG;AACH,MAAM,OAAO,eAAe;IAM1B,YAAY,IAAY,EAAE,IAAY,EAAE,SAAiB,EAAE,MAAc;QACvE,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,SAAS,GAAG,SAAS,CAAC;QAC3B,IAAI,CAAC,MAAM,GAAG,MAAM,CAAC;IACvB,CAAC;IAEO,MAAM,CAAC,aAAa,CAAC,KAAW,EAAE,CAAS;QACjD,MAAM,GAAG,GAAG,WAAW,CAAC,WAAW,EAAE,CAAC,eAAe,CAAC,KAAK,CAAC,SAAS,EAAE,CAAC,CAAC;QAEzE,IAAI,GAAG,KAAK,IAAI,EAAE;YAChB,OAAO,GAAG,CAAC;SACZ;QACD,MAAM,MAAM,GAAG,GAAG,CAAC,gBAAgB,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACjD,MAAM,MAAM,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC,KAAK,EAAE,MAAM,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAC7F,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC,CAAC;QAC5C,MAAM,MAAM,GAAG,MAAM,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC,CAAC;QAClE,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC;QACpC,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAChC,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,GAAG,KAAK,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC;QAC/D,MAAM,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC;QAC5B,MAAM,CAAC,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC;QACjC,MAAM,CAAC,EAAE,EAAE,IAAI,EAAE,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC;QACjC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACvC,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;QAE1C,IAAI,MAAM,KAAK,CAAC,EAAE;YAChB,OAAO,MAAM,CAAC;SACf;QACD,MAAM,MAAM,GAAG,IAAI,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC;QAE1C,OAAO,MAAM,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC,GAAG,MAAM,CAAC;IAC7C,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,MAAM,SAAS,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC;QACjC,MAAM,OAAO,GAAG,eAAe,CAAC,aAAa,CAAC,SAAS,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QAEtE,IAAI,OAAO,KAAK,CAAC,EAAE;YACjB,OAAO,QAAQ,CAAC,MAAM,CAAC;SACxB;QACD,MAAM,QAAQ,GAAG,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,QAAQ,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC;QAC3F,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,KAAK,CAAC,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;QACnF,MAAM,EAAE,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QAC5B,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,OAAO,GAAG,CAAC,CAAC,IAAI,CAAC,SAAS,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,EAAE,CAAC;QAEhF,OAAO,IAAI,CAAC,KAAK,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC;IACrC,CAAC;CACF"}

package/lib/force/EarthGravity.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 { J2000, Vector3D } from 'ootk-core';
+import { Force } from './Force';
+export declare class EarthGravity implements Force {
+    degree: number;
+    order: number;
+    _asphericalFlag: boolean;
+    constructor(degree: number, order: number);
+    _spherical(state: J2000): Vector3D;
+    _aspherical(state: J2000): Vector3D;
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/EarthGravity.js

@@ -0,0 +1,115 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { DataHandler, Earth, ITRF, Vector3D } from 'ootk-core';
+// / tesseral geopotential perturbations.
+export class EarthGravity {
+    constructor(degree, order) {
+        this.degree = Math.min(Math.max(degree, 0), 36);
+        this.order = Math.min(Math.max(order, 0), 36);
+        this._asphericalFlag = degree >= 2;
+    }
+    _spherical(state) {
+        const rMag = state.position.magnitude();
+        return state.position.scale(-Earth.mu / (rMag * rMag * rMag));
+    }
+    // eslint-disable-next-line max-statements
+    _aspherical(state) {
+        const posEcef = state.toITRF().position;
+        const ri = 1.0 / posEcef.magnitude();
+        const xor = posEcef.x * ri;
+        const yor = posEcef.y * ri;
+        const zor = posEcef.z * ri;
+        const ep = zor;
+        const reor = Earth.radiusEquator * ri;
+        let reorn = reor;
+        const muor2 = Earth.mu * ri * ri;
+        let sumH = 0.0;
+        let sumGm = 0.0;
+        let sumJ = 0.0;
+        let sumK = 0.0;
+        const cTil = new Float64Array(this.order + 4);
+        const sTil = new Float64Array(this.order + 4);
+        const pN = new Float64Array(this.order + 4);
+        const pNm1 = new Float64Array(this.order + 4);
+        const pNm2 = new Float64Array(this.order + 4);
+        pNm2[0] = 1.0;
+        pNm1[0] = ep;
+        pNm1[1] = 1.0;
+        cTil[0] = 1.0;
+        cTil[1] = xor;
+        sTil[1] = yor;
+        const dh = DataHandler.getInstance();
+        for (let n = 2, nm1 = 1, nm2 = 0, np1 = 3; n <= this.degree; nm2++, nm1++, n++, np1++) {
+            const twonm1 = 2.0 * n - 1.0;
+            reorn *= reor;
+            const cN0 = dh.getEgm96Coeffs(n, 0)[2];
+            pN[0] = (twonm1 * ep * pNm1[0] - nm1 * pNm2[0]) / n;
+            pN[1] = pNm2[1] + twonm1 * pNm1[0];
+            pN[2] = pNm2[2] + twonm1 * pNm1[1];
+            let sumHn = pN[1] * cN0;
+            let sumGmn = pN[0] * cN0 * np1;
+            if (this.order > 0) {
+                let sumJn = 0.0;
+                let sumKn = 0.0;
+                cTil[n] = cTil[1] * cTil[nm1] - sTil[1] * sTil[nm1];
+                sTil[n] = sTil[1] * cTil[nm1] + cTil[1] * sTil[nm1];
+                const lim = n < this.order ? n : this.order;
+                for (let m = 1, mm1 = 0, mm2 = -1, mp1 = 2, mp2 = 3; m <= lim; mm2++, mm1++, m++, mp1++, mp2++) {
+                    pN[mp1] = pNm2[mp1] + twonm1 * pNm1[m];
+                    const dm = m;
+                    const npmp1 = n + mp1;
+                    const pNm = pN[m];
+                    const pNmp1 = pN[mp1];
+                    const coefs = dh.getEgm96Coeffs(n, m);
+                    const cNm = coefs[2];
+                    const sNm = coefs[3];
+                    const mxPnm = dm * pNm;
+                    const bNmtil = cNm * cTil[m] + sNm * sTil[m];
+                    const pNmBnm = pNm * bNmtil;
+                    const bNmtm1 = cNm * cTil[mm1] + sNm * sTil[mm1];
+                    const aNmtm1 = cNm * sTil[mm1] - sNm * cTil[mm1];
+                    sumHn += pNmp1 * bNmtil;
+                    sumGmn += npmp1 * pNmBnm;
+                    sumJn += mxPnm * bNmtm1;
+                    sumKn -= mxPnm * aNmtm1;
+                }
+                sumJ += reorn * sumJn;
+                sumK += reorn * sumKn;
+            }
+            sumH += reorn * sumHn;
+            sumGm += reorn * sumGmn;
+            if (n < this.degree) {
+                for (let i = 0; i <= n; i++) {
+                    pNm2[i] = pNm1[i];
+                    pNm1[i] = pN[i];
+                }
+            }
+        }
+        const lambda = sumGm + ep * sumH;
+        const g = new Vector3D(-muor2 * (lambda * xor - sumJ), -muor2 * (lambda * yor - sumK), -muor2 * (lambda * zor - sumH));
+        return new ITRF(state.epoch, g, Vector3D.origin).toJ2000().position;
+    }
+    acceleration(state) {
+        let accVec = this._spherical(state);
+        if (this._asphericalFlag) {
+            accVec = accVec.add(this._aspherical(state));
+        }
+        return accVec;
+    }
+}
+//# sourceMappingURL=EarthGravity.js.map

package/lib/force/EarthGravity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"EarthGravity.js","sourceRoot":"","sources":["../../src/force/EarthGravity.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,WAAW,EAAE,KAAK,EAAE,IAAI,EAAS,QAAQ,EAAE,MAAM,WAAW,CAAC;AAMtE,yCAAyC;AACzC,MAAM,OAAO,YAAY;IAKvB,YAAY,MAAc,EAAE,KAAa;QACvC,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAChD,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC9C,IAAI,CAAC,eAAe,GAAG,MAAM,IAAI,CAAC,CAAC;IACrC,CAAC;IAED,UAAU,CAAC,KAAY;QACrB,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAExC,OAAO,KAAK,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC;IAChE,CAAC;IAED,0CAA0C;IAC1C,WAAW,CAAC,KAAY;QACtB,MAAM,OAAO,GAAG,KAAK,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC;QACxC,MAAM,EAAE,GAAG,GAAG,GAAG,OAAO,CAAC,SAAS,EAAE,CAAC;QACrC,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,OAAO,CAAC,CAAC,GAAG,EAAE,CAAC;QAE3B,MAAM,EAAE,GAAG,GAAG,CAAC;QACf,MAAM,IAAI,GAAG,KAAK,CAAC,aAAa,GAAG,EAAE,CAAC;QACtC,IAAI,KAAK,GAAG,IAAI,CAAC;QACjB,MAAM,KAAK,GAAG,KAAK,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;QAEjC,IAAI,IAAI,GAAG,GAAG,CAAC;QACf,IAAI,KAAK,GAAG,GAAG,CAAC;QAChB,IAAI,IAAI,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,GAAG,GAAG,CAAC;QAEf,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAE9C,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC5C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAC9C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QAE9C,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QACd,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;QAEd,MAAM,EAAE,GAAG,WAAW,CAAC,WAAW,EAAE,CAAC;QAErC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,MAAM,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE;YACrF,MAAM,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,GAAG,CAAC;YAE7B,KAAK,IAAI,IAAI,CAAC;YACd,MAAM,GAAG,GAAG,EAAE,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAEvC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACpD,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACnC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YAEnC,IAAI,KAAK,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YACxB,IAAI,MAAM,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;YAE/B,IAAI,IAAI,CAAC,KAAK,GAAG,CAAC,EAAE;gBAClB,IAAI,KAAK,GAAG,GAAG,CAAC;gBAChB,IAAI,KAAK,GAAG,GAAG,CAAC;gBAEhB,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBACpD,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;gBAEpD,MAAM,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC;gBAE5C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE;oBAC9F,EAAE,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAEvC,MAAM,EAAE,GAAG,CAAC,CAAC;oBACb,MAAM,KAAK,GAAG,CAAC,GAAG,GAAG,CAAC;oBAEtB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;oBAClB,MAAM,KAAK,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;oBAEtB,MAAM,KAAK,GAAG,EAAE,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;oBACtC,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBACrB,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAErB,MAAM,KAAK,GAAG,EAAE,GAAG,GAAG,CAAC;oBACvB,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAC7C,MAAM,MAAM,GAAG,GAAG,GAAG,MAAM,CAAC;oBAC5B,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;oBACjD,MAAM,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;oBAEjD,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;oBACxB,MAAM,IAAI,KAAK,GAAG,MAAM,CAAC;oBACzB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;oBACxB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;iBACzB;gBAED,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;gBACtB,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;aACvB;YAED,IAAI,IAAI,KAAK,GAAG,KAAK,CAAC;YACtB,KAAK,IAAI,KAAK,GAAG,MAAM,CAAC;YAExB,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE;gBACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,EAAE,EAAE;oBAC3B,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBAClB,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;iBACjB;aACF;SACF;QAED,MAAM,MAAM,GAAG,KAAK,GAAG,EAAE,GAAG,IAAI,CAAC;QACjC,MAAM,CAAC,GAAG,IAAI,QAAQ,CACpB,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,EAC9B,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,EAC9B,CAAC,KAAK,GAAG,CAAC,MAAM,GAAG,GAAG,GAAG,IAAI,CAAC,CAC/B,CAAC;QAEF,OAAO,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC;IACtE,CAAC;IAED,YAAY,CAAC,KAAY;QACvB,IAAI,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,KAAK,CAAC,CAAC;QAEpC,IAAI,IAAI,CAAC,eAAe,EAAE;YACxB,MAAM,GAAG,MAAM,CAAC,GAAG,CAAC,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC,CAAC;SAC9C;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}

package/{dist/ootk.d.ts → lib/force/Force.d.ts}

@@ -1,10 +1,8 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
+ * @author @thkruz Theodore Kruczek
  *
  * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 Theodore Kruczek
+ * @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
@@ -17,9 +15,11 @@
  * 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/>.
  */
-export { Sgp4 } from './sgp4/sgp4';
-export { Transforms } from './transforms/transforms';
-export { Utils } from './utils/utils';
-export { Tle } from './tle/tle';
-export * from './objects/index';
-export * from './types/types';
+import { J2000, Vector3D } from 'ootk-core';
+export declare abstract class Force {
+    /**
+     * Calculate the acceleration due to the perturbing force on a given
+     * state vector.
+     */
+    acceleration(state: J2000): Vector3D;
+}

package/lib/force/Force.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/>.
+ */
+/* eslint-disable @typescript-eslint/no-unused-vars */
+/* eslint-disable class-methods-use-this */
+// / Base class for perturbation forces.
+export class Force {
+    /**
+     * Calculate the acceleration due to the perturbing force on a given
+     * state vector.
+     */
+    acceleration(state) {
+        throw Error('Not implemented');
+    }
+}
+//# sourceMappingURL=Force.js.map

package/lib/force/Force.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Force.js","sourceRoot":"","sources":["../../src/force/Force.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAGH,sDAAsD;AACtD,2CAA2C;AAE3C,wCAAwC;AACxC,MAAM,OAAgB,KAAK;IACzB;;;OAGG;IACH,YAAY,CAAC,KAAY;QACvB,MAAM,KAAK,CAAC,iBAAiB,CAAC,CAAC;IACjC,CAAC;CACF"}