ootk 2.0.1 → 4.0.0-0

package/lib/optimize/SimpleLinearRegression.js

@@ -0,0 +1,107 @@
+/**
+ * @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/>.
+ */
+// / Simple linear regression _(y = mx + b)_.
+export class SimpleLinearRegression {
+    /**
+     * Create a new [SimpleLinearRegression] object from lists of x and y
+     * values.
+     */
+    constructor(xs, ys) {
+        this.xs = xs;
+        this.ys = ys;
+        // / Line slope
+        this._slope = 0.0;
+        // / Y-axis intercept
+        this._intercept = 0.0;
+        this._error = 0.0;
+        this.update();
+    }
+    // / Line slope
+    get slope() {
+        return this._slope;
+    }
+    // / Y-axis intercept
+    get intercept() {
+        return this._intercept;
+    }
+    // / Linear regression standard deviation
+    get error() {
+        return this._error;
+    }
+    // / Data length
+    get length() {
+        return Math.min(this.xs.length, this.ys.length);
+    }
+    _calcError() {
+        let total = 0.0;
+        for (let i = 0; i < this.length; i++) {
+            const delta = this.ys[i] - this.evaluate(this.xs[i]);
+            total += delta * delta;
+        }
+        this._error = Math.sqrt(total / (this.length - 1));
+    }
+    // / Update the linear fit with this object's current [xs] and [ys] values.
+    update() {
+        const n = Math.min(this.xs.length, this.ys.length);
+        let xMu = 0.0;
+        let yMu = 0.0;
+        for (let i = 0; i < n; i++) {
+            xMu += this.xs[i];
+            yMu += this.ys[i];
+        }
+        xMu /= n;
+        yMu /= n;
+        let pa = 0.0;
+        let xSig = 0.0;
+        let ySig = 0.0;
+        for (let i = 0; i < n; i++) {
+            const xd = this.xs[i] - xMu;
+            const yd = this.ys[i] - yMu;
+            pa += xd * yd;
+            xSig += xd * xd;
+            ySig += yd * yd;
+        }
+        const p = pa / (Math.sqrt(xSig) * Math.sqrt(ySig));
+        xSig = Math.sqrt(xSig / (n - 1));
+        ySig = Math.sqrt(ySig / (n - 1));
+        this._slope = p * (ySig / xSig);
+        this._intercept = yMu - this._slope * xMu;
+        this._calcError();
+    }
+    // / Evaluate this linear fit for y, given an [x] value.
+    evaluate(x) {
+        return this._slope * x + this._intercept;
+    }
+    /**
+     * Create a new [SimpleLinearRegression] object with outliers above the
+     * provided standard deviation [sigma] value removed.
+     */
+    filterOutliers(sigma = 1.0) {
+        const limit = this.error * sigma;
+        const xsOut = [];
+        const ysOut = [];
+        for (let i = 0; i < this.length; i++) {
+            if (Math.abs(this.ys[i] - this.evaluate(this.xs[i])) < limit) {
+                xsOut.push(this.xs[i]);
+                ysOut.push(this.ys[i]);
+            }
+        }
+        return new SimpleLinearRegression(xsOut, ysOut);
+    }
+}
+//# sourceMappingURL=SimpleLinearRegression.js.map

