@woosh/meep-engine 2.108.4 → 2.109.1

package/src/core/binary/split_by_2.spec.js

@@ -0,0 +1,22 @@
+import { int32_to_binary_string } from "./int32_to_binary_string.js";
+import { split_by_2 } from "./split_by_2.js";
+
+test("zero", () => {
+
+    expect(split_by_2(0)).toEqual(0);
+
+});
+
+test("max value", () => {
+
+    expect(
+        int32_to_binary_string(split_by_2(0xFFFFFFFF))
+    ).toEqual("01010101 01010101 01010101 01010101");
+
+});
+
+test("max 16bit value", () => {
+    expect(
+        int32_to_binary_string(split_by_2(0b1111111111111111))
+    ).toEqual("01010101 01010101 01010101 01010101");
+});

package/src/core/bvh2/bvh3/ebvh_build_for_geometry_morton.d.ts

@@ -4,6 +4,9 @@
  * @param {BVH} bvh
  * @param {number[]|Float32Array} index_array
  * @param {number[]|Float32Array} position_array
+ * @param {Uint32Array} [morton_codes]
+ * @param {AABB3} [bounds]
  */
-export function ebvh_build_for_geometry_morton(bvh: BVH, index_array: number[] | Float32Array, position_array: number[] | Float32Array): void;
+export function ebvh_build_for_geometry_morton(bvh: BVH, index_array: number[] | Float32Array, position_array: number[] | Float32Array, morton_codes?: Uint32Array, bounds?: AABB3): void;
+import { AABB3 } from "../../geom/3d/aabb/AABB3.js";
 //# sourceMappingURL=ebvh_build_for_geometry_morton.d.ts.map

package/src/core/bvh2/bvh3/ebvh_build_for_geometry_morton.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"ebvh_build_for_geometry_morton.d.ts","sourceRoot":"","sources":["../../../../../src/core/bvh2/bvh3/ebvh_build_for_geometry_morton.js"],"names":[],"mappings":"AAWA;;;;;;GAMG;AACH,sEAHW,MAAM,EAAE,GAAC,YAAY,kBACrB,MAAM,EAAE,GAAC,YAAY,QAoF/B"}
+{"version":3,"file":"ebvh_build_for_geometry_morton.d.ts","sourceRoot":"","sources":["../../../../../src/core/bvh2/bvh3/ebvh_build_for_geometry_morton.js"],"names":[],"mappings":"AAWA;;;;;;;;GAQG;AACH,sEALW,MAAM,EAAE,GAAC,YAAY,kBACrB,MAAM,EAAE,GAAC,YAAY,iBACrB,WAAW,WACX,KAAK,QA2Ff;sBA1GqB,6BAA6B"}

package/src/core/bvh2/bvh3/ebvh_build_for_geometry_morton.js

@@ -15,8 +15,16 @@ import { ebvh_build_hierarchy } from "./ebvh_build_hierarchy.js";
  * @param {BVH} bvh
  * @param {number[]|Float32Array} index_array
  * @param {number[]|Float32Array} position_array
+ * @param {Uint32Array} [morton_codes]
+ * @param {AABB3} [bounds]
  */
