@meimu/mutracer 0.2.0 → 0.3.0

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

@@ -1,56 +1,56 @@
 const displacementCommon = /*wgsl*/ `
 
 
-struct kc {
+struct jI {
 barycentric: vec3f,
 t: f32,
 normal: vec3f,
 
 }
 
-fn wO(
+fn xp(
 uu: f32,
-vA: vec3f, vB: vec3f,
+vK: vec3f, vL: vec3f,
 eAC: vec3f, eBD: vec3f,
-je: vec3f,
-qT: f32,
-qV: ptr<function, f32>
+iM: vec3f,
+qK: f32,
+qM: ptr<function, f32>
 ) -> bool {
 if (uu >= 0.0 && uu <= 1.0) {
-let pAB = mix(vA, vB, uu);
+let pAB = mix(vK, vL, uu);
 let pCD = mix(eAC, eBD, uu);
-var n = cross(je, pCD);
-let xf = 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) * xf;
-let vv = dot(n, je) * xf;
-if (vv >= 0.0 && vv <= 1.0 && tt > qT && tt < *qV) {
-*qV = 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 hx(
-qF: vec3f, je: vec3f, qT: f32, uM: f32,
+fn hc(
+qx: vec3f, iM: vec3f, qK: f32, uR: f32,
 pA: vec3f, pB: vec3f, pC: vec3f, pD: vec3f,
-qV: 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 vA = pA - qF;
-let vB = pB - qF;
+let vK = pA - qx;
+let vL = pB - qx;
 
 var u1: f32;
 var u2: f32;
 {
 let qN = cross(eAB, -eCD);
-let a = dot(qN, je);
-var b = dot(cross(vB, je), eBD);
-let c = dot(cross(vA, je), 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 xF = sqrt(det);
-let yo = -0.5 * (b + sign(b) * xF);
-u1 = yo / a;
-u2 = c / yo;
+let yg = sqrt(det);
+let yN = -0.5 * (b + sign(b) * yg);
+u1 = yN / a;
+u2 = c / yN;
 }
 }
 
-*qV = uM;
+*qM = uR;
 
 var hit = false;
 
-hit = wO(u1, vA, vB, eAC, eBD, je, qT, qV) || 
-wO(u2, vA, vB, eAC, eBD, je, qT, qV);
+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 tM(
-qF: vec3f, je: vec3f, qT: f32, uM: f32,
+fn tO(
+qx: vec3f, iM: vec3f, qK: f32, uR: f32,
 pA: vec3f, pB: vec3f, pC: vec3f,
 pD: vec3f, pE: vec3f, pF: vec3f,
 
-bX: ptr<function, f32>,
-bW: ptr<function, f32>
+bO: ptr<function, f32>,
+bN: ptr<function, f32>
 ) -> bool {
-*bX = INF;
-*bW = -INF;
+*bO = INF;
+*bN = -INF;
 
-let xV = maxDimension(abs(je));
-let gY = nk(je, xV);
+let yv = maxDimension(abs(iM));
+let gH = mU(iM, yv);
 var tt: f32;
 
-tt = ga(
-qF, je,
-fP(pC, pB, pA), xV, gY, INF
+tt = fM(
+qx, iM,
+fw(pC, pB, pA), yv, gH, INF
 ).t;
 if (tt > 0.0) {
-*bX = min(*bX, tt);
-*bW = max(*bW, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-tt = ga(
-qF, je,
-fP(pD, pE, pF), xV, gY, INF
+tt = fM(
+qx, iM,
+fw(pD, pE, pF), yv, gH, INF
 ).t;
 if (tt > 0.0) {
-*bX = min(*bX, tt);
-*bW = max(*bW, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hx(
-qF, je, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pA, pB, pD, pE,
 &tt)) {
-*bX = min(*bX, tt);
-*bW = max(*bW, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hx(
-qF, je, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pB, pC, pE, pF,
 &tt)) {
-*bX = min(*bX, tt);
-*bW = max(*bW, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
-if (hx(
-qF, je, -INF, INF,
+if (hc(
+qx, iM, -INF, INF,
 pC, pA, pF, pD,
 &tt)) {
-*bX = min(*bX, tt);
-*bW = max(*bW, tt);
+*bO = min(*bO, tt);
+*bN = max(*bN, tt);
 }
 
-*bX = max(*bX, qT);
-*bW = min(*bW, uM);
+*bO = max(*bO, qK);
+*bN = min(*bN, uR);
 
-return *bX <= *bW && *bW > 0.0;
+return *bO <= *bN && *bN > 0.0;
 }
 `;
 export default /* wgsl */ `
 
 ${displacementCommon}
 
-fn mz(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 vs = 1.0 / (d00 * d11 - d01 * d01);
-let v = clamp((d11 * d20 - d01 * d21) * vs, 0.0, 1.0);
-let w = clamp((d00 * d21 - d01 * d20) * vs, 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 pY(
-rs: vec3f, rS: vec3f, yt: f32,
+fn pS(
+rl: vec3f, rP: vec3f, yR: f32,
 _pA: vec3f, _pB: vec3f, _pC: vec3f,
 _nA: vec3f, _nB: vec3f, _nC: vec3f,
-xm: vec2f, xn: vec2f, xo: vec2f,
+xQ: vec2f, xR: vec2f, xS: vec2f,
 
-ku: i32,
-xl: vec2f,
+jW: i32,
+xP: vec2f,
 
-yd: f32, yc: f32,
+yE: f32, yD: f32,
 hScale: f32,
-wk: f32,
-vH: f32,
+wG: f32,
+vR: f32,
 uvProjection: i32,
 uvScale: vec2f,
-xM: bool,
-sk: bool
-) -> kc {
+ym: bool,
+se: bool
+) -> jI {
 
 
 
-let DT = select(max(abs(hScale) / 32.0, 1e-12), vH, vH > 0.0);
-let dt = DT / length(rS);
+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: kc;
+var result: jI;
 result.t = INF;
 
-var bX = INF;
-var bW = -INF;
+var bO = INF;
+var bN = -INF;
 
-let sg = sign(hScale);
-let wC = abs(hScale);
+let sb = sign(hScale);
+let wZ = abs(hScale);
 
-let uC = sg * normalize(cross(_pB - _pA, _pC - _pA));
+let uG = sb * normalize(cross(_pB - _pA, _pC - _pA));
 
 
-let dA = dot(_nA, uC);
-let dB = dot(_nB, uC);
-let dC = dot(_nC, uC);
+let dA = dot(_nA, uG);
+let dB = dot(_nB, uG);
+let dC = dot(_nC, uG);
 
-var minH = yd * hScale * min(min(dA, dB), dC);
-var maxH = yc * 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 vC = abs(maxH) - abs(minH);
+let vM = abs(maxH) - abs(minH);
 
 
 minH = sign(minH) * max(abs(minH) - DT * 2.0, 0.0);
-if (vC < 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 (!tM(
-rs, rS, 0.0, INF,
+if (!tO(
+rl, rP, 0.0, INF,
 pA, pB, pC,
 pD, pE, pF,
 
-&bX, &bW
+&bO, &bN
 )) {
 return result;
 }
 
-var yp = false;
+var yO = false;
 
-if (bW == bX) {
-bX = 0.0;
+if (bN == bO) {
+bO = 0.0;
 }
 
-var t = bX + dt * eC(2u);
-var s = rs + rS * t;
-var vq = dot(s - _pA, uC);
-var wV = select(0.0, INF, sg < 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 * vq / dA;
-var c1 = _pB + _nB * vq / dB;
-var c2 = _pC + _nC * vq / dC;
+var c0 = _pA + _nA * vE / dA;
+var c1 = _pB + _nB * vE / dB;
+var c2 = _pC + _nC * vE / dC;
 
-var uV: vec3f;
+var vc: vec3f;
 
-vq = abs(vq);
+vE = abs(vE);
 
-var sf = vq - wV;
-var qU = sf;
+var hDiff = vE - hSurf;
+var qL = hDiff;
 var step = 0;
 
 while (true) {
 
-let b = mz(s, c0, c1, c2);
+let b = mg(s, c0, c1, c2);
 
-uV = _pA * b.x + _pB * b.y + _pC * b.z;
+vc = _pA * b.x + _pB * b.y + _pC * b.z;
 
-vq = distance(uV, s);
+vE = distance(vc, s);
 
 var uv: vec2f;
 if (uvProjection > 0) {
 let n = normalize(_nA * b.x + _nB * b.y + _nC * b.z);
-uv = gU(uvProjection, uV - (_nA * b.x + _nB * b.y + _nC * b.z) * wk, n) * uvScale;
+uv = gE(uvProjection, vc - (_nA * b.x + _nB * b.y + _nC * b.z) * wG, n) * uvScale;
 }
 else {
-uv = xm * b.x + xn * b.y + xo * b.z;
+uv = xQ * b.x + xR * b.y + xS * b.z;
 }
-let h = sampleMap(ku, uv).x;
+let h = sampleMap(jW, uv).x;
 
-wV = max(wC * h, DT);
+hSurf = max(wZ * h, DT);
 
 
-let dh = dot(uC, s - c0);
+let dh = dot(uG, s - c0);
 c0 = c0 + _nA * dh / dA;
 c1 = c1 + _nB * dh / dB;
 c2 = c2 + _nC * dh / dC;
 
-qU = sf;
-sf = vq - wV;
+qL = hDiff;
+hDiff = vE - hSurf;
 
-if (!(select(sf * sg, sf * qU, sk) > 0.0
-&& t < bW + 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 + rS * dt;
+s = s + rP * dt;
 }
 
-if (select(sf * sg, sf * qU, sk) <= 0.0) {
-let dtt = dt * (abs(sf) / (abs(qU) + abs(sf) + EPS));
+if (select(hDiff * sb, hDiff * qL, se) <= 0.0) {
+let dtt = dt * (abs(hDiff) / (abs(qL) + abs(hDiff) + EPS));
 t = t - dtt;
 
-let yr = rs + rS * t;
-let b = mz(yr, 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 xX: vec2f;
-var xY: vec2f;
-var xZ: 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));  
-xX = gU(uvProjection, p_a - (_nA * b.x + _nB * b.y + _nC * b.z) * wk, NA) * uvScale;
-xY = gU(uvProjection, p_b - (_nA * (b.x + dbx) + _nB * b.y + _nC * (b.z - dbx)) * wk, NB) * uvScale;
-xZ = gU(uvProjection, p_c - (_nA * b.x + _nB * (b.y + dby) + _nC * (b.z - dby)) * wk, 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 {
-xX = xm * b.x + xn * b.y + xo * b.z;
-xY = xm * (b.x + dbx) + xn * b.y + xo * (b.z - dbx);
-xZ = xm * b.x + xn * (b.y + dby) + xo * (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 + sg * N * sampleMap(ku, xX).x * hScale;
-let s_b = p_b + sg * N * sampleMap(ku, xY).x * hScale;
-let s_c = p_c + sg * N * sampleMap(ku, xZ).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 = -sg * 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 - uC + sg * 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 gW(p: vec2u, frame: i32) {
+fn fi(p: vec2u, frame: i32) {
 
 s0 = vec4u(p, u32(frame), u32(p.x) + u32(p.y));
 }
 
 
-fn xd() {
+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 {
-xd(); return f32(s0.x)/f32(0xffffffffu);
+xF(); return f32(s0.x)/f32(0xffffffffu);
 }
 
 `;

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

