ootk 3.0.0 → 4.0.0-1

package/{lib/objects/space-object.d.ts → dist/orbit_determination/BatchLeastSquaresResult.js}

@@ -1,13 +1,7 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
- * with satellites and other orbital objects.
- *
- * @file The SpaceObject class is meant to provide properties to satellites/stars that
- * do not apply to ground based 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
@@ -20,19 +14,23 @@
  * 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 { BaseObject } from './base-object';
-import { SpaceObjectType } from '../types/types';
-interface ObjectInfo {
-    name?: string;
-    type?: SpaceObjectType;
-    rcs?: number;
-    vmag?: number;
-}
-export declare class SpaceObject extends BaseObject {
-    rcs: number;
-    vmag: number;
-    constructor(info: ObjectInfo);
-    getRcs(): number;
-    getVmag(): number;
+// / Batch least squares orbit determination result.
+export class BatchLeastSquaresResult {
+    state;
+    covariance;
+    rms;
+    /**
+     * Create a new [BatchLeastSquaresResult] object, containing the solved
+     * [state], [covariance], and root-mean-squared error [rms].
+     * @param state The solved state.
+     * @param covariance The solved covariance.
+     * @param rms The root-mean-squared error.
+     */
+    constructor(state, covariance, rms) {
+        this.state = state;
+        this.covariance = covariance;
+        this.rms = rms;
+        // Nothing to do here.
+    }
 }
-export {};
+//# sourceMappingURL=BatchLeastSquaresResult.js.map

package/dist/orbit_determination/BatchLeastSquaresResult.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BatchLeastSquaresResult.js","sourceRoot":"","sources":["../../src/orbit_determination/BatchLeastSquaresResult.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAKH,oDAAoD;AACpD,MAAM,OAAO,uBAAuB;IAQf;IAAqB;IAAoC;IAP5E;;;;;;OAMG;IACH,YAAmB,KAAY,EAAS,UAA2B,EAAS,GAAW;QAApE,UAAK,GAAL,KAAK,CAAO;QAAS,eAAU,GAAV,UAAU,CAAiB;QAAS,QAAG,GAAH,GAAG,CAAQ;QACrF,sBAAsB;IACxB,CAAC;CACF"}

package/dist/orbit_determination/GibbsIOD.d.ts

@@ -0,0 +1,38 @@
+/**
+ * @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, Kilometers, Vector3D } from 'ootk-core';
+/**
+ * Gibbs 3-position inital orbit determination.
+ */
+export declare class GibbsIOD {
+    mu: number;
+    constructor(mu?: number);
+    /** Abort solve if position plane exceeds this value. */
+    private static readonly coplanarThreshold_;
+    /**
+     * Attempt to create a state estimate from three inertial position vectors.
+     *
+     * Throws an error if the positions are not coplanar.
+     * @param r1 Position vector 1.
+     * @param r2 Position vector 2.
+     * @param r3 Position vector 3.
+     * @param t2 Time of position 2.
+     * @param t3 Time of position 3.
+     * @returns State estimate at time t2.
+     */
+    solve(r1: Vector3D<Kilometers>, r2: Vector3D<Kilometers>, r3: Vector3D<Kilometers>, t2: EpochUTC, t3: EpochUTC): J2000;
+}

package/dist/orbit_determination/GibbsIOD.js

@@ -0,0 +1,76 @@
+/**
+ * @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 { DEG2RAD, Earth, halfPi, J2000 } from 'ootk-core';
+import { ForceModel } from './../force/ForceModel';
+import { RungeKutta89Propagator } from './../propagator/RungeKutta89Propagator';
+/**
+ * Gibbs 3-position inital orbit determination.
+ */
+export class GibbsIOD {
+    mu;
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+        // Nothing to do here.
+    }
+    /** Abort solve if position plane exceeds this value. */
+    static coplanarThreshold_ = (5.0 * DEG2RAD);
+    /**
+     * Attempt to create a state estimate from three inertial position vectors.
+     *
+     * Throws an error if the positions are not coplanar.
+     * @param r1 Position vector 1.
+     * @param r2 Position vector 2.
+     * @param r3 Position vector 3.
+     * @param t2 Time of position 2.
+     * @param t3 Time of position 3.
+     * @returns State estimate at time t2.
+     */
+    solve(r1, r2, r3, t2, t3) {
+        const num = r1.normalize().dot(r2.normalize().cross(r3.normalize()));
+        const alpha = halfPi - Math.acos(num);
+        if (Math.abs(alpha) > GibbsIOD.coplanarThreshold_) {
+            throw new Error('Orbits are not coplanar.');
+        }
+        const r1m = r1.magnitude();
+        const r2m = r2.magnitude();
+        const r3m = r3.magnitude();
+        const d = r1.cross(r2).add(r2.cross(r3).add(r3.cross(r1)));
+        const n = r2.cross(r3).scale(r1m).add(r3.cross(r1).scale(r2m)).add(r1.cross(r2).scale(r3m));
+        const b = d.cross(r2);
+        const s = r1.scale(r2m - r3m).add(r2.scale(r3m - r1m).add(r3.scale(r1m - r2m)));
+        const nm = n.magnitude();
+        const dm = d.magnitude();
+        const vm = Math.sqrt(this.mu / (nm * dm));
+        const vlEci = b.scale((vm / r2m)).add(s.scale(vm));
+        const pv = new J2000(t2, r2, vlEci);
+        const forceModel = new ForceModel().setGravity(this.mu);
+        const orbit = new RungeKutta89Propagator(pv, forceModel);
+        const pv2 = new J2000(t2, r2, vlEci.negate());
+        const orbit2 = new RungeKutta89Propagator(pv2, forceModel);
+        const estP3 = orbit.propagate(t3).position;
+        const dist = estP3.subtract(r3).magnitude();
+        const estP3b = orbit2.propagate(t3).position;
+        const dist2 = estP3b.subtract(r3).magnitude();
+        if (dist <= dist2) {
+            orbit.reset();
+            return orbit.propagate(t2);
+        }
+        orbit2.reset();
+        return orbit2.propagate(t2);
+    }
+}
+//# sourceMappingURL=GibbsIOD.js.map

