@woosh/meep-engine 2.109.6 → 2.109.7

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

@@ -1,86 +0,0 @@
-/**
- * Math.sqrt(3)
- * @type {number}
- */
-const SQRT_3 = 1.7320508075688772;
-
-/**
- *
- * @param {number[]|Float32Array} 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]
-    ];
-}
-
-/**
- * This projects a vec3 to SH2/k space (i.e. we premultiply by 1/k)
- * these are second band SH factors, multiplied by 1/K ( K = Math.sqrt(15) / (2*Math.sqrt(Math.PI)) )
- * @param {number[]|Float32Array} output
- * @param {number} output_offset
- * @param {number} x
- * @param {number} y
- * @param {number} z
- */
-function project_v3_to_sh2k(output, output_offset, x, y, z) {
-    output[output_offset] = (y * x)
-    output[output_offset + 1] = -(y * z)
-    output[output_offset + 2] = 0.28867513459481287 * (3 * z * z - 1)
-    output[output_offset + 3] = -(z * x)
-    output[output_offset + 4] = 0.5 * (x * x - y * y)
-}
-
-/**
- *
- * @param {number[]} band2 vec5
- * @param {number[]} M mat3
- * @return {number[]} vec5
- */
-export function sh_rotate_band2(band2, M) {
-    const n = Math.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, 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 = new Float32Array(25);
-
-    project_v3_to_sh2k(ROverK, 0, M[0], M[1], M[2])                  // M * N0
-    project_v3_to_sh2k(ROverK, 5, M[6], M[7], M[8])                  // M * N1
-    project_v3_to_sh2k(ROverK, 10, n * (M[0] + M[3]), n * (M[1] + M[4]), n * (M[2] + M[5]))     // M * N2
-    project_v3_to_sh2k(ROverK, 15, n * (M[0] + M[6]), n * (M[1] + M[7]), n * (M[2] + M[8]))     // M * N3
-    project_v3_to_sh2k(ROverK, 20, 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;
-}

package/src/engine/graphics/sh3/SH3VisualisationMaterial.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"SH3VisualisationMaterial.d.ts","sourceRoot":"","sources":["../../../../../src/engine/graphics/sh3/SH3VisualisationMaterial.js"],"names":[],"mappings":"AA0DA;IAEI,cAoBC;CAEJ;+BAlFuC,OAAO"}

package/src/engine/graphics/sh3/SH3VisualisationMaterial.js

@@ -1,83 +0,0 @@
-import { ShaderMaterial, Vector3 } from "three";
-
-const shader_vx = `
-varying vec3 vNormal;
-
-void main() {
-    
-    vNormal = normalize( normalMatrix * normal );
-    
-    gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
-
-}
-`;
-
-const shader_fg = `
-#define RECIPROCAL_PI 0.318309886
-
-vec3 inverseTransformDirection( in vec3 normal, in mat4 matrix ) {
-	// matrix is assumed to be orthogonal
-	return normalize( ( vec4( normal, 0.0 ) * matrix ).xyz );
-}
-
-// source: https://graphics.stanford.edu/papers/envmap/envmap.pdf
-vec3 shGetIrradianceAt( in vec3 normal, in vec3 shCoefficients[ 9 ] ) {
-	// normal is assumed to have unit length
-	float x = normal.x, y = normal.y, z = normal.z;
-	
-	// band 0
-	vec3 result = shCoefficients[ 0 ] * 0.8862269254527579;
-	
-	// band 1
-	result -= shCoefficients[ 1 ] * 1.0233267079464885 * y;
-	result += shCoefficients[ 2 ] * 1.0233267079464885 * z;
-	result -= shCoefficients[ 3 ] * 1.0233267079464885 * x;
-	
-	// band 2
-	result += shCoefficients[ 4 ] * 0.8580855308097834 * x * y;
-	result -= shCoefficients[ 5 ] * 0.8580855308097834 * y * z;
-	result += shCoefficients[ 6 ] * ( 0.7431238683011272 * z * z - 0.24770795610037571 );
-	result -= shCoefficients[ 7 ] * 0.8580855308097834 * x * z;
-	result += shCoefficients[ 8 ] * 0.4290427654048917 * ( x * x - y * y );
-	
-	return result;
-}
-
-uniform vec3 sh[ 9 ]; // sh coefficients
-uniform float intensity; // light probe intensity
-varying vec3 vNormal;
-
-void main() {
-	vec3 normal = normalize( vNormal );
-	vec3 worldNormal = normalize(inverseTransformDirection( normal, viewMatrix ));
-	vec3 irradiance = shGetIrradianceAt( worldNormal, sh );
-	vec3 outgoingLight = 1.0  * irradiance * intensity;
-	gl_FragColor = linearToOutputTexel( vec4( outgoingLight, 1.0 ) );
-}
-`;
-
-export class SH3VisualisationMaterial extends ShaderMaterial {
-
-    constructor() {
-        const sh = [];
-
-        for (let i = 0; i < 9; i++) {
-            sh.push(new Vector3());
-        }
-
-        super({
-            fragmentShader: shader_fg,
-            vertexShader: shader_vx,
-            lights: false,
-            uniforms: {
-                intensity: {
-                    value: 1,
-                },
-                sh: {
-                    value: sh
-                }
-            }
-        })
-    }
-
-}

