ootk 3.0.0 → 4.0.0-1

package/dist/interpolator/VerletBlendInterpolator.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 { copySign, Earth, EpochWindow, J2000, } from 'ootk-core';
+import { CubicSplineInterpolator } from './CubicSplineInterpolator';
+import { LagrangeInterpolator } from './LagrangeInterpolator';
+import { StateInterpolator } from './StateInterpolator';
+/**
+ * Two-body Velocity Verlet Blend interpolator.
+ *
+ * The [VerletBlendInterpolator] retains the original ephemerides, so the
+ * original _"truth"_ states can be retrieved if needed without imparting any
+ * additional error, so this can be used to build other interpolator types.
+ * The implementation is simple and very tolerant when working with sparse
+ * ephemerides.
+ */
+export class VerletBlendInterpolator extends StateInterpolator {
+    ephemeris;
+    constructor(ephemeris) {
+        super();
+        this.ephemeris = ephemeris;
+    }
+    get sizeBytes() {
+        return (64 * 7 * this.ephemeris.length) / 8;
+    }
+    window() {
+        return new EpochWindow(this.ephemeris[0].epoch, this.ephemeris[this.ephemeris.length - 1].epoch);
+    }
+    static getClosest_(target, s1, s2) {
+        return target - s1.epoch.posix >= s2.epoch.posix - target ? s2 : s1;
+    }
+    matchState_(epoch) {
+        const target = epoch.posix;
+        if (target <= this.ephemeris[0].epoch.posix) {
+            return this.ephemeris[0];
+        }
+        if (target >= this.ephemeris[this.ephemeris.length - 1].epoch.posix) {
+            return this.ephemeris[this.ephemeris.length - 1];
+        }
+        let i = 0;
+        let j = this.ephemeris.length;
+        let mid = 0;
+        while (i < j) {
+            mid = (i + j) >> 1;
+            if (this.ephemeris[mid].epoch.posix === target) {
+                return this.ephemeris[mid];
+            }
+            if (target < this.ephemeris[mid].epoch.posix) {
+                if (mid > 0 && target > this.ephemeris[mid - 1].epoch.posix) {
+                    return VerletBlendInterpolator.getClosest_(target, this.ephemeris[mid - 1], this.ephemeris[mid]);
+                }
+                j = mid;
+            }
+            else {
+                if (mid < this.ephemeris.length - 1 && target < this.ephemeris[mid + 1].epoch.posix) {
+                    return VerletBlendInterpolator.getClosest_(target, this.ephemeris[mid], this.ephemeris[mid + 1]);
+                }
+                i = mid + 1;
+            }
+        }
+        return this.ephemeris[mid];
+    }
+    static _gravity(position) {
+        const r = position.magnitude();
+        return position.scale(-Earth.mu / (r * r * r));
+    }
+    static integrate_(state, step) {
+        const x0 = state.position;
+        const a0 = VerletBlendInterpolator._gravity(x0);
+        const v0 = state.velocity;
+        const x1 = x0
+            .add(v0.scale(step))
+            .add(a0.scale(0.5 * step * step));
+        const a1 = VerletBlendInterpolator._gravity(x1);
+        const v1 = v0
+            .add(a0.add(a1).scale(0.5 * step));
+        return new J2000(state.epoch.roll(step), x1, v1);
+    }
+    interpolate(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        let state = this.matchState_(epoch);
+        while (state.epoch.posix !== epoch.posix) {
+            const delta = epoch.posix - state.epoch.posix;
+            const stepMag = Math.min(5.0, Math.abs(delta));
+            const stepSize = copySign(stepMag, delta);
+            state = VerletBlendInterpolator.integrate_(state, stepSize);
+        }
+        return state;
+    }
+    getCachedState(epoch) {
+        if (!this.inWindow(epoch)) {
+            return null;
+        }
+        return this.matchState_(epoch);
+    }
+    toCubicSpline() {
+        return CubicSplineInterpolator.fromEphemeris(this.ephemeris);
+    }
+    toLagrange(order = 10) {
+        return LagrangeInterpolator.fromEphemeris(this.ephemeris, order);
+    }
+}
+//# sourceMappingURL=VerletBlendInterpolator.js.map

