ootk 6.0.0-rc.4 → 6.0.0

package/src/interpolator/__tests__/LagrangeInterpolator.test.ts

@@ -1,121 +0,0 @@
-import { EpochUTC, J2000, Kilometers, KilometersPerSecond, Seconds, Vector3D } from '../../main';
-import { LagrangeInterpolator } from '../LagrangeInterpolator';
-
-describe('LagrangeInterpolator', () => {
-  let testEphemeris: J2000[];
-  let interpolator: LagrangeInterpolator;
-
-  beforeEach(() => {
-    // Create sample ephemeris data for testing
-    testEphemeris = [
-      new J2000(new EpochUTC(0 as Seconds), new Vector3D(1000, 2000, 3000) as Vector3D<Kilometers>, new Vector3D(1, 2, 3) as Vector3D<KilometersPerSecond>),
-      new J2000(new EpochUTC(10 as Seconds), new Vector3D(1100, 2100, 3100) as Vector3D<Kilometers>, new Vector3D(1.1, 2.1, 3.1) as Vector3D<KilometersPerSecond>),
-      new J2000(new EpochUTC(20 as Seconds), new Vector3D(1200, 2200, 3200) as Vector3D<Kilometers>, new Vector3D(1.2, 2.2, 3.2) as Vector3D<KilometersPerSecond>),
-      new J2000(new EpochUTC(30 as Seconds), new Vector3D(1300, 2300, 3300) as Vector3D<Kilometers>, new Vector3D(1.3, 2.3, 3.3) as Vector3D<KilometersPerSecond>),
-      new J2000(new EpochUTC(40 as Seconds), new Vector3D(1400, 2400, 3400) as Vector3D<Kilometers>, new Vector3D(1.4, 2.4, 3.4) as Vector3D<KilometersPerSecond>),
-      new J2000(new EpochUTC(50 as Seconds), new Vector3D(1500, 2500, 3500) as Vector3D<Kilometers>, new Vector3D(1.5, 2.5, 3.5) as Vector3D<KilometersPerSecond>),
-    ];
-
-    interpolator = LagrangeInterpolator.fromEphemeris(testEphemeris, 6);
-  });
-
-  describe('constructor', () => {
-    it('should create instance with provided arrays', () => {
-      const t = new Float64Array([0, 10, 20]);
-      const x = new Float64Array([1000, 1100, 1200]);
-      const y = new Float64Array([2000, 2100, 2200]);
-      const z = new Float64Array([3000, 3100, 3200]);
-
-      const interp = new LagrangeInterpolator(t, x, y, z, 3);
-
-      expect(interp).toBeInstanceOf(LagrangeInterpolator);
-    });
-  });
-
-  describe('fromEphemeris', () => {
-    it('should create interpolator from ephemeris array', () => {
-      expect(interpolator).toBeInstanceOf(LagrangeInterpolator);
-    });
-
-    it('should use default order of 10', () => {
-      const interp = LagrangeInterpolator.fromEphemeris(testEphemeris);
-
-      expect(interp).toBeInstanceOf(LagrangeInterpolator);
-    });
-  });
-
-  describe('sizeBytes', () => {
-    it('should calculate correct size in bytes', () => {
-      const expectedSize = (64 * 4 * testEphemeris.length) / 8;
-
-      expect(interpolator.sizeBytes).toBe(expectedSize);
-    });
-  });
-
-  describe('interpolate', () => {
-    it('should interpolate position at given epoch', () => {
-      const epoch = new EpochUTC(15 as Seconds);
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).not.toBeNull();
-      expect(result?.epoch).toBe(epoch);
-      expect(result?.position.x).toBeCloseTo(1150, 0);
-      expect(result?.position.y).toBeCloseTo(2150, 0);
-      expect(result?.position.z).toBeCloseTo(3150, 0);
-    });
-
-    it('should return null for epoch outside window', () => {
-      const epoch = new EpochUTC(100 as Seconds);
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).toBeNull();
-    });
-
-    it('should interpolate at exact ephemeris point', () => {
-      const epoch = new EpochUTC(20 as Seconds);
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).not.toBeNull();
-      expect(result?.position.x).toBeCloseTo(1200, 0);
-    });
-  });
-
-  describe('window', () => {
-    it('should return correct epoch window', () => {
-      const window = interpolator.window();
-
-      expect(window.start.posix).toBe(0);
-      expect(window.end.posix).toBe(50);
-    });
-  });
-
-  describe('edge cases', () => {
-    it('should handle interpolation at start of window', () => {
-      const epoch = new EpochUTC(0 as Seconds);
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).not.toBeNull();
-      expect(result?.position.x).toBeCloseTo(1000, 0);
-    });
-
-    it('should handle interpolation at end of window', () => {
-      const epoch = new EpochUTC(50 as Seconds);
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).not.toBeNull();
-      expect(result?.position.x).toBeCloseTo(1500, 0);
-    });
-
-    it('should handle small ephemeris arrays', () => {
-      const smallEphemeris = [
-        new J2000(new EpochUTC(0 as Seconds), new Vector3D(1000, 2000, 3000) as Vector3D<Kilometers>, new Vector3D(1, 2, 3) as Vector3D<KilometersPerSecond>),
-        new J2000(new EpochUTC(10 as Seconds), new Vector3D(1100, 2100, 3100) as Vector3D<Kilometers>, new Vector3D(1.1, 2.1, 3.1) as Vector3D<KilometersPerSecond>),
-      ];
-
-      const smallInterpolator = LagrangeInterpolator.fromEphemeris(smallEphemeris, 2);
-      const result = smallInterpolator.interpolate(new EpochUTC(5 as Seconds));
-
-      expect(result).not.toBeNull();
-    });
-  });
-});

