ootk 3.0.0 → 4.0.0-1

package/dist/coordinate/Hill.js

@@ -0,0 +1,222 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Earth, J2000, Matrix, RelativeState, Vector3D, } from 'ootk-core';
+import { Thrust } from './../force/Thrust';
+// / Hill Modified Equidistant Cyllindrical _(EQCM)_ coordinates.
+export class Hill {
+    epoch;
+    position;
+    velocity;
+    semimajorAxis_;
+    meanMotion_;
+    constructor(epoch, position, velocity, semimajorAxis) {
+        this.epoch = epoch;
+        this.position = position;
+        this.velocity = velocity;
+        this.semimajorAxis_ = semimajorAxis;
+        this.meanMotion_ = Earth.smaToMeanMotion(this.semimajorAxis_);
+    }
+    static fromState(origin, radialPosition, intrackPosition, nodeVelocity, nodeOffsetTime) {
+        const a = origin.semimajorAxis;
+        const n = Earth.smaToMeanMotion(a);
+        const yDot = -3.0 * radialPosition * n * 0.5;
+        const z = (nodeVelocity / n) * Math.sin(n * -nodeOffsetTime);
+        const zDot = nodeVelocity * Math.cos(n * -nodeOffsetTime);
+        const r = new Vector3D(radialPosition, intrackPosition, z);
+        const v = new Vector3D(0.0, yDot, zDot);
+        return new Hill(origin.epoch, r, v, a);
+    }
+    static fromNmc(origin, majorAxisRange, nodeVelocity, nodeOffsetTime, translation = 0.0) {
+        const a = origin.semimajorAxis;
+        const n = Earth.smaToMeanMotion(a);
+        const xDot = majorAxisRange * n * 0.5;
+        const z = (nodeVelocity / n) * Math.sin(n * -nodeOffsetTime);
+        const zDot = nodeVelocity * Math.cos(n * -nodeOffsetTime);
+        const r = new Vector3D(0.0, majorAxisRange + translation, z);
+        const v = new Vector3D(xDot, 0.0, zDot);
+        return new Hill(origin.epoch, r, v, a);
+    }
+    static fromPerch(origin, perchRange, nodeVelocity, nodeOffsetTime) {
+        const a = origin.semimajorAxis;
+        const n = Earth.smaToMeanMotion(a);
+        const z = (nodeVelocity / n) * Math.sin(n * -nodeOffsetTime);
+        const zDot = nodeVelocity * Math.cos(n * -nodeOffsetTime);
+        const r = new Vector3D(0.0, perchRange, z);
+        const v = new Vector3D(0.0, 0.0, zDot);
+        return new Hill(origin.epoch, r, v, a);
+    }
+    get semimajorAxis() {
+        return this.semimajorAxis_;
+    }
+    set semimajorAxis(sma) {
+        this.semimajorAxis_ = sma;
+        this.meanMotion_ = Earth.smaToMeanMotion(this.semimajorAxis_);
+    }
+    get meanMotion() {
+        return this.meanMotion_;
+    }
+    toJ2000Matrix(origin, transform) {
+        const magrtgt = origin.position.magnitude();
+        const magrint = magrtgt + this.position.x;
+        const vtgtrsw = transform.multiplyVector3D(origin.velocity);
+        const lambdadottgt = vtgtrsw.y / magrtgt;
+        const lambdaint = this.position.y / magrtgt;
+        const phiint = this.position.z / magrtgt;
+        const sinphiint = Math.sin(phiint);
+        const cosphiint = Math.cos(phiint);
+        const sinlambdaint = Math.sin(lambdaint);
+        const coslambdaint = Math.cos(lambdaint);
+        const rotRswSez = new Matrix([
+            [sinphiint * coslambdaint, sinphiint * sinlambdaint, -cosphiint],
+            [-sinlambdaint, coslambdaint, 0],
+            [cosphiint * coslambdaint, cosphiint * sinlambdaint, sinphiint],
+        ]);
+        const rdotint = this.velocity.x + vtgtrsw.x;
+        const lambdadotint = this.velocity.y / magrtgt + lambdadottgt;
+        const phidotint = this.velocity.z / magrtgt;
+        const vintsez = new Vector3D(-magrint * phidotint, magrint * lambdadotint * cosphiint, rdotint);
+        const vintrsw = rotRswSez.transpose().multiplyVector3D(vintsez);
+        const vinteci = transform.transpose().multiplyVector3D(vintrsw);
+        const rintrsw = new Vector3D(cosphiint * magrint * coslambdaint, cosphiint * magrint * sinlambdaint, sinphiint * magrint);
+        const rinteci = transform.transpose().multiplyVector3D(rintrsw);
+        return new J2000(origin.epoch, rinteci, vinteci);
+    }
+    toJ2000(origin) {
+        return this.toJ2000Matrix(origin, RelativeState.createMatrix(origin.position, origin.velocity));
+    }
+    static transitionMatrix(t, meanMotion) {
+        const n = meanMotion;
+        const sn = Math.sin(n * t);
+        const cs = Math.cos(n * t);
+        return new Matrix([
+            [4.0 - 3.0 * cs, 0.0, 0.0, sn / n, (2.0 * (1.0 - cs)) / n, 0.0],
+            [6.0 * (sn - n * t), 1.0, 0.0, (-2.0 * (1.0 - cs)) / n, (4.0 * sn - 3.0 * n * t) / n, 0.0],
+            [0.0, 0.0, cs, 0.0, 0.0, sn / n],
+            [3.0 * n * sn, 0.0, 0.0, cs, 2.0 * sn, 0.0],
+            [-6.0 * n * (1.0 - cs), 0.0, 0.0, -2.0 * sn, 4.0 * cs - 3.0, 0.0],
+            [0.0, 0.0, -n * sn, 0.0, 0.0, cs],
+        ]);
+    }
+    transition(t) {
+        const sysMat = Hill.transitionMatrix(t, this.meanMotion_);
+        const res = sysMat.multiplyVector(this.position.join(this.velocity)).elements;
+        return new Hill(this.epoch.roll(t), new Vector3D(res[0], res[1], res[2]), new Vector3D(res[3], res[4], res[5]), this.semimajorAxis_);
+    }
+    transitionWithMatrix(stm, t) {
+        const res = stm.multiplyVector(this.position.join(this.velocity)).elements;
+        return new Hill(this.epoch.roll(t), new Vector3D(res[0], res[1], res[2]), new Vector3D(res[3], res[4], res[5]), this.semimajorAxis_);
+    }
+    propagate(newEpoch) {
+        return this.transition(newEpoch.difference(this.epoch));
+    }
+    propagateWithMatrix(stm, newEpoch) {
+        return this.transitionWithMatrix(stm, newEpoch.difference(this.epoch));
+    }
+    maneuver(maneuver) {
+        const state = this.propagate(maneuver.center);
+        return new Hill(state.epoch, state.position, state.velocity.add(maneuver.deltaV), state.semimajorAxis_);
+    }
+    ephemeris(start, stop, step = 60.0) {
+        const output = [];
+        let current = start;
+        while (stop >= current) {
+            output.push(this.propagate(current));
+            current = current.roll(step);
+        }
+        return output;
+    }
+    get period() {
+        return (2 * Math.PI) / this.meanMotion_;
+    }
+    nextRadialTangent() {
+        const x = this.position.x;
+        const xDot = this.velocity.x;
+        const yDot = this.velocity.y;
+        let t = Math.atan(-xDot / (3.0 * this.meanMotion_ * x + 2.0 * yDot)) / this.meanMotion_;
+        if (t <= 0) {
+            t = t + 0.5 * this.period;
+        }
+        else if (isNaN(t)) {
+            t = 0.5 * this.period;
+        }
+        return this.propagate(this.epoch.roll(t));
+    }
+    solveManeuver(waypoint, ignoreCrosstrack = false) {
+        const t = waypoint.epoch.difference(this.epoch);
+        const w = waypoint.relativePosition;
+        const sysMat = Hill.transitionMatrix(t, this.meanMotion_);
+        const posEquationMat = new Matrix([
+            [sysMat.elements[0][0], sysMat.elements[0][1], sysMat.elements[0][2]],
+            [sysMat.elements[1][0], sysMat.elements[1][1], sysMat.elements[1][2]],
+            [sysMat.elements[2][0], sysMat.elements[2][1], sysMat.elements[2][2]],
+        ]);
+        const solnVector = w
+            .subtract(posEquationMat.multiplyVector3D(this.position));
+        const velEquationMat = new Matrix([
+            [sysMat.elements[0][3], sysMat.elements[0][4], sysMat.elements[0][5]],
+            [sysMat.elements[1][3], sysMat.elements[1][4], sysMat.elements[1][5]],
+            [sysMat.elements[2][3], sysMat.elements[2][4], sysMat.elements[2][5]],
+        ]);
+        let result = velEquationMat
+            .inverse()
+            .multiplyVector3D(solnVector)
+            .subtract(this.velocity);
+        if (ignoreCrosstrack) {
+            result = new Vector3D(result.x, result.y, 0);
+        }
+        return new Thrust(this.epoch, result.x * 1000, result.y * 1000, result.z * 1000);
+    }
+    maneuverSequence(pivot, waypoints, preManeuvers = [], postManeuvers = []) {
+        let state = new Hill(this.epoch, this.position, this.velocity, this.semimajorAxis_);
+        preManeuvers = preManeuvers.slice();
+        postManeuvers = postManeuvers.slice();
+        let output = preManeuvers;
+        // Note difference was once compareTo
+        output.sort((a, b) => a.center.difference(b.center));
+        output = output.filter((mvr) => mvr.center >= this.epoch && mvr.center >= pivot);
+        for (const mvr of output) {
+            state = state.maneuver(mvr);
+        }
+        state = state.propagate(pivot);
+        for (const wpt of waypoints) {
+            const mvr = state.solveManeuver(wpt);
+            state = state.maneuver(mvr);
+            output.push(mvr);
+        }
+        output.push(...postManeuvers);
+        return output;
+    }
+    maneuverOrigin(maneuver) {
+        const state = this.propagate(maneuver.center);
+        const vInit = Math.sqrt(Earth.mu / this.semimajorAxis_);
+        const vFinal = vInit - maneuver.intrack * 1e-3;
+        const aFinal = Earth.mu / (vFinal * vFinal);
+        return new Hill(state.epoch, state.position, state.velocity.subtract(maneuver.deltaV), aFinal);
+    }
+    get name() {
+        return 'Hill';
+    }
+    toString() {
+        return [
+            `[${this.name}]`,
+            `  Epoch: ${this.epoch}`,
+            `  Position: ${this.position.toString(6)} km`,
+            `  Velocity: ${this.velocity.toString(9)} km/s`,
+        ].join('\n');
+    }
+}
+//# sourceMappingURL=Hill.js.map

