@meimu/mutracer 0.1.5 → 0.3.0

package/lib/renderer/wgsl/aovDisplay.wgsl.js

@@ -0,0 +1,70 @@
+// Independent of beauty composition and denoising; shared by WebGPU and native.
+export function aovDisplayShader(format) {
+    return /* wgsl */ `
+struct DisplayAovAccum { geometry: vec4f, position: vec4f, diffuse: vec4f, specular: vec4f, }
+struct DisplayAovCenter { identity: vec4u, motion: vec4f, }
+struct DisplayAovParams { sizeKind: vec4u, mapping: vec4f, flags: vec4u, }
+@group(0) @binding(0) var<storage, read> accumulation: array<DisplayAovAccum>;
+@group(0) @binding(1) var<storage, read> center: array<DisplayAovCenter>;
+@group(0) @binding(2) var image: texture_storage_2d<${format}, write>;
+@group(0) @binding(3) var<uniform> params: DisplayAovParams;
+
+fn tG(input: u32) -> u32 {
+var v = input;
+v = (v ^ (v >> 16u)) * 0x7feb352du;
+v = (v ^ (v >> 15u)) * 0x846ca68bu;
+return v ^ (v >> 16u);
+}
+fn tH(rgb: vec3f) -> vec3f {
+return select(1.055 * pow(max(rgb, vec3f(0.0)), vec3f(1.0 / 2.4)) - 0.055,
+12.92 * rgb, rgb <= vec3f(0.0031308));
+}
+fn sA(rgb: vec3f) -> vec3f {
+return select(pow((rgb + 0.055) / 1.055, vec3f(2.4)), rgb / 12.92,
+rgb <= vec3f(0.04045));
+}
+@compute @workgroup_size(8, 8)
+fn main(@builtin(global_invocation_id) id: vec3u) {
+let size = params.sizeKind.xy;
+if (any(id.xy >= size)) { return; }
+let index = id.x + (size.y - 1u - id.y) * size.x;
+let kind = params.sizeKind.w;
+let pixel = accumulation[index];
+let hit = pixel.position.w;
+var color = vec3f(0.0);
+var alpha = hit / f32(max(params.sizeKind.z, 1u));
+let inv = 1.0 / max(hit, 1.0);
+if (kind == 0u) {
+color = vec3f(clamp((pixel.geometry.x * inv - params.mapping.x) /
+(params.mapping.y - params.mapping.x), 0.0, 1.0));
+} else if (kind == 1u) {
+color = clamp(pixel.position.xyz * inv / params.mapping.z * 0.5 + 0.5,
+vec3f(0.0), vec3f(1.0));
+} else if (kind == 2u) {
+color = clamp(pixel.geometry.yzw * inv, vec3f(0.0), vec3f(1.0));
+} else if (kind == 3u || kind == 4u) {
+let r = select(pixel.diffuse, pixel.specular, kind == 4u);
+alpha = r.w / f32(max(params.sizeKind.z, 1u));
+color = tH(clamp(r.xyz / max(r.w, 1.0), vec3f(0.0), vec3f(1.0)));
+} else if (kind == 5u) {
+let identity = center[index].identity;
+alpha = select(1.0, 0.0, identity.x == 0xffffffffu);
+let hash = tG(identity[params.flags.x]);
+color = 0.2 + 0.8 * vec3f(f32(hash & 255u), f32((hash >> 8u) & 255u),
+f32((hash >> 16u) & 255u)) / 255.0;
+} else {
+let m = center[index].motion;
+alpha = select(0.0, 1.0, m.w == 3.0);
+let gz = length(m.xy);
+var hue = 0.0;
+if (gz > 0.0) { hue = fract(atan2(m.y, m.x) / 6.28318530718 + 1.0); }
+let p = abs(fract(vec3f(hue) + vec3f(0.0, 2.0 / 3.0, 1.0 / 3.0)) * 6.0 - 3.0);
+color = clamp(p - 1.0, vec3f(0.0), vec3f(1.0)) * clamp(gz / params.mapping.w, 0.0, 1.0);
+}
+if (alpha <= 0.0) { color = vec3f(0.0); }
+if (params.flags.y != 0u) { color = sA(color); }
+textureStore(image, vec2i(id.xy), vec4f(color, clamp(alpha, 0.0, 1.0)));
+}
+`;
+}
+//# sourceMappingURL=aovDisplay.wgsl.js.map

package/lib/renderer/wgsl/aovDisplay.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovDisplay.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovDisplay.wgsl.ts"],"names":[],"mappings":"AAAA,gFAAgF;AAChF,MAAM,UAAU,gBAAgB,CAAC,MAAoC;IACnE,OAAO,UAAU,CAAC;;;;;;sDAMkC,MAAM;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA2D3D,CAAC;AACF,CAAC","sourcesContent":["// Independent of beauty composition and denoising; shared by WebGPU and native.\nexport function aovDisplayShader(format: 'rgba16float' | 'rgba8unorm') {\n  return /* wgsl */ `\nstruct DisplayAovAccum { geometry: vec4f, position: vec4f, diffuse: vec4f, specular: vec4f, }\nstruct DisplayAovCenter { identity: vec4u, motion: vec4f, }\nstruct DisplayAovParams { sizeKind: vec4u, mapping: vec4f, flags: vec4u, }\n@group(0) @binding(0) var<storage, read> accumulation: array<DisplayAovAccum>;\n@group(0) @binding(1) var<storage, read> center: array<DisplayAovCenter>;\n@group(0) @binding(2) var image: texture_storage_2d<${format}, write>;\n@group(0) @binding(3) var<uniform> params: DisplayAovParams;\n\nfn displayAovHash(input: u32) -> u32 {\n  var v = input;\n  v = (v ^ (v >> 16u)) * 0x7feb352du;\n  v = (v ^ (v >> 15u)) * 0x846ca68bu;\n  return v ^ (v >> 16u);\n}\nfn displayAovSrgb(rgb: vec3f) -> vec3f {\n  return select(1.055 * pow(max(rgb, vec3f(0.0)), vec3f(1.0 / 2.4)) - 0.055,\n    12.92 * rgb, rgb <= vec3f(0.0031308));\n}\nfn displayAovLinear(rgb: vec3f) -> vec3f {\n  return select(pow((rgb + 0.055) / 1.055, vec3f(2.4)), rgb / 12.92,\n    rgb <= vec3f(0.04045));\n}\n@compute @workgroup_size(8, 8)\nfn main(@builtin(global_invocation_id) id: vec3u) {\n  let size = params.sizeKind.xy;\n  if (any(id.xy >= size)) { return; }\n  let index = id.x + (size.y - 1u - id.y) * size.x;\n  let kind = params.sizeKind.w;\n  let pixel = accumulation[index];\n  let hit = pixel.position.w;\n  var color = vec3f(0.0);\n  var alpha = hit / f32(max(params.sizeKind.z, 1u));\n  let inv = 1.0 / max(hit, 1.0);\n  if (kind == 0u) {\n    color = vec3f(clamp((pixel.geometry.x * inv - params.mapping.x) /\n      (params.mapping.y - params.mapping.x), 0.0, 1.0));\n  } else if (kind == 1u) {\n    color = clamp(pixel.position.xyz * inv / params.mapping.z * 0.5 + 0.5,\n      vec3f(0.0), vec3f(1.0));\n  } else if (kind == 2u) {\n    color = clamp(pixel.geometry.yzw * inv, vec3f(0.0), vec3f(1.0));\n  } else if (kind == 3u || kind == 4u) {\n    let r = select(pixel.diffuse, pixel.specular, kind == 4u);\n    alpha = r.w / f32(max(params.sizeKind.z, 1u));\n    color = displayAovSrgb(clamp(r.xyz / max(r.w, 1.0), vec3f(0.0), vec3f(1.0)));\n  } else if (kind == 5u) {\n    let identity = center[index].identity;\n    alpha = select(1.0, 0.0, identity.x == 0xffffffffu);\n    let hash = displayAovHash(identity[params.flags.x]);\n    color = 0.2 + 0.8 * vec3f(f32(hash & 255u), f32((hash >> 8u) & 255u),\n      f32((hash >> 16u) & 255u)) / 255.0;\n  } else {\n    let m = center[index].motion;\n    alpha = select(0.0, 1.0, m.w == 3.0);\n    let magnitude = length(m.xy);\n    var hue = 0.0;\n    if (magnitude > 0.0) { hue = fract(atan2(m.y, m.x) / 6.28318530718 + 1.0); }\n    let p = abs(fract(vec3f(hue) + vec3f(0.0, 2.0 / 3.0, 1.0 / 3.0)) * 6.0 - 3.0);\n    color = clamp(p - 1.0, vec3f(0.0), vec3f(1.0)) * clamp(magnitude / params.mapping.w, 0.0, 1.0);\n  }\n  if (alpha <= 0.0) { color = vec3f(0.0); }\n  if (params.flags.y != 0u) { color = displayAovLinear(color); }\n  textureStore(image, vec2i(id.xy), vec4f(color, clamp(alpha, 0.0, 1.0)));\n}\n`;\n}\n"]}

package/lib/renderer/wgsl/aovReflectance.wgsl.d.ts

@@ -0,0 +1,2 @@
+declare const _default: string;
+export default _default;

package/lib/renderer/wgsl/aovReflectance.wgsl.js

