@woosh/meep-engine 2.109.16 → 2.109.19

package/src/engine/graphics/sh3/lpv/depth/octahedral/v2/visualise.frag.glsl

@@ -1,5 +1,3 @@
-
-
 vec3 inverseTransformDirection(in vec3 normal, in mat4 matrix) {
     // matrix is assumed to be orthogonal
     return normalize((vec4(normal, 0.0) * matrix).xyz);
@@ -25,11 +23,22 @@ ivec2 lpv_index_to_256_coordinate(uint index) {
     return ivec2(int(pixel_x), int(pixel_y));
 }
 
+float sign_not_zero(float x) {
+    return x >= 0.0 ? 1.0 : -1.0;
+}
+
+vec2 sign_not_zero(vec2 x) {
+    return vec2(
+        sign_not_zero(x.x),
+        sign_not_zero(x.y)
+    );
+}
+
 vec3 sign_not_zero(vec3 x) {
     return vec3(
-        x.x >= 0.0 ? 1.0 : -1.0,
-        x.y >= 0.0 ? 1.0 : -1.0,
-        x.z >= 0.0 ? 1.0 : -1.0
+        sign_not_zero(x.x),
+        sign_not_zero(x.y),
+        sign_not_zero(x.z)
     );
 }
 
@@ -56,23 +65,45 @@ vec4 quadBlendWieghts(vec2 coords)
     return res;
 }
 
-
-vec2 VecToSphereOct(vec3 direction)
+//vec2 VecToSphereOct(vec3 v)
+//{
+//    float l1norm = abs(v.x) + abs(v.y) + abs(v.z);
+//
+//    vec2 result = v.xz / l1norm;
+//
+//    if (v.y < 0.0) {
+//        result = (1.0 - abs(result.yx)) * sign_not_zero(result.xy);
+//    }
+//
+//    return result;
+//}
+
+vec3 OctaSphereEnc(vec2 coord)
 {
-    vec3 octant = sign_not_zero(direction);
-
-    //  |x| + |y| + |z| = 1
-    float sum = dot(direction, octant);
-    vec3 octahedron = direction / sum;
+    coord = (coord - 0.5) * 2.0;
+    vec3 position = vec3(coord.x, 0.0, coord.y);
+    vec2 absolute = abs(position.xz);
+    position.y = 1.0 - absolute.x - absolute.y;
 
-    if (direction.y < 0.0)
+    if (position.y < 0.0)
     {
-        vec3 absolute = abs(octahedron);
+        position.xz = sign(position.xz) * vec2(1.0 - absolute.y, 1.0 - absolute.x);
+    }
+
+    return position;
+}
+
+vec2 VecToSphereOct(vec3 v)
+{
+    float l1norm = abs(v.x) + abs(v.y) + abs(v.z);
+
+    vec2 result = v.xz / l1norm;
 
-        octahedron.xz = octant.xz * vec2(1.0f - absolute.z, 1.0f - absolute.x);
+    if (v.y < 0.0) {
+        result = (1.0 - abs(result.yx)) * sign_not_zero(result.xy);
     }
 
-    return octahedron.xz;
+    return result;
 }
 
 float SampleBlended(sampler2D tex, vec2 uv0, vec2 uv1, vec2 uv2, vec4 weights) {
@@ -129,22 +160,245 @@ float lpv_probe_getDepth(uint probe_index, vec3 direction) {
     );
 }
 