package/dist/coordinate/Hill.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Hill.js","sourceRoot":"","sources":["../../src/coordinate/Hill.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EACL,KAAK,EAEL,KAAK,EAGL,MAAM,EAGN,aAAa,EAEb,QAAQ,GACT,MAAM,WAAW,CAAC;AACnB,OAAO,EAAE,MAAM,EAAE,MAAM,mBAAmB,CAAC;AAG3C,iEAAiE;AACjE,MAAM,OAAO,IAAI;IAKN;IACA;IACA;IAND,cAAc,CAAa;IAC3B,WAAW,CAAmB;IAEtC,YACS,KAAe,EACf,QAA8B,EAC9B,QAAuC,EAC9C,aAAyB;QAHlB,UAAK,GAAL,KAAK,CAAU;QACf,aAAQ,GAAR,QAAQ,CAAsB;QAC9B,aAAQ,GAAR,QAAQ,CAA+B;QAG9C,IAAI,CAAC,cAAc,GAAG,aAAa,CAAC;QACpC,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,eAAe,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;IAChE,CAAC;IAED,MAAM,CAAC,SAAS,CACd,MAAa,EACb,cAA0B,EAC1B,eAA2B,EAC3B,YAAiC,EACjC,cAAuB;QAEvB,MAAM,CAAC,GAAG,MAAM,CAAC,aAAa,CAAC;QAC/B,MAAM,CAAC,GAAG,KAAK,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,CAAC,GAAG,GAAG,cAAc,GAAG,CAAC,GAAG,GAA0B,CAAC;QACpE,MAAM,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAe,CAAC;QAC3E,MAAM,IAAI,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAwB,CAAC;QACjF,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,cAAc,EAAE,eAAe,EAAE,CAAC,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,GAA0B,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAE/D,OAAO,IAAI,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IAED,MAAM,CAAC,OAAO,CACZ,MAAa,EACb,cAA0B,EAC1B,YAAiC,EACjC,cAAuB,EACvB,WAAW,GAAG,GAAG;QAEjB,MAAM,CAAC,GAAG,MAAM,CAAC,aAAa,CAAC;QAC/B,MAAM,CAAC,GAAG,KAAK,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,IAAI,GAAG,cAAc,GAAG,CAAC,GAAG,GAA0B,CAAC;QAC7D,MAAM,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAe,CAAC;QAC3E,MAAM,IAAI,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAwB,CAAC;QACjF,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,GAAiB,EAAE,cAAc,GAAG,WAAyB,EAAE,CAAC,CAAC,CAAC;QACzF,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,IAAI,EAAE,GAA0B,EAAE,IAAI,CAAC,CAAC;QAE/D,OAAO,IAAI,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IAED,MAAM,CAAC,SAAS,CACd,MAAa,EACb,UAAsB,EACtB,YAAiC,EACjC,cAAuB;QAEvB,MAAM,CAAC,GAAG,MAAM,CAAC,aAAa,CAAC;QAC/B,MAAM,CAAC,GAAG,KAAK,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,CAAC,YAAY,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAe,CAAC;QAC3E,MAAM,IAAI,GAAG,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,cAAc,CAAwB,CAAC;QACjF,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,GAAiB,EAAE,UAAU,EAAE,CAAC,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,IAAI,QAAQ,CAAC,GAA0B,EAAE,GAA0B,EAAE,IAAI,CAAC,CAAC;QAErF,OAAO,IAAI,IAAI,CAAC,MAAM,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IAED,IAAI,aAAa;QACf,OAAO,IAAI,CAAC,cAAc,CAAC;IAC7B,CAAC;IAED,IAAI,aAAa,CAAC,GAAe;QAC/B,IAAI,CAAC,cAAc,GAAG,GAAG,CAAC;QAC1B,IAAI,CAAC,WAAW,GAAG,KAAK,CAAC,eAAe,CAAC,IAAI,CAAC,cAAc,CAAC,CAAC;IAChE,CAAC;IAED,IAAI,UAAU;QACZ,OAAO,IAAI,CAAC,WAAW,CAAC;IAC1B,CAAC;IAED,aAAa,CAAC,MAAa,EAAE,SAAiB;QAC5C,MAAM,OAAO,GAAG,MAAM,CAAC,QAAQ,CAAC,SAAS,EAAE,CAAC;QAC5C,MAAM,OAAO,GAAG,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC1C,MAAM,OAAO,GAAG,SAAS,CAAC,gBAAgB,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC;QAE5D,MAAM,YAAY,GAAG,OAAO,CAAC,CAAC,GAAG,OAAO,CAAC;QACzC,MAAM,SAAS,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,CAAC;QAC5C,MAAM,MAAM,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,CAAC;QACzC,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC;QACnC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QACzC,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC;QAEzC,MAAM,SAAS,GAAG,IAAI,MAAM,CAAC;YAC3B,CAAC,SAAS,GAAG,YAAY,EAAE,SAAS,GAAG,YAAY,EAAE,CAAC,SAAS,CAAC;YAChE,CAAC,CAAC,YAAY,EAAE,YAAY,EAAE,CAAC,CAAC;YAChC,CAAC,SAAS,GAAG,YAAY,EAAE,SAAS,GAAG,YAAY,EAAE,SAAS,CAAC;SAChE,CAAC,CAAC;QAEH,MAAM,OAAO,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC;QAC5C,MAAM,YAAY,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,GAAG,YAAY,CAAC;QAC9D,MAAM,SAAS,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,GAAG,OAAO,CAAC;QAC5C,MAAM,OAAO,GAAG,IAAI,QAAQ,CAAC,CAAC,OAAO,GAAG,SAAS,EAAE,OAAO,GAAG,YAAY,GAAG,SAAS,EAAE,OAAO,CAAC,CAAC;QAChG,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,EAAE,CAAC,gBAAgB,CAAC,OAAO,CAAC,CAAC;QAChE,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,EAAE,CAAC,gBAAgB,CAAC,OAAO,CAAkC,CAAC;QAEjG,MAAM,OAAO,GAAG,IAAI,QAAQ,CAC1B,SAAS,GAAG,OAAO,GAAG,YAAY,EAClC,SAAS,GAAG,OAAO,GAAG,YAAY,EAClC,SAAS,GAAG,OAAO,CACpB,CAAC;QAEF,MAAM,OAAO,GAAG,SAAS,CAAC,SAAS,EAAE,CAAC,gBAAgB,CAAC,OAAO,CAAyB,CAAC;QAExF,OAAO,IAAI,KAAK,CAAC,MAAM,CAAC,KAAK,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;IACnD,CAAC;IAED,OAAO,CAAC,MAAa;QACnB,OAAO,IAAI,CAAC,aAAa,CAAC,MAAM,EAAE,aAAa,CAAC,YAAY,CAAC,MAAM,CAAC,QAAQ,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;IAClG,CAAC;IAED,MAAM,CAAC,gBAAgB,CAAC,CAAS,EAAE,UAAkB;QACnD,MAAM,CAAC,GAAG,UAAU,CAAC;QACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAE3B,OAAO,IAAI,MAAM,CAAC;YAChB,CAAC,GAAG,GAAG,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC;YAC/D,CAAC,GAAG,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC;YAC1F,CAAC,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,CAAC,CAAC;YAChC,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,CAAC;YAC3C,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,GAAG,GAAG,EAAE,GAAG,GAAG,EAAE,GAAG,CAAC;YACjE,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;SAClC,CAAC,CAAC;IACL,CAAC;IAED,UAAU,CAAC,CAAU;QACnB,MAAM,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;QAC1D,MAAM,GAAG,GAAG,MAAM,CAAC,cAAc,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC;QAE9E,OAAO,IAAI,IAAI,CACb,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAClB,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAe,EAAE,GAAG,CAAC,CAAC,CAAe,EAAE,GAAG,CAAC,CAAC,CAAe,CAAC,EAC9E,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAwB,EAAE,GAAG,CAAC,CAAC,CAAwB,EAAE,GAAG,CAAC,CAAC,CAAwB,CAAC,EACzG,IAAI,CAAC,cAAc,CACpB,CAAC;IACJ,CAAC;IAED,oBAAoB,CAAC,GAAW,EAAE,CAAU;QAC1C,MAAM,GAAG,GAAG,GAAG,CAAC,cAAc,CAAC,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,QAAQ,CAAC;QAE3E,OAAO,IAAI,IAAI,CACb,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,EAClB,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAe,EAAE,GAAG,CAAC,CAAC,CAAe,EAAE,GAAG,CAAC,CAAC,CAAe,CAAC,EAC9E,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAwB,EAAE,GAAG,CAAC,CAAC,CAAwB,EAAE,GAAG,CAAC,CAAC,CAAwB,CAAC,EACzG,IAAI,CAAC,cAAc,CACpB,CAAC;IACJ,CAAC;IAED,SAAS,CAAC,QAAkB;QAC1B,OAAO,IAAI,CAAC,UAAU,CAAC,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;IAC1D,CAAC;IAED,mBAAmB,CAAC,GAAW,EAAE,QAAkB;QACjD,OAAO,IAAI,CAAC,oBAAoB,CAAC,GAAG,EAAE,QAAQ,CAAC,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC;IACzE,CAAC;IAED,QAAQ,CAAC,QAAgB;QACvB,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QAE9C,OAAO,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,KAAK,CAAC,cAAc,CAAC,CAAC;IAC1G,CAAC;IAED,SAAS,CAAC,KAAe,EAAE,IAAc,EAAE,OAAO,IAAe;QAC/D,MAAM,MAAM,GAAW,EAAE,CAAC;QAC1B,IAAI,OAAO,GAAG,KAAK,CAAC;QAEpB,OAAO,IAAI,IAAI,OAAO,EAAE;YACtB,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC;YACrC,OAAO,GAAG,OAAO,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;SAC9B;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,IAAI,MAAM;QACR,OAAO,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,WAAsB,CAAC;IACrD,CAAC;IAED,iBAAiB;QACf,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC1B,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC7B,MAAM,IAAI,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC7B,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,WAAW,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,WAAsB,CAAC;QAEnG,IAAI,CAAC,IAAI,CAAC,EAAE;YACV,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,MAAiB,CAAC;SACtC;aAAM,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE;YACnB,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,MAAiB,CAAC;SAClC;QAED,OAAO,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5C,CAAC;IAED,aAAa,CAAC,QAAkB,EAAE,gBAAgB,GAAG,KAAK;QACxD,MAAM,CAAC,GAAG,QAAQ,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QAChD,MAAM,CAAC,GAAG,QAAQ,CAAC,gBAAgB,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,gBAAgB,CAAC,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;QAC1D,MAAM,cAAc,GAAG,IAAI,MAAM,CAAC;YAChC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SACtE,CAAC,CAAC;QACH,MAAM,UAAU,GAAG,CAAC;aACjB,QAAQ,CAAC,cAAc,CAAC,gBAAgB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAA6C,CAAC;QACxG,MAAM,cAAc,GAAG,IAAI,MAAM,CAAC;YAChC,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;SACtE,CAAC,CAAC;QACH,IAAI,MAAM,GAAG,cAAc;aACxB,OAAO,EAAE;aACT,gBAAgB,CAAC,UAAU,CAAC;aAC5B,QAAQ,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAE3B,IAAI,gBAAgB,EAAE;YACpB,MAAM,GAAG,IAAI,QAAQ,CACnB,MAAM,CAAC,CAAC,EACR,MAAM,CAAC,CAAC,EACR,CAAwB,CACzB,CAAC;SACH;QAED,OAAO,IAAI,MAAM,CACf,IAAI,CAAC,KAAK,EACV,MAAM,CAAC,CAAC,GAAG,IAAuB,EAClC,MAAM,CAAC,CAAC,GAAG,IAAuB,EAClC,MAAM,CAAC,CAAC,GAAG,IAAuB,CACnC,CAAC;IACJ,CAAC;IAED,gBAAgB,CACd,KAAe,EACf,SAAqB,EACrB,eAAyB,EAAE,EAC3B,gBAA0B,EAAE;QAE5B,IAAI,KAAK,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,KAAK,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,cAAc,CAAC,CAAC;QAEpF,YAAY,GAAG,YAAY,CAAC,KAAK,EAAE,CAAC;QACpC,aAAa,GAAG,aAAa,CAAC,KAAK,EAAE,CAAC;QACtC,IAAI,MAAM,GAAG,YAAY,CAAC;QAE1B,qCAAqC;QACrC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,CAAC,UAAU,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;QACrD,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,MAAM,IAAI,IAAI,CAAC,KAAK,IAAI,GAAG,CAAC,MAAM,IAAI,KAAK,CAAC,CAAC;QACjF,KAAK,MAAM,GAAG,IAAI,MAAM,EAAE;YACxB,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;SAC7B;QACD,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;QAC/B,KAAK,MAAM,GAAG,IAAI,SAAS,EAAE;YAC3B,MAAM,GAAG,GAAG,KAAK,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;YAErC,KAAK,GAAG,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;YAC5B,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;SAClB;QACD,MAAM,CAAC,IAAI,CAAC,GAAG,aAAa,CAAC,CAAC;QAE9B,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,cAAc,CAAC,QAAgB;QAC7B,MAAM,KAAK,GAAG,IAAI,CAAC,SAAS,CAAC,QAAQ,CAAC,MAAM,CAAC,CAAC;QAC9C,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC;QACxD,MAAM,MAAM,GAAG,KAAK,GAAG,QAAQ,CAAC,OAAO,GAAG,IAAI,CAAC;QAC/C,MAAM,MAAM,GAAG,KAAK,CAAC,EAAE,GAAG,CAAC,MAAM,GAAG,MAAM,CAAe,CAAC;QAE1D,OAAO,IAAI,IAAI,CAAC,KAAK,CAAC,KAAK,EAAE,KAAK,CAAC,QAAQ,EAAE,KAAK,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,MAAM,CAAC,EAAE,MAAM,CAAC,CAAC;IACjG,CAAC;IAED,IAAI,IAAI;QACN,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,QAAQ;QACN,OAAO;YACL,IAAI,IAAI,CAAC,IAAI,GAAG;YAChB,YAAY,IAAI,CAAC,KAAK,EAAE;YACxB,eAAe,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK;YAC7C,eAAe,IAAI,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC,CAAC,OAAO;SAChD,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;IACf,CAAC;CACF"}