package/src/engine/graphics/sh3/gi/material/build_fragment_shader.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"build_fragment_shader.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/gi/material/build_fragment_shader.js"],"names":[],"mappings":"AAsDA;;;GAGG;AACH,8CAFW,MAAM,UAchB"}

package/src/engine/graphics/sh3/gi/material/build_fragment_shader.js

@@ -1,71 +0,0 @@
-import { insert_after } from "../../../../../core/primitives/strings/insert_after.js";
-
-const PREAMBLE = `
-
-        uniform float lpv_f_intensity;
-        
-        varying vec3[9] lpv_values;
-        
-        // vec3 shGetIrradianceAt( in vec3 normal, in vec3 shCoefficients[ 9 ] ) {
-        //     // normal is assumed to have unit length
-        //     float x = normal.x, y = normal.y, z = normal.z;
-        //    
-        //     // band 0
-        //     vec3 result = shCoefficients[ 0 ] * 0.886227;
-        //    
-        //     // band 1
-        //     result -= shCoefficients[ 1 ] * 2.0 * 0.511664 * y;
-        //     result += shCoefficients[ 2 ] * 2.0 * 0.511664 * z;
-        //     result -= shCoefficients[ 3 ] * 2.0 * 0.511664 * x;
-        //    
-        //     // band 2
-        //     result += shCoefficients[ 4 ] * 2.0 * 0.429043 * x * y;
-        //     result -= shCoefficients[ 5 ] * 2.0 * 0.429043 * y * z;
-        //     result += shCoefficients[ 6 ] * ( 0.743125 * z * z - 0.247708 );
-        //     result -= shCoefficients[ 7 ] * 2.0 * 0.429043 * x * z;
-        //     result += shCoefficients[ 8 ] * 0.429043 * ( x * x - y * y );
-        //    
-        //     return result;
-        // }
-        //
-        // vec3 inverseTransformDirection( in vec3 normal, in mat4 matrix ) {
-        //
-        //     // matrix is assumed to be orthogonal
-        //    
-        //     return normalize( ( vec4( normal, 0.0 ) * matrix ).xyz );
-        //
-        // }
-        //
-        // vec3 lpv_getLightProbeIrradiance( const in vec3 lightProbe[ 9 ], const in vec3 normal ) {
-        //
-        //     vec3 worldNormal = inverseTransformDirection( normal, viewMatrix );
-        //
-        //     vec3 irradiance = shGetIrradianceAt( worldNormal, lightProbe );
-        //
-        //     return irradiance;
-        //
-        // }
-`;
-
-const ACCUMULATION = `
-    irradiance += getLightProbeIrradiance(lpv_values, geometryNormal )* lpv_f_intensity;
-`;
-
-
-/**
- *
- * @param {string} source
- */
-export function build_fragment_shader(source) {
-
-    let result = PREAMBLE + source;
-
-
-    result = insert_after(
-        result,
-        '#include <lights_fragment_begin>',
-        ACCUMULATION
-    );
-
-    return result;
-}

