ootk 4.0.9 → 5.1.0

package/src/operations/Vector3D.ts

@@ -0,0 +1,289 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * 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 { Kilometers, KilometersPerSecond, Radians, linearDistance } from '../main.js';
+import { Matrix } from './Matrix.js';
+import { Vector } from './Vector.js';
+
+// / 3-dimensional vector.
+export class Vector3D<T extends number = number> {
+  constructor(
+    public x: T,
+    public y: T,
+    public z: T,
+  ) {
+    // Nothing to do here.
+  }
+
+  /**
+   * Create a new Vector3D object from the first three elements of a Vector
+   * object.
+   * @param v The Vector object to convert.
+   * @returns A new Vector3D object.
+   */
+  static fromVector<U extends number>(v: Vector<U>): Vector3D<U> {
+    return new Vector3D<U>(v.x as U, v.y as U, v.z as U);
+  }
+
+  // / Origin vector.
+  static readonly origin = new Vector3D<number>(0, 0, 0);
+
+  // / X-axis unit vector.
+  static readonly xAxis = new Vector3D<number>(1, 0, 0);
+
+  // / Y-axis unit vector.
+  static readonly yAxis = new Vector3D<number>(0, 1, 0);
+
+  // / Z-axis unit vector.
+  static readonly zAxis = new Vector3D<number>(0, 0, 1);
+
+  // / Negative x-axis unit vector.
+  static readonly xAxisNeg = new Vector3D<number>(-1, 0, 0);
+
+  // / Negative y-axis unit vector.
+  static readonly yAxisNeg = new Vector3D<number>(0, -1, 0);
+
+  // / Negative z-axis unit vector.
+  static readonly zAxisNeg = new Vector3D<number>(0, 0, -1);
+
+  // / Convert this to a [List] of doubles.
+  toList() {
+    return [this.x, this.y, this.z];
+  }
+
+  // / Convert this to a [Float64List] object.
+  toArray() {
+    return new Float64Array([this.x, this.y, this.z]);
+  }
+
+  /**
+   * Return the Vector3D element at the provided index.
+   * @deprecated don't do this
+   * @param index The index of the element to return.
+   * @returns The element at the provided index.
+   */
+  public getElement(index: number): number {
+    switch (index) {
+      case 0:
+        return this.x;
+      case 1:
+        return this.y;
+      case 2:
+        return this.z;
+      default:
+        throw new Error(`Index ${index} outside 3D vector bounds.`);
+    }
+  }
+
+  // / Convert this to a [Vector] object.
+  toVector() {
+    return new Vector(this.toList());
+  }
+
+  toString(fixed = -1): string {
+    if (fixed < 0) {
+      return `[${this.toList().join(', ')}]`;
+    }
+    const output = this.toList().map((e) => e.toFixed(fixed));
+
+    return `[${output.join(', ')}]`;
+  }
+
+  // / Return the magnitude of this vector.
+  magnitude(): T {
+    return Math.sqrt(this.x * this.x + this.y * this.y + this.z * this.z) as T;
+  }
+
+  // / Return the result of adding this to another [Vector3D].
+  add(v: Vector3D<T>): Vector3D<T> {
+    return new Vector3D<T>((this.x + v.x) as T, (this.y + v.y) as T, (this.z + v.z) as T);
+  }
+
+  // / Return the result of subtracting this and another [Vector3D].
+  subtract(v: Vector3D<T>): Vector3D<T> {
+    return new Vector3D<T>((this.x - v.x) as T, (this.y - v.y) as T, (this.z - v.z) as T);
+  }
+
+  // / Return a copy of this [Vector3D] scaled by [n];
+  scale<U extends number>(n: U): Vector3D<U> {
+    return new Vector3D<U>(this.x * n as U, this.y * n as U, this.z * n as U);
+  }
+
+  // / Return a copy of this [Vector3D] with the elements negated.
+  negate(): Vector3D<T> {
+    return this.scale(-1 as T);
+  }
+
+  /**
+   * Return the Euclidean distance between this and another Vector3D.
+   * @param v The other Vector3D.
+   * @returns The distance between this and the other Vector3D.
+   */
+  distance(v: Vector3D<T>): T {
+    return linearDistance(this, v);
+  }
+
+  /**
+   * Convert this to a unit Vector3D.
+   * @returns A unit Vector3D.
+   */
+  normalize(): Vector3D<T> {
+    const m = this.magnitude();
+
+    if (m === 0) {
+      return Vector3D.origin as Vector3D<T>;
+    }
+
+    return new Vector3D<T>(this.x / m as T, this.y / m as T, this.z / m as T);
+  }
+
+  // Calculate the dot product of this and another [Vector3D].
+  dot<T extends number>(v: Vector3D<T>): T {
+    return this.x * v.x + this.y * v.y + this.z * v.z as T;
+  }
+
+  // Calculate the outer product between this and another [Vector3D].
+  outer(v: Vector3D): Matrix {
+    return new Matrix([
+      [this.x * v.x, this.x * v.y, this.x * v.z],
+      [this.y * v.x, this.y * v.y, this.y * v.z],
+      [this.z * v.x, this.z * v.y, this.z * v.z],
+    ]);
+  }
+
+  // Calculate the cross product of this and another [Vector3D].
+  cross<U extends number>(v: Vector3D<U>): Vector3D<U> {
+    return new Vector3D<U>(
+      (this.y * v.z - this.z * v.y) as U,
+      (this.z * v.x - this.x * v.z) as U,
+      (this.x * v.y - this.y * v.x) as U,
+    );
+  }
+
+  // Calculate the skew-symmetric matrix for this [Vector3D].
+  skewSymmetric(): Matrix {
+    return new Matrix([
+      [0, -this.z, this.y],
+      [this.z, 0, -this.x],
+      [-this.y, this.x, 0],
+    ]);
+  }
+
+  /*
+   * Create a copy of this [Vector3D] rotated in the x-axis by angle [theta]
+   * _(rad)_.
+   */
+  rotX(theta: number): Vector3D {
+    const cosT = Math.cos(theta);
+    const sinT = Math.sin(theta);
+
+    return new Vector3D(this.x, cosT * this.y + sinT * this.z, -sinT * this.y + cosT * this.z);
+  }
+
+  /*
+   * Create a copy of this [Vector3D] rotated in the y-axis by angle [theta]
+   * _(rad)_.
+   */
+  rotY(theta: Radians): Vector3D<T> {
+    const cosT = Math.cos(theta);
+    const sinT = Math.sin(theta);
+
+    return new Vector3D<T>((cosT * this.x + -sinT * this.z) as T, this.y, (sinT * this.x + cosT * this.z) as T);
+  }
+
+  /*
+   * Create a copy of this [Vector3D] rotated in the z-axis by angle [theta]
+   * _(rad)_.
+   */
+  rotZ(theta: Radians): Vector3D<T> {
+    const cosT = Math.cos(theta);
+    const sinT = Math.sin(theta);
+
+    return new Vector3D((cosT * this.x + sinT * this.y) as T, (-sinT * this.x + cosT * this.y) as T, this.z);
+  }
+
+  // Calculate the angle _(rad)_ between this and another [Vector3D].
+  angle<U extends number>(v: Vector3D<U>): Radians {
+    const theta = Math.atan2(this.cross(v).magnitude(), this.dot(v)) as Radians;
+
+    return isNaN(theta) ? 0 as Radians : theta;
+  }
+
+  // Calculate the angle _(°)_ between this and another [Vector3D].
+  angleDegrees(v: Vector3D<T>): number {
+    return this.angle(v) * (180 / Math.PI);
+  }
+
+  /*
+   * Return `true` if line-of-sight exists between this and another [Vector3D]
+   * with a central body of the given [radius].
+   */
+  sight(v: Vector3D<KilometersPerSecond>, radius: Kilometers): boolean {
+    const r1Mag2 = this.magnitude() ** 2;
+    const r2Mag2 = v.magnitude() ** 2;
+    const rDot = this.dot(v);
+    const tMin = (r1Mag2 - rDot) / (r1Mag2 + r2Mag2 - 2 * rDot);
+    let los = false;
+
+    if (tMin < 0 || tMin > 1) {
+      los = true;
+    } else {
+      const c = (1 - tMin) * r1Mag2 + rDot * tMin;
+
+      if (c >= radius * radius) {
+        los = true;
+      }
+    }
+
+    return los;
+  }
+
+  // / Return the unit vector that bisects this and another [Vector3D].
+  bisect(v: Vector3D<T>): Vector3D<T> {
+    return this.scale(v.magnitude()).add(v.scale(this.magnitude())).normalize();
+  }
+
+  // / Convert this [Vector3D] into a row [Matrix].
+  row(): Matrix {
+    return new Matrix([[this.x, this.y, this.z]]);
+  }
+
+  // / Convert this [Vector3D] into a column [Matrix].
+  column(): Matrix {
+    return new Matrix([[this.x], [this.y], [this.z]]);
+  }
+
+  // / Join this and another [Vector3D] into a new [Vector] object.
+  join(v: Vector3D): Vector {
+    const output = new Float64Array(6);
+
+    output[0] = this.x;
+    output[1] = this.y;
+    output[2] = this.z;
+    output[3] = v.x;
+    output[4] = v.y;
+    output[5] = v.z;
+
+    return new Vector(output);
+  }
+}