package/lib/optimize/SimpleLinearRegression.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SimpleLinearRegression.js","sourceRoot":"","sources":["../../src/optimize/SimpleLinearRegression.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,6CAA6C;AAC7C,MAAM,OAAO,sBAAsB;IACjC;;;OAGG;IACH,YAAmB,EAAY,EAAS,EAAY;QAAjC,OAAE,GAAF,EAAE,CAAU;QAAS,OAAE,GAAF,EAAE,CAAU;QAIpD,eAAe;QACP,WAAM,GAAG,GAAG,CAAC;QAErB,qBAAqB;QACb,eAAU,GAAG,GAAG,CAAC;QAEjB,WAAM,GAAG,GAAG,CAAC;QATnB,IAAI,CAAC,MAAM,EAAE,CAAC;IAChB,CAAC;IAUD,eAAe;IACf,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,qBAAqB;IACrB,IAAI,SAAS;QACX,OAAO,IAAI,CAAC,UAAU,CAAC;IACzB,CAAC;IAED,yCAAyC;IACzC,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,MAAM,CAAC;IACrB,CAAC;IAED,gBAAgB;IAChB,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAClD,CAAC;IAEO,UAAU;QAChB,IAAI,KAAK,GAAG,GAAG,CAAC;QAEhB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACpC,MAAM,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YAErD,KAAK,IAAI,KAAK,GAAG,KAAK,CAAC;SACxB;QACD,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;IACrD,CAAC;IAED,2EAA2E;IAC3E,MAAM;QACJ,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,EAAE,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;QACnD,IAAI,GAAG,GAAG,GAAG,CAAC;QACd,IAAI,GAAG,GAAG,GAAG,CAAC;QAEd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,GAAG,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAClB,GAAG,IAAI,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;SACnB;QACD,GAAG,IAAI,CAAC,CAAC;QACT,GAAG,IAAI,CAAC,CAAC;QACT,IAAI,EAAE,GAAG,GAAG,CAAC;QACb,IAAI,IAAI,GAAG,GAAG,CAAC;QACf,IAAI,IAAI,GAAG,GAAG,CAAC;QAEf,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAC5B,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;YAE5B,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC;YACd,IAAI,IAAI,EAAE,GAAG,EAAE,CAAC;YAChB,IAAI,IAAI,EAAE,GAAG,EAAE,CAAC;SACjB;QACD,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC;QAEnD,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACjC,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC;QAChC,IAAI,CAAC,UAAU,GAAG,GAAG,GAAG,IAAI,CAAC,MAAM,GAAG,GAAG,CAAC;QAC1C,IAAI,CAAC,UAAU,EAAE,CAAC;IACpB,CAAC;IAED,wDAAwD;IACxD,QAAQ,CAAC,CAAS;QAChB,OAAO,IAAI,CAAC,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,UAAU,CAAC;IAC3C,CAAC;IAED;;;OAGG;IACH,cAAc,CAAC,KAAK,GAAG,GAAG;QACxB,MAAM,KAAK,GAAG,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACjC,MAAM,KAAK,GAAa,EAAE,CAAC;QAC3B,MAAM,KAAK,GAAa,EAAE,CAAC;QAE3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACpC,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,EAAE;gBAC5D,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvB,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;aACxB;SACF;QAED,OAAO,IAAI,sBAAsB,CAAC,KAAK,EAAE,KAAK,CAAC,CAAC;IAClD,CAAC;CACF"}

package/lib/optimize/SimplexEntry.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/>.
+ */
+/**
+ * Optimization cost function that takes in parameter array [xs] and returns
+ * a score.
+ */
+export type CostFunction = (points: Float64Array) => number;
+export declare class SimplexEntry {
+    private _f;
+    points: Float64Array;
+    score: number;
+    private _x;
+    constructor(_f: CostFunction, points: Float64Array);
+    getPoints(): Float64Array;
+    getScore(): number;
+    modify(n: number, xa: SimplexEntry, xb: SimplexEntry): SimplexEntry;
+    distance(se: SimplexEntry): number;
+}

package/lib/optimize/SimplexEntry.js

