@meimu/mutracer 0.2.0 → 0.3.0

package/lib/renderer/causticEmission.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"causticEmission.js","sourceRoot":"","sources":["../../src/renderer/causticEmission.ts"],"names":[],"mappings":"AACA,OAAO,EAAC,oBAAoB,EAAC,cAAc,EAAC,0BAA0B,EAAuB,MAAM,oBAAoB,CAAC;AACxH,OAAO,EAAC,gBAAgB,EAAC,oBAAoB,EAAC,MAAM,iBAAiB,CAAC;AAEtE,gEAAgE;AAChE,MAAM,UAAU,oBAAoB,CAAC,KAAwB,EAAC,MAAM,GAAC,CAAC;IACpE,IAAG,KAAK,CAAC,IAAI,KAAG,kBAAkB;QAAC,OAAO,cAAc,CAAC,KAAK,CAAC,KAAK,CAAC,GAAC,KAAK,CAAC,SAAS,GAAC,IAAI,CAAC,EAAE,GAAC,MAAM,GAAC,MAAM,CAAC;IAC5G,IAAG,KAAK,CAAC,IAAI,KAAG,WAAW;QAAC,OAAO,0BAA0B,CAAC,KAAK,CAAC,GAAC,CAAC,KAAK,CAAC,MAAM,IAAE,IAAI,CAAA,CAAC,CAAA,CAAC,GAAC,IAAI,CAAC,EAAE,CAAA,CAAC,CAAA,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAC,EAAE,GAAC,KAAK,CAAC,MAAM,CAAC,IAAE,CAAC,CAAC,CAAC;IACzI,IAAG,KAAK,CAAC,IAAI,KAAG,WAAW;QAAC,OAAO,0BAA0B,CAAC,KAAK,CAAC,CAAC;IACrE,MAAM,KAAK,GAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,GAAC,IAAI,CAAC,EAAE,GAAC,GAAG,CAAC,EAAC,KAAK,GAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,UAAU,EAAC,KAAK,CAAC,UAAU,CAAC,GAAC,IAAI,CAAC,EAAE,GAAC,GAAG,CAAC,CAAC;IAC3H,MAAM,IAAI,GAAC,GAAG,EAAC,OAAO,GAAC,GAAG,CAAC;IAAA,IAAI,QAAQ,GAAC,CAAC,CAAC;IAC1C,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,IAAI,EAAC,CAAC,EAAE;QAAC,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,OAAO,EAAC,CAAC,EAAE,EAAC,CAAC;YAChD,MAAM,MAAM,GAAC,CAAC,GAAC,CAAC,CAAC,GAAC,KAAK,CAAC,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,IAAI,CAAC;YACrC,MAAM,CAAC,GAAC,KAAK,GAAC,KAAK,CAAA,CAAC,CAAA,IAAI,CAAC,GAAG,CAAC,CAAC,EAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAC,CAAC,MAAM,GAAC,KAAK,CAAC,GAAC,CAAC,KAAK,GAAC,KAAK,CAAC,CAAC,CAAC,CAAA,CAAC,CAAA,CAAC,CAAC;YAC3E,MAAM,GAAG,GAAC,KAAK,CAAC,UAAU,EAAE,WAAW,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,EAAC,CAAC,GAAC,IAAI,CAAC,EAAE,GAAC,CAAC,CAAC,GAAC,EAAE,CAAC,GAAC,OAAO,GAAC,IAAI,CAAC,EAAE,CAAC,IAAE,CAAC,CAAC;YAC/F,QAAQ,IAAE,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAC,GAAG,CAAC;QAC5B,CAAC;IACD,QAAQ,IAAE,CAAC,GAAC,IAAI,CAAC,EAAE,GAAC,CAAC,CAAC,GAAC,KAAK,CAAC,GAAC,CAAC,IAAI,GAAC,OAAO,CAAC,CAAC;IAC7C,yEAAyE;IACzE,MAAM,SAAS,GAAC,cAAc,CAAC,KAAK,CAAC,KAAK,CAAC,GAAC,CAAC,KAAK,CAAC,MAAM,KAAG,CAAC,CAAA,CAAC,CAAA,KAAK,CAAC,SAAS,CAAA,CAAC,CAAA,oBAAoB,CAAC,KAAK,CAAC,GAAC,gBAAgB,CAAC,KAAK,CAAC,GAAC,CAAC,CAAC,CAAC;IACrI,OAAO,SAAS,GAAC,QAAQ,CAAC;AAC5B,CAAC;AACD,+EAA+E;AAC/E,MAAM,UAAU,oBAAoB,CAAC,MAAoC,EAAC,SAAqC,EAAC,MAAM,GAAC,CAAC,EAAC,gBAAgB,GAAC,CAAC;IACzI,MAAM,OAAO,GAAC,MAAM,CAAC,GAAG,CAAC,CAAC,KAAK,EAAC,WAAW,EAAC,EAAE,CAAA,CAAC,EAAC,UAAU,EAAC,CAAC,EAAC,WAAW,EAAC,MAAM,EAAC,oBAAoB,CAAC,KAAK,EAAC,MAAM,CAAC,EAAC,CAAC,CAAC,CAAC;IACtH,SAAS,CAAC,OAAO,CAAC,CAAC,QAAQ,EAAC,WAAW,EAAC,EAAE,CAAA,OAAO,CAAC,IAAI,CAAC,EAAC,UAAU,EAAC,CAAC,EAAC,WAAW,EAAC,MAAM,EAAC,IAAI,CAAC,EAAE,GAAC,QAAQ,CAAC,IAAI,GAAC,QAAQ,CAAC,iBAAiB,EAAC,CAAC,CAAC,CAAC;IAC5I,IAAG,gBAAgB,GAAC,CAAC;QAAC,OAAO,CAAC,IAAI,CAAC,EAAC,UAAU,EAAC,CAAC,EAAC,WAAW,EAAC,CAAC,EAAC,MAAM,EAAC,gBAAgB,GAAC,IAAI,CAAC,EAAE,GAAC,MAAM,GAAC,MAAM,EAAC,CAAC,CAAC;IAC/G,MAAM,YAAY,GAAC,oBAAoB,CAAC,OAAO,CAAC,EAAC,KAAK,GAAC,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC;IAChF,MAAM,IAAI,GAAC,IAAI,YAAY,CAAC,KAAK,GAAC,CAAC,CAAC,EAAC,KAAK,GAAC,IAAI,WAAW,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IACxE,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,KAAK,EAAC,CAAC,EAAE,EAAC,CAAC;QAAA,KAAK,CAAC,GAAG,CAAC,YAAY,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,EAAC,CAAC,GAAC,CAAC,CAAC,CAAC;QAAA,IAAI,CAAC,CAAC,GAAC,CAAC,GAAC,CAAC,CAAC,GAAC,YAAY,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAAA,CAAC;IACnH,OAAO,EAAC,IAAI,EAAC,KAAK,EAAC,UAAU,EAAC,YAAY,CAAC,WAAW,EAAC,CAAC;AAC1D,CAAC","sourcesContent":["import type {AnalyticLightState} from '../scene/lightState.js';\nimport {buildLightAliasTable,colorLuminance,estimateAnalyticLightPower,type EmissiveTriangle} from './lightSampling.js';\nimport {getSpotLightArea,getSpotLightRadiance} from './lightsData.js';\n\n/** Emitted flux, independent of the NEE selection heuristic. */\nexport function causticEmissionPower(light:AnalyticLightState,radius=1) {\n  if(light.type==='DirectionalLight')return colorLuminance(light.color)*light.intensity*Math.PI*radius*radius;\n  if(light.type==='RectLight')return estimateAnalyticLightPower(light)*(light.spread<=.001?1/Math.PI:Math.sin(Math.PI*.5*light.spread)**2);\n  if(light.type!=='SpotLight')return estimateAnalyticLightPower(light);\n  const outer=Math.cos(light.outerAngle*Math.PI/180),inner=Math.cos(Math.min(light.innerAngle,light.outerAngle)*Math.PI/180);\n  const rows=128,columns=128;let integral=0;\n  for(let y=0;y<rows;y++)for(let x=0;x<columns;x++){\n    const cosine=1-(1-outer)*(y+.5)/rows;\n    const t=inner>outer?Math.min(1,Math.max(0,(cosine-outer)/(inner-outer))):1;\n    const ies=light.iesProfile?.sampleBaked(Math.acos(cosine),2*Math.PI*(x+.5)/columns-Math.PI)??1;\n    integral+=t*t*(3-2*t)*ies;\n  }\n  integral*=2*Math.PI*(1-outer)/(rows*columns);\n  // A sphere projects to one quarter of its surface area in any direction.\n  const intensity=colorLuminance(light.color)*(light.radius===0?light.intensity:getSpotLightRadiance(light)*getSpotLightArea(light)/4);\n  return intensity*integral;\n}\n/** Two texels per alias choice: probability/alias/type/index, then its PMF. */\nexport function causticEmissionTable(lights:readonly AnalyticLightState[],triangles:readonly EmissiveTriangle[],radius=1,environmentPower=0) {\n  const entries=lights.map((light,sourceIndex)=>({sourceType:2,sourceIndex,weight:causticEmissionPower(light,radius)}));\n  triangles.forEach((triangle,sourceIndex)=>entries.push({sourceType:3,sourceIndex,weight:Math.PI*triangle.area*triangle.samplingLuminance}));\n  if(environmentPower>0)entries.push({sourceType:1,sourceIndex:0,weight:environmentPower*Math.PI*radius*radius});\n  const distribution=buildLightAliasTable(entries),count=distribution.pmfs.length;\n  const data=new Float32Array(count*8),words=new Uint32Array(data.buffer);\n  for(let i=0;i<count;i++){words.set(distribution.choices.subarray(i*4,i*4+4),i*8);data[i*8+4]=distribution.pmfs[i];}\n  return {data,count,totalPower:distribution.totalWeight};\n}\n"]}