package/dist/interpolator/VerletBlendInterpolator.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"VerletBlendInterpolator.js","sourceRoot":"","sources":["../../src/interpolator/VerletBlendInterpolator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EACL,QAAQ,EACR,KAAK,EAEL,WAAW,EACX,KAAK,GAKN,MAAM,WAAW,CAAC;AACnB,OAAO,EAAE,uBAAuB,EAAE,MAAM,2BAA2B,CAAC;AACpE,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC;AAC9D,OAAO,EAAE,iBAAiB,EAAE,MAAM,qBAAqB,CAAC;AAExD;;;;;;;;GAQG;AACH,MAAM,OAAO,uBAAwB,SAAQ,iBAAiB;IACzC;IAAnB,YAAmB,SAAkB;QACnC,KAAK,EAAE,CAAC;QADS,cAAS,GAAT,SAAS,CAAS;IAErC,CAAC;IAED,IAAI,SAAS;QACX,OAAO,CAAC,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;IAC9C,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,WAAW,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IACnG,CAAC;IAEO,MAAM,CAAC,WAAW,CAAC,MAAc,EAAE,EAAS,EAAE,EAAS;QAC7D,OAAO,MAAM,GAAG,EAAE,CAAC,KAAK,CAAC,KAAK,IAAI,EAAE,CAAC,KAAK,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;IACtE,CAAC;IAEO,WAAW,CAAC,KAAe;QACjC,MAAM,MAAM,GAAG,KAAK,CAAC,KAAK,CAAC;QAE3B,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;YAC3C,OAAO,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;SAC1B;QACD,IAAI,MAAM,IAAI,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;YACnE,OAAO,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;SAClD;QAED,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC;QAC9B,IAAI,GAAG,GAAG,CAAC,CAAC;QAEZ,OAAO,CAAC,GAAG,CAAC,EAAE;YACZ,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;YACnB,IAAI,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,KAAK,MAAM,EAAE;gBAC9C,OAAO,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;aAC5B;YACD,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;gBAC5C,IAAI,GAAG,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;oBAC3D,OAAO,uBAAuB,CAAC,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC,CAAC;iBAClG;gBACD,CAAC,GAAG,GAAG,CAAC;aACT;iBAAM;gBACL,IAAI,GAAG,GAAG,IAAI,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,IAAI,MAAM,GAAG,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,KAAK,CAAC,KAAK,EAAE;oBACnF,OAAO,uBAAuB,CAAC,WAAW,CAAC,MAAM,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,SAAS,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;iBAClG;gBACD,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;aACb;SACF;QAED,OAAO,IAAI,CAAC,SAAS,CAAC,GAAG,CAAC,CAAC;IAC7B,CAAC;IAEO,MAAM,CAAC,QAAQ,CAAC,QAAkB;QACxC,MAAM,CAAC,GAAG,QAAQ,CAAC,SAAS,EAAE,CAAC;QAE/B,OAAO,QAAQ,CAAC,KAAK,CAAC,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACjD,CAAC;IAEO,MAAM,CAAC,UAAU,CAAC,KAAY,EAAE,IAAa;QACnD,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,MAAM,EAAE,GAAG,uBAAuB,CAAC,QAAQ,CAAC,EAAE,CAAkC,CAAC;QACjF,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC;QAC1B,MAAM,EAAE,GAAG,EAAE;aACV,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,CAAoC,CAAC;aACtD,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,GAAG,IAAI,CAAoC,CAAC,CAAC;QACvE,MAAM,EAAE,GAAG,uBAAuB,CAAC,QAAQ,CAAC,EAAE,CAAkC,CAAC;QACjF,MAAM,EAAE,GAAG,EAAE;aACV,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,GAAG,GAAG,IAAI,CAAkC,CAAC,CAAC;QAEtE,OAAO,IAAI,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACnD,CAAC;IAED,WAAW,CAAC,KAAe;QACzB,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QACD,IAAI,KAAK,GAAG,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;QAEpC,OAAO,KAAK,CAAC,KAAK,CAAC,KAAK,KAAK,KAAK,CAAC,KAAK,EAAE;YACxC,MAAM,KAAK,GAAG,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,KAAK,CAAC,KAAK,CAAC;YAC9C,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;YAC/C,MAAM,QAAQ,GAAG,QAAQ,CAAC,OAAO,EAAE,KAAK,CAAY,CAAC;YAErD,KAAK,GAAG,uBAAuB,CAAC,UAAU,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;SAC7D;QAED,OAAO,KAAK,CAAC;IACf,CAAC;IAED,cAAc,CAAC,KAAe;QAC5B,IAAI,CAAC,IAAI,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE;YACzB,OAAO,IAAI,CAAC;SACb;QAED,OAAO,IAAI,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;IACjC,CAAC;IAED,aAAa;QACX,OAAO,uBAAuB,CAAC,aAAa,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;IAC/D,CAAC;IAED,UAAU,CAAC,KAAK,GAAG,EAAE;QACnB,OAAO,oBAAoB,CAAC,aAAa,CAAC,IAAI,CAAC,SAAS,EAAE,KAAK,CAAC,CAAC;IACnE,CAAC;CACF"}

package/{lib/ootk.js → dist/main.d.ts}

@@ -1,10 +1,12 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
+ * @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) 2020-2023 Theodore Kruczek
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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).
  *
  * 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,10 +19,11 @@
  * 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 { Sgp4 } from './sgp4/sgp4';
