@woosh/meep-engine 2.110.8 → 2.110.10

package/src/core/geom/3d/tetrahedra/delaunay/compute_delaunay_tetrahedral_mesh.js

@@ -1,15 +1,15 @@
 //
 
+import { assert } from "../../../../assert.js";
+import { BitSet } from "../../../../binary/BitSet.js";
 import { compute_bounding_simplex_3d } from "../compute_bounding_simplex_3d.js";
 import { TetrahedralMesh } from "../TetrahedralMesh.js";
+import { Cavity } from "./Cavity.js";
+import { fill_in_a_cavity } from "./fill_in_a_cavity.js";
+import { tetrahedral_mesh_compute_cavity } from "./tetrahedral_mesh_compute_cavity.js";
 import {
     tetrahedral_mesh_walk_towards_containing_tetrahedron
 } from "./tetrahedral_mesh_walk_towards_containing_tetrahedron.js";
-import { Cavity } from "./Cavity.js";
-import { tetrahedral_mesh_compute_cavity } from "./tetrahedral_mesh_compute_cavity.js";
-import { fill_in_a_cavity } from "./fill_in_a_cavity.js";
-import { BitSet } from "../../../../binary/BitSet.js";
-import { assert } from "../../../../assert.js";
 
 
 /**
@@ -18,7 +18,7 @@ import { assert } from "../../../../assert.js";
  * @param {TetrahedralMesh} mesh
  * @param {number[]} input serialized set of 3d points (x,y,z)
  * @param {number} [n] number of points in the input
- * @returns {boolean}
+ * @returns {boolean} whether mesh was created or not. Returns false on failure
  */
 export function compute_delaunay_tetrahedral_mesh(mesh, input, n = input.length / 3) {
     assert.notNull(mesh, 'mesh');
@@ -43,7 +43,7 @@ export function compute_delaunay_tetrahedral_mesh(mesh, input, n = input.length
 
     // start with a fairly low capacity for best-case scenario, rely on internal resizing after that
 
-    mesh.growCapacity(n + 1);
+    mesh.ensureCapacity(n + 1);
 
     const CAVITY = new Cavity();
 

package/src/core/geom/3d/tetrahedra/delaunay/tetrahedral_mesh_build_from_grid.spec.d.ts

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

package/src/core/geom/3d/tetrahedra/delaunay/tetrahedral_mesh_build_from_grid.spec.d.ts.map

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

package/src/core/geom/3d/tetrahedra/delaunay/tetrahedral_mesh_build_from_grid.spec.js

@@ -0,0 +1,61 @@
+import { AABB3 } from "../../aabb/AABB3.js";
+import { get_tetrahedron_volume } from "../get_tetrahedron_volume.js";
+import { TetrahedralMesh } from "../TetrahedralMesh.js";
+import { validate_tetrahedral_mesh } from "../validate_tetrahedral_mesh.js";
+import { tetrahedral_mesh_build_from_grid } from "./tetrahedral_mesh_build_from_grid.js";
+
+test("Single cube", () => {
+
+    const mesh = new TetrahedralMesh();
+
+    /**
+     *
+     * @type {number[]}
+     */
+    const points = [];
+
+    tetrahedral_mesh_build_from_grid(
+        mesh,
+        points,
+        new AABB3(0, 0, 0, 1, 1, 1),
+        2, 2, 2
+    );
+
+    // number of points
+    expect(points.length).toBe(8 * 3);
+
+    expect(mesh.count).toBeGreaterThanOrEqual(5);
+    expect(mesh.count).toBeLessThanOrEqual(6);
+
+    // mesh should only reference created points
+    const tetrahedrons = mesh.getLive();
+
+    for (let i = 0; i < tetrahedrons.length; i++) {
+        const tet_index = tetrahedrons[i];
+
+        for (let vertex_index = 0; vertex_index < 4; vertex_index++) {
+
+            const vert = mesh.getVertexIndex(tet_index, vertex_index);
+
+            expect(vert).toBeGreaterThanOrEqual(0);
+            expect(vert).toBeLessThan(8);
+        }
+
+        // tetrahedrons should have positive volume
+        const volume = get_tetrahedron_volume(mesh, points, tet_index);
+
+        expect(volume).toBeGreaterThan(0);
+    }
+
+    // do full validation
+    const problems = [];
+
+    const is_valid = validate_tetrahedral_mesh(mesh, points, (problem) => problems.push(problem));
+
+    if (!is_valid) {
+        console.error(problems.join('\n'));
+    }
+
+    expect(is_valid).toBe(true);
+
+});

package/src/core/geom/3d/tetrahedra/get_tetrahedron_volume.d.ts

@@ -0,0 +1,9 @@
+/**
+ *
+ * @param {TetrahedralMesh} mesh
+ * @param {ArrayLike<number>|number[]} points
+ * @param {number} tet_index
+ * @return {number}
+ */
+export function get_tetrahedron_volume(mesh: TetrahedralMesh, points: ArrayLike<number> | number[], tet_index: number): number;
+//# sourceMappingURL=get_tetrahedron_volume.d.ts.map

package/src/core/geom/3d/tetrahedra/get_tetrahedron_volume.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"get_tetrahedron_volume.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/tetrahedra/get_tetrahedron_volume.js"],"names":[],"mappings":"AAEA;;;;;;GAMG;AACH,sEAJW,UAAU,MAAM,CAAC,GAAC,MAAM,EAAE,aAC1B,MAAM,GACL,MAAM,CAUjB"}

