@woosh/meep-engine 2.114.0 → 2.115.0

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

@@ -13,9 +13,9 @@ import { intersect_ray_2d } from "../intersect_ray_2d.js";
  * @param {number[]|Float32Array} output result will be written here
  * @param {number} output_offset offset into output array where to write the result
  * @param {number} output_point_count number of points to comprise bounding polygon
- * @param {number[]} input_points Points to be bounded, must be an ordered loop
+ * @param {number[]} input_points Points to be bounded, must be an ordered loop. Must form a convex hull
  * @param {number} input_point_count number of points in the input to consider
- * @returns {boolean} whether or not a valid result could be found
+ * @returns {number} number of points added to the output
  *
  * @see 2014 "Implementation of linear minimum area enclosing triangle algorithm" Ovidiu Pârvu & David Gilbert
  * @see https://github.com/MagnusTiberius/humus3/blob/9d43412d6e6be869e224817f0a966fe48a1af40f/Util/ConvexHull.cpp#L281
@@ -24,9 +24,11 @@ export function fixed_convex_hull_humus(
     output, output_offset, output_point_count,
     input_points, input_point_count
 ) {
+
     if (output_point_count >= input_point_count) {
         // special case, input polygon is smaller or equal to in cardinality to what's desired
         array_copy(input_points, 0, output, output_offset, input_point_count * 2);
+
         return input_point_count;
     }
 
@@ -61,6 +63,7 @@ export function fixed_convex_hull_humus(
         combos++;
 
         for (let i = 0; i < k; i++) {
+            // pick point indices to represent a line
             const i0 = subset[i];
             const i1 = subset[(i + 1) % k];
 
@@ -143,7 +146,11 @@ export function fixed_convex_hull_humus(
         process_combination(s);
     }
 
-    console.log(`combinations:  ${combos}, valid: ${combos_valid} (${(combos_valid * 100 / combos).toFixed(2)}%)`);
+    // console.log(`combinations:  ${combos}, valid: ${combos_valid} (${(combos_valid * 100 / combos).toFixed(2)}%)`);
+    if(combos_valid === 0){
+        // no valid combinations
+        return 0;
+    }
 
-    return combos_valid > 0;
+    return output_point_count;
 }

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

@@ -1,7 +1,7 @@
 import { array_copy } from "../../../collection/array/array_copy.js";
 import { Uint32Heap } from "../../../collection/heap/Uint32Heap.js";
 import { intersect_ray_2d } from "../intersect_ray_2d.js";
-import { compute_triangle_area_2d } from "../compute_triangle_area_2d.js";
+import { triangle2d_compute_area } from "../triangle2d_compute_area.js";
 
 
 const v2_temp = new Float32Array(2);
@@ -158,7 +158,7 @@ function compute_collapse_cost(edge, vertex_buffer) {
     const v1x = vertex_buffer[edge.v1 * 2];
     const v1y = vertex_buffer[edge.v1 * 2 + 1];
 
-    const area = compute_triangle_area_2d(
+    const area = triangle2d_compute_area(
         v0x, v0y, v2_temp[0], v2_temp[1], v1x, v1y
     );
 

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/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 { compareArrays } from "../../collection/array/compareArrays.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";
 
 /**
  *

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/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/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"}