@meimu/mutracer 0.1.5 → 0.3.0

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

@@ -1,56 +1,56 @@
 const displacementCommon = /*wgsl*/ `
 
 
-struct jD {
+struct jI {
 barycentric: vec3f,
 t: f32,
 normal: vec3f,
 
 }
 
-fn uX(
+fn xp(
 uu: f32,
-tX: vec3f, tY: vec3f,
+vK: vec3f, vL: vec3f,
 eAC: vec3f, eBD: vec3f,
-iO: vec3f,
-pX: f32,
-pZ: ptr<function, f32>
+iM: vec3f,
+qK: f32,
+qM: ptr<function, f32>
 ) -> bool {
 if (uu >= 0.0 && uu <= 1.0) {
-let pAB = mix(tX, tY, uu);
+let pAB = mix(vK, vL, uu);
 let pCD = mix(eAC, eBD, uu);
-var n = cross(iO, pCD);
-let vA = 1.0 / dot(n, n);
+var n = cross(iM, pCD);
+let xH = 1.0 / dot(n, n);
 n = cross(n, pAB);
-let tt = dot(n, pCD) * vA;
-let vv = dot(n, iO) * vA;
-if (vv >= 0.0 && vv <= 1.0 && tt > pX && tt < *pZ) {
-*pZ = tt;
+let tt = dot(n, pCD) * xH;
+let vv = dot(n, iM) * xH;
+if (vv >= 0.0 && vv <= 1.0 && tt > qK && tt < *qM) {
+*qM = tt;
 return true;
 }
 }
 return false;
 }
 
-fn hj(
-pP: vec3f, iO: vec3f, pX: f32, tn: f32,
+fn hc(
+qx: vec3f, iM: vec3f, qK: f32, uR: f32,
 pA: vec3f, pB: vec3f, pC: vec3f, pD: vec3f,
-pZ: ptr<function, f32>
+qM: ptr<function, f32>
 ) -> bool {
 let eAB = pB - pA;
 let eAC = pC - pA;
 let eBD = pD - pB;
 let eCD = pD - pC;
-let tX = pA - pP;
-let tY = pB - pP;
+let vK = pA - qx;
+let vL = pB - qx;
 
 var u1: f32;
 var u2: f32;
 {
 let qN = cross(eAB, -eCD);
-let a = dot(qN, iO);
-var b = dot(cross(tY, iO), eBD);
-let c = dot(cross(tX, iO), eAC);
+let a = dot(qN, iM);
+var b = dot(cross(vL, iM), eBD);
+let c = dot(cross(vK, iM), eAC);
 b -= a + c;
 
 let det = b * b - 4.0 * a * c;
@@ -62,19 +62,19 @@ if (a == 0.0) {
 u1 = -c / b;
 u2 = -1.0;
 } else {
-let vQ = sqrt(det);
-let wL = -0.5 * (b + sign(b) * vQ);
-u1 = wL / a;
-u2 = c / wL;
+let yg = sqrt(det);
+let yN = -0.5 * (b + sign(b) * yg);
+u1 = yN / a;
+u2 = c / yN;
 }
 }
 
-*pZ = tn;
+*qM = uR;
 
 var hit = false;
 
-hit = uX(u1, tX, tY, eAC, eBD, iO, pX, pZ) || 
-uX(u2, tX, tY, eAC, eBD, iO, pX, pZ);
+hit = xp(u1, vK, vL, eAC, eBD, iM, qK, qM) || 
+xp(u2, vK, vL, eAC, eBD, iM, qK, qM);
 if (!hit) {
 return false;
 }
@@ -82,70 +82,70 @@ return false;
 return true;
 }
 
-fn sv(
-pP: vec3f, iO: vec3f, pX: f32, tn: f32,
+fn tO(
+qx: vec3f, iM: vec3f, qK: f32, uR: f32,
 pA: vec3f, pB: vec3f, pC: vec3f,
 pD: vec3f, pE: vec3f, pF: vec3f,
 
-bT: ptr<function, f32>,
-bS: ptr<function, f32>
+bO: ptr<function, f32>,
+bN: ptr<function, f32>
 ) -> bool {
-*bT = INF;
-*bS = -INF;
+*bO = INF;
+*bN = -INF;
 
-let wd = maxDimension(abs(iO));
-let gQ = mu(iO, wd);
+let yv = maxDimension(abs(iM));
+let gH = mU(iM, yv);
 var tt: f32;
 
 tt = fM(
-pP, iO,
-fB(pC, pB, pA), wd, gQ, INF
+qx, iM,
+fw(pC, pB, pA), yv, gH, INF
 ).t;
 if (tt > 0.0) {
-*bT = min(*bT, tt);
-*bS = max(*bS, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
 tt = fM(
-pP, iO,
-fB(pD, pE, pF), wd, gQ, INF
+qx, iM,
+fw(pD, pE, pF), yv, gH, INF
 ).t;
 if (tt > 0.0) {
-*bT = min(*bT, tt);
-*bS = max(*bS, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hj(
-pP, iO, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pA, pB, pD, pE,
 &tt)) {
-*bT = min(*bT, tt);
-*bS = max(*bS, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hj(
-pP, iO, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pB, pC, pE, pF,
 &tt)) {
-*bT = min(*bT, tt);
-*bS = max(*bS, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hj(
-pP, iO, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pC, pA, pF, pD,
 &tt)) {
-*bT = min(*bT, tt);
-*bS = max(*bS, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
 
-*bT = max(*bT, pX);
-*bS = min(*bS, tn);
+*bO = max(*bO, qK);
+*bN = min(*bN, uR);
 
-return *bT <= *bS && *bS > 0.0;
+return *bO <= *bN && *bN > 0.0;
 }
 `;
 export default /* wgsl */ `
 
 ${displacementCommon}
 
-fn lP(p: vec3f, p0: vec3f, p1: vec3f, p2: vec3f) -> vec3f {
+fn mg(p: vec3f, p0: vec3f, p1: vec3f, p2: vec3f) -> vec3f {
 let v0 = p1 - p0;
 let v1 = p2 - p0;
 let v2 = p - p0;
@@ -156,9 +156,9 @@ let d11 = dot(v1, v1);
 let d20 = dot(v2, v0);
 let d21 = dot(v2, v1);
 
-let tT = 1.0 / (d00 * d11 - d01 * d01);
-let v = clamp((d11 * d20 - d01 * d21) * tT, 0.0, 1.0);
-let w = clamp((d00 * d21 - d01 * d20) * tT, 0.0, 1.0);
+let vF = 1.0 / (d00 * d11 - d01 * d01);
+let v = clamp((d11 * d20 - d01 * d21) * vF, 0.0, 1.0);
+let w = clamp((d00 * d21 - d01 * d20) * vF, 0.0, 1.0);
 let u = 1.0 - v - w;
 
 return vec3f(u, v, w);
@@ -166,59 +166,59 @@ return vec3f(u, v, w);
 
 
 
-fn pc(
-qr: vec3f, qM: vec3f, wQ: f32,
+fn pS(
+rl: vec3f, rP: vec3f, yR: f32,
 _pA: vec3f, _pB: vec3f, _pC: vec3f,
 _nA: vec3f, _nB: vec3f, _nC: vec3f,
-vG: vec2f, vH: vec2f, vI: vec2f,
+xQ: vec2f, xR: vec2f, xS: vec2f,
 
-jQ: i32,
-vF: vec2f,
+jW: i32,
+xP: vec2f,
 
-wk: f32, wj: f32,
+yE: f32, yD: f32,
 hScale: f32,
-uJ: f32,
-uf: f32,
+wG: f32,
+vR: f32,
 uvProjection: i32,
 uvScale: vec2f,
-vV: bool,
-rc: bool
-) -> jD {
+ym: bool,
+se: bool
+) -> jI {
 
 
 
-let DT = select(max(abs(hScale) / 32.0, 1e-12), uf, uf > 0.0);
-let dt = DT / length(qM);
+let DT = select(max(abs(hScale) / 32.0, 1e-12), vR, vR > 0.0);
+let dt = DT / length(rP);
 let dbx = 0.001;
 let dby = 0.001;
 
-var result: jD;
+var result: jI;
 result.t = INF;
 
-var bT = INF;
-var bS = -INF;
+var bO = INF;
+var bN = -INF;
 
-let ra = sign(hScale);
-let uP = abs(hScale);
+let sb = sign(hScale);
+let wZ = abs(hScale);
 
-let te = ra * normalize(cross(_pB - _pA, _pC - _pA));
+let uG = sb * normalize(cross(_pB - _pA, _pC - _pA));
 
 
-let dA = dot(_nA, te);
-let dB = dot(_nB, te);
-let dC = dot(_nC, te);
+let dA = dot(_nA, uG);
+let dB = dot(_nB, uG);
+let dC = dot(_nC, uG);
 
-var minH = wk * hScale * min(min(dA, dB), dC);
-var maxH = wj * hScale * max(max(dA, dB), dC);
+var minH = yE * hScale * min(min(dA, dB), dC);
+var maxH = yD * hScale * max(max(dA, dB), dC);
 
 
 
 
-let tZ = abs(maxH) - abs(minH);
+let vM = abs(maxH) - abs(minH);
 
 
 minH = sign(minH) * max(abs(minH) - DT * 2.0, 0.0);
-if (tZ < DT * 2.0) {
+if (vM < DT * 2.0) {
 maxH = sign(maxH) * min(abs(maxH) + DT * 2.0, abs(hScale));
 }
 
@@ -230,70 +230,70 @@ let pD = _pA + _nA * maxH / dA;
 let pE = _pB + _nB * maxH / dB;
 let pF = _pC + _nC * maxH / dC;
 
-if (!sv(
-qr, qM, 0.0, INF,
+if (!tO(
+rl, rP, 0.0, INF,
 pA, pB, pC,
 pD, pE, pF,
 
-&bT, &bS
+&bO, &bN
 )) {
 return result;
 }
 
-var wM = false;
+var yO = false;
 
-if (bS == bT) {
-bT = 0.0;
+if (bN == bO) {
+bO = 0.0;
 }
 
-var t = bT + dt * eo(2u);
-var s = qr + qM * t;
-var tS = dot(s - _pA, te);
-var vd = select(0.0, INF, ra < 0.0);
+var t = bO + dt * ek(2u);
+var s = rl + rP * t;
+var vE = dot(s - _pA, uG);
+var hSurf = select(0.0, INF, sb < 0.0);
 
-var c0 = _pA + _nA * tS / dA;
-var c1 = _pB + _nB * tS / dB;
-var c2 = _pC + _nC * tS / dC;
+var c0 = _pA + _nA * vE / dA;
+var c1 = _pB + _nB * vE / dB;
+var c2 = _pC + _nC * vE / dC;
 
-var tF: vec3f;
+var vc: vec3f;
 
-tS = abs(tS);
+vE = abs(vE);
 
-var qZ = tS - vd;
-var pY = qZ;
+var hDiff = vE - hSurf;
+var qL = hDiff;
 var step = 0;
 
 while (true) {
 
-let b = lP(s, c0, c1, c2);
+let b = mg(s, c0, c1, c2);
 
-tF = _pA * b.x + _pB * b.y + _pC * b.z;
+vc = _pA * b.x + _pB * b.y + _pC * b.z;
 
-tS = distance(tF, s);
+vE = distance(vc, s);
 
 var uv: vec2f;
 if (uvProjection > 0) {
 let n = normalize(_nA * b.x + _nB * b.y + _nC * b.z);
-uv = gI(uvProjection, tF - (_nA * b.x + _nB * b.y + _nC * b.z) * uJ, n) * uvScale;
+uv = gE(uvProjection, vc - (_nA * b.x + _nB * b.y + _nC * b.z) * wG, n) * uvScale;
 }
 else {
-uv = vG * b.x + vH * b.y + vI * b.z;
+uv = xQ * b.x + xR * b.y + xS * b.z;
 }
-let h = sampleMap(jQ, uv).x;
+let h = sampleMap(jW, uv).x;
 
-vd = max(uP * h, DT);
+hSurf = max(wZ * h, DT);
 
 
-let dh = dot(te, s - c0);
+let dh = dot(uG, s - c0);
 c0 = c0 + _nA * dh / dA;
 c1 = c1 + _nB * dh / dB;
 c2 = c2 + _nC * dh / dC;
 
-pY = qZ;
-qZ = tS - vd;
+qL = hDiff;
+hDiff = vE - hSurf;
 
-if (!(select(qZ * ra, qZ * pY, rc) > 0.0
-&& t < bS + dt) 
+if (!(select(hDiff * sb, hDiff * qL, se) > 0.0
+&& t < bN + dt) 
 
 && step > 0) {
 break;
@@ -301,15 +301,15 @@ break;
 
 step++;
 t = t + dt;
-s = s + qM * dt;
+s = s + rP * dt;
 }
 
-if (select(qZ * ra, qZ * pY, rc) <= 0.0) {
-let dtt = dt * (abs(qZ) / (abs(pY) + abs(qZ) + EPS));
+if (select(hDiff * sb, hDiff * qL, se) <= 0.0) {
+let dtt = dt * (abs(hDiff) / (abs(qL) + abs(hDiff) + EPS));
 t = t - dtt;
 
-let wO = qr + qM * t;
-let b = lP(wO, c0, c1, c2);
+let yP = rl + rP * t;
+let b = mg(yP, c0, c1, c2);
 
 let N = normalize(_nA * b.x + _nB * b.y + _nC * b.z);
 
@@ -317,33 +317,33 @@ let p_a = _pA * b.x + _pB * b.y + _pC * b.z;
 let p_b = _pA * (b.x + dbx) + _pB * b.y + _pC * (b.z - dbx);
 let p_c = _pA * b.x + _pB * (b.y + dby) + _pC * (b.z - dby);
 
-var wg: vec2f;
-var wh: vec2f;
-var wi: vec2f;
+var yA: vec2f;
+var yB: vec2f;
+var yC: vec2f;
 
 if (uvProjection > 0) {
 let NA = normalize(_nA * b.x + _nB * b.y + _nC * b.z);
 let NB = normalize(_nA * (b.x + dbx) + _nB * b.y + _nC * (b.z - dbx));
 let NC = normalize(_nA * b.x + _nB * (b.y + dby) + _nC * (b.z - dby));  
-wg = gI(uvProjection, p_a - (_nA * b.x + _nB * b.y + _nC * b.z) * uJ, NA) * uvScale;
-wh = gI(uvProjection, p_b - (_nA * (b.x + dbx) + _nB * b.y + _nC * (b.z - dbx)) * uJ, NB) * uvScale;
-wi = gI(uvProjection, p_c - (_nA * b.x + _nB * (b.y + dby) + _nC * (b.z - dby)) * uJ, NC) * uvScale;
+yA = gE(uvProjection, p_a - (_nA * b.x + _nB * b.y + _nC * b.z) * wG, NA) * uvScale;
+yB = gE(uvProjection, p_b - (_nA * (b.x + dbx) + _nB * b.y + _nC * (b.z - dbx)) * wG, NB) * uvScale;
+yC = gE(uvProjection, p_c - (_nA * b.x + _nB * (b.y + dby) + _nC * (b.z - dby)) * wG, NC) * uvScale;
 }
 else {
-wg = vG * b.x + vH * b.y + vI * b.z;
-wh = vG * (b.x + dbx) + vH * b.y + vI * (b.z - dbx);
-wi = vG * b.x + vH * (b.y + dby) + vI * (b.z - dby);
+yA = xQ * b.x + xR * b.y + xS * b.z;
+yB = xQ * (b.x + dbx) + xR * b.y + xS * (b.z - dbx);
+yC = xQ * b.x + xR * (b.y + dby) + xS * (b.z - dby);
 }
 
-let s_a = p_a + ra * N * sampleMap(jQ, wg).x * hScale;
-let s_b = p_b + ra * N * sampleMap(jQ, wh).x * hScale;
-let s_c = p_c + ra * N * sampleMap(jQ, wi).x * hScale;
+let s_a = p_a + sb * N * sampleMap(jW, yA).x * hScale;
+let s_b = p_b + sb * N * sampleMap(jW, yB).x * hScale;
+let s_c = p_c + sb * N * sampleMap(jW, yC).x * hScale;
 
-let N_s = -ra * normalize(cross(s_c - s_a, s_b - s_a));
+let N_s = -sb * normalize(cross(s_c - s_a, s_b - s_a));
 
 result.t = t;
 result.barycentric = b;
-result.normal = normalize(N_s - te + ra * N);
+result.normal = normalize(N_s - uG + sb * N);
 }
 
 return result;

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

@@ -4,13 +4,13 @@
 export default /* wgsl */ `
 var<private> s0: vec4u;
 