package/src/engine/graphics/sh3/gi/material/build_vertex_shader.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"build_vertex_shader.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/gi/material/build_vertex_shader.js"],"names":[],"mappings":"AA+OA;;;GAGG;AACH,4CAFW,MAAM,UAchB"}

package/src/engine/graphics/sh3/gi/material/build_vertex_shader.js

@@ -1,256 +0,0 @@
-import { insert_after } from "../../../../../core/primitives/strings/insert_after.js";
-
-const PREAMBLE = `
-    #ifndef LPV_SHADER_CHUNK_PREAMBLE
-        #define LPV_SHADER_CHUNK_PREAMBLE
-        
-        precision highp usampler2D;
-        precision highp usampler3D;
-        
-        uniform usampler2D lpv_t_mesh_vertices;
-        uniform usampler2D lpv_t_mesh_neighbours;
-        uniform usampler3D lpv_t_mesh_lookup;
-        
-        uniform sampler2D lpv_t_probe_positions;
-        uniform sampler2D lpv_t_probe_data;
-        
-        uniform uint lpv_u_mesh_tet_count;
-        uniform vec3 lpv_v3_bounds_min;
-        uniform vec3 lpv_v3_bounds_max;
-        
-        varying vec3[9] lpv_values;
-        
-        #define SEARCH_STEP_LIMIT 64u
-        #define INVALID_TET 1073741823u
-        
-        ivec2 lpv_index_to_256_coordinate(uint index){
-        
-            uint pixel_x = index % 256u;
-            uint pixel_y = index / 256u; 
-            
-            return ivec2(int(pixel_x), int(pixel_y));
-        }
-        
-        uvec4 lpv_mesh_getVertices(uint tet_index){
-            ivec2 p = lpv_index_to_256_coordinate(tet_index);
-                    
-            return texelFetch(lpv_t_mesh_vertices, p, 0);
-            
-        }
-        
-        uvec4 lpv_mesh_getNeighbours(uint tet_index){
-            ivec2 p = lpv_index_to_256_coordinate(tet_index);
-            
-            return texelFetch(lpv_t_mesh_neighbours, p, 0);
-        }
-        
-        vec3[9] lpv_probe_getData(uint probe_index){
-            
-            
-            int slot = int( probe_index % 256u );
-            int column = int(probe_index / 256u); 
-            
-            int offset_x = int(slot * 9);
-            
-            vec3[9] result;
-            
-            for(int i=0; i<9; i++){
-                result[i] = texelFetch(lpv_t_probe_data, ivec2(offset_x + i, column),0).rgb;
-            }
-            
-            return result;
-        }
-        
-        vec3 lpv_probe_getPosition(uint probe_index){
-            return texelFetch(lpv_t_probe_positions, lpv_index_to_256_coordinate(probe_index), 0).rgb;
-        }
-        
-        mat3 lpv_mesh_makeMatrix(vec3 p0, vec3 p1, vec3 p2, vec3 p3){
-        
-            return inverse(
-                mat3(
-                    p0 - p3,
-                    p1 - p3,
-                    p2 - p3
-                )
-            );
-        
-        }
-        
-        vec4 lpv_mesh_getBarycentricCoordinates(uint tet_index, vec3 position){
-            uvec4 vertices = lpv_mesh_getVertices(tet_index);
-               
-            vec3 p0 = lpv_probe_getPosition(vertices[0]);
-            vec3 p1 = lpv_probe_getPosition(vertices[1]);
-            vec3 p2 = lpv_probe_getPosition(vertices[2]);
-            vec3 p3 = lpv_probe_getPosition(vertices[3]);
-        
-            mat3 matrix = lpv_mesh_makeMatrix(p0, p1, p2, p3);
-            
-            vec3 mult = matrix * (position - p3);
-            
-            return vec4(mult, 1.0 - mult.x - mult.y - mult.z);
-        }
-        
-        void lpv_getTetrahedron(
-            in vec3 position ,
-            in uint tet_guess,
-            out uint tet_index,
-            out vec4 weights
-        ){
-            uint came_from = 0u;
-            
-            tet_index = tet_guess;
-                    
-            for(uint i=0u; i < SEARCH_STEP_LIMIT; i++){
-            
-                weights = lpv_mesh_getBarycentricCoordinates(tet_index, position);
-            
-                if(weights.x >= 0.0 && weights.y >= 0.0 && weights.z >= 0.0 && weights.w >= 0.0){
-                    return;
-                }
-                
-                uvec4 neighbors = lpv_mesh_getNeighbours(tet_index);
-                
-                uint next_tet;
-                
-                // Otherwise find the smallest barycentric coord and move in that direction
-                if (weights.x < weights.y && weights.x < weights.z && weights.x < weights.w) {
-                    next_tet = neighbors[0];
-                } else if (weights.y < weights.z && weights.y < weights.w) {
-                    next_tet = neighbors[1];
-                } else if (weights.z < weights.w) {
-                    next_tet = neighbors[2];
-                } else {
-                    next_tet = neighbors[3];
-                }
-                
-                if(next_tet == came_from){
-                    // numerical instability caused us to oscillate on the boundary
-                    return;
-                }
-                
-                came_from = tet_index;  
-                tet_index = next_tet;
-            }
-        
-            tet_index = INVALID_TET;
-        }
-        
-        vec3[9] lpv_interpolate_probes(vec4 weights, uint tet_index){
-            
-            uvec4 vertices = lpv_mesh_getVertices(tet_index);
-            
-            vec3[9] probe0 = lpv_probe_getData(vertices[0]);
-            vec3[9] probe1 = lpv_probe_getData(vertices[1]);
-            vec3[9] probe2 = lpv_probe_getData(vertices[2]);
-            vec3[9] probe3 = lpv_probe_getData(vertices[3]);
-        
-            vec3[9] result;
-            
-            for(int i=0; i< 9; i++){
-                
-                result[i] = probe0[i]* weights[0]
-                            + probe1[i]* weights[1]
-                            + probe2[i]* weights[2]
-                            + probe3[i]* weights[3];
-            
-            }
-            
-            return result;
-        }
-        
-        vec4 lpv_bias_weight_by_normal(const in vec3 position, const in vec3 normal, const in uint tet_index, const in vec4 weights){
-            
-            uvec4 vertices = lpv_mesh_getVertices(tet_index);
-            
-            vec4 result = weights;
-            
-            for(int i = 0; i< 4; i++){
-                vec3 p = lpv_probe_getPosition(vertices[i]) - position;
-                
-                if(dot(normal, p) <= 0.0){
-                    // behind the surface, bias to 0
-                    result[i] = 0.0;
-                }
-            }
-            
-            float total_weight = result.x + result.y + result.z + result.w;
-            
-            // needs to add up to 1
-            result = result / total_weight;
-            
-            return result;
-        }
-        
-        uint lpv_guess_initial_tet(vec3 position){
-            
-            vec3 lpv_mesh_bounds_min = lpv_v3_bounds_min;
-            vec3 lpv_mesh_bounds_max = lpv_v3_bounds_max;
-            
-            vec3 lookup_coordinates = (position - lpv_mesh_bounds_min) / (lpv_mesh_bounds_max - lpv_mesh_bounds_min) ;
-           
-            return textureLod( lpv_t_mesh_lookup, lookup_coordinates, 0.0 ).r;
-        
-        }
-        
-    #endif
-`;
-
-
-const SHADER_CHUNK_ACCUMULATION = `
-    {
-        // lookup nearby tet
-        vec3 lpv_mesh_bounds_min = lpv_v3_bounds_min;
-        vec3 lpv_mesh_bounds_max = lpv_v3_bounds_max;
-        vec3 lookup_coordinates = (worldPosition.xyz - lpv_mesh_bounds_min) / (lpv_mesh_bounds_max - lpv_mesh_bounds_min) ;
-        
-        uint nearest_tet = lpv_guess_initial_tet( worldPosition.xyz );
-    
-        uint tet;
-        vec4 weights;
-        lpv_getTetrahedron(worldPosition.xyz, nearest_tet, tet, weights);
-        
-        
-        if(tet == INVALID_TET){
-            lpv_values = vec3[](
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0),
-                
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0),
-            
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0),
-                vec3(0.0, 0.0, 0.0)
-            );
-        }else{
-            
-	        vec3 worldNormal = normalize(inverseTransformDirection( normalize(normal), viewMatrix ));
-            
-            //weights = lpv_bias_weight_by_normal(worldPosition.xyz, normal, tet, weights);
-        
-            lpv_values = lpv_interpolate_probes(weights, tet);
-        }
-    }
-`
-
-/**
- *
- * @param {string} source
- */
-export function build_vertex_shader(source) {
-    let result = PREAMBLE + source;
-
-    result = insert_after(
-        result,
-        '#include <fog_vertex>',
-        SHADER_CHUNK_ACCUMULATION
-    );
-
-    // console.log(result)
-
-    return result;
-}

