@meimu/mutracer 0.1.4 → 0.2.0

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

@@ -6,31 +6,31 @@ import { INSTANCE_VEC4_STRIDE } from '../../bvhAccel.js';
 import { createShaderMixin, glsl, Shader } from '../../shader.js';
 function getBVHTraverseCode(handleHitCode, skipTriangleCode = '', initialTMax = 'RAY_MAX_DISTANCE') {
     const intersectTriangles = glsl `
-ivec3 vertexOffset = ivec3(floatBitsToInt(bv(cf * ${INSTANCE_VEC4_STRIDE}).y));
-for (int k = 0; k < bh; k++) {
+ivec3 vertexOffset = ivec3(floatBitsToInt(bF(bY * ${INSTANCE_VEC4_STRIDE}).y));
+for (int k = 0; k < primitiveCount; k++) {
 int triangleIndex = primitiveOffset + k;
 ${skipTriangleCode}
-vec4 r3 = dp(triangleIndex);
+vec4 r3 = dJ(triangleIndex);
 ivec3 index = floatBitsToInt(r3.xyz) + vertexOffset;
-fA tri = fA(
-kG(positionBuffer, index.x, vertexColumns).xyz,
-kG(positionBuffer, index.y, vertexColumns).xyz,
-kG(positionBuffer, index.z, vertexColumns).xyz
+fP tri = fP(
+lm(positionBuffer, index.x, vertexColumns).xyz,
+lm(positionBuffer, index.y, vertexColumns).xyz,
+lm(positionBuffer, index.z, vertexColumns).xyz
 );
-eC hit = fL(o, d, tri, ey, gP, tMax);
+eK hit = ga(o, d, tri, eq, gY, 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 kk = int(ae(
-floatBitsToInt(bv(cf * ${INSTANCE_VEC4_STRIDE}).z)
+int kP = int(ah(
+floatBitsToInt(bF(bY * ${INSTANCE_VEC4_STRIDE}).z)
 ).x);
-bool iN = fetchMatData(
-kk * MAT_DATA_SIZE + MAT_ROW_ALPHA
+bool jd = fetchMatData(
+kP * MAT_DATA_SIZE + MAT_ROW_ALPHA
 ).z > 0.5;
-if (!doubleSided && !iN && dot(faceNormal, d) > 0.0) {
+if (!doubleSided && !jd && dot(faceNormal, d) > 0.0) {
 continue;
 }
 
@@ -41,262 +41,262 @@ ${handleHitCode}
     return glsl `
 ivec2 stack[STACK_SIZE];
 int top = -1;
-int sO = -1;
+int ud = -1;
 float tMax = ${initialTMax};
 
 vec3 o = ray.o;
 vec3 d = ray.d;
 
-bool rb = false;
-vec3 hB;
-vec3 aB;
-ivec3 hw;
+bool sj = false;
+vec3 hS;
+vec3 aC;
+ivec3 hK;
 
-fA lC;
+fP mk;
 
 
 top = 0;
 stack[0] = ivec2(0, -1);
 
-int wr = 0;
-while(top >= 0 && wr < 5000) {
-wr++;
+int yk = 0;
+while(top >= 0 && yk < 5000) {
+yk++;
 ivec2 node = stack[top--];
 
-int cf = node.y;
+int bY = node.y;
 
-if (cf != sO) {
-if (cf < 0) {
+if (bY != ud) {
+if (bY < 0) {
 o = ray.o;
 d = ray.d;
 }
 else {
 mat4 w2l = mat4(
-bv(cf * ${INSTANCE_VEC4_STRIDE} + 1),
-bv(cf * ${INSTANCE_VEC4_STRIDE} + 2),
-bv(cf * ${INSTANCE_VEC4_STRIDE} + 3),
-bv(cf * ${INSTANCE_VEC4_STRIDE} + 4)
+bF(bY * ${INSTANCE_VEC4_STRIDE} + 1),
+bF(bY * ${INSTANCE_VEC4_STRIDE} + 2),
+bF(bY * ${INSTANCE_VEC4_STRIDE} + 3),
+bF(bY * ${INSTANCE_VEC4_STRIDE} + 4)
 );
 o = (w2l * vec4(ray.o, 1.0)).xyz;
 d = (w2l * vec4(ray.d, 0.0)).xyz;
 }
-ey = maxDimension(abs(d));
-gP = mv(d, ey);
-sO = cf;
+eq = maxDimension(abs(d));
+gY = nk(d, eq);
+ud = bY;
 }
 
-int bo = cf < 0 ? 0 : blasOffset;
-uint fJ = cf < 0 ? 0u : floatBitsToUint(bv(cf * ${INSTANCE_VEC4_STRIDE}).x);
-int gR = int(fJ >> 30u);
-int eW = int(fJ & 0x3fffffffu);
-int bh = 0;
+int bw = bY < 0 ? 0 : blasOffset;
+uint fY = bY < 0 ? 0u : floatBitsToUint(bF(bY * ${INSTANCE_VEC4_STRIDE}).x);
+int he = int(fY >> 30u);
+int fg = int(fY & 0x3fffffffu);
+int primitiveCount = 0;
 int primitiveOffset = 0;
 
-if (gR == 3) {
-uint hW = uint(node.x);
-if ((hW & 0x80000000u) != 0u) {
-bh = int((hW >> 23u) & 0xffu);
-primitiveOffset = bo + eW + int(hW & 0x007fffffu);
+if (he == 3) {
+uint iw = uint(node.x);
+if ((iw & 0x80000000u) != 0u) {
+primitiveCount = int((iw >> 23u) & 0xffu);
+primitiveOffset = bw + fg + int(iw & 0x007fffffu);
 }
 else {
-int i = node.x + bo;
-vec4 jA = dp(i);
-vec4 fF = dp(i + 1);
-vec4 mV = dp(i + 2);
-vec4 mW = dp(i + 3);
-vec4 mX = dp(i + 4);
-uint header = floatBitsToUint(jA.w);
-int rG = int(header << 24u) >> 24;
-int rH = int((header >> 8u) << 24u) >> 24;
-int rI = int((header >> 16u) << 24u) >> 24;
-vec3 hJ = vec3(
-uintBitsToFloat(uint(rG + 127) << 23u),
-uintBitsToFloat(uint(rH + 127) << 23u),
-uintBitsToFloat(uint(rI + 127) << 23u)
+int i = node.x + bw;
+vec4 jZ = dJ(i);
+vec4 fU = dJ(i + 1);
+vec4 nS = dJ(i + 2);
+vec4 nT = dJ(i + 3);
+vec4 nU = dJ(i + 4);
+uint header = floatBitsToUint(jZ.w);
+int sH = int(header << 24u) >> 24;
+int sI = int((header >> 8u) << 24u) >> 24;
+int sJ = int((header >> 16u) << 24u) >> 24;
+vec3 hZ = vec3(
+uintBitsToFloat(uint(sH + 127) << 23u),
+uintBitsToFloat(uint(sI + 127) << 23u),
+uintBitsToFloat(uint(sJ + 127) << 23u)
 );
-uint childBase = floatBitsToUint(fF.x);
-uint oQ = floatBitsToUint(fF.y);
+uint childBase = floatBitsToUint(fU.x);
+uint pM = floatBitsToUint(fU.y);
 uint internalMask = header >> 24u;
-uint ut = (d.x < 0.0 ? 4u : 0u) |
+uint vU = (d.x < 0.0 ? 4u : 0u) |
 (d.y < 0.0 ? 2u : 0u) |
 (d.z < 0.0 ? 1u : 0u);
-uint gI[12] = uint[12](
-floatBitsToUint(mV.x), floatBitsToUint(mV.y),
-floatBitsToUint(mV.z), floatBitsToUint(mV.w),
-floatBitsToUint(mW.x), floatBitsToUint(mW.y),
-floatBitsToUint(mW.z), floatBitsToUint(mW.w),
-floatBitsToUint(mX.x), floatBitsToUint(mX.y),
-floatBitsToUint(mX.z), floatBitsToUint(mX.w)
+uint gV[12] = uint[12](
+floatBitsToUint(nS.x), floatBitsToUint(nS.y),
+floatBitsToUint(nS.z), floatBitsToUint(nS.w),
+floatBitsToUint(nT.x), floatBitsToUint(nT.y),
+floatBitsToUint(nT.z), floatBitsToUint(nT.w),
+floatBitsToUint(nU.x), floatBitsToUint(nU.y),
+floatBitsToUint(nU.z), floatBitsToUint(nU.w)
 );
 
 for (int order = 7; order >= 0; order--) {
-uint slot = uint(order) ^ ut;
-uint sE = slot >= 4u
-? floatBitsToUint(fF.w)
-: floatBitsToUint(fF.z);
-uint meta = (sE >> ((slot & 3u) * 8u)) & 0xffu;
+uint slot = uint(order) ^ vU;
+uint tT = slot >= 4u
+? floatBitsToUint(fU.w)
+: floatBitsToUint(fU.z);
+uint meta = (tT >> ((slot & 3u) * 8u)) & 0xffu;
 if (meta == 0u) continue;
-uint lY = (slot & 3u) * 8u;
-uint kV = slot >> 2u;
-vec3 tH = vec3(
-float((gI[kV] >> lY) & 0xffu),
-float((gI[2u + kV] >> lY) & 0xffu),
-float((gI[4u + kV] >> lY) & 0xffu)
+uint mI = (slot & 3u) * 8u;
+uint lF = slot >> 2u;
+vec3 uX = vec3(
+float((gV[lF] >> mI) & 0xffu),
+float((gV[2u + lF] >> mI) & 0xffu),
+float((gV[4u + lF] >> mI) & 0xffu)
 );
-vec3 tG = vec3(
-float((gI[6u + kV] >> lY) & 0xffu),
-float((gI[8u + kV] >> lY) & 0xffu),
-float((gI[10u + kV] >> lY) & 0xffu)
+vec3 uW = vec3(
+float((gV[6u + lF] >> mI) & 0xffu),
+float((gV[8u + lF] >> mI) & 0xffu),
+float((gV[10u + lF] >> mI) & 0xffu)
 );
-vec3 uW = jA.xyz + tH * hJ;
-vec3 uV = jA.xyz + tG * hJ;
-if (cM(o, d, uW, uV, tMax) <= 0.0) continue;
-uint jC = meta & 0x1fu;
-if (jC >= 24u) {
-uint ns = internalMask & ((1u << slot) - 1u);
-uint sX = childBase + qT(ns);
+vec3 wN = jZ.xyz + uX * hZ;
+vec3 wM = jZ.xyz + uW * hZ;
+if (cU(o, d, wN, wM, tMax) <= 0.0) continue;
+uint kb = meta & 0x1fu;
+if (kb >= 24u) {
+uint om = internalMask & ((1u << slot) - 1u);
+uint un = childBase + rY(om);
 stack[++top] = ivec2(
-eW + int(sX * 5u), cf
+fg + int(un * 5u), bY
 );
 }
 else {
-uint triangleCount = qT(meta >> 5u);
-uint rS = oQ + jC;
-uint vh = 0x80000000u |
+uint triangleCount = rY(meta >> 5u);
+uint sV = pM + kb;
+uint wZ = 0x80000000u |
 (triangleCount << 23u) |
-(rS & 0x007fffffu);
-stack[++top] = ivec2(int(vh), cf);
+(sV & 0x007fffffu);
+stack[++top] = ivec2(int(wZ), bY);
 }
 }
 }
 }
-else if (gR == 2) {
-int i = node.x + bo;
-vec4 uY = dp(i);
-vec4 uZ = dp(i + 1);
-vec4 va = dp(i + 2);
-vec4 vb = dp(i + 3);
-bh = int(floatBitsToUint(va.w));
-if (bh > 0) {
-primitiveOffset = bo + eW + int(floatBitsToUint(vb.w));
+else if (he == 2) {
+int i = node.x + bw;
+vec4 wP = dJ(i);
+vec4 wQ = dJ(i + 1);
+vec4 wR = dJ(i + 2);
+vec4 wS = dJ(i + 3);
+primitiveCount = int(floatBitsToUint(wR.w));
+if (primitiveCount > 0) {
+primitiveOffset = bw + fg + int(floatBitsToUint(wS.w));
 }
 else {
-int left = eW + int(floatBitsToUint(uY.w) * 4u);
-int right = eW + int(floatBitsToUint(uZ.w) * 4u);
-float d0 = cM(o, d, uY.xyz, uZ.xyz, tMax);
-float d1 = cM(o, d, va.xyz, vb.xyz, tMax);
+int left = fg + int(floatBitsToUint(wP.w) * 4u);
+int right = fg + int(floatBitsToUint(wQ.w) * 4u);
+float d0 = cU(o, d, wP.xyz, wQ.xyz, tMax);
+float d1 = cU(o, d, wR.xyz, wS.xyz, tMax);
 if (d0 > d1) {
-if (d0 > 0.0) stack[++top] = ivec2(left, cf);
-if (d1 > 0.0) stack[++top] = ivec2(right, cf);
+if (d0 > 0.0) stack[++top] = ivec2(left, bY);
+if (d1 > 0.0) stack[++top] = ivec2(right, bY);
 }
 else {
-if (d1 > 0.0) stack[++top] = ivec2(right, cf);
-if (d0 > 0.0) stack[++top] = ivec2(left, cf);
+if (d1 > 0.0) stack[++top] = ivec2(right, bY);
+if (d0 > 0.0) stack[++top] = ivec2(left, bY);
 }
 }
 }
-else if (gR == 1) {
-uint hW = uint(node.x);
-if ((hW & 0x80000000u) != 0u) {
-bh = int((hW >> 23u) & 0xffu);
-primitiveOffset = bo + eW + int(hW & 0x007fffffu);
+else if (he == 1) {
+uint iw = uint(node.x);
+if ((iw & 0x80000000u) != 0u) {
+primitiveCount = int((iw >> 23u) & 0xffu);
+primitiveOffset = bw + fg + int(iw & 0x007fffffu);
 }
 else {
-int i = node.x + bo;
-vec4 ui = dp(i);
-vec4 ro = dp(i + 1);
-vec4 rp = dp(i + 2);
-vec4 rq = dp(i + 3);
-vec3 nodeMin = ui.xyz;
-vec3 ec = vec3(ui.w, ro.x, ro.y);
-uint mh[4] = uint[4](
-floatBitsToUint(ro.z), floatBitsToUint(ro.w),
-floatBitsToUint(rp.x), floatBitsToUint(rp.y)
+int i = node.x + bw;
+vec4 vJ = dJ(i);
+vec4 sz = dJ(i + 1);
+vec4 qR = dJ(i + 2);
+vec4 oQ = dJ(i + 3);
+vec3 nodeMin = vJ.xyz;
+vec3 er = vec3(vJ.w, sz.x, sz.y);
+uint mT[4] = uint[4](
+floatBitsToUint(sz.z), floatBitsToUint(sz.w),
+floatBitsToUint(qR.x), floatBitsToUint(qR.y)
 );
-vec4 rw = fc(floatBitsToUint(rp.z));
-vec4 rA = fc(floatBitsToUint(rp.w));
-vec4 rB = fc(floatBitsToUint(rq.x));
-vec4 rt = fc(floatBitsToUint(rq.y));
-vec4 ru = fc(floatBitsToUint(rq.z));
-vec4 rv = fc(floatBitsToUint(rq.w));
-vec3 gC = fb(d);
-vec3 bJ = eZ(gC, o);
-vec4 tx0 = (nodeMin.x + rw * ec.x - o.x) * gC.x;
-vec4 tx1 = (nodeMin.x + rt * ec.x - o.x) * gC.x;
-vec4 ty0 = (nodeMin.y + rA * ec.y - o.y) * gC.y;
-vec4 ty1 = (nodeMin.y + ru * ec.y - o.y) * gC.y;
-vec4 tz0 = (nodeMin.z + rB * ec.z - o.z) * gC.z;
-vec4 tz1 = (nodeMin.z + rv * ec.z - o.z) * gC.z;
-vec4 mL = max(max((min(tx0, tx1) - vec4(bJ.x)), (min(ty0, ty1) - vec4(bJ.y))), (min(tz0, tz1) - vec4(bJ.z)));
-vec4 rU = min(min((max(tx0, tx1) + vec4(bJ.x)), (max(ty0, ty1) + vec4(bJ.y))), (max(tz0, tz1) + vec4(bJ.z)));
-vec4 distance = mix(rU, mL, greaterThan(mL, vec4(0.0)));
-distance = mix(distance, vec4(-1.0), greaterThan(mL, min(rU, vec4(tMax))));
-float fd[4] = float[4](distance.x, distance.y, distance.z, distance.w);
+vec4 wx = fp(floatBitsToUint(qR.z));
+vec4 wy = fp(floatBitsToUint(qR.w));
+vec4 wB = fp(floatBitsToUint(oQ.x));
+vec4 wu = fp(floatBitsToUint(oQ.y));
+vec4 wv = fp(floatBitsToUint(oQ.z));
+vec4 ww = fp(floatBitsToUint(oQ.w));
+vec3 gr = fo(d);
+vec3 padding = fk(gr, o);
+vec4 tx0 = (nodeMin.x + wx * er.x - o.x) * gr.x;
+vec4 tx1 = (nodeMin.x + wu * er.x - o.x) * gr.x;
+vec4 ty0 = (nodeMin.y + wy * er.y - o.y) * gr.y;
+vec4 ty1 = (nodeMin.y + wv * er.y - o.y) * gr.y;
+vec4 tz0 = (nodeMin.z + wB * er.z - o.z) * gr.z;
+vec4 tz1 = (nodeMin.z + ww * er.z - o.z) * gr.z;
+vec4 nJ = max(max((min(tx0, tx1) - vec4(padding.x)), (min(ty0, ty1) - vec4(padding.y))), (min(tz0, tz1) - vec4(padding.z)));
+vec4 sY = min(min((max(tx0, tx1) + vec4(padding.x)), (max(ty0, ty1) + vec4(padding.y))), (max(tz0, tz1) + vec4(padding.z)));
+vec4 distance = mix(sY, nJ, greaterThan(nJ, vec4(0.0)));
+distance = mix(distance, vec4(-1.0), greaterThan(nJ, min(sY, vec4(tMax))));
+float fq[4] = float[4](distance.x, distance.y, distance.z, distance.w);
 
 
 for (int order = 0; order < 4; order++) {
-float dX = -1.0;
-int el = -1;
+float em = -1.0;
+int eB = -1;
 for (int child = 0; child < 4; child++) {
-if (mh[child] != 0u && fd[child] > dX) {
-dX = fd[child];
-el = child;
+if (mT[child] != 0u && fq[child] > em) {
+em = fq[child];
+eB = child;
 }
 }
-if (el >= 0 && dX > 0.0) {
-uint nD = mh[el];
-int jI = int(nD);
-if ((nD & 0x80000000u) == 0u) {
-jI += eW;
+if (eB >= 0 && em > 0.0) {
+uint oy = mT[eB];
+int ki = int(oy);
+if ((oy & 0x80000000u) == 0u) {
+ki += fg;
 }
-stack[++top] = ivec2(jI, cf);
-fd[el] = -1.0;
+stack[++top] = ivec2(ki, bY);
+fq[eB] = -1.0;
 }
 }
 }
 }
 else {
-int i = node.x + bo;
-vec4 r1 = dp(i);
-int bb = floatBitsToInt(r1.w);
+int i = node.x + bw;
+vec4 r1 = dJ(i);
+int bl = floatBitsToInt(r1.w);
 
-if (bb >= 0) {
+if (bl >= 0) {
 int left = node.x + 2;
-int right = bb + eW;
+int right = bl + fg;
 
-float d0 = cM(o, d, dp(left + bo).xyz, dp(left + 1 + bo).xyz, tMax);
-float d1 = cM(o, d, dp(right + bo).xyz, dp(right + 1 + bo).xyz, tMax);
+float d0 = cU(o, d, dJ(left + bw).xyz, dJ(left + 1 + bw).xyz, tMax);
+float d1 = cU(o, d, dJ(right + bw).xyz, dJ(right + 1 + bw).xyz, tMax);
 if (d0 > d1) {
 if (d0 > 0.0) {
-stack[++top] = ivec2(left, cf);
+stack[++top] = ivec2(left, bY);
 }
 if (d1 > 0.0) {
-stack[++top] = ivec2(right, cf);
+stack[++top] = ivec2(right, bY);
 }
 }
 else {
 if (d1 > 0.0) {
-stack[++top] = ivec2(right, cf);
+stack[++top] = ivec2(right, bY);
 }
 if (d0 > 0.0) {
-stack[++top] = ivec2(left, cf);
+stack[++top] = ivec2(left, bY);
 }
 }
 }
 else {
-if (cf < 0) {
-for (int k = 0; k < -bb; k++) {
-vec4 r3 = dp(i + 2 + k);
+if (bY < 0) {
+for (int k = 0; k < -bl; k++) {
+vec4 r3 = dJ(i + 2 + k);
 stack[++top] = ivec2(floatBitsToInt(r3.y), floatBitsToInt(r3.x));
 }
 }
 else {
-bh = -bb;
+primitiveCount = -bl;
 primitiveOffset = i + 2;
 }
 }
 }
-if (bh > 0) {
+if (primitiveCount > 0) {
 ${intersectTriangles}
 }
 }
@@ -321,13 +321,13 @@ export default createShaderMixin({
     },
     main: glsl `
 
-struct fA {
+struct fP {
 vec3 p0;
 vec3 p1;
 vec3 p2;
 };
 #ifdef NUM_LIGHTS
-float hx(float rad, vec3 pos, Ray r) {
+float hL(float rad, vec3 pos, Ray r) {
 vec3 op = pos - r.o;
 float eps = 0.001;
 float b = dot(op, r.d);
@@ -347,7 +347,7 @@ return t2;
 return INF;
 }
 
-float ni(vec3 pos, vec3 u, vec3 v, vec3 n, Ray r) {
+float og(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;
 
@@ -366,7 +366,7 @@ return t;
 return INF;
 }
 
-float mz(Ray r, Light light, out LightSampleRec lightSample) {
+float nr(Ray r, Light light, out LightSampleRec lightSample) {
 
 
 
@@ -383,15 +383,15 @@ vec3 normal = normalize(cross(u, v));
 if (dot(normal, r.d) > 0.) 
 return INF;
 
-d = ni(pos, u, v, normal, r);
+d = og(pos, u, v, normal, r);
 if (d < 0.) {
 d = INF;
 }
 else {
-float cc = dot(-r.d, normal);
+float ci = dot(-r.d, normal);
 float spread = light.spread;
-if (cc >= cos(spread * PI / 2.0)) {
-lightSample.pdf = (d * d) / (light.area * cc);
+if (ci >= cos(spread * PI / 2.0)) {
+lightSample.pdf = (d * d) / (light.area * ci);
 lightSample.emission = light.emission;
 }
 }
@@ -400,10 +400,10 @@ lightSample.emission = light.emission;
 #ifdef HAS_SPHERE_LIGHTS
 if (type == SPHERE_LIGHT) {
 if (light.radius <= 0.0) {
-vec3 ko = light.position - r.o;
-float projection = dot(ko, r.d);
-vec3 ba = ko - projection * r.d;
-if (projection <= EPS || dot(ba, ba) > 1e-8) return INF;
+vec3 kT = light.position - r.o;
+float projection = dot(kT, r.d);
+vec3 bk = kT - projection * r.d;
+if (projection <= EPS || dot(bk, bk) > 1e-8) return INF;
 d = projection;
 lightSample.direction = r.d;
 lightSample.normal = -r.d;
@@ -411,19 +411,19 @@ lightSample.dist = d;
 lightSample.pdf = 1.0;
 lightSample.isDelta = true;
 } else {
-d = hx(light.radius, light.position, r);
+d = hL(light.radius, light.position, r);
 lightSample.isDelta = false;
 }
 if (d < 0.) {
 d = INF;
 }
 else {
-float aj = distance(r.o, light.position);
-float cosThetaMax = aj <= light.radius
+float ak = distance(r.o, light.position);
+float cosThetaMax = ak <= light.radius
 ? -1.0
-: sqrt(max(0.0, 1.0 - light.radius * light.radius / (aj * aj)));
+: sqrt(max(0.0, 1.0 - light.radius * light.radius / (ak * ak)));
 if (light.radius > 0.0) {
-lightSample.pdf = aj <= light.radius
+lightSample.pdf = ak <= light.radius
 ? INV4PI
 : 1.0 / max(TWOPI * (1.0 - cosThetaMax), 1e-8);
 }
@@ -436,10 +436,10 @@ lightSample.emission = light.radius <= 0.0
 #ifdef HAS_SPOT_LIGHTS
 if (type == SPOT_LIGHT) {
 if (light.radius <= 0.0) {
-vec3 ko = light.position - r.o;
-float projection = dot(ko, r.d);
-vec3 ba = ko - projection * r.d;
-if (projection <= EPS || dot(ba, ba) > 1e-8) return INF;
+vec3 kT = light.position - r.o;
+float projection = dot(kT, r.d);
+vec3 bk = kT - projection * r.d;
+if (projection <= EPS || dot(bk, bk) > 1e-8) return INF;
 d = projection;
 lightSample.direction = r.d;
 lightSample.normal = -r.d;
@@ -447,23 +447,23 @@ lightSample.dist = d;
 lightSample.pdf = 1.0;
 lightSample.isDelta = true;
 } else {
-d = hx(light.radius, light.position, r);
+d = hL(light.radius, light.position, r);
 lightSample.isDelta = false;
 }
 if (d != INF) {
-vec3 ve = r.o + d * r.d;
-float mC = bL(light, normalize(r.o - ve));
-if (mC <= 0.0) return INF;
-float aj = distance(r.o, light.position);
-float cosThetaMax = aj <= light.radius
+vec3 wW = r.o + d * r.d;
+float nv = bP(light, normalize(r.o - wW));
+if (nv <= 0.0) return INF;
+float ak = distance(r.o, light.position);
+float cosThetaMax = ak <= light.radius
 ? -1.0
-: sqrt(max(0.0, 1.0 - light.radius * light.radius / (aj * aj)));
+: sqrt(max(0.0, 1.0 - light.radius * light.radius / (ak * ak)));
 if (light.radius > 0.0) {
-lightSample.pdf = aj <= light.radius
+lightSample.pdf = ak <= light.radius
 ? INV4PI
 : 1.0 / max(TWOPI * (1.0 - cosThetaMax), 1e-8);
 }
-lightSample.emission = light.emission * mC * (light.radius <= 0.0
+lightSample.emission = light.emission * nv * (light.radius <= 0.0
 ? 1.0 / max(d * d, 1e-8)
 : 1.0);
 }
@@ -474,96 +474,96 @@ return d;
 }
 #endif
 
-int fm(ivec4 m0, ivec4 m1, ivec4 m2, vec3 eK) {
-vec2 m00 = ae(m0.x);
-vec2 m01 = ae(m0.y);
-vec2 m02 = ae(m0.z);
-vec2 m10 = ae(m1.x);
-vec2 m11 = ae(m1.y);
-vec2 m12 = ae(m1.z);
-vec2 m20 = ae(m2.x);
-vec2 m21 = ae(m2.y);
-vec2 m22 = ae(m2.z);
-float w00 = m00.y * eK.x;
-float w01 = m01.y * eK.x;
-float w02 = m02.y * eK.x;
-float w10 = m10.y * eK.y;
-float w11 = m11.y * eK.y;
-float w12 = m12.y * eK.y;
-float w20 = m20.y * eK.z;
-float w21 = m21.y * eK.z;
-float w22 = m22.y * eK.z;
-
-float 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;
-
-float rand = eo(0u);
-float nj;
+int fF(ivec4 m0, ivec4 m1, ivec4 m2, vec3 eU) {
+vec2 m00 = ah(m0.x);
+vec2 m01 = ah(m0.y);
+vec2 m02 = ah(m0.z);
+vec2 m10 = ah(m1.x);
+vec2 m11 = ah(m1.y);
+vec2 m12 = ah(m1.z);
+vec2 m20 = ah(m2.x);
+vec2 m21 = ah(m2.y);
+vec2 m22 = ah(m2.z);
+float w00 = m00.y * eU.x;
+float w01 = m01.y * eU.x;
+float w02 = m02.y * eU.x;
+float w10 = m10.y * eU.y;
+float w11 = m11.y * eU.y;
+float w12 = m12.y * eU.y;
+float w20 = m20.y * eU.z;
+float w21 = m21.y * eU.z;
+float w22 = m22.y * eU.z;
+
+float 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;
+
+float rand = eC(0u);
+float oh;
 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 int(nj);
+return int(oh);
 }
 
 
-vec4 bv(int i);
+vec4 bF(int i);
 
 
-mat2x3 eE(vec3 normal, vec3 dp1, vec3 dp2, vec2 od, vec2 oe, vec4 aT) {
-float determinant = od.x * oe.y - od.y * oe.x;
-vec3 t = aT.xyz;
+mat2x3 eP(vec3 normal, vec3 dp1, vec3 dp2, vec2 oW, vec2 oX, vec4 aX) {
+float determinant = oW.x * oX.y - oW.y * oX.x;
+vec3 t = aX.xyz;
 vec3 b = vec3(0.0);
 
-float qC = abs(aT.w) > 1.5 ? 0.0 : aT.w;
+float rF = abs(aX.w) > 1.5 ? 0.0 : aX.w;
 if (abs(determinant) > 1e-12) {
-vec3 dt = (dp1 * oe.y - dp2 * od.y) / determinant;
-vec3 db = (dp2 * od.x - dp1 * oe.x) / determinant;
+vec3 dt = (dp1 * oX.y - dp2 * oW.y) / determinant;
+vec3 db = (dp2 * oW.x - dp1 * oX.x) / determinant;
 if (dot(t, t) < 1e-8) { t = dt; }
-if (qC == 0.0) { qC = (dot(cross(normal, dt), db) < 0.0 ? -1.0 : 1.0); }
+if (rF == 0.0) { rF = (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 dW = -1.0 / (ik + normal.z);
-float xy = normal.x * normal.y * dW;
-t = vec3(1.0 + ik * normal.x * normal.x * dW, ik * xy, -ik * normal.x);
-b = vec3(xy, ik + normal.y * normal.y * dW, -normal.y);
+float iI = (normal.z >= 0.0 ? 1.0 : -1.0);
+float cL = -1.0 / (iI + normal.z);
+float xy = normal.x * normal.y * cL;
+t = vec3(1.0 + iI * normal.x * normal.x * cL, iI * xy, -iI * normal.x);
+b = vec3(xy, iI + normal.y * normal.y * cL, -normal.y);
 }
 else { t = normalize(t); b = normalize(cross(normal, t)); }
-return mat2x3(t, b * (qC < 0.0 ? -1.0 : 1.0));
+return mat2x3(t, b * (rF < 0.0 ? -1.0 : 1.0));
 }
 
-void ht(inout ad si, fA tri, vec3 barycentric, ivec3 index) {
+void hH(inout SurfaceInteraction si, fP tri, vec3 barycentric, ivec3 index) {
 si.hit = true;
 si.position = barycentric.x * tri.p0 + barycentric.y * tri.p1 + barycentric.z * tri.p2;
-ivec2 i0 = gQ(index.x, vertexColumns);
-ivec2 i1 = gQ(index.y, vertexColumns);
-ivec2 i2 = gQ(index.z, vertexColumns);
+ivec2 i0 = hd(index.x, vertexColumns);
+ivec2 i1 = hd(index.y, vertexColumns);
+ivec2 i2 = hd(index.z, vertexColumns);
 si.vertexColor = vec3(1.0);
 #ifdef USE_VERTEX_COLOR
 si.vertexColor = barycentric.x * texelFetch(colorBuffer, i0, 0).rgb
@@ -582,14 +582,14 @@ si.vertexNormal = normal;
 
 
 vec3 geometricNormal = normalize(cross(tri.p1 - tri.p0, tri.p2 - tri.p0));
-vec3 km = barycentric.x * n0 + barycentric.y * n1 + barycentric.z * n2;
+vec3 kR = 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(km, geometricNormal) > 0.0) {
+if (dot(kR, 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);
@@ -607,7 +607,7 @@ 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;
 
 
 ivec4 m0 = ivec4(texelFetch(materialIndexBuffer, i0, 0));
@@ -615,7 +615,7 @@ ivec4 m1 = ivec4(texelFetch(materialIndexBuffer, i1, 0));
 ivec4 m2 = ivec4(texelFetch(materialIndexBuffer, i2, 0));
 
 if (m0.x > 0 || m1.x > 0 || m2.x > 0) {
-si.materialIndex = fm(m0, m1, m2, barycentric);
+si.materialIndex = fF(m0, m1, m2, barycentric);
 }
 
 vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
@@ -626,78 +626,80 @@ 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 cX = barycentric.x * secondaryUv0
+vec2 dj = barycentric.x * secondaryUv0
 + barycentric.y * secondaryUv1
 + barycentric.z * secondaryUv2;
 
-vec4 lM = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_UV);
-int uvProjection = int(lM.r);
+vec4 mw = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_UV);
+int uvProjection = int(mw.r);
 if (uvProjection > 0) {
-uv = gH(uvProjection, si.position, normal, lM.zw);
-cX = uv;
+uv = gU(uvProjection, si.position, normal, mw.zw);
+dj = uv;
 } else {
-uv *= lM.zw;
-cX *= lM.zw;
+uv *= mw.zw;
+dj *= mw.zw;
 }
 #else
 vec2 uv = vec2(0.0);
-vec2 cX = vec2(0.0);
+vec2 dj = vec2(0.0);
 #endif
 
-#ifdef NUM_COAT_NORMAL_MAPS
+#ifdef USE_COAT_NORMAL_FRAMES
 si.coatNormal = normal;
-float lw = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_LAYER_MAPS).z;
+#ifdef NUM_COAT_NORMAL_MAPS
+float pQ = 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 tm = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_FILM_TEXTURE).x;
-if (directionMaterial.x >= 0.0 || directionMaterial.y != 0.0 || tm != 0.0
-#ifdef NUM_COAT_NORMAL_MAPS
-|| lw >= 0.0
+float uL = fetchMatData(si.materialIndex * MAT_DATA_SIZE + MAT_ROW_FILM_TEXTURE).x;
+if (directionMaterial.x >= 0.0 || directionMaterial.y != 0.0 || uL != 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 aT = vec4(0.0);
+vec4 aX = vec4(0.0);
 #ifdef USE_TANGENT
-aT = barycentric.x * texelFetch(tangentBuffer, i0, 0)
+aX = barycentric.x * texelFetch(tangentBuffer, i0, 0)
 + barycentric.y * texelFetch(tangentBuffer, i1, 0)
 + barycentric.z * texelFetch(tangentBuffer, i2, 0);
 #endif
-mat2x3 eS = eE(normal, dp1, dp2, uv0-uv2, uv1-uv2, aT);
-si.geometryTangent = eS[0];
-si.geometryBitangent = eS[1];
-mat2x3 fY = eS;
+mat2x3 fc = eP(normal, dp1, dp2, uv0-uv2, uv1-uv2, aX);
+si.geometryTangent = fc[0];
+si.geometryBitangent = fc[1];
+mat2x3 gm = fc;
 #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;
-fY = eE(normal, dp1, dp2, a-c, b-c, (aT.w == 0.0 ? vec4(0.0) : aT));
+gm = eP(normal, dp1, dp2, a-c, b-c, (aX.w == 0.0 ? vec4(0.0) : aX));
 }
 #endif
-si.directionTangent = fY[0];
-si.directionBitangent = fY[1];
-#ifdef NUM_COAT_NORMAL_MAPS
-if (lw >= 0.0) {
-int kC = int(lw);
-mat2x3 mi = eS;
-if (mapUvChannel[kC] > 0.5) {
+si.directionTangent = gm[0];
+si.directionBitangent = gm[1];
+#if defined(USE_COAT_NORMAL_FRAMES) && defined(NUM_COAT_NORMAL_MAPS)
+if (pQ >= 0.0) {
+int lg = int(pQ);
+mat2x3 mU = fc;
+if (mapUvChannel[lg] > 0.5) {
 vec2 a = texelFetch(uv1Buffer, i0, 0).xy;
 vec2 b = texelFetch(uv1Buffer, i1, 0).xy;
 vec2 c = texelFetch(uv1Buffer, i2, 0).xy;
-mi = eE(normal, dp1, dp2, a-c, b-c,
-aT.w == 0.0 ? vec4(0.0) : aT);
+mU = eP(normal, dp1, dp2, a-c, b-c,
+aX.w == 0.0 ? vec4(0.0) : aX);
 }
-vec3 decoded = sampleMap(kC, getMapUvSet(kC, uv, cX)).rgb * 2.0 - 1.0;
-vec3 mapped = mi[0] * decoded.x * directionMaterial.z
-+ mi[1] * decoded.y * directionMaterial.w + normal * decoded.z;
+vec3 decoded = sampleMap(lg, getMapUvSet(lg, uv, dj)).rgb * 2.0 - 1.0;
+vec3 mapped = mU[0] * decoded.x * directionMaterial.z
++ mU[1] * decoded.y * directionMaterial.w + normal * decoded.z;
 if (dot(mapped, mapped) > 1e-12) { si.coatNormal = normalize(mapped); }
 }
 #endif
 }
 
 si.uv = uv;
-si.uv1 = cX;
+si.uv1 = dj;
 
 #ifdef USE_TANGENT
 vec3 t0 = texelFetch(tangentBuffer, i0, 0).xyz;
@@ -708,52 +710,52 @@ si.tangent = normalize(barycentric.x * t0 + barycentric.y * t1 + barycentric.z *
 #if defined(NUM_NORMAL_MAPS)
 vec3 dp1 = tri.p0 - tri.p2;
 vec3 dp2 = tri.p1 - tri.p2;
-vec2 od = uv0 - uv2;
-vec2 oe = uv1 - uv2;
-si.normal = om(
-si.materialIndex, uv, normal, dp1, dp2, od, oe
+vec2 oW = uv0 - uv2;
+vec2 oX = uv1 - uv2;
+si.normal = pd(
+si.materialIndex, uv, normal, dp1, dp2, oW, oX
 );
 #else
 si.normal = normal;
 #endif
 }
 
-struct eC {
+struct eK {
 float t;
 vec3 barycentric;
 };
 
-vec4 dp(int i) {
-return kG(bvhBuffer, i, bvhColumns);
+vec4 dJ(int i) {
+return lm(bvhBuffer, i, bvhColumns);
 }
 
-vec4 bv(int i) {
-return kG(instancesBuffer, i, instancesColumns);
+vec4 bF(int i) {
+return lm(instancesBuffer, i, instancesColumns);
 }
 
 
 
-uint qT(uint value) {
+uint rY(uint value) {
 value -= (value >> 1u) & 0x55555555u;
 value = (value & 0x33333333u) + ((value >> 2u) & 0x33333333u);
 value = (value + (value >> 4u)) & 0x0f0f0f0fu;
 return (value * 0x01010101u) >> 24u;
 }
 
-vec4 fc(uint lV) {
+vec4 fp(uint mE) {
 return vec4(
-float(lV & 0xffu),
-float((lV >> 8u) & 0xffu),
-float((lV >> 16u) & 0xffu),
-float(lV >> 24u)
+float(mE & 0xffu),
+float((mE >> 8u) & 0xffu),
+float((mE >> 16u) & 0xffu),
+float(mE >> 24u)
 ) / 255.0;
 }
 
 
 
 
-eC fL(vec3 o, vec3 d, fA tri, int ey, vec3 gP, float tMax) {
-eC ti;
+eK ga(vec3 o, vec3 d, fP tri, int eq, vec3 gY, float tMax) {
+eK ti;
 
 
 vec3 p0t = tri.p0 - o;
@@ -761,20 +763,20 @@ vec3 p1t = tri.p1 - o;
 vec3 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;
 }
 
 
-p0t.xy += gP.xy * p0t.z;
-p1t.xy += gP.xy * p1t.z;
-p2t.xy += gP.xy * p2t.z;
+p0t.xy += gY.xy * p0t.z;
+p1t.xy += gY.xy * p1t.z;
+p2t.xy += gY.xy * p2t.z;
 
 
 vec3 e = vec3(
@@ -784,9 +786,9 @@ p0t.x * p1t.y - p0t.y * p1t.x
 );
 
 #ifdef DEPTH_QUERY_MODE
-vec3 rE = abs(vec3(p1t.x * p2t.y, p2t.x * p0t.y, p0t.x * p1t.y))
+vec3 sF = 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), rE * 1e-6));
+e = mix(e, vec3(0.0), lessThanEqual(abs(e), sF * 1e-6));
 #endif
 
 
@@ -799,19 +801,19 @@ float det = e.x + e.y + e.z;
 if (det == 0.0) return ti;
 #endif
 
-p0t.z *= gP.z;
-p1t.z *= gP.z;
-p2t.z *= gP.z;
-float 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;
+float 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;
 }
 
-float sw = 1. / det;
-ti.t = rn * sw;
-ti.barycentric = e * sw;
+float tN = 1. / det;
+ti.t = sw * tN;
+ti.barycentric = e * tN;
 
 return ti;
 }
@@ -822,18 +824,18 @@ vec3 max;
 };
 
 
-vec3 fb(vec3 d) {
+vec3 fo(vec3 d) {
 #ifdef DEPTH_QUERY_MODE
-vec3 ap = mix(d, mix(vec3(1e-20), vec3(-1e-20), lessThan(d, vec3(0.0))), lessThan(abs(d), vec3(1e-20)));
-return 1.0 / ap;
+vec3 aq = mix(d, mix(vec3(1e-20), vec3(-1e-20), lessThan(d, vec3(0.0))), lessThan(abs(d), vec3(1e-20)));
+return 1.0 / aq;
 #else
 return 1.0 / d;
 #endif
 }
 
-vec3 eZ(vec3 gC, vec3 o) {
+vec3 fk(vec3 gr, vec3 o) {
 #ifdef DEPTH_QUERY_MODE
-return abs(gC) * max(vec3(1.0), abs(o)) * 1e-6;
+return abs(gr) * max(vec3(1.0), abs(o)) * 1e-6;
 #else
 return vec3(0.0);
 #endif
@@ -841,21 +843,21 @@ return vec3(0.0);
 
 
 
-float cM(vec3 o, vec3 d, vec3 vK, vec3 vJ, float tMax) {
+float cU(vec3 o, vec3 d, vec3 xq, vec3 xp, float tMax) {
 #ifdef DEPTH_QUERY_MODE
-vec3 gC = fb(d);
-vec3 uR = (vK - o) * gC;
-vec3 uS = (vJ - o) * gC;
-vec3 bJ = eZ(gC, o);
+vec3 gr = fo(d);
+vec3 wG = (xq - o) * gr;
+vec3 wI = (xp - o) * gr;
+vec3 padding = fk(gr, o);
 #else
-vec3 uR = (vK - o) / d;
-vec3 uS = (vJ - o) / d;
-vec3 bJ = vec3(0.0);
+vec3 wG = (xq - o) / d;
+vec3 wI = (xp - o) / d;
+vec3 padding = vec3(0.0);
 #endif
-vec3 mL = min(uR, uS) - bJ;
-vec3 rU = max(uR, uS) + bJ;
-float t0 = max(mL.x, max(mL.y, mL.z));
-float t1 = min(rU.x, min(rU.y, rU.z));
+vec3 nJ = min(wG, wI) - padding;
+vec3 sY = max(wG, wI) + padding;
+float t0 = max(nJ.x, max(nJ.y, nJ.z));
+float t1 = min(sY.x, min(sY.y, sY.z));
 
 return (t0 > t1 || t0 > tMax) ? -1.0 : (t0 > 0.0 ? t0 : t1);
 }
@@ -864,12 +866,12 @@ int maxDimension(vec3 v) {
 return v.x > v.y ? (v.x > v.z ? 0 : 2) : (v.y > v.z ? 1 : 2);
 }
 
-vec3 mv(vec3 d, int ey) {
+vec3 nk(vec3 d, int eq) {
 
 
-if (ey == 0) {
+if (eq == 0) {
 return vec3(-d.y, -d.z, 1.0) / d.x;
-} else if (ey == 1) {
+} else if (eq == 1) {
 return vec3(-d.z, -d.x, 1.0) / d.y;
 } else {
 return vec3(-d.x, -d.y, 1.0) / d.z;
@@ -877,50 +879,50 @@ return vec3(-d.x, -d.y, 1.0) / d.z;
 }
 
 
-void cZ(in Ray ray, inout ad si, bool doubleSided) {
+void intersectScene(in Ray ray, inout SurfaceInteraction si, bool doubleSided) {
 si.hit = false;
 
 if (isEmptyScene) {
 return;
 }
 
-int ey = maxDimension(abs(ray.d));
-vec3 gP = mv(ray.d, ey);
-int ao;
+int eq = maxDimension(abs(ray.d));
+vec3 gY = nk(ray.d, eq);
+int ap;
 
 ${getBVHTraverseCode(glsl `
 si.triangleIndex = triangleIndex;
-lC = tri;
-hB = faceNormal;
-aB = hit.barycentric;
-hw = index;
-rb = true;
-ao = cf;
+mk = tri;
+hS = faceNormal;
+aC = hit.barycentric;
+hK = index;
+sj = true;
+ap = bY;
 `)}
 
-if (rb) {
-si.instanceIndex = ao;
-si.materialIndex = int(ae(floatBitsToInt(bv(ao * ${INSTANCE_VEC4_STRIDE}).z)).x);
+if (sj) {
+si.instanceIndex = ap;
+si.materialIndex = int(ah(floatBitsToInt(bF(ap * ${INSTANCE_VEC4_STRIDE}).z)).x);
 
-ht(si, lC, aB, hw);
+hH(si, mk, aC, hK);
 mat4 l2w = mat4(
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 5),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 6),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 7),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 8)
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 5),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 6),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 7),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 8)
 );
 mat4 w2l = mat4(
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 1),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 2),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 3),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 4)
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 1),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 2),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 3),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 4)
 );
-mat3 aD = transpose(mat3(w2l));
-si.faceNormal = normalize(aD * hB);
-si.vertexNormal = normalize(aD * si.vertexNormal);
-si.normal = normalize(aD * si.normal);
-#ifdef NUM_COAT_NORMAL_MAPS
-si.coatNormal = normalize(aD * si.coatNormal);
+mat3 aH = transpose(mat3(w2l));
+si.faceNormal = normalize(aH * hS);
+si.vertexNormal = normalize(aH * si.vertexNormal);
+si.normal = normalize(aH * si.normal);
+#ifdef USE_COAT_NORMAL_FRAMES
+si.coatNormal = normalize(aH * si.coatNormal);
 #endif
 #ifdef USE_TANGENT
 if (dot(si.tangent, si.tangent) > 0.0) {
@@ -937,7 +939,7 @@ si.position = (l2w * vec4(si.position, 1.0)).xyz;
 }
 
 
-void iM(in Ray ray, inout ad si, int instanceIndex, float maxDistance) {
+void intersectSceneBevel(in Ray ray, inout SurfaceInteraction si, int instanceIndex, float maxDistance) {
 si.hit = false;
 
 if (isEmptyScene) {
@@ -945,46 +947,46 @@ return;
 }
 
 bool doubleSided = true;
-int ey = maxDimension(abs(ray.d));
-vec3 gP = mv(ray.d, ey);
-int ao;
+int eq = maxDimension(abs(ray.d));
+vec3 gY = nk(ray.d, eq);
+int ap;
 
 ${getBVHTraverseCode(glsl `
 si.triangleIndex = triangleIndex;
-lC = tri;
-hB = faceNormal;
-aB = hit.barycentric;
-hw = index;
-rb = true;
-ao = cf;
+mk = tri;
+hS = faceNormal;
+aC = hit.barycentric;
+hK = index;
+sj = true;
+ap = bY;
 `, glsl `
-if (cf != instanceIndex) {
+if (bY != instanceIndex) {
 continue;
 }
 `, 'maxDistance')}
 
-if (rb) {
-si.instanceIndex = ao;
-si.materialIndex = int(ae(floatBitsToInt(bv(ao * ${INSTANCE_VEC4_STRIDE}).z)).x);
-ht(si, lC, aB, hw);
+if (sj) {
+si.instanceIndex = ap;
+si.materialIndex = int(ah(floatBitsToInt(bF(ap * ${INSTANCE_VEC4_STRIDE}).z)).x);
+hH(si, mk, aC, hK);
 mat4 l2w = mat4(
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 5),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 6),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 7),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 8)
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 5),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 6),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 7),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 8)
 );
 mat4 w2l = mat4(
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 1),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 2),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 3),
-bv(ao * ${INSTANCE_VEC4_STRIDE} + 4)
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 1),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 2),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 3),
+bF(ap * ${INSTANCE_VEC4_STRIDE} + 4)
 );
-mat3 aD = transpose(mat3(w2l));
-si.faceNormal = normalize(aD * hB);
-si.vertexNormal = normalize(aD * si.vertexNormal);
-si.normal = normalize(aD * si.normal);
-#ifdef NUM_COAT_NORMAL_MAPS
-si.coatNormal = normalize(aD * si.coatNormal);
+mat3 aH = transpose(mat3(w2l));
+si.faceNormal = normalize(aH * hS);
+si.vertexNormal = normalize(aH * si.vertexNormal);
+si.normal = normalize(aH * si.normal);
+#ifdef USE_COAT_NORMAL_FRAMES
+si.coatNormal = normalize(aH * si.coatNormal);
 #endif
 #ifdef USE_TANGENT
 if (dot(si.tangent, si.tangent) > 0.0) {
@@ -1000,35 +1002,35 @@ si.position = (l2w * vec4(si.position, 1.0)).xyz;
 }
 }
 
-bool ex(in Ray ray, float lightDistance, int jL) {
+bool intersectSceneShadow(in Ray ray, float lightDistance, int kn) {
 
 if (isEmptyScene) {
 return false;
 }
 
 bool doubleSided = true;
-int ey = maxDimension(abs(ray.d));
-vec3 gP = mv(ray.d, ey);
+int eq = maxDimension(abs(ray.d));
+vec3 gY = nk(ray.d, eq);
 
 ${getBVHTraverseCode(glsl `
 
 
-bool jt = true;
+bool jS = true;
 vec3 barycentric = hit.barycentric;
 vec3 position = barycentric.x * tri.p0 + barycentric.y * tri.p1 + barycentric.z * tri.p2;
 
-if (jt) {
+if (jS) {
 #ifdef HAS_ALPHA
-int materialIndex = int(ae(floatBitsToInt(bv(cf * ${INSTANCE_VEC4_STRIDE}).z)).x);
-ivec2 i0 = gQ(index.x, vertexColumns);
-ivec2 i1 = gQ(index.y, vertexColumns);
-ivec2 i2 = gQ(index.z, vertexColumns);
+int materialIndex = int(ah(floatBitsToInt(bF(bY * ${INSTANCE_VEC4_STRIDE}).z)).x);
+ivec2 i0 = hd(index.x, vertexColumns);
+ivec2 i1 = hd(index.y, vertexColumns);
+ivec2 i2 = hd(index.z, vertexColumns);
 ivec4 m0 = ivec4(texelFetch(materialIndexBuffer, i0, 0));
 ivec4 m1 = ivec4(texelFetch(materialIndexBuffer, i1, 0));
 ivec4 m2 = ivec4(texelFetch(materialIndexBuffer, i2, 0));
 
 if (m0.x > 0 || m1.x > 0 || m2.x > 0) {
-materialIndex = fm(m0, m1, m2, barycentric);
+materialIndex = fF(m0, m1, m2, barycentric);
 }
 
 int ii = materialIndex * MAT_DATA_SIZE;
@@ -1042,12 +1044,12 @@ vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
 vec2 uv1 = texelFetch(uvBuffer, i1, 0).xy;
 vec2 uv2 = texelFetch(uvBuffer, i2, 0).xy;
 vec2 uv = barycentric.x * uv0 + barycentric.y * uv1 + barycentric.z * uv2;
-vec4 lM = fetchMatData(ii + MAT_ROW_UV);
-int uvProjection = int(lM.r);
+vec4 mw = fetchMatData(ii + MAT_ROW_UV);
+int uvProjection = int(mw.r);
 if (uvProjection > 0) {
-uv = gH(uvProjection, position, hB, lM.zw);
+uv = gU(uvProjection, position, hS, mw.zw);
 } else {
-uv *= lM.zw;
+uv *= mw.zw;
 }
 
 alpha *= sampleMap(diffuseMapIndex, uv).a;
@@ -1061,32 +1063,32 @@ vec2 uv0 = texelFetch(uvBuffer, i0, 0).xy;
 vec2 uv1 = texelFetch(uvBuffer, i1, 0).xy;
 vec2 uv2 = texelFetch(uvBuffer, i2, 0).xy;
 vec2 uv = barycentric.x * uv0 + barycentric.y * uv1 + barycentric.z * uv2;
-vec4 lM = fetchMatData(ii + MAT_ROW_UV);
-int uvProjection = int(lM.r);
+vec4 mw = fetchMatData(ii + MAT_ROW_UV);
+int uvProjection = int(mw.r);
 if (uvProjection > 0) {
-uv = gH(uvProjection, position, hB, lM.zw);
+uv = gU(uvProjection, position, hS, mw.zw);
 } else {
-uv *= lM.zw;
+uv *= mw.zw;
 }
 alpha *= sampleMap(alphaMapIndex, uv).g;
 }
 #endif
 float alphaCutoff = fetchMatData(ii + MAT_ROW_STANDALONE_MAPS).a;
-alpha = lr(
+alpha = mb(
 alpha, int(fetchMatData(ii + MAT_ROW_ALPHA).w), alphaCutoff
 );
-if (eo(0u) > alpha) {
+if (eC(0u) > alpha) {
 
 tMax = min(max(lightDistance, 0.0), RAY_MAX_DISTANCE);
-jt = false;
+jS = false;
 }
 #endif
 }
-if (jt) {
+if (jS) {
 return true;
 }
 `, glsl `
-if (triangleIndex == jL) {
+if (triangleIndex == kn) {
 continue;
 }
 `, 'min(max(lightDistance, 0.0), RAY_MAX_DISTANCE)')}