-export function ebvh_build_for_geometry_morton(bvh, index_array, position_array) {
+export function ebvh_build_for_geometry_morton(
+    bvh,
+    index_array,
+    position_array,
+    morton_codes = new Uint32Array(index_array.length / 3),
+    bounds = new AABB3()
+) {
     assert.defined(bvh, 'bvh');
 
     assert.defined(index_array, 'index_array');
@@ -25,16 +33,19 @@ export function ebvh_build_for_geometry_morton(bvh, index_array, position_array)
     assert.defined(position_array, 'position_array');
     assert.isArrayLike(position_array, 'position_array');
 
+    assert.isArrayLike(morton_codes, 'morton_codes');
+
     // clear out existing BVH
     bvh.release_all();
 
-    const aabb3 = new AABB3();
-
     // get bounds for the entire geometry
-    aabb3_from_v3_array(aabb3, position_array, position_array.length);
+    aabb3_from_v3_array(bounds, position_array, position_array.length);
 
     // allocate nodes
     const tri_count = index_array.length / 3;
+
+    assert.greaterThanOrEqual(morton_codes.length, tri_count);
+
     const node_leaf_count = tri_count;
     const node_bin_count = max2(0, node_leaf_count - 1);
 
@@ -59,9 +70,7 @@ export function ebvh_build_for_geometry_morton(bvh, index_array, position_array)
         sorted_triangle_order[i] = i;
     }
 
-    const morton_codes = new Uint32Array(tri_count);
-
-    build_triangle_morton_codes(morton_codes, tri_count, index_array, position_array, aabb3.x0, aabb3.y0, aabb3.z0, aabb3.x1, aabb3.y1, aabb3.z1);
+    build_triangle_morton_codes(morton_codes, tri_count, index_array, position_array, bounds.x0, bounds.y0, bounds.z0, bounds.x1, bounds.y1, bounds.z1);
 
     // sort leaves by morton codes
     array_quick_sort_by_lookup(sorted_triangle_order, morton_codes, 0, tri_count - 1);

package/src/core/bvh2/bvh3/ebvh_build_hierarchy.js

@@ -35,8 +35,8 @@ export function ebvh_build_hierarchy(
         let cursor = 0;
 
         while (cursor + 1 < unprocessed_node_count) {
-            const child_1 = unprocessed_nodes[cursor++];
             const child_2 = unprocessed_nodes[cursor++];
+            const child_1 = unprocessed_nodes[cursor++];
 
             const parent = node_pool[used_index++];
 

package/src/core/bvh2/bvh3/query/bvh_query_user_data_ray.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bvh_query_user_data_ray.d.ts","sourceRoot":"","sources":["../../../../../../src/core/bvh2/bvh3/query/bvh_query_user_data_ray.js"],"names":[],"mappings":"AAMA;;;;;;;;;;;;;GAaG;AACH,wDAXW,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,GACJ,MAAM,CAwElB"}
+{"version":3,"file":"bvh_query_user_data_ray.d.ts","sourceRoot":"","sources":["../../../../../../src/core/bvh2/bvh3/query/bvh_query_user_data_ray.js"],"names":[],"mappings":"AAMA;;;;;;;;;;;;;GAaG;AACH,wDAXW,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,GACJ,MAAM,CA6ElB"}

package/src/core/bvh2/bvh3/query/bvh_query_user_data_ray.js

@@ -29,14 +29,19 @@ export function bvh_query_user_data_ray(
         return 0;
     }
 
+    /**
+     * Move stack pointer to local variable scope to avoid de-referencing inside the loop
+     * @type {number}
+     */
+    let pointer = stack.pointer;
 
     /**
      *
      * @type {number}
      */
-    const stack_top = stack.pointer++;
+    const stack_top = pointer;
 
-    stack[stack_top] = root;
+    stack[pointer++] = root;
 
     let result_cursor = result_offset;
 
@@ -47,13 +52,13 @@ export function bvh_query_user_data_ray(
     const uint32 = bvh.__data_uint32;
 
     do {
-        stack.pointer--;
+        --pointer;
 
         /**
          *
          * @type {number}
          */
-        const node = stack[stack.pointer];
+        const node = stack[pointer];
 
         const address = node * ELEMENT_WORD_COUNT;
 
@@ -77,15 +82,15 @@ export function bvh_query_user_data_ray(
             // this is not a leaf node, push children onto traversal stack
             const child_2 = uint32[address + COLUMN_CHILD_2];
 
-            stack[stack.pointer++] = child_2;
-            stack[stack.pointer++] = child_1;
+            stack[pointer++] = child_2;
+            stack[pointer++] = child_1;
 
         } else {
             // leaf node
 
             result[result_cursor++] = uint32[address + COLUMN_USER_DATA];
         }
-    } while (stack.pointer > stack_top);
+    } while (pointer > stack_top);
 
     return result_cursor - result_offset;
 }

package/src/core/geom/2d/circle/circle_circumference.d.ts

@@ -0,0 +1,7 @@
+/**
+ *
+ * @param {number} radius
+ * @returns {number}
+ */
+export function circle_circumference(radius: number): number;
+//# sourceMappingURL=circle_circumference.d.ts.map

package/src/core/geom/2d/circle/circle_circumference.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"circle_circumference.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/2d/circle/circle_circumference.js"],"names":[],"mappings":"AAEA;;;;GAIG;AACH,6CAHW,MAAM,GACJ,MAAM,CAMlB"}

