@woosh/meep-engine 2.109.1 → 2.109.3

package/src/core/geom/Quaternion.js

@@ -9,6 +9,7 @@ import { computeHashFloat } from "../primitives/numbers/computeHashFloat.js";
 import { quat_decode_from_uint32 } from "./3d/quaternion/quat_decode_from_uint32.js";
 import { quat_encode_to_uint32 } from "./3d/quaternion/quat_encode_to_uint32.js";
 import { v3_dot } from "./vec3/v3_dot.js";
+import { v4_length } from "./vec4/v4_length.js";
 import Vector3 from "./Vector3.js";
 
 
@@ -451,7 +452,7 @@ class Quaternion {
         const z = this.z;
         const w = this.w;
 
-        return Math.sqrt(x * x + y * y + z * z + w * w);
+        return v4_length(x, y, z, w);
     }
 
     /**
@@ -463,48 +464,20 @@ class Quaternion {
         Quaternion.rotateTowards(this, this, other, max_delta);
     }
 
-    /**
-     *
-     * @param {Vector3} direction
-     * @param {Vector3} [up]
-     * @param {Vector3} [forward]
-     */
-    alignToDirection(direction, up = Vector3.up, forward = Vector3.forward) {
-
-        const dot = up.dot(direction);
-
-
-        if (Math.abs(dot + 1) < 0.000001) {
-            // vector a and b point exactly in the opposite direction,
-            // so it is a 180 degrees turn around the up-axis
-            this.set(forward.x, forward.y, forward.z, Math.PI);
-        } else if (Math.abs(dot - 1) < 0.000001) {
-            // vector a and b point exactly in the same direction
-            // so we return the identity quaternion
-            this.set(0, 0, 0, 1);
-        }
-
-        const angle = Math.acos(dot);
-
-        axis.crossVectors(up, direction);
-        axis.normalize();
-
-        this.fromAxisAngle(axis, angle);
-    }
-
     /**
      *
      * @param {Vector3} source
      * @param {Vector3} target
+     * @param {Vector3} [up]
      */