package/dist/coordinate/index.d.ts

@@ -0,0 +1 @@
+export { Hill } from './Hill';

package/dist/coordinate/index.js

@@ -0,0 +1,2 @@
+export { Hill } from './Hill';
+//# sourceMappingURL=index.js.map

package/dist/coordinate/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/coordinate/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,IAAI,EAAE,MAAM,QAAQ,CAAC"}

package/dist/covariance/CovarianceSample.d.ts

@@ -0,0 +1,43 @@
+/**
+ * @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 matrix_;
+    /**
+     * 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 The origin state.
+     * @param covariance The covariance.
+     * @param originForceModel The force model for the origin state.
+     * @param sampleForceModel The force model for the samples.
+     */
+    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/dist/covariance/CovarianceSample.js

@@ -0,0 +1,147 @@
+/**
+ * @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 {
+    _origin;
+    _samples = [];
+    matrix_ = Matrix.zero(6, 12);
+    /**
+     * 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 The origin state.
+     * @param covariance The covariance.
+     * @param originForceModel The force model for the origin state.
+     * @param sampleForceModel The force model for the samples.
+     */
+    constructor(state, covariance, originForceModel, sampleForceModel) {
+        originForceModel ??= new ForceModel().setGravity();
+        sampleForceModel ??= new ForceModel().setGravity();
+        this._origin = new RungeKutta89Propagator(state, originForceModel);
+        const s = covariance.matrix.cholesky().elements;
+        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).elements;
+        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(matrix) {
+        const pts = matrix.elements;
+        const c = 1.0 / 12.0;
+        const yu = new Vector([0, 0, 0, 0, 0, 0]).elements;
+        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.elements[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.matrix_.elements[0][i] = state.position.x;
+            this.matrix_.elements[1][i] = state.position.y;
+            this.matrix_.elements[2][i] = state.position.z;
+            this.matrix_.elements[3][i] = state.velocity.x;
+            this.matrix_.elements[4][i] = state.velocity.y;
+            this.matrix_.elements[5][i] = state.velocity.z;
+        }
+        const matrix = this._rebuildCovariance(this.matrix_);
+        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.matrix_.elements[0][i] = state.position.x;
+            this.matrix_.elements[1][i] = state.position.y;
+            this.matrix_.elements[2][i] = state.position.z;
+            this.matrix_.elements[3][i] = state.velocity.x;
+            this.matrix_.elements[4][i] = state.velocity.y;
+            this.matrix_.elements[5][i] = state.velocity.z;
+        }
+        const matrix = this._rebuildCovariance(this.matrix_);
+        return new StateCovariance(matrix, CovarianceFrame.RIC);
+    }
+}
+//# sourceMappingURL=CovarianceSample.js.map