package/src/engine/graphics/sh3/lpv/sh3_bake_depth.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"sh3_bake_depth.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/lpv/sh3_bake_depth.js"],"names":[],"mappings":"AAQA;;;;;;;;;;GAUG;AACH,uCATW,MAAM,EAAE,iBACR,MAAM,sCAEN,MAAM,cACN,MAAM,cACN,MAAM,cACN,MAAM,eACN,MAAM,QAsDhB"}

package/src/engine/graphics/sh3/path_tracer/build_bvh_nodes_distances_along_ray.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"build_bvh_nodes_distances_along_ray.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/path_tracer/build_bvh_nodes_distances_along_ray.js"],"names":[],"mappings":"AAKA;;;;;;;;;;;;GAYG;AACH,4DAXW,MAAM,EAAE,iBACR,MAAM,uBAEN,MAAM,EAAE,oBACR,MAAM,kBACN,MAAM,kBACN,MAAM,kBACN,MAAM,kBACN,MAAM,gBACN,MAAM,QA4BhB"}

package/src/engine/graphics/sh3/path_tracer/build_bvh_nodes_distances_sqr_to_point.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"build_bvh_nodes_distances_sqr_to_point.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/path_tracer/build_bvh_nodes_distances_sqr_to_point.js"],"names":[],"mappings":"AAKA;;;;;;;;;;;GAWG;AACH,+DAVW,MAAM,EAAE,iBACR,MAAM,uBAEN,MAAM,EAAE,oBACR,MAAM,kBACN,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,QA2BhB"}