@@ -0,0 +1,98 @@
+import proposals from './aovReflectanceSamples.wgsl.js';
+// Directional-hemispherical quadrature of the supported unlayered OpenPBR
+// reflection lobes. Keep the original 256 cosine/VNDF proposals and weights;
+// evaluate diffuse and reflection together without traversing the full layered
+// BSDF twice. Callers reject coat, fuzz, anisotropy, SSS and thin film first.
+export default proposals + /* wgsl */ `
+struct nQ {
+diffuse: vec3f,
+specular: vec3f,
+}
+fn hu(prepared: OpenPbrPreparedBsdf, qH: bool, kE: bool) -> nQ {
+var result = nQ();
+let view = prepared.viewLocal;
+if (view.z <= 1e-6) { return result; }
+let material = prepared.material;
+let alpha = vec2f(material.alphaX, material.alphaY);
+let lN = aD(view, alpha);
+let iR = prepared.weights[OPENPBR_LOBE_DIFFUSE];
+let fN = prepared.weights[OPENPBR_LOBE_METAL];
+let specularWeight = prepared.weights[OPENPBR_LOBE_SPECULAR];
+var dS = 1.0;
+if (kE && prepared.dielectricMsProbability > 0.0) {
+dS = max(aY(prepared.refractionEtaTOverEtaI, prepared.specularAlpha), 1e-12);
+}
+let ce = normalize(vec3f(alpha, 1.0) * view);
+let qG = min(ce.z, 1.0);
+let jN = 1.0 + qG;
+let count = 256u * min(rtUniform.loopUnrollGuard, 4u) / 4u;
+for (var i = 0u; i < count; i++) {
+let wn = i / 2u;
+let u = vec2f((f32(wn) + 0.5) / f32(count / 2u),
+fract(f32(reverseBits(wn)) * 2.3283064365386963e-10 + 0.37));
+let proposal = gV[wn];
+
+
+let hv = prepared.specularAlpha < 0.001 || count != 256u;
+var light = vec3f(proposal.xy, sqrt(max(0.0, 1.0 - dot(proposal.xy, proposal.xy))));
+if (hv) { light = aM(u); }
+if ((i & 1u) != 0u) {
+let kf = (1.0 - u.y) * jN;
+let dr = sqrt(max(0.0, (kf + (1.0 - qG)) * (u.y * jN)));
+var half = normalize(vec3f(alpha, 1.0) * vec3f(ce.xy + dr * proposal.zw, kf));
+if (hv) { half = aA(view, alpha, u); }
+light = reflect(-view, half);
+}
+
+if (light.z <= 0.0) { continue; }
+let half = normalize(view + light);
+let specularPdf = lN * sampling_ggx_ndf(half, alpha) / (4.0 * view.z);
+let pdf = 0.5 * (light.z * INVPI + specularPdf);
+
+
+let au = openpbr_to_local(prepared.basis, openpbr_to_world(prepared.basis, light));
+if (au.z <= 0.0) { continue; }
+let scale = abs(au.z) / (pdf * f32(count));
+if (qH && iR > 0.0) {
+let compensation = prepared.diffuseSpecularEnergyCompensation * ar(
+prepared.specularEtaTOverEtaI, prepared.specularAlpha, abs(au.z));
+var unusedPdf = 0.0;
+result.diffuse += iR * bw(material, prepared.diffuseColor,
+view, au, compensation, &unusedPdf) * scale;
+}
+if (kE) {
+let evaluatedHalf = normalize(view + au);
+let viewHalf = abs(dot(view, evaluatedHalf));
+var fresnel = vec3f(0.0);
+if (specularWeight > 0.0) {
+var dielectric = vec3f(ae(viewHalf, material.relativeIor));
+if (material.thinWalled) { dielectric = aq(dielectric, material.transmissionWeight); }
+fresnel += specularWeight * material.specularColor * dielectric;
+}
+if (fN > 0.0) {
+fresnel += fN * material.specularWeight * openpbr_metal_fresnel(viewHalf,
+material.baseWeight * material.baseColor, material.specularColor);
+}
+let rE = sampling_ggx_ndf(evaluatedHalf, alpha) * lN
+* aD(au, alpha) / (4.0 * max(au.z * view.z, 1e-8));
+var reflection = fresnel * rE;
+if (fN > 0.0 && prepared.metalMsFactor > 0.0) {
+let cv = min(prepared.metalMsFactor * aR(prepared.specularAlpha, au.z), 1.0 / max(au.z, 1e-8));
+reflection += fN * prepared.metalMsScale * cv * INVPI;
+}
+if (prepared.dielectricMsProbability > 0.0) {
+let oP = prepared.dielectricMsViewLoss * aV(
+prepared.refractionEtaTOverEtaI, prepared.specularAlpha, abs(au.z)) / dS;
+let cv = min(oP * prepared.dielectricMsReflectionRatio, 1.0 / abs(au.z));
+reflection += prepared.dielectricMsReflectionScale * cv * INVPI;
+}
+result.specular += reflection * scale;
+}
+}
+return result;
+}
+fn wl(prepared: OpenPbrPreparedBsdf) -> nQ {
+return hu(prepared, true, true);
+}
+`;
+//# sourceMappingURL=aovReflectance.wgsl.js.map

package/lib/renderer/wgsl/aovReflectance.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovReflectance.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovReflectance.wgsl.ts"],"names":[],"mappings":"AAAA,OAAO,SAAS,MAAM,iCAAiC,CAAC;AAExD,0EAA0E;AAC1E,6EAA6E;AAC7E,+EAA+E;AAC/E,8EAA8E;AAC9E,eAAe,SAAS,GAAG,UAAU,CAAC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA2FrC,CAAC","sourcesContent":["import proposals from './aovReflectanceSamples.wgsl.js';\n\n// Directional-hemispherical quadrature of the supported unlayered OpenPBR\n// reflection lobes. Keep the original 256 cosine/VNDF proposals and weights;\n// evaluate diffuse and reflection together without traversing the full layered\n// BSDF twice. Callers reject coat, fuzz, anisotropy, SSS and thin film first.\nexport default proposals + /* wgsl */ `\nstruct MtReflectance {\n  diffuse: vec3f,\n  specular: vec3f,\n}\nfn mtFrameReflectanceComponents(prepared: OpenPbrPreparedBsdf, wantDiffuse: bool, wantSpecular: bool) -> MtReflectance {\n  var result = MtReflectance();\n  let view = prepared.viewLocal;\n  if (view.z <= 1e-6) { return result; }\n  let material = prepared.material;\n  let alpha = vec2f(material.alphaX, material.alphaY);\n  let viewMasking = sampling_ggx_smith_g1(view, alpha);\n  let diffuseWeight = prepared.weights[OPENPBR_LOBE_DIFFUSE];\n  let metalWeight = prepared.weights[OPENPBR_LOBE_METAL];\n  let specularWeight = prepared.weights[OPENPBR_LOBE_SPECULAR];\n  var dielectricAverageLoss = 1.0;\n  if (wantSpecular && prepared.dielectricMsProbability > 0.0) {\n    dielectricAverageLoss = max(openpbr_ideal_dielectric_average_loss(prepared.refractionEtaTOverEtaI, prepared.specularAlpha), 1e-12);\n  }\n  let sphericalView = normalize(vec3f(alpha, 1.0) * view);\n  let viewZ = min(sphericalView.z, 1.0);\n  let capSpan = 1.0 + viewZ;\n  let count = 256u * min(rtUniform.loopUnrollGuard, 4u) / 4u;\n  for (var i = 0u; i < count; i++) {\n    let pair = i / 2u;\n    let u = vec2f((f32(pair) + 0.5) / f32(count / 2u),\n      fract(f32(reverseBits(pair)) * 2.3283064365386963e-10 + 0.37));\n    let proposal = mtReflectanceProposals[pair];\n    // Retain runtime proposals near delta lobes: tiny device-specific sin/cos\n    // rounding differences can otherwise be magnified by their narrow PDFs.\n    let runtimeProposal = prepared.specularAlpha < 0.001 || count != 256u;\n    var light = vec3f(proposal.xy, sqrt(max(0.0, 1.0 - dot(proposal.xy, proposal.xy))));\n    if (runtimeProposal) { light = sampling_cosine_hemisphere(u); }\n    if ((i & 1u) != 0u) {\n      let normalZ = (1.0 - u.y) * capSpan;\n      let radial = sqrt(max(0.0, (normalZ + (1.0 - viewZ)) * (u.y * capSpan)));\n      var half = normalize(vec3f(alpha, 1.0) * vec3f(sphericalView.xy + radial * proposal.zw, normalZ));\n      if (runtimeProposal) { half = sampling_ggx_visible_normal(view, alpha, u); }\n      light = reflect(-view, half);\n    }\n    // Rejected VNDF proposals retain their weight in the denominator.\n    if (light.z <= 0.0) { continue; }\n    let half = normalize(view + light);\n    let specularPdf = viewMasking * sampling_ggx_ndf(half, alpha) / (4.0 * view.z);\n    let pdf = 0.5 * (light.z * INVPI + specularPdf);\n    // Preserve the production evaluator's frame round trip, including FP32\n    // rounding near narrow highlights. Proposal PDFs use the original light.\n    let evaluatedLight = openpbr_to_local(prepared.basis, openpbr_to_world(prepared.basis, light));\n    if (evaluatedLight.z <= 0.0) { continue; }\n    let scale = abs(evaluatedLight.z) / (pdf * f32(count));\n    if (wantDiffuse && diffuseWeight > 0.0) {\n      let compensation = prepared.diffuseSpecularEnergyCompensation * openpbr_opaque_energy_complement(\n        prepared.specularEtaTOverEtaI, prepared.specularAlpha, abs(evaluatedLight.z));\n      var unusedPdf = 0.0;\n      result.diffuse += diffuseWeight * openpbr_eval_diffuse(material, prepared.diffuseColor,\n        view, evaluatedLight, compensation, &unusedPdf) * scale;\n    }\n    if (wantSpecular) {\n      let evaluatedHalf = normalize(view + evaluatedLight);\n      let viewHalf = abs(dot(view, evaluatedHalf));\n      var fresnel = vec3f(0.0);\n      if (specularWeight > 0.0) {\n        var dielectric = vec3f(openpbr_dielectric_fresnel(viewHalf, material.relativeIor));\n        if (material.thinWalled) { dielectric = openpbr_thin_wall_reflection_coefficient(dielectric, material.transmissionWeight); }\n        fresnel += specularWeight * material.specularColor * dielectric;\n      }\n      if (metalWeight > 0.0) {\n        fresnel += metalWeight * material.specularWeight * openpbr_metal_fresnel(viewHalf,\n          material.baseWeight * material.baseColor, material.specularColor);\n      }\n      let microfacet = sampling_ggx_ndf(evaluatedHalf, alpha) * viewMasking\n        * sampling_ggx_smith_g1(evaluatedLight, alpha) / (4.0 * max(evaluatedLight.z * view.z, 1e-8));\n      var reflection = fresnel * microfacet;\n      if (metalWeight > 0.0 && prepared.metalMsFactor > 0.0) {\n        let factor = min(prepared.metalMsFactor * openpbr_metal_energy_complement(prepared.specularAlpha, evaluatedLight.z), 1.0 / max(evaluatedLight.z, 1e-8));\n        reflection += metalWeight * prepared.metalMsScale * factor * INVPI;\n      }\n      if (prepared.dielectricMsProbability > 0.0) {\n        let factors = prepared.dielectricMsViewLoss * openpbr_ideal_dielectric_loss(\n          prepared.refractionEtaTOverEtaI, prepared.specularAlpha, abs(evaluatedLight.z)) / dielectricAverageLoss;\n        let factor = min(factors * prepared.dielectricMsReflectionRatio, 1.0 / abs(evaluatedLight.z));\n        reflection += prepared.dielectricMsReflectionScale * factor * INVPI;\n      }\n      result.specular += reflection * scale;\n    }\n  }\n  return result;\n}\nfn mtFrameReflectance(prepared: OpenPbrPreparedBsdf) -> MtReflectance {\n  return mtFrameReflectanceComponents(prepared, true, true);\n}\n`;\n"]}

package/lib/renderer/wgsl/aovReflectanceSamples.wgsl.d.ts