package/src/core/geom/3d/tetrahedra/get_tetrahedron_volume.js

@@ -0,0 +1,18 @@
+import { tetrahedron_compute_signed_volume } from "./tetrahedron_compute_signed_volume.js";
+
+/**
+ *
+ * @param {TetrahedralMesh} mesh
+ * @param {ArrayLike<number>|number[]} points
+ * @param {number} tet_index
+ * @return {number}
+ */
+export function get_tetrahedron_volume(mesh, points, tet_index) {
+
+    const a = mesh.getVertexIndex(tet_index, 0);
+    const b = mesh.getVertexIndex(tet_index, 1);
+    const c = mesh.getVertexIndex(tet_index, 2);
+    const d = mesh.getVertexIndex(tet_index, 3);
+
+    return tetrahedron_compute_signed_volume(points, a, b, c, d);
+}

package/src/core/geom/3d/tetrahedra/is_tetrahedron_degenerate.d.ts

@@ -0,0 +1,9 @@
+/**
+ * Vertices must be unique inside a single tetrahedron, that is - no repetitions are allowed
+ * If vertices are repeated - volume will be 0 and tetrahedron will be degenerate
+ * @param {TetrahedralMesh} mesh
+ * @param {number} tet
+ * @returns {boolean}
+ */
+export function is_tetrahedron_degenerate(mesh: TetrahedralMesh, tet: number): boolean;
+//# sourceMappingURL=is_tetrahedron_degenerate.d.ts.map

package/src/core/geom/3d/tetrahedra/is_tetrahedron_degenerate.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"is_tetrahedron_degenerate.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/tetrahedra/is_tetrahedron_degenerate.js"],"names":[],"mappings":"AAAA;;;;;;GAMG;AACH,sEAHW,MAAM,GACJ,OAAO,CAenB"}

package/src/core/geom/3d/tetrahedra/is_tetrahedron_degenerate.js

@@ -0,0 +1,21 @@
+/**
+ * Vertices must be unique inside a single tetrahedron, that is - no repetitions are allowed
+ * If vertices are repeated - volume will be 0 and tetrahedron will be degenerate
+ * @param {TetrahedralMesh} mesh
+ * @param {number} tet
+ * @returns {boolean}
+ */
+export function is_tetrahedron_degenerate(mesh, tet) {
+
+    for (let i0 = 0; i0 < 4; i0++) {
+        const v0 = mesh.getVertexIndex(tet, i0);
+        for (let i1 = i0 + 1; i1 < 4; i1++) {
+            if (v0 === mesh.getVertexIndex(tet, i1)) {
+                // identical vertices
+                return true;
+            }
+        }
+    }
+
+    return false;
+}

package/src/core/geom/3d/tetrahedra/validate_tetrahedral_mesh.d.ts

@@ -1,11 +1,3 @@
-/**
- *
- * @param {TetrahedralMesh} mesh
- * @param {ArrayLike<number>|number[]} points
- * @param {number} tet
- * @return {number}
- */
-export function get_tetrahedron_volume(mesh: TetrahedralMesh, points: ArrayLike<number> | number[], tet: number): number;
 /**
  *
  * @param {TetrahedralMesh} mesh
@@ -25,10 +17,10 @@ export function validate_tetrahedron_neighbourhood(mesh: TetrahedralMesh, tet: n
 /**
  *
  * @param {TetrahedralMesh} mesh
- * @param {ArrayLike<number>} points
+ * @param {ArrayLike<number>|number[]|Float32Array} points
  * @param {function(problem:string):*} [consumer]
  * @returns {boolean}
  */