-export { Transforms } from './transforms/transforms';
-export { Utils } from './utils/utils';
-export { Tle } from './tle/tle';
-export * from './objects/index';
+export * from 'ootk-core';
+export * from './coordinate';
+export * from './objects';
+export * from './operations';
+export * from './transforms';
+export * from './interfaces';
 export * from './types/types';
-//# sourceMappingURL=ootk.js.map
+export * from './utils';

package/dist/{objects/space-object.d.ts → main.js}

@@ -1,13 +1,12 @@
 /**
- * @author Theodore Kruczek.
- * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
+ * @author @thkruz 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.
- *
  * @license AGPL-3.0-or-later
- * @Copyright (c) 2020-2023 Theodore Kruczek
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * 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).
  *
  * 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 +19,14 @@
  * 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;
-}
-export {};
+export * from 'ootk-core';
+export * from './coordinate';
+export * from './objects';
+export * from './operations';
+export * from './transforms';
+export * from './interfaces';
+// eslint-disable-next-line @typescript-eslint/ban-ts-comment
+// @ts-ignore
+export * from './types/types';
+export * from './utils';
+//# sourceMappingURL=main.js.map

package/dist/main.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"main.js","sourceRoot":"","sources":["../src/main.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;GAoBG;AAEH,cAAc,WAAW,CAAC;AAC1B,cAAc,cAAc,CAAC;AAC7B,cAAc,WAAW,CAAC;AAC1B,cAAc,cAAc,CAAC;AAC7B,cAAc,cAAc,CAAC;AAC7B,cAAc,cAAc,CAAC;AAC7B,6DAA6D;AAC7D,aAAa;AACb,cAAc,eAAe,CAAC;AAC9B,cAAc,SAAS,CAAC"}

package/dist/maneuver/TwoBurnOrbitTransfer.d.ts

@@ -0,0 +1,41 @@
+/**
+ * @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, Seconds, SecondsPerMeterPerSecond } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+export declare class TwoBurnOrbitTransfer {
+    vInit: number;
+    vFinal: number;
+    vTransA: number;
+    vTransB: number;
+    tTrans: Seconds;
+    constructor(vInit: number, vFinal: number, vTransA: number, vTransB: number, tTrans: Seconds);
+    /**
+     * 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;
+    get deltaV(): number;
+    /**
+     * 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?: SecondsPerMeterPerSecond): [Thrust, Thrust];
+}

package/dist/maneuver/TwoBurnOrbitTransfer.js

@@ -0,0 +1,65 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ * @license AGPL-3.0-or-later
+ * @copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Earth } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+// / Container for a two-burn orbit transfer.
+export class TwoBurnOrbitTransfer {
+    vInit;
+    vFinal;
+    vTransA;
+    vTransB;
+    tTrans;
+    // / Create a new [TwoBurnOrbitTransfer] object.
+    constructor(vInit, vFinal, vTransA, vTransB, tTrans) {
+        this.vInit = vInit;
+        this.vFinal = vFinal;
+        this.vTransA = vTransA;
+        this.vTransB = vTransB;
+        this.tTrans = tTrans;
+        // 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, rFinal) {
+        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));
+        return new TwoBurnOrbitTransfer(vInit, vFinal, vTransA, vTransB, tTrans);
+    }
+    // / Return the total delta-velocity magnitude for both maneuvers _(km/s)_.
+    get deltaV() {
+        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, durationRate = 0.0) {
+        const mA = new Thrust(epoch, 0.0, this.vTransA * 1000.0, 0.0, durationRate);
+        const mB = new Thrust(epoch.roll(this.tTrans), 0.0, this.vTransB * 1000.0, 0.0, durationRate);
+        return [mA, mB];
+    }
+}
+//# sourceMappingURL=TwoBurnOrbitTransfer.js.map

package/dist/maneuver/TwoBurnOrbitTransfer.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"TwoBurnOrbitTransfer.js","sourceRoot":"","sources":["../../src/maneuver/TwoBurnOrbitTransfer.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAAE,KAAK,EAAgE,MAAM,WAAW,CAAC;AAChG,OAAO,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAEzC,6CAA6C;AAC7C,MAAM,OAAO,oBAAoB;IAGtB;IACA;IACA;IACA;IACA;IANT,gDAAgD;IAChD,YACS,KAAa,EACb,MAAc,EACd,OAAe,EACf,OAAe,EACf,MAAe;QAJf,UAAK,GAAL,KAAK,CAAQ;QACb,WAAM,GAAN,MAAM,CAAQ;QACd,YAAO,GAAP,OAAO,CAAQ;QACf,YAAO,GAAP,OAAO,CAAQ;QACf,WAAM,GAAN,MAAM,CAAS;QAEtB,sBAAsB;IACxB,CAAC;IAED;;;;;OAKG;IACH,MAAM,CAAC,eAAe,CAAC,KAAa,EAAE,MAAc;QAClD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,KAAK,CAAC,CAAC;QAC1C,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,GAAG,MAAM,CAAC,CAAC;QAC5C,MAAM,OAAO,GAAG,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,MAAM,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAC7E,MAAM,OAAO,GAAG,MAAM,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAC3E,MAAM,MAAM,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,KAAK,CAAC,EAAE,CAAC,CAAY,CAAC;QAExF,OAAO,IAAI,oBAAoB,CAAC,KAAK,EAAE,MAAM,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,CAAC,CAAC;IAC3E,CAAC;IAED,2EAA2E;IAC3E,IAAI,MAAM;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IACzD,CAAC;IAED;;;;;OAKG;IACH,WAAW,CAAC,KAAe,EAAE,eAAe,GAA+B;QACzE,MAAM,EAAE,GAAG,IAAI,MAAM,CACnB,KAAK,EACL,GAAsB,EACtB,IAAI,CAAC,OAAO,GAAG,MAAyB,EACxC,GAAsB,EACtB,YAAY,CACb,CAAC;QACF,MAAM,EAAE,GAAG,IAAI,MAAM,CACnB,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,EACvB,GAAsB,EACtB,IAAI,CAAC,OAAO,GAAG,MAAyB,EACxC,GAAsB,EACtB,YAAY,CACb,CAAC;QAEF,OAAO,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;IAClB,CAAC;CACF"}

package/dist/maneuver/Waypoint.d.ts

@@ -0,0 +1,91 @@
+/**
+ * @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';
+import { Thrust } from '../force/Thrust';
+import { StateInterpolator } from '../interpolator/StateInterpolator';
+import { ForceModel } from './../force/ForceModel';
+export declare class Waypoint {
+    epoch: EpochUTC;
+    relativePosition: Vector3D<Kilometers>;
+    constructor(epoch: EpochUTC, relativePosition: Vector3D<Kilometers>);
+    /**
+     * Return the perturbed error in a [maneuver] when compared against the
+     * target [waypoint] given an initial [state], [forceModel],
+     * [target] interpolator, and speculative relative maneuver
+     * [components] _(m/s)_.
+     * @param waypoint The waypoint to target.
+     * @param maneuver The maneuver to perturb.
+     * @param state The initial state of the interceptor.
+     * @param forceModel The force model to use for propagation.
+     * @param target The target interpolator.
+     * @param components The speculative maneuver components.
+     * @returns The perturbed error in the maneuver.
+     */
+    static _error(waypoint: Waypoint, maneuver: Thrust, state: J2000, forceModel: ForceModel, target: StateInterpolator, components: Float64Array): number;
+    /**
+     * Generate a score function for refining perturbed [waypoint] maneuvers.
+     *
+     * The score function takes an array of speculative radial, intrack, and
+     * crosstrack components _(m/s)_ and returns the propagated error from the
+     * desired waypoint target.
+     * @param waypoint The waypoint to target.
+     * @param maneuver The maneuver to perturb.
+     * @param state The initial state of the interceptor.
+     * @param forceModel The force model to use for propagation.
+     * @param target The target interpolator.
+     * @returns A score function for refining maneuvers.
+     */
+    static _refineManeuverScore(waypoint: Waypoint, maneuver: Thrust, state: J2000, forceModel: ForceModel, target: StateInterpolator): (components: Float64Array) => number;
+    /**
+     * Convert an array of [waypoints] into a maneuver sequence given an
+     * [interceptor] state, [pivot] epoch for the first burn to arrive at the
+     * first waypoint, and [target] ephemeris interpolator.
+     *
+     * Optional arguments are as follows:
+     *  - `preManeuvers`: maneuvers to execute before the pivot burn
+     * - `postManeuvers`: maneuvers to execute after the last pivot burn
+     * - `durationRate`: thruster duration rate _(s/m/s)_
+     * - `forceModel`: interceptor force model, defaults to two-body
+     * - `refine`: refine maneuvers to account for perturbations if `true`
+     * - `maxIter`: maximum refinement iterations per maneuver
+     * - `printIter`: print debug information on each refinement iteration
+     * @param interceptor The interceptor state.
+     * @param pivot The epoch of the first burn.
+     * @param waypoints The waypoints to target.
+     * @param target The target interpolator.
+     * @param preManeuvers The maneuvers to execute before the pivot burn.
+     * @param postManeuvers The maneuvers to execute after the last pivot burn.
+     * @param root0 The optional arguments.
+     * @param root0.durationRate The thruster duration rate.
+     * @param root0.forceModel The interceptor force model.
+     * @param root0.refine Whether to refine maneuvers to account for perturbations.
+     * @param root0.maxIter The maximum refinement iterations per maneuver.
+     * @param root0.printIter Whether to print debug information on each refinement iteration.
+     * @returns An array of maneuvers.
+     */
+    static toManeuvers(interceptor: J2000, pivot: EpochUTC, waypoints: Waypoint[], target: StateInterpolator, preManeuvers: Thrust[] | null, postManeuvers: Thrust[] | null, { durationRate, forceModel, refine, maxIter, printIter, }?: {
+        durationRate?: number;
+        forceModel?: ForceModel;
+        refine?: boolean;
+        maxIter?: number;
+        printIter?: boolean;
+    }): Thrust[];
+    static _refineManeuvers(waypoints: Waypoint[], maneuvers: Thrust[], interceptor: J2000, forceModel: ForceModel, target: StateInterpolator, { maxIter, printIter, }?: {
+        maxIter?: number;
+        printIter?: boolean;
+    }): Thrust[];
+}