@@ -0,0 +1,39 @@
+/**
+ * @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 { Vector } from 'ootk-core';
+export class SimplexEntry {
+    constructor(_f, points) {
+        this._f = _f;
+        this.points = points;
+        this._x = new Vector(points);
+        this.score = this._f(points);
+    }
+    getPoints() {
+        return this._x.toArray();
+    }
+    getScore() {
+        return this.score;
+    }
+    modify(n, xa, xb) {
+        return new SimplexEntry(this._f, this._x.add(xa._x.subtract(xb._x).scale(n)).toArray());
+    }
+    distance(se) {
+        return this._x.distance(se._x);
+    }
+}
+//# sourceMappingURL=SimplexEntry.js.map

package/lib/optimize/SimplexEntry.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"SimplexEntry.js","sourceRoot":"","sources":["../../src/optimize/SimplexEntry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAQnC,MAAM,OAAO,YAAY;IAIvB,YAAoB,EAAgB,EAAS,MAAoB;QAA7C,OAAE,GAAF,EAAE,CAAc;QAAS,WAAM,GAAN,MAAM,CAAc;QAC/D,IAAI,CAAC,EAAE,GAAG,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;QAC7B,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAC/B,CAAC;IAED,SAAS;QACP,OAAO,IAAI,CAAC,EAAE,CAAC,OAAO,EAAE,CAAC;IAC3B,CAAC;IAED,QAAQ;QACN,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;IAED,MAAM,CAAC,CAAS,EAAE,EAAgB,EAAE,EAAgB;QAClD,OAAO,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC,CAAC;IAC1F,CAAC;IAED,QAAQ,CAAC,EAAgB;QACvB,OAAO,IAAI,CAAC,EAAE,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC;IACjC,CAAC;CACF"}

package/lib/orbit_determination/BatchLeastSquaresOD.d.ts

@@ -0,0 +1,54 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { J2000 } from 'ootk-core';
+import { Observation } from '../observation/Observation';
+import { ForceModel } from './../force/ForceModel';
+import { BatchLeastSquaresResult } from './BatchLeastSquaresResult';
+/**
+ * Batch least squares orbit determination.
+ */
+export declare class BatchLeastSquaresOD {
+    private observations_;
+    private apriori_;
+    private forceModel_?;
+    private posStep_;
+    private velStep_;
+    private fastDerivatives_;
+    private propPairs_;
+    private propagator_;
+    private nominal_;
+    private start_;
+    /**
+     * Create a new [BatchLeastSquaresOD] object from a list of [Observation]
+     * objects, an [apriori] state estimate, and an optional
+     * spacecraft [forceModel].
+     */
+    constructor(observations_: Observation[], apriori_: J2000, forceModel_?: ForceModel, posStep_?: number, velStep_?: number, fastDerivatives_?: boolean);
+    private buildPropagator_;
+    private static stateToX0_;
+    private setPropagatorPairs_;
+    /**
+     * Attempt to solve a state estimate with the given root-mean-squared delta
+     * [tolerance].
+     */
+    solve({ tolerance, maxIter, printIter, }?: {
+        tolerance?: number;
+        maxIter?: number;
+        printIter?: boolean;
+    }): BatchLeastSquaresResult;
+}

package/lib/orbit_determination/BatchLeastSquaresOD.js