-vec3 lpv_probe_getPosition(uint probe_index) {
-    return texelFetch(lpv_t_probe_positions, lpv_index_to_256_coordinate(probe_index), 0).rgb;
+float lpv_bilinear_lerp(float v00, float v01, float v10, float v11, vec2 fraction) {
+
+    float x0 = mix(v00, v01, fraction.x);
+    float x1 = mix(v10, v11, fraction.x);
+
+    return mix(x0, x1, fraction.y);
+}
+
+vec2 lpv_sample_bilinear(sampler2D tex, ivec2 texel_position, vec2 fraction) {
+
+    float texel_00 = texelFetch(tex, texel_position, 0).r;
+    float texel_01 = texelFetch(tex, texel_position + ivec2(1, 0), 0).r;
+    float texel_10 = texelFetch(tex, texel_position + ivec2(0, 1), 0).r;
+    float texel_11 = texelFetch(tex, texel_position + ivec2(1, 1), 0).r;
+
+    return vec2(
+        lpv_bilinear_lerp(
+            texel_00, texel_01,
+            texel_10, texel_11,
+            fraction
+        ),
+        lpv_bilinear_lerp(
+            texel_00 * texel_00, texel_01 * texel_01,
+            texel_10 * texel_10, texel_11 * texel_11,
+            fraction
+        )
+    );
+}
+
+ivec2 wrapOctahedralTexelCoordinatesNone(in ivec2 texelCoord, in int texture_size) {
+    return texelCoord;
+}
+
+ivec2 wrapOctahedralTexelCoordinatesV0(in ivec2 texelCoord, in int texture_size) {
+    //    return texelCoord;
+    // Decrement the texture size by 1 to get the maximum valid index
+    int maxIndex = texture_size - 1;
+
+    // Wrap the texel coordinates
+    ivec2 wrappedCoord = texelCoord;
+
+    // Wrap the X coordinate
+    if (wrappedCoord.x > maxIndex) {
+        wrappedCoord.x = maxIndex;
+        wrappedCoord.y += 1;
+    }
+
+    // Wrap the Y coordinate
+    if (wrappedCoord.y > maxIndex) {
+        wrappedCoord.y = maxIndex - wrappedCoord.y; // Wrap from top to bottom
+    }
+
+    return wrappedCoord;
+}
+
+ivec2 wrapOctahedralTexelCoordinates(const in ivec2 texel, const in int texture_size) {
+    ivec2 wrapped = ((texel % texture_size) + texture_size) % texture_size;
+
+    int fx = (abs(texel.x / texture_size) + int(texel.x < 0));
+    int fy = (abs(texel.y / texture_size) + int(texel.y < 0));
+
+    if (((fx ^ fy) & 1) != 0) {
+        return (texture_size - (wrapped + ivec2(1)));
+        return wrapped;
+    } else {
+        return wrapped;
+    }
+}
+
+
+vec2 lpv_probe_getDepthBilinear(uint probe_index, vec3 direction) {
+    // get offset
+    int depth_tile_resolution = int(lpv_u_probe_depth_resolution);
+    ivec2 tile_resolution = ivec2(depth_tile_resolution);
+    const ivec2 atlas_size = ivec2(4096);
+
+    int tiles_per_row = atlas_size.x / depth_tile_resolution;
+
+    int tile_x = int(probe_index) % tiles_per_row;
+    int tile_y = int(probe_index) / tiles_per_row;
+
+    ivec2 tile_offset = ivec2(
+        tile_x * depth_tile_resolution,
+        tile_y * depth_tile_resolution
+    );
+
+    // convert direction to UV
+    vec2 octahedral_uv = clamp(VecToSphereOct(direction) * 0.5 + 0.5, 0.0, 1.0);
+    vec2 grid = octahedral_uv * vec2(depth_tile_resolution) - 0.5;
+
+    vec2 gridFrac = fract(grid);
+
+    ivec2 texel_position = ivec2(floor(grid));
+
+    ivec2 tile_p_00;
+    ivec2 tile_p_01;
+    ivec2 tile_p_10;
+    ivec2 tile_p_11;
+
+    tile_p_00 = wrapOctahedralTexelCoordinates(texel_position , depth_tile_resolution);
+    tile_p_01 = wrapOctahedralTexelCoordinates(texel_position + ivec2(1, 0), depth_tile_resolution);
+    tile_p_10 = wrapOctahedralTexelCoordinates(texel_position + ivec2(0, 1), depth_tile_resolution);
+    tile_p_11 = wrapOctahedralTexelCoordinates(texel_position + ivec2(1, 1), depth_tile_resolution);
+
+    float texel_00 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_00, 0).r;
+    float texel_01 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_01, 0).r;
+    float texel_10 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_10, 0).r;
+    float texel_11 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_11, 0).r;
+
+    return vec2(
+        lpv_bilinear_lerp(
+            texel_00, texel_01,
+            texel_10, texel_11,
+            gridFrac
+        ),
+        lpv_bilinear_lerp(
+            texel_00 * texel_00, texel_01 * texel_01,
+            texel_10 * texel_10, texel_11 * texel_11,
+            gridFrac
+        )
+    );
+}
+
+vec2 lpv_probe_getDepthTriangular(uint probe_index, vec3 direction) {
+    // get offset
+    uint depth_tile_resolution = lpv_u_probe_depth_resolution;
+    uvec2 atlas_size = uvec2(4096u);
+
+    uint tiles_per_row = atlas_size.x / depth_tile_resolution;
+
+    uint tile_x = probe_index % tiles_per_row;
+    uint tile_y = probe_index / tiles_per_row;
+
+    vec2 tile_offset = vec2(
+        tile_x * depth_tile_resolution,
+        tile_y * depth_tile_resolution
+    );
+
+    // convert direction to UV
+    vec2 octahedral_uv = clamp(VecToSphereOct(direction) * 0.5 + 0.5, 0.0, 1.0);
+    vec2 grid = octahedral_uv * vec2(depth_tile_resolution - 1u);
+
+    vec2 gridFrac = fract(grid);
+    vec2 gridFloor = floor(grid);
+
+    vec4 weights = quadBlendWieghts(gridFrac);
+
+    //3 nearest frames
+    vec2 frame0 = gridFloor;
+    vec2 frame1 = gridFloor + mix(vec2(0, 1), vec2(1, 0), weights.w);
+    vec2 frame2 = gridFloor + vec2(1.0, 1.0);
+
+    // move frames to atlas space
+    frame0 += tile_offset;
+    frame1 += tile_offset;
+    frame2 += tile_offset;
+
+    float samp0 = texelFetch(lpv_t_probe_depth, ivec2(frame0), 0).r;
+    float samp1 = texelFetch(lpv_t_probe_depth, ivec2(frame1), 0).r;
+    float samp2 = texelFetch(lpv_t_probe_depth, ivec2(frame2), 0).r;
+
+    float d0 = samp0 * weights.x;
+    float d1 = samp1 * weights.y;
+    float d2 = samp2 * weights.z;
+
+    float mean = d0 + d1 + d2;
+    float mean2 = samp0 * d0 + samp1 * d1 + samp2 * d2;
+
+    return vec2(mean, mean2);
 }