package/src/core/geom/2d/circle/circle_circumference.js

@@ -0,0 +1,12 @@
+import { assert } from "../../../assert.js";
+
+/**
+ *
+ * @param {number} radius
+ * @returns {number}
+ */
+export function circle_circumference(radius) {
+    assert.isNumber(radius, 'radius');
+
+    return 2 * Math.PI * radius;
+}

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

@@ -1 +1 @@
-{"version":3,"file":"compute_polygon_area_2d.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/compute_polygon_area_2d.js"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,gDAJW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,eAClC,MAAM,GACJ,MAAM,CA2BlB"}
+{"version":3,"file":"compute_polygon_area_2d.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/2d/compute_polygon_area_2d.js"],"names":[],"mappings":"AAAA;;;;;GAKG;AACH,gDAJW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,eAClC,MAAM,GACJ,MAAM,CA6BlB"}

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

@@ -8,11 +8,13 @@ export function compute_polygon_area_2d(points, point_count) {
     let sum = 0;
 
     for (let i = 0; i < point_count - 1; i++) {
-        const x0 = points[(i) * 2];
-        const y0 = points[(i) * 2 + 1];
+        const i2 = i * 2;
 
-        const x1 = points[(i + 1) * 2];
-        const y1 = points[(i + 1) * 2 + 1];
+        const x0 = points[i2];
+        const y0 = points[i2 + 1];
+
+        const x1 = points[i2 + 2];
+        const y1 = points[i2 + 3];
 
         sum += x0 * y1 - y0 * x1;
     }

package/src/core/geom/2d/lt-grid/LooseTightGrid.js

@@ -1,5 +1,5 @@
 import { assert } from "../../../assert.js";
-import { de_interleave_2_bits } from "../../../binary/de_interleave_bits_by_2.js";
+import { de_interleave_2_bits } from "../../../binary/de_interleave_2_bits.js";
 import { split_by_2 } from "../../../binary/split_by_2.js";
 import { UINT32_MAX } from "../../../binary/UINT32_MAX.js";
 import { SCRATCH_UINT32_TRAVERSAL_STACK } from "../../../collection/SCRATCH_UINT32_TRAVERSAL_STACK.js";

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

@@ -8,8 +8,8 @@
  * @param {number} by
  * @param {number} cx
  * @param {number} cy
- * @param {number} px
- * @param {number} py
+ * @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;
 //# sourceMappingURL=triangle2d_get_barycentric.d.ts.map

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

@@ -10,8 +10,8 @@ import { v2_dot } from "../vec2/v2_dot.js";
  * @param {number} by
  * @param {number} cx
  * @param {number} cy
- * @param {number} px
- * @param {number} py
+ * @param {number} px reference point X
+ * @param {number} py reference point Y
  */
 export function triangle2d_get_barycentric(
     output, output_offset,

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

@@ -0,0 +1,19 @@
+/**
+ * 3D ray
+ */
+export class Ray3 {
+    origin: Vector3;
+    direction: Vector3;
+    /**
+     *
+     * @param {number[]|mat4|Float32Array} m4
+     */
+    applyMatrix4(m4: number[] | mat4 | Float32Array): void;
+    /**
+     *
+     * @param {Ray3} other
+     */
+    copy(other: Ray3): void;
+}
+import Vector3 from "../Vector3.js";
+//# sourceMappingURL=Ray3.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"Ray3.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/3d/Ray3.js"],"names":[],"mappings":"AAEA;;GAEG;AACH;IACI,gBAAsB;IACtB,mBAAyB;IAEzB;;;OAGG;IACH,iBAFW,MAAM,EAAE,UAAM,YAAY,QAKpC;IAED;;;OAGG;IACH,YAFW,IAAI,QAKd;CACJ;oBA1BmB,eAAe"}

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

@@ -0,0 +1,27 @@
+import Vector3 from "../Vector3.js";
+
+/**
+ * 3D ray
+ */
+export class Ray3 {
+    origin = new Vector3()
+    direction = new Vector3()
+
+    /**
+     *
+     * @param {number[]|mat4|Float32Array} m4
+     */
+    applyMatrix4(m4) {
+        this.origin.applyMatrix4(m4);
+        this.direction.applyDirectionMatrix4(m4);
+    }
+
+    /**
+     *
+     * @param {Ray3} other
+     */
+    copy(other) {
+        this.origin.copy(other.origin);
+        this.direction.copy(other.direction);
+    }
+}

package/src/core/geom/3d/ray/ray3_shift_origin_along_direction.d.ts

@@ -0,0 +1,10 @@
+/**
+ *
+ * @param {number[]|ArrayLike<number>|Float32Array} output 6 component vector, [origin_x, origin_y, origin_z, direction_x, direction_y, direction_z]
+ * @param {number} output_offset
+ * @param {number[]|ArrayLike<number>|Float32Array} input 6 component vector, [origin_x, origin_y, origin_z, direction_x, direction_y, direction_z]
+ * @param {number} input_offset
+ * @param {number} distance
+ */
+export function ray3_shift_origin_along_direction(output: number[] | ArrayLike<number> | Float32Array, output_offset: number, input: number[] | ArrayLike<number> | Float32Array, input_offset: number, distance: number): void;
+//# sourceMappingURL=ray3_shift_origin_along_direction.d.ts.map

package/src/core/geom/3d/ray/ray3_shift_origin_along_direction.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ray3_shift_origin_along_direction.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/ray/ray3_shift_origin_along_direction.js"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,0DANW,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,iBACvC,MAAM,SACN,MAAM,EAAE,GAAC,UAAU,MAAM,CAAC,GAAC,YAAY,gBACvC,MAAM,YACN,MAAM,QAyBhB"}

package/src/core/geom/3d/ray/ray3_shift_origin_along_direction.js

@@ -0,0 +1,35 @@
+import { v3_shift_along_direction } from "../../vec3/v3_shift_along_direction.js";
+
+/**
+ *
+ * @param {number[]|ArrayLike<number>|Float32Array} output 6 component vector, [origin_x, origin_y, origin_z, direction_x, direction_y, direction_z]
+ * @param {number} output_offset
+ * @param {number[]|ArrayLike<number>|Float32Array} input 6 component vector, [origin_x, origin_y, origin_z, direction_x, direction_y, direction_z]
+ * @param {number} input_offset
+ * @param {number} distance
+ */
+export function ray3_shift_origin_along_direction(
+    output, output_offset,
+    input, input_offset,
+    distance
+) {
+    const origin_x = input[input_offset + 0];
+    const origin_y = input[input_offset + 1];
+    const origin_z = input[input_offset + 2];
+
+    const direction_x = input[input_offset + 3];
+    const direction_y = input[input_offset + 4];
+    const direction_z = input[input_offset + 5];
+
+    v3_shift_along_direction(
+        output, 0,
+        origin_x, origin_y, origin_z,
+        direction_x, direction_y, direction_z,
+        distance
+    );
+
+    output[output_offset + 3] = direction_x;
+    output[output_offset + 4] = direction_y;
+    output[output_offset + 5] = direction_z;
+}
+

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

@@ -0,0 +1 @@
+//# sourceMappingURL=compute_legendre_polynomial.d.ts.map

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

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

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

@@ -0,0 +1,9 @@
+/**
+ * Performs convolution of 3-band spherical harmonic coefficients with a cosine kernel
+ * @param {number[]} out
+ * @param {number} out_offset
+ * @param {number[]} input
+ * @param {number} input_offset
+ */
+export function sh3_convolve_with_cosine_kernel(out: number[], out_offset: number, input: number[], input_offset: number): void;
+//# sourceMappingURL=sh3_convolve_with_cosine_kernel.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"sh3_convolve_with_cosine_kernel.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_convolve_with_cosine_kernel.js"],"names":[],"mappings":"AAKA;;;;;;GAMG;AACH,qDALW,MAAM,EAAE,cACR,MAAM,SACN,MAAM,EAAE,gBACR,MAAM,QAchB"}

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

@@ -0,0 +1,26 @@
+// Cosine kernel for SH
+const CosineA0 = Math.PI;
+const CosineA1 = (2.0 * Math.PI) / 3.0;
+const CosineA2 = Math.PI / 4.0;
+
+/**
+ * Performs convolution of 3-band spherical harmonic coefficients with a cosine kernel
+ * @param {number[]} out
+ * @param {number} out_offset
+ * @param {number[]} input
+ * @param {number} input_offset
+ */
+export function sh3_convolve_with_cosine_kernel(out, out_offset, input, input_offset) {
+    out[out_offset + 0] = input[input_offset + 0] * CosineA0;
+
+    out[out_offset + 1] = input[input_offset + 1] * CosineA1;
+    out[out_offset + 2] = input[input_offset + 2] * CosineA1;
+    out[out_offset + 3] = input[input_offset + 3] * CosineA1;
+
+    out[out_offset + 4] = input[input_offset + 4] * CosineA2;
+    out[out_offset + 5] = input[input_offset + 5] * CosineA2;
+    out[out_offset + 6] = input[input_offset + 6] * CosineA2;
+    out[out_offset + 7] = input[input_offset + 7] * CosineA2;
+    out[out_offset + 8] = input[input_offset + 8] * CosineA2;
+}
+

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":"AAkSA;;;;;;;GAOG;AACH,qDANW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,iBAClC,MAAM,aACN,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,oBAClC,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":"AAgSA;;;;;;;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

@@ -3,7 +3,7 @@ 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 { sh3_rotate_bl } from "./sh3_rotate_bl.js";
 import { sh_index } from "./sh_index.js";
 
 /*
@@ -74,11 +74,9 @@ function sh3_min(input_sh3) {
     array_copy(y_axis, 0, M, 3, 3);
     array_copy(dir, 0, M, 6, 3);
 
-    mat3.transpose(M, M);
-
     const f = new Float32Array(9);
 
-    sh3_rotate(f, input_sh3, M);
+    sh3_rotate_bl(f, input_sh3, M);
     // here we're guaranteed to have normalize(float3{ -f[3], -f[1], f[2] }) == { 0, 0, 1 }
 
 

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

@@ -0,0 +1,8 @@
+/**
+ *
+ * @param {number[]|Float32Array} result sh3
+ * @param {number[]} sh input
+ * @param {number[]} rotation 3x3 rotation matrix
+ */
+export function sh3_rotate_bl(result: number[] | Float32Array, sh: number[], rotation: number[]): number[] | Float32Array;
+//# sourceMappingURL=sh3_rotate_bl.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"sh3_rotate_bl.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_rotate_bl.js"],"names":[],"mappings":"AAEA;;;;;GAKG;AACH,sCAJW,MAAM,EAAE,GAAC,YAAY,MACrB,MAAM,EAAE,YACR,MAAM,EAAE,2BAuFlB"}

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

@@ -0,0 +1,94 @@
+const scratch_mat_5x5 = new Float64Array(25);
+
+/**
+ *
+ * @param {number[]|Float32Array} result sh3
+ * @param {number[]} sh input
+ * @param {number[]} rotation 3x3 rotation matrix
+ */
+export function sh3_rotate_bl(result, sh, rotation) {
+    // Implementation adapted from BakerLab by https://github.com/TheRealMJP
+
+    const r00 = rotation[0];
+    const r10 = rotation[1];
+    const r20 = rotation[2];
+
+    const r01 = rotation[3];
+    const r11 = rotation[4];
+    const r21 = rotation[5];
+
+    const r02 = rotation[6];
+    const r12 = rotation[7];
+    const r22 = rotation[8];
+
+
+    // Constant
+    result[0] = sh[0];
+
+    // Linear
+    result[1] = r11 * sh[1] - r12 * sh[2] + r10 * sh[3];
+    result[2] = -r21 * sh[1] + r22 * sh[2] - r20 * sh[3];
+    result[3] = r01 * sh[1] - r02 * sh[2] + r00 * sh[3];
+
+    // Quadratic
+    const t41 = r01 * r00;
+    const t43 = r11 * r10;
+    const t48 = r11 * r12;
+    const t50 = r01 * r02;
+    const t55 = r02 * r02;
+    const t57 = r22 * r22;
+    const t58 = r12 * r12;
+    const t61 = r00 * r02;
+    const t63 = r10 * r12;
+    const t68 = r10 * r10;
+    const t70 = r01 * r01;
+    const t72 = r11 * r11;
+    const t74 = r00 * r00;
+    const t76 = r21 * r21;
+    const t78 = r20 * r20;
+
+    const v173 = 0.1732050808e1;
+    const v577 = 0.5773502693e0;
+    const v115 = 0.1154700539e1;
+    const v288 = 0.2886751347e0;
+    const v866 = 0.8660254040e0;
+
+    scratch_mat_5x5[0] = r11 * r00 + r01 * r10;
+    scratch_mat_5x5[1] = -r01 * r12 - r11 * r02;
+    scratch_mat_5x5[2] = v173 * r02 * r12;
+    scratch_mat_5x5[3] = -r10 * r02 - r00 * r12;
+    scratch_mat_5x5[4] = r00 * r10 - r01 * r11;
+    scratch_mat_5x5[5] = -r11 * r20 - r21 * r10;
+    scratch_mat_5x5[6] = r11 * r22 + r21 * r12;
+    scratch_mat_5x5[7] = -v173 * r22 * r12;
+    scratch_mat_5x5[8] = r20 * r12 + r10 * r22;
+    scratch_mat_5x5[9] = -r10 * r20 + r11 * r21;
+    scratch_mat_5x5[10] = -v577 * (t41 + t43) + v115 * r21 * r20;
+    scratch_mat_5x5[11] = v577 * (t48 + t50) - v115 * r21 * r22;
+    scratch_mat_5x5[12] = -0.5 * (t55 + t58) + t57;
+    scratch_mat_5x5[13] = v577 * (t61 + t63) - v115 * r20 * r22;
+    scratch_mat_5x5[14] = v288 * (t70 - t68 + t72 - t74) - v577 * (t76 - t78);
+    scratch_mat_5x5[15] = -r01 * r20 - r21 * r00;
+    scratch_mat_5x5[16] = r01 * r22 + r21 * r02;
+    scratch_mat_5x5[17] = -v173 * r22 * r02;
+    scratch_mat_5x5[18] = r00 * r22 + r20 * r02;
+    scratch_mat_5x5[19] = -r00 * r20 + r01 * r21;
+    scratch_mat_5x5[20] = t41 - t43;
+    scratch_mat_5x5[21] = -t50 + t48;
+    scratch_mat_5x5[22] = v866 * (t55 - t58);
+    scratch_mat_5x5[23] = t63 - t61;
+    scratch_mat_5x5[24] = 0.5 * (t74 - t68 - t70 + t72);
+
+    for (let i = 0; i < 5; ++i) {
+        const base = i * 5;
+
+        result[4 + i] = scratch_mat_5x5[base + 0] * sh[4]
+            + scratch_mat_5x5[base + 1] * sh[5]
+            + scratch_mat_5x5[base + 2] * sh[6]
+            + scratch_mat_5x5[base + 3] * sh[7]
+            + scratch_mat_5x5[base + 4] * sh[8]
+        ;
+    }
+
+    return result;
+}

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

@@ -0,0 +1,8 @@
+/**
+ *
+ * @param {number[]|Float32Array} result sh3
+ * @param {number[]} sh input
+ * @param {number[]} rotation 3x3 rotation matrix
+ */
+export function sh3_rotate_filament(result: number[] | Float32Array, sh: number[], rotation: number[]): void;
+//# sourceMappingURL=sh3_rotate_filament.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"sh3_rotate_filament.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_rotate_filament.js"],"names":[],"mappings":"AAGA;;;;;GAKG;AACH,4CAJW,MAAM,EAAE,GAAC,YAAY,MACrB,MAAM,EAAE,YACR,MAAM,EAAE,QAuBlB"}

package/src/core/geom/3d/sphere/harmonics/{sh3_rotate.js → sh3_rotate_filament.js}

@@ -7,7 +7,9 @@ import { sh_rotate_band2 } from "./sh_rotate_band2.js";
  * @param {number[]} sh input
  * @param {number[]} rotation 3x3 rotation matrix
  */
-export function sh3_rotate(result, sh, rotation) {
+export function sh3_rotate_filament(result, sh, rotation) {
+
+    // Adapted from filament https://github.com/google/filament/blob/e2da13f81734ccd6d6f9935061b771801682619e/libs/ibl/src/CubemapSH.cpp
 
     const b0 = sh[0];
     const band1 = [sh[1], sh[2], sh[3]];
@@ -26,4 +28,5 @@ export function sh3_rotate(result, sh, rotation) {
     result[6] = b2[2];
     result[7] = b2[3];
     result[8] = b2[4];
-}
+}
+

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

@@ -0,0 +1,12 @@
+/**
+ *
+ * @param {number[]|Float32Array|Float64Array} out
+ * @param {number} out_offset
+ * @param {number[]|Float32Array|Float64Array} basis
+ * @param {number} basis_offset
+ * @param {number[]|Float32Array|Float64Array} vector
+ * @param {number} vector_offset
+ * @param {number} dimensions
+ */
+export function sh3_vector_accumulate(out: number[] | Float32Array | Float64Array, out_offset: number, basis: number[] | Float32Array | Float64Array, basis_offset: number, vector: number[] | Float32Array | Float64Array, vector_offset: number, dimensions: number): void;
+//# sourceMappingURL=sh3_vector_accumulate.d.ts.map

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

@@ -0,0 +1 @@
+{"version":3,"file":"sh3_vector_accumulate.d.ts","sourceRoot":"","sources":["../../../../../../../src/core/geom/3d/sphere/harmonics/sh3_vector_accumulate.js"],"names":[],"mappings":"AAEA;;;;;;;;;GASG;AACH,2CARW,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,cAClC,MAAM,SACN,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,gBAClC,MAAM,UACN,MAAM,EAAE,GAAC,YAAY,GAAC,YAAY,iBAClC,MAAM,cACN,MAAM,QA6BhB"}

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

@@ -0,0 +1,40 @@
+import { assert } from "../../../../assert.js";
+
+/**
+ *
+ * @param {number[]|Float32Array|Float64Array} out
+ * @param {number} out_offset
+ * @param {number[]|Float32Array|Float64Array} basis
+ * @param {number} basis_offset
+ * @param {number[]|Float32Array|Float64Array} vector
+ * @param {number} vector_offset
+ * @param {number} dimensions
+ */
+export function sh3_vector_accumulate(
+    out, out_offset,
+    basis, basis_offset,
+    vector, vector_offset, dimensions
+) {
+
+    assert.isArrayLike(out, 'out');
+    assert.isNonNegativeInteger(out_offset, 'out_offset');
+
+    assert.isArrayLike(basis, 'basis');
+    assert.isNonNegativeInteger(basis_offset, 'basis_offset');
+
+    assert.isNonNegativeInteger(dimensions, 'dimensions');
+
+    for (let j = 0; j < 9; j++) {
+
+        const out_coefficient_offset = out_offset + j * dimensions;
+
+        const basis_value = basis[basis_offset + j];
+
+        for (let i = 0; i < dimensions; i++) {
+            const scalar = vector[vector_offset + i];
+
+            out[out_coefficient_offset + i] += basis_value * scalar;
+        }
+
+    }
+}

package/src/core/geom/3d/sphere/harmonics/sh3_vector_accumulate.spec.d.ts

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

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

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