@meimu/mutracer 0.1.4 → 0.2.0

package/lib/renderer/wgsl/chunks/intersect.wgsl.js

@@ -1,42 +1,42 @@
 // Portions adapted from ray-tracing-renderer@0.10.15
 // (https://github.com/hoverinc/ray-tracing-renderer).
 // Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.
-import { uvProjectionConstants } from '../../MaterialTable.js';
+import { uvProjectionConstants } from '../../materialConstants.js';
 import { INSTANCE_VEC4_STRIDE } from '../../bvhAccel.js';
 import pdm_displacementWgsl from './pdm_displacement.wgsl.js';
 import { createIntersectionAlphaWgsl } from './intersectionAlpha.wgsl.js';
 const displacementWgsl = pdm_displacementWgsl;
 const flatIntersectCode = function () {
     return /* wgsl */ `
-let hit = fL(o, d, tri, ey, gP, tMax);
+let hit = ga(o, d, tri, eq, gY, tMax);
 if (hit.t > 0.0) {
 let e1 = tri.p1 - tri.p0;
 let e2 = tri.p2 - tri.p0;
-eD = normalize(cross(e1, e2));
-let kk = i32(ae(bitcast<i32>(ou.z)).x);
-let iN = materialsBuffer[kk].doubleSided > 0.5;
-if (!doubleSided && !iN && dot(eD, d) > 0.0) {
-eD = vec3f(0.0);
+eL = normalize(cross(e1, e2));
+let kP = i32(ah(bitcast<i32>(pi.z)).x);
+let jd = materialsBuffer[kP].doubleSided > 0.5;
+if (!doubleSided && !jd && dot(eL, d) > 0.0) {
+eL = vec3f(0.0);
 continue;
 }
 
 tMax = hit.t;
-rb = true;
-ox = true;
-eY = o + tMax * d;
-aB = hit.barycentric;
+sj = true;
+pl = true;
+fi = o + tMax * d;
+aC = hit.barycentric;
 }
 `;
 };
 function getPDMDisplacmentCode(anyHit) {
     return /*wgsl */ `
-let kX = pc(
+let lH = pY(
 o, d, tMax,
 tri.p0, tri.p1, tri.p2,
 n0, n1, n2,
 uv0, uv1, uv2,
 displacementMapIndex,
-kJ(displacementMapIndex),
+jO(displacementMapIndex),
 minH, maxH,
 hScale,
 matData.displacementBias,
@@ -44,17 +44,17 @@ matData.displacementStepSize,
 i32(matData.displacementUvProjection),
 displacementUvProjectionScale,
 ${anyHit ? 'true' : 'false'},
-lD
+jA
 );
 `;
 }
 const displacementIntersectCode = function (anyHit, hasAnisotropyDirections) {
     return /* wgsl */ `
 
-let hScale = ou.w;
-let tA = bitcast<i32>(r3.w);
+let hScale = pi.w;
+let uQ = bitcast<i32>(r3.w);
 
-let materialIndex = i32(ae(bitcast<i32>(ou.z)).x);
+let materialIndex = i32(ah(bitcast<i32>(pi.z)).x);
 let matData = materialsBuffer[materialIndex];
 let n0 = normalize(attrBuffers[index.x * ATTR_STRIDE + NORMAL_OFF]${hasAnisotropyDirections ? '.xyz' : ''});
 let n1 = normalize(attrBuffers[index.y * ATTR_STRIDE + NORMAL_OFF]${hasAnisotropyDirections ? '.xyz' : ''});
@@ -64,9 +64,9 @@ tri.p0 += n0 * matData.displacementBias;
 tri.p1 += n1 * matData.displacementBias;
 tri.p2 += n2 * matData.displacementBias;
 }
-if (abs(hScale) > 0.0 && tA > 0) {
-let minH = f32(tA >> 8) / 255.0;
-let maxH = f32(tA & 0xff) / 255.0;
+if (abs(hScale) > 0.0 && uQ > 0) {
+let minH = f32(uQ >> 8) / 255.0;
+let maxH = f32(uQ & 0xff) / 255.0;
 let displacementUvProjectionScale = vec2f(matData.displacementUvProjectionScale);
 let displacementMapIndex = i32(matData.displacementMapIndex);
 
@@ -77,18 +77,18 @@ var uv2 = attrBuffers[index.z * ATTR_STRIDE + UV_OFF].xy;
 
 ${getPDMDisplacmentCode(anyHit)}
 
-if (kX.t < tMax) {
-eD = kX.normal;
+if (lH.t < tMax) {
+eL = lH.normal;
 
 
-tMax = kX.t;
+tMax = lH.t;
 
-aB = kX.barycentric;
-eY = o + tMax * d;
+aC = lH.barycentric;
+fi = o + tMax * d;
 
-bU = true;
-rb = true;
-ox = true;
+bZ = true;
+sj = true;
+pl = true;
 }
 }
 else {
@@ -97,7 +97,7 @@ ${flatIntersectCode()}
 `;
 };
 const cwbvhTraversalWgsl = /* wgsl */ `
-fn bp(value: u32) -> vec4f {
+fn bA(value: u32) -> vec4f {
 return vec4f(
 f32(value & 0xffu),
 f32((value >> 8u) & 0xffu),
@@ -112,118 +112,118 @@ f32(value >> 24u)
 
 
 
-fn kx(value: u32) -> u32 {
+fn lb(value: u32) -> u32 {
 return (bitcast<u32>(f32(value)) >> 23u) - 127u;
 }
 
-fn nZ(value: u32) -> u32 {
+fn oT(value: u32) -> u32 {
 return (bitcast<u32>(f32(value)) >> 23u) - 127u;
 }
 
-fn ov(
+fn pj(
 node: i32,
-tk: i32,
+uJ: i32,
 n0: vec4f,
 o: vec3f,
-gC: vec3f,
+gr: vec3f,
 tMax: f32
 ) -> u32 {
-let i = node + tk;
+let i = node + uJ;
 let n2 = bvhBuffer[i + 2];
 let n3 = bvhBuffer[i + 3];
 let n4 = bvhBuffer[i + 4];
 let header = bitcast<u32>(n0.w);
-let rG = bitcast<i32>((header & 0xffu) << 24u) >> 24;
-let rH = bitcast<i32>(((header >> 8u) & 0xffu) << 24u) >> 24;
-let rI = bitcast<i32>(((header >> 16u) & 0xffu) << 24u) >> 24;
-let hJ = vec3f(
-bitcast<f32>(u32(rG + 127) << 23u),
-bitcast<f32>(u32(rH + 127) << 23u),
-bitcast<f32>(u32(rI + 127) << 23u)
+let sH = bitcast<i32>((header & 0xffu) << 24u) >> 24;
+let sI = bitcast<i32>(((header >> 8u) & 0xffu) << 24u) >> 24;
+let sJ = bitcast<i32>(((header >> 16u) & 0xffu) << 24u) >> 24;
+let hZ = vec3f(
+bitcast<f32>(u32(sH + 127) << 23u),
+bitcast<f32>(u32(sI + 127) << 23u),
+bitcast<f32>(u32(sJ + 127) << 23u)
 );
-let bJ = eZ(gC, o);
-let eI = (n0.xyz - o) * gC;
-let dc = hJ * gC;
-let wy = bp(bitcast<u32>(n2.x));
-let wB = bp(bitcast<u32>(n2.y));
-let wC = bp(bitcast<u32>(n2.z));
-let wE = bp(bitcast<u32>(n2.w));
-let wF = bp(bitcast<u32>(n3.x));
-let wG = bp(bitcast<u32>(n3.y));
-let ws = bp(bitcast<u32>(n3.z));
-let wt = bp(bitcast<u32>(n3.w));
-let wu = bp(bitcast<u32>(n4.x));
-let wv = bp(bitcast<u32>(n4.y));
-let ww = bp(bitcast<u32>(n4.z));
-let wx = bp(bitcast<u32>(n4.w));
-
-let vp = wy * dc.x + eI.x;
-let vs = wC * dc.y + eI.y;
-let vu = wF * dc.z + eI.z;
-let vj = ws * dc.x + eI.x;
-let vl = wu * dc.y + eI.y;
-let vn = ww * dc.z + eI.z;
-let wI = max(
-max((min(vp, vj) - vec4f(bJ.x)), (min(vs, vl) - vec4f(bJ.y))),
-(min(vu, vn) - vec4f(bJ.z))
+let padding = fk(gr, o);
+let eS = (n0.xyz - o) * gr;
+let dq = hZ * gr;
+
+
+
+let vc = bitcast<vec4u>(n2);
+let vd = bitcast<vec4u>(n3);
+let ve = bitcast<vec4u>(n4);
+let negative = gr < vec3f(0.0);
+let wn = select(vc.xy, vd.zw, negative.x);
+let xt = select(vd.zw, vc.xy, negative.x);
+let wo = select(vc.zw, ve.xy, negative.y);
+let xu = select(ve.xy, vc.zw, negative.y);
+let wp = select(vd.xy, ve.zw, negative.z);
+let xv = select(ve.zw, vd.xy, negative.z);
+
+let yl = max(
+max(
+bA(wn.x) * dq.x + eS.x - vec4f(padding.x),
+bA(wo.x) * dq.y + eS.y - vec4f(padding.y)
+),
+bA(wp.x) * dq.z + eS.z - vec4f(padding.z)
 );
-let wn = min(
-min((max(vp, vj) + vec4f(bJ.x)), (max(vs, vl) + vec4f(bJ.y))),
-(max(vu, vn) + vec4f(bJ.z))
+let yg = min(
+min(
+bA(xt.x) * dq.x + eS.x + vec4f(padding.x),
+bA(xu.x) * dq.y + eS.y + vec4f(padding.y)
+),
+bA(xv.x) * dq.z + eS.z + vec4f(padding.z)
 );
-
-let vq = wB * dc.x + eI.x;
-let vt = wE * dc.y + eI.y;
-let vw = wG * dc.z + eI.z;
-let vk = wt * dc.x + eI.x;
-let vm = wv * dc.y + eI.y;
-let vo = wx * dc.z + eI.z;
-let wJ = max(
-max((min(vq, vk) - vec4f(bJ.x)), (min(vt, vm) - vec4f(bJ.y))),
-(min(vw, vo) - vec4f(bJ.z))
+let ym = max(
+max(
+bA(wn.y) * dq.x + eS.x - vec4f(padding.x),
+bA(wo.y) * dq.y + eS.y - vec4f(padding.y)
+),
+bA(wp.y) * dq.z + eS.z - vec4f(padding.z)
 );
-let wo = min(
-min((max(vq, vk) + vec4f(bJ.x)), (max(vt, vm) + vec4f(bJ.y))),
-(max(vw, vo) + vec4f(bJ.z))
+let yh = min(
+min(
+bA(xt.y) * dq.x + eS.x + vec4f(padding.x),
+bA(xu.y) * dq.y + eS.y + vec4f(padding.y)
+),
+bA(xv.y) * dq.z + eS.z + vec4f(padding.z)
 );
 
-let uK = (wI <= min(wn, vec4f(tMax))) & (wn > vec4f(0.0));
-let uL = (wJ <= min(wo, vec4f(tMax))) & (wo > vec4f(0.0));
-return select(0u, 1u, uK.x)
-| select(0u, 2u, uK.y)
-| select(0u, 4u, uK.z)
-| select(0u, 8u, uK.w)
-| select(0u, 16u, uL.x)
-| select(0u, 32u, uL.y)
-| select(0u, 64u, uL.z)
-| select(0u, 128u, uL.w);
+let wl = (yl <= min(yg, vec4f(tMax))) & (yg > vec4f(0.0));
+let wm = (ym <= min(yh, vec4f(tMax))) & (yh > vec4f(0.0));
+return select(0u, 1u, wl.x)
+| select(0u, 2u, wl.y)
+| select(0u, 4u, wl.z)
+| select(0u, 8u, wl.w)
+| select(0u, 16u, wm.x)
+| select(0u, 32u, wm.y)
+| select(0u, 64u, wm.z)
+| select(0u, 128u, wm.w);
 }
 
 `;
 function getBVHTraverseCode(handleHitCode, hasDisplacement, anyHit, blasLayoutMask, hasAnisotropyDirections, skipTriangleCode = '') {
     const intersectTriangles = /* wgsl */ `
-let ou = instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE}];
-let vertexOffset = vec3i(bitcast<i32>(ou.y));
-for (var k = 0; k < bh; k++) {
+let pi = instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE}];
+let vertexOffset = vec3i(bitcast<i32>(pi.y));
+for (var k = 0; k < primitiveCount; k++) {
 let triangleIndex = primitiveOffset + k;
 ${skipTriangleCode}
-if (bounces > 0 && triangleIndex == hO) {
+if (bounces > 0 && triangleIndex == fV) {
 continue;
 }
 
 let r3 = bvhBuffer[triangleIndex];
 let index = bitcast<vec3i>(r3.xyz) + vertexOffset;
-var tri = fA(
+var tri = fP(
 attrBuffers[index.x * ATTR_STRIDE + POSITION_OFF]${hasAnisotropyDirections ? '.xyz' : ''},
 attrBuffers[index.y * ATTR_STRIDE + POSITION_OFF]${hasAnisotropyDirections ? '.xyz' : ''},
 attrBuffers[index.z * ATTR_STRIDE + POSITION_OFF]${hasAnisotropyDirections ? '.xyz' : ''}
 );
-var ox = false;
+var pl = false;
 
 ${hasDisplacement
         ? displacementIntersectCode(anyHit, hasAnisotropyDirections)
         : flatIntersectCode()}
-if (ox) {
+if (pl) {
 ${handleHitCode}
 }
 }
@@ -231,285 +231,285 @@ ${handleHitCode}
     return /* wgsl */ `
 var d = ray.d;
 var o = ray.o;
-var mb: array<i32, TLAS_STACK_SIZE>;
-var jh = -1;
-var gC = fb(d);
-var bU = false;
-var rb = false;
-var eY: vec3f;
-var eD: vec3f;
-var aB: vec3f;
-
-
-if (cM(o, bvhBuffer[0].xyz, bvhBuffer[1].xyz, gC, tMax) > 0.0) {
-jh = 0;
-mb[0] = 0;
-}
-while (jh >= 0) {
-let sV = mb[jh];
-jh--;
+var jE = -1;
+var gr = fo(d);
+var bZ = false;
+var sj = false;
+var fi: vec3f;
+var eL: vec3f;
+var aC: vec3f;
+
+
+if (cU(o, bvhBuffer[0].xyz, bvhBuffer[1].xyz, gr, tMax) > 0.0) {
+jE = 0;
+mK[0] = 0;
+}
+while (jE >= 0) {
+let ul = mK[jE];
+jE--;
 o = ray.o;
 d = ray.d;
-gC = fb(d);
-
-let wf = bvhBuffer[sV];
-let eV = bitcast<i32>(wf.w);
-if (eV >= 0) {
-let left = sV + 2;
-let right = eV;
-let d0 = cM(
+gr = fo(d);
+
+let xW = bvhBuffer[ul];
+let ff = bitcast<i32>(xW.w);
+if (ff >= 0) {
+let left = ul + 2;
+let right = ff;
+let d0 = cU(
 o,
 bvhBuffer[left].xyz,
 bvhBuffer[left + 1].xyz,
-gC,
+gr,
 tMax
 );
-let d1 = cM(
+let d1 = cU(
 o,
 bvhBuffer[right].xyz,
 bvhBuffer[right + 1].xyz,
-gC,
+gr,
 tMax
 );
 
 if (d0 > d1) {
-if (d0 > 0.0) { jh++; mb[jh] = left; }
-if (d1 > 0.0) { jh++; mb[jh] = right; }
+if (d0 > 0.0) { jE++; mK[jE] = left; }
+if (d1 > 0.0) { jE++; mK[jE] = right; }
 }
 else {
-if (d1 > 0.0) { jh++; mb[jh] = right; }
-if (d0 > 0.0) { jh++; mb[jh] = left; }
+if (d1 > 0.0) { jE++; mK[jE] = right; }
+if (d0 > 0.0) { jE++; mK[jE] = left; }
 }
 }
 else {
-let qJ = -eV;
+let rQ = -ff;
 
 
 for (
-var mx = qJ - 1;
-mx >= 0;
-mx--
+var nq = rQ - 1;
+nq >= 0;
+nq--
 ) {
-let pI = bvhBuffer[sV + 2 + mx];
-let cf = bitcast<i32>(pI.x);
-let uU = bitcast<i32>(pI.y);
-let ou = instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE}];
+let qv = bvhBuffer[ul + 2 + nq];
+let bY = bitcast<i32>(qv.x);
+let wL = bitcast<i32>(qv.y);
+let pi = instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE}];
 let w2l = mat4x4<f32>(
-instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE} + 1],
-instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE} + 2],
-instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE} + 3],
-instancesBuffer[cf * ${INSTANCE_VEC4_STRIDE} + 4]
+instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE} + 1],
+instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE} + 2],
+instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE} + 3],
+instancesBuffer[bY * ${INSTANCE_VEC4_STRIDE} + 4]
 );
 o = (w2l * vec4f(ray.o, 1.0)).xyz;
 d = (w2l * vec4f(ray.d, 0.0)).xyz;
