ootk 3.0.0 → 4.0.0-0

package/lib/operations/BoxMuller.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 { Random, Vector } from 'ootk-core';
+export declare class BoxMuller {
+    private _index;
+    private _cache;
+    mu: number;
+    sigma: number;
+    rand: Random;
+    /**
+     * Create a new [BoxMuller] object with mean [mu], standard deviation
+     * [sigma], and [seed] number.
+     */
+    constructor(mu: number, sigma: number, seed?: number);
+    _generate(): void;
+    /**
+     * Generate a gaussian number, with mean [mu] and standard
+     * deviation [sigma].
+     */
+    nextGauss(): number;
+    /**
+     * Generate a [Vector] of gaussian numbers, with mean [mu] and standard
+     * deviation [sigma].
+     */
+    gaussVector(n: number): Vector;
+}

package/lib/operations/BoxMuller.js

@@ -0,0 +1,66 @@
+/**
+ * @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 { Random, TAU, Vector } from 'ootk-core';
+// / Box-Muller random Gaussian number generator.
+export class BoxMuller {
+    /**
+     * Create a new [BoxMuller] object with mean [mu], standard deviation
+     * [sigma], and [seed] number.
+     */
+    constructor(mu, sigma, seed = 0) {
+        this._index = 0;
+        this._cache = new Float64Array(2);
+        this.mu = mu;
+        this.sigma = sigma;
+        this.rand = new Random(seed);
+        this._generate();
+    }
+    // / Refill the cache with random Gaussian numbers.
+    _generate() {
+        this._index = 0;
+        const u1 = this.rand.nextFloat();
+        const u2 = this.rand.nextFloat();
+        const mag = this.sigma * Math.sqrt(-2.0 * Math.log(u1));
+        this._cache[0] = mag * Math.cos(TAU * u2) + this.mu;
+        this._cache[1] = mag * Math.sin(TAU * u2) + this.mu;
+    }
+    /**
+     * Generate a gaussian number, with mean [mu] and standard
+     * deviation [sigma].
+     */
+    nextGauss() {
+        if (this._index > 1) {
+            this._generate();
+        }
+        const result = this._cache[this._index];
+        this._index++;
+        return result;
+    }
+    /**
+     * Generate a [Vector] of gaussian numbers, with mean [mu] and standard
+     * deviation [sigma].
+     */
+    gaussVector(n) {
+        const result = new Float64Array(n);
+        for (let i = 0; i < n; i++) {
+            result[i] = this.nextGauss();
+        }
+        return new Vector(result);
+    }
+}
+//# sourceMappingURL=BoxMuller.js.map

package/lib/operations/BoxMuller.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"BoxMuller.js","sourceRoot":"","sources":["../../src/operations/BoxMuller.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,WAAW,CAAC;AAEhD,iDAAiD;AACjD,MAAM,OAAO,SAAS;IAUpB;;;OAGG;IACH,YAAY,EAAU,EAAE,KAAa,EAAE,IAAI,GAAG,CAAC;QAbvC,WAAM,GAAG,CAAC,CAAC;QACX,WAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAanC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC;QACb,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,IAAI,GAAG,IAAI,MAAM,CAAC,IAAI,CAAC,CAAC;QAC7B,IAAI,CAAC,SAAS,EAAE,CAAC;IACnB,CAAC;IAED,mDAAmD;IACnD,SAAS;QACP,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC;QAChB,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC;QACjC,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;QAExD,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACpD,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACtD,CAAC;IAED;;;OAGG;IACH,SAAS;QACP,IAAI,IAAI,CAAC,MAAM,GAAG,CAAC,EAAE;YACnB,IAAI,CAAC,SAAS,EAAE,CAAC;SAClB;QACD,MAAM,MAAM,GAAG,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;QAExC,IAAI,CAAC,MAAM,EAAE,CAAC;QAEd,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;OAGG;IACH,WAAW,CAAC,CAAS;QACnB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;QAEnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;SAC9B;QAED,OAAO,IAAI,MAAM,CAAC,MAAM,CAAC,CAAC;IAC5B,CAAC;CACF"}

package/lib/operations/EulerAngles.d.ts