@@ -0,0 +1,2 @@
+declare const _default: "\nconst mtReflectanceProposals: array<vec4f, 128> = array<vec4f, 128>(\n  vec4f(-0.971247852, -0.202764958, 0.999698818, 0.0245412290),\n  vec4f(-0.808660269, 0.547487736, 0.997290432, 0.0735645667),\n  vec4f(-0.942509234, 0.187289074, 0.992479563, 0.122410677),\n  vec4f(-0.763038993, -0.558020830, 0.985277653, 0.170961902),\n  vec4f(-0.929654300, -0.00785388052, 0.975702107, 0.219101235),\n  vec4f(-0.656618595, 0.740912974, 0.963776052, 0.266712785),\n  vec4f(-0.817265987, 0.373183846, 0.949528158, 0.313681751),\n  vec4f(-0.793490529, -0.386950523, 0.932992816, 0.359895051),\n  vec4f(-0.860713601, -0.105764769, 0.914209723, 0.405241340),\n  vec4f(-0.604139447, 0.619064212, 0.893224299, 0.449611336),\n  vec4f(-0.787263036, 0.281084418, 0.870086968, 0.492898226),\n  vec4f(-0.673318744, -0.468566477, 0.844853580, 0.534997642),\n  vec4f(-0.799623907, 0.0901312605, 0.817584813, 0.575808227),\n  vec4f(-0.570307195, -0.676738262, 0.788346410, 0.615231633),\n  vec4f(-0.627450705, 0.452229142, 0.757208824, 0.653172851),\n  vec4f(-0.698281944, -0.294418186, 0.724247098, 0.689540565),\n  vec4f(-0.726036906, -0.153989345, 0.689540505, 0.724247098),\n  vec4f(-0.523753047, 0.608020604, 0.653172791, 0.757208884),\n  vec4f(-0.671607554, 0.233185455, 0.615231574, 0.788346410),\n  vec4f(-0.491990626, -0.494420350, 0.575808167, 0.817584813),\n  vec4f(-0.678437173, 0.0412081257, 0.534997642, 0.844853580),\n  vec4f(-0.495611340, -0.807542920, 0.492898196, 0.870086968),\n  vec4f(-0.508582473, 0.402262419, 0.449611306, 0.893224299),\n  vec4f(-0.537693322, -0.333672792, 0.405241281, 0.914209783),\n  vec4f(-0.614562631, -0.0568606183, 0.359894961, 0.932992816),\n  vec4f(-0.452296555, 0.809743226, 0.313681662, 0.949528217),\n  vec4f(-0.492342144, 0.317682207, 0.266712755, 0.963776052),\n  vec4f(-0.403601944, -0.406031638, 0.219101220, 0.975702107),\n  vec4f(-0.537260950, 0.137945175, 0.170961857, 0.985277653),\n  vec4f(-0.404928178, -0.715918541, 0.122410625, 0.992479563),\n  vec4f(-0.386338592, 0.556550741, 0.0735644922, 0.997290432),\n  vec4f(-0.445942074, -0.242918074, 0.0245411359, 0.999698818),\n  vec4f(-0.459759980, -0.175696611, -0.0245412234, 0.999698818),\n  vec4f(-0.359771132, 0.682195961, -0.0735645741, 0.997290432),\n  vec4f(-0.412575692, 0.205535084, -0.122410707, 0.992479503),\n  vec4f(-0.336475104, -0.646421731, -0.170961946, 0.985277653),\n  vec4f(-0.429363400, 0.0166855231, -0.219101310, 0.975702107),\n  vec4f(-0.319253087, -0.925218344, -0.266712844, 0.963776052),\n  vec4f(-0.294470191, 0.430103362, -0.313681722, 0.949528158),\n  vec4f(-0.281283021, -0.387403429, -0.359895051, 0.932992816),\n  vec4f(-0.358183086, -0.0808178857, -0.405241340, 0.914209723),\n  vec4f(-0.271769434, 0.854052365, -0.449611366, 0.893224299),\n  vec4f(-0.244775921, 0.311606824, -0.492898166, 0.870087028),\n  vec4f(-0.244355842, -0.552090824, -0.534997702, 0.844853520),\n  vec4f(-0.283291280, 0.112162955, -0.575808167, 0.817584813),\n  vec4f(-0.224362880, -0.823741674, -0.615231693, 0.788346350),\n  vec4f(-0.210826442, 0.610893846, -0.653172851, 0.757208824),\n  vec4f(-0.191608176, -0.297364026, -0.689540684, 0.724246979),\n  vec4f(-0.207833603, -0.124338157, -0.724247158, 0.689540505),\n  vec4f(-0.176593989, 0.814833522, -0.757208824, 0.653172851),\n  vec4f(-0.158949688, 0.293469012, -0.788346469, 0.615231514),\n  vec4f(-0.152046368, -0.648668587, -0.817584813, 0.575808167),\n  vec4f(-0.167782038, 0.0643178672, -0.844853640, 0.534997463),\n  vec4f(-0.128430068, -0.907204390, -0.870087028, 0.492898136),\n  vec4f(-0.115809105, 0.572147131, -0.893224299, 0.449611366),\n  vec4f(-0.103222348, -0.403248310, -0.914209783, 0.405241221),\n  vec4f(-0.112737618, -0.0319865458, -0.932992816, 0.359895051),\n  vec4f(-0.0796749964, 0.955433667, -0.949528217, 0.313681602),\n  vec4f(-0.0672519580, 0.453793734, -0.963776052, 0.266712725),\n  vec4f(-0.0551275201, -0.538435221, -0.975702167, 0.219101071),\n  vec4f(-0.0426490456, 0.204346836, -0.985277653, 0.170961812),\n  vec4f(-0.0306719355, -0.790655255, -0.992479503, 0.122410700),\n  vec4f(-0.0184057076, 0.708510995, -0.997290492, 0.0735644475),\n  vec4f(-0.00613544276, -0.291185349, -0.999698818, 0.0245412104),\n  vec4f(0.00613530306, -0.244297966, -0.999698818, -0.0245413855),\n  vec4f(0.0184059385, 0.755400777, -0.997290432, -0.0735646188),\n  vec4f(0.0306059923, 0.253786176, -0.992479503, -0.122410871),\n  vec4f(0.0429276451, -0.743136168, -0.985277653, -0.170961991),\n  vec4f(0.0701737553, 0.00441496586, -0.975702107, -0.219101235),\n  vec4f(0.0674433559, -0.992085159, -0.963776052, -0.266712904),\n  vec4f(0.0794365108, 0.499346048, -0.949528158, -0.313681781),\n  vec4f(0.0915081054, -0.485832185, -0.932992756, -0.359895200),\n  vec4f(0.101071075, -0.0861614645, -0.914209723, -0.405241400),\n  vec4f(0.116242222, 0.872919321, -0.893224180, -0.449611515),\n  vec4f(0.126407236, 0.359021723, -0.870086908, -0.492898285),\n  vec4f(0.139908314, -0.603366435, -0.844853580, -0.534997642),\n  vec4f(0.168981716, 0.0979395211, -0.817584693, -0.575808346),\n  vec4f(0.164669126, -0.853637457, -0.788346410, -0.615231633),\n  vec4f(0.175589874, 0.605686426, -0.757208705, -0.653172970),\n  vec4f(0.181917489, -0.321471483, -0.724247038, -0.689540625),\n  vec4f(0.219242185, -0.135647133, -0.689540565, -0.724247098),\n  vec4f(0.212299466, 0.790099680, -0.653172910, -0.757208765),\n  vec4f(0.208243251, 0.240495890, -0.615231574, -0.788346469),\n  vec4f(0.234419063, -0.640313387, -0.575808048, -0.817584932),\n  vec4f(0.315854102, 0.0532567836, -0.534997344, -0.844853759),\n  vec4f(0.260333687, -0.894772291, -0.492898196, -0.870086968),\n  vec4f(0.265374005, 0.502337217, -0.449611217, -0.893224359),\n  vec4f(0.267429054, -0.339113116, -0.405241072, -0.914209843),\n  vec4f(0.380209267, -0.0445682108, -0.359895110, -0.932992756),\n  vec4f(0.307221651, 0.891683638, -0.313681692, -0.949528217),\n  vec4f(0.296788067, 0.329054117, -0.266712576, -0.963776112),\n  vec4f(0.317195147, -0.457708418, -0.219100907, -0.975702226),\n  vec4f(0.419729054, 0.148764044, -0.170961887, -0.985277653),\n  vec4f(0.349995196, -0.727014899, -0.122410536, -0.992479563),\n  vec4f(0.356363803, 0.594497383, -0.0735642835, -0.997290492),\n  vec4f(0.434345096, -0.231501341, -0.0245412868, -0.999698818),\n  vec4f(0.480473697, -0.164373189, 0.0245413091, -0.999698818),\n  vec4f(0.392657876, 0.681903303, 0.0735647827, -0.997290432),\n  vec4f(0.492657542, 0.218807966, 0.122411035, -0.992479503),\n  vec4f(0.409054011, -0.584142148, 0.170961916, -0.985277653),\n  vec4f(0.569577277, 0.0289491210, 0.219101399, -0.975702107),\n  vec4f(0.442882210, -0.855257273, 0.266713053, -0.963775992),\n  vec4f(0.459734946, 0.387970626, 0.313681692, -0.949528217),\n  vec4f(0.513070345, -0.343044043, 0.359895140, -0.932992756),\n  vec4f(0.629031241, -0.0690753013, 0.405241549, -0.914209664),\n  vec4f(0.483215511, 0.773316741, 0.449611664, -0.893224120),\n  vec4f(0.586818159, 0.310843110, 0.492898226, -0.870086968),\n  vec4f(0.528679490, -0.427168757, 0.534997761, -0.844853461),\n  vec4f(0.683721542, 0.126429096, 0.575808465, -0.817584634),\n  vec4f(0.517971814, -0.658907175, 0.615231574, -0.788346410),\n  vec4f(0.548397303, 0.476606160, 0.653172910, -0.757208765),\n  vec4f(0.695181608, -0.259932309, 0.689540744, -0.724246919),\n  vec4f(0.748597383, -0.117820688, 0.724247038, -0.689540625),\n  vec4f(0.556980729, 0.641256869, 0.757208884, -0.653172791),\n  vec4f(0.741829991, 0.268901259, 0.788346529, -0.615231454),\n  vec4f(0.647817433, -0.477341115, 0.817584991, -0.575807929),\n  vec4f(0.816602290, 0.0779312477, 0.844853580, -0.534997642),\n  vec4f(0.599919200, -0.671730816, 0.870087087, -0.492898077),\n  vec4f(0.724726558, 0.447135508, 0.893224418, -0.449611068),\n  vec4f(0.810727239, -0.307791352, 0.914209723, -0.405241370),\n  vec4f(0.882583082, -0.0201240778, 0.932992816, -0.359894961),\n  vec4f(0.645550549, 0.729843259, 0.949528277, -0.313681543),\n  vec4f(0.839177251, 0.362342149, 0.963776171, -0.266712397),\n  vec4f(0.839526892, -0.399391353, 0.975702107, -0.219101220),\n  vec4f(0.928934336, 0.175204813, 0.985277653, -0.170961723),\n  vec4f(0.774228692, -0.569184363, 0.992479563, -0.122410372),\n  vec4f(0.815476179, 0.537282884, 0.997290432, -0.0735645965),\n  vec4f(0.968665719, -0.214761987, 0.999698818, -0.0245411228)\n);\n";
+export default _default;

package/lib/renderer/wgsl/aovReflectanceSamples.wgsl.js