package/dist/covariance/CovarianceSample.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CovarianceSample.js","sourceRoot":"","sources":["../../src/covariance/CovarianceSample.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,2CAA2C;AAC3C,OAAO,EAGL,KAAK,EAGL,MAAM,EACN,aAAa,EACb,GAAG,EACH,MAAM,EACN,QAAQ,GACT,MAAM,WAAW,CAAC;AACnB,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;IACnB,OAAO,CAAyB;IAChC,QAAQ,GAA6B,EAAE,CAAC;IACxC,OAAO,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IAErC;;;;;;;;;OASG;IACH,YAAY,KAAY,EAAE,UAA2B,EAAE,gBAA6B,EAAE,gBAA6B;QACjH,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,QAAQ,CAAC;QAChD,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,QAAQ,CAAC;QAE7C,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,CAC1B,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAe,EAC5B,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAe,EAC5B,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAe,CAC7B,CAAC;YACF,MAAM,OAAO,GAAG,IAAI,QAAQ,CAC1B,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAwB,EACrC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAwB,EACrC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAwB,CACtC,CAAC;YAEF,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,MAAc;QACvC,MAAM,GAAG,GAAG,MAAM,CAAC,QAAQ,CAAC;QAC5B,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,QAAQ,CAAC;QACnD,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,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;aACtC;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,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;SAChD;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAErD,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,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;YAC/C,IAAI,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC,CAAC;SAChD;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,kBAAkB,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;QAErD,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;IAC1D,CAAC;CACF"}

