@meimu/mutracer 0.3.0 → 0.3.1

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

@@ -12,25 +12,25 @@ int triangleIndex = primitiveOffset + k;
 ${skipTriangleCode}
 vec4 r3 = fetchBVH(triangleIndex);
 ivec3 index = floatBitsToInt(r3.xyz) + vertexOffset;
-fw tri = fw(
+fL tri = fL(
 fetchData(positionBuffer, index.x, vertexColumns).xyz,
 fetchData(positionBuffer, index.y, vertexColumns).xyz,
 fetchData(positionBuffer, index.z, vertexColumns).xyz
 );
-ev hit = fM(o, d, tri, dZ, gH, tMax);
+eF hit = fY(o, d, tri, ei, gZ, tMax);
 
 if (hit.t > 0.0) {
 tMax = hit.t;
 vec3 e1 = tri.p1 - tri.p0;
 vec3 e2 = tri.p2 - tri.p0;
 vec3 faceNormal = normalize(cross(e1, e2));
-int ku = int(af(
+int lc = int(af(
 floatBitsToInt(fetchIns(bK * ${INSTANCE_VEC4_STRIDE}).z)
 ).x);
-bool iL = fetchMatData(
-ku * MAT_DATA_SIZE + MAT_ROW_ALPHA
+bool jm = fetchMatData(
+lc * MAT_DATA_SIZE + MAT_ROW_ALPHA
 ).z > 0.5;
-if (!doubleSided && !iL && dot(faceNormal, d) > 0.0) {
+if (!doubleSided && !jm && dot(faceNormal, d) > 0.0) {
 continue;
 }
 
@@ -41,31 +41,31 @@ ${handleHitCode}
     return glsl `
 ivec2 stack[STACK_SIZE];
 int top = -1;
-int ug = -1;
+int vQ = -1;
 float tMax = ${initialTMax};
 
 vec3 o = ray.o;
 vec3 d = ray.d;
 
-bool sd = false;
-vec3 hz;
+bool tD = false;
+vec3 hU;
 vec3 aw;
-ivec3 hr;
+ivec3 hL;
 
-fw lR;
+fL mL;
 
 
 top = 0;
 stack[0] = ivec2(0, -1);
 
-int yK = 0;
-while(top >= 0 && yK < 5000) {
-yK++;
+int AI = 0;
+while(top >= 0 && AI < 5000) {
+AI++;
 ivec2 node = stack[top--];
 
 int bK = node.y;
 
-if (bK != ug) {
+if (bK != vQ) {
 if (bK < 0) {
 o = ray.o;
 d = ray.d;
@@ -80,111 +80,111 @@ fetchIns(bK * ${INSTANCE_VEC4_STRIDE} + 4)
 o = (w2l * vec4(ray.o, 1.0)).xyz;
 d = (w2l * vec4(ray.d, 0.0)).xyz;
 }
-dZ = maxDimension(abs(d));
-gH = mU(d, dZ);
-ug = bK;
+ei = maxDimension(abs(d));
+gZ = nV(d, ei);
+vQ = bK;
 }
 
 int bufferOff = bK < 0 ? 0 : blasOffset;
-uint fK = bK < 0 ? 0u : floatBitsToUint(fetchIns(bK * ${INSTANCE_VEC4_STRIDE}).x);
-int gN = int(fK >> 30u);
-int eQ = int(fK & 0x3fffffffu);
+uint fW = bK < 0 ? 0u : floatBitsToUint(fetchIns(bK * ${INSTANCE_VEC4_STRIDE}).x);
+int hj = int(fW >> 30u);
+int fc = int(fW & 0x3fffffffu);
 int primitiveCount = 0;
 int primitiveOffset = 0;
 
-if (gN == 3) {
-uint hT = uint(node.x);
-if ((hT & 0x80000000u) != 0u) {
-primitiveCount = int((hT >> 23u) & 0xffu);
-primitiveOffset = bufferOff + eQ + int(hT & 0x007fffffu);
+if (hj == 3) {
+uint iE = uint(node.x);
+if ((iE & 0x80000000u) != 0u) {
+primitiveCount = int((iE >> 23u) & 0xffu);
+primitiveOffset = bufferOff + fc + int(iE & 0x007fffffu);
 }
 else {
 int i = node.x + bufferOff;
-vec4 jF = fetchBVH(i);
-vec4 fF = fetchBVH(i + 1);
-vec4 nI = fetchBVH(i + 2);
-vec4 nJ = fetchBVH(i + 3);
-vec4 nK = fetchBVH(i + 4);
-uint header = floatBitsToUint(jF.w);
-int sD = int(header << 24u) >> 24;
-int sE = int((header >> 8u) << 24u) >> 24;
-int sF = int((header >> 16u) << 24u) >> 24;
-vec3 hG = vec3(
-uintBitsToFloat(uint(sD + 127) << 23u),
-uintBitsToFloat(uint(sE + 127) << 23u),
-uintBitsToFloat(uint(sF + 127) << 23u)
+vec4 kj = fetchBVH(i);
+vec4 fR = fetchBVH(i + 1);
+vec4 oG = fetchBVH(i + 2);
+vec4 oH = fetchBVH(i + 3);
+vec4 oI = fetchBVH(i + 4);
+uint header = floatBitsToUint(kj.w);
+int uc = int(header << 24u) >> 24;
+int ud = int((header >> 8u) << 24u) >> 24;
+int ue = int((header >> 16u) << 24u) >> 24;
+vec3 ih = vec3(
+uintBitsToFloat(uint(uc + 127) << 23u),
+uintBitsToFloat(uint(ud + 127) << 23u),
+uintBitsToFloat(uint(ue + 127) << 23u)
 );
-uint childBase = floatBitsToUint(fF.x);
-uint ps = floatBitsToUint(fF.y);
+uint childBase = floatBitsToUint(fR.x);
+uint qE = floatBitsToUint(fR.y);
 uint internalMask = header >> 24u;
-uint wm = (d.x < 0.0 ? 4u : 0u) |
+uint yl = (d.x < 0.0 ? 4u : 0u) |
 (d.y < 0.0 ? 2u : 0u) |
 (d.z < 0.0 ? 1u : 0u);
-uint gF[12] = uint[12](
-floatBitsToUint(nI.x), floatBitsToUint(nI.y),
-floatBitsToUint(nI.z), floatBitsToUint(nI.w),
-floatBitsToUint(nJ.x), floatBitsToUint(nJ.y),
-floatBitsToUint(nJ.z), floatBitsToUint(nJ.w),
-floatBitsToUint(nK.x), floatBitsToUint(nK.y),
-floatBitsToUint(nK.z), floatBitsToUint(nK.w)
+uint gX[12] = uint[12](
+floatBitsToUint(oG.x), floatBitsToUint(oG.y),
+floatBitsToUint(oG.z), floatBitsToUint(oG.w),
+floatBitsToUint(oH.x), floatBitsToUint(oH.y),
+floatBitsToUint(oH.z), floatBitsToUint(oH.w),
+floatBitsToUint(oI.x), floatBitsToUint(oI.y),
+floatBitsToUint(oI.z), floatBitsToUint(oI.w)
 );
 
 for (int order = 7; order >= 0; order--) {
-uint slot = uint(order) ^ wm;
-uint tW = slot >= 4u
-? floatBitsToUint(fF.w)
-: floatBitsToUint(fF.z);
-uint meta = (tW >> ((slot & 3u) * 8u)) & 0xffu;
+uint slot = uint(order) ^ yl;
+uint vF = slot >= 4u
+? floatBitsToUint(fR.w)
+: floatBitsToUint(fR.z);
+uint meta = (vF >> ((slot & 3u) * 8u)) & 0xffu;
 if (meta == 0u) continue;
-uint mt = (slot & 3u) * 8u;
-uint lh = slot >> 2u;
-vec3 ve = vec3(
-float((gF[lh] >> mt) & 0xffu),
-float((gF[2u + lh] >> mt) & 0xffu),
-float((gF[4u + lh] >> mt) & 0xffu)
+uint nr = (slot & 3u) * 8u;
+uint me = slot >> 2u;
+vec3 wS = vec3(
+float((gX[me] >> nr) & 0xffu),
+float((gX[2u + me] >> nr) & 0xffu),
+float((gX[4u + me] >> nr) & 0xffu)
 );
-vec3 vd = vec3(
-float((gF[6u + lh] >> mt) & 0xffu),
-float((gF[8u + lh] >> mt) & 0xffu),
-float((gF[10u + lh] >> mt) & 0xffu)
+vec3 wR = vec3(
+float((gX[6u + me] >> nr) & 0xffu),
+float((gX[8u + me] >> nr) & 0xffu),
+float((gX[10u + me] >> nr) & 0xffu)
 );
-vec3 xn = jF.xyz + ve * hG;
-vec3 xm = jF.xyz + vd * hG;
-if (cH(o, d, xn, xm, tMax) <= 0.0) continue;
-uint jH = meta & 0x1fu;
-if (jH >= 24u) {
-uint oa = internalMask & ((1u << slot) - 1u);
-uint ur = childBase + rT(oa);
+vec3 zp = kj.xyz + wS * ih;
+vec3 zo = kj.xyz + wR * ih;
+if (bP(o, d, zp, zo, tMax) <= 0.0) continue;
+uint kl = meta & 0x1fu;
+if (kl >= 24u) {
+uint oZ = internalMask & ((1u << slot) - 1u);
+uint wb = childBase + tl(oZ);
 stack[++top] = ivec2(
-eQ + int(ur * 5u), bK
+fc + int(wb * 5u), bK
 );
 }
 else {
-uint triangleCount = rT(meta >> 5u);
-uint sT = ps + jH;
-uint xC = 0x80000000u |
+uint triangleCount = tl(meta >> 5u);
+uint ut = qE + kl;
+uint zC = 0x80000000u |
 (triangleCount << 23u) |
-(sT & 0x007fffffu);
-stack[++top] = ivec2(int(xC), bK);
+(ut & 0x007fffffu);
+stack[++top] = ivec2(int(zC), bK);
 }
 }
 }
 }
-else if (gN == 2) {
+else if (hj == 2) {
 int i = node.x + bufferOff;
-vec4 xq = fetchBVH(i);
-vec4 xr = fetchBVH(i + 1);
-vec4 xs = fetchBVH(i + 2);
-vec4 xt = fetchBVH(i + 3);
-primitiveCount = int(floatBitsToUint(xs.w));
+vec4 zr = fetchBVH(i);
+vec4 zs = fetchBVH(i + 1);
+vec4 zt = fetchBVH(i + 2);
+vec4 zu = fetchBVH(i + 3);
+primitiveCount = int(floatBitsToUint(zt.w));
 if (primitiveCount > 0) {
-primitiveOffset = bufferOff + eQ + int(floatBitsToUint(xt.w));
+primitiveOffset = bufferOff + fc + int(floatBitsToUint(zu.w));
 }
 else {
-int left = eQ + int(floatBitsToUint(xq.w) * 4u);
-int right = eQ + int(floatBitsToUint(xr.w) * 4u);
-float d0 = cH(o, d, xq.xyz, xr.xyz, tMax);
-float d1 = cH(o, d, xs.xyz, xt.xyz, tMax);
+int left = fc + int(floatBitsToUint(zr.w) * 4u);
+int right = fc + int(floatBitsToUint(zs.w) * 4u);
+float d0 = bP(o, d, zr.xyz, zs.xyz, tMax);
+float d1 = bP(o, d, zt.xyz, zu.xyz, tMax);
 if (d0 > d1) {
 if (d0 > 0.0) stack[++top] = ivec2(left, bK);
 if (d1 > 0.0) stack[++top] = ivec2(right, bK);
@@ -195,62 +195,62 @@ if (d0 > 0.0) stack[++top] = ivec2(left, bK);
 }
 }
 }
-else if (gN == 1) {
-uint hT = uint(node.x);
-if ((hT & 0x80000000u) != 0u) {
-primitiveCount = int((hT >> 23u) & 0xffu);
-primitiveOffset = bufferOff + eQ + int(hT & 0x007fffffu);
+else if (hj == 1) {
+uint iE = uint(node.x);
+if ((iE & 0x80000000u) != 0u) {
+primitiveCount = int((iE >> 23u) & 0xffu);
+primitiveOffset = bufferOff + fc + int(iE & 0x007fffffu);
 }
 else {
 int i = node.x + bufferOff;
-vec4 vW = fetchBVH(i);
-vec4 sq = fetchBVH(i + 1);
-vec4 qI = fetchBVH(i + 2);
-vec4 oF = fetchBVH(i + 3);
-vec3 nodeMin = vW.xyz;
-vec3 ea = vec3(vW.w, sq.x, sq.y);
-uint mF[4] = uint[4](
-floatBitsToUint(sq.z), floatBitsToUint(sq.w),
-floatBitsToUint(qI.x), floatBitsToUint(qI.y)
+vec4 xS = fetchBVH(i);
+vec4 tR = fetchBVH(i + 1);
+vec4 rW = fetchBVH(i + 2);
+vec4 pL = fetchBVH(i + 3);
+vec3 nodeMin = xS.xyz;
+vec3 ej = vec3(xS.w, tR.x, tR.y);
+uint nH[4] = uint[4](
+floatBitsToUint(tR.z), floatBitsToUint(tR.w),
+floatBitsToUint(rW.x), floatBitsToUint(rW.y)
 );
-vec4 wV = eX(floatBitsToUint(qI.z));
-vec4 wW = eX(floatBitsToUint(qI.w));
-vec4 wX = eX(floatBitsToUint(oF.x));
-vec4 wR = eX(floatBitsToUint(oF.y));
-vec4 wS = eX(floatBitsToUint(oF.z));
-vec4 wU = eX(floatBitsToUint(oF.w));
-vec3 gb = eW(d);
-vec3 padding = eU(gb, o);
-vec4 tx0 = (nodeMin.x + wV * ea.x - o.x) * gb.x;
-vec4 tx1 = (nodeMin.x + wR * ea.x - o.x) * gb.x;
-vec4 ty0 = (nodeMin.y + wW * ea.y - o.y) * gb.y;
-vec4 ty1 = (nodeMin.y + wS * ea.y - o.y) * gb.y;
-vec4 tz0 = (nodeMin.z + wX * ea.z - o.z) * gb.z;
-vec4 tz1 = (nodeMin.z + wU * ea.z - o.z) * gb.z;
-vec4 ns = max(max((min(tx0, tx1) - vec4(padding.x)), (min(ty0, ty1) - vec4(padding.y))), (min(tz0, tz1) - vec4(padding.z)));
-vec4 sW = min(min((max(tx0, tx1) + vec4(padding.x)), (max(ty0, ty1) + vec4(padding.y))), (max(tz0, tz1) + vec4(padding.z)));
-vec4 distance = mix(sW, ns, greaterThan(ns, vec4(0.0)));
-distance = mix(distance, vec4(-1.0), greaterThan(ns, min(sW, vec4(tMax))));
-float eY[4] = float[4](distance.x, distance.y, distance.z, distance.w);
+vec4 yV = fh(floatBitsToUint(rW.z));
+vec4 yW = fh(floatBitsToUint(rW.w));
+vec4 yX = fh(floatBitsToUint(pL.x));
+vec4 yS = fh(floatBitsToUint(pL.y));
+vec4 yT = fh(floatBitsToUint(pL.z));
+vec4 yU = fh(floatBitsToUint(pL.w));
+vec3 invD = traversalInverse(d);
+vec3 padding = ff(invD, o);
+vec4 tx0 = (nodeMin.x + yV * ej.x - o.x) * invD.x;
+vec4 tx1 = (nodeMin.x + yS * ej.x - o.x) * invD.x;
+vec4 ty0 = (nodeMin.y + yW * ej.y - o.y) * invD.y;
+vec4 ty1 = (nodeMin.y + yT * ej.y - o.y) * invD.y;
+vec4 tz0 = (nodeMin.z + yX * ej.z - o.z) * invD.z;
+vec4 tz1 = (nodeMin.z + yU * ej.z - o.z) * invD.z;
+vec4 ot = max(max((min(tx0, tx1) - vec4(padding.x)), (min(ty0, ty1) - vec4(padding.y))), (min(tz0, tz1) - vec4(padding.z)));
+vec4 uA = min(min((max(tx0, tx1) + vec4(padding.x)), (max(ty0, ty1) + vec4(padding.y))), (max(tz0, tz1) + vec4(padding.z)));
+vec4 distance = mix(uA, ot, greaterThan(ot, vec4(0.0)));
+distance = mix(distance, vec4(-1.0), greaterThan(ot, min(uA, vec4(tMax))));
+float fi[4] = float[4](distance.x, distance.y, distance.z, distance.w);
 
 
 for (int order = 0; order < 4; order++) {
-float dW = -1.0;
-int ej = -1;
+float ee = -1.0;
+int eq = -1;
 for (int child = 0; child < 4; child++) {
-if (mF[child] != 0u && eY[child] > dW) {
-dW = eY[child];
-ej = child;
+if (nH[child] != 0u && fi[child] > ee) {
+ee = fi[child];
+eq = child;
 }
 }
-if (ej >= 0 && dW > 0.0) {
-uint ol = mF[ej];
-int jO = int(ol);
-if ((ol & 0x80000000u) == 0u) {
-jO += eQ;
+if (eq >= 0 && ee > 0.0) {
+uint ph = nH[eq];
+int kv = int(ph);
+if ((ph & 0x80000000u) == 0u) {
+kv += fc;
 }
-stack[++top] = ivec2(jO, bK);
-eY[ej] = -1.0;
+stack[++top] = ivec2(kv, bK);
+fi[eq] = -1.0;
 }
 }
 }
@@ -258,14 +258,14 @@ eY[ej] = -1.0;
 else {
 int i = node.x + bufferOff;
 vec4 r1 = fetchBVH(i);
-int bf = floatBitsToInt(r1.w);
+int bh = floatBitsToInt(r1.w);
 
-if (bf >= 0) {
+if (bh >= 0) {
 int left = node.x + 2;
-int right = bf + eQ;
+int right = bh + fc;
 
-float d0 = cH(o, d, fetchBVH(left + bufferOff).xyz, fetchBVH(left + 1 + bufferOff).xyz, tMax);
-float d1 = cH(o, d, fetchBVH(right + bufferOff).xyz, fetchBVH(right + 1 + bufferOff).xyz, tMax);
+float d0 = bP(o, d, fetchBVH(left + bufferOff).xyz, fetchBVH(left + 1 + bufferOff).xyz, tMax);
+float d1 = bP(o, d, fetchBVH(right + bufferOff).xyz, fetchBVH(right + 1 + bufferOff).xyz, tMax);
 if (d0 > d1) {
 if (d0 > 0.0) {
 stack[++top] = ivec2(left, bK);
@@ -285,13 +285,16 @@ stack[++top] = ivec2(left, bK);
 }
 else {
 if (bK < 0) {
-for (int k = 0; k < -bf; k++) {
+for (int k = 0; k < -bh; k++) {
 vec4 r3 = fetchBVH(i + 2 + k);
+int box = floatBitsToInt(r3.w);
+if (bP(ray.o, ray.d, fetchBVH(box).xyz, fetchBVH(box + 1).xyz, tMax) > 0.0) {
 stack[++top] = ivec2(floatBitsToInt(r3.y), floatBitsToInt(r3.x));
 }
 }
+}
 else {
-primitiveCount = -bf;
+primitiveCount = -bh;
 primitiveOffset = i + 2;
 }
 }
@@ -303,31 +306,31 @@ ${intersectTriangles}
 `;
 }
 export default createShaderMixin({
+    defines: {
+        VERTEX_ATTRIBUTE_STRIDE: 2
+    },
     uniforms: {
         positionBuffer: Shader.uniform('sampler2D'),
-        normalBuffer: Shader.uniform('sampler2D'),
-        tangentBuffer: Shader.uniform('sampler2D'),
-        uvBuffer: Shader.uniform('sampler2D'),
-        uv1Buffer: Shader.uniform('sampler2D'),
-        colorBuffer: Shader.uniform('sampler2D'),
-        materialIndexBuffer: Shader.uniform('sampler2D'),
+        // Interleaved shading attributes; layout in webglVertexAttributes.ts.
+        vertexAttributes: Shader.uniform('sampler2D'),
+        // TLAS, BLASes and then the instance records.
         bvhBuffer: Shader.uniform('sampler2D'),
-        instancesBuffer: Shader.uniform('sampler2D'),
         bvhColumns: Shader.uniform('int', 1),
         vertexColumns: Shader.uniform('int', 1),
-        instancesColumns: Shader.uniform('int', 1),
+        vertexAttributeColumns: Shader.uniform('int', 1),
+        instancesOffset: Shader.uniform('int', 0),
         blasOffset: Shader.uniform('int'),
         isEmptyScene: Shader.uniform('bool')
     },
     main: glsl `
 
-struct fw {
+struct fL {
 vec3 p0;
 vec3 p1;
 vec3 p2;
 };
 #ifdef NUM_LIGHTS
-float hs(float rad, vec3 pos, Ray r) {
+float hM(float rad, vec3 pos, Ray r) {
 vec3 op = pos - r.o;
 float eps = 0.001;
 float b = dot(op, r.d);
@@ -347,7 +350,7 @@ return t2;
 return INF;
 }
 
-float nV(vec3 pos, vec3 u, vec3 v, vec3 n, Ray r) {
+float oU(vec3 pos, vec3 u, vec3 v, vec3 n, Ray r) {
 float dt = dot(r.d, n);
 float t = (dot(n, pos) - dot(n, r.o)) / dt;
 
@@ -376,7 +379,7 @@ float oneMinusCosThetaMax = sinThetaMaxSquared / (1.0 + cosThetaMax);
 return 1.0 / max(TWOPI * oneMinusCosThetaMax, 1e-8);
 }
 
-float mX(Ray r, Light light, out LightSampleRec lightSample) {
+float iB(Ray r, Light light, out LightSampleRec lightSample) {
 
 
 
@@ -393,15 +396,15 @@ vec3 normal = normalize(cross(u, v));
 if (dot(normal, r.d) > 0.) 
 return INF;
 
-d = nV(pos, u, v, normal, r);
+d = oU(pos, u, v, normal, r);
 if (d < 0.) {
 d = INF;
 }
 else {
-float ct = dot(-r.d, normal);
+float cv = dot(-r.d, normal);
 float spread = light.spread;
-if (ct >= cos(spread * PI / 2.0)) {
-lightSample.pdf = (d * d) / (light.area * ct);
+if (cv >= cos(spread * PI / 2.0)) {
+lightSample.pdf = (d * d) / (light.area * cv);
 lightSample.emission = light.emission;
 }
 }
@@ -410,10 +413,10 @@ lightSample.emission = light.emission;
 #ifdef HAS_SPHERE_LIGHTS
 if (type == SPHERE_LIGHT) {
 if (light.radius <= 0.0) {
-vec3 kz = light.position - r.o;
-float projection = dot(kz, r.d);
-vec3 bd = kz - projection * r.d;
-if (projection <= EPS || dot(bd, bd) > 1e-8) return INF;
+vec3 lj = light.position - r.o;
+float projection = dot(lj, r.d);
+vec3 bf = lj - projection * r.d;
+if (projection <= EPS || dot(bf, bf) > 1e-8) return INF;
 d = projection;
 lightSample.direction = r.d;
 lightSample.normal = -r.d;
@@ -421,7 +424,7 @@ lightSample.dist = d;
 lightSample.pdf = 1.0;
 lightSample.isDelta = true;
 } else {
-d = hs(light.radius, light.position, r);
+d = hM(light.radius, light.position, r);
 lightSample.isDelta = false;
 }
 if (d < 0.) {
@@ -440,10 +443,10 @@ lightSample.emission = light.radius <= 0.0
 #ifdef HAS_SPOT_LIGHTS
 if (type == SPOT_LIGHT) {
 if (light.radius <= 0.0) {
-vec3 kz = light.position - r.o;
-float projection = dot(kz, r.d);
-vec3 bd = kz - projection * r.d;
-if (projection <= EPS || dot(bd, bd) > 1e-8) return INF;
+vec3 lj = light.position - r.o;
+float projection = dot(lj, r.d);
+vec3 bf = lj - projection * r.d;
+if (projection <= EPS || dot(bf, bf) > 1e-8) return INF;
 d = projection;
 lightSample.direction = r.d;
 lightSample.normal = -r.d;
@@ -451,17 +454,17 @@ lightSample.dist = d;
 lightSample.pdf = 1.0;
 lightSample.isDelta = true;
 } else {
-d = hs(light.radius, light.position, r);
+d = hM(light.radius, light.position, r);
 lightSample.isDelta = false;
 }
 if (d != INF) {
-vec3 xv = r.o + d * r.d;
-float ne = evaluateSpotEmissionMultiplier(light, normalize(r.o - xv));
-if (ne <= 0.0) return INF;
+vec3 zA = r.o + d * r.d;
+float od = evaluateSpotEmissionMultiplier(light, normalize(r.o - zA));
+if (od <= 0.0) return INF;
 if (light.radius > 0.0) {
 lightSample.pdf = sphereSolidAnglePdf(r.o, light.position, light.radius);
 }
-lightSample.emission = light.emission * ne * (light.radius <= 0.0
+lightSample.emission = light.emission * od * (light.radius <= 0.0
 ? 1.0 / max(d * d, 1e-8)
 : 1.0);
 }
@@ -472,7 +475,7 @@ return d;
 }
 #endif
 
-int fh(ivec4 m0, ivec4 m1, ivec4 m2, vec3 bary) {
+int fv(ivec4 m0, ivec4 m1, ivec4 m2, vec3 bary) {
 vec2 m00 = af(m0.x);
 vec2 m01 = af(m0.y);
 vec2 m02 = af(m0.z);
@@ -492,86 +495,102 @@ float w20 = m20.y * bary.z;
 float w21 = m21.y * bary.z;
 float w22 = m22.y * bary.z;
 
-float hd = w00 + w01 + w02 + w10 + w11 + w12 + w20 + w21 + w22;
+float hB = w00 + w01 + w02 + w10 + w11 + w12 + w20 + w21 + w22;
 
-w00 = w00 / hd;
-w01 = w01 / hd + w00;
-w02 = w02 / hd + w01;
-w10 = w10 / hd + w02;
-w11 = w11 / hd + w10;
-w12 = w12 / hd + w11;
-w20 = w20 / hd + w12;
-w21 = w21 / hd + w20;
-w22 = w22 / hd + w21;
+w00 = w00 / hB;
+w01 = w01 / hB + w00;
+w02 = w02 / hB + w01;
+w10 = w10 / hB + w02;
+w11 = w11 / hB + w10;
+w12 = w12 / hB + w11;
+w20 = w20 / hB + w12;
+w21 = w21 / hB + w20;
+w22 = w22 / hB + w21;
 
-float rand = ek(0u);
-float nW;
+float rand = cm(0u);
+float oV;
 if (rand < w00) {
-nW = m00.x;
+oV = m00.x;
 } else if (rand < w01) {
-nW = m01.x;
+oV = m01.x;
 } else if (rand < w02) {
-nW = m02.x;
+oV = m02.x;
 } else if (rand < w10) {
-nW = m10.x;
+oV = m10.x;
 } else if (rand < w11) {
-nW = m11.x;
+oV = m11.x;
 } else if (rand < w12) {
-nW = m12.x;
+oV = m12.x;
 } else if (rand < w20) {
-nW = m20.x;
+oV = m20.x;
 } else if (rand < w21) {
-nW = m21.x;
+oV = m21.x;
 } else {
-nW = m22.x;
+oV = m22.x;
 }
-return int(nW);
+return int(oV);
 }
 
 
 vec4 fetchIns(int i);
 
+vec4 bs(int vertex, int slot) {
+return fetchData(vertexAttributes, vertex * VERTEX_ATTRIBUTE_STRIDE + slot, vertexAttributeColumns);
+}
+vec3 fetchVertexNormal(int vertex) { return bs(vertex, 0).xyz; }
+vec2 er(int vertex) { return bs(vertex, 1).xy; }
+ivec4 dT(int vertex) {
+return ivec4(bs(vertex, 0).w, bs(vertex, 1).zw, 0.0);
+}
+vec2 dK(int vertex) {
+#ifdef USE_UV1
+return bs(vertex, VERTEX_UV1_SLOT).xy;
+#else
+return vec2(0.0);
+#endif
+}
+#ifdef USE_TANGENT
+vec4 ef(int vertex) { return bs(vertex, VERTEX_TANGENT_SLOT); }
+#endif
 
-mat2x3 ez(vec3 normal, vec3 dp1, vec3 dp2, vec2 oM, vec2 oN, vec4 aS) {
-float determinant = oM.x * oN.y - oM.y * oN.x;
-vec3 t = aS.xyz;
+
+mat2x3 eI(vec3 normal, vec3 dp1, vec3 dp2, vec2 pU, vec2 pV, vec4 aT) {
+float determinant = pU.x * pV.y - pU.y * pV.x;
+vec3 t = aT.xyz;
 vec3 b = vec3(0.0);
 
-float rw = abs(aS.w) > 1.5 ? 0.0 : aS.w;
+float sO = abs(aT.w) > 1.5 ? 0.0 : aT.w;
 if (abs(determinant) > 1e-12) {
-vec3 dt = (dp1 * oN.y - dp2 * oM.y) / determinant;
-vec3 db = (dp2 * oM.x - dp1 * oN.x) / determinant;
+vec3 dt = (dp1 * pV.y - dp2 * pU.y) / determinant;
+vec3 db = (dp2 * pU.x - dp1 * pV.x) / determinant;
 if (dot(t, t) < 1e-8) { t = dt; }
-if (rw == 0.0) { rw = (dot(cross(normal, dt), db) < 0.0 ? -1.0 : 1.0); }
+if (sO == 0.0) { sO = (dot(cross(normal, dt), db) < 0.0 ? -1.0 : 1.0); }
 }
 t -= normal * dot(normal, t);
 if (dot(t, t) < 1e-8) {
-float ik = (normal.z >= 0.0 ? 1.0 : -1.0);
-float cv = -1.0 / (ik + normal.z);
-float xy = normal.x * normal.y * cv;
-t = vec3(1.0 + ik * normal.x * normal.x * cv, ik * xy, -ik * normal.x);
-b = vec3(xy, ik + normal.y * normal.y * cv, -normal.y);
+float iP = (normal.z >= 0.0 ? 1.0 : -1.0);
+float cA = -1.0 / (iP + normal.z);
+float xy = normal.x * normal.y * cA;
+t = vec3(1.0 + iP * normal.x * normal.x * cA, iP * xy, -iP * normal.x);
+b = vec3(xy, iP + normal.y * normal.y * cA, -normal.y);
 }
 else { t = normalize(t); b = normalize(cross(normal, t)); }
-return mat2x3(t, b * (rw < 0.0 ? -1.0 : 1.0));
+return mat2x3(t, b * (sO < 0.0 ? -1.0 : 1.0));
 }
 
-void hn(inout SurfaceInteraction si, fw tri, vec3 barycentric, ivec3 index) {
+void hI(inout SurfaceInteraction si, fL tri, vec3 barycentric, ivec3 index) {
 si.hit = true;
 si.position = barycentric.x * tri.p0 + barycentric.y * tri.p1 + barycentric.z * tri.p2;
-ivec2 i0 = gM(index.x, vertexColumns);
-ivec2 i1 = gM(index.y, vertexColumns);
-ivec2 i2 = gM(index.z, vertexColumns);
 si.vertexColor = vec3(1.0);
 #ifdef USE_VERTEX_COLOR
-si.vertexColor = barycentric.x * texelFetch(colorBuffer, i0, 0).rgb
-+ barycentric.y * texelFetch(colorBuffer, i1, 0).rgb
-+ barycentric.z * texelFetch(colorBuffer, i2, 0).rgb;
+si.vertexColor = barycentric.x * bs(index.x, VERTEX_COLOR_SLOT).rgb
++ barycentric.y * bs(index.y, VERTEX_COLOR_SLOT).rgb
++ barycentric.z * bs(index.z, VERTEX_COLOR_SLOT).rgb;
 #endif
 
-vec3 n0 = texelFetch(normalBuffer, i0, 0).xyz;
-vec3 n1 = texelFetch(normalBuffer, i1, 0).xyz;
-vec3 n2 = texelFetch(normalBuffer, i2, 0).xyz;
+vec3 n0 = fetchVertexNormal(index.x);
+vec3 n1 = fetchVertexNormal(index.y);
+vec3 n2 = fetchVertexNormal(index.z);
 vec3 normal = normalize(barycentric.x * n0 + barycentric.y * n1 + barycentric.z * n2);
 si.vertexNormal = normal;
 
@@ -580,14 +599,14 @@ si.vertexNormal = normal;
 
 
 vec3 geometricNormal = normalize(cross(tri.p1 - tri.p0, tri.p2 - tri.p0));
-vec3 kw = barycentric.x * n0 + barycentric.y * n1 + barycentric.z * n2;
+vec3 le = barycentric.x * n0 + barycentric.y * n1 + barycentric.z * n2;
 float a = dot(n2 - n0, tri.p0 - tri.p2);
 float b = dot(n2 - n1, tri.p1 - tri.p2);
 float c = dot(n1 - n0, tri.p1 - tri.p0);
 float h = a * barycentric.x * (barycentric.x - 1.0)
 + (a + b + c) * barycentric.x * barycentric.y
 + b * barycentric.y * (barycentric.y - 1.0);
-if (dot(kw, geometricNormal) > 0.0) {
+if (dot(le, geometricNormal) > 0.0) {
 float h0 = max(max(dot(tri.p1 - tri.p0, n0), dot(tri.p2 - tri.p0, n0)), 0.0);
 float h1 = max(max(dot(tri.p0 - tri.p1, n1), dot(tri.p2 - tri.p1, n1)), 0.0);
 float h2 = max(max(dot(tri.p0 - tri.p2, n2), dot(tri.p1 - tri.p2, n2)), 0.0);
@@ -605,120 +624,119 @@ 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 = kw * h;
-
+si.shadowOffset = le * h;
 
-ivec4 m0 = ivec4(texelFetch(materialIndexBuffer, i0, 0));
-ivec4 m1 = ivec4(texelFetch(materialIndexBuffer, i1, 0));
-ivec4 m2 = ivec4(texelFetch(materialIndexBuffer, i2, 0));
+ivec4 m0 = dT(index.x);
+ivec4 m1 = dT(index.y);
+ivec4 m2 = dT(index.z);
 
 if (m0.x > 0 || m1.x > 0 || m2.x > 0) {
-si.materialIndex = fh(m0, m1, m2, barycentric);
+si.materialIndex = fv(m0, m1, m2, barycentric);
 }
 
-vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
-vec2 uv1 = texelFetch(uvBuffer, i1, 0).xy;
-vec2 uv2 = texelFetch(uvBuffer, i2, 0).xy;
+vec2 uv0 = er(index.x);
+vec2 uv1 = er(index.y);
+vec2 uv2 = er(index.z);
 #if defined(NUM_DIFFUSE_MAPS) || defined(NUM_NORMAL_MAPS) || defined(NUM_METALLIC_ROUGHNESS_MAPS) || defined(NUM_ROUGHNESS_MAPS) || defined(NUM_METALNESS_MAPS) || defined(NUM_EMISSIVE_MAPS) || defined(NUM_THIN_FILM_MAPS) || defined(NUM_ANISOTROPY_MAPS) || defined(NUM_COAT_NORMAL_MAPS)
 vec2 uv = barycentric.x * uv0 + barycentric.y * uv1 + barycentric.z * uv2;
-vec2 secondaryUv0 = texelFetch(uv1Buffer, i0, 0).xy;
-vec2 secondaryUv1 = texelFetch(uv1Buffer, i1, 0).xy;
-vec2 secondaryUv2 = texelFetch(uv1Buffer, i2, 0).xy;
-vec2 cQ = barycentric.x * secondaryUv0
+vec2 secondaryUv0 = dK(index.x);
+vec2 secondaryUv1 = dK(index.y);
+vec2 secondaryUv2 = dK(index.z);
+vec2 cT = barycentric.x * secondaryUv0
 + barycentric.y * secondaryUv1
 + barycentric.z * secondaryUv2;
 
-vec4 mb = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_UV);
-int uvProjection = int(mb.r);
+vec4 mU = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_UV);
+int uvProjection = int(mU.r);
 if (uvProjection > 0) {
-uv = gE(uvProjection, si.position, normal, mb.zw);
-cQ = uv;
+uv = gW(uvProjection, si.position, normal, mU.zw);
+cT = uv;
 } else {
-uv *= mb.zw;
-cQ *= mb.zw;
+uv *= mU.zw;
+cT *= mU.zw;
 }
 #else
 vec2 uv = vec2(0.0);
-vec2 cQ = vec2(0.0);
+vec2 cT = vec2(0.0);
 #endif
 
 #ifdef USE_COAT_NORMAL_FRAMES
 si.coatNormal = normal;
 #ifdef NUM_COAT_NORMAL_MAPS
-float pH = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_LAYER_MAPS).z;
+float qK = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_LAYER_MAPS).z;
 #endif
 #endif
 vec4 directionMaterial = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_DIRECTION_NORMAL);
-float uQ = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_FILM_TEXTURE).x;
-if (directionMaterial.x >= 0.0 || directionMaterial.y != 0.0 || uQ != 0.0
+float wB = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_FILM_TEXTURE).x;
+if (directionMaterial.x >= 0.0 || directionMaterial.y != 0.0 || wB != 0.0
 #ifdef USE_COAT_NORMAL_FRAMES
 || fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_COAT).w > 0.0
 #endif
 ) {
 vec3 dp1 = tri.p0 - tri.p2;
 vec3 dp2 = tri.p1 - tri.p2;
-vec4 aS = vec4(0.0);
+vec4 aT = vec4(0.0);
 #ifdef USE_TANGENT
-aS = barycentric.x * texelFetch(tangentBuffer, i0, 0)
-+ barycentric.y * texelFetch(tangentBuffer, i1, 0)
-+ barycentric.z * texelFetch(tangentBuffer, i2, 0);
+aT = barycentric.x * ef(index.x)
++ barycentric.y * ef(index.y)
++ barycentric.z * ef(index.z);
 #endif
-mat2x3 eL = ez(normal, dp1, dp2, uv0-uv2, uv1-uv2, aS);
-si.geometryTangent = eL[0];
-si.geometryBitangent = eL[1];
-mat2x3 fX = eL;
+mat2x3 eV = eI(normal, dp1, dp2, uv0-uv2, uv1-uv2, aT);
+si.geometryTangent = eV[0];
+si.geometryBitangent = eV[1];
+mat2x3 gp = eV;
 #ifdef NUM_ANISOTROPY_MAPS
 if (directionMaterial.x >= 0.0 && mapUvChannel[int(directionMaterial.x)] > 0.5) {
-vec2 a = texelFetch(uv1Buffer, i0, 0).xy;
-vec2 b = texelFetch(uv1Buffer, i1, 0).xy;
-vec2 c = texelFetch(uv1Buffer, i2, 0).xy;
-fX = ez(normal, dp1, dp2, a-c, b-c, (aS.w == 0.0 ? vec4(0.0) : aS));
+vec2 a = dK(index.x);
+vec2 b = dK(index.y);
+vec2 c = dK(index.z);
+gp = eI(normal, dp1, dp2, a-c, b-c, (aT.w == 0.0 ? vec4(0.0) : aT));
 }
 #endif
-si.directionTangent = fX[0];
-si.directionBitangent = fX[1];
+si.directionTangent = gp[0];
+si.directionBitangent = gp[1];
 #if defined(USE_COAT_NORMAL_FRAMES) && defined(NUM_COAT_NORMAL_MAPS)
-if (pH >= 0.0) {
-int kP = int(pH);
-mat2x3 mG = eL;
-if (mapUvChannel[kP] > 0.5) {
-vec2 a = texelFetch(uv1Buffer, i0, 0).xy;
-vec2 b = texelFetch(uv1Buffer, i1, 0).xy;
-vec2 c = texelFetch(uv1Buffer, i2, 0).xy;
-mG = ez(normal, dp1, dp2, a-c, b-c,
-aS.w == 0.0 ? vec4(0.0) : aS);
-}
-vec3 decoded = sampleMap(kP, getMapUvSet(kP, uv, cQ)).rgb * 2.0 - 1.0;
-vec3 mapped = mG[0] * decoded.x * directionMaterial.z
-+ mG[1] * decoded.y * directionMaterial.w + normal * decoded.z;
+if (qK >= 0.0) {
+int lI = int(qK);
+mat2x3 nI = eV;
+if (mapUvChannel[lI] > 0.5) {
+vec2 a = dK(index.x);
+vec2 b = dK(index.y);
+vec2 c = dK(index.z);
+nI = eI(normal, dp1, dp2, a-c, b-c,
+aT.w == 0.0 ? vec4(0.0) : aT);
+}
+vec3 decoded = sampleMap(lI, getMapUvSet(lI, uv, cT)).rgb * 2.0 - 1.0;
+vec3 mapped = nI[0] * decoded.x * directionMaterial.z
++ nI[1] * decoded.y * directionMaterial.w + normal * decoded.z;
 if (dot(mapped, mapped) > 1e-12) { si.coatNormal = normalize(mapped); }
 }
 #endif
 }
 
 si.uv = uv;
-si.uv1 = cQ;
+si.uv1 = cT;
 
 #ifdef USE_TANGENT
-vec3 t0 = texelFetch(tangentBuffer, i0, 0).xyz;
-vec3 t1 = texelFetch(tangentBuffer, i1, 0).xyz;
-vec3 t2 = texelFetch(tangentBuffer, i2, 0).xyz;
+vec3 t0 = ef(index.x).xyz;
+vec3 t1 = ef(index.y).xyz;
+vec3 t2 = ef(index.z).xyz;
 si.tangent = normalize(barycentric.x * t0 + barycentric.y * t1 + barycentric.z * t2);
 #endif
 #if defined(NUM_NORMAL_MAPS)
 vec3 dp1 = tri.p0 - tri.p2;
 vec3 dp2 = tri.p1 - tri.p2;
-vec2 oM = uv0 - uv2;
-vec2 oN = uv1 - uv2;
-si.normal = oS(
-si.materialIndex, uv, normal, dp1, dp2, oM, oN
+vec2 pU = uv0 - uv2;
+vec2 pV = uv1 - uv2;
+si.normal = qa(
+si.materialIndex, uv, normal, dp1, dp2, pU, pV
 );
 #else
 si.normal = normal;
 #endif
 }
 
-struct ev {
+struct eF {
 float t;
 vec3 barycentric;
 };
@@ -728,32 +746,32 @@ return fetchData(bvhBuffer, i, bvhColumns);
 }
 
 vec4 fetchIns(int i) {
-return fetchData(instancesBuffer, i, instancesColumns);
+return fetchData(bvhBuffer, instancesOffset + i, bvhColumns);
 }
 
 
 
-uint rT(uint value) {
+uint tl(uint value) {
 value -= (value >> 1u) & 0x55555555u;
 value = (value & 0x33333333u) + ((value >> 2u) & 0x33333333u);
 value = (value + (value >> 4u)) & 0x0f0f0f0fu;
 return (value * 0x01010101u) >> 24u;
 }
 
-vec4 eX(uint mo) {
+vec4 fh(uint nk) {
 return vec4(
-float(mo & 0xffu),
-float((mo >> 8u) & 0xffu),
-float((mo >> 16u) & 0xffu),
-float(mo >> 24u)
+float(nk & 0xffu),
+float((nk >> 8u) & 0xffu),
+float((nk >> 16u) & 0xffu),
+float(nk >> 24u)
 ) / 255.0;
 }
 
 
 
 
-ev fM(vec3 o, vec3 d, fw tri, int dZ, vec3 gH, float tMax) {
-ev ti;
+eF fY(vec3 o, vec3 d, fL tri, int ei, vec3 gZ, float tMax) {
+eF ti;
 
 
 vec3 p0t = tri.p0 - o;
@@ -761,20 +779,20 @@ vec3 p1t = tri.p1 - o;
 vec3 p2t = tri.p2 - o;
 
 
-if (dZ == 0) {
+if (ei == 0) {
 p0t = p0t.yzx;
 p1t = p1t.yzx;
 p2t = p2t.yzx;
-} else if (dZ == 1) {
+} else if (ei == 1) {
 p0t = p0t.zxy;
 p1t = p1t.zxy;
 p2t = p2t.zxy;
 }
 
 
-p0t.xy += gH.xy * p0t.z;
-p1t.xy += gH.xy * p1t.z;
-p2t.xy += gH.xy * p2t.z;
+p0t.xy += gZ.xy * p0t.z;
+p1t.xy += gZ.xy * p1t.z;
+p2t.xy += gZ.xy * p2t.z;
 
 
 vec3 e = vec3(
@@ -784,9 +802,9 @@ p0t.x * p1t.y - p0t.y * p1t.x
 );
 
 #ifdef DEPTH_QUERY_MODE
-vec3 sB = abs(vec3(p1t.x * p2t.y, p2t.x * p0t.y, p0t.x * p1t.y))
+vec3 ua = abs(vec3(p1t.x * p2t.y, p2t.x * p0t.y, p0t.x * p1t.y))
 + abs(vec3(p1t.y * p2t.x, p2t.y * p0t.x, p0t.y * p1t.x));
-e = mix(e, vec3(0.0), lessThanEqual(abs(e), sB * 1e-6));
+e = mix(e, vec3(0.0), lessThanEqual(abs(e), ua * 1e-6));
 #endif
 
 
@@ -799,19 +817,19 @@ float det = e.x + e.y + e.z;
 if (det == 0.0) return ti;
 #endif
 
-p0t.z *= gH.z;
-p1t.z *= gH.z;
-p2t.z *= gH.z;
-float sp = (e.x * p0t.z + e.y * p1t.z + e.z * p2t.z);
+p0t.z *= gZ.z;
+p1t.z *= gZ.z;
+p2t.z *= gZ.z;
+float tQ = (e.x * p0t.z + e.y * p1t.z + e.z * p2t.z);
 
 
-if (abs(sp) > abs(tMax * det)) {
+if (abs(tQ) > abs(tMax * det)) {
 return ti;
 }
 
-float tP = 1. / det;
-ti.t = sp * tP;
-ti.barycentric = e * tP;
+float vu = 1. / det;
+ti.t = tQ * vu;
+ti.barycentric = e * vu;
 
 return ti;
 }
@@ -822,7 +840,7 @@ vec3 max;
 };
 
 
-vec3 eW(vec3 d) {
+vec3 traversalInverse(vec3 d) {
 #ifdef DEPTH_QUERY_MODE
 vec3 ah = mix(d, mix(vec3(1e-20), vec3(-1e-20), lessThan(d, vec3(0.0))), lessThan(abs(d), vec3(1e-20)));
 return 1.0 / ah;
@@ -831,9 +849,9 @@ return 1.0 / d;
 #endif
 }
 
-vec3 eU(vec3 gb, vec3 o) {
+vec3 ff(vec3 invD, vec3 o) {
 #ifdef DEPTH_QUERY_MODE
-return abs(gb) * max(vec3(1.0), abs(o)) * 1e-6;
+return abs(invD) * max(vec3(1.0), abs(o)) * 1e-6;
 #else
 return vec3(0.0);
 #endif
@@ -841,21 +859,21 @@ return vec3(0.0);
 
 
 
-float cH(vec3 o, vec3 d, vec3 xU, vec3 xT, float tMax) {
+float bP(vec3 o, vec3 d, vec3 zT, vec3 zS, float tMax) {
 #ifdef DEPTH_QUERY_MODE
-vec3 gb = eW(d);
-vec3 xg = (xU - o) * gb;
-vec3 xh = (xT - o) * gb;
-vec3 padding = eU(gb, o);
+vec3 invD = traversalInverse(d);
+vec3 zi = (zT - o) * invD;
+vec3 zj = (zS - o) * invD;
+vec3 padding = ff(invD, o);
 #else
-vec3 xg = (xU - o) / d;
-vec3 xh = (xT - o) / d;
+vec3 zi = (zT - o) / d;
+vec3 zj = (zS - o) / d;
 vec3 padding = vec3(0.0);
 #endif
-vec3 ns = min(xg, xh) - padding;
-vec3 sW = max(xg, xh) + padding;
-float t0 = max(ns.x, max(ns.y, ns.z));
-float t1 = min(sW.x, min(sW.y, sW.z));
+vec3 ot = min(zi, zj) - padding;
+vec3 uA = max(zi, zj) + padding;
+float t0 = max(ot.x, max(ot.y, ot.z));
+float t1 = min(uA.x, min(uA.y, uA.z));
 
 return (t0 > t1 || t0 > tMax) ? -1.0 : (t0 > 0.0 ? t0 : t1);
 }
@@ -864,12 +882,12 @@ int maxDimension(vec3 v) {
 return v.x > v.y ? (v.x > v.z ? 0 : 2) : (v.y > v.z ? 1 : 2);
 }
 
-vec3 mU(vec3 d, int dZ) {
+vec3 nV(vec3 d, int ei) {
 
 
-if (dZ == 0) {
+if (ei == 0) {
 return vec3(-d.y, -d.z, 1.0) / d.x;
-} else if (dZ == 1) {
+} else if (ei == 1) {
 return vec3(-d.z, -d.x, 1.0) / d.y;
 } else {
 return vec3(-d.x, -d.y, 1.0) / d.z;
@@ -884,25 +902,25 @@ if (isEmptyScene) {
 return;
 }
 
-int dZ = maxDimension(abs(ray.d));
-vec3 gH = mU(ray.d, dZ);
+int ei = maxDimension(abs(ray.d));
+vec3 gZ = nV(ray.d, ei);
 int am;
 
 ${getBVHTraverseCode(glsl `
 si.triangleIndex = triangleIndex;
-lR = tri;
-hz = faceNormal;
+mL = tri;
+hU = faceNormal;
 aw = hit.barycentric;
-hr = index;
-sd = true;
+hL = index;
+tD = true;
 am = bK;
 `)}
 
-if (sd) {
+if (tD) {
 si.instanceIndex = am;
 si.materialIndex = int(af(floatBitsToInt(fetchIns(am * ${INSTANCE_VEC4_STRIDE}).z)).x);
 
-hn(si, lR, aw, hr);
+hI(si, mL, aw, hL);
 mat4 l2w = mat4(
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 5),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 6),
@@ -915,12 +933,12 @@ fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 2),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 3),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 4)
 );
-mat3 aG = transpose(mat3(w2l));
-si.faceNormal = normalize(aG * hz);
-si.vertexNormal = normalize(aG * si.vertexNormal);
-si.normal = normalize(aG * si.normal);
+mat3 aH = transpose(mat3(w2l));
+si.faceNormal = normalize(aH * hU);
+si.vertexNormal = normalize(aH * si.vertexNormal);
+si.normal = normalize(aH * si.normal);
 #ifdef USE_COAT_NORMAL_FRAMES
-si.coatNormal = normalize(aG * si.coatNormal);
+si.coatNormal = normalize(aH * si.coatNormal);
 #endif
 #ifdef USE_TANGENT
 if (dot(si.tangent, si.tangent) > 0.0) {
@@ -945,17 +963,17 @@ return;
 }
 
 bool doubleSided = true;
-int dZ = maxDimension(abs(ray.d));
-vec3 gH = mU(ray.d, dZ);
+int ei = maxDimension(abs(ray.d));
+vec3 gZ = nV(ray.d, ei);
 int am;
 
 ${getBVHTraverseCode(glsl `
 si.triangleIndex = triangleIndex;
-lR = tri;
-hz = faceNormal;
+mL = tri;
+hU = faceNormal;
 aw = hit.barycentric;
-hr = index;
-sd = true;
+hL = index;
+tD = true;
 am = bK;
 `, glsl `
 if (bK != instanceIndex) {
@@ -963,10 +981,10 @@ continue;
 }
 `, 'maxDistance')}
 
-if (sd) {
+if (tD) {
 si.instanceIndex = am;
 si.materialIndex = int(af(floatBitsToInt(fetchIns(am * ${INSTANCE_VEC4_STRIDE}).z)).x);
-hn(si, lR, aw, hr);
+hI(si, mL, aw, hL);
 mat4 l2w = mat4(
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 5),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 6),
@@ -979,12 +997,12 @@ fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 2),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 3),
 fetchIns(am * ${INSTANCE_VEC4_STRIDE} + 4)
 );
-mat3 aG = transpose(mat3(w2l));
-si.faceNormal = normalize(aG * hz);
-si.vertexNormal = normalize(aG * si.vertexNormal);
-si.normal = normalize(aG * si.normal);
+mat3 aH = transpose(mat3(w2l));
+si.faceNormal = normalize(aH * hU);
+si.vertexNormal = normalize(aH * si.vertexNormal);
+si.normal = normalize(aH * si.normal);
 #ifdef USE_COAT_NORMAL_FRAMES
-si.coatNormal = normalize(aG * si.coatNormal);
+si.coatNormal = normalize(aH * si.coatNormal);
 #endif
 #ifdef USE_TANGENT
 if (dot(si.tangent, si.tangent) > 0.0) {
@@ -1000,35 +1018,32 @@ si.position = (l2w * vec4(si.position, 1.0)).xyz;
 }
 }
 
-bool intersectSceneShadow(in Ray ray, float lightDistance, int jR) {
+bool intersectSceneShadow(in Ray ray, float lightDistance, int kx) {
 
 if (isEmptyScene) {
 return false;
 }
 
 bool doubleSided = true;
-int dZ = maxDimension(abs(ray.d));
-vec3 gH = mU(ray.d, dZ);
+int ei = maxDimension(abs(ray.d));
+vec3 gZ = nV(ray.d, ei);
 
 ${getBVHTraverseCode(glsl `
 
 
-bool jA = true;
+bool kc = true;
 vec3 barycentric = hit.barycentric;
 vec3 position = barycentric.x * tri.p0 + barycentric.y * tri.p1 + barycentric.z * tri.p2;
 
-if (jA) {
+if (kc) {
 #ifdef HAS_ALPHA
 int materialIndex = int(af(floatBitsToInt(fetchIns(bK * ${INSTANCE_VEC4_STRIDE}).z)).x);
-ivec2 i0 = gM(index.x, vertexColumns);
-ivec2 i1 = gM(index.y, vertexColumns);
-ivec2 i2 = gM(index.z, vertexColumns);
-ivec4 m0 = ivec4(texelFetch(materialIndexBuffer, i0, 0));
-ivec4 m1 = ivec4(texelFetch(materialIndexBuffer, i1, 0));
-ivec4 m2 = ivec4(texelFetch(materialIndexBuffer, i2, 0));
+ivec4 m0 = dT(index.x);
+ivec4 m1 = dT(index.y);
+ivec4 m2 = dT(index.z);
 
 if (m0.x > 0 || m1.x > 0 || m2.x > 0) {
-materialIndex = fh(m0, m1, m2, barycentric);
+materialIndex = fv(m0, m1, m2, barycentric);
 }
 
 int ii = materialIndex * MAT_DATA_SIZE;
@@ -1038,16 +1053,16 @@ float alpha = fetchMatData(ii + MAT_ROW_ALPHA_FILM).r;
 if (mapAlphaUsage == ${mapAlphaUsageConstants.alpha + ''}) {
 int diffuseMapIndex = int(fetchMatData(ii + MAT_ROW_METALLIC_MAPS).z);
 if (diffuseMapIndex >= 0) {
-vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
-vec2 uv1 = texelFetch(uvBuffer, i1, 0).xy;
-vec2 uv2 = texelFetch(uvBuffer, i2, 0).xy;
+vec2 uv0 = er(index.x);
+vec2 uv1 = er(index.y);
+vec2 uv2 = er(index.z);
 vec2 uv = barycentric.x * uv0 + barycentric.y * uv1 + barycentric.z * uv2;
-vec4 mb = fetchMatData(ii + MAT_ROW_UV);
-int uvProjection = int(mb.r);
+vec4 mU = fetchMatData(ii + MAT_ROW_UV);
+int uvProjection = int(mU.r);
 if (uvProjection > 0) {
-uv = gE(uvProjection, position, hz, mb.zw);
+uv = gW(uvProjection, position, hU, mU.zw);
 } else {
-uv *= mb.zw;
+uv *= mU.zw;
 }
 
 alpha *= sampleMap(diffuseMapIndex, uv).a;
@@ -1057,36 +1072,36 @@ alpha *= sampleMap(diffuseMapIndex, uv).a;
 #ifdef NUM_ALPHA_MAPS
 int alphaMapIndex = int(fetchMatData(ii + MAT_ROW_STANDALONE_MAPS).z);
 if (alphaMapIndex >= 0) {
-vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
-vec2 uv1 = texelFetch(uvBuffer, i1, 0).xy;
-vec2 uv2 = texelFetch(uvBuffer, i2, 0).xy;
+vec2 uv0 = er(index.x);
+vec2 uv1 = er(index.y);
+vec2 uv2 = er(index.z);
 vec2 uv = barycentric.x * uv0 + barycentric.y * uv1 + barycentric.z * uv2;
-vec4 mb = fetchMatData(ii + MAT_ROW_UV);
-int uvProjection = int(mb.r);
+vec4 mU = fetchMatData(ii + MAT_ROW_UV);
+int uvProjection = int(mU.r);
 if (uvProjection > 0) {
-uv = gE(uvProjection, position, hz, mb.zw);
+uv = gW(uvProjection, position, hU, mU.zw);
 } else {
-uv *= mb.zw;
+uv *= mU.zw;
 }
 alpha *= sampleMap(alphaMapIndex, uv).g;
 }
 #endif
 float alphaCutoff = fetchMatData(ii + MAT_ROW_STANDALONE_MAPS).a;
-alpha = lJ(
+alpha = mB(
 alpha, int(fetchMatData(ii + MAT_ROW_ALPHA).w), alphaCutoff
 );
-if (ek(0u) > alpha) {
+if (cm(0u) > alpha) {
 
 tMax = min(max(lightDistance, 0.0), RAY_MAX_DISTANCE);
-jA = false;
+kc = false;
 }
 #endif
 }
-if (jA) {
+if (kc) {
 return true;
 }
 `, glsl `
-if (triangleIndex == jR) {
+if (triangleIndex == kx) {
 continue;
 }
 `, 'min(max(lightDistance, 0.0), RAY_MAX_DISTANCE)')}