@plasius/gpu-renderer 0.2.5 → 0.2.7

package/dist/index.cjs

@@ -24,9 +24,11 @@ __export(index_exports, {
   createGpuRenderer: () => createGpuRenderer,
   createRayTracingRenderPlan: () => createRayTracingRenderPlan,
   createRendererDebugHooks: () => createRendererDebugHooks,
+  createWavefrontAdaptiveSamplingLevels: () => createWavefrontAdaptiveSamplingLevels,
   createWavefrontBvhBuildLevels: () => createWavefrontBvhBuildLevels,
   createWavefrontBvhSortStages: () => createWavefrontBvhSortStages,
   createWavefrontEmissiveTriangleIndexSource: () => createWavefrontEmissiveTriangleIndexSource,
+  createWavefrontGpuMaterialSource: () => createWavefrontGpuMaterialSource,
   createWavefrontGpuMeshSource: () => createWavefrontGpuMeshSource,
   createWavefrontMeshAcceleration: () => createWavefrontMeshAcceleration,
   createWavefrontPathTracingComputeConfig: () => createWavefrontPathTracingComputeConfig,
@@ -70,36 +72,179 @@ __export(index_exports, {
 });
 module.exports = __toCommonJS(index_exports);
 
+// src/wavefront-frame-runtime.js
+function createGpuParallelismCounters() {
+  return {
+    directDispatches: 0,
+    directWorkgroups: 0,
+    directShaderInvocations: 0,
+    multiWorkgroupDispatches: 0,
+    largestDirectWorkgroupsPerDispatch: 0,
+    indirectDispatches: 0,
+    estimatedIndirectWorkgroupsUpperBound: 0,
+    estimatedIndirectShaderInvocationsUpperBound: 0,
+    indirectDispatchesWithMultiWorkgroupCapacity: 0,
+    largestEstimatedIndirectWorkgroupsPerDispatch: 0
+  };
+}
+function countDispatchWorkgroups(groups) {
+  return groups.reduce((product, value) => {
+    const numeric = Number(value ?? 1);
+    const count = Number.isFinite(numeric) ? Math.max(1, Math.trunc(numeric)) : 1;
+    return product * count;
+  }, 1);
+}
+function recordDirectDispatch(parallelism, groups, invocationsPerWorkgroup = 1) {
+  const workgroups = countDispatchWorkgroups(groups);
+  parallelism.directDispatches += 1;
+  parallelism.directWorkgroups += workgroups;
+  parallelism.directShaderInvocations += workgroups * invocationsPerWorkgroup;
+  parallelism.largestDirectWorkgroupsPerDispatch = Math.max(
+    parallelism.largestDirectWorkgroupsPerDispatch,
+    workgroups
+  );
+  if (workgroups > 1) {
+    parallelism.multiWorkgroupDispatches += 1;
+  }
+}
+function recordIndirectDispatch(parallelism, estimatedWorkgroupsUpperBound, invocationsPerWorkgroup = 1) {
+  const workgroups = Math.max(1, Math.trunc(Number(estimatedWorkgroupsUpperBound) || 1));
+  parallelism.indirectDispatches += 1;
+  parallelism.estimatedIndirectWorkgroupsUpperBound += workgroups;
+  parallelism.estimatedIndirectShaderInvocationsUpperBound += workgroups * invocationsPerWorkgroup;
+  parallelism.largestEstimatedIndirectWorkgroupsPerDispatch = Math.max(
+    parallelism.largestEstimatedIndirectWorkgroupsPerDispatch,
+    workgroups
+  );
+  if (workgroups > 1) {
+    parallelism.indirectDispatchesWithMultiWorkgroupCapacity += 1;
+  }
+}
+function createGpuParallelismDiagnostics(adapterDiagnostics, counters) {
+  const totalEstimatedWorkgroupsUpperBound = counters.directWorkgroups + counters.estimatedIndirectWorkgroupsUpperBound;
+  const totalEstimatedShaderInvocationsUpperBound = counters.directShaderInvocations + counters.estimatedIndirectShaderInvocationsUpperBound;
+  const exposesMultiWorkgroupParallelism = counters.multiWorkgroupDispatches > 0 || counters.indirectDispatchesWithMultiWorkgroupCapacity > 0;
+  return Object.freeze({
+    ...adapterDiagnostics,
+    directDispatches: counters.directDispatches,
+    directWorkgroups: counters.directWorkgroups,
+    directShaderInvocations: counters.directShaderInvocations,
+    multiWorkgroupDispatches: counters.multiWorkgroupDispatches,
+    largestDirectWorkgroupsPerDispatch: counters.largestDirectWorkgroupsPerDispatch,
+    indirectDispatches: counters.indirectDispatches,
+    estimatedIndirectWorkgroupsUpperBound: counters.estimatedIndirectWorkgroupsUpperBound,
+    estimatedIndirectShaderInvocationsUpperBound: counters.estimatedIndirectShaderInvocationsUpperBound,
+    indirectDispatchesWithMultiWorkgroupCapacity: counters.indirectDispatchesWithMultiWorkgroupCapacity,
+    largestEstimatedIndirectWorkgroupsPerDispatch: counters.largestEstimatedIndirectWorkgroupsPerDispatch,
+    totalEstimatedWorkgroupsUpperBound,
+    totalEstimatedShaderInvocationsUpperBound,
+    exposesMultiWorkgroupParallelism,
+    likelyUsesMoreThanOnePhysicalGpuCore: null,
+    coreUtilizationStatus: "not-exposed-by-webgpu"
+  });
+}
+function createGpuWorkerJobDiagnostics(parallelism, commandSubmissions, frameTimeMs, awaitedGpuCompletion) {
+  const directDispatchesCompleted = Math.max(0, Number(parallelism?.directDispatches ?? 0));
+  const indirectDispatchesCompleted = Math.max(
+    0,
+    Number(parallelism?.indirectDispatches ?? 0)
+  );
+  const completedPerFrame = directDispatchesCompleted + indirectDispatchesCompleted;
+  const completedPerSubmission = commandSubmissions > 0 ? completedPerFrame / commandSubmissions : completedPerFrame;
+  const completedPerSecond = awaitedGpuCompletion && frameTimeMs > 0 ? completedPerFrame * 1e3 / frameTimeMs : null;
+  return Object.freeze({
+    completedPerFrame,
+    completedPerSecond,
+    completedPerSubmission,
+    directDispatchesCompleted,
+    indirectDispatchesCompleted,
+    frameTimeMs,
+    awaitedGpuCompletion
+  });
+}
+function createGpuSubmissionBatcher({
+  device,
+  frameIndex,
+  maxFramePassesPerSubmission,
+  startingSubmissionCount = 0,
+  labelPrefix = "plasius.wavefront.frame"
+}) {
+  let encodedFramePasses = 0;
+  let submissionCount = 0;
+  let encoder = createCommandEncoder();
+  function createCommandEncoder() {
+    return device.createCommandEncoder({
+      label: `${labelPrefix}.${frameIndex}.batched.${startingSubmissionCount + submissionCount + 1}`
+    });
+  }
+  function submitCurrentEncoder() {
+    if (encodedFramePasses <= 0) {
+      return false;
+    }
+    device.queue.submit([encoder.finish()]);
+    submissionCount += 1;
+    encodedFramePasses = 0;
+    encoder = createCommandEncoder();
+    return true;
+  }
+  return Object.freeze({
+    reserve(passCount = 1) {
+      if (encodedFramePasses > 0 && encodedFramePasses + passCount > maxFramePassesPerSubmission) {
+        submitCurrentEncoder();
+      }
+      encodedFramePasses += passCount;
+      return encoder;
+    },
+    flush() {
+      submitCurrentEncoder();
+      return submissionCount;
+    },
+    getSubmissionCount() {
+      return submissionCount;
+    }
+  });
+}
+
 // src/wavefront-compute.js
 var DEFAULT_WIDTH = 1280;
 var DEFAULT_HEIGHT = 720;
 var DEFAULT_MAX_DEPTH = 6;
 var DEFAULT_TILE_SIZE = 128;
 var DEFAULT_SAMPLES_PER_PIXEL = 1;
+var MAX_SAMPLES_PER_PIXEL = 256;
+var DEFAULT_BRDF_LUT_SIZE = 128;
+var DEFAULT_BRDF_LUT_SAMPLE_COUNT = 256;
 var DEFAULT_MAX_FRAME_PASSES_PER_SUBMISSION = 256;
 var DEFAULT_SCENE_OBJECT_CAPACITY = 128;
 var DEFAULT_ENVIRONMENT_PORTAL_CAPACITY = 32;
+var DEFAULT_MEDIUM_PHASE_MODEL = 0;
 var WORKGROUP_SIZE = 64;
 var rendererWavefrontComputeMode = "webgpu-compute";
 var rendererWavefrontComputeWorkgroupSize = WORKGROUP_SIZE;
 var rendererWavefrontComputeStatsStride = 8;
 var RAY_RECORD_BYTES = 80;
-var HIT_RECORD_BYTES = 208;
-var SCENE_OBJECT_RECORD_BYTES = 96;
+var HIT_RECORD_BYTES = 256;
+var SCENE_OBJECT_RECORD_BYTES = 160;
 var MESH_VERTEX_RECORD_BYTES = 48;
-var MESH_RANGE_RECORD_BYTES = 96;
-var TRIANGLE_RECORD_BYTES = 208;
+var MESH_RANGE_RECORD_BYTES = 240;
+var TRIANGLE_RECORD_BYTES = 352;
+var GPU_MATERIAL_RECORD_BYTES = 192;
 var BVH_NODE_RECORD_BYTES = 48;
 var BVH_LEAF_REF_RECORD_BYTES = 16;
 var EMISSIVE_TRIANGLE_INDEX_BYTES = 4;
 var ENVIRONMENT_PORTAL_RECORD_BYTES = 96;
+var MEDIUM_TABLE_ROWS = 2;
 var ACCUMULATION_RECORD_BYTES = 16;
 var PATH_VERTEX_RECORD_BYTES = 16;
-var CONFIG_BUFFER_BYTES = 304;
+var GPU_SUBMITTED_WORK_TIMEOUT_MS = 5e3;
+var GPU_READBACK_COMPLETION_TIMEOUT_MS = 6e4;
+var GPU_MAX_SUBMITTED_WORK_TIMEOUT_MS = 6e4;
+var CONFIG_BUFFER_BYTES = 320;
 var COUNTER_DISPATCH_ARGS_OFFSET = 16;
 var INDIRECT_DISPATCH_ARGS_BYTES = 12;
 var COUNTER_BUFFER_BYTES = 32;
 var TRACE_STORAGE_BUFFER_BINDINGS = 10;
+var BRDF_LUT_UPLOAD_CACHE = /* @__PURE__ */ new Map();
 var MATERIAL_DIFFUSE = 0;
 var MATERIAL_METAL = 1;
 var MATERIAL_DIELECTRIC = 2;
@@ -134,6 +279,7 @@ var wavefrontPathTracingComputeLimits = Object.freeze({
   meshVertexRecordBytes: MESH_VERTEX_RECORD_BYTES,
   meshRangeRecordBytes: MESH_RANGE_RECORD_BYTES,
   triangleRecordBytes: TRIANGLE_RECORD_BYTES,
+  materialRecordBytes: GPU_MATERIAL_RECORD_BYTES,
   bvhNodeRecordBytes: BVH_NODE_RECORD_BYTES,
   bvhLeafReferenceRecordBytes: BVH_LEAF_REF_RECORD_BYTES,
   emissiveTriangleIndexBytes: EMISSIVE_TRIANGLE_INDEX_BYTES,
@@ -217,6 +363,32 @@ function asColor(value, fallback = [1, 1, 1, 1]) {
     clamp(readFiniteNumber("color[3]", value[3], fallback[3] ?? 1), 0, 1)
   ];
 }
+function deriveLegacySheenColor(baseColor, sheen, sheenTint) {
+  const sheenStrength = clamp(Number(sheen) || 0, 0, 1);
+  if (sheenStrength <= 0) {
+    return [0, 0, 0, 1];
+  }
+  const tint = clamp(Number(sheenTint) || 0, 0, 1);
+  const base = asColor(baseColor, [1, 1, 1, 1]);
+  return [
+    clamp((1 - tint) * sheenStrength + base[0] * tint * sheenStrength, 0, 1),
+    clamp((1 - tint) * sheenStrength + base[1] * tint * sheenStrength, 0, 1),
+    clamp((1 - tint) * sheenStrength + base[2] * tint * sheenStrength, 0, 1),
+    1
+  ];
+}
+function resolveSheenColor(input, fallbackBaseColor) {
+  if (input?.sheenColor || input?.material?.sheenColor) {
+    return asColor(input.sheenColor ?? input.material?.sheenColor, [0, 0, 0, 1]).map(
+      (value, index) => index < 3 ? clamp(value, 0, 1) : 1
+    );
+  }
+  return deriveLegacySheenColor(
+    fallbackBaseColor,
+    input?.sheen ?? input?.material?.sheen,
+    input?.sheenTint ?? input?.material?.sheenTint
+  );
+}
 function resolveEnvironmentMap(input = null) {
   const source = input && typeof input === "object" ? input : null;
   const hasTexture = Boolean(source?.view || source?.texture || source?.data);
@@ -226,6 +398,11 @@ function resolveEnvironmentMap(input = null) {
     enabled: hasTexture && source?.enabled !== false,
     width,
     height,
+    mipLevelCount: readPositiveInteger(
+      "environmentMap.mipLevelCount",
+      source?.mipLevelCount,
+      1
+    ),
     format: typeof source?.format === "string" ? source.format : "rgba16float",
     projection: typeof source?.projection === "string" ? source.projection : "equirectangular",
     texture: source?.texture ?? null,
@@ -237,7 +414,8 @@ function resolveEnvironmentMap(input = null) {
     ambientStrength: Math.max(
       0,
       readFiniteNumber("environmentMap.ambientStrength", source?.ambientStrength, 0.32)
-    )
+    ),
+    hasImportanceData: source?.hasImportanceData === true
   });
 }
 function resolveDeferredPathResolve(options = {}) {
@@ -402,6 +580,156 @@ function deriveBounds(input) {
   }
   return null;
 }
+function deriveBeerLambertAbsorptionFromAttenuationColor(attenuationColor, attenuationDistance, density = 1) {
+  const distance = Number(attenuationDistance);
+  const densityScale = Math.max(0, Number(density) || 0);
+  if (!Number.isFinite(distance) || distance <= 0 || densityScale <= 0) {
+    return [0, 0, 0];
+  }
+  return attenuationColor.slice(0, 3).map((channel) => {
+    const clamped = clamp(Number(channel) || 0, 1e-4, 1);
+    return Math.max(0, -Math.log(clamped) / distance * densityScale);
+  });
+}
+function readMediumPhaseModel(value) {
+  if (typeof value === "number" && Number.isFinite(value)) {
+    return Math.max(0, Math.trunc(value));
+  }
+  switch (String(value ?? "").trim().toLowerCase()) {
+    case "isotropic":
+    default:
+      return DEFAULT_MEDIUM_PHASE_MODEL;
+  }
+}
+function resolveWavefrontVolumeInput(input) {
+  return input?.volume ?? input?.material?.volume ?? null;
+}
+function normalizeWavefrontThickness(input, label) {
+  const volume = resolveWavefrontVolumeInput(input);
+  return Math.max(
+    0,
+    readFiniteNumber(
+      label,
+      input?.thickness ?? volume?.thickness ?? input?.material?.thickness,
+      0
+    )
+  );
+}
+function resolveWavefrontMediumId(input, fallbackId = 1) {
+  return input?.mediumRefId ?? input?.mediumId ?? input?.material?.mediumId ?? input?.materialRefId ?? input?.material?.id ?? input?.materialId ?? input?.id ?? fallbackId;
+}
+function deriveWavefrontTransportMedium(input, fallbackId = 1) {
+  const resolvedId = resolveWavefrontMediumId(input, fallbackId);
+  if (input?.medium) {
+    return normalizeWavefrontMedium(
+      {
+        ...input.medium,
+        id: input.medium.id ?? input.medium.mediumId ?? resolvedId
+      },
+      fallbackId
+    );
+  }
+  const volume = resolveWavefrontVolumeInput(input);
+  if (!volume) {
+    return null;
+  }
+  return normalizeWavefrontMedium(
+    {
+      id: resolvedId,
+      phaseModel: volume.phaseModel,
+      density: volume.density,
+      attenuationColor: volume.attenuationColor,
+      attenuationDistance: volume.attenuationDistance,
+      absorption: volume.absorption,
+      scattering: volume.scattering
+    },
+    fallbackId
+  );
+}
+function normalizeWavefrontMedium(input = {}, index = 0) {
+  const id = readNonNegativeInteger("medium id", input.id ?? input.mediumId, index);
+  const density = Math.max(0, readFiniteNumber("medium density", input.density, 1));
+  const attenuationColor = asColor(
+    input.attenuationColor ?? input.color ?? input.medium?.attenuationColor,
+    [1, 1, 1, 1]
+  );
+  const attenuationDistance = readFiniteNumber(
+    "medium attenuationDistance",
+    input.attenuationDistance ?? input.distance ?? input.medium?.attenuationDistance,
+    0
+  );
+  const absorption = Array.isArray(input.absorption) || Array.isArray(input.medium?.absorption) ? asVec3(input.absorption ?? input.medium?.absorption, [0, 0, 0]).map(
+    (value) => Math.max(0, Number(value) || 0)
+  ) : deriveBeerLambertAbsorptionFromAttenuationColor(
+    attenuationColor,
+    attenuationDistance,
+    density
+  );
+  const scattering = asVec3(
+    input.scattering ?? input.medium?.scattering,
+    [0, 0, 0]
+  ).map((value) => Math.max(0, Number(value) || 0));
+  return Object.freeze({
+    id,
+    phaseModel: readMediumPhaseModel(input.phaseModel ?? input.medium?.phaseModel),
+    density,
+    attenuationColor: Object.freeze(attenuationColor),
+    attenuationDistance,
+    absorption: Object.freeze(absorption),
+    scattering: Object.freeze(scattering)
+  });
+}
+function collectWavefrontMediums(options, meshes, sceneObjects = []) {
+  const mediumsById = /* @__PURE__ */ new Map();
+  mediumsById.set(
+    0,
+    Object.freeze({
+      id: 0,
+      phaseModel: DEFAULT_MEDIUM_PHASE_MODEL,
+      density: 0,
+      attenuationColor: Object.freeze([1, 1, 1, 1]),
+      attenuationDistance: 0,
+      absorption: Object.freeze([0, 0, 0]),
+      scattering: Object.freeze([0, 0, 0])
+    })
+  );
+  const register = (input, fallbackId = mediumsById.size) => {
+    if (!input) {
+      return;
+    }
+    const normalized = normalizeWavefrontMedium(
+      typeof input === "object" ? { id: fallbackId, ...input } : { id: fallbackId },
+      fallbackId
+    );
+    const existing = mediumsById.get(normalized.id);
+    if (existing && JSON.stringify(existing) !== JSON.stringify(normalized)) {
+      throw new Error(`Medium id ${normalized.id} is defined more than once with different values.`);
+    }
+    mediumsById.set(normalized.id, normalized);
+  };
+  for (const medium of options.mediums ?? []) {
+    register(medium);
+  }
+  for (const mesh of meshes) {
+    register(mesh.medium, mesh.mediumRefId ?? mesh.medium?.id ?? 0);
+  }
+  for (const mesh of meshes) {
+    if ((mesh.mediumRefId ?? 0) > 0 && !mediumsById.has(mesh.mediumRefId)) {
+      register({ id: mesh.mediumRefId });
+    }
+  }
+  for (const object of sceneObjects) {
+    register(object.medium, object.mediumRefId ?? object.medium?.id ?? 0);
+  }
+  for (const object of sceneObjects) {
+    if ((object.mediumRefId ?? 0) > 0 && !mediumsById.has(object.mediumRefId)) {
+      register({ id: object.mediumRefId });
+    }
+  }
+  return Object.freeze(
+    Array.from(mediumsById.values()).sort((left, right) => left.id - right.id)
+  );
+}
 function normalizeWavefrontSceneObject(input = {}, index = 0) {
   const bounds = deriveBounds(input);
   const kind = readObjectKind(input.kind ?? input.type ?? (bounds ? "box" : "sphere"));
@@ -411,7 +739,8 @@ function normalizeWavefrontSceneObject(input = {}, index = 0) {
     input.halfExtent ?? input.halfExtents ?? input.extents ?? bounds?.halfExtent,
     [0.5, 0.5, 0.5]
   ).map((value) => Math.max(value, 1e-3));
-  const materialKind = readMaterialKind(input.materialKind ?? input.material?.kind);
+  const materialKindInput = input.materialKind ?? input.material?.kind;
+  const materialKind = readMaterialKind(materialKindInput);
   const color = asColor(
     input.color ?? input.baseColor ?? input.albedo ?? input.material?.color ?? input.material?.baseColor,
     [0.72, 0.72, 0.68, 1]
@@ -420,19 +749,55 @@ function normalizeWavefrontSceneObject(input = {}, index = 0) {
     input.emission ?? input.emissive ?? input.material?.emission ?? input.material?.emissive,
     [0, 0, 0, 1]
   );
+  const opacity = clamp(readFiniteNumber("opacity", input.opacity ?? input.material?.opacity, color[3] ?? 1), 0, 1);
+  const transmission = clamp(
+    readFiniteNumber("transmission", input.transmission ?? input.material?.transmission, 0),
+    0,
+    1
+  );
+  const sheenColor = resolveSheenColor(input, color);
+  const specularColor = asColor(
+    input.specularColor ?? input.material?.specularColor,
+    [1, 1, 1, 1]
+  ).map((value, componentIndex) => componentIndex < 3 ? clamp(value, 0, 1) : 1);
+  const medium = deriveWavefrontTransportMedium(input, index + 1);
+  const resolvedMaterialKind = emission[0] > 0 || emission[1] > 0 || emission[2] > 0 ? MATERIAL_EMISSIVE : materialKindInput === void 0 || materialKindInput === null ? transmission > 1e-3 || opacity < 0.999 ? MATERIAL_TRANSPARENT : materialKind : materialKind;
   return Object.freeze({
     id: readNonNegativeInteger("id", input.id, index + 1),
     kind,
-    materialKind: emission[0] > 0 || emission[1] > 0 || emission[2] > 0 ? MATERIAL_EMISSIVE : materialKind,
+    materialKind: resolvedMaterialKind,
     flags: readNonNegativeInteger("flags", input.flags, 0),
+    mediumRefId: readNonNegativeInteger(
+      "mediumRefId",
+      input.mediumRefId ?? medium?.id ?? input.medium?.id ?? input.mediumId,
+      0
+    ),
+    medium,
     center: Object.freeze(center),
     halfExtent: Object.freeze(halfExtent),
     color: Object.freeze(color),
     emission: Object.freeze(emission),
     roughness: clamp(readFiniteNumber("roughness", input.roughness ?? input.material?.roughness, 0.72), 0, 1),
     metallic: clamp(readFiniteNumber("metallic", input.metallic ?? input.material?.metallic, 0), 0, 1),
-    opacity: clamp(readFiniteNumber("opacity", input.opacity ?? input.material?.opacity, color[3] ?? 1), 0, 1),
-    ior: clamp(readFiniteNumber("ior", input.ior ?? input.material?.ior, 1.45), 1, 3)
+    opacity,
+    ior: clamp(readFiniteNumber("ior", input.ior ?? input.material?.ior, 1.45), 1, 3),
+    sheen: clamp(readFiniteNumber("sheen", input.sheen ?? input.material?.sheen, 0), 0, 1),
+    sheenTint: clamp(readFiniteNumber("sheenTint", input.sheenTint ?? input.material?.sheenTint, 0), 0, 1),
+    sheenColor: Object.freeze(sheenColor),
+    clearcoat: clamp(readFiniteNumber("clearcoat", input.clearcoat ?? input.material?.clearcoat, 0), 0, 1),
+    clearcoatRoughness: clamp(
+      readFiniteNumber(
+        "clearcoatRoughness",
+        input.clearcoatRoughness ?? input.material?.clearcoatRoughness,
+        0.08
+      ),
+      0,
+      1
+    ),
+    specular: clamp(readFiniteNumber("specular", input.specular ?? input.material?.specular, 1), 0, 1),
+    specularColor: Object.freeze(specularColor),
+    thickness: normalizeWavefrontThickness(input, "thickness"),
+    transmission
   });
 }
 function createDefaultWavefrontSceneObjects() {
@@ -504,7 +869,8 @@ function normalizeWavefrontMesh(input = {}, meshIndex = 0) {
     input.uvs ?? input.texcoords ?? input.uv,
     (value) => readFiniteNumber("mesh uv", value, 0)
   ) : null;
-  const materialKind = readMaterialKind(input.materialKind ?? input.material?.kind);
+  const materialKindInput = input.materialKind ?? input.material?.kind;
+  const materialKind = readMaterialKind(materialKindInput);
   const color = asColor(
     input.color ?? input.baseColor ?? input.albedo ?? input.material?.color ?? input.material?.baseColor,
     [0.72, 0.72, 0.68, 1]
@@ -513,13 +879,26 @@ function normalizeWavefrontMesh(input = {}, meshIndex = 0) {
     input.emission ?? input.emissive ?? input.material?.emission ?? input.material?.emissive,
     [0, 0, 0, 1]
   );
+  const opacity = clamp(readFiniteNumber("opacity", input.opacity ?? input.material?.opacity, color[3] ?? 1), 0, 1);
+  const transmission = clamp(
+    readFiniteNumber("transmission", input.transmission ?? input.material?.transmission, 0),
+    0,
+    1
+  );
+  const sheenColor = resolveSheenColor(input, color);
+  const specularColor = asColor(
+    input.specularColor ?? input.material?.specularColor,
+    [1, 1, 1, 1]
+  ).map((value, componentIndex) => componentIndex < 3 ? clamp(value, 0, 1) : 1);
+  const medium = deriveWavefrontTransportMedium(input, meshIndex + 1);
+  const resolvedMaterialKind = emission[0] > 0 || emission[1] > 0 || emission[2] > 0 ? MATERIAL_EMISSIVE : materialKindInput === void 0 || materialKindInput === null ? transmission > 1e-3 || opacity < 0.999 ? MATERIAL_TRANSPARENT : materialKind : materialKind;
   return Object.freeze({
     id: readNonNegativeInteger("mesh id", input.id, meshIndex + 1),
     positions: Object.freeze(Array.from(positions, (value) => readFiniteNumber("mesh position", value, 0))),
     indices: Object.freeze(indices),
     normals: normals ? Object.freeze(normals) : null,
     uvs: uvs ? Object.freeze(uvs) : null,
-    materialKind: emission[0] > 0 || emission[1] > 0 || emission[2] > 0 ? MATERIAL_EMISSIVE : materialKind,
+    materialKind: resolvedMaterialKind,
     flags: readNonNegativeInteger("mesh flags", input.flags, 0),
     materialRefId: readNonNegativeInteger(
       "mesh materialRefId",
@@ -528,17 +907,176 @@ function normalizeWavefrontMesh(input = {}, meshIndex = 0) {
     ),
     mediumRefId: readNonNegativeInteger(
       "mesh mediumRefId",
-      input.mediumRefId ?? input.medium?.id ?? input.mediumId,
+      input.mediumRefId ?? medium?.id ?? input.medium?.id ?? input.mediumId ?? input.material?.mediumId,
       0
     ),
+    medium,
     color: Object.freeze(color),
     emission: Object.freeze(emission),
     roughness: clamp(readFiniteNumber("roughness", input.roughness ?? input.material?.roughness, 0.72), 0, 1),
     metallic: clamp(readFiniteNumber("metallic", input.metallic ?? input.material?.metallic, 0), 0, 1),
-    opacity: clamp(readFiniteNumber("opacity", input.opacity ?? input.material?.opacity, color[3] ?? 1), 0, 1),
-    ior: clamp(readFiniteNumber("ior", input.ior ?? input.material?.ior, 1.45), 1, 3)
+    opacity,
+    ior: clamp(readFiniteNumber("ior", input.ior ?? input.material?.ior, 1.45), 1, 3),
+    sheen: clamp(readFiniteNumber("sheen", input.sheen ?? input.material?.sheen, 0), 0, 1),
+    sheenTint: clamp(readFiniteNumber("sheenTint", input.sheenTint ?? input.material?.sheenTint, 0), 0, 1),
+    sheenColor: Object.freeze(sheenColor),
+    clearcoat: clamp(readFiniteNumber("clearcoat", input.clearcoat ?? input.material?.clearcoat, 0), 0, 1),
+    clearcoatRoughness: clamp(
+      readFiniteNumber(
+        "clearcoatRoughness",
+        input.clearcoatRoughness ?? input.material?.clearcoatRoughness,
+        0.08
+      ),
+      0,
+      1
+    ),
+    specular: clamp(readFiniteNumber("specular", input.specular ?? input.material?.specular, 1), 0, 1),
+    specularColor: Object.freeze(specularColor),
+    thickness: normalizeWavefrontThickness(input, "mesh thickness"),
+    transmission,
+    baseColorTexture: input.baseColorTexture ?? input.material?.baseColorTexture ?? null,
+    metallicRoughnessTexture: input.metallicRoughnessTexture ?? input.material?.metallicRoughnessTexture ?? null,
+    normalTexture: input.normalTexture ?? input.material?.normalTexture ?? null,
+    occlusionTexture: input.occlusionTexture ?? input.material?.occlusionTexture ?? null,
+    emissiveTexture: input.emissiveTexture ?? input.material?.emissiveTexture ?? null
   });
 }
+function clampUnit(value) {
+  return clamp(Number(value) || 0, 0, 1);
+}
+function srgbToLinear(value) {
+  const channel = clampUnit(value);
+  if (channel <= 0.04045) {
+    return channel / 12.92;
+  }
+  return ((channel + 0.055) / 1.055) ** 2.4;
+}
+function sampleTextureRgba(texture, uv = [0, 0], colorSpace = "linear") {
+  if (!texture || !Number.isFinite(texture.width) || !Number.isFinite(texture.height) || !texture.data || texture.width <= 0 || texture.height <= 0) {
+    return [1, 1, 1, 1];
+  }
+  const u = (uv[0] % 1 + 1) % 1;
+  const v = (uv[1] % 1 + 1) % 1;
+  const x = Math.min(texture.width - 1, Math.max(0, Math.round(u * (texture.width - 1))));
+  const y = Math.min(texture.height - 1, Math.max(0, Math.round((1 - v) * (texture.height - 1))));
+  const offset = (y * texture.width + x) * 4;
+  const data = texture.data;
+  const scale2 = resolveTextureSampleScale(data);
+  const defaultChannel = scale2 === 1 ? 1 : Math.round(1 / scale2);
+  const color = [
+    (data[offset] ?? defaultChannel) * scale2,
+    (data[offset + 1] ?? defaultChannel) * scale2,
+    (data[offset + 2] ?? defaultChannel) * scale2,
+    (data[offset + 3] ?? defaultChannel) * scale2
+  ];
+  if (colorSpace === "srgb") {
+    return [srgbToLinear(color[0]), srgbToLinear(color[1]), srgbToLinear(color[2]), color[3]];
+  }
+  return color;
+}
+function resolveTextureSampleScale(data) {
+  if (data instanceof Uint8Array || data instanceof Uint8ClampedArray) {
+    return 1 / 255;
+  }
+  if (data instanceof Uint16Array) {
+    return 1 / 65535;
+  }
+  if (Array.isArray(data) && data.some((value) => Number(value) > 1)) {
+    return 1 / 255;
+  }
+  return 1;
+}
+function normalizeVectorOrFallback(vector, fallback) {
+  return normalize(Array.isArray(vector) ? vector : fallback, fallback);
+}
+function buildTriangleTangentBasis(v0, v1, v2, uv0, uv1, uv2, fallbackNormal) {
+  const edge1 = subtract(v1, v0);
+  const edge2 = subtract(v2, v0);
+  const deltaUv1 = [uv1[0] - uv0[0], uv1[1] - uv0[1]];
+  const deltaUv2 = [uv2[0] - uv0[0], uv2[1] - uv0[1]];
+  const determinant = deltaUv1[0] * deltaUv2[1] - deltaUv1[1] * deltaUv2[0];
+  if (Math.abs(determinant) < 1e-6) {
+    const tangentFallback = Math.abs(fallbackNormal[1]) < 0.999 ? [0, 1, 0] : [1, 0, 0];
+    const tangent2 = normalize(cross(tangentFallback, fallbackNormal), [1, 0, 0]);
+    const bitangent2 = normalize(cross(fallbackNormal, tangent2), [0, 0, 1]);
+    return { tangent: tangent2, bitangent: bitangent2 };
+  }
+  const inverse = 1 / determinant;
+  const tangent = normalize(
+    [
+      inverse * (edge1[0] * deltaUv2[1] - edge2[0] * deltaUv1[1]),
+      inverse * (edge1[1] * deltaUv2[1] - edge2[1] * deltaUv1[1]),
+      inverse * (edge1[2] * deltaUv2[1] - edge2[2] * deltaUv1[1])
+    ],
+    [1, 0, 0]
+  );
+  const bitangent = normalize(
+    [
+      inverse * (-edge1[0] * deltaUv2[0] + edge2[0] * deltaUv1[0]),
+      inverse * (-edge1[1] * deltaUv2[0] + edge2[1] * deltaUv1[0]),
+      inverse * (-edge1[2] * deltaUv2[0] + edge2[2] * deltaUv1[0])
+    ],
+    [0, 0, 1]
+  );
+  return { tangent, bitangent };
+}
+function applyNormalMap(normal, tangent, bitangent, normalTexture, uv) {
+  if (!normalTexture) {
+    return normalizeVectorOrFallback(normal, [0, 1, 0]);
+  }
+  const sample = sampleTextureRgba(normalTexture, uv, "linear");
+  const strength = clampUnit(normalTexture.scale ?? 1);
+  const tangentNormal = normalize(
+    [
+      (sample[0] * 2 - 1) * strength,
+      (sample[1] * 2 - 1) * strength,
+      1 + (sample[2] * 2 - 1 - 1) * strength
+    ],
+    [0, 0, 1]
+  );
+  return normalize(
+    [
+      tangent[0] * tangentNormal[0] + bitangent[0] * tangentNormal[1] + normal[0] * tangentNormal[2],
+      tangent[1] * tangentNormal[0] + bitangent[1] * tangentNormal[1] + normal[1] * tangentNormal[2],
+      tangent[2] * tangentNormal[0] + bitangent[2] * tangentNormal[1] + normal[2] * tangentNormal[2]
+    ],
+    normal
+  );
+}
+function sampleBaseColor(mesh, uv) {
+  const sample = mesh.baseColorTexture ? sampleTextureRgba(mesh.baseColorTexture, uv, "srgb") : [1, 1, 1, 1];
+  return [
+    clampUnit(mesh.color[0] * sample[0]),
+    clampUnit(mesh.color[1] * sample[1]),
+    clampUnit(mesh.color[2] * sample[2]),
+    clampUnit((mesh.color[3] ?? 1) * sample[3])
+  ];
+}
+function sampleSurfaceMaterial(mesh, uv) {
+  const textureSample = mesh.metallicRoughnessTexture ? sampleTextureRgba(mesh.metallicRoughnessTexture, uv, "linear") : [1, 1, 1, 1];
+  return {
+    roughness: clamp(mesh.roughness * textureSample[1], 0, 1),
+    metallic: clamp(mesh.metallic * textureSample[2], 0, 1)
+  };
+}
+function averageColors(colors) {
+  const count = Math.max(colors.length, 1);
+  return colors.reduce(
+    (accumulator, color) => [
+      accumulator[0] + color[0] / count,
+      accumulator[1] + color[1] / count,
+      accumulator[2] + color[2] / count,
+      accumulator[3] + color[3] / count
+    ],
+    [0, 0, 0, 0]
+  );
+}
+function averageNumbers(values, fallback = 0) {
+  if (!Array.isArray(values) || values.length === 0) {
+    return fallback;
+  }
+  return values.reduce((sum, value) => sum + value, 0) / values.length;
+}
 function createMeshTriangleRecords(meshes) {
   const source = Array.isArray(meshes) ? meshes : [];
   let nextTriangleId = 0;
@@ -559,6 +1097,16 @@ function createMeshTriangleRecords(meshes) {
       const uv0 = mesh.uvs ? readVector2(mesh.uvs, a) : [0, 0];
       const uv1 = mesh.uvs ? readVector2(mesh.uvs, b) : [0, 0];
       const uv2 = mesh.uvs ? readVector2(mesh.uvs, c) : [0, 0];
+      const tangentBasis = buildTriangleTangentBasis(v0, v1, v2, uv0, uv1, uv2, faceNormal);
+      const shadedN0 = applyNormalMap(n0, tangentBasis.tangent, tangentBasis.bitangent, mesh.normalTexture, uv0);
+      const shadedN1 = applyNormalMap(n1, tangentBasis.tangent, tangentBasis.bitangent, mesh.normalTexture, uv1);
+      const shadedN2 = applyNormalMap(n2, tangentBasis.tangent, tangentBasis.bitangent, mesh.normalTexture, uv2);
+      const sampledColors = [sampleBaseColor(mesh, uv0), sampleBaseColor(mesh, uv1), sampleBaseColor(mesh, uv2)];
+      const sampledMaterials = [
+        sampleSurfaceMaterial(mesh, uv0),
+        sampleSurfaceMaterial(mesh, uv1),
+        sampleSurfaceMaterial(mesh, uv2)
+      ];
       const bounds = triangleBounds(v0, v1, v2);
       triangles.push(
         Object.freeze({
@@ -568,18 +1116,42 @@ function createMeshTriangleRecords(meshes) {
           flags: mesh.flags,
           materialRefId: mesh.materialRefId,
           mediumRefId: mesh.mediumRefId,
+          materialSlot: meshIndex,
           v0: Object.freeze(v0),
           v1: Object.freeze(v1),
           v2: Object.freeze(v2),
-          n0: Object.freeze(n0),
-          n1: Object.freeze(n1),
-          n2: Object.freeze(n2),
+          n0: Object.freeze(shadedN0),
+          n1: Object.freeze(shadedN1),
+          n2: Object.freeze(shadedN2),
           uv0: Object.freeze(uv0),
           uv1: Object.freeze(uv1),
           uv2: Object.freeze(uv2),
-          color: mesh.color,
+          color: Object.freeze(averageColors(sampledColors)),
           emission: mesh.emission,
-          material: Object.freeze([mesh.roughness, mesh.metallic, mesh.opacity, mesh.ior]),
+          material: Object.freeze([
+            averageNumbers(sampledMaterials.map((sample) => sample.roughness), mesh.roughness),
+            averageNumbers(sampledMaterials.map((sample) => sample.metallic), mesh.metallic),
+            mesh.opacity,
+            mesh.ior
+          ]),
+          materialResponse: Object.freeze([
+            mesh.sheenColor[0] ?? 0,
+            mesh.sheenColor[1] ?? 0,
+            mesh.sheenColor[2] ?? 0,
+            mesh.clearcoat
+          ]),
+          materialExtension: Object.freeze([
+            mesh.clearcoatRoughness,
+            mesh.specular,
+            mesh.transmission,
+            mesh.thickness
+          ]),
+          specularColor: Object.freeze([
+            mesh.specularColor[0] ?? 1,
+            mesh.specularColor[1] ?? 1,
+            mesh.specularColor[2] ?? 1,
+            1
+          ]),
           bounds: Object.freeze({
             min: Object.freeze(bounds.min),
             max: Object.freeze(bounds.max)
@@ -688,50 +1260,264 @@ function nextPowerOfTwo(value) {
   }
   return 2 ** Math.ceil(Math.log2(value));
 }
-function estimateBvhLeafSortCapacity(triangleCount) {
-  return triangleCount <= 0 ? 0 : nextPowerOfTwo(triangleCount);
-}
-function createWavefrontBvhSortStages(itemCountInput) {
-  const itemCount = readNonNegativeInteger("itemCount", itemCountInput, 0);
-  const sortCount = estimateBvhLeafSortCapacity(itemCount);
-  if (sortCount <= 1) {
-    return Object.freeze([]);
+function textureComponentToByte(value, fallback) {
+  const numeric = Number(value);
+  if (!Number.isFinite(numeric)) {
+    return fallback;
   }
-  const stages = [];
-  for (let sequenceSize = 2; sequenceSize <= sortCount; sequenceSize *= 2) {
-    for (let compareDistance = sequenceSize / 2; compareDistance >= 1; compareDistance /= 2) {
-      stages.push(
-        Object.freeze({
-          compareDistance,
-          sequenceSize
-        })
-      );
-    }
+  if (numeric >= 0 && numeric <= 1) {
+    return Math.max(0, Math.min(255, Math.round(numeric * 255)));
   }
-  return Object.freeze(stages);
+  return Math.max(0, Math.min(255, Math.round(numeric)));
 }
-function createWavefrontBvhBuildLevels(triangleCountInput) {
-  const triangleCount = readNonNegativeInteger("triangleCount", triangleCountInput, 0);
-  const internalCount = Math.max(0, triangleCount - 1);
-  if (internalCount === 0) {
-    return Object.freeze([]);
+function createSolidTextureSample(width, height, rgba) {
+  const data = new Uint8Array(width * height * 4);
+  for (let offset = 0; offset < data.length; offset += 4) {
+    data[offset] = rgba[0];
+    data[offset + 1] = rgba[1];
+    data[offset + 2] = rgba[2];
+    data[offset + 3] = rgba[3];
   }
-  const levels = [];
-  let depth = 0;
-  while (Math.pow(2, depth) - 1 < internalCount) {
-    depth += 1;
+  return Object.freeze({
+    width,
+    height,
+    data
+  });
+}
+function normalizeTextureSampleInput(texture, fallbackColor) {
+  if (!texture || !Number.isFinite(texture.width) || !Number.isFinite(texture.height) || texture.width <= 0 || texture.height <= 0) {
+    return createSolidTextureSample(1, 1, fallbackColor);
   }
-  for (let level = depth - 1; level >= 0; level -= 1) {
-    const start = Math.pow(2, level) - 1;
-    const end = Math.min(Math.pow(2, level + 1) - 2, internalCount - 1);
-    if (end >= start) {
-      levels.push(
-        Object.freeze({
-          start,
-          count: end - start + 1
-        })
-      );
-    }
+  const pixelCount = Math.trunc(texture.width) * Math.trunc(texture.height) * 4;
+  const source = ArrayBuffer.isView(texture.data) || Array.isArray(texture.data) ? texture.data : null;
+  if (!source || source.length < pixelCount) {
+    return createSolidTextureSample(1, 1, fallbackColor);
+  }
+  const data = new Uint8Array(pixelCount);
+  for (let index = 0; index < pixelCount; index += 1) {
+    data[index] = textureComponentToByte(source[index], fallbackColor[index % 4]);
+  }
+  return Object.freeze({
+    width: Math.trunc(texture.width),
+    height: Math.trunc(texture.height),
+    data
+  });
+}
+function buildTextureAtlas(textures, fallbackColor) {
+  const padding = 1;
+  const defaultTexture = createSolidTextureSample(1, 1, fallbackColor);
+  const uniqueEntries = [{ source: null, texture: defaultTexture }];
+  const bySource = /* @__PURE__ */ new Map();
+  for (const texture of Array.isArray(textures) ? textures : []) {
+    if (!texture || bySource.has(texture)) {
+      continue;
+    }
+    const normalized = normalizeTextureSampleInput(texture, fallbackColor);
+    bySource.set(texture, uniqueEntries.length);
+    uniqueEntries.push({ source: texture, texture: normalized });
+  }
+  const totalArea = uniqueEntries.reduce((sum, entry) => {
+    return sum + (entry.texture.width + padding * 2) * (entry.texture.height + padding * 2);
+  }, 0);
+  const maxTileWidth = uniqueEntries.reduce((maxWidth, entry) => {
+    return Math.max(maxWidth, entry.texture.width + padding * 2);
+  }, 1);
+  const targetWidth = Math.max(
+    maxTileWidth,
+    nextPowerOfTwo(Math.max(maxTileWidth, Math.ceil(Math.sqrt(totalArea))))
+  );
+  let cursorX = 0;
+  let cursorY = 0;
+  let rowHeight = 0;
+  let atlasWidth = 0;
+  const placements = uniqueEntries.map((entry) => {
+    const tileWidth = entry.texture.width + padding * 2;
+    const tileHeight = entry.texture.height + padding * 2;
+    if (cursorX > 0 && cursorX + tileWidth > targetWidth) {
+      cursorX = 0;
+      cursorY += rowHeight;
+      rowHeight = 0;
+    }
+    const placement = Object.freeze({
+      x: cursorX,
+      y: cursorY,
+      tileWidth,
+      tileHeight,
+      width: entry.texture.width,
+      height: entry.texture.height
+    });
+    cursorX += tileWidth;
+    atlasWidth = Math.max(atlasWidth, cursorX);
+    rowHeight = Math.max(rowHeight, tileHeight);
+    return placement;
+  });
+  const atlasHeight = Math.max(1, cursorY + rowHeight);
+  const atlasData = new Uint8Array(Math.max(1, atlasWidth * atlasHeight * 4));
+  const writePixel = (x, y, rgba) => {
+    const offset = (y * atlasWidth + x) * 4;
+    atlasData[offset] = rgba[0];
+    atlasData[offset + 1] = rgba[1];
+    atlasData[offset + 2] = rgba[2];
+    atlasData[offset + 3] = rgba[3];
+  };
+  const rects = placements.map((placement, entryIndex) => {
+    const { texture } = uniqueEntries[entryIndex];
+    for (let y = 0; y < placement.tileHeight; y += 1) {
+      for (let x = 0; x < placement.tileWidth; x += 1) {
+        const sampleX = Math.max(0, Math.min(texture.width - 1, x - padding));
+        const sampleY = Math.max(0, Math.min(texture.height - 1, y - padding));
+        const sourceOffset = (sampleY * texture.width + sampleX) * 4;
+        writePixel(placement.x + x, placement.y + y, texture.data.slice(sourceOffset, sourceOffset + 4));
+      }
+    }
+    return Object.freeze([
+      (placement.x + padding) / Math.max(1, atlasWidth),
+      (placement.y + padding) / Math.max(1, atlasHeight),
+      placement.width / Math.max(1, atlasWidth),
+      placement.height / Math.max(1, atlasHeight)
+    ]);
+  });
+  const rectBySource = /* @__PURE__ */ new Map();
+  uniqueEntries.forEach((entry, index) => {
+    if (entry.source) {
+      rectBySource.set(entry.source, rects[index]);
+    }
+  });
+  return Object.freeze({
+    width: Math.max(1, atlasWidth),
+    height: Math.max(1, atlasHeight),
+    data: atlasData,
+    defaultRect: rects[0],
+    resolveRect(texture) {
+      return rectBySource.get(texture) ?? rects[0];
+    }
+  });
+}
+function createWavefrontGpuMaterialSource(meshes = []) {
+  const source = Array.isArray(meshes) ? meshes : [meshes];
+  const normalized = source.map((meshInput, meshIndex) => normalizeWavefrontMesh(meshInput, meshIndex));
+  const baseColorAtlas = buildTextureAtlas(
+    normalized.map((mesh) => mesh.baseColorTexture),
+    [255, 255, 255, 255]
+  );
+  const metallicRoughnessAtlas = buildTextureAtlas(
+    normalized.map((mesh) => mesh.metallicRoughnessTexture),
+    [255, 255, 255, 255]
+  );
+  const normalAtlas = buildTextureAtlas(
+    normalized.map((mesh) => mesh.normalTexture),
+    [128, 128, 255, 255]
+  );
+  const occlusionAtlas = buildTextureAtlas(
+    normalized.map((mesh) => mesh.occlusionTexture),
+    [255, 255, 255, 255]
+  );
+  const emissiveAtlas = buildTextureAtlas(
+    normalized.map((mesh) => mesh.emissiveTexture),
+    [255, 255, 255, 255]
+  );
+  const bytes = new ArrayBuffer(Math.max(1, normalized.length) * GPU_MATERIAL_RECORD_BYTES);
+  const floatView = new Float32Array(bytes);
+  normalized.forEach((mesh, meshIndex) => {
+    const byteOffset = meshIndex * GPU_MATERIAL_RECORD_BYTES;
+    writeVec4(floatView, byteOffset, mesh.color);
+    writeVec4(floatView, byteOffset + 16, mesh.emission);
+    writeVec4(floatView, byteOffset + 32, [
+      mesh.roughness,
+      mesh.metallic,
+      mesh.opacity,
+      mesh.ior
+    ]);
+    writeVec4(floatView, byteOffset + 48, [
+      mesh.sheenColor[0] ?? 0,
+      mesh.sheenColor[1] ?? 0,
+      mesh.sheenColor[2] ?? 0,
+      mesh.clearcoat
+    ]);
+    writeVec4(floatView, byteOffset + 64, [
+      mesh.clearcoatRoughness,
+      mesh.specular,
+      mesh.transmission,
+      mesh.thickness
+    ]);
+    writeVec4(floatView, byteOffset + 80, [
+      mesh.specularColor[0] ?? 1,
+      mesh.specularColor[1] ?? 1,
+      mesh.specularColor[2] ?? 1,
+      1
+    ]);
+    writeVec4(floatView, byteOffset + 96, baseColorAtlas.resolveRect(mesh.baseColorTexture));
+    writeVec4(
+      floatView,
+      byteOffset + 112,
+      metallicRoughnessAtlas.resolveRect(mesh.metallicRoughnessTexture)
+    );
+    writeVec4(floatView, byteOffset + 128, normalAtlas.resolveRect(mesh.normalTexture));
+    writeVec4(floatView, byteOffset + 144, occlusionAtlas.resolveRect(mesh.occlusionTexture));
+    writeVec4(floatView, byteOffset + 160, emissiveAtlas.resolveRect(mesh.emissiveTexture));
+    writeVec4(floatView, byteOffset + 176, [
+      clampUnit(mesh.normalTexture?.scale ?? mesh.normalTexture?.strength ?? 1),
+      clampUnit(mesh.occlusionTexture?.strength ?? 1),
+      clampUnit(mesh.emissiveTexture?.strength ?? 1),
+      0
+    ]);
+  });
+  return Object.freeze({
+    buffer: bytes,
+    count: normalized.length,
+    recordBytes: GPU_MATERIAL_RECORD_BYTES,
+    records: Object.freeze(normalized),
+    baseColorAtlas,
+    metallicRoughnessAtlas,
+    normalAtlas,
+    occlusionAtlas,
+    emissiveAtlas
+  });
+}
+function estimateBvhLeafSortCapacity(triangleCount) {
+  return triangleCount <= 0 ? 0 : nextPowerOfTwo(triangleCount);
+}
+function createWavefrontBvhSortStages(itemCountInput) {
+  const itemCount = readNonNegativeInteger("itemCount", itemCountInput, 0);
+  const sortCount = estimateBvhLeafSortCapacity(itemCount);
+  if (sortCount <= 1) {
+    return Object.freeze([]);
+  }
+  const stages = [];
+  for (let sequenceSize = 2; sequenceSize <= sortCount; sequenceSize *= 2) {
+    for (let compareDistance = sequenceSize / 2; compareDistance >= 1; compareDistance /= 2) {
+      stages.push(
+        Object.freeze({
+          compareDistance,
+          sequenceSize
+        })
+      );
+    }
+  }
+  return Object.freeze(stages);
+}
+function createWavefrontBvhBuildLevels(triangleCountInput) {
+  const triangleCount = readNonNegativeInteger("triangleCount", triangleCountInput, 0);
+  const internalCount = Math.max(0, triangleCount - 1);
+  if (internalCount === 0) {
+    return Object.freeze([]);
+  }
+  const levels = [];
+  let depth = 0;
+  while (Math.pow(2, depth) - 1 < internalCount) {
+    depth += 1;
+  }
+  for (let level = depth - 1; level >= 0; level -= 1) {
+    const start = Math.pow(2, level) - 1;
+    const end = Math.min(Math.pow(2, level + 1) - 2, internalCount - 1);
+    if (end >= start) {
+      levels.push(
+        Object.freeze({
+          start,
+          count: end - start + 1
+        })
+      );
+    }
   }
   return Object.freeze(levels);
 }
@@ -745,9 +1531,10 @@ function resolveAccelerationBuildMode(options = {}) {
   }
   return mode;
 }
-function createWavefrontGpuMeshSource(meshes = []) {
+function createWavefrontGpuMeshSource(meshes = [], gpuMaterialSourceInput = null) {
   const source = Array.isArray(meshes) ? meshes : [meshes];
   const normalized = source.map((meshInput, meshIndex) => normalizeWavefrontMesh(meshInput, meshIndex));
+  const gpuMaterialSource = gpuMaterialSourceInput ?? createWavefrontGpuMaterialSource(normalized);
   const vertexCount = normalized.reduce((count, mesh) => count + mesh.positions.length / 3, 0);
   const indexCount = normalized.reduce((count, mesh) => count + mesh.indices.length, 0);
   const triangleCount = Math.floor(indexCount / 3);
@@ -799,7 +1586,7 @@ function createWavefrontGpuMeshSource(meshes = []) {
     meshUints[meshOffset + 8] = mesh.indices.length / 3;
     meshUints[meshOffset + 9] = meshVertexBase;
     meshUints[meshOffset + 10] = meshVertexCount;
-    meshUints[meshOffset + 11] = 0;
+    meshUints[meshOffset + 11] = meshIndex;
     const floatOffset = meshOffset;
     writeVec4(meshFloats, floatOffset * 4 + 48, mesh.color);
     writeVec4(meshFloats, floatOffset * 4 + 64, mesh.emission);
@@ -809,6 +1596,55 @@ function createWavefrontGpuMeshSource(meshes = []) {
       mesh.opacity,
       mesh.ior
     ]);
+    writeVec4(meshFloats, floatOffset * 4 + 96, [
+      mesh.sheenColor[0] ?? 0,
+      mesh.sheenColor[1] ?? 0,
+      mesh.sheenColor[2] ?? 0,
+      mesh.clearcoat
+    ]);
+    writeVec4(meshFloats, floatOffset * 4 + 112, [
+      mesh.clearcoatRoughness,
+      mesh.specular,
+      mesh.transmission,
+      mesh.thickness
+    ]);
+    writeVec4(meshFloats, floatOffset * 4 + 128, [
+      mesh.specularColor[0] ?? 1,
+      mesh.specularColor[1] ?? 1,
+      mesh.specularColor[2] ?? 1,
+      1
+    ]);
+    writeVec4(
+      meshFloats,
+      floatOffset * 4 + 144,
+      gpuMaterialSource.baseColorAtlas.resolveRect(mesh.baseColorTexture)
+    );
+    writeVec4(
+      meshFloats,
+      floatOffset * 4 + 160,
+      gpuMaterialSource.metallicRoughnessAtlas.resolveRect(mesh.metallicRoughnessTexture)
+    );
+    writeVec4(
+      meshFloats,
+      floatOffset * 4 + 176,
+      gpuMaterialSource.normalAtlas.resolveRect(mesh.normalTexture)
+    );
+    writeVec4(
+      meshFloats,
+      floatOffset * 4 + 192,
+      gpuMaterialSource.occlusionAtlas.resolveRect(mesh.occlusionTexture)
+    );
+    writeVec4(
+      meshFloats,
+      floatOffset * 4 + 208,
+      gpuMaterialSource.emissiveAtlas.resolveRect(mesh.emissiveTexture)
+    );
+    writeVec4(meshFloats, floatOffset * 4 + 224, [
+      clampUnit(mesh.normalTexture?.scale ?? mesh.normalTexture?.strength ?? 1),
+      clampUnit(mesh.occlusionTexture?.strength ?? 1),
+      clampUnit(mesh.emissiveTexture?.strength ?? 1),
+      0
+    ]);
     vertexCursor += meshVertexCount;
     indexCursor += mesh.indices.length;
     triangleCursor += mesh.indices.length / 3;
@@ -871,12 +1707,16 @@ function normalizeSceneObjects(sceneObjects, useDefaultScene = true) {
   const source = Array.isArray(sceneObjects) && sceneObjects.length > 0 ? sceneObjects : useDefaultScene ? createDefaultWavefrontSceneObjects() : [];
   return source.map((object, index) => normalizeWavefrontSceneObject(object, index));
 }
+function normalizeWavefrontMeshes(meshes) {
+  const source = Array.isArray(meshes) ? meshes : [];
+  return source.map((mesh, index) => normalizeWavefrontMesh(mesh, index));
+}
 function normalizeMeshes(options = {}) {
   if (Array.isArray(options.meshes)) {
-    return options.meshes;
+    return normalizeWavefrontMeshes(options.meshes);
   }
   if (options.mesh) {
-    return [options.mesh];
+    return normalizeWavefrontMeshes([options.mesh]);
   }
   return [];
 }
@@ -1111,12 +1951,14 @@ function estimateWavefrontPathTracingMemory(options = {}) {
     options.environmentPortalCapacity,
     0
   );
+  const materialCapacity = readNonNegativeInteger("materialCapacity", options.materialCapacity, 0);
   const queueBytes = tilePixelCapacity * RAY_RECORD_BYTES;
   const hitBytes = tilePixelCapacity * HIT_RECORD_BYTES;
   const accumulationBytes = tilePixelCapacity * ACCUMULATION_RECORD_BYTES;
   const pathVertexBytes = tilePixelCapacity * (maxDepth + 1) * PATH_VERTEX_RECORD_BYTES;
   const sceneObjectBytes = sceneObjectCapacity * SCENE_OBJECT_RECORD_BYTES;
   const triangleBytes = triangleCapacity * TRIANGLE_RECORD_BYTES;
+  const materialTableBytes = materialCapacity * GPU_MATERIAL_RECORD_BYTES;
   const bvhNodeBytes = bvhNodeCapacity * BVH_NODE_RECORD_BYTES;
   const bvhLeafReferenceBytes = bvhLeafSortCapacity * BVH_LEAF_REF_RECORD_BYTES;
   const emissiveTriangleMetadataBytes = emissiveTriangleCapacity * BVH_NODE_RECORD_BYTES;
@@ -1129,6 +1971,7 @@ function estimateWavefrontPathTracingMemory(options = {}) {
     pathVertexBytes,
     sceneObjectBytes,
     triangleBytes,
+    materialTableBytes,
     bvhNodeBytes,
     bvhLeafReferenceBytes,
     emissiveTriangleMetadataBytes,
@@ -1136,7 +1979,7 @@ function estimateWavefrontPathTracingMemory(options = {}) {
     configBytes: CONFIG_BUFFER_BYTES,
     counterBytes: COUNTER_BUFFER_BYTES,
     indirectDispatchBytes: INDIRECT_DISPATCH_ARGS_BYTES,
-    totalHotBufferBytes: queueBytes * 2 + hitBytes + accumulationBytes + pathVertexBytes + sceneObjectBytes + triangleBytes + bvhNodeBytes + bvhLeafReferenceBytes + emissiveTriangleMetadataBytes + environmentPortalBytes + CONFIG_BUFFER_BYTES + COUNTER_BUFFER_BYTES + INDIRECT_DISPATCH_ARGS_BYTES
+    totalHotBufferBytes: queueBytes * 2 + hitBytes + accumulationBytes + pathVertexBytes + sceneObjectBytes + triangleBytes + materialTableBytes + bvhNodeBytes + bvhLeafReferenceBytes + emissiveTriangleMetadataBytes + environmentPortalBytes + CONFIG_BUFFER_BYTES + COUNTER_BUFFER_BYTES + INDIRECT_DISPATCH_ARGS_BYTES
   });
 }
 function createWavefrontPathTracingComputeConfig(options = {}) {
@@ -1150,7 +1993,7 @@ function createWavefrontPathTracingComputeConfig(options = {}) {
   const samplesPerPixel = clamp(
     readPositiveInteger("samplesPerPixel", options.samplesPerPixel, DEFAULT_SAMPLES_PER_PIXEL),
     1,
-    64
+    MAX_SAMPLES_PER_PIXEL
   );
   const maxFramePassesPerSubmission = clamp(
     readPositiveInteger(
@@ -1168,7 +2011,8 @@ function createWavefrontPathTracingComputeConfig(options = {}) {
   );
   const meshes = normalizeMeshes(options);
   const meshSourceShape = estimateMeshSourceShape(meshes);
-  const gpuMeshSource = meshes.length > 0 ? createWavefrontGpuMeshSource(meshes) : createWavefrontGpuMeshSource([]);
+  const gpuMaterialSource = meshes.length > 0 ? createWavefrontGpuMaterialSource(meshes) : createWavefrontGpuMaterialSource([]);
+  const gpuMeshSource = meshes.length > 0 ? createWavefrontGpuMeshSource(meshes, gpuMaterialSource) : createWavefrontGpuMeshSource([]);
   const meshAcceleration = accelerationBuildMode === "cpu-debug" ? createWavefrontMeshAcceleration(meshes) : Object.freeze({ nodes: Object.freeze([]), triangles: Object.freeze([]) });
   const emissiveTriangleIndices = createWavefrontEmissiveTriangleIndexSource(
     meshes,
@@ -1179,6 +2023,7 @@ function createWavefrontPathTracingComputeConfig(options = {}) {
   const sceneObjects = Object.freeze(
     normalizeSceneObjects(options.sceneObjects, meshes.length === 0)
   );
+  const mediums = collectWavefrontMediums(options, meshes, sceneObjects);
   const sceneObjectCapacity = Math.max(
     sceneObjects.length,
     readPositiveInteger("sceneObjectCapacity", options.sceneObjectCapacity, DEFAULT_SCENE_OBJECT_CAPACITY)
@@ -1237,9 +2082,12 @@ function createWavefrontPathTracingComputeConfig(options = {}) {
     sceneObjects,
     sceneObjectCount: sceneObjects.length,
     sceneObjectCapacity,
+    mediums,
+    mediumCount: mediums.length,
     accelerationBuildMode,
     gpuAccelerationBuildRequired: accelerationBuildMode === "gpu" && triangleCount > 0,
     gpuMeshSource,
+    gpuMaterialSource,
     meshAcceleration,
     emissiveTriangleIndices,
     emissiveTriangleCount: emissiveTriangleIndices.count,
@@ -1270,6 +2118,7 @@ function createWavefrontPathTracingComputeConfig(options = {}) {
       maxDepth,
       sceneObjectCapacity,
       triangleCapacity,
+      materialCapacity: gpuMaterialSource.count,
       bvhNodeCapacity,
       bvhLeafSortCapacity,
       emissiveTriangleCapacity: emissiveTriangleIndices.capacity,
@@ -1336,16 +2185,35 @@ function packWavefrontSceneObjects(sceneObjects, capacity = sceneObjects.length)
     uintView[u32 + 1] = object.id;
     uintView[u32 + 2] = object.materialKind;
     uintView[u32 + 3] = object.flags;
-    writeVec4(floatView, byteOffset + 16, [...object.center, 0]);
-    writeVec4(floatView, byteOffset + 32, [...object.halfExtent, 0]);
-    writeVec4(floatView, byteOffset + 48, object.color);
-    writeVec4(floatView, byteOffset + 64, object.emission);
-    writeVec4(floatView, byteOffset + 80, [
+    uintView[u32 + 4] = object.mediumRefId;
+    writeVec4(floatView, byteOffset + 32, [...object.center, 0]);
+    writeVec4(floatView, byteOffset + 48, [...object.halfExtent, 0]);
+    writeVec4(floatView, byteOffset + 64, object.color);
+    writeVec4(floatView, byteOffset + 80, object.emission);
+    writeVec4(floatView, byteOffset + 96, [
       object.roughness,
       object.metallic,
       object.opacity,
       object.ior
     ]);
+    writeVec4(floatView, byteOffset + 112, [
+      object.sheenColor[0] ?? 0,
+      object.sheenColor[1] ?? 0,
+      object.sheenColor[2] ?? 0,
+      object.clearcoat
+    ]);
+    writeVec4(floatView, byteOffset + 128, [
+      object.clearcoatRoughness,
+      object.specular,
+      object.transmission,
+      object.thickness
+    ]);
+    writeVec4(floatView, byteOffset + 144, [
+      object.specularColor[0] ?? 1,
+      object.specularColor[1] ?? 1,
+      object.specularColor[2] ?? 1,
+      1
+    ]);
   });
   return Object.freeze({
     buffer: bytes,
@@ -1370,7 +2238,7 @@ function packWavefrontTriangles(triangles, capacity = triangles.length) {
     uintView[u32 + 3] = triangle.flags;
     uintView[u32 + 4] = triangle.materialRefId;
     uintView[u32 + 5] = triangle.mediumRefId;
-    uintView[u32 + 6] = 0;
+    uintView[u32 + 6] = triangle.materialSlot ?? 0;
     uintView[u32 + 7] = 0;
     writeVec4(floatView, byteOffset + 32, [...triangle.v0, 0]);
     writeVec4(floatView, byteOffset + 48, [...triangle.v1, 0]);
@@ -1383,6 +2251,15 @@ function packWavefrontTriangles(triangles, capacity = triangles.length) {
     writeVec4(floatView, byteOffset + 160, triangle.color);
     writeVec4(floatView, byteOffset + 176, triangle.emission);
     writeVec4(floatView, byteOffset + 192, triangle.material);
+    writeVec4(floatView, byteOffset + 208, triangle.materialResponse);
+    writeVec4(floatView, byteOffset + 224, triangle.materialExtension ?? [0.08, 1, 0, 0]);
+    writeVec4(floatView, byteOffset + 240, triangle.specularColor ?? [1, 1, 1, 1]);
+    writeVec4(floatView, byteOffset + 256, triangle.baseColorAtlas ?? [0, 0, 1, 1]);
+    writeVec4(floatView, byteOffset + 272, triangle.metallicRoughnessAtlas ?? [0, 0, 1, 1]);
+    writeVec4(floatView, byteOffset + 288, triangle.normalAtlas ?? [0, 0, 1, 1]);
+    writeVec4(floatView, byteOffset + 304, triangle.occlusionAtlas ?? [0, 0, 1, 1]);
+    writeVec4(floatView, byteOffset + 320, triangle.emissiveAtlas ?? [0, 0, 1, 1]);
+    writeVec4(floatView, byteOffset + 336, triangle.textureSettings ?? [1, 1, 1, 0]);
   });
   return Object.freeze({
     buffer: bytes,
@@ -1476,6 +2353,12 @@ function createConfigPayload(config, tile, frameIndex, buildRange = {}) {
     0,
     0
   ]);
+  writeVec4(floatView, 304, [
+    config.environmentMap.width ?? 1,
+    config.environmentMap.height ?? 1,
+    config.environmentMap.mipLevelCount ?? 1,
+    config.environmentMap.hasImportanceData ? 1 : 0
+  ]);
   return bytes;
 }
 function createTiles(width, height, tileSize) {
@@ -1649,7 +2532,8 @@ function intersectWavefrontReferenceTriangle(ray, triangle, options = {}) {
     position: Object.freeze(position),
     color: triangle.color,
     emission: triangle.emission,
-    material: triangle.material
+    material: triangle.material,
+    materialResponse: triangle.materialResponse
   });
 }
 function createWavefrontReferenceEnvironmentHit(config, ray) {
@@ -1675,7 +2559,8 @@ function createWavefrontReferenceEnvironmentHit(config, ray) {
     position: Object.freeze(add(ray.origin, scale(ray.direction, 1e3))),
     color: Object.freeze([0, 0, 0, 0]),
     emission: radiance,
-    material: Object.freeze([1, 0, 1, 1])
+    material: Object.freeze([1, 0, 1, 1]),
+    materialResponse: Object.freeze([0, 0, 0, 0])
   });
 }
 function traceWavefrontReferenceTriangles(config, ray, triangles, options = {}) {
@@ -1754,6 +2639,32 @@ function environmentMapIntegerScale(data) {
   }
   return 1;
 }
+function environmentMapHasSamplingData(environmentMap) {
+  if (!environmentMap || !environmentMap.data) {
+    return false;
+  }
+  const width = Math.max(1, environmentMap.width ?? 1);
+  const height = Math.max(1, environmentMap.height ?? 1);
+  return environmentMap.data.length >= width * height * 4;
+}
+function createRgba8TextureUpload(source) {
+  const width = Math.max(1, Math.trunc(source.width));
+  const height = Math.max(1, Math.trunc(source.height));
+  const bytesPerRow = alignTo(width * 4, 256);
+  const bytes = new Uint8Array(bytesPerRow * height);
+  const data = source.data instanceof Uint8Array ? source.data : new Uint8Array(source.data);
+  for (let y = 0; y < height; y += 1) {
+    const sourceOffset = y * width * 4;
+    const targetOffset = y * bytesPerRow;
+    bytes.set(data.subarray(sourceOffset, sourceOffset + width * 4), targetOffset);
+  }
+  return Object.freeze({
+    bytes,
+    bytesPerRow,
+    width,
+    height
+  });
+}
 function readEnvironmentMapComponent(data, index, fallback, integerScale = 1) {
   if (!data || index >= data.length) {
     return fallback;
@@ -1761,39 +2672,311 @@ function readEnvironmentMapComponent(data, index, fallback, integerScale = 1) {
   const value = Number(data[index]);
   return Number.isFinite(value) ? Math.max(0, value) * integerScale : fallback;
 }
-function createEnvironmentMapUploadBytes(environmentMap, fallbackColor) {
-  const width = Math.max(1, environmentMap.width);
-  const height = Math.max(1, environmentMap.height);
+function buildOrthonormalBasis(normal) {
+  const tangentFallback = Math.abs(normal[1]) < 0.999 ? [0, 1, 0] : [1, 0, 0];
+  const tangent = normalize(cross(tangentFallback, normal), [1, 0, 0]);
+  const bitangent = normalize(cross(normal, tangent), [0, 0, 1]);
+  return { tangent, bitangent };
+}
+function localToWorld(local, normal) {
+  const basis = buildOrthonormalBasis(normal);
+  return normalize(
+    add(
+      add(scale(basis.tangent, local[0]), scale(basis.bitangent, local[1])),
+      scale(normal, local[2])
+    ),
+    normal
+  );
+}
+function radicalInverseVdc(bits) {
+  let value = bits >>> 0;
+  value = (value << 16 | value >>> 16) >>> 0;
+  value = ((value & 1431655765) << 1 | (value & 2863311530) >>> 1) >>> 0;
+  value = ((value & 858993459) << 2 | (value & 3435973836) >>> 2) >>> 0;
+  value = ((value & 252645135) << 4 | (value & 4042322160) >>> 4) >>> 0;
+  value = ((value & 16711935) << 8 | (value & 4278255360) >>> 8) >>> 0;
+  return value * 23283064365386963e-26;
+}
+function hammersley(index, count) {
+  return [index / Math.max(count, 1), radicalInverseVdc(index)];
+}
+function importanceSampleGgx(sample, roughness, normal) {
+  const alpha = Math.max(roughness * roughness, 1e-4);
+  const phi = 2 * Math.PI * sample[0];
+  const cosTheta = Math.sqrt((1 - sample[1]) / (1 + (alpha * alpha - 1) * sample[1]));
+  const sinTheta = Math.sqrt(Math.max(0, 1 - cosTheta * cosTheta));
+  const halfVector = localToWorld(
+    [Math.cos(phi) * sinTheta, Math.sin(phi) * sinTheta, cosTheta],
+    normal
+  );
+  return normalize(halfVector, normal);
+}
+function geometrySchlickGgx(nDotV, roughness) {
+  const k = (roughness + 1) * (roughness + 1) / 8;
+  return nDotV / Math.max(nDotV * (1 - k) + k, 1e-6);
+}
+function geometrySmith(nDotV, nDotL, roughness) {
+  return geometrySchlickGgx(nDotV, roughness) * geometrySchlickGgx(nDotL, roughness);
+}
+function integrateBrdfSample(nDotV, roughness, sampleCount) {
+  const viewDirection = [Math.sqrt(Math.max(0, 1 - nDotV * nDotV)), 0, nDotV];
+  const normal = [0, 0, 1];
+  let scaleTerm = 0;
+  let biasTerm = 0;
+  for (let index = 0; index < sampleCount; index += 1) {
+    const xi = hammersley(index, sampleCount);
+    const halfVector = importanceSampleGgx(xi, roughness, normal);
+    const vDotH = Math.max(dot(viewDirection, halfVector), 0);
+    const lightDirection = normalize(
+      subtract(scale(halfVector, 2 * vDotH), viewDirection),
+      normal
+    );
+    const nDotL = Math.max(lightDirection[2], 0);
+    const nDotH = Math.max(halfVector[2], 0);
+    if (nDotL <= 0 || nDotH <= 0 || vDotH <= 0) {
+      continue;
+    }
+    const geometry = geometrySmith(nDotV, nDotL, roughness);
+    const visibility = geometry * vDotH / Math.max(nDotH * nDotV, 1e-6);
+    const fresnel = (1 - vDotH) ** 5;
+    scaleTerm += (1 - fresnel) * visibility;
+    biasTerm += fresnel * visibility;
+  }
+  return [scaleTerm / sampleCount, biasTerm / sampleCount];
+}
+function createBrdfLutUploadBytes(size = DEFAULT_BRDF_LUT_SIZE, sampleCount = DEFAULT_BRDF_LUT_SAMPLE_COUNT) {
+  const cacheKey = `${Math.max(1, Math.trunc(size))}:${Math.max(1, Math.trunc(sampleCount))}`;
+  const cached = BRDF_LUT_UPLOAD_CACHE.get(cacheKey);
+  if (cached) {
+    return cached;
+  }
+  const width = Math.max(1, Math.trunc(size));
+  const height = Math.max(1, Math.trunc(size));
   const rowBytes = width * 8;
   const bytesPerRow = alignTo(rowBytes, 256);
   const bytes = new Uint8Array(bytesPerRow * height);
+  const view = new DataView(bytes.buffer);
+  for (let y = 0; y < height; y += 1) {
+    const roughness = (y + 0.5) / height;
+    for (let x = 0; x < width; x += 1) {
+      const nDotV = Math.max((x + 0.5) / width, 1e-4);
+      const [scaleTerm, biasTerm] = integrateBrdfSample(nDotV, roughness, sampleCount);
+      const offset = y * bytesPerRow + x * 8;
+      view.setUint16(offset, float32ToFloat16Bits(scaleTerm), true);
+      view.setUint16(offset + 2, float32ToFloat16Bits(biasTerm), true);
+      view.setUint16(offset + 4, float32ToFloat16Bits(0), true);
+      view.setUint16(offset + 6, float32ToFloat16Bits(1), true);
+    }
+  }
+  const upload = Object.freeze({ bytes, bytesPerRow, width, height });
+  BRDF_LUT_UPLOAD_CACHE.set(cacheKey, upload);
+  return upload;
+}
+function createLinearEnvironmentPixels(environmentMap, fallbackColor) {
+  const width = Math.max(1, environmentMap.width);
+  const height = Math.max(1, environmentMap.height);
+  const pixels = new Float32Array(width * height * 4);
   const data = environmentMap.data;
   const integerScale = environmentMapIntegerScale(data);
-  const view = new DataView(bytes.buffer);
-  const writeComponent = (targetOffset, sourceOffset, fallback) => {
-    view.setUint16(
-      targetOffset,
-      float32ToFloat16Bits(
-        readEnvironmentMapComponent(data, sourceOffset, fallback, integerScale)
-      ),
-      true
-    );
+  for (let index = 0; index < width * height; index += 1) {
+    const sourceOffset = index * 4;
+    const targetOffset = index * 4;
+    pixels[targetOffset] = readEnvironmentMapComponent(data, sourceOffset, fallbackColor[0], integerScale);
+    pixels[targetOffset + 1] = readEnvironmentMapComponent(data, sourceOffset + 1, fallbackColor[1], integerScale);
+    pixels[targetOffset + 2] = readEnvironmentMapComponent(data, sourceOffset + 2, fallbackColor[2], integerScale);
+    pixels[targetOffset + 3] = readEnvironmentMapComponent(data, sourceOffset + 3, fallbackColor[3] ?? 1, integerScale);
+  }
+  return pixels;
+}
+function environmentUvToDirection(u, v, rotationRadians = 0) {
+  const angle = (u - rotationRadians / (2 * Math.PI) - 0.5) * 2 * Math.PI;
+  const theta = v * Math.PI;
+  const sinTheta = Math.sin(theta);
+  return [
+    Math.cos(angle) * sinTheta,
+    Math.cos(theta),
+    Math.sin(angle) * sinTheta
+  ];
+}
+function sampleEnvironmentPixelsBilinear(pixels, width, height, u, v) {
+  const wrappedU = (u % 1 + 1) % 1;
+  const clampedV = clamp(v, 0, 1);
+  const x = wrappedU * width - 0.5;
+  const y = clampedV * height - 0.5;
+  const x0 = (Math.floor(x) % width + width) % width;
+  const y0 = clamp(Math.floor(y), 0, height - 1);
+  const x1 = (x0 + 1) % width;
+  const y1 = clamp(y0 + 1, 0, height - 1);
+  const tx = x - Math.floor(x);
+  const ty = y - Math.floor(y);
+  const read = (px, py) => {
+    const offset = (py * width + px) * 4;
+    return [pixels[offset], pixels[offset + 1], pixels[offset + 2], pixels[offset + 3]];
   };
+  const a = read(x0, y0);
+  const b = read(x1, y0);
+  const c = read(x0, y1);
+  const d = read(x1, y1);
+  const mixPair = (first, second, factor) => first * (1 - factor) + second * factor;
+  return [
+    mixPair(mixPair(a[0], b[0], tx), mixPair(c[0], d[0], tx), ty),
+    mixPair(mixPair(a[1], b[1], tx), mixPair(c[1], d[1], tx), ty),
+    mixPair(mixPair(a[2], b[2], tx), mixPair(c[2], d[2], tx), ty),
+    mixPair(mixPair(a[3], b[3], tx), mixPair(c[3], d[3], tx), ty)
+  ];
+}
+function directionToEnvironmentUv(direction, rotationRadians = 0) {
+  const unitDirection = normalize(direction, [0, 1, 0]);
+  const rotationTurns = rotationRadians / (2 * Math.PI);
+  const u = ((Math.atan2(unitDirection[2], unitDirection[0]) / (2 * Math.PI) + 0.5 + rotationTurns) % 1 + 1) % 1;
+  const v = Math.acos(clamp(unitDirection[1], -1, 1)) / Math.PI;
+  return [u, clamp(v, 0, 1)];
+}
+function sampleEnvironmentRadiance(pixels, width, height, direction, rotationRadians = 0) {
+  const [u, v] = directionToEnvironmentUv(direction, rotationRadians);
+  return sampleEnvironmentPixelsBilinear(pixels, width, height, u, v);
+}
+function createFloat16RgbaUploadFromLevels(levels) {
+  return levels.map((level) => {
+    const rowBytes = level.width * 8;
+    const bytesPerRow = alignTo(rowBytes, 256);
+    const bytes = new Uint8Array(bytesPerRow * level.height);
+    const view = new DataView(bytes.buffer);
+    for (let y = 0; y < level.height; y += 1) {
+      for (let x = 0; x < level.width; x += 1) {
+        const sourceOffset = (y * level.width + x) * 4;
+        const targetOffset = y * bytesPerRow + x * 8;
+        view.setUint16(targetOffset, float32ToFloat16Bits(level.data[sourceOffset]), true);
+        view.setUint16(targetOffset + 2, float32ToFloat16Bits(level.data[sourceOffset + 1]), true);
+        view.setUint16(targetOffset + 4, float32ToFloat16Bits(level.data[sourceOffset + 2]), true);
+        view.setUint16(targetOffset + 6, float32ToFloat16Bits(level.data[sourceOffset + 3]), true);
+      }
+    }
+    return Object.freeze({ bytes, bytesPerRow, width: level.width, height: level.height });
+  });
+}
+function createPrefilteredEnvironmentLevels(environmentMap, fallbackColor) {
+  const sourcePixels = createLinearEnvironmentPixels(environmentMap, fallbackColor);
+  const sourceWidth = Math.max(1, environmentMap.width);
+  const sourceHeight = Math.max(1, environmentMap.height);
+  const mipLevelCount = Math.max(1, Math.floor(Math.log2(Math.max(sourceWidth, sourceHeight))) + 1);
+  const levels = [
+    Object.freeze({
+      width: sourceWidth,
+      height: sourceHeight,
+      data: sourcePixels
+    })
+  ];
+  for (let mipLevel = 1; mipLevel < mipLevelCount; mipLevel += 1) {
+    const width = Math.max(1, sourceWidth >> mipLevel);
+    const height = Math.max(1, sourceHeight >> mipLevel);
+    const roughness = mipLevelCount <= 1 ? 0 : mipLevel / (mipLevelCount - 1);
+    const data = new Float32Array(width * height * 4);
+    const sampleCount = roughness < 0.25 ? 64 : roughness < 0.6 ? 96 : 128;
+    for (let y = 0; y < height; y += 1) {
+      for (let x = 0; x < width; x += 1) {
+        const direction = environmentUvToDirection((x + 0.5) / width, (y + 0.5) / height, environmentMap.rotationRadians);
+        const normal = normalize(direction, [0, 1, 0]);
+        const viewDirection = normal;
+        let totalWeight = 0;
+        const accum = [0, 0, 0];
+        for (let sampleIndex = 0; sampleIndex < sampleCount; sampleIndex += 1) {
+          const xi = hammersley(sampleIndex, sampleCount);
+          const halfVector = importanceSampleGgx(xi, roughness, normal);
+          const viewDotHalf = Math.max(dot(viewDirection, halfVector), 0);
+          const lightDirection = normalize(
+            subtract(scale(halfVector, 2 * viewDotHalf), viewDirection),
+            normal
+          );
+          const nDotL = Math.max(dot(normal, lightDirection), 0);
+          if (nDotL <= 1e-6) {
+            continue;
+          }
+          const radiance = sampleEnvironmentRadiance(
+            sourcePixels,
+            sourceWidth,
+            sourceHeight,
+            lightDirection,
+            environmentMap.rotationRadians
+          );
+          accum[0] += radiance[0] * nDotL;
+          accum[1] += radiance[1] * nDotL;
+          accum[2] += radiance[2] * nDotL;
+          totalWeight += nDotL;
+        }
+        const offset = (y * width + x) * 4;
+        data[offset] = accum[0] / Math.max(totalWeight, 1e-6);
+        data[offset + 1] = accum[1] / Math.max(totalWeight, 1e-6);
+        data[offset + 2] = accum[2] / Math.max(totalWeight, 1e-6);
+        data[offset + 3] = 1;
+      }
+    }
+    levels.push(Object.freeze({ width, height, data }));
+  }
+  return Object.freeze({
+    levels,
+    mipLevelCount,
+    width: sourceWidth,
+    height: sourceHeight
+  });
+}
+function createEnvironmentSamplingTables(environmentMap, fallbackColor) {
+  if (!environmentMapHasSamplingData(environmentMap)) {
+    return Object.freeze({
+      width: 1,
+      height: 1,
+      pdf: new Float32Array([1]),
+      marginalCdf: new Float32Array([1]),
+      conditionalCdf: new Float32Array([1]),
+      hasImportanceData: false
+    });
+  }
+  const pixels = createLinearEnvironmentPixels(environmentMap, fallbackColor);
+  const width = Math.max(1, environmentMap.width);
+  const height = Math.max(1, environmentMap.height);
+  const pdf = new Float32Array(width * height);
+  const marginalCdf = new Float32Array(height);
+  const conditionalCdf = new Float32Array(width * height);
+  const rowSums = new Float32Array(height);
+  let totalWeight = 0;
   for (let y = 0; y < height; y += 1) {
+    const theta = (y + 0.5) / height * Math.PI;
+    const sinTheta = Math.max(Math.sin(theta), 1e-4);
+    let rowWeight = 0;
     for (let x = 0; x < width; x += 1) {
-      const sourceOffset = (y * width + x) * 4;
-      const targetOffset = y * bytesPerRow + x * 8;
-      writeComponent(targetOffset, sourceOffset, fallbackColor[0]);
-      writeComponent(targetOffset + 2, sourceOffset + 1, fallbackColor[1]);
-      writeComponent(targetOffset + 4, sourceOffset + 2, fallbackColor[2]);
-      writeComponent(targetOffset + 6, sourceOffset + 3, fallbackColor[3] ?? 1);
+      const offset = (y * width + x) * 4;
+      const luminance = pixels[offset] * 0.2126 + pixels[offset + 1] * 0.7152 + pixels[offset + 2] * 0.0722;
+      const weight = Math.max(luminance * sinTheta, 1e-6);
+      pdf[y * width + x] = weight;
+      rowWeight += weight;
+      conditionalCdf[y * width + x] = rowWeight;
+    }
+    rowSums[y] = rowWeight;
+    totalWeight += rowWeight;
+    if (rowWeight > 0) {
+      for (let x = 0; x < width; x += 1) {
+        conditionalCdf[y * width + x] /= rowWeight;
+      }
+    } else {
+      for (let x = 0; x < width; x += 1) {
+        conditionalCdf[y * width + x] = (x + 1) / width;
+      }
     }
+    marginalCdf[y] = totalWeight;
+  }
+  for (let y = 0; y < height; y += 1) {
+    marginalCdf[y] /= Math.max(totalWeight, 1e-6);
+  }
+  for (let index = 0; index < pdf.length; index += 1) {
+    pdf[index] /= Math.max(totalWeight, 1e-6);
   }
   return Object.freeze({
-    bytes,
-    bytesPerRow,
     width,
-    height
+    height,
+    pdf,
+    marginalCdf,
+    conditionalCdf,
+    hasImportanceData: true
   });
 }
 function createEnvironmentMapResource(device, constants, environmentMap, fallbackColor) {
@@ -1805,10 +2988,14 @@ function createEnvironmentMapResource(device, constants, environmentMap, fallbac
         addressModeU: "repeat",
         addressModeV: "clamp-to-edge",
         magFilter: "linear",
-        minFilter: "linear"
+        minFilter: "linear",
+        mipmapFilter: "linear"
       }),
       texture: null,
-      ownsTexture: false
+      ownsTexture: false,
+      width: Math.max(1, environmentMap.width),
+      height: Math.max(1, environmentMap.height),
+      mipLevelCount: Math.max(1, environmentMap.mipLevelCount ?? 1)
     });
   }
   if (environmentMap.texture && typeof environmentMap.texture.createView === "function") {
@@ -1819,15 +3006,91 @@ function createEnvironmentMapResource(device, constants, environmentMap, fallbac
         addressModeU: "repeat",
         addressModeV: "clamp-to-edge",
         magFilter: "linear",
-        minFilter: "linear"
+        minFilter: "linear",
+        mipmapFilter: "linear"
       }),
       texture: environmentMap.texture,
-      ownsTexture: false
+      ownsTexture: false,
+      width: Math.max(1, environmentMap.width),
+      height: Math.max(1, environmentMap.height),
+      mipLevelCount: Math.max(1, environmentMap.mipLevelCount ?? 1)
     });
   }
-  const upload = createEnvironmentMapUploadBytes(environmentMap, fallbackColor);
+  const prefiltered = createPrefilteredEnvironmentLevels(environmentMap, fallbackColor);
+  const uploads = createFloat16RgbaUploadFromLevels(prefiltered.levels);
   const texture = device.createTexture({
     label: environmentMap.enabled ? "plasius.wavefront.environmentMap" : "plasius.wavefront.environmentMapFallback",
+    size: { width: prefiltered.width, height: prefiltered.height },
+    format: "rgba16float",
+    mipLevelCount: prefiltered.mipLevelCount,
+    usage: constants.texture.TEXTURE_BINDING | constants.texture.COPY_DST
+  });
+  uploads.forEach((upload, mipLevel) => {
+    device.queue.writeTexture(
+      { texture, mipLevel },
+      upload.bytes,
+      { bytesPerRow: upload.bytesPerRow, rowsPerImage: upload.height },
+      { width: upload.width, height: upload.height, depthOrArrayLayers: 1 }
+    );
+  });
+  return Object.freeze({
+    view: texture.createView(),
+    sampler: environmentMap.sampler ?? device.createSampler({
+      label: "plasius.wavefront.environmentMapSampler",
+      addressModeU: "repeat",
+      addressModeV: "clamp-to-edge",
+      magFilter: "linear",
+      minFilter: "linear",
+      mipmapFilter: "linear"
+    }),
+    texture,
+    ownsTexture: true,
+    width: prefiltered.width,
+    height: prefiltered.height,
+    mipLevelCount: prefiltered.mipLevelCount
+  });
+}
+function createEnvironmentSamplingTextureResource(device, constants, environmentMap, fallbackColor) {
+  const tables = createEnvironmentSamplingTables(environmentMap, fallbackColor);
+  const rowBytes = tables.width * 8;
+  const bytesPerRow = alignTo(rowBytes, 256);
+  const bytes = new Uint8Array(bytesPerRow * tables.height);
+  const view = new DataView(bytes.buffer);
+  for (let y = 0; y < tables.height; y += 1) {
+    for (let x = 0; x < tables.width; x += 1) {
+      const probability = tables.pdf[y * tables.width + x];
+      const conditional = tables.conditionalCdf[y * tables.width + x];
+      const marginal = tables.marginalCdf[y];
+      const offset = y * bytesPerRow + x * 8;
+      view.setUint16(offset, float32ToFloat16Bits(probability), true);
+      view.setUint16(offset + 2, float32ToFloat16Bits(conditional), true);
+      view.setUint16(offset + 4, float32ToFloat16Bits(marginal), true);
+      view.setUint16(offset + 6, float32ToFloat16Bits(1), true);
+    }
+  }
+  const texture = device.createTexture({
+    label: "plasius.wavefront.environmentSampling",
+    size: { width: tables.width, height: tables.height },
+    format: "rgba16float",
+    usage: constants.texture.TEXTURE_BINDING | constants.texture.COPY_DST
+  });
+  device.queue.writeTexture(
+    { texture },
+    bytes,
+    { bytesPerRow, rowsPerImage: tables.height },
+    { width: tables.width, height: tables.height, depthOrArrayLayers: 1 }
+  );
+  return Object.freeze({
+    view: texture.createView(),
+    texture,
+    ownsTexture: true,
+    hasImportanceData: tables.hasImportanceData
+  });
+}
+function createBrdfLutResource(device, constants, size = DEFAULT_BRDF_LUT_SIZE) {
+  const upload = createBrdfLutUploadBytes(size);
+  const texture = device.createTexture({
+    label: "plasius.wavefront.brdfLut",
     size: { width: upload.width, height: upload.height },
     format: "rgba16float",
     usage: constants.texture.TEXTURE_BINDING | constants.texture.COPY_DST
@@ -1840,14 +3103,110 @@ function createEnvironmentMapResource(device, constants, environmentMap, fallbac
   );
   return Object.freeze({
     view: texture.createView(),
-    sampler: environmentMap.sampler ?? device.createSampler({
-      label: "plasius.wavefront.environmentMapSampler",
-      addressModeU: "repeat",
+    sampler: device.createSampler({
+      label: "plasius.wavefront.brdfLutSampler",
+      addressModeU: "clamp-to-edge",
       addressModeV: "clamp-to-edge",
       magFilter: "linear",
       minFilter: "linear"
     }),
     texture,
+    ownsTexture: true,
+    width: upload.width,
+    height: upload.height
+  });
+}
+function createMediumTextureResource(device, constants, mediums) {
+  const normalized = Array.isArray(mediums) && mediums.length > 0 ? mediums : [{ id: 0 }];
+  const width = Math.max(
+    1,
+    normalized.reduce((maximum, medium) => Math.max(maximum, medium.id ?? 0), 0) + 1
+  );
+  const level = {
+    width,
+    height: MEDIUM_TABLE_ROWS,
+    data: new Float32Array(width * MEDIUM_TABLE_ROWS * 4)
+  };
+  for (const medium of normalized) {
+    const mediumId = Math.max(0, Math.trunc(Number(medium.id) || 0));
+    const absorptionOffset = mediumId * 4;
+    level.data[absorptionOffset] = Math.max(0, medium.absorption?.[0] ?? 0);
+    level.data[absorptionOffset + 1] = Math.max(0, medium.absorption?.[1] ?? 0);
+    level.data[absorptionOffset + 2] = Math.max(0, medium.absorption?.[2] ?? 0);
+    level.data[absorptionOffset + 3] = Math.max(0, medium.phaseModel ?? 0);
+    const scatteringOffset = (width + mediumId) * 4;
+    level.data[scatteringOffset] = Math.max(0, medium.scattering?.[0] ?? 0);
+    level.data[scatteringOffset + 1] = Math.max(0, medium.scattering?.[1] ?? 0);
+    level.data[scatteringOffset + 2] = Math.max(0, medium.scattering?.[2] ?? 0);
+    level.data[scatteringOffset + 3] = Math.max(0, medium.density ?? 0);
+  }
+  const upload = createFloat16RgbaUploadFromLevels([level])[0];
+  const texture = device.createTexture({
+    label: "plasius.wavefront.mediumTable",
+    size: { width, height: MEDIUM_TABLE_ROWS },
+    format: "rgba16float",
+    usage: constants.texture.TEXTURE_BINDING | constants.texture.COPY_DST
+  });
+  device.queue.writeTexture(
+    { texture },
+    upload.bytes,
+    { bytesPerRow: upload.bytesPerRow, rowsPerImage: upload.height },
+    { width, height: MEDIUM_TABLE_ROWS, depthOrArrayLayers: 1 }
+  );
+  return Object.freeze({
+    texture,
+    view: texture.createView(),
+    ownsTexture: true,
+    count: normalized.length,
+    width
+  });
+}
+function mediumTablesEqual(left, right) {
+  const leftMediums = Array.isArray(left) ? left : [];
+  const rightMediums = Array.isArray(right) ? right : [];
+  if (leftMediums.length !== rightMediums.length) {
+    return false;
+  }
+  for (let index = 0; index < leftMediums.length; index += 1) {
+    const leftMedium = leftMediums[index];
+    const rightMedium = rightMediums[index];
+    if ((leftMedium?.id ?? 0) !== (rightMedium?.id ?? 0)) {
+      return false;
+    }
+    if ((leftMedium?.phaseModel ?? 0) !== (rightMedium?.phaseModel ?? 0)) {
+      return false;
+    }
+    if ((leftMedium?.density ?? 0) !== (rightMedium?.density ?? 0)) {
+      return false;
+    }
+    for (let component = 0; component < 3; component += 1) {
+      if ((leftMedium?.absorption?.[component] ?? 0) !== (rightMedium?.absorption?.[component] ?? 0)) {
+        return false;
+      }
+      if ((leftMedium?.scattering?.[component] ?? 0) !== (rightMedium?.scattering?.[component] ?? 0)) {
+        return false;
+      }
+    }
+  }
+  return true;
+}
+function createAtlasTextureResource(device, constants, atlas, label) {
+  const upload = createRgba8TextureUpload(atlas);
+  const texture = device.createTexture({
+    label,
+    size: { width: upload.width, height: upload.height },
+    format: "rgba8unorm",
+    usage: constants.texture.TEXTURE_BINDING | constants.texture.COPY_DST
+  });
+  device.queue.writeTexture(
+    { texture },
+    upload.bytes,
+    { bytesPerRow: upload.bytesPerRow, rowsPerImage: upload.height },
+    { width: upload.width, height: upload.height, depthOrArrayLayers: 1 }
+  );
+  return Object.freeze({
+    texture,
+    view: texture.createView(),
     ownsTexture: true
   });
 }
@@ -1893,6 +3252,24 @@ ${diagnostics}` : "";
     });
   }
 }
+async function assertShaderModuleCompiles(shaderModule, label) {
+  if (typeof shaderModule?.compilationInfo !== "function") {
+    return;
+  }
+  const info = await shaderModule.compilationInfo();
+  const messages = Array.isArray(info?.messages) ? info.messages : [];
+  const errors = messages.filter((message) => message?.type === "error");
+  if (errors.length <= 0) {
+    return;
+  }
+  const diagnostics = errors.map((message) => {
+    const line = Number.isFinite(message.lineNum) ? message.lineNum : "?";
+    const column = Number.isFinite(message.linePos) ? message.linePos : "?";
+    return `line ${line}:${column} ${message.message}`;
+  }).join("\n");
+  throw new Error(`WGSL compilation preflight failed for ${label}:
+${diagnostics}`);
+}
 async function createRenderPipeline(device, descriptor) {
   if (typeof device.createRenderPipelineAsync === "function") {
     return device.createRenderPipelineAsync(descriptor);
@@ -1901,6 +3278,7 @@ async function createRenderPipeline(device, descriptor) {
 }
 var WAVEFRONT_COMPUTE_WGSL = `
 const RAY_FLAG_GUIDED_EMISSIVE: u32 = 1u;
+const RAY_FLAG_DELTA_SAMPLE: u32 = 2u;
 
 struct RayRecord {
   rayId: u32,
@@ -1926,11 +3304,12 @@ struct HitRecord {
   primitiveId: u32,
   materialRefId: u32,
   mediumRefId: u32,
+  materialSlot: u32,
   pad0: u32,
   pad1: u32,
-  pad2: u32,
   distance: f32,
-  pad3: vec3<f32>,
+  occlusion: f32,
+  pad2: vec2<f32>,
   position: vec4<f32>,
   geometricNormal: vec4<f32>,
   shadingNormal: vec4<f32>,
@@ -1939,6 +3318,9 @@ struct HitRecord {
   color: vec4<f32>,
   emission: vec4<f32>,
   material: vec4<f32>,
+  materialResponse: vec4<f32>,
+  materialExtension: vec4<f32>,
+  specularColor: vec4<f32>,
 };
 
 struct SceneObject {
@@ -1946,11 +3328,18 @@ struct SceneObject {
   objectId: u32,
   materialKind: u32,
   flags: u32,
+  mediumRefId: u32,
+  pad0: u32,
+  pad1: u32,
+  pad2: u32,
   center: vec4<f32>,
   halfExtent: vec4<f32>,
   color: vec4<f32>,
   emission: vec4<f32>,
   material: vec4<f32>,
+  materialResponse: vec4<f32>,
+  materialExtension: vec4<f32>,
+  specularColor: vec4<f32>,
 };
 
 struct TriangleRecord {
@@ -1960,7 +3349,7 @@ struct TriangleRecord {
   flags: u32,
   materialRefId: u32,
   mediumRefId: u32,
-  pad0: u32,
+  materialSlot: u32,
   pad1: u32,
   v0: vec4<f32>,
   v1: vec4<f32>,
@@ -1973,6 +3362,15 @@ struct TriangleRecord {
   color: vec4<f32>,
   emission: vec4<f32>,
   material: vec4<f32>,
+  materialResponse: vec4<f32>,
+  materialExtension: vec4<f32>,
+  specularColor: vec4<f32>,
+  baseColorAtlas: vec4<f32>,
+  metallicRoughnessAtlas: vec4<f32>,
+  normalAtlas: vec4<f32>,
+  occlusionAtlas: vec4<f32>,
+  emissiveAtlas: vec4<f32>,
+  textureSettings: vec4<f32>,
 };
 
 struct BvhNode {
@@ -1993,10 +3391,10 @@ struct BvhLeafRef {
 
 struct ScatterResult {
   direction: vec4<f32>,
+  pdf: f32,
+  mediumRefId: u32,
   flags: u32,
   pad0: u32,
-  pad1: u32,
-  pad2: u32,
 };
 
 struct MeshVertex {
@@ -2017,10 +3415,19 @@ struct MeshRange {
   triangleCount: u32,
   firstVertex: u32,
   vertexCount: u32,
-  pad0: u32,
+  materialSlot: u32,
   color: vec4<f32>,
   emission: vec4<f32>,
   material: vec4<f32>,
+  materialResponse: vec4<f32>,
+  materialExtension: vec4<f32>,
+  specularColor: vec4<f32>,
+  baseColorAtlas: vec4<f32>,
+  metallicRoughnessAtlas: vec4<f32>,
+  normalAtlas: vec4<f32>,
+  occlusionAtlas: vec4<f32>,
+  emissiveAtlas: vec4<f32>,
+  textureSettings: vec4<f32>,
 };
 
 struct FrameConfig {
@@ -2061,6 +3468,7 @@ struct FrameConfig {
   _portalPad1: u32,
   environmentMapSettings: vec4<f32>,
   pathResolveSettings: vec4<f32>,
+  environmentMapMeta: vec4<f32>,
 };
 
 struct Counters {
@@ -2123,6 +3531,16 @@ struct EnvironmentPortal {
 @group(0) @binding(20) var environmentMapTexture: texture_2d<f32>;
 @group(0) @binding(21) var environmentMapSampler: sampler;
 @group(0) @binding(22) var<storage, read_write> pathVertices: array<vec4<f32>>;
+@group(0) @binding(23) var baseColorAtlasTexture: texture_2d<f32>;
+@group(0) @binding(24) var metallicRoughnessAtlasTexture: texture_2d<f32>;
+@group(0) @binding(25) var normalAtlasTexture: texture_2d<f32>;
+@group(0) @binding(26) var occlusionAtlasTexture: texture_2d<f32>;
+@group(0) @binding(27) var emissiveAtlasTexture: texture_2d<f32>;
+@group(0) @binding(28) var materialAtlasSampler: sampler;
+@group(0) @binding(29) var brdfLutTexture: texture_2d<f32>;
+@group(0) @binding(30) var brdfLutSampler: sampler;
+@group(0) @binding(31) var environmentSamplingTexture: texture_2d<f32>;
+@group(0) @binding(32) var mediumTableTexture: texture_2d<f32>;
 
 fn hash_u32(value: u32) -> u32 {
   var x = value;
@@ -2159,6 +3577,146 @@ fn safe_normalize(value: vec3<f32>, fallback: vec3<f32>) -> vec3<f32> {
   return value / len;
 }
 
+struct TangentBasis {
+  tangent: vec3<f32>,
+  bitangent: vec3<f32>,
+};
+
+struct SurfaceMaterialSample {
+  color: vec4<f32>,
+  emission: vec4<f32>,
+  material: vec4<f32>,
+  materialResponse: vec4<f32>,
+  materialExtension: vec4<f32>,
+  specularColor: vec4<f32>,
+  shadingNormal: vec3<f32>,
+  occlusion: f32,
+};
+
+fn srgb_to_linear_channel(value: f32) -> f32 {
+  if (value <= 0.04045) {
+    return value / 12.92;
+  }
+  return pow((value + 0.055) / 1.055, 2.4);
+}
+
+fn srgb_to_linear_vec3(value: vec3<f32>) -> vec3<f32> {
+  return vec3<f32>(
+    srgb_to_linear_channel(value.x),
+    srgb_to_linear_channel(value.y),
+    srgb_to_linear_channel(value.z)
+  );
+}
+
+fn wrap_uv(uv: vec2<f32>) -> vec2<f32> {
+  return fract(fract(uv) + vec2<f32>(1.0));
+}
+
+fn atlas_sample_uv(rect: vec4<f32>, uv: vec2<f32>) -> vec2<f32> {
+  let local = wrap_uv(uv);
+  let clamped = clamp(local, vec2<f32>(0.001), vec2<f32>(0.999));
+  return rect.xy + clamped * rect.zw;
+}
+
+fn sample_atlas(textureRef: texture_2d<f32>, rect: vec4<f32>, uv: vec2<f32>) -> vec4<f32> {
+  return textureSampleLevel(textureRef, materialAtlasSampler, atlas_sample_uv(rect, uv), 0.0);
+}
+
+fn build_triangle_tangent_basis(
+  triangle: TriangleRecord,
+  fallbackNormal: vec3<f32>
+) -> TangentBasis {
+  let edge1 = triangle.v1.xyz - triangle.v0.xyz;
+  let edge2 = triangle.v2.xyz - triangle.v0.xyz;
+  let uv0 = triangle.uv0uv1.xy;
+  let uv1 = triangle.uv0uv1.zw;
+  let uv2 = triangle.uv2Pad.xy;
+  let deltaUv1 = uv1 - uv0;
+  let deltaUv2 = uv2 - uv0;
+  let determinant = deltaUv1.x * deltaUv2.y - deltaUv1.y * deltaUv2.x;
+  if (abs(determinant) <= 0.000001) {
+    let tangentFallback = select(vec3<f32>(0.0, 1.0, 0.0), vec3<f32>(1.0, 0.0, 0.0), abs(fallbackNormal.y) >= 0.999);
+    let tangent = safe_normalize(cross(tangentFallback, fallbackNormal), vec3<f32>(1.0, 0.0, 0.0));
+    let bitangent = safe_normalize(cross(fallbackNormal, tangent), vec3<f32>(0.0, 0.0, 1.0));
+    return TangentBasis(tangent, bitangent);
+  }
+  let inverse = 1.0 / determinant;
+  let tangent = safe_normalize(
+    inverse * (edge1 * deltaUv2.y - edge2 * deltaUv1.y),
+    vec3<f32>(1.0, 0.0, 0.0)
+  );
+  let bitangent = safe_normalize(
+    inverse * (-edge1 * deltaUv2.x + edge2 * deltaUv1.x),
+    vec3<f32>(0.0, 0.0, 1.0)
+  );
+  return TangentBasis(tangent, bitangent);
+}
+
+fn sample_surface_material(
+  triangle: TriangleRecord,
+  uv: vec2<f32>,
+  geometricNormal: vec3<f32>,
+  shadingNormal: vec3<f32>
+) -> SurfaceMaterialSample {
+  let baseColorTexel = sample_atlas(baseColorAtlasTexture, triangle.baseColorAtlas, uv);
+  let baseColor = vec4<f32>(
+    clamp(triangle.color.rgb * srgb_to_linear_vec3(baseColorTexel.rgb), vec3<f32>(0.0), vec3<f32>(1.0)),
+    clamp(triangle.color.a * baseColorTexel.a, 0.0, 1.0)
+  );
+  let metallicRoughnessTexel = sample_atlas(
+    metallicRoughnessAtlasTexture,
+    triangle.metallicRoughnessAtlas,
+    uv
+  );
+  let normalTexel = sample_atlas(normalAtlasTexture, triangle.normalAtlas, uv);
+  let occlusionTexel = sample_atlas(occlusionAtlasTexture, triangle.occlusionAtlas, uv);
+  let emissiveTexel = sample_atlas(emissiveAtlasTexture, triangle.emissiveAtlas, uv);
+  let normalScale = clamp(triangle.textureSettings.x, 0.0, 1.0);
+  let tangentBasis = build_triangle_tangent_basis(triangle, geometricNormal);
+  let tangentNormal = safe_normalize(
+    vec3<f32>(
+      (normalTexel.x * 2.0 - 1.0) * normalScale,
+      (normalTexel.y * 2.0 - 1.0) * normalScale,
+      1.0 + ((normalTexel.z * 2.0 - 1.0) - 1.0) * normalScale
+    ),
+    vec3<f32>(0.0, 0.0, 1.0)
+  );
+  let mappedNormal = safe_normalize(
+    tangentBasis.tangent * tangentNormal.x +
+      tangentBasis.bitangent * tangentNormal.y +
+      shadingNormal * tangentNormal.z,
+    shadingNormal
+  );
+  let emission = vec4<f32>(
+    max(
+      triangle.emission.rgb *
+        srgb_to_linear_vec3(emissiveTexel.rgb) *
+        max(triangle.textureSettings.z, 0.0),
+      vec3<f32>(0.0)
+    ),
+    clamp(triangle.emission.a * emissiveTexel.a, 0.0, 1.0)
+  );
+  return SurfaceMaterialSample(
+    baseColor,
+    emission,
+    vec4<f32>(
+      clamp(triangle.material.x * metallicRoughnessTexel.y, 0.0, 1.0),
+      clamp(triangle.material.y * metallicRoughnessTexel.z, 0.0, 1.0),
+      clamp(triangle.material.z * baseColor.a, 0.0, 1.0),
+      clamp(triangle.material.w, 1.0, 3.0)
+    ),
+    triangle.materialResponse,
+    triangle.materialExtension,
+    triangle.specularColor,
+    repair_shading_normal(geometricNormal, mappedNormal),
+    clamp(
+      mix(1.0, occlusionTexel.x, clamp(triangle.textureSettings.y, 0.0, 1.0)),
+      0.0,
+      1.0
+    )
+  );
+}
+
 fn saturate(value: f32) -> f32 {
   return clamp(value, 0.0, 1.0);
 }
@@ -2167,6 +3725,10 @@ fn max_component(value: vec3<f32>) -> f32 {
   return max(max(value.x, value.y), value.z);
 }
 
+fn radiance_luminance(value: vec3<f32>) -> f32 {
+  return dot(value, vec3<f32>(0.2126, 0.7152, 0.0722));
+}
+
 fn environment_map_enabled() -> bool {
   return config.environmentMapSettings.x > 0.5;
 }
@@ -2287,12 +3849,349 @@ fn environment_portal_radiance_scale(origin: vec3<f32>, direction: vec3<f32>) ->
   return scale;
 }
 
-fn environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f32> {
-  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
-  let portalScale = environment_portal_radiance_scale(origin, rayDirection);
-  let portalHit = max_component(portalScale) > 0.0001;
-  return base_environment_radiance(rayDirection) *
-    select(vec3<f32>(1.0), portalScale, portalHit);
+fn environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f32> {
+  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
+  let portalScale = environment_portal_radiance_scale(origin, rayDirection);
+  let portalHit = max_component(portalScale) > 0.0001;
+  return base_environment_radiance(rayDirection) *
+    select(vec3<f32>(1.0), portalScale, portalHit);
+}
+
+fn direct_environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f32> {
+  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
+  let portalScale = environment_portal_radiance_scale(origin, rayDirection);
+  let portalHit = max_component(portalScale) > 0.0001;
+  if (
+    config.environmentPortalCount > 0u &&
+    config.environmentPortalMode == 2u &&
+    !portalHit
+  ) {
+    return vec3<f32>(0.0);
+  }
+  return base_environment_radiance(rayDirection) *
+    select(vec3<f32>(1.0), portalScale, portalHit);
+}
+
+fn radical_inverse_vdc(bitsValue: u32) -> f32 {
+  var bits = bitsValue;
+  bits = (bits << 16u) | (bits >> 16u);
+  bits = ((bits & 0x55555555u) << 1u) | ((bits & 0xaaaaaaaau) >> 1u);
+  bits = ((bits & 0x33333333u) << 2u) | ((bits & 0xccccccccu) >> 2u);
+  bits = ((bits & 0x0f0f0f0fu) << 4u) | ((bits & 0xf0f0f0f0u) >> 4u);
+  bits = ((bits & 0x00ff00ffu) << 8u) | ((bits & 0xff00ff00u) >> 8u);
+  return f32(bits) * 2.3283064365386963e-10;
+}
+
+fn hammersley_2d(index: u32, count: u32) -> vec2<f32> {
+  return vec2<f32>(f32(index) / max(f32(count), 1.0), radical_inverse_vdc(index));
+}
+
+fn build_basis_tangent(normal: vec3<f32>) -> vec3<f32> {
+  let tangentFallback = select(vec3<f32>(0.0, 1.0, 0.0), vec3<f32>(1.0, 0.0, 0.0), abs(normal.y) >= 0.999);
+  return safe_normalize(cross(tangentFallback, normal), vec3<f32>(1.0, 0.0, 0.0));
+}
+
+fn local_to_world(local: vec3<f32>, normal: vec3<f32>) -> vec3<f32> {
+  let tangent = build_basis_tangent(normal);
+  let bitangent = safe_normalize(cross(normal, tangent), vec3<f32>(0.0, 0.0, 1.0));
+  return safe_normalize(tangent * local.x + bitangent * local.y + normal * local.z, normal);
+}
+
+fn cosine_sample_hemisphere(sample: vec2<f32>, normal: vec3<f32>) -> vec3<f32> {
+  let phi = 6.28318530718 * sample.x;
+  let radius = sqrt(sample.y);
+  let x = cos(phi) * radius;
+  let y = sin(phi) * radius;
+  let z = sqrt(max(0.0, 1.0 - sample.y));
+  return local_to_world(vec3<f32>(x, y, z), normal);
+}
+
+fn importance_sample_ggx(sample: vec2<f32>, roughness: f32, normal: vec3<f32>) -> vec3<f32> {
+  let alpha = max(roughness * roughness, 0.0001);
+  let phi = 6.28318530718 * sample.x;
+  let cosTheta = sqrt((1.0 - sample.y) / max(1.0 + (alpha * alpha - 1.0) * sample.y, 0.0001));
+  let sinTheta = sqrt(max(0.0, 1.0 - cosTheta * cosTheta));
+  let localHalf = vec3<f32>(cos(phi) * sinTheta, sin(phi) * sinTheta, cosTheta);
+  return local_to_world(localHalf, normal);
+}
+
+fn distribution_ggx(normal: vec3<f32>, halfVector: vec3<f32>, roughness: f32) -> f32 {
+  let alpha = max(roughness * roughness, 0.0001);
+  let alpha2 = alpha * alpha;
+  let nDotH = saturate(dot(normal, halfVector));
+  let denominator = nDotH * nDotH * (alpha2 - 1.0) + 1.0;
+  return alpha2 / max(3.14159265359 * denominator * denominator, 0.000001);
+}
+
+fn geometry_schlick_ggx(nDotValue: f32, roughness: f32) -> f32 {
+  let k = ((roughness + 1.0) * (roughness + 1.0)) / 8.0;
+  return nDotValue / max(nDotValue * (1.0 - k) + k, 0.000001);
+}
+
+fn geometry_smith(normal: vec3<f32>, viewDirection: vec3<f32>, lightDirection: vec3<f32>, roughness: f32) -> f32 {
+  let nDotV = saturate(dot(normal, viewDirection));
+  let nDotL = saturate(dot(normal, lightDirection));
+  return geometry_schlick_ggx(nDotV, roughness) * geometry_schlick_ggx(nDotL, roughness);
+}
+
+fn fresnel_schlick(cosine: f32, f0: vec3<f32>) -> vec3<f32> {
+  return f0 + (vec3<f32>(1.0) - f0) * pow(1.0 - cosine, 5.0);
+}
+
+fn sample_brdf_lut(nDotV: f32, roughness: f32) -> vec2<f32> {
+  let uv = vec2<f32>(clamp(nDotV, 0.0, 1.0), clamp(roughness, 0.0, 1.0));
+  return textureSampleLevel(brdfLutTexture, brdfLutSampler, uv, 0.0).xy;
+}
+
+fn prefiltered_environment_radiance(direction: vec3<f32>, roughness: f32) -> vec3<f32> {
+  let uv = environment_map_uv(direction);
+  let maxLevel = max(config.environmentMapMeta.z - 1.0, 0.0);
+  let lod = clamp(roughness, 0.0, 1.0) * maxLevel;
+  let texel = max(textureSampleLevel(environmentMapTexture, environmentMapSampler, uv, lod).rgb, vec3<f32>(0.0));
+  return texel * max(config.environmentMapSettings.y, 0.0);
+}
+
+fn environment_pdf_dimensions() -> vec2<u32> {
+  return vec2<u32>(
+    max(u32(config.environmentMapMeta.x), 1u),
+    max(u32(config.environmentMapMeta.y), 1u)
+  );
+}
+
+fn environment_importance_sampling_enabled() -> bool {
+  return config.environmentMapMeta.w > 0.5;
+}
+
+fn uniform_sphere_pdf() -> f32 {
+  return 1.0 / (4.0 * 3.14159265359);
+}
+
+fn sample_uniform_sphere_direction(sample: vec2<f32>) -> vec3<f32> {
+  let z = 1.0 - 2.0 * sample.y;
+  let radial = sqrt(max(1.0 - z * z, 0.0));
+  let phi = sample.x * 6.28318530718;
+  return vec3<f32>(cos(phi) * radial, z, sin(phi) * radial);
+}
+
+fn environment_sampling_texel(x: u32, y: u32) -> vec4<f32> {
+  return textureLoad(environmentSamplingTexture, vec2<i32>(i32(x), i32(y)), 0);
+}
+
+fn environment_pdf_texel(x: u32, y: u32) -> f32 {
+  return environment_sampling_texel(x, y).x;
+}
+
+fn environment_row_cdf_texel(y: u32) -> f32 {
+  return environment_sampling_texel(0u, y).z;
+}
+
+fn environment_column_cdf_texel(x: u32, y: u32) -> f32 {
+  return environment_sampling_texel(x, y).y;
+}
+
+fn environment_direction_pdf(direction: vec3<f32>) -> f32 {
+  if (!environment_importance_sampling_enabled()) {
+    return uniform_sphere_pdf();
+  }
+  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
+  let uv = environment_map_uv(rayDirection);
+  let dimensions = environment_pdf_dimensions();
+  let width = max(f32(dimensions.x), 1.0);
+  let height = max(f32(dimensions.y), 1.0);
+  let x = min(u32(uv.x * width), dimensions.x - 1u);
+  let y = min(u32(uv.y * height), dimensions.y - 1u);
+  let discretePdf = max(environment_pdf_texel(x, y), 0.0);
+  let sinTheta = sqrt(max(1.0 - rayDirection.y * rayDirection.y, 0.0));
+  let solidAngle = max((2.0 * 3.14159265359 * 3.14159265359 * sinTheta) / (width * height), 0.000001);
+  return discretePdf / solidAngle;
+}
+
+fn sample_row_cdf(count: u32, sampleValue: f32) -> u32 {
+  if (count == 0u) {
+    return 0u;
+  }
+  var low = 0u;
+  var high = count - 1u;
+  loop {
+    if (low >= high) {
+      break;
+    }
+    let mid = (low + high) / 2u;
+    let cdfValue = environment_row_cdf_texel(mid);
+    if (sampleValue <= cdfValue) {
+      high = mid;
+    } else {
+      low = mid + 1u;
+    }
+  }
+  return min(low, count - 1u);
+}
+
+fn sample_column_cdf(row: u32, count: u32, sampleValue: f32) -> u32 {
+  if (count == 0u) {
+    return 0u;
+  }
+  var low = 0u;
+  var high = count - 1u;
+  loop {
+    if (low >= high) {
+      break;
+    }
+    let mid = (low + high) / 2u;
+    let cdfValue = environment_column_cdf_texel(mid, row);
+    if (sampleValue <= cdfValue) {
+      high = mid;
+    } else {
+      low = mid + 1u;
+    }
+  }
+  return min(low, count - 1u);
+}
+
+struct EnvironmentSample {
+  direction: vec3<f32>,
+  radiance: vec3<f32>,
+  pdf: f32,
+};
+
+fn sample_environment_importance(sample: vec2<f32>) -> EnvironmentSample {
+  if (!environment_importance_sampling_enabled()) {
+    let direction = sample_uniform_sphere_direction(sample);
+    return EnvironmentSample(direction, base_environment_radiance(direction), uniform_sphere_pdf());
+  }
+  let dimensions = environment_pdf_dimensions();
+  let row = sample_row_cdf(dimensions.y, sample.y);
+  let column = sample_column_cdf(row, dimensions.x, sample.x);
+  let uv = vec2<f32>(
+    (f32(column) + 0.5) / max(f32(dimensions.x), 1.0),
+    (f32(row) + 0.5) / max(f32(dimensions.y), 1.0)
+  );
+  let theta = uv.y * 3.14159265359;
+  let phi = (uv.x - 0.5 - config.environmentMapSettings.z / 6.28318530718) * 6.28318530718;
+  let sinTheta = sin(theta);
+  let direction = vec3<f32>(cos(phi) * sinTheta, cos(theta), sin(phi) * sinTheta);
+  let pdf = environment_direction_pdf(direction);
+  return EnvironmentSample(direction, base_environment_radiance(direction), pdf);
+}
+
+fn power_heuristic(pdfA: f32, pdfB: f32) -> f32 {
+  let a2 = pdfA * pdfA;
+  let b2 = pdfB * pdfB;
+  return a2 / max(a2 + b2, 0.000001);
+}
+
+fn visible_environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f32> {
+  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
+  let visible = !scene_visibility_blocked(origin, rayDirection, 1000000.0);
+  return select(vec3<f32>(0.0), direct_environment_radiance(origin, rayDirection), visible);
+}
+
+fn glossy_environment_direction(
+  incidentDirection: vec3<f32>,
+  normal: vec3<f32>,
+  roughness: f32,
+  normalBlendScale: f32
+) -> vec3<f32> {
+  let reflectionDirection = reflect(incidentDirection, normal);
+  let blend = clamp(roughness * roughness * normalBlendScale, 0.0, 0.92);
+  return safe_normalize(mix(reflectionDirection, normal, blend), normal);
+}
+
+fn surface_glossiness(hit: HitRecord) -> f32 {
+  let roughness = clamp(hit.material.x, 0.0, 1.0);
+  let metallic = clamp(hit.material.y, 0.0, 1.0);
+  let sheen = clamp(max_component(hit.materialResponse.xyz), 0.0, 1.0);
+  let clearcoat = clamp(hit.materialResponse.w, 0.0, 1.0);
+  let specularWeight = clamp(hit.materialExtension.y, 0.0, 1.0);
+  let transmission = clamp(hit.materialExtension.z, 0.0, 1.0);
+  let baseGloss =
+    max(
+      clearcoat,
+      max(sheen * 0.72, max(specularWeight * (0.38 + metallic * 0.62), transmission))
+    );
+  return clamp(baseGloss * (1.0 - roughness * 0.72) + metallic * (1.0 - roughness) * 0.35, 0.0, 1.0);
+}
+
+fn surface_specular_f0(hit: HitRecord, surfaceColor: vec3<f32>) -> vec3<f32> {
+  let metallic = clamp(hit.material.y, 0.0, 1.0);
+  let specularWeight = clamp(hit.materialExtension.y, 0.0, 1.0);
+  let specularColor = clamp(hit.specularColor.xyz, vec3<f32>(0.0), vec3<f32>(1.0));
+  let dielectricF0 = vec3<f32>(0.04) * specularWeight * specularColor;
+  return mix(dielectricF0, surfaceColor, metallic);
+}
+
+fn surface_bsdf_sampling_weights(hit: HitRecord) -> vec3<f32> {
+  let metallic = clamp(hit.material.y, 0.0, 1.0);
+  let clearcoat = clamp(hit.materialResponse.w, 0.0, 1.0);
+  let specularWeight = clamp(hit.materialExtension.y, 0.0, 1.0);
+  let diffuseWeight = clamp(
+    (1.0 - metallic) * max(1.0 - specularWeight * 0.5 - clearcoat * 0.25, 0.15),
+    0.0,
+    1.0
+  );
+  let specWeight = clamp(max(metallic, specularWeight * 0.75) * (1.0 - clearcoat * 0.5), 0.0, 1.0);
+  let clearcoatWeight = clamp(clearcoat, 0.0, 1.0);
+  let totalWeight = max(diffuseWeight + specWeight + clearcoatWeight, 0.000001);
+  return vec3<f32>(
+    diffuseWeight / totalWeight,
+    specWeight / totalWeight,
+    clearcoatWeight / totalWeight
+  );
+}
+
+fn evaluate_surface_bsdf(hit: HitRecord, viewDirection: vec3<f32>, lightDirection: vec3<f32>) -> vec3<f32> {
+  let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
+  let surfaceColor = clamp(max(hit.color.xyz, config.ambientColor.xyz * 0.35), vec3<f32>(0.0), vec3<f32>(1.0));
+  let roughness = clamp(hit.material.x, 0.0, 1.0);
+  let metallic = clamp(hit.material.y, 0.0, 1.0);
+  let clearcoat = clamp(hit.materialResponse.w, 0.0, 1.0);
+  let clearcoatRoughness = clamp(hit.materialExtension.x, 0.0, 1.0);
+  let occlusion = clamp(hit.occlusion, 0.0, 1.0);
+  let nDotV = saturate(dot(normal, viewDirection));
+  let nDotL = saturate(dot(normal, lightDirection));
+  if (nDotV <= 0.0 || nDotL <= 0.0) {
+    return vec3<f32>(0.0);
+  }
+  let halfVector = safe_normalize(viewDirection + lightDirection, normal);
+  let vDotH = saturate(dot(viewDirection, halfVector));
+  let f0 = surface_specular_f0(hit, surfaceColor);
+  let fresnel = fresnel_schlick(vDotH, f0);
+  let distribution = distribution_ggx(normal, halfVector, roughness);
+  let geometry = geometry_smith(normal, viewDirection, lightDirection, roughness);
+  let specular = (distribution * geometry * fresnel) / max(4.0 * nDotV * nDotL, 0.000001);
+  let diffuseWeight = (1.0 - metallic) * (1.0 - clearcoat * 0.24) * (1.0 - clamp(max_component(fresnel), 0.0, 0.98));
+  let diffuse = surfaceColor * diffuseWeight / 3.14159265359;
+  let clearcoatHalf = safe_normalize(viewDirection + lightDirection, normal);
+  let clearcoatDistribution = distribution_ggx(normal, clearcoatHalf, max(clearcoatRoughness, 0.02));
+  let clearcoatGeometry = geometry_smith(normal, viewDirection, lightDirection, max(clearcoatRoughness, 0.02));
+  let clearcoatFresnel = fresnel_schlick(saturate(dot(viewDirection, clearcoatHalf)), vec3<f32>(0.04));
+  let clearcoatTerm =
+    (clearcoatDistribution * clearcoatGeometry * clearcoatFresnel) /
+    max(4.0 * nDotV * nDotL, 0.000001) *
+    clearcoat;
+  return (diffuse + specular + clearcoatTerm) * mix(0.42, 1.0, occlusion);
+}
+
+fn diffuse_pdf(normal: vec3<f32>, lightDirection: vec3<f32>) -> f32 {
+  return saturate(dot(normal, lightDirection)) / 3.14159265359;
+}
+
+fn ggx_pdf(normal: vec3<f32>, viewDirection: vec3<f32>, lightDirection: vec3<f32>, roughness: f32) -> f32 {
+  let halfVector = safe_normalize(viewDirection + lightDirection, normal);
+  let nDotH = saturate(dot(normal, halfVector));
+  let vDotH = saturate(dot(viewDirection, halfVector));
+  let distribution = distribution_ggx(normal, halfVector, roughness);
+  return (distribution * nDotH) / max(4.0 * vDotH, 0.000001);
+}
+
+fn evaluate_surface_bsdf_pdf(hit: HitRecord, viewDirection: vec3<f32>, lightDirection: vec3<f32>) -> f32 {
+  let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
+  let roughness = clamp(hit.material.x, 0.0, 1.0);
+  let weights = surface_bsdf_sampling_weights(hit);
+  let diffuseTerm = diffuse_pdf(normal, lightDirection);
+  let specTerm = ggx_pdf(normal, viewDirection, lightDirection, max(roughness, 0.02));
+  let clearcoatTerm = ggx_pdf(normal, viewDirection, lightDirection, max(clamp(hit.materialExtension.x, 0.0, 1.0), 0.02));
+  return weights.x * diffuseTerm + weights.y * specTerm + weights.z * clearcoatTerm;
 }
 
 fn gated_environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f32> {
@@ -2307,12 +4206,102 @@ fn gated_environment_radiance(origin: vec3<f32>, direction: vec3<f32>) -> vec3<f
   return environment_radiance(origin, direction);
 }
 
+fn medium_dimensions() -> vec2<u32> {
+  return textureDimensions(mediumTableTexture);
+}
+
+fn medium_valid(mediumRefId: u32) -> bool {
+  let dimensions = medium_dimensions();
+  return mediumRefId > 0u && mediumRefId < dimensions.x;
+}
+
+fn medium_absorption(mediumRefId: u32) -> vec3<f32> {
+  if (!medium_valid(mediumRefId)) {
+    return vec3<f32>(0.0);
+  }
+  return max(
+    textureLoad(mediumTableTexture, vec2<i32>(i32(mediumRefId), 0), 0).xyz,
+    vec3<f32>(0.0)
+  );
+}
+
+fn medium_scattering(mediumRefId: u32) -> vec3<f32> {
+  if (!medium_valid(mediumRefId)) {
+    return vec3<f32>(0.0);
+  }
+  return max(
+    textureLoad(mediumTableTexture, vec2<i32>(i32(mediumRefId), 1), 0).xyz,
+    vec3<f32>(0.0)
+  );
+}
+
+fn medium_transmittance(mediumRefId: u32, distance: f32) -> vec3<f32> {
+  if (!medium_valid(mediumRefId) || distance <= 0.000001) {
+    return vec3<f32>(1.0);
+  }
+  let extinction = medium_absorption(mediumRefId) + medium_scattering(mediumRefId);
+  return vec3<f32>(
+    exp(-extinction.x * distance),
+    exp(-extinction.y * distance),
+    exp(-extinction.z * distance)
+  );
+}
+
+fn transmitted_medium_ref_id(ray: RayRecord, hit: HitRecord) -> u32 {
+  if (hit.mediumRefId == 0u) {
+    return ray.mediumRefId;
+  }
+  if (hit.frontFace == 1u) {
+    return hit.mediumRefId;
+  }
+  if (ray.mediumRefId == hit.mediumRefId) {
+    return 0u;
+  }
+  return ray.mediumRefId;
+}
+
 fn surface_path_response(hit: HitRecord) -> vec3<f32> {
   let color = clamp(hit.color.xyz, vec3<f32>(0.0), vec3<f32>(1.0));
   let opacity = clamp(hit.material.z, 0.0, 1.0);
+  let occlusion = clamp(hit.occlusion, 0.0, 1.0);
   let materialEnergy = select(0.68, 0.92, hit.materialKind == 1u || hit.materialKind == 2u);
   let transparentEnergy = select(materialEnergy, 0.9, hit.hitType == 3u);
-  return mix(vec3<f32>(1.0), color, max(opacity, 0.18)) * transparentEnergy;
+  return mix(vec3<f32>(1.0), color, max(opacity, 0.18)) * transparentEnergy * mix(0.55, 1.0, occlusion);
+}
+
+fn bounded_path_response_luminance(ray: RayRecord, hit: HitRecord) -> f32 {
+  let daylightFloor = max(config.pathResolveSettings.y, 0.0) * 0.08;
+  let hdriFloor = max(config.environmentMapSettings.w, 0.0) * 0.02;
+  let sceneFloor = max(daylightFloor, hdriFloor);
+  if (sceneFloor <= 0.000001) {
+    return 0.0;
+  }
+  let bounceRatio = select(
+    0.0,
+    f32(ray.bounce) / max(f32(config.maxDepth - 1u), 1.0),
+    config.maxDepth > 1u
+  );
+  let bounceScale = 1.0 - bounceRatio * 0.55;
+  let materialScale = select(1.0, 0.34, hit.materialKind == 1u || hit.materialKind == 2u);
+  let transparentScale = select(materialScale, 0.58, hit.hitType == 3u);
+  let opacityScale = mix(0.55, 1.0, clamp(hit.material.z, 0.0, 1.0));
+  return sceneFloor * bounceScale * transparentScale * opacityScale;
+}
+
+fn stabilize_surface_path_response(ray: RayRecord, hit: HitRecord, response: vec3<f32>) -> vec3<f32> {
+  let minimumLuminance = bounded_path_response_luminance(ray, hit);
+  let responseLuminance = radiance_luminance(response);
+  if (minimumLuminance <= 0.000001 || responseLuminance >= minimumLuminance) {
+    return response;
+  }
+  let tintBase = max(response, max(hit.color.xyz * 0.65, config.ambientColor.xyz * 0.35));
+  let tint = tintBase / max(max_component(tintBase), 0.0001);
+  let lifted = select(
+    tint * minimumLuminance,
+    response * (minimumLuminance / max(responseLuminance, 0.0001)),
+    responseLuminance > 0.0001
+  );
+  return clamp(lifted, vec3<f32>(0.0), vec3<f32>(0.98));
 }
 
 fn sunlit_baseline_radiance(normal: vec3<f32>) -> vec3<f32> {
@@ -2331,12 +4320,24 @@ fn sunlit_baseline_radiance(normal: vec3<f32>) -> vec3<f32> {
   return clamp_sample_radiance(sunTint * baseline * skyFacing * directionalWeight * 0.04);
 }
 
-fn terminal_surface_environment_source(hit: HitRecord) -> vec3<f32> {
+fn terminal_surface_environment_source(ray: RayRecord, hit: HitRecord) -> vec3<f32> {
   let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
-  let normalEnvironment = gated_environment_radiance(
-    hit.position.xyz + normal * 0.003,
-    normal
+  let origin = hit.position.xyz + normal * 0.003;
+  let roughness = clamp(hit.material.x, 0.0, 1.0);
+  let glossiness = surface_glossiness(hit);
+  let normalEnvironment = gated_environment_radiance(origin, normal);
+  let viewDirection = safe_normalize(-ray.direction.xyz, normal);
+  let reflectionDirection = glossy_environment_direction(
+    ray.direction.xyz,
+    normal,
+    roughness,
+    mix(0.88, 0.38, glossiness)
   );
+  let reflectionEnvironment = prefiltered_environment_radiance(reflectionDirection, roughness);
+  let surfaceColor = clamp(max(hit.color.xyz, config.ambientColor.xyz * 0.35), vec3<f32>(0.0), vec3<f32>(1.0));
+  let f0 = surface_specular_f0(hit, surfaceColor);
+  let brdfTerm = sample_brdf_lut(saturate(dot(normal, viewDirection)), roughness);
+  let specularEnvironment = reflectionEnvironment * (f0 * brdfTerm.x + vec3<f32>(brdfTerm.y));
   let sunlitFloor = sunlit_baseline_radiance(normal);
   let ambientFloor = select(
     max(config.ambientColor.xyz, sunlitFloor * 0.82),
@@ -2348,17 +4349,27 @@ fn terminal_surface_environment_source(hit: HitRecord) -> vec3<f32> {
     max(config.environmentMapSettings.w, max(0.12, config.pathResolveSettings.y * 0.42)),
     environment_map_enabled()
   );
-  let environmentFloor = max(ambientFloor, max(sunlitFloor, normalEnvironment * environmentInfluence));
+  let glossyEnvironment = max(
+    normalEnvironment,
+    max(reflectionEnvironment * mix(0.24, 0.92, glossiness), specularEnvironment)
+  );
+  let environmentFloor = max(ambientFloor, max(sunlitFloor, glossyEnvironment * environmentInfluence));
   let materialFloor = select(0.7, 1.0, hit.materialKind == 0u || hit.materialKind == 3u);
   return clamp_sample_radiance(environmentFloor * materialFloor);
 }
 
-fn terminal_surface_environment_contribution(ray: RayRecord, hit: HitRecord) -> vec3<f32> {
+fn terminal_surface_environment_contribution(
+  ray: RayRecord,
+  throughput: vec3<f32>,
+  hit: HitRecord
+) -> vec3<f32> {
   let surfaceColor = max(hit.color.xyz, config.ambientColor.xyz);
+  let occlusion = mix(0.75, 1.0, clamp(hit.occlusion, 0.0, 1.0));
   return clamp_sample_radiance(
-    ray.throughput.xyz *
+    throughput *
     surfaceColor *
-    terminal_surface_environment_source(hit)
+    terminal_surface_environment_source(ray, hit) *
+    occlusion
   );
 }
 
@@ -2391,6 +4402,10 @@ fn direct_environment_portal_irradiance(origin: vec3<f32>, normal: vec3<f32>) ->
     );
     let area = max(portal.position.w, 0.0001);
     let distanceFalloff = clamp(area / max(distanceSquared, area * 0.25), 0.0, 2.5);
+    let traceDistance = max(sqrt(distanceSquared) - 0.01, 0.01);
+    if (scene_visibility_blocked(origin, direction, traceDistance)) {
+      continue;
+    }
     irradiance = irradiance +
       portal.color.rgb *
       portal.normal.w *
@@ -2402,48 +4417,79 @@ fn direct_environment_portal_irradiance(origin: vec3<f32>, normal: vec3<f32>) ->
   return irradiance;
 }
 
-fn surface_direct_environment_contribution(ray: RayRecord, hit: HitRecord) -> vec3<f32> {
-  let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
-  let origin = hit.position.xyz + normal * 0.003;
-  let viewDirection = safe_normalize(-ray.direction.xyz, normal);
-  let surfaceColor = clamp(max(hit.color.xyz, config.ambientColor.xyz * 0.35), vec3<f32>(0.0), vec3<f32>(1.0));
-  let roughness = clamp(hit.material.x, 0.0, 1.0);
-  let metallic = clamp(hit.material.y, 0.0, 1.0);
-
-  let normalEnvironment = gated_environment_radiance(origin, normal);
-  let skyVisibility = 0.35 + saturate(normal.y * 0.5 + 0.5) * 0.45;
-  let sunlitFloor = sunlit_baseline_radiance(normal);
-  let ambientIrradiance = max(
-    select(config.ambientColor.xyz * 0.72, config.ambientColor.xyz * 0.28, environment_map_enabled()),
-    sunlitFloor * select(0.72, 0.45, environment_map_enabled())
-  );
-  let environmentIrradianceScale = select(
-    max(0.16, config.pathResolveSettings.y * 0.45),
-    max(config.environmentMapSettings.w, max(0.16, config.pathResolveSettings.y * 0.45)),
-    environment_map_enabled()
+fn visibility_test_ray(origin: vec3<f32>, direction: vec3<f32>) -> RayRecord {
+  let rayDirection = safe_normalize(direction, vec3<f32>(0.0, 1.0, 0.0));
+  return RayRecord(
+    0u,
+    0u,
+    0u,
+    0u,
+    0u,
+    0u,
+    0u,
+    0u,
+    vec4<f32>(origin, 1.0),
+    vec4<f32>(rayDirection, 0.0),
+    vec4<f32>(1.0, 1.0, 1.0, 1.0)
   );
-  let skyIrradiance = max(ambientIrradiance, normalEnvironment * skyVisibility * environmentIrradianceScale);
+}
 
-  let sunDirection = safe_normalize(
-    config.environmentSunDirectionIntensity.xyz,
-    vec3<f32>(0.0, 1.0, 0.0)
-  );
-  let sunFacing = saturate(dot(normal, sunDirection));
-  let sunRadiance = gated_environment_radiance(origin, sunDirection);
-  let sunIrradiance = sunRadiance * sunFacing * 0.2;
-  let portalIrradiance = direct_environment_portal_irradiance(origin, normal);
+fn scene_visibility_blocked(origin: vec3<f32>, direction: vec3<f32>, maxDistance: f32) -> bool {
+  let testRay = visibility_test_ray(origin, direction);
+  let nearest = max(maxDistance, 0.001);
 
-  let diffuseWeight = select(1.0 - metallic * 0.65, 0.22, hit.materialKind == 1u);
-  let diffuse = surfaceColor * (skyIrradiance + sunIrradiance + portalIrradiance) * diffuseWeight;
+  for (var objectIndex = 0u; objectIndex < config.sceneObjectCount; objectIndex = objectIndex + 1u) {
+    let object = sceneObjects[objectIndex];
+    var current = no_candidate();
+    if (object.kind == 1u) {
+      current = intersect_sphere(testRay, object);
+    } else if (object.kind == 2u) {
+      current = intersect_box(testRay, object);
+    }
+    if (current.hit == 1u && current.distance < nearest) {
+      return true;
+    }
+  }
 
-  let halfVector = safe_normalize(sunDirection + viewDirection, normal);
-  let specularPower = 8.0 + (1.0 - roughness) * 96.0;
-  let specularFacing = pow(saturate(dot(normal, halfVector)), specularPower) * sunFacing;
-  let specularTint = mix(vec3<f32>(0.04), surfaceColor, metallic);
-  let specular = specularTint * sunRadiance * specularFacing * select(0.16, 0.48, hit.materialKind == 1u || hit.materialKind == 2u);
+  let meshCandidate = intersect_bvh(testRay, nearest);
+  return meshCandidate.hit == 1u && meshCandidate.distance < nearest;
+}
 
-  let bounceWeight = select(1.0, 0.38, ray.bounce > 0u);
-  return clamp_sample_radiance(ray.throughput.xyz * (diffuse + specular) * bounceWeight);
+fn surface_direct_environment_contribution(ray: RayRecord, hit: HitRecord) -> vec3<f32> {
+  let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
+  let origin = hit.position.xyz + normal * 0.003;
+  let viewDirection = safe_normalize(-ray.direction.xyz, normal);
+  let lightSample = sample_environment_importance(vec2<f32>(
+    random01(mix_seed(ray.sourcePixelId, ray.sampleId, ray.bounce, config.frameIndex, 41u)),
+    random01(mix_seed(ray.sourcePixelId, ray.sampleId, ray.bounce, config.frameIndex, 43u))
+  ));
+  if (lightSample.pdf <= 0.000001) {
+    return vec3<f32>(0.0);
+  }
+  let lightDirection = safe_normalize(lightSample.direction, normal);
+  let nDotL = saturate(dot(normal, lightDirection));
+  if (nDotL <= 0.000001) {
+    return vec3<f32>(0.0);
+  }
+  if (scene_visibility_blocked(origin, lightDirection, 1000000.0)) {
+    return vec3<f32>(0.0);
+  }
+  let incidentRadiance = direct_environment_radiance(origin, lightDirection);
+  if (max_component(incidentRadiance) <= 0.000001) {
+    return vec3<f32>(0.0);
+  }
+  let bsdf = evaluate_surface_bsdf(hit, viewDirection, lightDirection);
+  if (max_component(bsdf) <= 0.000001) {
+    return vec3<f32>(0.0);
+  }
+  let bsdfPdf = evaluate_surface_bsdf_pdf(hit, viewDirection, lightDirection);
+  let misWeight = power_heuristic(lightSample.pdf, bsdfPdf);
+  let contribution =
+    ray.throughput.xyz *
+    bsdf *
+    incidentRadiance *
+    (nDotL * misWeight / max(lightSample.pdf, 0.000001));
+  return clamp_sample_radiance(contribution);
 }
 
 fn default_mesh_range() -> MeshRange {
@@ -2462,7 +4508,16 @@ fn default_mesh_range() -> MeshRange {
     0u,
     vec4<f32>(0.72, 0.72, 0.68, 1.0),
     vec4<f32>(0.0),
-    vec4<f32>(0.72, 0.0, 1.0, 1.45)
+    vec4<f32>(0.72, 0.0, 1.0, 1.45),
+    vec4<f32>(0.0, 0.0, 0.0, 0.08),
+    vec4<f32>(1.0, 1.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(1.0, 1.0, 1.0, 0.0)
   );
 }
 
@@ -2558,7 +4613,7 @@ fn prepareMeshTrianglesAndLeaves(@builtin(global_invocation_id) globalId: vec3<u
     mesh.flags,
     mesh.materialRefId,
     mesh.mediumRefId,
-    0u,
+    mesh.materialSlot,
     0u,
     vec4<f32>(vertex0.position.xyz, 0.0),
     vec4<f32>(vertex1.position.xyz, 0.0),
@@ -2570,7 +4625,16 @@ fn prepareMeshTrianglesAndLeaves(@builtin(global_invocation_id) globalId: vec3<u
     vec4<f32>(uv2, 0.0, 0.0),
     mesh.color,
     mesh.emission,
-    mesh.material
+    mesh.material,
+    mesh.materialResponse,
+    mesh.materialExtension,
+    mesh.specularColor,
+    mesh.baseColorAtlas,
+    mesh.metallicRoughnessAtlas,
+    mesh.normalAtlas,
+    mesh.occlusionAtlas,
+    mesh.emissiveAtlas,
+    mesh.textureSettings
   );
 
   let leafBase = config.triangleCount - 1u;
@@ -2729,7 +4793,8 @@ fn make_miss(ray: RayRecord) -> HitRecord {
     0u,
     0u,
     -1.0,
-    vec3<f32>(0.0),
+    1.0,
+    vec2<f32>(0.0),
     vec4<f32>(ray.origin.xyz + ray.direction.xyz * 1000.0, 1.0),
     vec4<f32>(-ray.direction.xyz, 0.0),
     vec4<f32>(-ray.direction.xyz, 0.0),
@@ -2737,7 +4802,10 @@ fn make_miss(ray: RayRecord) -> HitRecord {
     vec4<f32>(0.0),
     vec4<f32>(radiance, 1.0),
     vec4<f32>(0.0),
-    vec4<f32>(1.0, 0.0, 1.0, 1.0)
+    vec4<f32>(1.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 0.0, 0.08),
+    vec4<f32>(0.08, 1.0, 0.0, 0.0),
+    vec4<f32>(1.0, 1.0, 1.0, 1.0)
   );
 }
 
@@ -2772,7 +4840,7 @@ fn intersect_sphere(ray: RayRecord, object: SceneObject) -> Candidate {
     0xffffffffu,
     object.objectId,
     object.objectId,
-    0u
+    object.mediumRefId
   );
 }
 
@@ -2824,7 +4892,7 @@ fn intersect_box(ray: RayRecord, object: SceneObject) -> Candidate {
     0xffffffffu,
     object.objectId,
     object.objectId,
-    0u
+    object.mediumRefId
   );
 }
 
@@ -3032,6 +5100,19 @@ fn denoise_range_space(value: vec3<f32>) -> vec3<f32> {
   return value / (vec3<f32>(1.0) + value);
 }
 
+fn denoise_sample_count() -> f32 {
+  return clamp(1.0 / max(config.projectionAndSampling.z, 0.000001), 1.0, 256.0);
+}
+
+fn denoise_strength() -> f32 {
+  let spp = denoise_sample_count();
+  return clamp(0.44 / sqrt(spp), 0.08, 0.44);
+}
+
+fn denoise_kernel_radius() -> i32 {
+  return select(1i, 2i, denoise_sample_count() < 2.5);
+}
+
 @compute @workgroup_size(64)
 fn generatePrimaryRays(@builtin(global_invocation_id) globalId: vec3<u32>) {
   let index = globalId.x;
@@ -3062,6 +5143,10 @@ fn intersectActiveQueue(@builtin(global_invocation_id) globalId: vec3<u32>) {
   let ray = activeQueue[index];
   var nearest = 1000000.0;
   var hitObject = SceneObject(
+    0u,
+    0u,
+    0u,
+    0u,
     0u,
     0u,
     0u,
@@ -3070,7 +5155,10 @@ fn intersectActiveQueue(@builtin(global_invocation_id) globalId: vec3<u32>) {
     vec4<f32>(0.0),
     vec4<f32>(0.0),
     vec4<f32>(0.0),
-    vec4<f32>(1.0, 0.0, 1.0, 1.0)
+    vec4<f32>(1.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 0.0, 0.08),
+    vec4<f32>(0.08, 1.0, 0.0, 0.0),
+    vec4<f32>(1.0, 1.0, 1.0, 1.0)
   );
   var candidate = no_candidate();
   var hitTriangle = TriangleRecord(
@@ -3092,7 +5180,16 @@ fn intersectActiveQueue(@builtin(global_invocation_id) globalId: vec3<u32>) {
     vec4<f32>(0.0),
     vec4<f32>(0.0),
     vec4<f32>(0.0),
-    vec4<f32>(1.0, 0.0, 1.0, 1.0)
+    vec4<f32>(1.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 0.0, 0.08),
+    vec4<f32>(0.08, 1.0, 0.0, 0.0),
+    vec4<f32>(1.0, 1.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(0.0, 0.0, 1.0, 1.0),
+    vec4<f32>(1.0, 1.0, 1.0, 0.0)
   );
 
   for (var objectIndex = 0u; objectIndex < config.sceneObjectCount; objectIndex = objectIndex + 1u) {
@@ -3125,16 +5222,28 @@ fn intersectActiveQueue(@builtin(global_invocation_id) globalId: vec3<u32>) {
   let position = ray.origin.xyz + ray.direction.xyz * candidate.distance;
   let hitMaterialKind = select(hitObject.materialKind, hitTriangle.materialKind, candidate.triangleIndex != 0xffffffffu);
   let hitObjectId = select(hitObject.objectId, hitTriangle.meshId, candidate.triangleIndex != 0xffffffffu);
-  let hitColor = select(hitObject.color, hitTriangle.color, candidate.triangleIndex != 0xffffffffu);
-  let hitEmission = select(hitObject.emission, hitTriangle.emission, candidate.triangleIndex != 0xffffffffu);
-  let hitMaterial = select(hitObject.material, hitTriangle.material, candidate.triangleIndex != 0xffffffffu);
+  let meshSurface = sample_surface_material(
+    hitTriangle,
+    candidate.uv,
+    candidate.geometricNormal,
+    candidate.shadingNormal
+  );
+  let hitColor = select(hitObject.color, meshSurface.color, candidate.triangleIndex != 0xffffffffu);
+  let hitEmission = select(hitObject.emission, meshSurface.emission, candidate.triangleIndex != 0xffffffffu);
+  let hitMaterial = select(hitObject.material, meshSurface.material, candidate.triangleIndex != 0xffffffffu);
+  let hitMaterialResponse = select(hitObject.materialResponse, meshSurface.materialResponse, candidate.triangleIndex != 0xffffffffu);
+  let hitMaterialExtension = select(hitObject.materialExtension, meshSurface.materialExtension, candidate.triangleIndex != 0xffffffffu);
+  let hitSpecularColor = select(hitObject.specularColor, meshSurface.specularColor, candidate.triangleIndex != 0xffffffffu);
+  let hitShadingNormal = select(candidate.shadingNormal, meshSurface.shadingNormal, candidate.triangleIndex != 0xffffffffu);
   let hitPrimitiveId = select(candidate.primitiveId, hitTriangle.triangleId, candidate.triangleIndex != 0xffffffffu);
   let hitMaterialRefId = select(candidate.materialRefId, hitTriangle.materialRefId, candidate.triangleIndex != 0xffffffffu);
   let hitMediumRefId = select(candidate.mediumRefId, hitTriangle.mediumRefId, candidate.triangleIndex != 0xffffffffu);
+  let hitMaterialSlot = select(0u, hitTriangle.materialSlot, candidate.triangleIndex != 0xffffffffu);
+  let hitOcclusion = select(1.0, meshSurface.occlusion, candidate.triangleIndex != 0xffffffffu);
   var hitType = 0u;
   if (hitMaterialKind == 4u || emission_power(hitEmission) > 0.0001) {
     hitType = 1u;
-  } else if (hitMaterialKind == 3u || hitMaterial.z < 0.999) {
+  } else if (hitMaterialKind == 3u || hitMaterial.z < 0.999 || hitMaterialExtension.z > 0.001) {
     hitType = 3u;
   }
   atomicAdd(&counters.hitCount, 1u);
@@ -3148,19 +5257,23 @@ fn intersectActiveQueue(@builtin(global_invocation_id) globalId: vec3<u32>) {
     hitPrimitiveId,
     hitMaterialRefId,
     hitMediumRefId,
-    0u,
+    hitMaterialSlot,
     0u,
     0u,
     candidate.distance,
-    vec3<f32>(0.0),
+    hitOcclusion,
+    vec2<f32>(0.0),
     vec4<f32>(position, 1.0),
     vec4<f32>(candidate.geometricNormal, 0.0),
-    vec4<f32>(candidate.shadingNormal, 0.0),
+    vec4<f32>(hitShadingNormal, 0.0),
     vec4<f32>(candidate.barycentric, 0.0),
     vec4<f32>(candidate.uv, 0.0, 0.0),
     hitColor,
     hitEmission,
-    hitMaterial
+    hitMaterial,
+    hitMaterialResponse,
+    hitMaterialExtension,
+    hitSpecularColor
   );
 }
 
@@ -3229,60 +5342,106 @@ fn sample_environment_portal_direction(hit: HitRecord, seed: u32, fallback: vec3
 }
 
 fn scatter_direction(ray: RayRecord, hit: HitRecord, seed: u32) -> ScatterResult {
+  let normal = safe_normalize(hit.shadingNormal.xyz, vec3<f32>(0.0, 1.0, 0.0));
+  let viewDirection = safe_normalize(-ray.direction.xyz, normal);
   let roughness = clamp(hit.material.x, 0.0, 1.0);
-  if (hit.materialKind == 1u) {
+  let transmission = clamp(hit.materialExtension.z, 0.0, 1.0);
+  if (hit.materialKind == 1u && roughness <= 0.02) {
     return ScatterResult(
-      vec4<f32>(
-        safe_normalize(
-          reflect(ray.direction.xyz, hit.shadingNormal.xyz) + random_unit_vector(seed) * roughness,
-          hit.shadingNormal.xyz
-        ),
-        0.0
-      ),
-      0u,
-      0u,
+      vec4<f32>(reflect(ray.direction.xyz, normal), 0.0),
+      1.0,
+      ray.mediumRefId,
+      RAY_FLAG_DELTA_SAMPLE,
       0u,
-      0u
     );
   }
 
-  if (hit.materialKind == 2u || hit.materialKind == 3u) {
+  if (hit.materialKind == 2u || hit.materialKind == 3u || transmission > 0.001) {
     let ior = max(hit.material.w, 1.01);
     let etaRatio = select(ior, 1.0 / ior, hit.frontFace == 1u);
-    let cosTheta = min(dot(-ray.direction.xyz, hit.shadingNormal.xyz), 1.0);
+    let cosTheta = min(dot(-ray.direction.xyz, normal), 1.0);
     let sinTheta = sqrt(max(0.0, 1.0 - cosTheta * cosTheta));
     let cannotRefract = etaRatio * sinTheta > 1.0;
     let reflectChance = schlick(cosTheta, etaRatio);
-    if (cannotRefract || random01(seed + 23u) < reflectChance) {
-      return ScatterResult(vec4<f32>(reflect(ray.direction.xyz, hit.shadingNormal.xyz), 0.0), 0u, 0u, 0u, 0u);
+    let transmissionReflectChance = select(
+      reflectChance,
+      max(reflectChance, 1.0 - transmission),
+      transmission > 0.001
+    );
+    if (cannotRefract || random01(seed + 23u) < transmissionReflectChance) {
+      return ScatterResult(
+        vec4<f32>(reflect(ray.direction.xyz, normal), 0.0),
+        1.0,
+        ray.mediumRefId,
+        RAY_FLAG_DELTA_SAMPLE,
+        0u,
+      );
     }
-    return ScatterResult(vec4<f32>(refract_direction(ray.direction.xyz, hit.shadingNormal.xyz, etaRatio), 0.0), 0u, 0u, 0u, 0u);
-  }
-
-  let randomDiffuse = safe_normalize(
-    hit.shadingNormal.xyz + random_unit_vector(seed),
-    hit.shadingNormal.xyz
-  );
-  let guidedLight = sample_emissive_triangle_direction(hit, seed, randomDiffuse);
-  let canSampleLight = dot(hit.shadingNormal.xyz, guidedLight) > -0.04;
-  let guideProbability = select(0.38, 0.72, ray.bounce == 0u);
-  let useGuidedLight = canSampleLight && random01(seed + 37u) < guideProbability;
-  let guidedPortal = sample_environment_portal_direction(hit, seed, randomDiffuse);
-  let canSamplePortal = dot(hit.shadingNormal.xyz, guidedPortal) > -0.04;
-  let useGuidedPortal =
-    !useGuidedLight &&
-    canSamplePortal &&
-    config.environmentPortalCount > 0u &&
-    config.environmentPortalMode > 0u &&
-    random01(seed + 89u) < 0.58;
-  let guidedDirection = select(randomDiffuse, guidedPortal, useGuidedPortal);
-  return ScatterResult(
-    vec4<f32>(select(guidedDirection, guidedLight, useGuidedLight), 0.0),
-    select(0u, RAY_FLAG_GUIDED_EMISSIVE, useGuidedLight),
-    0u,
-    0u,
-    0u
-  );
+    return ScatterResult(
+      vec4<f32>(refract_direction(ray.direction.xyz, normal, etaRatio), 0.0),
+      1.0,
+      transmitted_medium_ref_id(ray, hit),
+      RAY_FLAG_DELTA_SAMPLE,
+      0u,
+    );
+  }
+
+  let guidedEmissiveAvailable = config.emissiveTriangleCount > 0u;
+  let guidedPortalAvailable =
+    config.environmentPortalCount > 0u && config.environmentPortalMode != 0u;
+  let guidedSelector = random01(seed + 17u);
+  if (guidedEmissiveAvailable && guidedSelector < 0.18) {
+    let guidedDirection = sample_emissive_triangle_direction(hit, seed + 101u, normal);
+    if (dot(normal, guidedDirection) > 0.000001) {
+      let guidedPdf = max(evaluate_surface_bsdf_pdf(hit, viewDirection, guidedDirection), 0.000001);
+      return ScatterResult(
+        vec4<f32>(guidedDirection, 0.0),
+        guidedPdf,
+        ray.mediumRefId,
+        RAY_FLAG_GUIDED_EMISSIVE,
+        0u,
+      );
+    }
+  }
+  if (guidedPortalAvailable && guidedSelector < 0.32) {
+    let guidedDirection = sample_environment_portal_direction(hit, seed + 131u, normal);
+    if (dot(normal, guidedDirection) > 0.000001) {
+      let guidedPdf = max(evaluate_surface_bsdf_pdf(hit, viewDirection, guidedDirection), 0.000001);
+      return ScatterResult(vec4<f32>(guidedDirection, 0.0), guidedPdf, ray.mediumRefId, 0u, 0u);
+    }
+  }
+
+  let weights = surface_bsdf_sampling_weights(hit);
+  let selector = random01(seed + 31u);
+  var lightDirection = normal;
+  if (selector < weights.x) {
+    lightDirection = cosine_sample_hemisphere(
+      vec2<f32>(random01(seed + 37u), random01(seed + 41u)),
+      normal
+    );
+  } else if (selector < weights.x + weights.y) {
+    let halfVector = importance_sample_ggx(
+      vec2<f32>(random01(seed + 47u), random01(seed + 53u)),
+      max(roughness, 0.02),
+      normal
+    );
+    lightDirection = safe_normalize(reflect(-viewDirection, halfVector), normal);
+  } else {
+    let halfVector = importance_sample_ggx(
+      vec2<f32>(random01(seed + 59u), random01(seed + 61u)),
+      max(clamp(hit.materialExtension.x, 0.0, 1.0), 0.02),
+      normal
+    );
+    lightDirection = safe_normalize(reflect(-viewDirection, halfVector), normal);
+  }
+  if (dot(normal, lightDirection) <= 0.000001) {
+    lightDirection = cosine_sample_hemisphere(
+      vec2<f32>(random01(seed + 67u), random01(seed + 71u)),
+      normal
+    );
+  }
+  let pdf = max(evaluate_surface_bsdf_pdf(hit, viewDirection, lightDirection), 0.000001);
+  return ScatterResult(vec4<f32>(lightDirection, 0.0), pdf, ray.mediumRefId, 0u, 0u);
 }
 
 @compute @workgroup_size(64)
@@ -3295,15 +5454,17 @@ fn resolveSurfaceRecords(@builtin(global_invocation_id) globalId: vec3<u32>) {
 
   let ray = activeQueue[index];
   let hit = hits[index];
+  let segmentTransmittance = medium_transmittance(ray.mediumRefId, hit.distance);
+  let arrivingThroughput = ray.throughput.xyz * segmentTransmittance;
   var contribution = vec3<f32>(0.0);
 
   if (hit.hitType == 1u) {
     let guidedLightWeight = select(1.0, 0.24, (ray.flags & RAY_FLAG_GUIDED_EMISSIVE) != 0u);
     let sourceRadiance = max(hit.emission.xyz, hit.color.xyz) * guidedLightWeight;
     if (deferred_path_resolve_enabled()) {
-      record_deferred_terminal_source(ray, sourceRadiance);
+      record_deferred_terminal_source(ray, sourceRadiance * segmentTransmittance);
     } else {
-      contribution = clamp_sample_radiance(ray.throughput.xyz * sourceRadiance);
+      contribution = clamp_sample_radiance(arrivingThroughput * sourceRadiance);
       accumulation[ray.rayId] =
         accumulation[ray.rayId] + vec4<f32>(contribution * sample_weight(), 1.0);
     }
@@ -3312,10 +5473,17 @@ fn resolveSurfaceRecords(@builtin(global_invocation_id) globalId: vec3<u32>) {
   }
 
   if (hit.hitType == 2u) {
+    var sourceRadiance = hit.color.xyz;
+    if ((ray.flags & RAY_FLAG_DELTA_SAMPLE) == 0u) {
+      let bsdfPdf = max(ray.throughput.w, 0.000001);
+      let lightPdf = environment_direction_pdf(ray.direction.xyz);
+      let misWeight = power_heuristic(bsdfPdf, lightPdf);
+      sourceRadiance = sourceRadiance * misWeight;
+    }
     if (deferred_path_resolve_enabled()) {
-      record_deferred_terminal_source(ray, hit.color.xyz);
+      record_deferred_terminal_source(ray, sourceRadiance * segmentTransmittance);
     } else {
-      contribution = clamp_sample_radiance(ray.throughput.xyz * max(hit.color.xyz, config.ambientColor.xyz));
+      contribution = clamp_sample_radiance(arrivingThroughput * sourceRadiance);
       accumulation[ray.rayId] =
         accumulation[ray.rayId] + vec4<f32>(contribution * sample_weight(), 1.0);
     }
@@ -3323,24 +5491,47 @@ fn resolveSurfaceRecords(@builtin(global_invocation_id) globalId: vec3<u32>) {
     return;
   }
 
-  let response = surface_path_response(hit);
+  let response = stabilize_surface_path_response(
+    ray,
+    hit,
+    surface_path_response(hit) * segmentTransmittance
+  );
   record_deferred_path_response(ray, response);
 
   let shouldEstimateDirectEnvironment =
-    !deferred_path_resolve_enabled() &&
     (hit.materialKind == 0u || hit.materialKind == 1u) &&
-    hit.material.z >= 0.95;
+    hit.material.z >= 0.95 &&
+    ray.bounce < 2u;
   if (shouldEstimateDirectEnvironment) {
-    let directEnvironment = surface_direct_environment_contribution(ray, hit);
+    let directEnvironment = surface_direct_environment_contribution(
+      RayRecord(
+        ray.rayId,
+        ray.parentRayId,
+        ray.sourcePixelId,
+        ray.sampleId,
+        ray.bounce,
+        ray.mediumRefId,
+        ray.flags,
+        0u,
+        ray.origin,
+        ray.direction,
+        vec4<f32>(arrivingThroughput, ray.throughput.w)
+      ),
+      hit
+    );
     accumulation[ray.rayId] =
       accumulation[ray.rayId] + vec4<f32>(directEnvironment * sample_weight(), 0.0);
   }
 
   if (ray.bounce + 1u >= config.maxDepth) {
     if (deferred_path_resolve_enabled()) {
-      record_deferred_terminal_source(ray, terminal_surface_environment_source(hit));
+      record_deferred_terminal_source(ray, terminal_surface_environment_source(ray, hit));
     } else {
-      let terminalEnvironment = terminal_surface_environment_contribution(ray, hit);
+      let terminalEnvironment = terminal_surface_environment_contribution(
+        ray,
+        arrivingThroughput,
+        hit
+      );
       accumulation[ray.rayId] =
         accumulation[ray.rayId] + vec4<f32>(terminalEnvironment * sample_weight(), 1.0);
     }
@@ -3353,9 +5544,13 @@ fn resolveSurfaceRecords(@builtin(global_invocation_id) globalId: vec3<u32>) {
   let nextIndex = atomicAdd(&counters.nextCount, 1u);
   if (nextIndex >= config.tilePixelCount) {
     if (deferred_path_resolve_enabled()) {
-      record_deferred_terminal_source(ray, terminal_surface_environment_source(hit));
+      record_deferred_terminal_source(ray, terminal_surface_environment_source(ray, hit));
     } else {
-      let overflowEnvironment = terminal_surface_environment_contribution(ray, hit);
+      let overflowEnvironment = terminal_surface_environment_contribution(
+        ray,
+        arrivingThroughput,
+        hit
+      );
       accumulation[ray.rayId] =
         accumulation[ray.rayId] + vec4<f32>(overflowEnvironment * sample_weight(), 1.0);
     }
@@ -3369,12 +5564,12 @@ fn resolveSurfaceRecords(@builtin(global_invocation_id) globalId: vec3<u32>) {
     ray.sourcePixelId,
     ray.sampleId,
     ray.bounce + 1u,
-    ray.mediumRefId,
+    scatter.mediumRefId,
     scatter.flags,
     0u,
     vec4<f32>(offset_origin(hit.position.xyz, hit.shadingNormal.xyz), 1.0),
     scatter.direction,
-    vec4<f32>(throughput, ray.throughput.w)
+    vec4<f32>(throughput, scatter.pdf)
   );
 }
 
@@ -3443,8 +5638,11 @@ fn denoiseLinearRadiance(@builtin(global_invocation_id) globalId: vec3<u32>) {
 
   let pixel = vec2<i32>(i32(x), i32(y));
   let center = textureLoad(denoiseInputRadiance, pixel, 0).xyz;
-  var sum = center * 1.4;
-  var totalWeight = 1.4;
+  let strength = denoise_strength();
+  let kernelRadius = denoise_kernel_radius();
+  let centerWeight = 1.7 - strength * 0.35;
+  var sum = center * centerWeight;
+  var totalWeight = centerWeight;
   let centerRange = denoise_range_space(center);
 
   for (var oy = -2i; oy <= 2i; oy = oy + 1i) {
@@ -3452,13 +5650,16 @@ fn denoiseLinearRadiance(@builtin(global_invocation_id) globalId: vec3<u32>) {
       if (ox == 0i && oy == 0i) {
         continue;
       }
+      if (abs(ox) > kernelRadius || abs(oy) > kernelRadius) {
+        continue;
+      }
       let sx = clamp(i32(x) + ox, 0i, i32(config.canvasWidth) - 1i);
       let sy = clamp(i32(y) + oy, 0i, i32(config.canvasHeight) - 1i);
       let sampleColor = textureLoad(denoiseInputRadiance, vec2<i32>(sx, sy), 0).xyz;
       let colorDistance = length(denoise_range_space(sampleColor) - centerRange);
-      let rangeWeight = 1.0 / (1.0 + colorDistance * 7.0);
-      let distanceWeight = 1.0 / (1.0 + f32(ox * ox + oy * oy) * 0.24);
-      let diagonalWeight = select(1.0, 0.78, abs(ox) + abs(oy) > 2i);
+      let rangeWeight = 1.0 / (1.0 + colorDistance * (11.0 + strength * 6.0));
+      let distanceWeight = 1.0 / (1.0 + f32(ox * ox + oy * oy) * (0.62 + strength * 0.24));
+      let diagonalWeight = select(1.0, 0.92, abs(ox) + abs(oy) > 1i);
       let weight = rangeWeight * diagonalWeight * distanceWeight;
       sum = sum + sampleColor * weight;
       totalWeight = totalWeight + weight;
@@ -3466,8 +5667,9 @@ fn denoiseLinearRadiance(@builtin(global_invocation_id) globalId: vec3<u32>) {
   }
 
   let filtered = sum / max(totalWeight, 0.0001);
-  let outlier = saturate(length(denoise_range_space(center) - denoise_range_space(filtered)) * 2.4);
-  let color = min(mix(center, filtered, 0.52 + outlier * 0.18), vec3<f32>(16.0));
+  let outlier = saturate(length(denoise_range_space(center) - denoise_range_space(filtered)) * 2.1);
+  let blend = min(0.3, strength * (0.62 + outlier * 0.12));
+  let color = min(mix(center, filtered, blend), vec3<f32>(16.0));
   textureStore(denoisedRadianceImage, pixel, vec4<f32>(color, 1.0));
 }
 
@@ -3481,8 +5683,10 @@ fn resolveDenoisedOutputImage(@builtin(global_invocation_id) globalId: vec3<u32>
 
   let pixel = vec2<i32>(i32(x), i32(y));
   let center = textureLoad(finalDenoiseInputRadiance, pixel, 0).xyz;
-  var sum = center * 1.25;
-  var totalWeight = 1.25;
+  let strength = denoise_strength();
+  let centerWeight = 1.35 - strength * 0.25;
+  var sum = center * centerWeight;
+  var totalWeight = centerWeight;
   let centerRange = denoise_range_space(center);
 
   for (var oy = -1i; oy <= 1i; oy = oy + 1i) {
@@ -3494,8 +5698,8 @@ fn resolveDenoisedOutputImage(@builtin(global_invocation_id) globalId: vec3<u32>
       let sy = clamp(i32(y) + oy, 0i, i32(config.canvasHeight) - 1i);
       let sampleColor = textureLoad(finalDenoiseInputRadiance, vec2<i32>(sx, sy), 0).xyz;
       let colorDistance = length(denoise_range_space(sampleColor) - centerRange);
-      let rangeWeight = 1.0 / (1.0 + colorDistance * 9.0);
-      let distanceWeight = 1.0 / (1.0 + f32(ox * ox + oy * oy) * 0.4);
+      let rangeWeight = 1.0 / (1.0 + colorDistance * (12.0 + strength * 8.0));
+      let distanceWeight = 1.0 / (1.0 + f32(ox * ox + oy * oy) * (0.82 + strength * 0.28));
       let weight = rangeWeight * distanceWeight;
       sum = sum + sampleColor * weight;
       totalWeight = totalWeight + weight;
@@ -3503,8 +5707,9 @@ fn resolveDenoisedOutputImage(@builtin(global_invocation_id) globalId: vec3<u32>
   }
 
   let filtered = sum / max(totalWeight, 0.0001);
-  let outlier = saturate(length(denoise_range_space(center) - denoise_range_space(filtered)) * 2.8);
-  let radiance = min(mix(center, filtered, 0.28 + outlier * 0.12), vec3<f32>(16.0));
+  let outlier = saturate(length(denoise_range_space(center) - denoise_range_space(filtered)) * 2.2);
+  let blend = min(0.18, strength * (0.42 + outlier * 0.08));
+  let radiance = min(mix(center, filtered, blend), vec3<f32>(16.0));
   textureStore(denoisedOutputImage, pixel, vec4<f32>(tone_map_radiance(radiance), 1.0));
 }
 `;
@@ -3580,95 +5785,25 @@ function createAdapterInfoSnapshot(adapter) {
     vendor: typeof info.vendor === "string" ? info.vendor : "",
     architecture: typeof info.architecture === "string" ? info.architecture : "",
     device: typeof info.device === "string" ? info.device : "",
-    description: typeof info.description === "string" ? info.description : ""
-  });
-}
-function createGpuAdapterParallelismDiagnostics(adapter, device) {
-  return Object.freeze({
-    physicalCoreCount: null,
-    physicalCoreCountAvailable: false,
-    physicalCoreCountUnavailableReason: "WebGPU does not expose physical GPU core counts.",
-    adapterInfo: createAdapterInfoSnapshot(adapter),
-    adapterLimits: Object.freeze({
-      maxComputeInvocationsPerWorkgroup: readGpuLimit(adapter, device, "maxComputeInvocationsPerWorkgroup"),
-      maxComputeWorkgroupSizeX: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeX"),
-      maxComputeWorkgroupSizeY: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeY"),
-      maxComputeWorkgroupSizeZ: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeZ"),
-      maxComputeWorkgroupsPerDimension: readGpuLimit(adapter, device, "maxComputeWorkgroupsPerDimension"),
-      maxStorageBuffersPerShaderStage: readGpuLimit(adapter, device, "maxStorageBuffersPerShaderStage"),
-      maxStorageBufferBindingSize: readGpuLimit(adapter, device, "maxStorageBufferBindingSize")
-    }),
-    configuredWorkgroupSize: WORKGROUP_SIZE
-  });
-}
-function createGpuParallelismCounters() {
-  return {
-    directDispatches: 0,
-    directWorkgroups: 0,
-    directShaderInvocations: 0,
-    multiWorkgroupDispatches: 0,
-    largestDirectWorkgroupsPerDispatch: 0,
-    indirectDispatches: 0,
-    estimatedIndirectWorkgroupsUpperBound: 0,
-    estimatedIndirectShaderInvocationsUpperBound: 0,
-    indirectDispatchesWithMultiWorkgroupCapacity: 0,
-    largestEstimatedIndirectWorkgroupsPerDispatch: 0
-  };
-}
-function countDispatchWorkgroups(groups) {
-  return groups.reduce((product, value) => {
-    const numeric = Number(value ?? 1);
-    const count = Number.isFinite(numeric) ? Math.max(1, Math.trunc(numeric)) : 1;
-    return product * count;
-  }, 1);
-}
-function recordDirectDispatch(parallelism, groups, invocationsPerWorkgroup = WORKGROUP_SIZE) {
-  const workgroups = countDispatchWorkgroups(groups);
-  parallelism.directDispatches += 1;
-  parallelism.directWorkgroups += workgroups;
-  parallelism.directShaderInvocations += workgroups * invocationsPerWorkgroup;
-  parallelism.largestDirectWorkgroupsPerDispatch = Math.max(
-    parallelism.largestDirectWorkgroupsPerDispatch,
-    workgroups
-  );
-  if (workgroups > 1) {
-    parallelism.multiWorkgroupDispatches += 1;
-  }
-}
-function recordIndirectDispatch(parallelism, estimatedWorkgroupsUpperBound, invocationsPerWorkgroup = WORKGROUP_SIZE) {
-  const workgroups = Math.max(1, Math.trunc(Number(estimatedWorkgroupsUpperBound) || 1));
-  parallelism.indirectDispatches += 1;
-  parallelism.estimatedIndirectWorkgroupsUpperBound += workgroups;
-  parallelism.estimatedIndirectShaderInvocationsUpperBound += workgroups * invocationsPerWorkgroup;
-  parallelism.largestEstimatedIndirectWorkgroupsPerDispatch = Math.max(
-    parallelism.largestEstimatedIndirectWorkgroupsPerDispatch,
-    workgroups
-  );
-  if (workgroups > 1) {
-    parallelism.indirectDispatchesWithMultiWorkgroupCapacity += 1;
-  }
-}
-function createGpuParallelismDiagnostics(adapterDiagnostics, counters) {
-  const totalEstimatedWorkgroupsUpperBound = counters.directWorkgroups + counters.estimatedIndirectWorkgroupsUpperBound;
-  const totalEstimatedShaderInvocationsUpperBound = counters.directShaderInvocations + counters.estimatedIndirectShaderInvocationsUpperBound;
-  const exposesMultiWorkgroupParallelism = counters.multiWorkgroupDispatches > 0 || counters.indirectDispatchesWithMultiWorkgroupCapacity > 0;
-  return Object.freeze({
-    ...adapterDiagnostics,
-    directDispatches: counters.directDispatches,
-    directWorkgroups: counters.directWorkgroups,
-    directShaderInvocations: counters.directShaderInvocations,
-    multiWorkgroupDispatches: counters.multiWorkgroupDispatches,
-    largestDirectWorkgroupsPerDispatch: counters.largestDirectWorkgroupsPerDispatch,
-    indirectDispatches: counters.indirectDispatches,
-    estimatedIndirectWorkgroupsUpperBound: counters.estimatedIndirectWorkgroupsUpperBound,
-    estimatedIndirectShaderInvocationsUpperBound: counters.estimatedIndirectShaderInvocationsUpperBound,
-    indirectDispatchesWithMultiWorkgroupCapacity: counters.indirectDispatchesWithMultiWorkgroupCapacity,
-    largestEstimatedIndirectWorkgroupsPerDispatch: counters.largestEstimatedIndirectWorkgroupsPerDispatch,
-    totalEstimatedWorkgroupsUpperBound,
-    totalEstimatedShaderInvocationsUpperBound,
-    exposesMultiWorkgroupParallelism,
-    likelyUsesMoreThanOnePhysicalGpuCore: null,
-    coreUtilizationStatus: "not-exposed-by-webgpu"
+    description: typeof info.description === "string" ? info.description : ""
+  });
+}
+function createGpuAdapterParallelismDiagnostics(adapter, device) {
+  return Object.freeze({
+    physicalCoreCount: null,
+    physicalCoreCountAvailable: false,
+    physicalCoreCountUnavailableReason: "WebGPU does not expose physical GPU core counts.",
+    adapterInfo: createAdapterInfoSnapshot(adapter),
+    adapterLimits: Object.freeze({
+      maxComputeInvocationsPerWorkgroup: readGpuLimit(adapter, device, "maxComputeInvocationsPerWorkgroup"),
+      maxComputeWorkgroupSizeX: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeX"),
+      maxComputeWorkgroupSizeY: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeY"),
+      maxComputeWorkgroupSizeZ: readGpuLimit(adapter, device, "maxComputeWorkgroupSizeZ"),
+      maxComputeWorkgroupsPerDimension: readGpuLimit(adapter, device, "maxComputeWorkgroupsPerDimension"),
+      maxStorageBuffersPerShaderStage: readGpuLimit(adapter, device, "maxStorageBuffersPerShaderStage"),
+      maxStorageBufferBindingSize: readGpuLimit(adapter, device, "maxStorageBufferBindingSize")
+    }),
+    configuredWorkgroupSize: WORKGROUP_SIZE
   });
 }
 function createEnvironmentMapSnapshot(environmentMap) {
@@ -3676,12 +5811,38 @@ function createEnvironmentMapSnapshot(environmentMap) {
     enabled: environmentMap.enabled,
     width: environmentMap.width,
     height: environmentMap.height,
+    mipLevelCount: environmentMap.mipLevelCount ?? 1,
     projection: environmentMap.projection,
     intensity: environmentMap.intensity,
     rotationRadians: environmentMap.rotationRadians,
-    ambientStrength: environmentMap.ambientStrength
+    ambientStrength: environmentMap.ambientStrength,
+    hasImportanceData: environmentMap.hasImportanceData === true
   });
 }
+function nowMs() {
+  if (typeof performance !== "undefined" && typeof performance.now === "function") {
+    return performance.now();
+  }
+  return Date.now();
+}
+function estimateSubmittedGpuWorkTimeoutMs(config, tileCount, overrideTimeoutMs = null) {
+  if (Number.isFinite(overrideTimeoutMs)) {
+    return Math.max(1, Math.trunc(Number(overrideTimeoutMs)));
+  }
+  const samplesPerPixel = Math.max(
+    1,
+    Number(config?.renderedSamplesPerPixel ?? config?.samplesPerPixel ?? 1)
+  );
+  const maxDepth = Math.max(1, Number(config?.maxDepth ?? 1));
+  const deferredResolvePasses = config?.deferredPathResolve ? 1 : 0;
+  const denoisePasses = config?.denoise ? samplesPerPixel < 4 ? 2 : 1 : 0;
+  const tiles = Math.max(1, Number(tileCount ?? 1));
+  const estimatedPasses = tiles * (samplesPerPixel * (maxDepth + 1 + deferredResolvePasses) + denoisePasses + 1);
+  return Math.min(
+    GPU_MAX_SUBMITTED_WORK_TIMEOUT_MS,
+    GPU_SUBMITTED_WORK_TIMEOUT_MS + estimatedPasses * 5
+  );
+}
 async function createWavefrontPathTracingComputeRenderer(options = {}) {
   assertAnalyticDisplayQualityPolicy(options);
   const constants = getGpuUsageConstants();
@@ -3886,6 +6047,65 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     config.environmentMap,
     config.environmentColor
   );
+  const environmentSamplingResource = createEnvironmentSamplingTextureResource(
+    device,
+    constants,
+    config.environmentMap,
+    config.environmentColor
+  );
+  let mediumTextureResource = createMediumTextureResource(
+    device,
+    constants,
+    config.mediums
+  );
+  config = Object.freeze({
+    ...config,
+    environmentMap: Object.freeze({
+      ...config.environmentMap,
+      width: environmentMapResource.width,
+      height: environmentMapResource.height,
+      mipLevelCount: environmentMapResource.mipLevelCount,
+      hasImportanceData: environmentSamplingResource.hasImportanceData
+    })
+  });
+  const brdfLutResource = createBrdfLutResource(device, constants);
+  const baseColorAtlasResource = createAtlasTextureResource(
+    device,
+    constants,
+    config.gpuMaterialSource.baseColorAtlas,
+    "plasius.wavefront.materialAtlas.baseColor"
+  );
+  const metallicRoughnessAtlasResource = createAtlasTextureResource(
+    device,
+    constants,
+    config.gpuMaterialSource.metallicRoughnessAtlas,
+    "plasius.wavefront.materialAtlas.metallicRoughness"
+  );
+  const normalAtlasResource = createAtlasTextureResource(
+    device,
+    constants,
+    config.gpuMaterialSource.normalAtlas,
+    "plasius.wavefront.materialAtlas.normal"
+  );
+  const occlusionAtlasResource = createAtlasTextureResource(
+    device,
+    constants,
+    config.gpuMaterialSource.occlusionAtlas,
+    "plasius.wavefront.materialAtlas.occlusion"
+  );
+  const emissiveAtlasResource = createAtlasTextureResource(
+    device,
+    constants,
+    config.gpuMaterialSource.emissiveAtlas,
+    "plasius.wavefront.materialAtlas.emissive"
+  );
+  const materialAtlasSampler = device.createSampler({
+    label: "plasius.wavefront.materialAtlasSampler",
+    addressModeU: "clamp-to-edge",
+    addressModeV: "clamp-to-edge",
+    magFilter: "linear",
+    minFilter: "linear"
+  });
   const traceBindGroupLayout = device.createBindGroupLayout({
     label: "plasius.wavefront.traceBindGroupLayout",
     entries: [
@@ -3915,7 +6135,17 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       { binding: 19, visibility: constants.shader.COMPUTE, buffer: { type: "read-only-storage" } },
       { binding: 20, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
       { binding: 21, visibility: constants.shader.COMPUTE, sampler: { type: "filtering" } },
-      { binding: 22, visibility: constants.shader.COMPUTE, buffer: { type: "storage" } }
+      { binding: 22, visibility: constants.shader.COMPUTE, buffer: { type: "storage" } },
+      { binding: 23, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 24, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 25, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 26, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 27, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 28, visibility: constants.shader.COMPUTE, sampler: { type: "filtering" } },
+      { binding: 29, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 30, visibility: constants.shader.COMPUTE, sampler: { type: "filtering" } },
+      { binding: 31, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } },
+      { binding: 32, visibility: constants.shader.COMPUTE, texture: { sampleType: "float" } }
     ]
   });
   const accelerationBindGroupLayout = device.createBindGroupLayout({
@@ -3994,6 +6224,7 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     label: "plasius.wavefront.computeShader",
     code: WAVEFRONT_COMPUTE_WGSL
   });
+  await assertShaderModuleCompiles(computeShader, "plasius.wavefront.computeShader");
   const pipelines = {
     prepareMeshTrianglesAndLeaves: await createComputePipeline(
       device,
@@ -4092,14 +6323,27 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
         { binding: 19, resource: { buffer: environmentPortalBuffer } },
         { binding: 20, resource: environmentMapResource.view },
         { binding: 21, resource: environmentMapResource.sampler },
-        { binding: 22, resource: { buffer: pathVertexBuffer } }
+        { binding: 22, resource: { buffer: pathVertexBuffer } },
+        { binding: 23, resource: baseColorAtlasResource.view },
+        { binding: 24, resource: metallicRoughnessAtlasResource.view },
+        { binding: 25, resource: normalAtlasResource.view },
+        { binding: 26, resource: occlusionAtlasResource.view },
+        { binding: 27, resource: emissiveAtlasResource.view },
+        { binding: 28, resource: materialAtlasSampler },
+        { binding: 29, resource: brdfLutResource.view },
+        { binding: 30, resource: brdfLutResource.sampler },
+        { binding: 31, resource: environmentSamplingResource.view },
+        { binding: 32, resource: mediumTextureResource.view }
       ]
     });
   }
-  const bindGroups = [
-    createTraceBindGroup(activeQueue, nextQueue, "plasius.wavefront.bind.activeNext"),
-    createTraceBindGroup(nextQueue, activeQueue, "plasius.wavefront.bind.nextActive")
-  ];
+  function createTraceBindGroups() {
+    return [
+      createTraceBindGroup(activeQueue, nextQueue, "plasius.wavefront.bind.activeNext"),
+      createTraceBindGroup(nextQueue, activeQueue, "plasius.wavefront.bind.nextActive")
+    ];
+  }
+  let bindGroups = createTraceBindGroups();
   const bvhBuildBindGroup = device.createBindGroup({
     label: "plasius.wavefront.bind.bvhBuild",
     layout: accelerationBindGroupLayout,
@@ -4145,6 +6389,11 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     outputView,
     "plasius.wavefront.bind.denoise.scratchToOutput"
   );
+  const denoiseDirectResolveBindGroup = createDenoiseResolveBindGroup(
+    radianceView,
+    outputView,
+    "plasius.wavefront.bind.denoise.radianceToOutput"
+  );
   const presentBindGroupLayout = device.createBindGroupLayout({
     label: "plasius.wavefront.presentBindGroupLayout",
     entries: [
@@ -4182,23 +6431,129 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
   let accelerationBuilt = !config.gpuAccelerationBuildRequired;
   let accelerationBuildCount = 0;
   let activeCameraOptions = options.camera ?? DEFAULT_CAMERA;
+  let lastCompletedFrameTimeMs = null;
+  let lastCompletedSamplesPerPixel = Math.max(1, config.samplesPerPixel);
   let lastGpuParallelism = createGpuParallelismDiagnostics(
     gpuAdapterParallelism,
     createGpuParallelismCounters()
   );
+  function resolveRenderedSamplesPerPixel(renderOptions = {}, awaitGPUCompletion = true) {
+    const targetSamplesPerPixel = clamp(
+      readPositiveInteger(
+        "samplesPerPixel",
+        renderOptions.samplesPerPixel,
+        config.samplesPerPixel
+      ),
+      1,
+      config.samplesPerPixel
+    );
+    const frameTimeBudgetMs = Number.isFinite(renderOptions.frameTimeBudgetMs) ? Math.max(0, Number(renderOptions.frameTimeBudgetMs)) : null;
+    const minimumSamplesPerPixel = clamp(
+      readPositiveInteger(
+        "minimumSamplesPerPixel",
+        renderOptions.minimumSamplesPerPixel,
+        frameTimeBudgetMs !== null && targetSamplesPerPixel > 1 ? 1 : targetSamplesPerPixel
+      ),
+      1,
+      targetSamplesPerPixel
+    );
+    if (frameTimeBudgetMs === null || !awaitGPUCompletion || targetSamplesPerPixel <= minimumSamplesPerPixel) {
+      return Object.freeze({
+        renderedSamplesPerPixel: targetSamplesPerPixel,
+        targetSamplesPerPixel,
+        minimumSamplesPerPixel,
+        frameTimeBudgetMs,
+        budgetConstrained: false
+      });
+    }
+    const estimatedSampleTimeMs = Number.isFinite(lastCompletedFrameTimeMs) && lastCompletedFrameTimeMs > 0 ? lastCompletedFrameTimeMs / Math.max(1, lastCompletedSamplesPerPixel) : null;
+    if (!Number.isFinite(estimatedSampleTimeMs) || estimatedSampleTimeMs <= 0) {
+      return Object.freeze({
+        renderedSamplesPerPixel: minimumSamplesPerPixel,
+        targetSamplesPerPixel,
+        minimumSamplesPerPixel,
+        frameTimeBudgetMs,
+        budgetConstrained: minimumSamplesPerPixel < targetSamplesPerPixel
+      });
+    }
+    const budgetLimitedSamples = clamp(
+      Math.floor(frameTimeBudgetMs / estimatedSampleTimeMs),
+      minimumSamplesPerPixel,
+      targetSamplesPerPixel
+    );
+    return Object.freeze({
+      renderedSamplesPerPixel: budgetLimitedSamples,
+      targetSamplesPerPixel,
+      minimumSamplesPerPixel,
+      frameTimeBudgetMs,
+      budgetConstrained: budgetLimitedSamples < targetSamplesPerPixel
+    });
+  }
+  function createFrameStats({
+    frameIndex,
+    accelerationBuildSubmitted,
+    frameSubmissionCount,
+    parallelismCounters,
+    renderedSamplesPerPixel,
+    targetSamplesPerPixel,
+    frameTimeBudgetMs,
+    budgetConstrained
+  }) {
+    lastGpuParallelism = createGpuParallelismDiagnostics(gpuAdapterParallelism, parallelismCounters);
+    const commandSubmissions = frameSubmissionCount + (accelerationBuildSubmitted ? 1 : 0);
+    return Object.freeze({
+      frame,
+      frameIndex,
+      width: config.width,
+      height: config.height,
+      maxDepth: config.maxDepth,
+      tiles: tiles.length,
+      tileSize: config.tileSize,
+      samplesPerPixel: targetSamplesPerPixel,
+      renderedSamplesPerPixel,
+      frameTimeBudgetMs,
+      budgetConstrained,
+      maxFramePassesPerSubmission: config.maxFramePassesPerSubmission,
+      screenRays: config.width * config.height,
+      primaryRays: config.width * config.height * renderedSamplesPerPixel,
+      sceneObjectCount: config.sceneObjectCount,
+      triangleCount: config.triangleCount,
+      emissiveTriangleCount: config.emissiveTriangleCount,
+      environmentPortalCount: config.environmentPortalCount,
+      environmentPortalMode: config.environmentPortalMode,
+      mediumCount: config.mediumCount,
+      environmentMap: createEnvironmentMapSnapshot(config.environmentMap),
+      deferredPathResolve: config.deferredPathResolve,
+      bvhNodeCount: config.bvhNodeCount,
+      displayQuality: config.displayQuality,
+      accelerationBuildMode: config.accelerationBuildMode,
+      gpuAccelerationBuildRequired: config.gpuAccelerationBuildRequired,
+      accelerationBuildSubmitted,
+      accelerationBuilt,
+      accelerationBuildCount,
+      commandSubmissions,
+      frameConfigSlots: frameConfigSlotCount,
+      gpuParallelism: lastGpuParallelism,
+      memory: config.memory
+    });
+  }
+  function writeFrameConfigSlot(slot, tile, frameIndex, buildRange = {}) {
+    if (slot >= frameConfigSlotCount) {
+      throw new Error("Wavefront frame config slot capacity exceeded.");
+    }
+    const offset = slot * configBufferStride;
+    device.queue.writeBuffer(
+      configBuffer,
+      offset,
+      createConfigPayload(config, tile, frameIndex, buildRange)
+    );
+    return offset;
+  }
   function createFrameConfigWriter(frameIndex) {
     let slot = 0;
     return (tile, buildRange = {}) => {
-      if (slot >= frameConfigSlotCount) {
-        throw new Error("Wavefront frame config slot capacity exceeded.");
-      }
-      const offset = slot * configBufferStride;
+      const offset = writeFrameConfigSlot(slot, tile, frameIndex, buildRange);
       slot += 1;
-      device.queue.writeBuffer(
-        configBuffer,
-        offset,
-        createConfigPayload(config, tile, frameIndex, buildRange)
-      );
       return offset;
     };
   }
@@ -4243,7 +6598,7 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     passEncoder.setPipeline(pipelines.prepareMeshTrianglesAndLeaves);
     const prepareWorkgroups = Math.ceil(config.bvhLeafSortCapacity / WORKGROUP_SIZE);
     passEncoder.dispatchWorkgroups(prepareWorkgroups);
-    recordDirectDispatch(parallelism, [prepareWorkgroups]);
+    recordDirectDispatch(parallelism, [prepareWorkgroups], WORKGROUP_SIZE);
     passEncoder.setPipeline(pipelines.sortBvhLeafRefs);
     for (let stageIndex = 0; stageIndex < config.bvhSortStages.length; stageIndex += 1) {
       passEncoder.setBindGroup(0, bvhBuildBindGroup, [
@@ -4251,13 +6606,13 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       ]);
       const sortWorkgroups = Math.ceil(config.bvhLeafSortCapacity / WORKGROUP_SIZE);
       passEncoder.dispatchWorkgroups(sortWorkgroups);
-      recordDirectDispatch(parallelism, [sortWorkgroups]);
+      recordDirectDispatch(parallelism, [sortWorkgroups], WORKGROUP_SIZE);
     }
     passEncoder.setBindGroup(0, bvhBuildBindGroup, [0]);
     passEncoder.setPipeline(pipelines.writeSortedBvhLeaves);
     const leafWriteWorkgroups = Math.ceil(config.triangleCount / WORKGROUP_SIZE);
     passEncoder.dispatchWorkgroups(leafWriteWorkgroups);
-    recordDirectDispatch(parallelism, [leafWriteWorkgroups]);
+    recordDirectDispatch(parallelism, [leafWriteWorkgroups], WORKGROUP_SIZE);
     passEncoder.setPipeline(pipelines.buildBvhInternalLevel);
     for (let levelIndex = 0; levelIndex < config.bvhBuildLevels.length; levelIndex += 1) {
       const buildLevel = config.bvhBuildLevels[levelIndex];
@@ -4266,7 +6621,7 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       ]);
       const levelWorkgroups = Math.ceil(buildLevel.count / WORKGROUP_SIZE);
       passEncoder.dispatchWorkgroups(levelWorkgroups);
-      recordDirectDispatch(parallelism, [levelWorkgroups]);
+      recordDirectDispatch(parallelism, [levelWorkgroups], WORKGROUP_SIZE);
     }
     passEncoder.end();
     device.queue.submit([encoder.finish()]);
@@ -4282,7 +6637,7 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     generatePass.setBindGroup(0, bindGroups[0], [configOffset]);
     generatePass.setPipeline(pipelines.generatePrimaryRays);
     generatePass.dispatchWorkgroups(tileWorkgroups);
-    recordDirectDispatch(parallelism, [tileWorkgroups]);
+    recordDirectDispatch(parallelism, [tileWorkgroups], WORKGROUP_SIZE);
     generatePass.end();
     for (let bounceIndex = 0; bounceIndex < config.maxDepth; bounceIndex += 1) {
       encoder.copyBufferToBuffer(
@@ -4298,10 +6653,10 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       passEncoder.setBindGroup(0, bindGroups[bounceIndex % 2], [configOffset]);
       passEncoder.setPipeline(pipelines.intersectActiveQueue);
       passEncoder.dispatchWorkgroupsIndirect(activeDispatchBuffer, 0);
-      recordIndirectDispatch(parallelism, tileWorkgroups);
+      recordIndirectDispatch(parallelism, tileWorkgroups, WORKGROUP_SIZE);
       passEncoder.setPipeline(pipelines.resolveSurfaceRecords);
       passEncoder.dispatchWorkgroupsIndirect(activeDispatchBuffer, 0);
-      recordIndirectDispatch(parallelism, tileWorkgroups);
+      recordIndirectDispatch(parallelism, tileWorkgroups, WORKGROUP_SIZE);
       passEncoder.setPipeline(pipelines.compactAndSwapQueues);
       passEncoder.dispatchWorkgroups(1);
       recordDirectDispatch(parallelism, [1], 1);
@@ -4316,30 +6671,45 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     passEncoder.setBindGroup(0, bindGroups[0], [configOffset]);
     passEncoder.setPipeline(pipelines.accumulateTerminalRadiance);
     passEncoder.dispatchWorkgroups(tileWorkgroups);
-    recordDirectDispatch(parallelism, [tileWorkgroups]);
+    recordDirectDispatch(parallelism, [tileWorkgroups], WORKGROUP_SIZE);
     passEncoder.end();
   }
-  function encodeDenoise(encoder, configOffset, parallelism) {
+  function encodeDenoise(encoder, configOffset, parallelism, renderedSamplesPerPixel = config.samplesPerPixel) {
     if (!config.denoise) {
       return;
     }
     const denoiseWorkgroupsX = Math.ceil(config.width / 8);
     const denoiseWorkgroupsY = Math.ceil(config.height / 8);
-    const radiancePass = encoder.beginComputePass({
-      label: "plasius.wavefront.denoiseRadiancePass"
-    });
-    radiancePass.setBindGroup(0, denoiseRadianceBindGroup, [configOffset]);
-    radiancePass.setPipeline(pipelines.denoiseLinearRadiance);
-    radiancePass.dispatchWorkgroups(denoiseWorkgroupsX, denoiseWorkgroupsY);
-    recordDirectDispatch(parallelism, [denoiseWorkgroupsX, denoiseWorkgroupsY]);
-    radiancePass.end();
+    const useTwoPassDenoise = renderedSamplesPerPixel < 4;
+    if (useTwoPassDenoise) {
+      const radiancePass = encoder.beginComputePass({
+        label: "plasius.wavefront.denoiseRadiancePass"
+      });
+      radiancePass.setBindGroup(0, denoiseRadianceBindGroup, [configOffset]);
+      radiancePass.setPipeline(pipelines.denoiseLinearRadiance);
+      radiancePass.dispatchWorkgroups(denoiseWorkgroupsX, denoiseWorkgroupsY);
+      recordDirectDispatch(
+        parallelism,
+        [denoiseWorkgroupsX, denoiseWorkgroupsY],
+        WORKGROUP_SIZE
+      );
+      radiancePass.end();
+    }
     const resolvePass = encoder.beginComputePass({
       label: "plasius.wavefront.denoiseResolvePass"
     });
-    resolvePass.setBindGroup(0, denoiseResolveBindGroup, [configOffset]);
+    resolvePass.setBindGroup(
+      0,
+      useTwoPassDenoise ? denoiseResolveBindGroup : denoiseDirectResolveBindGroup,
+      [configOffset]
+    );
     resolvePass.setPipeline(pipelines.resolveDenoisedOutputImage);
     resolvePass.dispatchWorkgroups(denoiseWorkgroupsX, denoiseWorkgroupsY);
-    recordDirectDispatch(parallelism, [denoiseWorkgroupsX, denoiseWorkgroupsY]);
+    recordDirectDispatch(
+      parallelism,
+      [denoiseWorkgroupsX, denoiseWorkgroupsY],
+      WORKGROUP_SIZE
+    );
     resolvePass.end();
   }
   function encodePresent(encoder) {
@@ -4360,98 +6730,213 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     passEncoder.draw(3);
     passEncoder.end();
   }
-  function dispatchFrame(frameIndex, parallelism) {
+  function dispatchFrame(frameIndex, parallelism, renderedSamplesPerPixel = config.samplesPerPixel) {
     const writeFrameConfig = createFrameConfigWriter(frameIndex);
-    let submissionCount = 0;
-    let encodedFramePasses = 0;
-    let encoder = device.createCommandEncoder({
-      label: `plasius.wavefront.frame.${frameIndex}.batched.${submissionCount + 1}`
+    const batch = createGpuSubmissionBatcher({
+      device,
+      frameIndex,
+      maxFramePassesPerSubmission: config.maxFramePassesPerSubmission
     });
-    function submitCurrentEncoder() {
-      if (encodedFramePasses <= 0) {
-        return;
-      }
-      device.queue.submit([encoder.finish()]);
-      submissionCount += 1;
-      encodedFramePasses = 0;
-      encoder = device.createCommandEncoder({
-        label: `plasius.wavefront.frame.${frameIndex}.batched.${submissionCount + 1}`
-      });
-    }
-    function reserveEncoder(passCount = 1) {
-      if (encodedFramePasses > 0 && encodedFramePasses + passCount > config.maxFramePassesPerSubmission) {
-        submitCurrentEncoder();
-      }
-      encodedFramePasses += passCount;
-      return encoder;
-    }
     for (const tile of tiles) {
-      for (let sampleIndex = 0; sampleIndex < config.samplesPerPixel; sampleIndex += 1) {
+      for (let sampleIndex = 0; sampleIndex < renderedSamplesPerPixel; sampleIndex += 1) {
         const configOffset = writeFrameConfig(tile, {
           sampleIndex,
-          sampleWeight: 1 / config.samplesPerPixel
+          sampleWeight: 1 / renderedSamplesPerPixel
         });
-        encodeTileSample(reserveEncoder(), tile, configOffset, parallelism);
+        encodeTileSample(
+          batch.reserve(config.maxDepth + 1),
+          tile,
+          configOffset,
+          parallelism
+        );
         if (config.deferredPathResolve) {
-          encodeTileOutput(reserveEncoder(), tile, configOffset, parallelism);
+          encodeTileOutput(batch.reserve(1), tile, configOffset, parallelism);
         }
       }
       if (!config.deferredPathResolve) {
         const outputConfigOffset = writeFrameConfig(tile, {
           sampleIndex: 0,
-          sampleWeight: 1 / config.samplesPerPixel
+          sampleWeight: 1 / renderedSamplesPerPixel
         });
-        encodeTileOutput(reserveEncoder(), tile, outputConfigOffset, parallelism);
+        encodeTileOutput(batch.reserve(1), tile, outputConfigOffset, parallelism);
       }
     }
     if (config.denoise) {
       const denoiseConfigOffset = writeFrameConfig(
         { x: 0, y: 0, width: config.width, height: config.height },
-        { sampleIndex: 0, sampleWeight: 1 / config.samplesPerPixel }
+        { sampleIndex: 0, sampleWeight: 1 / renderedSamplesPerPixel }
+      );
+      const denoisePassCount = renderedSamplesPerPixel < 4 ? 2 : 1;
+      encodeDenoise(
+        batch.reserve(denoisePassCount),
+        denoiseConfigOffset,
+        parallelism,
+        renderedSamplesPerPixel
       );
-      encodeDenoise(reserveEncoder(), denoiseConfigOffset, parallelism);
     }
-    encodePresent(reserveEncoder());
-    submitCurrentEncoder();
-    return submissionCount;
+    encodePresent(batch.reserve(1));
+    return batch.flush();
   }
-  function renderOnce() {
+  function renderOnce(renderOptions = {}, resolvedSamplingPlan = null) {
+    const frameStartTimeMs = nowMs();
     frame += 1;
     const frameIndex = frame + config.frameIndex;
+    const samplingPlan = resolvedSamplingPlan ?? resolveRenderedSamplesPerPixel(renderOptions, false);
     const parallelismCounters = createGpuParallelismCounters();
     const accelerationBuildSubmitted = dispatchGpuAccelerationBuild(frameIndex, parallelismCounters);
-    const frameSubmissionCount = dispatchFrame(frameIndex, parallelismCounters);
-    lastGpuParallelism = createGpuParallelismDiagnostics(gpuAdapterParallelism, parallelismCounters);
+    const frameSubmissionCount = dispatchFrame(
+      frameIndex,
+      parallelismCounters,
+      samplingPlan.renderedSamplesPerPixel
+    );
+    const frameTimeMs = Math.max(0, nowMs() - frameStartTimeMs);
     return Object.freeze({
-      frame,
-      width: config.width,
-      height: config.height,
-      maxDepth: config.maxDepth,
-      tiles: tiles.length,
-      tileSize: config.tileSize,
-      samplesPerPixel: config.samplesPerPixel,
+      ...createFrameStats({
+        frameIndex,
+        accelerationBuildSubmitted,
+        frameSubmissionCount,
+        parallelismCounters,
+        renderedSamplesPerPixel: samplingPlan.renderedSamplesPerPixel,
+        targetSamplesPerPixel: samplingPlan.targetSamplesPerPixel,
+        frameTimeBudgetMs: samplingPlan.frameTimeBudgetMs,
+        budgetConstrained: samplingPlan.budgetConstrained
+      }),
+      gpuWorkerJobs: createGpuWorkerJobDiagnostics(
+        lastGpuParallelism,
+        frameSubmissionCount + (accelerationBuildSubmitted ? 1 : 0),
+        frameTimeMs,
+        false
+      )
+    });
+  }
+  async function waitForSubmittedGpuWork(options2 = {}) {
+    if (typeof device.queue.onSubmittedWorkDone !== "function") {
+      return true;
+    }
+    const timeoutMs = Math.max(
+      1,
+      Number.isFinite(options2.timeoutMs) ? Number(options2.timeoutMs) : GPU_SUBMITTED_WORK_TIMEOUT_MS
+    );
+    const allowTimeout = options2.allowTimeout !== false;
+    const completionPromise = device.queue.onSubmittedWorkDone().then(
+      () => ({ status: "done" }),
+      (error) => {
+        throw error;
+      }
+    );
+    const lossPromise = typeof device.lost?.then === "function" ? device.lost.then((info) => {
+      throw new Error(
+        `WebGPU device lost while waiting for submitted work (${info?.reason ?? "unknown"}).`
+      );
+    }) : null;
+    let timeoutHandle = null;
+    let resolveTimeoutPromise = null;
+    let timeoutSettled = false;
+    const settleTimeoutPromise = (value) => {
+      if (timeoutSettled) {
+        return;
+      }
+      timeoutSettled = true;
+      resolveTimeoutPromise?.(value);
+    };
+    const timeoutPromise = new Promise((resolve) => {
+      resolveTimeoutPromise = resolve;
+      timeoutHandle = setTimeout(() => settleTimeoutPromise({ status: "timeout" }), timeoutMs);
+    });
+    let result;
+    try {
+      result = await Promise.race(
+        [completionPromise, timeoutPromise, lossPromise].filter(Boolean)
+      );
+    } finally {
+      if (timeoutHandle !== null) {
+        clearTimeout(timeoutHandle);
+        settleTimeoutPromise({ status: "cancelled" });
+      }
+    }
+    if (result?.status === "timeout") {
+      if (!allowTimeout) {
+        throw new Error(`Timed out after ${timeoutMs} ms waiting for submitted GPU work.`);
+      }
+      console.warn(
+        `[plasius.wavefront] Submitted GPU work did not report completion within ${timeoutMs} ms; continuing.`
+      );
+      return false;
+    }
+    return true;
+  }
+  function dispatchFrameAwaitingGpu(frameIndex, parallelism, renderedSamplesPerPixel = config.samplesPerPixel) {
+    const samplePassesPerSample = config.maxDepth + 1 + (config.deferredPathResolve ? 1 : 0);
+    const denoisePassCount = config.denoise ? renderedSamplesPerPixel < 4 ? 2 : 1 : 0;
+    const tailPassCount = denoisePassCount + 1;
+    const sampleBatchSize = Math.max(
+      1,
+      Math.floor(
+        Math.max(config.maxFramePassesPerSubmission - tailPassCount, 1) / Math.max(samplePassesPerSample, 1)
+      )
+    );
+    let submissionCount = 0;
+    for (const tile of tiles) {
+      for (let sampleStart = 0; sampleStart < renderedSamplesPerPixel; sampleStart += sampleBatchSize) {
+        const sampleEnd = Math.min(renderedSamplesPerPixel, sampleStart + sampleBatchSize);
+        const batch = createGpuSubmissionBatcher({
+          device,
+          frameIndex,
+          maxFramePassesPerSubmission: config.maxFramePassesPerSubmission,
+          startingSubmissionCount: submissionCount
+        });
+        let slot = 0;
+        for (let sampleIndex = sampleStart; sampleIndex < sampleEnd; sampleIndex += 1) {
+          const configOffset = writeFrameConfigSlot(slot, tile, frameIndex, {
+            sampleIndex,
+            sampleWeight: 1 / renderedSamplesPerPixel
+          });
+          slot += 1;
+          encodeTileSample(
+            batch.reserve(config.maxDepth + 1),
+            tile,
+            configOffset,
+            parallelism
+          );
+          if (config.deferredPathResolve) {
+            encodeTileOutput(batch.reserve(1), tile, configOffset, parallelism);
+          }
+        }
+        if (!config.deferredPathResolve && sampleEnd >= renderedSamplesPerPixel) {
+          const outputConfigOffset = writeFrameConfigSlot(slot, tile, frameIndex, {
+            sampleIndex: 0,
+            sampleWeight: 1 / renderedSamplesPerPixel
+          });
+          encodeTileOutput(batch.reserve(1), tile, outputConfigOffset, parallelism);
+        }
+        batch.flush();
+        submissionCount += batch.getSubmissionCount();
+      }
+    }
+    const tail = createGpuSubmissionBatcher({
+      device,
+      frameIndex,
       maxFramePassesPerSubmission: config.maxFramePassesPerSubmission,
-      screenRays: config.width * config.height,
-      primaryRays: config.width * config.height * config.samplesPerPixel,
-      sceneObjectCount: config.sceneObjectCount,
-      triangleCount: config.triangleCount,
-      emissiveTriangleCount: config.emissiveTriangleCount,
-      environmentPortalCount: config.environmentPortalCount,
-      environmentPortalMode: config.environmentPortalMode,
-      environmentMap: createEnvironmentMapSnapshot(config.environmentMap),
-      deferredPathResolve: config.deferredPathResolve,
-      bvhNodeCount: config.bvhNodeCount,
-      displayQuality: config.displayQuality,
-      accelerationBuildMode: config.accelerationBuildMode,
-      gpuAccelerationBuildRequired: config.gpuAccelerationBuildRequired,
-      accelerationBuildSubmitted,
-      accelerationBuilt,
-      accelerationBuildCount,
-      commandSubmissions: frameSubmissionCount + (accelerationBuildSubmitted ? 1 : 0),
-      frameConfigSlots: frameConfigSlotCount,
-      gpuParallelism: lastGpuParallelism,
-      memory: config.memory
+      startingSubmissionCount: submissionCount
     });
+    if (config.denoise) {
+      const denoiseConfigOffset = writeFrameConfigSlot(
+        0,
+        { x: 0, y: 0, width: config.width, height: config.height },
+        frameIndex,
+        { sampleIndex: 0, sampleWeight: 1 / renderedSamplesPerPixel }
+      );
+      encodeDenoise(
+        tail.reserve(denoisePassCount),
+        denoiseConfigOffset,
+        parallelism,
+        renderedSamplesPerPixel
+      );
+    }
+    encodePresent(tail.reserve(1));
+    tail.flush();
+    submissionCount += tail.getSubmissionCount();
+    return submissionCount;
   }
   async function readOutputProbe(optionsForProbe = {}) {
     const mapMode = constants.map;
@@ -4465,6 +6950,10 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       size: 256,
       usage: constants.buffer.COPY_DST | constants.buffer.MAP_READ
     });
+    await waitForSubmittedGpuWork({
+      timeoutMs: GPU_READBACK_COMPLETION_TIMEOUT_MS,
+      allowTimeout: false
+    });
     const encoder = device.createCommandEncoder({
       label: "plasius.wavefront.outputProbe.copy"
     });
@@ -4474,6 +6963,10 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       { width: 1, height: 1, depthOrArrayLayers: 1 }
     );
     device.queue.submit([encoder.finish()]);
+    await waitForSubmittedGpuWork({
+      timeoutMs: GPU_READBACK_COMPLETION_TIMEOUT_MS,
+      allowTimeout: false
+    });
     await readback.mapAsync(mapMode.READ);
     const bytes = new Uint8Array(readback.getMappedRange()).slice(0, 4);
     readback.unmap();
@@ -4486,7 +6979,57 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     });
   }
   async function renderFrame(renderOptions = {}) {
-    const frameStats = renderOnce();
+    const awaitGPUCompletion = renderOptions.awaitGPUCompletion !== false;
+    const samplingPlan = resolveRenderedSamplesPerPixel(renderOptions, awaitGPUCompletion);
+    const useThrottledHighSamplePath = awaitGPUCompletion && samplingPlan.renderedSamplesPerPixel >= 8;
+    const submittedWorkTimeoutMs = estimateSubmittedGpuWorkTimeoutMs(
+      { ...config, renderedSamplesPerPixel: samplingPlan.renderedSamplesPerPixel },
+      tiles.length,
+      renderOptions.submittedWorkTimeoutMs
+    );
+    const frameStartTimeMs = nowMs();
+    const submissionWaitOptions = awaitGPUCompletion ? { timeoutMs: submittedWorkTimeoutMs, allowTimeout: false } : { timeoutMs: submittedWorkTimeoutMs };
+    let frameStats;
+    if (useThrottledHighSamplePath) {
+      frame += 1;
+      const frameIndex = frame + config.frameIndex;
+      const parallelismCounters = createGpuParallelismCounters();
+      const accelerationBuildSubmitted = dispatchGpuAccelerationBuild(frameIndex, parallelismCounters);
+      const frameSubmissionCount = dispatchFrameAwaitingGpu(
+        frameIndex,
+        parallelismCounters,
+        samplingPlan.renderedSamplesPerPixel
+      );
+      frameStats = createFrameStats({
+        frameIndex,
+        accelerationBuildSubmitted,
+        frameSubmissionCount,
+        parallelismCounters,
+        renderedSamplesPerPixel: samplingPlan.renderedSamplesPerPixel,
+        targetSamplesPerPixel: samplingPlan.targetSamplesPerPixel,
+        frameTimeBudgetMs: samplingPlan.frameTimeBudgetMs,
+        budgetConstrained: samplingPlan.budgetConstrained
+      });
+    } else {
+      frameStats = renderOnce(renderOptions, samplingPlan);
+    }
+    if (awaitGPUCompletion) {
+      await waitForSubmittedGpuWork(submissionWaitOptions);
+    }
+    const frameTimeMs = Math.max(0, nowMs() - frameStartTimeMs);
+    if (awaitGPUCompletion) {
+      lastCompletedFrameTimeMs = frameTimeMs;
+      lastCompletedSamplesPerPixel = frameStats.renderedSamplesPerPixel ?? frameStats.samplesPerPixel;
+    }
+    frameStats = Object.freeze({
+      ...frameStats,
+      gpuWorkerJobs: createGpuWorkerJobDiagnostics(
+        frameStats.gpuParallelism,
+        frameStats.commandSubmissions,
+        frameTimeMs,
+        awaitGPUCompletion
+      )
+    });
     const probe = renderOptions.readOutputProbe === false ? null : await readOutputProbe(renderOptions.probe);
     const maxChannel = probe ? Math.max(...probe.rgba.slice(0, 3)) : 0;
     return Object.freeze({
@@ -4507,10 +7050,8 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       queueOverflow: 0
     });
   }
-  function updateSceneObjects(sceneObjects) {
-    const nextPackedScene = packWavefrontSceneObjects(sceneObjects, config.sceneObjectCapacity);
-    packedScene = nextPackedScene;
-    config = createWavefrontPathTracingComputeConfig({
+  function rebuildLiveConfig(overrides = {}) {
+    return createWavefrontPathTracingComputeConfig({
       ...options,
       canvas,
       width: config.width,
@@ -4521,26 +7062,35 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       sceneObjectCapacity: config.sceneObjectCapacity,
       sceneObjects: packedScene.objects,
       camera: activeCameraOptions,
-      frameIndex: config.frameIndex
+      environmentMap: {
+        ...config.environmentMap
+      },
+      frameIndex: config.frameIndex,
+      ...overrides
     });
+  }
+  function rebuildMediumResources(nextConfig) {
+    const previousMediumTextureResource = mediumTextureResource;
+    mediumTextureResource = createMediumTextureResource(device, constants, nextConfig.mediums);
+    bindGroups = createTraceBindGroups();
+    if (previousMediumTextureResource?.ownsTexture) {
+      previousMediumTextureResource.texture?.destroy?.();
+    }
+  }
+  function updateSceneObjects(sceneObjects) {
+    const nextPackedScene = packWavefrontSceneObjects(sceneObjects, config.sceneObjectCapacity);
+    packedScene = nextPackedScene;
+    const nextConfig = rebuildLiveConfig();
+    if (!mediumTablesEqual(config.mediums, nextConfig.mediums)) {
+      rebuildMediumResources(nextConfig);
+    }
+    config = nextConfig;
     device.queue.writeBuffer(sceneObjectBuffer, 0, packedScene.buffer);
     return config;
   }
   function updateCamera(cameraOptions = {}) {
     activeCameraOptions = cameraOptions;
-    config = createWavefrontPathTracingComputeConfig({
-      ...options,
-      canvas,
-      width: config.width,
-      height: config.height,
-      maxDepth: config.maxDepth,
-      tileSize: config.tileSize,
-      samplesPerPixel: config.samplesPerPixel,
-      sceneObjectCapacity: config.sceneObjectCapacity,
-      sceneObjects: packedScene.objects,
-      camera: activeCameraOptions,
-      frameIndex: config.frameIndex
-    });
+    config = rebuildLiveConfig();
     return config;
   }
   function getSnapshot() {
@@ -4558,6 +7108,7 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
       emissiveTriangleCount: config.emissiveTriangleCount,
       environmentPortalCount: config.environmentPortalCount,
       environmentPortalMode: config.environmentPortalMode,
+      mediumCount: config.mediumCount,
       environmentMap: createEnvironmentMapSnapshot(config.environmentMap),
       deferredPathResolve: config.deferredPathResolve,
       bvhNodeCount: config.bvhNodeCount,
@@ -4595,6 +7146,28 @@ async function createWavefrontPathTracingComputeRenderer(options = {}) {
     if (environmentMapResource.ownsTexture) {
       environmentMapResource.texture?.destroy?.();
     }
+    if (environmentSamplingResource.ownsTexture) {
+      environmentSamplingResource.texture?.destroy?.();
+    }
+    if (mediumTextureResource.ownsTexture) {
+      mediumTextureResource.texture?.destroy?.();
+    }
+    brdfLutResource.texture?.destroy?.();
+    if (baseColorAtlasResource.ownsTexture) {
+      baseColorAtlasResource.texture?.destroy?.();
+    }
+    if (metallicRoughnessAtlasResource.ownsTexture) {
+      metallicRoughnessAtlasResource.texture?.destroy?.();
+    }
+    if (normalAtlasResource.ownsTexture) {
+      normalAtlasResource.texture?.destroy?.();
+    }
+    if (occlusionAtlasResource.ownsTexture) {
+      occlusionAtlasResource.texture?.destroy?.();
+    }
+    if (emissiveAtlasResource.ownsTexture) {
+      emissiveAtlasResource.texture?.destroy?.();
+    }
     context.unconfigure?.();
   }
   return Object.freeze({
@@ -4747,6 +7320,48 @@ var rendererAccelerationStructurePolicies = Object.freeze(
     })
   )
 );
+function clampWavefrontAdaptiveSamplesPerPixel(value) {
+  if (!Number.isFinite(value)) {
+    return 1;
+  }
+  return Math.max(1, Math.min(256, Math.round(value)));
+}
+function createWavefrontAdaptiveSamplingLevels(options = {}) {
+  const requestedSamplesPerPixel = clampWavefrontAdaptiveSamplesPerPixel(
+    options.samplesPerPixel ?? 1
+  );
+  const minimumSamplesPerPixel = Math.min(
+    requestedSamplesPerPixel,
+    clampWavefrontAdaptiveSamplesPerPixel(options.minimumSamplesPerPixel ?? 1)
+  );
+  const frameTimeBudgetMs = Number.isFinite(options.frameTimeBudgetMs) ? Math.max(0, Number(options.frameTimeBudgetMs)) : 0;
+  const levels = /* @__PURE__ */ new Set([minimumSamplesPerPixel, requestedSamplesPerPixel]);
+  let currentSamplesPerPixel = minimumSamplesPerPixel;
+  while (currentSamplesPerPixel < requestedSamplesPerPixel) {
+    levels.add(currentSamplesPerPixel);
+    currentSamplesPerPixel *= 2;
+  }
+  levels.add(Math.min(currentSamplesPerPixel, requestedSamplesPerPixel));
+  return Object.freeze({
+    requestedSamplesPerPixel,
+    minimumSamplesPerPixel,
+    frameTimeBudgetMs,
+    levels: Object.freeze(
+      [...levels].sort((left, right) => left - right).map(
+        (samplesPerPixel) => Object.freeze({
+          id: `${samplesPerPixel}spp`,
+          label: `${samplesPerPixel} spp`,
+          estimatedCostMs: samplesPerPixel,
+          config: Object.freeze({
+            samplesPerPixel,
+            frameTimeBudgetMs,
+            minimumSamplesPerPixel
+          })
+        })
+      )
+    )
+  });
+}
 function createWavefrontField(name, type, description) {
   return Object.freeze({
     name,
@@ -5924,9 +8539,11 @@ var defaultRendererClearColor = DEFAULT_CLEAR_COLOR;
   createGpuRenderer,
   createRayTracingRenderPlan,
   createRendererDebugHooks,
+  createWavefrontAdaptiveSamplingLevels,
   createWavefrontBvhBuildLevels,
   createWavefrontBvhSortStages,
   createWavefrontEmissiveTriangleIndexSource,
+  createWavefrontGpuMaterialSource,
   createWavefrontGpuMeshSource,
   createWavefrontMeshAcceleration,
   createWavefrontPathTracingComputeConfig,