package/dist/orbit_determination/GibbsIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"GibbsIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/GibbsIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,OAAO,EAAE,KAAK,EAAY,MAAM,EAAE,KAAK,EAAsD,MAAM,WAAW,CAAC;AACxH,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,sBAAsB,EAAE,MAAM,wCAAwC,CAAC;AAEhF;;GAEG;AACH,MAAM,OAAO,QAAQ;IACA;IAAnB,YAAmB,KAAa,KAAK,CAAC,EAAE;QAArB,OAAE,GAAF,EAAE,CAAmB;QACtC,sBAAsB;IACxB,CAAC;IAED,wDAAwD;IAChD,MAAM,CAAU,kBAAkB,GAAY,CAAC,GAAG,GAAG,OAAO,CAAY,CAAC;IAEjF;;;;;;;;;;OAUG;IACH,KAAK,CACH,EAAwB,EACxB,EAAwB,EACxB,EAAwB,EACxB,EAAY,EACZ,EAAY;QAEZ,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,SAAS,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACrE,MAAM,KAAK,GAAG,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAEtC,IAAI,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,QAAQ,CAAC,kBAAkB,EAAE;YACjD,MAAM,IAAI,KAAK,CAAC,0BAA0B,CAAC,CAAC;SAC7C;QAED,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAE3B,MAAM,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC3D,MAAM,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5F,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QACtB,MAAM,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAEhF,MAAM,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;QAEzB,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAwB,CAAC;QACjE,MAAM,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,GAAG,CAAwB,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;QAE1E,MAAM,EAAE,GAAG,IAAI,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,CAAC;QAEpC,MAAM,UAAU,GAAG,IAAI,UAAU,EAAE,CAAC,UAAU,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACxD,MAAM,KAAK,GAAG,IAAI,sBAAsB,CAAC,EAAE,EAAE,UAAU,CAAC,CAAC;QAEzD,MAAM,GAAG,GAAG,IAAI,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC;QAC9C,MAAM,MAAM,GAAG,IAAI,sBAAsB,CAAC,GAAG,EAAE,UAAU,CAAC,CAAC;QAE3D,MAAM,KAAK,GAAG,KAAK,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC;QAC3C,MAAM,IAAI,GAAG,KAAK,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,SAAS,EAAE,CAAC;QAC5C,MAAM,MAAM,GAAG,MAAM,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC;QAC7C,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,SAAS,EAAE,CAAC;QAE9C,IAAI,IAAI,IAAI,KAAK,EAAE;YACjB,KAAK,CAAC,KAAK,EAAE,CAAC;YAEd,OAAO,KAAK,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;SAC5B;QACD,MAAM,CAAC,KAAK,EAAE,CAAC;QAEf,OAAO,MAAM,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC;IAC9B,CAAC"}

