@meimu/mutracer 0.1.5 → 0.3.0

package/lib/renderer/shaderCapabilities.d.ts

@@ -0,0 +1,23 @@
+export declare const MATERIAL_SHADER_CAPABILITIES: {
+    readonly transmission: 1;
+    readonly metal: 2;
+};
+export declare const LIGHT_SHADER_CAPABILITIES: {
+    readonly SphereLight: 1;
+    readonly RectLight: 2;
+    readonly SpotLight: 4;
+    readonly DirectionalLight: 8;
+    readonly emissive: 16;
+};
+export declare const ALL_MATERIAL_SHADER_CAPABILITIES = 3;
+export declare const ALL_LIGHT_SHADER_CAPABILITIES = 31;
+export declare function lightShaderCapabilities(types: readonly string[], emissiveCount: number): number;
+export declare function updateSceneShaderCapabilities(constants: Record<string, boolean | number>): void;
+export interface ShaderCapabilities {
+    SCENE_CAPABILITIES: number;
+    MATERIAL_CAPABILITIES: number;
+    LIGHT_CAPABILITIES: number;
+}
+export declare function shaderCapabilitiesInclude(supported: ShaderCapabilities, required: ShaderCapabilities): boolean;
+/** Retain support, not the scene's active flags. Commit only after compilation succeeds. */
+export declare function retainedShaderCapabilities(constants: Record<string, boolean | number>, previous?: ShaderCapabilities): ShaderCapabilities;

package/lib/renderer/shaderCapabilities.js

@@ -0,0 +1,55 @@
+import { wgpuSceneFeatureFlags, WGPU_SCENE_FEATURE_BITS } from './WGPUSceneFeatures.js';
+export const MATERIAL_SHADER_CAPABILITIES = { transmission: 1, metal: 2 };
+export const LIGHT_SHADER_CAPABILITIES = {
+    SphereLight: 1, RectLight: 2, SpotLight: 4, DirectionalLight: 8, emissive: 16
+};
+export const ALL_MATERIAL_SHADER_CAPABILITIES = 3;
+export const ALL_LIGHT_SHADER_CAPABILITIES = 31;
+export function lightShaderCapabilities(types, emissiveCount) {
+    let mask = emissiveCount > 0 ? LIGHT_SHADER_CAPABILITIES.emissive : 0;
+    for (const type of types) {
+        switch (type) {
+            case 'SphereLight':
+                mask |= LIGHT_SHADER_CAPABILITIES.SphereLight;
+                break;
+            case 'RectLight':
+                mask |= LIGHT_SHADER_CAPABILITIES.RectLight;
+                break;
+            case 'SpotLight':
+                mask |= LIGHT_SHADER_CAPABILITIES.SpotLight;
+                break;
+            case 'DirectionalLight':
+                mask |= LIGHT_SHADER_CAPABILITIES.DirectionalLight;
+                break;
+        }
+    }
+    return mask;
+}
+export function updateSceneShaderCapabilities(constants) {
+    constants.SCENE_CAPABILITIES = wgpuSceneFeatureFlags(constants);
+}
+export function shaderCapabilitiesInclude(supported, required) {
+    return (supported.SCENE_CAPABILITIES & required.SCENE_CAPABILITIES) === required.SCENE_CAPABILITIES &&
+        (supported.MATERIAL_CAPABILITIES & required.MATERIAL_CAPABILITIES) === required.MATERIAL_CAPABILITIES &&
+        (supported.LIGHT_CAPABILITIES & required.LIGHT_CAPABILITIES) === required.LIGHT_CAPABILITIES;
+}
+/** Retain support, not the scene's active flags. Commit only after compilation succeeds. */
+export function retainedShaderCapabilities(constants, previous) {
+    let scene = wgpuSceneFeatureFlags(constants) | (previous?.SCENE_CAPABILITIES ?? 0);
+    // Metalness is a continuous base-material edit, including the first change
+    // from a dielectric to a metal. Keep that support in every base variant.
+    const material = MATERIAL_SHADER_CAPABILITIES.metal |
+        Number(constants.MATERIAL_CAPABILITIES ?? 0) | (previous?.MATERIAL_CAPABILITIES ?? 0);
+    const commonSurface = WGPU_SCENE_FEATURE_BITS.HAS_COAT | WGPU_SCENE_FEATURE_BITS.HAS_FUZZ;
+    // Activate this small group on first use, keeping a plain scene's initial
+    // pipeline lean. Thin film, media, dispersion and alpha remain independent.
+    if ((scene & commonSurface) !== 0) {
+        scene |= commonSurface;
+    }
+    return {
+        SCENE_CAPABILITIES: scene,
+        MATERIAL_CAPABILITIES: material,
+        LIGHT_CAPABILITIES: Number(constants.LIGHT_CAPABILITIES ?? 0) | (previous?.LIGHT_CAPABILITIES ?? 0)
+    };
+}
+//# sourceMappingURL=shaderCapabilities.js.map

package/lib/renderer/shaderCapabilities.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"shaderCapabilities.js","sourceRoot":"","sources":["../../src/renderer/shaderCapabilities.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,qBAAqB,EAAE,uBAAuB,EAAC,MAAM,wBAAwB,CAAC;AAEtF,MAAM,CAAC,MAAM,4BAA4B,GAAG,EAAC,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAU,CAAC;AACjF,MAAM,CAAC,MAAM,yBAAyB,GAAG;IACvC,WAAW,EAAE,CAAC,EAAE,SAAS,EAAE,CAAC,EAAE,SAAS,EAAE,CAAC,EAAE,gBAAgB,EAAE,CAAC,EAAE,QAAQ,EAAE,EAAE;CACrE,CAAC;AACX,MAAM,CAAC,MAAM,gCAAgC,GAAG,CAAC,CAAC;AAClD,MAAM,CAAC,MAAM,6BAA6B,GAAG,EAAE,CAAC;AAEhD,MAAM,UAAU,uBAAuB,CAAC,KAAwB,EAAE,aAAqB;IACrF,IAAI,IAAI,GAAG,aAAa,GAAG,CAAC,CAAC,CAAC,CAAC,yBAAyB,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;IACtE,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,QAAQ,IAAI,EAAE,CAAC;YACb,KAAK,aAAa;gBAAE,IAAI,IAAI,yBAAyB,CAAC,WAAW,CAAC;gBAAC,MAAM;YACzE,KAAK,WAAW;gBAAE,IAAI,IAAI,yBAAyB,CAAC,SAAS,CAAC;gBAAC,MAAM;YACrE,KAAK,WAAW;gBAAE,IAAI,IAAI,yBAAyB,CAAC,SAAS,CAAC;gBAAC,MAAM;YACrE,KAAK,kBAAkB;gBAAE,IAAI,IAAI,yBAAyB,CAAC,gBAAgB,CAAC;gBAAC,MAAM;QACrF,CAAC;IACH,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,MAAM,UAAU,6BAA6B,CAAC,SAA2C;IACvF,SAAS,CAAC,kBAAkB,GAAG,qBAAqB,CAAC,SAAS,CAAC,CAAC;AAClE,CAAC;AAQD,MAAM,UAAU,yBAAyB,CAAC,SAA6B,EAAE,QAA4B;IACnG,OAAO,CAAC,SAAS,CAAC,kBAAkB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,KAAK,QAAQ,CAAC,kBAAkB;QACjG,CAAC,SAAS,CAAC,qBAAqB,GAAG,QAAQ,CAAC,qBAAqB,CAAC,KAAK,QAAQ,CAAC,qBAAqB;QACrG,CAAC,SAAS,CAAC,kBAAkB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,KAAK,QAAQ,CAAC,kBAAkB,CAAC;AACjG,CAAC;AAED,4FAA4F;AAC5F,MAAM,UAAU,0BAA0B,CAAC,SAA2C,EAAE,QAA6B;IACnH,IAAI,KAAK,GAAG,qBAAqB,CAAC,SAAS,CAAC,GAAG,CAAC,QAAQ,EAAE,kBAAkB,IAAI,CAAC,CAAC,CAAC;IACnF,2EAA2E;IAC3E,yEAAyE;IACzE,MAAM,QAAQ,GAAG,4BAA4B,CAAC,KAAK;QACjD,MAAM,CAAC,SAAS,CAAC,qBAAqB,IAAI,CAAC,CAAC,GAAG,CAAC,QAAQ,EAAE,qBAAqB,IAAI,CAAC,CAAC,CAAC;IACxF,MAAM,aAAa,GAAG,uBAAuB,CAAC,QAAQ,GAAG,uBAAuB,CAAC,QAAQ,CAAC;IAC1F,0EAA0E;IAC1E,4EAA4E;IAC5E,IAAI,CAAC,KAAK,GAAG,aAAa,CAAC,KAAK,CAAC,EAAE,CAAC;QAClC,KAAK,IAAI,aAAa,CAAC;IACzB,CAAC;IACD,OAAO;QACL,kBAAkB,EAAE,KAAK;QACzB,qBAAqB,EAAE,QAAQ;QAC/B,kBAAkB,EAAE,MAAM,CAAC,SAAS,CAAC,kBAAkB,IAAI,CAAC,CAAC,GAAG,CAAC,QAAQ,EAAE,kBAAkB,IAAI,CAAC,CAAC;KACpG,CAAC;AACJ,CAAC","sourcesContent":["import {wgpuSceneFeatureFlags, WGPU_SCENE_FEATURE_BITS} from './WGPUSceneFeatures.js';\n\nexport const MATERIAL_SHADER_CAPABILITIES = {transmission: 1, metal: 2} as const;\nexport const LIGHT_SHADER_CAPABILITIES = {\n  SphereLight: 1, RectLight: 2, SpotLight: 4, DirectionalLight: 8, emissive: 16\n} as const;\nexport const ALL_MATERIAL_SHADER_CAPABILITIES = 3;\nexport const ALL_LIGHT_SHADER_CAPABILITIES = 31;\n\nexport function lightShaderCapabilities(types: readonly string[], emissiveCount: number) {\n  let mask = emissiveCount > 0 ? LIGHT_SHADER_CAPABILITIES.emissive : 0;\n  for (const type of types) {\n    switch (type) {\n      case 'SphereLight': mask |= LIGHT_SHADER_CAPABILITIES.SphereLight; break;\n      case 'RectLight': mask |= LIGHT_SHADER_CAPABILITIES.RectLight; break;\n      case 'SpotLight': mask |= LIGHT_SHADER_CAPABILITIES.SpotLight; break;\n      case 'DirectionalLight': mask |= LIGHT_SHADER_CAPABILITIES.DirectionalLight; break;\n    }\n  }\n  return mask;\n}\n\nexport function updateSceneShaderCapabilities(constants: Record<string, boolean | number>) {\n  constants.SCENE_CAPABILITIES = wgpuSceneFeatureFlags(constants);\n}\n\nexport interface ShaderCapabilities {\n  SCENE_CAPABILITIES: number;\n  MATERIAL_CAPABILITIES: number;\n  LIGHT_CAPABILITIES: number;\n}\n\nexport function shaderCapabilitiesInclude(supported: ShaderCapabilities, required: ShaderCapabilities) {\n  return (supported.SCENE_CAPABILITIES & required.SCENE_CAPABILITIES) === required.SCENE_CAPABILITIES &&\n    (supported.MATERIAL_CAPABILITIES & required.MATERIAL_CAPABILITIES) === required.MATERIAL_CAPABILITIES &&\n    (supported.LIGHT_CAPABILITIES & required.LIGHT_CAPABILITIES) === required.LIGHT_CAPABILITIES;\n}\n\n/** Retain support, not the scene's active flags. Commit only after compilation succeeds. */\nexport function retainedShaderCapabilities(constants: Record<string, boolean | number>, previous?: ShaderCapabilities): ShaderCapabilities {\n  let scene = wgpuSceneFeatureFlags(constants) | (previous?.SCENE_CAPABILITIES ?? 0);\n  // Metalness is a continuous base-material edit, including the first change\n  // from a dielectric to a metal. Keep that support in every base variant.\n  const material = MATERIAL_SHADER_CAPABILITIES.metal |\n    Number(constants.MATERIAL_CAPABILITIES ?? 0) | (previous?.MATERIAL_CAPABILITIES ?? 0);\n  const commonSurface = WGPU_SCENE_FEATURE_BITS.HAS_COAT | WGPU_SCENE_FEATURE_BITS.HAS_FUZZ;\n  // Activate this small group on first use, keeping a plain scene's initial\n  // pipeline lean. Thin film, media, dispersion and alpha remain independent.\n  if ((scene & commonSurface) !== 0) {\n    scene |= commonSurface;\n  }\n  return {\n    SCENE_CAPABILITIES: scene,\n    MATERIAL_CAPABILITIES: material,\n    LIGHT_CAPABILITIES: Number(constants.LIGHT_CAPABILITIES ?? 0) | (previous?.LIGHT_CAPABILITIES ?? 0)\n  };\n}\n"]}