package/lib/renderer/causticUniform.d.ts

@@ -0,0 +1,12 @@
+export interface CausticLightUniformOptions {
+    casterSphere: readonly number[];
+    viewProjection: Float32Array | Float64Array | number[];
+    photons: number;
+    batch: number;
+    cameraScale: number;
+    casterCount?: number;
+    choiceCount?: number;
+    orthographic?: boolean;
+    experiments?: number;
+}
+export declare function packCausticLightUniform(options: CausticLightUniformOptions): ArrayBuffer;

package/lib/renderer/causticUniform.js

@@ -0,0 +1,22 @@
+import { researchTransport } from './researchTransport.js';
+import { matrixStorage } from '../math/storage.js';
+import { CAUSTIC_UNIFORM_BYTES } from './wgsl/causticLight.compute.js';
+export function packCausticLightUniform(options) {
+    const data = new ArrayBuffer(CAUSTIC_UNIFORM_BYTES), f = new Float32Array(data), u = new Uint32Array(data);
+    const matrix = matrixStorage(options.viewProjection);
+    for (let i = 0; i < 4; i++)
+        f[i] = options.casterSphere[i];
+    for (let i = 0; i < 16; i++)
+        f[i + 4] = matrix[i];
+    u[20] = options.photons;
+    u[21] = options.batch;
+    u[22] = 0;
+    f[23] = options.cameraScale;
+    u[24] = options.casterCount ?? 0;
+    u[25] = options.choiceCount ?? 0;
+    u[26] = +!!options.orthographic;
+    u[27] = options.experiments ?? 0;
+    researchTransport()?.packUniform(data, options);
+    return data;
+}
+//# sourceMappingURL=causticUniform.js.map

package/lib/renderer/causticUniform.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"causticUniform.js","sourceRoot":"","sources":["../../src/renderer/causticUniform.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,iBAAiB,EAAC,MAAM,wBAAwB,CAAC;AACzD,OAAO,EAAC,aAAa,EAAC,MAAM,oBAAoB,CAAC;AACjD,OAAO,EAAC,qBAAqB,EAAC,MAAM,gCAAgC,CAAC;AAErE,MAAM,UAAU,uBAAuB,CAAC,OAAkC;IACxE,MAAM,IAAI,GAAC,IAAI,WAAW,CAAC,qBAAqB,CAAC,EAAC,CAAC,GAAC,IAAI,YAAY,CAAC,IAAI,CAAC,EAAC,CAAC,GAAC,IAAI,WAAW,CAAC,IAAI,CAAC,CAAC;IACnG,MAAM,MAAM,GAAC,aAAa,CAAC,OAAO,CAAC,cAAc,CAAC,CAAC;IAAA,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,CAAC,EAAC,CAAC,EAAE;QAAC,CAAC,CAAC,CAAC,CAAC,GAAC,OAAO,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC;IAAA,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,EAAE,EAAC,CAAC,EAAE;QAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1I,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,OAAO,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,KAAK,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,CAAC,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,WAAW,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,WAAW,IAAE,CAAC,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,WAAW,IAAE,CAAC,CAAC;IAAA,CAAC,CAAC,EAAE,CAAC,GAAC,CAAC,CAAC,CAAC,OAAO,CAAC,YAAY,CAAC;IACpK,CAAC,CAAC,EAAE,CAAC,GAAC,OAAO,CAAC,WAAW,IAAE,CAAC,CAAC;IAC7B,iBAAiB,EAAE,EAAE,WAAW,CAAC,IAAI,EAAE,OAAO,CAAC,CAAC;IAChD,OAAO,IAAI,CAAC;AACd,CAAC","sourcesContent":["import {researchTransport} from './researchTransport.js';\nimport {matrixStorage} from '../math/storage.js';\nimport {CAUSTIC_UNIFORM_BYTES} from './wgsl/causticLight.compute.js';\nexport interface CausticLightUniformOptions {casterSphere:readonly number[];viewProjection:Float32Array|Float64Array|number[];photons:number;batch:number;cameraScale:number;casterCount?:number;choiceCount?:number;orthographic?:boolean;experiments?:number}\nexport function packCausticLightUniform(options:CausticLightUniformOptions) {\n  const data=new ArrayBuffer(CAUSTIC_UNIFORM_BYTES),f=new Float32Array(data),u=new Uint32Array(data);\n  const matrix=matrixStorage(options.viewProjection);for(let i=0;i<4;i++)f[i]=options.casterSphere[i];for(let i=0;i<16;i++)f[i+4]=matrix[i];\n  u[20]=options.photons;u[21]=options.batch;u[22]=0;f[23]=options.cameraScale;u[24]=options.casterCount??0;u[25]=options.choiceCount??0;u[26]=+!!options.orthographic;\n  u[27]=options.experiments??0;\n  researchTransport()?.packUniform(data, options);\n  return data;\n}\n"]}

package/lib/renderer/caustics.d.ts

@@ -0,0 +1,16 @@
+import type Mesh from '../scene/MeshState.js';
+import type OpenPBRMaterial from '../scene/MaterialState.js';
+export declare const CAUSTIC_ROUGHNESS_MAX = 0.1;
+/** Material-only classification, shared by caster collection and GPU estimators. */
+export declare function isCausticSpecular(m: Pick<OpenPBRMaterial, 'baseWeight' | 'baseMetalness' | 'transmissionWeight' | 'subsurfaceWeight' | 'fuzzWeight' | 'coatWeight' | 'specularRoughness'>): boolean;
+export type CasterSphere = [number, number, number, number];
+/** One conservative world-space sphere per smooth mesh instance. */
+export declare function causticCasterSpheres(meshes: readonly Mesh[]): CasterSphere[];
+/** Conservative sphere enclosing the caster mixture. */
+export declare function causticCasterSphere(meshes: readonly Mesh[]): CasterSphere;
+/** CPU reference for the solid-angle mixture, including overlapping cones. */
+export declare function casterDirectionPdf(origin: readonly number[], direction: readonly number[], spheres: readonly CasterSphere[], weights?: readonly number[]): number;
+/** Stochastic smooth coverage currently disagrees with PT energy.
+ * Keep those scenes on PT until transport and endpoint weights are verified.
+ */
+export declare function hasStochasticCausticCaster(meshes: readonly Mesh[]): boolean;

package/lib/renderer/caustics.js