package/dist/maneuver/Waypoint.js

@@ -0,0 +1,179 @@
+/**
+ * @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 { RIC, Vector3D, } from 'ootk-core';
+import { Thrust } from '../force/Thrust';
+import { ForceModel } from './../force/ForceModel';
+import { DownhillSimplex } from './../optimize/DownhillSimplex';
+import { LambertIOD } from './../orbit_determination/LambertIOD';
+import { RungeKutta89Propagator } from './../propagator/RungeKutta89Propagator';
+// / Relative waypoint targeting.
+export class Waypoint {
+    epoch;
+    relativePosition;
+    // / Create a new [Waypoint] object.
+    constructor(epoch, relativePosition) {
+        this.epoch = epoch;
+        this.relativePosition = relativePosition;
+    }
+    /**
+     * Return the perturbed error in a [maneuver] when compared against the
+     * target [waypoint] given an initial [state], [forceModel],
+     * [target] interpolator, and speculative relative maneuver
+     * [components] _(m/s)_.
+     * @param waypoint The waypoint to target.
+     * @param maneuver The maneuver to perturb.
+     * @param state The initial state of the interceptor.
+     * @param forceModel The force model to use for propagation.
+     * @param target The target interpolator.
+     * @param components The speculative maneuver components.
+     * @returns The perturbed error in the maneuver.
+     */
+    static _error(waypoint, maneuver, state, forceModel, target, components) {
+        const testManeuver = new Thrust(maneuver.center, components[0], components[1], components[2], maneuver.durationRate);
+        const propagator = new RungeKutta89Propagator(state, forceModel);
+        const maneuverSteps = propagator.maneuver(testManeuver);
+        const postManeuver = maneuverSteps[maneuverSteps.length - 1];
+        const interceptor = new RungeKutta89Propagator(postManeuver, forceModel).propagate(waypoint.epoch);
+        const targetState = target.interpolate(waypoint.epoch);
+        if (targetState === null) {
+            throw new Error('Error calculation failed; epoch is outside the target interpolator ephemeris window.');
+        }
+        const expected = waypoint.relativePosition;
+        const actual = RIC.fromJ2000(interceptor, targetState);
+        return actual.position.distance(expected);
+    }
+    /**
+     * Generate a score function for refining perturbed [waypoint] maneuvers.
+     *
+     * The score function takes an array of speculative radial, intrack, and
+     * crosstrack components _(m/s)_ and returns the propagated error from the
+     * desired waypoint target.
+     * @param waypoint The waypoint to target.
+     * @param maneuver The maneuver to perturb.
+     * @param state The initial state of the interceptor.
+     * @param forceModel The force model to use for propagation.
+     * @param target The target interpolator.
+     * @returns A score function for refining maneuvers.
+     */
+    static _refineManeuverScore(waypoint, maneuver, state, forceModel, target) {
+        return (components) => Waypoint._error(waypoint, maneuver, state, forceModel, target, components);
+    }
+    /**
+     * Convert an array of [waypoints] into a maneuver sequence given an
+     * [interceptor] state, [pivot] epoch for the first burn to arrive at the
+     * first waypoint, and [target] ephemeris interpolator.
+     *
+     * Optional arguments are as follows:
+     *  - `preManeuvers`: maneuvers to execute before the pivot burn
+     * - `postManeuvers`: maneuvers to execute after the last pivot burn
+     * - `durationRate`: thruster duration rate _(s/m/s)_
+     * - `forceModel`: interceptor force model, defaults to two-body
+     * - `refine`: refine maneuvers to account for perturbations if `true`
+     * - `maxIter`: maximum refinement iterations per maneuver
+     * - `printIter`: print debug information on each refinement iteration
+     * @param interceptor The interceptor state.
+     * @param pivot The epoch of the first burn.
+     * @param waypoints The waypoints to target.
+     * @param target The target interpolator.
+     * @param preManeuvers The maneuvers to execute before the pivot burn.
+     * @param postManeuvers The maneuvers to execute after the last pivot burn.
+     * @param root0 The optional arguments.
+     * @param root0.durationRate The thruster duration rate.
+     * @param root0.forceModel The interceptor force model.
+     * @param root0.refine Whether to refine maneuvers to account for perturbations.
+     * @param root0.maxIter The maximum refinement iterations per maneuver.
+     * @param root0.printIter Whether to print debug information on each refinement iteration.
+     * @returns An array of maneuvers.
+     */
+    static toManeuvers(interceptor, pivot, waypoints, target, preManeuvers, postManeuvers, { durationRate = 0.0, forceModel, refine = false, maxIter = 500, printIter = false, } = {}) {
+        const preMnv = preManeuvers ?? [];
+        const postMnv = postManeuvers ?? [];
+        let state = interceptor;
+        if (preMnv.length > 0) {
+            for (const maneuver of preMnv) {
+                const mvrStep = new RungeKutta89Propagator(state, forceModel).maneuver(maneuver);
+                state = mvrStep[mvrStep.length - 1];
+            }
+        }
+        state = new RungeKutta89Propagator(state, forceModel).propagate(pivot);
+        const pivotState = state;
+        let waypointManeuvers = [];
+        for (const wp of waypoints) {
+            const targetWp = new RIC(wp.relativePosition, Vector3D.origin);
+            const targetState = target.interpolate(wp.epoch);
+            if (targetState === null) {
+                throw new Error('Waypoint outside target interpolator window.');
+            }
+            const wpState = targetWp.toJ2000(targetState);
+            const tof = wp.epoch.difference(state.epoch);
+            const revs = Math.floor(tof / state.period);
+            const shortPath = LambertIOD.useShortPath(state, targetState);
+            const lambert = new LambertIOD();
+            const components = lambert.estimate(state.position, wpState.position, state.epoch, wp.epoch, {
+                posigrade: shortPath,
+                nRev: revs,
+            });
+            if (components === null) {
+                throw new Error('Lambert solve result is null.');
+            }
+            const componentsRel = RIC.fromJ2000(components, state).velocity.scale(1e3);
+            const maneuver = new Thrust(state.epoch, componentsRel.x, componentsRel.y, componentsRel.z, durationRate);
+            const tempProp = new RungeKutta89Propagator(state, forceModel);
+            tempProp.maneuver(maneuver);
+            state = tempProp.propagate(wp.epoch);
+            waypointManeuvers.push(maneuver);
+        }
+        if (refine) {
+            forceModel ??= new ForceModel().setGravity();
+            waypointManeuvers = this._refineManeuvers(waypoints, waypointManeuvers, pivotState, forceModel, target, {
+                maxIter,
+                printIter,
+            });
+        }
+        const output = [];
+        output.push(...preMnv);
+        output.push(...waypointManeuvers);
+        output.push(...postMnv);
+        return output;
+    }
+    static _refineManeuvers(waypoints, maneuvers, interceptor, forceModel, target, { maxIter = 500, printIter = false, } = {}) {
+        let state = interceptor;
+        const output = [];
+        for (let i = 0; i < maneuvers.length; i++) {
+            const maneuver = maneuvers[i];
+            const waypoint = waypoints[i];
+            const guess = maneuver.deltaV.scale(1e3).toArray();
+            const simplex = DownhillSimplex.generateSimplex(guess, 1e-1);
+            const scoreFn = this._refineManeuverScore(waypoint, maneuver, state, forceModel, target);
+            const results = DownhillSimplex.solveSimplex(scoreFn, simplex, {
+                maxIter,
+                xTolerance: 1e-6,
+                fTolerance: 1e-6,
+                printIter,
+            });
+            const tR = results[0];
+            const tI = results[1];
+            const tC = results[2];
+            const newManeuver = new Thrust(maneuver.center, tR, tI, tC, maneuver.durationRate);
+            output.push(newManeuver);
+            const mvrStep = new RungeKutta89Propagator(state, forceModel).maneuver(newManeuver);
+            state = mvrStep[mvrStep.length - 1];
+        }
+        return output;
+    }
+}
+//# sourceMappingURL=Waypoint.js.map