-gC = fb(d);
-ey = maxDimension(abs(d));
-gP = mv(d, ey);
+gr = fo(d);
+eq = maxDimension(abs(d));
+gY = nk(d, eq);
 
-let fJ = bitcast<u32>(ou.x);
-let gR = fJ >> 30u;
-let blasOffset = i32(fJ & 0x3fffffffu);
-let bo = rtUniform.blasOffset;
+let fY = bitcast<u32>(pi.x);
+let he = fY >> 30u;
+let blasOffset = i32(fY & 0x3fffffffu);
+let bw = rtUniform.blasOffset;
 
 ${(blasLayoutMask & 8) !== 0 ? /* wgsl */ `
-if (gR == 3u) {
+if (he == 3u) {
 
 
 
 
-var sm: array<vec2u, BLAS_STACK_SIZE>;
-var pn = -1;
-var rX = 0u;
-var kZ = 0u;
+var qj = -1;
+var td = 0u;
+var lJ = 0u;
 
-var ep = 0x80000000u;
-let ut = select(0u, 4u, d.x < 0.0)
+var eD = 0x80000000u;
+let vU = select(0u, 4u, d.x < 0.0)
 | select(0u, 2u, d.y < 0.0)
 | select(0u, 1u, d.z < 0.0);
-let ql = 7u - ut;
+let rn = 7u - vU;
 
 
 
 
-while (rX < 4096u) {
-rX++;
+while (td < 4096u) {
+td++;
 
-let internalMask = ep & 0xffu;
-let qh =
-kx(ep >> 24u) + 24u;
-ep &= ~(1u << qh);
-if ((ep & 0xff000000u) != 0u) {
-pn++;
+let internalMask = eD & 0xffu;
+let ri =
+lb(eD >> 24u) + 24u;
+eD &= ~(1u << ri);
+if ((eD & 0xff000000u) != 0u) {
+qj++;
 
-sm[pn] = vec2u(kZ, ep);
+tD[qj] = vec2u(lJ, eD);
 }
 
-let slot = (qh - 24u) ^ ql;
-let ns = (1u << slot) - 1u;
-let relativeIndex = countOneBits(internalMask & ns);
-let sk = kZ + relativeIndex;
-let node = blasOffset + i32(sk * 5u);
-let jA = bvhBuffer[node + bo];
-let ue = ov(
-node, bo, jA, o, gC, tMax
+let slot = (ri - 24u) ^ rn;
+let om = (1u << slot) - 1u;
+let relativeIndex = countOneBits(internalMask & om);
+let tA = lJ + relativeIndex;
+let node = blasOffset + i32(tA * 5u);
+let jZ = bvhBuffer[node + bw];
+let vG = pj(
+node, bw, jZ, o, gr, tMax
 );
 
 
-ep = 0u;
-if (ue != 0u) {
-let fF = bvhBuffer[node + bo + 1];
-let rO = bitcast<u32>(jA.w) >> 24u;
-var tv = 0u;
-var kg = ue;
-while (kg != 0u) {
-let rC = kx(kg);
-kg &= ~(1u << rC);
-let sE = select(
-bitcast<u32>(fF.z),
-bitcast<u32>(fF.w),
-rC >= 4u
+eD = 0u;
+if (vG != 0u) {
+let fU = bvhBuffer[node + bw + 1];
+let sP = bitcast<u32>(jZ.w) >> 24u;
+var uO = 0u;
+var kM = vG;
+while (kM != 0u) {
+let sC = lb(kM);
+kM &= ~(1u << sC);
+let tT = select(
+bitcast<u32>(fU.z),
+bitcast<u32>(fU.w),
+sC >= 4u
 );
-let td =
-(sE >> ((rC & 3u) * 8u)) & 0xffu;
-let jC = td & 0x1fu;
-let uT = select(
-jC,
-jC ^ ql,
-jC >= 24u
+let uB =
+(tT >> ((sC & 3u) * 8u)) & 0xffu;
+let kb = uB & 0x1fu;
+let wK = select(
+kb,
+kb ^ rn,
+kb >= 24u
 );
-let uO = (td >> 5u) & 0x7u;
-tv |= uO << uT;
-}
-
-kZ = bitcast<u32>(fF.x);
-ep = (tv & 0xff000000u) | rO;
-let oQ = bitcast<u32>(fF.y);
-var oR = tv & 0x00ffffffu;
-while (oR != 0u) {
-let nT = nZ(oR);
-oR &= ~(1u << nT);
-let bh = 1;
-let primitiveOffset = bo + blasOffset
-+ i32(oQ + nT);
+let wt = (uB >> 5u) & 0x7u;
+uO |= wt << wK;
+}
+
+lJ = bitcast<u32>(fU.x);
+eD = (uO & 0xff000000u) | sP;
+let pM = bitcast<u32>(fU.y);
+var pN = uO & 0x00ffffffu;
+while (pN != 0u) {
+let oO = oT(pN);
+pN &= ~(1u << oO);
+let primitiveCount = 1;
+let primitiveOffset = bw + blasOffset
++ i32(pM + oO);
 ${intersectTriangles}
 }
 }
 
-if ((ep & 0xff000000u) != 0u) { continue; }
-if (pn >= 0) {
-let qm = sm[pn];
-kZ = qm.x;
-ep = qm.y;
-pn--;
+if ((eD & 0xff000000u) != 0u) { continue; }
+if (qj >= 0) {
+let ro = tD[qj];
+lJ = ro.x;
+eD = ro.y;
+qj--;
 }
 else { break; }
 }
 }
 ` : ''}
 
