@woosh/meep-engine 2.107.0 → 2.108.2

package/src/core/color/rgb_to_luminance.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"rgb_to_luminance.d.ts","sourceRoot":"","sources":["../../../../src/core/color/rgb_to_luminance.js"],"names":[],"mappings":"AAAA;;;;;;GAMG;AACH,oCALW,MAAM,KACN,MAAM,KACN,MAAM,GACL,MAAM,CAIjB"}

package/src/core/color/rgb_to_luminance.js

@@ -0,0 +1,10 @@
+/**
+ *
+ * @param {number} r
+ * @param {number} g
+ * @param {number} b
+ * @return {number}
+ */
+export function rgb_to_luminance(r, g, b) {
+    return 0.299 * r + 0.587 * g + 0.114 * b;
+}

package/src/core/color/sRGB/linear_to_sRGB.d.ts

@@ -1,3 +1,4 @@
+export function convert_channel_linear_to_sRGB(c: any): number;
 /**
  *
  * @param {number[]|ArrayLike<number>} output

package/src/core/color/sRGB/linear_to_sRGB.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"linear_to_sRGB.d.ts","sourceRoot":"","sources":["../../../../../src/core/color/sRGB/linear_to_sRGB.js"],"names":[],"mappings":"AAQA;;;;;;GAMG;AACH,uCALW,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,iBAC1B,MAAM,SACN,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,gBAC1B,MAAM,QAahB"}
+{"version":3,"file":"linear_to_sRGB.d.ts","sourceRoot":"","sources":["../../../../../src/core/color/sRGB/linear_to_sRGB.js"],"names":[],"mappings":"AAAA,+DAMC;AAED;;;;;;GAMG;AACH,uCALW,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,iBAC1B,MAAM,SACN,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,gBAC1B,MAAM,QAahB"}

package/src/core/color/sRGB/linear_to_sRGB.js

@@ -1,4 +1,4 @@
-function convert(c) {
+export function convert_channel_linear_to_sRGB(c) {
     if (c < 0.0031308) {
         return c * 12.92;
     } else {
@@ -21,7 +21,7 @@ export function linear_to_sRGB(
     const g = input[input_offset + 1];
     const b = input[input_offset + 2];
 
-    output[output_offset] = convert(r);
-    output[output_offset + 1] = convert(g);
-    output[output_offset + 2] = convert(b);
+    output[output_offset] = convert_channel_linear_to_sRGB(r);
+    output[output_offset + 1] = convert_channel_linear_to_sRGB(g);
+    output[output_offset + 2] = convert_channel_linear_to_sRGB(b);
 }

package/src/core/geom/3d/aabb/AABB3.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"AABB3.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/aabb/AABB3.js"],"names":[],"mappings":"AAeA;;;GAGG;AACH;IACI;;;;;;;;;OASG;IACH,iBARW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,EAYhB;IAyCD,gBAEC;IA7BD,aAEC;IA6BD,gBAEC;IA5BD,aAEC;IA4BD,gBAEC;IA5BD,aAEC;IA4BD,gBAEC;IA3BD,aAEC;IA2BD,gBAEC;IA1BD,aAEC;IA0BD,gBAEC;IA1BD,aAEC;IAGG,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAGf;;;;;;OAMG;IACH,iBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;;;;OAOG;IACH,8BANW,MAAM,KACN,MAAM,KACN,MAAM,aACN,MAAM,GACJ,OAAO,CAMnB;IAED,4BAKC;IAED;;;OAGG;IACH,sBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,kBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,iBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFW,KAAK,QAIf;IAED;;;;;;;;OAQG;IACH,wFAsCC;IAED;;;;OAIG;IACH,cAHW,KAAK,GACH,OAAO,CAInB;IAED;;;;;;;;;OASG;IACH,iFAFa,OAAO,CASnB;IAED;;;;;;;;OAQG;IACH,iGA6BC;IAED,oCAEC;IAED,0BAEC;IAED;;;;;OAKG;IACH,cAJW,MAAM,KACN,MAAM,KACN,MAAM,QAOhB;IAED;;;;;;OAMG;IACH,oBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,MAAM,CAQlB;IAED;;;;OAIG;IACH,mBAHW,KAAK,GACH,MAAM,CAIlB;IAED;;;;;;;;;;OAUG;IACH,mBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAclB;IAED;;;;OAIG;IACH,wBAHW,KAAK,GACH,MAAM,CAIlB;IAED;;;;;;;;;;OAUG;IACH,sBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAsClB;IAED;;;;;;OAMG;IACH,qBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,OAAO,CA6BnB;IAED;;;;OAIG;IACH,iBAHW,KAAK,GACH,OAAO,CAInB;IAED;;;;;;;;;OASG;IACH,iBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,OAAO,CA6BnB;IAED;;;;;;;;;OASG;IACH,iBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,iBAHW,KAAK,GACH,OAAO,CAInB;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,oBAEC;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,qBAEC;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,oBAEC;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,kCAMC;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,iCAKC;IAGD;;;;;;;;OAQG;IACH,4EAEC;IAED,kGAEC;IAED;;;;OAIG;IACH,kCAFa,OAAO,CAMnB;IAED;;;;OAIG;IACH,oDAmBC;IAED;;;OAGG;IACH,mBAFW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,QAa5C;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,WAClC,MAAM,QAYhB;IAED;;;;OAIG;IACH,sBAHW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,WAClC,MAAM,QAchB;IAED;;;OAGG;IACH,wBAHW,MAAM,KAAK,OAWrB;IAED;;;;;;;OAOG;IACH,oCANW,MAAM,YACN,MAAM,YACN,MAAM,UACN,MAAM,GACJ,MAAM,CAQlB;IAED;;;;;;;OAOG;IACH,wBANW,MAAM,YACN,MAAM,YACN,MAAM,UACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,4BAFY,OAAO,CAMlB;IAED;;;OAGG;IACH,+BAHW,OAAO,GACL,OAAO,CAanB;IAED;;;;OAIG;IACH,0CAFa,MAAM,CAoBlB;IAED;;;;OAIG;IACH,sCAHW,MAAM,EAAE,GACN,MAAM,CAMlB;IAED;;;;OAIG;IACH;;QAFa,OAAO,CAYnB;IAED;;;;OAIG;IACH,gCAHW,UAAU,MAAM,CAAC,GAAC,MAAM,EAAE,GAAC,YAAY,GACrC,OAAO,CAYnB;IAED;;;OAGG;IACH,qBAFW,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,QAWjD;IAED;;;OAGG;IACH,YAFW,MAAM,QAUhB;IAED;;;OAGG;IACH,SAFa,KAAK,CAQjB;IAED;;;;;;;aAEC;IAIL;;;OAGG;IACH,kBAFU,OAAO,CAEM;IAEvB;;;;OAIG;IACH,iBAFU,MAAM,CAEM;IAx4BlB,mDAUC;CAg3BJ"}
+{"version":3,"file":"AABB3.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/aabb/AABB3.js"],"names":[],"mappings":"AAgBA;;;GAGG;AACH;IACI;;;;;;;;;OASG;IACH,iBARW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,EAYhB;IAyCD,gBAEC;IA7BD,aAEC;IA6BD,gBAEC;IA5BD,aAEC;IA4BD,gBAEC;IA5BD,aAEC;IA4BD,gBAEC;IA3BD,aAEC;IA2BD,gBAEC;IA1BD,aAEC;IA0BD,gBAEC;IA1BD,aAEC;IAGG,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAIX,QAAW;IAGf;;;;;;OAMG;IACH,iBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;;;;OAOG;IACH,8BANW,MAAM,KACN,MAAM,KACN,MAAM,aACN,MAAM,GACJ,OAAO,CAMnB;IAED,4BAKC;IAED;;;OAGG;IACH,sBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,kBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,iBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFW,KAAK,QAIf;IAED;;;;;;;;OAQG;IACH,wFAsCC;IAED;;;;OAIG;IACH,cAHW,KAAK,GACH,OAAO,CAInB;IAED;;;;;;;;;OASG;IACH,iFAFa,OAAO,CASnB;IAED;;;;;;;;OAQG;IACH,iGA6BC;IAED,oCAEC;IAED,0BAEC;IAED;;;;;OAKG;IACH,cAJW,MAAM,KACN,MAAM,KACN,MAAM,QAOhB;IAED;;;;;;OAMG;IACH,oBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,MAAM,CAQlB;IAED;;;;OAIG;IACH,mBAHW,KAAK,GACH,MAAM,CAIlB;IAED;;;;;;;;;;OAUG;IACH,mBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAclB;IAED;;;;OAIG;IACH,wBAHW,KAAK,GACH,MAAM,CAIlB;IAED;;;;;;;;;;OAUG;IACH,sBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAsClB;IAED;;;;;;OAMG;IACH,qBALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,OAAO,CA6BnB;IAED;;;;OAIG;IACH,iBAHW,KAAK,GACH,OAAO,CAInB;IAED;;;;;;;;;OASG;IACH,iBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,OAAO,CA6BnB;IAED;;;;;;;;;OASG;IACH,iBARW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,iBAHW,KAAK,GACH,OAAO,CAInB;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,oBAEC;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,qBAEC;IAED;;;OAGG;IACH,eAFa,MAAM,CAIlB;IAED,oBAEC;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,mBAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,mBAFW,OAAO,QAQjB;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,cAFa,MAAM,CAIlB;IAED,sBAEC;IAED;;;OAGG;IACH,mBAFW,OAAO,WASjB;IAGD;;;;;;;;OAQG;IACH,4EAEC;IAED,kGAEC;IAED;;;;OAIG;IACH,kCAFa,OAAO,CAMnB;IAED;;;;OAIG;IACH,oDAmBC;IAED;;;OAGG;IACH,mBAFW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,QAa5C;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,WAClC,MAAM,QAYhB;IAED;;;;OAIG;IACH,sBAHW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,WAClC,MAAM,QAchB;IAED;;;OAGG;IACH,wBAHW,MAAM,KAAK,OAWrB;IAED;;;;;;;OAOG;IACH,oCANW,MAAM,YACN,MAAM,YACN,MAAM,UACN,MAAM,GACJ,MAAM,CAQlB;IAED;;;;;;;OAOG;IACH,wBANW,MAAM,YACN,MAAM,YACN,MAAM,UACN,MAAM,GACJ,OAAO,CAInB;IAED;;;;OAIG;IACH,4BAFY,OAAO,CAMlB;IAED;;;OAGG;IACH,+BAHW,OAAO,GACL,OAAO,CAanB;IAED;;;;OAIG;IACH,0CAFa,MAAM,CAoBlB;IAED;;;;OAIG;IACH,sCAHW,MAAM,EAAE,GACN,MAAM,CAMlB;IAED;;;;OAIG;IACH;;QAFa,OAAO,CAYnB;IAED;;;;OAIG;IACH,gCAHW,UAAU,MAAM,CAAC,GAAC,MAAM,EAAE,GAAC,YAAY,GACrC,OAAO,CAYnB;IAED;;;OAGG;IACH,qBAFW,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,QAWjD;IAED;;;OAGG;IACH,YAFW,MAAM,QAUhB;IAED;;;OAGG;IACH,SAFa,KAAK,CAQjB;IAED;;;;;;;aAEC;IAIL;;;OAGG;IACH,kBAFU,OAAO,CAEM;IAEvB;;;;OAIG;IACH,iBAFU,MAAM,CAEM;IA14BlB,mDAUC;CAk3BJ;oBAr6BmB,kBAAkB"}

package/src/core/geom/3d/aabb/AABB3.js

@@ -1,4 +1,5 @@
 import { assert } from "../../../assert.js";
+import Vector3 from "../../Vector3.js";
 import computeMortonCode from "../morton/Morton.js";
 import { aabb3_array_intersects_point } from "./aabb3_array_intersects_point.js";
 import { aabb3_build_corners } from "./aabb3_build_corners.js";
@@ -630,11 +631,13 @@ export class AABB3 {
      * Get center position of the box
      * @param {Vector3} target where to write result
      */