-fn kh(p: vec2u, frame: i32) {
+fn fi(p: vec2u, frame: i32) {
 
 s0 = vec4u(p, u32(frame), u32(p.x) + u32(p.y));
 }
 
 
-fn vx() {
+fn xF() {
 s0 = s0 * 1664525u + 1013904223u;
 s0.x += s0.y*s0.w; s0.y += s0.z*s0.x; s0.z += s0.x*s0.y; s0.w += s0.y*s0.z;
 s0 = s0 ^ (s0>>vec4(16u));
@@ -18,7 +18,7 @@ s0.x += s0.y*s0.w; s0.y += s0.z*s0.x; s0.z += s0.x*s0.y; s0.w += s0.y*s0.z;
 }
 
 fn rand() -> f32 {
-vx(); return f32(s0.x)/f32(0xffffffffu);
+xF(); return f32(s0.x)/f32(0xffffffffu);
 }
 
 `;

package/lib/renderer/wgsl/chunks/rayTraceCore.wgsl.d.ts

@@ -1,5 +1,7 @@
-export default function rayTraceCore({ hasAnisotropyDirections, hasCoatNormals, hasSceneAlpha }: {
+export default function rayTraceCore({ hasAnisotropyDirections, hasCoatNormals, hasSceneAlpha, caustics, lightCapabilities }: {
     hasAnisotropyDirections?: boolean;
     hasCoatNormals?: boolean;
     hasSceneAlpha?: boolean;
+    caustics?: boolean;
+    lightCapabilities?: number;
 }): string;

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

@@ -1,7 +1,8 @@
-import { materialBufferWgsl } from '../../materialLayout.js';
+import { materialBufferWgsl } from '../../materialLayoutSchema.js';
+import { ALL_LIGHT_SHADER_CAPABILITIES, LIGHT_SHADER_CAPABILITIES } from '../../shaderCapabilities.js';
 import { ShadowCatcherMaterial, UnlitMaterial } from '../../../constants.js';
-import { CONSTANT_DIRECTIONAL_LIGHT, CONSTANT_RECT_LIGHT, CONSTANT_SPHERE_LIGHT, CONSTANT_SPOT_LIGHT } from '../../lightsUtil.js';
-export default function rayTraceCore({ hasAnisotropyDirections = true, hasCoatNormals = false, hasSceneAlpha = true }) {
+import { CONSTANT_DIRECTIONAL_LIGHT, CONSTANT_RECT_LIGHT, CONSTANT_SPHERE_LIGHT, CONSTANT_SPOT_LIGHT } from '../../lightConstants.js';
+export default function rayTraceCore({ hasAnisotropyDirections = true, hasCoatNormals = false, hasSceneAlpha = true, caustics = false, lightCapabilities = ALL_LIGHT_SHADER_CAPABILITIES }) {
     return /*wgsl*/ `
 
 const luminance = vec3f(0.2126, 0.7152, 0.0722);
@@ -15,13 +16,18 @@ const RECT_LIGHT = ${CONSTANT_RECT_LIGHT};
 const SPHERE_LIGHT = ${CONSTANT_SPHERE_LIGHT};
 const DIRECTIONAL_LIGHT = ${CONSTANT_DIRECTIONAL_LIGHT};
 const SPOT_LIGHT = ${CONSTANT_SPOT_LIGHT};
+const HAS_SCENE_SPHERE_LIGHT = ${(lightCapabilities & LIGHT_SHADER_CAPABILITIES.SphereLight) !== 0};
+const HAS_SCENE_RECT_LIGHT = ${(lightCapabilities & LIGHT_SHADER_CAPABILITIES.RectLight) !== 0};
+const HAS_SCENE_SPOT_LIGHT = ${(lightCapabilities & LIGHT_SHADER_CAPABILITIES.SpotLight) !== 0};
+const HAS_SCENE_DIRECTIONAL_LIGHT = ${(lightCapabilities & LIGHT_SHADER_CAPABILITIES.DirectionalLight) !== 0};
+const HAS_SCENE_EMISSIVE_LIGHT = ${(lightCapabilities & LIGHT_SHADER_CAPABILITIES.emissive) !== 0};
 
 struct Ray {
 o: vec3f,
 d: vec3f,
 };
 
-struct ad {
+struct SurfaceInteraction {
 position: vec3f,
 
 
@@ -56,6 +62,8 @@ barycentrics: vec2f,
 
 wavelength: f32,
 
+spectralChannel: i32,
+
 mediumType: i32,
 mediumExtinction: vec3f,
 mediumScattering: vec3f,
@@ -66,7 +74,7 @@ mediumAnisotropy: f32,
 
 
 
-struct J {
+struct OpenPbrResolvedMaterial {
 baseColor: vec3f,
 transmissionWeight: f32,
 
@@ -75,7 +83,9 @@ specularAnisotropy: f32,
 
 baseMetalness: f32,
 specularIor: f32,
-specularRoughness: f32,
+specularRoughness: f32,${caustics ? /* wgsl */ `
+causticRoughness: f32,
+causticNormalMapped: bool,` : ''}
 
 transmissionDispersionScale: f32,
 dispersionAbbeNumber: f32,
@@ -147,7 +157,13 @@ fromEye: bool,
 abort: bool,
 scatteringPdf: f32,
 scatteringIsDelta: bool,
-isNextRefraction: bool,
+isNextRefraction: bool,${caustics ? /* wgsl */ `
+causticK: u32,
+causticSingleT: bool,
+causticReceiverPrimary: bool,
+causticHasReceiver: bool,
+causticPhotonReceiver:bool,
+causticMneeReceiver: bool,` : ''}
 };
 
 struct LightSampleRec {
@@ -203,7 +219,7 @@ secondaryUv1: vec2f,
 secondaryUv2: vec2f,
 };
 
-struct pe {
+struct pW {
 probability: f32,
 aliasIndex: i32,
 sourceType: i32,
@@ -250,8 +266,12 @@ sceneFeatureFlags: u32,
 sequenceSeed: u32,
 zSobolLogSPP: u32,
 sampleSequenceOffset: u32,
+
+shadowCatcherChoiceOffset: i32,
+shadowCatcherChoiceCount: i32,
 environmentRotation: mat4x4f,
 backgroundMapTransform: vec4f,
+backgroundMapRotation: mat4x4f,
 }
 
 
@@ -273,44 +293,44 @@ backgroundMapTransform: vec4f,
 @group(1) @binding(9) var openPbrEnergyLut: texture_2d<f32>;
 @group(1) @binding(10) var adaptiveSampleTargets: texture_2d<f32>;
 
-fn cL(dV: i32, ls: i32) -> f32 {
-let im = dV + ls / 4;
+fn cF(dV: i32, lK: i32) -> f32 {
+let im = dV + lK / 4;
 let width = rtUniform.lightSamplingTextureWidth;
 let packed = textureLoad(lightSamplingData, vec2i(im % width, im / width), 0);
-return packed[ls % 4];
+return packed[lK % 4];
 }
 
-fn kR(dV: i32, uv: vec2f) -> f32 {
+fn ld(dV: i32, uv: vec2f) -> f32 {
 if (dV < 0) { return 1.0; }
-let fQ = fract(uv) * 128.0 - vec2f(0.5);
-let x0 = i32(floor(fQ.x));
-let y0 = i32(floor(fQ.y));
-let fx = fract(fQ.x);
-let fy = fract(fQ.y);
+let pixel = fract(uv) * 128.0 - vec2f(0.5);
+let x0 = i32(floor(pixel.x));
+let y0 = i32(floor(pixel.y));
+let fx = fract(pixel.x);
+let fy = fract(pixel.y);
 let xa = select(x0, 127, x0 < 0);
 let xb = select(x0 + 1, 0, x0 >= 127);
 let ya = select(y0, 127, y0 < 0);
 let yb = select(y0 + 1, 0, y0 >= 127);
-let a = cL(dV, ya * 128 + xa);
-let b = cL(dV, ya * 128 + xb);
-let c = cL(dV, yb * 128 + xa);
-let d = cL(dV, yb * 128 + xb);
+let a = cF(dV, ya * 128 + xa);
+let b = cF(dV, ya * 128 + xb);
+let c = cF(dV, yb * 128 + xa);
+let d = cF(dV, yb * 128 + xb);
 return mix(mix(a, b, fx), mix(c, d, fx), fy);
 }
 
-fn bL(light: Light, dK: vec3f) -> f32 {
+fn evaluateSpotEmissionMultiplier(light: Light, dI: vec3f) -> f32 {
 let axis = normalize(light.u);
-let cc = dot(axis, dK);
-var tP = select(0.0, 1.0, cc >= light.v.x);
+let ct = dot(axis, dI);
+var vw = select(0.0, 1.0, ct >= light.v.x);
 if (light.v.y > light.v.x + 1e-6) {
-tP = smoothstep(light.v.x, light.v.y, cc);
+vw = smoothstep(light.v.x, light.v.y, ct);
 }
-if (tP <= 0.0 || light.textureIndex < 0.0) { return tP; }
-let uj = light.iesReference;
-let vD = cross(axis, uj);
-let vertical = acos(clamp(cc, -1.0, 1.0));
-let horizontal = atan2(dot(dK, vD), dot(dK, uj));
-return tP * kR(i32(light.textureIndex), vec2f(vertical / PI, horizontal / TWOPI + 0.5));
+if (vw <= 0.0 || light.textureIndex < 0.0) { return vw; }
+let vX = light.iesReference;
+let xN = cross(axis, vX);
+let vertical = acos(clamp(ct, -1.0, 1.0));
+let horizontal = atan2(dot(dI, xN), dot(dI, vX));
+return vw * ld(i32(light.textureIndex), vec2f(vertical / PI, horizontal / TWOPI + 0.5));
 }
 
 
@@ -322,7 +342,7 @@ return tP * kR(i32(light.textureIndex), vec2f(vertical / PI, horizontal / TWOPI
 @group(2) @binding(2) var<uniform> analyticLights: array<Light, NUM_UNIFORM_LIGHTS>;
 
 
-fn iC(ray: ptr<function, Ray>, origin: vec3f, direction: vec3f) {
+fn initRay(ray: ptr<function, Ray>, origin: vec3f, direction: vec3f) {
 (*ray).o = origin;
 (*ray).d = direction;
 }