-    lookAt(source, target) {
+    lookAt(source, target, up = Vector3.up) {
         const forward = new Vector3();
 
         forward.subVectors(target, source);
 
         forward.normalize();
 
-        this.alignToDirection(forward);
+        this.lookRotation(forward, up);
     }
 
     /**
@@ -887,22 +860,23 @@ class Quaternion {
         const qW_1 = qW + l0;
 
         //normalize result
-        const l1 = Math.sqrt(p0 + qW_1 * qW_1);
+        const l1 = 1 / Math.sqrt(p0 + qW_1 * qW_1);
 
-        const x = qX / l1;
-        const y = qY / l1;
-        const z = qZ / l1;
-        const w = qW / l1;
+        const x = qX * l1;
+        const y = qY * l1;
+        const z = qZ * l1;
+        const w = qW * l1;
 
         this.set(x, y, z, w);
     }
 
     /**
+     * NOTE: Vectors need to be normalized
      * Based on blog post: http://lolengine.net/blog/2013/09/18/beautiful-maths-quaternion-from-vectors
      * @param {Vector3} from
      * @param {Vector3} to
      */
-    fromUnitVectors2(from, to) {
+    fromUnitVectors(from, to) {
         assert.ok(from.isNormalized(), `from vector is not normalized, length = ${from.length()}`);
         assert.ok(to.isNormalized(), `to vector is not normalized, length = ${to.length()}`);
 
@@ -923,7 +897,7 @@ class Quaternion {
 
         const uv_dot = v3_dot(ax, ay, az, bx, by, bz);
 
-        if (uv_dot === -1) {
+        if (uv_dot < -0.9999999) {
             //to vector is opposite, produce a reversal quaternion
 
             tempvec3.crossVectors(Vector3.left, from);
@@ -967,51 +941,6 @@ class Quaternion {
         );
     }
 
-    /**
-     * Input vectors must be normalized
-     * @param {Vector3} from
-     * @param {Vector3} to
-     */
-    fromUnitVectors(from, to) {
-        assert.ok(from.isNormalized(), `from vector is not normalized, length = ${from.length()}`);
-        assert.ok(to.isNormalized(), `to vector is not normalized, length = ${to.length()}`);
-
-        const dot = from.dot(to);
-
-        let x, y, z, w;
-
-        if (dot < -0.9999999) {
-            const tempvec3 = new Vector3();
-
-            tempvec3.crossVectors(Vector3.left, from);
-
-            if (tempvec3.lengthSqr() < 0.00001) {
-                tempvec3.crossVectors(Vector3.up, from);
-            }
-
-            tempvec3.normalize();
-
-            this.fromAxisAngle(tempvec3, Math.PI);
-        } else if (dot > 0.9999999) {
-            x = 0;
-            y = 0;
-            z = 0;
-            w = 1;
-
-            this.set(x, y, z, w);
-        } else {
-            //cross product from x to
-
-            x = from.y * to.z - from.z * to.y;
-            y = from.z * to.x - from.x * to.z;
-            z = from.x * to.y - from.y * to.x;
-            w = 1 + dot;
-
-            this.set(x, y, z, w);
-            this.normalize();
-        }
-    }
-
     /**
      * @deprecated
      * @param {Matrix4} m
@@ -1216,17 +1145,6 @@ class Quaternion {
         this.slerpQuaternions(this, other, t);
     }
 
-    /**
-     * @see https://github.com/gareth-cross/quat/blob/master/include/quaternion.hpp
-     * TODO implement
-     * @param {Quaternion} a
-     * @param {Quaternion} b
-     * @param {number} t
-     */
-    quaternionsInterpolateRungeKutta4(a, b, t) {
-        throw new Error('Not Implemented');
-    }
-
     /**
      *
      * @param {function} handler
@@ -1557,13 +1475,6 @@ Quaternion.prototype.toArray = Quaternion.prototype.writeToArray;
 Quaternion.identity = Object.freeze(new Quaternion(0, 0, 0, 1));
 
 
-/**
- *
- * @type {Vector3}
- */
-const axis = new Vector3();
-
-
 const tempvec3 = new Vector3();
 
 export default Quaternion;

package/src/core/geom/Quaternion.spec.js

@@ -70,14 +70,14 @@ test('multiply quaternions', () => {
     });
 });
 
-test('fromUnitVectors2', () => {
+test('fromUnitVectors', () => {
     const q = new Quaternion();
 
     const v = new Vector3();
 
     function check(v0, v1) {
 
-        q.fromUnitVectors2(v0, v1);
+        q.fromUnitVectors(v0, v1);
 
         v.copy(v0);
 

package/src/core/geom/Vector1.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Vector1.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector1.js"],"names":[],"mappings":";AAWA;IA2PI;;;;;OAKG;IACH,kBAJW,OAAO,KACP,OAAO,GACN,MAAM,CAIjB;IAlQD;;;;;;OAMG;IACH,gBALW,MAAM,EAchB;IAHG,UAAU;IAEV,0DAA6B;IAWjC;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,iBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,OAHW,MAAM,GACJ,OAAO,CAgBnB;IAED;;;OAGG;IACH,aAFW,MAAM,QAOhB;IAED;;;OAGG;IACH,UAFY,OAAO,CAIlB;IAED,kBAEC;IAED,kBAEC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;OAGG;IACH,WAFW,OAAO,8BAAwB,WAIzC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;;OAKG;IACH,WAJW,MAAM,QACN,MAAM,GACJ,OAAO,CAInB;IAED;;;OAGG;IACH,YAFW,OAAO,8BAAwB,QAIzC;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;OAEG;IACH,QAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,4BAGC;IAED,iBAEC;IAED,yBAEC;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED,oBAEC;IAED;;;OAGG;IACH,2CAEC;IAED;;;OAGG;IACH,6CAIC;IAgBD,mBAEC;IAND,gBAEC;IAaL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAXvB,sDAIC;CACJ;;;cAYS,OAAO;aAMP,OAAO;;mBAzSE,4BAA4B"}
+{"version":3,"file":"Vector1.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector1.js"],"names":[],"mappings":";AAWA;IAmQI;;;;;OAKG;IACH,kBAJW,OAAO,KACP,OAAO,GACN,MAAM,CAIjB;IA1QD;;;;;;OAMG;IACH,gBALW,MAAM,EAchB;IAHG,UAAU;IAEV,0DAA6B;IAWjC;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,iBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,OAHW,MAAM,GACJ,OAAO,CAgBnB;IAED;;;OAGG;IACH,aAFW,MAAM,QAOhB;IAED;;;OAGG;IACH,UAFY,OAAO,CAIlB;IAED,kBAEC;IAED,kBAEC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;OAGG;IACH,WAFW,OAAO,8BAAwB,WAIzC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;;OAKG;IACH,WAJW,MAAM,QACN,MAAM,GACJ,OAAO,CAInB;IAED;;;OAGG;IACH,eAEC;IAED;;;OAGG;IACH,YAFW,OAAO,8BAAwB,QAIzC;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;OAEG;IACH,QAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,4BAGC;IAED,iBAEC;IAED,yBAEC;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED,oBAEC;IAED;;;OAGG;IACH,2CAEC;IAED;;;OAGG;IACH,6CAIC;IAgBD,mBAEC;IAND,gBAEC;IAaL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAXvB,sDAIC;CACJ;;;cAYS,OAAO;aAMP,OAAO;;mBAjTE,4BAA4B"}

package/src/core/geom/Vector1.js

@@ -20,7 +20,7 @@ class Vector1 extends Number {
     constructor(x = 0) {
         super();
 
-        assert.isNumber( x, 'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         this.x = x;
@@ -67,7 +67,7 @@ class Vector1 extends Number {
      * @returns {Vector1}
      */
     set(x) {
-        assert.isNumber( x,'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         const oldValue = this.x;
@@ -87,7 +87,7 @@ class Vector1 extends Number {
      * @param {number} x
      */
     setSilent(x) {
-        assert.isNumber( x,'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         this.x = x;
@@ -171,6 +171,14 @@ class Vector1 extends Number {
         return this.set(clamp(this.x, low, high));
     }
 
+    /**
+     * Negate sign of the value, if value is 7, it will become -7 and vice versa
+     * Zero is unaffected
+     */
+    negate() {
+        this.set(-this.x);
+    }
+
     /**
      *
      * @param {Vector1|Vector2|Vector3|Vector4} other
@@ -225,7 +233,7 @@ class Vector1 extends Number {
      * @param {number[]} array
      * @param {number} offset
      */
-    readFromArray(array, offset=0) {
+    readFromArray(array, offset = 0) {
         this.set(array[offset]);
     }
 
@@ -234,11 +242,11 @@ class Vector1 extends Number {
      * @param {number[]} array
      * @param {number} offset
      */
-    writeToArray(array, offset=0) {
+    writeToArray(array, offset = 0) {
         array[offset] = this.x;
     }
 
-    asArray(){
+    asArray() {
         return [this.x];
     }
 

package/src/core/geom/vec3/v3_array_displace_in_direction.d.ts

@@ -0,0 +1,14 @@
+/**
+ *
+ * @param {number[]|Float32Array} output
+ * @param {number} output_offset
+ * @param {number} origin_x
+ * @param {number} origin_y
+ * @param {number} origin_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
+ * @param {number} distance
+ */
+export function v3_array_displace_in_direction(output: number[] | Float32Array, output_offset: number, origin_x: number, origin_y: number, origin_z: number, direction_x: number, direction_y: number, direction_z: number, distance: number): void;
+//# sourceMappingURL=v3_array_displace_in_direction.d.ts.map

package/src/core/geom/vec3/v3_array_displace_in_direction.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v3_array_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_array_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,uDAVW,MAAM,EAAE,GAAC,YAAY,iBACrB,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,YACN,MAAM,QAYhB"}

package/src/core/geom/vec3/{v3_shift_along_direction.js → v3_array_displace_in_direction.js}

@@ -10,7 +10,7 @@
  * @param {number} direction_z
  * @param {number} distance
  */
-export function v3_shift_along_direction(
+export function v3_array_displace_in_direction(
     output, output_offset,
     origin_x, origin_y, origin_z,
     direction_x, direction_y, direction_z,

package/src/core/geom/vec3/v3_displace_in_direction.d.ts

@@ -2,13 +2,13 @@
  *
  * @param {Vector3} result
  * @param {number} distance
- * @param {number} sourceX
- * @param {number} sourceY
- * @param {number} sourceZ
- * @param {number} directionX
- * @param {number} directionY
- * @param {number} directionZ
+ * @param {number} origin_x
+ * @param {number} source_y
+ * @param {number} source_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
  * @returns {boolean} true if displaced successfully, false if direction had 0 magnitude
  */
-export function v3_displace_in_direction(result: Vector3, distance: number, sourceX: number, sourceY: number, sourceZ: number, directionX: number, directionY: number, directionZ: number): boolean;
+export function v3_displace_in_direction(result: Vector3, distance: number, origin_x: number, source_y: number, source_z: number, direction_x: number, direction_y: number, direction_z: number): boolean;
 //# sourceMappingURL=v3_displace_in_direction.d.ts.map

package/src/core/geom/vec3/v3_displace_in_direction.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"v3_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,oEATW,MAAM,WACN,MAAM,WACN,MAAM,WACN,MAAM,cACN,MAAM,cACN,MAAM,cACN,MAAM,GACJ,OAAO,CAgCnB"}
+{"version":3,"file":"v3_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,oEATW,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,GACJ,OAAO,CAmCnB"}

package/src/core/geom/vec3/v3_displace_in_direction.js

@@ -2,40 +2,43 @@
  *
  * @param {Vector3} result
  * @param {number} distance
- * @param {number} sourceX
- * @param {number} sourceY
- * @param {number} sourceZ
- * @param {number} directionX
- * @param {number} directionY
- * @param {number} directionZ
+ * @param {number} origin_x
+ * @param {number} source_y
+ * @param {number} source_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
  * @returns {boolean} true if displaced successfully, false if direction had 0 magnitude
  */
 export function v3_displace_in_direction(
     result,
     distance,
-    sourceX, sourceY, sourceZ,
-    directionX, directionY, directionZ
+    origin_x, source_y, source_z,
+    direction_x, direction_y, direction_z
 ) {
+    
     if (distance === 0) {
-        result.set(sourceX, sourceY, sourceZ);
+        result.set(origin_x, source_y, source_z);
         return true;
     }
 
-    //scale direction to have length equal to distance
+    // scale direction to have length equal to distance
 
-    const directionLength = Math.sqrt(directionX * directionX + directionY * directionY + directionZ * directionZ);
+    const direction_length_sqr = direction_x * direction_x + direction_y * direction_y + direction_z * direction_z;
 
-    if(directionLength === 0){
+    if (direction_length_sqr === 0) {
         // direction vector is zero length, avoid division by zero, clamp to source instead
-        result.set(sourceX, sourceY, sourceZ);
+        result.set(origin_x, source_y, source_z);
         return false;
     }
 
-    const m = distance / directionLength;
+    const direction_length = Math.sqrt(direction_length_sqr);
 
-    const x = sourceX + directionX * m;
-    const y = sourceY + directionY * m;
-    const z = sourceZ + directionZ * m;
+    const m = distance / direction_length;
+
+    const x = origin_x + direction_x * m;
+    const y = source_y + direction_y * m;
+    const z = source_z + direction_z * m;
 
     result.set(x, y, z);
 

package/src/core/geom/vec3/v3_displace_in_direction.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=v3_displace_in_direction.spec.d.ts.map

package/src/core/geom/vec3/v3_displace_in_direction.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v3_displace_in_direction.spec.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.spec.js"],"names":[],"mappings":""}

package/src/core/geom/vec3/v3_displace_in_direction.spec.js

@@ -0,0 +1,80 @@
+import Vector3 from "../Vector3.js";
+import { v3_displace_in_direction } from "./v3_displace_in_direction.js";
+
+test("zero distance", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, 0,
+        0, 0, 0,
+        0, 0, 1
+    )).toBe(true);
+
+});
+
+test("direction magnitude is 0", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, 1,
+        0, 0, 0,
+        0, 0, 0
+    )).toBe(false);
+
+});
+
+test("cardinal directions", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 0, 1
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 2, 0]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 0, -1
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 2, 6]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 1, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, -1, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, -1, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 5, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        1, 0, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([-2, 2, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        -1, 0, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([4, 2, 3]);
+
+});

package/src/core/geom/vec4/v4_length.d.ts

@@ -0,0 +1,10 @@
+/**
+ *
+ * @param {number} x
+ * @param {number} y
+ * @param {number} z
+ * @param {number} w
+ * @returns {number}
+ */
+export function v4_length(x: number, y: number, z: number, w: number): number;
+//# sourceMappingURL=v4_length.d.ts.map

package/src/core/geom/vec4/v4_length.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v4_length.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec4/v4_length.js"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,6BANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,MAAM,CAIlB"}

package/src/core/geom/vec4/v4_length.js

@@ -0,0 +1,13 @@
+import { v4_length_sqr } from "./v4_length_sqr.js";
+
+/**
+ *
+ * @param {number} x
+ * @param {number} y
+ * @param {number} z
+ * @param {number} w
+ * @returns {number}
+ */
+export function v4_length(x, y, z, w) {
+    return Math.sqrt(v4_length_sqr(x, y, z, w));
+}

package/src/core/geom/vec4/v4_length_sqr.js

@@ -8,4 +8,5 @@
  */
 export function v4_length_sqr(x, y, z, w) {
     return x * x + y * y + z * z + w * w;
-}
+}
+

package/src/core/graph/eigen/matrix_eigenvalues_in_place.spec.js

@@ -10,7 +10,22 @@ function roughly_equal_array_vectors(actual, expected) {
     }
 }
 
-test("sample", () => {
+test("2x2 sample A", () => {
+
+
+    const t = new SquareMatrix(2, BinaryDataType.Float32);
+    t.fromArray([
+        1,-2,
+        5,-7
+    ]);
+
+    const eigen = matrix_eigenvalues_in_place(t);
+
+
+    roughly_equal_array_vectors(eigen, [-5.44949, -0.55051]);
+});
+
+test("3x3 sample A", () => {
 
 
     const t = new SquareMatrix(3, BinaryDataType.Float32);
@@ -25,3 +40,18 @@ test("sample", () => {
 
     roughly_equal_array_vectors(eigen, [9.62192181, -9.12030391, 3.4983821]);
 });
+test("3x3 sample B", () => {
+
+
+    const t = new SquareMatrix(3, BinaryDataType.Float32);
+    t.fromArray([
+        1, 2, 3,
+        4, 5, 6,
+        7, 8, 9
+    ]);
+
+    const eigen = matrix_eigenvalues_in_place(t);
+
+
+    roughly_equal_array_vectors(eigen, [16.1168, -1.11684, 0]);
+});

package/src/core/math/FLT_EPSILON_64.d.ts

@@ -0,0 +1,6 @@
+/**
+ * Machine epsilon, 2^(-53)
+ * @type {number}
+ */
+export const FLT_EPSILON_64: number;
+//# sourceMappingURL=FLT_EPSILON_64.d.ts.map

package/src/core/math/FLT_EPSILON_64.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FLT_EPSILON_64.d.ts","sourceRoot":"","sources":["../../../../src/core/math/FLT_EPSILON_64.js"],"names":[],"mappings":"AAAA;;;GAGG;AACH,6BAFU,MAAM,CAEqC"}

package/src/core/math/FLT_EPSILON_64.js

@@ -0,0 +1,5 @@
+/**
+ * Machine epsilon, 2^(-53)
+ * @type {number}
+ */
+export const FLT_EPSILON_64 = 1.1102230246251565e-16;

package/src/core/math/gaussian.d.ts

@@ -1,7 +1,7 @@
 /**
  *
- * @param {number} sigma
- * @param {number} v
+ * @param {number} sigma distance from mean
+ * @param {number} v Variance
  * @returns {number}
  */
 export function gaussian(sigma: number, v: number): number;

package/src/core/math/gaussian.js

@@ -1,7 +1,7 @@
 /**
  *
- * @param {number} sigma
- * @param {number} v
+ * @param {number} sigma distance from mean
+ * @param {number} v Variance
  * @returns {number}
  */
 export function gaussian(sigma, v) {

package/src/core/math/linalg/README.md

@@ -0,0 +1 @@
+Linear Argebra package

package/src/core/math/linalg/lu_factor_linear_system.d.ts

@@ -0,0 +1,13 @@
+/**
+ * Factor linear equations Ax = b using LU decomposition A = LU where L is
+ * lower triangular matrix and U is upper triangular matrix. Input is
+ * square matrix A, integer array of pivot indices index[0->n-1], and size
+ * of square matrix n. Output factorization LU is in matrix A. If error is
+ * found, method returns 0.
+ * @param {number[]} A Square matrix, will be modified
+ * @param {number[]|Uint32Array} index will be written to
+ * @param {number} size Size of the matrix ( 4x4 matrix has size 4 )
+ * @returns {boolean}
+ */
+export function lu_factor_linear_system(A: number[], index: number[] | Uint32Array, size: number): boolean;
+//# sourceMappingURL=lu_factor_linear_system.d.ts.map

package/src/core/math/linalg/lu_factor_linear_system.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"lu_factor_linear_system.d.ts","sourceRoot":"","sources":["../../../../../src/core/math/linalg/lu_factor_linear_system.js"],"names":[],"mappings":"AAOA;;;;;;;;;;GAUG;AACH,2CALW,MAAM,EAAE,SACR,MAAM,EAAE,GAAC,WAAW,QACpB,MAAM,GACJ,OAAO,CA+FnB"}