+vec3 lpv_DEBUG(uint probe_index, vec3 direction) {
+    // get offset
+    int depth_tile_resolution = int(lpv_u_probe_depth_resolution);
+    const ivec2 atlas_size = ivec2(4096);
 
-float lpv_probe_getVisibilityMask(vec3 position, uint probe_index) {
+    int tiles_per_row = atlas_size.x / depth_tile_resolution;
 
-    vec3 probe_position = lpv_probe_getPosition(probe_index);
+    int tile_x = int(probe_index) % tiles_per_row;
+    int tile_y = int(probe_index) / tiles_per_row;
 
-    vec3 local_probe_offset = probe_position - position;
+    ivec2 tile_offset = ivec2(
+        tile_x * depth_tile_resolution,
+        tile_y * depth_tile_resolution
+    );
 
-    vec3 direction = normalize(local_probe_offset);
+    // convert direction to UV
+    vec2 octahedral_uv = clamp(VecToSphereOct(direction) * 0.5 + 0.5, 0.0, 1.0);
+    vec2 grid = octahedral_uv * vec2(depth_tile_resolution);
 
-    float depth = lpv_probe_getDepth(probe_index, -direction);
+    vec2 gridFrac = fract(grid);
+
+    ivec2 texel_position = ivec2(floor(grid));
+
+    ivec2 tile_p_01 = wrapOctahedralTexelCoordinates(texel_position + ivec2(1, 0), depth_tile_resolution);
+    ivec2 tile_p_10 = wrapOctahedralTexelCoordinates(texel_position + ivec2(0, 1), depth_tile_resolution);
+    ivec2 tile_p_11 = wrapOctahedralTexelCoordinates(texel_position + ivec2(1, 1), depth_tile_resolution);
+
+    return mix(vec3(
+                   tile_p_01.x == depth_tile_resolution ? 1.0 : 0.0,
+                   tile_p_10.y == depth_tile_resolution ? 1.0 : 0.0,
+                   0.0
+               ), vec3(octahedral_uv, 0.0), 0.0);
+
+    float texel_00 = texelFetch(lpv_t_probe_depth, tile_offset + texel_position, 0).r;
+    float texel_01 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_01, 0).r;
+    float texel_10 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_10, 0).r;
+    float texel_11 = texelFetch(lpv_t_probe_depth, tile_offset + tile_p_11, 0).r;
+
+
+    vec3 dir_00 = OctaSphereEnc(vec2(texel_position) / float(depth_tile_resolution));
+    vec3 dir_01 = OctaSphereEnc(vec2(tile_p_01) / float(depth_tile_resolution));
+    vec3 dir_10 = OctaSphereEnc(vec2(tile_p_10) / float(depth_tile_resolution));
+    vec3 dir_11 = OctaSphereEnc(vec2(tile_p_11) / float(depth_tile_resolution));
+
+
+    vec3 dir = vec3(
+        lpv_bilinear_lerp(
+            dir_00.x, dir_01.x,
+            dir_10.x, dir_11.x,
+            gridFrac
+        ), lpv_bilinear_lerp(
+            dir_00.y, dir_01.y,
+            dir_10.y, dir_11.y,
+            gridFrac
+        ), lpv_bilinear_lerp(
+            dir_00.z, dir_01.z,
+            dir_10.z, dir_11.z,
+            gridFrac
+        )
+    ) * 0.5 + 0.5;
+
+    return dir;
+
+    return vec3(octahedral_uv, 1.0);
+}
 