package/src/operations/operations.ts

@@ -0,0 +1,27 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @description Orbital Object ToolKit (ootk) is a collection of tools for working
+ * with satellites and other orbital objects.
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2025 Kruczek Labs LLC
+ *
+ * Many of the classes are based off of the work of @david-rc-dayton and his
+ * Pious Squid library (https://github.com/david-rc-dayton/pious_squid) which
+ * is licensed under the MIT license.
+ *
+ * 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/>.
+ */
+
+export { Matrix } from './Matrix.js';
+export { Random } from './Random.js';
+export { Vector } from './Vector.js';
+export { Vector3D } from './Vector3D.js';

package/src/optimize/ChebyshevCompressor.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { EpochUTC, Seconds, Vector3D } from 'ootk-core';
+import { EpochUTC, Seconds, Vector3D } from '../main.js';
 import { StateInterpolator } from '../interpolator/StateInterpolator.js';
 import { ChebyshevCoefficients } from './../interpolator/ChebyshevCoefficients.js';
 import { ChebyshevInterpolator } from './../interpolator/ChebyshevInterpolator.js';
@@ -23,7 +23,7 @@ import { ChebyshevInterpolator } from './../interpolator/ChebyshevInterpolator.j
 // / Ephemeris compressor.
 export class ChebyshevCompressor {
   // / Create a new [ChebyshevCompressor] object from an [Interpolator].
-  constructor(private _interpolator: StateInterpolator) {
+  constructor(private readonly interpolator_: StateInterpolator) {
     // Do nothing.
   }
 
@@ -41,7 +41,7 @@ export class ChebyshevCompressor {
     for (let i = 0; i < n; i++) {
       const x = ChebyshevCompressor._cosPi((i + h) / n);
       const seconds = x * (h * (b - a)) + h * (b + a) as Seconds;
-      const state = this._interpolator.interpolate(new EpochUTC(seconds))!;
+      const state = this.interpolator_.interpolate(new EpochUTC(seconds))!;
       const fx = state.position.x;
       const fy = state.position.y;
       const fz = state.position.z;
@@ -78,8 +78,8 @@ export class ChebyshevCompressor {
    * @returns A new [ChebyshevInterpolator] object.
    */
   compress(cpr = 21): ChebyshevInterpolator {
-    const { start, end } = this._interpolator.window();
-    const period = this._interpolator.interpolate(start)!.period;
+    const { start, end } = this.interpolator_.window();
+    const period = this.interpolator_.interpolate(start)!.period;
     const coefficients: ChebyshevCoefficients[] = [];
     let current = start;
 

package/src/optimize/DownhillSimplex.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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

package/src/optimize/GoldenSection.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -21,7 +21,7 @@ import { DifferentiableFunction } from '../main.js';
 export class GoldenSection {
   private static readonly _grInv: number = 1.0 / (0.5 * (Math.sqrt(5) + 1));
 
-  private static _check(fc: number, fd: number, solveMax: boolean): boolean {
+  private static check_(fc: number, fd: number, solveMax: boolean): boolean {
     return solveMax ? fc > fd : fc < fd;
   }
 
@@ -57,7 +57,7 @@ export class GoldenSection {
     let d = a + (b - a) * GoldenSection._grInv;
 
     while (Math.abs(b - a) > tolerance) {
-      if (GoldenSection._check(f(c), f(d), solveMax)) {
+      if (GoldenSection.check_(f(c), f(d), solveMax)) {
         b = d;
       } else {
         a = c;

package/src/optimize/PolynomialRegression.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
  * terms of the GNU Affero General Public License as published by the Free Software
@@ -17,7 +17,7 @@
 
 /* eslint-disable require-jsdoc */
 /* eslint-disable func-style */
-import { evalPoly } from 'ootk-core';
+import { evalPoly } from '../main.js';
 import { DownhillSimplex } from './DownhillSimplex.js';
 import { PolynomicalRegressionResult } from './PolynomicalRegressionResult.js';
 // / Polynomial regression optimizer.

package/src/optimize/PolynomicalRegressionResult.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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

package/src/optimize/SimpleLinearRegression.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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

package/src/optimize/SimplexEntry.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Vector } from 'ootk-core';
+import { Vector } from '../main.js';
 
 /**
  * Optimization cost function that takes in parameter array [xs] and returns
@@ -25,15 +25,15 @@ export type CostFunction = (points: Float64Array) => number;
 
 export class SimplexEntry {
   score: number;
-  private _x: Vector;
+  private readonly x_: Vector;
 
-  constructor(private _f: CostFunction, public points: Float64Array) {
-    this._x = new Vector(points);
-    this.score = this._f(points);
+  constructor(private readonly f_: CostFunction, public points: Float64Array) {
+    this.x_ = new Vector(points);
+    this.score = this.f_(points);
   }
 
   getPoints(): Float64Array {
-    return this._x.toArray();
+    return this.x_.toArray();
   }
 
   getScore(): number {
@@ -41,10 +41,10 @@ export class SimplexEntry {
   }
 
   modify(n: number, xa: SimplexEntry, xb: SimplexEntry): SimplexEntry {
-    return new SimplexEntry(this._f, this._x.add(xa._x.subtract(xb._x).scale(n)).toArray());
+    return new SimplexEntry(this.f_, this.x_.add(xa.x_.subtract(xb.x_).scale(n)).toArray());
   }
 
   distance(se: SimplexEntry): number {
-    return this._x.distance(se._x);
+    return this.x_.distance(se.x_);
   }
 }

package/src/orbit_determination/BatchLeastSquaresOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { concat, EpochUTC, J2000, Kilometers, KilometersPerSecond, Matrix, Vector, Vector3D } from 'ootk-core';
+import { concat, EpochUTC, J2000, Kilometers, KilometersPerSecond, Matrix, Vector, Vector3D } from '../main.js';
 import { Observation } from '../observation/Observation.js';
 import { PropagatorPairs } from '../observation/PropagatorPairs.js';
 import { Propagator } from '../propagator/Propagator.js';
@@ -30,13 +30,13 @@ import { BatchLeastSquaresResult } from './BatchLeastSquaresResult.js';
  */
 export class BatchLeastSquaresOD {
   /** Propagator pair cache, for generating observation Jacobians. */
-  private propPairs_: PropagatorPairs;
+  private readonly propPairs_: PropagatorPairs;
   /**  Nominal state propagator. */
   private propagator_: Propagator;
   /**  State estimate during solve. */
-  private nominal_: J2000;
+  private readonly nominal_: J2000;
   /**  Solve start epoch. */
-  private start_: EpochUTC;
+  private readonly start_: EpochUTC;
 
   /**
    * Create a new [BatchLeastSquaresOD] object from a list of [Observation]
@@ -51,12 +51,12 @@ export class BatchLeastSquaresOD {
    * @returns [BatchLeastSquaresOD] object.
    */
   constructor(
-    private observations_: Observation[],
-    private apriori_: J2000,
-    private forceModel_?: ForceModel,
-    private posStep_: number = 1e-5,
-    private velStep_: number = 1e-5,
-    private fastDerivatives_: boolean = false,
+    private readonly observations_: Observation[],
+    private readonly apriori_: J2000,
+    private readonly forceModel_?: ForceModel,
+    private readonly posStep_: number = 1e-5,
+    private readonly velStep_: number = 1e-5,
+    private readonly fastDerivatives_: boolean = false,
   ) {
     this.observations_.sort((a, b) => a.epoch.posix - b.epoch.posix);
     this.start_ = this.observations_[0].epoch;

package/src/orbit_determination/BatchLeastSquaresResult.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { J2000 } from 'ootk-core';
+import { J2000 } from '../main.js';
 import { StateCovariance } from './../covariance/StateCovariance.js';
 
 // / Batch least squares orbit determination result.

package/src/orbit_determination/GibbsIOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,9 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { DEG2RAD, Earth, EpochUTC, halfPi, J2000, Kilometers, KilometersPerSecond, Radians, Vector3D } from 'ootk-core';
+import {
+  DEG2RAD, Earth, EpochUTC, halfPi, J2000, Kilometers, KilometersPerSecond, Radians, Vector3D,
+} from '../main.js';
 import { ForceModel } from './../force/ForceModel.js';
 import { RungeKutta89Propagator } from './../propagator/RungeKutta89Propagator.js';
 

package/src/orbit_determination/GoodingIOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Earth, J2000, Kilometers, KilometersPerSecond, Seconds, TAU, Vector3D } from 'ootk-core';
+import { Earth, J2000, Kilometers, KilometersPerSecond, Seconds, TAU, Vector3D } from '../main.js';
 import { ForceModel } from '../force/ForceModel.js';
 import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator.js';
 import { ObservationOptical } from './../observation/ObservationOptical.js';

package/src/orbit_determination/HerrickGibbsIOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, Vector3D } from 'ootk-core';
+import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, Vector3D } from '../main.js';
 
 /**
  * Herrik-Gibbs 3-position initial orbit determination.

package/src/orbit_determination/LambertIOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, Radians, Vector3D } from 'ootk-core';
+import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, Radians, Vector3D } from '../main.js';
 
 // / Lambert two-position and time initial orbit determination.
 export class LambertIOD {

package/src/orbit_determination/ModifiedGoodingIOD.ts

@@ -1,7 +1,7 @@
 /**
  * @author @thkruz Theodore Kruczek
  * @license AGPL-3.0-or-later
- * @copyright (c) 2020-2024 Theodore Kruczek
+ * @copyright (c) 2025 Kruczek Labs LLC
  *
  * 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
@@ -15,7 +15,7 @@
  * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
  */
 
-import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, RadecTopocentric, Vector3D } from 'ootk-core';
+import { Earth, EpochUTC, J2000, Kilometers, KilometersPerSecond, RadecTopocentric, Vector3D } from '../main.js';
 import { ForceModel } from '../force/ForceModel.js';
 import { RungeKutta89Propagator } from '../propagator/RungeKutta89Propagator.js';
 import { ObservationOptical } from './../observation/ObservationOptical.js';
@@ -38,8 +38,8 @@ type SolveOptions = {
  * Used for orbit determination from multiple optical observations.
  */
 export class ModifiedGoodingIOD {
-  private observations_: ObservationOptical[];
-  private mu_: number;
+  private readonly observations_: ObservationOptical[];
+  private readonly mu_: number;
 
   constructor(observations: ObservationOptical[], mu: number = Earth.mu) {
     this.observations_ = observations;

package/src/orbit_determination/index.ts

@@ -0,0 +1,7 @@
+export { BatchLeastSquaresOD } from './BatchLeastSquaresOD.js';
+export { BatchLeastSquaresResult } from './BatchLeastSquaresResult.js';
+export { GibbsIOD } from './GibbsIOD.js';
+export { GoodingIOD } from './GoodingIOD.js';
+export { HerrickGibbsIOD } from './HerrickGibbsIOD.js';
+export { LambertIOD } from './LambertIOD.js';
+export { ModifiedGoodingIOD } from './ModifiedGoodingIOD.js';