package/src/interpolator/__tests__/StateInterpolator.test.ts

@@ -1,57 +0,0 @@
-import { EpochUTC, EpochWindow, J2000, Seconds } from '../../main';
-import { StateInterpolator } from '../StateInterpolator';
-
-// Concrete implementation for testing
-class TestStateInterpolator extends StateInterpolator {
-  private start_: EpochUTC;
-  private end_: EpochUTC;
-
-  constructor() {
-    super();
-    this.start_ = new EpochUTC(0 as Seconds);
-    this.end_ = new EpochUTC(86400 as Seconds);
-  }
-
-  override window(): EpochWindow {
-    return new EpochWindow(this.start_, this.end_);
-  }
-
-  override interpolate(epoch: EpochUTC): J2000 | null {
-    return null;
-  }
-}
-
-describe('StateInterpolator', () => {
-  let interpolator: TestStateInterpolator;
-
-  beforeEach(() => {
-    interpolator = new TestStateInterpolator();
-  });
-
-  describe('interpolate', () => {
-    it('should be an abstract method that must be implemented', () => {
-      const epoch = EpochUTC.fromDateTime(new Date());
-
-      expect(() => interpolator.interpolate(epoch)).not.toThrow();
-    });
-
-    it('should return null in test implementation', () => {
-      const epoch = EpochUTC.fromDateTime(new Date());
-      const result = interpolator.interpolate(epoch);
-
-      expect(result).toBeNull();
-    });
-  });
-
-  describe('sizeBytes', () => {
-    it('should throw "Not implemented" error', () => {
-      expect(() => interpolator.sizeBytes).toThrow('Not implemented.');
-    });
-  });
-
-  describe('inheritance', () => {
-    it('should extend Interpolator', () => {
-      expect(interpolator).toBeInstanceOf(StateInterpolator);
-    });
-  });
-});

package/src/interpolator/__tests__/VerletBlendInterpolator.test.ts