package/lib/renderer/shaderCapacity.d.ts

@@ -0,0 +1,2 @@
+/** Grow in coarse buckets; scene edits must not shrink a compiled stack. */
+export declare function retainedShaderCapacity(required: number, previous?: number, minimum?: number): number;

package/lib/renderer/shaderCapacity.js

@@ -0,0 +1,5 @@
+/** Grow in coarse buckets; scene edits must not shrink a compiled stack. */
+export function retainedShaderCapacity(required, previous = 0, minimum = 8) {
+    return Math.max(previous, minimum, 2 ** Math.ceil(Math.log2(Math.max(1, required))));
+}
+//# sourceMappingURL=shaderCapacity.js.map

package/lib/renderer/shaderCapacity.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"shaderCapacity.js","sourceRoot":"","sources":["../../src/renderer/shaderCapacity.ts"],"names":[],"mappings":"AAAA,4EAA4E;AAC5E,MAAM,UAAU,sBAAsB,CAAC,QAAgB,EAAE,QAAQ,GAAG,CAAC,EAAE,OAAO,GAAG,CAAC;IAChF,OAAO,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,OAAO,EAAE,CAAC,IAAI,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;AACvF,CAAC","sourcesContent":["/** Grow in coarse buckets; scene edits must not shrink a compiled stack. */\nexport function retainedShaderCapacity(required: number, previous = 0, minimum = 8) {\n  return Math.max(previous, minimum, 2 ** Math.ceil(Math.log2(Math.max(1, required))));\n}\n\n"]}

package/lib/renderer/skyModel.d.ts

@@ -0,0 +1,9 @@
+/** Apparent sun radius in degrees; larger than the real 0.27° so the disk reads at map resolution. */
+export declare const SUN_RADIUS = 0.4;
+/** Aerosol density relative to a very clean atmosphere (optical depth about 0.005). */
+export declare const TURBIDITY = 5;
+export declare const GROUND_ALBEDO: number[];
+/** Sun plus sky irradiance on a horizontal plane with the sun at 33°. */
+export declare const ILLUMINANCE: number;
+/** Horizontal illuminance of this atmosphere for unit solar irradiance and a 33° sun. */
+export declare const REFERENCE_ILLUMINANCE = 0.488;

package/lib/renderer/skyModel.js

@@ -0,0 +1,10 @@
+/** Apparent sun radius in degrees; larger than the real 0.27° so the disk reads at map resolution. */
+export const SUN_RADIUS = 0.4;
+/** Aerosol density relative to a very clean atmosphere (optical depth about 0.005). */
+export const TURBIDITY = 5;
+export const GROUND_ALBEDO = [0.2, 0.18, 0.15];
+/** Sun plus sky irradiance on a horizontal plane with the sun at 33°. */
+export const ILLUMINANCE = Math.PI;
+/** Horizontal illuminance of this atmosphere for unit solar irradiance and a 33° sun. */
+export const REFERENCE_ILLUMINANCE = 0.488;
+//# sourceMappingURL=skyModel.js.map

package/lib/renderer/skyModel.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"skyModel.js","sourceRoot":"","sources":["../../src/renderer/skyModel.ts"],"names":[],"mappings":"AAAA,sGAAsG;AACtG,MAAM,CAAC,MAAM,UAAU,GAAG,GAAG,CAAC;AAC9B,uFAAuF;AACvF,MAAM,CAAC,MAAM,SAAS,GAAG,CAAC,CAAC;AAC3B,MAAM,CAAC,MAAM,aAAa,GAAG,CAAC,GAAG,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;AAC/C,yEAAyE;AACzE,MAAM,CAAC,MAAM,WAAW,GAAG,IAAI,CAAC,EAAE,CAAC;AACnC,yFAAyF;AACzF,MAAM,CAAC,MAAM,qBAAqB,GAAG,KAAK,CAAC","sourcesContent":["/** Apparent sun radius in degrees; larger than the real 0.27° so the disk reads at map resolution. */\nexport const SUN_RADIUS = 0.4;\n/** Aerosol density relative to a very clean atmosphere (optical depth about 0.005). */\nexport const TURBIDITY = 5;\nexport const GROUND_ALBEDO = [0.2, 0.18, 0.15];\n/** Sun plus sky irradiance on a horizontal plane with the sun at 33°. */\nexport const ILLUMINANCE = Math.PI;\n/** Horizontal illuminance of this atmosphere for unit solar irradiance and a 33° sun. */\nexport const REFERENCE_ILLUMINANCE = 0.488;\n\n"]}

package/lib/renderer/spectrum.d.ts

@@ -0,0 +1,21 @@
+/**
+ * Continuous wavelength sampling shared by camera paths and caustic photons.
+ *
+ * CIE 1931 matching functions use the multi-lobe fit from Wyman, Sloan and
+ * Shirley, "Simple Analytic Approximations to the CIE XYZ Color Matching
+ * Functions", JCGT 2(2), 2013. Only the published formula and coefficients are
+ * used. Out-of-gamut negatives clamp to zero, and each channel is normalized so
+ * that a uniform wavelength keeps white light white.
+ */
+export declare const SPECTRUM_WAVELENGTH_MIN = 400;
+export declare const SPECTRUM_WAVELENGTH_RANGE = 300;
+export declare function spectrumRgb(wavelength: number): number[];
+export declare const SPECTRUM_PHOTON_DISTRIBUTION = 3;
+/** Maps a unit random number to a wavelength and its probability density. */
+export declare function spectrumSample(distribution: number, u: number): {
+    wavelength: number;
+    pdf: number;
+};
+export declare function spectrumIorWavelength(wavelength: number): number;
+export declare const spectrumWgsl: string;
+export declare const spectrumGlsl: string;