@@ -0,0 +1,50 @@
+/**
+ * @author @thkruz Theodore Kruczek
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2024 Theodore Kruczek
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
+ * without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
+ * See the GNU Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { Matrix, Vector3D } from 'ootk-core';
+export declare class EulerAngles {
+    roll: number;
+    pitch: number;
+    yaw: number;
+    /**
+     * Create a new [EulerAngles] object from [roll], [pitch], and [yaw]
+     * angles _(rad)_.
+     */
+    constructor(roll: number, pitch: number, yaw: number);
+    /**
+     * Create a new [EulerAngles] object from [roll], [pitch], and [yaw]
+     * angles _(deg)_.
+     */
+    static fromDegrees(rollDeg: number, pitchDeg: number, yawDeg: number): EulerAngles;
+    /**
+     * Create a new [EulerAngles] object from 3-2-1 ordered direction cosine
+     * matrix [c].
+     */
+    static fromDcm321(c: Matrix): EulerAngles;
+    get rollDegrees(): number;
+    get pitchDegrees(): number;
+    get yawDegrees(): number;
+    get phi(): number;
+    get theta(): number;
+    get psi(): number;
+    get phiDegrees(): number;
+    get thetaDegrees(): number;
+    get psiDegrees(): number;
+    toString(precision?: number): string;
+    dcm321(): Matrix;
+    rotateVector321(v: Vector3D): Vector3D;
+}

package/lib/operations/EulerAngles.js

@@ -0,0 +1,111 @@
+/**
+ * @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, Matrix, RAD2DEG } from 'ootk-core';
+// / Class containing Euler angles.
+export class EulerAngles {
+    /**
+     * Create a new [EulerAngles] object from [roll], [pitch], and [yaw]
+     * angles _(rad)_.
+     */
+    constructor(roll, pitch, yaw) {
+        this.roll = roll;
+        this.pitch = pitch;
+        this.yaw = yaw;
+    }
+    /**
+     * Create a new [EulerAngles] object from [roll], [pitch], and [yaw]
+     * angles _(deg)_.
+     */
+    static fromDegrees(rollDeg, pitchDeg, yawDeg) {
+        const roll = rollDeg * DEG2RAD;
+        const pitch = pitchDeg * DEG2RAD;
+        const yaw = yawDeg * DEG2RAD;
+        return new EulerAngles(roll, pitch, yaw);
+    }
+    /**
+     * Create a new [EulerAngles] object from 3-2-1 ordered direction cosine
+     * matrix [c].
+     */
+    static fromDcm321(c) {
+        const roll = Math.atan(c[1][2] / c[2][2]);
+        const pitch = -Math.asin(c[0][2]);
+        const yaw = Math.atan(c[0][1] / c[0][0]);
+        return new EulerAngles(roll, pitch, yaw);
+    }
+    // / Roll component _(deg)_.
+    get rollDegrees() {
+        return this.roll * RAD2DEG;
+    }
+    // / Pitch component _(deg)_.
+    get pitchDegrees() {
+        return this.pitch * RAD2DEG;
+    }
+    // / Yaw component _(deg)_.
+    get yawDegrees() {
+        return this.yaw * RAD2DEG;
+    }
+    // / Roll component alias _(rad)_.
+    get phi() {
+        return this.roll;
+    }
+    // / Pitch component alias _(rad)_.
+    get theta() {
+        return this.pitch;
+    }
+    // / Yaw component alias _(rad)_.
+    get psi() {
+        return this.yaw;
+    }
+    // / Roll component alias _(deg)_.
+    get phiDegrees() {
+        return this.phi * RAD2DEG;
+    }
+    // / Pitch component alias _(deg)_.
+    get thetaDegrees() {
+        return this.theta * RAD2DEG;
+    }
+    // / Yaw component alias _(deg)_.
+    get psiDegrees() {
+        return this.psi * RAD2DEG;
+    }
+    toString(precision = 6) {
+        const rollStr = this.rollDegrees.toFixed(precision);
+        const pitchStr = this.pitchDegrees.toFixed(precision);
+        const yawStr = this.yawDegrees.toFixed(precision);
+        return `Euler(roll: ${rollStr}°, pitch: ${pitchStr}°, yaw: ${yawStr}°)`;
+    }
+    // / Convert this to a 3-2-1 ordered direction cosine [Matrix].
+    dcm321() {
+        const sPhi = Math.sin(this.phi);
+        const cPhi = Math.cos(this.phi);
+        const sTheta = Math.sin(this.theta);
+        const cTheta = Math.cos(this.theta);
+        const sPsi = Math.sin(this.psi);
+        const cPsi = Math.cos(this.psi);
+        return new Matrix([
+            [cTheta * cPsi, cTheta * sPsi, -sTheta],
+            [sPhi * sTheta * cPsi - cPhi * sPsi, sPhi * sTheta * sPsi + cPhi * cPsi, sPhi * cTheta],
+            [cPhi * sTheta * cPsi + sPhi * sPsi, cPhi * sTheta * sPsi - sPhi * cPsi, cPhi * cTheta],
+        ]);
+    }
+    // / Perform a 3-2-1 ordered rotation on provided vector [v].
+    rotateVector321(v) {
+        return this.dcm321().multiplyVector3D(v);
+    }
+}
+//# sourceMappingURL=EulerAngles.js.map