@@ -0,0 +1,93 @@
+import { computeGeometryBounds } from './geometryBvhBounds.js';
+import { matrixStorage } from '../math/storage.js';
+import { read, write } from '../math/access.js';
+import { meshInstanceCount, meshInstanceTransform } from './meshInstances.js';
+import { mat4, vec3 } from '../math.js';
+import { resolveMaterialAlphaMode } from '../scene/MaterialState.js';
+export const CAUSTIC_ROUGHNESS_MAX = 0.1;
+/** Material-only classification, shared by caster collection and GPU estimators. */
+export function isCausticSpecular(m) {
+    return m.baseWeight * (1 - m.baseMetalness) * (1 - m.transmissionWeight) <= 1e-6 &&
+        m.subsurfaceWeight <= 0 && m.fuzzWeight <= 0 && m.coatWeight <= 0 && m.specularRoughness <= CAUSTIC_ROUGHNESS_MAX;
+}
+const casterBounds = new WeakMap();
+function localCasterBounds(geometry) {
+    const cached = casterBounds.get(geometry), positions = geometry.position.array;
+    if (cached && cached.version === geometry.version && cached.positions === positions)
+        return cached.bounds;
+    const value = computeGeometryBounds(geometry), bounds = [matrixStorage(value[0]), matrixStorage(value[1])];
+    casterBounds.set(geometry, { version: geometry.version, positions, bounds });
+    return bounds;
+}
+/** One conservative world-space sphere per smooth mesh instance. */
+export function causticCasterSpheres(meshes) {
+    const spheres = [];
+    const matrix = mat4.create(), point = vec3.create();
+    for (const mesh of meshes) {
+        if (!mesh.materials.some(isCausticSpecular))
+            continue;
+        // Geometry.boundingBox is optional caller metadata and starts at zero.
+        // Caster collection must also work for a raw public Geometry descriptor.
+        const bounds = localCasterBounds(mesh.geometry), minCorner = bounds[0], maxCorner = bounds[1];
+        if (!Number.isFinite(minCorner[0]) || !Number.isFinite(maxCorner[0]))
+            continue;
+        let expansion = 0;
+        for (const material of mesh.materials)
+            if (isCausticSpecular(material))
+                expansion = Math.max(expansion, Math.abs(material.displacementBias) + (material.displacementMap ? Math.abs(material.displacementScale) : 0));
+        for (let instance = 0; instance < meshInstanceCount(mesh); instance++) {
+            const min = [Infinity, Infinity, Infinity], max = [-Infinity, -Infinity, -Infinity];
+            meshInstanceTransform(mesh, instance, matrix);
+            for (let corner = 0; corner < 8; corner++) {
+                for (let axis = 0; axis < 3; axis++)
+                    write(point, axis, ((corner & (1 << axis)) ? maxCorner[axis] + expansion : minCorner[axis] - expansion));
+                vec3.transformMat4(point, point, matrix);
+                for (let axis = 0; axis < 3; axis++) {
+                    min[axis] = Math.min(min[axis], read(point, axis));
+                    max[axis] = Math.max(max[axis], read(point, axis));
+                }
+            }
+            spheres.push([(min[0] + max[0]) * .5, (min[1] + max[1]) * .5, (min[2] + max[2]) * .5, Math.hypot(max[0] - min[0], max[1] - min[1], max[2] - min[2]) * .501]);
+        }
+    }
+    return spheres;
+}
+/** Conservative sphere enclosing the caster mixture. */
+export function causticCasterSphere(meshes) {
+    const spheres = causticCasterSpheres(meshes);
+    if (!spheres.length)
+        return [0, 0, 0, 0];
+    const min = [Infinity, Infinity, Infinity], max = [-Infinity, -Infinity, -Infinity];
+    for (const s of spheres) {
+        const components = [s[0], s[1], s[2]];
+        for (let i = 0; i < 3; i++) {
+            min[i] = Math.min(min[i], components[i] - s[3]);
+            max[i] = Math.max(max[i], components[i] + s[3]);
+        }
+    }
+    const center = min.map((v, i) => (v + max[i]) * .5);
+    let radius = 0;
+    for (const s of spheres)
+        radius = Math.max(radius, Math.hypot(s[0] - center[0], s[1] - center[1], s[2] - center[2]) + s[3]);
+    return [center[0], center[1], center[2], radius];
+}
+/** CPU reference for the solid-angle mixture, including overlapping cones. */
+export function casterDirectionPdf(origin, direction, spheres, weights = []) {
+    let total = 0, overlaps = 0;
+    for (const [index, s] of spheres.entries()) {
+        const weight = weights[index] ?? 1;
+        const d = Math.hypot(s[0] - origin[0], s[1] - origin[1], s[2] - origin[2]);
+        const cosine = d > s[3] ? Math.sqrt(Math.max(0, 1 - (s[3] / d) ** 2)) : -1;
+        total += 2 * Math.PI * (1 - cosine) * weight;
+        if (d <= s[3] || ((s[0] - origin[0]) * direction[0] + (s[1] - origin[1]) * direction[1] + (s[2] - origin[2]) * direction[2]) / d >= cosine)
+            overlaps += weight;
+    }
+    return total > 0 ? overlaps / total : 0;
+}
+/** Stochastic smooth coverage currently disagrees with PT energy.
+ * Keep those scenes on PT until transport and endpoint weights are verified.
+ */
+export function hasStochasticCausticCaster(meshes) {
+    return meshes.some(mesh => mesh.materials.some(m => isCausticSpecular(m) && resolveMaterialAlphaMode(m) === 'blend'));
+}
+//# sourceMappingURL=caustics.js.map

