@zephyr3d/scene 0.9.27 → 0.10.0

package/dist/shaders/skin_brdf.js

@@ -0,0 +1,243 @@
+import { distributionGGX, visGGX, fresnelSchlick } from './pbr.js';
+import { SSSProfile, SKIN_TRANSMISSION_NO_DATA_ENCODING } from '../material/sssprofile.js';
+
+/** Opacity below which the dual-lobe specular fades back to a single lobe. */ const OPACITY_THRESHOLD_EPS = 0.1;
+/**
+ * Multiple-scattering energy terms for a GGX specular lobe.
+ *
+ * @param scope - Shader scope.
+ * @param roughness - Roughness of the lobe being corrected.
+ * @param NoV - Clamped dot(normal, viewDir).
+ * @param F0 - Fresnel reflectance at normal incidence.
+ * @returns `vec2(W, A)`: specular gain and reflected share.
+ *
+ * @internal
+ */ function skinSpecularEnergyTerms(scope, roughness, NoV, F0) {
+    const pb = scope.$builder;
+    const funcName = 'lib_skinSpecularEnergyTerms';
+    pb.func(funcName, [
+        pb.float('roughness'),
+        pb.float('NoV'),
+        pb.float('F0')
+    ], function() {
+        this.$l.r = pb.max(this.roughness, 1e-3);
+        this.$l.c = pb.max(this.NoV, 1e-3);
+        this.$l.E = pb.sub(1, pb.clamp(pb.mul(pb.pow(this.r, pb.div(this.c, this.r)), pb.div(pb.add(pb.mul(this.r, this.c), 0.0266916), pb.add(0.466495, this.c))), 0, 1));
+        this.E = pb.max(this.E, 1e-4);
+        this.$l.oneMinusC = pb.sub(1, this.c);
+        this.$l.pow5 = pb.mul(pb.mul(this.oneMinusC, this.oneMinusC), pb.mul(this.oneMinusC, this.oneMinusC), this.oneMinusC);
+        this.$l.Ef = pb.mul(this.pow5, pb.pow(pb.add(pb.mul(2.36651, pb.pow(this.c, pb.mul(4.7703, this.r))), 0.0387332), this.r));
+        this.$l.F90 = pb.clamp(pb.mul(50, this.F0), 0, 1);
+        this.$l.W = pb.add(1, pb.mul(this.F0, pb.div(pb.sub(1, this.E), this.E)));
+        this.$l.A = pb.mul(this.W, pb.add(pb.mul(this.E, this.F0), pb.mul(this.Ef, pb.sub(this.F90, this.F0))));
+        this.$return(pb.vec2(this.W, pb.clamp(this.A, 0, 1)));
+    });
+    return pb.getGlobalScope()[funcName](roughness, NoV, F0);
+}
+/**
+ * Burley diffuse BRDF.
+ *
+ * @param scope - Shader scope.
+ * @param NdotV - Clamped dot(normal, viewDir).
+ * @param NdotL - Clamped dot(normal, lightDir).
+ * @param VdotH - dot(viewDir, halfVector).
+ * @param roughness - Material roughness.
+ * @returns Scalar diffuse BRDF value; the caller applies albedo and `NoL`.
+ *
+ * @internal
+ */ function skinDiffuseBRDF(scope, NdotV, NdotL, VdotH, roughness) {
+    const pb = scope.$builder;
+    const funcName = 'lib_skinDiffuseBurley';
+    pb.func(funcName, [
+        pb.float('NdotV'),
+        pb.float('NdotL'),
+        pb.float('VdotH'),
+        pb.float('roughness')
+    ], function() {
+        this.$l.fd90 = pb.add(0.5, pb.mul(2, this.VdotH, this.VdotH, this.roughness));
+        this.$l.fd90m1 = pb.sub(this.fd90, 1);
+        this.$l.oneMinusNoV = pb.sub(1, this.NdotV);
+        this.$l.pow5V = pb.mul(pb.mul(this.oneMinusNoV, this.oneMinusNoV), pb.mul(this.oneMinusNoV, this.oneMinusNoV), this.oneMinusNoV);
+        this.$l.oneMinusNoL = pb.sub(1, this.NdotL);
+        this.$l.pow5L = pb.mul(pb.mul(this.oneMinusNoL, this.oneMinusNoL), pb.mul(this.oneMinusNoL, this.oneMinusNoL), this.oneMinusNoL);
+        this.$l.fdV = pb.add(1, pb.mul(this.fd90m1, this.pow5V));
+        this.$l.fdL = pb.add(1, pb.mul(this.fd90m1, this.pow5L));
+        this.$return(pb.mul(1 / Math.PI, this.fdV, this.fdL));
+    });
+    return pb.getGlobalScope()[funcName](NdotV, NdotL, VdotH, roughness);
+}
+/**
+ * Resolves the two specular lobe roughnesses from the profile.
+ *
+ * @param scope - Shader scope.
+ * @param roughness - Material roughness.
+ * @param opacity - Subsurface opacity, used to fade the effect out.
+ * @param lobe0Scale - Profile roughness multiplier for the narrow lobe.
+ * @param lobe1Scale - Profile roughness multiplier for the wide lobe.
+ * @returns `vec2(lobe0Roughness, lobe1Roughness)`.
+ *
+ * @internal
+ */ function skinDualSpecularRoughness(scope, roughness, opacity, lobe0Scale, lobe1Scale) {
+    const pb = scope.$builder;
+    const funcName = 'lib_skinDualSpecularRoughness';
+    pb.func(funcName, [
+        pb.float('roughness'),
+        pb.float('opacity'),
+        pb.float('lobe0Scale'),
+        pb.float('lobe1Scale')
+    ], function() {
+        this.$l.fade = pb.clamp(pb.mul(pb.sub(this.opacity, OPACITY_THRESHOLD_EPS), 10), 0, 1);
+        this.$l.scale0 = pb.mix(pb.float(1), this.lobe0Scale, this.fade);
+        this.$l.scale1 = pb.mix(pb.float(1), this.lobe1Scale, this.fade);
+        this.$return(pb.vec2(pb.max(pb.clamp(pb.mul(this.roughness, this.scale0), 0, 1), 0.02), pb.clamp(pb.mul(this.roughness, this.scale1), 0, 1)));
+    });
+    return pb.getGlobalScope()[funcName](roughness, opacity, lobe0Scale, lobe1Scale);
+}
+/**
+ * Dual-lobe GGX specular for skin.
+ *
+ * @param scope - Shader scope.
+ * @param NoH - dot(normal, halfVector).
+ * @param NoV - dot(normal, viewDir).
+ * @param NoL - dot(normal, lightDir).
+ * @param VoH - dot(viewDir, halfVector).
+ * @param lobeRoughness - `vec2` of the two lobe roughnesses.
+ * @param lobeMix - Blend between the narrow and wide lobe.
+ * @param F0 - Specular colour at normal incidence.
+ * @returns vec3 specular contribution; the caller applies light colour and `NoL`.
+ *
+ * @internal
+ */ function skinDualLobeSpecular(scope, NoH, NoV, NoL, VoH, lobeRoughness, lobeMix, F0) {
+    const pb = scope.$builder;
+    const funcName = 'lib_skinDualLobeSpecular';
+    pb.func(funcName, [
+        pb.float('NoH'),
+        pb.float('NoV'),
+        pb.float('NoL'),
+        pb.float('VoH'),
+        pb.vec2('lobeRoughness'),
+        pb.float('lobeMix'),
+        pb.float('F0')
+    ], function() {
+        this.$l.a0 = pb.mul(this.lobeRoughness.x, this.lobeRoughness.x);
+        this.$l.a1 = pb.mul(this.lobeRoughness.y, this.lobeRoughness.y);
+        this.$l.avgRoughness = pb.mix(this.lobeRoughness.x, this.lobeRoughness.y, this.lobeMix);
+        this.$l.avgAlpha = pb.mul(this.avgRoughness, this.avgRoughness);
+        this.$l.D = pb.mix(distributionGGX(this, this.NoH, this.a0), distributionGGX(this, this.NoH, this.a1), this.lobeMix);
+        this.$l.Vis = visGGX(this, this.NoV, this.NoL, this.avgAlpha);
+        this.$l.F = fresnelSchlick(this, this.VoH, pb.vec3(this.F0), pb.vec3(1));
+        this.$return(pb.mul(this.F, pb.mul(this.D, this.Vis)));
+    });
+    return pb.getGlobalScope()[funcName](NoH, NoV, NoL, VoH, lobeRoughness, lobeMix, F0);
+}
+/**
+ * Fetches one column of a profile row from the packed table.
+ *
+ * @remarks
+ * The table is `rgba32f`, hence `unfilterable-float` on WebGPU: the sampler must
+ * be non-filtering, so this is an exact texel fetch and anything wanting
+ * interpolation between columns has to lerp two fetches by hand.
+ *
+ * Emitted inline rather than behind a shader function: on WebGPU a texture
+ * parameter needs its sampler passed alongside it, and `textureSampleLevel` on
+ * the parameter alone compiles to nothing usable.
+ *
+ * @param scope - Shader scope.
+ * @param tex - The packed profile table, as declared in the caller's scope.
+ * @param texelSize - `vec2(1 / columns, 1 / rows)`.
+ * @param row - Precomputed `v` coordinate of the profile's row.
+ * @param column - Column to read, as a float.
+ * @returns The `vec4` stored at that column.
+ *
+ * @internal
+ */ function readSSSProfileColumn(scope, tex, texelSize, row, column) {
+    const pb = scope.$builder;
+    return pb.textureSampleLevel(tex, pb.vec2(pb.mul(pb.add(column, 0.5), texelSize.x), row), 0);
+}
+/**
+ * Back-lit subsurface transmission.
+ *
+ * @param scope - Shader scope.
+ * @param profileTex - The packed profile table.
+ * @param profileTexelSize - `vec2(1 / columns, 1 / rows)` of that table.
+ * @param profileId - Normalized profile id.
+ * @param thickness - Encoded optical depth from the thickness pass, where 1
+ *   means "nothing in the way".
+ * @param normal - Shading normal.
+ * @param viewVec - Unit vector from the surface towards the eye.
+ * @param lightDir - Unit vector from the surface towards the light.
+ * @returns vec3 transmitted radiance; the caller applies light colour,
+ *   attenuation and the transmission shadow.
+ *
+ * @internal
+ */ function skinTransmission(scope, profileTex, profileTexelSize, profileId, thickness, normal, viewVec, lightDir) {
+    const pb = scope.$builder;
+    const lutOffset = SSSProfile.transmissionLutOffset;
+    const lutSize = SSSProfile.transmissionLutSize;
+    const lastLutColumn = lutOffset + lutSize - 1;
+    // Table reads stay inline (see readSSSProfileColumn); the locals are prefixed
+    // so they cannot collide with the caller's. Rows are addressed by the profile
+    // id, which arrives normalized because it rides in an 8-bit channel.
+    scope.$l.zSkinTrRow = pb.mul(pb.add(pb.mul(pb.clamp(profileId, 0, 1), 255), 0.5), profileTexelSize.y);
+    // GetTransmissionProfile. The index is `opticalDepth / MAX * (size - 1)` and
+    // the decode is `opticalDepth = (1 - thickness) * MAX`, so the MAX cancels and
+    // the encoded thickness maps onto the table directly.
+    scope.$l.zSkinTrIndex = pb.mul(pb.clamp(pb.sub(1, thickness), 0, 1), lutSize - 1);
+    scope.$l.zSkinTrI0 = pb.floor(scope.zSkinTrIndex);
+    scope.$l.zSkinTrC0 = pb.add(scope.zSkinTrI0, lutOffset);
+    scope.$l.zSkinTrC1 = pb.min(pb.add(scope.zSkinTrC0, 1), lastLutColumn);
+    scope.$l.zSkinTrProfile = pb.mix(readSSSProfileColumn(scope, profileTex, profileTexelSize, scope.zSkinTrRow, scope.zSkinTrC0).rgb, readSSSProfileColumn(scope, profileTex, profileTexelSize, scope.zSkinTrRow, scope.zSkinTrC1).rgb, pb.sub(scope.zSkinTrIndex, scope.zSkinTrI0));
+    // (extinctionScale, normalScale, scatteringDistribution, 1 / ior)
+    scope.$l.zSkinTrParams = readSSSProfileColumn(scope, profileTex, profileTexelSize, scope.zSkinTrRow, pb.float(SSSProfile.transmissionParamColumn));
+    return skinTransmissionPhase(scope, scope.zSkinTrProfile, scope.zSkinTrParams, thickness, normal, viewVec, lightDir);
+}
+/**
+ * The part of the transmission BxDF that touches no textures.
+ *
+ * @remarks
+ * Split out so the math can live behind a shader function while the table reads
+ * stay inline at the call site, where the sampler is in scope.
+ *
+ * @param scope - Shader scope.
+ * @param profile - The interpolated transmission profile, `rgb`.
+ * @param params - `(extinctionScale, normalScale, scatteringDistribution, 1 / ior)`.
+ * @param thickness - The raw encoding, still needed to spot the "no data" sentinel.
+ * @param normal - Shading normal.
+ * @param viewVec - Unit vector from the surface towards the eye.
+ * @param lightDir - Unit vector from the surface towards the light.
+ * @returns vec3 transmitted radiance.
+ *
+ * @internal
+ */ function skinTransmissionPhase(scope, profile, params, thickness, normal, viewVec, lightDir) {
+    const pb = scope.$builder;
+    const funcName = 'lib_skinTransmissionPhase';
+    pb.func(funcName, [
+        pb.vec3('profile'),
+        pb.vec4('params'),
+        pb.float('thickness'),
+        pb.vec3('normal'),
+        pb.vec3('viewVec'),
+        pb.vec3('lightDir')
+    ], function() {
+        // "No data" must be rejected rather than decoded: it decodes to zero
+        // optical depth, the *most* transmissive entry, so letting it through lights
+        // the subject from behind with a light that was never behind it. The
+        // threshold sits in the gap between the largest encoding the pass can write
+        // and the sentinel, so it clips no real measurement.
+        this.$if(pb.greaterThan(this.thickness, 0.5 * (1 + SKIN_TRANSMISSION_NO_DATA_ENCODING)), function() {
+            this.$return(pb.vec3(0));
+        });
+        this.$l.refracV = pb.refract(this.viewVec, pb.neg(this.normal), this.params.w);
+        this.$l.cosJ = pb.dot(pb.neg(this.lightDir), this.refracV);
+        this.$l.g = this.params.z;
+        this.$l.g2 = pb.mul(this.g, this.g);
+        this.$l.gcos = pb.sub(1, pb.mul(this.g, this.cosJ));
+        this.$l.gcos2 = pb.mul(this.gcos, this.gcos);
+        this.$l.phase = pb.div(pb.mul(0.5, pb.sub(1, this.g2)), pb.max(this.gcos2, 1e-5));
+        this.$return(pb.mul(this.profile, this.phase));
+    });
+    return pb.getGlobalScope()[funcName](profile, params, thickness, normal, viewVec, lightDir);
+}
+
+export { skinDiffuseBRDF, skinDualLobeSpecular, skinDualSpecularRoughness, skinSpecularEnergyTerms, skinTransmission };
+//# sourceMappingURL=skin_brdf.js.map