package/lib/renderer/spectrum.js

@@ -0,0 +1,154 @@
+/**
+ * Continuous wavelength sampling shared by camera paths and caustic photons.
+ *
+ * CIE 1931 matching functions use the multi-lobe fit from Wyman, Sloan and
+ * Shirley, "Simple Analytic Approximations to the CIE XYZ Color Matching
+ * Functions", JCGT 2(2), 2013. Only the published formula and coefficients are
+ * used. Out-of-gamut negatives clamp to zero, and each channel is normalized so
+ * that a uniform wavelength keeps white light white.
+ */
+export const SPECTRUM_WAVELENGTH_MIN = 400;
+export const SPECTRUM_WAVELENGTH_RANGE = 300;
+// Amplitude, mean, sigma below the mean and sigma above it.
+const XYZ_LOBES = [
+    [[1.056, 599.8, 37.9, 31.0], [0.362, 442.0, 16.0, 26.7], [-0.065, 501.1, 20.4, 26.2]],
+    [[0.821, 568.8, 46.9, 40.5], [0.286, 530.9, 16.3, 31.1]],
+    [[1.217, 437.0, 11.8, 36.0], [0.681, 459.0, 26.0, 13.8]]
+];
+// Linear sRGB primaries, D65.
+const XYZ_TO_RGB = [
+    [3.2406, -1.5372, -0.4986],
+    [-0.9689, 1.8758, 0.0415],
+    [0.0557, -0.204, 1.057]
+];
+function unnormalizedRgb(wavelength) {
+    const xyz = XYZ_LOBES.map((lobes) => {
+        let sum = 0;
+        for (const lobe of lobes) {
+            const t = (wavelength - lobe[1]) / (wavelength < lobe[1] ? lobe[2] : lobe[3]);
+            sum += lobe[0] * Math.exp(-0.5 * t * t);
+        }
+        return sum;
+    });
+    return XYZ_TO_RGB.map((row) => Math.max(0, row[0] * xyz[0] + row[1] * xyz[1] + row[2] * xyz[2]));
+}
+function channelScale() {
+    const steps = 3000;
+    const sum = [0, 0, 0];
+    for (let i = 0; i < steps; i++) {
+        const rgb = unnormalizedRgb(SPECTRUM_WAVELENGTH_MIN + (SPECTRUM_WAVELENGTH_RANGE * (i + 0.5)) / steps);
+        for (let c = 0; c < 3; c++)
+            sum[c] += rgb[c] / steps;
+    }
+    // Rounded so JavaScript and native builds generate identical shader text.
+    return sum.map((mean) => Number((1 / mean).toFixed(6)));
+}
+const SCALE = channelScale();
+export function spectrumRgb(wavelength) {
+    return unnormalizedRgb(wavelength).map((value, c) => value * SCALE[c]);
+}
+// Inverse CDFs in 5 nm bins: one per RGB channel for camera paths, and r+g+b for
+// photons. Narrow bins keep each sampled response close to the fit.
+const SAMPLE_BINS = 60;
+const SAMPLE_BIN_WIDTH = SPECTRUM_WAVELENGTH_RANGE / SAMPLE_BINS;
+export const SPECTRUM_PHOTON_DISTRIBUTION = 3;
+function sampleCdfs() {
+    const substeps = 20;
+    const cdfs = [];
+    for (let distribution = 0; distribution <= SPECTRUM_PHOTON_DISTRIBUTION; distribution++) {
+        const cdf = [0];
+        for (let bin = 0; bin < SAMPLE_BINS; bin++) {
+            let mass = 0;
+            for (let i = 0; i < substeps; i++) {
+                const rgb = spectrumRgb(SPECTRUM_WAVELENGTH_MIN + (bin + (i + 0.5) / substeps) * SAMPLE_BIN_WIDTH);
+                mass += distribution === SPECTRUM_PHOTON_DISTRIBUTION ? rgb[0] + rgb[1] + rgb[2] : rgb[distribution];
+            }
+            cdf.push(cdf[bin] + mass);
+        }
+        for (const value of cdf)
+            cdfs.push(Number((value / cdf[SAMPLE_BINS]).toFixed(7)));
+    }
+    return cdfs;
+}
+const SAMPLE_CDFS = sampleCdfs();
+/** Maps a unit random number to a wavelength and its probability density. */
+export function spectrumSample(distribution, u) {
+    const first = distribution * (SAMPLE_BINS + 1);
+    let low = 0;
+    let high = SAMPLE_BINS;
+    while (high - low > 1) {
+        const middle = (low + high) >> 1;
+        if (SAMPLE_CDFS[first + middle] <= u)
+            low = middle;
+        else
+            high = middle;
+    }
+    const mass = SAMPLE_CDFS[first + low + 1] - SAMPLE_CDFS[first + low];
+    return {
+        wavelength: SPECTRUM_WAVELENGTH_MIN + SAMPLE_BIN_WIDTH * (low + (u - SAMPLE_CDFS[first + low]) / mass),
+        pdf: mass / SAMPLE_BIN_WIDTH
+    };
+}
+// Cauchy IOR grows quickly at short wavelengths and strong dispersion scales turn
+// violet light into total internal reflection. Below 500 nm the refracting
+// wavelength eases toward 475 nm, while color still follows the sampled wavelength.
+const IOR_KNEE = 500;
+const IOR_KNEE_WIDTH = 25;
+export function spectrumIorWavelength(wavelength) {
+    if (wavelength >= IOR_KNEE)
+        return wavelength;
+    return IOR_KNEE - IOR_KNEE_WIDTH * (1 - Math.exp((wavelength - IOR_KNEE) / IOR_KNEE_WIDTH));
+}
+const CDF_LITERALS = SAMPLE_CDFS.map((value) => value.toFixed(7)).join(',');
+const lobeSum = (lobes) => lobes.map((lobe) => `${lobe[0]}f*spectrumLobe(wavelength,${lobe[1]}f,${lobe[2]}f,${lobe[3]}f)`).join('+');
+export const spectrumWgsl = /* wgsl */ `
+fn spectrumLobe(wavelength:f32,mean:f32,below:f32,above:f32)->f32 {
+let t=(wavelength-mean)/select(above,below,wavelength<mean);
+return exp(-0.5*t*t);
+}
+fn spectrumRgb(wavelength:f32)->vec3f {
+let xyz=vec3f(${XYZ_LOBES.map(lobeSum).join(',')});
+let rgb=mat3x3f(${[0, 1, 2].map((c) => XYZ_TO_RGB.map((row) => `${row[c]}f`).join(',')).join(',')})*xyz;
+return max(rgb,vec3f(0.0))*vec3f(${SCALE.map((s) => `${s}f`).join(',')});
+}
+fn spectrumIorWavelength(wavelength:f32)->f32 {
+if(wavelength>=${IOR_KNEE}.0) { return wavelength; }
+return ${IOR_KNEE}.0-${IOR_KNEE_WIDTH}.0*(1.0-exp((wavelength-${IOR_KNEE}.0)/${IOR_KNEE_WIDTH}.0));
+}
+const SPECTRUM_CDF=array<f32,${(SAMPLE_BINS + 1) * 4}>(${CDF_LITERALS});
+fn spectrumSample(distribution:u32,u:f32,pdf:ptr<function,f32>)->f32 {
+let first=distribution*${SAMPLE_BINS + 1}u;
+var low=0u;
+var high=${SAMPLE_BINS}u;
+while(high-low>1u) {
+let middle=(low+high)/2u;
+if(SPECTRUM_CDF[first+middle]<=u) { low=middle; } else { high=middle; }
+}
+let mass=SPECTRUM_CDF[first+low+1u]-SPECTRUM_CDF[first+low];
+*pdf=mass/${SAMPLE_BIN_WIDTH}.0;
+return ${SPECTRUM_WAVELENGTH_MIN}.0+${SAMPLE_BIN_WIDTH}.0*(f32(low)+(u-SPECTRUM_CDF[first+low])/mass);
+}
+
+fn cU(channel:u32,u:f32)->f32 {
+var pdf=0.0;
+return spectrumIorWavelength(spectrumSample(channel,u,&pdf));
+}
+`;
+export const spectrumGlsl = /* glsl */ `
+const float SPECTRUM_CDF[${(SAMPLE_BINS + 1) * 4}]=float[${(SAMPLE_BINS + 1) * 4}](${CDF_LITERALS});
+
+float cU(int channel,float u) {
+int first=channel*${SAMPLE_BINS + 1};
+int low=0;
+int high=${SAMPLE_BINS};
+while(high-low>1) {
+int middle=(low+high)/2;
+if(SPECTRUM_CDF[first+middle]<=u) { low=middle; } else { high=middle; }
+}
+float start=SPECTRUM_CDF[first+low];
+float wavelength=${SPECTRUM_WAVELENGTH_MIN}.0+${SAMPLE_BIN_WIDTH}.0*(float(low)+(u-start)/(SPECTRUM_CDF[first+low+1]-start));
+if(wavelength>=${IOR_KNEE}.0) { return wavelength; }
+return ${IOR_KNEE}.0-${IOR_KNEE_WIDTH}.0*(1.0-exp((wavelength-${IOR_KNEE}.0)/${IOR_KNEE_WIDTH}.0));
+}
+`;
+//# sourceMappingURL=spectrum.js.map