package/dist/orbit_determination/GoodingIOD.d.ts

@@ -0,0 +1,51 @@
+/**
+ * @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, Kilometers } from 'ootk-core';
+import { ObservationOptical } from './../observation/ObservationOptical';
+/**
+ * Gooding angles-only initial orbit determination.
+ *
+ * Used for orbit determination from three optical observations.
+ */
+export declare class GoodingIOD {
+    private readonly _mu;
+    private readonly o1_;
+    private readonly o2_;
+    private readonly o3_;
+    private vObserverPosition1_;
+    private vObserverPosition2_;
+    private vObserverPosition3_;
+    private r_;
+    private v_;
+    private t_;
+    private r1_;
+    private r2_;
+    private r3_;
+    private rho1_;
+    private rho2_;
+    private rho3_;
+    private d1_;
+    private d3_;
+    private facFiniteDiff_;
+    private readonly _forceModel;
+    constructor(o1: ObservationOptical, o2: ObservationOptical, o3: ObservationOptical, mu?: number);
+    private _getPositionOnLoS2;
+    private _modifyIterate;
+    private _computeDerivatives;
+    solve(r1Init: Kilometers, r3Init: Kilometers, nRev?: number, direction?: boolean): J2000;
+    private _solveRangeProblem;
+}

package/dist/orbit_determination/GoodingIOD.js