-export function validate_tetrahedral_mesh(mesh: TetrahedralMesh, points: ArrayLike<number>, consumer?: typeof noop): boolean;
+export function validate_tetrahedral_mesh(mesh: TetrahedralMesh, points: ArrayLike<number> | number[] | Float32Array, consumer?: typeof noop): boolean;
 import { noop } from "../../../function/noop.js";
 //# sourceMappingURL=validate_tetrahedral_mesh.d.ts.map

package/src/core/geom/3d/tetrahedra/validate_tetrahedral_mesh.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"validate_tetrahedral_mesh.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/tetrahedra/validate_tetrahedral_mesh.js"],"names":[],"mappings":"AAIA;;;;;;GAMG;AACH,sEAJW,UAAU,MAAM,CAAC,GAAC,MAAM,EAAE,OAC1B,MAAM,GACL,MAAM,CAUjB;AAED;;;;;;;GAOG;AACH,+DALW,MAAM,gBACN,MAAM,mBACG,MAAM,WACd,OAAO,CA6DlB;AAGD;;;;;GAKG;AACH,+EAHW,MAAM,mCAehB;AAuBD;;;;;;GAMG;AACH,yEAJW,UAAU,MAAM,CAAC,2BAEf,OAAO,CA4BnB;qBArKoB,2BAA2B"}
+{"version":3,"file":"validate_tetrahedral_mesh.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/tetrahedra/validate_tetrahedral_mesh.js"],"names":[],"mappings":"AAKA;;;;;;;GAOG;AACH,+DALW,MAAM,gBACN,MAAM,mBACG,MAAM,WACd,OAAO,CA6DlB;AAGD;;;;;GAKG;AACH,+EAHW,MAAM,mCAehB;AAED;;;;;;GAMG;AACH,yEAJW,UAAU,MAAM,CAAC,GAAC,MAAM,EAAE,GAAC,YAAY,2BAErC,OAAO,CA4BnB;qBAhIoB,2BAA2B"}

package/src/core/geom/3d/tetrahedra/validate_tetrahedral_mesh.js

@@ -1,23 +1,7 @@
 import { noop } from "../../../function/noop.js";
+import { get_tetrahedron_volume } from "./get_tetrahedron_volume.js";
+import { is_tetrahedron_degenerate } from "./is_tetrahedron_degenerate.js";
 import { INVALID_NEIGHBOUR } from "./TetrahedralMesh.js";
-import { tetrahedron_compute_signed_volume } from "./tetrahedron_compute_signed_volume.js";
-
-/**
- *
- * @param {TetrahedralMesh} mesh
- * @param {ArrayLike<number>|number[]} points
- * @param {number} tet
- * @return {number}
- */
-export function get_tetrahedron_volume(mesh, points, tet) {
-
-    const a = mesh.getVertexIndex(tet, 0);
-    const b = mesh.getVertexIndex(tet, 1);
-    const c = mesh.getVertexIndex(tet, 2);
-    const d = mesh.getVertexIndex(tet, 3);
-
-    return tetrahedron_compute_signed_volume(points, a, b, c, d);
-}
 
 /**
  *
@@ -112,28 +96,7 @@ export function validate_tetrahedron_neighbourhood(mesh, tet, consumer = noop) {
 /**
  *
  * @param {TetrahedralMesh} mesh
- * @param {number} tet
- * @returns {boolean}
- */
-function is_tetrahedron_degenerate(mesh, tet) {
-
-    for (let i0 = 0; i0 < 4; i0++) {
-        const v0 = mesh.getVertexIndex(tet, i0);
-        for (let i1 = i0 + 1; i1 < 4; i1++) {
-            if (v0 === mesh.getVertexIndex(tet, i1)) {
-                // identical vertices
-                return true;
-            }
-        }
-    }
-
-    return false;
-}
-
-/**
- *
- * @param {TetrahedralMesh} mesh
- * @param {ArrayLike<number>} points
+ * @param {ArrayLike<number>|number[]|Float32Array} points
  * @param {function(problem:string):*} [consumer]
  * @returns {boolean}
  */

package/src/core/geom/3d/v3_compute_triangle_normal.js

@@ -54,4 +54,4 @@ export function v3_compute_triangle_normal(
     result[result_offset] = x;
     result[result_offset + 1] = y;
     result[result_offset + 2] = z;
-}
+}

package/src/core/geom/3d/v3_negate_array.d.ts

@@ -0,0 +1,9 @@
+/**
+ *
+ * @param {number[]|Float32Array} result
+ * @param {number} result_offset
+ * @param {number[]|Float32Array} input
+ * @param {number} input_offset
+ */
+export function v3_negate_array(result: number[] | Float32Array, result_offset: number, input: number[] | Float32Array, input_offset: number): void;
+//# sourceMappingURL=v3_negate_array.d.ts.map