@@ -0,0 +1,118 @@
+/**
+ * @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 { concat, J2000, Matrix, Vector, Vector3D } from 'ootk-core';
+import { PropagatorPairs } from '../observation/PropagatorPairs';
+import { CovarianceFrame, StateCovariance } from './../covariance/StateCovariance';
+import { ForceModel } from './../force/ForceModel';
+import { KeplerPropagator } from './../propagator/KeplerPropagator';
+import { RungeKutta89Propagator } from './../propagator/RungeKutta89Propagator';
+import { BatchLeastSquaresResult } from './BatchLeastSquaresResult';
+/**
+ * Batch least squares orbit determination.
+ */
+export class BatchLeastSquaresOD {
+    /**
+     * Create a new [BatchLeastSquaresOD] object from a list of [Observation]
+     * objects, an [apriori] state estimate, and an optional
+     * spacecraft [forceModel].
+     */
+    constructor(observations_, apriori_, forceModel_, posStep_ = 1e-5, velStep_ = 1e-5, fastDerivatives_ = false) {
+        this.observations_ = observations_;
+        this.apriori_ = apriori_;
+        this.forceModel_ = forceModel_;
+        this.posStep_ = posStep_;
+        this.velStep_ = velStep_;
+        this.fastDerivatives_ = fastDerivatives_;
+        this.observations_.sort((a, b) => a.epoch.posix - b.epoch.posix);
+        this.start_ = this.observations_[0].epoch;
+        this.propPairs_ = new PropagatorPairs(this.posStep_, this.velStep_);
+        this.forceModel_ ??= new ForceModel().setGravity();
+        this.propagator_ = new RungeKutta89Propagator(this.apriori_, this.forceModel_);
+        this.nominal_ = this.propagator_.propagate(this.start_);
+    }
+    buildPropagator_(x0, simple) {
+        const state = new J2000(this.nominal_.epoch, new Vector3D(x0[0], x0[1], x0[2]), new Vector3D(x0[3], x0[4], x0[5]));
+        if (simple) {
+            return new KeplerPropagator(state.toClassicalElements());
+        }
+        return new RungeKutta89Propagator(state, this.forceModel_);
+    }
+    static stateToX0_(state) {
+        return concat(state.position.toArray(), state.velocity.toArray());
+    }
+    setPropagatorPairs_(x0) {
+        const pl = this.buildPropagator_(x0, this.fastDerivatives_);
+        for (let i = 0; i < 6; i++) {
+            const step = this.propPairs_.step(i);
+            const xh = x0.slice();
+            xh[i] += step;
+            const ph = this.buildPropagator_(xh, this.fastDerivatives_);
+            this.propPairs_.set(i, ph, pl);
+        }
+    }
+    /**
+     * Attempt to solve a state estimate with the given root-mean-squared delta
+     * [tolerance].
+     */
+    solve({ tolerance = 1e-3, maxIter = 250, printIter = false, } = {}) {
+        let breakFlag = false;
+        const xNom = BatchLeastSquaresOD.stateToX0_(this.nominal_);
+        let weightedRms = Infinity;
+        const atwaMatInit = Matrix.zero(6, 6);
+        const atwbMatInit = Matrix.zero(6, 1);
+        let atwaMat = atwaMatInit;
+        for (let iter = 0; iter < maxIter; iter++) {
+            atwaMat = atwaMatInit;
+            let atwbMat = atwbMatInit;
+            this.propagator_ = this.buildPropagator_(xNom, false);
+            this.setPropagatorPairs_(xNom);
+            let rmsTotal = 0.0;
+            let measCount = 0;
+            for (const ob of this.observations_) {
+                const noise = ob.noise;
+                const aMat = ob.jacobian(this.propPairs_);
+                const aMatTN = aMat.transpose().multiply(noise);
+                const bMat = ob.residual(this.propagator_);
+                atwaMat = atwaMat.add(aMatTN.multiply(aMat));
+                atwbMat = atwbMat.add(aMatTN.multiply(bMat));
+                rmsTotal += bMat.transpose().multiply(noise).multiply(bMat)[0][0];
+                measCount += noise.rows;
+            }
+            const newWeightedRms = Math.sqrt(rmsTotal / measCount);
+            if (printIter) {
+                // eslint-disable-next-line no-console
+                console.log(`${iter + 1}: rms=${newWeightedRms} x=${new Vector(xNom)}`);
+            }
+            if (Math.abs((weightedRms - newWeightedRms) / weightedRms) <= tolerance) {
+                breakFlag = true;
+            }
+            weightedRms = newWeightedRms;
+            const dX = atwaMat.inverse().multiply(atwbMat);
+            for (let i = 0; i < 6; i++) {
+                xNom[i] += dX[i][0];
+            }
+            if (breakFlag) {
+                break;
+            }
+        }
+        const p = atwaMat.inverse();
+        const covariance = new StateCovariance(p, CovarianceFrame.ECI);
+        return new BatchLeastSquaresResult(this.buildPropagator_(xNom, false).propagate(this.start_), covariance, weightedRms);
+    }
+}
+//# sourceMappingURL=BatchLeastSquaresOD.js.map

