@meimu/mutracer 0.0.0-stage → 0.1.0

package/lib/renderer/glsl/chunks/openPbrBsdf.glsl.js

@@ -0,0 +1,1767 @@
+/*
+ * OpenPBR energy-compensation, EON and thin-film portions Copyright 2026 Adobe.
+ * Licensed under the Apache License, Version 2.0.
+ * https://www.apache.org/licenses/LICENSE-2.0
+ */
+/*
+ * Renderer-specific OpenPBR layered BSDF.
+ *
+ * The prepare/eval/sample/pdf organization and fixed-lobe composition follow:
+ * - https://github.com/adobe/openpbr-bsdf
+ * - https://github.com/portsmouth/OpenPBR-viewer
+ *
+ * This compact GPU implementation reuses the renderer's GGX VNDF and volume
+ * integrator. Adobe's opaque-dielectric energy tables are used for the base
+ * rough-diffuse/specular coupling; fuzz still uses an analytic approximation.
+ */
+import { glsl } from '../../shader.js';
+import { OPENPBR_METAL_ENERGY_LUT_ROW, OPENPBR_METAL_ENERGY_LUT_AVERAGE_ROW, OPENPBR_DIELECTRIC_ENERGY_LUT_ROW, OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW, OPENPBR_DIELECTRIC_REFLECTION_RATIO_LUT_ROW } from '../../openPbrEnergyLut.js';
+export default glsl `
+const int OPENPBR_LOBE_FUZZ = 0;
+const int OPENPBR_LOBE_COAT = 1;
+const int OPENPBR_LOBE_METAL = 2;
+const int OPENPBR_LOBE_SPECULAR = 3;
+const int OPENPBR_LOBE_TRANSMISSION = 4;
+const int OPENPBR_LOBE_DIFFUSE = 5;
+const int OPENPBR_LOBE_SUBSURFACE = 6;
+const int OPENPBR_LOBE_COUNT = 7;
+
+struct bj {
+vec3 tangent;
+vec3 bitangent;
+vec3 normal;
+};
+
+struct ag {
+J material;
+bj basis;
+bj coatBasis;
+bool directionFrames;
+#ifdef NUM_COAT_NORMAL_MAPS
+vec3 coatViewLocal;
+#endif
+vec3 viewLocal;
+vec3 diffuseColor;
+float specularEtaTOverEtaI;
+float refractionEtaTOverEtaI;
+float transmissionReflectionIor;
+float thinFilmExteriorIor;
+float thinFilmBaseIor;
+float specularAlpha;
+float diffuseSpecularEnergyCompensation;
+vec3 coatBaseScale;
+vec3 metalMsScale;
+float metalMsFactor;
+float metalMsProbability;
+
+vec3 dielectricMsReflectionScale;
+vec3 dielectricMsTransmissionScale;
+float dielectricMsViewLoss;
+float dielectricMsReflectionRatio;
+float dielectricMsProbability;
+
+float weights[OPENPBR_LOBE_COUNT];
+float probabilities[OPENPBR_LOBE_COUNT];
+};
+
+float aw(vec3 color) {
+return dot(color, vec3(0.212671, 0.715160, 0.072169));
+}
+
+float cV(float bu) {
+float value = (bu - 1.0) / (bu + 1.0);
+return value * value;
+}
+
+float aY(
+float bu,
+float specularWeight
+) {
+float sE = min(
+cV(bu) * specularWeight,
+0.9999
+);
+float sG = sqrt(sE);
+float kV = (1.0 + sG) / (1.0 - sG);
+return bu < 1.0 ? 1.0 / kV : kV;
+}
+
+float af(float cc, float relativeIor) {
+if (abs(relativeIor - 1.0) < 1e-6) {
+return 0.0;
+}
+return ct(cc, relativeIor);
+}
+
+float aC(float ior) {
+if (ior < 1.0) {
+return 15.0 - ((1.0 / ior - 1.0) / 1.5) * 15.0;
+}
+return 16.0 + ((ior - 1.0) / 1.5) * 15.0;
+}
+
+float F(int rX, int row) {
+return texelFetch(openPbrEnergyLut, ivec2(rX, row), 0).r;
+}
+
+float fY(
+float dI,
+float bR,
+float hO
+) {
+float bA = clamp(dI, 0.0, 31.0);
+float aG = clamp(bR, 0.0, 31.0);
+float bt = clamp(hO, 0.0, 31.0);
+int dP = int(floor(bA));
+int aM = int(floor(aG));
+int cI = int(floor(bt));
+int hB = min(dP + 1, 31);
+int cn = min(aM + 1, 31);
+int eO = min(cI + 1, 31);
+float cp = bA - float(dP);
+float cy = aG - float(aM);
+float ku = bt - float(cI);
+
+float uk = F(cI, dP * 32 + aM);
+float ul = F(eO, dP * 32 + aM);
+float um = F(cI, dP * 32 + cn);
+float un = F(eO, dP * 32 + cn);
+float uo = F(cI, hB * 32 + aM);
+float uq = F(eO, hB * 32 + aM);
+float ur = F(cI, hB * 32 + cn);
+float us = F(eO, hB * 32 + cn);
+float v00 = mix(uk, uo, cp);
+float v01 = mix(ul, uq, cp);
+float v10 = mix(um, ur, cp);
+float v11 = mix(un, us, cp);
+return mix(
+mix(v00, v10, cy),
+mix(v01, v11, cy),
+ku
+);
+}
+
+float gx(float dI, float bR) {
+float bA = clamp(dI, 0.0, 31.0);
+float aG = clamp(bR, 0.0, 31.0);
+int dP = int(floor(bA));
+int aM = int(floor(aG));
+int hB = min(dP + 1, 31);
+int cn = min(aM + 1, 31);
+float cp = bA - float(dP);
+float cy = aG - float(aM);
+float v00 = F(aM, 1024 + dP);
+float v01 = F(cn, 1024 + dP);
+float v10 = F(aM, 1024 + hB);
+float v11 = F(cn, 1024 + hB);
+return mix(
+mix(v00, v10, cp),
+mix(v01, v11, cp),
+cy
+);
+}
+
+
+
+float ef(
+float dI,
+float bR,
+float hO
+) {
+float bA = clamp(dI, 0.0, 31.0);
+float aG = clamp(bR, 0.0, 31.0);
+float bt = clamp(hO, 0.0, 31.0);
+int dP = int(floor(bA));
+int aM = int(floor(aG));
+int cI = int(floor(bt));
+int hB = min(dP + 1, 31);
+int cn = min(aM + 1, 31);
+int eO = min(cI + 1, 31);
+float cp = bA - float(dP);
+float cy = aG - float(aM);
+float ku = bt - float(cI);
+
+float uk = F(cI, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + dP * 32 + aM);
+float ul = F(eO, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + dP * 32 + aM);
+float um = F(cI, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + dP * 32 + cn);
+float un = F(eO, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + dP * 32 + cn);
+float uo = F(cI, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + hB * 32 + aM);
+float uq = F(eO, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + hB * 32 + aM);
+float ur = F(cI, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + hB * 32 + cn);
+float us = F(eO, ${OPENPBR_DIELECTRIC_ENERGY_LUT_ROW} + hB * 32 + cn);
+float v00 = mix(uk, uo, cp);
+float v01 = mix(ul, uq, cp);
+float v10 = mix(um, ur, cp);
+float v11 = mix(un, us, cp);
+return mix(
+mix(v00, v10, cy),
+mix(v01, v11, cy),
+ku
+);
+}
+
+float ee(float dI, float bR) {
+float bA = clamp(dI, 0.0, 31.0);
+float aG = clamp(bR, 0.0, 31.0);
+int dP = int(floor(bA));
+int aM = int(floor(aG));
+int hB = min(dP + 1, 31);
+int cn = min(aM + 1, 31);
+float cp = bA - float(dP);
+float cy = aG - float(aM);
+float v00 = F(aM, ${OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW} + dP);
+float v01 = F(cn, ${OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW} + dP);
+float v10 = F(aM, ${OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW} + hB);
+float v11 = F(cn, ${OPENPBR_DIELECTRIC_ENERGY_LUT_AVERAGE_ROW} + hB);
+return mix(
+mix(v00, v10, cp),
+mix(v01, v11, cp),
+cy
+);
+}
+
+float eF(float dI, float bR) {
+float bA = clamp(dI, 0.0, 31.0);
+float aG = clamp(bR, 0.0, 31.0);
+int dP = int(floor(bA));
+int aM = int(floor(aG));
+int hB = min(dP + 1, 31);
+int cn = min(aM + 1, 31);
+float cp = bA - float(dP);
+float cy = aG - float(aM);
+float v00 = F(aM, ${OPENPBR_DIELECTRIC_REFLECTION_RATIO_LUT_ROW} + dP);
+float v01 = F(cn, ${OPENPBR_DIELECTRIC_REFLECTION_RATIO_LUT_ROW} + dP);
+float v10 = F(aM, ${OPENPBR_DIELECTRIC_REFLECTION_RATIO_LUT_ROW} + hB);
+float v11 = F(cn, ${OPENPBR_DIELECTRIC_REFLECTION_RATIO_LUT_ROW} + hB);
+return mix(
+mix(v00, v10, cp),
+mix(v01, v11, cp),
+cy
+);
+}
+
+float bO(float ior, float alpha, float au) {
+return ef(aC(ior), sqrt(alpha) * 31.0, au * 31.0);
+}
+float cO(float ior, float alpha) {
+return ee(aC(ior), sqrt(alpha) * 31.0);
+}
+float ch(float ior, float alpha) {
+return eF(aC(ior), sqrt(alpha) * 31.0);
+}
+
+
+float bH(float alpha, float au) {
+float ia = clamp(sqrt(alpha) * 31.0, 0.0, 31.0);
+float bt = clamp(au * 31.0, 0.0, 31.0);
+int r0 = int(floor(ia));
+int c0 = int(floor(bt));
+int r1 = min(r0 + 1, 31);
+int c1 = min(c0 + 1, 31);
+return mix(
+mix(F(c0, ${OPENPBR_METAL_ENERGY_LUT_ROW} + r0), F(c1, ${OPENPBR_METAL_ENERGY_LUT_ROW} + r0), fract(bt)),
+mix(F(c0, ${OPENPBR_METAL_ENERGY_LUT_ROW} + r1), F(c1, ${OPENPBR_METAL_ENERGY_LUT_ROW} + r1), fract(bt)),
+fract(ia));
+}
+
+float cu(float alpha) {
+float index = clamp(sqrt(alpha) * 31.0, 0.0, 31.0);
+int i0 = int(floor(index));
+return mix(F(i0, ${OPENPBR_METAL_ENERGY_LUT_AVERAGE_ROW}),
+F(min(i0 + 1, 31), ${OPENPBR_METAL_ENERGY_LUT_AVERAGE_ROW}), fract(index));
+}
+
+float aX(float value, float ior) {
+if (ior > 2.5 || ior < 0.4) {
+float tI = cV(2.5);
+float progress = (cV(ior) - tI) / (1.0 - tI);
+return (1.0 - progress) * value;
+}
+return value;
+}
+
+float aI(
+float ior,
+float alpha,
+float au
+) {
+float value = fY(
+aC(ior),
+sqrt(alpha) * 31.0,
+au * 31.0
+);
+return aX(value, ior);
+}
+
+float cq(float ior, float alpha) {
+float value = gx(
+aC(ior),
+sqrt(alpha) * 31.0
+);
+return aX(value, ior);
+}
+
+vec3 openpbr_to_world(bj basis, vec3 direction) {
+return direction.x * basis.tangent
++ direction.y * basis.bitangent
++ direction.z * basis.normal;
+}
+
+vec3 openpbr_to_local(bj basis, vec3 direction) {
+return vec3(
+dot(direction, basis.tangent),
+dot(direction, basis.bitangent),
+dot(direction, basis.normal)
+);
+}
+
+bj fg(
+vec3 aU,
+vec3 normal,
+float anisotropy
+) {
+bj basis;
+basis.normal = normal;
+basis.tangent = vec3(0.0);
+#ifdef USE_TANGENT
+basis.tangent = aU;
+#endif
+if (dot(basis.tangent, basis.tangent) < 0.001 || anisotropy <= 0.0) {
+aK(normal, basis.tangent, basis.bitangent);
+}
+else {
+basis.bitangent = normalize(cross(normal, basis.tangent));
+basis.tangent = normalize(cross(basis.bitangent, normal));
+}
+return basis;
+}
+
+float jR(
+float av,
+float ak,
+float L,
+float S
+) {
+return (av * L - ak * S)
+/ (av * L + ak * S);
+}
+
+float jQ(
+float av,
+float ak,
+float L,
+float S
+) {
+return (ak * L - av * S)
+/ (av * S + ak * L);
+}
+
+float gX(
+float av,
+float ak,
+float L,
+float S
+) {
+return 2.0 * av * L
+/ (av * L + ak * S);
+}
+
+float gW(
+float av,
+float ak,
+float L,
+float S
+) {
+return 2.0 * av * L
+/ (av * S + ak * L);
+}
+
+vec2 aV(vec2 a, vec2 b) {
+return vec2(a.x * b.x - a.y * b.y, a.x * b.y + a.y * b.x);
+}
+
+vec2 aZ(vec2 a, vec2 b) {
+float ap = max(dot(b, b), 1e-9);
+return vec2(
+(a.x * b.x + a.y * b.y) / ap,
+(a.y * b.x - a.x * b.y) / ap
+);
+}
+
+vec2 hU(vec2 value) {
+float ns = length(value);
+float sC = sqrt(max(0.0, 0.5 * (ns + value.x)));
+float fc = sqrt(max(0.0, 0.5 * (ns - value.x)));
+return vec2(sC, value.y >= 0.0 ? fc : -fc);
+}
+
+vec2 bX(
+float L,
+float av,
+vec2 ak,
+bool qu
+) {
+float ob = sqrt(max(0.0, 1.0 - L * L));
+vec2 sj = aZ(vec2(av, 0.0), ak);
+vec2 hw = sj * ob;
+vec2 S = hU(
+vec2(1.0, 0.0) - aV(hw, hw)
+);
+if (aV(ak, S).y < 0.0) {
+S = -S;
+}
+
+vec2 fH = vec2(av * L, 0.0);
+vec2 eA = aV(ak, S);
+if (qu) {
+vec2 ci = ak * L;
+vec2 ce = S * av;
+return aZ(
+ci - ce,
+ci + ce
+);
+}
+return aZ(
+fH - eA,
+fH + eA
+);
+}
+
+float eq(
+float iD,
+float gG,
+float ov,
+float mA,
+vec2 jW,
+float phase
+) {
+vec2 mX = vec2(cos(phase), sin(phase));
+vec2 hX = aV(jW, mX);
+vec2 kC = hX * (gG * mA);
+vec2 ap = vec2(1.0, 0.0) - hX * ov;
+vec2 hv = vec2(iD, 0.0)
++ aZ(kC, ap);
+return dot(hv, hv);
+}
+
+float ar(
+float L,
+float wavelength,
+float cm,
+float aL,
+float ck,
+vec2 eu
+) {
+L = clamp(abs(L), 1e-5, 1.0);
+float hq = pow(aL / ck, 2.0)
+* (1.0 - L * L);
+if (hq >= 1.0) {
+return 1.0;
+}
+float cU = sqrt(1.0 - hq);
+float kP = jR(
+aL, ck, L, cU
+);
+float kO = jQ(
+aL, ck, L, cU
+);
+float jd = gX(
+aL, ck, L, cU
+);
+float jc = gW(
+aL, ck, L, cU
+);
+float gj = ck * cU
+/ max(aL * L, 1e-9);
+vec2 sP = bX(
+cU, ck, eu, false
+);
+vec2 sO = bX(
+cU, ck, eu, true
+);
+float phase = 2.0 * PI * (2.0 * ck * cm * cU) / wavelength;
+float tu = eq(
+kP,
+jd,
+-kP,
+jd * gj,
+sP,
+phase
+);
+float tr = eq(
+kO,
+jc,
+-kO,
+jc * gj,
+sO,
+phase
+);
+return clamp(0.5 * (tu + tr), 0.0, 1.0);
+}
+
+vec3 hD(
+float L,
+float cm,
+float ck,
+float eu,
+float aL
+) {
+return vec3(
+ar(
+L, 620.0, cm, aL, ck, vec2(eu, 0.0)
+),
+ar(
+L, 540.0, cm, aL, ck, vec2(eu, 0.0)
+),
+ar(
+L, 450.0, cm, aL, ck, vec2(eu, 0.0)
+)
+);
+}
+
+
+
+vec3 dT(float au, vec3 f0, vec3 tint) {
+vec3 qW = mix(f0, vec3(1.0), aE(au));
+if (all(equal(tint, vec3(1.0)))) { return qW; }
+float dj = 1.0 / 7.0;
+float lW = 1.0 - dj;
+vec3 fF = mix(f0, vec3(1.0), aE(dj));
+float dg = 1.0 - au;
+float qh = au * aE(au) * dg
+/ (dj * pow(lW, 6.0));
+return clamp(qW - qh * (fF - tint * fF), vec3(0.0), vec3(1.0));
+}
+
+vec2 cQ(float f0, float tb) {
+float dM = clamp(f0, 0.0, 0.99);
+float gF = sqrt(dM);
+float qt = (1.0 - dM) / (1.0 + dM);
+float qr = (1.0 + gF) / max(1.0 - gF, 1e-9);
+float ior = mix(qr, qt, clamp(tb, 0.0, 1.0));
+float jG = (
+(ior + 1.0) * (ior + 1.0) * dM
+- (ior - 1.0) * (ior - 1.0)
+) / max(1.0 - dM, 1e-9);
+return vec2(ior, sqrt(max(jG, 0.0)));
+}
+
+vec3 iS(
+float L,
+float cm,
+float ck,
+vec3 qS,
+vec3 jN,
+float aL
+) {
+return vec3(
+ar(
+L,
+620.0,
+cm,
+aL,
+ck,
+cQ(qS.r, jN.r)
+),
+ar(
+L,
+540.0,
+cm,
+aL,
+ck,
+cQ(qS.g, jN.g)
+),
+ar(
+L,
+450.0,
+cm,
+aL,
+ck,
+cQ(qS.b, jN.b)
+)
+);
+}
+
+float bN(float cm, float weight) {
+return clamp(weight, 0.0, 1.0) * smoothstep(1e-9, 30.0, cm);
+}
+
+vec3 bW(
+vec3 ei,
+float L,
+float cm,
+float ck,
+float eu,
+float weight,
+float aL
+) {
+vec3 pd = hD(
+L,
+cm,
+ck,
+eu,
+aL
+);
+return mix(
+ei,
+pd,
+bN(cm, weight)
+);
+}
+
+vec3 lc(
+vec3 ei,
+float L,
+float cm,
+float ck,
+vec3 qS,
+vec3 jN,
+float weight,
+float aL
+) {
+vec3 pd = iS(
+L,
+cm,
+ck,
+qS,
+jN,
+aL
+);
+return mix(
+ei,
+pd,
+bN(cm, weight)
+);
+}
+
+vec3 cD(
+J material,
+vec3 color,
+vec3 viewLocal,
+vec3 lightLocal,
+float fo,
+out float pdf
+) {
+pdf = 0.0;
+if (lightLocal.z <= 0.0) {
+return vec3(0.0);
+}
+
+float roughness = material.diffuseRoughness;
+float gQ = 0.5 - 2.0 / (3.0 * PI);
+float nU = 2.0 / 3.0 - 28.0 / (15.0 * PI);
+float om = lightLocal.z;
+float mn = viewLocal.z;
+float hM = dot(lightLocal, viewLocal) - om * mn;
+float gN = hM > 0.0
+? hM / max(max(om, mn), 1e-7)
+: hM;
+float ih = 1.0 / (1.0 + gQ * roughness);
+vec3 nf = color * INVPI * ih
+* (1.0 + roughness * gN);
+
+float dc = 1.0 - om;
+float cP = 1.0 - mn;
+float oc = dc * (
+0.0571085289 + dc * (
+0.491881867 + dc * (
+-0.332181442 + dc * 0.0714429953
+)
+)
+);
+float mW = cP * (
+0.0571085289 + cP * (
+0.491881867 + cP * (
+-0.332181442 + cP * 0.0714429953
+)
+)
+);
+float pi = (1.0 + roughness * oc)
+/ (1.0 + gQ * roughness);
+float oq = (1.0 + roughness * mW)
+/ (1.0 + gQ * roughness);
+float fQ = ih * (1.0 + nU * roughness);
+vec3 hC = color * color * fQ
+/ max(vec3(1e-7), vec3(1.0) - color * (1.0 - fQ));
+vec3 lp = hC * INVPI
+* max(1e-7, 1.0 - oq)
+* max(1e-7, 1.0 - pi)
+/ max(1e-7, 1.0 - fQ);
+pdf = lightLocal.z * INVPI;
+return (nf + lp) * fo;
+}
+
+float aO(float viewCosine, float roughness) {
+float x = clamp(viewCosine, 0.0, 1.0);
+float y = clamp(roughness, 0.01, 1.0);
+float scale = y * (0.0206607 + 1.58491 * y)
+/ (0.0379424 + y * (1.32227 + y));
+float mean = y * (-0.193854 + y * (-1.14885 + y * (1.7932 - 0.95943 * y * y)))
+/ (0.046391 + y);
+float offset = y * (0.000654023 + (-0.0207818 + 0.119681 * y) * y)
+/ (1.26264 + y * (-1.92021 + y));
+float gaussian = exp(-0.5 * pow((x - mean) / scale, 2.0))
+/ (scale * sqrt(2.0 * PI));
+return max(gaussian + offset, 0.0);
+}
+
+vec3 jP(
+J material,
+vec3 viewLocal,
+vec3 lightLocal,
+out float pdf
+) {
+pdf = 0.0;
+if (viewLocal.z <= 0.0 || lightLocal.z <= 0.0) {
+return vec3(0.0);
+}
+float albedo = aO(
+viewLocal.z,
+material.fuzzRoughness
+);
+pdf = lightLocal.z * INVPI;
+return material.fuzzColor * albedo * INVPI;
+}
+
+vec3 am(
+J material,
+vec3 viewLocal,
+vec3 lightLocal,
+vec3 V,
+vec3 cK,
+out float pdf
+) {
+pdf = 0.0;
+if (lightLocal.z <= 0.0) {
+return vec3(0.0);
+}
+
+vec2 alpha = vec2(material.alphaX, material.alphaY);
+float distribution = gE(V, alpha);
+float jM = bk(viewLocal, alpha);
+float ea = jM
+* bk(lightLocal, alpha);
+
+pdf = jM * distribution / (4.0 * max(viewLocal.z, 1e-8));
+return cK * distribution * ea
+/ (4.0 * max(lightLocal.z * viewLocal.z, 1e-8));
+}
+
+vec3 dk(
+J material,
+float relativeIor,
+vec3 tint,
+vec3 viewLocal,
+vec3 lightLocal,
+vec3 V,
+vec3 cK,
+out float pdf
+) {
+pdf = 0.0;
+if (lightLocal.z >= 0.0) {
+return vec3(0.0);
+}
+
+float kz = dot(lightLocal, V);
+float bQ = dot(viewLocal, V);
+
+
+if (bQ <= 0.0 || kz >= 0.0) {
+return vec3(0.0);
+}
+vec2 alpha = vec2(material.alphaX, material.alphaY);
+float distribution = gE(V, alpha);
+float jM = bk(viewLocal, alpha);
+float ea = jM
+* bk(lightLocal, alpha);
+float ap = kz + bQ * relativeIor;
+ap *= ap;
+float pn = abs(kz) / max(ap, 1e-8);
+float jm = relativeIor * relativeIor;
+
+pdf = jM * max(0.0, bQ) * distribution * pn
+/ max(viewLocal.z, 1e-8);
+return tint * (vec3(1.0) - cK) * distribution * ea
+* abs(bQ) * pn * jm
+/ max(abs(lightLocal.z * viewLocal.z), 1e-8);
+}
+
+vec3 bn(J material) {
+#ifdef HAS_MEDIUM
+if (material.mediumType > 0) {
+return vec3(1.0);
+}
+#endif
+return material.transmissionColor;
+}
+
+float cS(float ah) {
+return 2.0 * ah / (1.0 + ah);
+}
+
+vec3 aP(
+vec3 ah,
+float transmissionWeight
+) {
+vec3 gf = 2.0 * ah
+/ (vec3(1.0) + ah);
+return mix(
+ah,
+gf,
+clamp(transmissionWeight, 0.0, 1.0)
+);
+}
+
+vec3 bG(
+J material,
+float L,
+float bu
+) {
+vec3 tint = bn(material);
+if (!material.thinWalled) {
+return tint;
+}
+float iF = max(0.0, 1.0 - L * L);
+float dN = iF
+/ max(bu * bu, 1e-9);
+if (dN >= 1.0) {
+return vec3(0.0);
+}
+float S = sqrt(1.0 - dN);
+return pow(max(tint, vec3(0.0)), vec3(1.0 / max(S, 1e-5)));
+}
+
+vec3 dL(
+J material,
+vec3 viewLocal,
+vec3 lightLocal,
+out float pdf
+) {
+pdf = 0.0;
+if (viewLocal.z <= 0.0 || lightLocal.z >= 0.0) {
+return vec3(0.0);
+}
+float S = -lightLocal.z;
+float nR = ct(
+viewLocal.z,
+material.relativeIor
+);
+pdf = S * INVPI;
+return vec3(1.0 - nR) * INVPI;
+}
+
+void dQ(inout ag prepared) {
+float total = prepared.probabilities[OPENPBR_LOBE_FUZZ];
+for (int hR = 1; hR < OPENPBR_LOBE_COUNT; hR++) {
+total += prepared.probabilities[hR];
+}
+if (total <= 1e-8) {
+for (int hR = 0; hR < OPENPBR_LOBE_COUNT; hR++) {
+prepared.probabilities[hR] = 0.0;
+}
+return;
+}
+float og = 1.0 / total;
+for (int hR = 0; hR < OPENPBR_LOBE_COUNT; hR++) {
+prepared.probabilities[hR] *= og;
+}
+}
+
+
+
+float bM(float eta) {
+if (abs(eta - 1.0) < 1e-6) { return 0.0; }
+if (eta > 1.0) { return (eta - 1.0) / (4.08567 + 1.00071 * eta); }
+return 0.997118 + 0.1014 * eta - 0.965241 * eta * eta - 0.130607 * eta * eta * eta;
+}
+
+vec3 bP(vec3 f0, vec3 tint) {
+float dj = 1.0 / 7.0;
+vec3 fF = mix(f0, vec3(1.0), aE(dj));
+vec3 b = fF * (vec3(1.0) - tint) / (dj * pow(1.0 - dj, 6.0));
+return clamp(f0 + (vec3(1.0) - f0) / 21.0 - b / 126.0, vec3(0.0), vec3(1.0));
+}
+
+
+
+vec3 cC(J material, float au) {
+if (au <= 0.0) { return vec3(1.0); }
+float lr = sqrt(max(0.0, 1.0 - (1.0 - au * au) / (material.coatIor * material.coatIor)));
+vec3 tl = pow(max(material.coatColor, vec3(0.0)), vec3(0.5 / max(lr, 1e-5)));
+return mix(vec3(1.0), tl, material.coatWeight);
+}
+
+float gi(J material, float au) {
+if (material.coatRoughness <= ROUGHNESS_MIN) {
+return af(au, 1.0 / material.coatIor);
+}
+
+
+return clamp(1.0 - aI(material.coatIor,
+material.coatRoughness * material.coatRoughness, au), 0.0, 1.0);
+}
+
+vec3 gV(ag prepared) {
+J material = prepared.material;
+float dielectric = 1.0 - material.baseMetalness;
+float opaque = dielectric * (1.0 - material.transmissionWeight) * (1.0 - material.subsurfaceWeight);
+float sss = dielectric * (1.0 - material.transmissionWeight) * material.subsurfaceWeight;
+float tT = dielectric * material.transmissionWeight;
+float fL = bM(prepared.specularEtaTOverEtaI);
+float kI = bM(1.0 / prepared.transmissionReflectionIor);
+float pw = clamp(opaque * fL + (sss + tT) * kI + material.baseMetalness, 0.0, 1.0);
+float ji = mix(1.0, material.specularRoughness, pw);
+float ks = bM(material.coatIor);
+float kr = 1.0 - (1.0 - ks) / (material.coatIor * material.coatIor);
+float k = mix(ks, kr, ji);
+vec3 eb = material.baseMetalness * material.specularWeight * bP(material.baseWeight * material.baseColor, material.specularColor)
++ opaque * mix(material.baseWeight * material.baseColor, vec3(1.0), fL)
++ sss * material.subsurfaceColor + vec3(tT);
+vec3 qH = vec3(1.0 - k) / max(vec3(1.0) - eb * k, vec3(1e-6));
+return mix(vec3(1.0), qH, material.coatWeight * material.coatDarkening);
+}
+
+vec3 ds(ag prepared, float kY) {
+if (prepared.material.coatWeight <= 0.0) { return vec3(1.0); }
+if (prepared.material.inside) {
+
+return cC(prepared.material, -kY);
+}
+
+return prepared.coatBaseScale * cC(prepared.material, kY);
+}
+
+
+bj du(bj fallback, vec3 tangent, vec3 bitangent, vec2 direction) {
+vec3 t = tangent - fallback.normal * dot(fallback.normal, tangent);
+if (dot(t, t) < 1e-8) { t = fallback.tangent; }
+t = normalize(t);
+vec3 b = cross(fallback.normal, t) * (dot(cross(fallback.normal, t), bitangent) < 0.0 ? -1.0 : 1.0);
+t = t * direction.x + b * direction.y;
+if (dot(t, t) < 1e-8) { return fallback; }
+bj basis = fallback;
+basis.tangent = normalize(t);
+basis.bitangent = normalize(cross(basis.normal, basis.tangent));
+basis.tangent = normalize(cross(basis.bitangent, basis.normal));
+return basis;
+}
+
+ag eG(
+vec3 tangent,
+J material,
+vec3 bF,
+vec3 normal
+) {
+ag prepared;
+prepared.material = material;
+prepared.basis = fg(
+tangent,
+normal,
+max(material.specularAnisotropy, material.coatAnisotropy)
+);
+prepared.coatBasis = prepared.basis;
+prepared.directionFrames = false;
+#ifdef NUM_COAT_NORMAL_MAPS
+
+
+if (material.coatNormalEnabled) {
+prepared.coatBasis = fg(material.geometryTangent, material.coatNormal, material.coatAnisotropy);
+prepared.directionFrames = true;
+}
+#endif
+if (material.specularAnisotropy > 0.0 && material.anisotropyDirectionEnabled) {
+prepared.basis = du(prepared.basis, material.directionTangent, material.directionBitangent, material.anisotropyDirection);
+prepared.directionFrames = true;
+}
+if (material.coatAnisotropy > 0.0 && material.coatAnisotropyRotation != 0.0) {
+prepared.coatBasis = du(prepared.coatBasis, material.geometryTangent, material.geometryBitangent, vec2(cos(material.coatAnisotropyRotation), sin(material.coatAnisotropyRotation)));
+prepared.directionFrames = true;
+}
+prepared.viewLocal = openpbr_to_local(prepared.basis, bF);
+#ifdef NUM_COAT_NORMAL_MAPS
+prepared.coatViewLocal = openpbr_to_local(prepared.coatBasis, bF);
+#endif
+prepared.diffuseColor = material.baseColor;
+
+float cA = 1.0 - material.baseMetalness;
+float dK = cA
+* (1.0 - material.transmissionWeight);
+prepared.weights[OPENPBR_LOBE_FUZZ] = material.fuzzWeight;
+prepared.weights[OPENPBR_LOBE_COAT] = material.coatWeight;
+if (material.inside) { prepared.weights[OPENPBR_LOBE_COAT] = 0.0; }
+prepared.weights[OPENPBR_LOBE_METAL] = material.baseMetalness;
+prepared.weights[OPENPBR_LOBE_SPECULAR] = cA;
+prepared.weights[OPENPBR_LOBE_TRANSMISSION] = cA
+* material.transmissionWeight;
+prepared.weights[OPENPBR_LOBE_DIFFUSE] = dK
+* (1.0 - material.subsurfaceWeight) * material.baseWeight;
+prepared.weights[OPENPBR_LOBE_SUBSURFACE] = dK
+* material.subsurfaceWeight;
+
+prepared.specularEtaTOverEtaI = aY(
+1.0 / material.relativeIor,
+material.specularWeight
+);
+prepared.material.relativeIor = 1.0 / prepared.specularEtaTOverEtaI;
+prepared.refractionEtaTOverEtaI = prepared.specularEtaTOverEtaI;
+prepared.transmissionReflectionIor = prepared.material.relativeIor;
+
+
+if (material.coatWeight > 0.0) {
+float cj = 1.0;
+if (!material.inside) { cj = mix(1.0, material.coatIor, material.coatWeight); }
+float bq = cj * aY(
+material.specularIor / cj, material.specularWeight);
+float er = bq / cj;
+if (cj > bq && bq >= 1.0) {
+er = 1.0 / er;
+}
+prepared.specularEtaTOverEtaI = er;
+prepared.transmissionReflectionIor = 1.0 / er;
+prepared.refractionEtaTOverEtaI = bq;
+if (material.inside) {
+prepared.refractionEtaTOverEtaI = 1.0 / bq;
+prepared.transmissionReflectionIor = bq;
+}
+prepared.material.relativeIor = 1.0 / prepared.refractionEtaTOverEtaI;
+}
+
+
+
+prepared.thinFilmExteriorIor = 1.0;
+prepared.thinFilmBaseIor = prepared.specularEtaTOverEtaI;
+if (material.coatWeight > 0.0) {
+if (material.inside) {
+prepared.thinFilmExteriorIor = 1.0 / prepared.refractionEtaTOverEtaI;
+prepared.thinFilmBaseIor = 1.0;
+}
+else {
+prepared.thinFilmExteriorIor = mix(1.0, material.coatIor, material.coatWeight);
+prepared.thinFilmBaseIor = prepared.refractionEtaTOverEtaI;
+}
+}
+prepared.specularAlpha = material.specularRoughness
+* material.specularRoughness;
+float viewCosine = abs(prepared.viewLocal.z);
+float ah = af(
+viewCosine,
+1.0 / prepared.specularEtaTOverEtaI
+);
+vec3 ft = material.specularColor * (
+material.thinWalled
+? mix(
+ah,
+cS(ah),
+material.transmissionWeight
+)
+: ah
+);
+vec3 hS = material.specularWeight * dT(
+viewCosine,
+material.baseWeight * material.baseColor,
+material.specularColor
+);
+float eL = material.thinWalled
+? cS(af(viewCosine, prepared.transmissionReflectionIor))
+: af(viewCosine, prepared.transmissionReflectionIor);
+vec3 gH = bG(
+material,
+viewCosine,
+prepared.refractionEtaTOverEtaI
+);
+float im = aI(
+prepared.specularEtaTOverEtaI,
+prepared.specularAlpha,
+viewCosine
+);
+float gk = cq(
+prepared.specularEtaTOverEtaI,
+prepared.specularAlpha
+);
+prepared.diffuseSpecularEnergyCompensation = im
+/ max(gk, 1e-12);
+
+prepared.probabilities[OPENPBR_LOBE_FUZZ] = prepared.weights[OPENPBR_LOBE_FUZZ]
+* aw(material.fuzzColor)
+* max(aO(
+viewCosine,
+material.fuzzRoughness
+), 0.02);
+prepared.probabilities[OPENPBR_LOBE_COAT] = 0.0;
+prepared.coatBaseScale = vec3(1.0);
+if (material.coatWeight > 0.0 && !material.inside) {
+#ifdef NUM_COAT_NORMAL_MAPS
+float coatViewCosine = material.coatNormalEnabled ? max(prepared.coatViewLocal.z, 0.0) : viewCosine;
+float kp = gi(material, coatViewCosine);
+#else
+float kp = gi(material, viewCosine);
+#endif
+prepared.probabilities[OPENPBR_LOBE_COAT] = prepared.weights[OPENPBR_LOBE_COAT] * kp;
+prepared.coatBaseScale = cC(material,
+#ifdef NUM_COAT_NORMAL_MAPS
+coatViewCosine
+#else
+viewCosine
+#endif
+)
+* (1.0 - prepared.weights[OPENPBR_LOBE_COAT] * kp) * gV(prepared);
+}
+prepared.probabilities[OPENPBR_LOBE_METAL] = prepared.weights[OPENPBR_LOBE_METAL]
+* aw(hS);
+prepared.metalMsScale = vec3(0.0);
+prepared.metalMsFactor = 0.0;
+prepared.metalMsProbability = 0.0;
+if (material.specularRoughness >= 0.04 && prepared.weights[OPENPBR_LOBE_METAL] > 0.0) {
+vec3 hI = bP(material.baseWeight * material.baseColor, material.specularColor);
+prepared.metalMsScale = material.specularWeight * hI * hI;
+float pJ = bH(prepared.specularAlpha, viewCosine);
+prepared.metalMsFactor = pJ / max(cu(prepared.specularAlpha), 1e-12);
+float ff = prepared.weights[OPENPBR_LOBE_METAL] * aw(prepared.metalMsScale) * pJ;
+prepared.probabilities[OPENPBR_LOBE_METAL] += ff;
+prepared.metalMsProbability = ff / max(prepared.probabilities[OPENPBR_LOBE_METAL], 1e-12);
+}
+prepared.probabilities[OPENPBR_LOBE_SPECULAR] = prepared.weights[OPENPBR_LOBE_SPECULAR]
+* aw(ft);
+prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] = prepared.weights[OPENPBR_LOBE_TRANSMISSION]
+* (1.0 - eL)
+* aw(gH);
+
+
+
+
+
+
+prepared.dielectricMsReflectionScale = vec3(0.0);
+prepared.dielectricMsTransmissionScale = vec3(0.0);
+prepared.dielectricMsViewLoss = 0.0;
+prepared.dielectricMsReflectionRatio = 0.0;
+prepared.dielectricMsProbability = 0.0;
+if (!material.thinWalled && material.transmissionWeight > 0.0 &&
+prepared.specularAlpha >= 0.0016 && prepared.weights[OPENPBR_LOBE_TRANSMISSION] > 0.0 &&
+abs(prepared.refractionEtaTOverEtaI - 1.0) > 1e-6) {
+float lq = prepared.weights[OPENPBR_LOBE_TRANSMISSION] + prepared.weights[OPENPBR_LOBE_SUBSURFACE];
+prepared.dielectricMsReflectionScale = (material.inside ? vec3(1.0) : material.specularColor) * lq;
+prepared.dielectricMsTransmissionScale = gH * prepared.weights[OPENPBR_LOBE_TRANSMISSION]
++ vec3(prepared.weights[OPENPBR_LOBE_SUBSURFACE]);
+prepared.dielectricMsReflectionRatio = ch(prepared.refractionEtaTOverEtaI, prepared.specularAlpha);
+prepared.dielectricMsViewLoss = bO(prepared.refractionEtaTOverEtaI, prepared.specularAlpha, viewCosine);
+vec3 nI = mix(prepared.dielectricMsTransmissionScale, prepared.dielectricMsReflectionScale, prepared.dielectricMsReflectionRatio);
+float ff = aw(nI) * prepared.dielectricMsViewLoss;
+prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] += ff;
+prepared.dielectricMsProbability = ff / max(prepared.probabilities[OPENPBR_LOBE_TRANSMISSION], 1e-12);
+}
+prepared.probabilities[OPENPBR_LOBE_DIFFUSE] = prepared.weights[OPENPBR_LOBE_DIFFUSE]
+* aw(prepared.diffuseColor)
+* prepared.diffuseSpecularEnergyCompensation;
+prepared.probabilities[OPENPBR_LOBE_SUBSURFACE] = prepared.weights[OPENPBR_LOBE_SUBSURFACE]
+* (1.0 - eL)
+* aw(material.subsurfaceColor)
+* prepared.diffuseSpecularEnergyCompensation;
+
+
+if (material.coatWeight > 0.0 && material.thinFilmWeight > 0.0) {
+float kW = 0.02 * bN(material.thinFilmThickness, material.thinFilmWeight);
+prepared.probabilities[OPENPBR_LOBE_SPECULAR] = max(prepared.probabilities[OPENPBR_LOBE_SPECULAR], prepared.weights[OPENPBR_LOBE_SPECULAR] * kW);
+prepared.probabilities[OPENPBR_LOBE_METAL] = max(prepared.probabilities[OPENPBR_LOBE_METAL], prepared.weights[OPENPBR_LOBE_METAL] * material.specularWeight * kW);
+}
+if (material.coatWeight > 0.0 && !material.inside) {
+float gm = max(max(prepared.coatBaseScale.x, prepared.coatBaseScale.y), prepared.coatBaseScale.z);
+prepared.probabilities[OPENPBR_LOBE_METAL] *= gm;
+prepared.probabilities[OPENPBR_LOBE_SPECULAR] *= gm;
+prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] *= gm;
+prepared.probabilities[OPENPBR_LOBE_DIFFUSE] *= gm;
+prepared.probabilities[OPENPBR_LOBE_SUBSURFACE] *= gm;
+}
+dQ(prepared);
+return prepared;
+}
+
+float cR(ag prepared, float au, bool reflection) {
+if (prepared.material.coatWeight <= 0.0) {
+return af(au, prepared.material.relativeIor);
+}
+float kg = af(au, prepared.transmissionReflectionIor);
+if (!reflection) { return kg; }
+float mU = af(au, 1.0 / prepared.specularEtaTOverEtaI);
+float rC = prepared.material.transmissionWeight + (1.0 - prepared.material.transmissionWeight) * prepared.material.subsurfaceWeight;
+return mix(mU, kg, rC);
+}
+
+vec3 bI(
+ag prepared,
+vec3 lightDirection,
+out float pdf
+) {
+J material = prepared.material;
+vec3 viewLocal = prepared.viewLocal;
+vec3 lightLocal = openpbr_to_local(prepared.basis, lightDirection);
+bool reflection = lightLocal.z * viewLocal.z > 0.0;
+#ifdef NUM_COAT_NORMAL_MAPS
+bool independentCoat = material.coatNormalEnabled;
+bool geometricReflection = dot(material.geometricNormal, lightDirection) > 0.0;
+float coatLightCosine = dot(prepared.coatBasis.normal, lightDirection);
+bool coatReflection = geometricReflection && prepared.coatViewLocal.z > 0.0 && coatLightCosine > 0.0;
+#endif
+
+vec3 V;
+if (reflection) {
+V = normalize(lightLocal + viewLocal);
+}
+else {
+V = normalize(lightLocal + viewLocal * material.relativeIor);
+}
+if (V.z < 0.0) {
+V = -V;
+}
+
+vec3 response = vec3(0.0);
+pdf = 0.0;
+float aH = 0.0;
+float bQ = abs(dot(viewLocal, V));
+float diffuseSpecularEnergyCompensation =
+prepared.diffuseSpecularEnergyCompensation
+* aI(
+prepared.specularEtaTOverEtaI,
+prepared.specularAlpha,
+abs(lightLocal.z)
+);
+
+#ifdef NUM_COAT_NORMAL_MAPS
+if (!independentCoat || (lightLocal.z != 0.0 && reflection == geometricReflection)) {
+#endif
+if (reflection) {
+if (prepared.weights[OPENPBR_LOBE_DIFFUSE] > 0.0) {
+response += prepared.weights[OPENPBR_LOBE_DIFFUSE] * cD(
+material,
+prepared.diffuseColor,
+viewLocal,
+lightLocal,
+diffuseSpecularEnergyCompensation,
+aH
+);
+pdf += prepared.probabilities[OPENPBR_LOBE_DIFFUSE] * aH;
+}
+
+if (prepared.weights[OPENPBR_LOBE_SUBSURFACE] > 0.0 && material.thinWalled) {
+response += 0.5 * prepared.weights[OPENPBR_LOBE_SUBSURFACE] * cD(
+material,
+material.subsurfaceColor,
+viewLocal,
+lightLocal,
+diffuseSpecularEnergyCompensation,
+aH
+);
+pdf += prepared.probabilities[OPENPBR_LOBE_SUBSURFACE] * 0.5 * aH;
+}
+
+if (prepared.weights[OPENPBR_LOBE_SPECULAR] > 0.0) {
+vec3 cr = vec3(cR(prepared, bQ, true));
+if (material.thinFilmWeight > 0.0) {
+cr = bW(
+cr,
+bQ,
+material.thinFilmThickness,
+material.thinFilmIor,
+prepared.thinFilmBaseIor,
+material.thinFilmWeight,
+prepared.thinFilmExteriorIor
+);
+}
+if (material.thinWalled) {
+cr = aP(
+cr,
+material.transmissionWeight
+);
+}
+vec3 cK = material.specularColor * cr;
+response += prepared.weights[OPENPBR_LOBE_SPECULAR] * am(
+material,
+viewLocal,
+lightLocal,
+V,
+cK,
+aH
+);
+pdf += prepared.probabilities[OPENPBR_LOBE_SPECULAR] * aH;
+}
+
+if (prepared.weights[OPENPBR_LOBE_METAL] > 0.0) {
+vec3 cK = dT(
+bQ,
+material.baseWeight * material.baseColor,
+material.specularColor
+);
+if (material.thinFilmWeight > 0.0) {
+cK = lc(
+cK,
+bQ,
+material.thinFilmThickness,
+material.thinFilmIor,
+material.baseWeight * material.baseColor,
+material.specularColor,
+material.thinFilmWeight,
+prepared.thinFilmExteriorIor
+);
+}
+cK *= material.specularWeight;
+response += prepared.weights[OPENPBR_LOBE_METAL] * am(
+material,
+viewLocal,
+lightLocal,
+V,
+cK,
+aH
+);
+if (prepared.metalMsFactor > 0.0) {
+float dW = min(prepared.metalMsFactor * bH(prepared.specularAlpha, lightLocal.z), 1.0 / max(lightLocal.z, 1e-8));
+response += prepared.weights[OPENPBR_LOBE_METAL] * prepared.metalMsScale * dW * INVPI;
+aH = mix(aH, 0.5 * INVPI, prepared.metalMsProbability);
+}
+pdf += prepared.probabilities[OPENPBR_LOBE_METAL] * aH;
+}
+}
+else {
+vec3 cK = vec3(cR(prepared, bQ, false));
+if (material.thinFilmWeight > 0.0) {
+cK = bW(
+cK,
+bQ,
+material.thinFilmThickness,
+material.thinFilmIor,
+prepared.thinFilmBaseIor,
+material.thinFilmWeight,
+prepared.thinFilmExteriorIor
+);
+}
+if (prepared.weights[OPENPBR_LOBE_TRANSMISSION] > 0.0) {
+if (material.thinWalled) {
+vec3 kM = vec3(lightLocal.xy, -lightLocal.z);
+vec3 dd = normalize(kM + viewLocal);
+if (dd.z < 0.0) {
+dd = -dd;
+}
+float nD = abs(dot(viewLocal, dd));
+vec3 lv = vec3(cR(prepared, nD, false));
+if (material.thinFilmWeight > 0.0) {
+lv = bW(
+lv,
+nD,
+material.thinFilmThickness,
+material.thinFilmIor,
+prepared.thinFilmBaseIor,
+material.thinFilmWeight,
+prepared.thinFilmExteriorIor
+);
+}
+vec3 gf = aP(
+lv,
+1.0
+);
+vec3 tint = bG(
+material,
+abs(viewLocal.z),
+prepared.refractionEtaTOverEtaI
+);
+response += prepared.weights[OPENPBR_LOBE_TRANSMISSION]
+* am(
+material,
+viewLocal,
+kM,
+dd,
+tint * (vec3(1.0) - gf),
+aH
+);
+}
+else {
+response += prepared.weights[OPENPBR_LOBE_TRANSMISSION]
+* dk(
+material,
+material.relativeIor,
+bn(material),
+viewLocal,
+lightLocal,
+V,
+cK,
+aH
+);
+}
+if (prepared.dielectricMsProbability > 0.0) {
+aH *= 1.0 - prepared.dielectricMsProbability;
+}
+pdf += prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] * aH;
+}
+if (prepared.weights[OPENPBR_LOBE_SUBSURFACE] > 0.0) {
+if (material.thinWalled) {
+vec3 nT = vec3(lightLocal.xy, -lightLocal.z);
+response += 0.5 * prepared.weights[OPENPBR_LOBE_SUBSURFACE]
+* cD(
+material,
+material.subsurfaceColor,
+viewLocal,
+nT,
+diffuseSpecularEnergyCompensation,
+aH
+);
+aH *= 0.5;
+}
+else {
+response += prepared.weights[OPENPBR_LOBE_SUBSURFACE] * dL(
+material,
+viewLocal,
+lightLocal,
+aH
+);
+}
+pdf += prepared.probabilities[OPENPBR_LOBE_SUBSURFACE] * aH;
+}
+}
+
+
+
+
+if (prepared.dielectricMsProbability > 0.0) {
+float au = abs(lightLocal.z);
+if (au > 0.0) {
+float ior = (reflection ? prepared.refractionEtaTOverEtaI : 1.0 / prepared.refractionEtaTOverEtaI);
+float iU = (reflection ? prepared.dielectricMsReflectionRatio : 1.0 - prepared.dielectricMsReflectionRatio);
+float td = prepared.dielectricMsViewLoss * bO(ior, prepared.specularAlpha, au)
+/ max(cO(ior, prepared.specularAlpha), 1e-12);
+float dW = min(td * iU, 1.0 / au);
+vec3 scale = (reflection ? prepared.dielectricMsReflectionScale : prepared.dielectricMsTransmissionScale);
+response += scale * dW * INVPI;
+pdf += prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] * prepared.dielectricMsProbability * iU * 0.5 * INVPI;
+}
+}
+
+#ifdef NUM_COAT_NORMAL_MAPS
+} 
+#endif
+if (material.coatWeight > 0.0) {
+response *= ds(prepared,
+#ifdef NUM_COAT_NORMAL_MAPS
+independentCoat ? coatLightCosine : lightLocal.z
+#else
+lightLocal.z
+#endif
+);
+}
+#ifdef NUM_COAT_NORMAL_MAPS
+if (independentCoat) { response *= abs(lightLocal.z); }
+#endif
+
+if (
+#ifdef NUM_COAT_NORMAL_MAPS
+(independentCoat ? coatReflection : reflection)
+#else
+reflection
+#endif
+&& prepared.weights[OPENPBR_LOBE_COAT] > 0.0) {
+vec3 dH = viewLocal;
+vec3 kt = lightLocal;
+vec3 gq = V;
+float jD = bQ;
+if (prepared.directionFrames) {
+dH = openpbr_to_local(prepared.coatBasis, openpbr_to_world(prepared.basis, viewLocal));
+#ifdef NUM_COAT_NORMAL_MAPS
+if (independentCoat) { dH = prepared.coatViewLocal; }
+#endif
+kt = openpbr_to_local(prepared.coatBasis, lightDirection);
+gq = normalize(kt + dH);
+if (gq.z < 0.0) { gq = -gq; }
+jD = abs(dot(dH, gq));
+}
+
+vec3 kU = vec3(af(
+jD,
+1.0 / material.coatIor
+));
+J dG = material;
+vec2 jp = aJ(material.coatRoughness, material.coatAnisotropy);
+dG.alphaX = jp.x;
+dG.alphaY = jp.y;
+response += prepared.weights[OPENPBR_LOBE_COAT]
+#ifdef NUM_COAT_NORMAL_MAPS
+* (independentCoat ? max(coatLightCosine, 0.0) : 1.0)
+#endif
+* am(
+dG,
+dH,
+kt,
+gq,
+kU,
+aH
+);
+pdf += prepared.probabilities[OPENPBR_LOBE_COAT] * aH;
+}
+
+if (
+#ifdef NUM_COAT_NORMAL_MAPS
+(independentCoat ? geometricReflection : reflection)
+#else
+reflection
+#endif
+&& prepared.weights[OPENPBR_LOBE_FUZZ] > 0.0) {
+float qn = aO(
+abs(viewLocal.z),
+material.fuzzRoughness
+);
+response *= max(
+vec3(0.0),
+vec3(1.0 - prepared.weights[OPENPBR_LOBE_FUZZ] * qn)
+);
+response += prepared.weights[OPENPBR_LOBE_FUZZ]
+#ifdef NUM_COAT_NORMAL_MAPS
+* (independentCoat ? max(lightLocal.z, 0.0) : 1.0)
+#endif
+* jP(
+material,
+viewLocal,
+lightLocal,
+aH
+);
+pdf += prepared.probabilities[OPENPBR_LOBE_FUZZ] * aH;
+}
+
+#ifdef NUM_COAT_NORMAL_MAPS
+if (independentCoat) { return response; }
+#endif
+return response * abs(lightLocal.z);
+}
+
+vec3 qU(
+vec3 tangent,
+J material,
+vec3 bF,
+vec3 normal,
+vec3 lightDirection,
+out float pdf
+) {
+ag prepared = eG(
+tangent,
+material,
+bF,
+normal
+);
+return bI(prepared, lightDirection, pdf);
+}
+
+float vJ(
+vec3 tangent,
+J material,
+vec3 bF,
+vec3 normal,
+vec3 lightDirection
+) {
+float pdf;
+qU(
+tangent,
+material,
+bF,
+normal,
+lightDirection,
+pdf
+);
+return pdf;
+}
+
+vec3 dv(
+ag prepared,
+vec3 cl,
+out vec3 lightDirection,
+out float pdf,
+out bool isTransmissionEvent
+) {
+J material = prepared.material;
+vec3 viewLocal = prepared.viewLocal;
+float cx = cl.x;
+float fp = cl.y;
+float bm = cl.z;
+pdf = 0.0;
+isTransmissionEvent = false;
+lightDirection = prepared.basis.normal;
+
+
+
+if (!prepared.material.thinWalled &&
+abs(prepared.material.relativeIor - 1.0) <= 1e-6 &&
+prepared.probabilities[OPENPBR_LOBE_TRANSMISSION] >= 1.0 - 1e-6) {
+lightDirection = openpbr_to_world(prepared.basis, -viewLocal);
+pdf = 1.0;
+isTransmissionEvent = true;
+return (bn(prepared.material) * ds(prepared,
+#ifdef NUM_COAT_NORMAL_MAPS
+prepared.material.coatNormalEnabled ? dot(prepared.coatBasis.normal, lightDirection) : -prepared.viewLocal.z
+#else
+-prepared.viewLocal.z
+#endif
+));
+}
+
+float cdf[OPENPBR_LOBE_COUNT];
+cdf[OPENPBR_LOBE_FUZZ] = prepared.probabilities[OPENPBR_LOBE_FUZZ];
+for (int hR = 1; hR < OPENPBR_LOBE_COUNT; hR++) {
+cdf[hR] = cdf[hR - 1] + prepared.probabilities[hR];
+}
+if (cdf[OPENPBR_LOBE_SUBSURFACE] <= 0.0) {
+return vec3(0.0);
+}
+
+vec3 lightLocal;
+if (bm < cdf[OPENPBR_LOBE_FUZZ] || (
+bm >= cdf[OPENPBR_LOBE_TRANSMISSION] && bm < cdf[OPENPBR_LOBE_DIFFUSE]
+)) {
+lightLocal = bi(vec2(cx, fp));
+}
+else if (bm < cdf[OPENPBR_LOBE_COAT]) {
+vec3 dH = viewLocal;
+if (prepared.directionFrames) { dH = openpbr_to_local(prepared.coatBasis, openpbr_to_world(prepared.basis, viewLocal)); }
+#ifdef NUM_COAT_NORMAL_MAPS
+if (material.coatNormalEnabled) { dH = prepared.coatViewLocal; }
+#endif
+vec3 V = aQ(
+dH,
+
+aJ(material.coatRoughness, material.coatAnisotropy),
+vec2(cx, fp)
+);
+if (V.z < 0.0) {
+V = -V;
+}
+lightLocal = normalize(reflect(-dH, V));
+#ifdef NUM_COAT_NORMAL_MAPS
+if (material.coatNormalEnabled) {
+lightDirection = openpbr_to_world(prepared.coatBasis, lightLocal);
+if (lightLocal.z <= 0.0 || dot(material.geometricNormal, lightDirection) <= 0.0) { return vec3(0.0); }
+return bI(prepared, lightDirection, pdf);
+}
+#endif
+if (prepared.directionFrames) { lightLocal = openpbr_to_local(prepared.basis, openpbr_to_world(prepared.coatBasis, lightLocal)); }
+}
+else if (bm < cdf[OPENPBR_LOBE_SPECULAR]) {
+bool tg = bm < cdf[OPENPBR_LOBE_METAL];
+float rw = (bm - cdf[OPENPBR_LOBE_COAT]) / max(prepared.probabilities[OPENPBR_LOBE_METAL], 1e-12);
+if (tg && rw < prepared.metalMsProbability) {
+float z = cx;
+float fi = sqrt(max(0.0, 1.0 - z * z));
+float br = TWOPI * fp;
+lightLocal = vec3(fi * cos(br), fi * sin(br), z);
+}
+else {
+vec3 V = aQ(
+viewLocal,
+vec2(material.alphaX, material.alphaY),
+vec2(cx, fp)
+);
+if (V.z < 0.0) {
+V = -V;
+}
+lightLocal = normalize(reflect(-viewLocal, V));
+}
+}
+else if (bm < cdf[OPENPBR_LOBE_TRANSMISSION]) {
+isTransmissionEvent = true;
+
+
+float ue = (bm - cdf[OPENPBR_LOBE_SPECULAR]) / max(prepared.probabilities[OPENPBR_LOBE_TRANSMISSION], 1e-12);
+if (prepared.dielectricMsProbability > 0.0 && ue < prepared.dielectricMsProbability) {
+float iU = prepared.dielectricMsReflectionRatio;
+bool lM = cx < iU;
+float z;
+if (lM) { z = min(cx / iU, 0.99999994); }
+else { z = -min((cx - iU) / (1.0 - iU), 0.99999994); }
+float fi = sqrt(max(0.0, 1.0 - z * z));
+float br = TWOPI * fp;
+lightLocal = vec3(fi * cos(br), fi * sin(br), z);
+isTransmissionEvent = !lM;
+}
+else {
+vec3 V = aQ(
+viewLocal,
+vec2(material.alphaX, material.alphaY),
+vec2(cx, fp)
+);
+if (V.z < 0.0) {
+V = -V;
+}
+float bQ = abs(dot(viewLocal, V));
+vec3 cK = vec3(af(
+bQ,
+prepared.material.relativeIor
+));
+if (material.thinFilmWeight > 0.0) {
+cK = bW(
+cK,
+bQ,
+material.thinFilmThickness,
+material.thinFilmIor,
+prepared.thinFilmBaseIor,
+material.thinFilmWeight,
+prepared.thinFilmExteriorIor
+);
+}
+if (material.thinWalled) {
+vec3 hW = normalize(reflect(-viewLocal, V));
+lightLocal = vec3(hW.xy, -hW.z);
+}
+else {
+lightLocal = refract(
+-viewLocal,
+V,
+prepared.material.relativeIor
+);
+if (dot(lightLocal, lightLocal) <= 1e-8) {
+isTransmissionEvent = false;
+return vec3(0.0);
+}
+lightLocal = normalize(lightLocal);
+}
+}
+}
+else {
+isTransmissionEvent = true;
+if (material.thinWalled) {
+float ky = cx;
+bool pH = ky >= 0.5;
+cx = pH
+? (ky - 0.5) * 2.0
+: ky * 2.0;
+lightLocal = bi(vec2(cx, fp));
+if (pH) {
+lightLocal.z = -lightLocal.z;
+}
+else {
+isTransmissionEvent = false;
+}
+}
+else {
+lightLocal = bi(vec2(cx, fp));
+lightLocal.z = -lightLocal.z;
+}
+}
+
+lightDirection = openpbr_to_world(prepared.basis, lightLocal);
+
+if (lightLocal.z == 0.0 || (lightLocal.z < 0.0) != isTransmissionEvent) {
+return vec3(0.0);
+}
+#ifdef NUM_COAT_NORMAL_MAPS
+if (material.coatNormalEnabled &&
+((dot(material.geometricNormal, lightDirection) <= 0.0) != isTransmissionEvent)) { return vec3(0.0); }
+#endif
+return bI(prepared, lightDirection, pdf);
+}
+
+vec3 su(
+vec3 tangent,
+J material,
+vec3 cl,
+vec3 bF,
+vec3 normal,
+out vec3 lightDirection,
+out float pdf,
+out bool isTransmissionEvent
+) {
+ag prepared = eG(tangent, material, bF, normal);
+return dv(prepared, cl, lightDirection, pdf, isTransmissionEvent);
+}`;
+//# sourceMappingURL=openPbrBsdf.glsl.js.map