@@ -1,6 +1,7 @@
-export default function rayTraceCore({ hasAnisotropyDirections, hasCoatNormals, hasSceneAlpha, lightCapabilities }: {
+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

@@ -2,7 +2,7 @@ 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 '../../lightConstants.js';
-export default function rayTraceCore({ hasAnisotropyDirections = true, hasCoatNormals = false, hasSceneAlpha = true, lightCapabilities = ALL_LIGHT_SHADER_CAPABILITIES }) {
+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);
@@ -62,6 +62,8 @@ barycentrics: vec2f,
 
 wavelength: f32,
 
+spectralChannel: i32,
+
 mediumType: i32,
 mediumExtinction: vec3f,
 mediumScattering: vec3f,
@@ -81,7 +83,9 @@ specularAnisotropy: f32,
 
 baseMetalness: f32,
 specularIor: f32,
-specularRoughness: f32,
+specularRoughness: f32,${caustics ? /* wgsl */ `
+causticRoughness: f32,
+causticNormalMapped: bool,` : ''}
 
 transmissionDispersionScale: f32,
 dispersionAbbeNumber: f32,
@@ -153,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 {
@@ -209,7 +219,7 @@ secondaryUv1: vec2f,
 secondaryUv2: vec2f,
 };
 
-struct qb {
+struct pW {
 probability: f32,
 aliasIndex: i32,
 sourceType: i32,
@@ -283,15 +293,15 @@ backgroundMapRotation: mat4x4f,
 @group(1) @binding(9) var openPbrEnergyLut: texture_2d<f32>;
 @group(1) @binding(10) var adaptiveSampleTargets: texture_2d<f32>;
 
-fn cS(el: i32, mc: i32) -> f32 {
-let iK = el + mc / 4;
+fn cF(dV: i32, lK: i32) -> f32 {
+let im = dV + lK / 4;
 let width = rtUniform.lightSamplingTextureWidth;
-let packed = textureLoad(lightSamplingData, vec2i(iK % width, iK / width), 0);
-return packed[mc % 4];
+let packed = textureLoad(lightSamplingData, vec2i(im % width, im / width), 0);
+return packed[lK % 4];
 }
 
-fn lB(el: i32, uv: vec2f) -> f32 {
-if (el < 0) { return 1.0; }
+fn ld(dV: i32, uv: vec2f) -> f32 {
+if (dV < 0) { return 1.0; }
 let pixel = fract(uv) * 128.0 - vec2f(0.5);
 let x0 = i32(floor(pixel.x));
 let y0 = i32(floor(pixel.y));
@@ -301,26 +311,26 @@ 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 = cS(el, ya * 128 + xa);
-let b = cS(el, ya * 128 + xb);
-let c = cS(el, yb * 128 + xa);
-let d = cS(el, 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 bP(light: Light, dV: vec3f) -> f32 {
+fn evaluateSpotEmissionMultiplier(light: Light, dI: vec3f) -> f32 {
 let axis = normalize(light.u);
-let ci = dot(axis, dV);
-var vm = select(0.0, 1.0, ci >= 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) {
-vm = smoothstep(light.v.x, light.v.y, ci);
+vw = smoothstep(light.v.x, light.v.y, ct);
 }
-if (vm <= 0.0 || light.textureIndex < 0.0) { return vm; }
-let vK = light.iesReference;
-let xj = cross(axis, vK);
-let vertical = acos(clamp(ci, -1.0, 1.0));
-let horizontal = atan2(dot(dV, xj), dot(dV, vK));
-return vm * lB(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));
 }
 
 
@@ -332,7 +342,7 @@ return vm * lB(i32(light.textureIndex), vec2f(vertical / PI, horizontal / TWOPI
 @group(2) @binding(2) var<uniform> analyticLights: array<Light, NUM_UNIFORM_LIGHTS>;
 
 
-fn gC(ray: ptr<function, Ray>, origin: vec3f, direction: vec3f) {
+fn initRay(ray: ptr<function, Ray>, origin: vec3f, direction: vec3f) {
 (*ray).o = origin;
 (*ray).d = direction;
 }