package/lib/orbit_determination/BatchLeastSquaresOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BatchLeastSquaresOD.js","sourceRoot":"","sources":["../../src/orbit_determination/BatchLeastSquaresOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAY,KAAK,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,MAAM,WAAW,CAAC;AAE9E,OAAO,EAAE,eAAe,EAAE,MAAM,gCAAgC,CAAC;AAEjE,OAAO,EAAE,eAAe,EAAE,eAAe,EAAE,MAAM,iCAAiC,CAAC;AACnF,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,gBAAgB,EAAE,MAAM,kCAAkC,CAAC;AACpE,OAAO,EAAE,sBAAsB,EAAE,MAAM,wCAAwC,CAAC;AAChF,OAAO,EAAE,uBAAuB,EAAE,MAAM,2BAA2B,CAAC;AAEpE;;GAEG;AACH,MAAM,OAAO,mBAAmB;IAU9B;;;;OAIG;IACH,YACU,aAA4B,EAC5B,QAAe,EACf,WAAwB,EACxB,WAAmB,IAAI,EACvB,WAAmB,IAAI,EACvB,mBAA4B,KAAK;QALjC,kBAAa,GAAb,aAAa,CAAe;QAC5B,aAAQ,GAAR,QAAQ,CAAO;QACf,gBAAW,GAAX,WAAW,CAAa;QACxB,aAAQ,GAAR,QAAQ,CAAe;QACvB,aAAQ,GAAR,QAAQ,CAAe;QACvB,qBAAgB,GAAhB,gBAAgB,CAAiB;QAEzC,IAAI,CAAC,aAAa,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;QACjE,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC;QAC1C,IAAI,CAAC,UAAU,GAAG,IAAI,eAAe,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,QAAQ,CAAC,CAAC;QACpE,IAAI,CAAC,WAAW,KAAK,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;QACnD,IAAI,CAAC,WAAW,GAAG,IAAI,sBAAsB,CAAC,IAAI,CAAC,QAAQ,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;QAC/E,IAAI,CAAC,QAAQ,GAAG,IAAI,CAAC,WAAW,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC1D,CAAC;IAEO,gBAAgB,CAAC,EAAgB,EAAE,MAAe;QACxD,MAAM,KAAK,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,EAAE,IAAI,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,QAAQ,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEnH,IAAI,MAAM,EAAE;YACV,OAAO,IAAI,gBAAgB,CAAC,KAAK,CAAC,mBAAmB,EAAE,CAAC,CAAC;SAC1D;QAED,OAAO,IAAI,sBAAsB,CAAC,KAAK,EAAE,IAAI,CAAC,WAAW,CAAC,CAAC;IAC7D,CAAC;IAEO,MAAM,CAAC,UAAU,CAAC,KAAY;QACpC,OAAO,MAAM,CAAC,KAAK,CAAC,QAAQ,CAAC,OAAO,EAAE,EAAE,KAAK,CAAC,QAAQ,CAAC,OAAO,EAAE,CAAC,CAAC;IACpE,CAAC;IAEO,mBAAmB,CAAC,EAAgB;QAC1C,MAAM,EAAE,GAAG,IAAI,CAAC,gBAAgB,CAAC,EAAE,EAAE,IAAI,CAAC,gBAAgB,CAAC,CAAC;QAE5D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,IAAI,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;YACrC,MAAM,EAAE,GAAG,EAAE,CAAC,KAAK,EAAE,CAAC;YAEtB,EAAE,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC;YACd,MAAM,EAAE,GAAG,IAAI,CAAC,gBAAgB,CAAC,EAAE,EAAE,IAAI,CAAC,gBAAgB,CAAC,CAAC;YAE5D,IAAI,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;SAChC;IACH,CAAC;IAED;;;OAGG;IACH,KAAK,CAAC,EACJ,SAAS,GAAG,IAAI,EAChB,OAAO,GAAG,GAAG,EACb,SAAS,GAAG,KAAK,MAKf,EAAE;QACJ,IAAI,SAAS,GAAG,KAAK,CAAC;QACtB,MAAM,IAAI,GAAG,mBAAmB,CAAC,UAAU,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;QAC3D,IAAI,WAAW,GAAG,QAAQ,CAAC;QAC3B,MAAM,WAAW,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtC,MAAM,WAAW,GAAG,MAAM,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtC,IAAI,OAAO,GAAG,WAAW,CAAC;QAE1B,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,OAAO,EAAE,IAAI,EAAE,EAAE;YACzC,OAAO,GAAG,WAAW,CAAC;YACtB,IAAI,OAAO,GAAG,WAAW,CAAC;YAE1B,IAAI,CAAC,WAAW,GAAG,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC;YACtD,IAAI,CAAC,mBAAmB,CAAC,IAAI,CAAC,CAAC;YAC/B,IAAI,QAAQ,GAAG,GAAG,CAAC;YACnB,IAAI,SAAS,GAAG,CAAC,CAAC;YAElB,KAAK,MAAM,EAAE,IAAI,IAAI,CAAC,aAAa,EAAE;gBACnC,MAAM,KAAK,GAAG,EAAE,CAAC,KAAK,CAAC;gBACvB,MAAM,IAAI,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,UAAU,CAAC,CAAC;gBAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;gBAChD,MAAM,IAAI,GAAG,EAAE,CAAC,QAAQ,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;gBAE3C,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC;gBAC7C,OAAO,GAAG,OAAO,CAAC,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC;gBAC7C,QAAQ,IAAI,IAAI,CAAC,SAAS,EAAE,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClE,SAAS,IAAI,KAAK,CAAC,IAAI,CAAC;aACzB;YACD,MAAM,cAAc,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,SAAS,CAAC,CAAC;YAEvD,IAAI,SAAS,EAAE;gBACb,sCAAsC;gBACtC,OAAO,CAAC,GAAG,CAAC,GAAG,IAAI,GAAG,CAAC,SAAS,cAAc,MAAM,IAAI,MAAM,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;aACzE;YACD,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,WAAW,GAAG,cAAc,CAAC,GAAG,WAAW,CAAC,IAAI,SAAS,EAAE;gBACvE,SAAS,GAAG,IAAI,CAAC;aAClB;YACD,WAAW,GAAG,cAAc,CAAC;YAC7B,MAAM,EAAE,GAAG,OAAO,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC,OAAO,CAAC,CAAC;YAE/C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;gBAC1B,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;aACrB;YACD,IAAI,SAAS,EAAE;gBACb,MAAM;aACP;SACF;QACD,MAAM,CAAC,GAAG,OAAO,CAAC,OAAO,EAAE,CAAC;QAC5B,MAAM,UAAU,GAAG,IAAI,eAAe,CAAC,CAAC,EAAE,eAAe,CAAC,GAAG,CAAC,CAAC;QAE/D,OAAO,IAAI,uBAAuB,CAChC,IAAI,CAAC,gBAAgB,CAAC,IAAI,EAAE,KAAK,CAAC,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,CAAC,EACzD,UAAU,EACV,WAAW,CACZ,CAAC;IACJ,CAAC;CACF"}