package/src/core/geom/3d/v3_negate_array.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v3_negate_array.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/3d/v3_negate_array.js"],"names":[],"mappings":"AAAA;;;;;;GAMG;AACH,wCALW,MAAM,EAAE,GAAC,YAAY,iBACrB,MAAM,SACN,MAAM,EAAE,GAAC,YAAY,gBACrB,MAAM,QAUhB"}

package/src/core/geom/3d/v3_negate_array.js

@@ -0,0 +1,16 @@
+/**
+ *
+ * @param {number[]|Float32Array} result
+ * @param {number} result_offset
+ * @param {number[]|Float32Array} input
+ * @param {number} input_offset
+ */
+export function v3_negate_array(result, result_offset, input, input_offset) {
+    const x = input[input_offset];
+    const y = input[input_offset + 1];
+    const z = input[input_offset + 2];
+
+    result[result_offset] = -x;
+    result[result_offset + 1] = -y;
+    result[result_offset + 2] = -z;
+}

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

@@ -1 +1 @@
-{"version":3,"file":"Quaternion.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Quaternion.js"],"names":[],"mappings":";AAsBA;IA+2CI;;;;OAIG;IACH,kCAFa,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,0BALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,6BALW,UAAU,QACV,UAAU,MACV,UAAU,aACV,MAAM,QAkBhB;IAj6CD;;;;;;;OAOG;IACH,gBANW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,EA+BhB;IA1BG;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,oBAFU,OAAO,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAEnD;IAqBjC,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IA+BD;;;;;;;;OAQG;IACH,kBAPW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAiDhB;IAED;;;;;OAKG;IACH,uBAJW,OAAO,QACP,OAAO,QAMjB;IAED;;;;OAIG;IACH,WAHW,UAAU,GACT,MAAM,CAQjB;IAED;;;OAGG;IACH,mBAFW,UAAU,QAKpB;IAED;;OAEG;IACH,eAqBC;IAED;;;;OAIG;IACH,eAHW,UAAU,GACT,MAAM,CAoBjB;IAED;;;;OAIG;IACH,oBAHW,OAAO,SACP,MAAM,QAIhB;IAED;;;;;;OAMG;IACH,mBALW,MAAM,MACN,MAAM,MACN,MAAM,SACN,MAAM,QAwBhB;IAED;;;;;OAKG;IACH,2BAJW,OAAO,SACP,UAAU,SACV,UAAU,QA+BpB;IAED;;;;OAIG;IACH,wBAHW,OAAO,GACL,MAAM,CAYlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAkBlB;IAED,kBAWC;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACL,UAAU,CASrB;IAED;;OAEG;IACH,gBAFW,UAAU,QAIpB;IAED;;;;OAIG;IACH,2BAHW,UAAU,UACV,UAAU,QAcpB;IAED;;;;;;;;;;;OAWG;IACH,yBAVW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACL,UAAU,CAYrB;IAED;;;OAGG;IACH,UAFY,MAAM,CASjB;IAED;;;;OAIG;IACH,qBAHW,UAAU,aACV,MAAM,QAIhB;IAED;;;;;OAKG;IACH,eAJW,OAAO,UACP,OAAO,OACP,OAAO,QAUjB;IAED;;;OAGG;IACH,wBAFsB,MAAM,QAM3B;IAED;;;;;;;OAOG;IACH,kBANW,MAAM,KACN,MAAM,KACN,MAAM,mBAEJ,UAAU,CAmCtB;IAED;;;;OAIG;IACH,yBAFW,OAAO,QA4BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA2BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA0BjB;IAED;;;;;;OAMG;IACH,mBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAIhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;OAIG;IACH,yBAHW,OAAO,MACP,OAAO,QA+CjB;IAED;;;;;OAKG;IACH,sBAHW,OAAO,MACP,OAAO,QAiEjB;IAED;;;OAGG;IACH,wCASC;IAED;;;;;;;;;;;;;;OAcG;IACH,6BAXW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,GACJ,UAAU,CA4EtB;IAED;;;;OAIG;IACH,YAHW,UAAU,KACV,MAAM,QAMhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,UACV,UAAU,KACV,MAAM,QAgBhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,MACV,UAAU,KACV,MAAM,QAiDhB;IAGD;;;;OAIG;IACH,aAHW,UAAU,KACV,MAAM,QAIhB;IAED;;;OAGG;IACH,iCAIC;IAED;;;;OAIG;IACH,YAHW,UAAU,GACR,UAAU,CAItB;IAED;;;OAGG;IACH,SAFa,UAAU,CAQtB;IAED;;;;;;;OAOG;IACH,OANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CA+BtB;IAED;;;OAGG;IACH,aAFa,UAAU,CAItB;IAED;;;;;MAOC;IAED,yBAEC;IAED;;;OAGG;IACH,2CAKC;IAED;;;OAGG;IACH,6CAOC;IAED;;;OAGG;IACH,kDAKC;IAED;;;OAGG;IACH,oDAOC;IAED;;;;OAIG;IACH,wBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,kBAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAShB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAOhB;IAED,oBAEC;IAED;;;;OAIG;IACH,cAHW,UAAU,GACR,OAAO,CAQnB;IAED;;;OAGG;IACH,QAFa,MAAM,CAQlB;IAED;;;;;OAKG;IACH,qBAJW,UAAU,cACV,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;;;;OASG;IACH,uBAKC;IAED;;;;OAIG;IACH,sBAHsB,MAAM,GAChB,UAAU,CAkBrB;IAED,mBAEC;IAwDL,mBAjKe,MAAM,EAAE,WACR,MAAM,UAgKS;IAC9B,iBApJe,MAAM,EAAE,WACR,MAAM,UAmJO;IA71CxB;;OAEG;IACH,sDAOC;CAg1CJ;;kBAOS,UAAU;;mBA/7CD,4BAA4B;oBAW3B,cAAc"}
+{"version":3,"file":"Quaternion.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Quaternion.js"],"names":[],"mappings":";AAsBA;IAg0CI;;;;OAIG;IACH,kCAFa,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,0BALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,6BALW,UAAU,QACV,UAAU,MACV,UAAU,aACV,MAAM,QAkBhB;IAl3CD;;;;;;;OAOG;IACH,gBANW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,EA+BhB;IA1BG;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,oBAFU,OAAO,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAEnD;IAqBjC,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IA+BD;;;;;;;;OAQG;IACH,kBAPW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAiDhB;IAED;;;;;OAKG;IACH,uBAJW,OAAO,QACP,OAAO,QAMjB;IAED;;;;OAIG;IACH,WAHW,UAAU,GACT,MAAM,CAQjB;IAED;;;OAGG;IACH,mBAFW,UAAU,QAKpB;IAED;;OAEG;IACH,eAqBC;IAED;;;;OAIG;IACH,eAHW,UAAU,GACT,MAAM,CAoBjB;IAED;;;;OAIG;IACH,oBAHW,OAAO,SACP,MAAM,QAIhB;IAED;;;;;;OAMG;IACH,mBALW,MAAM,MACN,MAAM,MACN,MAAM,SACN,MAAM,QA6BhB;IAED;;;;;OAKG;IACH,2BAJW,OAAO,SACP,UAAU,SACV,UAAU,QA+BpB;IAED;;;;OAIG;IACH,wBAHW,OAAO,GACL,MAAM,CAYlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAkBlB;IAED,kBAWC;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACL,UAAU,CASrB;IAED;;OAEG;IACH,gBAFW,UAAU,QAIpB;IAED;;;;OAIG;IACH,2BAHW,UAAU,UACV,UAAU,QAcpB;IAED;;;;;;;;;;;OAWG;IACH,yBAVW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACL,UAAU,CAYrB;IAED;;;OAGG;IACH,UAFY,MAAM,CASjB;IAED;;;;OAIG;IACH,qBAHW,UAAU,aACV,MAAM,QAIhB;IAED;;;;;OAKG;IACH,eAJW,OAAO,UACP,OAAO,OACP,OAAO,QAUjB;IAED;;;OAGG;IACH,wBAFsB,MAAM,QAM3B;IAED;;;;;;;OAOG;IACH,kBANW,MAAM,KACN,MAAM,KACN,MAAM,mBAEJ,UAAU,CAmCtB;IAED;;;;OAIG;IACH,yBAFW,OAAO,QA4BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA2BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA0BjB;IAED;;;;;;OAMG;IACH,mBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAIhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;OAKG;IACH,sBAHW,OAAO,MACP,OAAO,QAiEjB;IAED;;;OAGG;IACH,wCASC;IAED;;;;;;;;;;;;;;OAcG;IACH,6BAXW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,GACJ,UAAU,CA4EtB;IAED;;;;OAIG;IACH,YAHW,UAAU,KACV,MAAM,QAMhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,UACV,UAAU,KACV,MAAM,QAgBhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,MACV,UAAU,KACV,MAAM,QAiDhB;IAGD;;;;OAIG;IACH,aAHW,UAAU,KACV,MAAM,QAIhB;IAED;;;OAGG;IACH,iCAIC;IAED;;;;OAIG;IACH,YAHW,UAAU,GACR,UAAU,CAItB;IAED;;;OAGG;IACH,SAFa,UAAU,CAQtB;IAED;;;;;;;OAOG;IACH,OANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CA+BtB;IAED;;;OAGG;IACH,aAFa,UAAU,CAItB;IAED;;;;;MAOC;IAED,yBAEC;IAED;;;OAGG;IACH,2CAKC;IAED;;;OAGG;IACH,6CAOC;IAED;;;OAGG;IACH,kDAKC;IAED;;;OAGG;IACH,oDAOC;IAED;;;;OAIG;IACH,wBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,kBAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAShB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAOhB;IAED,oBAEC;IAED;;;;OAIG;IACH,cAHW,UAAU,GACR,OAAO,CAQnB;IAED;;;OAGG;IACH,QAFa,MAAM,CAQlB;IAED;;;;;OAKG;IACH,qBAJW,UAAU,cACV,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;;;;OASG;IACH,uBAKC;IAED;;;;OAIG;IACH,sBAHsB,MAAM,GAChB,UAAU,CAkBrB;IAED,mBAEC;IAwDL,mBAjKe,MAAM,EAAE,WACR,MAAM,UAgKS;IAC9B,iBApJe,MAAM,EAAE,WACR,MAAM,UAmJO;IA9yCxB;;OAEG;IACH,sDAOC;CAiyCJ;;kBAOS,UAAU;;mBAh5CD,4BAA4B;oBAW3B,cAAc"}