-var gt: array<i32, BLAS_STACK_SIZE>;
-var dS = -1;
-if (gR != 3u) {
-dS = 0;
+var eh = -1;
+if (he != 3u) {
+eh = 0;
 
-gt[0] = uU;
+gJ[0] = wL;
 }
 
-while (dS >= 0) {
+while (eh >= 0) {
 
-let node = gt[dS];
-dS--;
-var bh = 0;
+let node = gJ[eh];
+eh--;
+var primitiveCount = 0;
 var primitiveOffset = 0;
 
 ${(blasLayoutMask & 4) !== 0 ? /* wgsl */ `
-if (gR == 2u) {
-let i = node + bo;
-let uY = bvhBuffer[i];
-let uZ = bvhBuffer[i + 1];
-let va = bvhBuffer[i + 2];
-let vb = bvhBuffer[i + 3];
-bh = i32(bitcast<u32>(va.w));
-if (bh > 0) {
-primitiveOffset = bo + blasOffset
-+ i32(bitcast<u32>(vb.w));
+if (he == 2u) {
+let i = node + bw;
+let wP = bvhBuffer[i];
+let wQ = bvhBuffer[i + 1];
+let wR = bvhBuffer[i + 2];
+let wS = bvhBuffer[i + 3];
+primitiveCount = i32(bitcast<u32>(wR.w));
+if (primitiveCount > 0) {
+primitiveOffset = bw + blasOffset
++ i32(bitcast<u32>(wS.w));
 }
 else {
-let left = blasOffset + i32(bitcast<u32>(uY.w) * 4u);
-let right = blasOffset + i32(bitcast<u32>(uZ.w) * 4u);
-let oy = cM(
-o, uY.xyz, uZ.xyz, gC, tMax
+let left = blasOffset + i32(bitcast<u32>(wP.w) * 4u);
+let right = blasOffset + i32(bitcast<u32>(wQ.w) * 4u);
+let pm = cU(
+o, wP.xyz, wQ.xyz, gr, tMax
 );
-let nu = cM(
-o, va.xyz, vb.xyz, gC, tMax
+let oq = cU(
+o, wR.xyz, wS.xyz, gr, tMax
 );
-if (oy > nu) {
-if (oy > 0.0) {
-dS++; gt[dS] = left;
+if (pm > oq) {
+if (pm > 0.0) {
+eh++; gJ[eh] = left;
 }
-if (nu > 0.0) {
-dS++; gt[dS] = right;
+if (oq > 0.0) {
+eh++; gJ[eh] = right;
 }
 }
 else {
-if (nu > 0.0) {
-dS++; gt[dS] = right;
+if (oq > 0.0) {
+eh++; gJ[eh] = right;
 }
-if (oy > 0.0) {
-dS++; gt[dS] = left;
+if (pm > 0.0) {
+eh++; gJ[eh] = left;
 }
 }
 }
 }
 ` : ''}
 ${(blasLayoutMask & 2) !== 0 ? /* wgsl */ `
-if (gR == 1u) {
-let hW = bitcast<u32>(node);
-if ((hW & 0x80000000u) != 0u) {
-bh = i32((hW >> 23u) & 0xffu);
-primitiveOffset = bo + blasOffset
-+ i32(hW & 0x007fffffu);
+if (he == 1u) {
+let iw = bitcast<u32>(node);
+if ((iw & 0x80000000u) != 0u) {
+primitiveCount = i32((iw >> 23u) & 0xffu);
+primitiveOffset = bw + blasOffset
++ i32(iw & 0x007fffffu);
 }
 else {
-let i = node + bo;
-let ui = bvhBuffer[i];
-let ro = bvhBuffer[i + 1];
-let rp = bvhBuffer[i + 2];
-let rq = bvhBuffer[i + 3];
-let nodeMin = ui.xyz;
-let ec = vec3f(ui.w, ro.x, ro.y);
-let mh = array<u32, 4>(
-bitcast<u32>(ro.z), bitcast<u32>(ro.w),
-bitcast<u32>(rp.x), bitcast<u32>(rp.y)
+let i = node + bw;
+let vJ = bvhBuffer[i];
+let sz = bvhBuffer[i + 1];
+let qR = bvhBuffer[i + 2];
+let oQ = bvhBuffer[i + 3];
+let nodeMin = vJ.xyz;
+let er = vec3f(vJ.w, sz.x, sz.y);
+let mT = array<u32, 4>(
+bitcast<u32>(sz.z), bitcast<u32>(sz.w),
+bitcast<u32>(qR.x), bitcast<u32>(qR.y)
 );
-let rw = unpack4x8unorm(bitcast<u32>(rp.z));
-let rA = unpack4x8unorm(bitcast<u32>(rp.w));
-let rB = unpack4x8unorm(bitcast<u32>(rq.x));
-let rt = unpack4x8unorm(bitcast<u32>(rq.y));
-let ru = unpack4x8unorm(bitcast<u32>(rq.z));
-let rv = unpack4x8unorm(bitcast<u32>(rq.w));
-let tx0 = (nodeMin.x + rw * ec.x - o.x) * gC.x;
-let tx1 = (nodeMin.x + rt * ec.x - o.x) * gC.x;
-let ty0 = (nodeMin.y + rA * ec.y - o.y) * gC.y;
-let ty1 = (nodeMin.y + ru * ec.y - o.y) * gC.y;
-let tz0 = (nodeMin.z + rB * ec.z - o.z) * gC.z;
-let tz1 = (nodeMin.z + rv * ec.z - o.z) * gC.z;
-let bJ = eZ(gC, o);
-let mL = max(
-max((min(tx0, tx1) - vec4f(bJ.x)), (min(ty0, ty1) - vec4f(bJ.y))),
-(min(tz0, tz1) - vec4f(bJ.z))
+
+
+let negative = gr < vec3f(0.0);
+let wn = unpack4x8unorm(bitcast<u32>(select(qR.z, oQ.y, negative.x)));
+let xt = unpack4x8unorm(bitcast<u32>(select(oQ.y, qR.z, negative.x)));
+let wo = unpack4x8unorm(bitcast<u32>(select(qR.w, oQ.z, negative.y)));
+let xu = unpack4x8unorm(bitcast<u32>(select(oQ.z, qR.w, negative.y)));
+let wp = unpack4x8unorm(bitcast<u32>(select(oQ.x, oQ.w, negative.z)));
+let xv = unpack4x8unorm(bitcast<u32>(select(oQ.w, oQ.x, negative.z)));
+let padding = fk(gr, o);
+let nJ = max(
+max(
+(nodeMin.x + wn * er.x - o.x) * gr.x - vec4f(padding.x),
+(nodeMin.y + wo * er.y - o.y) * gr.y - vec4f(padding.y)
+),
+(nodeMin.z + wp * er.z - o.z) * gr.z - vec4f(padding.z)
 );
-let rU = min(
-min((max(tx0, tx1) + vec4f(bJ.x)), (max(ty0, ty1) + vec4f(bJ.y))),
-(max(tz0, tz1) + vec4f(bJ.z))
+let sY = min(
+min(
+(nodeMin.x + xt * er.x - o.x) * gr.x + vec4f(padding.x),
+(nodeMin.y + xu * er.y - o.y) * gr.y + vec4f(padding.y)
+),
+(nodeMin.z + xv * er.z - o.z) * gr.z + vec4f(padding.z)
 );
-var distance = select(rU, mL, mL > vec4f(0.0));
+var distance = select(sY, nJ, nJ > vec4f(0.0));
 distance = select(
 distance,
 vec4f(-1.0),
-mL > min(rU, vec4f(tMax))
+nJ > min(sY, vec4f(tMax))
 );
-var fd = array<f32, 4>(
+var fq = array<f32, 4>(
 distance.x,
 distance.y,
 distance.z,
@@ -518,83 +518,83 @@ distance.w
 
 
 for (var order = 0; order < 4; order++) {
-var dX = -1.0;
-var el = -1;
+var em = -1.0;
+var eB = -1;
 for (var child = 0; child < 4; child++) {
 if (
-mh[child] != 0u &&
-fd[child] > dX
+mT[child] != 0u &&
+fq[child] > em
 ) {
-dX = fd[child];
-el = child;
+em = fq[child];
+eB = child;
 }
 }
-if (el >= 0 && dX > 0.0) {
-let nD = mh[el];
-var jI = bitcast<i32>(nD);
-if ((nD & 0x80000000u) == 0u) {
-jI += blasOffset;
+if (eB >= 0 && em > 0.0) {
+let oy = mT[eB];
+var ki = bitcast<i32>(oy);
+if ((oy & 0x80000000u) == 0u) {
+ki += blasOffset;
 }
-dS++;
-gt[dS] = jI;
-fd[el] = -1.0;
+eh++;
+gJ[eh] = ki;
+fq[eB] = -1.0;
 }
 }
 }
 }
 ` : ''}
 ${(blasLayoutMask & 1) !== 0 ? /* wgsl */ `
-if (gR == 0u) {
-let i = node + bo;
+if (he == 0u) {
+let i = node + bw;
 let r1 = bvhBuffer[i];
-let bb = bitcast<i32>(r1.w);
+let bl = bitcast<i32>(r1.w);
 
-if (bb >= 0) {
+if (bl >= 0) {
 let left = node + 2;
-let right = bb + blasOffset;
-let d0 = cM(
+let right = bl + blasOffset;
+let d0 = cU(
 o,
-bvhBuffer[left + bo].xyz,
-bvhBuffer[left + 1 + bo].xyz,
-gC,
+bvhBuffer[left + bw].xyz,
+bvhBuffer[left + 1 + bw].xyz,
+gr,
 tMax
 );
-let d1 = cM(
+let d1 = cU(
 o,
-bvhBuffer[right + bo].xyz,
-bvhBuffer[right + 1 + bo].xyz,
-gC,
+bvhBuffer[right + bw].xyz,
+bvhBuffer[right + 1 + bw].xyz,
+gr,
 tMax
 );
 
 if (d0 > d1) {
 if (d0 > 0.0) {
-dS++; gt[dS] = left;
+eh++; gJ[eh] = left;
 
 }
 if (d1 > 0.0) {
-dS++; gt[dS] = right;
+eh++; gJ[eh] = right;
 
 }
 }
 else {
 if (d1 > 0.0) {
-dS++; gt[dS] = right;
+eh++; gJ[eh] = right;
 
 }
 if (d0 > 0.0) {
-dS++; gt[dS] = left;
+eh++; gJ[eh] = left;
 
 }
 }
 }
 else {
-bh = -bb;
+primitiveCount = -bl;
 primitiveOffset = i + 2;
 }
 }
 ` : ''}
-if (bh > 0) {
+if (primitiveCount > 0) {
 ${intersectTriangles}
 }
 }
@@ -604,19 +604,221 @@ ${intersectTriangles}
 
 `;
 }
-export default function getIntersectCode({ hasDisplacement, hasAnisotropyDirections = true, hasCoatNormals = false, blasLayoutMask = 1 }) {
+// Raw first-hit contract for native diagnostics. Reuse the same traversal and
+// triangle tests without surface completion or any integrator/alpha sampling.
+export function createSceneHitQueryWgsl(hasAnisotropyDirections = false) {
+    return /* wgsl */ `
+struct rd {
+hit: bool, instance: i32, record: i32, t: f32, barycentric: vec3f,
+}
+fn vW(ray: Ray) -> rd {
+var result: rd;
+result.instance = -1; result.record = -1; result.t = RAY_MAX_DISTANCE;
+if (rtUniform.isEmpty > 0) { return result; }
+let bounces = 0;
+let doubleSided = false;
+let jA = false;
+let fV = -1;
+var tMax = RAY_MAX_DISTANCE;
+var eq: i32;
+var gY: vec3f;
+${getBVHTraverseCode(`
+result.hit = true; result.instance = bY; result.record = triangleIndex;
+result.t = tMax; result.barycentric = aC;
+`, false, false, 1, hasAnisotropyDirections)}
+return result;
+}
+`;
+}
+export function createSceneIntersectionsWgsl({ hasDisplacement, hasAnisotropyDirections = true, hasCoatNormals = false, blasLayoutMask = 1, closestName = 'intersectScene', shadowName = 'intersectSceneShadow', bevelName = 'intersectSceneBevel', closestHitHook = '' }) {
+    return /* wgsl */ `
+fn ${closestName}(ray: Ray, si: ptr<function, SurfaceInteraction>, bounces: i32, doubleSided: bool, jA: bool) {
+
+(*si).hit = false;
+
+if (rtUniform.isEmpty > 0) {
+return;
+}
+
+var tMax = RAY_MAX_DISTANCE;
+var eq: i32;
+var gY: vec3f;
+
+var hK: vec3i;
+var mk: fP;
+var ap: i32;
+var fV: i32 = (*si).triangleIndex;
+
+${getBVHTraverseCode(`
+(*si).triangleIndex = triangleIndex;
+mk = tri;
+hK = index;
+ap = bY;
+`, hasDisplacement, false, blasLayoutMask, hasAnisotropyDirections)}
+
+${closestHitHook}  if (sj) {
+mV(
+si,
+mk,
+aC,
+hK,
+bZ,
+ap,
+eL,
+fi,
+-1
+);
+}
+}
+
+
+fn ${bevelName}(ray: Ray, si: ptr<function, SurfaceInteraction>, instanceIndex: i32, maxDistance: f32) {
+(*si).hit = false;
+
+if (rtUniform.isEmpty > 0) {
+return;
+}
+
+let bounces = 0;
+let doubleSided = true;
+let jA = false;
+var tMax = maxDistance;
+var eq: i32;
+var gY: vec3f;
+var hK: vec3i;
+var mk: fP;
+var ap: i32;
+var fV: i32 = -1;
+
+${getBVHTraverseCode(`
+(*si).triangleIndex = triangleIndex;
+mk = tri;
+hK = index;
+ap = bY;
+`, hasDisplacement, false, blasLayoutMask, hasAnisotropyDirections, `
+if (bY != instanceIndex) {
+continue;
+}
+`)}
+
+if (sj) {
+(*si).instanceIndex = ap;
+(*si).normal = eL;
+(*si).position = fi;
+(*si).materialIndex = i32(ah(bitcast<i32>(instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE}].z)).x);
+hH(
+si,
+mk,
+aC,
+hK,
+bZ,
+-1
+);
+
+let l2w = mat4x4<f32>(
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 5],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 6],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 7],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 8]
+);
+let w2l = mat4x4<f32>(
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 1],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 2],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 3],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 4]
+);
+let aH = transpose(mat3x3<f32>(
+w2l[0].xyz,
+w2l[1].xyz,
+w2l[2].xyz
+));
+(*si).faceNormal = normalize(aH * eL);
+(*si).vertexNormal = normalize(aH * (*si).vertexNormal);
+(*si).normal = normalize(aH * (*si).normal);
+${hasCoatNormals ? /* wgsl */ `  (*si).coatNormal = normalize(aH * (*si).coatNormal);` : ''}
+if (USE_TANGENT && dot((*si).tangent, (*si).tangent) > 0.0) {
+(*si).tangent = normalize((l2w * vec4f((*si).tangent, 0.0)).xyz);
+}
+${hasAnisotropyDirections ? /* wgsl */ `
+(*si).directionTangent = (l2w * vec4f((*si).directionTangent, 0.0)).xyz;
+(*si).directionBitangent = (l2w * vec4f((*si).directionBitangent, 0.0)).xyz;
+(*si).geometryTangent = (l2w * vec4f((*si).geometryTangent, 0.0)).xyz;
+(*si).geometryBitangent = (l2w * vec4f((*si).geometryBitangent, 0.0)).xyz;
+` : /* wgsl */ ``}
+(*si).shadowOffset = (l2w * vec4f((*si).shadowOffset, 0.0)).xyz;
+(*si).position = (l2w * vec4((*si).position, 1.0)).xyz;
+}
+}
+
+fn ${shadowName}(ray: Ray, lightDistance: f32, kn: i32) -> bool {
+if (rtUniform.isEmpty > 0) {
+return false;
+}
+
+let doubleSided = true;
+let rq = min(max(lightDistance, 0.0), RAY_MAX_DISTANCE);
+var tMax = rq;
+var eq: i32;
+var gY: vec3f;
+
+let bounces = 1;
+var jA = false;
+var fV: i32 = kn;
+
+${getBVHTraverseCode(
+    /*wgsl */ `
+
+
+
+var jS = true;
+
+if (jS) {
+if (HAS_ALPHA) {
+let materialIndex = i32(ah(bitcast<i32>(pi.z)).x);
+if (!oC(
+materialIndex,
+index,
+aC,
+fi,
+eL
+)) {
+
+tMax = rq;
+jS = false;
+}
+}
+}
+if (jS) {
+return true;
+}
+`, hasDisplacement, true, blasLayoutMask, hasAnisotropyDirections)}
+
+return false;
+}
+
+`;
+}
+export default function getIntersectCode({ hasDisplacement, hasAnisotropyDirections = true, hasCoatNormals = false, blasLayoutMask = 1, sceneIntersections = '' }) {
     return /*wgsl*/ `
 