@@ -0,0 +1,247 @@
+/**
+ * @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, TAU, Vector3D } from 'ootk-core';
+import { ForceModel } from '../force/ForceModel';
+import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator';
+import { GibbsIOD } from './GibbsIOD';
+import { LambertIOD } from './LambertIOD';
+/**
+ * Gooding angles-only initial orbit determination.
+ *
+ * Used for orbit determination from three optical observations.
+ */
+export class GoodingIOD {
+    _mu;
+    o1_;
+    o2_;
+    o3_;
+    vObserverPosition1_ = Vector3D.origin;
+    vObserverPosition2_ = Vector3D.origin;
+    vObserverPosition3_ = Vector3D.origin;
+    r_ = 0.0;
+    v_ = 0.0;
+    t_ = 0.0;
+    r1_ = 0.0;
+    r2_ = 0.0;
+    r3_ = 0.0;
+    rho1_ = 0.0;
+    rho2_ = 0.0;
+    rho3_ = 0.0;
+    d1_ = 0.0;
+    d3_ = 0.0;
+    facFiniteDiff_ = 0.0;
+    _forceModel = new ForceModel().setGravity(1.0);
+    constructor(o1, o2, o3, mu = Earth.mu) {
+        this._mu = mu;
+        this.o1_ = o1;
+        this.o2_ = o2;
+        this.o3_ = o3;
+    }
+    _getPositionOnLoS2({ e1, r01, e3, r03, t13, t12, nRev, posigrade, }) {
+        const p1 = this.vObserverPosition1_.add(e1.scale(r01));
+        this.r1_ = p1.magnitude();
+        const p3 = this.vObserverPosition3_.add(e3.scale(r03));
+        this.r3_ = p3.magnitude();
+        const p13 = p1.cross(p3);
+        let th = Math.atan2(p13.magnitude(), p1.dot(p3));
+        if (!posigrade) {
+            th = TAU - th;
+        }
+        const v1 = new Float64Array(2);
+        const exitflag = LambertIOD.solve(this.r1_, this.r3_, th, t13, nRev, v1);
+        if (exitflag) {
+            const pn = p1.cross(p3);
+            const pt = pn.cross(p1);
+            let rt = pt.magnitude();
+            if (!posigrade) {
+                rt = -rt;
+            }
+            const vel1 = p1.scale(v1[0] / this.r1_).add(pt.scale(v1[1] / rt));
+            const p2 = new RungeKutta89Propagator(new J2000(this.o1_.epoch, p1, vel1), this._forceModel).propagate(this.o1_.epoch.roll(t12)).position;
+            return p2;
+        }
+        return null;
+    }
+    _modifyIterate(lineOfSight1, lineOfSight3) {
+        const r13 = this.vObserverPosition3_.subtract(this.vObserverPosition1_);
+        this.d1_ = r13.dot(lineOfSight1);
+        this.d3_ = r13.dot(lineOfSight3);
+        const d2 = lineOfSight1.dot(lineOfSight3);
+        const d4 = 1.0 - d2 * d2;
+        this.rho1_ = Math.max((this.d1_ - this.d3_ * d2) / d4, 0.0);
+        this.rho3_ = Math.max((this.d1_ * d2 - this.d3_) / d4, 0.0);
+    }
+    _computeDerivatives({ x, y, lineOfSight1, lineOfSight3, pin, ein, t13, t12, nrev, direction, fd, gd, }) {
+        const p = pin.normalize();
+        const en = ein.normalize();
+        const dx = this.facFiniteDiff_ * x;
+        const dy = this.facFiniteDiff_ * y;
+        const cm1 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x - dx,
+            e3: lineOfSight3,
+            r03: y,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this.vObserverPosition2_);
+        const fm1 = p.dot(cm1);
+        const gm1 = en.dot(cm1);
+        const cp1 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x + dx,
+            e3: lineOfSight3,
+            r03: y,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this.vObserverPosition2_);
+        const fp1 = p.dot(cp1);
+        const gp1 = en.dot(cp1);
+        const fx = (fp1 - fm1) / (2.0 * dx);
+        const gx = (gp1 - gm1) / (2.0 * dx);
+        const cm3 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x,
+            e3: lineOfSight3,
+            r03: y - dy,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this.vObserverPosition2_);
+        const fm3 = p.dot(cm3);
+        const gm3 = en.dot(cm3);
+        const cp3 = this._getPositionOnLoS2({
+            e1: lineOfSight1,
+            r01: x,
+            e3: lineOfSight3,
+            r03: y + dy,
+            t13,
+            t12,
+            nRev: nrev,
+            posigrade: direction,
+        }).subtract(this.vObserverPosition2_);
+        const fp3 = p.dot(cp3);
+        const gp3 = en.dot(cp3);
+        const fy = (fp3 - fm3) / (2.0 * dy);
+        const gy = (gp3 - gm3) / (2.0 * dy);
+        fd[0] = fx;
+        fd[1] = fy;
+        gd[0] = gx;
+        gd[1] = gy;
+    }
+    solve(r1Init, r3Init, nRev = 0, direction = true) {
+        const lineOfSight1 = this.o1_.observation.lineOfSight();
+        const lineOfSight2 = this.o2_.observation.lineOfSight();
+        const lineOfSight3 = this.o3_.observation.lineOfSight();
+        this.r_ = Math.max(r1Init, r3Init);
+        this.v_ = Math.sqrt(this._mu / this.r_);
+        this.t_ = this.r_ / this.v_;
+        this.vObserverPosition1_ = this.o1_.site.position.scale(1.0 / this.r_);
+        this.vObserverPosition2_ = this.o2_.site.position.scale(1.0 / this.r_);
+        this.vObserverPosition3_ = this.o3_.site.position.scale(1.0 / this.r_);
+        const maxiter = 100;
+        this._solveRangeProblem({
+            rho1init: r1Init / this.r_,
+            rho3init: r3Init / this.r_,
+            t13: this.o3_.epoch.difference(this.o1_.epoch) / this.t_,
+            t12: this.o2_.epoch.difference(this.o1_.epoch) / this.t_,
+            nrev: nRev,
+            direction,
+            lineOfSight1,
+            lineOfSight2,
+            lineOfSight3,
+            maxIterations: maxiter,
+        });
+        const gibbs = new GibbsIOD(this._mu);
+        const p1 = this.vObserverPosition1_.add(lineOfSight1.scale(this.rho1_)).scale(this.r_);
+        const p2 = this.vObserverPosition2_.add(lineOfSight2.scale(this.rho2_)).scale(this.r_);
+        const p3 = this.vObserverPosition3_.add(lineOfSight3.scale(this.rho3_)).scale(this.r_);
+        return gibbs.solve(p1, p2, p3, this.o2_.epoch, this.o3_.epoch);
+    }
+    _solveRangeProblem({ rho1init, rho3init, t13, t12, nrev, direction, lineOfSight1, lineOfSight2, lineOfSight3, maxIterations, }) {
+        const arbf = 1e-6;
+        const cvtol = 1e-14;
+        this.rho1_ = rho1init;
+        this.rho3_ = rho3init;
+        let iter = 0;
+        let stoppingCriterion = 10.0 * cvtol;
+        while (iter < maxIterations && Math.abs(stoppingCriterion) > cvtol) {
+            this.facFiniteDiff_ = arbf;
+            const p2 = this._getPositionOnLoS2({
+                e1: lineOfSight1,
+                r01: this.rho1_,
+                e3: lineOfSight3,
+                r03: this.rho3_,
+                t13,
+                t12,
+                nRev: nrev,
+                posigrade: direction,
+            });
+            if (p2 === null) {
+                this._modifyIterate(lineOfSight1, lineOfSight3);
+            }
+            else {
+                this.r2_ = p2.magnitude();
+                const c = p2.subtract(this.vObserverPosition2_);
+                this.rho2_ = c.magnitude();
+                const cr = lineOfSight2.dot(c);
+                const u = lineOfSight2.cross(c);
+                const p = u.cross(lineOfSight2).normalize();
+                const ent = lineOfSight2.cross(p);
+                const enr = ent.magnitude();
+                if (enr === 0.0) {
+                    return;
+                }
+                const en = ent.normalize();
+                const fc = p.dot(c);
+                const fd = new Float64Array(2);
+                const gd = new Float64Array(2);
+                this._computeDerivatives({
+                    x: this.rho1_,
+                    y: this.rho3_,
+                    lineOfSight1,
+                    lineOfSight3,
+                    pin: p,
+                    ein: en,
+                    t13,
+                    t12,
+                    nrev,
+                    direction,
+                    fd,
+                    gd,
+                });
+                const fr1 = fd[0];
+                const fr3 = fd[1];
+                const gr1 = gd[0];
+                const gr3 = gd[1];
+                const detj = fr1 * gr3 - fr3 * gr1;
+                this.d3_ = (-gr3 * fc) / detj;
+                this.d1_ = (gr1 * fc) / detj;
+                this.rho1_ = this.rho1_ + this.d3_;
+                this.rho3_ = this.rho3_ + this.d1_;
+                const den = Math.max(cr, this.r2_);
+                stoppingCriterion = fc / den;
+            }
+            ++iter;
+        }
+    }
+}
+//# sourceMappingURL=GoodingIOD.js.map