package/dist/covariance/StateCovariance.d.ts

@@ -0,0 +1,42 @@
+/**
+ * @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].
+     * @param matrix The covariance matrix.
+     * @param frame The covariance frame.
+     * @returns A new [StateCovariance] object.
+     */
+    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/dist/covariance/StateCovariance.js

@@ -0,0 +1,65 @@
+/**
+ * @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 {
+    matrix;
+    frame;
+    /**
+     * Create a new [StateCovariance] object given its covariance [matrix] and
+     * [CovarianceFrame].
+     * @param matrix The covariance matrix.
+     * @param frame The covariance frame.
+     * @returns A new [StateCovariance] object.
+     */
+    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.elements[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.elements[i][i];
+            result[i] = Math.sqrt(variance);
+        }
+        return new Vector(result);
+    }
+}
+//# sourceMappingURL=StateCovariance.js.map

package/dist/covariance/StateCovariance.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"StateCovariance.js","sourceRoot":"","sources":["../../src/covariance/StateCovariance.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;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;IAQP;IAAuB;IAP1C;;;;;;OAMG;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,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC;SAChE;QAED,OAAO,IAAI,eAAe,CAAC,MAAM,EAAE,KAAK,CAAC,CAAC;IAC5C,CAAC;IAED;;;OAGG;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,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAE5C,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/dist/force/AtmosphericDrag.d.ts

@@ -0,0 +1,31 @@
+/**
+ * @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;
+}