package/lib/renderer/spectrum.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"spectrum.js","sourceRoot":"","sources":["../../src/renderer/spectrum.ts"],"names":[],"mappings":"AAAA;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,uBAAuB,GAAG,GAAG,CAAC;AAC3C,MAAM,CAAC,MAAM,yBAAyB,GAAG,GAAG,CAAC;AAE7C,4DAA4D;AAC5D,MAAM,SAAS,GAAiB;IAC9B,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACrF,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;IACxD,CAAC,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,EAAE,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;CACzD,CAAC;AACF,8BAA8B;AAC9B,MAAM,UAAU,GAAe;IAC7B,CAAC,MAAM,EAAE,CAAC,MAAM,EAAE,CAAC,MAAM,CAAC;IAC1B,CAAC,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC;IACzB,CAAC,MAAM,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC;CACxB,CAAC;AAEF,SAAS,eAAe,CAAC,UAAkB;IACzC,MAAM,GAAG,GAAG,SAAS,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE;QAClC,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;YACzB,MAAM,CAAC,GAAG,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,UAAU,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAC9E,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC1C,CAAC;QACD,OAAO,GAAG,CAAC;IACb,CAAC,CAAC,CAAC;IACH,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AACnG,CAAC;AAED,SAAS,YAAY;IACnB,MAAM,KAAK,GAAG,IAAI,CAAC;IACnB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,MAAM,GAAG,GAAG,eAAe,CAAC,uBAAuB,GAAG,CAAC,yBAAyB,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACvG,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;YAAE,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;IACvD,CAAC;IACD,0EAA0E;IAC1E,OAAO,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1D,CAAC;AACD,MAAM,KAAK,GAAG,YAAY,EAAE,CAAC;AAE7B,MAAM,UAAU,WAAW,CAAC,UAAkB;IAC5C,OAAO,eAAe,CAAC,UAAU,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;AACzE,CAAC;AAED,iFAAiF;AACjF,oEAAoE;AACpE,MAAM,WAAW,GAAG,EAAE,CAAC;AACvB,MAAM,gBAAgB,GAAG,yBAAyB,GAAG,WAAW,CAAC;AACjE,MAAM,CAAC,MAAM,4BAA4B,GAAG,CAAC,CAAC;AAC9C,SAAS,UAAU;IACjB,MAAM,QAAQ,GAAG,EAAE,CAAC;IACpB,MAAM,IAAI,GAAa,EAAE,CAAC;IAC1B,KAAK,IAAI,YAAY,GAAG,CAAC,EAAE,YAAY,IAAI,4BAA4B,EAAE,YAAY,EAAE,EAAE,CAAC;QACxF,MAAM,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;QAChB,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,WAAW,EAAE,GAAG,EAAE,EAAE,CAAC;YAC3C,IAAI,IAAI,GAAG,CAAC,CAAC;YACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,CAAC,EAAE,EAAE,CAAC;gBAClC,MAAM,GAAG,GAAG,WAAW,CAAC,uBAAuB,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,QAAQ,CAAC,GAAG,gBAAgB,CAAC,CAAC;gBACnG,IAAI,IAAI,YAAY,KAAK,4BAA4B,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,YAAY,CAAC,CAAC;YACvG,CAAC;YACD,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,CAAC;QAC5B,CAAC;QACD,KAAK,MAAM,KAAK,IAAI,GAAG;YAAE,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,KAAK,GAAG,GAAG,CAAC,WAAW,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IACpF,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AACD,MAAM,WAAW,GAAG,UAAU,EAAE,CAAC;AAEjC,6EAA6E;AAC7E,MAAM,UAAU,cAAc,CAAC,YAAoB,EAAE,CAAS;IAC5D,MAAM,KAAK,GAAG,YAAY,GAAG,CAAC,WAAW,GAAG,CAAC,CAAC,CAAC;IAC/C,IAAI,GAAG,GAAG,CAAC,CAAC;IACZ,IAAI,IAAI,GAAG,WAAW,CAAC;IACvB,OAAO,IAAI,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC;QACtB,MAAM,MAAM,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;QACjC,IAAI,WAAW,CAAC,KAAK,GAAG,MAAM,CAAC,IAAI,CAAC;YAAE,GAAG,GAAG,MAAM,CAAC;;YAC9C,IAAI,GAAG,MAAM,CAAC;IACrB,CAAC;IACD,MAAM,IAAI,GAAG,WAAW,CAAC,KAAK,GAAG,GAAG,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC;IACrE,OAAO;QACL,UAAU,EAAE,uBAAuB,GAAG,gBAAgB,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,WAAW,CAAC,KAAK,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC;QACtG,GAAG,EAAE,IAAI,GAAG,gBAAgB;KAC7B,CAAC;AACJ,CAAC;AAED,kFAAkF;AAClF,2EAA2E;AAC3E,oFAAoF;AACpF,MAAM,QAAQ,GAAG,GAAG,CAAC;AACrB,MAAM,cAAc,GAAG,EAAE,CAAC;AAE1B,MAAM,UAAU,qBAAqB,CAAC,UAAkB;IACtD,IAAI,UAAU,IAAI,QAAQ;QAAE,OAAO,UAAU,CAAC;IAC9C,OAAO,QAAQ,GAAG,cAAc,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,UAAU,GAAG,QAAQ,CAAC,GAAG,cAAc,CAAC,CAAC,CAAC;AAC9F,CAAC;AAED,MAAM,YAAY,GAAG,WAAW,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;AAC5E,MAAM,OAAO,GAAG,CAAC,KAAiB,EAAE,EAAE,CACpC,KAAK,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,6BAA6B,IAAI,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;AAE5G,MAAM,CAAC,MAAM,YAAY,GAAG,UAAU,CAAC;;;;;;gBAMrB,SAAS,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;kBAC9B,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,UAAU,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;mCAC9D,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC;;;iBAGrD,QAAQ;SAChB,QAAQ,MAAM,cAAc,2BAA2B,QAAQ,OAAO,cAAc;;+BAEhE,CAAC,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,YAAY;;yBAE1C,WAAW,GAAG,CAAC;;WAE7B,WAAW;;;;;;YAMV,gBAAgB;SACnB,uBAAuB,MAAM,gBAAgB;;;;;;;CAOvD,CAAC;AAEF,MAAM,CAAC,MAAM,YAAY,GAAG,UAAU,CAAC;2BACZ,CAAC,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,WAAW,CAAC,WAAW,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,YAAY;;;oBAG3E,WAAW,GAAG,CAAC;;WAExB,WAAW;;;;;;mBAMH,uBAAuB,MAAM,gBAAgB;iBAC/C,QAAQ;SAChB,QAAQ,MAAM,cAAc,2BAA2B,QAAQ,OAAO,cAAc;;CAE9F,CAAC","sourcesContent":["/**\n * Continuous wavelength sampling shared by camera paths and caustic photons.\n *\n * CIE 1931 matching functions use the multi-lobe fit from Wyman, Sloan and\n * Shirley, \"Simple Analytic Approximations to the CIE XYZ Color Matching\n * Functions\", JCGT 2(2), 2013. Only the published formula and coefficients are\n * used. Out-of-gamut negatives clamp to zero, and each channel is normalized so\n * that a uniform wavelength keeps white light white.\n */\nexport const SPECTRUM_WAVELENGTH_MIN = 400;\nexport const SPECTRUM_WAVELENGTH_RANGE = 300;\n\n// Amplitude, mean, sigma below the mean and sigma above it.\nconst XYZ_LOBES: number[][][] = [\n  [[1.056, 599.8, 37.9, 31.0], [0.362, 442.0, 16.0, 26.7], [-0.065, 501.1, 20.4, 26.2]],\n  [[0.821, 568.8, 46.9, 40.5], [0.286, 530.9, 16.3, 31.1]],\n  [[1.217, 437.0, 11.8, 36.0], [0.681, 459.0, 26.0, 13.8]]\n];\n// Linear sRGB primaries, D65.\nconst XYZ_TO_RGB: number[][] = [\n  [3.2406, -1.5372, -0.4986],\n  [-0.9689, 1.8758, 0.0415],\n  [0.0557, -0.204, 1.057]\n];\n\nfunction unnormalizedRgb(wavelength: number): number[] {\n  const xyz = XYZ_LOBES.map((lobes) => {\n    let sum = 0;\n    for (const lobe of lobes) {\n      const t = (wavelength - lobe[1]) / (wavelength < lobe[1] ? lobe[2] : lobe[3]);\n      sum += lobe[0] * Math.exp(-0.5 * t * t);\n    }\n    return sum;\n  });\n  return XYZ_TO_RGB.map((row) => Math.max(0, row[0] * xyz[0] + row[1] * xyz[1] + row[2] * xyz[2]));\n}\n\nfunction channelScale(): number[] {\n  const steps = 3000;\n  const sum = [0, 0, 0];\n  for (let i = 0; i < steps; i++) {\n    const rgb = unnormalizedRgb(SPECTRUM_WAVELENGTH_MIN + (SPECTRUM_WAVELENGTH_RANGE * (i + 0.5)) / steps);\n    for (let c = 0; c < 3; c++) sum[c] += rgb[c] / steps;\n  }\n  // Rounded so JavaScript and native builds generate identical shader text.\n  return sum.map((mean) => Number((1 / mean).toFixed(6)));\n}\nconst SCALE = channelScale();\n\nexport function spectrumRgb(wavelength: number): number[] {\n  return unnormalizedRgb(wavelength).map((value, c) => value * SCALE[c]);\n}\n\n// Inverse CDFs in 5 nm bins: one per RGB channel for camera paths, and r+g+b for\n// photons. Narrow bins keep each sampled response close to the fit.\nconst SAMPLE_BINS = 60;\nconst SAMPLE_BIN_WIDTH = SPECTRUM_WAVELENGTH_RANGE / SAMPLE_BINS;\nexport const SPECTRUM_PHOTON_DISTRIBUTION = 3;\nfunction sampleCdfs(): number[] {\n  const substeps = 20;\n  const cdfs: number[] = [];\n  for (let distribution = 0; distribution <= SPECTRUM_PHOTON_DISTRIBUTION; distribution++) {\n    const cdf = [0];\n    for (let bin = 0; bin < SAMPLE_BINS; bin++) {\n      let mass = 0;\n      for (let i = 0; i < substeps; i++) {\n        const rgb = spectrumRgb(SPECTRUM_WAVELENGTH_MIN + (bin + (i + 0.5) / substeps) * SAMPLE_BIN_WIDTH);\n        mass += distribution === SPECTRUM_PHOTON_DISTRIBUTION ? rgb[0] + rgb[1] + rgb[2] : rgb[distribution];\n      }\n      cdf.push(cdf[bin] + mass);\n    }\n    for (const value of cdf) cdfs.push(Number((value / cdf[SAMPLE_BINS]).toFixed(7)));\n  }\n  return cdfs;\n}\nconst SAMPLE_CDFS = sampleCdfs();\n\n/** Maps a unit random number to a wavelength and its probability density. */\nexport function spectrumSample(distribution: number, u: number) {\n  const first = distribution * (SAMPLE_BINS + 1);\n  let low = 0;\n  let high = SAMPLE_BINS;\n  while (high - low > 1) {\n    const middle = (low + high) >> 1;\n    if (SAMPLE_CDFS[first + middle] <= u) low = middle;\n    else high = middle;\n  }\n  const mass = SAMPLE_CDFS[first + low + 1] - SAMPLE_CDFS[first + low];\n  return {\n    wavelength: SPECTRUM_WAVELENGTH_MIN + SAMPLE_BIN_WIDTH * (low + (u - SAMPLE_CDFS[first + low]) / mass),\n    pdf: mass / SAMPLE_BIN_WIDTH\n  };\n}\n\n// Cauchy IOR grows quickly at short wavelengths and strong dispersion scales turn\n// violet light into total internal reflection. Below 500 nm the refracting\n// wavelength eases toward 475 nm, while color still follows the sampled wavelength.\nconst IOR_KNEE = 500;\nconst IOR_KNEE_WIDTH = 25;\n\nexport function spectrumIorWavelength(wavelength: number) {\n  if (wavelength >= IOR_KNEE) return wavelength;\n  return IOR_KNEE - IOR_KNEE_WIDTH * (1 - Math.exp((wavelength - IOR_KNEE) / IOR_KNEE_WIDTH));\n}\n\nconst CDF_LITERALS = SAMPLE_CDFS.map((value) => value.toFixed(7)).join(',');\nconst lobeSum = (lobes: number[][]) =>\n  lobes.map((lobe) => `${lobe[0]}f*spectrumLobe(wavelength,${lobe[1]}f,${lobe[2]}f,${lobe[3]}f)`).join('+');\n\nexport const spectrumWgsl = /* wgsl */ `\nfn spectrumLobe(wavelength:f32,mean:f32,below:f32,above:f32)->f32 {\n  let t=(wavelength-mean)/select(above,below,wavelength<mean);\n  return exp(-0.5*t*t);\n}\nfn spectrumRgb(wavelength:f32)->vec3f {\n  let xyz=vec3f(${XYZ_LOBES.map(lobeSum).join(',')});\n  let rgb=mat3x3f(${[0, 1, 2].map((c) => XYZ_TO_RGB.map((row) => `${row[c]}f`).join(',')).join(',')})*xyz;\n  return max(rgb,vec3f(0.0))*vec3f(${SCALE.map((s) => `${s}f`).join(',')});\n}\nfn spectrumIorWavelength(wavelength:f32)->f32 {\n  if(wavelength>=${IOR_KNEE}.0) { return wavelength; }\n  return ${IOR_KNEE}.0-${IOR_KNEE_WIDTH}.0*(1.0-exp((wavelength-${IOR_KNEE}.0)/${IOR_KNEE_WIDTH}.0));\n}\nconst SPECTRUM_CDF=array<f32,${(SAMPLE_BINS + 1) * 4}>(${CDF_LITERALS});\nfn spectrumSample(distribution:u32,u:f32,pdf:ptr<function,f32>)->f32 {\n  let first=distribution*${SAMPLE_BINS + 1}u;\n  var low=0u;\n  var high=${SAMPLE_BINS}u;\n  while(high-low>1u) {\n    let middle=(low+high)/2u;\n    if(SPECTRUM_CDF[first+middle]<=u) { low=middle; } else { high=middle; }\n  }\n  let mass=SPECTRUM_CDF[first+low+1u]-SPECTRUM_CDF[first+low];\n  *pdf=mass/${SAMPLE_BIN_WIDTH}.0;\n  return ${SPECTRUM_WAVELENGTH_MIN}.0+${SAMPLE_BIN_WIDTH}.0*(f32(low)+(u-SPECTRUM_CDF[first+low])/mass);\n}\n// Refracting wavelength of a camera path that carries one RGB channel.\nfn spectrumCameraWavelength(channel:u32,u:f32)->f32 {\n  var pdf=0.0;\n  return spectrumIorWavelength(spectrumSample(channel,u,&pdf));\n}\n`;\n\nexport const spectrumGlsl = /* glsl */ `\nconst float SPECTRUM_CDF[${(SAMPLE_BINS + 1) * 4}]=float[${(SAMPLE_BINS + 1) * 4}](${CDF_LITERALS});\n// Refracting wavelength of a camera path that carries one RGB channel.\nfloat spectrumCameraWavelength(int channel,float u) {\n  int first=channel*${SAMPLE_BINS + 1};\n  int low=0;\n  int high=${SAMPLE_BINS};\n  while(high-low>1) {\n    int middle=(low+high)/2;\n    if(SPECTRUM_CDF[first+middle]<=u) { low=middle; } else { high=middle; }\n  }\n  float start=SPECTRUM_CDF[first+low];\n  float wavelength=${SPECTRUM_WAVELENGTH_MIN}.0+${SAMPLE_BIN_WIDTH}.0*(float(low)+(u-start)/(SPECTRUM_CDF[first+low+1]-start));\n  if(wavelength>=${IOR_KNEE}.0) { return wavelength; }\n  return ${IOR_KNEE}.0-${IOR_KNEE_WIDTH}.0*(1.0-exp((wavelength-${IOR_KNEE}.0)/${IOR_KNEE_WIDTH}.0));\n}\n`;\n"]}