@@ -0,0 +1,135 @@
+// Generated by scripts/generate-aov-reflectance-samples.mjs; do not edit.
+// xy: concentric cosine disk; zw: spherical-cap azimuth cosine/sine.
+export default /* wgsl */ `
+const gV: array<vec4f, 128> = array<vec4f, 128>(
+vec4f(-0.971247852, -0.202764958, 0.999698818, 0.0245412290),
+vec4f(-0.808660269, 0.547487736, 0.997290432, 0.0735645667),
+vec4f(-0.942509234, 0.187289074, 0.992479563, 0.122410677),
+vec4f(-0.763038993, -0.558020830, 0.985277653, 0.170961902),
+vec4f(-0.929654300, -0.00785388052, 0.975702107, 0.219101235),
+vec4f(-0.656618595, 0.740912974, 0.963776052, 0.266712785),
+vec4f(-0.817265987, 0.373183846, 0.949528158, 0.313681751),
+vec4f(-0.793490529, -0.386950523, 0.932992816, 0.359895051),
+vec4f(-0.860713601, -0.105764769, 0.914209723, 0.405241340),
+vec4f(-0.604139447, 0.619064212, 0.893224299, 0.449611336),
+vec4f(-0.787263036, 0.281084418, 0.870086968, 0.492898226),
+vec4f(-0.673318744, -0.468566477, 0.844853580, 0.534997642),
+vec4f(-0.799623907, 0.0901312605, 0.817584813, 0.575808227),
+vec4f(-0.570307195, -0.676738262, 0.788346410, 0.615231633),
+vec4f(-0.627450705, 0.452229142, 0.757208824, 0.653172851),
+vec4f(-0.698281944, -0.294418186, 0.724247098, 0.689540565),
+vec4f(-0.726036906, -0.153989345, 0.689540505, 0.724247098),
+vec4f(-0.523753047, 0.608020604, 0.653172791, 0.757208884),
+vec4f(-0.671607554, 0.233185455, 0.615231574, 0.788346410),
+vec4f(-0.491990626, -0.494420350, 0.575808167, 0.817584813),
+vec4f(-0.678437173, 0.0412081257, 0.534997642, 0.844853580),
+vec4f(-0.495611340, -0.807542920, 0.492898196, 0.870086968),
+vec4f(-0.508582473, 0.402262419, 0.449611306, 0.893224299),
+vec4f(-0.537693322, -0.333672792, 0.405241281, 0.914209783),
+vec4f(-0.614562631, -0.0568606183, 0.359894961, 0.932992816),
+vec4f(-0.452296555, 0.809743226, 0.313681662, 0.949528217),
+vec4f(-0.492342144, 0.317682207, 0.266712755, 0.963776052),
+vec4f(-0.403601944, -0.406031638, 0.219101220, 0.975702107),
+vec4f(-0.537260950, 0.137945175, 0.170961857, 0.985277653),
+vec4f(-0.404928178, -0.715918541, 0.122410625, 0.992479563),
+vec4f(-0.386338592, 0.556550741, 0.0735644922, 0.997290432),
+vec4f(-0.445942074, -0.242918074, 0.0245411359, 0.999698818),
+vec4f(-0.459759980, -0.175696611, -0.0245412234, 0.999698818),
+vec4f(-0.359771132, 0.682195961, -0.0735645741, 0.997290432),
+vec4f(-0.412575692, 0.205535084, -0.122410707, 0.992479503),
+vec4f(-0.336475104, -0.646421731, -0.170961946, 0.985277653),
+vec4f(-0.429363400, 0.0166855231, -0.219101310, 0.975702107),
+vec4f(-0.319253087, -0.925218344, -0.266712844, 0.963776052),
+vec4f(-0.294470191, 0.430103362, -0.313681722, 0.949528158),
+vec4f(-0.281283021, -0.387403429, -0.359895051, 0.932992816),
+vec4f(-0.358183086, -0.0808178857, -0.405241340, 0.914209723),
+vec4f(-0.271769434, 0.854052365, -0.449611366, 0.893224299),
+vec4f(-0.244775921, 0.311606824, -0.492898166, 0.870087028),
+vec4f(-0.244355842, -0.552090824, -0.534997702, 0.844853520),
+vec4f(-0.283291280, 0.112162955, -0.575808167, 0.817584813),
+vec4f(-0.224362880, -0.823741674, -0.615231693, 0.788346350),
+vec4f(-0.210826442, 0.610893846, -0.653172851, 0.757208824),
+vec4f(-0.191608176, -0.297364026, -0.689540684, 0.724246979),
+vec4f(-0.207833603, -0.124338157, -0.724247158, 0.689540505),
+vec4f(-0.176593989, 0.814833522, -0.757208824, 0.653172851),
+vec4f(-0.158949688, 0.293469012, -0.788346469, 0.615231514),
+vec4f(-0.152046368, -0.648668587, -0.817584813, 0.575808167),
+vec4f(-0.167782038, 0.0643178672, -0.844853640, 0.534997463),
+vec4f(-0.128430068, -0.907204390, -0.870087028, 0.492898136),
+vec4f(-0.115809105, 0.572147131, -0.893224299, 0.449611366),
+vec4f(-0.103222348, -0.403248310, -0.914209783, 0.405241221),
+vec4f(-0.112737618, -0.0319865458, -0.932992816, 0.359895051),
+vec4f(-0.0796749964, 0.955433667, -0.949528217, 0.313681602),
+vec4f(-0.0672519580, 0.453793734, -0.963776052, 0.266712725),
+vec4f(-0.0551275201, -0.538435221, -0.975702167, 0.219101071),
+vec4f(-0.0426490456, 0.204346836, -0.985277653, 0.170961812),
+vec4f(-0.0306719355, -0.790655255, -0.992479503, 0.122410700),
+vec4f(-0.0184057076, 0.708510995, -0.997290492, 0.0735644475),
+vec4f(-0.00613544276, -0.291185349, -0.999698818, 0.0245412104),
+vec4f(0.00613530306, -0.244297966, -0.999698818, -0.0245413855),
+vec4f(0.0184059385, 0.755400777, -0.997290432, -0.0735646188),
+vec4f(0.0306059923, 0.253786176, -0.992479503, -0.122410871),
+vec4f(0.0429276451, -0.743136168, -0.985277653, -0.170961991),
+vec4f(0.0701737553, 0.00441496586, -0.975702107, -0.219101235),
+vec4f(0.0674433559, -0.992085159, -0.963776052, -0.266712904),
+vec4f(0.0794365108, 0.499346048, -0.949528158, -0.313681781),
+vec4f(0.0915081054, -0.485832185, -0.932992756, -0.359895200),
+vec4f(0.101071075, -0.0861614645, -0.914209723, -0.405241400),
+vec4f(0.116242222, 0.872919321, -0.893224180, -0.449611515),
+vec4f(0.126407236, 0.359021723, -0.870086908, -0.492898285),
+vec4f(0.139908314, -0.603366435, -0.844853580, -0.534997642),
+vec4f(0.168981716, 0.0979395211, -0.817584693, -0.575808346),
+vec4f(0.164669126, -0.853637457, -0.788346410, -0.615231633),
+vec4f(0.175589874, 0.605686426, -0.757208705, -0.653172970),
+vec4f(0.181917489, -0.321471483, -0.724247038, -0.689540625),
+vec4f(0.219242185, -0.135647133, -0.689540565, -0.724247098),
+vec4f(0.212299466, 0.790099680, -0.653172910, -0.757208765),
+vec4f(0.208243251, 0.240495890, -0.615231574, -0.788346469),
+vec4f(0.234419063, -0.640313387, -0.575808048, -0.817584932),
+vec4f(0.315854102, 0.0532567836, -0.534997344, -0.844853759),
+vec4f(0.260333687, -0.894772291, -0.492898196, -0.870086968),
+vec4f(0.265374005, 0.502337217, -0.449611217, -0.893224359),
+vec4f(0.267429054, -0.339113116, -0.405241072, -0.914209843),
+vec4f(0.380209267, -0.0445682108, -0.359895110, -0.932992756),
+vec4f(0.307221651, 0.891683638, -0.313681692, -0.949528217),
+vec4f(0.296788067, 0.329054117, -0.266712576, -0.963776112),
+vec4f(0.317195147, -0.457708418, -0.219100907, -0.975702226),
+vec4f(0.419729054, 0.148764044, -0.170961887, -0.985277653),
+vec4f(0.349995196, -0.727014899, -0.122410536, -0.992479563),
+vec4f(0.356363803, 0.594497383, -0.0735642835, -0.997290492),
+vec4f(0.434345096, -0.231501341, -0.0245412868, -0.999698818),
+vec4f(0.480473697, -0.164373189, 0.0245413091, -0.999698818),
+vec4f(0.392657876, 0.681903303, 0.0735647827, -0.997290432),
+vec4f(0.492657542, 0.218807966, 0.122411035, -0.992479503),
+vec4f(0.409054011, -0.584142148, 0.170961916, -0.985277653),
+vec4f(0.569577277, 0.0289491210, 0.219101399, -0.975702107),
+vec4f(0.442882210, -0.855257273, 0.266713053, -0.963775992),
+vec4f(0.459734946, 0.387970626, 0.313681692, -0.949528217),
+vec4f(0.513070345, -0.343044043, 0.359895140, -0.932992756),
+vec4f(0.629031241, -0.0690753013, 0.405241549, -0.914209664),
+vec4f(0.483215511, 0.773316741, 0.449611664, -0.893224120),
+vec4f(0.586818159, 0.310843110, 0.492898226, -0.870086968),
+vec4f(0.528679490, -0.427168757, 0.534997761, -0.844853461),
+vec4f(0.683721542, 0.126429096, 0.575808465, -0.817584634),
+vec4f(0.517971814, -0.658907175, 0.615231574, -0.788346410),
+vec4f(0.548397303, 0.476606160, 0.653172910, -0.757208765),
+vec4f(0.695181608, -0.259932309, 0.689540744, -0.724246919),
+vec4f(0.748597383, -0.117820688, 0.724247038, -0.689540625),
+vec4f(0.556980729, 0.641256869, 0.757208884, -0.653172791),
+vec4f(0.741829991, 0.268901259, 0.788346529, -0.615231454),
+vec4f(0.647817433, -0.477341115, 0.817584991, -0.575807929),
+vec4f(0.816602290, 0.0779312477, 0.844853580, -0.534997642),
+vec4f(0.599919200, -0.671730816, 0.870087087, -0.492898077),
+vec4f(0.724726558, 0.447135508, 0.893224418, -0.449611068),
+vec4f(0.810727239, -0.307791352, 0.914209723, -0.405241370),
+vec4f(0.882583082, -0.0201240778, 0.932992816, -0.359894961),
+vec4f(0.645550549, 0.729843259, 0.949528277, -0.313681543),
+vec4f(0.839177251, 0.362342149, 0.963776171, -0.266712397),
+vec4f(0.839526892, -0.399391353, 0.975702107, -0.219101220),
+vec4f(0.928934336, 0.175204813, 0.985277653, -0.170961723),
+vec4f(0.774228692, -0.569184363, 0.992479563, -0.122410372),
+vec4f(0.815476179, 0.537282884, 0.997290432, -0.0735645965),
+vec4f(0.968665719, -0.214761987, 0.999698818, -0.0245411228)
+);
+`;
+//# sourceMappingURL=aovReflectanceSamples.wgsl.js.map