package/dist/orbit_determination/GoodingIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"GoodingIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/GoodingIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,KAAK,EAAE,KAAK,EAA4C,GAAG,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAClG,OAAO,EAAE,UAAU,EAAE,MAAM,qBAAqB,CAAC;AACjD,OAAO,EAAE,sBAAsB,EAAE,MAAM,sCAAsC,CAAC;AAE9E,OAAO,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AACtC,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAyC1C;;;;GAIG;AACH,MAAM,OAAO,UAAU;IACJ,GAAG,CAAS;IACZ,GAAG,CAAqB;IACxB,GAAG,CAAqB;IACxB,GAAG,CAAqB;IACjC,mBAAmB,GAAG,QAAQ,CAAC,MAA8B,CAAC;IAC9D,mBAAmB,GAAG,QAAQ,CAAC,MAA8B,CAAC;IAC9D,mBAAmB,GAAG,QAAQ,CAAC,MAA8B,CAAC;IAC9D,EAAE,GAAG,GAAiB,CAAC;IACvB,EAAE,GAAG,GAAG,CAAC;IACT,EAAE,GAAG,GAAG,CAAC;IACT,GAAG,GAAG,GAAG,CAAC;IACV,GAAG,GAAG,GAAG,CAAC;IACV,GAAG,GAAG,GAAG,CAAC;IACV,KAAK,GAAG,GAAiB,CAAC;IAC1B,KAAK,GAAG,GAAiB,CAAC;IAC1B,KAAK,GAAG,GAAiB,CAAC;IAC1B,GAAG,GAAG,GAAiB,CAAC;IACxB,GAAG,GAAG,GAAiB,CAAC;IACxB,cAAc,GAAG,GAAG,CAAC;IACZ,WAAW,GAAe,IAAI,UAAU,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC;IAE5E,YAAY,EAAsB,EAAE,EAAsB,EAAE,EAAsB,EAAE,KAAa,KAAK,CAAC,EAAE;QACvG,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;QACd,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC;IAChB,CAAC;IAEO,kBAAkB,CAAC,EACzB,EAAE,EACF,GAAG,EACH,EAAE,EACF,GAAG,EACH,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,GACe;QACxB,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvD,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAE1B,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAEvD,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAE1B,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;QAEzB,IAAI,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,SAAS,EAAE,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,IAAI,CAAC,SAAS,EAAE;YACd,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;SACf;QAED,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAC/B,MAAM,QAAQ,GAAG,UAAU,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC,CAAC;QAEzE,IAAI,QAAQ,EAAE;YACZ,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YACxB,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC;YAExB,IAAI,EAAE,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;YAExB,IAAI,CAAC,SAAS,EAAE;gBACd,EAAE,GAAG,CAAC,EAAgB,CAAC;aACxB;YAED,MAAM,IAAI,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAkC,CAAC;YAEnG,MAAM,EAAE,GAAG,IAAI,sBAAsB,CAAC,IAAI,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,EAAE,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC,SAAS,CACpG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CACzB,CAAC,QAAQ,CAAC;YAEX,OAAO,EAAE,CAAC;SACX;QAED,OAAO,IAAI,CAAC;IACd,CAAC;IAEO,cAAc,CAAC,YAAsB,EAAE,YAAsB;QACnE,MAAM,GAAG,GAAG,IAAI,CAAC,mBAAmB,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAExE,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,YAAY,CAAe,CAAC;QAC/C,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,YAAY,CAAe,CAAC;QAC/C,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;QAC1C,MAAM,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC;QAEzB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,EAAE,EAAE,GAAG,CAAe,CAAC;QAC1E,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,CAAe,CAAC;IAC5E,CAAC;IAEO,mBAAmB,CAAC,EAC1B,CAAC,EACD,CAAC,EACD,YAAY,EACZ,YAAY,EACZ,GAAG,EACH,GAAG,EACH,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,EACT,EAAE,EACF,EAAE,GACuB;QACzB,MAAM,CAAC,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAC1B,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;QAE3B,MAAM,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,CAAe,CAAC;QACjD,MAAM,EAAE,GAAG,IAAI,CAAC,cAAc,GAAG,CAAe,CAAC;QAEjD,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAgB;YACzB,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAgB;YACzB,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QAEpC,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAgB;YACzB,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,GAAG,GAAG,IAAI,CAAC,kBAAkB,CAAC;YAClC,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC;YACN,EAAE,EAAE,YAAY;YAChB,GAAG,EAAE,CAAC,GAAG,EAAgB;YACzB,GAAG;YACH,GAAG;YACH,IAAI,EAAE,IAAI;YACV,SAAS,EAAE,SAAS;SACrB,CAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;QAEvC,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAExB,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QACpC,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;QAEpC,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACX,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACb,CAAC;IAED,KAAK,CAAC,MAAkB,EAAE,MAAkB,EAAE,IAAI,GAAG,CAAC,EAAE,SAAS,GAAG,IAAI;QACtE,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAA0B,CAAC;QAChF,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAA0B,CAAC;QAChF,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,EAA0B,CAAC;QAEhF,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,MAAM,CAAe,CAAC;QACjD,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACxC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC;QAE5B,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAgB,CAAC,CAAC;QACrF,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAgB,CAAC,CAAC;QACrF,IAAI,CAAC,mBAAmB,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAC,EAAgB,CAAC,CAAC;QAErF,MAAM,OAAO,GAAG,GAAG,CAAC;QAEpB,IAAI,CAAC,kBAAkB,CAAC;YACtB,QAAQ,EAAE,MAAM,GAAG,IAAI,CAAC,EAAgB;YACxC,QAAQ,EAAE,MAAM,GAAG,IAAI,CAAC,EAAgB;YACxC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAa;YACnE,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,EAAa;YACnE,IAAI,EAAE,IAAI;YACV,SAAS;YACT,YAAY;YACZ,YAAY;YACZ,YAAY;YACZ,aAAa,EAAE,OAAO;SACvB,CAAC,CAAC;QACH,MAAM,KAAK,GAAG,IAAI,QAAQ,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrC,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,IAAI,CAAC,mBAAmB,CAAC,GAAG,CAAC,YAAY,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;QAEvF,OAAO,KAAK,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IACjE,CAAC;IAEO,kBAAkB,CAAC,EACzB,QAAQ,EACR,QAAQ,EACR,GAAG,EACH,GAAG,EACH,IAAI,EACJ,SAAS,EACT,YAAY,EACZ,YAAY,EACZ,YAAY,EACZ,aAAa,GACW;QACxB,MAAM,IAAI,GAAG,IAAI,CAAC;QAClB,MAAM,KAAK,GAAG,KAAK,CAAC;QAEpB,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC;QACtB,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC;QAEtB,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,iBAAiB,GAAG,IAAI,GAAG,KAAK,CAAC;QAErC,OAAO,IAAI,GAAG,aAAa,IAAI,IAAI,CAAC,GAAG,CAAC,iBAAiB,CAAC,GAAG,KAAK,EAAE;YAClE,IAAI,CAAC,cAAc,GAAG,IAAI,CAAC;YAE3B,MAAM,EAAE,GAAG,IAAI,CAAC,kBAAkB,CAAC;gBACjC,EAAE,EAAE,YAAY;gBAChB,GAAG,EAAE,IAAI,CAAC,KAAK;gBACf,EAAE,EAAE,YAAY;gBAChB,GAAG,EAAE,IAAI,CAAC,KAAK;gBACf,GAAG;gBACH,GAAG;gBACH,IAAI,EAAE,IAAI;gBACV,SAAS,EAAE,SAAS;aACrB,CAAC,CAAC;YAEH,IAAI,EAAE,KAAK,IAAI,EAAE;gBACf,IAAI,CAAC,cAAc,CAAC,YAAY,EAAE,YAAY,CAAC,CAAC;aACjD;iBAAM;gBACL,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;gBAC1B,MAAM,CAAC,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,mBAAmB,CAAC,CAAC;gBAEhD,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,SAAS,EAAE,CAAC;gBAC3B,MAAM,EAAE,GAAG,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAE/B,MAAM,CAAC,GAAG,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAChC,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,SAAS,EAAE,CAAC;gBAC5C,MAAM,GAAG,GAAG,YAAY,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBAElC,MAAM,GAAG,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;gBAE5B,IAAI,GAAG,KAAK,GAAG,EAAE;oBACf,OAAO;iBACR;gBAED,MAAM,EAAE,GAAG,GAAG,CAAC,SAAS,EAAE,CAAC;gBAE3B,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAEpB,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;gBAC/B,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;gBAE/B,IAAI,CAAC,mBAAmB,CAAC;oBACvB,CAAC,EAAE,IAAI,CAAC,KAAK;oBACb,CAAC,EAAE,IAAI,CAAC,KAAK;oBACb,YAAY;oBACZ,YAAY;oBACZ,GAAG,EAAE,CAAC;oBACN,GAAG,EAAE,EAAE;oBACP,GAAG;oBACH,GAAG;oBACH,IAAI;oBACJ,SAAS;oBACT,EAAE;oBACF,EAAE;iBACH,CAAC,CAAC;gBAEH,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBAClB,MAAM,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;gBAEnC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAkB,CAAC;gBAC5C,IAAI,CAAC,GAAG,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAkB,CAAC;gBAE3C,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAiB,CAAC;gBACjD,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAiB,CAAC;gBAEjD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;gBAEnC,iBAAiB,GAAG,EAAE,GAAG,GAAG,CAAC;aAC9B;YAED,EAAE,IAAI,CAAC;SACR;IACH,CAAC;CACF"}