package/src/core/geom/Quaternion.js

@@ -13,9 +13,9 @@ import { v4_length } from "./vec4/v4_length.js";
 import Vector3 from "./Vector3.js";
 
 
-const forward = new Vector3();
-const up = new Vector3();
-const right = new Vector3();
+const scratch_v3_a = new Vector3();
+const scratch_v3_b = new Vector3();
+const scratch_v3_c = new Vector3();
 
 const sin = Math.sin;
 const cos = Math.cos;
@@ -116,12 +116,12 @@ class Quaternion {
      * @param {number} uz up vector
      */
     _lookRotation(fx, fy, fz, ux, uy, uz) {
-        forward.set(fx, fy, fz);
-        forward.normalize();
+        scratch_v3_a.set(fx, fy, fz);
+        scratch_v3_a.normalize();
 
-        right._crossVectors(ux, uy, uz, forward.x, forward.y, forward.z);
+        scratch_v3_c._crossVectors(ux, uy, uz, scratch_v3_a.x, scratch_v3_a.y, scratch_v3_a.z);
 
-        if (right.lengthSq() === 0) {
+        if (scratch_v3_c.lengthSq() === 0) {
             // up and forward are parallel
 
             /*
@@ -130,32 +130,32 @@ class Quaternion {
              */
 
             if (Math.abs(uz) === 1) {
-                forward.x += 0.001;
+                scratch_v3_a.x += 0.001;
             } else {
-                forward.z += 0.001;
+                scratch_v3_a.z += 0.001;
             }
 
-            forward.normalize();
-            right._crossVectors(ux, uy, uz, forward.x, forward.y, forward.z);
+            scratch_v3_a.normalize();
+            scratch_v3_c._crossVectors(ux, uy, uz, scratch_v3_a.x, scratch_v3_a.y, scratch_v3_a.z);
 
         }
 
-        right.normalize();
+        scratch_v3_c.normalize();
 
-        up.crossVectors(forward, right);
+        scratch_v3_b.crossVectors(scratch_v3_a, scratch_v3_c);
 
         // construct partial transform matrix
-        const m00 = right.x;
-        const m10 = right.y;
-        const m20 = right.z;
+        const m00 = scratch_v3_c.x;
+        const m10 = scratch_v3_c.y;
+        const m20 = scratch_v3_c.z;
 
-        const m01 = up.x;
-        const m11 = up.y;
-        const m21 = up.z;
+        const m01 = scratch_v3_b.x;
+        const m11 = scratch_v3_b.y;
+        const m21 = scratch_v3_b.z;
 
-        const m02 = forward.x;
-        const m12 = forward.y;
-        const m22 = forward.z;
+        const m02 = scratch_v3_a.x;
+        const m12 = scratch_v3_a.y;
+        const m22 = scratch_v3_a.z;
 
         this.__setFromRotationMatrix(
             m00, m01, m02,
@@ -189,7 +189,7 @@ class Quaternion {
     }
 
     /**
-     *
+     * Makes this quaternion into an inverse of the other
      * @param {Quaternion} other
      */
     copyInverse(other) {
@@ -209,7 +209,7 @@ class Quaternion {
         const length_sqr = x * x + y * y + z * z + w * w;
 
         if (length_sqr === 0) {
-            this.set(0, 0, 0, 0);
+            this.set(0, 0, 0, 1);
             return;
         }
 
@@ -265,6 +265,11 @@ class Quaternion {
      * @param {number} angle
      */
     _fromAxisAngle(ax, ay, az, angle) {
+        assert.isNumber(ax, 'ax');
+        assert.isNumber(ay, 'ay');
+        assert.isNumber(az, 'az');
+        assert.isNumber(angle, 'angle');
+
         const halfAngle = angle * 0.5;
 
         const sinA2 = sin(halfAngle);
@@ -303,7 +308,7 @@ class Quaternion {
         // perform projection of rotation onto axis
         const d = v3_dot(x, y, z, axis.x, axis.y, axis.z);
 
-        const mag2 = (axis.x * axis.x + axis.y * axis.y + axis.z * axis.z);
+        const mag2 = axis.x * axis.x + axis.y * axis.y + axis.z * axis.z;
 
         const m = d / mag2;
 
@@ -818,58 +823,6 @@ class Quaternion {
         this.set(_x, _y, _z, _w);
     }
 
-    /**
-     * Based on blog post: http://lolengine.net/blog/2013/09/18/beautiful-maths-quaternion-from-vectors
-     * @param {Vector3} from
-     * @param {Vector3} to
-     */
-    fromNonUnitVectors(from, to) {
-        //quat quat::fromtwovectors(vec3 u, vec3 v)
-        // {
-        //     vec3 w = cross(u, v);
-        //     quat q = quat(dot(u, v), w.x, w.y, w.z);
-        //     q.w += length(q);
-        //     return normalize(q);
-        // }
-
-        const ax = from.x;
-        const ay = from.y;
-        const az = from.z;
-
-        const bx = to.x;
-        const by = to.y;
-        const bz = to.z;
-
-        //compute cross product
-        const crossX = ay * bz - az * by;
-        const crossY = az * bx - ax * bz;
-        const crossZ = ax * by - ay * bx;
-
-        const uv_dot = v3_dot(ax, ay, az, bx, by, bz);
-
-
-        const qX = crossX;
-        const qY = crossY;
-        const qZ = crossZ;
-        const qW = uv_dot;
-
-        const p0 = qX * qX + qY * qY + qZ * qZ;
-
-        const l0 = Math.sqrt(p0 + qW * qW);
-
-        const qW_1 = qW + l0;
-
-        //normalize result
-        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;
-
-        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

package/src/core/geom/Vector3.d.ts

@@ -112,6 +112,10 @@ export default class Vector3 implements Vector3Like {
 
     writeToArray(array: number[] | ArrayLike<number>, offset?: number): void
 
+    toArray(array: number[] | ArrayLike<number>, offset?: number): void
+
+    fromArray(array: number[] | ArrayLike<number>, offset?: number): void
+
     asArray(): number[]
 
     static distance(a: Vector3Like, b: Vector3Like): number

package/src/core/geom/Vector3.js

@@ -505,10 +505,10 @@ class Vector3 {
 
         // calculate quat * vector
 
-        var ix = qw * x + qy * z - qz * y;
-        var iy = qw * y + qz * x - qx * z;
-        var iz = qw * z + qx * y - qy * x;
-        var iw = -qx * x - qy * y - qz * z;
+        const ix = qw * x + qy * z - qz * y;
+        const iy = qw * y + qz * x - qx * z;
+        const iz = qw * z + qx * y - qy * x;
+        const iw = -qx * x - qy * y - qz * z;
 
         // calculate result * inverse quat
 

package/src/core/geom/mat3/m3_cm_invert.d.ts

@@ -0,0 +1,7 @@
+/**
+ *
+ * @param {Float32Array} out
+ * @param {Float32Array} a
+ */
+export function m3_cm_invert(out: Float32Array, a: Float32Array): boolean;
+//# sourceMappingURL=m3_cm_invert.d.ts.map

package/src/core/geom/mat3/m3_cm_invert.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"m3_cm_invert.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/mat3/m3_cm_invert.js"],"names":[],"mappings":"AAAA;;;;GAIG;AACH,kCAHW,YAAY,KACZ,YAAY,WAyCtB"}

package/src/core/geom/mat3/m3_cm_invert.js

@@ -0,0 +1,45 @@
+/**
+ *
+ * @param {Float32Array} out
+ * @param {Float32Array} a
+ */
+export function m3_cm_invert(out, a) {
+    // adapted from gl-matrix
+
+    const a00 = a[0],
+        a01 = a[1],
+        a02 = a[2];
+
+    const a10 = a[3],
+        a11 = a[4],
+        a12 = a[5];
+
+    const a20 = a[6],
+        a21 = a[7],
+        a22 = a[8];
+
+    const b01 = a22 * a11 - a12 * a21;
+    const b11 = -a22 * a10 + a12 * a20;
+    const b21 = a21 * a10 - a11 * a20;
+
+    // Calculate the determinant
+    let det = a00 * b01 + a01 * b11 + a02 * b21;
+
+    if (!det) {
+        return false;
+    }
+
+    det = 1.0 / det;
+
+    out[0] = b01 * det;
+    out[1] = (-a22 * a01 + a02 * a21) * det;
+    out[2] = (a12 * a01 - a02 * a11) * det;
+    out[3] = b11 * det;
+    out[4] = (a22 * a00 - a02 * a20) * det;
+    out[5] = (-a12 * a00 + a02 * a10) * det;
+    out[6] = b21 * det;
+    out[7] = (-a21 * a00 + a01 * a20) * det;
+    out[8] = (a11 * a00 - a01 * a10) * det;
+
+    return true;
+}

package/src/core/geom/vec2/v2_matrix3_cm_multiply.d.ts

@@ -0,0 +1,10 @@
+/**
+ * Column-major ( OpenGL style) matrix multiplication
+ * @param {Float32Array|number[]} out
+ * @param {number} out_offset
+ * @param {number} x
+ * @param {number} y
+ * @param {Float32Array|number[]} m3
+ */
+export function v2_matrix3_cm_multiply(out: Float32Array | number[], out_offset: number, x: number, y: number, m3: Float32Array | number[]): void;
+//# sourceMappingURL=v2_matrix3_cm_multiply.d.ts.map

package/src/core/geom/vec2/v2_matrix3_cm_multiply.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v2_matrix3_cm_multiply.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec2/v2_matrix3_cm_multiply.js"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AACH,4CANW,YAAY,GAAC,MAAM,EAAE,cACrB,MAAM,KACN,MAAM,KACN,MAAM,MACN,YAAY,GAAC,MAAM,EAAE,QAa/B"}

package/src/core/geom/vec2/v2_matrix3_cm_multiply.js

@@ -0,0 +1,20 @@
+/**
+ * Column-major ( OpenGL style) matrix multiplication
+ * @param {Float32Array|number[]} out
+ * @param {number} out_offset
+ * @param {number} x
+ * @param {number} y
+ * @param {Float32Array|number[]} m3
+ */
+export function v2_matrix3_cm_multiply(
+    out, out_offset,
+    x, y,
+    m3
+) {
+
+    const tx = m3[0] * x + m3[1] * y + m3[2];
+    const ty = m3[3] * x + m3[4] * y + m3[5];
+
+    out[out_offset] = tx;
+    out[out_offset + 1] = ty;
+}

package/src/core/geom/vec2/v2_matrix3_rm_multiply.d.ts

@@ -0,0 +1,10 @@
+/**
+ * Row-major ( OpenGL style) matrix multiplication
+ * @param {Float32Array|number[]} out
+ * @param {number} out_offset
+ * @param {number} x
+ * @param {number} y
+ * @param {Float32Array|number[]} m3
+ */
+export function v2_matrix3_rm_multiply(out: Float32Array | number[], out_offset: number, x: number, y: number, m3: Float32Array | number[]): void;
+//# sourceMappingURL=v2_matrix3_rm_multiply.d.ts.map

package/src/core/geom/vec2/v2_matrix3_rm_multiply.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v2_matrix3_rm_multiply.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec2/v2_matrix3_rm_multiply.js"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AACH,4CANW,YAAY,GAAC,MAAM,EAAE,cACrB,MAAM,KACN,MAAM,KACN,MAAM,MACN,YAAY,GAAC,MAAM,EAAE,QAa/B"}

package/src/core/geom/vec2/v2_matrix3_rm_multiply.js

@@ -0,0 +1,21 @@
+/**
+ * Row-major ( OpenGL style) matrix multiplication
+ * @param {Float32Array|number[]} out
+ * @param {number} out_offset
+ * @param {number} x
+ * @param {number} y
+ * @param {Float32Array|number[]} m3
+ */
+export function v2_matrix3_rm_multiply(
+    out, out_offset,
+    x, y,
+    m3
+) {
+
+    const tx = m3[0] * x + m3[3] * y + m3[6];
+    const ty = m3[1] * x + m3[4] * y + m3[7];
+
+    out[out_offset] = tx;
+    out[out_offset + 1] = ty;
+}
+