package/lib/renderer/wgsl/aovReflectanceSamples.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovReflectanceSamples.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovReflectanceSamples.wgsl.ts"],"names":[],"mappings":"AAAA,0EAA0E;AAC1E,qEAAqE;AACrE,eAAe,UAAU,CAAC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAmIzB,CAAC","sourcesContent":["// Generated by scripts/generate-aov-reflectance-samples.mjs; do not edit.\n// xy: concentric cosine disk; zw: spherical-cap azimuth cosine/sine.\nexport default /* wgsl */ `\nconst mtReflectanceProposals: array<vec4f, 128> = array<vec4f, 128>(\n  vec4f(-0.971247852, -0.202764958, 0.999698818, 0.0245412290),\n  vec4f(-0.808660269, 0.547487736, 0.997290432, 0.0735645667),\n  vec4f(-0.942509234, 0.187289074, 0.992479563, 0.122410677),\n  vec4f(-0.763038993, -0.558020830, 0.985277653, 0.170961902),\n  vec4f(-0.929654300, -0.00785388052, 0.975702107, 0.219101235),\n  vec4f(-0.656618595, 0.740912974, 0.963776052, 0.266712785),\n  vec4f(-0.817265987, 0.373183846, 0.949528158, 0.313681751),\n  vec4f(-0.793490529, -0.386950523, 0.932992816, 0.359895051),\n  vec4f(-0.860713601, -0.105764769, 0.914209723, 0.405241340),\n  vec4f(-0.604139447, 0.619064212, 0.893224299, 0.449611336),\n  vec4f(-0.787263036, 0.281084418, 0.870086968, 0.492898226),\n  vec4f(-0.673318744, -0.468566477, 0.844853580, 0.534997642),\n  vec4f(-0.799623907, 0.0901312605, 0.817584813, 0.575808227),\n  vec4f(-0.570307195, -0.676738262, 0.788346410, 0.615231633),\n  vec4f(-0.627450705, 0.452229142, 0.757208824, 0.653172851),\n  vec4f(-0.698281944, -0.294418186, 0.724247098, 0.689540565),\n  vec4f(-0.726036906, -0.153989345, 0.689540505, 0.724247098),\n  vec4f(-0.523753047, 0.608020604, 0.653172791, 0.757208884),\n  vec4f(-0.671607554, 0.233185455, 0.615231574, 0.788346410),\n  vec4f(-0.491990626, -0.494420350, 0.575808167, 0.817584813),\n  vec4f(-0.678437173, 0.0412081257, 0.534997642, 0.844853580),\n  vec4f(-0.495611340, -0.807542920, 0.492898196, 0.870086968),\n  vec4f(-0.508582473, 0.402262419, 0.449611306, 0.893224299),\n  vec4f(-0.537693322, -0.333672792, 0.405241281, 0.914209783),\n  vec4f(-0.614562631, -0.0568606183, 0.359894961, 0.932992816),\n  vec4f(-0.452296555, 0.809743226, 0.313681662, 0.949528217),\n  vec4f(-0.492342144, 0.317682207, 0.266712755, 0.963776052),\n  vec4f(-0.403601944, -0.406031638, 0.219101220, 0.975702107),\n  vec4f(-0.537260950, 0.137945175, 0.170961857, 0.985277653),\n  vec4f(-0.404928178, -0.715918541, 0.122410625, 0.992479563),\n  vec4f(-0.386338592, 0.556550741, 0.0735644922, 0.997290432),\n  vec4f(-0.445942074, -0.242918074, 0.0245411359, 0.999698818),\n  vec4f(-0.459759980, -0.175696611, -0.0245412234, 0.999698818),\n  vec4f(-0.359771132, 0.682195961, -0.0735645741, 0.997290432),\n  vec4f(-0.412575692, 0.205535084, -0.122410707, 0.992479503),\n  vec4f(-0.336475104, -0.646421731, -0.170961946, 0.985277653),\n  vec4f(-0.429363400, 0.0166855231, -0.219101310, 0.975702107),\n  vec4f(-0.319253087, -0.925218344, -0.266712844, 0.963776052),\n  vec4f(-0.294470191, 0.430103362, -0.313681722, 0.949528158),\n  vec4f(-0.281283021, -0.387403429, -0.359895051, 0.932992816),\n  vec4f(-0.358183086, -0.0808178857, -0.405241340, 0.914209723),\n  vec4f(-0.271769434, 0.854052365, -0.449611366, 0.893224299),\n  vec4f(-0.244775921, 0.311606824, -0.492898166, 0.870087028),\n  vec4f(-0.244355842, -0.552090824, -0.534997702, 0.844853520),\n  vec4f(-0.283291280, 0.112162955, -0.575808167, 0.817584813),\n  vec4f(-0.224362880, -0.823741674, -0.615231693, 0.788346350),\n  vec4f(-0.210826442, 0.610893846, -0.653172851, 0.757208824),\n  vec4f(-0.191608176, -0.297364026, -0.689540684, 0.724246979),\n  vec4f(-0.207833603, -0.124338157, -0.724247158, 0.689540505),\n  vec4f(-0.176593989, 0.814833522, -0.757208824, 0.653172851),\n  vec4f(-0.158949688, 0.293469012, -0.788346469, 0.615231514),\n  vec4f(-0.152046368, -0.648668587, -0.817584813, 0.575808167),\n  vec4f(-0.167782038, 0.0643178672, -0.844853640, 0.534997463),\n  vec4f(-0.128430068, -0.907204390, -0.870087028, 0.492898136),\n  vec4f(-0.115809105, 0.572147131, -0.893224299, 0.449611366),\n  vec4f(-0.103222348, -0.403248310, -0.914209783, 0.405241221),\n  vec4f(-0.112737618, -0.0319865458, -0.932992816, 0.359895051),\n  vec4f(-0.0796749964, 0.955433667, -0.949528217, 0.313681602),\n  vec4f(-0.0672519580, 0.453793734, -0.963776052, 0.266712725),\n  vec4f(-0.0551275201, -0.538435221, -0.975702167, 0.219101071),\n  vec4f(-0.0426490456, 0.204346836, -0.985277653, 0.170961812),\n  vec4f(-0.0306719355, -0.790655255, -0.992479503, 0.122410700),\n  vec4f(-0.0184057076, 0.708510995, -0.997290492, 0.0735644475),\n  vec4f(-0.00613544276, -0.291185349, -0.999698818, 0.0245412104),\n  vec4f(0.00613530306, -0.244297966, -0.999698818, -0.0245413855),\n  vec4f(0.0184059385, 0.755400777, -0.997290432, -0.0735646188),\n  vec4f(0.0306059923, 0.253786176, -0.992479503, -0.122410871),\n  vec4f(0.0429276451, -0.743136168, -0.985277653, -0.170961991),\n  vec4f(0.0701737553, 0.00441496586, -0.975702107, -0.219101235),\n  vec4f(0.0674433559, -0.992085159, -0.963776052, -0.266712904),\n  vec4f(0.0794365108, 0.499346048, -0.949528158, -0.313681781),\n  vec4f(0.0915081054, -0.485832185, -0.932992756, -0.359895200),\n  vec4f(0.101071075, -0.0861614645, -0.914209723, -0.405241400),\n  vec4f(0.116242222, 0.872919321, -0.893224180, -0.449611515),\n  vec4f(0.126407236, 0.359021723, -0.870086908, -0.492898285),\n  vec4f(0.139908314, -0.603366435, -0.844853580, -0.534997642),\n  vec4f(0.168981716, 0.0979395211, -0.817584693, -0.575808346),\n  vec4f(0.164669126, -0.853637457, -0.788346410, -0.615231633),\n  vec4f(0.175589874, 0.605686426, -0.757208705, -0.653172970),\n  vec4f(0.181917489, -0.321471483, -0.724247038, -0.689540625),\n  vec4f(0.219242185, -0.135647133, -0.689540565, -0.724247098),\n  vec4f(0.212299466, 0.790099680, -0.653172910, -0.757208765),\n  vec4f(0.208243251, 0.240495890, -0.615231574, -0.788346469),\n  vec4f(0.234419063, -0.640313387, -0.575808048, -0.817584932),\n  vec4f(0.315854102, 0.0532567836, -0.534997344, -0.844853759),\n  vec4f(0.260333687, -0.894772291, -0.492898196, -0.870086968),\n  vec4f(0.265374005, 0.502337217, -0.449611217, -0.893224359),\n  vec4f(0.267429054, -0.339113116, -0.405241072, -0.914209843),\n  vec4f(0.380209267, -0.0445682108, -0.359895110, -0.932992756),\n  vec4f(0.307221651, 0.891683638, -0.313681692, -0.949528217),\n  vec4f(0.296788067, 0.329054117, -0.266712576, -0.963776112),\n  vec4f(0.317195147, -0.457708418, -0.219100907, -0.975702226),\n  vec4f(0.419729054, 0.148764044, -0.170961887, -0.985277653),\n  vec4f(0.349995196, -0.727014899, -0.122410536, -0.992479563),\n  vec4f(0.356363803, 0.594497383, -0.0735642835, -0.997290492),\n  vec4f(0.434345096, -0.231501341, -0.0245412868, -0.999698818),\n  vec4f(0.480473697, -0.164373189, 0.0245413091, -0.999698818),\n  vec4f(0.392657876, 0.681903303, 0.0735647827, -0.997290432),\n  vec4f(0.492657542, 0.218807966, 0.122411035, -0.992479503),\n  vec4f(0.409054011, -0.584142148, 0.170961916, -0.985277653),\n  vec4f(0.569577277, 0.0289491210, 0.219101399, -0.975702107),\n  vec4f(0.442882210, -0.855257273, 0.266713053, -0.963775992),\n  vec4f(0.459734946, 0.387970626, 0.313681692, -0.949528217),\n  vec4f(0.513070345, -0.343044043, 0.359895140, -0.932992756),\n  vec4f(0.629031241, -0.0690753013, 0.405241549, -0.914209664),\n  vec4f(0.483215511, 0.773316741, 0.449611664, -0.893224120),\n  vec4f(0.586818159, 0.310843110, 0.492898226, -0.870086968),\n  vec4f(0.528679490, -0.427168757, 0.534997761, -0.844853461),\n  vec4f(0.683721542, 0.126429096, 0.575808465, -0.817584634),\n  vec4f(0.517971814, -0.658907175, 0.615231574, -0.788346410),\n  vec4f(0.548397303, 0.476606160, 0.653172910, -0.757208765),\n  vec4f(0.695181608, -0.259932309, 0.689540744, -0.724246919),\n  vec4f(0.748597383, -0.117820688, 0.724247038, -0.689540625),\n  vec4f(0.556980729, 0.641256869, 0.757208884, -0.653172791),\n  vec4f(0.741829991, 0.268901259, 0.788346529, -0.615231454),\n  vec4f(0.647817433, -0.477341115, 0.817584991, -0.575807929),\n  vec4f(0.816602290, 0.0779312477, 0.844853580, -0.534997642),\n  vec4f(0.599919200, -0.671730816, 0.870087087, -0.492898077),\n  vec4f(0.724726558, 0.447135508, 0.893224418, -0.449611068),\n  vec4f(0.810727239, -0.307791352, 0.914209723, -0.405241370),\n  vec4f(0.882583082, -0.0201240778, 0.932992816, -0.359894961),\n  vec4f(0.645550549, 0.729843259, 0.949528277, -0.313681543),\n  vec4f(0.839177251, 0.362342149, 0.963776171, -0.266712397),\n  vec4f(0.839526892, -0.399391353, 0.975702107, -0.219101220),\n  vec4f(0.928934336, 0.175204813, 0.985277653, -0.170961723),\n  vec4f(0.774228692, -0.569184363, 0.992479563, -0.122410372),\n  vec4f(0.815476179, 0.537282884, 0.997290432, -0.0735645965),\n  vec4f(0.968665719, -0.214761987, 0.999698818, -0.0245411228)\n);\n`;\n"]}

package/lib/renderer/wgsl/aovResolve.wgsl.d.ts

@@ -0,0 +1,5 @@
+/** Shared normalization for readback and GPU denoiser snapshots. */
+export declare const aovResolveCommon = "\nstruct AovResolvedAccumulation { geometry: vec4f, position: vec4f, diffuse: vec4f, specular: vec4f, }\nstruct AovResolvedCenter { identity: vec4u, motion: vec4f, }\n@group(0) @binding(0) var<storage,read> accumulation: array<AovResolvedAccumulation>;\n@group(0) @binding(1) var<storage,read> center: array<AovResolvedCenter>;\n@group(0) @binding(3) var<uniform> parameters: vec4u;\nfn aovResolvedValue(source: u32, kind: u32) -> vec4u {\n  if(kind==5u){return center[source].identity;}\n  var value=vec4f(0.0);\n  if(kind==6u){let motion=center[source].motion;value=vec4f(motion.xy,0.0,motion.w);}\n  else{\n    let pixel=accumulation[source];let hits=pixel.position.w;\n    let scale=select(0.0,1.0/max(hits,1.0),hits>0.0);let coverage=hits/f32(max(parameters.z,1u));\n    if(kind==0u){value=vec4f(pixel.geometry.x*scale,0.0,0.0,coverage);}\n    else if(kind==1u){value=vec4f(pixel.position.xyz*scale,coverage);}\n    else if(kind==2u){value=vec4f(pixel.geometry.yzw*scale,coverage);}\n    else{let reflectance=select(pixel.diffuse,pixel.specular,kind==4u);\n      let inv=select(0.0,1.0/max(reflectance.w,1.0),reflectance.w>0.0);\n      value=vec4f(reflectance.xyz*inv,reflectance.w/f32(max(parameters.z,1u)));\n    }\n  }\n  return bitcast<vec4u>(value);\n}\n";
+/** Bit 8 of kind requests top row first. Kept for the native single-output pass. */
+declare const _default: string;
+export default _default;