package/dist/orbit_determination/HerrickGibbsIOD.d.ts

@@ -0,0 +1,33 @@
+/**
+ * @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, Kilometers, Vector3D } from 'ootk-core';
+/**
+ * Herrik-Gibbs 3-position initial orbit determination.
+ *
+ * Possibly better than regular Gibbs IOD for closely spaced position
+ * vectors (less than 5°).
+ */
+export declare class HerrickGibbsIOD {
+    mu: number;
+    /**
+     * Create a new [HerrickGibbsIOD] object with optional
+     * gravitational parameter [mu].
+     * @param mu Gravitational parameter (default: Earth.mu)
+     */
+    constructor(mu?: number);
+    solve(r1: Vector3D, t1: EpochUTC, r2: Vector3D<Kilometers>, t2: EpochUTC, r3: Vector3D, t3: EpochUTC): J2000;
+}

package/dist/orbit_determination/HerrickGibbsIOD.js

@@ -0,0 +1,50 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Earth, J2000 } from 'ootk-core';
+/**
+ * Herrik-Gibbs 3-position initial orbit determination.
+ *
+ * Possibly better than regular Gibbs IOD for closely spaced position
+ * vectors (less than 5°).
+ */
+export class HerrickGibbsIOD {
+    mu;
+    /**
+     * Create a new [HerrickGibbsIOD] object with optional
+     * gravitational parameter [mu].
+     * @param mu Gravitational parameter (default: Earth.mu)
+     */
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+        // Nothing to do here.
+    }
+    // / Attempt to create a state estimate from three inertial position vectors.
+    solve(r1, t1, r2, t2, r3, t3) {
+        const dt31 = t3.difference(t1);
+        const dt32 = t3.difference(t2);
+        const dt21 = t2.difference(t1);
+        const r1m = r1.magnitude();
+        const r2m = r2.magnitude();
+        const r3m = r3.magnitude();
+        const vA = r1.scale(-dt32 * (1.0 / (dt21 * dt31) + this.mu / (12.0 * r1m * r1m * r1m)));
+        const vB = r2.scale((dt32 - dt21) * (1.0 / (dt21 * dt32) + this.mu / (12.0 * r2m * r2m * r2m)));
+        const vC = r3.scale(dt21 * (1.0 / (dt32 * dt31) + this.mu / (12.0 * r3m * r3m * r3m)));
+        const v2 = vA.add(vB).add(vC);
+        return new J2000(t2, r2, v2);
+    }
+}
+//# sourceMappingURL=HerrickGibbsIOD.js.map