package/lib/renderer/textureAtlas.d.ts

@@ -0,0 +1,60 @@
+import type TextureState from '../scene/TextureState.js';
+export declare const DEFAULT_MAX_TEXTURE_ARRAY_BYTES: number;
+export interface MaterialTableOptions {
+    maxTextureArrayBytes?: number;
+    maxTextureArrayLayers?: number;
+    maxTextureDimension2D?: number;
+}
+export interface TextureAtlasStats {
+    estimatedBytes: number;
+    height: number;
+    layers: number;
+    packTimeMs: number;
+    scale: number;
+    sourcePixels: number;
+    textureCount: number;
+    width: number;
+}
+export interface AtlasPlacement {
+    height: number;
+    index: number;
+    layer: number;
+    width: number;
+    x: number;
+    y: number;
+}
+interface AtlasShelf {
+    height: number;
+    x: number;
+    y: number;
+}
+interface AtlasPage {
+    nextY: number;
+    placements: AtlasPlacement[];
+    shelves: AtlasShelf[];
+}
+export declare function planTextureAtlas(textures: TextureState[], mapColorLinear: number[], options?: MaterialTableOptions, floatStorage?: boolean): {
+    pages: AtlasPage[];
+    placements: AtlasPlacement[];
+    width: number;
+    height: number;
+    floatStorage: boolean;
+    bytesPerPixel: number;
+    mapInfo: number[];
+    mapTransform: number[];
+    mapWrap: number[];
+    mapUvChannel: number[];
+    mapColorLinear: number[];
+    stats: {
+        estimatedBytes: number;
+        height: number;
+        layers: number;
+        packTimeMs: number;
+        scale: number;
+        sourcePixels: number;
+        textureCount: number;
+        width: number;
+    };
+};
+export type TextureAtlasPlan = ReturnType<typeof planTextureAtlas>;
+export {};