package/lib/renderer/caustics.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"caustics.js","sourceRoot":"","sources":["../../src/renderer/caustics.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,qBAAqB,EAAC,MAAM,wBAAwB,CAAC;AAE7D,OAAO,EAAC,aAAa,EAAC,MAAM,oBAAoB,CAAC;AACjD,OAAO,EAAC,IAAI,EAAC,KAAK,EAAC,MAAM,mBAAmB,CAAC;AAE7C,OAAO,EAAC,iBAAiB,EAAC,qBAAqB,EAAC,MAAM,oBAAoB,CAAC;AAC3E,OAAO,EAAC,IAAI,EAAC,IAAI,EAAC,MAAM,YAAY,CAAC;AAErC,OAAO,EAAC,wBAAwB,EAAC,MAAM,2BAA2B,CAAC;AAEnE,MAAM,CAAC,MAAM,qBAAqB,GAAG,GAAG,CAAC;AACzC,oFAAoF;AACpF,MAAM,UAAU,iBAAiB,CAAC,CAAwJ;IACxL,OAAO,CAAC,CAAC,UAAU,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,kBAAkB,CAAC,IAAI,IAAI;QAC9E,CAAC,CAAC,gBAAgB,IAAI,CAAC,IAAI,CAAC,CAAC,UAAU,IAAI,CAAC,IAAI,CAAC,CAAC,UAAU,IAAI,CAAC,IAAI,CAAC,CAAC,iBAAiB,IAAI,qBAAqB,CAAC;AACtH,CAAC;AAGD,MAAM,YAAY,GAAC,IAAI,OAAO,EAAuF,CAAC;AACtH,SAAS,iBAAiB,CAAC,QAAiB;IAC3C,MAAM,MAAM,GAAC,YAAY,CAAC,GAAG,CAAC,QAAQ,CAAC,EAAC,SAAS,GAAC,QAAQ,CAAC,QAAQ,CAAC,KAAK,CAAC;IAC1E,IAAG,MAAM,IAAE,MAAM,CAAC,OAAO,KAAG,QAAQ,CAAC,OAAO,IAAE,MAAM,CAAC,SAAS,KAAG,SAAS;QAAC,OAAO,MAAM,CAAC,MAAM,CAAC;IAChG,MAAM,KAAK,GAAC,qBAAqB,CAAC,QAAQ,CAAC,EAAC,MAAM,GAA6B,CAAC,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,EAAC,aAAa,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACjI,YAAY,CAAC,GAAG,CAAC,QAAQ,EAAC,EAAC,OAAO,EAAC,QAAQ,CAAC,OAAO,EAAC,SAAS,EAAC,MAAM,EAAC,CAAC,CAAC;IAAA,OAAO,MAAM,CAAC;AACtF,CAAC;AACD,oEAAoE;AACpE,MAAM,UAAU,oBAAoB,CAAC,MAAuB;IAC1D,MAAM,OAAO,GAAgB,EAAE,CAAC;IAChC,MAAM,MAAM,GAAC,IAAI,CAAC,MAAM,EAAE,EAAC,KAAK,GAAC,IAAI,CAAC,MAAM,EAAE,CAAC;IAC/C,KAAI,MAAM,IAAI,IAAI,MAAM,EAAE,CAAC;QACzB,IAAG,CAAC,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,iBAAiB,CAAC;YAAC,SAAS;QACpD,uEAAuE;QACvE,yEAAyE;QACzE,MAAM,MAAM,GAAC,iBAAiB,CAAC,IAAI,CAAC,QAAQ,CAAC,EAAC,SAAS,GAAC,MAAM,CAAC,CAAC,CAAC,EAAC,SAAS,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QACtF,IAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,IAAE,CAAC,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;YAAC,SAAS;QAC3E,IAAI,SAAS,GAAC,CAAC,CAAC;QAChB,KAAI,MAAM,QAAQ,IAAI,IAAI,CAAC,SAAS;YAAC,IAAG,iBAAiB,CAAC,QAAQ,CAAC;gBAAC,SAAS,GAAC,IAAI,CAAC,GAAG,CAAC,SAAS,EAC9F,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,gBAAgB,CAAC,GAAC,CAAC,QAAQ,CAAC,eAAe,CAAA,CAAC,CAAA,IAAI,CAAC,GAAG,CAAC,QAAQ,CAAC,iBAAiB,CAAC,CAAA,CAAC,CAAA,CAAC,CAAC,CAAC,CAAC;QACzG,KAAI,IAAI,QAAQ,GAAC,CAAC,EAAC,QAAQ,GAAC,iBAAiB,CAAC,IAAI,CAAC,EAAC,QAAQ,EAAE,EAAE,CAAC;YAC/D,MAAM,GAAG,GAAC,CAAC,QAAQ,EAAC,QAAQ,EAAC,QAAQ,CAAC,EAAC,GAAG,GAAC,CAAC,CAAC,QAAQ,EAAC,CAAC,QAAQ,EAAC,CAAC,QAAQ,CAAC,CAAC;YAC3E,qBAAqB,CAAC,IAAI,EAAC,QAAQ,EAAC,MAAM,CAAC,CAAC;YAC5C,KAAI,IAAI,MAAM,GAAC,CAAC,EAAC,MAAM,GAAC,CAAC,EAAC,MAAM,EAAE,EAAE,CAAC;gBACnC,KAAI,IAAI,IAAI,GAAC,CAAC,EAAC,IAAI,GAAC,CAAC,EAAC,IAAI,EAAE;oBAAC,KAAK,CAAC,KAAK,EAAC,IAAI,EAAC,CAAC,CAAC,MAAM,GAAC,CAAC,CAAC,IAAE,IAAI,CAAC,CAAC,CAAA,CAAC,CAAA,SAAS,CAAC,IAAI,CAAC,GAAC,SAAS,CAAA,CAAC,CAAA,SAAS,CAAC,IAAI,CAAC,GAAC,SAAS,CAAC,CAAC,CAAC;gBACxH,IAAI,CAAC,aAAa,CAAC,KAAK,EAAC,KAAK,EAAC,MAAM,CAAC,CAAC;gBACvC,KAAI,IAAI,IAAI,GAAC,CAAC,EAAC,IAAI,GAAC,CAAC,EAAC,IAAI,EAAE,EAAC,CAAC;oBAAA,GAAG,CAAC,IAAI,CAAC,GAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,EAAC,IAAI,CAAC,KAAK,EAAC,IAAI,CAAC,CAAC,CAAC;oBAAA,GAAG,CAAC,IAAI,CAAC,GAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,EAAC,IAAI,CAAC,KAAK,EAAC,IAAI,CAAC,CAAC,CAAC;gBAAA,CAAC;YAC/H,CAAC;YACD,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAC,EAAE,EAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAC,EAAE,EAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAC,EAAE,EAAC,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,EAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,EAAC,GAAG,CAAC,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAC,IAAI,CAAC,CAAC,CAAC;QACtI,CAAC;IACH,CAAC;IACD,OAAO,OAAO,CAAC;AACjB,CAAC;AACD,wDAAwD;AACxD,MAAM,UAAU,mBAAmB,CAAC,MAAsB;IACxD,MAAM,OAAO,GAAC,oBAAoB,CAAC,MAAM,CAAC,CAAC;IAAA,IAAG,CAAC,OAAO,CAAC,MAAM;QAAC,OAAO,CAAC,CAAC,EAAC,CAAC,EAAC,CAAC,EAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,GAAG,GAAC,CAAC,QAAQ,EAAC,QAAQ,EAAC,QAAQ,CAAC,EAAC,GAAG,GAAC,CAAC,CAAC,QAAQ,EAAC,CAAC,QAAQ,EAAC,CAAC,QAAQ,CAAC,CAAC;IAC3E,KAAI,MAAM,CAAC,IAAI,OAAO,EAAC,CAAC;QAAA,MAAM,UAAU,GAAU,CAAC,CAAC,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAAA,KAAI,IAAI,CAAC,GAAC,CAAC,EAAC,CAAC,GAAC,CAAC,EAAC,CAAC,EAAE,EAAC,CAAC;YAAA,GAAG,CAAC,CAAC,CAAC,GAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAC,UAAU,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAAA,GAAG,CAAC,CAAC,CAAC,GAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAC,UAAU,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAAA,CAAC;IAAA,CAAC;IAChL,MAAM,MAAM,GAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,EAAC,CAAC,EAAC,EAAE,CAAA,CAAC,CAAC,GAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAC,EAAE,CAAC,CAAC;IAC3C,IAAI,MAAM,GAAC,CAAC,CAAC;IAAA,KAAI,MAAM,CAAC,IAAI,OAAO;QAAC,MAAM,GAAC,IAAI,CAAC,GAAG,CAAC,MAAM,EAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC1H,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAC,MAAM,CAAC,CAAC,CAAC,EAAC,MAAM,CAAC,CAAC,CAAC,EAAC,MAAM,CAAC,CAAC;AAChD,CAAC;AACD,8EAA8E;AAC9E,MAAM,UAAU,kBAAkB,CAAC,MAAwB,EAAC,SAA2B,EAAC,OAA+B,EAAC,UAA0B,EAAE;IAClJ,IAAI,KAAK,GAAC,CAAC,EAAC,QAAQ,GAAC,CAAC,CAAC;IACvB,KAAI,MAAM,CAAC,KAAK,EAAC,CAAC,CAAC,IAAI,OAAO,CAAC,OAAO,EAAE,EAAC,CAAC;QAAA,MAAM,MAAM,GAAC,OAAO,CAAC,KAAK,CAAC,IAAE,CAAC,CAAC;QAAA,MAAM,CAAC,GAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,EAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;QACxI,MAAM,MAAM,GAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAA,CAAC,CAAA,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,IAAE,CAAC,CAAC,CAAC,CAAA,CAAC,CAAA,CAAC,CAAC,CAAC;QAC5D,KAAK,IAAE,CAAC,GAAC,IAAI,CAAC,EAAE,GAAC,CAAC,CAAC,GAAC,MAAM,CAAC,GAAC,MAAM,CAAC;QACnC,IAAG,CAAC,IAAE,CAAC,CAAC,CAAC,CAAC,IAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAC,SAAS,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAC,SAAS,CAAC,CAAC,CAAC,GAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAC,MAAM,CAAC,CAAC,CAAC,CAAC,GAAC,SAAS,CAAC,CAAC,CAAC,CAAC,GAAC,CAAC,IAAE,MAAM;YAAC,QAAQ,IAAE,MAAM,CAAC;IACrI,CAAC;IACD,OAAO,KAAK,GAAC,CAAC,CAAA,CAAC,CAAA,QAAQ,GAAC,KAAK,CAAA,CAAC,CAAA,CAAC,CAAC;AAClC,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,0BAA0B,CAAC,MAAsB;IAChE,OAAO,MAAM,CAAC,IAAI,CAAC,IAAI,CAAA,EAAE,CAAA,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,CAAC,CAAA,EAAE,CAAA,iBAAiB,CAAC,CAAC,CAAC,IAAE,wBAAwB,CAAC,CAAC,CAAC,KAAG,OAAO,CAAC,CAAC,CAAC;AAC/G,CAAC","sourcesContent":["import {computeGeometryBounds} from './geometryBvhBounds.js';\nimport type Geometry from '../scene/Geometry.js';\nimport {matrixStorage} from '../math/storage.js';\nimport {read,write} from '../math/access.js';\nimport type Mesh from '../scene/MeshState.js';\nimport {meshInstanceCount,meshInstanceTransform} from './meshInstances.js';\nimport {mat4,vec3} from '../math.js';\nimport type OpenPBRMaterial from '../scene/MaterialState.js';\nimport {resolveMaterialAlphaMode} from '../scene/MaterialState.js';\n\nexport const CAUSTIC_ROUGHNESS_MAX = 0.1;\n/** Material-only classification, shared by caster collection and GPU estimators. */\nexport function isCausticSpecular(m: Pick<OpenPBRMaterial, 'baseWeight' | 'baseMetalness' | 'transmissionWeight' | 'subsurfaceWeight' | 'fuzzWeight' | 'coatWeight' | 'specularRoughness'>) {\n  return m.baseWeight * (1 - m.baseMetalness) * (1 - m.transmissionWeight) <= 1e-6 &&\n    m.subsurfaceWeight <= 0 && m.fuzzWeight <= 0 && m.coatWeight <= 0 && m.specularRoughness <= CAUSTIC_ROUGHNESS_MAX;\n}\n\nexport type CasterSphere = [number,number,number,number];\nconst casterBounds=new WeakMap<Geometry,{version:number;positions:Float32Array;bounds:[Float64Array,Float64Array]}>();\nfunction localCasterBounds(geometry:Geometry){\n const cached=casterBounds.get(geometry),positions=geometry.position.array;\n if(cached&&cached.version===geometry.version&&cached.positions===positions)return cached.bounds;\n const value=computeGeometryBounds(geometry),bounds:[Float64Array,Float64Array]=[matrixStorage(value[0]),matrixStorage(value[1])];\n casterBounds.set(geometry,{version:geometry.version,positions,bounds});return bounds;\n}\n/** One conservative world-space sphere per smooth mesh instance. */\nexport function causticCasterSpheres(meshes: readonly Mesh[]): CasterSphere[] {\n  const spheres:CasterSphere[]=[];\n  const matrix=mat4.create(),point=vec3.create();\n  for(const mesh of meshes) {\n    if(!mesh.materials.some(isCausticSpecular))continue;\n    // Geometry.boundingBox is optional caller metadata and starts at zero.\n    // Caster collection must also work for a raw public Geometry descriptor.\n    const bounds=localCasterBounds(mesh.geometry),minCorner=bounds[0],maxCorner=bounds[1];\n    if(!Number.isFinite(minCorner[0])||!Number.isFinite(maxCorner[0]))continue;\n    let expansion=0;\n    for(const material of mesh.materials)if(isCausticSpecular(material))expansion=Math.max(expansion,\n      Math.abs(material.displacementBias)+(material.displacementMap?Math.abs(material.displacementScale):0));\n    for(let instance=0;instance<meshInstanceCount(mesh);instance++) {\n      const min=[Infinity,Infinity,Infinity],max=[-Infinity,-Infinity,-Infinity];\n      meshInstanceTransform(mesh,instance,matrix);\n      for(let corner=0;corner<8;corner++) {\n        for(let axis=0;axis<3;axis++)write(point,axis,((corner&(1<<axis))?maxCorner[axis]+expansion:minCorner[axis]-expansion));\n        vec3.transformMat4(point,point,matrix);\n        for(let axis=0;axis<3;axis++){min[axis]=Math.min(min[axis],read(point,axis));max[axis]=Math.max(max[axis],read(point,axis));}\n      }\n      spheres.push([(min[0]+max[0])*.5,(min[1]+max[1])*.5,(min[2]+max[2])*.5,Math.hypot(max[0]-min[0],max[1]-min[1],max[2]-min[2])*.501]);\n    }\n  }\n  return spheres;\n}\n/** Conservative sphere enclosing the caster mixture. */\nexport function causticCasterSphere(meshes:readonly Mesh[]):CasterSphere {\n  const spheres=causticCasterSpheres(meshes);if(!spheres.length)return [0,0,0,0];\n  const min=[Infinity,Infinity,Infinity],max=[-Infinity,-Infinity,-Infinity];\n  for(const s of spheres){const components:number[]=[s[0],s[1],s[2]];for(let i=0;i<3;i++){min[i]=Math.min(min[i],components[i]-s[3]);max[i]=Math.max(max[i],components[i]+s[3]);}}\n  const center=min.map((v,i)=>(v+max[i])*.5);\n  let radius=0;for(const s of spheres)radius=Math.max(radius,Math.hypot(s[0]-center[0],s[1]-center[1],s[2]-center[2])+s[3]);\n  return [center[0],center[1],center[2],radius];\n}\n/** CPU reference for the solid-angle mixture, including overlapping cones. */\nexport function casterDirectionPdf(origin:readonly number[],direction:readonly number[],spheres:readonly CasterSphere[],weights:readonly number[]=[]) {\n  let total=0,overlaps=0;\n  for(const [index,s] of spheres.entries()){const weight=weights[index]??1;const d=Math.hypot(s[0]-origin[0],s[1]-origin[1],s[2]-origin[2]);\n    const cosine=d>s[3]?Math.sqrt(Math.max(0,1-(s[3]/d)**2)):-1;\n    total+=2*Math.PI*(1-cosine)*weight;\n    if(d<=s[3]||((s[0]-origin[0])*direction[0]+(s[1]-origin[1])*direction[1]+(s[2]-origin[2])*direction[2])/d>=cosine)overlaps+=weight;\n  }\n  return total>0?overlaps/total:0;\n}\n\n/** Stochastic smooth coverage currently disagrees with PT energy.\n * Keep those scenes on PT until transport and endpoint weights are verified.\n */\nexport function hasStochasticCausticCaster(meshes:readonly Mesh[]){\n return meshes.some(mesh=>mesh.materials.some(m=>isCausticSpecular(m)&&resolveMaterialAlphaMode(m)==='blend'));\n}\n"]}

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

