@tinypoly/gpu-renderer 0.0.4 → 0.1.0

package/dist/subsurface.d.ts

@@ -0,0 +1,76 @@
+import type { SerializedMaterial } from "./upstream.js";
+/**
+ * How a material scatters light inside itself, from `userData.cpuRenderer.subsurface`. `true` takes every default.
+ *
+ * The material's own colour (with its map) is what the medium reflects in all once light has wandered through it, as
+ * in Blender's Principled BSDF, so a textured skin keeps its texture; `radius` says how far each colour gets before it
+ * scatters, which is what makes red bleed through an ear and blue stay at the surface.
+ */
+export type SubsurfaceInput = boolean | {
+    /** Share of the diffuse layer that scatters inside rather than at the surface, `0` to `1`. Defaults to `1`. */
+    weight?: number;
+    /**
+     * Mean free path of red, green and blue, as multiples of `scale`; one number for all three. Defaults to
+     * `[1, 0.2, 0.1]`, Blender's, which suits skin.
+     */
+    radius?: number | readonly number[];
+    /** World units `radius` is measured in. Defaults to `0.05`: five centimetres in a scene modelled in metres. */
+    scale?: number;
+    /** Henyey-Greenstein anisotropy of each scattering, `-1` (back) to `1` (forward). Defaults to `0`. */
+    anisotropy?: number;
+};
+/** A medium resolved to what the kernel reads. */
+export interface SubsurfaceMedium {
+    weight: number;
+    /** Mean free path per channel, in world units. */
+    radius: [number, number, number];
+    anisotropy: number;
+}
+/** A snapshot material, with the medium `serializeLive` attaches; a snapshot built elsewhere can carry one too. */
+export type SubsurfaceMaterial = SerializedMaterial & {
+    subsurface?: SubsurfaceInput;
+};
+/** The medium an input describes, with its defaults, or `null` when it describes none. */
+export declare function subsurfaceMedium(input: unknown): SubsurfaceMedium | null;
+/** The medium a snapshot material carries, if any. */
+export declare function subsurfaceOf(material: SerializedMaterial | undefined): SubsurfaceMedium | null;
+/**
+ * A mesh filled with a medium, from `userData.cpuRenderer.volume`: its surface is only the boundary, invisible, and
+ * light scatters and fades in the space it encloses. The mesh has to be closed, with its normals pointing out.
+ */
+export type VolumeInput = {
+    /** How much of the light the medium takes per scene unit: one over the mean distance a ray gets. Defaults to `1`. */
+    density?: number;
+    /** What it keeps at each scattering, `#rrggbb` or linear rgb: white for fog and clouds, dark for soot. */
+    color?: string | readonly number[];
+    /** Henyey-Greenstein anisotropy, `-1` (back) to `1` (forward). Defaults to `0`; haze and dust are forward, ~0.7. */
+    anisotropy?: number;
+};
+export interface VolumeMedium {
+    density: number;
+    /** Single scattering albedo, linear rgb. */
+    albedo: [number, number, number];
+    anisotropy: number;
+}
+/** A snapshot material that may carry a volume, as `serializeLive` attaches it. */
+export type VolumeMaterial = SerializedMaterial & {
+    volume?: VolumeInput;
+};
+/** The medium an input describes, with its defaults, or `null` when it describes none. */
+export declare function volumeMedium(input: unknown): VolumeMedium | null;
+/**
+ * Rounded edges at shading time, from `userData.cpuRenderer.bevel`: the radius, in scene units, over which the
+ * shading normal turns from one face to the next, as Cycles' bevel node rounds them. A number, or `{ radius }`.
+ */
+export type BevelInput = number | {
+    radius?: number;
+};
+/** A snapshot material that may carry a bevel, as `serializeLive` attaches it. */
+export type BevelMaterial = SerializedMaterial & {
+    bevel?: BevelInput;
+};
+/** The bevel radius an input asks for, or zero. */
+export declare function bevelRadius(input: unknown): number;
+export declare function bevelOf(material: SerializedMaterial | undefined): number;
+export declare function volumeOf(material: SerializedMaterial | undefined): VolumeMedium | null;
+//# sourceMappingURL=subsurface.d.ts.map

package/dist/subsurface.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"subsurface.d.ts","sourceRoot":"","sources":["../src/subsurface.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,kBAAkB,EAAE,MAAM,eAAe,CAAC;AAExD;;;;;;GAMG;AACH,MAAM,MAAM,eAAe,GAAG,OAAO,GAAG;IACtC,+GAA+G;IAC/G,MAAM,CAAC,EAAE,MAAM,CAAC;IAChB;;;OAGG;IACH,MAAM,CAAC,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,CAAC;IACpC,+GAA+G;IAC/G,KAAK,CAAC,EAAE,MAAM,CAAC;IACf,sGAAsG;IACtG,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB,CAAC;AAEF,kDAAkD;AAClD,MAAM,WAAW,gBAAgB;IAC/B,MAAM,EAAE,MAAM,CAAC;IACf,kDAAkD;IAClD,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,UAAU,EAAE,MAAM,CAAC;CACpB;AAED,mHAAmH;AACnH,MAAM,MAAM,kBAAkB,GAAG,kBAAkB,GAAG;IAAE,UAAU,CAAC,EAAE,eAAe,CAAA;CAAE,CAAC;AAQvF,0FAA0F;AAC1F,wBAAgB,gBAAgB,CAAC,KAAK,EAAE,OAAO,GAAG,gBAAgB,GAAG,IAAI,CAyBxE;AAED,sDAAsD;AACtD,wBAAgB,YAAY,CAAC,QAAQ,EAAE,kBAAkB,GAAG,SAAS,GAAG,gBAAgB,GAAG,IAAI,CAE9F;AAED;;;GAGG;AACH,MAAM,MAAM,WAAW,GAAG;IACxB,qHAAqH;IACrH,OAAO,CAAC,EAAE,MAAM,CAAC;IACjB,0GAA0G;IAC1G,KAAK,CAAC,EAAE,MAAM,GAAG,SAAS,MAAM,EAAE,CAAC;IACnC,oHAAoH;IACpH,UAAU,CAAC,EAAE,MAAM,CAAC;CACrB,CAAC;AAEF,MAAM,WAAW,YAAY;IAC3B,OAAO,EAAE,MAAM,CAAC;IAChB,4CAA4C;IAC5C,MAAM,EAAE,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC;IACjC,UAAU,EAAE,MAAM,CAAC;CACpB;AAED,mFAAmF;AACnF,MAAM,MAAM,cAAc,GAAG,kBAAkB,GAAG;IAAE,MAAM,CAAC,EAAE,WAAW,CAAA;CAAE,CAAC;AAmB3E,0FAA0F;AAC1F,wBAAgB,YAAY,CAAC,KAAK,EAAE,OAAO,GAAG,YAAY,GAAG,IAAI,CAchE;AAED;;;GAGG;AACH,MAAM,MAAM,UAAU,GAAG,MAAM,GAAG;IAAE,MAAM,CAAC,EAAE,MAAM,CAAA;CAAE,CAAC;AAEtD,kFAAkF;AAClF,MAAM,MAAM,aAAa,GAAG,kBAAkB,GAAG;IAAE,KAAK,CAAC,EAAE,UAAU,CAAA;CAAE,CAAC;AAExE,mDAAmD;AACnD,wBAAgB,WAAW,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM,CAIlD;AAED,wBAAgB,OAAO,CAAC,QAAQ,EAAE,kBAAkB,GAAG,SAAS,GAAG,MAAM,CAExE;AAED,wBAAgB,QAAQ,CAAC,QAAQ,EAAE,kBAAkB,GAAG,SAAS,GAAG,YAAY,GAAG,IAAI,CAEtF"}

package/dist/subsurface.js

@@ -0,0 +1,70 @@
+const DEFAULT_RADIUS = [1, 0.2, 0.1];
+const DEFAULT_SCALE = 0.05;
+const finite = (value, fallback) => typeof value === "number" && Number.isFinite(value) ? value : fallback;
+/** The medium an input describes, with its defaults, or `null` when it describes none. */
+export function subsurfaceMedium(input) {
+    if (input === true)
+        input = {};
+    if (!input || typeof input !== "object")
+        return null;
+    const value = input;
+    const weight = Math.min(1, Math.max(0, finite(value.weight, 1)));
+    if (weight <= 0)
+        return null;
+    const scale = Math.max(0, finite(value.scale, DEFAULT_SCALE));
+    const radius = typeof value.radius === "number"
+        ? [value.radius, value.radius, value.radius]
+        : Array.isArray(value.radius) && value.radius.length >= 3
+            ? value.radius.slice(0, 3)
+            : DEFAULT_RADIUS;
+    // A path that never travels would divide by zero; a hair of distance keeps it a very dense medium instead.
+    const paths = radius.map((channel, k) => Math.max(1e-5, finite(channel, DEFAULT_RADIUS[k]) * scale));
+    return {
+        weight,
+        radius: [paths[0], paths[1], paths[2]],
+        anisotropy: Math.min(0.95, Math.max(-0.95, finite(value.anisotropy, 0))),
+    };
+}
+/** The medium a snapshot material carries, if any. */
+export function subsurfaceOf(material) {
+    return subsurfaceMedium(material?.subsurface);
+}
+const srgbToLinear = (c) => (c <= 0.04045 ? c / 12.92 : ((c + 0.055) / 1.055) ** 2.4);
+function albedoOf(color) {
+    if (Array.isArray(color) && color.length >= 3) {
+        return [0, 1, 2].map(k => Math.min(1, Math.max(0, finite(color[k], 1))));
+    }
+    const match = typeof color === "string" ? /^#?([0-9a-f]{6})$/i.exec(color) : null;
+    if (!match)
+        return [1, 1, 1];
+    const value = parseInt(match[1], 16);
+    return [(value >> 16) & 255, (value >> 8) & 255, value & 255].map(c => srgbToLinear(c / 255));
+}
+/** The medium an input describes, with its defaults, or `null` when it describes none. */
+export function volumeMedium(input) {
+    if (input === true)
+        input = {};
+    if (!input || typeof input !== "object")
+        return null;
+    const value = input;
+    const density = Math.max(0, finite(value.density, 1));
+    if (density <= 0)
+        return null;
+    return {
+        density,
+        albedo: albedoOf(value.color),
+        anisotropy: Math.min(0.95, Math.max(-0.95, finite(value.anisotropy, 0))),
+    };
+}
+/** The bevel radius an input asks for, or zero. */
+export function bevelRadius(input) {
+    const radius = typeof input === "number" ? input : input?.radius;
+    return typeof radius === "number" && Number.isFinite(radius) && radius > 0 ? radius : 0;
+}
+export function bevelOf(material) {
+    return bevelRadius(material?.bevel);
+}
+export function volumeOf(material) {
+    return volumeMedium(material?.volume);
+}
+//# sourceMappingURL=subsurface.js.map

package/dist/subsurface.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"subsurface.js","sourceRoot":"","sources":["../src/subsurface.ts"],"names":[],"mappings":"AAkCA,MAAM,cAAc,GAA6B,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;AAC/D,MAAM,aAAa,GAAG,IAAI,CAAC;AAE3B,MAAM,MAAM,GAAG,CAAC,KAAc,EAAE,QAAgB,EAAE,EAAE,CAClD,OAAO,KAAK,KAAK,QAAQ,IAAI,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,QAAQ,CAAC;AAEzE,0FAA0F;AAC1F,MAAM,UAAU,gBAAgB,CAAC,KAAc;IAC7C,IAAI,KAAK,KAAK,IAAI;QAAE,KAAK,GAAG,EAAE,CAAC;IAC/B,IAAI,CAAC,KAAK,IAAI,OAAO,KAAK,KAAK,QAAQ;QAAE,OAAO,IAAI,CAAC;IAErD,MAAM,KAAK,GAAG,KAA0C,CAAC;IACzD,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;IAEjE,IAAI,MAAM,IAAI,CAAC;QAAE,OAAO,IAAI,CAAC;IAE7B,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,KAAK,EAAE,aAAa,CAAC,CAAC,CAAC;IAE9D,MAAM,MAAM,GAAG,OAAO,KAAK,CAAC,MAAM,KAAK,QAAQ;QAC7C,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,EAAE,KAAK,CAAC,MAAM,CAAC;QAC5C,CAAC,CAAC,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,MAAM,CAAC,IAAI,KAAK,CAAC,MAAM,CAAC,MAAM,IAAI,CAAC;YACvD,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC,CAAC;YAC1B,CAAC,CAAC,cAAc,CAAC;IAErB,2GAA2G;IAC3G,MAAM,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,MAAM,CAAC,OAAO,EAAE,cAAc,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;IAErG,OAAO;QACL,MAAM;QACN,MAAM,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;QACtC,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,MAAM,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC;KACzE,CAAC;AACJ,CAAC;AAED,sDAAsD;AACtD,MAAM,UAAU,YAAY,CAAC,QAAwC;IACnE,OAAO,gBAAgB,CAAE,QAA2C,EAAE,UAAU,CAAC,CAAC;AACpF,CAAC;AAyBD,MAAM,YAAY,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC,CAAC;AAE9F,SAAS,QAAQ,CAAC,KAAc;IAC9B,IAAI,KAAK,CAAC,OAAO,CAAC,KAAK,CAAC,IAAI,KAAK,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QAC9C,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAA6B,CAAC;IACvG,CAAC;IAED,MAAM,KAAK,GAAG,OAAO,KAAK,KAAK,QAAQ,CAAC,CAAC,CAAC,oBAAoB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC;IAElF,IAAI,CAAC,KAAK;QAAE,OAAO,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAE7B,MAAM,KAAK,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;IAErC,OAAO,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,KAAK,IAAI,CAAC,CAAC,GAAG,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,YAAY,CAAC,CAAC,GAAG,GAAG,CAAC,CAClE,CAAC;AAC7B,CAAC;AAED,0FAA0F;AAC1F,MAAM,UAAU,YAAY,CAAC,KAAc;IACzC,IAAI,KAAK,KAAK,IAAI;QAAE,KAAK,GAAG,EAAE,CAAC;IAC/B,IAAI,CAAC,KAAK,IAAI,OAAO,KAAK,KAAK,QAAQ;QAAE,OAAO,IAAI,CAAC;IAErD,MAAM,KAAK,GAAG,KAAoB,CAAC;IACnC,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,CAAC,KAAK,CAAC,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC;IAEtD,IAAI,OAAO,IAAI,CAAC;QAAE,OAAO,IAAI,CAAC;IAE9B,OAAO;QACL,OAAO;QACP,MAAM,EAAE,QAAQ,CAAC,KAAK,CAAC,KAAK,CAAC;QAC7B,UAAU,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,EAAE,MAAM,CAAC,KAAK,CAAC,UAAU,EAAE,CAAC,CAAC,CAAC,CAAC;KACzE,CAAC;AACJ,CAAC;AAWD,mDAAmD;AACnD,MAAM,UAAU,WAAW,CAAC,KAAc;IACxC,MAAM,MAAM,GAAG,OAAO,KAAK,KAAK,QAAQ,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAE,KAAqC,EAAE,MAAM,CAAC;IAElG,OAAO,OAAO,MAAM,KAAK,QAAQ,IAAI,MAAM,CAAC,QAAQ,CAAC,MAAM,CAAC,IAAI,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;AAC1F,CAAC;AAED,MAAM,UAAU,OAAO,CAAC,QAAwC;IAC9D,OAAO,WAAW,CAAE,QAAsC,EAAE,KAAK,CAAC,CAAC;AACrE,CAAC;AAED,MAAM,UAAU,QAAQ,CAAC,QAAwC;IAC/D,OAAO,YAAY,CAAE,QAAuC,EAAE,MAAM,CAAC,CAAC;AACxE,CAAC"}