-    //    return step(depth, length(local_probe_offset));
-    return step(depth, 20.0);
+
+vec3 lpv_probe_getPosition(uint probe_index) {
+    return texelFetch(lpv_t_probe_positions, lpv_index_to_256_coordinate(probe_index), 0).rgb;
 }
 
 varying vec3 vNormal;
@@ -154,8 +408,10 @@ void main() {
     vec3 normal = normalize(vNormal);
     vec3 worldNormal = normalize(inverseTransformDirection(normal, viewMatrix));
 
-    float depth = lpv_probe_getDepth(instance_index , worldNormal) / 7.0;
+    float depth = lpv_probe_getDepthBilinear(instance_index, worldNormal).r / 7.0;
+//    float depth = lpv_probe_getDepthTriangular(instance_index, worldNormal).r / 7.0;
 
+    //    gl_FragColor = vec4(lpv_DEBUG(instance_index, worldNormal), 1.0);
     gl_FragColor = vec4(vec3(depth), 1.0);
-//    gl_FragColor = vec4( 1.0);
+    //    gl_FragColor = vec4( 1.0);
 }

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

@@ -1 +1 @@
-{"version":3,"file":"probe_volume_to_textures.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/lpv/probe_volume_to_textures.js"],"names":[],"mappings":"AAGA;;;;;GAKG;AACH,6EAHW,MAAM,EAAE,eACR,MAAM,QAqChB;AAGD;;;;;;GAMG;AACH,2DALW,oBAAoB,UACpB,MAAM,EAAE,eACR,MAAM,eACN,MAAM,QAiBhB;qCApEoC,uCAAuC"}
+{"version":3,"file":"probe_volume_to_textures.d.ts","sourceRoot":"","sources":["../../../../../../src/engine/graphics/sh3/lpv/probe_volume_to_textures.js"],"names":[],"mappings":"AAGA;;;;;GAKG;AACH,6EAHW,MAAM,EAAE,eACR,MAAM,QAiChB;AAGD;;;;;;GAMG;AACH,2DALW,oBAAoB,UACpB,MAAM,EAAE,eACR,MAAM,eACN,MAAM,QAiBhB;qCAhEoC,uCAAuC"}

package/src/engine/graphics/sh3/lpv/probe_volume_to_textures.js

@@ -27,18 +27,14 @@ export function octahedral_depth_to_atlas(lpv, atlas_data, atlas_width) {
 
                 const tile_pixel_index = pixel_y * tile_resolution + pixel_x;
 
-                const source_pixel_index = source_offset + tile_pixel_index;
-                const source_pixel_address = source_pixel_index * 2;
+                const source_value = probe_depth[source_offset + tile_pixel_index];
 
                 const target_pixel_x = tile_x * tile_resolution + pixel_x;
                 const target_pixel_y = tile_y * tile_resolution + pixel_y;
 
                 const target_pixel_index = target_pixel_y * atlas_width + target_pixel_x;
 
-                const target_pixel_address = target_pixel_index * 2;
-
-                atlas_data[target_pixel_address] = probe_depth[source_pixel_address];
-                atlas_data[target_pixel_address + 1] = probe_depth[source_pixel_address + 1];
+                atlas_data[target_pixel_index] = source_value;
             }
         }
     }

