@woosh/meep-engine 2.92.23 → 2.93.1

package/build/meep.cjs

@@ -3145,6 +3145,9 @@ const forward = new Vector3$1();
 const up = new Vector3$1();
 const right = new Vector3$1();
 
+const sin$1 = Math.sin;
+const cos$1 = Math.cos;
+
 let Quaternion$1 = class Quaternion {
     /**
      *
@@ -3180,6 +3183,8 @@ let Quaternion$1 = class Quaternion {
         this.onChanged = new Signal();
     }
 
+    // Making Quaternion comply to array interface
+
     get 0() {
         return this.x;
     }
@@ -3212,6 +3217,19 @@ let Quaternion$1 = class Quaternion {
         this.w = v;
     }
 
+    /**
+     * Making quaternion iterable
+     */
+    * [Symbol.iterator]() {
+
+        yield this.x;
+        yield this.y;
+        yield this.z;
+        yield this.w;
+
+    }
+
+
     /**
      *
      * @param {number} fx forward vector
@@ -3252,58 +3270,22 @@ let Quaternion$1 = class Quaternion {
 
         // construct partial transform matrix
         const m00 = right.x;
-        const m01 = right.y;
-        const m02 = right.z;
-        const m10 = up.x;
-        const m11 = up.y;
-        const m12 = up.z;
-        const m20 = forward.x;
-        const m21 = forward.y;
-        const m22 = forward.z;
-
-
-        const num8 = (m00 + m11) + m22;
-
-        let _x, _y, _z, _w;
-
-        if (num8 > 0) {
-            let num = Math.sqrt(num8 + 1);
-            _w = num * 0.5;
-            num = 0.5 / num;
-            _x = (m12 - m21) * num;
-            _y = (m20 - m02) * num;
-            _z = (m01 - m10) * num;
-        } else if ((m00 >= m11) && (m00 >= m22)) {
-
-            const num7 = Math.sqrt(((1 + m00) - m11) - m22);
-            const num4 = 0.5 / num7;
-
-            _x = 0.5 * num7;
-            _y = (m01 + m10) * num4;
-            _z = (m02 + m20) * num4;
-            _w = (m12 - m21) * num4;
-
-        } else if (m11 > m22) {
+        const m10 = right.y;
+        const m20 = right.z;
 
-            const num6 = Math.sqrt(((1 + m11) - m00) - m22);
-            const num3 = 0.5 / num6;
-
-            _x = (m10 + m01) * num3;
-            _y = 0.5 * num6;
-            _z = (m21 + m12) * num3;
-            _w = (m20 - m02) * num3;
-        } else {
-
-            const num5 = Math.sqrt(((1 + m22) - m00) - m11);
-            const num2 = 0.5 / num5;
+        const m01 = up.x;
+        const m11 = up.y;
+        const m21 = up.z;
 
-            _x = (m20 + m02) * num2;
-            _y = (m21 + m12) * num2;
-            _z = 0.5 * num5;
-            _w = (m01 - m10) * num2;
-        }
+        const m02 = forward.x;
+        const m12 = forward.y;
+        const m22 = forward.z;
 
-        this.set(_x, _y, _z, _w);
+        this.__setFromRotationMatrix(
+            m00, m01, m02,
+            m10, m11, m12,
+            m20, m21, m22
+        );
     }
 
     /**
@@ -3311,13 +3293,24 @@ let Quaternion$1 = class Quaternion {
      * @param {Vector3} vForward
      * @param {Vector3} [vUp=Vector3.up]
      *
-     * @source http://answers.unity3d.com/questions/467614/what-is-the-source-code-of-quaternionlookrotation.html
      */
     lookRotation(vForward, vUp = Vector3$1.up) {
 
         this._lookRotation(vForward.x, vForward.y, vForward.z, vUp.x, vUp.y, vUp.z);
     }
 
+    /**
+     *
+     * @param {Quaternion} other
+     * @return {number}
+     */
+    dot(other) {
+        return this.x * other.x
+            + this.y * other.y
+            + this.z * other.z
+            + this.w * other.w
+            ;
+    }
 
     /**
      *
@@ -3337,14 +3330,14 @@ let Quaternion$1 = class Quaternion {
         const z = this.z;
         const w = this.w;
 
-        const dot = x * x + y * y + z * z + w * w;
+        const length_sqr = x * x + y * y + z * z + w * w;
 
-        if (dot === 0) {
+        if (length_sqr === 0) {
             this.set(0, 0, 0, 0);
             return;
         }
 
-        const invDot = 1.0 / dot;
+        const invDot = 1.0 / length_sqr;
 
         const _x = -x * invDot;
         const _y = -y * invDot;
@@ -3398,8 +3391,8 @@ let Quaternion$1 = class Quaternion {
     _fromAxisAngle(ax, ay, az, angle) {
         const halfAngle = angle * 0.5;
 
-        const sinA2 = Math.sin(halfAngle);
-        const cosA2 = Math.cos(halfAngle);
+        const sinA2 = sin$1(halfAngle);
+        const cosA2 = cos$1(halfAngle);
 
         const qx = ax * sinA2;
         const qy = ay * sinA2;
@@ -3480,7 +3473,7 @@ let Quaternion$1 = class Quaternion {
     toAxisAngle(axis) {
         const rad = Math.acos(this.w) * 2.0;
 
-        const s = Math.sin(rad * 0.5);
+        const s = sin$1(rad * 0.5);
 
         if (Math.abs(s) > EPSILON) {
             axis.set(
@@ -3523,7 +3516,6 @@ let Quaternion$1 = class Quaternion {
     }
 
     /**
-     * @see http://www.euclideanspace.com/maths/algebra/realNormedAlgebra/quaternions/code/index.htm
      * @param {Quaternion} other
      */
     multiply(other) {
@@ -3563,6 +3555,8 @@ let Quaternion$1 = class Quaternion {
      */
     _multiplyQuaternions(ax, ay, az, aw, bx, by, bz, bw) {
 
+        // see http://www.euclideanspace.com/maths/algebra/realNormedAlgebra/quaternions/code/index.htm
+
         const x = ax * bw + aw * bx + ay * bz - az * by;
         const y = ay * bw + aw * by + az * bx - ax * bz;
         const z = az * bw + aw * bz + ax * by - ay * bx;
@@ -3652,7 +3646,7 @@ let Quaternion$1 = class Quaternion {
      * @param {number} x
      * @param {number} y
      * @param {number} z
-     * @param {String} order a combination of capital letters X,Y,Z. Examples: XYZ, YXZ
+     * @param {String} [order] a combination of capital letters X,Y,Z. Examples: XYZ, YXZ
      * @returns {Quaternion}
      */
     __setFromEuler(x, y, z, order = 'XYZ') {
@@ -3683,6 +3677,8 @@ let Quaternion$1 = class Quaternion {
 
             this.fromEulerAnglesXZY(x, y, z);
 
+        } else {
+            throw new Error(`Invalid order '${order}', bust be 3 capital letters consisting of X,Y and Z`);
         }
 
         return this;
@@ -3807,16 +3803,13 @@ let Quaternion$1 = class Quaternion {
         const half_y = y * 0.5;
         const half_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
+        const s1 = sin$1(half_x);
+        const s2 = sin$1(half_y);
+        const s3 = sin$1(half_z);
 
-        const s1 = sin(half_x);
-        const s2 = sin(half_y);
-        const s3 = sin(half_z);
-
-        const c1 = cos(half_x);
-        const c2 = cos(half_y);
-        const c3 = cos(half_z);
+        const c1 = cos$1(half_x);
+        const c2 = cos$1(half_y);
+        const c3 = cos$1(half_z);
 
         const _x = s1 * c2 * c3 + c1 * s2 * s3;
         const _y = c1 * s2 * c3 - s1 * c2 * s3;
@@ -3840,16 +3833,13 @@ let Quaternion$1 = class Quaternion {
         const scaled_y = y * 0.5;
         const scaled_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
-
-        const s1 = sin(scaled_x);
-        const s2 = sin(scaled_y);
-        const s3 = sin(scaled_z);
+        const s1 = sin$1(scaled_x);
+        const s2 = sin$1(scaled_y);
+        const s3 = sin$1(scaled_z);
 
-        const c1 = cos(scaled_x);
-        const c2 = cos(scaled_y);
-        const c3 = cos(scaled_z);
+        const c1 = cos$1(scaled_x);
+        const c2 = cos$1(scaled_y);
+        const c3 = cos$1(scaled_z);
 
         const _x = s1 * c2 * c3 + c1 * s2 * s3;
         const _y = c1 * s2 * c3 - s1 * c2 * s3;
@@ -3873,16 +3863,13 @@ let Quaternion$1 = class Quaternion {
         const scaled_y = y * 0.5;
         const scaled_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
-
-        const s1 = sin(scaled_x);
-        const s2 = sin(scaled_y);
-        const s3 = sin(scaled_z);
+        const s1 = sin$1(scaled_x);
+        const s2 = sin$1(scaled_y);
+        const s3 = sin$1(scaled_z);
 
-        const c1 = cos(scaled_x);
-        const c2 = cos(scaled_y);
-        const c3 = cos(scaled_z);
+        const c1 = cos$1(scaled_x);
+        const c2 = cos$1(scaled_y);
+        const c3 = cos$1(scaled_z);
 
         const _x = s1 * c2 * c3 - c1 * s2 * s3;
         const _y = c1 * s2 * c3 + s1 * c2 * s3;
@@ -3906,16 +3893,13 @@ let Quaternion$1 = class Quaternion {
         const scaled_y = y * 0.5;
         const scaled_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
+        const s1 = sin$1(scaled_x);
+        const s2 = sin$1(scaled_y);
+        const s3 = sin$1(scaled_z);
 
-        const s1 = sin(scaled_x);
-        const s2 = sin(scaled_y);
-        const s3 = sin(scaled_z);
-
-        const c1 = cos(scaled_x);
-        const c2 = cos(scaled_y);
-        const c3 = cos(scaled_z);
+        const c1 = cos$1(scaled_x);
+        const c2 = cos$1(scaled_y);
+        const c3 = cos$1(scaled_z);
 
         const _x = s1 * c2 * c3 - c1 * s2 * s3;
         const _y = c1 * s2 * c3 + s1 * c2 * s3;
@@ -3939,16 +3923,13 @@ let Quaternion$1 = class Quaternion {
         const scaled_y = y * 0.5;
         const scaled_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
-
-        const s1 = sin(scaled_x);
-        const s2 = sin(scaled_y);
-        const s3 = sin(scaled_z);
+        const s1 = sin$1(scaled_x);
+        const s2 = sin$1(scaled_y);
+        const s3 = sin$1(scaled_z);
 
-        const c1 = cos(scaled_x);
-        const c2 = cos(scaled_y);
-        const c3 = cos(scaled_z);
+        const c1 = cos$1(scaled_x);
+        const c2 = cos$1(scaled_y);
+        const c3 = cos$1(scaled_z);
 
         const _x = s1 * c2 * c3 + c1 * s2 * s3;
         const _y = c1 * s2 * c3 + s1 * c2 * s3;
@@ -3972,16 +3953,13 @@ let Quaternion$1 = class Quaternion {
         const scaled_y = y * 0.5;
         const scaled_z = z * 0.5;
 
-        const sin = Math.sin;
-        const cos = Math.cos;
+        const s1 = sin$1(scaled_x);
+        const s2 = sin$1(scaled_y);
+        const s3 = sin$1(scaled_z);
 
-        const s1 = sin(scaled_x);
-        const s2 = sin(scaled_y);
-        const s3 = sin(scaled_z);
-
-        const c1 = cos(scaled_x);
-        const c2 = cos(scaled_y);
-        const c3 = cos(scaled_z);
+        const c1 = cos$1(scaled_x);
+        const c2 = cos$1(scaled_y);
+        const c3 = cos$1(scaled_z);
 
         const _x = s1 * c2 * c3 - c1 * s2 * s3;
         const _y = c1 * s2 * c3 - s1 * c2 * s3;
@@ -4155,7 +4133,7 @@ let Quaternion$1 = class Quaternion {
     }
 
     /**
-     *
+     * @deprecated
      * @param {Matrix4} m
      */
     setFromRotationMatrix(m) {
@@ -4311,9 +4289,9 @@ let Quaternion$1 = class Quaternion {
         if ((1.0 - cosom) > EPSILON) {
             // standard case (slerp)
             omega = Math.acos(cosom);
-            sinom = Math.sin(omega);
-            scale0 = Math.sin((1.0 - t) * omega) / sinom;
-            scale1 = Math.sin(t * omega) / sinom;
+            sinom = sin$1(omega);
+            scale0 = sin$1((1.0 - t) * omega) / sinom;
+            scale1 = sin$1(t * omega) / sinom;
         } else {
             // "from" and "to" quaternions are very close
             //  ... so we can do a linear interpolation
@@ -4569,15 +4547,6 @@ let Quaternion$1 = class Quaternion {
             && epsilonEquals(this.w, w, tolerance);
     }
 
-    * [Symbol.iterator]() {
-
-        yield this.x;
-        yield this.y;
-        yield this.z;
-        yield this.w;
-
-    }
-
     /**
      * Based on http://planning.cs.uiuc.edu/node198.html
      * @param {function():number} [random]
@@ -4594,10 +4563,10 @@ let Quaternion$1 = class Quaternion {
         const u3 = 2 * Math.PI * random();
 
         return this.set(
-            sqrt1u1 * Math.cos(u2),
-            sqrtu1 * Math.sin(u3),
-            sqrtu1 * Math.cos(u3),
-            sqrt1u1 * Math.sin(u2),
+            sqrt1u1 * cos$1(u2),
+            sqrtu1 * sin$1(u3),
+            sqrtu1 * cos$1(u3),
+            sqrt1u1 * sin$1(u2),
         );
     }
 
@@ -4645,6 +4614,7 @@ let Quaternion$1 = class Quaternion {
         const angle = from.angleTo(to);
 
         if (angle === 0) {
+            // we're already where we need tobe. Also - avoid division by 0
             result.copy(to);
         } else {
             // clamp to 1, to make sure we don't overshoot