package/lib/renderer/wgsl/aovResolve.wgsl.js

@@ -0,0 +1,38 @@
+/** Shared normalization for readback and GPU denoiser snapshots. */
+export const aovResolveCommon = /* wgsl */ `
+struct AovResolvedAccumulation { geometry: vec4f, position: vec4f, diffuse: vec4f, specular: vec4f, }
+struct AovResolvedCenter { identity: vec4u, motion: vec4f, }
+@group(0) @binding(0) var<storage,read> accumulation: array<AovResolvedAccumulation>;
+@group(0) @binding(1) var<storage,read> center: array<AovResolvedCenter>;
+@group(0) @binding(3) var<uniform> parameters: vec4u;
+fn oj(source: u32, kind: u32) -> vec4u {
+if(kind==5u){return center[source].identity;}
+var value=vec4f(0.0);
+if(kind==6u){let motion=center[source].motion;value=vec4f(motion.xy,0.0,motion.w);}
+else{
+let pixel=accumulation[source];let hits=pixel.position.w;
+let scale=select(0.0,1.0/max(hits,1.0),hits>0.0);let coverage=hits/f32(max(parameters.z,1u));
+if(kind==0u){value=vec4f(pixel.geometry.x*scale,0.0,0.0,coverage);}
+else if(kind==1u){value=vec4f(pixel.position.xyz*scale,coverage);}
+else if(kind==2u){value=vec4f(pixel.geometry.yzw*scale,coverage);}
+else{let reflectance=select(pixel.diffuse,pixel.specular,kind==4u);
+let inv=select(0.0,1.0/max(reflectance.w,1.0),reflectance.w>0.0);
+value=vec4f(reflectance.xyz*inv,reflectance.w/f32(max(parameters.z,1u)));
+}
+}
+return bitcast<vec4u>(value);
+}
+`;
+/** Bit 8 of kind requests top row first. Kept for the native single-output pass. */
+export default aovResolveCommon + /* wgsl */ `
+@group(0) @binding(2) var<storage,read_write> output: array<vec4u>;
+@compute @workgroup_size(8,8)
+fn main(@builtin(global_invocation_id) id: vec3u) {
+let width=parameters.x;let height=parameters.y;
+if(id.x>=width||id.y>=height){return;}
+let index=id.x+id.y*width;
+let source=id.x+select(id.y,height-1u-id.y,(parameters.w&256u)!=0u)*width;
+output[index]=oj(source,parameters.w&255u);
+}
+`;
+//# sourceMappingURL=aovResolve.wgsl.js.map

package/lib/renderer/wgsl/aovResolve.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovResolve.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovResolve.wgsl.ts"],"names":[],"mappings":"AAAA,oEAAoE;AACpE,MAAM,CAAC,MAAM,gBAAgB,GAAG,UAAU,CAAC;;;;;;;;;;;;;;;;;;;;;;;CAuB1C,CAAC;AACF,oFAAoF;AACpF,eAAe,gBAAgB,GAAG,UAAU,CAAC;;;;;;;;;;CAU5C,CAAC","sourcesContent":["/** Shared normalization for readback and GPU denoiser snapshots. */\nexport const aovResolveCommon = /* wgsl */ `\nstruct AovResolvedAccumulation { geometry: vec4f, position: vec4f, diffuse: vec4f, specular: vec4f, }\nstruct AovResolvedCenter { identity: vec4u, motion: vec4f, }\n@group(0) @binding(0) var<storage,read> accumulation: array<AovResolvedAccumulation>;\n@group(0) @binding(1) var<storage,read> center: array<AovResolvedCenter>;\n@group(0) @binding(3) var<uniform> parameters: vec4u;\nfn aovResolvedValue(source: u32, kind: u32) -> vec4u {\n  if(kind==5u){return center[source].identity;}\n  var value=vec4f(0.0);\n  if(kind==6u){let motion=center[source].motion;value=vec4f(motion.xy,0.0,motion.w);}\n  else{\n    let pixel=accumulation[source];let hits=pixel.position.w;\n    let scale=select(0.0,1.0/max(hits,1.0),hits>0.0);let coverage=hits/f32(max(parameters.z,1u));\n    if(kind==0u){value=vec4f(pixel.geometry.x*scale,0.0,0.0,coverage);}\n    else if(kind==1u){value=vec4f(pixel.position.xyz*scale,coverage);}\n    else if(kind==2u){value=vec4f(pixel.geometry.yzw*scale,coverage);}\n    else{let reflectance=select(pixel.diffuse,pixel.specular,kind==4u);\n      let inv=select(0.0,1.0/max(reflectance.w,1.0),reflectance.w>0.0);\n      value=vec4f(reflectance.xyz*inv,reflectance.w/f32(max(parameters.z,1u)));\n    }\n  }\n  return bitcast<vec4u>(value);\n}\n`;\n/** Bit 8 of kind requests top row first. Kept for the native single-output pass. */\nexport default aovResolveCommon + /* wgsl */ `\n@group(0) @binding(2) var<storage,read_write> output: array<vec4u>;\n@compute @workgroup_size(8,8)\nfn main(@builtin(global_invocation_id) id: vec3u) {\n  let width=parameters.x;let height=parameters.y;\n  if(id.x>=width||id.y>=height){return;}\n  let index=id.x+id.y*width;\n  let source=id.x+select(id.y,height-1u-id.y,(parameters.w&256u)!=0u)*width;\n  output[index]=aovResolvedValue(source,parameters.w&255u);\n}\n`;\n"]}

package/lib/renderer/wgsl/aovSnapshot.wgsl.d.ts

@@ -0,0 +1,4 @@
+import { type AovKind } from '../aov.js';
+export declare const AOV_SNAPSHOT_BATCH_SIZE = 6;
+export declare const aovSnapshotBinding: (index: number) => number;
+export declare function aovSnapshotShader(kinds: readonly AovKind[]): string;

package/lib/renderer/wgsl/aovSnapshot.wgsl.js

@@ -0,0 +1,21 @@
+import { AOV_KINDS } from '../aov.js';
+import { aovResolveCommon } from './aovResolve.wgsl.js';
+// Two input storage buffers leave six output slots on the baseline WebGPU
+// eight-storage-buffer limit. Binding 3 remains the shared uniform binding.
+export const AOV_SNAPSHOT_BATCH_SIZE = 6;
+export const aovSnapshotBinding = (index) => index === 0 ? 2 : index + 3;
+export function aovSnapshotShader(kinds) {
+    if (!kinds.length || kinds.length > AOV_SNAPSHOT_BATCH_SIZE)
+        throw new RangeError('Invalid AOV snapshot batch size');
+    return aovResolveCommon + kinds.map((_, i) => `@group(0) @binding(${aovSnapshotBinding(i)}) var<storage,read_write> output${i}: array<vec4u>;`).join('\n') + /* wgsl */ `
+@compute @workgroup_size(8,8)
+fn main(@builtin(global_invocation_id) id: vec3u) {
+let width=parameters.x;let height=parameters.y;
+if(id.x>=width||id.y>=height){return;}
+let index=id.x+id.y*width;
+let source=id.x+(height-1u-id.y)*width;
+${kinds.map((kind, i) => `output${i}[index]=oj(source,${AOV_KINDS.indexOf(kind)}u);`).join('\n  ')}
+}
+`;
+}
+//# sourceMappingURL=aovSnapshot.wgsl.js.map

package/lib/renderer/wgsl/aovSnapshot.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovSnapshot.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovSnapshot.wgsl.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,SAAS,EAAe,MAAM,WAAW,CAAC;AAClD,OAAO,EAAC,gBAAgB,EAAC,MAAM,sBAAsB,CAAC;AAEtD,0EAA0E;AAC1E,4EAA4E;AAC5E,MAAM,CAAC,MAAM,uBAAuB,GAAG,CAAC,CAAC;AACzC,MAAM,CAAC,MAAM,kBAAkB,GAAG,CAAC,KAAa,EAAE,EAAE,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;AAEjF,MAAM,UAAU,iBAAiB,CAAC,KAAyB;IACzD,IAAI,CAAC,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM,GAAG,uBAAuB;QAAE,MAAM,IAAI,UAAU,CAAC,iCAAiC,CAAC,CAAC;IACrH,OAAO,gBAAgB,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAC3C,sBAAsB,kBAAkB,CAAC,CAAC,CAAC,mCAAmC,CAAC,iBAAiB,CACjG,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,UAAU,CAAC;;;;;;;EAOxB,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC,SAAS,CAAC,qBAAmC,SAAS,CAAC,OAAO,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,CAAC,MAAM,CAAC;;CAEjH,CAAC;AACF,CAAC","sourcesContent":["import {AOV_KINDS, type AovKind} from '../aov.js';\nimport {aovResolveCommon} from './aovResolve.wgsl.js';\n\n// Two input storage buffers leave six output slots on the baseline WebGPU\n// eight-storage-buffer limit. Binding 3 remains the shared uniform binding.\nexport const AOV_SNAPSHOT_BATCH_SIZE = 6;\nexport const aovSnapshotBinding = (index: number) => index === 0 ? 2 : index + 3;\n\nexport function aovSnapshotShader(kinds: readonly AovKind[]): string {\n  if (!kinds.length || kinds.length > AOV_SNAPSHOT_BATCH_SIZE) throw new RangeError('Invalid AOV snapshot batch size');\n  return aovResolveCommon + kinds.map((_, i) =>\n    `@group(0) @binding(${aovSnapshotBinding(i)}) var<storage,read_write> output${i}: array<vec4u>;`\n  ).join('\\n') + /* wgsl */ `\n@compute @workgroup_size(8,8)\nfn main(@builtin(global_invocation_id) id: vec3u) {\n  let width=parameters.x;let height=parameters.y;\n  if(id.x>=width||id.y>=height){return;}\n  let index=id.x+id.y*width;\n  let source=id.x+(height-1u-id.y)*width;\n  ${kinds.map((kind, i) => `output${i}[index]=aovResolvedValue(source,${AOV_KINDS.indexOf(kind)}u);`).join('\\n  ')}\n}\n`;\n}\n"]}

package/lib/renderer/wgsl/aovTrace.wgsl.d.ts