package/dist/texturePack.d.ts

@@ -9,10 +9,10 @@ import type { SerializedScene, SerializedTexture } from "./upstream.js";
  */
 export interface PackedTextures {
     /**
-     * One buffer for all of it: `TEXTURE_WORDS` words per texture first, as `TextureInfo` in `wgsl/common.wgsl.ts`, then
-     * the texels each entry points at. A byte texture takes one word per texel, RGBA in the low to high bytes; a float
-     * texture takes four, each the bits of one linear channel. Keeping the two in one binding is what holds the shader
-     * inside the eight storage buffers a WebGPU device has to offer.
+     * One buffer for all of it: the lookup tables first (see `TEXTURE_TABLE`), then `TEXTURE_WORDS` words per texture, as
+     * `TextureInfo` in `wgsl/common.wgsl.ts`, then the texels each entry points at. A byte texture takes one word per
+     * texel, RGBA in the low to high bytes; a float texture takes four, each the bits of one linear channel. Keeping them
+     * in one binding is what holds the shader inside the eight storage buffers a WebGPU device has to offer.
      */
     data: Uint32Array;
     count: number;
@@ -28,6 +28,17 @@ export interface TextureExtras {
 }
 /** Sixteen words per entry: see `TextureInfo` in `wgsl/common.wgsl.ts`. */
 export declare const TEXTURE_WORDS = 16;
+/**
+ * Lookup tables that lead the buffer, at offsets the shader knows as constants: the second Sobol dimension one byte of
+ * the index at a time (four tables of 256), then the linear value of each of the 256 sRGB codes. Both replace work the
+ * kernel did on every read with the very numbers it computed, so nothing they feed changes.
+ */
+export declare const SOBOL_TABLE = 0;
+export declare const SRGB_TABLE: number;
+/** Where the texture table starts, after the lookup tables. */
+export declare const TEXTURE_TABLE: number;
+/** Where a texture's entry starts in the buffer. */
+export declare const textureEntry: (index: number) => number;
 /** What a material's map index means when the material has no such map. */
 export declare const NO_TEXTURE = 4294967295;
 /** Three's own wrap constants, which the shader compares against. */
@@ -42,6 +53,7 @@ export declare const UV_CHANNELS = 2;
 export interface TexturePlan {
     textures: SerializedTexture[];
     indexOf: Map<string, number>;
+    /** The lookup tables and the texture table together: where the texels begin. */
     tableWords: number;
     /** Where the words appended after the texels begin. */
     extraOffset: number;

package/dist/texturePack.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"texturePack.d.ts","sourceRoot":"","sources":["../src/texturePack.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,eAAe,EAAE,iBAAiB,EAAE,MAAM,eAAe,CAAC;AAExE;;;;;;;GAOG;AACH,MAAM,WAAW,cAAc;IAC7B;;;;;OAKG;IACH,IAAI,EAAE,WAAW,CAAC;IAClB,KAAK,EAAE,MAAM,CAAC;IACd,wFAAwF;IACxF,OAAO,EAAE,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7B,wGAAwG;IACxG,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,oHAAoH;AACpH,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,EAAE,SAAS,iBAAiB,EAAE,CAAC;IACxC,KAAK,CAAC,EAAE,WAAW,CAAC;CACrB;AAED,2EAA2E;AAC3E,eAAO,MAAM,aAAa,KAAK,CAAC;AAEhC,2EAA2E;AAC3E,eAAO,MAAM,UAAU,aAAa,CAAC;AASrC,qEAAqE;AACrE,QAAA,MAAM,MAAM,OAAO,CAAC;AACpB,QAAA,MAAM,eAAe,OAAO,CAAC;AAE7B,+FAA+F;AAC/F,eAAO,MAAM,WAAW,IAAI,CAAC;AAY7B;;;GAGG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,EAAE,iBAAiB,EAAE,CAAC;IAC9B,OAAO,EAAE,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7B,UAAU,EAAE,MAAM,CAAC;IACnB,uDAAuD;IACvD,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,eAAe,EAAE,aAAa,GAAE,SAAS,iBAAiB,EAAO,GAAG,WAAW,CAalH;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,eAAe,EAAE,MAAM,GAAE,aAAkB,GAAG,cAAc,CAE/F;AAED,oFAAoF;AACpF,wBAAgB,aAAa,CAAC,IAAI,EAAE,WAAW,EAAE,UAAU,EAAE,WAAW,GAAG,cAAc,CAwDxF;AAED,OAAO,EAAE,eAAe,EAAE,MAAM,EAAE,CAAC"}
+{"version":3,"file":"texturePack.d.ts","sourceRoot":"","sources":["../src/texturePack.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,eAAe,EAAE,iBAAiB,EAAE,MAAM,eAAe,CAAC;AAExE;;;;;;;GAOG;AACH,MAAM,WAAW,cAAc;IAC7B;;;;;OAKG;IACH,IAAI,EAAE,WAAW,CAAC;IAClB,KAAK,EAAE,MAAM,CAAC;IACd,wFAAwF;IACxF,OAAO,EAAE,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7B,wGAAwG;IACxG,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,oHAAoH;AACpH,MAAM,WAAW,aAAa;IAC5B,QAAQ,CAAC,EAAE,SAAS,iBAAiB,EAAE,CAAC;IACxC,KAAK,CAAC,EAAE,WAAW,CAAC;CACrB;AAED,2EAA2E;AAC3E,eAAO,MAAM,aAAa,KAAK,CAAC;AAEhC;;;;GAIG;AACH,eAAO,MAAM,WAAW,IAAI,CAAC;AAC7B,eAAO,MAAM,UAAU,QAAwB,CAAC;AAChD,+DAA+D;AAC/D,eAAO,MAAM,aAAa,QAAmB,CAAC;AA6C9C,oDAAoD;AACpD,eAAO,MAAM,YAAY,GAAI,OAAO,MAAM,WAA0C,CAAC;AAErF,2EAA2E;AAC3E,eAAO,MAAM,UAAU,aAAa,CAAC;AASrC,qEAAqE;AACrE,QAAA,MAAM,MAAM,OAAO,CAAC;AACpB,QAAA,MAAM,eAAe,OAAO,CAAC;AAE7B,+FAA+F;AAC/F,eAAO,MAAM,WAAW,IAAI,CAAC;AAe7B;;;GAGG;AACH,MAAM,WAAW,WAAW;IAC1B,QAAQ,EAAE,iBAAiB,EAAE,CAAC;IAC9B,OAAO,EAAE,GAAG,CAAC,MAAM,EAAE,MAAM,CAAC,CAAC;IAC7B,gFAAgF;IAChF,UAAU,EAAE,MAAM,CAAC;IACnB,uDAAuD;IACvD,WAAW,EAAE,MAAM,CAAC;CACrB;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,eAAe,EAAE,aAAa,GAAE,SAAS,iBAAiB,EAAO,GAAG,WAAW,CAalH;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,eAAe,EAAE,MAAM,GAAE,aAAkB,GAAG,cAAc,CAE/F;AAED,oFAAoF;AACpF,wBAAgB,aAAa,CAAC,IAAI,EAAE,WAAW,EAAE,UAAU,EAAE,WAAW,GAAG,cAAc,CA+DxF;AAED,OAAO,EAAE,eAAe,EAAE,MAAM,EAAE,CAAC"}

package/dist/texturePack.js

@@ -1,5 +1,49 @@
 /** Sixteen words per entry: see `TextureInfo` in `wgsl/common.wgsl.ts`. */
 export const TEXTURE_WORDS = 16;
+/**
+ * Lookup tables that lead the buffer, at offsets the shader knows as constants: the second Sobol dimension one byte of
+ * the index at a time (four tables of 256), then the linear value of each of the 256 sRGB codes. Both replace work the
+ * kernel did on every read with the very numbers it computed, so nothing they feed changes.
+ */
+export const SOBOL_TABLE = 0;
+export const SRGB_TABLE = SOBOL_TABLE + 4 * 256;
+/** Where the texture table starts, after the lookup tables. */
+export const TEXTURE_TABLE = SRGB_TABLE + 256;
+/** The four byte tables of the Sobol sequence's second dimension, whose generator matrix is Pascal's triangle mod 2. */
+function sobolTables() {
+    const directions = [];
+    let direction = 0x80000000;
+    for (let bit = 0; bit < 32; bit++) {
+        directions.push(direction >>> 0);
+        direction = (direction ^ (direction >>> 1)) >>> 0;
+    }
+    const tables = new Uint32Array(4 * 256);
+    for (let byte = 0; byte < 4; byte++) {
+        for (let value = 0; value < 256; value++) {
+            let xor = 0;
+            for (let bit = 0; bit < 8; bit++)
+                if (value & (1 << bit))
+                    xor ^= directions[byte * 8 + bit];
+            tables[byte * 256 + value] = xor >>> 0;
+        }
+    }
+    return tables;
+}
+/** The linear value of every sRGB code, worked out in double precision and rounded once. */
+function srgbTable() {
+    return Float32Array.from({ length: 256 }, (_, code) => {
+        const value = Math.fround(code / 255);
+        return value <= 0.04045 ? value / 12.92 : ((value + 0.055) / 1.055) ** 2.4;
+    });
+}
+const LOOKUP = (() => {
+    const words = new Uint32Array(TEXTURE_TABLE);
+    words.set(sobolTables(), SOBOL_TABLE);
+    words.set(new Uint32Array(srgbTable().buffer), SRGB_TABLE);
+    return words;
+})();
+/** Where a texture's entry starts in the buffer. */
+export const textureEntry = (index) => TEXTURE_TABLE + index * TEXTURE_WORDS;
 /** What a material's map index means when the material has no such map. */
 export const NO_TEXTURE = 0xffffffff;
 const TEXTURE_FLOAT = 1;
@@ -13,6 +57,8 @@ const REPEAT = 1000;
 const MIRRORED_REPEAT = 1002;
 /** How many UV sets the geometry carries, and so the highest channel a texture may ask for. */
 export const UV_CHANNELS = 2;
+/** Whether a word's first byte is its lowest, as it is on the device: then RGBA bytes are already RGBA words. */
+const LITTLE_ENDIAN = new Uint8Array(new Uint32Array([1]).buffer)[0] === 1;
 const texelCount = (texture) => Math.max(0, texture.width) * Math.max(0, texture.height);
 const isFloat = (texture) => texture.data instanceof Float32Array;
 /**
@@ -27,7 +73,7 @@ export function planTextures(scene, extraTextures = []) {
         .filter(texture => texelCount(texture) > 0 && texture.data.length >= texelCount(texture) * 4)
         .filter(texture => !seen.has(texture.id) && !!seen.add(texture.id));
     const indexOf = new Map(textures.map((texture, index) => [texture.id, index]));
-    const tableWords = Math.max(1, textures.length) * TEXTURE_WORDS;
+    const tableWords = TEXTURE_TABLE + Math.max(1, textures.length) * TEXTURE_WORDS;
     const texelWords = textures.reduce((sum, texture) => sum + texelCount(texture) * (isFloat(texture) ? 4 : 1), 0);
     return { textures, indexOf, tableWords, extraOffset: tableWords + texelWords };
 }
@@ -42,8 +88,9 @@ export function writeTextures(plan, extraWords) {
     const words = data;
     const floats = new Float32Array(data.buffer);
     let texelAt = tableWords;
+    data.set(LOOKUP, 0);
     textures.forEach((texture, index) => {
-        const at = index * TEXTURE_WORDS;
+        const at = textureEntry(index);
         const texels = texelCount(texture);
         const float = isFloat(texture);
         words[at] = texture.width;
@@ -69,12 +116,18 @@ export function writeTextures(plan, extraWords) {
         }
         else {
             const source = texture.data;
-            for (let texel = 0; texel < texels; texel++) {
-                const offset = texel * 4;
-                words[texelAt + texel] = source[offset]
-                    | (source[offset + 1] << 8)
-                    | (source[offset + 2] << 16)
-                    | (source[offset + 3] << 24);
+            if (LITTLE_ENDIAN) {
+                // The bytes already sit in the order the words want them: copied whole, as the GPU reads them.
+                new Uint8Array(data.buffer, texelAt * 4, texels * 4).set(source.subarray(0, texels * 4));
+            }
+            else {
+                for (let texel = 0; texel < texels; texel++) {
+                    const offset = texel * 4;
+                    words[texelAt + texel] = source[offset]
+                        | (source[offset + 1] << 8)
+                        | (source[offset + 2] << 16)
+                        | (source[offset + 3] << 24);
+                }
             }
             texelAt += texels;
         }

package/dist/texturePack.js.map

@@ -1 +1 @@
-{"version":3,"file":"texturePack.js","sourceRoot":"","sources":["../src/texturePack.ts"],"names":[],"mappings":"AA+BA,2EAA2E;AAC3E,MAAM,CAAC,MAAM,aAAa,GAAG,EAAE,CAAC;AAEhC,2EAA2E;AAC3E,MAAM,CAAC,MAAM,UAAU,GAAG,UAAU,CAAC;AAErC,MAAM,aAAa,GAAG,CAAC,CAAC;AACxB,MAAM,YAAY,GAAG,CAAC,CAAC;AACvB,MAAM,eAAe,GAAG,CAAC,CAAC;AAC1B,MAAM,cAAc,GAAG,CAAC,CAAC;AACzB,mFAAmF;AACnF,MAAM,kBAAkB,GAAG,EAAE,CAAC;AAE9B,qEAAqE;AACrE,MAAM,MAAM,GAAG,IAAI,CAAC;AACpB,MAAM,eAAe,GAAG,IAAI,CAAC;AAE7B,+FAA+F;AAC/F,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAC;AAE7B,MAAM,UAAU,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;AAE5G,MAAM,OAAO,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,OAAO,CAAC,IAAI,YAAY,YAAY,CAAC;AAErF;;;GAGG;AACH,MAAM,SAAS,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,UAAU,KAAK,MAAM,CAAC;AAcrG,MAAM,UAAU,YAAY,CAAC,KAAsB,EAAE,gBAA8C,EAAE;IACnG,MAAM,IAAI,GAAG,IAAI,GAAG,EAAU,CAAC;IAE/B,wEAAwE;IACxE,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,QAAQ,IAAI,EAAE,EAAE,GAAG,aAAa,CAAC;SACzD,MAAM,CAAC,OAAO,CAAC,EAAE,CAAC,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,IAAI,CAAC,MAAM,IAAI,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;SAC5F,MAAM,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC;IAEtE,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,aAAa,CAAC;IAChE,MAAM,UAAU,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,OAAO,EAAE,EAAE,CAAC,GAAG,GAAG,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAEhH,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,UAAU,GAAG,UAAU,EAAE,CAAC;AACjF,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,KAAsB,EAAE,SAAwB,EAAE;IAC7E,OAAO,aAAa,CAAC,YAAY,CAAC,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,MAAM,CAAC,KAAK,IAAI,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;AACjG,CAAC;AAED,oFAAoF;AACpF,MAAM,UAAU,aAAa,CAAC,IAAiB,EAAE,UAAuB;IACtE,MAAM,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,IAAI,CAAC;IAE/C,oFAAoF;IACpF,MAAM,IAAI,GAAG,IAAI,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC;IAChF,MAAM,KAAK,GAAG,IAAI,CAAC;IACnB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC7C,IAAI,OAAO,GAAG,UAAU,CAAC;IAEzB,QAAQ,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,KAAK,EAAE,EAAE;QAClC,MAAM,EAAE,GAAG,KAAK,GAAG,aAAa,CAAC;QACjC,MAAM,MAAM,GAAG,UAAU,CAAC,OAAO,CAAC,CAAC;QACnC,MAAM,KAAK,GAAG,OAAO,CAAC,OAAO,CAAC,CAAC;QAE/B,KAAK,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC1B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QAC/B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;QACxB,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;cACpC,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC9B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC9B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,OAAO,GAAG,CAAC,CAAC,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC;QAC5E,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAElB,mGAAmG;QACnG,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,EAAE,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE3F,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;QAE1D,IAAI,KAAK,EAAE,CAAC;YACV,yEAAyE;YACzE,MAAM,CAAC,GAAG,CAAE,OAAO,CAAC,IAAqB,CAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YAC5E,OAAO,IAAI,MAAM,GAAG,CAAC,CAAC;QACxB,CAAC;aAAM,CAAC;YACN,MAAM,MAAM,GAAG,OAAO,CAAC,IAAkB,CAAC;YAE1C,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC;gBAC5C,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,CAAC;gBAEzB,KAAK,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;sBACnC,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;sBACzB,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;sBAC1B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;YACjC,CAAC;YAED,OAAO,IAAI,MAAM,CAAC;QACpB,CAAC;IACH,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IAE9B,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,CAAC;AACzE,CAAC;AAED,OAAO,EAAE,eAAe,EAAE,MAAM,EAAE,CAAC"}
+{"version":3,"file":"texturePack.js","sourceRoot":"","sources":["../src/texturePack.ts"],"names":[],"mappings":"AA+BA,2EAA2E;AAC3E,MAAM,CAAC,MAAM,aAAa,GAAG,EAAE,CAAC;AAEhC;;;;GAIG;AACH,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAC;AAC7B,MAAM,CAAC,MAAM,UAAU,GAAG,WAAW,GAAG,CAAC,GAAG,GAAG,CAAC;AAChD,+DAA+D;AAC/D,MAAM,CAAC,MAAM,aAAa,GAAG,UAAU,GAAG,GAAG,CAAC;AAE9C,wHAAwH;AACxH,SAAS,WAAW;IAClB,MAAM,UAAU,GAAa,EAAE,CAAC;IAChC,IAAI,SAAS,GAAG,UAAU,CAAC;IAE3B,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,EAAE,CAAC;QAClC,UAAU,CAAC,IAAI,CAAC,SAAS,KAAK,CAAC,CAAC,CAAC;QACjC,SAAS,GAAG,CAAC,SAAS,GAAG,CAAC,SAAS,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;IACpD,CAAC;IAED,MAAM,MAAM,GAAG,IAAI,WAAW,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAExC,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,IAAI,GAAG,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;QACpC,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,GAAG,EAAE,KAAK,EAAE,EAAE,CAAC;YACzC,IAAI,GAAG,GAAG,CAAC,CAAC;YAEZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE;gBAAE,IAAI,KAAK,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC;oBAAE,GAAG,IAAI,UAAU,CAAC,IAAI,GAAG,CAAC,GAAG,GAAG,CAAC,CAAC;YAE5F,MAAM,CAAC,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC;QACzC,CAAC;IACH,CAAC;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AAED,4FAA4F;AAC5F,SAAS,SAAS;IAChB,OAAO,YAAY,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,GAAG,EAAE,EAAE,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE;QACpD,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC;QAEtC,OAAO,KAAK,IAAI,OAAO,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,KAAK,CAAC,IAAI,GAAG,CAAC;IAC7E,CAAC,CAAC,CAAC;AACL,CAAC;AAED,MAAM,MAAM,GAAG,CAAC,GAAG,EAAE;IACnB,MAAM,KAAK,GAAG,IAAI,WAAW,CAAC,aAAa,CAAC,CAAC;IAE7C,KAAK,CAAC,GAAG,CAAC,WAAW,EAAE,EAAE,WAAW,CAAC,CAAC;IACtC,KAAK,CAAC,GAAG,CAAC,IAAI,WAAW,CAAC,SAAS,EAAE,CAAC,MAAM,CAAC,EAAE,UAAU,CAAC,CAAC;IAE3D,OAAO,KAAK,CAAC;AACf,CAAC,CAAC,EAAE,CAAC;AAEL,oDAAoD;AACpD,MAAM,CAAC,MAAM,YAAY,GAAG,CAAC,KAAa,EAAE,EAAE,CAAC,aAAa,GAAG,KAAK,GAAG,aAAa,CAAC;AAErF,2EAA2E;AAC3E,MAAM,CAAC,MAAM,UAAU,GAAG,UAAU,CAAC;AAErC,MAAM,aAAa,GAAG,CAAC,CAAC;AACxB,MAAM,YAAY,GAAG,CAAC,CAAC;AACvB,MAAM,eAAe,GAAG,CAAC,CAAC;AAC1B,MAAM,cAAc,GAAG,CAAC,CAAC;AACzB,mFAAmF;AACnF,MAAM,kBAAkB,GAAG,EAAE,CAAC;AAE9B,qEAAqE;AACrE,MAAM,MAAM,GAAG,IAAI,CAAC;AACpB,MAAM,eAAe,GAAG,IAAI,CAAC;AAE7B,+FAA+F;AAC/F,MAAM,CAAC,MAAM,WAAW,GAAG,CAAC,CAAC;AAE7B,iHAAiH;AACjH,MAAM,aAAa,GAAG,IAAI,UAAU,CAAC,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC;AAE3E,MAAM,UAAU,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,MAAM,CAAC,CAAC;AAE5G,MAAM,OAAO,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,OAAO,CAAC,IAAI,YAAY,YAAY,CAAC;AAErF;;;GAGG;AACH,MAAM,SAAS,GAAG,CAAC,OAA0B,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,OAAO,CAAC,IAAI,OAAO,CAAC,UAAU,KAAK,MAAM,CAAC;AAerG,MAAM,UAAU,YAAY,CAAC,KAAsB,EAAE,gBAA8C,EAAE;IACnG,MAAM,IAAI,GAAG,IAAI,GAAG,EAAU,CAAC;IAE/B,wEAAwE;IACxE,MAAM,QAAQ,GAAG,CAAC,GAAG,KAAK,CAAC,QAAQ,IAAI,EAAE,EAAE,GAAG,aAAa,CAAC;SACzD,MAAM,CAAC,OAAO,CAAC,EAAE,CAAC,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,IAAI,OAAO,CAAC,IAAI,CAAC,MAAM,IAAI,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,CAAC;SAC5F,MAAM,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,EAAE,CAAC,CAAC,CAAC;IAEtE,MAAM,OAAO,GAAG,IAAI,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC,OAAO,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC;IAC/E,MAAM,UAAU,GAAG,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,MAAM,CAAC,GAAG,aAAa,CAAC;IAChF,MAAM,UAAU,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,OAAO,EAAE,EAAE,CAAC,GAAG,GAAG,UAAU,CAAC,OAAO,CAAC,GAAG,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;IAEhH,OAAO,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,WAAW,EAAE,UAAU,GAAG,UAAU,EAAE,CAAC;AACjF,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,KAAsB,EAAE,SAAwB,EAAE;IAC7E,OAAO,aAAa,CAAC,YAAY,CAAC,KAAK,EAAE,MAAM,CAAC,QAAQ,CAAC,EAAE,MAAM,CAAC,KAAK,IAAI,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;AACjG,CAAC;AAED,oFAAoF;AACpF,MAAM,UAAU,aAAa,CAAC,IAAiB,EAAE,UAAuB;IACtE,MAAM,EAAE,QAAQ,EAAE,OAAO,EAAE,UAAU,EAAE,GAAG,IAAI,CAAC;IAE/C,oFAAoF;IACpF,MAAM,IAAI,GAAG,IAAI,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,WAAW,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC;IAChF,MAAM,KAAK,GAAG,IAAI,CAAC;IACnB,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC;IAC7C,IAAI,OAAO,GAAG,UAAU,CAAC;IAEzB,IAAI,CAAC,GAAG,CAAC,MAAM,EAAE,CAAC,CAAC,CAAC;IAEpB,QAAQ,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,KAAK,EAAE,EAAE;QAClC,MAAM,EAAE,GAAG,YAAY,CAAC,KAAK,CAAC,CAAC;QAC/B,MAAM,MAAM,GAAG,UAAU,CAAC,OAAO,CAAC,CAAC;QACnC,MAAM,KAAK,GAAG,OAAO,CAAC,OAAO,CAAC,CAAC;QAE/B,KAAK,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC1B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,MAAM,CAAC;QAC/B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC;QACxB,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,SAAS,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,OAAO,CAAC,OAAO,CAAC,CAAC,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC;cACvC,CAAC,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,cAAc,CAAC,CAAC,CAAC,CAAC,CAAC;cACpC,CAAC,OAAO,CAAC,UAAU,CAAC,CAAC,CAAC,kBAAkB,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC9B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,OAAO,CAAC,KAAK,CAAC;QAC9B,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,OAAO,GAAG,CAAC,CAAC,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC;QAC5E,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAElB,mGAAmG;QACnG,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,EAAE,MAAM,KAAK,CAAC,CAAC,CAAC,CAAC,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAE3F,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;QAE1D,IAAI,KAAK,EAAE,CAAC;YACV,yEAAyE;YACzE,MAAM,CAAC,GAAG,CAAE,OAAO,CAAC,IAAqB,CAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,EAAE,OAAO,CAAC,CAAC;YAC5E,OAAO,IAAI,MAAM,GAAG,CAAC,CAAC;QACxB,CAAC;aAAM,CAAC;YACN,MAAM,MAAM,GAAG,OAAO,CAAC,IAAkB,CAAC;YAE1C,IAAI,aAAa,EAAE,CAAC;gBAClB,+FAA+F;gBAC/F,IAAI,UAAU,CAAC,IAAI,CAAC,MAAM,EAAE,OAAO,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,EAAE,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;YAC3F,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,KAAK,GAAG,CAAC,EAAE,KAAK,GAAG,MAAM,EAAE,KAAK,EAAE,EAAE,CAAC;oBAC5C,MAAM,MAAM,GAAG,KAAK,GAAG,CAAC,CAAC;oBAEzB,KAAK,CAAC,OAAO,GAAG,KAAK,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;0BACnC,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;0BACzB,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC;0BAC1B,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,CAAC;gBACjC,CAAC;YACH,CAAC;YAED,OAAO,IAAI,MAAM,CAAC;QACpB,CAAC;IACH,CAAC,CAAC,CAAC;IAEH,IAAI,CAAC,GAAG,CAAC,UAAU,EAAE,OAAO,CAAC,CAAC;IAE9B,OAAO,EAAE,IAAI,EAAE,KAAK,EAAE,QAAQ,CAAC,MAAM,EAAE,OAAO,EAAE,WAAW,EAAE,OAAO,EAAE,CAAC;AACzE,CAAC;AAED,OAAO,EAAE,eAAe,EAAE,MAAM,EAAE,CAAC"}

package/dist/volumeContainment.d.ts

@@ -0,0 +1,17 @@
+import type { PackedGeometry } from "./geometry.js";
+/** What the kernel reads as no material at all. */
+export declare const NO_MATERIAL = 4294967295;
+/** The triangles that bound a volume, each with its material's index. */
+export interface VolumeTriangles {
+    /** Nine floats a triangle. */
+    positions: Float32Array;
+    materials: Uint32Array;
+}
+/** The boundaries of every volume in a packed scene, or `null` when it has none. */
+export declare function volumeTriangles(geometry: PackedGeometry): VolumeTriangles | null;
+/**
+ * The volume a point stands inside, by the parity of the boundary crossings of a ray from it: an odd count of one
+ * volume's triangles means inside it. The camera usually stands in none, so the kernel only asks this of the camera.
+ */
+export declare function containingVolume(volumes: VolumeTriangles | null, point: readonly number[]): number;
+//# sourceMappingURL=volumeContainment.d.ts.map

package/dist/volumeContainment.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"volumeContainment.d.ts","sourceRoot":"","sources":["../src/volumeContainment.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,cAAc,EAAE,MAAM,eAAe,CAAC;AAGpD,mDAAmD;AACnD,eAAO,MAAM,WAAW,aAAa,CAAC;AAEtC,yEAAyE;AACzE,MAAM,WAAW,eAAe;IAC9B,8BAA8B;IAC9B,SAAS,EAAE,YAAY,CAAC;IACxB,SAAS,EAAE,WAAW,CAAC;CACxB;AAED,oFAAoF;AACpF,wBAAgB,eAAe,CAAC,QAAQ,EAAE,cAAc,GAAG,eAAe,GAAG,IAAI,CAmBhF;AAKD;;;GAGG;AACH,wBAAgB,gBAAgB,CAAC,OAAO,EAAE,eAAe,GAAG,IAAI,EAAE,KAAK,EAAE,SAAS,MAAM,EAAE,GAAG,MAAM,CAsClG"}

package/dist/volumeContainment.js

@@ -0,0 +1,61 @@
+import { MATERIAL_VOLUME } from "./materialLayout.js";
+/** What the kernel reads as no material at all. */
+export const NO_MATERIAL = 0xffffffff;
+/** The boundaries of every volume in a packed scene, or `null` when it has none. */
+export function volumeTriangles(geometry) {
+    const volume = geometry.materialProperties.map(properties => (properties.flags & MATERIAL_VOLUME) !== 0);
+    if (!volume.some(Boolean))
+        return null;
+    const picked = [];
+    for (let triangle = 0; triangle < geometry.triangleCount; triangle++)
+        if (volume[geometry.materials[triangle]])
+            picked.push(triangle);
+    const positions = new Float32Array(picked.length * 9);
+    const materials = new Uint32Array(picked.length);
+    picked.forEach((triangle, k) => {
+        positions.set(geometry.positions.subarray(triangle * 9, triangle * 9 + 9), k * 9);
+        materials[k] = geometry.materials[triangle];
+    });
+    return { positions, materials };
+}
+/** An irrational direction, so a ray from the camera does not run along an edge of an axis aligned box. */
+const PROBE = [0.4123, 0.7071, 0.5745].map((value, _, all) => value / Math.hypot(...all));
+/**
+ * The volume a point stands inside, by the parity of the boundary crossings of a ray from it: an odd count of one
+ * volume's triangles means inside it. The camera usually stands in none, so the kernel only asks this of the camera.
+ */
+export function containingVolume(volumes, point) {
+    if (!volumes)
+        return NO_MATERIAL;
+    const crossings = new Map();
+    const p = volumes.positions;
+    const [dx, dy, dz] = PROBE;
+    for (let k = 0; k < volumes.materials.length; k++) {
+        const at = k * 9;
+        const e1x = p[at + 3] - p[at], e1y = p[at + 4] - p[at + 1], e1z = p[at + 5] - p[at + 2];
+        const e2x = p[at + 6] - p[at], e2y = p[at + 7] - p[at + 1], e2z = p[at + 8] - p[at + 2];
+        // Moller and Trumbore.
+        const hx = dy * e2z - dz * e2y, hy = dz * e2x - dx * e2z, hz = dx * e2y - dy * e2x;
+        const det = e1x * hx + e1y * hy + e1z * hz;
+        if (Math.abs(det) < 1e-12)
+            continue;
+        const inverse = 1 / det;
+        const sx = point[0] - p[at], sy = point[1] - p[at + 1], sz = point[2] - p[at + 2];
+        const u = (sx * hx + sy * hy + sz * hz) * inverse;
+        if (u < 0 || u > 1)
+            continue;
+        const qx = sy * e1z - sz * e1y, qy = sz * e1x - sx * e1z, qz = sx * e1y - sy * e1x;
+        const v = (dx * qx + dy * qy + dz * qz) * inverse;
+        if (v < 0 || u + v > 1)
+            continue;
+        if ((e2x * qx + e2y * qy + e2z * qz) * inverse > 0) {
+            const material = volumes.materials[k];
+            crossings.set(material, (crossings.get(material) ?? 0) + 1);
+        }
+    }
+    for (const [material, count] of crossings)
+        if (count % 2 === 1)
+            return material;
+    return NO_MATERIAL;
+}
+//# sourceMappingURL=volumeContainment.js.map

package/dist/volumeContainment.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"volumeContainment.js","sourceRoot":"","sources":["../src/volumeContainment.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,eAAe,EAAE,MAAM,qBAAqB,CAAC;AAEtD,mDAAmD;AACnD,MAAM,CAAC,MAAM,WAAW,GAAG,UAAU,CAAC;AAStC,oFAAoF;AACpF,MAAM,UAAU,eAAe,CAAC,QAAwB;IACtD,MAAM,MAAM,GAAG,QAAQ,CAAC,kBAAkB,CAAC,GAAG,CAAC,UAAU,CAAC,EAAE,CAAC,CAAC,UAAU,CAAC,KAAK,GAAG,eAAe,CAAC,KAAK,CAAC,CAAC,CAAC;IAEzG,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,OAAO,CAAC;QAAE,OAAO,IAAI,CAAC;IAEvC,MAAM,MAAM,GAAa,EAAE,CAAC;IAE5B,KAAK,IAAI,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAAC,aAAa,EAAE,QAAQ,EAAE;QAClE,IAAI,MAAM,CAAC,QAAQ,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC;YAAE,MAAM,CAAC,IAAI,CAAC,QAAQ,CAAC,CAAC;IAElE,MAAM,SAAS,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACtD,MAAM,SAAS,GAAG,IAAI,WAAW,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC;IAEjD,MAAM,CAAC,OAAO,CAAC,CAAC,QAAQ,EAAE,CAAC,EAAE,EAAE;QAC7B,SAAS,CAAC,GAAG,CAAC,QAAQ,CAAC,SAAS,CAAC,QAAQ,CAAC,QAAQ,GAAG,CAAC,EAAE,QAAQ,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAClF,SAAS,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,SAAS,CAAC,QAAQ,CAAC,CAAC;IAC9C,CAAC,CAAC,CAAC;IAEH,OAAO,EAAE,SAAS,EAAE,SAAS,EAAE,CAAC;AAClC,CAAC;AAED,2GAA2G;AAC3G,MAAM,KAAK,GAAG,CAAC,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;AAE1F;;;GAGG;AACH,MAAM,UAAU,gBAAgB,CAAC,OAA+B,EAAE,KAAwB;IACxF,IAAI,CAAC,OAAO;QAAE,OAAO,WAAW,CAAC;IAEjC,MAAM,SAAS,GAAG,IAAI,GAAG,EAAkB,CAAC;IAC5C,MAAM,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC;IAC5B,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,KAAK,CAAC;IAE3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAClD,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACxF,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QACxF,uBAAuB;QACvB,MAAM,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;QACnF,MAAM,GAAG,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC;QAE3C,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,KAAK;YAAE,SAAS;QAEpC,MAAM,OAAO,GAAG,CAAC,GAAG,GAAG,CAAC;QACxB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAClF,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC;QAElD,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC;YAAE,SAAS;QAE7B,MAAM,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,CAAC;QACnF,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC;QAElD,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC;YAAE,SAAS;QAEjC,IAAI,CAAC,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE,CAAC,GAAG,OAAO,GAAG,CAAC,EAAE,CAAC;YACnD,MAAM,QAAQ,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC;YAEtC,SAAS,CAAC,GAAG,CAAC,QAAQ,EAAE,CAAC,SAAS,CAAC,GAAG,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9D,CAAC;IACH,CAAC;IAED,KAAK,MAAM,CAAC,QAAQ,EAAE,KAAK,CAAC,IAAI,SAAS;QAAE,IAAI,KAAK,GAAG,CAAC,KAAK,CAAC;YAAE,OAAO,QAAQ,CAAC;IAEhF,OAAO,WAAW,CAAC;AACrB,CAAC"}

package/dist/wgsl/bsdf.wgsl.d.ts

@@ -9,5 +9,5 @@
  * for the diffuse lobe, the visible normal distribution (Heitz 2018) for GGX, anisotropic or not, and the Blinn-Phong
  * lobe for Phong, mixed by how much each lobe reflects.
  */
-export declare const BSDF_WGSL = "\n/** How a surface answers a direction: what it reflects, and how likely sampling was to have chosen it. */\nstruct BsdfEval {\n  // The bidirectional distribution times the cosine at the surface, which is what a contribution is multiplied by.\n  value: vec3<f32>,\n  // Solid angle density of this direction under this renderer's own sampling, for the power heuristic.\n  pdf: f32,\n};\n\nstruct BsdfSample {\n  direction: vec3<f32>,\n  // The distribution over its own density: what the path's throughput is multiplied by.\n  weight: vec3<f32>,\n  pdf: f32,\n  valid: bool,\n  // The diffuse lobe chose the direction, rather than a specular one: which of the denoiser's two signals it feeds.\n  diffuse: bool,\n};\n\nfn luminanceOf(color: vec3<f32>) -> f32 {\n  return dot(color, vec3<f32>(0.2126, 0.7152, 0.0722));\n}\n\nfn max3(color: vec3<f32>) -> f32 {\n  return max(color.r, max(color.g, color.b));\n}\n\n/** Schlick's approximation, as in Three: the reflectance always reaches one at a grazing angle. */\nfn fresnelSchlick(f0: vec3<f32>, cosine: f32) -> vec3<f32> {\n  return schlick3(f0, 1.0, cosine);\n}\n\n/** Reflectance at normal incidence: the dielectric's own, turned towards the base colour as the surface turns metal. */\nfn surfaceF0(surface: Surface) -> vec3<f32> {\n  return mix(surface.specularColor, surface.albedo, surface.metalness);\n}\n\n/** GGX, with alpha the square of the roughness, as Three parameterises it. */\nfn ggxDistribution(alpha2: f32, nh: f32) -> f32 {\n  let denominator = nh * nh * (alpha2 - 1.0) + 1.0;\n\n  return alpha2 / (PI * denominator * denominator);\n}\n\n/** Smith's height correlated visibility, which already carries the 1/(4 nl nv) the lobe would otherwise need. */\nfn ggxVisibility(alpha2: f32, nl: f32, nv: f32) -> f32 {\n  return 0.5 / (nl * sqrt(nv * nv * (1.0 - alpha2) + alpha2) + nv * sqrt(nl * nl * (1.0 - alpha2) + alpha2) + 1e-7);\n}\n\n/** The separable masking term, which is the one the visible normal sampling is built on. */\nfn ggxMasking(alpha2: f32, nv: f32) -> f32 {\n  return 2.0 * nv / (nv + sqrt(alpha2 + (1.0 - alpha2) * nv * nv));\n}\n\n/** Karis' analytic fit of the GGX directional albedo (Three's DFGApprox): the scale and bias applied to F0 and F90. */\nfn dfgApprox(roughness: f32, nv: f32) -> vec2<f32> {\n  let r = roughness * vec4<f32>(-1.0, -0.0275, -0.572, 0.022) + vec4<f32>(1.0, 0.0425, 1.04, -0.04);\n  let a004 = min(r.x * r.x, exp2(-9.28 * nv)) * r.x + r.y;\n\n  return vec2<f32>(-1.04, 1.04) * a004 + r.zw;\n}\n\n/** Three's IBLSheenBRDF: the Charlie sheen lobe integrated over the hemisphere. */\nfn sheenAlbedo(cosine: f32, roughness: f32) -> f32 {\n  let r2 = roughness * roughness;\n  let inverse = 1.0 / (roughness + 0.1);\n  let a = -1.9362 + 1.0678 * roughness + 0.4573 * r2 - 0.8469 * inverse;\n  let b = -0.6014 + 0.5538 * roughness - 0.4670 * r2 - 0.1255 * inverse;\n\n  return clamp(exp(a * clamp(cosine, 0.0, 1.0) + b), 0.0, 1.0);\n}\n\n/** Three's D_Charlie times V_Neubelt. */\nfn sheenLobe(roughness: f32, nh: f32, nv: f32, nl: f32) -> f32 {\n  let inverse = 1.0 / (roughness * roughness);\n  let sin2h = max(1.0 - nh * nh, 0.0078125);\n  let distribution = (2.0 + inverse) * pow(sin2h, inverse * 0.5) / (2.0 * PI);\n\n  return distribution * clamp(1.0 / (4.0 * (nl + nv - nl * nv)), 0.0, 1.0);\n}\n\n/** The frame the specular lobe is defined in: Three's anisotropy direction when there is one, any frame otherwise. */\nfn specularBasis(surface: Surface) -> mat3x3<f32> {\n  let normal = surface.normal;\n\n  if (surface.anisotropy > 0.0) {\n    let t = surface.anisotropyT - normal * dot(normal, surface.anisotropyT);\n    let tangent = normalize(t);\n\n    if (all(tangent == tangent) && dot(t, t) > 1e-12) {\n      return mat3x3<f32>(tangent, cross(normal, tangent), normal);\n    }\n  }\n\n  return orthonormalBasis(normal);\n}\n\n/** The GGX roughness along the frame's two axes: Three's alphaT along the anisotropy, alpha across it. */\nfn specularAlphas(surface: Surface) -> vec2<f32> {\n  let alpha = surface.roughness * surface.roughness;\n\n  if (surface.anisotropy > 0.0) {\n    return vec2<f32>(max(surface.alphaT, 1e-4), max(alpha, 1e-4));\n  }\n\n  return vec2<f32>(max(alpha, 1e-4));\n}\n\n/** Anisotropic GGX in the lobe's own frame; with equal alphas it is exactly the isotropic distribution. */\nfn ggxAnisotropic(h: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let x = h.x / alphas.x;\n  let y = h.y / alphas.y;\n  let d = x * x + y * y + h.z * h.z;\n\n  return 1.0 / (PI * alphas.x * alphas.y * d * d);\n}\n\n/** Smith's masking of one direction, in the lobe's frame. */\nfn ggxMaskingAnisotropic(v: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let projected = alphas.x * alphas.x * v.x * v.x + alphas.y * alphas.y * v.y * v.y;\n\n  return 2.0 * v.z / (v.z + sqrt(v.z * v.z + projected));\n}\n\n/** Three's V_GGX_SmithCorrelated_Anisotropic, in the lobe's frame. */\nfn ggxVisibilityAnisotropic(l: vec3<f32>, v: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let gv = l.z * length(vec3<f32>(alphas.x * v.x, alphas.y * v.y, v.z));\n  let gl = v.z * length(vec3<f32>(alphas.x * l.x, alphas.y * l.y, l.z));\n\n  return 0.5 / max(gv + gl, 1e-7);\n}\n\n/**\n * A half vector drawn from the distribution of normals the view direction can actually see (Heitz 2018), which is what\n * keeps the specular lobe from wasting samples on microfacets facing away. The view is in the lobe's frame.\n */\nfn sampleVisibleNormal(view: vec3<f32>, alphas: vec2<f32>, u: vec2<f32>) -> vec3<f32> {\n  let stretched = normalize(vec3<f32>(alphas.x * view.x, alphas.y * view.y, view.z));\n  let lengthSquared = stretched.x * stretched.x + stretched.y * stretched.y;\n  let t1 = select(\n    vec3<f32>(1.0, 0.0, 0.0),\n    vec3<f32>(-stretched.y, stretched.x, 0.0) * inverseSqrt(max(lengthSquared, 1e-20)),\n    lengthSquared > 1e-20);\n  let t2 = cross(stretched, t1);\n  let radius = sqrt(u.x);\n  let angle = 2.0 * PI * u.y;\n  let px = radius * cos(angle);\n  let raw = radius * sin(angle);\n  let lift = 0.5 * (1.0 + stretched.z);\n  let py = (1.0 - lift) * sqrt(max(0.0, 1.0 - px * px)) + lift * raw;\n  let normal = t1 * px + t2 * py + stretched * sqrt(max(0.0, 1.0 - px * px - py * py));\n\n  return normalize(vec3<f32>(alphas.x * normal.x, alphas.y * normal.y, max(normal.z, 1e-6)));\n}\n\n/** What the clearcoat takes back from everything under it, seen from this view: Three's (1 - clearcoat * Fcc). */\nfn coatTransmittance(surface: Surface, view: vec3<f32>) -> f32 {\n  if (surface.clearcoat <= 0.0) {\n    return 1.0;\n  }\n\n  let nv = clamp(dot(surface.clearcoatNormal, view), 0.0, 1.0);\n\n  return 1.0 - surface.clearcoat * schlick3(vec3<f32>(0.04), 1.0, nv).x;\n}\n\n/** How often each lobe is chosen, as (diffuse, specular, clearcoat): by how much of the reflection each accounts for. */\nfn lobeWeights(surface: Surface, view: vec3<f32>) -> vec3<f32> {\n  let nv = max(dot(surface.normal, view), 1e-4);\n\n  if (surface.kind == KIND_PHONG) {\n    let diffuse = luminanceOf(surface.albedo) + 0.02;\n    let specular = luminanceOf(surface.specularColor) * surface.specularStrength * (surface.shininess * 0.125 + 1.0) + 0.02;\n\n    return vec3<f32>(diffuse, specular, 0.0) / (diffuse + specular);\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return vec3<f32>(1.0, 0.0, 0.0);\n  }\n\n  let base = coatTransmittance(surface, view);\n  // The floors keep every lobe reachable: the diffuse one also carries the sheen and the multiple scattering term.\n  let diffuse = (luminanceOf(surface.albedo * (1.0 - surface.metalness)) + max3(surface.sheenColor)) * base + 0.02;\n  let specular = luminanceOf(schlick3(surface.specularColor, surface.specularF90, nv)\n    * (1.0 - surface.metalness) + surface.albedo * surface.metalness) * base + 0.05;\n  let coat = select(0.0, (1.0 - base) + 0.02, surface.clearcoat > 0.0);\n  let total = diffuse + specular + coat;\n\n  return vec3<f32>(diffuse, specular, coat) / total;\n}\n\n/** Three's physical direct lighting for one light direction, times the cosine. */\nfn physicalValue(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> vec3<f32> {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let nv = max(dot(normal, view), 1e-4);\n  var value = vec3<f32>(0.0);\n\n  if (nl > 0.0) {\n    let basis = specularBasis(surface);\n    let alphas = specularAlphas(surface);\n    let half = normalize(light + view);\n    let vh = clamp(dot(view, half), 0.0, 1.0);\n    let f0 = surfaceF0(surface);\n    var fresnel = schlick3(f0, surface.specularF90, vh);\n\n    if (surface.iridescence > 0.0) {\n      fresnel = mix(fresnel, surface.iridescenceFresnel, surface.iridescence);\n    }\n\n    let localLight = vec3<f32>(dot(light, basis[0]), dot(light, basis[1]), nl);\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), nv);\n    let localHalf = vec3<f32>(dot(half, basis[0]), dot(half, basis[1]), max(dot(half, normal), 0.0));\n    let single = fresnel * (ggxVisibilityAnisotropic(localLight, localView, alphas) * ggxAnisotropic(localHalf, alphas));\n\n    // Three's BRDF_GGX_Multiscatter: the energy single scattering loses at high roughness, returned.\n    let dfgV = dfgApprox(surface.roughness, nv);\n    let dfgL = dfgApprox(surface.roughness, nl);\n    let essV = f0 * dfgV.x + surface.specularF90 * dfgV.y;\n    let essL = f0 * dfgL.x + surface.specularF90 * dfgL.y;\n    let emsV = 1.0 - (dfgV.x + dfgV.y);\n    let emsL = 1.0 - (dfgL.x + dfgL.y);\n    let average = f0 + (vec3<f32>(1.0) - f0) * 0.047619;\n    let multiple = essV * essL * average / (vec3<f32>(1.0) - emsV * emsL * average + 1e-7) * (emsV * emsL);\n\n    var diffuse = surface.albedo * (1.0 - surface.metalness) / PI;\n    var specular = single + multiple;\n\n    // Three's sheen: its lobe sits on top and takes back the energy it reflects from the layers below.\n    if (any(surface.sheenColor > vec3<f32>(0.0))) {\n      let energy = 1.0 - max3(surface.sheenColor)\n        * max(sheenAlbedo(nv, surface.sheenRoughness), sheenAlbedo(nl, surface.sheenRoughness));\n\n      diffuse = diffuse * energy;\n      specular = specular * energy\n        + surface.sheenColor * sheenLobe(surface.sheenRoughness, clamp(dot(normal, half), 0.0, 1.0), nv, nl);\n    }\n\n    value = (diffuse + specular) * nl;\n  }\n\n  if (surface.clearcoat > 0.0) {\n    let coatNormal = surface.clearcoatNormal;\n    value = value * coatTransmittance(surface, view);\n\n    let nlCoat = dot(coatNormal, light);\n    let nvCoat = max(dot(coatNormal, view), 1e-4);\n\n    if (nlCoat > 0.0) {\n      let half = normalize(light + view);\n      let alpha = surface.clearcoatRoughness * surface.clearcoatRoughness;\n      let alpha2 = alpha * alpha;\n      let fresnel = schlick3(vec3<f32>(0.04), 1.0, clamp(dot(view, half), 0.0, 1.0)).x;\n      let lobe = fresnel * ggxVisibility(alpha2, nlCoat, nvCoat) * ggxDistribution(alpha2, max(dot(coatNormal, half), 0.0));\n\n      value = value + vec3<f32>(surface.clearcoat * lobe * nlCoat);\n    }\n  }\n\n  return value;\n}\n\n/** Three's D_BlinnPhong times its implicit visibility of a quarter: the Phong highlight's distribution. */\nfn blinnPhong(shininess: f32, nh: f32) -> f32 {\n  return 0.25 * (shininess * 0.5 + 1.0) / PI * pow(max(nh, 0.0), shininess);\n}\n\n/** Density of the Blinn-Phong lobe's sampling, per solid angle of the light direction. */\nfn blinnPhongPdf(shininess: f32, nh: f32, vh: f32) -> f32 {\n  return (shininess + 1.0) / (2.0 * PI) * pow(max(nh, 0.0), shininess) / max(4.0 * vh, 1e-6);\n}\n\n/** Density with which the GGX lobe (VNDF) would have produced this light direction, in the lobe's frame. */\nfn ggxPdf(localLight: vec3<f32>, localView: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let half = normalize(localLight + localView);\n\n  if (half.z <= 0.0 || localView.z <= 0.0) {\n    return 0.0;\n  }\n\n  return ggxMaskingAnisotropic(localView, alphas) * ggxAnisotropic(half, alphas) / (4.0 * localView.z);\n}\n\n/** The mixture density of every lobe this surface samples with. */\nfn bsdfPdf(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> f32 {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let weights = lobeWeights(surface, view);\n  var pdf = weights.x * max(nl, 0.0) / PI;\n\n  if (surface.kind == KIND_PHONG) {\n    let half = normalize(light + view);\n\n    return pdf + weights.y * blinnPhongPdf(surface.shininess, dot(normal, half), max(dot(view, half), 0.0));\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return pdf;\n  }\n\n  let basis = specularBasis(surface);\n  let alphas = specularAlphas(surface);\n\n  pdf = pdf + weights.y * ggxPdf(\n    vec3<f32>(dot(light, basis[0]), dot(light, basis[1]), nl),\n    vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, normal)),\n    alphas);\n\n  if (weights.z > 0.0) {\n    let coatBasis = orthonormalBasis(surface.clearcoatNormal);\n    let alpha = max(surface.clearcoatRoughness * surface.clearcoatRoughness, 1e-4);\n\n    pdf = pdf + weights.z * ggxPdf(\n      vec3<f32>(dot(light, coatBasis[0]), dot(light, coatBasis[1]), dot(light, surface.clearcoatNormal)),\n      vec3<f32>(dot(view, coatBasis[0]), dot(view, coatBasis[1]), dot(view, surface.clearcoatNormal)),\n      vec2<f32>(alpha));\n  }\n\n  return pdf;\n}\n\n/**\n * What the surface sends towards the view direction, and the density of that direction. A Lambert material keeps\n * Three's Lambert term alone, which is all its own shading has; a toon material reads its ramp in place of the cosine.\n */\nfn evaluateBsdf(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> BsdfEval {\n  let normal = surface.normal;\n  var result: BsdfEval;\n  result.value = vec3<f32>(0.0);\n  result.pdf = 0.0;\n\n  let nl = dot(normal, light);\n  let nv = dot(normal, view);\n\n  if (nv <= 0.0) {\n    return result;\n  }\n\n  if (surface.kind == KIND_TOON) {\n    // Three's toon lighting reads its ramp across the whole sphere of directions, the side facing away included.\n    result.value = surface.albedo / PI * toonRamp(surface.gradientMap, nl);\n    result.pdf = max(nl, 0.0) / PI;\n\n    return result;\n  }\n\n  // A clearcoat may face a light its base turns away from, so only a surface without one stops here.\n  if (nl <= 0.0 && surface.clearcoat <= 0.0) {\n    return result;\n  }\n\n  if ((surface.flags & MATERIAL_PBR) != 0u) {\n    result.value = physicalValue(surface, view, light);\n  } else if (surface.kind == KIND_PHONG) {\n    let half = normalize(light + view);\n    let fresnel = schlick3(surface.specularColor, 1.0, clamp(dot(view, half), 0.0, 1.0));\n    let specular = fresnel * blinnPhong(surface.shininess, dot(normal, half)) * surface.specularStrength;\n\n    result.value = (surface.albedo / PI + specular) * max(nl, 0.0);\n  } else {\n    result.value = surface.albedo / PI * max(nl, 0.0);\n  }\n\n  result.pdf = bsdfPdf(surface, view, light);\n\n  return result;\n}\n\n/**\n * Reflectance of a smooth dielectric interface, from both polarisations rather than Schlick's fit: refraction needs the\n * exact cut-off past which nothing gets through. eta is the index the ray leaves over the index it enters.\n */\nfn fresnelDielectric(cosine: f32, eta: f32) -> f32 {\n  let c = clamp(cosine, 0.0, 1.0);\n  let sin2 = eta * eta * (1.0 - c * c);\n\n  if (sin2 >= 1.0) {\n    return 1.0;\n  }\n\n  let t = sqrt(1.0 - sin2);\n  let s = (eta * c - t) / (eta * c + t);\n  let p = (c - eta * t) / (c + eta * t);\n\n  return 0.5 * (s * s + p * p);\n}\n\n/** Where a transmissive surface sends the path: off it, or through it. */\nstruct DielectricSample {\n  direction: vec3<f32>,\n  weight: vec3<f32>,\n  // The path went through the surface instead of reflecting off it.\n  crossed: bool,\n  // It went through without bending, so the density that found the next hit is still the one before this surface.\n  straight: bool,\n  valid: bool,\n};\n\n/**\n * The transmissive share of a physical material. Reflection and refraction are chosen in proportion to the Fresnel\n * reflectance on a microfacet drawn from the visible normals, which leaves only the masking of the outgoing direction\n * in the weight. Light is tinted by the base colour once, on the way in, which is where Three applies it.\n *\n * A thin wall (zero thickness, Three's default) bends nothing: what is not reflected carries straight on, exactly as a\n * shadow ray crossing it does in thinCrossing, so the two estimators agree on a window. ior is the index this path\n * sees, which dispersion makes depend on the colour the path carries.\n */\nfn sampleDielectric(surface: Surface, view: vec3<f32>, ior: f32, lobe: f32, u: vec2<f32>) -> DielectricSample {\n  let normal = surface.normal;\n  var result: DielectricSample;\n  result.valid = false;\n  result.crossed = false;\n  result.straight = false;\n  result.weight = vec3<f32>(0.0);\n\n  let nv = dot(normal, view);\n\n  if (nv <= 0.0) {\n    return result;\n  }\n\n  let alpha = max(surface.roughness * surface.roughness, 1e-4);\n  let alpha2 = alpha * alpha;\n  let basis = orthonormalBasis(normal);\n  let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), nv);\n  let localHalf = sampleVisibleNormal(localView, vec2<f32>(alpha), u);\n  let half = basis[0] * localHalf.x + basis[1] * localHalf.y + normal * localHalf.z;\n\n  if (surface.thickness <= 0.0) {\n    let reflectance = fresnelDielectric(abs(dot(surface.facing, view)), 1.0 / ior);\n\n    if (lobe >= reflectance) {\n      result.direction = -view;\n      result.weight = surface.albedo;\n      result.crossed = true;\n      result.straight = true;\n      result.valid = true;\n\n      return result;\n    }\n\n    result.direction = reflect(-view, half);\n  } else {\n    let eta = select(ior, 1.0 / ior, surface.front);\n    let refracted = refract(-view, half, eta);\n    // Past the critical angle refract() returns nothing, and everything reflects.\n    let reflectance = select(fresnelDielectric(dot(view, half), eta), 1.0, all(refracted == vec3<f32>(0.0)));\n\n    if (lobe >= reflectance) {\n      result.direction = normalize(refracted);\n\n      let no = -dot(normal, result.direction);\n\n      if (no <= 0.0) {\n        return result;\n      }\n\n      result.weight = vec3<f32>(ggxMasking(alpha2, no)) * select(vec3<f32>(1.0), surface.albedo, surface.front);\n      result.crossed = true;\n      result.valid = true;\n\n      return result;\n    }\n\n    result.direction = reflect(-view, half);\n  }\n\n  let nl = dot(normal, result.direction);\n\n  if (nl <= 0.0) {\n    return result;\n  }\n\n  result.weight = vec3<f32>(ggxMasking(alpha2, nl));\n  result.valid = true;\n\n  return result;\n}\n\n/** A half vector around the normal distributed as cos^s, which is the Blinn-Phong lobe's own shape. */\nfn sampleBlinnPhongHalf(normal: vec3<f32>, shininess: f32, u: vec2<f32>) -> vec3<f32> {\n  let cosine = pow(u.x, 1.0 / (shininess + 1.0));\n  let sine = sqrt(max(0.0, 1.0 - cosine * cosine));\n  let angle = 2.0 * PI * u.y;\n\n  return normalize(orthonormalBasis(normal) * vec3<f32>(sine * cos(angle), sine * sin(angle), cosine));\n}\n\n/** Which way the path carries on, and what its throughput is multiplied by. */\nfn sampleBsdf(surface: Surface, view: vec3<f32>, lobe: f32, u: vec2<f32>) -> BsdfSample {\n  let normal = surface.normal;\n  var result: BsdfSample;\n  result.valid = false;\n  result.weight = vec3<f32>(0.0);\n  result.pdf = 0.0;\n  result.diffuse = true;\n\n  if (dot(normal, view) <= 0.0) {\n    return result;\n  }\n\n  // Lambert, and the bounce of a toon surface: Three lights a toon surface's indirect light with a plain Lambert term.\n  if ((surface.flags & MATERIAL_PBR) == 0u && surface.kind != KIND_PHONG) {\n    result.direction = cosineDirection(normal, u);\n    // The cosine weight cancels the Lambert cosine exactly, which leaves the albedo alone.\n    result.weight = surface.albedo;\n    result.pdf = max(dot(normal, result.direction), 0.0) / PI;\n    result.valid = result.pdf > 0.0;\n\n    return result;\n  }\n\n  let weights = lobeWeights(surface, view);\n\n  result.diffuse = lobe < weights.x;\n\n  if (lobe < weights.x) {\n    result.direction = cosineDirection(normal, u);\n  } else if (surface.kind == KIND_PHONG) {\n    result.direction = reflect(-view, sampleBlinnPhongHalf(normal, surface.shininess, u));\n  } else if (lobe < weights.x + weights.y) {\n    let basis = specularBasis(surface);\n    // The lobe is drawn in its own frame, so the view direction goes into that frame and the result comes back out.\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, normal));\n    let localHalf = sampleVisibleNormal(localView, specularAlphas(surface), u);\n    let half = basis * localHalf;\n\n    result.direction = reflect(-view, half);\n  } else {\n    let coatNormal = surface.clearcoatNormal;\n    let basis = orthonormalBasis(coatNormal);\n    let alpha = max(surface.clearcoatRoughness * surface.clearcoatRoughness, 1e-4);\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, coatNormal));\n    let half = basis * sampleVisibleNormal(localView, vec2<f32>(alpha), u);\n\n    result.direction = reflect(-view, half);\n  }\n\n  let evaluated = evaluateBsdf(surface, view, result.direction);\n\n  if (evaluated.pdf <= 0.0 || all(evaluated.value <= vec3<f32>(0.0))) {\n    return result;\n  }\n\n  result.pdf = evaluated.pdf;\n  result.weight = evaluated.value / evaluated.pdf;\n  result.valid = true;\n\n  return result;\n}\n\n/**\n * Only the diffuse part of evaluateBsdf, times the cosine: what of the light from this direction the denoiser counts\n * as diffuse. The rest of the value is specular.\n */\nfn diffuseValue(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> vec3<f32> {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let nv = max(dot(normal, view), 1e-4);\n\n  if (surface.kind == KIND_TOON) {\n    return surface.albedo / PI * toonRamp(surface.gradientMap, nl);\n  }\n\n  if (nl <= 0.0) {\n    return vec3<f32>(0.0);\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return surface.albedo / PI * nl;\n  }\n\n  var diffuse = surface.albedo * (1.0 - surface.metalness) / PI;\n\n  if (any(surface.sheenColor > vec3<f32>(0.0))) {\n    diffuse = diffuse * (1.0 - max3(surface.sheenColor)\n      * max(sheenAlbedo(nv, surface.sheenRoughness), sheenAlbedo(nl, surface.sheenRoughness)));\n  }\n\n  return diffuse * nl * coatTransmittance(surface, view);\n}\n\n/** How much of the light a surface's specular lobes reflect in all, seen from this view: what that signal is divided by. */\nfn specularAlbedo(surface: Surface, view: vec3<f32>) -> vec3<f32> {\n  let nv = clamp(dot(surface.normal, view), 0.0, 1.0);\n\n  if (surface.kind == KIND_PHONG) {\n    return surface.specularColor * surface.specularStrength;\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return vec3<f32>(0.0);\n  }\n\n  let dfg = dfgApprox(surface.roughness, nv);\n  var albedo = surfaceF0(surface) * dfg.x + vec3<f32>(surface.specularF90 * dfg.y);\n\n  if (surface.clearcoat > 0.0) {\n    let coatDfg = dfgApprox(surface.clearcoatRoughness, clamp(dot(surface.clearcoatNormal, view), 0.0, 1.0));\n    albedo = albedo * coatTransmittance(surface, view) + vec3<f32>(surface.clearcoat * (0.04 * coatDfg.x + coatDfg.y));\n  }\n\n  // What crosses a transmissive surface is the specular signal too, and it is not filtered by any colour of the surface.\n  return mix(albedo, vec3<f32>(1.0), clamp(surface.transmission, 0.0, 1.0));\n}\n";
+export declare const BSDF_WGSL = "\n/** How a surface answers a direction: what it reflects, and how likely sampling was to have chosen it. */\nstruct BsdfEval {\n  // The bidirectional distribution times the cosine at the surface, which is what a contribution is multiplied by.\n  value: vec3<f32>,\n  // Solid angle density of this direction under this renderer's own sampling, for the power heuristic.\n  pdf: f32,\n};\n\nstruct BsdfSample {\n  direction: vec3<f32>,\n  // The distribution over its own density: what the path's throughput is multiplied by.\n  weight: vec3<f32>,\n  pdf: f32,\n  valid: bool,\n  // The diffuse lobe chose the direction, rather than a specular one: which of the denoiser's two signals it feeds.\n  diffuse: bool,\n};\n\nfn luminanceOf(color: vec3<f32>) -> f32 {\n  return dot(color, vec3<f32>(0.2126, 0.7152, 0.0722));\n}\n\nfn max3(color: vec3<f32>) -> f32 {\n  return max(color.r, max(color.g, color.b));\n}\n\n/** Schlick's approximation, as in Three: the reflectance always reaches one at a grazing angle. */\nfn fresnelSchlick(f0: vec3<f32>, cosine: f32) -> vec3<f32> {\n  return schlick3(f0, 1.0, cosine);\n}\n\n/** Reflectance at normal incidence: the dielectric's own, turned towards the base colour as the surface turns metal. */\nfn surfaceF0(surface: Surface) -> vec3<f32> {\n  return mix(surface.specularColor, surface.albedo, surface.metalness);\n}\n\n/** GGX, with alpha the square of the roughness, as Three parameterises it. */\nfn ggxDistribution(alpha2: f32, nh: f32) -> f32 {\n  let denominator = nh * nh * (alpha2 - 1.0) + 1.0;\n\n  return alpha2 / (PI * denominator * denominator);\n}\n\n/** Smith's height correlated visibility, which already carries the 1/(4 nl nv) the lobe would otherwise need. */\nfn ggxVisibility(alpha2: f32, nl: f32, nv: f32) -> f32 {\n  return 0.5 / (nl * sqrt(nv * nv * (1.0 - alpha2) + alpha2) + nv * sqrt(nl * nl * (1.0 - alpha2) + alpha2) + 1e-7);\n}\n\n/** The separable masking term, which is the one the visible normal sampling is built on. */\nfn ggxMasking(alpha2: f32, nv: f32) -> f32 {\n  return 2.0 * nv / (nv + sqrt(alpha2 + (1.0 - alpha2) * nv * nv));\n}\n\n/** Karis' analytic fit of the GGX directional albedo (Three's DFGApprox): the scale and bias applied to F0 and F90. */\nfn dfgApprox(roughness: f32, nv: f32) -> vec2<f32> {\n  let r = roughness * vec4<f32>(-1.0, -0.0275, -0.572, 0.022) + vec4<f32>(1.0, 0.0425, 1.04, -0.04);\n  let a004 = min(r.x * r.x, exp2(-9.28 * nv)) * r.x + r.y;\n\n  return vec2<f32>(-1.04, 1.04) * a004 + r.zw;\n}\n\n/** Three's IBLSheenBRDF: the Charlie sheen lobe integrated over the hemisphere. */\nfn sheenAlbedo(cosine: f32, roughness: f32) -> f32 {\n  let r2 = roughness * roughness;\n  let inverse = 1.0 / (roughness + 0.1);\n  let a = -1.9362 + 1.0678 * roughness + 0.4573 * r2 - 0.8469 * inverse;\n  let b = -0.6014 + 0.5538 * roughness - 0.4670 * r2 - 0.1255 * inverse;\n\n  return clamp(exp(a * clamp(cosine, 0.0, 1.0) + b), 0.0, 1.0);\n}\n\n/** Three's D_Charlie times V_Neubelt. */\nfn sheenLobe(roughness: f32, nh: f32, nv: f32, nl: f32) -> f32 {\n  let inverse = 1.0 / (roughness * roughness);\n  let sin2h = max(1.0 - nh * nh, 0.0078125);\n  let distribution = (2.0 + inverse) * pow(sin2h, inverse * 0.5) / (2.0 * PI);\n\n  return distribution * clamp(1.0 / (4.0 * (nl + nv - nl * nv)), 0.0, 1.0);\n}\n\n/** The frame the specular lobe is defined in: Three's anisotropy direction when there is one, any frame otherwise. */\nfn specularBasis(surface: Surface) -> mat3x3<f32> {\n  let normal = surface.normal;\n\n  if (surface.anisotropy > 0.0) {\n    let t = surface.anisotropyT - normal * dot(normal, surface.anisotropyT);\n    let tangent = normalize(t);\n\n    if (all(tangent == tangent) && dot(t, t) > 1e-12) {\n      return mat3x3<f32>(tangent, cross(normal, tangent), normal);\n    }\n  }\n\n  return orthonormalBasis(normal);\n}\n\n/** The GGX roughness along the frame's two axes: Three's alphaT along the anisotropy, alpha across it. */\nfn specularAlphas(surface: Surface) -> vec2<f32> {\n  let alpha = surface.roughness * surface.roughness;\n\n  if (surface.anisotropy > 0.0) {\n    return vec2<f32>(max(surface.alphaT, 1e-4), max(alpha, 1e-4));\n  }\n\n  return vec2<f32>(max(alpha, 1e-4));\n}\n\n/** Anisotropic GGX in the lobe's own frame; with equal alphas it is exactly the isotropic distribution. */\nfn ggxAnisotropic(h: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let x = h.x / alphas.x;\n  let y = h.y / alphas.y;\n  let d = x * x + y * y + h.z * h.z;\n\n  return 1.0 / (PI * alphas.x * alphas.y * d * d);\n}\n\n/** Smith's masking of one direction, in the lobe's frame. */\nfn ggxMaskingAnisotropic(v: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let projected = alphas.x * alphas.x * v.x * v.x + alphas.y * alphas.y * v.y * v.y;\n\n  return 2.0 * v.z / (v.z + sqrt(v.z * v.z + projected));\n}\n\n/** Three's V_GGX_SmithCorrelated_Anisotropic, in the lobe's frame. */\nfn ggxVisibilityAnisotropic(l: vec3<f32>, v: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let gv = l.z * length(vec3<f32>(alphas.x * v.x, alphas.y * v.y, v.z));\n  let gl = v.z * length(vec3<f32>(alphas.x * l.x, alphas.y * l.y, l.z));\n\n  return 0.5 / max(gv + gl, 1e-7);\n}\n\n/**\n * A half vector drawn from the distribution of normals the view direction can actually see (Heitz 2018), which is what\n * keeps the specular lobe from wasting samples on microfacets facing away. The view is in the lobe's frame.\n */\nfn sampleVisibleNormal(view: vec3<f32>, alphas: vec2<f32>, u: vec2<f32>) -> vec3<f32> {\n  let stretched = normalize(vec3<f32>(alphas.x * view.x, alphas.y * view.y, view.z));\n  let lengthSquared = stretched.x * stretched.x + stretched.y * stretched.y;\n  let t1 = select(\n    vec3<f32>(1.0, 0.0, 0.0),\n    vec3<f32>(-stretched.y, stretched.x, 0.0) * inverseSqrt(max(lengthSquared, 1e-20)),\n    lengthSquared > 1e-20);\n  let t2 = cross(stretched, t1);\n  let radius = sqrt(u.x);\n  let angle = 2.0 * PI * u.y;\n  let px = radius * cos(angle);\n  let raw = radius * sin(angle);\n  let lift = 0.5 * (1.0 + stretched.z);\n  let py = (1.0 - lift) * sqrt(max(0.0, 1.0 - px * px)) + lift * raw;\n  let normal = t1 * px + t2 * py + stretched * sqrt(max(0.0, 1.0 - px * px - py * py));\n\n  return normalize(vec3<f32>(alphas.x * normal.x, alphas.y * normal.y, max(normal.z, 1e-6)));\n}\n\n/** What the clearcoat takes back from everything under it, seen from this view: Three's (1 - clearcoat * Fcc). */\nfn coatTransmittance(surface: Surface, view: vec3<f32>) -> f32 {\n  if (surface.clearcoat <= 0.0) {\n    return 1.0;\n  }\n\n  let nv = clamp(dot(surface.clearcoatNormal, view), 0.0, 1.0);\n\n  return 1.0 - surface.clearcoat * schlick3(vec3<f32>(0.04), 1.0, nv).x;\n}\n\n/**\n * A point inside a volume, where light scatters by the medium's phase function instead of off a surface. Only the\n * integrator makes one: its normal faces back along the ray, its albedo is the medium's, its anisotropy the phase's.\n */\nconst KIND_VOLUME: u32 = 100u;\n\n/** How often each lobe is chosen, as (diffuse, specular, clearcoat): by how much of the reflection each accounts for. */\nfn lobeWeights(surface: Surface, view: vec3<f32>) -> vec3<f32> {\n  let nv = max(dot(surface.normal, view), 1e-4);\n\n  if (surface.kind == KIND_PHONG) {\n    let diffuse = luminanceOf(surface.albedo) + 0.02;\n    let specular = luminanceOf(surface.specularColor) * surface.specularStrength * (surface.shininess * 0.125 + 1.0) + 0.02;\n\n    return vec3<f32>(diffuse, specular, 0.0) / (diffuse + specular);\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return vec3<f32>(1.0, 0.0, 0.0);\n  }\n\n  let base = coatTransmittance(surface, view);\n  // The floors keep every lobe reachable: the diffuse one also carries the sheen and the multiple scattering term.\n  let diffuse = (luminanceOf(surface.albedo * ((1.0 - surface.metalness) * (1.0 - surface.subsurface)))\n    + max3(surface.sheenColor)) * base + 0.02;\n  let specular = luminanceOf(schlick3(surface.specularColor, surface.specularF90, nv)\n    * (1.0 - surface.metalness) + surface.albedo * surface.metalness) * base + 0.05;\n  let coat = select(0.0, (1.0 - base) + 0.02, surface.clearcoat > 0.0);\n  let total = diffuse + specular + coat;\n\n  return vec3<f32>(diffuse, specular, coat) / total;\n}\n\n/** Three's physical direct lighting for one light direction, times the cosine. */\nfn physicalValue(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> vec3<f32> {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let nv = max(dot(normal, view), 1e-4);\n  var value = vec3<f32>(0.0);\n\n  if (nl > 0.0) {\n    let basis = specularBasis(surface);\n    let alphas = specularAlphas(surface);\n    let half = normalize(light + view);\n    let vh = clamp(dot(view, half), 0.0, 1.0);\n    let f0 = surfaceF0(surface);\n    var fresnel = schlick3(f0, surface.specularF90, vh);\n\n    if (surface.iridescence > 0.0) {\n      fresnel = mix(fresnel, surface.iridescenceFresnel, surface.iridescence);\n    }\n\n    let localLight = vec3<f32>(dot(light, basis[0]), dot(light, basis[1]), nl);\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), nv);\n    let localHalf = vec3<f32>(dot(half, basis[0]), dot(half, basis[1]), max(dot(half, normal), 0.0));\n    let single = fresnel * (ggxVisibilityAnisotropic(localLight, localView, alphas) * ggxAnisotropic(localHalf, alphas));\n\n    // Three's BRDF_GGX_Multiscatter: the energy single scattering loses at high roughness, returned.\n    let dfgV = dfgApprox(surface.roughness, nv);\n    let dfgL = dfgApprox(surface.roughness, nl);\n    let essV = f0 * dfgV.x + surface.specularF90 * dfgV.y;\n    let essL = f0 * dfgL.x + surface.specularF90 * dfgL.y;\n    let emsV = 1.0 - (dfgV.x + dfgV.y);\n    let emsL = 1.0 - (dfgL.x + dfgL.y);\n    let average = f0 + (vec3<f32>(1.0) - f0) * 0.047619;\n    let multiple = essV * essL * average / (vec3<f32>(1.0) - emsV * emsL * average + 1e-7) * (emsV * emsL);\n\n    // What scatters inside the material leaves it elsewhere, by the random walk, not from this point.\n    var diffuse = surface.albedo * ((1.0 - surface.metalness) * (1.0 - surface.subsurface)) / PI;\n    var specular = single + multiple;\n\n    // Three's sheen: its lobe sits on top and takes back the energy it reflects from the layers below.\n    if (any(surface.sheenColor > vec3<f32>(0.0))) {\n      let energy = 1.0 - max3(surface.sheenColor)\n        * max(sheenAlbedo(nv, surface.sheenRoughness), sheenAlbedo(nl, surface.sheenRoughness));\n\n      diffuse = diffuse * energy;\n      specular = specular * energy\n        + surface.sheenColor * sheenLobe(surface.sheenRoughness, clamp(dot(normal, half), 0.0, 1.0), nv, nl);\n    }\n\n    value = (diffuse + specular) * nl;\n  }\n\n  if (surface.clearcoat > 0.0) {\n    let coatNormal = surface.clearcoatNormal;\n    value = value * coatTransmittance(surface, view);\n\n    let nlCoat = dot(coatNormal, light);\n    let nvCoat = max(dot(coatNormal, view), 1e-4);\n\n    if (nlCoat > 0.0) {\n      let half = normalize(light + view);\n      let alpha = surface.clearcoatRoughness * surface.clearcoatRoughness;\n      let alpha2 = alpha * alpha;\n      let fresnel = schlick3(vec3<f32>(0.04), 1.0, clamp(dot(view, half), 0.0, 1.0)).x;\n      let lobe = fresnel * ggxVisibility(alpha2, nlCoat, nvCoat) * ggxDistribution(alpha2, max(dot(coatNormal, half), 0.0));\n\n      value = value + vec3<f32>(surface.clearcoat * lobe * nlCoat);\n    }\n  }\n\n  return value;\n}\n\n/** Three's D_BlinnPhong times its implicit visibility of a quarter: the Phong highlight's distribution. */\nfn blinnPhong(shininess: f32, nh: f32) -> f32 {\n  return 0.25 * (shininess * 0.5 + 1.0) / PI * pow(max(nh, 0.0), shininess);\n}\n\n/** Density of the Blinn-Phong lobe's sampling, per solid angle of the light direction. */\nfn blinnPhongPdf(shininess: f32, nh: f32, vh: f32) -> f32 {\n  return (shininess + 1.0) / (2.0 * PI) * pow(max(nh, 0.0), shininess) / max(4.0 * vh, 1e-6);\n}\n\n/** Density with which the GGX lobe (VNDF) would have produced this light direction, in the lobe's frame. */\nfn ggxPdf(localLight: vec3<f32>, localView: vec3<f32>, alphas: vec2<f32>) -> f32 {\n  let half = normalize(localLight + localView);\n\n  if (half.z <= 0.0 || localView.z <= 0.0) {\n    return 0.0;\n  }\n\n  return ggxMaskingAnisotropic(localView, alphas) * ggxAnisotropic(half, alphas) / (4.0 * localView.z);\n}\n\n/** The mixture density of every lobe this surface samples with. */\nfn bsdfPdf(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> f32 {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let weights = lobeWeights(surface, view);\n  var pdf = weights.x * max(nl, 0.0) / PI;\n\n  if (surface.kind == KIND_PHONG) {\n    let half = normalize(light + view);\n\n    return pdf + weights.y * blinnPhongPdf(surface.shininess, dot(normal, half), max(dot(view, half), 0.0));\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return pdf;\n  }\n\n  let basis = specularBasis(surface);\n  let alphas = specularAlphas(surface);\n\n  pdf = pdf + weights.y * ggxPdf(\n    vec3<f32>(dot(light, basis[0]), dot(light, basis[1]), nl),\n    vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, normal)),\n    alphas);\n\n  if (weights.z > 0.0) {\n    let coatBasis = orthonormalBasis(surface.clearcoatNormal);\n    let alpha = max(surface.clearcoatRoughness * surface.clearcoatRoughness, 1e-4);\n\n    pdf = pdf + weights.z * ggxPdf(\n      vec3<f32>(dot(light, coatBasis[0]), dot(light, coatBasis[1]), dot(light, surface.clearcoatNormal)),\n      vec3<f32>(dot(view, coatBasis[0]), dot(view, coatBasis[1]), dot(view, surface.clearcoatNormal)),\n      vec2<f32>(alpha));\n  }\n\n  return pdf;\n}\n\n/**\n * What the surface sends towards the view direction, and the density of that direction. A Lambert material keeps\n * Three's Lambert term alone, which is all its own shading has; a toon material reads its ramp in place of the cosine.\n */\nfn evaluateBsdf(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> BsdfEval {\n  let normal = surface.normal;\n  var result: BsdfEval;\n  result.value = vec3<f32>(0.0);\n  result.pdf = 0.0;\n\n  let nl = dot(normal, light);\n  let nv = dot(normal, view);\n\n  // In a medium every direction is open: the phase function stands for the lobes, and samples itself exactly.\n  if (surface.kind == KIND_VOLUME) {\n    let phase = henyeyGreenstein(dot(-view, light), surface.anisotropy);\n\n    result.value = surface.albedo * phase;\n    result.pdf = phase;\n\n    return result;\n  }\n\n  if (nv <= 0.0) {\n    return result;\n  }\n\n  if (surface.kind == KIND_TOON) {\n    // Three's toon lighting reads its ramp across the whole sphere of directions, the side facing away included.\n    result.value = surface.albedo / PI * toonRamp(surface.gradientMap, nl);\n    result.pdf = max(nl, 0.0) / PI;\n\n    return result;\n  }\n\n  // A clearcoat may face a light its base turns away from, so only a surface without one stops here.\n  if (nl <= 0.0 && surface.clearcoat <= 0.0) {\n    return result;\n  }\n\n  if ((surface.flags & MATERIAL_PBR) != 0u) {\n    result.value = physicalValue(surface, view, light);\n  } else if (surface.kind == KIND_PHONG) {\n    let half = normalize(light + view);\n    let fresnel = schlick3(surface.specularColor, 1.0, clamp(dot(view, half), 0.0, 1.0));\n    let specular = fresnel * blinnPhong(surface.shininess, dot(normal, half)) * surface.specularStrength;\n\n    result.value = (surface.albedo / PI + specular) * max(nl, 0.0);\n  } else {\n    result.value = surface.albedo / PI * max(nl, 0.0);\n  }\n\n  result.pdf = bsdfPdf(surface, view, light);\n\n  return result;\n}\n\n/**\n * Reflectance of a smooth dielectric interface, from both polarisations rather than Schlick's fit: refraction needs the\n * exact cut-off past which nothing gets through. eta is the index the ray leaves over the index it enters.\n */\nfn fresnelDielectric(cosine: f32, eta: f32) -> f32 {\n  let c = clamp(cosine, 0.0, 1.0);\n  let sin2 = eta * eta * (1.0 - c * c);\n\n  if (sin2 >= 1.0) {\n    return 1.0;\n  }\n\n  let t = sqrt(1.0 - sin2);\n  let s = (eta * c - t) / (eta * c + t);\n  let p = (c - eta * t) / (c + eta * t);\n\n  return 0.5 * (s * s + p * p);\n}\n\n/** Where a transmissive surface sends the path: off it, or through it. */\nstruct DielectricSample {\n  direction: vec3<f32>,\n  weight: vec3<f32>,\n  // The path went through the surface instead of reflecting off it.\n  crossed: bool,\n  // It went through without bending, so the density that found the next hit is still the one before this surface.\n  straight: bool,\n  valid: bool,\n};\n\n/**\n * The transmissive share of a physical material. Reflection and refraction are chosen in proportion to the Fresnel\n * reflectance on a microfacet drawn from the visible normals, which leaves only the masking of the outgoing direction\n * in the weight. Light is tinted by the base colour once, on the way in, which is where Three applies it.\n *\n * A thin wall (zero thickness, Three's default) bends nothing: what is not reflected carries straight on, exactly as a\n * shadow ray crossing it does in thinCrossing, so the two estimators agree on a window. ior is the index this path\n * sees, which dispersion makes depend on the colour the path carries.\n */\nfn sampleDielectric(surface: Surface, view: vec3<f32>, ior: f32, lobe: f32, u: vec2<f32>) -> DielectricSample {\n  let normal = surface.normal;\n  var result: DielectricSample;\n  result.valid = false;\n  result.crossed = false;\n  result.straight = false;\n  result.weight = vec3<f32>(0.0);\n\n  let nv = dot(normal, view);\n\n  if (nv <= 0.0) {\n    return result;\n  }\n\n  let alpha = max(surface.roughness * surface.roughness, 1e-4);\n  let alpha2 = alpha * alpha;\n  let basis = orthonormalBasis(normal);\n  let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), nv);\n  let localHalf = sampleVisibleNormal(localView, vec2<f32>(alpha), u);\n  let half = basis[0] * localHalf.x + basis[1] * localHalf.y + normal * localHalf.z;\n\n  if (surface.thickness <= 0.0) {\n    let reflectance = fresnelDielectric(abs(dot(surface.facing, view)), 1.0 / ior);\n\n    if (lobe >= reflectance) {\n      result.direction = -view;\n      result.weight = surface.albedo;\n      result.crossed = true;\n      result.straight = true;\n      result.valid = true;\n\n      return result;\n    }\n\n    result.direction = reflect(-view, half);\n  } else {\n    let eta = select(ior, 1.0 / ior, surface.front);\n    let refracted = refract(-view, half, eta);\n    // Past the critical angle refract() returns nothing, and everything reflects.\n    let reflectance = select(fresnelDielectric(dot(view, half), eta), 1.0, all(refracted == vec3<f32>(0.0)));\n\n    if (lobe >= reflectance) {\n      result.direction = normalize(refracted);\n\n      let no = -dot(normal, result.direction);\n\n      if (no <= 0.0) {\n        return result;\n      }\n\n      result.weight = vec3<f32>(ggxMasking(alpha2, no)) * select(vec3<f32>(1.0), surface.albedo, surface.front);\n      result.crossed = true;\n      result.valid = true;\n\n      return result;\n    }\n\n    result.direction = reflect(-view, half);\n  }\n\n  let nl = dot(normal, result.direction);\n\n  if (nl <= 0.0) {\n    return result;\n  }\n\n  result.weight = vec3<f32>(ggxMasking(alpha2, nl));\n  result.valid = true;\n\n  return result;\n}\n\n/** A half vector around the normal distributed as cos^s, which is the Blinn-Phong lobe's own shape. */\nfn sampleBlinnPhongHalf(normal: vec3<f32>, shininess: f32, u: vec2<f32>) -> vec3<f32> {\n  let cosine = pow(u.x, 1.0 / (shininess + 1.0));\n  let sine = sqrt(max(0.0, 1.0 - cosine * cosine));\n  let angle = 2.0 * PI * u.y;\n\n  return normalize(orthonormalBasis(normal) * vec3<f32>(sine * cos(angle), sine * sin(angle), cosine));\n}\n\n/** Which way the path carries on, and what its throughput is multiplied by. */\nfn sampleBsdf(surface: Surface, view: vec3<f32>, lobe: f32, u: vec2<f32>) -> BsdfSample {\n  let normal = surface.normal;\n  var result: BsdfSample;\n  result.valid = false;\n  result.weight = vec3<f32>(0.0);\n  result.pdf = 0.0;\n  result.diffuse = true;\n\n  if (surface.kind == KIND_VOLUME) {\n    result.direction = sampleHenyeyGreenstein(-view, surface.anisotropy, u);\n    result.weight = surface.albedo;\n    result.pdf = henyeyGreenstein(dot(-view, result.direction), surface.anisotropy);\n    result.valid = true;\n\n    return result;\n  }\n\n  if (dot(normal, view) <= 0.0) {\n    return result;\n  }\n\n  // Lambert, and the bounce of a toon surface: Three lights a toon surface's indirect light with a plain Lambert term.\n  if ((surface.flags & MATERIAL_PBR) == 0u && surface.kind != KIND_PHONG) {\n    result.direction = cosineDirection(normal, u);\n    // The cosine weight cancels the Lambert cosine exactly, which leaves the albedo alone.\n    result.weight = surface.albedo;\n    result.pdf = max(dot(normal, result.direction), 0.0) / PI;\n    result.valid = result.pdf > 0.0;\n\n    return result;\n  }\n\n  let weights = lobeWeights(surface, view);\n\n  result.diffuse = lobe < weights.x;\n\n  if (lobe < weights.x) {\n    result.direction = cosineDirection(normal, u);\n  } else if (surface.kind == KIND_PHONG) {\n    result.direction = reflect(-view, sampleBlinnPhongHalf(normal, surface.shininess, u));\n  } else if (lobe < weights.x + weights.y) {\n    let basis = specularBasis(surface);\n    // The lobe is drawn in its own frame, so the view direction goes into that frame and the result comes back out.\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, normal));\n    let localHalf = sampleVisibleNormal(localView, specularAlphas(surface), u);\n    let half = basis * localHalf;\n\n    result.direction = reflect(-view, half);\n  } else {\n    let coatNormal = surface.clearcoatNormal;\n    let basis = orthonormalBasis(coatNormal);\n    let alpha = max(surface.clearcoatRoughness * surface.clearcoatRoughness, 1e-4);\n    let localView = vec3<f32>(dot(view, basis[0]), dot(view, basis[1]), dot(view, coatNormal));\n    let half = basis * sampleVisibleNormal(localView, vec2<f32>(alpha), u);\n\n    result.direction = reflect(-view, half);\n  }\n\n  let evaluated = evaluateBsdf(surface, view, result.direction);\n\n  if (evaluated.pdf <= 0.0 || all(evaluated.value <= vec3<f32>(0.0))) {\n    return result;\n  }\n\n  result.pdf = evaluated.pdf;\n  result.weight = evaluated.value / evaluated.pdf;\n  result.valid = true;\n\n  return result;\n}\n\n/**\n * Only the diffuse part of evaluateBsdf, times the cosine: what of the light from this direction the denoiser counts\n * as diffuse. The rest of the value is specular.\n */\nfn diffuseValue(surface: Surface, view: vec3<f32>, light: vec3<f32>) -> vec3<f32> {\n  let normal = surface.normal;\n  let nl = dot(normal, light);\n  let nv = max(dot(normal, view), 1e-4);\n\n  if (surface.kind == KIND_TOON) {\n    return surface.albedo / PI * toonRamp(surface.gradientMap, nl);\n  }\n\n  // Light scattered in a medium is all diffuse, as the denoiser sees it.\n  if (surface.kind == KIND_VOLUME) {\n    return surface.albedo * henyeyGreenstein(dot(-view, light), surface.anisotropy);\n  }\n\n  if (nl <= 0.0) {\n    return vec3<f32>(0.0);\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return surface.albedo / PI * nl;\n  }\n\n  var diffuse = surface.albedo * ((1.0 - surface.metalness) * (1.0 - surface.subsurface)) / PI;\n\n  if (any(surface.sheenColor > vec3<f32>(0.0))) {\n    diffuse = diffuse * (1.0 - max3(surface.sheenColor)\n      * max(sheenAlbedo(nv, surface.sheenRoughness), sheenAlbedo(nl, surface.sheenRoughness)));\n  }\n\n  return diffuse * nl * coatTransmittance(surface, view);\n}\n\n/** How much of the light a surface's specular lobes reflect in all, seen from this view: what that signal is divided by. */\nfn specularAlbedo(surface: Surface, view: vec3<f32>) -> vec3<f32> {\n  let nv = clamp(dot(surface.normal, view), 0.0, 1.0);\n\n  if (surface.kind == KIND_PHONG) {\n    return surface.specularColor * surface.specularStrength;\n  }\n\n  if ((surface.flags & MATERIAL_PBR) == 0u) {\n    return vec3<f32>(0.0);\n  }\n\n  let dfg = dfgApprox(surface.roughness, nv);\n  var albedo = surfaceF0(surface) * dfg.x + vec3<f32>(surface.specularF90 * dfg.y);\n\n  if (surface.clearcoat > 0.0) {\n    let coatDfg = dfgApprox(surface.clearcoatRoughness, clamp(dot(surface.clearcoatNormal, view), 0.0, 1.0));\n    albedo = albedo * coatTransmittance(surface, view) + vec3<f32>(surface.clearcoat * (0.04 * coatDfg.x + coatDfg.y));\n  }\n\n  // What crosses a transmissive surface is the specular signal too, and it is not filtered by any colour of the surface.\n  return mix(albedo, vec3<f32>(1.0), clamp(surface.transmission, 0.0, 1.0));\n}\n";
 //# sourceMappingURL=bsdf.wgsl.d.ts.map

package/dist/wgsl/bsdf.wgsl.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"bsdf.wgsl.d.ts","sourceRoot":"","sources":["../../src/wgsl/bsdf.wgsl.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,eAAO,MAAM,SAAS,s8tBA8lBrB,CAAC"}
+{"version":3,"file":"bsdf.wgsl.d.ts","sourceRoot":"","sources":["../../src/wgsl/bsdf.wgsl.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;GAUG;AACH,eAAO,MAAM,SAAS,+twBA8nBrB,CAAC"}