package/lib/orbit_determination/BatchLeastSquaresResult.d.ts

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

package/lib/orbit_determination/BatchLeastSquaresResult.js

@@ -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/>.
+ */
+// / Batch least squares orbit determination result.
+export class BatchLeastSquaresResult {
+    /**
+     * Create a new [BatchLeastSquaresResult] object, containing the solved
+     * [state], [covariance], and root-mean-squared error [rms].
+     */
+    constructor(state, covariance, rms) {
+        this.state = state;
+        this.covariance = covariance;
+        this.rms = rms;
+        // Nothing to do here.
+    }
+}
+//# sourceMappingURL=BatchLeastSquaresResult.js.map

package/lib/orbit_determination/BatchLeastSquaresResult.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BatchLeastSquaresResult.js","sourceRoot":"","sources":["../../src/orbit_determination/BatchLeastSquaresResult.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAKH,oDAAoD;AACpD,MAAM,OAAO,uBAAuB;IAClC;;;OAGG;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/lib/orbit_determination/GibbsIOD.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, 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.
+     */
+    solve(r1: Vector3D, r2: Vector3D, r3: Vector3D, t2: EpochUTC, t3: EpochUTC): J2000;
+}

package/lib/orbit_determination/GibbsIOD.js

@@ -0,0 +1,70 @@
+/**
+ * @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 {
+    constructor(mu = Earth.mu) {
+        this.mu = mu;
+        // Nothing to do here.
+    }
+    /**
+     * Attempt to create a state estimate from three inertial position vectors.
+     *
+     * Throws an error if the positions are not coplanar.
+     */
+    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);
+    }
+}
+/** Abort solve if position plane exceeds this value. */
+GibbsIOD.coplanarThreshold_ = (5.0 * DEG2RAD);
+//# sourceMappingURL=GibbsIOD.js.map

package/lib/orbit_determination/GibbsIOD.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"GibbsIOD.js","sourceRoot":"","sources":["../../src/orbit_determination/GibbsIOD.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,OAAO,EAAE,KAAK,EAAY,MAAM,EAAE,KAAK,EAAqB,MAAM,WAAW,CAAC;AACvF,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,sBAAsB,EAAE,MAAM,wCAAwC,CAAC;AAEhF;;GAEG;AACH,MAAM,OAAO,QAAQ;IACnB,YAAmB,KAAa,KAAK,CAAC,EAAE;QAArB,OAAE,GAAF,EAAE,CAAmB;QACtC,sBAAsB;IACxB,CAAC;IAKD;;;;OAIG;IACH,KAAK,CAAC,EAAY,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY,EAAE,EAAY;QACxE,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,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;QAEjD,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;;AApDD,wDAAwD;AAChC,2BAAkB,GAAY,CAAC,GAAG,GAAG,OAAO,CAAY,CAAC"}

package/lib/orbit_determination/GoodingIOD.d.ts

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