+
+
+
+
+var<private> mK: array<i32, TLAS_STACK_SIZE>;
+var<private> gJ: array<i32, BLAS_STACK_SIZE>;
+${(blasLayoutMask & 8) !== 0 ? /* wgsl */ `var<private> tD: array<vec2u, BLAS_STACK_SIZE>;` : ''}
+
 ${hasDisplacement ? displacementWgsl : ''}
 ${(blasLayoutMask & 8) !== 0 ? cwbvhTraversalWgsl : ''}
 
-struct fA {
+struct fP {
 p0: vec3f,
 p1: vec3f,
 p2: vec3f,
 };
 
-fn hx(rad: f32, pos: vec3f, r: Ray) -> f32 {
+fn hL(rad: f32, pos: vec3f, r: Ray) -> f32 {
 let op: vec3f = pos - r.o;
 let eps: f32 = 0.001;
 let b: f32 = dot(op, r.d);
@@ -639,7 +841,7 @@ return t2;
 return INF;
 }
 
-fn ni(pos: vec3f, u: vec3f, v: vec3f, n: vec3f, r: Ray) -> f32 {
+fn og(pos: vec3f, u: vec3f, v: vec3f, n: vec3f, r: Ray) -> f32 {
 let dt: f32 = dot(r.d, n);
 let t: f32 = (dot(n, pos) - dot(n, r.o)) / dt;
 
@@ -659,14 +861,14 @@ return INF;
 }
 
 fn sphereSolidAnglePdf(origin: vec3f, center: vec3f, radius: f32) -> f32 {
-let aj = distance(origin, center);
-if (aj <= radius) { return INV4PI; }
-let jv = min(1.0, radius * radius / (aj * aj));
-let cosThetaMax = sqrt(max(0.0, 1.0 - jv));
+let ak = distance(origin, center);
+if (ak <= radius) { return INV4PI; }
+let jV = min(1.0, radius * radius / (ak * ak));
+let cosThetaMax = sqrt(max(0.0, 1.0 - jV));
 return 1.0 / max(TWOPI * (1.0 - cosThetaMax), 1e-8);
 }
 
-fn mz(r: Ray, light: Light, lightSample: ptr<function, LightSampleRec>) -> f32 {
+fn nr(r: Ray, light: Light, lightSample: ptr<function, LightSampleRec>) -> f32 {
 
 
 var d: f32 = INF;
@@ -678,9 +880,9 @@ var d: f32 = INF;
 (*lightSample).isDelta = false;
 var u = light.u;
 var v = light.v;
-var dP = i32(light.litType);
+var ec = i32(light.litType);
 
-if (dP == RECT_LIGHT) {
+if (HAS_SCENE_RECT_LIGHT && ec == RECT_LIGHT) {
 
 var pos = light.position - u * 0.5 - v * 0.5;
 var normal = normalize(cross(u, v));
@@ -688,18 +890,18 @@ if (dot(normal, r.d) > 0.0) {
 return INF;
 }
 
-d = ni(pos, u, v, normal, r);
+d = og(pos, u, v, normal, r);
 if (d < 0.0) {
 d = INF;
 }
 else {
-var cc: f32 = dot(-r.d, normal);
+var ci: f32 = dot(-r.d, normal);
 let spread = light.radiusOrSpread;
-if (cc >= cos(spread * PI / 2.0)) {
+if (ci >= cos(spread * PI / 2.0)) {
 (*lightSample).direction = r.d;
 (*lightSample).normal = normal;
 (*lightSample).dist = d;
-(*lightSample).pdf = (d * d) / (light.area * cc);
+(*lightSample).pdf = (d * d) / (light.area * ci);
 (*lightSample).emission = light.emission;
 (*lightSample).isDelta = false;
 }
@@ -709,12 +911,12 @@ d = INF;
 }
 }
 
-else if (dP == SPHERE_LIGHT) {
+else if (HAS_SCENE_SPHERE_LIGHT && ec == SPHERE_LIGHT) {
 if (light.radiusOrSpread <= 0.0) {
-let ko = light.position - r.o;
-let projection = dot(ko, r.d);
-let ba = ko - projection * r.d;
-if (projection <= EPS || dot(ba, ba) > 1e-8) {
+let kT = light.position - r.o;
+let projection = dot(kT, r.d);
+let bk = kT - projection * r.d;
+if (projection <= EPS || dot(bk, bk) > 1e-8) {
 return INF;
 }
 d = projection;
@@ -725,7 +927,7 @@ d = projection;
 (*lightSample).isDelta = true;
 }
 else {
-d = hx(light.radiusOrSpread, light.position, r);
+d = hL(light.radiusOrSpread, light.position, r);
 (*lightSample).isDelta = false;
 if (d > 0.0 && d != INF) {
 (*lightSample).direction = r.d;
@@ -754,12 +956,12 @@ light.radiusOrSpread <= 0.0
 }
 }
 
-else if (dP == SPOT_LIGHT) {
+else if (HAS_SCENE_SPOT_LIGHT && ec == SPOT_LIGHT) {
 if (light.radiusOrSpread <= 0.0) {
-let ko = light.position - r.o;
-let projection = dot(ko, r.d);
-let ba = ko - projection * r.d;
-if (projection <= EPS || dot(ba, ba) > 1e-8) {
+let kT = light.position - r.o;
+let projection = dot(kT, r.d);
+let bk = kT - projection * r.d;
+if (projection <= EPS || dot(bk, bk) > 1e-8) {
 return INF;
 }
 d = projection;
@@ -770,13 +972,13 @@ d = projection;
 (*lightSample).isDelta = true;
 }
 else {
-d = hx(light.radiusOrSpread, light.position, r);
+d = hL(light.radiusOrSpread, light.position, r);
 (*lightSample).isDelta = false;
 }
 if (d != INF) {
-let ve = r.o + d * r.d;
-let mC = bL(light, normalize(r.o - ve));
-if (mC <= 0.0) { return INF; }
+let wW = r.o + d * r.d;
+let nv = bP(light, normalize(r.o - wW));
+if (nv <= 0.0) { return INF; }
 if (light.radiusOrSpread > 0.0) {
 (*lightSample).direction = r.d;
 (*lightSample).dist = d;
@@ -785,7 +987,7 @@ r.o + r.d * d - light.position
 );
 (*lightSample).pdf = sphereSolidAnglePdf(r.o, light.position, light.radiusOrSpread);
 }
-(*lightSample).emission = light.emission * mC * select(
+(*lightSample).emission = light.emission * nv * select(
 1.0,
 1.0 / max(d * d, 1e-8),
 light.radiusOrSpread <= 0.0
@@ -793,7 +995,7 @@ light.radiusOrSpread <= 0.0
 }
 }
 
-else if (dP != RECT_LIGHT && dP != SPHERE_LIGHT && dP != SPOT_LIGHT) {
+else if (ec != RECT_LIGHT && ec != SPHERE_LIGHT && ec != SPOT_LIGHT) {
 
 d = INF;
 }
@@ -801,111 +1003,111 @@ d = INF;
 return d;
 }
 
-fn fm(m0: vec3i, m1: vec3i, m2: vec3i, eK: vec3f) -> i32 {
-let m00 = ae(m0.x);
-let m01 = ae(m0.y);
-let m02 = ae(m0.z);
-let m10 = ae(m1.x);
-let m11 = ae(m1.y);
-let m12 = ae(m1.z);
-let m20 = ae(m2.x);
-let m21 = ae(m2.y);
-let m22 = ae(m2.z);
-var w00 = m00.y * eK.x;
-var w01 = m01.y * eK.x;
-var w02 = m02.y * eK.x;
-var w10 = m10.y * eK.y;
-var w11 = m11.y * eK.y;
-var w12 = m12.y * eK.y;
-var w20 = m20.y * eK.z;
-var w21 = m21.y * eK.z;
-var w22 = m22.y * eK.z;
-
-var hk = w00 + w01 + w02 + w10 + w11 + w12 + w20 + w21 + w22;
-
-w00 = w00 / hk;
-w01 = w01 / hk + w00;
-w02 = w02 / hk + w01;
-w10 = w10 / hk + w02;
-w11 = w11 / hk + w10;
-w12 = w12 / hk + w11;
-w20 = w20 / hk + w12;
-w21 = w21 / hk + w20;
-w22 = w22 / hk + w21;
-
-let rand = eo(0u);
-var nj: f32;
+fn fF(m0: vec3i, m1: vec3i, m2: vec3i, eU: vec3f) -> i32 {
+let m00 = ah(m0.x);
+let m01 = ah(m0.y);
+let m02 = ah(m0.z);
+let m10 = ah(m1.x);
+let m11 = ah(m1.y);
+let m12 = ah(m1.z);
+let m20 = ah(m2.x);
+let m21 = ah(m2.y);
+let m22 = ah(m2.z);
+var w00 = m00.y * eU.x;
+var w01 = m01.y * eU.x;
+var w02 = m02.y * eU.x;
+var w10 = m10.y * eU.y;
+var w11 = m11.y * eU.y;
+var w12 = m12.y * eU.y;
+var w20 = m20.y * eU.z;
+var w21 = m21.y * eU.z;
+var w22 = m22.y * eU.z;
+
+var hy = w00 + w01 + w02 + w10 + w11 + w12 + w20 + w21 + w22;
+
+w00 = w00 / hy;
+w01 = w01 / hy + w00;
+w02 = w02 / hy + w01;
+w10 = w10 / hy + w02;
+w11 = w11 / hy + w10;
+w12 = w12 / hy + w11;
+w20 = w20 / hy + w12;
+w21 = w21 / hy + w20;
+w22 = w22 / hy + w21;
+
+let rand = eC(0u);
+var oh: f32;
 if (rand < w00) {
-nj = m00.x;
+oh = m00.x;
 } else if (rand < w01) {
-nj = m01.x;
+oh = m01.x;
 } else if (rand < w02) {
-nj = m02.x;
+oh = m02.x;
 } else if (rand < w10) {
-nj = m10.x;
+oh = m10.x;
 } else if (rand < w11) {
-nj = m11.x;
+oh = m11.x;
 } else if (rand < w12) {
-nj = m12.x;
+oh = m12.x;
 } else if (rand < w20) {
-nj = m20.x;
+oh = m20.x;
 } else if (rand < w21) {
-nj = m21.x;
+oh = m21.x;
 } else {
-nj = m22.x;
+oh = m22.x;
 }
-return i32(nj);
+return i32(oh);
 }
 
-fn lW(
-pa: i32,
+fn mF(
+pW: i32,
 index: vec3i,
-eK: vec3f
+eU: vec3f
 ) -> i32 {
 let m0 = vec3i(attrBuffers[index.x * ATTR_STRIDE + MATERIALINDEX_OFF]${hasAnisotropyDirections ? '.xyz' : ''});
 let m1 = vec3i(attrBuffers[index.y * ATTR_STRIDE + MATERIALINDEX_OFF]${hasAnisotropyDirections ? '.xyz' : ''});
 let m2 = vec3i(attrBuffers[index.z * ATTR_STRIDE + MATERIALINDEX_OFF]${hasAnisotropyDirections ? '.xyz' : ''});
 if (m0.x > 0 || m1.x > 0 || m2.x > 0) {
-return fm(m0, m1, m2, eK);
+return fF(m0, m1, m2, eU);
 }
-return pa;
+return pW;
 }
 
 
 
 ${hasAnisotropyDirections ? /* wgsl */ `
-fn eE(normal: vec3f, dp1: vec3f, dp2: vec3f, od: vec2f, oe: vec2f, input: vec4f) -> mat2x3f {
-let determinant = od.x * oe.y - od.y * oe.x;
+fn eP(normal: vec3f, dp1: vec3f, dp2: vec3f, oW: vec2f, oX: vec2f, input: vec4f) -> mat2x3f {
+let determinant = oW.x * oX.y - oW.y * oX.x;
 var t = input.xyz;
 var b = vec3f(0.0);
 
 var sign = select(input.w, 0.0, abs(input.w) > 1.5);
 if (abs(determinant) > 1e-12) {
-let dt = (dp1 * oe.y - dp2 * od.y) / determinant;
-let db = (dp2 * od.x - dp1 * oe.x) / determinant;
+let dt = (dp1 * oX.y - dp2 * oW.y) / determinant;
+let db = (dp2 * oW.x - dp1 * oX.x) / determinant;
 if (dot(t, t) < 1e-8) { t = dt; }
 if (sign == 0.0) { sign = select(1.0, -1.0, dot(cross(normal, dt), db) < 0.0); }
 }
 t -= normal * dot(normal, t);
 if (dot(t, t) < 1e-8) {
-let ik = select(-1.0, 1.0, normal.z >= 0.0);
-let dW = -1.0 / (ik + normal.z);
-let xy = normal.x * normal.y * dW;
-t = vec3f(1.0 + ik * normal.x * normal.x * dW, ik * xy, -ik * normal.x);
-b = vec3f(xy, ik + normal.y * normal.y * dW, -normal.y);
+let iI = select(-1.0, 1.0, normal.z >= 0.0);
+let cL = -1.0 / (iI + normal.z);
+let xy = normal.x * normal.y * cL;
+t = vec3f(1.0 + iI * normal.x * normal.x * cL, iI * xy, -iI * normal.x);
+b = vec3f(xy, iI + normal.y * normal.y * cL, -normal.y);
 }
 else { t = normalize(t); b = normalize(cross(normal, t)); }
 return mat2x3f(t, b * select(1.0, -1.0, sign < 0.0));
 }
 
 ` : /* wgsl */ ``}
-fn ht(
-si: ptr<function, ad>,
-tri: fA,
-eK: vec3f,
+fn hH(
+si: ptr<function, SurfaceInteraction>,
+tri: fP,
+eU: vec3f,
 index: vec3i,
-bU: bool,
-fv: i32
+bZ: bool,
+fM: i32
 ) {
 (*si).hit = true;
 
@@ -917,26 +1119,26 @@ let i1 = index.y;
 let i2 = index.z;
 (*si).vertexColor = vec3f(1.0);
 if (USE_VERTEX_COLOR) {
-(*si).vertexColor = eK.x * attrBuffers[i0 * ATTR_STRIDE + COLOR_OFF].xyz
-+ eK.y * attrBuffers[i1 * ATTR_STRIDE + COLOR_OFF].xyz
-+ eK.z * attrBuffers[i2 * ATTR_STRIDE + COLOR_OFF].xyz;
+(*si).vertexColor = eU.x * attrBuffers[i0 * ATTR_STRIDE + COLOR_OFF].xyz
++ eU.y * attrBuffers[i1 * ATTR_STRIDE + COLOR_OFF].xyz
++ eU.z * attrBuffers[i2 * ATTR_STRIDE + COLOR_OFF].xyz;
 }
-let oY = uvProjection == ${uvProjectionConstants.spherical};
+let pU = uvProjection == ${uvProjectionConstants.spherical};
 
 var normal: vec3f = (*si).normal;
-var gZ: vec3f;
+var hq: vec3f;
 var n0 = vec3f(0.0);
 var n1 = vec3f(0.0);
 var n2 = vec3f(0.0);
 
 
-if (!bU || oY) {
+if (!bZ || pU) {
 n0 = attrBuffers[i0 * ATTR_STRIDE + NORMAL_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
 n1 = attrBuffers[i1 * ATTR_STRIDE + NORMAL_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
 n2 = attrBuffers[i2 * ATTR_STRIDE + NORMAL_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
-gZ = normalize(eK.x * n0 + eK.y * n1 + eK.z * n2);
-if (!bU) {
-normal = gZ;
+hq = normalize(eU.x * n0 + eU.y * n1 + eU.z * n2);
+if (!bZ) {
+normal = hq;
 }
 }
 
@@ -944,19 +1146,19 @@ var uv: vec2f;
 var uv0: vec2f;
 var uv1: vec2f;
 var uv2: vec2f;
-var cX = vec2f(0.0);
+var dj = vec2f(0.0);
 if (USE_UV) {
 uv0 = attrBuffers[i0 * ATTR_STRIDE + UV_OFF].xy;
 uv1 = attrBuffers[i1 * ATTR_STRIDE + UV_OFF].xy;
 uv2 = attrBuffers[i2 * ATTR_STRIDE + UV_OFF].xy;
-uv = eK.x * uv0 + eK.y * uv1 + eK.z * uv2;
+uv = eU.x * uv0 + eU.y * uv1 + eU.z * uv2;
 if (USE_UV1) {
 let secondaryUv0 = attrBuffers[i0 * ATTR_STRIDE + UV1_OFF].xy;
 let secondaryUv1 = attrBuffers[i1 * ATTR_STRIDE + UV1_OFF].xy;
 let secondaryUv2 = attrBuffers[i2 * ATTR_STRIDE + UV1_OFF].xy;
-cX = eK.x * secondaryUv0
-+ eK.y * secondaryUv1
-+ eK.z * secondaryUv2;
+dj = eU.x * secondaryUv0
++ eU.y * secondaryUv1
++ eU.z * secondaryUv2;
 }
 }
 
@@ -964,28 +1166,28 @@ if (USE_TANGENT) {
 let t0 = attrBuffers[i0 * ATTR_STRIDE + TANGENT_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
 let t1 = attrBuffers[i1 * ATTR_STRIDE + TANGENT_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
 let t2 = attrBuffers[i2 * ATTR_STRIDE + TANGENT_OFF]${hasAnisotropyDirections ? '.xyz' : ''};
-(*si).tangent = normalize(eK.x * t0 + eK.y * t1 + eK.z * t2);
+(*si).tangent = normalize(eU.x * t0 + eU.y * t1 + eU.z * t2);
 }
 
-if (fv >= 0) {
-(*si).materialIndex = fv;
+if (fM >= 0) {
+(*si).materialIndex = fM;
 }
 else {
-(*si).materialIndex = lW(
+(*si).materialIndex = mF(
 (*si).materialIndex,
 index,
-eK
+eU
 );
 }
 
 if (uvProjection > 0) {
 
-uv = gH(uvProjection, (*si).position, select(normal, gZ, oY));
-cX = uv;
+uv = gU(uvProjection, (*si).position, select(normal, hq, pU));
+dj = uv;
 }
 
 uv = uv * material.uvScale;
-cX = cX * material.uvScale;
+dj = dj * material.uvScale;
 
 
 if (uvProjection == ${uvProjectionConstants.triplanar}) {
@@ -995,44 +1197,44 @@ uv = fract(uv);
 ${hasAnisotropyDirections ? /* wgsl */ `
 ${hasCoatNormals ? /* wgsl */ `  (*si).coatNormal = normal;` : ''}
 let directionMaterial = materialsBuffer[(*si).materialIndex];
-if (directionMaterial.anisotropyMapIndex >= 0.0 || directionMaterial.specularAnisotropyRotation != 0.0 || directionMaterial.coatAnisotropyRotation != 0.0${hasCoatNormals ? " || (USE_COAT_NORMAL_MAPS && directionMaterial.coatNormalMapIndex >= 0.0)" : ""}) {
+if (directionMaterial.anisotropyMapIndex >= 0.0 || directionMaterial.specularAnisotropyRotation != 0.0 || directionMaterial.coatAnisotropyRotation != 0.0${hasCoatNormals ? " || directionMaterial.coatWeight > 0.0" : ""}) {
 let dp1 = tri.p0 - tri.p2;
 let dp2 = tri.p1 - tri.p2;
-var aT = vec4f(0.0);
+var aX = vec4f(0.0);
 if (USE_TANGENT) {
-aT = eK.x * attrBuffers[i0 * ATTR_STRIDE + TANGENT_OFF]
-+ eK.y * attrBuffers[i1 * ATTR_STRIDE + TANGENT_OFF]
-+ eK.z * attrBuffers[i2 * ATTR_STRIDE + TANGENT_OFF];
-}
-let eS = eE(normal, dp1, dp2, uv0-uv2, uv1-uv2, aT);
-(*si).geometryTangent = eS[0];
-(*si).geometryBitangent = eS[1];
-var fY = eS;
+aX = eU.x * attrBuffers[i0 * ATTR_STRIDE + TANGENT_OFF]
++ eU.y * attrBuffers[i1 * ATTR_STRIDE + TANGENT_OFF]
++ eU.z * attrBuffers[i2 * ATTR_STRIDE + TANGENT_OFF];
+}
+let fc = eP(normal, dp1, dp2, uv0-uv2, uv1-uv2, aX);
+(*si).geometryTangent = fc[0];
+(*si).geometryBitangent = fc[1];
+var gm = fc;
 if (USE_UV1 && directionMaterial.anisotropyMapIndex >= 0.0 && mapInfo[i32(directionMaterial.anisotropyMapIndex) + NUM_MAPS * 2].x > 0.5) {
 let a = attrBuffers[i0 * ATTR_STRIDE + UV1_OFF].xy;
 let b = attrBuffers[i1 * ATTR_STRIDE + UV1_OFF].xy;
 let c = attrBuffers[i2 * ATTR_STRIDE + UV1_OFF].xy;
-fY = eE(normal, dp1, dp2, a-c, b-c, select(aT, vec4f(0.0), aT.w == 0.0));
+gm = eP(normal, dp1, dp2, a-c, b-c, select(aX, vec4f(0.0), aX.w == 0.0));
 }
-(*si).directionTangent = fY[0];
-(*si).directionBitangent = fY[1];
+(*si).directionTangent = gm[0];
+(*si).directionBitangent = gm[1];
 ${hasCoatNormals ? /* wgsl */ `
 if (USE_COAT_NORMAL_MAPS && directionMaterial.coatNormalMapIndex >= 0.0) {
-let kC = i32(directionMaterial.coatNormalMapIndex);
-var mi = eS;
-if (USE_UV1 && mapInfo[kC + NUM_MAPS * 2].x > 0.5) {
+let lg = i32(directionMaterial.coatNormalMapIndex);
+var mU = fc;
+if (USE_UV1 && mapInfo[lg + NUM_MAPS * 2].x > 0.5) {
 let a = attrBuffers[i0 * ATTR_STRIDE + UV1_OFF].xy;
 let b = attrBuffers[i1 * ATTR_STRIDE + UV1_OFF].xy;
 let c = attrBuffers[i2 * ATTR_STRIDE + UV1_OFF].xy;
-mi = eE(normal, dp1, dp2, a-c, b-c,
-select(aT, vec4f(0.0), aT.w == 0.0));
+mU = eP(normal, dp1, dp2, a-c, b-c,
+select(aX, vec4f(0.0), aX.w == 0.0));
 }
 
 
-let decoded = sampleMap(kC, getMapUvSet(kC, uv,
-select(cX, uv, uvProjection > 0))).rgb * 2.0 - 1.0;
-let mapped = mi[0] * decoded.x * directionMaterial.coatNormalScale.x
-+ mi[1] * decoded.y * directionMaterial.coatNormalScale.y + normal * decoded.z;
+let decoded = sampleMap(lg, getMapUvSet(lg, uv,
+select(dj, uv, uvProjection > 0))).rgb * 2.0 - 1.0;
+let mapped = mU[0] * decoded.x * directionMaterial.coatNormalScale.x
++ mU[1] * decoded.y * directionMaterial.coatNormalScale.y + normal * decoded.z;
 if (dot(mapped, mapped) > 1e-12) { (*si).coatNormal = normalize(mapped); }
 }
 ` : ''}
@@ -1042,19 +1244,19 @@ if (dot(mapped, mapped) > 1e-12) { (*si).coatNormal = normalize(mapped); }
 (*si).normal = normal;
 (*si).vertexNormal = normal;
 (*si).shadowOffset = vec3f(0.0);
-if (!bU) {
+if (!bZ) {
 
 
 
 let geometricNormal = normalize(cross(tri.p1 - tri.p0, tri.p2 - tri.p0));
-let km = eK.x * n0 + eK.y * n1 + eK.z * n2;
+let kR = eU.x * n0 + eU.y * n1 + eU.z * n2;
 let a = dot(n2 - n0, tri.p0 - tri.p2);
 let b = dot(n2 - n1, tri.p1 - tri.p2);
 let c = dot(n1 - n0, tri.p1 - tri.p0);
-var h = a * eK.x * (eK.x - 1.0)
-+ (a + b + c) * eK.x * eK.y
-+ b * eK.y * (eK.y - 1.0);
-if (dot(km, geometricNormal) > 0.0) {
+var h = a * eU.x * (eU.x - 1.0)
++ (a + b + c) * eU.x * eU.y
++ b * eU.y * (eU.y - 1.0);
+if (dot(kR, geometricNormal) > 0.0) {
 var h0 = max(max(dot(tri.p1 - tri.p0, n0), dot(tri.p2 - tri.p0, n0)), 0.0);
 var h1 = max(max(dot(tri.p0 - tri.p1, n1), dot(tri.p2 - tri.p1, n1)), 0.0);
 var h2 = max(max(dot(tri.p0 - tri.p2, n2), dot(tri.p1 - tri.p2, n2)), 0.0);
@@ -1072,22 +1274,22 @@ h1 = max(dot((*si).position - tri.p1, n1) + h1, 0.0);
 h2 = max(dot((*si).position - tri.p2, n2) + h2, 0.0);
 h = min(-min(min(h0, h1), h2), h * 0.5);
 }
-(*si).shadowOffset = km * h;
+(*si).shadowOffset = kR * h;
 }
 (*si).uv = uv;
-(*si).uv1 = cX;
+(*si).uv1 = dj;
 if (USE_NORMAL_MAPS) {
 let dp1 = tri.p0 - tri.p2;
 let dp2 = tri.p1 - tri.p2;
-let od = uv0 - uv2;
-let oe = uv1 - uv2;
-(*si).normal = om(
-(*si).materialIndex, uv, normal, dp1, dp2, od, oe
+let oW = uv0 - uv2;
+let oX = uv1 - uv2;
+(*si).normal = pd(
+(*si).materialIndex, uv, normal, dp1, dp2, oW, oX
 );
 }
 }
 
-struct eC {
+struct eK {
 t: f32,
 barycentric: vec3f,
 };
@@ -1096,8 +1298,8 @@ barycentric: vec3f,
 
 
 
-fn fL(o: vec3f, d: vec3f, tri: fA, ey: i32, gP: vec3f, tMax: f32) -> eC {
-var ti: eC;
+fn ga(o: vec3f, d: vec3f, tri: fP, eq: i32, gY: vec3f, tMax: f32) -> eK {
+var ti: eK;
 
 
 var p0t = tri.p0 - o;
@@ -1105,11 +1307,11 @@ var p1t = tri.p1 - o;
 var p2t = tri.p2 - o;
 
 
-if (ey == 0) {
+if (eq == 0) {
 p0t = p0t.yzx;
 p1t = p1t.yzx;
 p2t = p2t.yzx;
-} else if (ey == 1) {
+} else if (eq == 1) {
 p0t = p0t.zxy;
 p1t = p1t.zxy;
 p2t = p2t.zxy;
@@ -1117,9 +1319,9 @@ p2t = p2t.zxy;
 
 
 
-p0t = vec3f(p0t.xy + gP.xy * p0t.z, p0t.z);
-p1t = vec3f(p1t.xy + gP.xy * p1t.z, p1t.z);
-p2t = vec3f(p2t.xy + gP.xy * p2t.z, p2t.z);
+p0t = vec3f(p0t.xy + gY.xy * p0t.z, p0t.z);
+p1t = vec3f(p1t.xy + gY.xy * p1t.z, p1t.z);
+p2t = vec3f(p2t.xy + gY.xy * p2t.z, p2t.z);
 
 
 var e = vec3f(
@@ -1131,9 +1333,9 @@ p0t.x * p1t.y - p0t.y * p1t.x
 if (DEPTH_QUERY_MODE) {
 
 
-let rE = abs(vec3f(p1t.x * p2t.y, p2t.x * p0t.y, p0t.x * p1t.y))
+let sF = abs(vec3f(p1t.x * p2t.y, p2t.x * p0t.y, p0t.x * p1t.y))
 + abs(vec3f(p1t.y * p2t.x, p2t.y * p0t.x, p0t.y * p1t.x));
-e = select(e, vec3f(0.0), abs(e) <= rE * 1e-6);
+e = select(e, vec3f(0.0), abs(e) <= sF * 1e-6);
 }
 
 
@@ -1144,19 +1346,19 @@ return ti;
 let det = e.x + e.y + e.z;
 if (DEPTH_QUERY_MODE && det == 0.0) { return ti; }
 
-p0t.z *= gP.z;
-p1t.z *= gP.z;
-p2t.z *= gP.z;
-let rn = (e.x * p0t.z + e.y * p1t.z + e.z * p2t.z);
+p0t.z *= gY.z;
+p1t.z *= gY.z;
+p2t.z *= gY.z;
+let sw = (e.x * p0t.z + e.y * p1t.z + e.z * p2t.z);
 
 
-if (abs(rn) > abs(tMax * det)) {
+if (abs(sw) > abs(tMax * det)) {
 return ti;
 }
 
-let sw = 1. / det;
-ti.t = rn * sw;
-ti.barycentric = e * sw;
+let tN = 1. / det;
+ti.t = sw * tN;
+ti.barycentric = e * tN;
 
 return ti;
 }
@@ -1166,26 +1368,26 @@ return ti;
 
 override DEPTH_QUERY_MODE: bool = false;
 
-fn fb(d: vec3f) -> vec3f {
+fn fo(d: vec3f) -> vec3f {
 if (!DEPTH_QUERY_MODE) { return 1.0 / d; }
 return 1.0 / select(d, select(vec3f(1e-20), vec3f(-1e-20), d < vec3f(0.0)), abs(d) < vec3f(1e-20));
 }
 
-fn eZ(gC: vec3f, o: vec3f) -> vec3f {
+fn fk(gr: vec3f, o: vec3f) -> vec3f {
 if (!DEPTH_QUERY_MODE) { return vec3f(0.0); }
-return abs(gC) * max(vec3f(1.0), abs(o)) * 1e-6;
+return abs(gr) * max(vec3f(1.0), abs(o)) * 1e-6;
 }
 
 
 
-fn cM(o: vec3f, vK: vec3f, vJ: vec3f, gC: vec3f, tMax: f32) -> f32 {
-let uR = (vK - o) * gC;
-let uS = (vJ - o) * gC;
-let bJ = eZ(gC, o);
-let mL = min(uR, uS) - bJ;
-let rU = max(uR, uS) + bJ;
-let t0 = max(mL.x, max(mL.y, mL.z));
-let t1 = min(rU.x, min(rU.y, rU.z));
+fn cU(o: vec3f, xq: vec3f, xp: vec3f, gr: vec3f, tMax: f32) -> f32 {
+let wG = (xq - o) * gr;
+let wI = (xp - o) * gr;
+let padding = fk(gr, o);
+let nJ = min(wG, wI) - padding;
+let sY = max(wG, wI) + padding;
+let t0 = max(nJ.x, max(nJ.y, nJ.z));
+let t1 = min(sY.x, min(sY.y, sY.z));
 return select(select(t1, t0, t0 > 0.0), -1.0, (t0 > t1 || t0 > tMax));
 }
 
@@ -1193,69 +1395,69 @@ fn maxDimension(v: vec3f) -> i32 {
 return select(select(2, 1, v.y > v.z), select(2, 0, v.x > v.z), v.x > v.y);
 }
 
-fn mv(d: vec3f, ey: i32) -> vec3f {
+fn nk(d: vec3f, eq: i32) -> vec3f {
 
 
-if (ey == 0) {
+if (eq == 0) {
 return vec3f(-d.y, -d.z, 1.0) / d.x;
-} else if (ey == 1) {
+} else if (eq == 1) {
 return vec3f(-d.z, -d.x, 1.0) / d.y;
 } else {
 return vec3f(-d.x, -d.y, 1.0) / d.z;
 }
 }
 
-fn mj(
-si: ptr<function, ad>,
-lC: fA,
-aB: vec3f,
-hw: vec3i,
-bU: bool,
-ao: i32,
-eD: vec3f,
-eY: vec3f,
-fv: i32
+fn mV(
+si: ptr<function, SurfaceInteraction>,
+mk: fP,
+aC: vec3f,
+hK: vec3i,
+bZ: bool,
+ap: i32,
+eL: vec3f,
+fi: vec3f,
+fM: i32
 ) {
-(*si).instanceIndex = ao;
-(*si).barycentrics = aB.xy;
-(*si).normal = eD;
-(*si).position = eY;
-(*si).materialIndex = i32(ae(
-bitcast<i32>(instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE}].z)
+(*si).instanceIndex = ap;
+(*si).barycentrics = aC.xy;
+(*si).normal = eL;
+(*si).position = fi;
+(*si).materialIndex = i32(ah(
+bitcast<i32>(instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE}].z)
 ).x);
 
 
-ht(
+hH(
 si,
-lC,
-aB,
-hw,
-bU,
-fv
+mk,
+aC,
+hK,
+bZ,
+fM
 );
 
 let l2w = mat4x4<f32>(
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 5],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 6],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 7],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 8]
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 5],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 6],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 7],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 8]
 );
 let w2l = mat4x4<f32>(
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 1],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 2],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 3],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 4]
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 1],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 2],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 3],
+instancesBuffer[ap * ${INSTANCE_VEC4_STRIDE} + 4]
 );
-let aD = transpose(mat3x3<f32>(
+let aH = transpose(mat3x3<f32>(
 w2l[0].xyz,
 w2l[1].xyz,
 w2l[2].xyz
 ));
 
-(*si).faceNormal = normalize(aD * eD);
-(*si).vertexNormal = normalize(aD * (*si).vertexNormal);
-(*si).normal = normalize(aD * (*si).normal);
-${hasCoatNormals ? /* wgsl */ `  (*si).coatNormal = normalize(aD * (*si).coatNormal);` : ''}
+(*si).faceNormal = normalize(aH * eL);
+(*si).vertexNormal = normalize(aH * (*si).vertexNormal);
+(*si).normal = normalize(aH * (*si).normal);
+${hasCoatNormals ? /* wgsl */ `  (*si).coatNormal = normalize(aH * (*si).coatNormal);` : ''}
 if (USE_TANGENT && dot((*si).tangent, (*si).tangent) > 0.0) {
 (*si).tangent = normalize((l2w * vec4f((*si).tangent, 0.0)).xyz);
 }
@@ -1274,171 +1476,6 @@ ${createIntersectionAlphaWgsl(hasAnisotropyDirections)}
 
 
 
-
-fn cZ(ray: Ray, si: ptr<function, ad>, bounces: i32, doubleSided: bool, lD: bool) {
-
-(*si).hit = false;
-
-if (rtUniform.isEmpty > 0) {
-return;
-}
-
-var tMax = RAY_MAX_DISTANCE;
-var ey: i32;
-var gP: vec3f;
-
-var hw: vec3i;
-var lC: fA;
-var ao: i32;
-var hO: i32 = (*si).triangleIndex;
-
-${getBVHTraverseCode(`
-(*si).triangleIndex = triangleIndex;
-lC = tri;
-hw = index;
-ao = cf;
-`, hasDisplacement, false, blasLayoutMask, hasAnisotropyDirections)}
-
-if (rb) {
-mj(
-si,
-lC,
-aB,
-hw,
-bU,
-ao,
-eD,
-eY,
--1
-);
-}
-}
-
-
-fn iM(ray: Ray, si: ptr<function, ad>, instanceIndex: i32, maxDistance: f32) {
-(*si).hit = false;
-
-if (rtUniform.isEmpty > 0) {
-return;
-}
-
-let bounces = 0;
-let doubleSided = true;
-let lD = false;
-var tMax = maxDistance;
-var ey: i32;
-var gP: vec3f;
-var hw: vec3i;
-var lC: fA;
-var ao: i32;
-var hO: i32 = -1;
-
-${getBVHTraverseCode(`
-(*si).triangleIndex = triangleIndex;
-lC = tri;
-hw = index;
-ao = cf;
-`, hasDisplacement, false, blasLayoutMask, hasAnisotropyDirections, `
-if (cf != instanceIndex) {
-continue;
-}
-`)}
-
-if (rb) {
-(*si).instanceIndex = ao;
-(*si).normal = eD;
-(*si).position = eY;
-(*si).materialIndex = i32(ae(bitcast<i32>(instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE}].z)).x);
-ht(
-si,
-lC,
-aB,
-hw,
-bU,
--1
-);
-
-let l2w = mat4x4<f32>(
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 5],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 6],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 7],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 8]
-);
-let w2l = mat4x4<f32>(
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 1],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 2],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 3],
-instancesBuffer[ao * ${INSTANCE_VEC4_STRIDE} + 4]
-);
-let aD = transpose(mat3x3<f32>(
-w2l[0].xyz,
-w2l[1].xyz,
-w2l[2].xyz
-));
-(*si).faceNormal = normalize(aD * eD);
-(*si).vertexNormal = normalize(aD * (*si).vertexNormal);
-(*si).normal = normalize(aD * (*si).normal);
-${hasCoatNormals ? /* wgsl */ `  (*si).coatNormal = normalize(aD * (*si).coatNormal);` : ''}
-if (USE_TANGENT && dot((*si).tangent, (*si).tangent) > 0.0) {
-(*si).tangent = normalize((l2w * vec4f((*si).tangent, 0.0)).xyz);
-}
-${hasAnisotropyDirections ? /* wgsl */ `
-(*si).directionTangent = (l2w * vec4f((*si).directionTangent, 0.0)).xyz;
-(*si).directionBitangent = (l2w * vec4f((*si).directionBitangent, 0.0)).xyz;
-(*si).geometryTangent = (l2w * vec4f((*si).geometryTangent, 0.0)).xyz;
-(*si).geometryBitangent = (l2w * vec4f((*si).geometryBitangent, 0.0)).xyz;
-` : /* wgsl */ ``}
-(*si).shadowOffset = (l2w * vec4f((*si).shadowOffset, 0.0)).xyz;
-(*si).position = (l2w * vec4((*si).position, 1.0)).xyz;
-}
-}
-
-fn ex(ray: Ray, lightDistance: f32, jL: i32) -> bool {
-if (rtUniform.isEmpty > 0) {
-return false;
-}
-
-let doubleSided = true;
-let qo = min(max(lightDistance, 0.0), RAY_MAX_DISTANCE);
-var tMax = qo;
-var ey: i32;
-var gP: vec3f;
-
-let bounces = 1;
-var lD = false;
-var hO: i32 = jL;
-
-${getBVHTraverseCode(
-    /*wgsl */ `
-
-
-
-var jt = true;
-
-if (jt) {
-if (HAS_ALPHA) {
-let materialIndex = i32(ae(bitcast<i32>(ou.z)).x);
-if (!nG(
-materialIndex,
-index,
-aB,
-eY,
-eD
-)) {
-
-tMax = qo;
-jt = false;
-}
-}
-}
-if (jt) {
-return true;
-}
-`, hasDisplacement, true, blasLayoutMask, hasAnisotropyDirections)}
-
-return false;
-}
-
-`;
+${sceneIntersections || createSceneIntersectionsWgsl({ hasDisplacement, hasAnisotropyDirections, hasCoatNormals, blasLayoutMask })}`;
 }
 //# sourceMappingURL=intersect.wgsl.js.map