package/dist/shaders/skin_brdf.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"skin_brdf.js","sources":["../../src/shaders/skin_brdf.ts"],"sourcesContent":["import type { PBInsideFunctionScope, PBShaderExp } from '@zephyr3d/device';\r\nimport { distributionGGX, fresnelSchlick, visGGX } from './pbr';\r\nimport { SKIN_TRANSMISSION_NO_DATA_ENCODING, SSSProfile } from '../material/sssprofile';\r\n\r\n/** Opacity below which the dual-lobe specular fades back to a single lobe. */\r\nconst OPACITY_THRESHOLD_EPS = 0.1;\r\n\r\n/**\r\n * Multiple-scattering energy terms for a GGX specular lobe.\r\n *\r\n * @param scope - Shader scope.\r\n * @param roughness - Roughness of the lobe being corrected.\r\n * @param NoV - Clamped dot(normal, viewDir).\r\n * @param F0 - Fresnel reflectance at normal incidence.\r\n * @returns `vec2(W, A)`: specular gain and reflected share.\r\n *\r\n * @internal\r\n */\r\nexport function skinSpecularEnergyTerms(\r\n  scope: PBInsideFunctionScope,\r\n  roughness: PBShaderExp,\r\n  NoV: PBShaderExp,\r\n  F0: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const funcName = 'lib_skinSpecularEnergyTerms';\r\n  pb.func(funcName, [pb.float('roughness'), pb.float('NoV'), pb.float('F0')], function () {\r\n    this.$l.r = pb.max(this.roughness, 1e-3);\r\n    this.$l.c = pb.max(this.NoV, 1e-3);\r\n    this.$l.E = pb.sub(\r\n      1,\r\n      pb.clamp(\r\n        pb.mul(\r\n          pb.pow(this.r, pb.div(this.c, this.r)),\r\n          pb.div(pb.add(pb.mul(this.r, this.c), 0.0266916), pb.add(0.466495, this.c))\r\n        ),\r\n        0,\r\n        1\r\n      )\r\n    );\r\n    this.E = pb.max(this.E, 1e-4);\r\n    this.$l.oneMinusC = pb.sub(1, this.c);\r\n    this.$l.pow5 = pb.mul(\r\n      pb.mul(this.oneMinusC, this.oneMinusC),\r\n      pb.mul(this.oneMinusC, this.oneMinusC),\r\n      this.oneMinusC\r\n    );\r\n    this.$l.Ef = pb.mul(\r\n      this.pow5,\r\n      pb.pow(pb.add(pb.mul(2.36651, pb.pow(this.c, pb.mul(4.7703, this.r))), 0.0387332), this.r)\r\n    );\r\n    this.$l.F90 = pb.clamp(pb.mul(50, this.F0), 0, 1);\r\n    this.$l.W = pb.add(1, pb.mul(this.F0, pb.div(pb.sub(1, this.E), this.E)));\r\n    this.$l.A = pb.mul(this.W, pb.add(pb.mul(this.E, this.F0), pb.mul(this.Ef, pb.sub(this.F90, this.F0))));\r\n    this.$return(pb.vec2(this.W, pb.clamp(this.A, 0, 1)));\r\n  });\r\n  return pb.getGlobalScope()[funcName](roughness, NoV, F0) as PBShaderExp;\r\n}\r\n\r\n/**\r\n * Burley diffuse BRDF.\r\n *\r\n * @param scope - Shader scope.\r\n * @param NdotV - Clamped dot(normal, viewDir).\r\n * @param NdotL - Clamped dot(normal, lightDir).\r\n * @param VdotH - dot(viewDir, halfVector).\r\n * @param roughness - Material roughness.\r\n * @returns Scalar diffuse BRDF value; the caller applies albedo and `NoL`.\r\n *\r\n * @internal\r\n */\r\nexport function skinDiffuseBRDF(\r\n  scope: PBInsideFunctionScope,\r\n  NdotV: PBShaderExp,\r\n  NdotL: PBShaderExp,\r\n  VdotH: PBShaderExp,\r\n  roughness: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const funcName = 'lib_skinDiffuseBurley';\r\n  pb.func(\r\n    funcName,\r\n    [pb.float('NdotV'), pb.float('NdotL'), pb.float('VdotH'), pb.float('roughness')],\r\n    function () {\r\n      this.$l.fd90 = pb.add(0.5, pb.mul(2, this.VdotH, this.VdotH, this.roughness));\r\n      this.$l.fd90m1 = pb.sub(this.fd90, 1);\r\n      this.$l.oneMinusNoV = pb.sub(1, this.NdotV);\r\n      this.$l.pow5V = pb.mul(\r\n        pb.mul(this.oneMinusNoV, this.oneMinusNoV),\r\n        pb.mul(this.oneMinusNoV, this.oneMinusNoV),\r\n        this.oneMinusNoV\r\n      );\r\n      this.$l.oneMinusNoL = pb.sub(1, this.NdotL);\r\n      this.$l.pow5L = pb.mul(\r\n        pb.mul(this.oneMinusNoL, this.oneMinusNoL),\r\n        pb.mul(this.oneMinusNoL, this.oneMinusNoL),\r\n        this.oneMinusNoL\r\n      );\r\n      this.$l.fdV = pb.add(1, pb.mul(this.fd90m1, this.pow5V));\r\n      this.$l.fdL = pb.add(1, pb.mul(this.fd90m1, this.pow5L));\r\n      this.$return(pb.mul(1 / Math.PI, this.fdV, this.fdL));\r\n    }\r\n  );\r\n  return pb.getGlobalScope()[funcName](NdotV, NdotL, VdotH, roughness) as PBShaderExp;\r\n}\r\n\r\n/**\r\n * Resolves the two specular lobe roughnesses from the profile.\r\n *\r\n * @param scope - Shader scope.\r\n * @param roughness - Material roughness.\r\n * @param opacity - Subsurface opacity, used to fade the effect out.\r\n * @param lobe0Scale - Profile roughness multiplier for the narrow lobe.\r\n * @param lobe1Scale - Profile roughness multiplier for the wide lobe.\r\n * @returns `vec2(lobe0Roughness, lobe1Roughness)`.\r\n *\r\n * @internal\r\n */\r\nexport function skinDualSpecularRoughness(\r\n  scope: PBInsideFunctionScope,\r\n  roughness: PBShaderExp,\r\n  opacity: PBShaderExp,\r\n  lobe0Scale: PBShaderExp,\r\n  lobe1Scale: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const funcName = 'lib_skinDualSpecularRoughness';\r\n  pb.func(\r\n    funcName,\r\n    [pb.float('roughness'), pb.float('opacity'), pb.float('lobe0Scale'), pb.float('lobe1Scale')],\r\n    function () {\r\n      this.$l.fade = pb.clamp(pb.mul(pb.sub(this.opacity, OPACITY_THRESHOLD_EPS), 10), 0, 1);\r\n      this.$l.scale0 = pb.mix(pb.float(1), this.lobe0Scale, this.fade);\r\n      this.$l.scale1 = pb.mix(pb.float(1), this.lobe1Scale, this.fade);\r\n      this.$return(\r\n        pb.vec2(\r\n          pb.max(pb.clamp(pb.mul(this.roughness, this.scale0), 0, 1), 0.02),\r\n          pb.clamp(pb.mul(this.roughness, this.scale1), 0, 1)\r\n        )\r\n      );\r\n    }\r\n  );\r\n  return pb.getGlobalScope()[funcName](roughness, opacity, lobe0Scale, lobe1Scale) as PBShaderExp;\r\n}\r\n\r\n/**\r\n * Dual-lobe GGX specular for skin.\r\n *\r\n * @param scope - Shader scope.\r\n * @param NoH - dot(normal, halfVector).\r\n * @param NoV - dot(normal, viewDir).\r\n * @param NoL - dot(normal, lightDir).\r\n * @param VoH - dot(viewDir, halfVector).\r\n * @param lobeRoughness - `vec2` of the two lobe roughnesses.\r\n * @param lobeMix - Blend between the narrow and wide lobe.\r\n * @param F0 - Specular colour at normal incidence.\r\n * @returns vec3 specular contribution; the caller applies light colour and `NoL`.\r\n *\r\n * @internal\r\n */\r\nexport function skinDualLobeSpecular(\r\n  scope: PBInsideFunctionScope,\r\n  NoH: PBShaderExp,\r\n  NoV: PBShaderExp,\r\n  NoL: PBShaderExp,\r\n  VoH: PBShaderExp,\r\n  lobeRoughness: PBShaderExp,\r\n  lobeMix: PBShaderExp,\r\n  F0: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const funcName = 'lib_skinDualLobeSpecular';\r\n  pb.func(\r\n    funcName,\r\n    [\r\n      pb.float('NoH'),\r\n      pb.float('NoV'),\r\n      pb.float('NoL'),\r\n      pb.float('VoH'),\r\n      pb.vec2('lobeRoughness'),\r\n      pb.float('lobeMix'),\r\n      pb.float('F0')\r\n    ],\r\n    function () {\r\n      this.$l.a0 = pb.mul(this.lobeRoughness.x, this.lobeRoughness.x);\r\n      this.$l.a1 = pb.mul(this.lobeRoughness.y, this.lobeRoughness.y);\r\n      this.$l.avgRoughness = pb.mix(this.lobeRoughness.x, this.lobeRoughness.y, this.lobeMix);\r\n      this.$l.avgAlpha = pb.mul(this.avgRoughness, this.avgRoughness);\r\n      this.$l.D = pb.mix(\r\n        distributionGGX(this, this.NoH, this.a0),\r\n        distributionGGX(this, this.NoH, this.a1),\r\n        this.lobeMix\r\n      );\r\n      this.$l.Vis = visGGX(this, this.NoV, this.NoL, this.avgAlpha);\r\n      this.$l.F = fresnelSchlick(this, this.VoH, pb.vec3(this.F0), pb.vec3(1));\r\n      this.$return(pb.mul(this.F, pb.mul(this.D, this.Vis)));\r\n    }\r\n  );\r\n  return pb.getGlobalScope()[funcName](NoH, NoV, NoL, VoH, lobeRoughness, lobeMix, F0) as PBShaderExp;\r\n}\r\n\r\n/**\r\n * Fetches one column of a profile row from the packed table.\r\n *\r\n * @remarks\r\n * The table is `rgba32f`, hence `unfilterable-float` on WebGPU: the sampler must\r\n * be non-filtering, so this is an exact texel fetch and anything wanting\r\n * interpolation between columns has to lerp two fetches by hand.\r\n *\r\n * Emitted inline rather than behind a shader function: on WebGPU a texture\r\n * parameter needs its sampler passed alongside it, and `textureSampleLevel` on\r\n * the parameter alone compiles to nothing usable.\r\n *\r\n * @param scope - Shader scope.\r\n * @param tex - The packed profile table, as declared in the caller's scope.\r\n * @param texelSize - `vec2(1 / columns, 1 / rows)`.\r\n * @param row - Precomputed `v` coordinate of the profile's row.\r\n * @param column - Column to read, as a float.\r\n * @returns The `vec4` stored at that column.\r\n *\r\n * @internal\r\n */\r\nfunction readSSSProfileColumn(\r\n  scope: PBInsideFunctionScope,\r\n  tex: PBShaderExp,\r\n  texelSize: PBShaderExp,\r\n  row: PBShaderExp,\r\n  column: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  return pb.textureSampleLevel(tex, pb.vec2(pb.mul(pb.add(column, 0.5), texelSize.x), row), 0);\r\n}\r\n\r\n/**\r\n * Back-lit subsurface transmission.\r\n *\r\n * @param scope - Shader scope.\r\n * @param profileTex - The packed profile table.\r\n * @param profileTexelSize - `vec2(1 / columns, 1 / rows)` of that table.\r\n * @param profileId - Normalized profile id.\r\n * @param thickness - Encoded optical depth from the thickness pass, where 1\r\n *   means \"nothing in the way\".\r\n * @param normal - Shading normal.\r\n * @param viewVec - Unit vector from the surface towards the eye.\r\n * @param lightDir - Unit vector from the surface towards the light.\r\n * @returns vec3 transmitted radiance; the caller applies light colour,\r\n *   attenuation and the transmission shadow.\r\n *\r\n * @internal\r\n */\r\nexport function skinTransmission(\r\n  scope: PBInsideFunctionScope,\r\n  profileTex: PBShaderExp,\r\n  profileTexelSize: PBShaderExp,\r\n  profileId: PBShaderExp,\r\n  thickness: PBShaderExp,\r\n  normal: PBShaderExp,\r\n  viewVec: PBShaderExp,\r\n  lightDir: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const lutOffset = SSSProfile.transmissionLutOffset;\r\n  const lutSize = SSSProfile.transmissionLutSize;\r\n  const lastLutColumn = lutOffset + lutSize - 1;\r\n  // Table reads stay inline (see readSSSProfileColumn); the locals are prefixed\r\n  // so they cannot collide with the caller's. Rows are addressed by the profile\r\n  // id, which arrives normalized because it rides in an 8-bit channel.\r\n  scope.$l.zSkinTrRow = pb.mul(pb.add(pb.mul(pb.clamp(profileId, 0, 1), 255), 0.5), profileTexelSize.y);\r\n  // GetTransmissionProfile. The index is `opticalDepth / MAX * (size - 1)` and\r\n  // the decode is `opticalDepth = (1 - thickness) * MAX`, so the MAX cancels and\r\n  // the encoded thickness maps onto the table directly.\r\n  scope.$l.zSkinTrIndex = pb.mul(pb.clamp(pb.sub(1, thickness), 0, 1), lutSize - 1);\r\n  scope.$l.zSkinTrI0 = pb.floor(scope.zSkinTrIndex);\r\n  scope.$l.zSkinTrC0 = pb.add(scope.zSkinTrI0, lutOffset);\r\n  scope.$l.zSkinTrC1 = pb.min(pb.add(scope.zSkinTrC0, 1), lastLutColumn);\r\n  scope.$l.zSkinTrProfile = pb.mix(\r\n    readSSSProfileColumn(scope, profileTex, profileTexelSize, scope.zSkinTrRow, scope.zSkinTrC0).rgb,\r\n    readSSSProfileColumn(scope, profileTex, profileTexelSize, scope.zSkinTrRow, scope.zSkinTrC1).rgb,\r\n    pb.sub(scope.zSkinTrIndex, scope.zSkinTrI0)\r\n  );\r\n  // (extinctionScale, normalScale, scatteringDistribution, 1 / ior)\r\n  scope.$l.zSkinTrParams = readSSSProfileColumn(\r\n    scope,\r\n    profileTex,\r\n    profileTexelSize,\r\n    scope.zSkinTrRow,\r\n    pb.float(SSSProfile.transmissionParamColumn)\r\n  );\r\n  return skinTransmissionPhase(\r\n    scope,\r\n    scope.zSkinTrProfile,\r\n    scope.zSkinTrParams,\r\n    thickness,\r\n    normal,\r\n    viewVec,\r\n    lightDir\r\n  );\r\n}\r\n\r\n/**\r\n * The part of the transmission BxDF that touches no textures.\r\n *\r\n * @remarks\r\n * Split out so the math can live behind a shader function while the table reads\r\n * stay inline at the call site, where the sampler is in scope.\r\n *\r\n * @param scope - Shader scope.\r\n * @param profile - The interpolated transmission profile, `rgb`.\r\n * @param params - `(extinctionScale, normalScale, scatteringDistribution, 1 / ior)`.\r\n * @param thickness - The raw encoding, still needed to spot the \"no data\" sentinel.\r\n * @param normal - Shading normal.\r\n * @param viewVec - Unit vector from the surface towards the eye.\r\n * @param lightDir - Unit vector from the surface towards the light.\r\n * @returns vec3 transmitted radiance.\r\n *\r\n * @internal\r\n */\r\nfunction skinTransmissionPhase(\r\n  scope: PBInsideFunctionScope,\r\n  profile: PBShaderExp,\r\n  params: PBShaderExp,\r\n  thickness: PBShaderExp,\r\n  normal: PBShaderExp,\r\n  viewVec: PBShaderExp,\r\n  lightDir: PBShaderExp\r\n): PBShaderExp {\r\n  const pb = scope.$builder;\r\n  const funcName = 'lib_skinTransmissionPhase';\r\n  pb.func(\r\n    funcName,\r\n    [\r\n      pb.vec3('profile'),\r\n      pb.vec4('params'),\r\n      pb.float('thickness'),\r\n      pb.vec3('normal'),\r\n      pb.vec3('viewVec'),\r\n      pb.vec3('lightDir')\r\n    ],\r\n    function () {\r\n      // \"No data\" must be rejected rather than decoded: it decodes to zero\r\n      // optical depth, the *most* transmissive entry, so letting it through lights\r\n      // the subject from behind with a light that was never behind it. The\r\n      // threshold sits in the gap between the largest encoding the pass can write\r\n      // and the sentinel, so it clips no real measurement.\r\n      this.$if(pb.greaterThan(this.thickness, 0.5 * (1 + SKIN_TRANSMISSION_NO_DATA_ENCODING)), function () {\r\n        this.$return(pb.vec3(0));\r\n      });\r\n      this.$l.refracV = pb.refract(this.viewVec, pb.neg(this.normal), this.params.w);\r\n      this.$l.cosJ = pb.dot(pb.neg(this.lightDir), this.refracV);\r\n      this.$l.g = this.params.z;\r\n      this.$l.g2 = pb.mul(this.g, this.g);\r\n      this.$l.gcos = pb.sub(1, pb.mul(this.g, this.cosJ));\r\n      this.$l.gcos2 = pb.mul(this.gcos, this.gcos);\r\n      this.$l.phase = pb.div(pb.mul(0.5, pb.sub(1, this.g2)), pb.max(this.gcos2, 1e-5));\r\n      this.$return(pb.mul(this.profile, this.phase));\r\n    }\r\n  );\r\n  return pb.getGlobalScope()[funcName](profile, params, thickness, normal, viewVec, lightDir) as PBShaderExp;\r\n}\r\n"],"names":["OPACITY_THRESHOLD_EPS","skinSpecularEnergyTerms","scope","roughness","NoV","F0","pb","$builder","funcName","func","float","$l","r","max","c","E","sub","clamp","mul","pow","div","add","oneMinusC","pow5","Ef","F90","W","A","$return","vec2","getGlobalScope","skinDiffuseBRDF","NdotV","NdotL","VdotH","fd90","fd90m1","oneMinusNoV","pow5V","oneMinusNoL","pow5L","fdV","fdL","Math","PI","skinDualSpecularRoughness","opacity","lobe0Scale","lobe1Scale","fade","scale0","mix","scale1","skinDualLobeSpecular","NoH","NoL","VoH","lobeRoughness","lobeMix","a0","x","a1","y","avgRoughness","avgAlpha","D","distributionGGX","Vis","visGGX","F","fresnelSchlick","vec3","readSSSProfileColumn","tex","texelSize","row","column","textureSampleLevel","skinTransmission","profileTex","profileTexelSize","profileId","thickness","normal","viewVec","lightDir","lutOffset","SSSProfile","transmissionLutOffset","lutSize","transmissionLutSize","lastLutColumn","zSkinTrRow","zSkinTrIndex","zSkinTrI0","floor","zSkinTrC0","zSkinTrC1","min","zSkinTrProfile","rgb","zSkinTrParams","transmissionParamColumn","skinTransmissionPhase","profile","params","vec4","$if","greaterThan","SKIN_TRANSMISSION_NO_DATA_ENCODING","refracV","refract","neg","w","cosJ","dot","g","z","g2","gcos","gcos2","phase"],"mappings":";;;AAIA,+EACA,MAAMA,qBAAwB,GAAA,GAAA;AAE9B;;;;;;;;;;IAWO,SAASC,uBACdC,CAAAA,KAA4B,EAC5BC,SAAsB,EACtBC,GAAgB,EAChBC,EAAe,EAAA;IAEf,MAAMC,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,6BAAA;IACjBF,EAAGG,CAAAA,IAAI,CAACD,QAAU,EAAA;AAACF,QAAAA,EAAAA,CAAGI,KAAK,CAAC,WAAA,CAAA;AAAcJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,KAAA,CAAA;AAAQJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,IAAA;KAAM,EAAE,WAAA;QAC1E,IAAI,CAACC,EAAE,CAACC,CAAC,GAAGN,EAAGO,CAAAA,GAAG,CAAC,IAAI,CAACV,SAAS,EAAE,IAAA,CAAA;QACnC,IAAI,CAACQ,EAAE,CAACG,CAAC,GAAGR,EAAGO,CAAAA,GAAG,CAAC,IAAI,CAACT,GAAG,EAAE,IAAA,CAAA;QAC7B,IAAI,CAACO,EAAE,CAACI,CAAC,GAAGT,EAAGU,CAAAA,GAAG,CAChB,CAAA,EACAV,EAAGW,CAAAA,KAAK,CACNX,EAAAA,CAAGY,GAAG,CACJZ,EAAGa,CAAAA,GAAG,CAAC,IAAI,CAACP,CAAC,EAAEN,EAAAA,CAAGc,GAAG,CAAC,IAAI,CAACN,CAAC,EAAE,IAAI,CAACF,CAAC,CAAA,CAAA,EACpCN,EAAGc,CAAAA,GAAG,CAACd,EAAAA,CAAGe,GAAG,CAACf,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACN,CAAC,EAAE,IAAI,CAACE,CAAC,CAAG,EAAA,SAAA,CAAA,EAAYR,EAAGe,CAAAA,GAAG,CAAC,QAAA,EAAU,IAAI,CAACP,CAAC,CAAA,CAAA,CAAA,EAE3E,CACA,EAAA,CAAA,CAAA,CAAA;QAGJ,IAAI,CAACC,CAAC,GAAGT,EAAAA,CAAGO,GAAG,CAAC,IAAI,CAACE,CAAC,EAAE,IAAA,CAAA;QACxB,IAAI,CAACJ,EAAE,CAACW,SAAS,GAAGhB,EAAGU,CAAAA,GAAG,CAAC,CAAA,EAAG,IAAI,CAACF,CAAC,CAAA;AACpC,QAAA,IAAI,CAACH,EAAE,CAACY,IAAI,GAAGjB,EAAGY,CAAAA,GAAG,CACnBZ,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACI,SAAS,EAAE,IAAI,CAACA,SAAS,CACrChB,EAAAA,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACI,SAAS,EAAE,IAAI,CAACA,SAAS,CACrC,EAAA,IAAI,CAACA,SAAS,CAAA;AAEhB,QAAA,IAAI,CAACX,EAAE,CAACa,EAAE,GAAGlB,EAAAA,CAAGY,GAAG,CACjB,IAAI,CAACK,IAAI,EACTjB,EAAGa,CAAAA,GAAG,CAACb,EAAGe,CAAAA,GAAG,CAACf,EAAGY,CAAAA,GAAG,CAAC,OAAA,EAASZ,GAAGa,GAAG,CAAC,IAAI,CAACL,CAAC,EAAER,EAAGY,CAAAA,GAAG,CAAC,MAAA,EAAQ,IAAI,CAACN,CAAC,KAAK,SAAY,CAAA,EAAA,IAAI,CAACA,CAAC,CAAA,CAAA;AAE3F,QAAA,IAAI,CAACD,EAAE,CAACc,GAAG,GAAGnB,GAAGW,KAAK,CAACX,EAAGY,CAAAA,GAAG,CAAC,EAAI,EAAA,IAAI,CAACb,EAAE,GAAG,CAAG,EAAA,CAAA,CAAA;AAC/C,QAAA,IAAI,CAACM,EAAE,CAACe,CAAC,GAAGpB,EAAAA,CAAGe,GAAG,CAAC,CAAGf,EAAAA,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACb,EAAE,EAAEC,EAAGc,CAAAA,GAAG,CAACd,EAAAA,CAAGU,GAAG,CAAC,CAAG,EAAA,IAAI,CAACD,CAAC,CAAG,EAAA,IAAI,CAACA,CAAC,CAAA,CAAA,CAAA;AACtE,QAAA,IAAI,CAACJ,EAAE,CAACgB,CAAC,GAAGrB,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACQ,CAAC,EAAEpB,GAAGe,GAAG,CAACf,GAAGY,GAAG,CAAC,IAAI,CAACH,CAAC,EAAE,IAAI,CAACV,EAAE,CAAA,EAAGC,GAAGY,GAAG,CAAC,IAAI,CAACM,EAAE,EAAElB,EAAGU,CAAAA,GAAG,CAAC,IAAI,CAACS,GAAG,EAAE,IAAI,CAACpB,EAAE,CAAA,CAAA,CAAA,CAAA;AACnG,QAAA,IAAI,CAACuB,OAAO,CAACtB,EAAGuB,CAAAA,IAAI,CAAC,IAAI,CAACH,CAAC,EAAEpB,GAAGW,KAAK,CAAC,IAAI,CAACU,CAAC,EAAE,CAAG,EAAA,CAAA,CAAA,CAAA,CAAA;AACnD,KAAA,CAAA;AACA,IAAA,OAAOrB,GAAGwB,cAAc,EAAE,CAACtB,QAAS,CAAA,CAACL,WAAWC,GAAKC,EAAAA,EAAAA,CAAAA;AACvD;AAEA;;;;;;;;;;;IAYO,SAAS0B,eAAAA,CACd7B,KAA4B,EAC5B8B,KAAkB,EAClBC,KAAkB,EAClBC,KAAkB,EAClB/B,SAAsB,EAAA;IAEtB,MAAMG,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,uBAAA;IACjBF,EAAGG,CAAAA,IAAI,CACLD,QACA,EAAA;AAACF,QAAAA,EAAAA,CAAGI,KAAK,CAAC,OAAA,CAAA;AAAUJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,OAAA,CAAA;AAAUJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,OAAA,CAAA;AAAUJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,WAAA;KAAa,EAChF,WAAA;QACE,IAAI,CAACC,EAAE,CAACwB,IAAI,GAAG7B,EAAGe,CAAAA,GAAG,CAAC,GAAA,EAAKf,EAAGY,CAAAA,GAAG,CAAC,CAAG,EAAA,IAAI,CAACgB,KAAK,EAAE,IAAI,CAACA,KAAK,EAAE,IAAI,CAAC/B,SAAS,CAAA,CAAA;QAC3E,IAAI,CAACQ,EAAE,CAACyB,MAAM,GAAG9B,EAAGU,CAAAA,GAAG,CAAC,IAAI,CAACmB,IAAI,EAAE,CAAA,CAAA;QACnC,IAAI,CAACxB,EAAE,CAAC0B,WAAW,GAAG/B,EAAGU,CAAAA,GAAG,CAAC,CAAA,EAAG,IAAI,CAACgB,KAAK,CAAA;AAC1C,QAAA,IAAI,CAACrB,EAAE,CAAC2B,KAAK,GAAGhC,EAAGY,CAAAA,GAAG,CACpBZ,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACmB,WAAW,EAAE,IAAI,CAACA,WAAW,CACzC/B,EAAAA,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACmB,WAAW,EAAE,IAAI,CAACA,WAAW,CACzC,EAAA,IAAI,CAACA,WAAW,CAAA;QAElB,IAAI,CAAC1B,EAAE,CAAC4B,WAAW,GAAGjC,EAAGU,CAAAA,GAAG,CAAC,CAAA,EAAG,IAAI,CAACiB,KAAK,CAAA;AAC1C,QAAA,IAAI,CAACtB,EAAE,CAAC6B,KAAK,GAAGlC,EAAGY,CAAAA,GAAG,CACpBZ,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACqB,WAAW,EAAE,IAAI,CAACA,WAAW,CACzCjC,EAAAA,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACqB,WAAW,EAAE,IAAI,CAACA,WAAW,CACzC,EAAA,IAAI,CAACA,WAAW,CAAA;AAElB,QAAA,IAAI,CAAC5B,EAAE,CAAC8B,GAAG,GAAGnC,EAAAA,CAAGe,GAAG,CAAC,CAAA,EAAGf,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACkB,MAAM,EAAE,IAAI,CAACE,KAAK,CAAA,CAAA;AACtD,QAAA,IAAI,CAAC3B,EAAE,CAAC+B,GAAG,GAAGpC,EAAAA,CAAGe,GAAG,CAAC,CAAA,EAAGf,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACkB,MAAM,EAAE,IAAI,CAACI,KAAK,CAAA,CAAA;AACtD,QAAA,IAAI,CAACZ,OAAO,CAACtB,EAAGY,CAAAA,GAAG,CAAC,CAAIyB,GAAAA,IAAAA,CAAKC,EAAE,EAAE,IAAI,CAACH,GAAG,EAAE,IAAI,CAACC,GAAG,CAAA,CAAA;AACrD,KAAA,CAAA;IAEF,OAAOpC,EAAAA,CAAGwB,cAAc,EAAE,CAACtB,SAAS,CAACwB,KAAAA,EAAOC,OAAOC,KAAO/B,EAAAA,SAAAA,CAAAA;AAC5D;AAEA;;;;;;;;;;;IAYO,SAAS0C,yBAAAA,CACd3C,KAA4B,EAC5BC,SAAsB,EACtB2C,OAAoB,EACpBC,UAAuB,EACvBC,UAAuB,EAAA;IAEvB,MAAM1C,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,+BAAA;IACjBF,EAAGG,CAAAA,IAAI,CACLD,QACA,EAAA;AAACF,QAAAA,EAAAA,CAAGI,KAAK,CAAC,WAAA,CAAA;AAAcJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,SAAA,CAAA;AAAYJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,YAAA,CAAA;AAAeJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,YAAA;KAAc,EAC5F,WAAA;QACE,IAAI,CAACC,EAAE,CAACsC,IAAI,GAAG3C,EAAGW,CAAAA,KAAK,CAACX,EAAGY,CAAAA,GAAG,CAACZ,EAAGU,CAAAA,GAAG,CAAC,IAAI,CAAC8B,OAAO,EAAE9C,qBAAAA,CAAAA,EAAwB,KAAK,CAAG,EAAA,CAAA,CAAA;AACpF,QAAA,IAAI,CAACW,EAAE,CAACuC,MAAM,GAAG5C,EAAAA,CAAG6C,GAAG,CAAC7C,EAAAA,CAAGI,KAAK,CAAC,IAAI,IAAI,CAACqC,UAAU,EAAE,IAAI,CAACE,IAAI,CAAA;AAC/D,QAAA,IAAI,CAACtC,EAAE,CAACyC,MAAM,GAAG9C,EAAAA,CAAG6C,GAAG,CAAC7C,EAAAA,CAAGI,KAAK,CAAC,IAAI,IAAI,CAACsC,UAAU,EAAE,IAAI,CAACC,IAAI,CAAA;AAC/D,QAAA,IAAI,CAACrB,OAAO,CACVtB,GAAGuB,IAAI,CACLvB,GAAGO,GAAG,CAACP,EAAGW,CAAAA,KAAK,CAACX,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACf,SAAS,EAAE,IAAI,CAAC+C,MAAM,GAAG,CAAG,EAAA,CAAA,CAAA,EAAI,OAC5D5C,EAAGW,CAAAA,KAAK,CAACX,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACf,SAAS,EAAE,IAAI,CAACiD,MAAM,GAAG,CAAG,EAAA,CAAA,CAAA,CAAA,CAAA;AAGvD,KAAA,CAAA;IAEF,OAAO9C,EAAAA,CAAGwB,cAAc,EAAE,CAACtB,SAAS,CAACL,SAAAA,EAAW2C,SAASC,UAAYC,EAAAA,UAAAA,CAAAA;AACvE;AAEA;;;;;;;;;;;;;;AAcC,IACM,SAASK,oBAAAA,CACdnD,KAA4B,EAC5BoD,GAAgB,EAChBlD,GAAgB,EAChBmD,GAAgB,EAChBC,GAAgB,EAChBC,aAA0B,EAC1BC,OAAoB,EACpBrD,EAAe,EAAA;IAEf,MAAMC,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,0BAAA;IACjBF,EAAGG,CAAAA,IAAI,CACLD,QACA,EAAA;AACEF,QAAAA,EAAAA,CAAGI,KAAK,CAAC,KAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,KAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,KAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,KAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGuB,IAAI,CAAC,eAAA,CAAA;AACRvB,QAAAA,EAAAA,CAAGI,KAAK,CAAC,SAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGI,KAAK,CAAC,IAAA;KACV,EACD,WAAA;AACE,QAAA,IAAI,CAACC,EAAE,CAACgD,EAAE,GAAGrD,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACuC,aAAa,CAACG,CAAC,EAAE,IAAI,CAACH,aAAa,CAACG,CAAC,CAAA;AAC9D,QAAA,IAAI,CAACjD,EAAE,CAACkD,EAAE,GAAGvD,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACuC,aAAa,CAACK,CAAC,EAAE,IAAI,CAACL,aAAa,CAACK,CAAC,CAAA;QAC9D,IAAI,CAACnD,EAAE,CAACoD,YAAY,GAAGzD,EAAG6C,CAAAA,GAAG,CAAC,IAAI,CAACM,aAAa,CAACG,CAAC,EAAE,IAAI,CAACH,aAAa,CAACK,CAAC,EAAE,IAAI,CAACJ,OAAO,CAAA;AACtF,QAAA,IAAI,CAAC/C,EAAE,CAACqD,QAAQ,GAAG1D,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAAC6C,YAAY,EAAE,IAAI,CAACA,YAAY,CAAA;AAC9D,QAAA,IAAI,CAACpD,EAAE,CAACsD,CAAC,GAAG3D,EAAG6C,CAAAA,GAAG,CAChBe,eAAAA,CAAgB,IAAI,EAAE,IAAI,CAACZ,GAAG,EAAE,IAAI,CAACK,EAAE,CACvCO,EAAAA,eAAAA,CAAgB,IAAI,EAAE,IAAI,CAACZ,GAAG,EAAE,IAAI,CAACO,EAAE,CACvC,EAAA,IAAI,CAACH,OAAO,CAAA;AAEd,QAAA,IAAI,CAAC/C,EAAE,CAACwD,GAAG,GAAGC,MAAAA,CAAO,IAAI,EAAE,IAAI,CAAChE,GAAG,EAAE,IAAI,CAACmD,GAAG,EAAE,IAAI,CAACS,QAAQ,CAAA;QAC5D,IAAI,CAACrD,EAAE,CAAC0D,CAAC,GAAGC,cAAe,CAAA,IAAI,EAAE,IAAI,CAACd,GAAG,EAAElD,EAAAA,CAAGiE,IAAI,CAAC,IAAI,CAAClE,EAAE,CAAA,EAAGC,EAAGiE,CAAAA,IAAI,CAAC,CAAA,CAAA,CAAA;AACrE,QAAA,IAAI,CAAC3C,OAAO,CAACtB,GAAGY,GAAG,CAAC,IAAI,CAACmD,CAAC,EAAE/D,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAAC+C,CAAC,EAAE,IAAI,CAACE,GAAG,CAAA,CAAA,CAAA;AACrD,KAAA,CAAA;IAEF,OAAO7D,EAAAA,CAAGwB,cAAc,EAAE,CAACtB,QAAAA,CAAS,CAAC8C,GAAAA,EAAKlD,GAAKmD,EAAAA,GAAAA,EAAKC,GAAKC,EAAAA,aAAAA,EAAeC,OAASrD,EAAAA,EAAAA,CAAAA;AACnF;AAEA;;;;;;;;;;;;;;;;;;;;IAqBA,SAASmE,oBACPtE,CAAAA,KAA4B,EAC5BuE,GAAgB,EAChBC,SAAsB,EACtBC,GAAgB,EAChBC,MAAmB,EAAA;IAEnB,MAAMtE,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,OAAOD,GAAGuE,kBAAkB,CAACJ,KAAKnE,EAAGuB,CAAAA,IAAI,CAACvB,EAAGY,CAAAA,GAAG,CAACZ,EAAAA,CAAGe,GAAG,CAACuD,MAAAA,EAAQ,MAAMF,SAAUd,CAAAA,CAAC,GAAGe,GAAM,CAAA,EAAA,CAAA,CAAA;AAC5F;AAEA;;;;;;;;;;;;;;;;AAgBC,IACM,SAASG,gBAAAA,CACd5E,KAA4B,EAC5B6E,UAAuB,EACvBC,gBAA6B,EAC7BC,SAAsB,EACtBC,SAAsB,EACtBC,MAAmB,EACnBC,OAAoB,EACpBC,QAAqB,EAAA;IAErB,MAAM/E,EAAAA,GAAKJ,MAAMK,QAAQ;IACzB,MAAM+E,SAAAA,GAAYC,WAAWC,qBAAqB;IAClD,MAAMC,OAAAA,GAAUF,WAAWG,mBAAmB;IAC9C,MAAMC,aAAAA,GAAgBL,YAAYG,OAAU,GAAA,CAAA;;;;IAI5CvF,KAAMS,CAAAA,EAAE,CAACiF,UAAU,GAAGtF,GAAGY,GAAG,CAACZ,EAAGe,CAAAA,GAAG,CAACf,EAAAA,CAAGY,GAAG,CAACZ,EAAAA,CAAGW,KAAK,CAACgE,SAAAA,EAAW,GAAG,CAAI,CAAA,EAAA,GAAA,CAAA,EAAM,GAAMD,CAAAA,EAAAA,gBAAAA,CAAiBlB,CAAC,CAAA;;;;AAIpG5D,IAAAA,KAAAA,CAAMS,EAAE,CAACkF,YAAY,GAAGvF,EAAAA,CAAGY,GAAG,CAACZ,EAAAA,CAAGW,KAAK,CAACX,GAAGU,GAAG,CAAC,GAAGkE,SAAY,CAAA,EAAA,CAAA,EAAG,IAAIO,OAAU,GAAA,CAAA,CAAA;IAC/EvF,KAAMS,CAAAA,EAAE,CAACmF,SAAS,GAAGxF,GAAGyF,KAAK,CAAC7F,MAAM2F,YAAY,CAAA;IAChD3F,KAAMS,CAAAA,EAAE,CAACqF,SAAS,GAAG1F,GAAGe,GAAG,CAACnB,KAAM4F,CAAAA,SAAS,EAAER,SAAAA,CAAAA;AAC7CpF,IAAAA,KAAAA,CAAMS,EAAE,CAACsF,SAAS,GAAG3F,EAAG4F,CAAAA,GAAG,CAAC5F,EAAAA,CAAGe,GAAG,CAACnB,KAAM8F,CAAAA,SAAS,EAAE,CAAIL,CAAAA,EAAAA,aAAAA,CAAAA;AACxDzF,IAAAA,KAAAA,CAAMS,EAAE,CAACwF,cAAc,GAAG7F,EAAAA,CAAG6C,GAAG,CAC9BqB,oBAAAA,CAAqBtE,KAAO6E,EAAAA,UAAAA,EAAYC,kBAAkB9E,KAAM0F,CAAAA,UAAU,EAAE1F,KAAM8F,CAAAA,SAAS,EAAEI,GAAG,EAChG5B,oBAAqBtE,CAAAA,KAAAA,EAAO6E,YAAYC,gBAAkB9E,EAAAA,KAAAA,CAAM0F,UAAU,EAAE1F,MAAM+F,SAAS,CAAA,CAAEG,GAAG,EAChG9F,GAAGU,GAAG,CAACd,MAAM2F,YAAY,EAAE3F,MAAM4F,SAAS,CAAA,CAAA;;AAG5C5F,IAAAA,KAAAA,CAAMS,EAAE,CAAC0F,aAAa,GAAG7B,qBACvBtE,KACA6E,EAAAA,UAAAA,EACAC,gBACA9E,EAAAA,KAAAA,CAAM0F,UAAU,EAChBtF,EAAAA,CAAGI,KAAK,CAAC6E,WAAWe,uBAAuB,CAAA,CAAA;IAE7C,OAAOC,qBAAAA,CACLrG,KACAA,EAAAA,KAAAA,CAAMiG,cAAc,EACpBjG,MAAMmG,aAAa,EACnBnB,SACAC,EAAAA,MAAAA,EACAC,OACAC,EAAAA,QAAAA,CAAAA;AAEJ;AAEA;;;;;;;;;;;;;;;;;AAiBC,IACD,SAASkB,qBAAAA,CACPrG,KAA4B,EAC5BsG,OAAoB,EACpBC,MAAmB,EACnBvB,SAAsB,EACtBC,MAAmB,EACnBC,OAAoB,EACpBC,QAAqB,EAAA;IAErB,MAAM/E,EAAAA,GAAKJ,MAAMK,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,2BAAA;IACjBF,EAAGG,CAAAA,IAAI,CACLD,QACA,EAAA;AACEF,QAAAA,EAAAA,CAAGiE,IAAI,CAAC,SAAA,CAAA;AACRjE,QAAAA,EAAAA,CAAGoG,IAAI,CAAC,QAAA,CAAA;AACRpG,QAAAA,EAAAA,CAAGI,KAAK,CAAC,WAAA,CAAA;AACTJ,QAAAA,EAAAA,CAAGiE,IAAI,CAAC,QAAA,CAAA;AACRjE,QAAAA,EAAAA,CAAGiE,IAAI,CAAC,SAAA,CAAA;AACRjE,QAAAA,EAAAA,CAAGiE,IAAI,CAAC,UAAA;KACT,EACD,WAAA;;;;;;AAME,QAAA,IAAI,CAACoC,GAAG,CAACrG,EAAAA,CAAGsG,WAAW,CAAC,IAAI,CAAC1B,SAAS,EAAE,GAAA,IAAO,CAAA,GAAI2B,kCAAiC,CAAK,CAAA,EAAA,WAAA;AACvF,YAAA,IAAI,CAACjF,OAAO,CAACtB,EAAAA,CAAGiE,IAAI,CAAC,CAAA,CAAA,CAAA;AACvB,SAAA,CAAA;QACA,IAAI,CAAC5D,EAAE,CAACmG,OAAO,GAAGxG,EAAGyG,CAAAA,OAAO,CAAC,IAAI,CAAC3B,OAAO,EAAE9E,EAAG0G,CAAAA,GAAG,CAAC,IAAI,CAAC7B,MAAM,GAAG,IAAI,CAACsB,MAAM,CAACQ,CAAC,CAAA;AAC7E,QAAA,IAAI,CAACtG,EAAE,CAACuG,IAAI,GAAG5G,GAAG6G,GAAG,CAAC7G,EAAG0G,CAAAA,GAAG,CAAC,IAAI,CAAC3B,QAAQ,CAAG,EAAA,IAAI,CAACyB,OAAO,CAAA;QACzD,IAAI,CAACnG,EAAE,CAACyG,CAAC,GAAG,IAAI,CAACX,MAAM,CAACY,CAAC;AACzB,QAAA,IAAI,CAAC1G,EAAE,CAAC2G,EAAE,GAAGhH,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACkG,CAAC,EAAE,IAAI,CAACA,CAAC,CAAA;AAClC,QAAA,IAAI,CAACzG,EAAE,CAAC4G,IAAI,GAAGjH,EAAAA,CAAGU,GAAG,CAAC,CAAA,EAAGV,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACkG,CAAC,EAAE,IAAI,CAACF,IAAI,CAAA,CAAA;AACjD,QAAA,IAAI,CAACvG,EAAE,CAAC6G,KAAK,GAAGlH,EAAGY,CAAAA,GAAG,CAAC,IAAI,CAACqG,IAAI,EAAE,IAAI,CAACA,IAAI,CAAA;QAC3C,IAAI,CAAC5G,EAAE,CAAC8G,KAAK,GAAGnH,EAAGc,CAAAA,GAAG,CAACd,EAAAA,CAAGY,GAAG,CAAC,GAAKZ,EAAAA,EAAAA,CAAGU,GAAG,CAAC,CAAA,EAAG,IAAI,CAACsG,EAAE,CAAA,CAAA,EAAIhH,EAAGO,CAAAA,GAAG,CAAC,IAAI,CAAC2G,KAAK,EAAE,IAAA,CAAA,CAAA;AAC3E,QAAA,IAAI,CAAC5F,OAAO,CAACtB,EAAAA,CAAGY,GAAG,CAAC,IAAI,CAACsF,OAAO,EAAE,IAAI,CAACiB,KAAK,CAAA,CAAA;AAC9C,KAAA,CAAA;IAEF,OAAOnH,EAAAA,CAAGwB,cAAc,EAAE,CAACtB,QAAAA,CAAS,CAACgG,OAASC,EAAAA,MAAAA,EAAQvB,SAAWC,EAAAA,MAAAA,EAAQC,OAASC,EAAAA,QAAAA,CAAAA;AACpF;;;;"}