@@ -0,0 +1,2 @@
+declare const _default: "\noverride AOV_REFLECTANCE: bool = true;\nstruct AovAccumulation {\n  depthRoughness: vec4f,\n  position: vec4f,\n  diffuse: vec4f,\n  specular: vec4f,\n}\nstruct AovCenter { identity: vec4u, motion: vec4f, }\n@group(3) @binding(5) var<storage, read_write> aovAccum: array<AovAccumulation>;\n@group(3) @binding(6) var<storage, read_write> aovCenter: array<AovCenter>;\n@group(3) @binding(7) var<storage, read> aovScene: array<u32>;\n@group(3) @binding(9) var<storage, read> aovPrevious: array<AovCenter>;\n@group(3) @binding(8) var<uniform> aovParameters: vec4u;\n\nfn aovMatrix(offset: u32) -> mat4x4f {\n  return mat4x4f(bitcast<vec4f>(vec4u(aovScene[offset], aovScene[offset+1u], aovScene[offset+2u], aovScene[offset+3u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+4u], aovScene[offset+5u], aovScene[offset+6u], aovScene[offset+7u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+8u], aovScene[offset+9u], aovScene[offset+10u], aovScene[offset+11u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+12u], aovScene[offset+13u], aovScene[offset+14u], aovScene[offset+15u])));\n}\nfn aovVector(offset:u32)->vec3f {return bitcast<vec3f>(vec3u(aovScene[offset],aovScene[offset+1u],aovScene[offset+2u]));}\nfn aovSurface(rayInput: Ray, si: ptr<function, SurfaceInteraction>, material: ptr<function, OpenPbrResolvedMaterial>) -> bool {\n  var ray = rayInput;\n  (*si).triangleIndex = -1;\n  for (var layer = 0; layer < 32; layer++) {\n    intersectScene(ray, si, 0, false, false);\n    if (!(*si).hit) { return false; }\n    getMaterial(material, (*si).materialIndex, (*si).uv, (*si).uv1, (*si).instanceIndex, (*si).vertexColor);\n    if ((*material).materialType != SHADOW_CATCHER && (*material).alpha >= 0.5) { return true; }\n    ray.o = (*si).position + ray.d * (EPS * 0.02);\n  }\n  return false;\n}\n@compute @workgroup_size(8, 8)\nfn aovTrace(@builtin(global_invocation_id) id: vec3u) {\n  let size = aovParameters.xy;\n  if (any(id.xy >= size)) { return; }\n  let index = id.x + id.y * size.x;\n  let sample = aovParameters.z;\n  if (sample == 0u) { aovAccum[index] = AovAccumulation(); }\n  samplingPixel = id.xy; pixelSampleIndex = sample + rtUniform.sampleSequenceOffset;\n  samplingBounce = 0u; hashedCounter = 0u;\n  rng_initialize(id.xy, i32(pixelSampleIndex ^ sampleHash(rtUniform.sequenceSeed)));\n  let uv = (vec2f(id.xy) + sample2D(CAMERA_PIXEL)) / vec2f(size);\n  var ray = initCameraRay(uv);\n  if (cameraUniform.apertureRadius > 0.0) {\n    let eye = cameraUniform.viewInv[3].xyz;\n    let forward = -normalize(cameraUniform.viewInv[2].xyz);\n    let focusPoint = eye + ray.d * (cameraUniform.focusDistance / dot(ray.d, forward));\n    let lens = cameraUniform.apertureRadius * sampleCircle(sample2D(CAMERA_LENS));\n    let origin = (cameraUniform.viewInv * vec4f(lens, 0.0, 1.0)).xyz;\n    initRay(&ray, origin, normalize(focusPoint-origin));\n  }\n  var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;\n  if (!aovSurface(ray, &si, &material)) { return; }\n  let depth = max(dot(si.position-cameraUniform.viewInv[3].xyz, -normalize(cameraUniform.viewInv[2].xyz)), 0.0);\n  let unsupported = (material.materialType != STANDARD) | (material.alpha < 1.0) |\n      (material.coatWeight > 0.0) | (material.fuzzWeight > 0.0) |\n      (material.specularAnisotropy > 0.0) | (material.subsurfaceWeight > 0.0) |\n      (dot(material.mediumScattering, material.mediumScattering) > 0.0) |\n      (material.bevel > 0.0) | (material.thinFilmWeight > 0.0) |\n      (material.transmissionDispersionScale > 0.0);\n  aovAccum[index].depthRoughness += vec4f(depth, material.specularRoughness, material.baseMetalness, select(0.0,1.0,unsupported));\n  if (AOV_REFLECTANCE && (aovParameters.w & 24u) != 0u && !unsupported) {\n    var normal=si.normal; var face=si.faceNormal;\n    var cosine=-dot(normal,ray.d); var faceCosine=-dot(face,ray.d);\n    if (material.thinWalled && (material.transmissionWeight>0.0 || material.subsurfaceWeight>0.0)) {\n      if(cosine<0.0){normal=-normal;cosine=-cosine;}\n      if(faceCosine<0.0){face=-face;faceCosine=-faceCosine;}\n    }\n    if(faceCosine>=0.0 && cosine<0.0){normal=face;cosine=faceCosine;}\n    normal=select(-normal,normal,cosine>=0.0);\n    material.inside=(cosine<=0.0) && !material.thinWalled;\n    material.relativeIor=select(material.specularIor,1.0/material.specularIor,cosine>0.0);\n    let prepared=openpbr_prepare(si.tangent,material,-ray.d,normal);\n    let reflectance=mtFrameReflectanceComponents(prepared, (aovParameters.w & 8u) != 0u, (aovParameters.w & 16u) != 0u);\n    aovAccum[index].diffuse+=vec4f(reflectance.diffuse,1.0);\n    aovAccum[index].specular+=vec4f(reflectance.specular,1.0);\n  }\n  aovAccum[index].position += vec4f(si.position, 1.0);\n}\n@compute @workgroup_size(8, 8)\nfn aovTraceCenter(@builtin(global_invocation_id) id: vec3u) {\n  let size = aovParameters.xy;\n  if (any(id.xy >= size)) { return; }\n  let index = id.x + id.y * size.x;\n  aovCenter[index] = AovCenter(vec4u(0xffffffffu), vec4f(0.0));\n  let uv = (vec2f(id.xy) + 0.5) / vec2f(size);\n  let ray = initCameraRay(uv);\n  var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;\n  if (!aovSurface(ray, &si, &material)) { return; }\n  let offset = 52u + u32(si.instanceIndex) * 24u;\n  let materialOffset = aovScene[16u];\n  let primitiveOffset = aovScene[17u];\n  aovCenter[index].identity = vec4u(aovScene[offset], aovScene[offset+1u],\n    aovScene[materialOffset + u32(si.materialIndex)], aovScene[primitiveOffset + u32(si.triangleIndex)]);\n  aovCenter[index].motion.z = dot(si.position-aovVector(36u),aovVector(40u));\n  if (aovScene[offset+2u] == 0u) { return; }\n  let oldWorld = aovMatrix(offset+4u) * vec4f(si.position,1.0);\n  let oldClip = aovMatrix(0u) * oldWorld;\n  if (oldClip.w <= 0.0) { return; }\n  let oldNdc=oldClip.xyz/oldClip.w;\n  let oldPixel=(oldNdc.xy*0.5+0.5)*vec2f(size);\n  if (any(oldPixel < vec2f(0.0)) || any(oldPixel >= vec2f(size)) || abs(oldNdc.z)>1.0) { return; }\n  let previous=aovPrevious[u32(oldPixel.x)+u32(oldPixel.y)*size.x];\n  let oldDepth=dot(oldWorld.xyz-aovVector(44u),aovVector(48u));\n  if (any(previous.identity.xy != aovCenter[index].identity.xy) || abs(previous.motion.z-oldDepth)>max(1e-4,abs(oldDepth)*1e-3)) { return; }\n  let motion=(oldPixel-(vec2f(id.xy)+0.5))*vec2f(1.0,-1.0);\n  aovCenter[index].motion=vec4f(motion,aovCenter[index].motion.z,3.0);\n}\n";
+export default _default;

package/lib/renderer/wgsl/aovTrace.wgsl.js

@@ -0,0 +1,115 @@
+// These entry points use their own output bindings and never mutate beauty RNG
+// state or accumulation. Storage usage remains four scene + three AOV buffers.
+export default /* wgsl */ `
+override AOV_REFLECTANCE: bool = true;
+struct AovAccumulation {
+depthRoughness: vec4f,
+position: vec4f,
+diffuse: vec4f,
+specular: vec4f,
+}
+struct AovCenter { identity: vec4u, motion: vec4f, }
+@group(3) @binding(5) var<storage, read_write> aovAccum: array<AovAccumulation>;
+@group(3) @binding(6) var<storage, read_write> aovCenter: array<AovCenter>;
+@group(3) @binding(7) var<storage, read> aovScene: array<u32>;
+@group(3) @binding(9) var<storage, read> aovPrevious: array<AovCenter>;
+@group(3) @binding(8) var<uniform> aovParameters: vec4u;
+
+fn oi(offset: u32) -> mat4x4f {
+return mat4x4f(bitcast<vec4f>(vec4u(aovScene[offset], aovScene[offset+1u], aovScene[offset+2u], aovScene[offset+3u])),
+bitcast<vec4f>(vec4u(aovScene[offset+4u], aovScene[offset+5u], aovScene[offset+6u], aovScene[offset+7u])),
+bitcast<vec4f>(vec4u(aovScene[offset+8u], aovScene[offset+9u], aovScene[offset+10u], aovScene[offset+11u])),
+bitcast<vec4f>(vec4u(aovScene[offset+12u], aovScene[offset+13u], aovScene[offset+14u], aovScene[offset+15u])));
+}
+fn iu(offset:u32)->vec3f {return bitcast<vec3f>(vec3u(aovScene[offset],aovScene[offset+1u],aovScene[offset+2u]));}
+fn mC(xG: Ray, si: ptr<function, SurfaceInteraction>, material: ptr<function, OpenPbrResolvedMaterial>) -> bool {
+var ray = xG;
+(*si).triangleIndex = -1;
+for (var layer = 0; layer < 32; layer++) {
+intersectScene(ray, si, 0, false, false);
+if (!(*si).hit) { return false; }
+getMaterial(material, (*si).materialIndex, (*si).uv, (*si).uv1, (*si).instanceIndex, (*si).vertexColor);
+if ((*material).materialType != SHADOW_CATCHER && (*material).alpha >= 0.5) { return true; }
+ray.o = (*si).position + ray.d * (EPS * 0.02);
+}
+return false;
+}
+@compute @workgroup_size(8, 8)
+fn aovTrace(@builtin(global_invocation_id) id: vec3u) {
+let size = aovParameters.xy;
+if (any(id.xy >= size)) { return; }
+let index = id.x + id.y * size.x;
+let sample = aovParameters.z;
+if (sample == 0u) { aovAccum[index] = AovAccumulation(); }
+aI = id.xy; aC = sample + rtUniform.sampleSequenceOffset;
+samplingBounce = 0u; ca = 0u;
+fi(id.xy, i32(aC ^ cB(rtUniform.sequenceSeed)));
+let uv = (vec2f(id.xy) + sample2D(CAMERA_PIXEL)) / vec2f(size);
+var ray = eT(uv);
+if (cameraUniform.apertureRadius > 0.0) {
+let eye = cameraUniform.viewInv[3].xyz;
+let forward = -normalize(cameraUniform.viewInv[2].xyz);
+let ju = eye + ray.d * (cameraUniform.focusDistance / dot(ray.d, forward));
+let uV = cameraUniform.apertureRadius * kq(sample2D(CAMERA_LENS));
+let origin = (cameraUniform.viewInv * vec4f(uV, 0.0, 1.0)).xyz;
+initRay(&ray, origin, normalize(ju-origin));
+}
+var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;
+if (!mC(ray, &si, &material)) { return; }
+let depth = max(dot(si.position-cameraUniform.viewInv[3].xyz, -normalize(cameraUniform.viewInv[2].xyz)), 0.0);
+let unsupported = (material.materialType != STANDARD) | (material.alpha < 1.0) |
+(material.coatWeight > 0.0) | (material.fuzzWeight > 0.0) |
+(material.specularAnisotropy > 0.0) | (material.subsurfaceWeight > 0.0) |
+(dot(material.mediumScattering, material.mediumScattering) > 0.0) |
+(material.bevel > 0.0) | (material.thinFilmWeight > 0.0) |
+(material.transmissionDispersionScale > 0.0);
+aovAccum[index].depthRoughness += vec4f(depth, material.specularRoughness, material.baseMetalness, select(0.0,1.0,unsupported));
+if (AOV_REFLECTANCE && (aovParameters.w & 24u) != 0u && !unsupported) {
+var normal=si.normal; var face=si.faceNormal;
+var cosine=-dot(normal,ray.d); var fa=-dot(face,ray.d);
+if (material.thinWalled && (material.transmissionWeight>0.0 || material.subsurfaceWeight>0.0)) {
+if(cosine<0.0){normal=-normal;cosine=-cosine;}
+if(fa<0.0){face=-face;fa=-fa;}
+}
+if(fa>=0.0 && cosine<0.0){normal=face;cosine=fa;}
+normal=select(-normal,normal,cosine>=0.0);
+material.inside=(cosine<=0.0) && !material.thinWalled;
+material.relativeIor=select(material.specularIor,1.0/material.specularIor,cosine>0.0);
+let prepared=bV(si.tangent,material,-ray.d,normal);
+let reflectance=hu(prepared, (aovParameters.w & 8u) != 0u, (aovParameters.w & 16u) != 0u);
+aovAccum[index].diffuse+=vec4f(reflectance.diffuse,1.0);
+aovAccum[index].specular+=vec4f(reflectance.specular,1.0);
+}
+aovAccum[index].position += vec4f(si.position, 1.0);
+}
+@compute @workgroup_size(8, 8)
+fn aovTraceCenter(@builtin(global_invocation_id) id: vec3u) {
+let size = aovParameters.xy;
+if (any(id.xy >= size)) { return; }
+let index = id.x + id.y * size.x;
+aovCenter[index] = AovCenter(vec4u(0xffffffffu), vec4f(0.0));
+let uv = (vec2f(id.xy) + 0.5) / vec2f(size);
+let ray = eT(uv);
+var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;
+if (!mC(ray, &si, &material)) { return; }
+let offset = 52u + u32(si.instanceIndex) * 24u;
+let materialOffset = aovScene[16u];
+let primitiveOffset = aovScene[17u];
+aovCenter[index].identity = vec4u(aovScene[offset], aovScene[offset+1u],
+aovScene[materialOffset + u32(si.materialIndex)], aovScene[primitiveOffset + u32(si.triangleIndex)]);
+aovCenter[index].motion.z = dot(si.position-iu(36u),iu(40u));
+if (aovScene[offset+2u] == 0u) { return; }
+let qs = oi(offset+4u) * vec4f(si.position,1.0);
+let pg = oi(0u) * qs;
+if (pg.w <= 0.0) { return; }
+let ua=pg.xyz/pg.w;
+let hU=(ua.xy*0.5+0.5)*vec2f(size);
+if (any(hU < vec2f(0.0)) || any(hU >= vec2f(size)) || abs(ua.z)>1.0) { return; }
+let previous=aovPrevious[u32(hU.x)+u32(hU.y)*size.x];
+let qr=dot(qs.xyz-iu(44u),iu(48u));
+if (any(previous.identity.xy != aovCenter[index].identity.xy) || abs(previous.motion.z-qr)>max(1e-4,abs(qr)*1e-3)) { return; }
+let motion=(hU-(vec2f(id.xy)+0.5))*vec2f(1.0,-1.0);
+aovCenter[index].motion=vec4f(motion,aovCenter[index].motion.z,3.0);
+}
+`;
+//# sourceMappingURL=aovTrace.wgsl.js.map