package/lib/renderer/textureAtlas.js

@@ -0,0 +1,145 @@
+import { textureWrapIndex } from '../scene/TextureState.js';
+export const DEFAULT_MAX_TEXTURE_ARRAY_BYTES = 384 * 1024 * 1024;
+/**
+ * Packs material textures into a small 2D-array atlas.
+ *
+ * The input rectangles are sorted from tallest/largest to smallest, then
+ * placed with a shelf allocator. Each atlas layer contains horizontal shelves;
+ * a texture is inserted into the first shelf with enough remaining width, or a
+ * new shelf is opened below the existing ones. A new array layer is allocated
+ * only when none of the current layers has enough room.
+ *
+ * `makeTextureArray` first tries this packing at full resolution. If the result
+ * exceeds the adapter layer limit or the configured memory budget, it performs
+ * a binary search for the largest uniform scale that fits. This preserves full
+ * resolution whenever packing alone is sufficient and degrades predictably on
+ * constrained hardware.
+ *
+ * The returned placements keep the original logical texture indices. Shaders
+ * use per-texture atlas metadata instead of treating the logical index as the
+ * physical array layer.
+ */
+function packTextureAtlas(textures, pageSize, scale) {
+    const pages = [];
+    const placements = new Array(textures.length);
+    const items = textures
+        .map((texture, index) => ({
+        height: Math.max(1, Math.floor(texture.image.height * scale)),
+        index,
+        width: Math.max(1, Math.floor(texture.image.width * scale))
+    }))
+        .sort((a, b) => b.height - a.height || b.width - a.width);
+    for (const item of items) {
+        let placement;
+        for (let layer = 0; layer < pages.length && !placement; layer++) {
+            const page = pages[layer];
+            for (const shelf of page.shelves) {
+                if (item.height <= shelf.height &&
+                    shelf.x + item.width <= pageSize) {
+                    placement = {
+                        height: item.height, index: item.index, width: item.width,
+                        layer,
+                        x: shelf.x,
+                        y: shelf.y
+                    };
+                    shelf.x += item.width;
+                    break;
+                }
+            }
+            if (!placement && page.nextY + item.height <= pageSize) {
+                placement = {
+                    height: item.height, index: item.index, width: item.width,
+                    layer,
+                    x: 0,
+                    y: page.nextY
+                };
+                page.shelves.push({
+                    height: item.height,
+                    x: item.width,
+                    y: page.nextY
+                });
+                page.nextY += item.height;
+            }
+        }
+        if (!placement) {
+            const layer = pages.length;
+            placement = { height: item.height, index: item.index, width: item.width, layer, x: 0, y: 0 };
+            pages.push({
+                nextY: item.height,
+                placements: [],
+                shelves: [{ height: item.height, x: item.width, y: 0 }]
+            });
+        }
+        pages[placement.layer].placements.push(placement);
+        placements[item.index] = placement;
+    }
+    return { pages, placements };
+}
+export function planTextureAtlas(textures, mapColorLinear, options = {}, floatStorage = mapColorLinear.some(value => !!value)) {
+    const bytesPerPixel = floatStorage ? 8 : 4;
+    if (!textures.length)
+        return { pages: [], placements: [],
+            width: 1, height: 1, floatStorage, bytesPerPixel, mapInfo: [],
+            mapTransform: [], mapWrap: [], mapUvChannel: [], mapColorLinear,
+            stats: { estimatedBytes: 0, height: 1, layers: 0, packTimeMs: 0, scale: 1,
+                sourcePixels: 0, textureCount: 0, width: 1 } };
+    const maxTextureArrayLayers = options.maxTextureArrayLayers ?? Infinity;
+    const maxTextureDimension2D = options.maxTextureDimension2D ?? Infinity;
+    const maxTextureArrayBytes = options.maxTextureArrayBytes ?? DEFAULT_MAX_TEXTURE_ARRAY_BYTES;
+    const sourcePixels = textures.reduce((sum, texture) => sum + texture.image.width * texture.image.height, 0);
+    const maxSourceDimension = textures.reduce((max, texture) => Math.max(max, texture.image.width, texture.image.height), 1);
+    const maxPageSizeFromBudget = Math.max(1, Math.floor(Math.sqrt(maxTextureArrayBytes / bytesPerPixel)));
+    const pageSize = Math.max(1, Math.floor(Math.min(maxSourceDimension, maxTextureDimension2D, maxPageSizeFromBudget)));
+    const pageBytes = pageSize * pageSize * bytesPerPixel;
+    const maxPageCount = Math.max(1, Math.min(maxTextureArrayLayers, Math.floor(maxTextureArrayBytes / pageBytes)));
+    const maxScale = Math.min(1, pageSize / maxSourceDimension);
+    let packed = packTextureAtlas(textures, pageSize, maxScale);
+    let scale = maxScale;
+    if (packed.pages.length > maxPageCount) {
+        let low = 0;
+        let high = maxScale;
+        let bestPacked;
+        let bestScale = 0;
+        for (let i = 0; i < 14; i++) {
+            const candidateScale = (low + high) / 2;
+            const candidate = packTextureAtlas(textures, pageSize, candidateScale);
+            if (candidate.pages.length <= maxPageCount) {
+                low = candidateScale;
+                bestScale = candidateScale;
+                bestPacked = candidate;
+            }
+            else {
+                high = candidateScale;
+            }
+        }
+        if (!bestPacked) {
+            throw new Error(`Unable to fit ${textures.length} textures into ${maxPageCount} atlas layers.`);
+        }
+        packed = bestPacked;
+        scale = bestScale;
+    }
+    const atlasWidth = packed.placements.reduce((max, placement) => Math.max(max, placement.x + placement.width), 1);
+    const atlasHeight = packed.placements.reduce((max, placement) => Math.max(max, placement.y + placement.height), 1);
+    // mapInfo = [uvScale.x, uvScale.y, layer + uvOffset.x, uvOffset.y].
+    // Packing layer and X offset into one float keeps the shader metadata to one
+    // vec4 per logical texture; floor(z) restores the layer and fract(z) the X.
+    const mapInfo = packed.placements.flatMap((placement) => {
+        const offsetX = (placement.x + 0.5) / atlasWidth;
+        const offsetY = (placement.y + 0.5) / atlasHeight;
+        return [
+            Math.max(placement.width - 1, 0) / atlasWidth,
+            Math.max(placement.height - 1, 0) / atlasHeight,
+            placement.layer + offsetX,
+            offsetY
+        ];
+    });
+    return { pages: packed.pages, placements: packed.placements, width: atlasWidth, height: atlasHeight,
+        floatStorage, bytesPerPixel, mapInfo,
+        mapTransform: textures.flatMap(texture => [texture.uvScale[0], texture.uvScale[1], texture.uvOffset[0], texture.uvOffset[1]]),
+        mapWrap: textures.flatMap(texture => [textureWrapIndex(texture.wrapS), textureWrapIndex(texture.wrapT)]),
+        mapUvChannel: textures.map(texture => texture.uvChannel), mapColorLinear,
+        stats: { estimatedBytes: atlasWidth * atlasHeight * packed.pages.length * bytesPerPixel,
+            height: atlasHeight, layers: packed.pages.length, packTimeMs: 0, scale, sourcePixels,
+            textureCount: textures.length, width: atlasWidth } };
+}
+//# sourceMappingURL=textureAtlas.js.map