package/dist/orbit_determination/HerrickGibbsIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"HerrickGibbsIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/HerrickGibbsIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,KAAK,EAAY,KAAK,EAA6C,MAAM,WAAW,CAAC;AAE9F;;;;;GAKG;AACH,MAAM,OAAO,eAAe;IAMP;IALnB;;;;OAIG;IACH,YAAmB,KAAa,KAAK,CAAC,EAAE;QAArB,OAAE,GAAF,EAAE,CAAmB;QACtC,sBAAsB;IACxB,CAAC;IAED,6EAA6E;IAC7E,KAAK,CAAC,EAAY,EAAE,EAAY,EAAE,EAAwB,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY;QAClG,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,IAAI,GAAG,EAAE,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC;QAC/B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,GAAG,GAAG,EAAE,CAAC,SAAS,EAAE,CAAC;QAC3B,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACxF,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAChG,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAkC,CAAC;QAE/D,OAAO,IAAI,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IAC/B,CAAC;CACF"}

package/dist/orbit_determination/LambertIOD.d.ts

@@ -0,0 +1,59 @@
+/**
+ * @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, Kilometers, Vector3D } from 'ootk-core';
+export declare class LambertIOD {
+    mu: number;
+    constructor(mu?: number);
+    /**
+     * Try to guess the short path argument given an [interceptor] and
+     * [target] state.
+     * @param interceptor Interceptor
+     * @param target Target
+     * @returns True if the short path should be used, false otherwise.
+     */
+    static useShortPath(interceptor: J2000, target: J2000): boolean;
+    private static timeOfFlight_;
+    /**
+     * Attempt to solve output velocity [v1] _(km/s)_ given radii [r1] and
+     * [r2] _(canonical)_, sweep angle [dth] _(rad)_, time of flight [tau]
+     * _(canonical)_, and number of revolutions _(mRev)_.
+     * @param r1 Radius 1
+     * @param r2 Radius 2
+     * @param dth Sweep angle
+     * @param tau Time of flight
+     * @param mRev Number of revolutions
+     * @param v1 Output velocity
+     * @returns True if successful, false otherwise.
+     */
+    static solve(r1: number, r2: number, dth: number, tau: number, mRev: number, v1: Float64Array): boolean;
+    /**
+     * Estimate a state vector for inertial position [p1] _(km)_ given the
+     * two epoch and positions.
+     * @param p1 Position vector 1
+     * @param p2 Position vector 2
+     * @param t1 Epoch 1
+     * @param t2 Epoch 2
+     * @param root0 Optional parameters
+     * @param root0.posigrade If true, use the positive root (default: true)
+     * @param root0.nRev Number of revolutions (default: 0)
+     * @returns A [J2000] object with the estimated state vector.
+     */
+    estimate(p1: Vector3D<Kilometers>, p2: Vector3D<Kilometers>, t1: EpochUTC, t2: EpochUTC, { posigrade, nRev }?: {
+        posigrade?: boolean;
+        nRev?: number;
+    }): J2000 | null;
+}