package/lib/operations/EulerAngles.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"EulerAngles.js","sourceRoot":"","sources":["../../src/operations/EulerAngles.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,OAAO,EAAE,MAAM,EAAE,OAAO,EAAY,MAAM,WAAW,CAAC;AAE/D,mCAAmC;AACnC,MAAM,OAAO,WAAW;IAUtB;;;OAGG;IACH,YAAY,IAAY,EAAE,KAAa,EAAE,GAAW;QAClD,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC;QACjB,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC;QACnB,IAAI,CAAC,GAAG,GAAG,GAAG,CAAC;IACjB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,WAAW,CAAC,OAAe,EAAE,QAAgB,EAAE,MAAc;QAClE,MAAM,IAAI,GAAG,OAAO,GAAG,OAAO,CAAC;QAC/B,MAAM,KAAK,GAAG,QAAQ,GAAG,OAAO,CAAC;QACjC,MAAM,GAAG,GAAG,MAAM,GAAG,OAAO,CAAC;QAE7B,OAAO,IAAI,WAAW,CAAC,IAAI,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,UAAU,CAAC,CAAS;QACzB,MAAM,IAAI,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC1C,MAAM,KAAK,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClC,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEzC,OAAO,IAAI,WAAW,CAAC,IAAI,EAAE,KAAK,EAAE,GAAG,CAAC,CAAC;IAC3C,CAAC;IAED,4BAA4B;IAC5B,IAAI,WAAW;QACb,OAAO,IAAI,CAAC,IAAI,GAAG,OAAO,CAAC;IAC7B,CAAC;IAED,6BAA6B;IAC7B,IAAI,YAAY;QACd,OAAO,IAAI,CAAC,KAAK,GAAG,OAAO,CAAC;IAC9B,CAAC;IAED,2BAA2B;IAC3B,IAAI,UAAU;QACZ,OAAO,IAAI,CAAC,GAAG,GAAG,OAAO,CAAC;IAC5B,CAAC;IAED,kCAAkC;IAClC,IAAI,GAAG;QACL,OAAO,IAAI,CAAC,IAAI,CAAC;IACnB,CAAC;IAED,mCAAmC;IACnC,IAAI,KAAK;QACP,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;IAED,iCAAiC;IACjC,IAAI,GAAG;QACL,OAAO,IAAI,CAAC,GAAG,CAAC;IAClB,CAAC;IAED,kCAAkC;IAClC,IAAI,UAAU;QACZ,OAAO,IAAI,CAAC,GAAG,GAAG,OAAO,CAAC;IAC5B,CAAC;IAED,mCAAmC;IACnC,IAAI,YAAY;QACd,OAAO,IAAI,CAAC,KAAK,GAAG,OAAO,CAAC;IAC9B,CAAC;IAED,iCAAiC;IACjC,IAAI,UAAU;QACZ,OAAO,IAAI,CAAC,GAAG,GAAG,OAAO,CAAC;IAC5B,CAAC;IAED,QAAQ,CAAC,SAAS,GAAG,CAAC;QACpB,MAAM,OAAO,GAAG,IAAI,CAAC,WAAW,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACpD,MAAM,QAAQ,GAAG,IAAI,CAAC,YAAY,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACtD,MAAM,MAAM,GAAG,IAAI,CAAC,UAAU,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QAElD,OAAO,eAAe,OAAO,aAAa,QAAQ,WAAW,MAAM,IAAI,CAAC;IAC1E,CAAC;IAED,+DAA+D;IAC/D,MAAM;QACJ,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAChC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAChC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACpC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACpC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAChC,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QAEhC,OAAO,IAAI,MAAM,CAAC;YAChB,CAAC,MAAM,GAAG,IAAI,EAAE,MAAM,GAAG,IAAI,EAAE,CAAC,MAAM,CAAC;YACvC,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,EAAE,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,EAAE,IAAI,GAAG,MAAM,CAAC;YACvF,CAAC,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,EAAE,IAAI,GAAG,MAAM,GAAG,IAAI,GAAG,IAAI,GAAG,IAAI,EAAE,IAAI,GAAG,MAAM,CAAC;SACxF,CAAC,CAAC;IACL,CAAC;IAED,6DAA6D;IAC7D,eAAe,CAAC,CAAW;QACzB,OAAO,IAAI,CAAC,MAAM,EAAE,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;IAC3C,CAAC;CACF"}