package/lib/renderer/textureAtlas.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"textureAtlas.js","sourceRoot":"","sources":["../../src/renderer/textureAtlas.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,gBAAgB,EAAE,MAAM,0BAA0B,CAAC;AAE5D,MAAM,CAAC,MAAM,+BAA+B,GAAG,GAAG,GAAG,IAAI,GAAG,IAAI,CAAC;AAwCjE;;;;;;;;;;;;;;;;;;GAkBG;AACH,SAAS,gBAAgB,CACvB,QAAwB,EACxB,QAAgB,EAChB,KAAa;IAEb,MAAM,KAAK,GAAgB,EAAE,CAAC;IAC9B,MAAM,UAAU,GAAG,IAAI,KAAK,CAAiB,QAAQ,CAAC,MAAM,CAAC,CAAC;IAC9D,MAAM,KAAK,GAAG,QAAQ;SACnB,GAAG,CAAC,CAAC,OAAO,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;QACxB,MAAM,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC;QAC7D,KAAK;QACL,KAAK,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC;KAC5D,CAAC,CAAC;SACF,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,MAAM,IAAI,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;IAE5D,KAAK,MAAM,IAAI,IAAI,KAAK,EAAE,CAAC;QACzB,IAAI,SAAqC,CAAC;QAC1C,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,SAAS,EAAE,KAAK,EAAE,EAAE,CAAC;YAChE,MAAM,IAAI,GAAG,KAAK,CAAC,KAAK,CAAC,CAAC;YAC1B,KAAK,MAAM,KAAK,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC;gBACjC,IACE,IAAI,CAAC,MAAM,IAAI,KAAK,CAAC,MAAM;oBAC3B,KAAK,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,IAAI,QAAQ,EAChC,CAAC;oBACD,SAAS,GAAG;wBACV,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK;wBACzD,KAAK;wBACL,CAAC,EAAE,KAAK,CAAC,CAAC;wBACV,CAAC,EAAE,KAAK,CAAC,CAAC;qBACX,CAAC;oBACF,KAAK,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC;oBACtB,MAAM;gBACR,CAAC;YACH,CAAC;YACD,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,IAAI,QAAQ,EAAE,CAAC;gBACvD,SAAS,GAAG;oBACV,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK;oBACzD,KAAK;oBACL,CAAC,EAAE,CAAC;oBACJ,CAAC,EAAE,IAAI,CAAC,KAAK;iBACd,CAAC;gBACF,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC;oBAChB,MAAM,EAAE,IAAI,CAAC,MAAM;oBACnB,CAAC,EAAE,IAAI,CAAC,KAAK;oBACb,CAAC,EAAE,IAAI,CAAC,KAAK;iBACd,CAAC,CAAC;gBACH,IAAI,CAAC,KAAK,IAAI,IAAI,CAAC,MAAM,CAAC;YAC5B,CAAC;QACH,CAAC;QAED,IAAI,CAAC,SAAS,EAAE,CAAC;YACf,MAAM,KAAK,GAAG,KAAK,CAAC,MAAM,CAAC;YAC3B,SAAS,GAAG,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,IAAI,CAAC,KAAK,EAAE,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;YAC7F,KAAK,CAAC,IAAI,CAAC;gBACT,KAAK,EAAE,IAAI,CAAC,MAAM;gBAClB,UAAU,EAAE,EAAE;gBACd,OAAO,EAAE,CAAC,EAAE,MAAM,EAAE,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE,IAAI,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC;aACxD,CAAC,CAAC;QACL,CAAC;QACD,KAAK,CAAC,SAAS,CAAC,KAAK,CAAC,CAAC,UAAU,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC;QAClD,UAAU,CAAC,IAAI,CAAC,KAAK,CAAC,GAAG,SAAS,CAAC;IACrC,CAAC;IAED,OAAO,EAAE,KAAK,EAAE,UAAU,EAAE,CAAC;AAC/B,CAAC;AAGD,MAAM,UAAU,gBAAgB,CAAC,QAAwB,EAAE,cAAwB,EAAE,UAAgC,EAAE,EACrH,YAAY,GAAG,cAAc,CAAC,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC,KAAK,CAAC;IACpD,MAAM,aAAa,GAAG,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC3C,IAAI,CAAC,QAAQ,CAAC,MAAM;QAAE,OAAO,EAAE,KAAK,EAAE,EAAiB,EAAE,UAAU,EAAE,EAAsB;YACzF,KAAK,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,aAAa,EAAE,OAAO,EAAE,EAAc;YACzE,YAAY,EAAE,EAAc,EAAE,OAAO,EAAE,EAAc,EAAE,YAAY,EAAE,EAAc,EAAE,cAAc;YACnG,KAAK,EAAE,EAAC,cAAc,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,UAAU,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC;gBACtE,YAAY,EAAE,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EAAC,EAAE,CAAC;IAClD,MAAM,qBAAqB,GAAG,OAAO,CAAC,qBAAqB,IAAI,QAAQ,CAAC;IACxE,MAAM,qBAAqB,GAAG,OAAO,CAAC,qBAAqB,IAAI,QAAQ,CAAC;IACxE,MAAM,oBAAoB,GACxB,OAAO,CAAC,oBAAoB,IAAI,+BAA+B,CAAC;IAClE,MAAM,YAAY,GAAG,QAAQ,CAAC,MAAM,CAClC,CAAC,GAAG,EAAE,OAAO,EAAE,EAAE,CAAC,GAAG,GAAG,OAAO,CAAC,KAAK,CAAC,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC,MAAM,EAClE,CAAC,CACF,CAAC;IACF,MAAM,kBAAkB,GAAG,QAAQ,CAAC,MAAM,CACxC,CAAC,GAAG,EAAE,OAAO,EAAE,EAAE,CACf,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,OAAO,CAAC,KAAK,CAAC,KAAK,EAAE,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,EAC1D,CAAC,CACF,CAAC;IACF,MAAM,qBAAqB,GAAG,IAAI,CAAC,GAAG,CACpC,CAAC,EACD,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,oBAAoB,GAAG,aAAa,CAAC,CAAC,CAC5D,CAAC;IACF,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CACvB,CAAC,EACD,IAAI,CAAC,KAAK,CACR,IAAI,CAAC,GAAG,CACN,kBAAkB,EAClB,qBAAqB,EACrB,qBAAqB,CACtB,CACF,CACF,CAAC;IACF,MAAM,SAAS,GAAG,QAAQ,GAAG,QAAQ,GAAG,aAAa,CAAC;IACtD,MAAM,YAAY,GAAG,IAAI,CAAC,GAAG,CAC3B,CAAC,EACD,IAAI,CAAC,GAAG,CACN,qBAAqB,EACrB,IAAI,CAAC,KAAK,CAAC,oBAAoB,GAAG,SAAS,CAAC,CAC7C,CACF,CAAC;IACF,MAAM,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CAAC;IAC5D,IAAI,MAAM,GAAG,gBAAgB,CAAC,QAAQ,EAAE,QAAQ,EAAE,QAAQ,CAAC,CAAC;IAC5D,IAAI,KAAK,GAAG,QAAQ,CAAC;IAErB,IAAI,MAAM,CAAC,KAAK,CAAC,MAAM,GAAG,YAAY,EAAE,CAAC;QACvC,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,IAAI,GAAG,QAAQ,CAAC;QACpB,IAAI,UAA2D,CAAC;QAChE,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;YAC5B,MAAM,cAAc,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC;YACxC,MAAM,SAAS,GAAG,gBAAgB,CAChC,QAAQ,EACR,QAAQ,EACR,cAAc,CACf,CAAC;YACF,IAAI,SAAS,CAAC,KAAK,CAAC,MAAM,IAAI,YAAY,EAAE,CAAC;gBAC3C,GAAG,GAAG,cAAc,CAAC;gBACrB,SAAS,GAAG,cAAc,CAAC;gBAC3B,UAAU,GAAG,SAAS,CAAC;YACzB,CAAC;iBAAM,CAAC;gBACN,IAAI,GAAG,cAAc,CAAC;YACxB,CAAC;QACH,CAAC;QACD,IAAI,CAAC,UAAU,EAAE,CAAC;YAChB,MAAM,IAAI,KAAK,CACb,iBAAiB,QAAQ,CAAC,MAAM,kBAAkB,YAAY,gBAAgB,CAC/E,CAAC;QACJ,CAAC;QACD,MAAM,GAAG,UAAU,CAAC;QACpB,KAAK,GAAG,SAAS,CAAC;IACpB,CAAC;IAED,MAAM,UAAU,GAAG,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,EAAE,CAAC,CAAC,CAAC;IACjH,MAAM,WAAW,GAAG,MAAM,CAAC,UAAU,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,SAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,SAAS,CAAC,CAAC,GAAG,SAAS,CAAC,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC;IACnH,oEAAoE;IACpE,6EAA6E;IAC7E,4EAA4E;IAC5E,MAAM,OAAO,GAAG,MAAM,CAAC,UAAU,CAAC,OAAO,CAAC,CAAC,SAAS,EAAE,EAAE;QACtD,MAAM,OAAO,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC;QACjD,MAAM,OAAO,GAAG,CAAC,SAAS,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,WAAW,CAAC;QAClD,OAAO;YACL,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,UAAU;YAC7C,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,WAAW;YAC/C,SAAS,CAAC,KAAK,GAAG,OAAO;YACzB,OAAO;SACR,CAAC;IACJ,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,KAAK,EAAE,MAAM,CAAC,KAAK,EAAE,UAAU,EAAE,MAAM,CAAC,UAAU,EAAE,KAAK,EAAE,UAAU,EAAE,MAAM,EAAE,WAAW;QACjG,YAAY,EAAE,aAAa,EAAE,OAAO;QACpC,YAAY,EAAE,QAAQ,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,OAAO,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7H,OAAO,EAAE,QAAQ,CAAC,OAAO,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,gBAAgB,CAAC,OAAO,CAAC,KAAK,CAAC,EAAE,gBAAgB,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC;QACxG,YAAY,EAAE,QAAQ,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,OAAO,CAAC,SAAS,CAAC,EAAE,cAAc;QACxE,KAAK,EAAE,EAAE,cAAc,EAAE,UAAU,GAAG,WAAW,GAAG,MAAM,CAAC,KAAK,CAAC,MAAM,GAAG,aAAa;YACrF,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,MAAM,CAAC,KAAK,CAAC,MAAM,EAAE,UAAU,EAAE,CAAC,EAAE,KAAK,EAAE,YAAY;YACpF,YAAY,EAAE,QAAQ,CAAC,MAAM,EAAE,KAAK,EAAE,UAAU,EAAE,EAAE,CAAC;AAC3D,CAAC","sourcesContent":["import type TextureState from '../scene/TextureState.js';\nimport { textureWrapIndex } from '../scene/TextureState.js';\n\nexport const DEFAULT_MAX_TEXTURE_ARRAY_BYTES = 384 * 1024 * 1024;\n\nexport interface MaterialTableOptions {\n  maxTextureArrayBytes?: number;\n  maxTextureArrayLayers?: number;\n  maxTextureDimension2D?: number;\n}\n\nexport interface TextureAtlasStats {\n  estimatedBytes: number;\n  height: number;\n  layers: number;\n  packTimeMs: number;\n  scale: number;\n  sourcePixels: number;\n  textureCount: number;\n  width: number;\n}\n\nexport interface AtlasPlacement {\n  height: number;\n  index: number;\n  layer: number;\n  width: number;\n  x: number;\n  y: number;\n}\n\ninterface AtlasShelf {\n  height: number;\n  x: number;\n  y: number;\n}\n\ninterface AtlasPage {\n  nextY: number;\n  placements: AtlasPlacement[];\n  shelves: AtlasShelf[];\n}\n\n/**\n * Packs material textures into a small 2D-array atlas.\n *\n * The input rectangles are sorted from tallest/largest to smallest, then\n * placed with a shelf allocator. Each atlas layer contains horizontal shelves;\n * a texture is inserted into the first shelf with enough remaining width, or a\n * new shelf is opened below the existing ones. A new array layer is allocated\n * only when none of the current layers has enough room.\n *\n * `makeTextureArray` first tries this packing at full resolution. If the result\n * exceeds the adapter layer limit or the configured memory budget, it performs\n * a binary search for the largest uniform scale that fits. This preserves full\n * resolution whenever packing alone is sufficient and degrades predictably on\n * constrained hardware.\n *\n * The returned placements keep the original logical texture indices. Shaders\n * use per-texture atlas metadata instead of treating the logical index as the\n * physical array layer.\n */\nfunction packTextureAtlas(\n  textures: TextureState[],\n  pageSize: number,\n  scale: number\n) {\n  const pages: AtlasPage[] = [];\n  const placements = new Array<AtlasPlacement>(textures.length);\n  const items = textures\n    .map((texture, index) => ({\n      height: Math.max(1, Math.floor(texture.image.height * scale)),\n      index,\n      width: Math.max(1, Math.floor(texture.image.width * scale))\n    }))\n    .sort((a, b) => b.height - a.height || b.width - a.width);\n\n  for (const item of items) {\n    let placement: AtlasPlacement | undefined;\n    for (let layer = 0; layer < pages.length && !placement; layer++) {\n      const page = pages[layer];\n      for (const shelf of page.shelves) {\n        if (\n          item.height <= shelf.height &&\n          shelf.x + item.width <= pageSize\n        ) {\n          placement = {\n            height: item.height, index: item.index, width: item.width,\n            layer,\n            x: shelf.x,\n            y: shelf.y\n          };\n          shelf.x += item.width;\n          break;\n        }\n      }\n      if (!placement && page.nextY + item.height <= pageSize) {\n        placement = {\n          height: item.height, index: item.index, width: item.width,\n          layer,\n          x: 0,\n          y: page.nextY\n        };\n        page.shelves.push({\n          height: item.height,\n          x: item.width,\n          y: page.nextY\n        });\n        page.nextY += item.height;\n      }\n    }\n\n    if (!placement) {\n      const layer = pages.length;\n      placement = { height: item.height, index: item.index, width: item.width, layer, x: 0, y: 0 };\n      pages.push({\n        nextY: item.height,\n        placements: [],\n        shelves: [{ height: item.height, x: item.width, y: 0 }]\n      });\n    }\n    pages[placement.layer].placements.push(placement);\n    placements[item.index] = placement;\n  }\n\n  return { pages, placements };\n}\n\n\nexport function planTextureAtlas(textures: TextureState[], mapColorLinear: number[], options: MaterialTableOptions = {},\n  floatStorage = mapColorLinear.some(value => !!value)) {\n  const bytesPerPixel = floatStorage ? 8 : 4;\n  if (!textures.length) return { pages: [] as AtlasPage[], placements: [] as AtlasPlacement[],\n    width: 1, height: 1, floatStorage, bytesPerPixel, mapInfo: [] as number[],\n    mapTransform: [] as number[], mapWrap: [] as number[], mapUvChannel: [] as number[], mapColorLinear,\n    stats: {estimatedBytes: 0, height: 1, layers: 0, packTimeMs: 0, scale: 1,\n      sourcePixels: 0, textureCount: 0, width: 1} };\n  const maxTextureArrayLayers = options.maxTextureArrayLayers ?? Infinity;\n  const maxTextureDimension2D = options.maxTextureDimension2D ?? Infinity;\n  const maxTextureArrayBytes =\n    options.maxTextureArrayBytes ?? DEFAULT_MAX_TEXTURE_ARRAY_BYTES;\n  const sourcePixels = textures.reduce(\n    (sum, texture) => sum + texture.image.width * texture.image.height,\n    0\n  );\n  const maxSourceDimension = textures.reduce(\n    (max, texture) =>\n      Math.max(max, texture.image.width, texture.image.height),\n    1\n  );\n  const maxPageSizeFromBudget = Math.max(\n    1,\n    Math.floor(Math.sqrt(maxTextureArrayBytes / bytesPerPixel))\n  );\n  const pageSize = Math.max(\n    1,\n    Math.floor(\n      Math.min(\n        maxSourceDimension,\n        maxTextureDimension2D,\n        maxPageSizeFromBudget\n      )\n    )\n  );\n  const pageBytes = pageSize * pageSize * bytesPerPixel;\n  const maxPageCount = Math.max(\n    1,\n    Math.min(\n      maxTextureArrayLayers,\n      Math.floor(maxTextureArrayBytes / pageBytes)\n    )\n  );\n  const maxScale = Math.min(1, pageSize / maxSourceDimension);\n  let packed = packTextureAtlas(textures, pageSize, maxScale);\n  let scale = maxScale;\n\n  if (packed.pages.length > maxPageCount) {\n    let low = 0;\n    let high = maxScale;\n    let bestPacked: ReturnType<typeof packTextureAtlas> | undefined;\n    let bestScale = 0;\n    for (let i = 0; i < 14; i++) {\n      const candidateScale = (low + high) / 2;\n      const candidate = packTextureAtlas(\n        textures,\n        pageSize,\n        candidateScale\n      );\n      if (candidate.pages.length <= maxPageCount) {\n        low = candidateScale;\n        bestScale = candidateScale;\n        bestPacked = candidate;\n      } else {\n        high = candidateScale;\n      }\n    }\n    if (!bestPacked) {\n      throw new Error(\n        `Unable to fit ${textures.length} textures into ${maxPageCount} atlas layers.`\n      );\n    }\n    packed = bestPacked;\n    scale = bestScale;\n  }\n\n  const atlasWidth = packed.placements.reduce((max, placement) => Math.max(max, placement.x + placement.width), 1);\n  const atlasHeight = packed.placements.reduce((max, placement) => Math.max(max, placement.y + placement.height), 1);\n  // mapInfo = [uvScale.x, uvScale.y, layer + uvOffset.x, uvOffset.y].\n  // Packing layer and X offset into one float keeps the shader metadata to one\n  // vec4 per logical texture; floor(z) restores the layer and fract(z) the X.\n  const mapInfo = packed.placements.flatMap((placement) => {\n    const offsetX = (placement.x + 0.5) / atlasWidth;\n    const offsetY = (placement.y + 0.5) / atlasHeight;\n    return [\n      Math.max(placement.width - 1, 0) / atlasWidth,\n      Math.max(placement.height - 1, 0) / atlasHeight,\n      placement.layer + offsetX,\n      offsetY\n    ];\n  });\n\n  return { pages: packed.pages, placements: packed.placements, width: atlasWidth, height: atlasHeight,\n    floatStorage, bytesPerPixel, mapInfo,\n    mapTransform: textures.flatMap(texture => [texture.uvScale[0], texture.uvScale[1], texture.uvOffset[0], texture.uvOffset[1]]),\n    mapWrap: textures.flatMap(texture => [textureWrapIndex(texture.wrapS), textureWrapIndex(texture.wrapT)]),\n    mapUvChannel: textures.map(texture => texture.uvChannel), mapColorLinear,\n    stats: { estimatedBytes: atlasWidth * atlasHeight * packed.pages.length * bytesPerPixel,\n      height: atlasHeight, layers: packed.pages.length, packTimeMs: 0, scale, sourcePixels,\n      textureCount: textures.length, width: atlasWidth } };\n}\nexport type TextureAtlasPlan = ReturnType<typeof planTextureAtlas>;\n"]}

package/lib/renderer/textureAtlasPixels.d.ts

@@ -0,0 +1,7 @@
+import type { TextureDataSource } from '../scene/textureData.js';
+import type { AtlasPlacement } from './textureAtlas.js';
+export declare function writeAtlasPixel(data: Uint16Array | Uint8Array, index: number, value: number): void;
+export declare function writeRawAtlasPlacement(data: Uint16Array | Uint8Array, width: number, placement: AtlasPlacement, texture: {
+    image: TextureDataSource;
+    flipY: boolean;
+}, floatStorage: boolean): void;