package/lib/renderer/wgsl/aovTrace.wgsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"aovTrace.wgsl.js","sourceRoot":"","sources":["../../../src/renderer/wgsl/aovTrace.wgsl.ts"],"names":[],"mappings":"AAAA,+EAA+E;AAC/E,+EAA+E;AAC/E,eAAe,UAAU,CAAC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+GzB,CAAC","sourcesContent":["// These entry points use their own output bindings and never mutate beauty RNG\n// state or accumulation. Storage usage remains four scene + three AOV buffers.\nexport default /* wgsl */ `\noverride AOV_REFLECTANCE: bool = true;\nstruct AovAccumulation {\n  depthRoughness: vec4f,\n  position: vec4f,\n  diffuse: vec4f,\n  specular: vec4f,\n}\nstruct AovCenter { identity: vec4u, motion: vec4f, }\n@group(3) @binding(5) var<storage, read_write> aovAccum: array<AovAccumulation>;\n@group(3) @binding(6) var<storage, read_write> aovCenter: array<AovCenter>;\n@group(3) @binding(7) var<storage, read> aovScene: array<u32>;\n@group(3) @binding(9) var<storage, read> aovPrevious: array<AovCenter>;\n@group(3) @binding(8) var<uniform> aovParameters: vec4u;\n\nfn aovMatrix(offset: u32) -> mat4x4f {\n  return mat4x4f(bitcast<vec4f>(vec4u(aovScene[offset], aovScene[offset+1u], aovScene[offset+2u], aovScene[offset+3u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+4u], aovScene[offset+5u], aovScene[offset+6u], aovScene[offset+7u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+8u], aovScene[offset+9u], aovScene[offset+10u], aovScene[offset+11u])),\n    bitcast<vec4f>(vec4u(aovScene[offset+12u], aovScene[offset+13u], aovScene[offset+14u], aovScene[offset+15u])));\n}\nfn aovVector(offset:u32)->vec3f {return bitcast<vec3f>(vec3u(aovScene[offset],aovScene[offset+1u],aovScene[offset+2u]));}\nfn aovSurface(rayInput: Ray, si: ptr<function, SurfaceInteraction>, material: ptr<function, OpenPbrResolvedMaterial>) -> bool {\n  var ray = rayInput;\n  (*si).triangleIndex = -1;\n  for (var layer = 0; layer < 32; layer++) {\n    intersectScene(ray, si, 0, false, false);\n    if (!(*si).hit) { return false; }\n    getMaterial(material, (*si).materialIndex, (*si).uv, (*si).uv1, (*si).instanceIndex, (*si).vertexColor);\n    if ((*material).materialType != SHADOW_CATCHER && (*material).alpha >= 0.5) { return true; }\n    ray.o = (*si).position + ray.d * (EPS * 0.02);\n  }\n  return false;\n}\n@compute @workgroup_size(8, 8)\nfn aovTrace(@builtin(global_invocation_id) id: vec3u) {\n  let size = aovParameters.xy;\n  if (any(id.xy >= size)) { return; }\n  let index = id.x + id.y * size.x;\n  let sample = aovParameters.z;\n  if (sample == 0u) { aovAccum[index] = AovAccumulation(); }\n  samplingPixel = id.xy; pixelSampleIndex = sample + rtUniform.sampleSequenceOffset;\n  samplingBounce = 0u; hashedCounter = 0u;\n  rng_initialize(id.xy, i32(pixelSampleIndex ^ sampleHash(rtUniform.sequenceSeed)));\n  let uv = (vec2f(id.xy) + sample2D(CAMERA_PIXEL)) / vec2f(size);\n  var ray = initCameraRay(uv);\n  if (cameraUniform.apertureRadius > 0.0) {\n    let eye = cameraUniform.viewInv[3].xyz;\n    let forward = -normalize(cameraUniform.viewInv[2].xyz);\n    let focusPoint = eye + ray.d * (cameraUniform.focusDistance / dot(ray.d, forward));\n    let lens = cameraUniform.apertureRadius * sampleCircle(sample2D(CAMERA_LENS));\n    let origin = (cameraUniform.viewInv * vec4f(lens, 0.0, 1.0)).xyz;\n    initRay(&ray, origin, normalize(focusPoint-origin));\n  }\n  var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;\n  if (!aovSurface(ray, &si, &material)) { return; }\n  let depth = max(dot(si.position-cameraUniform.viewInv[3].xyz, -normalize(cameraUniform.viewInv[2].xyz)), 0.0);\n  let unsupported = (material.materialType != STANDARD) | (material.alpha < 1.0) |\n      (material.coatWeight > 0.0) | (material.fuzzWeight > 0.0) |\n      (material.specularAnisotropy > 0.0) | (material.subsurfaceWeight > 0.0) |\n      (dot(material.mediumScattering, material.mediumScattering) > 0.0) |\n      (material.bevel > 0.0) | (material.thinFilmWeight > 0.0) |\n      (material.transmissionDispersionScale > 0.0);\n  aovAccum[index].depthRoughness += vec4f(depth, material.specularRoughness, material.baseMetalness, select(0.0,1.0,unsupported));\n  if (AOV_REFLECTANCE && (aovParameters.w & 24u) != 0u && !unsupported) {\n    var normal=si.normal; var face=si.faceNormal;\n    var cosine=-dot(normal,ray.d); var faceCosine=-dot(face,ray.d);\n    if (material.thinWalled && (material.transmissionWeight>0.0 || material.subsurfaceWeight>0.0)) {\n      if(cosine<0.0){normal=-normal;cosine=-cosine;}\n      if(faceCosine<0.0){face=-face;faceCosine=-faceCosine;}\n    }\n    if(faceCosine>=0.0 && cosine<0.0){normal=face;cosine=faceCosine;}\n    normal=select(-normal,normal,cosine>=0.0);\n    material.inside=(cosine<=0.0) && !material.thinWalled;\n    material.relativeIor=select(material.specularIor,1.0/material.specularIor,cosine>0.0);\n    let prepared=openpbr_prepare(si.tangent,material,-ray.d,normal);\n    let reflectance=mtFrameReflectanceComponents(prepared, (aovParameters.w & 8u) != 0u, (aovParameters.w & 16u) != 0u);\n    aovAccum[index].diffuse+=vec4f(reflectance.diffuse,1.0);\n    aovAccum[index].specular+=vec4f(reflectance.specular,1.0);\n  }\n  aovAccum[index].position += vec4f(si.position, 1.0);\n}\n@compute @workgroup_size(8, 8)\nfn aovTraceCenter(@builtin(global_invocation_id) id: vec3u) {\n  let size = aovParameters.xy;\n  if (any(id.xy >= size)) { return; }\n  let index = id.x + id.y * size.x;\n  aovCenter[index] = AovCenter(vec4u(0xffffffffu), vec4f(0.0));\n  let uv = (vec2f(id.xy) + 0.5) / vec2f(size);\n  let ray = initCameraRay(uv);\n  var si: SurfaceInteraction; var material: OpenPbrResolvedMaterial;\n  if (!aovSurface(ray, &si, &material)) { return; }\n  let offset = 52u + u32(si.instanceIndex) * 24u;\n  let materialOffset = aovScene[16u];\n  let primitiveOffset = aovScene[17u];\n  aovCenter[index].identity = vec4u(aovScene[offset], aovScene[offset+1u],\n    aovScene[materialOffset + u32(si.materialIndex)], aovScene[primitiveOffset + u32(si.triangleIndex)]);\n  aovCenter[index].motion.z = dot(si.position-aovVector(36u),aovVector(40u));\n  if (aovScene[offset+2u] == 0u) { return; }\n  let oldWorld = aovMatrix(offset+4u) * vec4f(si.position,1.0);\n  let oldClip = aovMatrix(0u) * oldWorld;\n  if (oldClip.w <= 0.0) { return; }\n  let oldNdc=oldClip.xyz/oldClip.w;\n  let oldPixel=(oldNdc.xy*0.5+0.5)*vec2f(size);\n  if (any(oldPixel < vec2f(0.0)) || any(oldPixel >= vec2f(size)) || abs(oldNdc.z)>1.0) { return; }\n  let previous=aovPrevious[u32(oldPixel.x)+u32(oldPixel.y)*size.x];\n  let oldDepth=dot(oldWorld.xyz-aovVector(44u),aovVector(48u));\n  if (any(previous.identity.xy != aovCenter[index].identity.xy) || abs(previous.motion.z-oldDepth)>max(1e-4,abs(oldDepth)*1e-3)) { return; }\n  let motion=(oldPixel-(vec2f(id.xy)+0.5))*vec2f(1.0,-1.0);\n  aovCenter[index].motion=vec4f(motion,aovCenter[index].motion.z,3.0);\n}\n`;\n"]}