package/dist/shaders/ssgi.js

@@ -0,0 +1,63 @@
+/**
+ * Animated interleaved-gradient-noise pair used to decorrelate SSGI rays.
+ *
+ * Produces two low-discrepancy random numbers in [0, 1) for a given fragment,
+ * advanced per sample/frame through golden-ratio sequences so the temporal
+ * resolve pass converges the Monte-Carlo estimate.
+ *
+ * @param scope - The shader scope
+ * @param fragCoord - Fragment coordinate in pixels
+ * @param seed - Sequence index (frameIndex * sampleCount + sampleIndex)
+ * @returns A vec2 of pseudo random numbers in [0, 1)
+ *
+ * @internal
+ */ function SSGI_rand2(scope, fragCoord, seed) {
+    const pb = scope.$builder;
+    pb.func('SSGI_ign', [
+        pb.vec2('p')
+    ], function() {
+        this.$return(pb.fract(pb.mul(52.9829189, pb.fract(pb.add(pb.mul(0.06711056, this.p.x), pb.mul(0.00583715, this.p.y))))));
+    });
+    pb.func('SSGI_rand2', [
+        pb.vec2('fragCoord'),
+        pb.float('seed')
+    ], function() {
+        this.$l.n1 = this.SSGI_ign(this.fragCoord);
+        this.$l.n2 = this.SSGI_ign(pb.add(this.fragCoord, pb.vec2(5.588238, 5.588238)));
+        // Golden-ratio (and its square) sequences advance the pair every sample
+        this.$return(pb.fract(pb.add(pb.vec2(this.n1, this.n2), pb.mul(pb.vec2(0.61803398875, 0.38196601125), this.seed))));
+    });
+    return scope.SSGI_rand2(fragCoord, seed);
+}
+/**
+ * Cosine-weighted hemisphere sampling around a normal.
+ *
+ * With cosine-weighted sampling the Monte-Carlo estimator of the rendering
+ * equation's diffuse term reduces to the plain average of the sampled
+ * radiances, so the SSGI trace pass can simply accumulate hit colors.
+ *
+ * @param scope - The shader scope
+ * @param normal - Unit normal the hemisphere is oriented around
+ * @param rand - Two random numbers in [0, 1)
+ * @returns A unit direction inside the hemisphere
+ *
+ * @internal
+ */ function SSGI_cosineSampleHemisphere(scope, normal, rand) {
+    const pb = scope.$builder;
+    pb.func('SSGI_cosineSampleHemisphere', [
+        pb.vec3('n'),
+        pb.vec2('u')
+    ], function() {
+        this.$l.up = this.$choice(pb.lessThan(pb.abs(this.n.z), 0.999), pb.vec3(0, 0, 1), pb.vec3(1, 0, 0));
+        this.$l.tangent = pb.normalize(pb.cross(this.up, this.n));
+        this.$l.bitangent = pb.cross(this.n, this.tangent);
+        this.$l.r = pb.sqrt(this.u.x);
+        this.$l.phi = pb.mul(this.u.y, Math.PI * 2);
+        this.$l.dir = pb.add(pb.mul(this.tangent, pb.mul(this.r, pb.cos(this.phi))), pb.mul(this.bitangent, pb.mul(this.r, pb.sin(this.phi))), pb.mul(this.n, pb.sqrt(pb.max(0, pb.sub(1, this.u.x)))));
+        this.$return(pb.normalize(this.dir));
+    });
+    return scope.SSGI_cosineSampleHemisphere(normal, rand);
+}
+
+export { SSGI_cosineSampleHemisphere, SSGI_rand2 };
+//# sourceMappingURL=ssgi.js.map