@@ -1,116 +0,0 @@
-import { EpochUTC, J2000, Kilometers, KilometersPerSecond, Vector3D } from '../../main';
-import { VerletBlendInterpolator } from '../VerletBlendInterpolator';
-
-describe('VerletBlendInterpolator', () => {
-  let sampleEphemeris: J2000[];
-  let interpolator: VerletBlendInterpolator;
-
-  beforeEach(() => {
-    const epoch1 = EpochUTC.fromDateTimeString('2024-01-01T00:00:00.000Z');
-    const epoch2 = EpochUTC.fromDateTimeString('2024-01-01T00:01:00.000Z');
-    const epoch3 = EpochUTC.fromDateTimeString('2024-01-01T00:02:00.000Z');
-
-    sampleEphemeris = [
-      new J2000(epoch1, new Vector3D(7000, 0, 0) as Vector3D<Kilometers>, new Vector3D(0, 7.5, 0) as Vector3D<KilometersPerSecond>),
-      new J2000(epoch2, new Vector3D(7100, 450, 0) as Vector3D<Kilometers>, new Vector3D(-0.5, 7.4, 0) as Vector3D<KilometersPerSecond>),
-      new J2000(epoch3, new Vector3D(7150, 890, 0) as Vector3D<Kilometers>, new Vector3D(-1.0, 7.3, 0) as Vector3D<KilometersPerSecond>),
-    ];
-
-    interpolator = new VerletBlendInterpolator(sampleEphemeris);
-  });
-
-  describe('constructor', () => {
-    it('should create an interpolator with ephemeris', () => {
-      expect(interpolator).toBeDefined();
-      expect(interpolator.ephemeris).toEqual(sampleEphemeris);
-    });
-  });
-
-  describe('sizeBytes', () => {
-    it('should calculate memory size correctly', () => {
-      const expectedSize = (64 * 7 * sampleEphemeris.length) / 8;
-
-      expect(interpolator.sizeBytes).toBe(expectedSize);
-    });
-  });
-
-  describe('window', () => {
-    it('should return the correct epoch window', () => {
-      const window = interpolator.window();
-
-      expect(window.start).toEqual(sampleEphemeris[0].epoch);
-      expect(window.end).toEqual(sampleEphemeris[2].epoch);
-    });
-  });
-
-  describe('getCachedState', () => {
-    it('should return exact match for cached epoch', () => {
-      const result = interpolator.getCachedState(sampleEphemeris[1].epoch);
-
-      expect(result).toEqual(sampleEphemeris[1]);
-    });
-
-    it('should return closest state for non-exact epoch', () => {
-      const testEpoch = EpochUTC.fromDateTimeString('2024-01-01T00:00:30.000Z');
-      const result = interpolator.getCachedState(testEpoch);
-
-      expect(result).toBeDefined();
-      expect([sampleEphemeris[0], sampleEphemeris[1]]).toContain(result);
-    });
-
-    it('should return null for epoch outside window', () => {
-      const beforeWindow = EpochUTC.fromDateTimeString('2023-12-31T23:59:00.000Z');
-      const afterWindow = EpochUTC.fromDateTimeString('2024-01-01T00:03:00.000Z');
-
-      expect(interpolator.getCachedState(beforeWindow)).toBeNull();
-      expect(interpolator.getCachedState(afterWindow)).toBeNull();
-    });
-  });
-
-  describe('interpolate', () => {
-    it('should return exact state for cached epoch', () => {
-      const result = interpolator.interpolate(sampleEphemeris[0].epoch);
-
-      expect(result).toBeDefined();
-      expect(result?.epoch).toEqual(sampleEphemeris[0].epoch);
-    });
-
-    it('should interpolate state for epoch within window', () => {
-      const testEpoch = EpochUTC.fromDateTimeString('2024-01-01T00:00:30.000Z');
-      const result = interpolator.interpolate(testEpoch);
-
-      expect(result).toBeDefined();
-      expect(result?.epoch.posix).toBeCloseTo(testEpoch.posix, 6);
-    });
-
-    it('should return null for epoch outside window', () => {
-      const beforeWindow = EpochUTC.fromDateTimeString('2023-12-31T23:59:00.000Z');
-      const afterWindow = EpochUTC.fromDateTimeString('2024-01-01T00:03:00.000Z');
-
-      expect(interpolator.interpolate(beforeWindow)).toBeNull();
-      expect(interpolator.interpolate(afterWindow)).toBeNull();
-    });
-  });
-
-  describe('toCubicSpline', () => {
-    it('should create a CubicSplineInterpolator', () => {
-      const spline = interpolator.toCubicSpline();
-
-      expect(spline).toBeDefined();
-    });
-  });
-
-  describe('toLagrange', () => {
-    it('should create a LagrangeInterpolator with default order', () => {
-      const lagrange = interpolator.toLagrange();
-
-      expect(lagrange).toBeDefined();
-    });
-
-    it('should create a LagrangeInterpolator with custom order', () => {
-      const lagrange = interpolator.toLagrange(5);
-
-      expect(lagrange).toBeDefined();
-    });
-  });
-});

package/src/interpolator/index.ts