package/src/engine/graphics/sh3/path_tracer/ray_reflect.d.ts

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

package/src/engine/graphics/sh3/path_tracer/ray_reflect.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"ray_reflect.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/path_tracer/ray_reflect.js"],"names":[],"mappings":""}

package/src/engine/graphics/sh3/path_tracer/ray_reflect.js

@@ -1,27 +0,0 @@
-import { v3_dot } from "../../../../core/geom/vec3/v3_dot.js";
-
-/**
- *
- * @param {number[]} out
- * @param {number} output_offset
- * @param {number[]} v
- * @param {number} input_offset
- * @param {number[]} n
- * @param {number} normal_offset
- */
-function ray_reflect(out, output_offset, v, input_offset, n, normal_offset) {
-    const input_x = v[0 + input_offset];
-    const input_y = v[1 + input_offset];
-    const input_z = v[2 + input_offset];
-
-    const nx = n[normal_offset + 0];
-    const ny = n[normal_offset + 1];
-    const nz = n[normal_offset + 2];
-
-    const dot = v3_dot(input_x, input_y, input_z, nx, ny, nz);
-
-    out[output_offset + 0] = input_x - 2 * dot * nx;
-    out[output_offset + 1] = input_y - 2 * dot * ny;
-    out[output_offset + 2] = input_z - 2 * dot * nz;
-
-}

package/src/engine/graphics/sh3/path_tracer/sort_bvh_nodes_by_distance_to_point.d.ts.map

@@ -1 +0,0 @@
-{"version":3,"file":"sort_bvh_nodes_by_distance_to_point.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/path_tracer/sort_bvh_nodes_by_distance_to_point.js"],"names":[],"mappings":"AAMA;;;;;;;;;GASG;AACH,yEAPW,MAAM,EAAE,GAAC,WAAW,oBACpB,MAAM,kBACN,MAAM,SACN,MAAM,SACN,MAAM,SACN,MAAM,QA+BhB"}