package/dist/shaders/ssgi.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ssgi.js","sources":["../../src/shaders/ssgi.ts"],"sourcesContent":["import type { PBInsideFunctionScope, PBShaderExp } from '@zephyr3d/device';\n\n/**\n * Animated interleaved-gradient-noise pair used to decorrelate SSGI rays.\n *\n * Produces two low-discrepancy random numbers in [0, 1) for a given fragment,\n * advanced per sample/frame through golden-ratio sequences so the temporal\n * resolve pass converges the Monte-Carlo estimate.\n *\n * @param scope - The shader scope\n * @param fragCoord - Fragment coordinate in pixels\n * @param seed - Sequence index (frameIndex * sampleCount + sampleIndex)\n * @returns A vec2 of pseudo random numbers in [0, 1)\n *\n * @internal\n */\nexport function SSGI_rand2(scope: PBInsideFunctionScope, fragCoord: PBShaderExp, seed: PBShaderExp) {\n  const pb = scope.$builder;\n  pb.func('SSGI_ign', [pb.vec2('p')], function () {\n    this.$return(\n      pb.fract(\n        pb.mul(52.9829189, pb.fract(pb.add(pb.mul(0.06711056, this.p.x), pb.mul(0.00583715, this.p.y))))\n      )\n    );\n  });\n  pb.func('SSGI_rand2', [pb.vec2('fragCoord'), pb.float('seed')], function () {\n    this.$l.n1 = this.SSGI_ign(this.fragCoord);\n    this.$l.n2 = this.SSGI_ign(pb.add(this.fragCoord, pb.vec2(5.588238, 5.588238)));\n    // Golden-ratio (and its square) sequences advance the pair every sample\n    this.$return(\n      pb.fract(pb.add(pb.vec2(this.n1, this.n2), pb.mul(pb.vec2(0.61803398875, 0.38196601125), this.seed)))\n    );\n  });\n  return scope.SSGI_rand2(fragCoord, seed);\n}\n\n/**\n * Cosine-weighted hemisphere sampling around a normal.\n *\n * With cosine-weighted sampling the Monte-Carlo estimator of the rendering\n * equation's diffuse term reduces to the plain average of the sampled\n * radiances, so the SSGI trace pass can simply accumulate hit colors.\n *\n * @param scope - The shader scope\n * @param normal - Unit normal the hemisphere is oriented around\n * @param rand - Two random numbers in [0, 1)\n * @returns A unit direction inside the hemisphere\n *\n * @internal\n */\nexport function SSGI_cosineSampleHemisphere(\n  scope: PBInsideFunctionScope,\n  normal: PBShaderExp,\n  rand: PBShaderExp\n) {\n  const pb = scope.$builder;\n  pb.func('SSGI_cosineSampleHemisphere', [pb.vec3('n'), pb.vec2('u')], function () {\n    this.$l.up = this.$choice(pb.lessThan(pb.abs(this.n.z), 0.999), pb.vec3(0, 0, 1), pb.vec3(1, 0, 0));\n    this.$l.tangent = pb.normalize(pb.cross(this.up, this.n));\n    this.$l.bitangent = pb.cross(this.n, this.tangent);\n    this.$l.r = pb.sqrt(this.u.x);\n    this.$l.phi = pb.mul(this.u.y, Math.PI * 2);\n    this.$l.dir = pb.add(\n      pb.mul(this.tangent, pb.mul(this.r, pb.cos(this.phi))),\n      pb.mul(this.bitangent, pb.mul(this.r, pb.sin(this.phi))),\n      pb.mul(this.n, pb.sqrt(pb.max(0, pb.sub(1, this.u.x))))\n    );\n    this.$return(pb.normalize(this.dir));\n  });\n  return scope.SSGI_cosineSampleHemisphere(normal, rand);\n}\n"],"names":["SSGI_rand2","scope","fragCoord","seed","pb","$builder","func","vec2","$return","fract","mul","add","p","x","y","float","$l","n1","SSGI_ign","n2","SSGI_cosineSampleHemisphere","normal","rand","vec3","up","$choice","lessThan","abs","n","z","tangent","normalize","cross","bitangent","r","sqrt","u","phi","Math","PI","dir","cos","sin","max","sub"],"mappings":"AAEA;;;;;;;;;;;;;AAaC,IACM,SAASA,UAAAA,CAAWC,KAA4B,EAAEC,SAAsB,EAAEC,IAAiB,EAAA;IAChG,MAAMC,EAAAA,GAAKH,MAAMI,QAAQ;IACzBD,EAAGE,CAAAA,IAAI,CAAC,UAAY,EAAA;AAACF,QAAAA,EAAAA,CAAGG,IAAI,CAAC,GAAA;KAAK,EAAE,WAAA;AAClC,QAAA,IAAI,CAACC,OAAO,CACVJ,EAAAA,CAAGK,KAAK,CACNL,EAAAA,CAAGM,GAAG,CAAC,YAAYN,EAAGK,CAAAA,KAAK,CAACL,EAAAA,CAAGO,GAAG,CAACP,EAAAA,CAAGM,GAAG,CAAC,UAAY,EAAA,IAAI,CAACE,CAAC,CAACC,CAAC,CAAA,EAAGT,EAAGM,CAAAA,GAAG,CAAC,UAAY,EAAA,IAAI,CAACE,CAAC,CAACE,CAAC,CAAA,CAAA,CAAA,CAAA,CAAA,CAAA;AAGlG,KAAA,CAAA;IACAV,EAAGE,CAAAA,IAAI,CAAC,YAAc,EAAA;AAACF,QAAAA,EAAAA,CAAGG,IAAI,CAAC,WAAA,CAAA;AAAcH,QAAAA,EAAAA,CAAGW,KAAK,CAAC,MAAA;KAAQ,EAAE,WAAA;QAC9D,IAAI,CAACC,EAAE,CAACC,EAAE,GAAG,IAAI,CAACC,QAAQ,CAAC,IAAI,CAAChB,SAAS,CAAA;AACzC,QAAA,IAAI,CAACc,EAAE,CAACG,EAAE,GAAG,IAAI,CAACD,QAAQ,CAACd,GAAGO,GAAG,CAAC,IAAI,CAACT,SAAS,EAAEE,EAAGG,CAAAA,IAAI,CAAC,QAAU,EAAA,QAAA,CAAA,CAAA,CAAA;;AAEpE,QAAA,IAAI,CAACC,OAAO,CACVJ,EAAAA,CAAGK,KAAK,CAACL,EAAGO,CAAAA,GAAG,CAACP,EAAAA,CAAGG,IAAI,CAAC,IAAI,CAACU,EAAE,EAAE,IAAI,CAACE,EAAE,CAAGf,EAAAA,EAAAA,CAAGM,GAAG,CAACN,EAAGG,CAAAA,IAAI,CAAC,aAAA,EAAe,aAAgB,CAAA,EAAA,IAAI,CAACJ,IAAI,CAAA,CAAA,CAAA,CAAA;AAEtG,KAAA,CAAA;IACA,OAAOF,KAAAA,CAAMD,UAAU,CAACE,SAAWC,EAAAA,IAAAA,CAAAA;AACrC;AAEA;;;;;;;;;;;;;AAaC,IACM,SAASiB,2BAAAA,CACdnB,KAA4B,EAC5BoB,MAAmB,EACnBC,IAAiB,EAAA;IAEjB,MAAMlB,EAAAA,GAAKH,MAAMI,QAAQ;IACzBD,EAAGE,CAAAA,IAAI,CAAC,6BAA+B,EAAA;AAACF,QAAAA,EAAAA,CAAGmB,IAAI,CAAC,GAAA,CAAA;AAAMnB,QAAAA,EAAAA,CAAGG,IAAI,CAAC,GAAA;KAAK,EAAE,WAAA;AACnE,QAAA,IAAI,CAACS,EAAE,CAACQ,EAAE,GAAG,IAAI,CAACC,OAAO,CAACrB,EAAGsB,CAAAA,QAAQ,CAACtB,EAAAA,CAAGuB,GAAG,CAAC,IAAI,CAACC,CAAC,CAACC,CAAC,CAAG,EAAA,KAAA,CAAA,EAAQzB,GAAGmB,IAAI,CAAC,CAAG,EAAA,CAAA,EAAG,CAAInB,CAAAA,EAAAA,EAAAA,CAAGmB,IAAI,CAAC,GAAG,CAAG,EAAA,CAAA,CAAA,CAAA;AAChG,QAAA,IAAI,CAACP,EAAE,CAACc,OAAO,GAAG1B,GAAG2B,SAAS,CAAC3B,EAAG4B,CAAAA,KAAK,CAAC,IAAI,CAACR,EAAE,EAAE,IAAI,CAACI,CAAC,CAAA,CAAA;AACvD,QAAA,IAAI,CAACZ,EAAE,CAACiB,SAAS,GAAG7B,EAAG4B,CAAAA,KAAK,CAAC,IAAI,CAACJ,CAAC,EAAE,IAAI,CAACE,OAAO,CAAA;AACjD,QAAA,IAAI,CAACd,EAAE,CAACkB,CAAC,GAAG9B,EAAAA,CAAG+B,IAAI,CAAC,IAAI,CAACC,CAAC,CAACvB,CAAC,CAAA;AAC5B,QAAA,IAAI,CAACG,EAAE,CAACqB,GAAG,GAAGjC,GAAGM,GAAG,CAAC,IAAI,CAAC0B,CAAC,CAACtB,CAAC,EAAEwB,IAAAA,CAAKC,EAAE,GAAG,CAAA,CAAA;AACzC,QAAA,IAAI,CAACvB,EAAE,CAACwB,GAAG,GAAGpC,EAAGO,CAAAA,GAAG,CAClBP,EAAAA,CAAGM,GAAG,CAAC,IAAI,CAACoB,OAAO,EAAE1B,EAAGM,CAAAA,GAAG,CAAC,IAAI,CAACwB,CAAC,EAAE9B,EAAAA,CAAGqC,GAAG,CAAC,IAAI,CAACJ,GAAG,KACnDjC,EAAGM,CAAAA,GAAG,CAAC,IAAI,CAACuB,SAAS,EAAE7B,EAAGM,CAAAA,GAAG,CAAC,IAAI,CAACwB,CAAC,EAAE9B,EAAGsC,CAAAA,GAAG,CAAC,IAAI,CAACL,GAAG,CAAA,CAAA,CAAA,EACrDjC,EAAGM,CAAAA,GAAG,CAAC,IAAI,CAACkB,CAAC,EAAExB,EAAG+B,CAAAA,IAAI,CAAC/B,EAAAA,CAAGuC,GAAG,CAAC,CAAA,EAAGvC,EAAGwC,CAAAA,GAAG,CAAC,CAAG,EAAA,IAAI,CAACR,CAAC,CAACvB,CAAC,CAAA,CAAA,CAAA,CAAA,CAAA;QAErD,IAAI,CAACL,OAAO,CAACJ,EAAAA,CAAG2B,SAAS,CAAC,IAAI,CAACS,GAAG,CAAA,CAAA;AACpC,KAAA,CAAA;IACA,OAAOvC,KAAAA,CAAMmB,2BAA2B,CAACC,MAAQC,EAAAA,IAAAA,CAAAA;AACnD;;;;"}

package/dist/shaders/ssgi_reprojection.js

@@ -0,0 +1,118 @@
+import { Matrix4x4 } from '@zephyr3d/base';
+
+/**
+ * Shared helpers for validating SSGI temporal history across frames.
+ *
+ * SSGI stores the surface depth of a frame as that frame's view-space Z (the
+ * distance along the camera forward axis, in world units). Comparing such a
+ * value against the *current* frame's view-space Z is only meaningful while the
+ * camera is still: rotating the camera turns the forward axis, so a static world
+ * point changes its view-space Z by roughly
+ *
+ *   dz ~ dTheta * viewZ * tan(halfFovX) * |ndc.x|
+ *
+ * which is zero at the screen center and grows with both depth and screen
+ * eccentricity. At a plain mouse-look rate that already exceeds any reasonable
+ * absolute tolerance for distant off-center pixels, so history gets rejected in
+ * screen-edge bands even though it is perfectly valid. Every consumer therefore
+ * has to convert the current depth into the *previous* frame's view space before
+ * comparing.
+ *
+ * @internal
+ */ /**
+ * Relative share of the compared depth added to the absolute tolerance. 2% keeps
+ * a 10m surface within 0.2m, comfortably above rgba16f depth quantization while
+ * still separating surfaces that a bounce should not share.
+ * @internal
+ */ const SSGI_DEPTH_TOLERANCE_SLOPE = 0.02;
+/**
+ * `prevVP * inverse(currentVP)`: maps a current-frame clip position onto the
+ * previous frame's clip space, where `w` is the previous view-space Z.
+ *
+ * The current half uses the jittered inverse VP because every consumer
+ * reconstructs from a depth buffer that was rasterized with jitter applied.
+ * Jitter only shifts clip x/y, so pairing it with the unjittered `prevVPMatrix`
+ * still produces the correct `w`.
+ *
+ * Falls back to the current VP before any previous transform exists, which
+ * degenerates into comparing the current depth against itself - harmless,
+ * because the history textures are empty on that frame anyway.
+ *
+ * @param ctx - Current draw context
+ * @returns The reprojection matrix
+ * @internal
+ */ function ssgiCurrentToPrevClipMatrix(ctx) {
+    const camera = ctx.camera;
+    const prevVP = camera.prevVPMatrix ?? camera.viewProjectionMatrix;
+    return Matrix4x4.multiply(prevVP, camera.jitteredInvVPMatrix);
+}
+/**
+ * Converts the current frame's linear depth at `uv` into the view-space Z it
+ * would have had in the previous frame, so it can be compared against the depth
+ * stored in the SSGI surface history.
+ *
+ * `currentToPrevClipMatrix` must be `prevViewProjection * inverse(currentViewProjection)`,
+ * mapping a current-frame clip position onto the previous frame's clip space. In
+ * a perspective projection the resulting `w` is exactly the previous view-space
+ * Z in world units, which is the quantity the surface pass wrote.
+ *
+ * @param scope - Current shader scope
+ * @param currentToPrevClipMatrix - `prevVP * invCurrentVP`, as a mat4
+ * @param uv - Screen uv of the pixel being validated
+ * @param linearDepth - Normalized linear depth (viewZ / far) at `uv`
+ * @param nearFar - Camera near and far planes
+ * @returns Previous-frame view-space Z, in world units
+ * @internal
+ */ function ssgiPreviousViewDepth(scope, currentToPrevClipMatrix, uv, linearDepth, nearFar) {
+    const pb = scope.$builder;
+    const funcName = 'Z_ssgiPreviousViewDepth';
+    pb.func(funcName, [
+        pb.mat4('toPrevClip'),
+        pb.vec2('uv'),
+        pb.float('linearDepth'),
+        pb.vec2('nearFar')
+    ], function() {
+        // Invert the perspective depth mapping. `linearDepth` is the normalized
+        // viewZ / far that the depth prepass stored, so the non-linear depth is
+        // (near / linearDepth - far) / (near - far) - identical to the
+        // reconstruction in SSGI_getPosition. Operating on the normalized value
+        // directly matters: converting to world-unit viewZ first and reusing the
+        // same expression silently drops a factor of `far`, which pins every
+        // pixel at the far plane and rejects all history.
+        this.$l.safeLinearDepth = pb.max(this.linearDepth, 1e-6);
+        this.$l.ndcZ = pb.div(pb.sub(pb.div(this.nearFar.x, this.safeLinearDepth), this.nearFar.y), pb.sub(this.nearFar.x, this.nearFar.y));
+        this.$l.clipPos = pb.vec4(pb.sub(pb.mul(this.uv, 2), pb.vec2(1)), pb.sub(pb.mul(pb.clamp(this.ndcZ, 0, 1), 2), 1), 1);
+        // clipPos is already normalized (w = 1 after the NDC divide), so the
+        // reprojected w is directly the previous view-space Z.
+        this.$return(pb.mul(this.toPrevClip, this.clipPos).w);
+    });
+    return scope[funcName](currentToPrevClipMatrix, uv, linearDepth, nearFar);
+}
+/**
+ * Depth tolerance for a history comparison at a given depth.
+ *
+ * A single absolute tolerance cannot serve the whole depth range: it is far too
+ * loose next to the camera, where a few centimetres separate distinct surfaces,
+ * and too tight in the distance, where depth precision and residual reprojection
+ * error both scale with viewZ. Scaling the floor by depth keeps the test at a
+ * roughly constant relative strictness.
+ *
+ * @param scope - Current shader scope
+ * @param baseTolerance - Configured absolute tolerance, in world units
+ * @param viewDepth - View-space Z the comparison happens at, in world units
+ * @returns Effective tolerance, in world units
+ * @internal
+ */ function ssgiDepthTolerance(scope, baseTolerance, viewDepth) {
+    const pb = scope.$builder;
+    const funcName = 'Z_ssgiDepthTolerance';
+    pb.func(funcName, [
+        pb.float('baseTolerance'),
+        pb.float('viewDepth')
+    ], function() {
+        this.$return(pb.max(pb.max(this.baseTolerance, 1e-4), pb.mul(pb.abs(this.viewDepth), SSGI_DEPTH_TOLERANCE_SLOPE)));
+    });
+    return scope[funcName](baseTolerance, viewDepth);
+}
+
+export { SSGI_DEPTH_TOLERANCE_SLOPE, ssgiCurrentToPrevClipMatrix, ssgiDepthTolerance, ssgiPreviousViewDepth };
+//# sourceMappingURL=ssgi_reprojection.js.map

package/dist/shaders/ssgi_reprojection.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ssgi_reprojection.js","sources":["../../src/shaders/ssgi_reprojection.ts"],"sourcesContent":["import { Matrix4x4 } from '@zephyr3d/base';\nimport type { PBInsideFunctionScope, PBShaderExp } from '@zephyr3d/device';\nimport type { DrawContext } from '../render/drawable';\n\n/**\n * Shared helpers for validating SSGI temporal history across frames.\n *\n * SSGI stores the surface depth of a frame as that frame's view-space Z (the\n * distance along the camera forward axis, in world units). Comparing such a\n * value against the *current* frame's view-space Z is only meaningful while the\n * camera is still: rotating the camera turns the forward axis, so a static world\n * point changes its view-space Z by roughly\n *\n *   dz ~ dTheta * viewZ * tan(halfFovX) * |ndc.x|\n *\n * which is zero at the screen center and grows with both depth and screen\n * eccentricity. At a plain mouse-look rate that already exceeds any reasonable\n * absolute tolerance for distant off-center pixels, so history gets rejected in\n * screen-edge bands even though it is perfectly valid. Every consumer therefore\n * has to convert the current depth into the *previous* frame's view space before\n * comparing.\n *\n * @internal\n */\n\n/**\n * Relative share of the compared depth added to the absolute tolerance. 2% keeps\n * a 10m surface within 0.2m, comfortably above rgba16f depth quantization while\n * still separating surfaces that a bounce should not share.\n * @internal\n */\nexport const SSGI_DEPTH_TOLERANCE_SLOPE = 0.02;\n\n/**\n * `prevVP * inverse(currentVP)`: maps a current-frame clip position onto the\n * previous frame's clip space, where `w` is the previous view-space Z.\n *\n * The current half uses the jittered inverse VP because every consumer\n * reconstructs from a depth buffer that was rasterized with jitter applied.\n * Jitter only shifts clip x/y, so pairing it with the unjittered `prevVPMatrix`\n * still produces the correct `w`.\n *\n * Falls back to the current VP before any previous transform exists, which\n * degenerates into comparing the current depth against itself - harmless,\n * because the history textures are empty on that frame anyway.\n *\n * @param ctx - Current draw context\n * @returns The reprojection matrix\n * @internal\n */\nexport function ssgiCurrentToPrevClipMatrix(ctx: DrawContext): Matrix4x4 {\n  const camera = ctx.camera;\n  const prevVP = camera.prevVPMatrix ?? camera.viewProjectionMatrix;\n  return Matrix4x4.multiply(prevVP, camera.jitteredInvVPMatrix);\n}\n\n/**\n * Converts the current frame's linear depth at `uv` into the view-space Z it\n * would have had in the previous frame, so it can be compared against the depth\n * stored in the SSGI surface history.\n *\n * `currentToPrevClipMatrix` must be `prevViewProjection * inverse(currentViewProjection)`,\n * mapping a current-frame clip position onto the previous frame's clip space. In\n * a perspective projection the resulting `w` is exactly the previous view-space\n * Z in world units, which is the quantity the surface pass wrote.\n *\n * @param scope - Current shader scope\n * @param currentToPrevClipMatrix - `prevVP * invCurrentVP`, as a mat4\n * @param uv - Screen uv of the pixel being validated\n * @param linearDepth - Normalized linear depth (viewZ / far) at `uv`\n * @param nearFar - Camera near and far planes\n * @returns Previous-frame view-space Z, in world units\n * @internal\n */\nexport function ssgiPreviousViewDepth(\n  scope: PBInsideFunctionScope,\n  currentToPrevClipMatrix: PBShaderExp,\n  uv: PBShaderExp,\n  linearDepth: PBShaderExp,\n  nearFar: PBShaderExp\n): PBShaderExp {\n  const pb = scope.$builder;\n  const funcName = 'Z_ssgiPreviousViewDepth';\n  pb.func(\n    funcName,\n    [pb.mat4('toPrevClip'), pb.vec2('uv'), pb.float('linearDepth'), pb.vec2('nearFar')],\n    function () {\n      // Invert the perspective depth mapping. `linearDepth` is the normalized\n      // viewZ / far that the depth prepass stored, so the non-linear depth is\n      // (near / linearDepth - far) / (near - far) - identical to the\n      // reconstruction in SSGI_getPosition. Operating on the normalized value\n      // directly matters: converting to world-unit viewZ first and reusing the\n      // same expression silently drops a factor of `far`, which pins every\n      // pixel at the far plane and rejects all history.\n      this.$l.safeLinearDepth = pb.max(this.linearDepth, 1e-6);\n      this.$l.ndcZ = pb.div(\n        pb.sub(pb.div(this.nearFar.x, this.safeLinearDepth), this.nearFar.y),\n        pb.sub(this.nearFar.x, this.nearFar.y)\n      );\n      this.$l.clipPos = pb.vec4(\n        pb.sub(pb.mul(this.uv, 2), pb.vec2(1)),\n        pb.sub(pb.mul(pb.clamp(this.ndcZ, 0, 1), 2), 1),\n        1\n      );\n      // clipPos is already normalized (w = 1 after the NDC divide), so the\n      // reprojected w is directly the previous view-space Z.\n      this.$return(pb.mul(this.toPrevClip, this.clipPos).w);\n    }\n  );\n  return scope[funcName](currentToPrevClipMatrix, uv, linearDepth, nearFar) as PBShaderExp;\n}\n\n/**\n * Depth tolerance for a history comparison at a given depth.\n *\n * A single absolute tolerance cannot serve the whole depth range: it is far too\n * loose next to the camera, where a few centimetres separate distinct surfaces,\n * and too tight in the distance, where depth precision and residual reprojection\n * error both scale with viewZ. Scaling the floor by depth keeps the test at a\n * roughly constant relative strictness.\n *\n * @param scope - Current shader scope\n * @param baseTolerance - Configured absolute tolerance, in world units\n * @param viewDepth - View-space Z the comparison happens at, in world units\n * @returns Effective tolerance, in world units\n * @internal\n */\nexport function ssgiDepthTolerance(\n  scope: PBInsideFunctionScope,\n  baseTolerance: PBShaderExp,\n  viewDepth: PBShaderExp\n): PBShaderExp {\n  const pb = scope.$builder;\n  const funcName = 'Z_ssgiDepthTolerance';\n  pb.func(funcName, [pb.float('baseTolerance'), pb.float('viewDepth')], function () {\n    this.$return(\n      pb.max(pb.max(this.baseTolerance, 1e-4), pb.mul(pb.abs(this.viewDepth), SSGI_DEPTH_TOLERANCE_SLOPE))\n    );\n  });\n  return scope[funcName](baseTolerance, viewDepth) as PBShaderExp;\n}\n"],"names":["SSGI_DEPTH_TOLERANCE_SLOPE","ssgiCurrentToPrevClipMatrix","ctx","camera","prevVP","prevVPMatrix","viewProjectionMatrix","Matrix4x4","multiply","jitteredInvVPMatrix","ssgiPreviousViewDepth","scope","currentToPrevClipMatrix","uv","linearDepth","nearFar","pb","$builder","funcName","func","mat4","vec2","float","$l","safeLinearDepth","max","ndcZ","div","sub","x","y","clipPos","vec4","mul","clamp","$return","toPrevClip","w","ssgiDepthTolerance","baseTolerance","viewDepth","abs"],"mappings":";;AAIA;;;;;;;;;;;;;;;;;;;;;;;;IA2BaA,MAAAA,0BAAAA,GAA6B;AAE1C;;;;;;;;;;;;;;;;IAiBO,SAASC,2BAAAA,CAA4BC,GAAgB,EAAA;IAC1D,MAAMC,MAAAA,GAASD,IAAIC,MAAM;AACzB,IAAA,MAAMC,MAASD,GAAAA,MAAAA,CAAOE,YAAY,IAAIF,OAAOG,oBAAoB;AACjE,IAAA,OAAOC,SAAUC,CAAAA,QAAQ,CAACJ,MAAAA,EAAQD,OAAOM,mBAAmB,CAAA;AAC9D;AAEA;;;;;;;;;;;;;;;;;IAkBO,SAASC,qBAAAA,CACdC,KAA4B,EAC5BC,uBAAoC,EACpCC,EAAe,EACfC,WAAwB,EACxBC,OAAoB,EAAA;IAEpB,MAAMC,EAAAA,GAAKL,MAAMM,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,yBAAA;IACjBF,EAAGG,CAAAA,IAAI,CACLD,QACA,EAAA;AAACF,QAAAA,EAAAA,CAAGI,IAAI,CAAC,YAAA,CAAA;AAAeJ,QAAAA,EAAAA,CAAGK,IAAI,CAAC,IAAA,CAAA;AAAOL,QAAAA,EAAAA,CAAGM,KAAK,CAAC,aAAA,CAAA;AAAgBN,QAAAA,EAAAA,CAAGK,IAAI,CAAC,SAAA;KAAW,EACnF,WAAA;;;;;;;;QAQE,IAAI,CAACE,EAAE,CAACC,eAAe,GAAGR,EAAGS,CAAAA,GAAG,CAAC,IAAI,CAACX,WAAW,EAAE,IAAA,CAAA;AACnD,QAAA,IAAI,CAACS,EAAE,CAACG,IAAI,GAAGV,EAAAA,CAAGW,GAAG,CACnBX,EAAAA,CAAGY,GAAG,CAACZ,EAAAA,CAAGW,GAAG,CAAC,IAAI,CAACZ,OAAO,CAACc,CAAC,EAAE,IAAI,CAACL,eAAe,CAAA,EAAG,IAAI,CAACT,OAAO,CAACe,CAAC,CAAA,EACnEd,GAAGY,GAAG,CAAC,IAAI,CAACb,OAAO,CAACc,CAAC,EAAE,IAAI,CAACd,OAAO,CAACe,CAAC,CAAA,CAAA;AAEvC,QAAA,IAAI,CAACP,EAAE,CAACQ,OAAO,GAAGf,GAAGgB,IAAI,CACvBhB,EAAGY,CAAAA,GAAG,CAACZ,EAAGiB,CAAAA,GAAG,CAAC,IAAI,CAACpB,EAAE,EAAE,CAAA,CAAA,EAAIG,EAAGK,CAAAA,IAAI,CAAC,CACnCL,CAAAA,CAAAA,EAAAA,EAAAA,CAAGY,GAAG,CAACZ,EAAAA,CAAGiB,GAAG,CAACjB,EAAAA,CAAGkB,KAAK,CAAC,IAAI,CAACR,IAAI,EAAE,CAAG,EAAA,CAAA,CAAA,EAAI,IAAI,CAC7C,CAAA,EAAA,CAAA,CAAA;;;AAIF,QAAA,IAAI,CAACS,OAAO,CAACnB,EAAAA,CAAGiB,GAAG,CAAC,IAAI,CAACG,UAAU,EAAE,IAAI,CAACL,OAAO,EAAEM,CAAC,CAAA;AACtD,KAAA,CAAA;AAEF,IAAA,OAAO1B,KAAK,CAACO,QAAAA,CAAS,CAACN,uBAAAA,EAAyBC,IAAIC,WAAaC,EAAAA,OAAAA,CAAAA;AACnE;AAEA;;;;;;;;;;;;;;AAcC,IACM,SAASuB,kBAAAA,CACd3B,KAA4B,EAC5B4B,aAA0B,EAC1BC,SAAsB,EAAA;IAEtB,MAAMxB,EAAAA,GAAKL,MAAMM,QAAQ;AACzB,IAAA,MAAMC,QAAW,GAAA,sBAAA;IACjBF,EAAGG,CAAAA,IAAI,CAACD,QAAU,EAAA;AAACF,QAAAA,EAAAA,CAAGM,KAAK,CAAC,eAAA,CAAA;AAAkBN,QAAAA,EAAAA,CAAGM,KAAK,CAAC,WAAA;KAAa,EAAE,WAAA;QACpE,IAAI,CAACa,OAAO,CACVnB,EAAAA,CAAGS,GAAG,CAACT,EAAAA,CAAGS,GAAG,CAAC,IAAI,CAACc,aAAa,EAAE,IAAA,CAAA,EAAOvB,EAAGiB,CAAAA,GAAG,CAACjB,EAAAA,CAAGyB,GAAG,CAAC,IAAI,CAACD,SAAS,CAAGxC,EAAAA,0BAAAA,CAAAA,CAAAA,CAAAA;AAE5E,KAAA,CAAA;AACA,IAAA,OAAOW,KAAK,CAACO,QAAS,CAAA,CAACqB,aAAeC,EAAAA,SAAAA,CAAAA;AACxC;;;;"}