package/dist/maneuver/Waypoint.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Waypoint.js","sourceRoot":"","sources":["../../src/maneuver/Waypoint.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;GAeG;AAEH,OAAO,EAML,GAAG,EAEH,QAAQ,GACT,MAAM,WAAW,CAAC;AACnB,OAAO,EAAE,MAAM,EAAE,MAAM,iBAAiB,CAAC;AAEzC,OAAO,EAAE,UAAU,EAAE,MAAM,uBAAuB,CAAC;AACnD,OAAO,EAAE,eAAe,EAAE,MAAM,+BAA+B,CAAC;AAChE,OAAO,EAAE,UAAU,EAAE,MAAM,qCAAqC,CAAC;AACjE,OAAO,EAAE,sBAAsB,EAAE,MAAM,wCAAwC,CAAC;AAEhF,iCAAiC;AACjC,MAAM,OAAO,QAAQ;IACnB,KAAK,CAAW;IAChB,gBAAgB,CAAuB;IAEvC,oCAAoC;IACpC,YAAY,KAAe,EAAE,gBAAsC;QACjE,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,gBAAgB,GAAG,gBAAgB,CAAC;IAC3C,CAAC;IAED;;;;;;;;;;;;OAYG;IACH,MAAM,CAAC,MAAM,CACX,QAAkB,EAClB,QAAgB,EAChB,KAAY,EACZ,UAAsB,EACtB,MAAyB,EACzB,UAAwB;QAExB,MAAM,YAAY,GAAG,IAAI,MAAM,CAC7B,QAAQ,CAAC,MAAM,EACf,UAAU,CAAC,CAAC,CAAoB,EAChC,UAAU,CAAC,CAAC,CAAoB,EAChC,UAAU,CAAC,CAAC,CAAoB,EAChC,QAAQ,CAAC,YAAY,CACtB,CAAC;QACF,MAAM,UAAU,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;QACjE,MAAM,aAAa,GAAG,UAAU,CAAC,QAAQ,CAAC,YAAY,CAAC,CAAC;QACxD,MAAM,YAAY,GAAG,aAAa,CAAC,aAAa,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;QAC7D,MAAM,WAAW,GAAG,IAAI,sBAAsB,CAAC,YAAY,EAAE,UAAU,CAAC,CAAC,SAAS,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QACnG,MAAM,WAAW,GAAG,MAAM,CAAC,WAAW,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC;QAEvD,IAAI,WAAW,KAAK,IAAI,EAAE;YACxB,MAAM,IAAI,KAAK,CAAC,sFAAsF,CAAC,CAAC;SACzG;QACD,MAAM,QAAQ,GAAG,QAAQ,CAAC,gBAAgB,CAAC;QAC3C,MAAM,MAAM,GAAG,GAAG,CAAC,SAAS,CAAC,WAAW,EAAE,WAAW,CAAC,CAAC;QAEvD,OAAO,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;IAC5C,CAAC;IAED;;;;;;;;;;;;OAYG;IACH,MAAM,CAAC,oBAAoB,CACzB,QAAkB,EAClB,QAAgB,EAChB,KAAY,EACZ,UAAsB,EACtB,MAAyB;QAEzB,OAAO,CAAC,UAAwB,EAAE,EAAE,CAAC,QAAQ,CAAC,MAAM,CAAC,QAAQ,EAAE,QAAQ,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,EAAE,UAAU,CAAC,CAAC;IAClH,CAAC;IAED;;;;;;;;;;;;;;;;;;;;;;;;;;OA0BG;IACH,MAAM,CAAC,WAAW,CAChB,WAAkB,EAClB,KAAe,EACf,SAAqB,EACrB,MAAyB,EACzB,YAA6B,EAC7B,aAA8B,EAC9B,EACE,YAAY,GAAG,GAAG,EAClB,UAAU,EACV,MAAM,GAAG,KAAK,EACd,OAAO,GAAG,GAAG,EACb,SAAS,GAAG,KAAK,MAOf,EAAE;QAEN,MAAM,MAAM,GAAG,YAAY,IAAI,EAAE,CAAC;QAClC,MAAM,OAAO,GAAG,aAAa,IAAI,EAAE,CAAC;QACpC,IAAI,KAAK,GAAG,WAAW,CAAC;QAExB,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE;YACrB,KAAK,MAAM,QAAQ,IAAI,MAAM,EAAE;gBAC7B,MAAM,OAAO,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;gBAEjF,KAAK,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;aACrC;SACF;QACD,KAAK,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC;QACvE,MAAM,UAAU,GAAG,KAAK,CAAC;QACzB,IAAI,iBAAiB,GAAa,EAAE,CAAC;QAErC,KAAK,MAAM,EAAE,IAAI,SAAS,EAAE;YAC1B,MAAM,QAAQ,GAAG,IAAI,GAAG,CAAC,EAAE,CAAC,gBAAgB,EAAE,QAAQ,CAAC,MAAuC,CAAC,CAAC;YAChG,MAAM,WAAW,GAAG,MAAM,CAAC,WAAW,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;YAEjD,IAAI,WAAW,KAAK,IAAI,EAAE;gBACxB,MAAM,IAAI,KAAK,CAAC,8CAA8C,CAAC,CAAC;aACjE;YACD,MAAM,OAAO,GAAG,QAAQ,CAAC,OAAO,CAAC,WAAW,CAAC,CAAC;YAC9C,MAAM,GAAG,GAAG,EAAE,CAAC,KAAK,CAAC,UAAU,CAAC,KAAK,CAAC,KAAK,CAAC,CAAC;YAC7C,MAAM,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,KAAK,CAAC,MAAM,CAAC,CAAC;YAC5C,MAAM,SAAS,GAAG,UAAU,CAAC,YAAY,CAAC,KAAK,EAAE,WAAW,CAAC,CAAC;YAC9D,MAAM,OAAO,GAAG,IAAI,UAAU,EAAE,CAAC;YACjC,MAAM,UAAU,GAAG,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAC,QAAQ,EAAE,OAAO,CAAC,QAAQ,EAAE,KAAK,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,EAAE;gBAC3F,SAAS,EAAE,SAAS;gBACpB,IAAI,EAAE,IAAI;aACX,CAAC,CAAC;YAEH,IAAI,UAAU,KAAK,IAAI,EAAE;gBACvB,MAAM,IAAI,KAAK,CAAC,+BAA+B,CAAC,CAAC;aAClD;YACD,MAAM,aAAa,GAAG,GAAG,CAAC,SAAS,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,GAAG,CAA8B,CAAC;YACxG,MAAM,QAAQ,GAAG,IAAI,MAAM,CACzB,KAAK,CAAC,KAAK,EACX,aAAa,CAAC,CAAC,EACf,aAAa,CAAC,CAAC,EACf,aAAa,CAAC,CAAC,EACf,YAAwC,CACzC,CAAC;YACF,MAAM,QAAQ,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC;YAE/D,QAAQ,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;YAC5B,KAAK,GAAG,QAAQ,CAAC,SAAS,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC;YACrC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;SAClC;QACD,IAAI,MAAM,EAAE;YACV,UAAU,KAAK,IAAI,UAAU,EAAE,CAAC,UAAU,EAAE,CAAC;YAC7C,iBAAiB,GAAG,IAAI,CAAC,gBAAgB,CAAC,SAAS,EAAE,iBAAiB,EAAE,UAAU,EAAE,UAAU,EAAE,MAAM,EAAE;gBACtG,OAAO;gBACP,SAAS;aACV,CAAC,CAAC;SACJ;QACD,MAAM,MAAM,GAAa,EAAE,CAAC;QAE5B,MAAM,CAAC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAC;QACvB,MAAM,CAAC,IAAI,CAAC,GAAG,iBAAiB,CAAC,CAAC;QAClC,MAAM,CAAC,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC;QAExB,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,MAAM,CAAC,gBAAgB,CACrB,SAAqB,EACrB,SAAmB,EACnB,WAAkB,EAClB,UAAsB,EACtB,MAAyB,EACzB,EACE,OAAO,GAAG,GAAG,EACb,SAAS,GAAG,KAAK,MAIf,EAAE;QAEN,IAAI,KAAK,GAAG,WAAW,CAAC;QACxB,MAAM,MAAM,GAAa,EAAE,CAAC;QAE5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;YACzC,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,QAAQ,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,KAAK,GAAG,QAAQ,CAAC,MAAM,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,CAAC;YACnD,MAAM,OAAO,GAAG,eAAe,CAAC,eAAe,CAAC,KAAK,EAAE,IAAI,CAAC,CAAC;YAC7D,MAAM,OAAO,GAAG,IAAI,CAAC,oBAAoB,CAAC,QAAQ,EAAE,QAAQ,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,CAAC,CAAC;YACzF,MAAM,OAAO,GAAG,eAAe,CAAC,YAAY,CAAC,OAAO,EAAE,OAAO,EAAE;gBAC7D,OAAO;gBACP,UAAU,EAAE,IAAI;gBAChB,UAAU,EAAE,IAAI;gBAChB,SAAS;aACV,CAAC,CAAC;YACH,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAoB,CAAC;YACzC,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAoB,CAAC;YACzC,MAAM,EAAE,GAAG,OAAO,CAAC,CAAC,CAAoB,CAAC;YACzC,MAAM,WAAW,GAAG,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,QAAQ,CAAC,YAAY,CAAC,CAAC;YAEnF,MAAM,CAAC,IAAI,CAAC,WAAW,CAAC,CAAC;YACzB,MAAM,OAAO,GAAG,IAAI,sBAAsB,CAAC,KAAK,EAAE,UAAU,CAAC,CAAC,QAAQ,CAAC,WAAW,CAAC,CAAC;YAEpF,KAAK,GAAG,OAAO,CAAC,OAAO,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;SACrC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;CACF"}