@woosh/meep-engine 2.114.0 → 2.116.0

package/src/core/geom/2d/convex-hull/orientation3_array.js

@@ -1,4 +1,4 @@
-import { orientation3 } from "./orientation3.js";
+import { orient2d } from "robust-predicates";
 
 /**
  *
@@ -18,5 +18,5 @@ export function orientation3_array(points, a, b, c) {
     const cx = points[c * 2];
     const cy = points[c * 2 + 1];
 
-    return orientation3(ax, ay, bx, by, cx, cy);
+    return orient2d(ax, ay, bx, by, cx, cy);
 }

package/src/core/geom/2d/intersect_ray_2d.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"intersect_ray_2d.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/intersect_ray_2d.js"],"names":[],"mappings":"AAEA;;;;;;;;;;;;;GAaG;AACH,sCAZW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,cAClC,MAAM,cACN,MAAM,cACN,MAAM,iBACN,MAAM,iBACN,MAAM,cACN,MAAM,cACN,MAAM,iBACN,MAAM,iBACN,MAAM,GACJ,OAAO,CAkCnB"}
+{"version":3,"file":"intersect_ray_2d.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/intersect_ray_2d.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;GAaG;AACH,sCAZW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,cAClC,MAAM,cACN,MAAM,cACN,MAAM,iBACN,MAAM,iBACN,MAAM,cACN,MAAM,cACN,MAAM,iBACN,MAAM,iBACN,MAAM,GACJ,OAAO,CAoCnB"}

package/src/core/geom/2d/intersect_ray_2d.js

@@ -1,5 +1,3 @@
-import { fabsf } from "../../math/fabsf.js";
-
 /**
  * NOTE: that direction must not be normalized for correct results and instead derived from line length
  * @param {number[]|Float32Array|Float64Array} out
@@ -20,15 +18,17 @@ export function intersect_ray_2d(
     origin_b_x, origin_b_y, direction_b_x, direction_b_y
 ) {
 
-    const d = (direction_a_x * direction_b_y - direction_a_y * direction_b_x);
+    const denominator = (direction_a_x * direction_b_y - direction_a_y * direction_b_x);
 
-    if (fabsf(d) < 0.000000000001) // Parallel lines
+    if (Math.abs(denominator) < 1e-7){
+        // Parallel lines
         return false;
+    }
 
     const oab_y = origin_a_y - origin_b_y;
     const oab_x = origin_a_x - origin_b_x;
 
-    const t = (direction_b_x * oab_y - direction_b_y * oab_x) / d;
+    const t = (direction_b_x * oab_y - direction_b_y * oab_x) / denominator;
 
     if (t < 0.5) // Intersects on the wrong side
         return false;

package/src/core/geom/2d/{compute_triangle_area_2d.d.ts → triangle2d_compute_area.d.ts}

@@ -8,5 +8,5 @@
  * @param {number} y2
  * @returns {number}
  */
-export function compute_triangle_area_2d(x0: number, y0: number, x1: number, y1: number, x2: number, y2: number): number;
-//# sourceMappingURL=compute_triangle_area_2d.d.ts.map
+export function triangle2d_compute_area(x0: number, y0: number, x1: number, y1: number, x2: number, y2: number): number;
+//# sourceMappingURL=triangle2d_compute_area.d.ts.map

package/src/core/geom/2d/triangle2d_compute_area.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"triangle2d_compute_area.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/triangle2d_compute_area.js"],"names":[],"mappings":"AAAA;;;;;;;;;GASG;AACH,4CARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAclB"}

package/src/core/geom/2d/{compute_triangle_area_2d.js → triangle2d_compute_area.js}

@@ -8,7 +8,7 @@
  * @param {number} y2
  * @returns {number}
  */
-export function compute_triangle_area_2d(
+export function triangle2d_compute_area(
     x0, y0,
     x1, y1,
     x2, y2

package/src/core/geom/2d/triangle2d_contains_point.d.ts

@@ -0,0 +1,14 @@
+/**
+ *
+ * @param {number} ax
+ * @param {number} ay
+ * @param {number} bx
+ * @param {number} by
+ * @param {number} cx
+ * @param {number} cy
+ * @param {number} px reference point X
+ * @param {number} py reference point Y
+ * @returns {boolean}
+ */
+export function triangle2d_contains_point(ax: number, ay: number, bx: number, by: number, cx: number, cy: number, px: number, py: number): boolean;
+//# sourceMappingURL=triangle2d_contains_point.d.ts.map

package/src/core/geom/2d/triangle2d_contains_point.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"triangle2d_contains_point.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/triangle2d_contains_point.js"],"names":[],"mappings":"AAIA;;;;;;;;;;;GAWG;AACH,8CAVW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,OAAO,CAwBnB"}