package/src/engine/graphics/sh3/lpv/serialization/LightProbeVolumeSerializationAdapter.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"LightProbeVolumeSerializationAdapter.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/lpv/serialization/LightProbeVolumeSerializationAdapter.js"],"names":[],"mappings":"AAGA;IAEI,+BAAwB;IAGxB;;;;OAIG;IACH,8BAFW,gBAAgB,QAiB1B;IAED;;;;OAIG;IACH,gCAFW,gBAAgB,QAiB1B;CACJ;gDAnD+C,mEAAmE;iCAClF,wBAAwB"}
+{"version":3,"file":"LightProbeVolumeSerializationAdapter.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/lpv/serialization/LightProbeVolumeSerializationAdapter.js"],"names":[],"mappings":"AAGA;IAEI,+BAAwB;IAGxB;;;;OAIG;IACH,8BAFW,gBAAgB,QAiB1B;IAED;;;;OAIG;IACH,gCAFW,gBAAgB,QAkB1B;CACJ;gDApD+C,mEAAmE;iCAClF,wBAAwB"}

package/src/engine/graphics/sh3/lpv/serialization/LightProbeVolumeSerializationAdapter.js

@@ -45,6 +45,7 @@ export class LightProbeVolumeSerializationAdapter extends BinaryClassSerializati
 
         // depth
         const depth_resolution = buffer.readUint32();
+        value.depth_map_resolution = depth_resolution;
         buffer.readFloat32Array(value.depth, 0, depth_resolution * depth_resolution * count);
 
         value.mesh.deserialize(buffer);

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

@@ -464,7 +464,7 @@ function* render(target, pt, camera, scene, progress = { current: 0, total: 0 })
 
                     for (let i = 0; i < pixel_sample_count; i++) {
 
-                        ray.projectFromCamera(
+                        ray.projectFromPerspectiveCamera(
                             _x + random() * pixel_scale_x * 2,
                             _y + random() * pixel_scale_y * 2,
                             0.5,

package/src/engine/graphics/sh3/path_tracer/sampling/getBiasedNormalSample.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"getBiasedNormalSample.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/getBiasedNormalSample.js"],"names":[],"mappings":"AAEA;;;;;;;;;GASG;AACH,2CAPW,MAAM,EAAE,cACR,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,SACN,MAAM,gBACK,MAAM,QAsD3B"}
+{"version":3,"file":"getBiasedNormalSample.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/getBiasedNormalSample.js"],"names":[],"mappings":"AAGA;;;;;;;;;GASG;AACH,2CAPW,MAAM,EAAE,cACR,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,SACN,MAAM,gBACK,MAAM,QAsD3B"}

package/src/engine/graphics/sh3/path_tracer/sampling/getBiasedNormalSample.js