-    getCenter(target) {
+    getCenter(target = new Vector3()) {
         const x = this.getCenterX();
         const y = this.getCenterY();
         const z = this.getCenterZ();
         target.set(x, y, z);
+
+        return target;
     }
 
 

package/src/core/geom/3d/aabb/aabb3_intersects_ray.js

@@ -35,7 +35,7 @@ export function aabb3_intersects_ray(
     const diff_x = origin_x - center_x;
 
 
-    if (diff_x * direction_x >= 0.0 && abs(diff_x) > extents_x) {
+    if (diff_x * direction_x >= 0 && abs(diff_x) > extents_x) {
         return false;
     }
 
@@ -47,7 +47,7 @@ export function aabb3_intersects_ray(
     const diff_y = origin_y - center_y;
 
 
-    if (diff_y * direction_y >= 0.0 && abs(diff_y) > extents_y) {
+    if (diff_y * direction_y >= 0 && abs(diff_y) > extents_y) {
         return false;
     }
 
@@ -59,7 +59,7 @@ export function aabb3_intersects_ray(
     const diff_z = origin_z - center_z;
 
 
-    if (diff_z * direction_z >= 0.0 && abs(diff_z) > extents_z) {
+    if (diff_z * direction_z >= 0 && abs(diff_z) > extents_z) {
         return false;
     }
 

package/src/core/geom/3d/shape/UnitSphereShape3D.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"UnitSphereShape3D.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/shape/UnitSphereShape3D.js"],"names":[],"mappings":"AAKA;;GAEG;AACH;IACI,qBAEC;IAED,2BAEC;IAED,wCAQC;IAED,4DAYC;IAED,mEAEC;IAED,6CAMC;IAED,oCAMC;IAED,kFAEC;IAED,eAEC;CACJ;;;;gCApE+B,sBAAsB"}
+{"version":3,"file":"UnitSphereShape3D.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/shape/UnitSphereShape3D.js"],"names":[],"mappings":"AAKA;;GAEG;AACH;IACI,qBAEC;IAED,2BAEC;IAED,wCAQC;IAED,4DAYC;IAED,mEAEC;IAED,6CAMC;IAED,oCAMC;IAED,kFAEC;IAED,eAEC;CACJ;;;;gCAlE+B,sBAAsB"}

package/src/core/geom/3d/shape/UnitSphereShape3D.js

@@ -1,28 +1,28 @@
-import { AbstractShape3D } from "./AbstractShape3D.js";
 import { randomPointInSphere } from "../../random/randomPointInSphere.js";
 import { v3_length } from "../../vec3/v3_length.js";
+import { AbstractShape3D } from "./AbstractShape3D.js";
 import { compute_signed_distance_gradient_by_sampling } from "./util/compute_signed_distance_gradient_by_sampling.js";
 
 /**
- * Sphere with diameter of 1 (radius = 0.5)
+ * Sphere with diameter of 2 (radius = 1)
  */
 export class UnitSphereShape3D extends AbstractShape3D {
     get volume() {
-        return 0.523598776;
+        return (4 / 3) * Math.PI;
     }
 
     get surface_area() {
-        return 3.14159265;
+        return Math.PI * 4;
     }
 
     compute_bounding_box(result) {
-        result[0] = -0.5;
-        result[1] = -0.5;
-        result[2] = -0.5;
+        result[0] = -1;
+        result[1] = -1;
+        result[2] = -1;
 
-        result[3] = 0.5;
-        result[4] = 0.5;
-        result[5] = 0.5;
+        result[3] = 1;
+        result[4] = 1;
+        result[5] = 1;
     }
 
     nearest_point_on_surface(result, reference) {
@@ -31,8 +31,8 @@ export class UnitSphereShape3D extends AbstractShape3D {
         const r_y = reference[1];
         const r_z = reference[2];
 
-        // normalize vector to 0.5 radius
-        const d = 0.5 / Math.hypot(r_x, r_y, r_z);
+        // normalize vector to 1 radius
+        const d = 1 / Math.hypot(r_x, r_y, r_z);
 
         result[0] = r_x * d;
         result[1] = r_y * d;
@@ -48,7 +48,7 @@ export class UnitSphereShape3D extends AbstractShape3D {
             point[0],
             point[1],
             point[2]
-        ) - 0.5;
+        ) - 1;
     }
 
     contains_point(point) {
@@ -56,7 +56,7 @@ export class UnitSphereShape3D extends AbstractShape3D {
         const y = point[1];
         const z = point[2];
 
-        return (x * x + y * y + z * z) < 0.25;
+        return (x * x + y * y + z * z) < 1;
     }
 
     sample_random_point_in_volume(result, result_offset, random) {

package/src/core/geom/3d/sphere/harmonics/README.md

@@ -8,8 +8,22 @@ There are basically 3 approaches,
 * Optimization of the polynomial (generally to remove negative values)
 * Piece-wise reconstruction
 
+## Coefficients of SH basis
+see [3] (Stupid Spherical Harmonics)
+
+        1.0 / (2.0 * M_SQRT_PI),    // 0: 0  0
+        -M_SQRT_3 / (2.0 * M_SQRT_PI),    // 1: 1 -1
+        M_SQRT_3 / (2.0 * M_SQRT_PI),    // 2: 1  0
+        -M_SQRT_3 / (2.0 * M_SQRT_PI),    // 3: 1  1
+        M_SQRT_15 / (2.0 * M_SQRT_PI),    // 4: 2 -2
+        -M_SQRT_15 / (2.0 * M_SQRT_PI),    // 5: 2 -1
+        M_SQRT_5 / (4.0 * M_SQRT_PI),    // 6: 2  0
+        -M_SQRT_15 / (2.0 * M_SQRT_PI),    // 7: 2  1
+        M_SQRT_15 / (4.0 * M_SQRT_PI)     // 8: 2  2
+
 ### References:
 
 ---
-* 2008 - GDC - Stupid Spherical Harmonics (SH) Tricks by Peter-Pike Sloan
-* [Computing product of 2 sphecial harmonics](https://github.com/ruba/cortex-vfx/blob/c256b7274516a02fa64a0260351b86f1d6ebe8c2/include/IECore/SphericalHarmonicsAlgo.inl)
+* [1] 2008 - GDC - Stupid Spherical Harmonics (SH) Tricks by Peter-Pike Sloan
+* [2] [Computing product of 2 sphecial harmonics](https://github.com/ruba/cortex-vfx/blob/c256b7274516a02fa64a0260351b86f1d6ebe8c2/include/IECore/SphericalHarmonicsAlgo.inl)
+* [3] Stupid Spherical Harmonics

package/src/core/geom/3d/sphere/harmonics/SH3_COEFFICIENTS.d.ts

@@ -0,0 +1,6 @@
+/**
+ * @readonly
+ * @type {Float32Array}
+ */
+export const SH3_COEFFICIENTS: Float32Array;
+//# sourceMappingURL=SH3_COEFFICIENTS.d.ts.map

package/src/core/geom/3d/sphere/harmonics/SH3_COEFFICIENTS.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"SH3_COEFFICIENTS.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/SH3_COEFFICIENTS.js"],"names":[],"mappings":"AAAA;;;GAGG;AACH,+BAFU,YAAY,CAcnB"}

package/src/core/geom/3d/sphere/harmonics/SH3_COEFFICIENTS.js

@@ -0,0 +1,17 @@
+/**
+ * @readonly
+ * @type {Float32Array}
+ */
+export const SH3_COEFFICIENTS = new Float32Array([
+    0.28209479177387814,    // 0: 0  0
+
+    -0.4886025119029199,    // 1: 1 -1
+    0.4886025119029199,     // 2: 1  0
+    -0.4886025119029199,    // 3: 1  1
+
+    1.0925484305920792,     // 4: 2 -2
+    -1.0925484305920792,    // 5: 2 -1
+    0.31539156525252005,    // 6: 2  0
+    -1.0925484305920792,    // 7: 2  1
+    0.5462742152960396      // 8: 2  2
+]);

package/src/core/geom/3d/sphere/harmonics/sh3_basis_at.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"sh3_basis_at.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_basis_at.js"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AACH,gCALW,MAAM,KACN,MAAM,KACN,MAAM,WACN,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,GAAC,YAAY,QAsB9D"}
+{"version":3,"file":"sh3_basis_at.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_basis_at.js"],"names":[],"mappings":"AAGA;;;;;;;;GAQG;AACH,gCALW,MAAM,KACN,MAAM,KACN,MAAM,WACN,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,GAAC,YAAY,QA+B9D"}

package/src/core/geom/3d/sphere/harmonics/sh3_basis_at.js

@@ -1,3 +1,6 @@
+import { assert } from "../../../../assert.js";
+
+
 /**
  * evaluate the basis functions
  * shBasis is an Array[ 9 ]
@@ -9,22 +12,31 @@
  */
 export function sh3_basis_at(x, y, z, shBasis) {
 
+    assert.isNumber(x,'x');
+    assert.notNaN(x,'x');
+
+    assert.isNumber(y,'y');
+    assert.notNaN(y,'y');
+
+    assert.isNumber(z,'z');
+    assert.notNaN(z,'z');
+
     // normal is assumed to be unit length
 
 
     // band 0
-    shBasis[0] = 0.282095;
+    shBasis[0] = 0.28209479177387814;
 
     // band 1
-    shBasis[1] = 0.488603 * y;
-    shBasis[2] = 0.488603 * z;
-    shBasis[3] = 0.488603 * x;
+    shBasis[1] = -0.4886025119029199 * y;
+    shBasis[2] = 0.4886025119029199 * z;
+    shBasis[3] = -0.4886025119029199 * x;
 
     // band 2
-    shBasis[4] = 1.092548 * x * y;
-    shBasis[5] = 1.092548 * y * z;
-    shBasis[6] = 0.315392 * (3 * z * z - 1);
-    shBasis[7] = 1.092548 * x * z;
-    shBasis[8] = 0.546274 * (x * x - y * y);
+    shBasis[4] = 1.0925484305920792 * x * y;
+    shBasis[5] = -1.0925484305920792 * y * z;
+    shBasis[6] = 0.31539156525252005 * (3 * z * z - 1);
+    shBasis[7] = -1.0925484305920792 * x * z;
+    shBasis[8] = 0.5462742152960396 * (x * x - y * y);
 
 }

package/src/core/geom/3d/sphere/harmonics/sh3_dering_optimize_positive.d.ts

@@ -1,10 +1,10 @@
 /**
  *
- * @param {number[]|Float32Array} result
+ * @param {number[]|Float32Array|Float64Array} result
  * @param {number} result_offset
- * @param {number[]|Float32Array} harmonics
+ * @param {number[]|Float32Array|Float64Array} harmonics
  * @param {number} harmonics_offset
  * @param {number} dimension_count
  */
-export function sh3_dering_optimize_positive(result: number[] | Float32Array, result_offset: number, harmonics: number[] | Float32Array, harmonics_offset: number, dimension_count?: number): void;
+export function sh3_dering_optimize_positive(result: number[] | Float32Array | Float64Array, result_offset: number, harmonics: number[] | Float32Array | Float64Array, harmonics_offset: number, dimension_count?: number): void;
 //# sourceMappingURL=sh3_dering_optimize_positive.d.ts.map

package/src/core/geom/3d/sphere/harmonics/sh3_dering_optimize_positive.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"sh3_dering_optimize_positive.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_dering_optimize_positive.js"],"names":[],"mappings":"AA6cA;;;;;;;GAOG;AACH,qDANW,MAAM,EAAE,GAAC,YAAY,iBACrB,MAAM,aACN,MAAM,EAAE,GAAC,YAAY,oBACrB,MAAM,oBACN,MAAM,QAIhB"}
+{"version":3,"file":"sh3_dering_optimize_positive.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_dering_optimize_positive.js"],"names":[],"mappings":"AAkSA;;;;;;;GAOG;AACH,qDANW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,iBAClC,MAAM,aACN,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,oBAClC,MAAM,oBACN,MAAM,QAIhB"}

package/src/core/geom/3d/sphere/harmonics/sh3_dering_optimize_positive.js

@@ -2,172 +2,27 @@ import { mat3, vec3 } from "gl-matrix";
 import { assert } from "../../../../assert.js";
 import { array_copy } from "../../../../collection/array/array_copy.js";
 import { min2 } from "../../../../math/min2.js";
+import { SH3_COEFFICIENTS } from "./SH3_COEFFICIENTS.js";
+import { sh3_rotate } from "./sh3_rotate.js";
+import { sh_index } from "./sh_index.js";
 
 /*
 @see https://github.com/Bestmaker602/olacziy/blob/212b64ea5f1856b390cdf7629801243f76a4466d/libs/ibl/src/CubemapSH.cpp
  */
 
 
-const F_PI = 3.14159265358979323846264338327950288;
-const F_SQRT1_2 = 0.707106781186547524400844362104849039;
-const M_SQRT_3 = 1.7320508076;
-
-const M_SQRT_PI = 1.7724538509;
-const M_SQRT_5 = 2.2360679775;
-const M_SQRT_15 = 3.8729833462;
-
-/**
- *
- * @param {number[]} M 5x5 matrix
- * @param {number[]} x vec5
- * @return {number[]} vec5
- */
-function multiply_5d(M, x) {
-    return [
-        M[0] * x[0] + M[5] * x[1] + M[10] * x[2] + M[15] * x[3] + M[20] * x[4],
-        M[1] * x[0] + M[6] * x[1] + M[11] * x[2] + M[16] * x[3] + M[21] * x[4],
-        M[2] * x[0] + M[7] * x[1] + M[12] * x[2] + M[17] * x[3] + M[22] * x[4],
-        M[3] * x[0] + M[8] * x[1] + M[13] * x[2] + M[18] * x[3] + M[23] * x[4],
-        M[4] * x[0] + M[9] * x[1] + M[14] * x[2] + M[19] * x[3] + M[24] * x[4]
-    ];
-}
-
-/**
- * compute Index of spherical harmonics coefficient
- * @param {number} m signed offset from zonal harmonic
- * @param {number} l band
- * @return {number}
- */
-function sh_index(m, l) {
-    return l * (l + 1) + m;
-}
-
-
-/*
- * utilities to rotate very low order spherical harmonics (up to 3rd band)
- * @param {number[]} band1
- * @param {number[]} M 3x3 matrix
- * @returns {number[]}
- */
-function rotateShericalHarmonicBand1(band1, M) {
-
-    // inverse() is not constexpr -- so we pre-calculate it in mathematica
-    //
-    //    constexpr float3 N0{ 1, 0, 0 };
-    //    constexpr float3 N1{ 0, 1, 0 };
-    //    constexpr float3 N2{ 0, 0, 1 };
-    //
-    //    constexpr mat3f A1 = { // this is the projection of N0, N1, N2 to SH space
-    //            float3{ -N0.y, N0.z, -N0.x },
-    //            float3{ -N1.y, N1.z, -N1.x },
-    //            float3{ -N2.y, N2.z, -N2.x }
-    //    };
-    //
-    //    const mat3f invA1 = inverse(A1);
-
-    const invA1TimesK = [
-        0, -1, 0,
-        0, 0, 1,
-        -1, 0, 0
-    ];
-
-    const R1OverK = [
-        -M[1], M[2], -M[0],
-        -M[4], M[5], -M[3],
-        -M[7], M[8], -M[6]
-    ];
-
-
-    mat3.scale(invA1TimesK, invA1TimesK, band1);
-
-
-    mat3.multiply(R1OverK, R1OverK, invA1TimesK);
-
-    return R1OverK;
-}
-
-/**
- *  This projects a vec3 to SH2/k space (i.e. we premultiply by 1/k)
- *  below can't be constexpr
- * @param {number} x
- * @param {number} y
- * @param {number} z
- * @return {number[]} vec5
- */
-function project_v3_to_sh(x, y, z) {
-    return [
-        (y * x),
-        -(y * z),
-        1 / (2 * M_SQRT_3) * ((3 * z * z - 1)),
-        -(z * x),
-        0.5 * ((x * x - y * y))
-    ];
-}
-
-/**
- *
- * @param {number[]} band2 vec5
- * @param {number[]} M mat3
- * @return {number[]} vec5
- */
-function rotateShericalHarmonicBand2(band2, M) {
-    const n = F_SQRT1_2;
-
-    //  Below we precompute (with help of Mathematica):
-    //    constexpr float3 N0{ 1, 0, 0 };
-    //    constexpr float3 N1{ 0, 0, 1 };
-    //    constexpr float3 N2{ n, n, 0 };
-    //    constexpr float3 N3{ n, 0, n };
-    //    constexpr float3 N4{ 0, n, n };
-    //    constexpr float M_SQRT_PI = 1.7724538509f;
-    //    constexpr float M_SQRT_15 = 3.8729833462f;
-    //    constexpr float k = M_SQRT_15 / (2.0f * M_SQRT_PI);
-    //    --> k * inverse(mat5{project(N0), project(N1), project(N2), project(N3), project(N4)})
-    const invATimesK = [
-        0, 1, 2, 0, 0,
-        -1, 0, 0, 0, -2,
-        0, M_SQRT_3, 0, 0, 0,
-        1, 1, 0, -2, 0,
-        2, 1, 0, 0, 0
-    ];
-
-    // this is: invA * k * band2
-    // 5x5 matrix by vec5 (this a lot of zeroes and constants, which the compiler should eliminate)
-    const invATimesKTimesBand2 = multiply_5d(invATimesK, band2);
-
-    // this is: mat5{project(N0), project(N1), project(N2), project(N3), project(N4)} / k
-    // (the 1/k comes from project(), see above)
-    const ROverK =
-        project_v3_to_sh(M[0], M[1], M[2])                  // M * N0
-            .concat(project_v3_to_sh(M[6], M[7], M[8]))                  // M * N1
-            .concat(project_v3_to_sh(n * (M[0] + M[3]), n * (M[1] + M[4]), n * (M[2] + M[5])))     // M * N2
-            .concat(project_v3_to_sh(n * (M[0] + M[6]), n * (M[1] + M[7]), n * (M[2] + M[8])))     // M * N3
-            .concat(project_v3_to_sh(n * (M[3] + M[6]), n * (M[4] + M[7]), n * (M[5] + M[8])))      // M * N4
-    ;
-
-    // notice how "k" disappears
-    // this is: (R / k) * (invA * k) * band2 == R * invA * band2
-    const result = multiply_5d(ROverK, invATimesKTimesBand2);
-
-    return result;
-}
-
 /**
  * SH from environment with high dynamic range (or high frequencies -- high dynamic range creates
  * high frequencies) exhibit "ringing" and negative values when reconstructed.
  * To mitigate this, we need to low-pass the input image -- or equivalently window the SH by
  * coefficient that tapper towards zero with the band.
  *
- * We use ideas and techniques from
- *    Stupid Spherical Harmonics (SH)
- *    Deringing Spherical Harmonics
- * by Peter-Pike Sloan
- * https://www.ppsloan.org/publications/shdering.pdf
- * @param {number} w
+ * We use ideas and techniques from Stupid Spherical Harmonics (SH) Deringing Spherical Harmonics by Peter-Pike Sloan ( see https://www.ppsloan.org/publications/shdering.pdf )
+ * @param {number} w cutoff
  * @param {number} l
  * @returns {number}
  */
-function sincWindow(l, w) {
+function sinc_window(l, w) {
     if (l === 0) {
         return 1.0;
     } else if (l >= w) {
@@ -176,7 +31,7 @@ function sincWindow(l, w) {
 
     // we use a sinc window scaled to the desired window size in bands units
     // a sinc window only has zonal harmonics
-    let x = ((F_PI) * l) / w;
+    let x = (Math.PI * l) / w;
 
     x = Math.sin(x) / x;
 
@@ -188,52 +43,15 @@ function sincWindow(l, w) {
 }
 
 /**
- *
- * @param {number[]|Float32Array} result sh3
- * @param {number[]} sh input
- * @param {number[]} M mat3
- */
-function rotate_sh3_bands(result, sh, M) {
-
-    const b0 = sh[0];
-    const band1 = [sh[1], sh[2], sh[3]];
-    const b1 = rotateShericalHarmonicBand1(band1, M);
-    const band2 = [sh[4], sh[5], sh[6], sh[7], sh[8]];
-    const b2 = rotateShericalHarmonicBand2(band2, M);
-
-    result[0] = b0;
-
-    result[1] = b1[0];
-    result[2] = b1[1];
-    result[3] = b1[2];
-
-    result[4] = b2[0];
-    result[5] = b2[1];
-    result[6] = b2[2];
-    result[7] = b2[3];
-    result[8] = b2[4];
-}
-
-/**
- *
- * @param {number[]} input_sh3 sh3
+ * Minimum values encoded on the spherical harmonics
+ * Lets us find ringing artifacts, which are cases where harmonics encode a negative value
+ * @param {number[]|Float32Array} input_sh3 sh3
  * @returns {number}
  */
-function shmin(input_sh3) {
+function sh3_min(input_sh3) {
     // See "Deringing Spherical Harmonics" by Peter-Pike Sloan
     // https://www.ppsloan.org/publications/shdering.pdf
 
-    const A = [
-        1.0 / (2.0 * M_SQRT_PI),    // 0: 0  0
-        -M_SQRT_3 / (2.0 * M_SQRT_PI),    // 1: 1 -1
-        M_SQRT_3 / (2.0 * M_SQRT_PI),    // 2: 1  0
-        -M_SQRT_3 / (2.0 * M_SQRT_PI),    // 3: 1  1
-        M_SQRT_15 / (2.0 * M_SQRT_PI),    // 4: 2 -2
-        -M_SQRT_15 / (2.0 * M_SQRT_PI),    // 5: 2 -1
-        M_SQRT_5 / (4.0 * M_SQRT_PI),    // 6: 2  0
-        -M_SQRT_15 / (2.0 * M_SQRT_PI),    // 7: 2  1
-        M_SQRT_15 / (4.0 * M_SQRT_PI)     // 8: 2  2
-    ];
 
     // first this to do is to rotate the SH to align Z with the optimal linear direction
     const dir = vec3.fromValues(-input_sh3[3], -input_sh3[1], input_sh3[2]);
@@ -259,7 +77,7 @@ function shmin(input_sh3) {
 
     const f = new Float32Array(9);
 
-    rotate_sh3_bands(f, input_sh3, M);
+    sh3_rotate(f, input_sh3, M);
     // here we're guaranteed to have normalize(float3{ -f[3], -f[1], f[2] }) == { 0, 0, 1 }
 
 
@@ -270,7 +88,7 @@ function shmin(input_sh3) {
     // of z such as:  m2min = -m2max * (1 - z^2) =  m2max * z^2 - m2max
     //      with m2max = A[8] * std::sqrt(f[8] * f[8] + f[4] * f[4]);
     // We can therefore include this in the ZH min computation (which is function of z^2 as well)
-    const m2max = A[8] * Math.sqrt(f[8] * f[8] + f[4] * f[4]);
+    const m2max = SH3_COEFFICIENTS[8] * Math.sqrt(f[8] * f[8] + f[4] * f[4]);
 
     // Find the min of the zonal harmonics
     // -----------------------------------
@@ -286,9 +104,9 @@ function shmin(input_sh3) {
     //
     // We also needs to check that -1 < z < 1, otherwise the min is either in z = -1 or 1
     //
-    const a = 3 * A[6] * f[6] + m2max;
-    const b = A[2] * f[2];
-    const c = A[0] * f[0] - A[6] * f[6] - m2max;
+    const a = 3 * SH3_COEFFICIENTS[6] * f[6] + m2max;
+    const b = SH3_COEFFICIENTS[2] * f[2];
+    const c = SH3_COEFFICIENTS[0] * f[0] - SH3_COEFFICIENTS[6] * f[6] - m2max;
 
     const zmin = -b / (2.0 * a);
     const m0min_z = a * zmin * zmin + b * zmin + c;
@@ -304,7 +122,7 @@ function shmin(input_sh3) {
     // The function considered is:
     //        Y(x, y, z) = A[5] * f[5] * s.y * s.z
     //                   + A[7] * f[7] * s.z * s.x
-    const d = A[4] * Math.sqrt(f[5] * f[5] + f[7] * f[7]);
+    const d = SH3_COEFFICIENTS[4] * Math.sqrt(f[5] * f[5] + f[7] * f[7]);
 
     // the |m|=1 function is minimal in -0.5 -- use that to skip the Newton's loop when possible
     let minimum = m0min - 0.5 * d;
@@ -315,7 +133,7 @@ function shmin(input_sh3) {
 
 
         let dz;
-        let z = -F_SQRT1_2;   // we start guessing at the min of |m|=1 function
+        let z = -Math.SQRT1_2;   // we start guessing at the min of |m|=1 function
         do {
             minimum = func_(z, a, b, c, d); // evaluate our function
             dz = increment_(z, a, b, d); // refine our guess by this amount
@@ -366,14 +184,22 @@ function increment_(x, a, b, d) {
  * @param {number} numBands
  */
 function windowing(f, cutoff, numBands) {
+
     for (let l = 0; l < numBands; l++) {
-        const w = sincWindow(l, cutoff);
+
+        const w = sinc_window(l, cutoff);
+
         f[sh_index(0, l)] *= w;
+
         for (let m = 1; m <= l; m++) {
+
             f[sh_index(-m, l)] *= w;
             f[sh_index(m, l)] *= w;
+
         }
+
     }
+
     return f;
 }
 
@@ -390,7 +216,8 @@ function windowing(f, cutoff, numBands) {
 function windowSH(
     output, output_offset,
     input, input_offset,
-    numBands, channel_count, cutoff = 0
+    numBands, channel_count,
+    cutoff = 0
 ) {
     assert.isNonNegativeInteger(channel_count, 'channel_count');
     assert.greaterThan(channel_count, 0, 'channel_count must be greater than 0');
@@ -424,7 +251,9 @@ function windowSH(
 
                 const m = 0.5 * (l + r);
 
-                if (shmin(windowing(SH, m, numBands)) < 0) {
+                windowing(SH, m, numBands);
+
+                if (sh3_min(SH) < 0) {
                     r = m;
                 } else {
                     l = m;
@@ -442,7 +271,7 @@ function windowSH(
 
     for (let l = 0; l < numBands; l++) {
         // scale each band individually
-        let w = sincWindow(l, cutoff);
+        const w = sinc_window(l, cutoff);
 
         // write zonal harmonic
         for (let i = 0; i < channel_count; i++) {
@@ -461,9 +290,9 @@ function windowSH(
 
 /**
  *
- * @param {number[]|Float32Array} result
+ * @param {number[]|Float32Array|Float64Array} result
  * @param {number} result_offset
- * @param {number[]|Float32Array} harmonics
+ * @param {number[]|Float32Array|Float64Array} harmonics
  * @param {number} harmonics_offset
  * @param {number} dimension_count
  */