package/lib/operations/Quaternion.d.ts

@@ -0,0 +1,57 @@
+/**
+ * @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 { Matrix, Vector, Vector3D } from 'ootk-core';
+export declare class Quaternion {
+    x: number;
+    y: number;
+    z: number;
+    w: number;
+    constructor(x: number, y: number, z: number, w: number);
+    static zero: Quaternion;
+    static one: Quaternion;
+    static xAxis: Quaternion;
+    static yAxis: Quaternion;
+    static zAxis: Quaternion;
+    toString(precision?: number): string;
+    positivePolar(): Quaternion;
+    magnitudeSquared(): number;
+    magnitude(): number;
+    scale(n: number): Quaternion;
+    negate(): Quaternion;
+    normalize(): Quaternion;
+    conjugate(): Quaternion;
+    inverse(): Quaternion;
+    add(q: Quaternion): Quaternion;
+    subtract(q: Quaternion): Quaternion;
+    addReal(n: number): Quaternion;
+    multiply(q: Quaternion): Quaternion;
+    dot(q: Quaternion): number;
+    rotateVector(v: Vector): Vector;
+    rotateVector3D(v: Vector3D): Vector3D;
+    lerp(q: Quaternion, t: number): Quaternion;
+    slerp(q: Quaternion, t: number): Quaternion;
+    toVector3D(): Vector3D;
+    angle(q: Quaternion): number;
+    geodesicAngle(q: Quaternion): number;
+    distance(q: Quaternion): number;
+    delta(qTo: Quaternion): Quaternion;
+    toDirectionCosineMatrix(): Matrix;
+    toRotationMatrix(): Matrix;
+    vectorAngle(observer: Vector3D, target: Vector3D, forward: Vector3D): number;
+    kinematics(angularVelocity: Vector3D): Quaternion;
+}

package/lib/operations/Quaternion.js

@@ -0,0 +1,164 @@
+/**
+ * @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 { Matrix, Vector, Vector3D, wrapAngle } from 'ootk-core';
+export class Quaternion {
+    constructor(x, y, z, w) {
+        this.x = x;
+        this.y = y;
+        this.z = z;
+        this.w = w;
+    }
+    toString(precision = 8) {
+        const xStr = this.x.toFixed(precision);
+        const yStr = this.y.toFixed(precision);
+        const zStr = this.z.toFixed(precision);
+        const wStr = this.w.toFixed(precision);
+        return `Q(x: ${xStr}, y: ${yStr}, z: ${zStr}, w: ${wStr})`;
+    }
+    positivePolar() {
+        return this.w >= 0 ? this.normalize() : this.negate().normalize();
+    }
+    magnitudeSquared() {
+        return this.w * this.w + this.x * this.x + this.y * this.y + this.z * this.z;
+    }
+    magnitude() {
+        return Math.sqrt(this.magnitudeSquared());
+    }
+    scale(n) {
+        return new Quaternion(n * this.x, n * this.y, n * this.z, n * this.w);
+    }
+    negate() {
+        return this.scale(-1);
+    }
+    normalize() {
+        const m = this.magnitude();
+        if (m === 0) {
+            return Quaternion.zero;
+        }
+        return this.scale(1 / m);
+    }
+    conjugate() {
+        return new Quaternion(-this.x, -this.y, -this.z, this.w);
+    }
+    inverse() {
+        return this.conjugate().scale(1 / this.magnitudeSquared());
+    }
+    add(q) {
+        return new Quaternion(this.x + q.x, this.y + q.y, this.z + q.z, this.w + q.w);
+    }
+    subtract(q) {
+        return new Quaternion(this.x - q.x, this.y - q.y, this.z - q.z, this.w - q.w);
+    }
+    addReal(n) {
+        return new Quaternion(this.x, this.y, this.z, this.w + n);
+    }
+    multiply(q) {
+        const mx = this.w * q.x + this.x * q.w + this.y * q.z - this.z * q.y;
+        const my = this.w * q.y - this.x * q.z + this.y * q.w + this.z * q.x;
+        const mz = this.w * q.z + this.x * q.y - this.y * q.x + this.z * q.w;
+        const mw = this.w * q.w - this.x * q.x - this.y * q.y - this.z * q.z;
+        return new Quaternion(mx, my, mz, mw);
+    }
+    dot(q) {
+        return this.x * q.x + this.y * q.y + this.z * q.z + this.w * q.w;
+    }
+    rotateVector(v) {
+        const q = this.multiply(new Quaternion(v.x, v.y, v.z, 0)).multiply(this.conjugate());
+        return new Vector([q.x, q.y, q.z]);
+    }
+    rotateVector3D(v) {
+        const q = this.multiply(new Quaternion(v.x, v.y, v.z, 0)).multiply(this.conjugate());
+        return new Vector3D(q.x, q.y, q.z);
+    }
+    lerp(q, t) {
+        const f = 1.0 - t;
+        return new Quaternion(f * this.x + t * q.x, f * this.y + t * q.y, f * this.z + t * q.z, f * this.w + t * q.w).positivePolar();
+    }
+    slerp(q, t) {
+        let qp = q;
+        let dotP = this.dot(qp);
+        if (dotP < 0) {
+            dotP = -dotP;
+            qp = qp.negate();
+        }
+        if (dotP > 0.9995) {
+            return this.lerp(qp, t);
+        }
+        const theta = Math.acos(dotP);
+        const sinTheta = Math.sin(theta);
+        const f1 = Math.sin((1.0 - t) * theta) / sinTheta;
+        const f2 = Math.sin(t * theta) / sinTheta;
+        return new Quaternion(f1 * this.x + f2 * qp.x, f1 * this.y + f2 * qp.y, f1 * this.z + f2 * qp.z, f1 * this.w + f2 * qp.w).positivePolar();
+    }
+    toVector3D() {
+        return new Vector3D(this.x, this.y, this.z);
+    }
+    angle(q) {
+        const c = this.multiply(q.conjugate()).normalize();
+        return 2 * Math.atan2(c.toVector3D().magnitude(), c.w);
+    }
+    geodesicAngle(q) {
+        const p = this.dot(q);
+        return wrapAngle(Math.acos(2 * p * p - 1.0));
+    }
+    distance(q) {
+        const m01 = this.subtract(q).magnitude();
+        const p01 = this.add(q).magnitude();
+        return m01 < p01 ? m01 : p01;
+    }
+    delta(qTo) {
+        return this.inverse().multiply(qTo);
+    }
+    toDirectionCosineMatrix() {
+        const w2 = this.w * this.w;
+        const x2 = this.x * this.x;
+        const y2 = this.y * this.y;
+        const z2 = this.z * this.z;
+        const m = [
+            [w2 + x2 - y2 - z2, 2.0 * (this.x * this.y + this.z * this.w), 2.0 * (this.x * this.z - this.y * this.w)],
+            [2.0 * (this.x * this.y - this.z * this.w), w2 - x2 + y2 - z2, 2.0 * (this.y * this.z + this.x * this.w)],
+            [2.0 * (this.x * this.z + this.y * this.w), 2.0 * (this.y * this.z - this.x * this.w), w2 - x2 - y2 + z2],
+        ];
+        return new Matrix(m);
+    }
+    toRotationMatrix() {
+        return this.toDirectionCosineMatrix().transpose();
+    }
+    vectorAngle(observer, target, forward) {
+        const delta = target.subtract(observer);
+        const transform = this.toDirectionCosineMatrix().multiplyVector3D(delta);
+        return forward.angle(transform);
+    }
+    kinematics(angularVelocity) {
+        const wPrime = new Vector([0, angularVelocity[0], angularVelocity[1], angularVelocity[2]]);
+        const qMat = new Matrix([
+            [this.x, this.w, -this.z, this.y],
+            [this.y, this.z, this.w, -this.x],
+            [this.z, -this.y, this.x, this.w],
+            [this.w, -this.x, -this.y, -this.z],
+        ]);
+        const result = qMat.multiplyVector(wPrime).scale(0.5);
+        return new Quaternion(result[0], result[1], result[2], result[3]);
+    }
+}
+Quaternion.zero = new Quaternion(0, 0, 0, 0);
+Quaternion.one = new Quaternion(0, 0, 0, 1);
+Quaternion.xAxis = new Quaternion(1, 0, 0, 0);
+Quaternion.yAxis = new Quaternion(0, 1, 0, 0);
+Quaternion.zAxis = new Quaternion(0, 0, 1, 0);
+//# sourceMappingURL=Quaternion.js.map

package/lib/operations/Quaternion.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"Quaternion.js","sourceRoot":"","sources":["../../src/operations/Quaternion.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,EAAE,QAAQ,EAAE,SAAS,EAAE,MAAM,WAAW,CAAC;AAEhE,MAAM,OAAO,UAAU;IAMrB,YAAY,CAAS,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QACpD,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;QACX,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC;IACb,CAAC;IAQD,QAAQ,CAAC,SAAS,GAAG,CAAC;QACpB,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QACvC,MAAM,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,OAAO,CAAC,SAAS,CAAC,CAAC;QAEvC,OAAO,QAAQ,IAAI,QAAQ,IAAI,QAAQ,IAAI,QAAQ,IAAI,GAAG,CAAC;IAC7D,CAAC;IAED,aAAa;QACX,OAAO,IAAI,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,SAAS,EAAE,CAAC;IACpE,CAAC;IAED,gBAAgB;QACd,OAAO,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;IAC/E,CAAC;IAED,SAAS;QACP,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,gBAAgB,EAAE,CAAC,CAAC;IAC5C,CAAC;IAED,KAAK,CAAC,CAAS;QACb,OAAO,IAAI,UAAU,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACxE,CAAC;IAED,MAAM;QACJ,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;IACxB,CAAC;IAED,SAAS;QACP,MAAM,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;QAE3B,IAAI,CAAC,KAAK,CAAC,EAAE;YACX,OAAO,UAAU,CAAC,IAAI,CAAC;SACxB;QAED,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC3B,CAAC;IAED,SAAS;QACP,OAAO,IAAI,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC;IAED,OAAO;QACL,OAAO,IAAI,CAAC,SAAS,EAAE,CAAC,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,gBAAgB,EAAE,CAAC,CAAC;IAC7D,CAAC;IAED,GAAG,CAAC,CAAa;QACf,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC;IAED,QAAQ,CAAC,CAAa;QACpB,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChF,CAAC;IAED,OAAO,CAAC,CAAS;QACf,OAAO,IAAI,UAAU,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC5D,CAAC;IAED,QAAQ,CAAC,CAAa;QACpB,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACrE,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACrE,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACrE,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAErE,OAAO,IAAI,UAAU,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;IACxC,CAAC;IAED,GAAG,CAAC,CAAa;QACf,OAAO,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACnE,CAAC;IAED,YAAY,CAAC,CAAS;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC;QAErF,OAAO,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IAED,cAAc,CAAC,CAAW;QACxB,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,IAAI,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC;QAErF,OAAO,IAAI,QAAQ,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACrC,CAAC;IAED,IAAI,CAAC,CAAa,EAAE,CAAS;QAC3B,MAAM,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QAElB,OAAO,IAAI,UAAU,CACnB,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EACpB,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EACpB,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EACpB,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CACrB,CAAC,aAAa,EAAE,CAAC;IACpB,CAAC;IAED,KAAK,CAAC,CAAa,EAAE,CAAS;QAC5B,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC;QAExB,IAAI,IAAI,GAAG,CAAC,EAAE;YACZ,IAAI,GAAG,CAAC,IAAI,CAAC;YACb,EAAE,GAAG,EAAE,CAAC,MAAM,EAAE,CAAC;SAClB;QACD,IAAI,IAAI,GAAG,MAAM,EAAE;YACjB,OAAO,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;SACzB;QACD,MAAM,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC;QAC9B,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,QAAQ,CAAC;QAClD,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,QAAQ,CAAC;QAE1C,OAAO,IAAI,UAAU,CACnB,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EACvB,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EACvB,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,EACvB,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CACxB,CAAC,aAAa,EAAE,CAAC;IACpB,CAAC;IAED,UAAU;QACR,OAAO,IAAI,QAAQ,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9C,CAAC;IAED,KAAK,CAAC,CAAa;QACjB,MAAM,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC,SAAS,EAAE,CAAC;QAEnD,OAAO,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,UAAU,EAAE,CAAC,SAAS,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IACzD,CAAC;IAED,aAAa,CAAC,CAAa;QACzB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAEtB,OAAO,SAAS,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC/C,CAAC;IAED,QAAQ,CAAC,CAAa;QACpB,MAAM,GAAG,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QACzC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC;QAEpC,OAAO,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC;IAC/B,CAAC;IAED,KAAK,CAAC,GAAe;QACnB,OAAO,IAAI,CAAC,OAAO,EAAE,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC;IACtC,CAAC;IAED,uBAAuB;QACrB,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC;QAC3B,MAAM,CAAC,GAAG;YACR,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACzG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;YACzG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;SAC1G,CAAC;QAEF,OAAO,IAAI,MAAM,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED,gBAAgB;QACd,OAAO,IAAI,CAAC,uBAAuB,EAAE,CAAC,SAAS,EAAE,CAAC;IACpD,CAAC;IAED,WAAW,CAAC,QAAkB,EAAE,MAAgB,EAAE,OAAiB;QACjE,MAAM,KAAK,GAAG,MAAM,CAAC,QAAQ,CAAC,QAAQ,CAAC,CAAC;QACxC,MAAM,SAAS,GAAG,IAAI,CAAC,uBAAuB,EAAE,CAAC,gBAAgB,CAAC,KAAK,CAAC,CAAC;QAEzE,OAAO,OAAO,CAAC,KAAK,CAAC,SAAS,CAAC,CAAC;IAClC,CAAC;IAED,UAAU,CAAC,eAAyB;QAClC,MAAM,MAAM,GAAG,IAAI,MAAM,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3F,MAAM,IAAI,GAAG,IAAI,MAAM,CAAC;YACtB,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;YACjC,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;YACjC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC;YACjC,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC;SACpC,CAAC,CAAC;QACH,MAAM,MAAM,GAAG,IAAI,CAAC,cAAc,CAAC,MAAM,CAAC,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;QAEtD,OAAO,IAAI,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACpE,CAAC;;AA5LM,eAAI,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AAClC,cAAG,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACjC,gBAAK,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACnC,gBAAK,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;AACnC,gBAAK,GAAG,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC"}

package/lib/operations/RandomGaussianSource.d.ts

@@ -0,0 +1,25 @@
+/**
+ * @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, Vector3D } from 'ootk-core';
+export declare class RandomGaussianSource {
+    private _boxMuller;
+    constructor(seed?: number);
+    nextGauss(): number;
+    gaussVector(n: number): Vector;
+    gaussSphere(radius?: number): Vector3D;
+}

package/lib/operations/RandomGaussianSource.js

@@ -0,0 +1,43 @@
+/**
+ * @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';
+import { BoxMuller } from './BoxMuller';
+export class RandomGaussianSource {
+    constructor(seed = 0) {
+        this._boxMuller = new BoxMuller(0, 1, seed);
+    }
+    nextGauss() {
+        return this._boxMuller.nextGauss();
+    }
+    gaussVector(n) {
+        if (n < 1) {
+            throw new Error('n must be greater than 0');
+        }
+        const result = new Vector([this.nextGauss()]);
+        for (let i = 0; i < n; i++) {
+            if (i > 0) {
+                result.add(new Vector[this.nextGauss()]());
+            }
+        }
+        return result;
+    }
+    gaussSphere(radius = 1.0) {
+        return this.gaussVector(3).toVector3D(0).normalize().scale(radius);
+    }
+}
+//# sourceMappingURL=RandomGaussianSource.js.map

package/lib/operations/RandomGaussianSource.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"RandomGaussianSource.js","sourceRoot":"","sources":["../../src/operations/RandomGaussianSource.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;GAgBG;AAEH,OAAO,EAAE,MAAM,EAAY,MAAM,WAAW,CAAC;AAC7C,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AAExC,MAAM,OAAO,oBAAoB;IAG/B,YAAY,IAAI,GAAG,CAAC;QAClB,IAAI,CAAC,UAAU,GAAG,IAAI,SAAS,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;IAC9C,CAAC;IAED,SAAS;QACP,OAAO,IAAI,CAAC,UAAU,CAAC,SAAS,EAAE,CAAC;IACrC,CAAC;IAED,WAAW,CAAC,CAAS;QACnB,IAAI,CAAC,GAAG,CAAC,EAAE;YACT,MAAM,IAAI,KAAK,CAAC,0BAA0B,CAAC,CAAC;SAC7C;QAED,MAAM,MAAM,GAAG,IAAI,MAAM,CAAC,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QAE9C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;YAC1B,IAAI,CAAC,GAAG,CAAC,EAAE;gBACT,MAAM,CAAC,GAAG,CAAC,IAAI,MAAM,CAAC,IAAI,CAAC,SAAS,EAAE,CAAC,EAAE,CAAC,CAAC;aAC5C;SACF;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,WAAW,CAAC,MAAM,GAAG,GAAG;QACtB,OAAO,IAAI,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;IACrE,CAAC;CACF"}

package/lib/operations/index.d.ts

@@ -0,0 +1,4 @@
+export { BoxMuller } from './BoxMuller';
+export { EulerAngles } from './EulerAngles';
+export { Quaternion } from './Quaternion';
+export { RandomGaussianSource } from './RandomGaussianSource';

package/lib/operations/index.js

@@ -0,0 +1,5 @@
+export { BoxMuller } from './BoxMuller';
+export { EulerAngles } from './EulerAngles';
+export { Quaternion } from './Quaternion';
+export { RandomGaussianSource } from './RandomGaussianSource';
+//# sourceMappingURL=index.js.map

package/lib/operations/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/operations/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,SAAS,EAAE,MAAM,aAAa,CAAC;AACxC,OAAO,EAAE,WAAW,EAAE,MAAM,eAAe,CAAC;AAC5C,OAAO,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAC1C,OAAO,EAAE,oBAAoB,EAAE,MAAM,wBAAwB,CAAC"}

package/lib/optimize/ChebyshevCompressor.d.ts

@@ -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/>.
+ */
+import { StateInterpolator } from '../interpolator/StateInterpolator';
+import { ChebyshevInterpolator } from './../interpolator/ChebyshevInterpolator';
+export declare class ChebyshevCompressor {
+    private _interpolator;
+    constructor(_interpolator: StateInterpolator);
+    private static _cosPi;
+    private _fitCoefficient;
+    private _fitWindow;
+    /**
+     * Compress this object's interpolater, using the provided coefficients
+     * per revolution [cpr].
+     */
+    compress(cpr?: number): ChebyshevInterpolator;
+}