@@ -1,3 +1,4 @@
+import { v3_length } from "../../../../../core/geom/vec3/v3_length.js";
 import { PI2 } from "../../../../../core/math/PI2.js";
 
 /**
@@ -37,7 +38,7 @@ export function getBiasedNormalSample(
     }
 
     // normalize orthonormal vector
-    const o1_norm = 1 / Math.hypot(o1_x, o1_y, o1_z);
+    const o1_norm = 1 / v3_length(o1_x, o1_y, o1_z);
 
     o1_x *= o1_norm;
     o1_y *= o1_norm;

package/src/engine/graphics/sh3/path_tracer/sampling/getCosineDirection.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"getCosineDirection.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/getCosineDirection.js"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,wCANW,MAAM,EAAE,cACR,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,gBACK,MAAM,QAgC3B"}
+{"version":3,"file":"getCosineDirection.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/getCosineDirection.js"],"names":[],"mappings":"AAGA;;;;;;;GAOG;AACH,wCANW,MAAM,EAAE,cACR,MAAM,UACN,MAAM,EAAE,iBACR,MAAM,gBACK,MAAM,QAgC3B"}

package/src/engine/graphics/sh3/path_tracer/sampling/getCosineDirection.js

@@ -1,3 +1,4 @@
+import { v3_length } from "../../../../../core/geom/vec3/v3_length.js";
 import { PI2 } from "../../../../../core/math/PI2.js";
 
 /**
@@ -33,7 +34,7 @@ export function getCosineDirection(
     let out_dir_y = dir_y + k1;
     let out_dir_z = dir_z + r_y;
 
-    const len2 = 1 / Math.hypot(out_dir_x, out_dir_y, out_dir_z);
+    const len2 = 1 / v3_length(out_dir_x, out_dir_y, out_dir_z);
 
     out[out_offset] = out_dir_x * len2;
     out[out_offset + 1] = out_dir_y * len2;

package/src/engine/graphics/sh3/path_tracer/sampling/v3_orthonormal_matrix_from_normal.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"v3_orthonormal_matrix_from_normal.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/v3_orthonormal_matrix_from_normal.js"],"names":[],"mappings":"AAAA;;;;;;;GAOG;AACH,uDANW,MAAM,EAAE,cACR,MAAM,SACN,MAAM,SACN,MAAM,SACN,MAAM,QAsChB"}
+{"version":3,"file":"v3_orthonormal_matrix_from_normal.d.ts","sourceRoot":"","sources":["../../../../../../../src/engine/graphics/sh3/path_tracer/sampling/v3_orthonormal_matrix_from_normal.js"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,uDANW,MAAM,EAAE,cACR,MAAM,SACN,MAAM,SACN,MAAM,SACN,MAAM,QAsChB"}

package/src/engine/graphics/sh3/path_tracer/sampling/v3_orthonormal_matrix_from_normal.js

@@ -1,3 +1,5 @@
+import { v3_length } from "../../../../../core/geom/vec3/v3_length.js";
+
 /**
  * Build orthonormal matrix from direction vector
  * @param {number[]} out
@@ -21,7 +23,7 @@ export function v3_orthonormal_matrix_from_normal(out, out_offset, dir_x, dir_y,
     }
 
     // normalize orthonormal vector
-    const o1_norm = 1 / Math.hypot(o1_x, o1_y, o1_z);
+    const o1_norm = 1 / v3_length(o1_x, o1_y, o1_z);
 
     o1_x *= o1_norm;
     o1_y *= o1_norm;

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

@@ -22,6 +22,7 @@ import { kelvin_to_rgb } from "../../../core/color/kelvin/kelvin_to_rgb.js";
 import { compose_matrix4_array } from "../../../core/geom/3d/compose_matrix4_array.js";
 import { visualize_tetrahedral_mesh } from "../../../core/geom/3d/tetrahedra/visualize_tetrahedral_mesh.js";
 import Quaternion from "../../../core/geom/Quaternion.js";
+import { v3_length } from "../../../core/geom/vec3/v3_length.js";
 import Vector3 from "../../../core/geom/Vector3.js";
 import { max2 } from "../../../core/math/max2.js";
 import { min2 } from "../../../core/math/min2.js";
@@ -245,7 +246,7 @@ async function getVolume({
         engine,
         ecd,
         bounds: mesh_bounds,
-        density: 500
+        density: 5000
     });
 
     const buffer = new BinaryBuffer();
@@ -383,13 +384,13 @@ async function main(engine) {
     // const path = 'data/models/samples/cyberpunk_bike/scene.gltf';
     // const path = 'data/models/LowPolyTownshipSet/Town_Hall/model.gltf';
     // const path = 'data/models/sibenik/3-window-less/model.gltf';
-    // const path = 'data/models/samples/gi_box_01/model.glb';
-    // const path = 'data/models/sponza-pbr/gltf/sponza.glb';
-    const path = 'data/models/pica_pica/pica_pica.gltf';
+    const path = 'data/models/sponza-pbr/gltf/sponza.glb';
     // const path = 'data/models/samples/susanne.glb';
     // const path = 'data/models/samples/teapot.gltf';
     // const path = 'data/models/samples/salle_de_bain/model.glb';
     // const path = 'data/models/samples/conference/model-no-curtains.glb';
+    // const path = 'data/models/pica_pica/pica_pica.gltf';
+    // const path = 'data/models/samples/gi_box_01/model.glb';
 
     const mesh_asset = await engine.assetManager.promise(path, 'model/gltf+json');
     const gltf = mesh_asset.gltf;
@@ -445,7 +446,7 @@ async function main(engine) {
         graphics: engine.graphics
     });
 
-    const model_footprint = Math.hypot(mesh_bounds.getExtentsX(), mesh_bounds.getExtentsY(), mesh_bounds.getExtentsZ());
+    const model_footprint = v3_length(mesh_bounds.getExtentsX(), mesh_bounds.getExtentsY(), mesh_bounds.getExtentsZ());
 
 
     const camera = ecd.getAnyComponent(Camera);
@@ -499,7 +500,7 @@ async function main(engine) {
 
 
     const probe_viz = volume.visualize_probes({ size: model_footprint * 0.002 });
-    probe_viz.build(ecd);
+    // probe_viz.build(ecd);
 
     const mesh_viz = visualize_tetrahedral_mesh({
         mesh: volume.mesh,
@@ -508,6 +509,7 @@ async function main(engine) {
     });
 
     const depth_viz = VLPDepthMapVisualisation.from(volume).build();
+    depth_viz.build(ecd);
 
     // volume.visualize_mesh({ ecd, opacity: 0.5 });