@meimu/mutracer 0.3.0 → 0.3.1

package/lib/renderer/glsl/rayTrace.frag.js.map

@@ -1 +1 @@
-{"version":3,"file":"rayTrace.frag.js","sourceRoot":"","sources":["../../../src/renderer/glsl/rayTrace.frag.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,iBAAiB,EAAC,MAAM,yBAAyB,CAAC;AAC1D,OAAO,EAAE,eAAe,EAAE,MAAM,kBAAkB,CAAC;AACnD,OAAO,aAAa,MAAM,4BAA4B,CAAC;AACvD,qDAAqD;AACrD,sDAAsD;AACtD,kFAAkF;AAElF,OAAO,UAAU,MAAM,yBAAyB,CAAC;AACjD,OAAO,EAAC,sBAAsB,EAAC,MAAM,+BAA+B,CAAC;AACrE,OAAO,aAAa,MAAM,gCAAgC,CAAC;AAC3D,OAAO,EAAC,wBAAwB,EAAC,MAAM,iCAAiC,CAAC;AACzE,OAAO,SAAS,MAAM,4BAA4B,CAAC;AACnD,OAAO,QAAQ,MAAM,uBAAuB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,EAAC,wBAAwB,EAAC,MAAM,iCAAiC,CAAC;AACzE,OAAO,eAAe,MAAM,8BAA8B,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAC5D,OAAO,eAAe,MAAM,8BAA8B,CAAC;AAC3D,OAAO,YAAY,MAAM,2BAA2B,CAAC;AACrD,OAAO,EAAC,YAAY,EAAC,MAAM,gBAAgB,CAAC;AAC5C,OAAO,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AACxC,OAAO,gBAAgB,MAAM,+BAA+B,CAAC;AAE7D,MAAM,2BAA2B,GAAG,IAAI,CAAA;;;;;;;CAOvC,CAAC;AAEF,0BAA0B;AAC1B,iEAAiE;AACjE,kFAAkF;AAClF,4CAA4C;AAC5C,MAAM,oBAAoB,GAAG,IAAI,CAAA;;;;;;;;CAQhC,CAAC;AACF,MAAM,iBAAiB,GAAG,IAAI,CAAA;;;0CAGY,eAAe;;;;;;;;;CASxD,CAAC;AAEF,MAAM,CAAC,OAAO,UAAU,4BAA4B,CAAC,GAAG,OAAkB;IACxE,MAAM,QAAQ,GAAG,iBAAiB,EAAE,EAAE,QAAQ,CAAC,OAAO,CAAC,CAAC;IACxD,MAAM,OAAO,GAAG,QAAQ,EAAE,OAAO,IAAI,KAAK,CAAC;IAC3C,MAAM,SAAS,GAAG,QAAQ,EAAE,SAAS,CAAC;IACtC,MAAM,kBAAkB,GAAG,wBAAwB,CAAC,SAAS,CAAC,CAAC;IAC/D,MAAM,YAAY,GAAG,sBAAsB,CAAC,SAAS,CAAC,CAAC;IACvD,MAAM,kBAAkB,GAAG,wBAAwB,CAAC,SAAS,CAAC,CAAC;IAC/D,MAAM,QAAQ,GAAG,UAAU,EAAE,CAAC;IAC9B,OAAO,IAAI,cAAc,CAAC;QAC1B,IAAI,EAAE,cAAc;QACpB,OAAO,EAAE,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,SAAS,CAAC;QACjD,OAAO,EAAE;YACP,OAAO,EAAE,CAAC;SACX;QACD,QAAQ,EAAE;YACR,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,WAAW,CAAC,EAAE,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC;YAC3F,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC;YACjI,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAC,CAAC,CAAC;YAC9C,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC;YAChD,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC;YAEhD,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YACzC,OAAO,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YAC/B,0BAA0B,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,CAAC;YACnD,gBAAgB,EAAE,MAAM,CAAC,OAAO,CAAC,WAAW,CAAC;YAE7C,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YAEpC,eAAe,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;YACzC,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC;YAC1C,qBAAqB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;YAC/C,WAAW,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;SACtC;QACD,QAAQ,EAAE,CAAC,YAAY,EAAE,kBAAkB,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACrH,IAAI,EAAE,IAAI,CAAA;EACV,UAAU;;EAEV,QAAQ;EACR,MAAM,CAAC,IAAI;;EAEX,aAAa;;EAEb,YAAY,CAAC,IAAI;;EAEjB,gBAAgB;;EAEhB,kBAAkB,CAAC,IAAI;EACvB,SAAS,CAAC,IAAI;EACd,MAAM,CAAC,IAAI;;EAEX,YAAY;EACZ,YAAY;EACZ,eAAe;;;;;;;EAOf,eAAe;EACf,QAAQ,EAAE,YAAY,IAAI,EAAE,GAAG,kBAAkB;;EAEjD,aAAa;;;;;;;;;;;;;;;;EAgBb,QAAQ,CAAC,CAAC,CAAC,2BAA2B,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EA0LrC,OAAO,CAAC,CAAC,CAAC,gCAAgC,CAAC,CAAC,CAAC,EAAE,2DAA2D,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;sBAyJ7G,OAAO,CAAC,CAAC,CAAC,UAAU,CAAC;;sGAE2D,CAAC,CAAC,CAAC,EAAE;;;;;;EAM3G,QAAQ,CAAC,CAAC,CAAC,oBAAoB,CAAC,CAAC,CAAC,iBAAiB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAwGnD,SAAS,EAAE,SAAS,IAAI,EAAE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;cAkDZ,QAAQ,EAAE,WAAW,IAAI,YAAY;;;;;;;;;;EAUnD;KACC,CAAC,CAAC;AACL,CAAC","sourcesContent":["import {researchTransport} from '../researchTransport.js';\nimport { MAX_PATH_EVENTS } from '../pathLimits.js';\nimport cameraRayGlsl from './chunks/cameraRay.glsl.js';\n// Portions adapted from ray-tracing-renderer@0.10.15\n// (https://github.com/hoverinc/ray-tracing-renderer).\n// Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.\n\nimport commonGlsl from './chunks/common.glsl.js';\nimport {createRayTraceCoreGlsl} from './chunks/rayTraceCore.glsl.js';\nimport textureLinear from './chunks/textureLinear.glsl.js';\nimport {createMaterialBufferGlsl} from './chunks/materialBuffer.glsl.js';\nimport intersect from './chunks/intersect.glsl.js';\nimport randGlsl from './chunks/rand.glsl.js';\nimport random from './chunks/random.glsl.js';\nimport envMap from './chunks/envMap.glsl.js';\nimport {createSampleMaterialGlsl} from './chunks/sampleMaterial.glsl.js';\nimport sampleLightGlsl from './chunks/sampleLight.glsl.js';\nimport { FragmentShader, glsl, Shader } from '../shader.js';\nimport openPbrBsdfGlsl from './chunks/openPbrBsdf.glsl.js';\nimport samplingGlsl from './chunks/sampling.glsl.js';\nimport {spectrumGlsl} from '../spectrum.js';\nimport { isIntelMac } from '../util.js';\nimport uvProjectionGlsl from './chunks/uvProjection.glsl.js';\n\nconst doubleTransmissionCheckCode = glsl`\n#if defined(HAS_TRANSMISSION) || defined(HAS_ALPHA)\n  if (!si.isTransparency && BOUNCES < 16 && depth * 2 > BOUNCES) {\n    path.abort = true;\n    return;\n  }\n#endif\n`;\n\n// bounce fix on intel mac\n// The intel mac only have one bounce if the loop is not unrolled\n// But windows will have extremely slow shader compilation if the loop is unrolled\n// So here the code are injected separately.\nconst bounceCodeOnIntelMac = glsl`\n// Mark loop can be unrolled\nfor(int _idx_ = 0; _idx_ < BOUNCES; _idx_++) {{\n  if (path.abort) {\n    return vec4(path.li, path.alpha);\n  }\n  bounce(path, depth, si, mat, outAux);\n}}\n`;\nconst bounceCodeOnOther = glsl`\nbool isOpaque = !si.isTransparency;\nint activeBounceLimit = min(BOUNCES, bounceLimit);\nfor(int idx = 0; idx < min(BOUNCES * 3, ${MAX_PATH_EVENTS}); idx++) {\n  if (isOpaque && idx > activeBounceLimit) {\n    path.abort = true;\n  }\n  if (path.abort) {\n    return vec4(path.li, path.alpha);\n  }\n  bounce(path, depth, si, mat, outAux);\n}\n`;\n\nexport default function createRayTraceFragmentShader(...options: unknown[]) {\n  const research = researchTransport()?.fragment(options);\n  const enabled = research?.enabled ?? false;\n  const fragments = research?.fragments;\n  const materialBufferGlsl = createMaterialBufferGlsl(fragments);\n  const rayTraceCore = createRayTraceCoreGlsl(fragments);\n  const sampleMaterialGlsl = createSampleMaterialGlsl(fragments);\n  const intelMac = isIntelMac();\n  return new FragmentShader({\n  name: 'rayTraceFrag',\n  outputs: ['color', 'albedo', 'normal', 'moments'],\n  defines: {\n    BOUNCES: 2\n  },\n  uniforms: {\n    emissiveLightData: Shader.uniform('sampler2D'), emissiveLightColumns: Shader.uniform('int'),\n    emissiveChoiceOffset: Shader.uniform('int'), emissiveChoiceCount: Shader.uniform('int'), emissivePmfOffset: Shader.uniform('int'),\n    emissiveWeightTotal: Shader.uniform('float',0),\n    cameraApertureRadius: Shader.uniform('float', 0),\n    cameraFocusDistance: Shader.uniform('float', 10),\n\n    viewProjectionInv: Shader.uniform('mat4'),\n    viewInv: Shader.uniform('mat4'),\n    transmissiveOutputAlphaMin: Shader.uniform('float'),\n    openPbrEnergyLut: Shader.uniform('sampler2D'),\n\n    viewportSize: Shader.uniform('vec2'),\n\n    spectralChannel: Shader.uniform('int', 0),\n    bevelEnabled: Shader.uniform('bool', true),\n    dispersionWavelengths: Shader.uniform('int', 3),\n    bounceLimit: Shader.uniform('int', 8)\n  },\n  includes: [rayTraceCore, materialBufferGlsl, intersect, random, envMap, ...(research?.mixin ? [research.mixin] : [])],\n  main: glsl`\n${commonGlsl}\n\n${randGlsl}\n${random.main}\n\n${textureLinear}\n\n${rayTraceCore.main}\n\n${uvProjectionGlsl}\n\n${materialBufferGlsl.main}\n${intersect.main}\n${envMap.main}\n\n${samplingGlsl}\n${spectrumGlsl}\n${openPbrBsdfGlsl}\n\nfloat heroWavelength;\nfloat pathWavelength;\nint heroChannel;\nint pathChannel;\n\n${sampleLightGlsl}\n${research?.declarations ?? ''}${sampleMaterialGlsl}\n\n${cameraRayGlsl}\nstruct RayTraceAuxiliary {\n  // Denoising guides and auxiliary pixel classification.\n  vec3 albedo;\n  vec3 normal;\n  bool isBlank;\n  // The primary ray hit a shadow catcher, whose sample the catcher pass owns.\n  bool isShadowCatcher;\n};\n\nvoid bounce(inout Path path, inout int depth, inout SurfaceInteraction si, inout OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux) {\n  samplingBounce = uint(max(depth, 0));\n  // The bounce will be doubled on the transmission material.\n  // For other materials. We will quite early\n  // TODO It is only used on Mac with AMD card. If I enabled it on windows.\n  // It will cause shader compilation too slow.\n${intelMac ? doubleTransmissionCheckCode : ''}\n\n  bool isShadowCatcher = mat.materialType == SHADOW_CATCHER;\n  if (isShadowCatcher) {\n    if (path.fromEye) {\n      // Match the composed shadow with neutral receiver guides, while keeping\n      // background coverage and radiance independent of the invisible ground.\n      outAux.albedo = vec3(1.0);\n      outAux.normal = dot(si.normal, path.ray.d) > 0.0 ? -si.normal : si.normal;\n      si.hit = false;\n      outAux.isShadowCatcher = true;\n    }\n    else {\n      // The catcher belongs only to the shadow-mask pass. Keep its default\n      // white material out of glossy and transmission reflections.\n#ifdef USE_ENVMAP\n      if (path.isNextRefraction) {\n        // The invisible catcher can sit between a transmitted camera ray and\n        // the grounded skybox. Use the displayed background here instead of\n        // terminating early with the raw lighting environment.\n        vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n        bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n        float backgroundMisWeight = path.scatteringIsDelta ? 1.0\n          : sampling_power_heuristic(path.scatteringPdf,\n              envMapPdf(environmentUv(path.ray.d)) * envLightSelectionPmf);\n        path.li += backgroundMisWeight * path.beta * bg.rgb * bg.a;\n        path.alpha = max(path.alpha, bg.a);\n      }\n      else {\n        vec2 envUv = environmentUv(path.ray.d);\n        vec3 emission = textureLinear(envMap, envUv).rgb;\n        emission = sampling_spectral_rgb_filter(emission, pathChannel);\n        path.li += path.beta * emission *\n          sampling_power_heuristic(\n            path.scatteringPdf,\n            envMapPdf(envUv) * envLightSelectionPmf\n          );\n      }\n#else\n      vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n      bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n      path.li += path.beta * bg.rgb * bg.a;\n#endif\n      path.abort = true;\n      return;\n    }\n  }\n\n  if (!si.hit) {\n    // Sample from background if first ray not hit anything.\n    vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n    // Dispersion renders one additive pass per RGB wavelength. Background\n    // radiance has to be split in the same way as lights and emissive\n    // surfaces, otherwise the full skybox is accumulated three times.\n    bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n    if (path.fromEye) {\n      path.li = bg.rgb * bg.a;\n      path.alpha = bg.a;\n      // Only treat as blank pixels when the background is transparent\n      outAux.isBlank = bg.a < 0.01;\n    }\n    else {\n      Ray cam = initCameraRay(v_Texcoord + pixelJitter());\n      bool towardsCam = dot(cam.d, path.ray.d) > 0.0;\n      if (transmissiveOutputAlphaMin < 1.0) {\n        if (towardsCam) {\n          path.alpha = max(path.alpha, bg.a);\n        } else {\n          path.alpha = 1.0;\n        }\n      }\n#if !defined(USE_ENVMAP)\n      // TODO treat background as a light source\n      if (towardsCam) {\n        path.li += path.beta * sampling_power_heuristic(path.scatteringPdf, 1.0) * bg.rgb;\n      }\n#endif\n    }\n    path.abort = true;\n    return;\n  }\n\n  bool sampleIgnored = false;\n  if ((mat.materialType == STANDARD || mat.materialType == UNLIT)) {\n    if (transmissiveOutputAlphaMin < 1.0) {\n      // TODO increase the alpha when there are more bounces.\n      float alpha = max(1.0 - mat.transmissionWeight, 1.0 - pow(1.0 - transmissiveOutputAlphaMin, float(depth)));\n      path.alpha = max(alpha, path.alpha);\n    }\n    float alpha = mat.alpha;\n    if (alpha < 1.0) {\n      // Alpha is sampled as stochastic coverage. An accepted surface branch\n      // must use full radiance weight; multiplying by alpha again is biased.\n      path.alpha = max(path.alpha, alpha);\n      if (hashedSample(3u) > alpha) {\n        sampleIgnored = true;\n      }\n    }\n    if (!sampleIgnored) {\n      SampleMaterial(si, mat, path);\n    }\n    else {\n      path.ray.o = si.position + path.ray.d * EPS;\n    }\n  }\n\n  if (!si.albedoSampled && !sampleIgnored) {\n    // Select denoising guides.\n    if ((mat.materialType == STANDARD || mat.materialType == UNLIT) &&\n      !(mat.transmissionWeight > 0.5 ||\n        (mat.specularRoughness <= 0.02 && mat.baseMetalness > 0.98))) {\n      // Alpha-blended decals can be skipped stochastically. In that case keep\n      // looking for the first useful rough-surface guide behind the decal,\n      // even though it is reached after depth zero.\n      outAux.albedo = mat.baseColor;\n      outAux.normal = dot(si.normal, si.position - vec3(viewInv[3])) > 0.0\n        ? -si.normal\n        : si.normal;\n      si.albedoSampled = true;\n    }\n    if (depth > 4) {\n      si.albedoSampled = true;\n    }\n  }\n\n\n  // Russian Roulette sampling\n#ifdef USE_RR\n  if (depth >= RR_MAX && !sampleIgnored) {\n    // Max-component survival (pbrt-v4) avoids over-amplifying saturated paths.\n    float maxComponent = max(max(path.beta.x, path.beta.y), path.beta.z);\n    float q = max(0.0, 1.0 - maxComponent);\n    if (sample1D(bounceDim(RUSSIAN_ROULETTE)) < q) {\n      // TODO we force use alpha 1.0 here.\n      // To avoid the glass pixels is transparent even if it's on an opaque object.\n      path.alpha = 1.0;\n      path.abort = true;\n    }\n    else {\n      path.beta /= 1.0 - q;\n    }\n  }\n#endif\n\n  if (path.abort) {\n    return;\n  }\n\n  // Next ray hit\n  intersectScene(path.ray, si, true);\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  si.isTransparency = mat.transmissionWeight > 0.001 || mat.subsurfaceWeight > 0.001 || mat.alpha < 0.999;\n\n  // Dynamic preview only becomes spectral when the primary visible surface\n  // actually disperses transmitted light. Direct background and ordinary\n  // opaque materials stay full RGB and therefore remain noise-free.\n  if (dispersionWavelengths == 1 && si.hit && mat.transmissionWeight > 0.001 && mat.transmissionDispersionScale > 0.001) {\n    pathWavelength = heroWavelength;\n    pathChannel = heroChannel;\n  }\n  if (mat.transmissionWeight > 0.001 && pathWavelength > 0.0) {\n    mat.specularIor = sampling_cauchy_ior(\n      pathWavelength,\n      mat.specularIor,\n      mat.transmissionDispersionScale,\n      mat.dispersionAbbeNumber\n    );\n  }\n\n  if (sampleIgnored) {\n    return;\n  }\n\n  bool isAlphaOccluded = hashedSample(3u) < mat.alpha;\n#ifdef NUM_LIGHTS\n  // Use light if hit the light\n  LightSampleRec lightSample;\n  for (int i = 0; i < NUM_LIGHTS; i++) {\n    // Fetch light Data\n    Light light = fetchLight(i);\n    float d = intersectLight(path.ray, light, lightSample);\n    bool occluded = si.hit && d > distance(si.position, path.ray.o);\n    if (d != INF && !(occluded && isAlphaOccluded)) {\n      lightSample.pdf *= light.selectionPmf;\n      float misWeight = sampling_power_heuristic(path.scatteringPdf, lightSample.pdf);\n      lightSample.emission = sampling_spectral_rgb_filter(lightSample.emission, pathChannel);\n      ${enabled ? 'if (!causticHitOwned(path)) { ' : ''}path.li += misWeight * path.beta * lightSample.emission;${enabled ? ' }' : ''}\n      // TODO set alpha to 1.0 will cause some sharp edges\n      path.alpha = 1.0;\n      path.abort = true;\n    }\n  }\n#endif\n\n#ifdef USE_ENVMAP\n  if (!(si.hit && isAlphaOccluded)) {\n    if (path.isNextRefraction) {\n      // Camera-visible transmission must use the displayed background. The\n      // raw environment remains the lighting source for non-transmission rays.\n      vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n      bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n      float backgroundMisWeight = path.scatteringIsDelta ? 1.0\n          : sampling_power_heuristic(path.scatteringPdf,\n              envMapPdf(environmentUv(path.ray.d)) * envLightSelectionPmf);\n        path.li += backgroundMisWeight * path.beta * bg.rgb * bg.a;\n      path.alpha = max(path.alpha, bg.a);\n      path.abort = true;\n      return;\n    }\n    vec2 uv = environmentUv(path.ray.d);\n    vec3 emission = textureLinear(envMap, uv).rgb;\n    float misWeight = sampling_power_heuristic(\n      path.scatteringPdf,\n      envMapPdf(uv) * envLightSelectionPmf\n    );\n    emission = sampling_spectral_rgb_filter(emission, pathChannel);\n    path.li += path.beta * emission * misWeight;\n\n  // TODO force the min alpha to be 1 when showing envmap as background\n  #if !defined(USE_BACKGROUND)\n    path.alpha = 1.0;\n  #endif\n    path.abort = true;\n    return;\n  }\n#endif\n\n  depth++;\n}\n\n\nvec2 bevelCubicSample(float radius, float xi) {\n  float x = 0.25;\n  for (int iteration = 0; iteration < 10; iteration++) {\n    float x2 = x * x;\n    float x3 = x2 * x;\n    float oneMinusX = 1.0 - x;\n    float value = 10.0 * x2 - 20.0 * x3 + 15.0 * x2 * x2 -\n      4.0 * x2 * x3 - xi;\n    float derivative = 20.0 * x * oneMinusX * oneMinusX * oneMinusX;\n    if (abs(value) < 1e-6 || derivative == 0.0) break;\n    x = clamp(x - value / derivative, 0.0, 1.0);\n  }\n  float diskRadius = x * radius;\n  return vec2(\n    diskRadius,\n    sqrt(max(radius * radius - diskRadius * diskRadius, 0.0))\n  );\n}\n\nfloat bevelCubicPdf(float radius, float distanceToSource) {\n  if (distanceToSource >= radius) return 0.0;\n  float falloff = radius - distanceToSource;\n  return falloff * falloff * falloff;\n}\n\n// Shader bevel. Usage: docs/materials.md (Renderer extensions)\nvec3 getBevelNormal(SurfaceInteraction source, float radius) {\n  vec3 geometricNormal = source.faceNormal;\n  if (dot(geometricNormal, source.vertexNormal) < 0.0) geometricNormal = -geometricNormal;\n  vec3 tangent;\n  vec3 bitangent;\n  sampling_orthonormal_basis(geometricNormal, tangent, bitangent);\n  vec3 normalSum = vec3(0.0);\n\n  for (int sampleIndex = 0; sampleIndex < BEVEL_SAMPLES; sampleIndex++) {\n    vec2 bevelRnd = sample2D(CAMERA_BEVEL_BASE + uint(sampleIndex));\n    float axisSample = bevelRnd.x;\n    float radialSample = bevelRnd.y;\n    vec3 diskNormal = geometricNormal;\n    vec3 diskTangent = tangent;\n    vec3 diskBitangent = bitangent;\n    float pickPdfNormal = 0.5;\n    float pickPdfTangent = 0.25;\n    float pickPdfBitangent = 0.25;\n    if (axisSample < 0.5) {\n      axisSample *= 2.0;\n    }\n    else if (axisSample < 0.75) {\n      vec3 swapAxis = diskNormal;\n      diskNormal = diskTangent;\n      diskTangent = swapAxis;\n      pickPdfNormal = 0.25;\n      pickPdfTangent = 0.5;\n      axisSample = (axisSample - 0.5) * 4.0;\n    }\n    else {\n      vec3 swapAxis = diskNormal;\n      diskNormal = diskBitangent;\n      diskBitangent = swapAxis;\n      pickPdfNormal = 0.25;\n      pickPdfBitangent = 0.5;\n      axisSample = (axisSample - 0.75) * 4.0;\n    }\n\n    vec2 diskSample = bevelCubicSample(radius, radialSample);\n    float phi = TWOPI * axisSample;\n    vec3 diskOffset = diskSample.x *\n      (cos(phi) * diskTangent + sin(phi) * diskBitangent);\n    Ray ray;\n    initRay(\n      ray,\n      source.position + diskNormal * diskSample.y + diskOffset,\n      -diskNormal\n    );\n    SurfaceInteraction hit;\n    intersectSceneBevel(\n      ray,\n      hit,\n      source.instanceIndex,\n      max(2.0 * diskSample.y, EPS)\n    );\n    if (!hit.hit) continue;\n\n    float pdfNormal = pickPdfNormal * abs(dot(diskNormal, hit.faceNormal));\n    float pdfTangent = pickPdfTangent * abs(dot(diskTangent, hit.faceNormal));\n    float pdfBitangent = pickPdfBitangent * abs(dot(diskBitangent, hit.faceNormal));\n    float weight = pdfNormal / max(\n      pdfNormal * pdfNormal + pdfTangent * pdfTangent +\n        pdfBitangent * pdfBitangent,\n      1e-8\n    );\n    weight *= bevelCubicPdf(radius, distance(hit.position, source.position)) /\n      max(bevelCubicPdf(radius, diskSample.x), 1e-8);\n    normalSum += weight * hit.vertexNormal;\n  }\n\n  if (dot(normalSum, normalSum) <= 1e-12) return source.normal;\n  vec3 bevelVertexNormal = normalize(normalSum);\n  vec3 result = source.normal + bevelVertexNormal - source.vertexNormal;\n  return dot(result, result) > 1e-12 ? normalize(result) : source.normal;\n}\n\nvec4 integrateOnePathSample(SurfaceInteraction si, Ray ray, OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux) {\n  Path path;\n  path.ray = ray;\n  path.li = vec3(0);\n  path.alpha = transmissiveOutputAlphaMin < 1.0 ? 0.0 : 1.0;\n  path.beta = vec3(1.0);\n  path.fromEye = true;${enabled ? /* glsl */ `\n  path.causticK = 0; path.causticSingleT = false;\n  path.causticReceiverPrimary = false; path.causticHasReceiver = false; path.causticMneeReceiver=false;` : ''}\n  path.scatteringPdf = 0.0;\n  path.scatteringIsDelta = false;\n  path.isNextRefraction = false;\n  path.abort = false;\n  int depth = 0;\n${intelMac ? bounceCodeOnIntelMac : bounceCodeOnOther}\n  return vec4(path.li, path.alpha);\n}\n\nvec4 integrateOnePath(SurfaceInteraction si, Ray ray, OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux) {\n  float coverage = clamp(mat.alpha, 0.0, 1.0);\n  if (!si.hit || (mat.materialType != STANDARD && mat.materialType != UNLIT) || coverage >= 1.0) {\n    return integrateOnePathSample(si, ray, mat, outAux);\n  }\n\n  // glTF BLEND is source-over opacity, not a physical glass boundary. Split\n  // the primary ray so reflected radiance is attenuated by alpha every sample\n  // while the complementary branch continues straight to the scene behind it.\n  // Secondary alpha hits remain stochastic to avoid exponential path growth.\n  Ray backRay = ray;\n  backRay.o = si.position + ray.d * (EPS * 0.02);\n  SurfaceInteraction backSi = si;\n  backSi.albedoSampled = false;\n  intersectScene(backRay, backSi, false);\n  OpenPbrResolvedMaterial backMat;\n  getMaterial(backMat, backSi.materialIndex, backSi.uv, backSi.uv1, backSi.instanceIndex, backSi.vertexColor);\n  RayTraceAuxiliary backAux = outAux;\n  if (backSi.hit && (backMat.materialType == STANDARD || backMat.materialType == UNLIT)) {\n    backAux.albedo = backMat.transmissionWeight > 0.5\n      ? vec3(1.0)\n      : backMat.baseColor;\n    backAux.normal = dot(\n      backSi.normal,\n      backSi.position - vec3(viewInv[3])\n    ) > 0.0\n      ? -backSi.normal\n      : backSi.normal;\n    backSi.albedoSampled = backMat.transmissionWeight > 0.5;\n  }\n  vec4 back = integrateOnePathSample(backSi, backRay, backMat, backAux);\n  if (coverage <= 0.0) {\n    outAux = backAux;\n    return back;\n  }\n\n  SurfaceInteraction frontSi = si;\n  OpenPbrResolvedMaterial frontMat = mat;\n  frontMat.alpha = 1.0;\n  RayTraceAuxiliary frontAux = outAux;\n  vec4 front = integrateOnePathSample(frontSi, ray, frontMat, frontAux);\n  outAux.albedo = mix(backAux.albedo, frontAux.albedo, coverage);\n  vec3 blendedNormal = mix(backAux.normal, frontAux.normal, coverage);\n  outAux.normal = dot(blendedNormal, blendedNormal) > 1e-12\n    ? normalize(blendedNormal)\n    : vec3(0.0);\n  outAux.isBlank = backAux.isBlank && frontAux.isBlank;\n  return mix(back, front, coverage);\n}\n\n// Path tracing integrator as described in\n// http://www.pbr-book.org/3ed-2018/Light_Transport_I_Surface_Reflection/Path_Tracing.html#\nvec4 integrator(Ray ray, out RayTraceAuxiliary outAux) {\n\n  SurfaceInteraction si;\n  OpenPbrResolvedMaterial mat;\n\n  si.albedoSampled = false;\n#ifdef HAS_MEDIUM\n  si.mediumType = 0;\n  si.mediumExtinction = vec3(0.0);\n  si.mediumScattering = vec3(0.0);\n  si.mediumAnisotropy = 0.0;\n#endif\n\n  intersectScene(ray, si, false);\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  si.isTransparency = mat.transmissionWeight > 0.001 || mat.subsurfaceWeight > 0.001 || mat.alpha < 0.999;\n\n#ifdef HAS_BEVEL\n  if (bevelEnabled && mat.bevel > 0.0) {\n    si.normal = getBevelNormal(si, mat.bevel);\n  }\n#endif\n\n  // Keep a clean first-hit fallback for denoising. Glass retains this stable\n  // guide; other delta paths may replace it at their first non-delta hit.\n  if (si.hit && (mat.materialType == STANDARD || mat.materialType == UNLIT)) {\n    outAux.albedo = mat.transmissionWeight > 0.5 ? vec3(1.0) : mat.baseColor;\n    outAux.normal = dot(si.normal, si.position - vec3(viewInv[3])) > 0.0\n      ? -si.normal\n      : si.normal;\n    // A refracted/reflected continuation is stochastic and therefore not a\n    // clean OIDN feature. Keep the stable primary glass guide instead.\n    si.albedoSampled = mat.transmissionWeight > 0.5;\n  }\n\n  if (mat.transmissionDispersionScale > 0.001) {\n    mat.specularIor = sampling_cauchy_ior(\n      pathWavelength,\n      mat.specularIor,\n      mat.transmissionDispersionScale,\n      mat.dispersionAbbeNumber\n    );\n  }\n  return integrateOnePath(si, ray, mat, outAux);\n}\n\nvoid main() {\n\n${fragments?.debugInit ?? ''}  initRandom();\n  rng_initialize(v_Texcoord * viewportSize, int(uint(sampleCount) ^ sampleHash(uint(sequenceSeed))));\n\n\n  pathWavelength = 0.0;\n  pathChannel = 0;\n  if (spectralChannel > 0) {\n    float wavelengthSample = sample1D(CAMERA_WAVELENGTH);\n    if (dispersionWavelengths == 1) {\n      heroChannel = min(int(wavelengthSample * 3.0), 2) + 1;\n      heroWavelength = spectrumCameraWavelength(heroChannel - 1, fract(wavelengthSample * 3.0));\n    } else {\n      pathChannel = spectralChannel;\n      pathWavelength = spectrumCameraWavelength(spectralChannel - 1, wavelengthSample);\n    }\n  }\n\n  // Thin lens model with depth-of-field\n  // http://www.pbr-book.org/3ed-2018/Camera_Models/Projective_Camera_Models.html#TheThinLensModelandDepthofField\n  Ray cam = primaryCameraRay();\n\n  RayTraceAuxiliary outAux;\n  outAux.albedo = vec3(0.0);\n  outAux.normal = vec3(0.0);\n  outAux.isBlank = false;\n  outAux.isShadowCatcher = false;\n\n  vec4 liAndAlpha = integrator(cam, outAux);\n\n  // Clamp at Integration Level to reduce fireflies.\n  // TODO Path level configuration?\n#ifdef CLAMPING\n  float clamping = float(CLAMPING);\n  liAndAlpha.rgb = clamp(liAndAlpha.rgb, -clamping, clamping);\n#endif\n\n  // Each sample will be accumulated three times.\n  if (pathWavelength > 0.0) {\n    // RGB will be normalized using accumated hit in out_albedo.a\n    liAndAlpha.rgb *= 3.0;\n    if (dispersionWavelengths != 1) {\n      // Full spectral samples use three additive passes. Preview only submits\n      // one hero wavelength, so its alpha must not be divided.\n      liAndAlpha.a /= 3.0;\n    }\n  }\n  if (!(liAndAlpha.x < INF && liAndAlpha.x > -EPS)) {\n    liAndAlpha = vec4(0, 0, 0, 1);\n  }\n\n  out_color = ${research?.debugOutput || 'liAndAlpha'};\n  out_albedo = vec4(outAux.albedo, outAux.isBlank ? 0.0 : 1.0);\n  out_normal = vec4(outAux.normal, 1.0);\n  // Moments and count of the built-in denoiser's guide value. The catcher\n  // pass accumulates its own samples (see shadowCatcher.frag).\n  float guideValue = dot(liAndAlpha.rgb, vec3(0.2126, 0.7152, 0.0722));\n  out_moments = outAux.isShadowCatcher\n    ? vec4(0.0)\n    : vec4(guideValue, guideValue * guideValue, 1.0, 0.0);\n\n}`\n  });\n}\n"]}
+{"version":3,"file":"rayTrace.frag.js","sourceRoot":"","sources":["../../../src/renderer/glsl/rayTrace.frag.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,kBAAkB,EAAE,gBAAgB,EAAE,UAAU,EAAE,aAAa,EAAE,oBAAoB,EAAsB,MAAM,iCAAiC,CAAC;AAC3J,OAAO,EAAC,iBAAiB,EAAC,MAAM,yBAAyB,CAAC;AAC1D,OAAO,EAAE,qBAAqB,EAAE,MAAM,kBAAkB,CAAC;AACzD,OAAO,aAAa,MAAM,4BAA4B,CAAC;AACvD,qDAAqD;AACrD,sDAAsD;AACtD,kFAAkF;AAElF,OAAO,UAAU,MAAM,yBAAyB,CAAC;AACjD,OAAO,EAAC,sBAAsB,EAAC,MAAM,+BAA+B,CAAC;AACrE,OAAO,aAAa,MAAM,gCAAgC,CAAC;AAC3D,OAAO,EAAC,wBAAwB,EAAC,MAAM,iCAAiC,CAAC;AACzE,OAAO,SAAS,MAAM,4BAA4B,CAAC;AACnD,OAAO,QAAQ,MAAM,uBAAuB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,EAAC,wBAAwB,EAAC,MAAM,iCAAiC,CAAC;AACzE,OAAO,eAAe,MAAM,8BAA8B,CAAC;AAC3D,OAAO,EAAE,cAAc,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAC5D,OAAO,eAAe,MAAM,8BAA8B,CAAC;AAC3D,OAAO,YAAY,MAAM,2BAA2B,CAAC;AACrD,OAAO,EAAC,YAAY,EAAC,MAAM,gBAAgB,CAAC;AAC5C,OAAO,EAAE,UAAU,EAAE,MAAM,YAAY,CAAC;AACxC,OAAO,gBAAgB,MAAM,+BAA+B,CAAC;AAE7D,MAAM,2BAA2B,GAAG,IAAI,CAAA;;;;;;;CAOvC,CAAC;AAEF,0BAA0B;AAC1B,iEAAiE;AACjE,kFAAkF;AAClF,4CAA4C;AAC5C,MAAM,oBAAoB,GAAG,IAAI,CAAA;;;;;;;;CAQhC,CAAC;AACF,MAAM,iBAAiB,GAAG,IAAI,CAAA;;;0CAGY,qBAAqB;;;;;;;;;CAS9D,CAAC;AAEF,MAAM,CAAC,OAAO,UAAU,4BAA4B,CAAC,GAAG,OAAkB;IACxE,MAAM,uBAAuB,GAAI,OAAO,CAAC,CAAC,CAAqD,EAAE,uBAAuB,IAAI,KAAK,CAAC;IAClI,MAAM,WAAW,GAAI,OAAO,CAAC,CAAC,CAAiD,EAAE,mBAAmB,IAAI,KAAK,CAAC;IAC9G,MAAM,aAAa,GAAmB;QACpC,IAAI,EAAE,CAAC,MAAM,EAAE,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,oBAAoB,EAAE,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,yBAAyB,EAAE,KAAK,CAAC;QACjG,GAAG,EAAE,CAAC,mBAAmB,EAAE,QAAQ,CAAC,EAAE,KAAK,EAAE,CAAC,KAAK,EAAE,OAAO,CAAC,EAAE,KAAK,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC;QACpF,OAAO,EAAE,CAAC,MAAM,EAAE,eAAe,CAAC,EAAE,MAAM,EAAE,CAAC,MAAM,EAAE,UAAU,CAAC,EAAE,GAAG,oBAAoB;KAC1F,CAAC;IACF,MAAM,QAAQ,GAAG,iBAAiB,EAAE,EAAE,QAAQ,CAAC,OAAO,OAAO,CAAC,CAAC,CAAC,KAAK,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC;IAC9F,MAAM,OAAO,GAAG,QAAQ,EAAE,OAAO,IAAI,KAAK,CAAC;IAC3C,MAAM,SAAS,GAAG,QAAQ,EAAE,SAAS,CAAC;IACtC,MAAM,kBAAkB,GAAG,wBAAwB,CAAC,SAAS,CAAC,CAAC;IAC/D,MAAM,YAAY,GAAG,sBAAsB,CAAC,SAAS,CAAC,CAAC;IACvD,MAAM,kBAAkB,GAAG,wBAAwB,CAAC,SAAS,EAAE,uBAAuB,CAAC,CAAC;IACxF,MAAM,QAAQ,GAAG,UAAU,EAAE,CAAC;IAC9B,OAAO,IAAI,cAAc,CAAC;QAC1B,IAAI,EAAE,cAAc;QACpB,OAAO,EAAE,CAAC,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,SAAS,CAAC;QACjD,OAAO,EAAE;YACP,OAAO,EAAE,CAAC;SACX;QACD,QAAQ,EAAE;YACR,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,WAAW,CAAC,EAAE,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC;YAC3F,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,CAAC;YACjI,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAC,CAAC,CAAC;YAC9C,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC;YAChD,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC;YAEhD,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YACzC,aAAa,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YACrC,OAAO,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YAC/B,0BAA0B,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,CAAC;YACnD,gBAAgB,EAAE,MAAM,CAAC,OAAO,CAAC,WAAW,CAAC;YAE7C,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;YAEpC,eAAe,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;YACzC,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,EAAE,IAAI,CAAC;YAC1C,qBAAqB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;YAC/C,WAAW,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;SACtC;QACD,QAAQ,EAAE,CAAC,YAAY,EAAE,kBAAkB,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,EAAE,GAAG,CAAC,QAAQ,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACrH,IAAI,EAAE,IAAI,CAAA;EACV,UAAU;;EAEV,QAAQ;EACR,MAAM,CAAC,IAAI;;EAEX,aAAa;;EAEb,YAAY,CAAC,IAAI;;EAEjB,gBAAgB;;EAEhB,kBAAkB,CAAC,IAAI;EACvB,SAAS,CAAC,IAAI;EACd,MAAM,CAAC,IAAI;;EAEX,YAAY;EACZ,YAAY;EACZ,eAAe;;;;;;;EAOf,eAAe;EACf,QAAQ,EAAE,YAAY,IAAI,EAAE,GAAG,kBAAkB;;EAEjD,aAAa;;;;;;;;;;EAUb,WAAW,CAAC,CAAC,CAAC,kBAAkB,GAAG,gBAAgB,CAAC,WAAW,EAAE,aAAa,CAAC,CAAC,CAAC,CAAC,EAAE,QAAQ,WAAW,CAAC,CAAC,CAAC,oBAAoB,CAAC,CAAC,CAAC,QAAQ,oIAAoI,WAAW,CAAC,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,EAAE,MAAM,WAAW,CAAC,CAAC,CAAC,4BAA4B,CAAC,CAAC,CAAC,EAAE;;;;;;EAM/W,QAAQ,CAAC,CAAC,CAAC,2BAA2B,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAkF3C,WAAW,CAAC,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,+BAA+B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAmE1E,WAAW,CAAC,CAAC,CAAC,UAAU,CAAC;;;CAG1B,CAAC,CAAC,CAAC,EAAE;EACJ,WAAW,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,UAAU,CAAC;;;CAG/B;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAkCO,OAAO,CAAC,CAAC,CAAC,gCAAgC,CAAC,CAAC,CAAC,EAAE,2DAA2D,OAAO,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;sBAyJ7G,OAAO,CAAC,CAAC,CAAC,UAAU,CAAC;;sGAE2D,CAAC,CAAC,CAAC,EAAE;;;;;;EAM3G,WAAW,CAAC,CAAC,CAAC,UAAU,CAAC;;;;;;;EAOrB,aAAa,CAAC,WAAW,EAAE,aAAa,CAAC;;;;iBAI5B,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,SAAS,WAAW,QAAQ,CAAC,CAAC,CAAC,SAAS,CAAC,CAAC,CAAC,oBAAoB,qBAAqB,GAAG;EACzH,QAAQ,CAAC,CAAC,CAAC,kCAAkC,CAAC,CAAC,CAAC,iEAAiE;;;;;;;;;;;;;EAazG,UAAU,CAAC,WAAW,EAAE,aAAa,CAAC;;;;;;4DAMc,CAAC,CAAC,CAAC,UAAU,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,oBAAoB,CAAC,CAAC,CAAC,iBAAiB;kCAC9F;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAoElC,WAAW,CAAC,CAAC,CAAC,2BAA2B,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;;;;;;;;;;;;EAsB9C,WAAW,CAAC,CAAC,CAAC,UAAU,CAAC;;;;;aAKZ,OAAO,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,sEAAsE;;;CAGxG,CAAC,CAAC,CAAC,EAAE;;;;;;;;;;;OAWC,WAAW,CAAC,CAAC,CAAC,yBAAyB,CAAC,CAAC,CAAC,MAAM,IAAI,WAAW,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,EAAE;;EAEzF,SAAS,EAAE,SAAS,IAAI,EAAE;;;;EAI1B,WAAW,CAAC,CAAC,CAAC,0BAA0B,CAAC,CAAC,CAAC,EAAE;;;;;;;;gBAQzB,WAAW,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,iBAAiB;sBACjB,WAAW,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,iBAAiB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;cAqC/E,QAAQ,EAAE,WAAW,IAAI,YAAY;;;;;;;;;;GAUlD,WAAW,CAAC,CAAC,CAAC,UAAU,CAAC;;;;;;;;;;;CAW3B,CAAC,CAAC,CAAC,EAAE,EAAE;KACL,CAAC,CAAC;AACL,CAAC","sourcesContent":["import {spectralPrefixGlsl, prefixCheckpoint, prefixSave, prefixRestore, prefixRandomBindings, type PrefixBindings} from './chunks/spectralPrefix.glsl.js';\nimport {researchTransport} from '../researchTransport.js';\nimport { MAX_PATH_EVENTS_WEBGL } from '../pathLimits.js';\nimport cameraRayGlsl from './chunks/cameraRay.glsl.js';\n// Portions adapted from ray-tracing-renderer@0.10.15\n// (https://github.com/hoverinc/ray-tracing-renderer).\n// Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.\n\nimport commonGlsl from './chunks/common.glsl.js';\nimport {createRayTraceCoreGlsl} from './chunks/rayTraceCore.glsl.js';\nimport textureLinear from './chunks/textureLinear.glsl.js';\nimport {createMaterialBufferGlsl} from './chunks/materialBuffer.glsl.js';\nimport intersect from './chunks/intersect.glsl.js';\nimport randGlsl from './chunks/rand.glsl.js';\nimport random from './chunks/random.glsl.js';\nimport envMap from './chunks/envMap.glsl.js';\nimport {createSampleMaterialGlsl} from './chunks/sampleMaterial.glsl.js';\nimport sampleLightGlsl from './chunks/sampleLight.glsl.js';\nimport { FragmentShader, glsl, Shader } from '../shader.js';\nimport openPbrBsdfGlsl from './chunks/openPbrBsdf.glsl.js';\nimport samplingGlsl from './chunks/sampling.glsl.js';\nimport {spectrumGlsl} from '../spectrum.js';\nimport { isIntelMac } from '../util.js';\nimport uvProjectionGlsl from './chunks/uvProjection.glsl.js';\n\nconst doubleTransmissionCheckCode = glsl`\n#if defined(HAS_TRANSMISSION) || defined(HAS_ALPHA)\n  if (!si.isTransparency && BOUNCES < 16 && depth * 2 > BOUNCES) {\n    path.abort = true;\n    return;\n  }\n#endif\n`;\n\n// bounce fix on intel mac\n// The intel mac only have one bounce if the loop is not unrolled\n// But windows will have extremely slow shader compilation if the loop is unrolled\n// So here the code are injected separately.\nconst bounceCodeOnIntelMac = glsl`\n// Mark loop can be unrolled\nfor(int _idx_ = 0; _idx_ < BOUNCES; _idx_++) {{\n  if (path.abort) {\n    return vec4(path.li, path.alpha);\n  }\n  bounce(path, depth, si, mat, outAux);\n}}\n`;\nconst bounceCodeOnOther = glsl`\nbool isOpaque = !si.isTransparency;\nint activeBounceLimit = min(BOUNCES, bounceLimit);\nfor(int idx = 0; idx < min(BOUNCES * 3, ${MAX_PATH_EVENTS_WEBGL}); idx++) {\n  if (isOpaque && idx > activeBounceLimit) {\n    path.abort = true;\n  }\n  if (path.abort) {\n    return vec4(path.li, path.alpha);\n  }\n  bounce(path, depth, si, mat, outAux);\n}\n`;\n\nexport default function createRayTraceFragmentShader(...options: unknown[]) {\n  const channelIndependentLight = (options[0] as {channelIndependentLight?: boolean} | undefined)?.channelIndependentLight ?? false;\n  const sharePrefix = (options[0] as {shareSpectralPrefix?: boolean} | undefined)?.shareSpectralPrefix ?? false;\n  const surfacePrefix: PrefixBindings = {\n    path: ['Path', 'path'], si: ['SurfaceInteraction', 'si'], mat: ['OpenPbrResolvedMaterial', 'mat'],\n    aux: ['RayTraceAuxiliary', 'outAux'], depth: ['int', 'depth'], event: ['int', 'idx'],\n    ignored: ['bool', 'sampleIgnored'], opaque: ['bool', 'isOpaque'], ...prefixRandomBindings\n  };\n  const research = researchTransport()?.fragment(typeof options[0] === 'object' ? [] : options);\n  const enabled = research?.enabled ?? false;\n  const fragments = research?.fragments;\n  const materialBufferGlsl = createMaterialBufferGlsl(fragments);\n  const rayTraceCore = createRayTraceCoreGlsl(fragments);\n  const sampleMaterialGlsl = createSampleMaterialGlsl(fragments, channelIndependentLight);\n  const intelMac = isIntelMac();\n  return new FragmentShader({\n  name: 'rayTraceFrag',\n  outputs: ['color', 'albedo', 'normal', 'moments'],\n  defines: {\n    BOUNCES: 2\n  },\n  uniforms: {\n    emissiveLightData: Shader.uniform('sampler2D'), emissiveLightColumns: Shader.uniform('int'),\n    emissiveChoiceOffset: Shader.uniform('int'), emissiveChoiceCount: Shader.uniform('int'), emissivePmfOffset: Shader.uniform('int'),\n    emissiveWeightTotal: Shader.uniform('float',0),\n    cameraApertureRadius: Shader.uniform('float', 0),\n    cameraFocusDistance: Shader.uniform('float', 10),\n\n    viewProjectionInv: Shader.uniform('mat4'),\n    projectionInv: Shader.uniform('mat4'),\n    viewInv: Shader.uniform('mat4'),\n    transmissiveOutputAlphaMin: Shader.uniform('float'),\n    openPbrEnergyLut: Shader.uniform('sampler2D'),\n\n    viewportSize: Shader.uniform('vec2'),\n\n    spectralChannel: Shader.uniform('int', 0),\n    bevelEnabled: Shader.uniform('bool', true),\n    dispersionWavelengths: Shader.uniform('int', 3),\n    bounceLimit: Shader.uniform('int', 8)\n  },\n  includes: [rayTraceCore, materialBufferGlsl, intersect, random, envMap, ...(research?.mixin ? [research.mixin] : [])],\n  main: glsl`\n${commonGlsl}\n\n${randGlsl}\n${random.main}\n\n${textureLinear}\n\n${rayTraceCore.main}\n\n${uvProjectionGlsl}\n\n${materialBufferGlsl.main}\n${intersect.main}\n${envMap.main}\n\n${samplingGlsl}\n${spectrumGlsl}\n${openPbrBsdfGlsl}\n\nfloat heroWavelength;\nfloat pathWavelength;\nint heroChannel;\nint pathChannel;\n\n${sampleLightGlsl}\n${research?.declarations ?? ''}${sampleMaterialGlsl}\n\n${cameraRayGlsl}\nstruct RayTraceAuxiliary {\n  // Denoising guides and auxiliary pixel classification.\n  vec3 albedo;\n  vec3 normal;\n  bool isBlank;\n  // The primary ray hit a shadow catcher, whose sample the catcher pass owns.\n  bool isShadowCatcher;\n};\n\n${sharePrefix ? spectralPrefixGlsl + prefixCheckpoint('spSurface', surfacePrefix) : ''}void ${sharePrefix ? 'shadeSurfaceVertex' : 'bounce'}(inout Path path, inout int depth, inout SurfaceInteraction si, inout OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux${sharePrefix ? ', out bool sampleIgnored' : ''}) {${sharePrefix ? '\\n  sampleIgnored = false;' : ''}\n  samplingBounce = uint(max(depth, 0));\n  // The bounce will be doubled on the transmission material.\n  // For other materials. We will quite early\n  // TODO It is only used on Mac with AMD card. If I enabled it on windows.\n  // It will cause shader compilation too slow.\n${intelMac ? doubleTransmissionCheckCode : ''}\n\n  bool isShadowCatcher = mat.materialType == SHADOW_CATCHER;\n  if (isShadowCatcher) {\n    if (path.fromEye) {\n      // Match the composed shadow with neutral receiver guides, while keeping\n      // background coverage and radiance independent of the invisible ground.\n      outAux.albedo = vec3(1.0);\n      outAux.normal = dot(si.normal, path.ray.d) > 0.0 ? -si.normal : si.normal;\n      si.hit = false;\n      outAux.isShadowCatcher = true;\n    }\n    else {\n      // The catcher belongs only to the shadow-mask pass. Keep its default\n      // white material out of glossy and transmission reflections.\n#ifdef USE_ENVMAP\n      if (path.isNextRefraction) {\n        // The invisible catcher can sit between a transmitted camera ray and\n        // the grounded skybox. Use the displayed background here instead of\n        // terminating early with the raw lighting environment.\n        vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n        bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n        float backgroundMisWeight = path.scatteringIsDelta ? 1.0\n          : sampling_power_heuristic(path.scatteringPdf,\n              envMapPdf(environmentUv(path.ray.d)) * envLightSelectionPmf);\n        path.li += backgroundMisWeight * path.beta * bg.rgb * bg.a;\n        path.alpha = max(path.alpha, bg.a);\n      }\n      else {\n        vec2 envUv = environmentUv(path.ray.d);\n        vec3 emission = textureLinear(envMap, envUv).rgb;\n        emission = sampling_spectral_rgb_filter(emission, pathChannel);\n        path.li += path.beta * emission *\n          sampling_power_heuristic(\n            path.scatteringPdf,\n            envMapPdf(envUv) * envLightSelectionPmf\n          );\n      }\n#else\n      vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n      bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n      path.li += path.beta * bg.rgb * bg.a;\n#endif\n      path.abort = true;\n      return;\n    }\n  }\n\n  if (!si.hit) {\n    // Sample from background if first ray not hit anything.\n    vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n    // Dispersion renders one additive pass per RGB wavelength. Background\n    // radiance has to be split in the same way as lights and emissive\n    // surfaces, otherwise the full skybox is accumulated three times.\n    bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n    if (path.fromEye) {\n      path.li = bg.rgb * bg.a;\n      path.alpha = bg.a;\n      // Only treat as blank pixels when the background is transparent\n      outAux.isBlank = bg.a < 0.01;\n    }\n    else {\n      Ray cam = initCameraRay(v_Texcoord + pixelJitter());\n      bool towardsCam = dot(cam.d, path.ray.d) > 0.0;\n      if (transmissiveOutputAlphaMin < 1.0) {\n        if (towardsCam) {\n          path.alpha = max(path.alpha, bg.a);\n        } else {\n          path.alpha = 1.0;\n        }\n      }\n#if !defined(USE_ENVMAP)\n      // TODO treat background as a light source\n      if (towardsCam) {\n        path.li += path.beta * sampling_power_heuristic(path.scatteringPdf, 1.0) * bg.rgb;\n      }\n#endif\n    }\n    path.abort = true;\n    return;\n  }\n\n${sharePrefix ? '  sampleIgnored = false;' : '  bool sampleIgnored = false;'}\n  if ((mat.materialType == STANDARD || mat.materialType == UNLIT)) {\n    if (transmissiveOutputAlphaMin < 1.0) {\n      // TODO increase the alpha when there are more bounces.\n      float alpha = max(1.0 - mat.transmissionWeight, 1.0 - pow(1.0 - transmissiveOutputAlphaMin, float(depth)));\n      path.alpha = max(alpha, path.alpha);\n    }\n    float alpha = mat.alpha;\n    if (alpha < 1.0) {\n      // Alpha is sampled as stochastic coverage. An accepted surface branch\n      // must use full radiance weight; multiplying by alpha again is biased.\n      path.alpha = max(path.alpha, alpha);\n      if (hashedSample(3u) > alpha) {\n        sampleIgnored = true;\n      }\n    }\n    if (!sampleIgnored) {\n      SampleMaterial(si, mat, path);\n    }\n    else {\n      path.ray.o = si.position + path.ray.d * EPS;\n    }\n  }\n\n  if (!si.albedoSampled && !sampleIgnored) {\n    // Select denoising guides.\n    if ((mat.materialType == STANDARD || mat.materialType == UNLIT) &&\n      !(mat.transmissionWeight > 0.5 ||\n        (mat.specularRoughness <= 0.02 && mat.baseMetalness > 0.98))) {\n      // Alpha-blended decals can be skipped stochastically. In that case keep\n      // looking for the first useful rough-surface guide behind the decal,\n      // even though it is reached after depth zero.\n      outAux.albedo = mat.baseColor;\n      outAux.normal = dot(si.normal, si.position - vec3(viewInv[3])) > 0.0\n        ? -si.normal\n        : si.normal;\n      si.albedoSampled = true;\n    }\n    if (depth > 4) {\n      si.albedoSampled = true;\n    }\n  }\n\n\n  // Russian Roulette sampling\n#ifdef USE_RR\n  int rouletteStart = path.scatteringMediumEntered ? RR_MAX : 2;\n  if (depth >= rouletteStart && !sampleIgnored) {\n    // Max-component survival (pbrt-v4) avoids over-amplifying saturated paths.\n    float maxComponent = max(max(path.beta.x, path.beta.y), path.beta.z);\n    float q = max(0.0, 1.0 - maxComponent);\n    if (sample1D(bounceDim(RUSSIAN_ROULETTE)) < q) {\n      // TODO we force use alpha 1.0 here.\n      // To avoid the glass pixels is transparent even if it's on an opaque object.\n      path.alpha = 1.0;\n      path.abort = true;\n    }\n    else {\n      path.beta /= 1.0 - q;\n    }\n  }\n#endif\n\n  if (path.abort) {\n    return;\n  }\n\n${sharePrefix ? /* glsl */ `  return;\n}\nvoid completeSurfaceVertex(inout Path path, inout int depth, inout SurfaceInteraction si, inout OpenPbrResolvedMaterial mat, bool sampleIgnored) {\n` : ''}  // Next ray hit\n${sharePrefix ? '' : /* glsl */ `  intersectScene(path.ray, si, true);\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  si.isTransparency = mat.transmissionWeight > 0.001 || mat.subsurfaceWeight > 0.001 || mat.alpha < 0.999;\n`}\n  // Dynamic preview only becomes spectral when the primary visible surface\n  // actually disperses transmitted light. Direct background and ordinary\n  // opaque materials stay full RGB and therefore remain noise-free.\n  if (dispersionWavelengths == 1 && si.hit && mat.transmissionWeight > 0.001 && mat.transmissionDispersionScale > 0.001) {\n    pathWavelength = heroWavelength;\n    pathChannel = heroChannel;\n  }\n  if (mat.transmissionWeight > 0.001 && pathWavelength > 0.0) {\n    mat.specularIor = sampling_cauchy_ior(\n      pathWavelength,\n      mat.specularIor,\n      mat.transmissionDispersionScale,\n      mat.dispersionAbbeNumber\n    );\n  }\n\n  if (sampleIgnored) {\n    return;\n  }\n\n  bool isAlphaOccluded = hashedSample(3u) < mat.alpha;\n#ifdef NUM_LIGHTS\n  // Use light if hit the light\n  LightSampleRec lightSample;\n  for (int i = 0; i < NUM_LIGHTS; i++) {\n    // Fetch light Data\n    Light light = fetchLight(i);\n    float d = intersectLight(path.ray, light, lightSample);\n    bool occluded = si.hit && d > distance(si.position, path.ray.o);\n    if (d != INF && !(occluded && isAlphaOccluded)) {\n      lightSample.pdf *= light.selectionPmf;\n      float misWeight = sampling_power_heuristic(path.scatteringPdf, lightSample.pdf);\n      lightSample.emission = sampling_spectral_rgb_filter(lightSample.emission, pathChannel);\n      ${enabled ? 'if (!causticHitOwned(path)) { ' : ''}path.li += misWeight * path.beta * lightSample.emission;${enabled ? ' }' : ''}\n      // TODO set alpha to 1.0 will cause some sharp edges\n      path.alpha = 1.0;\n      path.abort = true;\n    }\n  }\n#endif\n\n#ifdef USE_ENVMAP\n  if (!(si.hit && isAlphaOccluded)) {\n    if (path.isNextRefraction) {\n      // Camera-visible transmission must use the displayed background. The\n      // raw environment remains the lighting source for non-transmission rays.\n      vec4 bg = sampleBackgroundFromRay(path.ray.o, path.ray.d);\n      bg.rgb = sampling_spectral_rgb_filter(bg.rgb, pathChannel);\n      float backgroundMisWeight = path.scatteringIsDelta ? 1.0\n          : sampling_power_heuristic(path.scatteringPdf,\n              envMapPdf(environmentUv(path.ray.d)) * envLightSelectionPmf);\n        path.li += backgroundMisWeight * path.beta * bg.rgb * bg.a;\n      path.alpha = max(path.alpha, bg.a);\n      path.abort = true;\n      return;\n    }\n    vec2 uv = environmentUv(path.ray.d);\n    vec3 emission = textureLinear(envMap, uv).rgb;\n    float misWeight = sampling_power_heuristic(\n      path.scatteringPdf,\n      envMapPdf(uv) * envLightSelectionPmf\n    );\n    emission = sampling_spectral_rgb_filter(emission, pathChannel);\n    path.li += path.beta * emission * misWeight;\n\n  // TODO force the min alpha to be 1 when showing envmap as background\n  #if !defined(USE_BACKGROUND)\n    path.alpha = 1.0;\n  #endif\n    path.abort = true;\n    return;\n  }\n#endif\n\n  depth++;\n}\n\n\nvec2 bevelCubicSample(float radius, float xi) {\n  float x = 0.25;\n  for (int iteration = 0; iteration < 10; iteration++) {\n    float x2 = x * x;\n    float x3 = x2 * x;\n    float oneMinusX = 1.0 - x;\n    float value = 10.0 * x2 - 20.0 * x3 + 15.0 * x2 * x2 -\n      4.0 * x2 * x3 - xi;\n    float derivative = 20.0 * x * oneMinusX * oneMinusX * oneMinusX;\n    if (abs(value) < 1e-6 || derivative == 0.0) break;\n    x = clamp(x - value / derivative, 0.0, 1.0);\n  }\n  float diskRadius = x * radius;\n  return vec2(\n    diskRadius,\n    sqrt(max(radius * radius - diskRadius * diskRadius, 0.0))\n  );\n}\n\nfloat bevelCubicPdf(float radius, float distanceToSource) {\n  if (distanceToSource >= radius) return 0.0;\n  float falloff = radius - distanceToSource;\n  return falloff * falloff * falloff;\n}\n\n// Shader bevel. Usage: docs/materials.md (Renderer extensions)\nvec3 getBevelNormal(SurfaceInteraction source, float radius) {\n  vec3 geometricNormal = source.faceNormal;\n  if (dot(geometricNormal, source.vertexNormal) < 0.0) geometricNormal = -geometricNormal;\n  vec3 tangent;\n  vec3 bitangent;\n  sampling_orthonormal_basis(geometricNormal, tangent, bitangent);\n  vec3 normalSum = vec3(0.0);\n\n  for (int sampleIndex = 0; sampleIndex < BEVEL_SAMPLES; sampleIndex++) {\n    vec2 bevelRnd = sample2D(CAMERA_BEVEL_BASE + uint(sampleIndex));\n    float axisSample = bevelRnd.x;\n    float radialSample = bevelRnd.y;\n    vec3 diskNormal = geometricNormal;\n    vec3 diskTangent = tangent;\n    vec3 diskBitangent = bitangent;\n    float pickPdfNormal = 0.5;\n    float pickPdfTangent = 0.25;\n    float pickPdfBitangent = 0.25;\n    if (axisSample < 0.5) {\n      axisSample *= 2.0;\n    }\n    else if (axisSample < 0.75) {\n      vec3 swapAxis = diskNormal;\n      diskNormal = diskTangent;\n      diskTangent = swapAxis;\n      pickPdfNormal = 0.25;\n      pickPdfTangent = 0.5;\n      axisSample = (axisSample - 0.5) * 4.0;\n    }\n    else {\n      vec3 swapAxis = diskNormal;\n      diskNormal = diskBitangent;\n      diskBitangent = swapAxis;\n      pickPdfNormal = 0.25;\n      pickPdfBitangent = 0.5;\n      axisSample = (axisSample - 0.75) * 4.0;\n    }\n\n    vec2 diskSample = bevelCubicSample(radius, radialSample);\n    float phi = TWOPI * axisSample;\n    vec3 diskOffset = diskSample.x *\n      (cos(phi) * diskTangent + sin(phi) * diskBitangent);\n    Ray ray;\n    initRay(\n      ray,\n      source.position + diskNormal * diskSample.y + diskOffset,\n      -diskNormal\n    );\n    SurfaceInteraction hit;\n    intersectSceneBevel(\n      ray,\n      hit,\n      source.instanceIndex,\n      max(2.0 * diskSample.y, EPS)\n    );\n    if (!hit.hit) continue;\n\n    float pdfNormal = pickPdfNormal * abs(dot(diskNormal, hit.faceNormal));\n    float pdfTangent = pickPdfTangent * abs(dot(diskTangent, hit.faceNormal));\n    float pdfBitangent = pickPdfBitangent * abs(dot(diskBitangent, hit.faceNormal));\n    float weight = pdfNormal / max(\n      pdfNormal * pdfNormal + pdfTangent * pdfTangent +\n        pdfBitangent * pdfBitangent,\n      1e-8\n    );\n    weight *= bevelCubicPdf(radius, distance(hit.position, source.position)) /\n      max(bevelCubicPdf(radius, diskSample.x), 1e-8);\n    normalSum += weight * hit.vertexNormal;\n  }\n\n  if (dot(normalSum, normalSum) <= 1e-12) return source.normal;\n  vec3 bevelVertexNormal = normalize(normalSum);\n  vec3 result = source.normal + bevelVertexNormal - source.vertexNormal;\n  return dot(result, result) > 1e-12 ? normalize(result) : source.normal;\n}\n\nvec4 integrateOnePathSample(SurfaceInteraction si, Ray ray, OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux) {\n  Path path; path.scatteringMediumEntered = false;\n  path.ray = ray;\n  path.li = vec3(0);\n  path.alpha = transmissiveOutputAlphaMin < 1.0 ? 0.0 : 1.0;\n  path.beta = vec3(1.0);\n  path.fromEye = true;${enabled ? /* glsl */ `\n  path.causticK = 0; path.causticSingleT = false;\n  path.causticReceiverPrimary = false; path.causticHasReceiver = false; path.causticMneeReceiver=false;` : ''}\n  path.scatteringPdf = 0.0;\n  path.scatteringIsDelta = false;\n  path.isNextRefraction = false;\n  path.abort = false;\n  int depth = 0;\n${sharePrefix ? /* glsl */ `  int spSlot = spCall++;\n  if (spActive && spChannel > 1 && spCompleted[spSlot].ready) return spReplay(spSlot, outAux);\n  bool isOpaque = !si.isTransparency;\n  bool sampleIgnored = false, spResume = false;\n  int spEvent = 0;\n  if (spActive && spChannel > 1 && spForkCall == spSlot) {\n    int idx = 0;\n    ${prefixRestore('spSurface', surfacePrefix)}\n    spEvent = idx; spResume = true; spFork(path);\n  }\n  int activeBounceLimit = min(BOUNCES, bounceLimit);\n  for (int idx = ${intelMac ? '0' : 'spEvent'}; idx < ${intelMac ? 'BOUNCES' : `min(BOUNCES * 3, ${MAX_PATH_EVENTS_WEBGL})`}; idx++) {\n${intelMac ? '    if (idx < spEvent) continue;' : '    if (isOpaque && idx > activeBounceLimit) path.abort = true;'}\n    if (path.abort) break;\n    if (!spResume) {\n      shadeSurfaceVertex(path, depth, si, mat, outAux, sampleIgnored);\n      if (!path.abort) {\n        intersectScene(path.ray, si, true);\n        getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n        si.isTransparency = mat.transmissionWeight > 0.001 || mat.subsurfaceWeight > 0.001 || mat.alpha < 0.999;\n      }\n    }\n    spResume = false;\n    if (!path.abort) {\n      if (spActive && mat.transmissionWeight > 0.001 && spDisperses(mat)) {\n        ${prefixSave('spSurface', surfacePrefix)}\n        spForkCall = spSlot; spFork(path);\n      }\n      completeSurfaceVertex(path, depth, si, mat, sampleIgnored);\n    }\n  }\n  return spFinish(spSlot, vec4(path.li, path.alpha), outAux);` : /* glsl */ `${intelMac ? bounceCodeOnIntelMac : bounceCodeOnOther}\n  return vec4(path.li, path.alpha);`}\n}\n\nvec4 integrateOnePath(SurfaceInteraction si, Ray ray, OpenPbrResolvedMaterial mat, inout RayTraceAuxiliary outAux) {\n  float coverage = clamp(mat.alpha, 0.0, 1.0);\n  if (!si.hit || (mat.materialType != STANDARD && mat.materialType != UNLIT) || coverage >= 1.0) {\n    return integrateOnePathSample(si, ray, mat, outAux);\n  }\n\n  // glTF BLEND is source-over opacity, not a physical glass boundary. Split\n  // the primary ray so reflected radiance is attenuated by alpha every sample\n  // while the complementary branch continues straight to the scene behind it.\n  // Secondary alpha hits remain stochastic to avoid exponential path growth.\n  Ray backRay = ray;\n  backRay.o = si.position + ray.d * (EPS * 0.02);\n  SurfaceInteraction backSi = si;\n  backSi.albedoSampled = false;\n  intersectScene(backRay, backSi, false);\n  OpenPbrResolvedMaterial backMat;\n  getMaterial(backMat, backSi.materialIndex, backSi.uv, backSi.uv1, backSi.instanceIndex, backSi.vertexColor);\n  RayTraceAuxiliary backAux = outAux;\n  if (backSi.hit && (backMat.materialType == STANDARD || backMat.materialType == UNLIT)) {\n    backAux.albedo = backMat.transmissionWeight > 0.5\n      ? vec3(1.0)\n      : backMat.baseColor;\n    backAux.normal = dot(\n      backSi.normal,\n      backSi.position - vec3(viewInv[3])\n    ) > 0.0\n      ? -backSi.normal\n      : backSi.normal;\n    backSi.albedoSampled = backMat.transmissionWeight > 0.5;\n  }\n  vec4 back = integrateOnePathSample(backSi, backRay, backMat, backAux);\n  if (coverage <= 0.0) {\n    outAux = backAux;\n    return back;\n  }\n\n  SurfaceInteraction frontSi = si;\n  OpenPbrResolvedMaterial frontMat = mat;\n  frontMat.alpha = 1.0;\n  RayTraceAuxiliary frontAux = outAux;\n  vec4 front = integrateOnePathSample(frontSi, ray, frontMat, frontAux);\n  outAux.albedo = mix(backAux.albedo, frontAux.albedo, coverage);\n  vec3 blendedNormal = mix(backAux.normal, frontAux.normal, coverage);\n  outAux.normal = dot(blendedNormal, blendedNormal) > 1e-12\n    ? normalize(blendedNormal)\n    : vec3(0.0);\n  outAux.isBlank = backAux.isBlank && frontAux.isBlank;\n  return mix(back, front, coverage);\n}\n\n// Path tracing integrator as described in\n// http://www.pbr-book.org/3ed-2018/Light_Transport_I_Surface_Reflection/Path_Tracing.html#\nvec4 integrator(Ray ray, out RayTraceAuxiliary outAux) {\n\n  SurfaceInteraction si;\n  OpenPbrResolvedMaterial mat;\n\n  si.albedoSampled = false;\n#ifdef HAS_MEDIUM\n  si.mediumType = 0;\n  si.mediumExtinction = vec3(0.0);\n  si.mediumScattering = vec3(0.0);\n  si.mediumAnisotropy = 0.0;\n#endif\n\n${sharePrefix ? '  if (spChannel == 1) {\\n' : ''}  intersectScene(ray, si, false);\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  si.isTransparency = mat.transmissionWeight > 0.001 || mat.subsurfaceWeight > 0.001 || mat.alpha < 0.999;\n\n#ifdef HAS_BEVEL\n  if (bevelEnabled && mat.bevel > 0.0) {\n    si.normal = getBevelNormal(si, mat.bevel);\n  }\n#endif\n\n  // Keep a clean first-hit fallback for denoising. Glass retains this stable\n  // guide; other delta paths may replace it at their first non-delta hit.\n  if (si.hit && (mat.materialType == STANDARD || mat.materialType == UNLIT)) {\n    outAux.albedo = mat.transmissionWeight > 0.5 ? vec3(1.0) : mat.baseColor;\n    outAux.normal = dot(si.normal, si.position - vec3(viewInv[3])) > 0.0\n      ? -si.normal\n      : si.normal;\n    // A refracted/reflected continuation is stochastic and therefore not a\n    // clean OIDN feature. Keep the stable primary glass guide instead.\n    si.albedoSampled = mat.transmissionWeight > 0.5;\n  }\n\n${sharePrefix ? /* glsl */ `    spPrimary = SpectralPrimary(si, mat, outAux, s0, hashedCounter, samplingBounce);\n  } else {\n    si = spPrimary.si; mat = spPrimary.mat; outAux = spPrimary.aux;\n    s0 = spPrimary.rng; hashedCounter = spPrimary.counter; samplingBounce = spPrimary.bounce;\n  }\n  spActive = ${enabled ? 'false' : '!(mat.transmissionDispersionScale > 0.001 && mat.specularIor != 1.0)'};\n  pathChannel = spActive ? 0 : spChannel;\n  spCall = 0;\n` : ''}  if (mat.transmissionDispersionScale > 0.001) {\n    mat.specularIor = sampling_cauchy_ior(\n      pathWavelength,\n      mat.specularIor,\n      mat.transmissionDispersionScale,\n      mat.dispersionAbbeNumber\n    );\n  }\n  return integrateOnePath(si, ray, mat, outAux);\n}\n\nvoid ${sharePrefix ? 'traceSpectralWavelength' : 'main'}(${sharePrefix ? 'int channel' : ''}) {\n\n${fragments?.debugInit ?? ''}  initRandom();\n  rng_initialize(v_Texcoord * viewportSize, int(uint(sampleCount) ^ sampleHash(uint(sequenceSeed))));\n\n\n${sharePrefix ? '  spChannel = channel;\\n' : ''}  pathWavelength = 0.0;\n  pathChannel = 0;\n  if (spectralChannel > 0) {\n    float wavelengthSample = sample1D(CAMERA_WAVELENGTH);\n    if (dispersionWavelengths == 1) {\n      heroChannel = min(int(wavelengthSample * 3.0), 2) + 1;\n      heroWavelength = spectrumCameraWavelength(heroChannel - 1, fract(wavelengthSample * 3.0));\n    } else {\n      pathChannel = ${sharePrefix ? 'spChannel' : 'spectralChannel'};\n      pathWavelength = spectrumCameraWavelength(${sharePrefix ? 'spChannel' : 'spectralChannel'} - 1, wavelengthSample);\n    }\n  }\n\n  // Thin lens model with depth-of-field\n  // http://www.pbr-book.org/3ed-2018/Camera_Models/Projective_Camera_Models.html#TheThinLensModelandDepthofField\n  Ray cam = primaryCameraRay();\n\n  RayTraceAuxiliary outAux;\n  outAux.albedo = vec3(0.0);\n  outAux.normal = vec3(0.0);\n  outAux.isBlank = false;\n  outAux.isShadowCatcher = false;\n\n  vec4 liAndAlpha = integrator(cam, outAux);\n\n  // Clamp at Integration Level to reduce fireflies.\n  // TODO Path level configuration?\n#ifdef CLAMPING\n  float clamping = float(CLAMPING);\n  liAndAlpha.rgb = clamp(liAndAlpha.rgb, -clamping, clamping);\n#endif\n\n  // Each sample will be accumulated three times.\n  if (pathWavelength > 0.0) {\n    // RGB will be normalized using accumated hit in out_albedo.a\n    liAndAlpha.rgb *= 3.0;\n    if (dispersionWavelengths != 1) {\n      // Full spectral samples use three additive passes. Preview only submits\n      // one hero wavelength, so its alpha must not be divided.\n      liAndAlpha.a /= 3.0;\n    }\n  }\n  if (!(liAndAlpha.x < INF && liAndAlpha.x > -EPS)) {\n    liAndAlpha = vec4(0, 0, 0, 1);\n  }\n\n  out_color = ${research?.debugOutput || 'liAndAlpha'};\n  out_albedo = vec4(outAux.albedo, outAux.isBlank ? 0.0 : 1.0);\n  out_normal = vec4(outAux.normal, 1.0);\n  // Moments and count of the built-in denoiser's guide value. The catcher\n  // pass accumulates its own samples (see shadowCatcher.frag).\n  float guideValue = dot(liAndAlpha.rgb, vec3(0.2126, 0.7152, 0.0722));\n  out_moments = outAux.isShadowCatcher\n    ? vec4(0.0)\n    : vec4(guideValue, guideValue * guideValue, 1.0, 0.0);\n\n}${sharePrefix ? /* glsl */ `\nvoid main() {\n  spCompleted[0].ready = false; spCompleted[1].ready = false;\n  vec4 color = vec4(0.0), albedo = vec4(0.0), normal = vec4(0.0), moments = vec4(0.0);\n  // A uniform bound keeps the driver from cloning the complete optical graph.\n  for (int channel = 1; channel <= dispersionWavelengths; channel++) {\n    traceSpectralWavelength(channel);\n    color += out_color; albedo += out_albedo; normal += out_normal; moments += out_moments;\n  }\n  out_color = color; out_albedo = albedo; out_normal = normal; out_moments = moments;\n}\n` : ''}`\n  });\n}\n"]}

package/lib/renderer/glsl/shadowCatcher.frag.js

@@ -26,6 +26,7 @@ export default new FragmentShader({
         spectralChannel: Shader.uniform('int', 0),
         viewportSize: Shader.uniform('vec2'),
         viewProjectionInv: Shader.uniform('mat4'),
+        projectionInv: Shader.uniform('mat4'),
         viewInv: Shader.uniform('mat4'),
         cameraApertureRadius: Shader.uniform('float', 0),
         cameraFocusDistance: Shader.uniform('float', 10)
@@ -52,10 +53,10 @@ ${cameraRayGlsl}
 
 
 
-vec4 hR(inout Ray ray, out vec4 background, out vec2 shadowCounts) {
+vec4 hq(inout Ray ray, out vec4 background, out vec2 shadowCounts) {
 background = vec4(0.0);
 shadowCounts = vec2(0.0);
-Path path;
+Path path; path.scatteringMediumEntered = false;
 path.ray = ray;
 path.li = vec3(0);
 path.alpha = 1.0;
@@ -82,7 +83,7 @@ if (mat.materialType != SHADOW_CATCHER) {
 return vec4(0.0, 0.0, 0.0, 1.0);
 }
 
-vec4 bg = cn(ray.o, ray.d);
+vec4 bg = cq(ray.o, ray.d);
 background = vec4(bg.rgb * bg.a, bg.a);
 if (!path.abort) {
 shadowCounts = sampleShadowCatcher(si, 0, path);
@@ -91,14 +92,14 @@ return vec4(1.0, bg.a, bg.a >= 0.01 ? 1.0 : 0.0, 0.0);
 }
 
 void main() {
-rA();
-fi(v_Texcoord * viewportSize, int(uint(sampleCount) ^ cB(uint(sequenceSeed))));
+sP();
+fA(v_Texcoord * viewportSize, int(uint(sampleCount) ^ cD(uint(sequenceSeed))));
 
 if (spectralChannel > 0) { sample1D(CAMERA_WAVELENGTH); }
-Ray cam = ow();
+Ray cam = pq();
 vec4 background;
 vec2 shadowCounts;
-vec4 liAndAlpha = hR(cam, background, shadowCounts);
+vec4 liAndAlpha = hq(cam, background, shadowCounts);
 
 if (!(liAndAlpha.x < INF && liAndAlpha.x > -EPS)) {
 liAndAlpha = vec4(0.0);
@@ -112,11 +113,11 @@ out_catcherBackground = background;
 
 
 float visibility = 1.0 - (shadowCounts.y > 0.0 ? shadowCounts.x / shadowCounts.y : 0.0);
-float jw = visibility * (background.a >= 0.01
+float jZ = visibility * (background.a >= 0.01
 ? dot(background.rgb / background.a, vec3(0.2126, 0.7152, 0.0722))
 : 1.0);
 out_catcherMoments = liAndAlpha.r > 0.0
-? vec4(jw, jw * jw, 1.0, shadowCounts.y)
+? vec4(jZ, jZ * jZ, 1.0, shadowCounts.y)
 : vec4(0.0);
 
 }`

package/lib/renderer/glsl/shadowCatcher.frag.js.map

@@ -1 +1 @@
-{"version":3,"file":"shadowCatcher.frag.js","sourceRoot":"","sources":["../../../src/renderer/glsl/shadowCatcher.frag.ts"],"names":[],"mappings":"AAAA,qDAAqD;AACrD,sDAAsD;AACtD,kFAAkF;AAElF,OAAO,aAAa,MAAM,4BAA4B,CAAC;AACvD,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,YAAY,MAAM,+BAA+B,CAAC;AACzD,OAAO,aAAa,MAAM,gCAAgC,CAAC;AAC3D,OAAO,kBAAkB,MAAM,iCAAiC,CAAC;AACjE,OAAO,SAAS,MAAM,4BAA4B,CAAC;AACnD,OAAO,QAAQ,MAAM,uBAAuB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,mBAAmB,MAAM,sCAAsC,CAAC;AACvE,OAAO,WAAW,MAAM,8BAA8B,CAAC;AACvD,OAAO,EAAE,cAAc,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAC5D,OAAO,YAAY,MAAM,2BAA2B,CAAC;AACrD,OAAO,gBAAgB,MAAM,+BAA+B,CAAC;AAE7D,eAAe,IAAI,cAAc,CAAC;IAChC,IAAI,EAAE,eAAe;IACrB,OAAO,EAAE,CAAC,eAAe,EAAE,mBAAmB,EAAE,gBAAgB,CAAC;IACjE,OAAO,EAAE;QACP,OAAO,EAAE,CAAC;KACX;IACD,QAAQ,EAAE;QACR,qBAAqB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;QAC/C,eAAe,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;QACzC,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QACpC,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QACzC,OAAO,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QAC/B,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC;QAChD,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC;KACjD;IACD,QAAQ,EAAE,CAAC,YAAY,EAAE,kBAAkB,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC;IACvE,IAAI,EAAE,IAAI,CAAA;EACV,MAAM;EACN,QAAQ;EACR,MAAM,CAAC,IAAI;;EAEX,aAAa;;EAEb,YAAY,CAAC,IAAI;EACjB,gBAAgB;;EAEhB,kBAAkB,CAAC,IAAI;EACvB,SAAS,CAAC,IAAI;EACd,MAAM,CAAC,IAAI;;EAEX,YAAY;EACZ,WAAW;EACX,mBAAmB;EACnB,aAAa;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAuEb;CACD,CAAC,CAAC","sourcesContent":["// Portions adapted from ray-tracing-renderer@0.10.15\n// (https://github.com/hoverinc/ray-tracing-renderer).\n// Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.\n\nimport cameraRayGlsl from './chunks/cameraRay.glsl.js';\nimport common from './chunks/common.glsl.js';\nimport rayTraceCore from './chunks/rayTraceCore.glsl.js';\nimport textureLinear from './chunks/textureLinear.glsl.js';\nimport materialBufferGlsl from './chunks/materialBuffer.glsl.js';\nimport intersect from './chunks/intersect.glsl.js';\nimport randGlsl from './chunks/rand.glsl.js';\nimport random from './chunks/random.glsl.js';\nimport envMap from './chunks/envMap.glsl.js';\nimport sampleShadowCatcher from './chunks/sampleShadowCatcher.glsl.js';\nimport sampleLight from './chunks/sampleLight.glsl.js';\nimport { FragmentShader, glsl, Shader } from '../shader.js';\nimport samplingGlsl from './chunks/sampling.glsl.js';\nimport uvProjectionGlsl from './chunks/uvProjection.glsl.js';\n\nexport default new FragmentShader({\n  name: 'shadowCatcher',\n  outputs: ['shadowCatcher', 'catcherBackground', 'catcherMoments'],\n  defines: {\n    BOUNCES: 2\n  },\n  uniforms: {\n    dispersionWavelengths: Shader.uniform('int', 3),\n    spectralChannel: Shader.uniform('int', 0),\n    viewportSize: Shader.uniform('vec2'),\n    viewProjectionInv: Shader.uniform('mat4'),\n    viewInv: Shader.uniform('mat4'),\n    cameraApertureRadius: Shader.uniform('float', 0),\n    cameraFocusDistance: Shader.uniform('float', 10)\n  },\n  includes: [rayTraceCore, materialBufferGlsl, intersect, random, envMap],\n  main: glsl`\n${common}\n${randGlsl}\n${random.main}\n\n${textureLinear}\n\n${rayTraceCore.main}\n${uvProjectionGlsl}\n\n${materialBufferGlsl.main}\n${intersect.main}\n${envMap.main}\n\n${samplingGlsl}\n${sampleLight}\n${sampleShadowCatcher}\n${cameraRayGlsl}\n\n// Path tracing integrator as described in\n// http://www.pbr-book.org/3ed-2018/Light_Transport_I_Surface_Reflection/Path_Tracing.html#\nvec4 integrator(inout Ray ray, out vec4 background, out vec2 shadowCounts) {\n  background = vec4(0.0);\n  shadowCounts = vec2(0.0);\n  Path path;\n  path.ray = ray;\n  path.li = vec3(0);\n  path.alpha = 1.0;\n  path.beta = vec3(1.0);\n  path.fromEye = true;\n  path.scatteringPdf = 0.0;\n  path.abort = false;\n\n  SurfaceInteraction si;\n\n  intersectScene(path.ray, si, false);\n  if (!si.hit) {\n    // Red stores catcher coverage. Keep it zero for pixels where the primary\n    // ray did not hit the catcher so denoised shadows cannot bleed onto the\n    // model silhouette during composition.\n    return vec4(0.0, 0.0, 0.0, 1.0);\n  }\n\n  OpenPbrResolvedMaterial mat;\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  if (mat.materialType != SHADOW_CATCHER) {\n    // The main path tracer already renders regular geometry. Only evaluate\n    // shadow opacity when the primary ray actually hits the catcher.\n    return vec4(0.0, 0.0, 0.0, 1.0);\n  }\n\n  vec4 bg = sampleBackgroundFromRay(ray.o, ray.d);\n  background = vec4(bg.rgb * bg.a, bg.a);\n  if (!path.abort) {\n    shadowCounts = sampleShadowCatcher(si, 0, path);\n  }\n  return vec4(1.0, bg.a, bg.a >= 0.01 ? 1.0 : 0.0, 0.0);\n}\n\nvoid main() {\n  initRandom();\n  rng_initialize(v_Texcoord * viewportSize, int(uint(sampleCount) ^ sampleHash(uint(sequenceSeed))));\n  // Match the beauty pass's random stream in spectral renders.\n  if (spectralChannel > 0) { sample1D(CAMERA_WAVELENGTH); }\n  Ray cam = primaryCameraRay();\n  vec4 background;\n  vec2 shadowCounts;\n  vec4 liAndAlpha = integrator(cam, background, shadowCounts);\n\n  if (!(liAndAlpha.x < INF && liAndAlpha.x > -EPS)) {\n    liAndAlpha = vec4(0.0);\n    shadowCounts = vec2(0.0);\n  }\n\n  // Alpha and the moments' w sum the occluded and total light sample counts;\n  // the output pass resolves their ratio into shadow visibility.\n  out_shadowCatcher = vec4(liAndAlpha.rgb, shadowCounts.x);\n  out_catcherBackground = background;\n  // The composed catcher is the background darkened by the shadow. Over a\n  // transparent background the shadow is only its opacity.\n  float visibility = 1.0 - (shadowCounts.y > 0.0 ? shadowCounts.x / shadowCounts.y : 0.0);\n  float guideValue = visibility * (background.a >= 0.01\n    ? dot(background.rgb / background.a, vec3(0.2126, 0.7152, 0.0722))\n    : 1.0);\n  out_catcherMoments = liAndAlpha.r > 0.0\n    ? vec4(guideValue, guideValue * guideValue, 1.0, shadowCounts.y)\n    : vec4(0.0);\n\n}`\n});\n"]}
+{"version":3,"file":"shadowCatcher.frag.js","sourceRoot":"","sources":["../../../src/renderer/glsl/shadowCatcher.frag.ts"],"names":[],"mappings":"AAAA,qDAAqD;AACrD,sDAAsD;AACtD,kFAAkF;AAElF,OAAO,aAAa,MAAM,4BAA4B,CAAC;AACvD,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,YAAY,MAAM,+BAA+B,CAAC;AACzD,OAAO,aAAa,MAAM,gCAAgC,CAAC;AAC3D,OAAO,kBAAkB,MAAM,iCAAiC,CAAC;AACjE,OAAO,SAAS,MAAM,4BAA4B,CAAC;AACnD,OAAO,QAAQ,MAAM,uBAAuB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,MAAM,MAAM,yBAAyB,CAAC;AAC7C,OAAO,mBAAmB,MAAM,sCAAsC,CAAC;AACvE,OAAO,WAAW,MAAM,8BAA8B,CAAC;AACvD,OAAO,EAAE,cAAc,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,cAAc,CAAC;AAC5D,OAAO,YAAY,MAAM,2BAA2B,CAAC;AACrD,OAAO,gBAAgB,MAAM,+BAA+B,CAAC;AAE7D,eAAe,IAAI,cAAc,CAAC;IAChC,IAAI,EAAE,eAAe;IACrB,OAAO,EAAE,CAAC,eAAe,EAAE,mBAAmB,EAAE,gBAAgB,CAAC;IACjE,OAAO,EAAE;QACP,OAAO,EAAE,CAAC;KACX;IACD,QAAQ,EAAE;QACR,qBAAqB,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;QAC/C,eAAe,EAAE,MAAM,CAAC,OAAO,CAAC,KAAK,EAAE,CAAC,CAAC;QACzC,YAAY,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QACpC,iBAAiB,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QACzC,aAAa,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QACrC,OAAO,EAAE,MAAM,CAAC,OAAO,CAAC,MAAM,CAAC;QAC/B,oBAAoB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,CAAC,CAAC;QAChD,mBAAmB,EAAE,MAAM,CAAC,OAAO,CAAC,OAAO,EAAE,EAAE,CAAC;KACjD;IACD,QAAQ,EAAE,CAAC,YAAY,EAAE,kBAAkB,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,CAAC;IACvE,IAAI,EAAE,IAAI,CAAA;EACV,MAAM;EACN,QAAQ;EACR,MAAM,CAAC,IAAI;;EAEX,aAAa;;EAEb,YAAY,CAAC,IAAI;EACjB,gBAAgB;;EAEhB,kBAAkB,CAAC,IAAI;EACvB,SAAS,CAAC,IAAI;EACd,MAAM,CAAC,IAAI;;EAEX,YAAY;EACZ,WAAW;EACX,mBAAmB;EACnB,aAAa;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EAuEb;CACD,CAAC,CAAC","sourcesContent":["// Portions adapted from ray-tracing-renderer@0.10.15\n// (https://github.com/hoverinc/ray-tracing-renderer).\n// Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.\n\nimport cameraRayGlsl from './chunks/cameraRay.glsl.js';\nimport common from './chunks/common.glsl.js';\nimport rayTraceCore from './chunks/rayTraceCore.glsl.js';\nimport textureLinear from './chunks/textureLinear.glsl.js';\nimport materialBufferGlsl from './chunks/materialBuffer.glsl.js';\nimport intersect from './chunks/intersect.glsl.js';\nimport randGlsl from './chunks/rand.glsl.js';\nimport random from './chunks/random.glsl.js';\nimport envMap from './chunks/envMap.glsl.js';\nimport sampleShadowCatcher from './chunks/sampleShadowCatcher.glsl.js';\nimport sampleLight from './chunks/sampleLight.glsl.js';\nimport { FragmentShader, glsl, Shader } from '../shader.js';\nimport samplingGlsl from './chunks/sampling.glsl.js';\nimport uvProjectionGlsl from './chunks/uvProjection.glsl.js';\n\nexport default new FragmentShader({\n  name: 'shadowCatcher',\n  outputs: ['shadowCatcher', 'catcherBackground', 'catcherMoments'],\n  defines: {\n    BOUNCES: 2\n  },\n  uniforms: {\n    dispersionWavelengths: Shader.uniform('int', 3),\n    spectralChannel: Shader.uniform('int', 0),\n    viewportSize: Shader.uniform('vec2'),\n    viewProjectionInv: Shader.uniform('mat4'),\n    projectionInv: Shader.uniform('mat4'),\n    viewInv: Shader.uniform('mat4'),\n    cameraApertureRadius: Shader.uniform('float', 0),\n    cameraFocusDistance: Shader.uniform('float', 10)\n  },\n  includes: [rayTraceCore, materialBufferGlsl, intersect, random, envMap],\n  main: glsl`\n${common}\n${randGlsl}\n${random.main}\n\n${textureLinear}\n\n${rayTraceCore.main}\n${uvProjectionGlsl}\n\n${materialBufferGlsl.main}\n${intersect.main}\n${envMap.main}\n\n${samplingGlsl}\n${sampleLight}\n${sampleShadowCatcher}\n${cameraRayGlsl}\n\n// Path tracing integrator as described in\n// http://www.pbr-book.org/3ed-2018/Light_Transport_I_Surface_Reflection/Path_Tracing.html#\nvec4 integrator(inout Ray ray, out vec4 background, out vec2 shadowCounts) {\n  background = vec4(0.0);\n  shadowCounts = vec2(0.0);\n  Path path; path.scatteringMediumEntered = false;\n  path.ray = ray;\n  path.li = vec3(0);\n  path.alpha = 1.0;\n  path.beta = vec3(1.0);\n  path.fromEye = true;\n  path.scatteringPdf = 0.0;\n  path.abort = false;\n\n  SurfaceInteraction si;\n\n  intersectScene(path.ray, si, false);\n  if (!si.hit) {\n    // Red stores catcher coverage. Keep it zero for pixels where the primary\n    // ray did not hit the catcher so denoised shadows cannot bleed onto the\n    // model silhouette during composition.\n    return vec4(0.0, 0.0, 0.0, 1.0);\n  }\n\n  OpenPbrResolvedMaterial mat;\n  getMaterial(mat, si.materialIndex, si.uv, si.uv1, si.instanceIndex, si.vertexColor);\n  if (mat.materialType != SHADOW_CATCHER) {\n    // The main path tracer already renders regular geometry. Only evaluate\n    // shadow opacity when the primary ray actually hits the catcher.\n    return vec4(0.0, 0.0, 0.0, 1.0);\n  }\n\n  vec4 bg = sampleBackgroundFromRay(ray.o, ray.d);\n  background = vec4(bg.rgb * bg.a, bg.a);\n  if (!path.abort) {\n    shadowCounts = sampleShadowCatcher(si, 0, path);\n  }\n  return vec4(1.0, bg.a, bg.a >= 0.01 ? 1.0 : 0.0, 0.0);\n}\n\nvoid main() {\n  initRandom();\n  rng_initialize(v_Texcoord * viewportSize, int(uint(sampleCount) ^ sampleHash(uint(sequenceSeed))));\n  // Match the beauty pass's random stream in spectral renders.\n  if (spectralChannel > 0) { sample1D(CAMERA_WAVELENGTH); }\n  Ray cam = primaryCameraRay();\n  vec4 background;\n  vec2 shadowCounts;\n  vec4 liAndAlpha = integrator(cam, background, shadowCounts);\n\n  if (!(liAndAlpha.x < INF && liAndAlpha.x > -EPS)) {\n    liAndAlpha = vec4(0.0);\n    shadowCounts = vec2(0.0);\n  }\n\n  // Alpha and the moments' w sum the occluded and total light sample counts;\n  // the output pass resolves their ratio into shadow visibility.\n  out_shadowCatcher = vec4(liAndAlpha.rgb, shadowCounts.x);\n  out_catcherBackground = background;\n  // The composed catcher is the background darkened by the shadow. Over a\n  // transparent background the shadow is only its opacity.\n  float visibility = 1.0 - (shadowCounts.y > 0.0 ? shadowCounts.x / shadowCounts.y : 0.0);\n  float guideValue = visibility * (background.a >= 0.01\n    ? dot(background.rgb / background.a, vec3(0.2126, 0.7152, 0.0722))\n    : 1.0);\n  out_catcherMoments = liAndAlpha.r > 0.0\n    ? vec4(guideValue, guideValue * guideValue, 1.0, shadowCounts.y)\n    : vec4(0.0);\n\n}`\n});\n"]}

package/lib/renderer/glsl/webglPipeline.frag.js

@@ -7,22 +7,22 @@ export const brightFragment = new FragmentShader({
         scale: Shader.uniform('float', 1)
     },
     main: glsl `
-const vec3 xc = vec3(0.2125, 0.7154, 0.0721);
-vec4 tV(vec4 a, vec4 b, vec4 c) {
+const vec3 ze = vec3(0.2125, 0.7154, 0.0721);
+vec4 vE(vec4 a, vec4 b, vec4 c) {
 return a + b + c - min(min(a, b), c) - max(max(a, b), c);
 }
 void main() {
 vec2 d = 1.0 / vec2(colorTexRegion());
-vec4 cR = sample_colorTex(v_Texcoord);
+vec4 texel = sample_colorTex(v_Texcoord);
 vec4 s1 = sample_colorTex(v_Texcoord - vec2(d.x, 0.0));
 vec4 s2 = sample_colorTex(v_Texcoord + vec2(d.x, 0.0));
 vec4 s3 = sample_colorTex(v_Texcoord - vec2(0.0, d.y));
 vec4 s4 = sample_colorTex(v_Texcoord + vec2(0.0, d.y));
-cR = tV(tV(cR, s1, s2), s3, s4);
-float lum = dot(cR.rgb, xc);
-out_color = lum > threshold && cR.a > 0.0
-? vec4(cR.rgb * scale, cR.a)
-: vec4(0.0, 0.0, 0.0, cR.a);
+texel = vE(vE(texel, s1, s2), s3, s4);
+float lum = dot(texel.rgb, ze);
+out_color = lum > threshold && texel.a > 0.0
+? vec4(texel.rgb * scale, texel.a)
+: vec4(0.0, 0.0, 0.0, texel.a);
 }`
 });
 export const composeFragment = new FragmentShader({
@@ -48,29 +48,29 @@ export const composeFragment = new FragmentShader({
         sRGBOutput: Shader.uniform('bool', 1)
     },
     main: glsl `
-vec3 nT(vec3 color) {
+vec3 oS(vec3 color) {
 color = max(color, 0.0);
 return mix(1.055 * pow(color, vec3(1.0 / 2.4)) - 0.055,
 color * 12.92, lessThanEqual(color, vec3(0.0031308)));
 }
-vec3 kb(vec3 color) {
+vec3 kH(vec3 color) {
 return clamp(color, 0.0, 1.0);
 }
-vec3 iN(vec3 color) {
+vec3 jp(vec3 color) {
 return clamp(color / (vec3(1.0) + color), 0.0, 1.0);
 }
-vec3 jP(vec3 color) {
+vec3 kw(vec3 color) {
 color = max(vec3(0.0), color - 0.004);
 return pow((color * (6.2 * color + 0.5)) /
 (color * (6.2 * color + 1.7) + 0.06), vec3(2.2));
 }
-vec3 pm(vec3 value) {
-vec3 iF = value * (value + 0.0245786) - 0.000090537;
+vec3 qw(vec3 value) {
+vec3 je = value * (value + 0.0245786) - 0.000090537;
 vec3 ah = value * (0.983729 * value + 0.4329510) + 0.238081;
-return iF / ah;
+return je / ah;
 }
-vec3 lO(vec3 color) {
-const mat3 qk = mat3(
+vec3 mH(vec3 color) {
+const mat3 rt = mat3(
 vec3(0.59719, 0.07600, 0.02840),
 vec3(0.35458, 0.90834, 0.13383),
 vec3(0.04823, 0.01566, 0.83777)
@@ -80,8 +80,8 @@ vec3(1.60475, -0.10208, -0.00327),
 vec3(-0.53108, 1.10813, -0.07276),
 vec3(-0.07367, -0.00605, 1.07602)
 );
-color = qk * (color / 0.6);
-return clamp(outputMatrix * pm(color), 0.0, 1.0);
+color = rt * (color / 0.6);
+return clamp(outputMatrix * qw(color), 0.0, 1.0);
 }
 const mat3 LINEAR_REC2020_TO_LINEAR_SRGB = mat3(
 vec3(1.6605, -0.1246, -0.0182),
@@ -93,44 +93,44 @@ vec3(0.6274, 0.0691, 0.0164),
 vec3(0.3293, 0.9195, 0.0880),
 vec3(0.0433, 0.0113, 0.8956)
 );
-vec3 fT(vec3 value) {
-vec3 nB = value * value;
-vec3 rO = nB * nB;
-return 15.5 * rO * nB - 40.14 * rO * value
-+ 31.96 * rO - 6.868 * nB * value
-+ 0.4298 * nB + 0.1191 * value - 0.00232;
-}
-vec3 my(vec3 color) {
-const mat3 qm = mat3(
+vec3 gj(vec3 value) {
+vec3 ox = value * value;
+vec3 tf = ox * ox;
+return 15.5 * tf * ox - 40.14 * tf * value
++ 31.96 * tf - 6.868 * ox * value
++ 0.4298 * ox + 0.1191 * value - 0.00232;
+}
+vec3 nC(vec3 color) {
+const mat3 rv = mat3(
 vec3(0.856627153315983, 0.137318972929847, 0.11189821299995),
 vec3(0.0951212405381588, 0.761241990602591, 0.0767994186031903),
 vec3(0.0482516061458583, 0.101439036467562, 0.811302368396859)
 );
-const mat3 ph = mat3(
+const mat3 qr = mat3(
 vec3(1.1271005818144368, -0.1413297634984383, -0.14132976349843826),
 vec3(-0.11060664309660323, 1.157823702216272, -0.11060664309660294),
 vec3(-0.016493938717834573, -0.016493938717834257, 1.2519364065950405)
 );
-color = qm * (LINEAR_SRGB_TO_LINEAR_REC2020 * color);
+color = rv * (LINEAR_SRGB_TO_LINEAR_REC2020 * color);
 color = clamp((log2(max(color, 1e-10)) + 12.47393) / 16.5, 0.0, 1.0);
-color = ph * fT(color);
+color = qr * gj(color);
 color = pow(max(color, 0.0), vec3(2.2));
 return clamp(LINEAR_REC2020_TO_LINEAR_SRGB * color, 0.0, 1.0);
 }
-vec3 bZ(vec3 color) {
-if (toneMappingMode == 1) return lO(color);
-if (toneMappingMode == 2) return jP(color);
-if (toneMappingMode == 3) return iN(color);
-if (toneMappingMode == 4) return kb(color);
-return my(color);
+vec3 bH(vec3 color) {
+if (toneMappingMode == 1) return mH(color);
+if (toneMappingMode == 2) return kw(color);
+if (toneMappingMode == 3) return jp(color);
+if (toneMappingMode == 4) return kH(color);
+return nC(color);
 }
-vec3 ga(vec3 color) {
-vec3 uL = vec3(
+vec3 fm(vec3 color) {
+vec3 wu = vec3(
 1.0 + temperature * 0.16 + tint * 0.04,
 1.0 - tint * 0.12,
 1.0 - temperature * 0.16 + tint * 0.04
 );
-color = max(color * uL, 0.0);
+color = max(color * wu, 0.0);
 color = max((color - 0.18) * contrast + 0.18, 0.0);
 float luminance = dot(color, vec3(0.2126, 0.7152, 0.0722));
 return max(mix(vec3(luminance), color, saturation), 0.0);
@@ -146,12 +146,12 @@ return;
 vec4 color = sample_colorTex(uv);
 vec4 bloom = sample_bloomTex(uv) * bloomIntensity;
 color.rgb += bloom.rgb;
-color.rgb = ga(color.rgb * exposure);
+color.rgb = fm(color.rgb * exposure);
 if (sRGBOutput) {
 float edge = smoothstep(0.2, 0.72, distance(uv, vec2(0.5)));
 color.rgb *= 1.0 - vignette * edge;
-color.rgb = bZ(color.rgb);
-color.rgb = nT(color.rgb);
+color.rgb = bH(color.rgb);
+color.rgb = oS(color.rgb);
 }
 color.a = max(color.a, min(dot(bloom.rgb, vec3(0.2125, 0.7154, 0.0721)), 1.0));
 
@@ -159,32 +159,32 @@ if (shadowCatcherEnabled > 0.5) {
 float coverage = clamp(shadow.r, 0.0, 1.0);
 
 
-float fR = clamp(
+float gd = clamp(
 (1.0 - shadow.a) / max(coverage, 0.000001), 0.0, 1.0
 );
-fR = max(fR - (1.0 - shadowCatcherAlpha), 0.0);
+gd = max(gd - (1.0 - shadowCatcherAlpha), 0.0);
 if (coverage > 0.0) {
 vec4 foreground = sample_foregroundTex(uv);
 vec4 background = sample_catcherBackgroundTex(uv);
-vec3 foregroundColor = ga((foreground.rgb + bloom.rgb) * exposure);
-vec3 backgroundColor = ga((background.rgb / max(background.a, 0.000001) + bloom.rgb) * exposure);
+vec3 foregroundColor = fm((foreground.rgb + bloom.rgb) * exposure);
+vec3 backgroundColor = fm((background.rgb / max(background.a, 0.000001) + bloom.rgb) * exposure);
 if (sRGBOutput) {
 float edge = smoothstep(0.2, 0.72, distance(uv, vec2(0.5)));
 foregroundColor *= 1.0 - vignette * edge;
 backgroundColor *= 1.0 - vignette * edge;
-foregroundColor = nT(bZ(foregroundColor));
-backgroundColor = nT(bZ(backgroundColor));
+foregroundColor = oS(bH(foregroundColor));
+backgroundColor = oS(bH(backgroundColor));
 }
 
-float opacity = fR;
-fR *= coverage;
-float mD = clamp(dot(bloom.rgb, vec3(0.2125, 0.7154, 0.0721)), 0.0, 1.0);
-float hq = max(foreground.a, mD * (1.0 - coverage));
-float jf = max(background.a, mD * coverage);
-float alpha = hq + jf * (1.0 - opacity) + fR;
-vec3 premultiplied = foregroundColor * hq
-+ backgroundColor * jf * (1.0 - opacity)
-+ shadowCatcherColor * fR;
+float opacity = gd;
+gd *= coverage;
+float nF = clamp(dot(bloom.rgb, vec3(0.2125, 0.7154, 0.0721)), 0.0, 1.0);
+float hK = max(foreground.a, nF * (1.0 - coverage));
+float jJ = max(background.a, nF * coverage);
+float alpha = hK + jJ * (1.0 - opacity) + gd;
+vec3 premultiplied = foregroundColor * hK
++ backgroundColor * jJ * (1.0 - opacity)
++ shadowCatcherColor * gd;
 color = vec4(premultiplied / max(alpha, 0.000001), alpha);
 }
 }

package/lib/renderer/gpuMemory.js

@@ -3,7 +3,7 @@ export function textureBytes(texture) {
     if (!texture)
         return 0;
     const bytes = texture.format === 'rgba32float' ? 16 :
-        texture.format === 'rgba16float' ? 8 :
+        texture.format === 'rgba16float' || texture.format === 'rg32float' ? 8 :
             texture.format === 'r32float' ? 4 : 4;
     return texture.width * texture.height * texture.depthOrArrayLayers * bytes;
 }

package/lib/renderer/gpuMemory.js.map

@@ -1 +1 @@
-{"version":3,"file":"gpuMemory.js","sourceRoot":"","sources":["../../src/renderer/gpuMemory.ts"],"names":[],"mappings":"AAAA,gFAAgF;AAChF,MAAM,UAAU,YAAY,CAAC,OAA+B;IAC1D,IAAI,CAAC,OAAO;QAAE,OAAO,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,KAAK,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACnD,OAAO,CAAC,MAAM,KAAK,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACtC,OAAO,CAAC,MAAM,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,OAAO,OAAO,CAAC,KAAK,GAAG,OAAO,CAAC,MAAM,GAAG,OAAO,CAAC,kBAAkB,GAAG,KAAK,CAAC;AAC7E,CAAC","sourcesContent":["/** Payload bytes of the uncompressed render-target formats used by MuTracer. */\nexport function textureBytes(texture: GPUTexture | undefined): number {\n  if (!texture) return 0;\n  const bytes = texture.format === 'rgba32float' ? 16 :\n    texture.format === 'rgba16float' ? 8 :\n    texture.format === 'r32float' ? 4 : 4;\n  return texture.width * texture.height * texture.depthOrArrayLayers * bytes;\n}\n"]}
+{"version":3,"file":"gpuMemory.js","sourceRoot":"","sources":["../../src/renderer/gpuMemory.ts"],"names":[],"mappings":"AAAA,gFAAgF;AAChF,MAAM,UAAU,YAAY,CAAC,OAA+B;IAC1D,IAAI,CAAC,OAAO;QAAE,OAAO,CAAC,CAAC;IACvB,MAAM,KAAK,GAAG,OAAO,CAAC,MAAM,KAAK,aAAa,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QACnD,OAAO,CAAC,MAAM,KAAK,aAAa,IAAI,OAAO,CAAC,MAAM,KAAK,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACxE,OAAO,CAAC,MAAM,KAAK,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACxC,OAAO,OAAO,CAAC,KAAK,GAAG,OAAO,CAAC,MAAM,GAAG,OAAO,CAAC,kBAAkB,GAAG,KAAK,CAAC;AAC7E,CAAC","sourcesContent":["/** Payload bytes of the uncompressed render-target formats used by MuTracer. */\nexport function textureBytes(texture: GPUTexture | undefined): number {\n  if (!texture) return 0;\n  const bytes = texture.format === 'rgba32float' ? 16 :\n    texture.format === 'rgba16float' || texture.format === 'rg32float' ? 8 :\n    texture.format === 'r32float' ? 4 : 4;\n  return texture.width * texture.height * texture.depthOrArrayLayers * bytes;\n}\n"]}

package/lib/renderer/openPbrEnergyLut.d.ts

@@ -8,7 +8,8 @@
  * This repository re-quantized the Adobe data to 8-bit row-normalized values.
  * Directional rows store Uint16 min/max followed by 32 Uint8 samples.
  * Metal and ideal-dielectric averages store exact little-endian Uint16 values.
- * They decode directly to Float32 without row quantization or repeated texels.
+ * Dielectric averages decode directly to Float32. Metal averages are normalized
+ * to the decoded piecewise-linear directional rows (see below).
  * Storage uses differences from the previous row: little-endian min/max
  * modulo 65536 and samples modulo 256. The first row is relative to zero;
  * differences continue across all concatenated energy tables.

package/lib/renderer/openPbrEnergyLut.js

@@ -8,7 +8,8 @@
  * This repository re-quantized the Adobe data to 8-bit row-normalized values.
  * Directional rows store Uint16 min/max followed by 32 Uint8 samples.
  * Metal and ideal-dielectric averages store exact little-endian Uint16 values.
- * They decode directly to Float32 without row quantization or repeated texels.
+ * Dielectric averages decode directly to Float32. Metal averages are normalized
+ * to the decoded piecewise-linear directional rows (see below).
  * Storage uses differences from the previous row: little-endian min/max
  * modulo 65536 and samples modulo 256. The first row is relative to zero;
  * differences continue across all concatenated energy tables.
@@ -963,6 +964,22 @@ export function getOpenPbrEnergyLutData() {
             : OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW * 32 + index - 32;
         data[offset] = value;
     }
+    // Integrating before interpolation in roughness commutes with that linear
+    // interpolation. Use the exact integral of 2*mu times each linear segment,
+    // rather than the separately tabulated u16 average: low-roughness rows have
+    // narrow grazing features that the 32 cosine samples cannot resolve exactly.
+    // This makes the compensation lobe integrate to its directional loss for any
+    // roughness without changing the directional table or the Fresnel model.
+    for (let row = 0; row < OPENPBR_ENERGY_LUT_SIZE; row++) {
+        const offset = (OPENPBR_METAL_ENERGY_LUT_ROW + row) * OPENPBR_ENERGY_LUT_WIDTH;
+        const h = 1 / (OPENPBR_ENERGY_LUT_SIZE - 1);
+        let average = 0;
+        for (let i = 0; i < OPENPBR_ENERGY_LUT_SIZE - 1; i++) {
+            const a = data[offset + i], b = data[offset + i + 1];
+            average += h * (i * h * (a + b) + h * (a + 2 * b) / 3);
+        }
+        data[OPENPBR_METAL_ENERGY_LUT_AVERAGE_ROW * OPENPBR_ENERGY_LUT_WIDTH + row] = average;
+    }
     cachedData = data;
     return data;
 }