@@ -1,32 +0,0 @@
-/**
- * @author 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 { Interpolator } from './Interpolator';
-export { FieldInterpolator } from './FieldInterpolator';
-export { StateInterpolator } from './StateInterpolator';
-export { CubicSpline } from './CubicSpline';
-export { CubicSplineInterpolator } from './CubicSplineInterpolator';
-export { LagrangeInterpolator } from './LagrangeInterpolator';
-export { ChebyshevCoefficients } from './ChebyshevCoefficients';
-export { ChebyshevInterpolator } from './ChebyshevInterpolator';
-export { VerletBlendInterpolator } from './VerletBlendInterpolator';

package/src/main.ts

@@ -1,87 +0,0 @@
-/**
- * @author 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 * from './errors/index';
-
-export * from './enums/index';
-
-export * from './types/types';
-
-export * from './interfaces/index';
-
-export * from './time/index';
-
-export * from './transforms/index';
-
-export * from './utils/index';
-
-export * from './operations/operations';
-
-export { BaseObject } from './objects/index';
-
-export { Earth } from './body/index';
-
-export * from './coordinate/index';
-
-export * from './observation/index';
-
-export * from './data/DataHandler';
-
-export * from './sgp4/index';
-
-export * from './objects/index';
-
-export * from './body/index';
-
-export * from './operations/index';
-
-export * from './force/index';
-
-export * from './propagator/index';
-
-export * from './orbit-determination/index';
-
-export * from './covariance/index';
-
-export * from './conjunction/index';
-
-export * from './sensor/index';
-
-export * from './parsers/index';
-
-export * from './comm/index';
-
-export * from './interpolator/index';
-
-export * from './maneuver/index';
-
-// Note: fetch module not exported - HorizonsAPI uses browser fetch which isn't available in Node.js
-// export * from './fetch/index';
-
-export * from './optimize/index';
-
-export * from './scheduling/index';
-
-export * from './constellation/index';
-
-export * from './orbit-design/index';

package/src/maneuver/TwoBurnOrbitTransfer.ts

@@ -1,79 +0,0 @@
-/**
- * @author Theodore Kruczek
- * @license AGPL-3.0-or-later
- * @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
- * 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 { Thrust } from '../force/Thrust';
-import { Earth, EpochUTC, MetersPerSecond, Seconds, SecondsPerMeterPerSecond } from '../main';
-
-// / Container for a two-burn orbit transfer.
-export class TwoBurnOrbitTransfer {
-  // / Create a new [TwoBurnOrbitTransfer] object.
-  constructor(
-    public vInit: number,
-    public vFinal: number,
-    public vTransA: number,
-    public vTransB: number,
-    public tTrans: Seconds,
-  ) {
-    // Nothing to do here.
-  }
-
-  /**
-   * Calculates the parameters for a Hohmann transfer orbit between two circular orbits.
-   * @param rInit The initial radius of the orbit. (km)
-   * @param rFinal The final radius of the orbit. (km)
-   * @returns An instance of TwoBurnOrbitTransfer containing the calculated parameters.
-   */
-  static hohmannTransfer(rInit: number, rFinal: number): TwoBurnOrbitTransfer {
-    const vInit = Math.sqrt(Earth.mu / rInit);
-    const vFinal = Math.sqrt(Earth.mu / rFinal);
-    const vTransA = vInit * (Math.sqrt((2.0 * rFinal) / (rInit + rFinal)) - 1.0);
-    const vTransB = vFinal * (1.0 - Math.sqrt((2 * rInit) / (rInit + rFinal)));
-    const tTrans = Math.PI * Math.sqrt((rInit + rFinal) ** 3 / (8.0 * Earth.mu)) as Seconds;
-
-    return new TwoBurnOrbitTransfer(vInit, vFinal, vTransA, vTransB, tTrans);
-  }
-
-  // / Return the total delta-velocity magnitude for both maneuvers _(km/s)_.
-  get deltaV(): number {
-    return Math.abs(this.vTransA) + Math.abs(this.vTransB);
-  }
-
-  /**
-   * Calculates the two maneuver thrusts required for a two-burn orbit transfer.
-   * @param epoch The epoch of the maneuver.
-   * @param durationRate The duration rate of the maneuver (s/m/s).
-   * @returns An array containing the two thrust objects representing the maneuvers.
-   */
-  toManeuvers(epoch: EpochUTC, durationRate = 0.0 as SecondsPerMeterPerSecond): [Thrust, Thrust] {
-    const mA = new Thrust(
-      epoch,
-      0.0 as MetersPerSecond,
-      this.vTransA * 1000.0 as MetersPerSecond,
-      0.0 as MetersPerSecond,
-      durationRate,
-    );
-    const mB = new Thrust(
-      epoch.roll(this.tTrans),
-      0.0 as MetersPerSecond,
-      this.vTransB * 1000.0 as MetersPerSecond,
-      0.0 as MetersPerSecond,
-      durationRate,
-    );
-
-    return [mA, mB];
-  }
-}