package/src/core/geom/2d/triangle2d_contains_point.js

@@ -0,0 +1,39 @@
+import { triangle2d_get_barycentric } from "./triangle2d_get_barycentric.js";
+
+const scratch = new Float32Array(3);
+
+/**
+ *
+ * @param {number} ax
+ * @param {number} ay
+ * @param {number} bx
+ * @param {number} by
+ * @param {number} cx
+ * @param {number} cy
+ * @param {number} px reference point X
+ * @param {number} py reference point Y
+ * @returns {boolean}
+ */
+export function triangle2d_contains_point(
+    ax, ay,
+    bx, by,
+    cx, cy,
+    px, py
+) {
+    triangle2d_get_barycentric(
+        scratch, 0,
+        ax, ay,
+        bx, by,
+        cx, cy,
+        px, py
+    );
+
+    const u = scratch[0];
+    const v = scratch[1];
+    const w = scratch[2];
+
+    return u >= 0
+        && v >= 0
+        && w >= 0
+        && Math.abs(1 - u - v - w) < 1e-7;
+}

package/src/core/geom/2d/triangle2d_get_barycentric.d.ts

@@ -1,6 +1,6 @@
 /**
  *
- * @param {number[]} output
+ * @param {number[]|TypedArray|Float32Array} output
  * @param {number} output_offset
  * @param {number} ax
  * @param {number} ay
@@ -11,5 +11,5 @@
  * @param {number} px reference point X
  * @param {number} py reference point Y
  */
-export function triangle2d_get_barycentric(output: number[], output_offset: number, ax: number, ay: number, bx: number, by: number, cx: number, cy: number, px: number, py: number): void;
+export function triangle2d_get_barycentric(output: number[] | TypedArray | Float32Array, output_offset: number, ax: number, ay: number, bx: number, by: number, cx: number, cy: number, px: number, py: number): void;
 //# sourceMappingURL=triangle2d_get_barycentric.d.ts.map

package/src/core/geom/2d/triangle2d_get_barycentric.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"triangle2d_get_barycentric.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/triangle2d_get_barycentric.js"],"names":[],"mappings":"AAEA;;;;;;;;;;;;GAYG;AACH,mDAXW,MAAM,EAAE,iBACR,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAkChB"}
+{"version":3,"file":"triangle2d_get_barycentric.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/triangle2d_get_barycentric.js"],"names":[],"mappings":"AAEA;;;;;;;;;;;;GAYG;AACH,mDAXW,MAAM,EAAE,gBAAY,YAAY,iBAChC,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAkChB"}

package/src/core/geom/2d/triangle2d_get_barycentric.js