@@ -9,7 +9,7 @@ export default new FragmentShader({
         mapping: Shader.uniform('vec4'), idComponent: Shader.uniform('int')
     },
     main: /* glsl */ `
-uint vn(uint v) {
+uint vB(uint v) {
 v = (v ^ (v >> 16u)) * 0x7feb352du;
 v = (v ^ (v >> 15u)) * 0x846ca68bu;
 return v ^ (v >> 16u);
@@ -36,16 +36,16 @@ color = srgb(clamp(r.xyz / max(r.w, 1.0), 0.0, 1.0));
 } else if (kind == 5) {
 uvec4 ids = texelFetch(identity, p, 0);
 alpha = ids.x == 0xffffffffu ? 0.0 : 1.0;
-uint h = vn(ids[idComponent]);
+uint h = vB(ids[idComponent]);
 color = .2 + .8 * vec3(float(h & 255u), float((h >> 8u) & 255u),
 float((h >> 16u) & 255u)) / 255.0;
 } else {
 vec4 m = texelFetch(motion, p, 0);
 alpha = m.w == 3.0 ? 1.0 : 0.0;
-float gR = length(m.xy), hue = 0.0;
-if (gR > 0.0) hue = fract(atan(m.y, m.x) / 6.28318530718 + 1.0);
+float gz = length(m.xy), hue = 0.0;
+if (gz > 0.0) hue = fract(atan(m.y, m.x) / 6.28318530718 + 1.0);
 vec3 v = abs(fract(vec3(hue) + vec3(0.0, 2.0 / 3.0, 1.0 / 3.0)) * 6.0 - 3.0);
-color = clamp(v - 1.0, 0.0, 1.0) * clamp(gR / mapping.w, 0.0, 1.0);
+color = clamp(v - 1.0, 0.0, 1.0) * clamp(gz / mapping.w, 0.0, 1.0);
 }
 if (alpha <= 0.0) color = vec3(0.0);
 out_color = vec4(color, clamp(alpha, 0.0, 1.0));

package/lib/renderer/glsl/aovReflectance.glsl.js

@@ -1,82 +1,82 @@
 import proposals from './aovReflectanceSamples.glsl.js';
 // GLSL mirror of the shared 256-proposal, fused OpenPBR reflection quadrature.
 export default proposals + /* glsl */ `
-void aovReflectance(ag prepared, bool qQ, bool kY, out vec3 diffuse, out vec3 specular) {
+void aovReflectance(OpenPbrPreparedBsdf prepared, bool qH, bool kE, out vec3 diffuse, out vec3 specular) {
 diffuse=vec3(0.0); specular=vec3(0.0);
 vec3 view = prepared.viewLocal;
 if (view.z <= 1e-6) { return; }
 OpenPbrResolvedMaterial material = prepared.material;
 vec2 alpha = vec2(material.alphaX, material.alphaY);
-float mg = aF(view, alpha);
-float jl = prepared.weights[OPENPBR_LOBE_DIFFUSE];
-float gb = prepared.weights[OPENPBR_LOBE_METAL];
+float lN = aD(view, alpha);
+float iR = prepared.weights[OPENPBR_LOBE_DIFFUSE];
+float fN = prepared.weights[OPENPBR_LOBE_METAL];
 float specularWeight = prepared.weights[OPENPBR_LOBE_SPECULAR];
-float ei = 1.0;
-if (kY && prepared.dielectricMsProbability > 0.0) {
-ei = max(bf(prepared.refractionEtaTOverEtaI, prepared.specularAlpha), 1e-12);
+float dS = 1.0;
+if (kE && prepared.dielectricMsProbability > 0.0) {
+dS = max(aY(prepared.refractionEtaTOverEtaI, prepared.specularAlpha), 1e-12);
 }
-vec3 cr = normalize(vec3(alpha, 1.0) * view);
-float qP = min(cr.z, 1.0);
-float kh = 1.0 + qP;
+vec3 ce = normalize(vec3(alpha, 1.0) * view);
+float qG = min(ce.z, 1.0);
+float jN = 1.0 + qG;
 uint count = 256u;
 for (uint i = 0u; i < uint(NO_UNROLL(256)); i++) {
-uint vV = i / 2u;
-vec2 u = vec2((float(vV) + 0.5) / float(count / 2u),
-fract(float(reverseBits(vV)) * 2.3283064365386963e-10 + 0.37));
-vec4 kK = hp[vV];
+uint wn = i / 2u;
+vec2 u = vec2((float(wn) + 0.5) / float(count / 2u),
+fract(float(reverseBits(wn)) * 2.3283064365386963e-10 + 0.37));
+vec4 proposal = gV[wn];
 
-bool hO = prepared.specularAlpha < 0.001;
-vec3 light = vec3(kK.xy, sqrt(max(0.0, 1.0 - dot(kK.xy, kK.xy))));
-if (hO) { light = aO(u); }
+bool hv = prepared.specularAlpha < 0.001;
+vec3 light = vec3(proposal.xy, sqrt(max(0.0, 1.0 - dot(proposal.xy, proposal.xy))));
+if (hv) { light = aM(u); }
 if ((i & 1u) != 0u) {
-float kC = (1.0 - u.y) * kh;
-float ex = sqrt(max(0.0, (kC + (1.0 - qP)) * (u.y * kh)));
-vec3 ad = normalize(vec3(alpha, 1.0) * vec3(cr.xy + ex * kK.zw, kC));
-if (hO) { ad = aD(view, alpha, u); }
-light = reflect(-view, ad);
+float kf = (1.0 - u.y) * jN;
+float dr = sqrt(max(0.0, (kf + (1.0 - qG)) * (u.y * jN)));
+vec3 aa = normalize(vec3(alpha, 1.0) * vec3(ce.xy + dr * proposal.zw, kf));
+if (hv) { aa = aA(view, alpha, u); }
+light = reflect(-view, aa);
 }
 
 if (light.z <= 0.0) { continue; }
-vec3 ad = normalize(view + light);
-float qH = mg * cW(ad, alpha) / (4.0 * view.z);
-float pdf = 0.5 * (light.z * INVPI + qH);
+vec3 aa = normalize(view + light);
+float specularPdf = lN * sampling_ggx_ndf(aa, alpha) / (4.0 * view.z);
+float pdf = 0.5 * (light.z * INVPI + specularPdf);
 
 
-vec3 aA = openpbr_to_local(prepared.basis, openpbr_to_world(prepared.basis, light));
-if (aA.z <= 0.0) { continue; }
-float scale = abs(aA.z) / (pdf * float(count));
-if (qQ && jl > 0.0) {
-float compensation = prepared.diffuseSpecularEnergyCompensation * aw(
-prepared.specularEtaTOverEtaI, prepared.specularAlpha, abs(aA.z));
-float sZ = 0.0;
-diffuse += jl * bJ(material, prepared.diffuseColor,
-view, aA, compensation, sZ) * scale;
+vec3 au = openpbr_to_local(prepared.basis, openpbr_to_world(prepared.basis, light));
+if (au.z <= 0.0) { continue; }
+float scale = abs(au.z) / (pdf * float(count));
+if (qH && iR > 0.0) {
+float compensation = prepared.diffuseSpecularEnergyCompensation * ar(
+prepared.specularEtaTOverEtaI, prepared.specularAlpha, abs(au.z));
+float unusedPdf = 0.0;
+diffuse += iR * bw(material, prepared.diffuseColor,
+view, au, compensation, unusedPdf) * scale;
 }
-if (kY) {
-vec3 jn = normalize(view + aA);
-float bm = abs(dot(view, jn));
-vec3 bG = vec3(0.0);
+if (kE) {
+vec3 evaluatedHalf = normalize(view + au);
+float viewHalf = abs(dot(view, evaluatedHalf));
+vec3 fresnel = vec3(0.0);
 if (specularWeight > 0.0) {
-vec3 dielectric = vec3(af(bm, material.relativeIor));
-if (material.thinWalled) { dielectric = av(dielectric, material.transmissionWeight); }
-bG += specularWeight * material.specularColor * dielectric;
+vec3 dielectric = vec3(ae(viewHalf, material.relativeIor));
+if (material.thinWalled) { dielectric = aq(dielectric, material.transmissionWeight); }
+fresnel += specularWeight * material.specularColor * dielectric;
 }
-if (gb > 0.0) {
-bG += gb * material.specularWeight * cw(bm,
+if (fN > 0.0) {
+fresnel += fN * material.specularWeight * openpbr_metal_fresnel(viewHalf,
 material.baseWeight * material.baseColor, material.specularColor);
 }
-float rJ = cW(jn, alpha) * mg
-* aF(aA, alpha) / (4.0 * max(aA.z * view.z, 1e-8));
-vec3 reflection = bG * rJ;
-if (gb > 0.0 && prepared.metalMsFactor > 0.0) {
-float cL = min(prepared.metalMsFactor * aV(prepared.specularAlpha, aA.z), 1.0 / max(aA.z, 1e-8));
-reflection += gb * prepared.metalMsScale * cL * INVPI;
+float rE = sampling_ggx_ndf(evaluatedHalf, alpha) * lN
+* aD(au, alpha) / (4.0 * max(au.z * view.z, 1e-8));
+vec3 reflection = fresnel * rE;
+if (fN > 0.0 && prepared.metalMsFactor > 0.0) {
+float cv = min(prepared.metalMsFactor * aR(prepared.specularAlpha, au.z), 1.0 / max(au.z, 1e-8));
+reflection += fN * prepared.metalMsScale * cv * INVPI;
 }
 if (prepared.dielectricMsProbability > 0.0) {
-float oZ = prepared.dielectricMsViewLoss * bb(
-prepared.refractionEtaTOverEtaI, prepared.specularAlpha, abs(aA.z)) / ei;
-float cL = min(oZ * prepared.dielectricMsReflectionRatio, 1.0 / abs(aA.z));
-reflection += prepared.dielectricMsReflectionScale * cL * INVPI;
+float oP = prepared.dielectricMsViewLoss * aV(
+prepared.refractionEtaTOverEtaI, prepared.specularAlpha, abs(au.z)) / dS;
+float cv = min(oP * prepared.dielectricMsReflectionRatio, 1.0 / abs(au.z));
+reflection += prepared.dielectricMsReflectionScale * cv * INVPI;
 }
 specular += reflection * scale;
 }

package/lib/renderer/glsl/aovReflectanceSamples.glsl.js

@@ -1,7 +1,7 @@
 // Generated by scripts/generate-aov-reflectance-samples.mjs; do not edit.
 // xy: concentric cosine disk; zw: spherical-cap azimuth cosine/sine.
 export default /* glsl */ `
-const vec4 hp[128] = vec4[128](
+const vec4 gV[128] = vec4[128](
 vec4(-0.971247852, -0.202764958, 0.999698818, 0.0245412290),
 vec4(-0.808660269, 0.547487736, 0.997290432, 0.0735645667),
 vec4(-0.942509234, 0.187289074, 0.992479563, 0.122410677),

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

@@ -11,19 +11,19 @@ import random from './chunks/random.glsl.js';
 import cameraRay from './chunks/cameraRay.glsl.js';
 import rand from './chunks/rand.glsl.js';
 const surface = /* glsl */ `
-uint ke(int index) { ivec2 size = textureSize(aovScene,0); uvec4 v = texelFetch(aovScene, ivec2((index/4)%size.x, (index/4)/size.x),0); return v[index%4]; }
-vec3 iN(int offset){return vec3(uintBitsToFloat(ke(offset)),uintBitsToFloat(ke(offset+1)),uintBitsToFloat(ke(offset+2)));}
-mat4 ow(int offset) { mat4 m; for (int i=0;i<4;i++) for(int j=0;j<4;j++) m[i][j]=uintBitsToFloat(ke(offset+i*4+j)); return m; }
-Ray wJ(vec2 uv) {
+uint jK(int index) { ivec2 size = textureSize(aovScene,0); uvec4 v = texelFetch(aovScene, ivec2((index/4)%size.x, (index/4)/size.x),0); return v[index%4]; }
+vec3 iu(int offset){return vec3(uintBitsToFloat(jK(offset)),uintBitsToFloat(jK(offset+1)),uintBitsToFloat(jK(offset+2)));}
+mat4 oi(int offset) { mat4 m; for (int i=0;i<4;i++) for(int j=0;j<4;j++) m[i][j]=uintBitsToFloat(jK(offset+i*4+j)); return m; }
+Ray xj(vec2 uv) {
 vec4 origin = viewProjectionInv*vec4(uv*2.0-1.0,-1.0,1.0);
 vec4 target = viewProjectionInv*vec4(uv*2.0-1.0,1.0,1.0);
-Ray ray; gC(ray, origin.xyz/origin.w, normalize(target.xyz*origin.w-origin.xyz*target.w)); return ray;
+Ray ray; initRay(ray, origin.xyz/origin.w, normalize(target.xyz*origin.w-origin.xyz*target.w)); return ray;
 }
-bool mQ(Ray wX, out SurfaceInteraction si, out OpenPbrResolvedMaterial material) {
-Ray ray=wX; si.triangleIndex=-1;
+bool mC(Ray xA, out SurfaceInteraction si, out OpenPbrResolvedMaterial material) {
+Ray ray=xA; si.triangleIndex=-1;
 for(int layer=0;layer<32;layer++) {
 intersectScene(ray,si,false); if(!si.hit) return false;
-dg(material,si.materialIndex,si.uv,si.uv1,si.instanceIndex,si.vertexColor);
+getMaterial(material,si.materialIndex,si.uv,si.uv1,si.instanceIndex,si.vertexColor);
 if(material.materialType!=SHADOW_CATCHER && material.alpha>=0.5) return true;
 ray.o=si.position+ray.d*(EPS*0.02);
 }
@@ -57,56 +57,56 @@ ${surface}
 void main() {
 ${center ? /* glsl */ `
 out_id=uvec4(0xffffffffu); out_motion=vec4(0.0);
-Ray ray=wJ((floor(gl_FragCoord.xy)+0.5)/viewportSize);
+Ray ray=xj((floor(gl_FragCoord.xy)+0.5)/viewportSize);
 SurfaceInteraction si; OpenPbrResolvedMaterial material;
-if(!mQ(ray,si,material)) return;
+if(!mC(ray,si,material)) return;
 int offset=52+si.instanceIndex*24;
-out_id=uvec4(ke(offset),ke(offset+1),ke(int(ke(16))+si.materialIndex),ke(int(ke(17))+si.triangleIndex));
-out_motion.z=dot(si.position-iN(36),iN(40));
-if(ke(offset+2)==0u)return;
-vec4 qA=ow(offset+4)*vec4(si.position,1.0);
-vec4 ps=ow(0)*qA;
-if(ps.w<=0.0)return;
-vec3 tX=ps.xyz/ps.w;vec2 iz=(tX.xy*0.5+0.5)*viewportSize;
-if(any(lessThan(iz,vec2(0.0)))||any(greaterThanEqual(iz,viewportSize))||abs(tX.z)>1.0)return;
-uvec4 identity=texelFetch(previousId,ivec2(iz),0);float qz=texelFetch(previousMotion,ivec2(iz),0).z;
-if(any(notEqual(identity.xy,out_id.xy))||abs(qz-dot(qA.xyz-iN(44),iN(48)))>max(1e-4,abs(qz)*1e-3))return;
-out_motion=vec4((iz-gl_FragCoord.xy)*vec2(1.0,-1.0),out_motion.z,3.0);
+out_id=uvec4(jK(offset),jK(offset+1),jK(int(jK(16))+si.materialIndex),jK(int(jK(17))+si.triangleIndex));
+out_motion.z=dot(si.position-iu(36),iu(40));
+if(jK(offset+2)==0u)return;
+vec4 qs=oi(offset+4)*vec4(si.position,1.0);
+vec4 pg=oi(0)*qs;
+if(pg.w<=0.0)return;
+vec3 ua=pg.xyz/pg.w;vec2 hU=(ua.xy*0.5+0.5)*viewportSize;
+if(any(lessThan(hU,vec2(0.0)))||any(greaterThanEqual(hU,viewportSize))||abs(ua.z)>1.0)return;
+uvec4 identity=texelFetch(previousId,ivec2(hU),0);float qr=texelFetch(previousMotion,ivec2(hU),0).z;
+if(any(notEqual(identity.xy,out_id.xy))||abs(qr-dot(qs.xyz-iu(44),iu(48)))>max(1e-4,abs(qr)*1e-3))return;
+out_motion=vec4((hU-gl_FragCoord.xy)*vec2(1.0,-1.0),out_motion.z,3.0);
 ` : /* glsl */ `
 ivec2 coord=ivec2(gl_FragCoord.xy);
 out_geometry=aovSample==0 ? vec4(0.0) : texelFetch(previousGeometry,coord,0);
 out_position=aovSample==0 ? vec4(0.0) : texelFetch(previousPosition,coord,0);
 out_diffuse=aovSample==0 ? vec4(0.0) : texelFetch(previousDiffuse,coord,0);
 out_specular=aovSample==0 ? vec4(0.0) : texelFetch(previousSpecular,coord,0);
-rG(); aZ=uvec2(floor(gl_FragCoord.xy));
-gW(floor(gl_FragCoord.xy),int(aQ ^ dG(uint(sequenceSeed))));
-Ray ray=wJ((floor(gl_FragCoord.xy)+dZ(CAMERA_PIXEL))/viewportSize);
+rA(); aI=uvec2(floor(gl_FragCoord.xy));
+fi(floor(gl_FragCoord.xy),int(aC ^ cB(uint(sequenceSeed))));
+Ray ray=xj((floor(gl_FragCoord.xy)+sample2D(CAMERA_PIXEL))/viewportSize);
 if(cameraApertureRadius>0.0) {
 vec3 eye=viewInv[3].xyz, forward=-normalize(viewInv[2].xyz);
-vec3 jN=eye+ray.d*(cameraFocusDistance/dot(ray.d,forward));
-vec2 xU=cameraApertureRadius*kN(dZ(CAMERA_LENS));
-vec3 origin=(viewInv*vec4(xU,0.0,1.0)).xyz;
-gC(ray,origin,normalize(jN-origin));
+vec3 ju=eye+ray.d*(cameraFocusDistance/dot(ray.d,forward));
+vec2 uV=cameraApertureRadius*kq(sample2D(CAMERA_LENS));
+vec3 origin=(viewInv*vec4(uV,0.0,1.0)).xyz;
+initRay(ray,origin,normalize(ju-origin));
 }
 SurfaceInteraction si; OpenPbrResolvedMaterial material;
-if(!mQ(ray,si,material)) return;
+if(!mC(ray,si,material)) return;
 float depth=max(dot(si.position-viewInv[3].xyz,-normalize(viewInv[2].xyz)),0.0);
 bool unsupported=material.materialType!=STANDARD || material.alpha<1.0 || material.coatWeight>0.0 || material.fuzzWeight>0.0 || material.specularAnisotropy>0.0 || material.subsurfaceWeight>0.0 || dot(material.mediumScattering,material.mediumScattering)>0.0 || material.bevel>0.0 || material.thinFilmWeight>0.0 || material.transmissionDispersionScale>0.0;
 out_geometry+=vec4(depth,material.specularRoughness,material.baseMetalness,unsupported?1.0:0.0);
 #ifdef AOV_REFLECTANCE
 if((aovMask&24)!=0 && !unsupported){
-vec3 normal=si.normal,face=si.faceNormal;float cosine=-dot(normal,ray.d),fu=-dot(face,ray.d);
+vec3 normal=si.normal,face=si.faceNormal;float cosine=-dot(normal,ray.d),fa=-dot(face,ray.d);
 if(material.thinWalled&&(material.transmissionWeight>0.0||material.subsurfaceWeight>0.0)){
-if(cosine<0.0){normal=-normal;cosine=-cosine;}if(fu<0.0){face=-face;fu=-fu;}
+if(cosine<0.0){normal=-normal;cosine=-cosine;}if(fa<0.0){face=-face;fa=-fa;}
 }
-if(fu>=0.0&&cosine<0.0){normal=face;cosine=fu;}
+if(fa>=0.0&&cosine<0.0){normal=face;cosine=fa;}
 if(cosine<0.0)normal=-normal;
 material.inside=cosine<=0.0&&!material.thinWalled;material.relativeIor=cosine>0.0?1.0/material.specularIor:material.specularIor;
 vec3 tangent=vec3(0.0);
 #ifdef USE_TANGENT
 tangent=si.tangent;
 #endif
-ag prepared=ds(tangent,material,-ray.d,normal);
+OpenPbrPreparedBsdf prepared=bV(tangent,material,-ray.d,normal);
 vec3 diffuse,specular;aovReflectance(prepared,(aovMask&8)!=0,(aovMask&16)!=0,diffuse,specular);
 out_diffuse+=vec4(diffuse,1.0);out_specular+=vec4(specular,1.0);
 }

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

@@ -3,16 +3,16 @@
 // Copyright (c) 2019 HOVER, MIT license; see LICENSES/RayTracingRenderer-MIT.txt.
 // Shared primary-ray construction keeps beauty and catcher coverage aligned.
 export default /* glsl */ `
-Ray fj(vec2 coord) {
+Ray eT(vec2 coord) {
 vec4 origin = viewProjectionInv * vec4(coord * 2.0 - 1.0, -1.0, 1.0);
 vec4 target = viewProjectionInv * vec4(coord * 2.0 - 1.0, 1.0, 1.0);
 Ray cam;
-gC(cam, origin.xyz / origin.w, normalize(target.xyz * origin.w - origin.xyz * target.w));
+initRay(cam, origin.xyz / origin.w, normalize(target.xyz * origin.w - origin.xyz * target.w));
 
 return cam;
 }
 
-vec2 kN(vec2 p) {
+vec2 kq(vec2 p) {
 p = 2.0 * p - 1.0;
 bool greater = abs(p.x) > abs(p.y);
 float r = greater ? p.x : p.y;
@@ -20,16 +20,16 @@ float theta = greater ? 0.25 * PI * p.y / p.x : PI * (0.5 - 0.25 * p.x / p.y);
 return r * vec2(cos(theta), sin(theta));
 }
 
-Ray oG() {
-Ray cam = fj(v_Texcoord + tg());
+Ray ow() {
+Ray cam = eT(v_Texcoord + te());
 if (cameraApertureRadius > 0.0) {
 
 vec3 eye = viewInv[3].xyz;
 vec3 forward = -normalize(viewInv[2].xyz);
-vec3 jN = eye + cam.d * (cameraFocusDistance / dot(cam.d, forward));
-vec2 sK = cameraApertureRadius * kN(dZ(CAMERA_LENS));
-vec3 origin = vec3(viewInv * vec4(sK, 0.0, 1.0));
-gC(cam, origin, normalize(jN - origin));
+vec3 ju = eye + cam.d * (cameraFocusDistance / dot(cam.d, forward));
+vec2 sH = cameraApertureRadius * kq(sample2D(CAMERA_LENS));
+vec3 origin = vec3(viewInv * vec4(sH, 0.0, 1.0));
+initRay(cam, origin, normalize(ju - origin));
 }
 return cam;
 }

package/lib/renderer/glsl/chunks/emissiveLight.glsl.d.ts

@@ -0,0 +1,2 @@
+declare const _default: "\nvec4 emissiveTexel(int i) { return fetchData(emissiveLightData,i,emissiveLightColumns); }\nvoid sampleEmissiveSceneLight(vec3 x, vec3 rnd, inout LightSampleRec ls) {\n  ls.pdf=0.0;ls.isDelta=false;ls.sourceType=0;ls.sourceIndex=-1;\n  if(emissiveChoiceCount<=0) return;\n  float scaled=rnd.z*float(emissiveChoiceCount);int column=min(int(scaled),emissiveChoiceCount-1);\n  vec4 first=emissiveTexel(emissiveChoiceOffset+column);\n  int chosen=fract(scaled)<first.x?column:floatBitsToInt(first.y);\n  vec4 choice=emissiveTexel(emissiveChoiceOffset+chosen);int kind=floatBitsToInt(choice.z),index=floatBitsToInt(choice.w);\n  float pmf=emissiveTexel(emissivePmfOffset+chosen).x;\n  if(kind==1) {\n#ifdef USE_ENVMAP\n    sampleEnvLight(rnd.xy,ls);ls.pdf*=pmf;\n#endif\n  } else if(kind==2) {\n#ifdef NUM_LIGHTS\n    Light light=fetchLight(index);\n    \n    ls.sourceType=2;ls.sourceIndex=index;\n    sampleAnalyticLight(light,x,rnd.xy,ls);ls.pdf*=pmf;\n#endif\n  } else if(kind==3) {\n    vec4 a=emissiveTexel(index*6),b=emissiveTexel(index*6+1),c=emissiveTexel(index*6+2);\n    vec4 d=emissiveTexel(index*6+3),e=emissiveTexel(index*6+4),f=emissiveTexel(index*6+5);\n    float u=sqrt(rnd.x);vec3 bary=vec3(1.0-u,u*(1.0-rnd.y),u*rnd.y);\n    vec3 p=a.xyz*bary.x+b.xyz*bary.y+c.xyz*bary.z;\n    ls.normal=normalize(cross(b.xyz-a.xyz,c.xyz-a.xyz));ls.direction=normalize(p-x);ls.dist=distance(x,p);\n    float cosine=dot(ls.normal,-ls.direction);if(cosine<=0.0) return;\n    OpenPbrResolvedMaterial m;\n    getMaterial(m,int(c.w),bary.x*d.xy+bary.y*d.zw+bary.z*e.xy,bary.x*e.zw+bary.y*f.xy+bary.z*f.zw,0,vec3(1.0));\n    ls.emission=m.emission;ls.pdf=pmf*ls.dist*ls.dist/max(b.w*cosine,1e-8);ls.sourceType=1;ls.sourceIndex=index;\n  }\n}\n";
+export default _default;

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

@@ -0,0 +1,35 @@
+export default /* glsl */ `
+vec4 ei(int i) { return fetchData(emissiveLightData,i,emissiveLightColumns); }
+void jb(vec3 x, vec3 rnd, inout LightSampleRec ls) {
+ls.pdf=0.0;ls.isDelta=false;ls.sourceType=0;ls.sourceIndex=-1;
+if(emissiveChoiceCount<=0) return;
+float scaled=rnd.z*float(emissiveChoiceCount);int column=min(int(scaled),emissiveChoiceCount-1);
+vec4 first=ei(emissiveChoiceOffset+column);
+int chosen=fract(scaled)<first.x?column:floatBitsToInt(first.y);
+vec4 choice=ei(emissiveChoiceOffset+chosen);int kind=floatBitsToInt(choice.z),index=floatBitsToInt(choice.w);
+float pmf=ei(emissivePmfOffset+chosen).x;
+if(kind==1) {
+#ifdef USE_ENVMAP
+ig(rnd.xy,ls);ls.pdf*=pmf;
+#endif
+} else if(kind==2) {
+#ifdef NUM_LIGHTS
+Light light=fetchLight(index);
+
+ls.sourceType=2;ls.sourceIndex=index;
+sampleAnalyticLight(light,x,rnd.xy,ls);ls.pdf*=pmf;
+#endif
+} else if(kind==3) {
+vec4 a=ei(index*6),b=ei(index*6+1),c=ei(index*6+2);
+vec4 d=ei(index*6+3),e=ei(index*6+4),f=ei(index*6+5);
+float u=sqrt(rnd.x);vec3 bary=vec3(1.0-u,u*(1.0-rnd.y),u*rnd.y);
+vec3 p=a.xyz*bary.x+b.xyz*bary.y+c.xyz*bary.z;
+ls.normal=normalize(cross(b.xyz-a.xyz,c.xyz-a.xyz));ls.direction=normalize(p-x);ls.dist=distance(x,p);
+float cosine=dot(ls.normal,-ls.direction);if(cosine<=0.0) return;
+OpenPbrResolvedMaterial m;
+getMaterial(m,int(c.w),bary.x*d.xy+bary.y*d.zw+bary.z*e.xy,bary.x*e.zw+bary.y*f.xy+bary.z*f.zw,0,vec3(1.0));
+ls.emission=m.emission;ls.pdf=pmf*ls.dist*ls.dist/max(b.w*cosine,1e-8);ls.sourceType=1;ls.sourceIndex=index;
+}
+}
+`;
+//# sourceMappingURL=emissiveLight.glsl.js.map

package/lib/renderer/glsl/chunks/emissiveLight.glsl.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"emissiveLight.glsl.js","sourceRoot":"","sources":["../../../../src/renderer/glsl/chunks/emissiveLight.glsl.ts"],"names":[],"mappings":"AAAA,eAAe,UAAU,CAAC;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CAiCzB,CAAC","sourcesContent":["export default /* glsl */ `\nvec4 emissiveTexel(int i) { return fetchData(emissiveLightData,i,emissiveLightColumns); }\nvoid sampleEmissiveSceneLight(vec3 x, vec3 rnd, inout LightSampleRec ls) {\n  ls.pdf=0.0;ls.isDelta=false;ls.sourceType=0;ls.sourceIndex=-1;\n  if(emissiveChoiceCount<=0) return;\n  float scaled=rnd.z*float(emissiveChoiceCount);int column=min(int(scaled),emissiveChoiceCount-1);\n  vec4 first=emissiveTexel(emissiveChoiceOffset+column);\n  int chosen=fract(scaled)<first.x?column:floatBitsToInt(first.y);\n  vec4 choice=emissiveTexel(emissiveChoiceOffset+chosen);int kind=floatBitsToInt(choice.z),index=floatBitsToInt(choice.w);\n  float pmf=emissiveTexel(emissivePmfOffset+chosen).x;\n  if(kind==1) {\n#ifdef USE_ENVMAP\n    sampleEnvLight(rnd.xy,ls);ls.pdf*=pmf;\n#endif\n  } else if(kind==2) {\n#ifdef NUM_LIGHTS\n    Light light=fetchLight(index);\n    \n    ls.sourceType=2;ls.sourceIndex=index;\n    sampleAnalyticLight(light,x,rnd.xy,ls);ls.pdf*=pmf;\n#endif\n  } else if(kind==3) {\n    vec4 a=emissiveTexel(index*6),b=emissiveTexel(index*6+1),c=emissiveTexel(index*6+2);\n    vec4 d=emissiveTexel(index*6+3),e=emissiveTexel(index*6+4),f=emissiveTexel(index*6+5);\n    float u=sqrt(rnd.x);vec3 bary=vec3(1.0-u,u*(1.0-rnd.y),u*rnd.y);\n    vec3 p=a.xyz*bary.x+b.xyz*bary.y+c.xyz*bary.z;\n    ls.normal=normalize(cross(b.xyz-a.xyz,c.xyz-a.xyz));ls.direction=normalize(p-x);ls.dist=distance(x,p);\n    float cosine=dot(ls.normal,-ls.direction);if(cosine<=0.0) return;\n    OpenPbrResolvedMaterial m;\n    getMaterial(m,int(c.w),bary.x*d.xy+bary.y*d.zw+bary.z*e.xy,bary.x*e.zw+bary.y*f.xy+bary.z*f.zw,0,vec3(1.0));\n    ls.emission=m.emission;ls.pdf=pmf*ls.dist*ls.dist/max(b.w*cosine,1e-8);ls.sourceType=1;ls.sourceIndex=index;\n  }\n}\n`;\n"]}