@@ -2,7 +2,7 @@ import { v2_dot } from "../vec2/v2_dot.js";
 
 /**
  *
- * @param {number[]} output
+ * @param {number[]|TypedArray|Float32Array} output
  * @param {number} output_offset
  * @param {number} ax
  * @param {number} ay

package/src/core/geom/3d/topology/expandConnectivityByLocality.js

@@ -1,5 +1,5 @@
 import { array_copy } from "../../../collection/array/array_copy.js";
-import { arraySetSortingDiff } from "../../../collection/array/arraySetSortingDiff.js";
+import { array_set_diff_sorting } from "../../../collection/array/array_set_diff_sorting.js";
 import { compareFaces } from "./compareFaces.js";
 import { computeTopoMeshVertexDuplicates } from "./computeTopoMeshVertexDuplicates.js";
 import { query_edge_other_vertex } from "./query/query_edge_other_vertex.js";
@@ -53,7 +53,7 @@ function connectEdges(a, b) {
     const faces_a = a.faces;
     const faces_b = b.faces;
 
-    const diff = arraySetSortingDiff(faces_a, faces_b, compareFaces);
+    const diff = array_set_diff_sorting(faces_a, faces_b, compareFaces);
 
     const unique_a = diff.uniqueA;
     const unique_b = diff.uniqueB;

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

@@ -1 +1 @@
-{"version":3,"file":"Vector2.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector2.js"],"names":[],"mappings":";AAmBA;IACI;;;;;OAKG;IACH,gBAJW,MAAM,MACN,MAAM,EAuBhB;IAbG;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV,0DAA6B;IAGjC;;;;OAIG;IACH,iBAHW,MAAM,EAAE,WACR,MAAM,QAOhB;IAED;;;;OAIG;IACH,eAHW,MAAM,EAAE,WACR,MAAM,QAKhB;IAED,iBAMC;IAED;;;;;OAKG;IACH,OAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAqBnB;IAED;;;;OAIG;IACH,aAHW,MAAM,KACN,MAAM,QAYhB;IAED;;;;OAIG;IACH,QAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,QAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,KACP,OAAO,QAIjB;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;OAGG;IACH,QAFa,OAAO,CAInB;IAED;;OAEG;IACH,cAIC;IAED;;;OAGG;IACH,OAFa,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;;OAKG;IACH,aAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAMnB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,YAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;OAGG;IACH,UAFa,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,gCAEC;IAED;;;OAGG;IACH,eAFW,MAAM,QAIhB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,6BAFa,OAAO,CAOnB;IAED;;;MAEC;IAED,0BAOC;IAED;;;OAGG;IACH,2CAGC;IAED;;;OAGG;IACH,6CAKC;IAED;;;OAGG;IACH,kDAGC;IAED;;;OAGG;IACH,oDAKC;IAED;;;OAGG;IACH,UAFa,OAAO,CAInB;IAED;;;;;;OAMG;IACH,uEAIC;IAED;;;;OAIG;IACH,8CAIC;IAED;;;;OAIG;IACH,iDAIC;IAED;;;;OAIG;IACH,qBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;;OAKG;IACH,eAJW,OAAO,KACP,OAAO,YACP,MAAM,QAOhB;IAED;;;OAGG;IACH,sBAFW,MAAM,EAAE,QAUlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;;OAKG;IACH,kBAJW,MAAM,KACN,MAAM,GACJ,MAAM,CAMlB;IAED;;;;;OAKG;IACH,eAJW,MAAM,KACN,MAAM,GACJ,MAAM,CAIlB;IAED;;;;OAIG;IACH,2BAHW,OAAO,GACL,MAAM,CAOlB;IAED;;OAEG;IACH,UAFa,MAAM,CAIlB;IAED;;OAEG;IACH,kBAWC;IAGD;;;OAGG;IACH,QAFa,MAAM,CAOlB;IAED;;;OAGG;IACH,cAFW,MAAM,QAahB;IAED;;;;;OAKG;IACH,6CAFa,OAAO,CAQnB;IAED,mBAEC;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;;OAKG;IACH,qBAJW,OAAO,cACP,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;OAMG;IACH,kBALW,MAAM,KACN,MAAM,cACN,MAAM,GACL,OAAO,CAKlB;IAUD,mBAEC;IAVD,gBAEC;IAUD,mBAEC;IAVD,gBAEC;IAsDL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAI3B;;OAEG;IACH,sBAznBe,MAAM,EAAE,WACR,MAAM,UAwnBS;IAE9B;;OAEG;IACH,uBA1oBe,MAAM,EAAE,WACR,MAAM,UAyoBU;IA3D3B,sDAKC;CACJ;;mBAIS,OAAO;qBAMP,OAAO;qBAMP,OAAO;sBAMP,OAAO;qBAMP,OAAO;oBAKP,OAAO;;;mBAlqBE,4BAA4B;4BAQnB,uBAAuB"}
+{"version":3,"file":"Vector2.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector2.js"],"names":[],"mappings":";AAmBA;IACI;;;;;OAKG;IACH,gBAJW,MAAM,MACN,MAAM,EAuBhB;IAbG;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV,0DAA6B;IAGjC;;;;OAIG;IACH,iBAHW,MAAM,EAAE,GAAC,YAAY,WACrB,MAAM,QAOhB;IAED;;;;OAIG;IACH,eAHW,MAAM,EAAE,GAAC,YAAY,WACrB,MAAM,QAKhB;IAED,iBAMC;IAED;;;;;OAKG;IACH,OAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAqBnB;IAED;;;;OAIG;IACH,aAHW,MAAM,KACN,MAAM,QAYhB;IAED;;;;OAIG;IACH,QAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,QAHW,MAAM,GACJ,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,KACP,OAAO,QAIjB;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;OAGG;IACH,QAFa,OAAO,CAInB;IAED;;OAEG;IACH,cAIC;IAED;;;OAGG;IACH,OAFa,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;;OAKG;IACH,aAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAMnB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,YAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;OAGG;IACH,UAFa,OAAO,CAInB;IAED;;;;;OAKG;IACH,QAJW,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,gCAEC;IAED;;;OAGG;IACH,eAFW,MAAM,QAIhB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,6BAFa,OAAO,CAOnB;IAED;;;MAEC;IAED,0BAOC;IAED;;;OAGG;IACH,2CAGC;IAED;;;OAGG;IACH,6CAKC;IAED;;;OAGG;IACH,kDAGC;IAED;;;OAGG;IACH,oDAKC;IAED;;;OAGG;IACH,UAFa,OAAO,CAInB;IAED;;;;;;OAMG;IACH,uEAIC;IAED;;;;OAIG;IACH,8CAIC;IAED;;;;OAIG;IACH,iDAIC;IAED;;;;OAIG;IACH,qBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;;OAKG;IACH,eAJW,OAAO,KACP,OAAO,YACP,MAAM,QAOhB;IAED;;;OAGG;IACH,sBAFW,MAAM,EAAE,QAUlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;;OAKG;IACH,kBAJW,MAAM,KACN,MAAM,GACJ,MAAM,CAMlB;IAED;;;;;OAKG;IACH,eAJW,MAAM,KACN,MAAM,GACJ,MAAM,CAIlB;IAED;;;;OAIG;IACH,2BAHW,OAAO,GACL,MAAM,CAOlB;IAED;;OAEG;IACH,UAFa,MAAM,CAIlB;IAED;;OAEG;IACH,kBAWC;IAGD;;;OAGG;IACH,QAFa,MAAM,CAOlB;IAED;;;OAGG;IACH,cAFW,MAAM,QAahB;IAED;;;;;OAKG;IACH,6CAFa,OAAO,CAQnB;IAED,mBAEC;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;;OAKG;IACH,qBAJW,OAAO,cACP,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;OAMG;IACH,kBALW,MAAM,KACN,MAAM,cACN,MAAM,GACL,OAAO,CAKlB;IAUD,mBAEC;IAVD,gBAEC;IAUD,mBAEC;IAVD,gBAEC;IAsDL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAI3B;;OAEG;IACH,sBAznBe,MAAM,EAAE,GAAC,YAAY,WACrB,MAAM,UAwnBS;IAE9B;;OAEG;IACH,uBA1oBe,MAAM,EAAE,GAAC,YAAY,WACrB,MAAM,UAyoBU;IA3D3B,sDAKC;CACJ;;mBAIS,OAAO;qBAMP,OAAO;qBAMP,OAAO;sBAMP,OAAO;qBAMP,OAAO;oBAKP,OAAO;;;mBAlqBE,4BAA4B;4BAQnB,uBAAuB"}

package/src/core/geom/Vector2.js

@@ -48,7 +48,7 @@ class Vector2 {
 
     /**
      *
-     * @param {number[]} array
+     * @param {number[]|Float32Array} array
      * @param {number} offset
      */
     fromArray(array, offset = 0) {
@@ -60,7 +60,7 @@ class Vector2 {
 
     /**
      *
-     * @param {number[]} array
+     * @param {number[]|Float32Array} array
      * @param {number} offset
      */
     toArray(array, offset = 0) {

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

@@ -1,5 +1,5 @@
 /**
- * Kaiser Bessel Function of third kind
+ * Kaiser Bessel Function of third kind (BesselJ(3,x))
  * NOTE: this is an approximation, if input value is sufficiently large - results will be way off
  * This is a Tylor series approximation
  * @see https://www.shadertoy.com/view/wlf3RH

package/src/core/math/bessel_3.js

@@ -1,5 +1,5 @@
 /**
- * Kaiser Bessel Function of third kind
+ * Kaiser Bessel Function of third kind (BesselJ(3,x))
  * NOTE: this is an approximation, if input value is sufficiently large - results will be way off
  * This is a Tylor series approximation
  * @see https://www.shadertoy.com/view/wlf3RH

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

@@ -1,5 +1,5 @@
 /**
- * Compute lowest positive integer divisor of a given "product" value
+ * Compute the largest positive integer divisor of a given "product" value
  * @param {number} product
  * @param {number} limit
  * @returns {number}

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

@@ -1 +1 @@
-{"version":3,"file":"computeWholeDivisorLow.d.ts","sourceRoot":"","sources":["../../../../src/core/math/computeWholeDivisorLow.js"],"names":[],"mappings":"AAEA;;;;;GAKG;AACH,gDAJW,MAAM,SACN,MAAM,GACJ,MAAM,CAuBlB"}
+{"version":3,"file":"computeWholeDivisorLow.d.ts","sourceRoot":"","sources":["../../../../src/core/math/computeWholeDivisorLow.js"],"names":[],"mappings":"AAEA;;;;;GAKG;AACH,gDAJW,MAAM,SACN,MAAM,GACJ,MAAM,CA0BlB"}

package/src/core/math/computeWholeDivisorLow.js

@@ -1,18 +1,21 @@
 import { assert } from "../assert.js";
 
 /**
- * Compute lowest positive integer divisor of a given "product" value
+ * Compute the largest positive integer divisor of a given "product" value
  * @param {number} product
  * @param {number} limit
  * @returns {number}
  */
 export function computeWholeDivisorLow(product, limit) {
-    assert.typeOf(product, 'number', 'product');
-    assert.typeOf(limit, 'number', 'limit');
 
+    assert.isNumber(product, 'product');
+    assert.isInteger(product, 'product');
 
-    assert.ok(Number.isInteger(product), `expected product to be an integer, instead got '${product}'`);
-    assert.ok(Number.isInteger(limit), `expected limit to be an integer, instead got '${limit}'`);
+    assert.isNumber(limit, 'limit');
+    assert.isInteger(limit, 'limit');
+
+    assert.isFiniteNumber(limit,'limit');
+    assert.isFiniteNumber(product,'product');
 
     let i = limit;
 

package/src/core/math/interval/separation1D.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"separation1D.d.ts","sourceRoot":"","sources":["../../../../../src/core/math/interval/separation1D.js"],"names":[],"mappings":"AAIA;;;;;;;;GAQG;AACH,iCANW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CASlB"}

package/src/core/math/{separation1D.js → interval/separation1D.js}

@@ -1,6 +1,6 @@
-import { assert } from "../assert.js";
-import { max2 } from "./max2.js";
-import { min2 } from "./min2.js";
+import { assert } from "../../assert.js";
+import { max2 } from "../max2.js";
+import { min2 } from "../min2.js";
 
 /**
  * Returns a number representing overlap distance between two 1D line segments.

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

@@ -1 +1 @@
-{"version":3,"file":"solveQuadratic.d.ts","sourceRoot":"","sources":["../../../../src/core/math/solveQuadratic.js"],"names":[],"mappings":"AAEA;;;;;;;;GAQG;AACH,uCAJW,MAAM,EAAE,iBACR,MAAM,KAJN,MAAM,KACN,MAAM,KACN,MAAM,GAGJ,MAAM,CAuClB"}
+{"version":3,"file":"solveQuadratic.d.ts","sourceRoot":"","sources":["../../../../src/core/math/solveQuadratic.js"],"names":[],"mappings":"AAGA;;;;;;;;GAQG;AACH,uCAJW,MAAM,EAAE,iBACR,MAAM,KAJN,MAAM,KACN,MAAM,KACN,MAAM,GAGJ,MAAM,CAgDlB"}

package/src/core/math/solveQuadratic.js

@@ -1,3 +1,4 @@
+import { assert } from "../assert.js";
 import { EPSILON } from "./EPSILON.js";
 
 /**
@@ -10,10 +11,17 @@ import { EPSILON } from "./EPSILON.js";
  * @returns {number} number of found solutions (roots)
  */
 export function solveQuadratic(result, result_offset, a, b, c) {
+    assert.isNumber(a, 'a');
+    assert.isNumber(b, 'b');
+    assert.isNumber(c, 'c');
+
 
     if (Math.abs(a) < EPSILON) {
+
         if (Math.abs(b) < EPSILON) {
+
             if (Math.abs(c) < EPSILON) {
+                // special case, all variables are close to 0
 
                 result[result_offset] = 0;
                 result[result_offset + 1] = 0;
@@ -21,12 +29,14 @@ export function solveQuadratic(result, result_offset, a, b, c) {
                 return 2;
 
             }
+
         } else {
             result[result_offset] = -c / b;
             result[result_offset + 1] = -c / b;
 
             return 2;
         }
+
     } else {
         // 2 root solution
 

package/src/core/model/object/compareValues.js

@@ -1,8 +1,8 @@
-import { compareStrings } from "../../primitives/strings/compareStrings.js";
-import { compareBooleans } from "../../primitives/boolean/compareBooleans.js";
+import { array_compare } from "../../collection/array/array_compare.js";
 import { extractFunctionBody } from "../../function/extractFunctionBody.js";
-import { compareArrays } from "../../primitives/array/compareArrays.js";
+import { compareBooleans } from "../../primitives/boolean/compareBooleans.js";
 import { number_compare_ascending } from "../../primitives/numbers/number_compare_ascending.js";
+import { compareStrings } from "../../primitives/strings/compareStrings.js";
 
 /**
  *
@@ -42,13 +42,13 @@ export function compareValues(a, b) {
         }
 
         if (Array.isArray(a) && Array.isArray(b)) {
-            return compareArrays(a, b, compareValues);
+            return array_compare(a, b, compareValues);
         }
 
         const aKeys = Object.keys(a);
         const bKeys = Object.keys(b);
 
-        const dKeys = compareArrays(aKeys, bKeys, compareStrings);
+        const dKeys = array_compare(aKeys, bKeys, compareStrings);
 
         if (dKeys !== 0) {
             return dKeys;

package/src/engine/EntityCreator.d.ts

@@ -1,14 +1,3 @@
-/**
- *
- * @param timeout
- * @param action
- * @returns {Entity}
- */
-export function createTimer({ timeout, action }: {
-    timeout: any;
-    action: any;
-}): Entity;
-import Entity from './ecs/Entity.js';
 /**
  *
  * @param {Vector3} [position]
@@ -21,5 +10,16 @@ import Entity from './ecs/Entity.js';
  * @returns {Entity}
  */
 export function createSound({ position, timeout, url, track, positioned, channel, volume }?: Vector3): Entity;
+/**
+ *
+ * @param timeout
+ * @param action
+ * @returns {Entity}
+ */
+export function createTimer({ timeout, action }: {
+    timeout: any;
+    action: any;
+}): Entity;
 import Vector3 from "../core/geom/Vector3.js";
+import Entity from './ecs/Entity.js';
 //# sourceMappingURL=EntityCreator.d.ts.map

package/src/engine/EntityCreator.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"EntityCreator.d.ts","sourceRoot":"","sources":["../../../src/engine/EntityCreator.js"],"names":[],"mappings":"AAwEA;;;;;GAKG;AACH;;;IAFa,MAAM,CAkBlB;mBA3FkB,iBAAiB;AASpC;;;;;;;;;;GAUG;AACH,6FATW,OAAO,GAOL,MAAM,CAiDlB;oBAtEmB,yBAAyB"}
+{"version":3,"file":"EntityCreator.d.ts","sourceRoot":"","sources":["../../../src/engine/EntityCreator.js"],"names":[],"mappings":"AAYA;;;;;;;;;;GAUG;AACH,6FATW,OAAO,GAOL,MAAM,CAiDlB;AAED;;;;;GAKG;AACH;;;IAFa,MAAM,CAkBlB;oBA9FmB,yBAAyB;mBAG1B,iBAAiB"}

package/src/engine/EntityCreator.js

@@ -21,7 +21,7 @@ import { SoundEmitterChannels } from "./sound/ecs/emitter/SoundEmitterSystem.js"
  * @param {number} [volume]
  * @returns {Entity}
  */
-function createSound(
+export function createSound(
     {
         position = Vector3.zero,
         timeout = 60,
@@ -92,6 +92,4 @@ export function createTimer({ timeout, action }) {
     }));
 
     return builder;
-}
-
-export { createSound };
+}

package/src/engine/ecs/EntityComponentDataset.js

@@ -1,7 +1,7 @@
 import { assert } from "../../core/assert.js";
 import { BitSet } from "../../core/binary/BitSet.js";
+import { array_set_diff } from "../../core/collection/array/array_set_diff.js";
 import { array_shrink_to_size } from "../../core/collection/array/array_shrink_to_size.js";
-import { arraySetDiff } from "../../core/collection/array/arraySetDiff.js";
 import { findSignalHandlerIndexByHandle } from "../../core/events/signal/findSignalHandlerIndexByHandle.js";
 import {
     findSignalHandlerIndexByHandleAndContext
@@ -463,7 +463,7 @@ export class EntityComponentDataset {
 
         const newComponentTypeCount = map.length;
 
-        const diff = arraySetDiff(map, this.componentTypeMap);
+        const diff = array_set_diff(map, this.componentTypeMap);
 
         const typesToAdd = diff.uniqueA;
         const typesToRemove = diff.uniqueB;
@@ -681,7 +681,7 @@ export class EntityComponentDataset {
      * @returns {boolean} false if no new classes were added, true if at least one new class was added
      */
     registerManyComponentTypes(types) {
-        const diff = arraySetDiff(types, this.componentTypeMap);
+        const diff = array_set_diff(types, this.componentTypeMap);
 
         if (diff.uniqueA.length === 0) {
             // all classes area already registered

package/src/engine/ecs/grid/HeightMap2NormalMap.js

@@ -1,7 +1,7 @@
 import ImageFilter from '../../graphics/filter/ImageFilter.js';
 import NormalMapShader from '../../graphics/shaders/NormalMapShader2.js';
 import { Sampler2D } from '../../graphics/texture/sampler/Sampler2D.js';
-import { sampler2DtoFloat32Texture } from "../../graphics/texture/sampler/sampler2DToFloat32Texture.js";
+import { sampler2d_to_float32_texture } from "../../graphics/texture/sampler/sampler2d_to_float32_texture.js";
 
 
 function convertChannel(v) {
@@ -51,7 +51,7 @@ function heightMap2NormalMap(renderer, sampler) {
     const width = sampler.width;
     const height = sampler.height;
 
-    const texture = sampler2DtoFloat32Texture(sampler);
+    const texture = sampler2d_to_float32_texture(sampler);
 
     //construct shader
     const shader = new NormalMapShader();

package/src/engine/ecs/grid/NormalMap2AOMap.js

@@ -3,7 +3,7 @@ import ImageFilter from '../../graphics/filter/ImageFilter.js';
 import AmbientOcclusionShader from '../../graphics/shaders/AmbientOcclusionShader.js';
 import { DenoiseShader } from "../../graphics/shaders/DenoiseShader.js";
 import { Sampler2D } from '../../graphics/texture/sampler/Sampler2D.js';
-import { sampler2DtoFloat32Texture } from "../../graphics/texture/sampler/sampler2DToFloat32Texture.js";
+import { sampler2d_to_float32_texture } from "../../graphics/texture/sampler/sampler2d_to_float32_texture.js";
 
 function filterResult2Texture(data, width, height) {
     const result = new DataTexture();
@@ -46,9 +46,9 @@ function normalMap2OcclusionMap(
     const width = resultSize.x;
     const height = resultSize.y;
     //
-    const normalTexture = sampler2DtoFloat32Texture(normalMap);
+    const normalTexture = sampler2d_to_float32_texture(normalMap);
 
-    const heightTexture = sampler2DtoFloat32Texture(heightMap);
+    const heightTexture = sampler2d_to_float32_texture(heightMap);
 
 
     //construct shader

package/src/engine/graphics/ecs/mesh/SkeletonUtils.js

@@ -1,5 +1,5 @@
 import { Bone } from "three";
-import { arrayPickBestElement } from "../../../../core/collection/array/arrayPickBestElement.js";
+import { array_pick_best_element } from "../../../../core/collection/array/array_pick_best_element.js";
 import { string_compute_similarity } from "../../../../core/primitives/strings/string_compute_similarity.js";
 import { extractName } from "../../../../extractName.js";
 import { BoneMapping } from "./skeleton/BoneMapping.js";
@@ -43,7 +43,7 @@ export function getSkeletonBone(component, boneName) {
     //bone not found
 
     //try to find similar bones
-    const bestMatch = arrayPickBestElement(bones, function (bone) {
+    const bestMatch = array_pick_best_element(bones, function (bone) {
         return string_compute_similarity(bone.name, boneName);
     });
 

package/src/engine/graphics/impostors/octahedral/util/build_cutout_from_atlas_by_alpha.d.ts

@@ -4,7 +4,7 @@
  * @param {ImpostorDescription} impostor
  * @param {THREE.WebGLRenderer} renderer
  * @param {number} vertex_limit how many vertices to use for the billboard, more vertices will provide better fit and lower area reducing overdraw at the cost of extra triangles. Prefer values between 4 and 16
- * @returns {number[]}
+ * @returns {Float32Array}
  */
-export function build_cutout_from_atlas_by_alpha({ impostor, renderer, vertex_limit }: ImpostorDescription): number[];
+export function build_cutout_from_atlas_by_alpha({ impostor, renderer, vertex_limit }: ImpostorDescription): Float32Array;
 //# sourceMappingURL=build_cutout_from_atlas_by_alpha.d.ts.map

package/src/engine/graphics/impostors/octahedral/util/build_cutout_from_atlas_by_alpha.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"build_cutout_from_atlas_by_alpha.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/impostors/octahedral/util/build_cutout_from_atlas_by_alpha.js"],"names":[],"mappings":"AAeA;;;;;;;GAOG;AACH,6GAFa,MAAM,EAAE,CAqHpB"}
+{"version":3,"file":"build_cutout_from_atlas_by_alpha.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/impostors/octahedral/util/build_cutout_from_atlas_by_alpha.js"],"names":[],"mappings":"AAQA;;;;;;;GAOG;AACH,6GAFa,YAAY,CA0ExB"}

package/src/engine/graphics/impostors/octahedral/util/build_cutout_from_atlas_by_alpha.js

@@ -1,17 +1,10 @@
 //
 
 import { assert } from "../../../../../core/assert.js";
-import { array_copy } from "../../../../../core/collection/array/array_copy.js";
-import { compute_polygon_area_2d } from "../../../../../core/geom/2d/compute_polygon_area_2d.js";
-import { convex_hull_jarvis_2d } from "../../../../../core/geom/2d/convex-hull/convex_hull_jarvis_2d.js";
-import { fixed_convex_hull_relaxation } from "../../../../../core/geom/2d/convex-hull/fixed_convex_hull_relaxation.js";
 import { max2 } from "../../../../../core/math/max2.js";
 import { renderSamplerFromTexture } from "../../../render/visibility/hiz/buildCanvasViewFromTexture.js";
 import { Sampler2D } from "../../../texture/sampler/Sampler2D.js";
-import {
-    make_edge_condition_channel_threshold
-} from "../../../texture/sampler/search/make_edge_condition_channel_threshold.js";
-import { sampler2d_find_pixels } from "../../../texture/sampler/search/sampler2d_find_pixels.js";
+import { sampler2d_build_cutout } from "./sampler2d_build_cutout.js";
 
 /**
  * TODO For this specific use-case with octahedral impostors we can get away tracing a circular shape in general.
@@ -19,7 +12,7 @@ import { sampler2d_find_pixels } from "../../../texture/sampler/search/sampler2d
  * @param {ImpostorDescription} impostor
  * @param {THREE.WebGLRenderer} renderer
  * @param {number} vertex_limit how many vertices to use for the billboard, more vertices will provide better fit and lower area reducing overdraw at the cost of extra triangles. Prefer values between 4 and 16
- * @returns {number[]}
+ * @returns {Float32Array}
  */
 export function build_cutout_from_atlas_by_alpha({
 
@@ -50,7 +43,7 @@ export function build_cutout_from_atlas_by_alpha({
     const frame_size = impostor.resolution / impostor.frame_count;
     const frame_size_i = Math.ceil(frame_size);
 
-    const accummulator = Sampler2D.uint8(1, frame_size_i, frame_size_i);
+    const accumulator = Sampler2D.uint8(1, frame_size_i, frame_size_i);
 
     for (let i = 0; i < impostor.frame_count; i++) {
 
@@ -76,44 +69,14 @@ export function build_cutout_from_atlas_by_alpha({
 
                         const local_index = y * frame_size_i + x;
 
-                        accummulator.data[local_index] = max2(v, accummulator.data[local_index]);
+                        accumulator.data[local_index] = max2(v, accumulator.data[local_index]);
                     }
                 }
             }
         }
     }
 
-
-    const edge_pixels = sampler2d_find_pixels(accummulator, make_edge_condition_channel_threshold(0, 5));
-
-    // filtered pixels
-    const filtered_edge_pixels = [];
-
-    const half_pixel = 0.5;
-
-    for (let i = 0; i < edge_pixels.length; i += 2) {
-        const x = edge_pixels[i];
-        const y = edge_pixels[i + 1];
-
-        filtered_edge_pixels.push(x - half_pixel, y - half_pixel);
-        filtered_edge_pixels.push(x + half_pixel, y - half_pixel);
-        filtered_edge_pixels.push(x + half_pixel, y + half_pixel);
-        filtered_edge_pixels.push(x - half_pixel, y + half_pixel);
-    }
-
-    const hull_point_indices = convex_hull_jarvis_2d(filtered_edge_pixels, filtered_edge_pixels.length / 2);
-
-    const hull_points = [];
-
-    for (let i = 0; i < hull_point_indices.length; i++) {
-        const hpi = hull_point_indices[i];
-
-        array_copy(filtered_edge_pixels, hpi * 2, hull_points, i * 2, 2);
-    }
-
-    const cutout_points = new Float32Array(vertex_limit * 2);
-
-    fixed_convex_hull_relaxation(cutout_points, 0, vertex_limit, hull_points, hull_point_indices.length);
+    const cutout_points = sampler2d_build_cutout(accumulator, vertex_limit, frame_size);
 
     // normalize points
     for (let i = 0; i < vertex_limit; i++) {
@@ -121,19 +84,6 @@ export function build_cutout_from_atlas_by_alpha({
         cutout_points[i * 2 + 1] /= frame_size;
     }
 
-    const area = compute_polygon_area_2d(cutout_points, vertex_limit);
-
-    console.log(`Cutout area: ${area}`);
-
-    if (area >= 0.99 && vertex_limit >= 4) {
-        // default back to aligned quad, produced result is a degradation
-        return new Float32Array([
-            0, 0,
-            1, 0,
-            1, 1,
-            0, 1
-        ]);
-    }
 
     return cutout_points;
 }

package/src/engine/graphics/impostors/octahedral/util/sampler2d_build_cutout.d.ts

@@ -0,0 +1,9 @@
+/**
+ *
+ * @param {Sampler2D} accumulator
+ * @param {number} vertex_limit
+ * @param {number} frame_size
+ * @returns {Float32Array}
+ */
+export function sampler2d_build_cutout(accumulator: Sampler2D, vertex_limit: number, frame_size: number): Float32Array;
+//# sourceMappingURL=sampler2d_build_cutout.d.ts.map

package/src/engine/graphics/impostors/octahedral/util/sampler2d_build_cutout.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"sampler2d_build_cutout.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/impostors/octahedral/util/sampler2d_build_cutout.js"],"names":[],"mappings":"AASA;;;;;;GAMG;AACH,6EAJW,MAAM,cACN,MAAM,GACJ,YAAY,CAiDxB"}