softbox 4.0.0

package/dist/ThreePathTracingService-Bv1CmdCK.cjs

@@ -0,0 +1,3977 @@
+"use strict";Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const l=require("three"),T=require("./index-C4Chx9Wl.cjs");function $e(s){const e=Object.create(null,{[Symbol.toStringTag]:{value:"Module"}});if(s){for(const t in s)if(t!=="default"){const r=Object.getOwnPropertyDescriptor(s,t);Object.defineProperty(e,t,r.get?r:{enumerable:!0,get:()=>s[t]})}}return e.default=s,Object.freeze(e)}const A=$e(l);function je(s){switch(s){case 1:return"R";case 2:return"RG";case 3:return"RGBA";case 4:return"RGBA"}throw new Error}function Ye(s){switch(s){case 1:return l.RedFormat;case 2:return l.RGFormat;case 3:return l.RGBAFormat;case 4:return l.RGBAFormat}}function de(s){switch(s){case 1:return l.RedIntegerFormat;case 2:return l.RGIntegerFormat;case 3:return l.RGBAIntegerFormat;case 4:return l.RGBAIntegerFormat}}class Oe extends l.DataTexture{constructor(){super(),this.minFilter=l.NearestFilter,this.magFilter=l.NearestFilter,this.generateMipmaps=!1,this.overrideItemSize=null,this._forcedType=null}updateFrom(e){const t=this.overrideItemSize,r=e.itemSize,i=e.count;if(t!==null){if(r*i%t!==0)throw new Error("VertexAttributeTexture: overrideItemSize must divide evenly into buffer length.");e.itemSize=t,e.count=i*r/t}const o=e.itemSize,a=e.count,c=e.normalized,h=e.array.constructor,d=h.BYTES_PER_ELEMENT;let m=this._forcedType,f=o;if(m===null)switch(h){case Float32Array:m=l.FloatType;break;case Uint8Array:case Uint16Array:case Uint32Array:m=l.UnsignedIntType;break;case Int8Array:case Int16Array:case Int32Array:m=l.IntType;break}let n,x,p,v,u=je(o);switch(m){case l.FloatType:p=1,x=Ye(o),c&&d===1?(v=h,u+="8",h===Uint8Array?n=l.UnsignedByteType:(n=l.ByteType,u+="_SNORM")):(v=Float32Array,u+="32F",n=l.FloatType);break;case l.IntType:u+=d*8+"I",p=c?Math.pow(2,h.BYTES_PER_ELEMENT*8-1):1,x=de(o),d===1?(v=Int8Array,n=l.ByteType):d===2?(v=Int16Array,n=l.ShortType):(v=Int32Array,n=l.IntType);break;case l.UnsignedIntType:u+=d*8+"UI",p=c?Math.pow(2,h.BYTES_PER_ELEMENT*8-1):1,x=de(o),d===1?(v=Uint8Array,n=l.UnsignedByteType):d===2?(v=Uint16Array,n=l.UnsignedShortType):(v=Uint32Array,n=l.UnsignedIntType);break}f===3&&(x===l.RGBAFormat||x===l.RGBAIntegerFormat)&&(f=4);const g=Math.ceil(Math.sqrt(a))||1,y=f*g*g,b=new v(y),w=e.normalized;e.normalized=!1;for(let R=0;R<a;R++){const _=f*R;b[_]=e.getX(R)/p,o>=2&&(b[_+1]=e.getY(R)/p),o>=3&&(b[_+2]=e.getZ(R)/p,f===4&&(b[_+3]=1)),o>=4&&(b[_+3]=e.getW(R)/p)}e.normalized=w,this.internalFormat=u,this.format=x,this.type=n,this.image.width=g,this.image.height=g,this.image.data=b,this.needsUpdate=!0,this.dispose(),e.itemSize=r,e.count=i}}class Be extends Oe{constructor(){super(),this._forcedType=l.UnsignedIntType}}class Ne extends Oe{constructor(){super(),this._forcedType=l.FloatType}}class Xe{constructor(){this.index=new Be,this.position=new Ne,this.bvhBounds=new l.DataTexture,this.bvhContents=new l.DataTexture,this._cachedIndexAttr=null,this.index.overrideItemSize=3}updateFrom(e){const{geometry:t}=e;if(Ke(e,this.bvhBounds,this.bvhContents),this.position.updateFrom(t.attributes.position),e.indirect){const r=e._indirectBuffer;if(this._cachedIndexAttr===null||this._cachedIndexAttr.count!==r.length)if(t.index)this._cachedIndexAttr=t.index.clone();else{const i=T.getIndexArray(T.getVertexCount(t));this._cachedIndexAttr=new l.BufferAttribute(i,1,!1)}Qe(t,r,this._cachedIndexAttr),this.index.updateFrom(this._cachedIndexAttr)}else this.index.updateFrom(t.index)}dispose(){const{index:e,position:t,bvhBounds:r,bvhContents:i}=this;e&&e.dispose(),t&&t.dispose(),r&&r.dispose(),i&&i.dispose()}}function Qe(s,e,t){const r=t.array,i=s.index?s.index.array:null;for(let o=0,a=e.length;o<a;o++){const c=3*o,h=3*e[o];for(let d=0;d<3;d++)r[c+d]=i?i[h+d]:h+d}}function Ke(s,e,t){const r=s._roots;if(r.length!==1)throw new Error("MeshBVHUniformStruct: Multi-root BVHs not supported.");const i=r[0],o=new Uint16Array(i),a=new Uint32Array(i),c=new Float32Array(i),h=i.byteLength/T.BYTES_PER_NODE,d=2*Math.ceil(Math.sqrt(h/2)),m=new Float32Array(4*d*d),f=Math.ceil(Math.sqrt(h)),n=new Uint32Array(2*f*f);for(let x=0;x<h;x++){const p=x*T.BYTES_PER_NODE/4,v=p*2,u=T.BOUNDING_DATA_INDEX(p);for(let g=0;g<3;g++)m[8*x+0+g]=c[u+0+g],m[8*x+4+g]=c[u+3+g];if(T.IS_LEAF(v,o)){const g=T.COUNT(v,o),y=T.OFFSET(p,a),b=T.LEAFNODE_MASK_32|g;n[x*2+0]=b,n[x*2+1]=y}else{const g=a[p+6],y=T.SPLIT_AXIS(p,a);n[x*2+0]=y,n[x*2+1]=g}}e.image.data=m,e.image.width=d,e.image.height=d,e.format=l.RGBAFormat,e.type=l.FloatType,e.internalFormat="RGBA32F",e.minFilter=l.NearestFilter,e.magFilter=l.NearestFilter,e.generateMipmaps=!1,e.needsUpdate=!0,e.dispose(),t.image.data=n,t.image.width=f,t.image.height=f,t.format=l.RGIntegerFormat,t.type=l.UnsignedIntType,t.internalFormat="RG32UI",t.minFilter=l.NearestFilter,t.magFilter=l.NearestFilter,t.generateMipmaps=!1,t.needsUpdate=!0,t.dispose()}const Ze=`
+
+// A stack of uint32 indices can can store the indices for
+// a perfectly balanced tree with a depth up to 31. Lower stack
+// depth gets higher performance.
+//
+// However not all trees are balanced. Best value to set this to
+// is the trees max depth.
+#ifndef BVH_STACK_DEPTH
+#define BVH_STACK_DEPTH 60
+#endif
+
+#ifndef INFINITY
+#define INFINITY 1e20
+#endif
+
+// Utilities
+uvec4 uTexelFetch1D( usampler2D tex, uint index ) {
+
+	uint width = uint( textureSize( tex, 0 ).x );
+	uvec2 uv;
+	uv.x = index % width;
+	uv.y = index / width;
+
+	return texelFetch( tex, ivec2( uv ), 0 );
+
+}
+
+ivec4 iTexelFetch1D( isampler2D tex, uint index ) {
+
+	uint width = uint( textureSize( tex, 0 ).x );
+	uvec2 uv;
+	uv.x = index % width;
+	uv.y = index / width;
+
+	return texelFetch( tex, ivec2( uv ), 0 );
+
+}
+
+vec4 texelFetch1D( sampler2D tex, uint index ) {
+
+	uint width = uint( textureSize( tex, 0 ).x );
+	uvec2 uv;
+	uv.x = index % width;
+	uv.y = index / width;
+
+	return texelFetch( tex, ivec2( uv ), 0 );
+
+}
+
+vec4 textureSampleBarycoord( sampler2D tex, vec3 barycoord, uvec3 faceIndices ) {
+
+	return
+		barycoord.x * texelFetch1D( tex, faceIndices.x ) +
+		barycoord.y * texelFetch1D( tex, faceIndices.y ) +
+		barycoord.z * texelFetch1D( tex, faceIndices.z );
+
+}
+
+void ndcToCameraRay(
+	vec2 coord, mat4 cameraWorld, mat4 invProjectionMatrix,
+	out vec3 rayOrigin, out vec3 rayDirection
+) {
+
+	// get camera look direction and near plane for camera clipping
+	vec4 lookDirection = cameraWorld * vec4( 0.0, 0.0, - 1.0, 0.0 );
+	vec4 nearVector = invProjectionMatrix * vec4( 0.0, 0.0, - 1.0, 1.0 );
+	float near = abs( nearVector.z / nearVector.w );
+
+	// get the camera direction and position from camera matrices
+	vec4 origin = cameraWorld * vec4( 0.0, 0.0, 0.0, 1.0 );
+	vec4 direction = invProjectionMatrix * vec4( coord, 0.5, 1.0 );
+	direction /= direction.w;
+	direction = cameraWorld * direction - origin;
+
+	// slide the origin along the ray until it sits at the near clip plane position
+	origin.xyz += direction.xyz * near / dot( direction, lookDirection );
+
+	rayOrigin = origin.xyz;
+	rayDirection = direction.xyz;
+
+}
+`,Je=`
+
+#ifndef TRI_INTERSECT_EPSILON
+#define TRI_INTERSECT_EPSILON 1e-5
+#endif
+
+// Raycasting
+bool intersectsBounds( vec3 rayOrigin, vec3 rayDirection, vec3 boundsMin, vec3 boundsMax, out float dist ) {
+
+	// https://www.reddit.com/r/opengl/comments/8ntzz5/fast_glsl_ray_box_intersection/
+	// https://tavianator.com/2011/ray_box.html
+	vec3 invDir = 1.0 / rayDirection;
+
+	// find intersection distances for each plane
+	vec3 tMinPlane = invDir * ( boundsMin - rayOrigin );
+	vec3 tMaxPlane = invDir * ( boundsMax - rayOrigin );
+
+	// get the min and max distances from each intersection
+	vec3 tMinHit = min( tMaxPlane, tMinPlane );
+	vec3 tMaxHit = max( tMaxPlane, tMinPlane );
+
+	// get the furthest hit distance
+	vec2 t = max( tMinHit.xx, tMinHit.yz );
+	float t0 = max( t.x, t.y );
+
+	// get the minimum hit distance
+	t = min( tMaxHit.xx, tMaxHit.yz );
+	float t1 = min( t.x, t.y );
+
+	// set distance to 0.0 if the ray starts inside the box
+	dist = max( t0, 0.0 );
+
+	return t1 >= dist;
+
+}
+
+bool intersectsTriangle(
+	vec3 rayOrigin, vec3 rayDirection, vec3 a, vec3 b, vec3 c,
+	out vec3 barycoord, out vec3 norm, out float dist, out float side
+) {
+
+	// https://stackoverflow.com/questions/42740765/intersection-between-line-and-triangle-in-3d
+	vec3 edge1 = b - a;
+	vec3 edge2 = c - a;
+	norm = cross( edge1, edge2 );
+
+	float det = - dot( rayDirection, norm );
+	float invdet = 1.0 / det;
+
+	vec3 AO = rayOrigin - a;
+	vec3 DAO = cross( AO, rayDirection );
+
+	vec4 uvt;
+	uvt.x = dot( edge2, DAO ) * invdet;
+	uvt.y = - dot( edge1, DAO ) * invdet;
+	uvt.z = dot( AO, norm ) * invdet;
+	uvt.w = 1.0 - uvt.x - uvt.y;
+
+	// set the hit information
+	barycoord = uvt.wxy; // arranged in A, B, C order
+	dist = uvt.z;
+	side = sign( det );
+	norm = side * normalize( norm );
+
+	// add an epsilon to avoid misses between triangles
+	uvt += vec4( TRI_INTERSECT_EPSILON );
+
+	return all( greaterThanEqual( uvt, vec4( 0.0 ) ) );
+
+}
+
+bool intersectTriangles(
+	// geometry info and triangle range
+	sampler2D positionAttr, usampler2D indexAttr, uint offset, uint count,
+
+	// ray
+	vec3 rayOrigin, vec3 rayDirection,
+
+	// outputs
+	inout float minDistance, inout uvec4 faceIndices, inout vec3 faceNormal, inout vec3 barycoord,
+	inout float side, inout float dist
+) {
+
+	bool found = false;
+	vec3 localBarycoord, localNormal;
+	float localDist, localSide;
+	for ( uint i = offset, l = offset + count; i < l; i ++ ) {
+
+		uvec3 indices = uTexelFetch1D( indexAttr, i ).xyz;
+		vec3 a = texelFetch1D( positionAttr, indices.x ).rgb;
+		vec3 b = texelFetch1D( positionAttr, indices.y ).rgb;
+		vec3 c = texelFetch1D( positionAttr, indices.z ).rgb;
+
+		if (
+			intersectsTriangle( rayOrigin, rayDirection, a, b, c, localBarycoord, localNormal, localDist, localSide )
+			&& localDist < minDistance
+		) {
+
+			found = true;
+			minDistance = localDist;
+
+			faceIndices = uvec4( indices.xyz, i );
+			faceNormal = localNormal;
+
+			side = localSide;
+			barycoord = localBarycoord;
+			dist = localDist;
+
+		}
+
+	}
+
+	return found;
+
+}
+
+bool intersectsBVHNodeBounds( vec3 rayOrigin, vec3 rayDirection, sampler2D bvhBounds, uint currNodeIndex, out float dist ) {
+
+	uint cni2 = currNodeIndex * 2u;
+	vec3 boundsMin = texelFetch1D( bvhBounds, cni2 ).xyz;
+	vec3 boundsMax = texelFetch1D( bvhBounds, cni2 + 1u ).xyz;
+	return intersectsBounds( rayOrigin, rayDirection, boundsMin, boundsMax, dist );
+
+}
+
+// use a macro to hide the fact that we need to expand the struct into separate fields
+#define	bvhIntersectFirstHit(		bvh,		rayOrigin, rayDirection, faceIndices, faceNormal, barycoord, side, dist	)	_bvhIntersectFirstHit(		bvh.position, bvh.index, bvh.bvhBounds, bvh.bvhContents,		rayOrigin, rayDirection, faceIndices, faceNormal, barycoord, side, dist	)
+
+bool _bvhIntersectFirstHit(
+	// bvh info
+	sampler2D bvh_position, usampler2D bvh_index, sampler2D bvh_bvhBounds, usampler2D bvh_bvhContents,
+
+	// ray
+	vec3 rayOrigin, vec3 rayDirection,
+
+	// output variables split into separate variables due to output precision
+	inout uvec4 faceIndices, inout vec3 faceNormal, inout vec3 barycoord,
+	inout float side, inout float dist
+) {
+
+	// stack needs to be twice as long as the deepest tree we expect because
+	// we push both the left and right child onto the stack every traversal
+	int pointer = 0;
+	uint stack[ BVH_STACK_DEPTH ];
+	stack[ 0 ] = 0u;
+
+	float triangleDistance = INFINITY;
+	bool found = false;
+	while ( pointer > - 1 && pointer < BVH_STACK_DEPTH ) {
+
+		uint currNodeIndex = stack[ pointer ];
+		pointer --;
+
+		// check if we intersect the current bounds
+		float boundsHitDistance;
+		if (
+			! intersectsBVHNodeBounds( rayOrigin, rayDirection, bvh_bvhBounds, currNodeIndex, boundsHitDistance )
+			|| boundsHitDistance > triangleDistance
+		) {
+
+			continue;
+
+		}
+
+		uvec2 boundsInfo = uTexelFetch1D( bvh_bvhContents, currNodeIndex ).xy;
+		bool isLeaf = bool( boundsInfo.x & 0xffff0000u );
+
+		if ( isLeaf ) {
+
+			uint count = boundsInfo.x & 0x0000ffffu;
+			uint offset = boundsInfo.y;
+
+			found = intersectTriangles(
+				bvh_position, bvh_index, offset, count,
+				rayOrigin, rayDirection, triangleDistance,
+				faceIndices, faceNormal, barycoord, side, dist
+			) || found;
+
+		} else {
+
+			uint leftIndex = currNodeIndex + 1u;
+			uint splitAxis = boundsInfo.x & 0x0000ffffu;
+			uint rightIndex = currNodeIndex + boundsInfo.y;
+
+			bool leftToRight = rayDirection[ splitAxis ] >= 0.0;
+			uint c1 = leftToRight ? leftIndex : rightIndex;
+			uint c2 = leftToRight ? rightIndex : leftIndex;
+
+			// set c2 in the stack so we traverse it later. We need to keep track of a pointer in
+			// the stack while we traverse. The second pointer added is the one that will be
+			// traversed first
+			pointer ++;
+			stack[ pointer ] = c2;
+
+			pointer ++;
+			stack[ pointer ] = c1;
+
+		}
+
+	}
+
+	return found;
+
+}
+`,et=`
+struct BVH {
+
+	usampler2D index;
+	sampler2D position;
+
+	sampler2D bvhBounds;
+	usampler2D bvhContents;
+
+};
+`;function ze(s,e,t=0){if(s.isInterleavedBufferAttribute){const r=s.itemSize;for(let i=0,o=s.count;i<o;i++){const a=i+t;e.setX(a,s.getX(i)),r>=2&&e.setY(a,s.getY(i)),r>=3&&e.setZ(a,s.getZ(i)),r>=4&&e.setW(a,s.getW(i))}}else{const r=e.array,i=r.constructor,o=r.BYTES_PER_ELEMENT*s.itemSize*t;new i(r.buffer,o,s.array.length).set(s.array)}}function H(s,e=null){const t=s.array.constructor,r=s.normalized,i=s.itemSize,o=e===null?s.count:e;return new l.BufferAttribute(new t(i*o),i,r)}function N(s,e){if(!s&&!e)return!0;if(!!s!=!!e)return!1;const t=s.count===e.count,r=s.normalized===e.normalized,i=s.array.constructor===e.array.constructor,o=s.itemSize===e.itemSize;return!(!t||!r||!i||!o)}function tt(s){const e=s[0].index!==null,t=new Set(Object.keys(s[0].attributes));if(!s[0].getAttribute("position"))throw new Error("StaticGeometryGenerator: position attribute is required.");for(let r=0;r<s.length;++r){const i=s[r];let o=0;if(e!==(i.index!==null))throw new Error("StaticGeometryGenerator: All geometries must have compatible attributes; make sure index attribute exists among all geometries, or in none of them.");for(const a in i.attributes){if(!t.has(a))throw new Error('StaticGeometryGenerator: All geometries must have compatible attributes; make sure "'+a+'" attribute exists among all geometries, or in none of them.');o++}if(o!==t.size)throw new Error("StaticGeometryGenerator: All geometries must have the same number of attributes.")}}function rt(s){let e=0;for(let t=0,r=s.length;t<r;t++)e+=s[t].getIndex().count;return e}function at(s){let e=0;for(let t=0,r=s.length;t<r;t++)e+=s[t].getAttribute("position").count;return e}function it(s,e,t){s.index&&s.index.count!==e&&s.setIndex(null);const r=s.attributes;for(const i in r)r[i].count!==t&&s.deleteAttribute(i)}function st(s,e={},t=new l.BufferGeometry){const{useGroups:r=!1,forceUpdate:i=!1,skipAssigningAttributes:o=[],overwriteIndex:a=!0}=e;tt(s);const c=s[0].index!==null,h=c?rt(s):-1,d=at(s);if(it(t,h,d),r){let f=0;for(let n=0,x=s.length;n<x;n++){const p=s[n];let v;c?v=p.getIndex().count:v=p.getAttribute("position").count,t.addGroup(f,v,n),f+=v}}if(c){let f=!1;if(t.index||(t.setIndex(new l.BufferAttribute(new Uint32Array(h),1,!1)),f=!0),f||a){let n=0,x=0;const p=t.getIndex();for(let v=0,u=s.length;v<u;v++){const g=s[v],y=g.getIndex();if(!(!i&&!f&&o[v]))for(let w=0;w<y.count;++w)p.setX(n+w,y.getX(w)+x);n+=y.count,x+=g.getAttribute("position").count}}}const m=Object.keys(s[0].attributes);for(let f=0,n=m.length;f<n;f++){let x=!1;const p=m[f];if(!t.getAttribute(p)){const g=s[0].getAttribute(p);t.setAttribute(p,H(g,d)),x=!0}let v=0;const u=t.getAttribute(p);for(let g=0,y=s.length;g<y;g++){const b=s[g],w=!i&&!x&&o[g],R=b.getAttribute(p);w||ze(R,u,v),v+=R.count}}}function ot(s,e,t){const r=s.index,o=s.attributes.position.count,a=r?r.count:o;let c=s.groups;c.length===0&&(c=[{count:a,start:0,materialIndex:0}]);let h=s.getAttribute("materialIndex");if(!h||h.count!==o){let m;t.length<=255?m=new Uint8Array(o):m=new Uint16Array(o),h=new l.BufferAttribute(m,1,!1),s.deleteAttribute("materialIndex"),s.setAttribute("materialIndex",h)}const d=h.array;for(let m=0;m<c.length;m++){const f=c[m],n=f.start,x=f.count,p=Math.min(x,a-n),v=Array.isArray(e)?e[f.materialIndex]:e,u=t.indexOf(v);for(let g=0;g<p;g++){let y=n+g;r&&(y=r.getX(y)),d[y]=u}}}function nt(s,e){if(!s.index){const t=s.attributes.position.count,r=new Array(t);for(let i=0;i<t;i++)r[i]=i;s.setIndex(r)}if(!s.attributes.normal&&e&&e.includes("normal")&&s.computeVertexNormals(),!s.attributes.uv&&e&&e.includes("uv")){const t=s.attributes.position.count;s.setAttribute("uv",new l.BufferAttribute(new Float32Array(t*2),2,!1))}if(!s.attributes.uv2&&e&&e.includes("uv2")){const t=s.attributes.position.count;s.setAttribute("uv2",new l.BufferAttribute(new Float32Array(t*2),2,!1))}if(!s.attributes.tangent&&e&&e.includes("tangent"))if(s.attributes.uv&&s.attributes.normal)s.computeTangents();else{const t=s.attributes.position.count;s.setAttribute("tangent",new l.BufferAttribute(new Float32Array(t*4),4,!1))}if(!s.attributes.color&&e&&e.includes("color")){const t=s.attributes.position.count,r=new Float32Array(t*4);r.fill(1),s.setAttribute("color",new l.BufferAttribute(r,4))}}function se(s){let e=0;if(s.byteLength!==0){const t=new Uint8Array(s);for(let r=0;r<s.byteLength;r++){const i=t[r];e=(e<<5)-e+i,e|=0}}return e}function me(s){let e=s.uuid;const t=Object.values(s.attributes);s.index&&(t.push(s.index),e+=`index|${s.index.version}`);const r=Object.keys(t).sort();for(const i of r){const o=t[i];e+=`${i}_${o.version}|`}return e}function pe(s){const e=s.skeleton;return e?(e.boneTexture||e.computeBoneTexture(),`${se(e.boneTexture.image.data.buffer)}_${e.boneTexture.uuid}`):null}class lt{constructor(e=null){this.matrixWorld=new l.Matrix4,this.geometryHash=null,this.skeletonHash=null,this.primitiveCount=-1,e!==null&&this.updateFrom(e)}updateFrom(e){const t=e.geometry,r=(t.index?t.index.count:t.attributes.position.count)/3;this.matrixWorld.copy(e.matrixWorld),this.geometryHash=me(t),this.primitiveCount=r,this.skeletonHash=pe(e)}didChange(e){const t=e.geometry,r=(t.index?t.index.count:t.attributes.position.count)/3;return!(this.matrixWorld.equals(e.matrixWorld)&&this.geometryHash===me(t)&&this.skeletonHash===pe(e)&&this.primitiveCount===r)}}const M=new l.Vector3,C=new l.Vector3,D=new l.Vector3,ge=new l.Vector4,U=new l.Vector3,J=new l.Vector3,ve=new l.Vector4,xe=new l.Vector4,G=new l.Matrix4,ye=new l.Matrix4;function be(s,e,t){const r=s.skeleton,i=s.geometry,o=r.bones,a=r.boneInverses;ve.fromBufferAttribute(i.attributes.skinIndex,e),xe.fromBufferAttribute(i.attributes.skinWeight,e),G.elements.fill(0);for(let c=0;c<4;c++){const h=xe.getComponent(c);if(h!==0){const d=ve.getComponent(c);ye.multiplyMatrices(o[d].matrixWorld,a[d]),ct(G,ye,h)}}return G.multiply(s.bindMatrix).premultiply(s.bindMatrixInverse),t.transformDirection(G),t}function ee(s,e,t,r,i){U.set(0,0,0);for(let o=0,a=s.length;o<a;o++){const c=e[o],h=s[o];c!==0&&(J.fromBufferAttribute(h,r),t?U.addScaledVector(J,c):U.addScaledVector(J.sub(i),c))}i.add(U)}function ct(s,e,t){const r=s.elements,i=e.elements;for(let o=0,a=i.length;o<a;o++)r[o]+=i[o]*t}function ut(s){const{index:e,attributes:t}=s;if(e)for(let r=0,i=e.count;r<i;r+=3){const o=e.getX(r),a=e.getX(r+2);e.setX(r,a),e.setX(r+2,o)}else for(const r in t){const i=t[r],o=i.itemSize;for(let a=0,c=i.count;a<c;a+=3)for(let h=0;h<o;h++){const d=i.getComponent(a,h),m=i.getComponent(a+2,h);i.setComponent(a,h,m),i.setComponent(a+2,h,d)}}return s}function ft(s,e={},t=new l.BufferGeometry){e={applyWorldTransforms:!0,attributes:[],...e};const r=s.geometry,i=e.applyWorldTransforms,o=e.attributes.includes("normal"),a=e.attributes.includes("tangent"),c=r.attributes,h=t.attributes;for(const y in t.attributes)(!e.attributes.includes(y)||!(y in r.attributes))&&t.deleteAttribute(y);!t.index&&r.index&&(t.index=r.index.clone()),h.position||t.setAttribute("position",H(c.position)),o&&!h.normal&&c.normal&&t.setAttribute("normal",H(c.normal)),a&&!h.tangent&&c.tangent&&t.setAttribute("tangent",H(c.tangent)),N(r.index,t.index),N(c.position,h.position),o&&N(c.normal,h.normal),a&&N(c.tangent,h.tangent);const d=c.position,m=o?c.normal:null,f=a?c.tangent:null,n=r.morphAttributes.position,x=r.morphAttributes.normal,p=r.morphAttributes.tangent,v=r.morphTargetsRelative,u=s.morphTargetInfluences,g=new l.Matrix3;g.getNormalMatrix(s.matrixWorld),r.index&&t.index.array.set(r.index.array);for(let y=0,b=c.position.count;y<b;y++)M.fromBufferAttribute(d,y),m&&C.fromBufferAttribute(m,y),f&&(ge.fromBufferAttribute(f,y),D.fromBufferAttribute(f,y)),u&&(n&&ee(n,u,v,y,M),x&&ee(x,u,v,y,C),p&&ee(p,u,v,y,D)),s.isSkinnedMesh&&(s.applyBoneTransform(y,M),m&&be(s,y,C),f&&be(s,y,D)),i&&M.applyMatrix4(s.matrixWorld),h.position.setXYZ(y,M.x,M.y,M.z),m&&(i&&C.applyNormalMatrix(g),h.normal.setXYZ(y,C.x,C.y,C.z)),f&&(i&&D.transformDirection(s.matrixWorld),h.tangent.setXYZW(y,D.x,D.y,D.z,ge.w));for(const y in e.attributes){const b=e.attributes[y];b==="position"||b==="tangent"||b==="normal"||!(b in c)||(h[b]||t.setAttribute(b,H(c[b])),N(c[b],h[b]),ze(c[b],h[b]))}return s.matrixWorld.determinant()<0&&ut(t),t}class ht extends l.BufferGeometry{constructor(){super(),this.version=0,this.hash=null,this._diff=new lt}isCompatible(e,t){const r=e.geometry;for(let i=0;i<t.length;i++){const o=t[i],a=r.attributes[o],c=this.attributes[o];if(a&&!N(a,c))return!1}return!0}updateFrom(e,t){const r=this._diff;return r.didChange(e)?(ft(e,t,this),r.updateFrom(e),this.version++,this.hash=`${this.uuid}_${this.version}`,!0):!1}}const ae=0,He=1,We=2;function dt(s,e){for(let t=0,r=s.length;t<r;t++)s[t].traverseVisible(o=>{o.isMesh&&e(o)})}function mt(s){const e=[];for(let t=0,r=s.length;t<r;t++){const i=s[t];Array.isArray(i.material)?e.push(...i.material):e.push(i.material)}return e}function pt(s,e,t){if(s.length===0){e.setIndex(null);const r=e.attributes;for(const i in r)e.deleteAttribute(i);for(const i in t.attributes)e.setAttribute(t.attributes[i],new l.BufferAttribute(new Float32Array(0),4,!1))}else st(s,t,e);for(const r in e.attributes)e.attributes[r].needsUpdate=!0}class gt{constructor(e){this.objects=null,this.useGroups=!0,this.applyWorldTransforms=!0,this.generateMissingAttributes=!0,this.overwriteIndex=!0,this.attributes=["position","normal","color","tangent","uv","uv2"],this._intermediateGeometry=new Map,this._geometryMergeSets=new WeakMap,this._mergeOrder=[],this._dummyMesh=null,this.setObjects(e||[])}_getDummyMesh(){if(!this._dummyMesh){const e=new l.MeshBasicMaterial,t=new l.BufferGeometry;t.setAttribute("position",new l.BufferAttribute(new Float32Array(9),3)),this._dummyMesh=new l.Mesh(t,e)}return this._dummyMesh}_getMeshes(){const e=[];return dt(this.objects,t=>{e.push(t)}),e.sort((t,r)=>t.uuid>r.uuid?1:t.uuid<r.uuid?-1:0),e.length===0&&e.push(this._getDummyMesh()),e}_updateIntermediateGeometries(){const{_intermediateGeometry:e}=this,t=this._getMeshes(),r=new Set(e.keys()),i={attributes:this.attributes,applyWorldTransforms:this.applyWorldTransforms};for(let o=0,a=t.length;o<a;o++){const c=t[o],h=c.uuid;r.delete(h);let d=e.get(h);(!d||!d.isCompatible(c,this.attributes))&&(d&&d.dispose(),d=new ht,e.set(h,d)),d.updateFrom(c,i)&&this.generateMissingAttributes&&nt(d,this.attributes)}r.forEach(o=>{e.delete(o)})}setObjects(e){Array.isArray(e)?this.objects=[...e]:this.objects=[e]}generate(e=new l.BufferGeometry){const{useGroups:t,overwriteIndex:r,_intermediateGeometry:i,_geometryMergeSets:o}=this,a=this._getMeshes(),c=[],h=[],d=o.get(e)||[];this._updateIntermediateGeometries();let m=!1;a.length!==d.length&&(m=!0);for(let n=0,x=a.length;n<x;n++){const p=a[n],v=i.get(p.uuid);h.push(v);const u=d[n];!u||u.uuid!==v.uuid?(c.push(!1),m=!0):u.version!==v.version?c.push(!1):c.push(!0)}pt(h,e,{useGroups:t,forceUpdate:m,skipAssigningAttributes:c,overwriteIndex:r}),m&&e.dispose(),o.set(e,h.map(n=>({version:n.version,uuid:n.uuid})));let f=ae;return m?f=We:c.includes(!1)&&(f=He),{changeType:f,materials:mt(a),geometry:e}}}function vt(s){const e=new Set;for(let t=0,r=s.length;t<r;t++){const i=s[t];for(const o in i){const a=i[o];a&&a.isTexture&&e.add(a)}}return Array.from(e)}function xt(s){const e=[],t=new Set;for(let i=0,o=s.length;i<o;i++)s[i].traverse(a=>{a.visible&&(a.isRectAreaLight||a.isSpotLight||a.isPointLight||a.isDirectionalLight)&&(e.push(a),a.iesMap&&t.add(a.iesMap))});const r=Array.from(t).sort((i,o)=>i.uuid<o.uuid?1:i.uuid>o.uuid?-1:0);return{lights:e,iesTextures:r}}class yt{get initialized(){return!!this.bvh}constructor(e){this.bvhOptions={},this.attributes=["position","normal","tangent","color","uv","uv2"],this.generateBVH=!0,this.bvh=null,this.geometry=new l.BufferGeometry,this.staticGeometryGenerator=new gt(e),this._bvhWorker=null,this._pendingGenerate=null,this._buildAsync=!1}setObjects(e){this.staticGeometryGenerator.setObjects(e)}setBVHWorker(e){this._bvhWorker=e}async generateAsync(e=null){if(!this._bvhWorker)throw new Error('PathTracingSceneGenerator: "setBVHWorker" must be called before "generateAsync" can be called.');if(this.bvh instanceof Promise)return this._pendingGenerate||(this._pendingGenerate=new Promise(async()=>(await this.bvh,this._pendingGenerate=null,this.generateAsync(e)))),this._pendingGenerate;{this._buildAsync=!0;const t=this.generate(e);return this._buildAsync=!1,t.bvh=this.bvh=await t.bvh,t}}generate(e=null){const{staticGeometryGenerator:t,geometry:r,attributes:i}=this,o=t.objects;t.attributes=i,o.forEach(f=>{f.traverse(n=>{n.isSkinnedMesh&&n.skeleton&&n.skeleton.update()})});const a=t.generate(r),c=a.materials,h=vt(c),{lights:d,iesTextures:m}=xt(o);if(a.changeType!==ae&&ot(r,c,c),this.generateBVH){if(this.bvh instanceof Promise)throw new Error("PathTracingSceneGenerator: BVH is already building asynchronously.");if(a.changeType===We){const f={strategy:T.SAH,maxLeafTris:1,indirect:!0,onProgress:e,...this.bvhOptions};this._buildAsync?this.bvh=this._bvhWorker.generate(r,f):this.bvh=new T.MeshBVH(r,f)}else a.changeType===He&&this.bvh.refit()}return{bvhChanged:a.changeType!==ae,bvh:this.bvh,lights:d,iesTextures:m,geometry:r,materials:c,textures:h,objects:o}}}const bt=new l.OrthographicCamera(-1,1,1,-1,0,1);class Tt extends l.BufferGeometry{constructor(){super(),this.setAttribute("position",new l.Float32BufferAttribute([-1,3,0,-1,-1,0,3,-1,0],3)),this.setAttribute("uv",new l.Float32BufferAttribute([0,2,0,0,2,0],2))}}const wt=new Tt;class z{constructor(e){this._mesh=new l.Mesh(wt,e)}dispose(){this._mesh.geometry.dispose()}render(e){e.render(this._mesh,bt)}get material(){return this._mesh.material}set material(e){this._mesh.material=e}}class oe extends l.ShaderMaterial{set needsUpdate(e){super.needsUpdate=!0,this.dispatchEvent({type:"recompilation"})}constructor(e){super(e);for(const t in this.uniforms)Object.defineProperty(this,t,{get(){return this.uniforms[t].value},set(r){this.uniforms[t].value=r}})}setDefine(e,t=void 0){if(t==null){if(e in this.defines)return delete this.defines[e],this.needsUpdate=!0,!0}else if(this.defines[e]!==t)return this.defines[e]=t,this.needsUpdate=!0,!0;return!1}}class Rt extends oe{constructor(e){super({blending:l.NoBlending,uniforms:{target1:{value:null},target2:{value:null},opacity:{value:1}},vertexShader:`
+
+				varying vec2 vUv;
+
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}`,fragmentShader:`
+
+				uniform float opacity;
+
+				uniform sampler2D target1;
+				uniform sampler2D target2;
+
+				varying vec2 vUv;
+
+				void main() {
+
+					vec4 color1 = texture2D( target1, vUv );
+					vec4 color2 = texture2D( target2, vUv );
+
+					float invOpacity = 1.0 - opacity;
+					float totalAlpha = color1.a * invOpacity + color2.a * opacity;
+
+					if ( color1.a != 0.0 || color2.a != 0.0 ) {
+
+						gl_FragColor.rgb = color1.rgb * ( invOpacity * color1.a / totalAlpha ) + color2.rgb * ( opacity * color2.a / totalAlpha );
+						gl_FragColor.a = totalAlpha;
+
+					} else {
+
+						gl_FragColor = vec4( 0.0 );
+
+					}
+
+				}`}),this.setValues(e)}}function V(s=1){let e="uint";return s>1&&(e="uvec"+s),`
+		${e} sobolReverseBits( ${e} x ) {
+
+			x = ( ( ( x & 0xaaaaaaaau ) >> 1 ) | ( ( x & 0x55555555u ) << 1 ) );
+			x = ( ( ( x & 0xccccccccu ) >> 2 ) | ( ( x & 0x33333333u ) << 2 ) );
+			x = ( ( ( x & 0xf0f0f0f0u ) >> 4 ) | ( ( x & 0x0f0f0f0fu ) << 4 ) );
+			x = ( ( ( x & 0xff00ff00u ) >> 8 ) | ( ( x & 0x00ff00ffu ) << 8 ) );
+			return ( ( x >> 16 ) | ( x << 16 ) );
+
+		}
+
+		${e} sobolHashCombine( uint seed, ${e} v ) {
+
+			return seed ^ ( v + ${e}( ( seed << 6 ) + ( seed >> 2 ) ) );
+
+		}
+
+		${e} sobolLaineKarrasPermutation( ${e} x, ${e} seed ) {
+
+			x += seed;
+			x ^= x * 0x6c50b47cu;
+			x ^= x * 0xb82f1e52u;
+			x ^= x * 0xc7afe638u;
+			x ^= x * 0x8d22f6e6u;
+			return x;
+
+		}
+
+		${e} nestedUniformScrambleBase2( ${e} x, ${e} seed ) {
+
+			x = sobolLaineKarrasPermutation( x, seed );
+			x = sobolReverseBits( x );
+			return x;
+
+		}
+	`}function q(s=1){let e="uint",t="float",r="",i=".r",o="1u";return s>1&&(e="uvec"+s,t="vec"+s,r=s+"",s===2?(i=".rg",o="uvec2( 1u, 2u )"):s===3?(i=".rgb",o="uvec3( 1u, 2u, 3u )"):(i="",o="uvec4( 1u, 2u, 3u, 4u )")),`
+
+		${t} sobol${r}( int effect ) {
+
+			uint seed = sobolGetSeed( sobolBounceIndex, uint( effect ) );
+			uint index = sobolPathIndex;
+
+			uint shuffle_seed = sobolHashCombine( seed, 0u );
+			uint shuffled_index = nestedUniformScrambleBase2( sobolReverseBits( index ), shuffle_seed );
+			${t} sobol_pt = sobolGetTexturePoint( shuffled_index )${i};
+			${e} result = ${e}( sobol_pt * 16777216.0 );
+
+			${e} seed2 = sobolHashCombine( seed, ${o} );
+			result = nestedUniformScrambleBase2( result, seed2 );
+
+			return SOBOL_FACTOR * ${t}( result >> 8 );
+
+		}
+	`}const Le=`
+
+	// Utils
+	const float SOBOL_FACTOR = 1.0 / 16777216.0;
+	const uint SOBOL_MAX_POINTS = 256u * 256u;
+
+	${V(1)}
+	${V(2)}
+	${V(3)}
+	${V(4)}
+
+	uint sobolHash( uint x ) {
+
+		// finalizer from murmurhash3
+		x ^= x >> 16;
+		x *= 0x85ebca6bu;
+		x ^= x >> 13;
+		x *= 0xc2b2ae35u;
+		x ^= x >> 16;
+		return x;
+
+	}
+
+`,It=`
+
+	const uint SOBOL_DIRECTIONS_1[ 32 ] = uint[ 32 ](
+		0x80000000u, 0xc0000000u, 0xa0000000u, 0xf0000000u,
+		0x88000000u, 0xcc000000u, 0xaa000000u, 0xff000000u,
+		0x80800000u, 0xc0c00000u, 0xa0a00000u, 0xf0f00000u,
+		0x88880000u, 0xcccc0000u, 0xaaaa0000u, 0xffff0000u,
+		0x80008000u, 0xc000c000u, 0xa000a000u, 0xf000f000u,
+		0x88008800u, 0xcc00cc00u, 0xaa00aa00u, 0xff00ff00u,
+		0x80808080u, 0xc0c0c0c0u, 0xa0a0a0a0u, 0xf0f0f0f0u,
+		0x88888888u, 0xccccccccu, 0xaaaaaaaau, 0xffffffffu
+	);
+
+	const uint SOBOL_DIRECTIONS_2[ 32 ] = uint[ 32 ](
+		0x80000000u, 0xc0000000u, 0x60000000u, 0x90000000u,
+		0xe8000000u, 0x5c000000u, 0x8e000000u, 0xc5000000u,
+		0x68800000u, 0x9cc00000u, 0xee600000u, 0x55900000u,
+		0x80680000u, 0xc09c0000u, 0x60ee0000u, 0x90550000u,
+		0xe8808000u, 0x5cc0c000u, 0x8e606000u, 0xc5909000u,
+		0x6868e800u, 0x9c9c5c00u, 0xeeee8e00u, 0x5555c500u,
+		0x8000e880u, 0xc0005cc0u, 0x60008e60u, 0x9000c590u,
+		0xe8006868u, 0x5c009c9cu, 0x8e00eeeeu, 0xc5005555u
+	);
+
+	const uint SOBOL_DIRECTIONS_3[ 32 ] = uint[ 32 ](
+		0x80000000u, 0xc0000000u, 0x20000000u, 0x50000000u,
+		0xf8000000u, 0x74000000u, 0xa2000000u, 0x93000000u,
+		0xd8800000u, 0x25400000u, 0x59e00000u, 0xe6d00000u,
+		0x78080000u, 0xb40c0000u, 0x82020000u, 0xc3050000u,
+		0x208f8000u, 0x51474000u, 0xfbea2000u, 0x75d93000u,
+		0xa0858800u, 0x914e5400u, 0xdbe79e00u, 0x25db6d00u,
+		0x58800080u, 0xe54000c0u, 0x79e00020u, 0xb6d00050u,
+		0x800800f8u, 0xc00c0074u, 0x200200a2u, 0x50050093u
+	);
+
+	const uint SOBOL_DIRECTIONS_4[ 32 ] = uint[ 32 ](
+		0x80000000u, 0x40000000u, 0x20000000u, 0xb0000000u,
+		0xf8000000u, 0xdc000000u, 0x7a000000u, 0x9d000000u,
+		0x5a800000u, 0x2fc00000u, 0xa1600000u, 0xf0b00000u,
+		0xda880000u, 0x6fc40000u, 0x81620000u, 0x40bb0000u,
+		0x22878000u, 0xb3c9c000u, 0xfb65a000u, 0xddb2d000u,
+		0x78022800u, 0x9c0b3c00u, 0x5a0fb600u, 0x2d0ddb00u,
+		0xa2878080u, 0xf3c9c040u, 0xdb65a020u, 0x6db2d0b0u,
+		0x800228f8u, 0x400b3cdcu, 0x200fb67au, 0xb00ddb9du
+	);
+
+	uint getMaskedSobol( uint index, uint directions[ 32 ] ) {
+
+		uint X = 0u;
+		for ( int bit = 0; bit < 32; bit ++ ) {
+
+			uint mask = ( index >> bit ) & 1u;
+			X ^= mask * directions[ bit ];
+
+		}
+		return X;
+
+	}
+
+	vec4 generateSobolPoint( uint index ) {
+
+		if ( index >= SOBOL_MAX_POINTS ) {
+
+			return vec4( 0.0 );
+
+		}
+
+		// NOTE: this sobol "direction" is also available but we can't write out 5 components
+		// uint x = index & 0x00ffffffu;
+		uint x = sobolReverseBits( getMaskedSobol( index, SOBOL_DIRECTIONS_1 ) ) & 0x00ffffffu;
+		uint y = sobolReverseBits( getMaskedSobol( index, SOBOL_DIRECTIONS_2 ) ) & 0x00ffffffu;
+		uint z = sobolReverseBits( getMaskedSobol( index, SOBOL_DIRECTIONS_3 ) ) & 0x00ffffffu;
+		uint w = sobolReverseBits( getMaskedSobol( index, SOBOL_DIRECTIONS_4 ) ) & 0x00ffffffu;
+
+		return vec4( x, y, z, w ) * SOBOL_FACTOR;
+
+	}
+
+`,St=`
+
+	// Seeds
+	uniform sampler2D sobolTexture;
+	uint sobolPixelIndex = 0u;
+	uint sobolPathIndex = 0u;
+	uint sobolBounceIndex = 0u;
+
+	uint sobolGetSeed( uint bounce, uint effect ) {
+
+		return sobolHash(
+			sobolHashCombine(
+				sobolHashCombine(
+					sobolHash( bounce ),
+					sobolPixelIndex
+				),
+				effect
+			)
+		);
+
+	}
+
+	vec4 sobolGetTexturePoint( uint index ) {
+
+		if ( index >= SOBOL_MAX_POINTS ) {
+
+			index = index % SOBOL_MAX_POINTS;
+
+		}
+
+		uvec2 dim = uvec2( textureSize( sobolTexture, 0 ).xy );
+		uint y = index / dim.x;
+		uint x = index - y * dim.x;
+		vec2 uv = vec2( x, y ) / vec2( dim );
+		return texture( sobolTexture, uv );
+
+	}
+
+	${q(1)}
+	${q(2)}
+	${q(3)}
+	${q(4)}
+
+`;class _t extends oe{constructor(){super({blending:l.NoBlending,uniforms:{resolution:{value:new l.Vector2}},vertexShader:`
+
+				varying vec2 vUv;
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}
+			`,fragmentShader:`
+
+				${Le}
+				${It}
+
+				varying vec2 vUv;
+				uniform vec2 resolution;
+				void main() {
+
+					uint index = uint( gl_FragCoord.y ) * uint( resolution.x ) + uint( gl_FragCoord.x );
+					gl_FragColor = generateSobolPoint( index );
+
+				}
+			`})}}class Ft{generate(e,t=256){const r=new l.WebGLRenderTarget(t,t,{type:l.FloatType,format:l.RGBAFormat,minFilter:l.NearestFilter,magFilter:l.NearestFilter,generateMipmaps:!1}),i=e.getRenderTarget();e.setRenderTarget(r);const o=new z(new _t);return o.material.resolution.set(t,t),o.render(e),e.setRenderTarget(i),o.dispose(),r}}class At extends l.PerspectiveCamera{set bokehSize(e){this.fStop=this.getFocalLength()/e}get bokehSize(){return this.getFocalLength()/this.fStop}constructor(...e){super(...e),this.fStop=1.4,this.apertureBlades=0,this.apertureRotation=0,this.focusDistance=25,this.anamorphicRatio=1}copy(e,t){return super.copy(e,t),this.fStop=e.fStop,this.apertureBlades=e.apertureBlades,this.apertureRotation=e.apertureRotation,this.focusDistance=e.focusDistance,this.anamorphicRatio=e.anamorphicRatio,this}}class Mt{constructor(){this.bokehSize=0,this.apertureBlades=0,this.apertureRotation=0,this.focusDistance=10,this.anamorphicRatio=1}updateFrom(e){e instanceof At?(this.bokehSize=e.bokehSize,this.apertureBlades=e.apertureBlades,this.apertureRotation=e.apertureRotation,this.focusDistance=e.focusDistance,this.anamorphicRatio=e.anamorphicRatio):(this.bokehSize=0,this.apertureRotation=0,this.apertureBlades=0,this.focusDistance=10,this.anamorphicRatio=1)}}function te(s){const e=new Uint16Array(s.length);for(let t=0,r=s.length;t<r;++t)e[t]=l.DataUtils.toHalfFloat(s[t]);return e}function Te(s,e,t=0,r=s.length){let i=t,o=t+r-1;for(;i<o;){const a=i+o>>1;s[a]<e?i=a+1:o=a}return i-t}function Ct(s,e,t){return .2126*s+.7152*e+.0722*t}function Dt(s,e=l.HalfFloatType){const t=s.clone();t.source=new l.Source({...t.image});const{width:r,height:i,data:o}=t.image;let a=o;if(t.type!==e){e===l.HalfFloatType?a=new Uint16Array(o.length):a=new Float32Array(o.length);let c;o instanceof Int8Array||o instanceof Int16Array||o instanceof Int32Array?c=2**(8*o.BYTES_PER_ELEMENT-1)-1:c=2**(8*o.BYTES_PER_ELEMENT)-1;for(let h=0,d=o.length;h<d;h++){let m=o[h];t.type===l.HalfFloatType&&(m=l.DataUtils.fromHalfFloat(o[h])),t.type!==l.FloatType&&t.type!==l.HalfFloatType&&(m/=c),e===l.HalfFloatType&&(a[h]=l.DataUtils.toHalfFloat(m))}t.image.data=a,t.type=e}if(t.flipY){const c=a;a=a.slice();for(let h=0;h<i;h++)for(let d=0;d<r;d++){const m=i-h-1,f=4*(h*r+d),n=4*(m*r+d);a[n+0]=c[f+0],a[n+1]=c[f+1],a[n+2]=c[f+2],a[n+3]=c[f+3]}t.flipY=!1,t.image.data=a}return t}class Pt{constructor(){const e=new l.DataTexture(te(new Float32Array([0,0,0,0])),1,1);e.type=l.HalfFloatType,e.format=l.RGBAFormat,e.minFilter=l.LinearFilter,e.magFilter=l.LinearFilter,e.wrapS=l.RepeatWrapping,e.wrapT=l.RepeatWrapping,e.generateMipmaps=!1,e.needsUpdate=!0;const t=new l.DataTexture(te(new Float32Array([0,1])),1,2);t.type=l.HalfFloatType,t.format=l.RedFormat,t.minFilter=l.LinearFilter,t.magFilter=l.LinearFilter,t.generateMipmaps=!1,t.needsUpdate=!0;const r=new l.DataTexture(te(new Float32Array([0,0,1,1])),2,2);r.type=l.HalfFloatType,r.format=l.RedFormat,r.minFilter=l.LinearFilter,r.magFilter=l.LinearFilter,r.generateMipmaps=!1,r.needsUpdate=!0,this.map=e,this.marginalWeights=t,this.conditionalWeights=r,this.totalSum=0}dispose(){this.marginalWeights.dispose(),this.conditionalWeights.dispose(),this.map.dispose()}updateFrom(e){const t=Dt(e);t.wrapS=l.RepeatWrapping,t.wrapT=l.ClampToEdgeWrapping;const{width:r,height:i,data:o}=t.image,a=new Float32Array(r*i),c=new Float32Array(r*i),h=new Float32Array(i),d=new Float32Array(i);let m=0,f=0;for(let u=0;u<i;u++){let g=0;for(let y=0;y<r;y++){const b=u*r+y,w=l.DataUtils.fromHalfFloat(o[4*b+0]),R=l.DataUtils.fromHalfFloat(o[4*b+1]),_=l.DataUtils.fromHalfFloat(o[4*b+2]),F=Ct(w,R,_);g+=F,m+=F,a[b]=F,c[b]=g}if(g!==0)for(let y=u*r,b=u*r+r;y<b;y++)a[y]/=g,c[y]/=g;f+=g,h[u]=g,d[u]=f}if(f!==0)for(let u=0,g=h.length;u<g;u++)h[u]/=f,d[u]/=f;const n=new Uint16Array(i),x=new Uint16Array(r*i);for(let u=0;u<i;u++){const g=(u+1)/i,y=Te(d,g);n[u]=l.DataUtils.toHalfFloat((y+.5)/i)}for(let u=0;u<i;u++)for(let g=0;g<r;g++){const y=u*r+g,b=(g+1)/r,w=Te(c,b,u*r,r);x[y]=l.DataUtils.toHalfFloat((w+.5)/r)}this.dispose();const{marginalWeights:p,conditionalWeights:v}=this;p.image={width:i,height:1,data:n},p.needsUpdate=!0,v.image={width:r,height:i,data:x},v.needsUpdate=!0,this.totalSum=m,this.map=t}}const re=6,Et=0,kt=1,Ot=2,Bt=3,Nt=4,S=new l.Vector3,I=new l.Vector3,we=new l.Matrix4,E=new l.Quaternion,Re=new l.Vector3,k=new l.Vector3,zt=new l.Vector3(0,1,0);class Ht{constructor(){const e=new l.DataTexture(new Float32Array(4),1,1);e.format=l.RGBAFormat,e.type=l.FloatType,e.wrapS=l.ClampToEdgeWrapping,e.wrapT=l.ClampToEdgeWrapping,e.generateMipmaps=!1,e.minFilter=l.NearestFilter,e.magFilter=l.NearestFilter,this.tex=e,this.count=0}updateFrom(e,t=[]){const r=this.tex,i=Math.max(e.length*re,1),o=Math.ceil(Math.sqrt(i));r.image.width!==o&&(r.dispose(),r.image.data=new Float32Array(o*o*4),r.image.width=o,r.image.height=o);const a=r.image.data;for(let h=0,d=e.length;h<d;h++){const m=e[h],f=h*re*4;let n=0;for(let p=0;p<re*4;p++)a[f+p]=0;m.getWorldPosition(I),a[f+n++]=I.x,a[f+n++]=I.y,a[f+n++]=I.z;let x=Et;if(m.isRectAreaLight&&m.isCircular?x=kt:m.isSpotLight?x=Ot:m.isDirectionalLight?x=Bt:m.isPointLight&&(x=Nt),a[f+n++]=x,a[f+n++]=m.color.r,a[f+n++]=m.color.g,a[f+n++]=m.color.b,a[f+n++]=m.intensity,m.getWorldQuaternion(E),m.isRectAreaLight)S.set(m.width,0,0).applyQuaternion(E),a[f+n++]=S.x,a[f+n++]=S.y,a[f+n++]=S.z,n++,I.set(0,m.height,0).applyQuaternion(E),a[f+n++]=I.x,a[f+n++]=I.y,a[f+n++]=I.z,a[f+n++]=S.cross(I).length()*(m.isCircular?Math.PI/4:1);else if(m.isSpotLight){const p=m.radius||0;Re.setFromMatrixPosition(m.matrixWorld),k.setFromMatrixPosition(m.target.matrixWorld),we.lookAt(Re,k,zt),E.setFromRotationMatrix(we),S.set(1,0,0).applyQuaternion(E),a[f+n++]=S.x,a[f+n++]=S.y,a[f+n++]=S.z,n++,I.set(0,1,0).applyQuaternion(E),a[f+n++]=I.x,a[f+n++]=I.y,a[f+n++]=I.z,a[f+n++]=Math.PI*p*p,a[f+n++]=p,a[f+n++]=m.decay,a[f+n++]=m.distance,a[f+n++]=Math.cos(m.angle),a[f+n++]=Math.cos(m.angle*(1-m.penumbra)),a[f+n++]=m.iesMap?t.indexOf(m.iesMap):-1}else if(m.isPointLight){const p=S.setFromMatrixPosition(m.matrixWorld);a[f+n++]=p.x,a[f+n++]=p.y,a[f+n++]=p.z,n++,n+=4,n+=1,a[f+n++]=m.decay,a[f+n++]=m.distance}else if(m.isDirectionalLight){const p=S.setFromMatrixPosition(m.matrixWorld),v=I.setFromMatrixPosition(m.target.matrixWorld);k.subVectors(p,v).normalize(),a[f+n++]=k.x,a[f+n++]=k.y,a[f+n++]=k.z}}this.count=e.length;const c=se(a.buffer);return this.hash!==c?(this.hash=c,r.needsUpdate=!0,!0):!1}}function Ie(s,e,t,r,i){if(e>r)throw new Error;const o=s.length/e,a=s.constructor.BYTES_PER_ELEMENT*8;let c=1;switch(s.constructor){case Uint8Array:case Uint16Array:case Uint32Array:c=2**a-1;break;case Int8Array:case Int16Array:case Int32Array:c=2**(a-1)-1;break}for(let h=0;h<o;h++){const d=4*h,m=e*h;for(let f=0;f<r;f++)t[i+d+f]=e>=f+1?s[m+f]/c:0}}class Wt extends l.DataArrayTexture{constructor(){super(),this._textures=[],this.type=l.FloatType,this.format=l.RGBAFormat,this.internalFormat="RGBA32F"}updateAttribute(e,t){const r=this._textures[e];r.updateFrom(t);const i=r.image,o=this.image;if(i.width!==o.width||i.height!==o.height)throw new Error("FloatAttributeTextureArray: Attribute must be the same dimensions when updating single layer.");const{width:a,height:c,data:h}=o,m=a*c*4*e;let f=t.itemSize;f===3&&(f=4),Ie(r.image.data,f,h,4,m),this.dispose(),this.needsUpdate=!0}setAttributes(e){const t=e[0].count,r=e.length;for(let f=0,n=r;f<n;f++)if(e[f].count!==t)throw new Error("FloatAttributeTextureArray: All attributes must have the same item count.");const i=this._textures;for(;i.length<r;){const f=new Ne;i.push(f)}for(;i.length>r;)i.pop();for(let f=0,n=r;f<n;f++)i[f].updateFrom(e[f]);const a=i[0].image,c=this.image;(a.width!==c.width||a.height!==c.height||a.depth!==r)&&(c.width=a.width,c.height=a.height,c.depth=r,c.data=new Float32Array(c.width*c.height*c.depth*4));const{data:h,width:d,height:m}=c;for(let f=0,n=r;f<n;f++){const x=i[f],v=d*m*4*f;let u=e[f].itemSize;u===3&&(u=4),Ie(x.image.data,u,h,4,v)}this.dispose(),this.needsUpdate=!0}}class Lt extends Wt{updateNormalAttribute(e){this.updateAttribute(0,e)}updateTangentAttribute(e){this.updateAttribute(1,e)}updateUvAttribute(e){this.updateAttribute(2,e)}updateColorAttribute(e){this.updateAttribute(3,e)}updateFrom(e,t,r,i){this.setAttributes([e,t,r,i])}}function ne(s,e){return s.uuid<e.uuid?1:s.uuid>e.uuid?-1:0}function ie(s){return`${s.source.uuid}:${s.colorSpace}`}function Ut(s){const e=new Set,t=[];for(let r=0,i=s.length;r<i;r++){const o=s[r],a=ie(o);e.has(a)||(e.add(a),t.push(o))}return t}function Gt(s){const e=s.map(r=>r.iesMap||null).filter(r=>r),t=new Set(e);return Array.from(t).sort(ne)}function Vt(s){const e=new Set;for(let r=0,i=s.length;r<i;r++){const o=s[r];for(const a in o){const c=o[a];c&&c.isTexture&&e.add(c)}}const t=Array.from(e);return Ut(t).sort(ne)}function qt(s){const e=[];return s.traverse(t=>{t.visible&&(t.isRectAreaLight||t.isSpotLight||t.isPointLight||t.isDirectionalLight)&&e.push(t)}),e.sort(ne)}const Ue=45,Se=Ue*4;class $t{constructor(){this._features={}}isUsed(e){return e in this._features}setUsed(e,t=!0){t===!1?delete this._features[e]:this._features[e]=!0}reset(){this._features={}}}class jt extends l.DataTexture{constructor(){super(new Float32Array(4),1,1),this.format=l.RGBAFormat,this.type=l.FloatType,this.wrapS=l.ClampToEdgeWrapping,this.wrapT=l.ClampToEdgeWrapping,this.minFilter=l.NearestFilter,this.magFilter=l.NearestFilter,this.generateMipmaps=!1,this.features=new $t}updateFrom(e,t){function r(p,v,u=-1){if(v in p&&p[v]){const g=ie(p[v]);return f[g]}else return u}function i(p,v,u){return v in p?p[v]:u}function o(p,v,u,g){const y=p[v]&&p[v].isTexture?p[v]:null;if(y){y.matrixAutoUpdate&&y.updateMatrix();const b=y.matrix.elements;let w=0;u[g+w++]=b[0],u[g+w++]=b[3],u[g+w++]=b[6],w++,u[g+w++]=b[1],u[g+w++]=b[4],u[g+w++]=b[7],w++}return 8}let a=0;const c=e.length*Ue,h=Math.ceil(Math.sqrt(c))||1,{image:d,features:m}=this,f={};for(let p=0,v=t.length;p<v;p++)f[ie(t[p])]=p;d.width!==h&&(this.dispose(),d.data=new Float32Array(h*h*4),d.width=h,d.height=h);const n=d.data;m.reset();for(let p=0,v=e.length;p<v;p++){const u=e[p];if(u.isFogVolumeMaterial){m.setUsed("FOG");for(let b=0;b<Se;b++)n[a+b]=0;n[a+0+0]=u.color.r,n[a+0+1]=u.color.g,n[a+0+2]=u.color.b,n[a+8+3]=i(u,"emissiveIntensity",0),n[a+12+0]=u.emissive.r,n[a+12+1]=u.emissive.g,n[a+12+2]=u.emissive.b,n[a+52+1]=u.density,n[a+52+3]=0,n[a+56+2]=4,a+=Se;continue}n[a++]=u.color.r,n[a++]=u.color.g,n[a++]=u.color.b,n[a++]=r(u,"map"),n[a++]=i(u,"metalness",0),n[a++]=r(u,"metalnessMap"),n[a++]=i(u,"roughness",0),n[a++]=r(u,"roughnessMap"),n[a++]=i(u,"ior",1.5),n[a++]=i(u,"transmission",0),n[a++]=r(u,"transmissionMap"),n[a++]=i(u,"emissiveIntensity",0),"emissive"in u?(n[a++]=u.emissive.r,n[a++]=u.emissive.g,n[a++]=u.emissive.b):(n[a++]=0,n[a++]=0,n[a++]=0),n[a++]=r(u,"emissiveMap"),n[a++]=r(u,"normalMap"),"normalScale"in u?(n[a++]=u.normalScale.x,n[a++]=u.normalScale.y):(n[a++]=1,n[a++]=1),n[a++]=i(u,"clearcoat",0),n[a++]=r(u,"clearcoatMap"),n[a++]=i(u,"clearcoatRoughness",0),n[a++]=r(u,"clearcoatRoughnessMap"),n[a++]=r(u,"clearcoatNormalMap"),"clearcoatNormalScale"in u?(n[a++]=u.clearcoatNormalScale.x,n[a++]=u.clearcoatNormalScale.y):(n[a++]=1,n[a++]=1),a++,n[a++]=i(u,"sheen",0),"sheenColor"in u?(n[a++]=u.sheenColor.r,n[a++]=u.sheenColor.g,n[a++]=u.sheenColor.b):(n[a++]=0,n[a++]=0,n[a++]=0),n[a++]=r(u,"sheenColorMap"),n[a++]=i(u,"sheenRoughness",0),n[a++]=r(u,"sheenRoughnessMap"),n[a++]=r(u,"iridescenceMap"),n[a++]=r(u,"iridescenceThicknessMap"),n[a++]=i(u,"iridescence",0),n[a++]=i(u,"iridescenceIOR",1.3);const g=i(u,"iridescenceThicknessRange",[100,400]);n[a++]=g[0],n[a++]=g[1],"specularColor"in u?(n[a++]=u.specularColor.r,n[a++]=u.specularColor.g,n[a++]=u.specularColor.b):(n[a++]=1,n[a++]=1,n[a++]=1),n[a++]=r(u,"specularColorMap"),n[a++]=i(u,"specularIntensity",1),n[a++]=r(u,"specularIntensityMap");const y=i(u,"thickness",0)===0&&i(u,"attenuationDistance",1/0)===1/0;if(n[a++]=Number(y),a++,"attenuationColor"in u?(n[a++]=u.attenuationColor.r,n[a++]=u.attenuationColor.g,n[a++]=u.attenuationColor.b):(n[a++]=1,n[a++]=1,n[a++]=1),n[a++]=i(u,"attenuationDistance",1/0),n[a++]=r(u,"alphaMap"),n[a++]=u.opacity,n[a++]=u.alphaTest,!y&&u.transmission>0)n[a++]=0;else switch(u.side){case l.FrontSide:n[a++]=1;break;case l.BackSide:n[a++]=-1;break;case l.DoubleSide:n[a++]=0;break}n[a++]=Number(i(u,"matte",!1)),n[a++]=Number(i(u,"castShadow",!0)),n[a++]=Number(u.vertexColors)|Number(u.flatShading)<<1,n[a++]=Number(u.transparent),a+=o(u,"map",n,a),a+=o(u,"metalnessMap",n,a),a+=o(u,"roughnessMap",n,a),a+=o(u,"transmissionMap",n,a),a+=o(u,"emissiveMap",n,a),a+=o(u,"normalMap",n,a),a+=o(u,"clearcoatMap",n,a),a+=o(u,"clearcoatNormalMap",n,a),a+=o(u,"clearcoatRoughnessMap",n,a),a+=o(u,"sheenColorMap",n,a),a+=o(u,"sheenRoughnessMap",n,a),a+=o(u,"iridescenceMap",n,a),a+=o(u,"iridescenceThicknessMap",n,a),a+=o(u,"specularColorMap",n,a),a+=o(u,"specularIntensityMap",n,a)}const x=se(n.buffer);return this.hash!==x?(this.hash=x,this.needsUpdate=!0,!0):!1}}const _e=new l.Color;function Yt(s){return s?`${s.uuid}:${s.version}`:null}function Xt(s,e){for(const t in e)t in s&&(s[t]=e[t])}class Fe extends l.WebGLArrayRenderTarget{constructor(e,t,r){const i={format:l.RGBAFormat,type:l.UnsignedByteType,minFilter:l.LinearFilter,magFilter:l.LinearFilter,wrapS:l.RepeatWrapping,wrapT:l.RepeatWrapping,generateMipmaps:!1,...r};super(e,t,1,i),Xt(this.texture,i),this.texture.setTextures=(...a)=>{this.setTextures(...a)},this.hashes=[null];const o=new z(new Qt);this.fsQuad=o}setTextures(e,t,r=this.width,i=this.height){const o=e.getRenderTarget(),a=e.toneMapping,c=e.getClearAlpha();e.getClearColor(_e);const h=t.length||1;(r!==this.width||i!==this.height||this.depth!==h)&&(this.setSize(r,i,h),this.hashes=new Array(h).fill(null)),e.setClearColor(0,0),e.toneMapping=l.NoToneMapping;const d=this.fsQuad,m=this.hashes;let f=!1;for(let n=0,x=h;n<x;n++){const p=t[n],v=Yt(p);p&&(m[n]!==v||p.isWebGLRenderTarget)&&(p.matrixAutoUpdate=!1,p.matrix.identity(),d.material.map=p,e.setRenderTarget(this,n),d.render(e),p.updateMatrix(),p.matrixAutoUpdate=!0,m[n]=v,f=!0)}return d.material.map=null,e.setClearColor(_e,c),e.setRenderTarget(o),e.toneMapping=a,f}dispose(){super.dispose(),this.fsQuad.dispose()}}class Qt extends l.ShaderMaterial{get map(){return this.uniforms.map.value}set map(e){this.uniforms.map.value=e}constructor(){super({uniforms:{map:{value:null}},vertexShader:`
+				varying vec2 vUv;
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}
+			`,fragmentShader:`
+				uniform sampler2D map;
+				varying vec2 vUv;
+				void main() {
+
+					gl_FragColor = texture2D( map, vUv );
+
+				}
+			`})}}function Kt(s,e=Math.random()){for(let t=s.length-1;t>0;t--){const r=Math.floor(e()*(t+1)),i=s[t];s[t]=s[r],s[r]=i}return s}class Zt{constructor(e,t,r=Math.random){const i=e**t,o=new Uint16Array(i);let a=i;for(let c=0;c<i;c++)o[c]=c;this.samples=new Float32Array(t),this.strataCount=e,this.reset=function(){for(let c=0;c<i;c++)o[c]=c;a=0},this.reshuffle=function(){a=0},this.next=function(){const{samples:c}=this;a>=o.length&&(Kt(o,r),this.reshuffle());let h=o[a++];for(let d=0;d<t;d++)c[d]=(h%e+r())/e,h=Math.floor(h/e);return c}}}class Jt{constructor(e,t,r=Math.random){let i=0;for(const h of t)i+=h;const o=new Float32Array(i),a=[];let c=0;for(const h of t){const d=new Zt(e,h,r);d.samples=new Float32Array(o.buffer,c,d.samples.length),c+=d.samples.length*4,a.push(d)}this.samples=o,this.strataCount=e,this.next=function(){for(const h of a)h.next();return o},this.reshuffle=function(){for(const h of a)h.reshuffle()},this.reset=function(){for(const h of a)h.reset()}}}class er{constructor(e=0){this.m=2147483648,this.a=1103515245,this.c=12345,this.seed=e}nextInt(){return this.seed=(this.a*this.seed+this.c)%this.m,this.seed}nextFloat(){return this.nextInt()/(this.m-1)}}class tr extends l.DataTexture{constructor(e=1,t=1,r=8){super(new Float32Array(1),1,1,l.RGBAFormat,l.FloatType),this.minFilter=l.NearestFilter,this.magFilter=l.NearestFilter,this.strata=r,this.sampler=null,this.generator=new er,this.stableNoise=!1,this.random=()=>this.stableNoise?this.generator.nextFloat():Math.random(),this.init(e,t,r)}init(e=this.image.height,t=this.image.width,r=this.strata){const{image:i}=this;if(i.width===t&&i.height===e&&this.sampler!==null)return;const o=new Array(e*t).fill(4),a=new Jt(r,o,this.random);i.width=t,i.height=e,i.data=a.samples,this.sampler=a,this.dispose(),this.next()}next(){this.sampler.next(),this.needsUpdate=!0}reset(){this.sampler.reset(),this.generator.seed=0}}function rr(s,e=Math.random){for(let t=s.length-1;t>0;t--){const r=~~((e()-1e-6)*t),i=s[t];s[t]=s[r],s[r]=i}}function ar(s,e){s.fill(0);for(let t=0;t<e;t++)s[t]=1}class Ae{constructor(e){this.count=0,this.size=-1,this.sigma=-1,this.radius=-1,this.lookupTable=null,this.score=null,this.binaryPattern=null,this.resize(e),this.setSigma(1.5)}findVoid(){const{score:e,binaryPattern:t}=this;let r=1/0,i=-1;for(let o=0,a=t.length;o<a;o++){if(t[o]!==0)continue;const c=e[o];c<r&&(r=c,i=o)}return i}findCluster(){const{score:e,binaryPattern:t}=this;let r=-1/0,i=-1;for(let o=0,a=t.length;o<a;o++){if(t[o]!==1)continue;const c=e[o];c>r&&(r=c,i=o)}return i}setSigma(e){if(e===this.sigma)return;const t=~~(Math.sqrt(20*e**2)+1),r=2*t+1,i=new Float32Array(r*r),o=e*e;for(let a=-t;a<=t;a++)for(let c=-t;c<=t;c++){const h=(t+c)*r+a+t,d=a*a+c*c;i[h]=Math.E**(-d/(2*o))}this.lookupTable=i,this.sigma=e,this.radius=t}resize(e){this.size!==e&&(this.size=e,this.score=new Float32Array(e*e),this.binaryPattern=new Uint8Array(e*e))}invert(){const{binaryPattern:e,score:t,size:r}=this;t.fill(0);for(let i=0,o=e.length;i<o;i++)if(e[i]===0){const a=~~(i/r),c=i-a*r;this.updateScore(c,a,1),e[i]=1}else e[i]=0}updateScore(e,t,r){const{size:i,score:o,lookupTable:a}=this,c=this.radius,h=2*c+1;for(let d=-c;d<=c;d++)for(let m=-c;m<=c;m++){const f=(c+m)*h+d+c,n=a[f];let x=e+d;x=x<0?i+x:x%i;let p=t+m;p=p<0?i+p:p%i;const v=p*i+x;o[v]+=r*n}}addPointIndex(e){this.binaryPattern[e]=1;const t=this.size,r=~~(e/t),i=e-r*t;this.updateScore(i,r,1),this.count++}removePointIndex(e){this.binaryPattern[e]=0;const t=this.size,r=~~(e/t),i=e-r*t;this.updateScore(i,r,-1),this.count--}copy(e){this.resize(e.size),this.score.set(e.score),this.binaryPattern.set(e.binaryPattern),this.setSigma(e.sigma),this.count=e.count}}class ir{constructor(){this.random=Math.random,this.sigma=1.5,this.size=64,this.majorityPointsRatio=.1,this.samples=new Ae(1),this.savedSamples=new Ae(1)}generate(){const{samples:e,savedSamples:t,sigma:r,majorityPointsRatio:i,size:o}=this;e.resize(o),e.setSigma(r);const a=Math.floor(o*o*i),c=e.binaryPattern;ar(c,a),rr(c,this.random);for(let f=0,n=c.length;f<n;f++)c[f]===1&&e.addPointIndex(f);for(;;){const f=e.findCluster();e.removePointIndex(f);const n=e.findVoid();if(f===n){e.addPointIndex(f);break}e.addPointIndex(n)}const h=new Uint32Array(o*o);t.copy(e);let d;for(d=e.count-1;d>=0;){const f=e.findCluster();e.removePointIndex(f),h[f]=d,d--}const m=o*o;for(d=t.count;d<m/2;){const f=t.findVoid();t.addPointIndex(f),h[f]=d,d++}for(t.invert();d<m;){const f=t.findCluster();t.removePointIndex(f),h[f]=d,d++}return{data:h,maxValue:m}}}function sr(s){return s>=3?4:s}function or(s){switch(s){case 1:return l.RedFormat;case 2:return l.RGFormat;default:return l.RGBAFormat}}class nr extends l.DataTexture{constructor(e=64,t=1){super(new Float32Array(4),1,1,l.RGBAFormat,l.FloatType),this.minFilter=l.NearestFilter,this.magFilter=l.NearestFilter,this.size=e,this.channels=t,this.update()}update(){const e=this.channels,t=this.size,r=new ir;r.channels=e,r.size=t;const i=sr(e),o=or(i);(this.image.width!==t||o!==this.format)&&(this.image.width=t,this.image.height=t,this.image.data=new Float32Array(t**2*i),this.format=o,this.dispose());const a=this.image.data;for(let c=0,h=e;c<h;c++){const d=r.generate(),m=d.data,f=d.maxValue;for(let n=0,x=m.length;n<x;n++){const p=m[n]/f;a[n*i+c]=p}}this.needsUpdate=!0}}const lr=`
+
+	struct PhysicalCamera {
+
+		float focusDistance;
+		float anamorphicRatio;
+		float bokehSize;
+		int apertureBlades;
+		float apertureRotation;
+
+	};
+
+`,cr=`
+
+	struct EquirectHdrInfo {
+
+		sampler2D marginalWeights;
+		sampler2D conditionalWeights;
+		sampler2D map;
+
+		float totalSum;
+
+	};
+
+`,ur=`
+
+	#define RECT_AREA_LIGHT_TYPE 0
+	#define CIRC_AREA_LIGHT_TYPE 1
+	#define SPOT_LIGHT_TYPE 2
+	#define DIR_LIGHT_TYPE 3
+	#define POINT_LIGHT_TYPE 4
+
+	struct LightsInfo {
+
+		sampler2D tex;
+		uint count;
+
+	};
+
+	struct Light {
+
+		vec3 position;
+		int type;
+
+		vec3 color;
+		float intensity;
+
+		vec3 u;
+		vec3 v;
+		float area;
+
+		// spot light fields
+		float radius;
+		float near;
+		float decay;
+		float distance;
+		float coneCos;
+		float penumbraCos;
+		int iesProfile;
+
+	};
+
+	Light readLightInfo( sampler2D tex, uint index ) {
+
+		uint i = index * 6u;
+
+		vec4 s0 = texelFetch1D( tex, i + 0u );
+		vec4 s1 = texelFetch1D( tex, i + 1u );
+		vec4 s2 = texelFetch1D( tex, i + 2u );
+		vec4 s3 = texelFetch1D( tex, i + 3u );
+
+		Light l;
+		l.position = s0.rgb;
+		l.type = int( round( s0.a ) );
+
+		l.color = s1.rgb;
+		l.intensity = s1.a;
+
+		l.u = s2.rgb;
+		l.v = s3.rgb;
+		l.area = s3.a;
+
+		if ( l.type == SPOT_LIGHT_TYPE || l.type == POINT_LIGHT_TYPE ) {
+
+			vec4 s4 = texelFetch1D( tex, i + 4u );
+			vec4 s5 = texelFetch1D( tex, i + 5u );
+			l.radius = s4.r;
+			l.decay = s4.g;
+			l.distance = s4.b;
+			l.coneCos = s4.a;
+
+			l.penumbraCos = s5.r;
+			l.iesProfile = int( round( s5.g ) );
+
+		} else {
+
+			l.radius = 0.0;
+			l.decay = 0.0;
+			l.distance = 0.0;
+
+			l.coneCos = 0.0;
+			l.penumbraCos = 0.0;
+			l.iesProfile = - 1;
+
+		}
+
+		return l;
+
+	}
+
+`,fr=`
+
+	struct Material {
+
+		vec3 color;
+		int map;
+
+		float metalness;
+		int metalnessMap;
+
+		float roughness;
+		int roughnessMap;
+
+		float ior;
+		float transmission;
+		int transmissionMap;
+
+		float emissiveIntensity;
+		vec3 emissive;
+		int emissiveMap;
+
+		int normalMap;
+		vec2 normalScale;
+
+		float clearcoat;
+		int clearcoatMap;
+		int clearcoatNormalMap;
+		vec2 clearcoatNormalScale;
+		float clearcoatRoughness;
+		int clearcoatRoughnessMap;
+
+		int iridescenceMap;
+		int iridescenceThicknessMap;
+		float iridescence;
+		float iridescenceIor;
+		float iridescenceThicknessMinimum;
+		float iridescenceThicknessMaximum;
+
+		vec3 specularColor;
+		int specularColorMap;
+
+		float specularIntensity;
+		int specularIntensityMap;
+		bool thinFilm;
+
+		vec3 attenuationColor;
+		float attenuationDistance;
+
+		int alphaMap;
+
+		bool castShadow;
+		float opacity;
+		float alphaTest;
+
+		float side;
+		bool matte;
+
+		float sheen;
+		vec3 sheenColor;
+		int sheenColorMap;
+		float sheenRoughness;
+		int sheenRoughnessMap;
+
+		bool vertexColors;
+		bool flatShading;
+		bool transparent;
+		bool fogVolume;
+
+		mat3 mapTransform;
+		mat3 metalnessMapTransform;
+		mat3 roughnessMapTransform;
+		mat3 transmissionMapTransform;
+		mat3 emissiveMapTransform;
+		mat3 normalMapTransform;
+		mat3 clearcoatMapTransform;
+		mat3 clearcoatNormalMapTransform;
+		mat3 clearcoatRoughnessMapTransform;
+		mat3 sheenColorMapTransform;
+		mat3 sheenRoughnessMapTransform;
+		mat3 iridescenceMapTransform;
+		mat3 iridescenceThicknessMapTransform;
+		mat3 specularColorMapTransform;
+		mat3 specularIntensityMapTransform;
+
+	};
+
+	mat3 readTextureTransform( sampler2D tex, uint index ) {
+
+		mat3 textureTransform;
+
+		vec4 row1 = texelFetch1D( tex, index );
+		vec4 row2 = texelFetch1D( tex, index + 1u );
+
+		textureTransform[0] = vec3(row1.r, row2.r, 0.0);
+		textureTransform[1] = vec3(row1.g, row2.g, 0.0);
+		textureTransform[2] = vec3(row1.b, row2.b, 1.0);
+
+		return textureTransform;
+
+	}
+
+	Material readMaterialInfo( sampler2D tex, uint index ) {
+
+		uint i = index * 45u;
+
+		vec4 s0 = texelFetch1D( tex, i + 0u );
+		vec4 s1 = texelFetch1D( tex, i + 1u );
+		vec4 s2 = texelFetch1D( tex, i + 2u );
+		vec4 s3 = texelFetch1D( tex, i + 3u );
+		vec4 s4 = texelFetch1D( tex, i + 4u );
+		vec4 s5 = texelFetch1D( tex, i + 5u );
+		vec4 s6 = texelFetch1D( tex, i + 6u );
+		vec4 s7 = texelFetch1D( tex, i + 7u );
+		vec4 s8 = texelFetch1D( tex, i + 8u );
+		vec4 s9 = texelFetch1D( tex, i + 9u );
+		vec4 s10 = texelFetch1D( tex, i + 10u );
+		vec4 s11 = texelFetch1D( tex, i + 11u );
+		vec4 s12 = texelFetch1D( tex, i + 12u );
+		vec4 s13 = texelFetch1D( tex, i + 13u );
+		vec4 s14 = texelFetch1D( tex, i + 14u );
+
+		Material m;
+		m.color = s0.rgb;
+		m.map = int( round( s0.a ) );
+
+		m.metalness = s1.r;
+		m.metalnessMap = int( round( s1.g ) );
+		m.roughness = s1.b;
+		m.roughnessMap = int( round( s1.a ) );
+
+		m.ior = s2.r;
+		m.transmission = s2.g;
+		m.transmissionMap = int( round( s2.b ) );
+		m.emissiveIntensity = s2.a;
+
+		m.emissive = s3.rgb;
+		m.emissiveMap = int( round( s3.a ) );
+
+		m.normalMap = int( round( s4.r ) );
+		m.normalScale = s4.gb;
+
+		m.clearcoat = s4.a;
+		m.clearcoatMap = int( round( s5.r ) );
+		m.clearcoatRoughness = s5.g;
+		m.clearcoatRoughnessMap = int( round( s5.b ) );
+		m.clearcoatNormalMap = int( round( s5.a ) );
+		m.clearcoatNormalScale = s6.rg;
+
+		m.sheen = s6.a;
+		m.sheenColor = s7.rgb;
+		m.sheenColorMap = int( round( s7.a ) );
+		m.sheenRoughness = s8.r;
+		m.sheenRoughnessMap = int( round( s8.g ) );
+
+		m.iridescenceMap = int( round( s8.b ) );
+		m.iridescenceThicknessMap = int( round( s8.a ) );
+		m.iridescence = s9.r;
+		m.iridescenceIor = s9.g;
+		m.iridescenceThicknessMinimum = s9.b;
+		m.iridescenceThicknessMaximum = s9.a;
+
+		m.specularColor = s10.rgb;
+		m.specularColorMap = int( round( s10.a ) );
+
+		m.specularIntensity = s11.r;
+		m.specularIntensityMap = int( round( s11.g ) );
+		m.thinFilm = bool( s11.b );
+
+		m.attenuationColor = s12.rgb;
+		m.attenuationDistance = s12.a;
+
+		m.alphaMap = int( round( s13.r ) );
+
+		m.opacity = s13.g;
+		m.alphaTest = s13.b;
+		m.side = s13.a;
+
+		m.matte = bool( s14.r );
+		m.castShadow = bool( s14.g );
+		m.vertexColors = bool( int( s14.b ) & 1 );
+		m.flatShading = bool( int( s14.b ) & 2 );
+		m.fogVolume = bool( int( s14.b ) & 4 );
+		m.transparent = bool( s14.a );
+
+		uint firstTextureTransformIdx = i + 15u;
+
+		// mat3( 1.0 ) is an identity matrix
+		m.mapTransform = m.map == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx );
+		m.metalnessMapTransform = m.metalnessMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 2u );
+		m.roughnessMapTransform = m.roughnessMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 4u );
+		m.transmissionMapTransform = m.transmissionMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 6u );
+		m.emissiveMapTransform = m.emissiveMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 8u );
+		m.normalMapTransform = m.normalMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 10u );
+		m.clearcoatMapTransform = m.clearcoatMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 12u );
+		m.clearcoatNormalMapTransform = m.clearcoatNormalMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 14u );
+		m.clearcoatRoughnessMapTransform = m.clearcoatRoughnessMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 16u );
+		m.sheenColorMapTransform = m.sheenColorMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 18u );
+		m.sheenRoughnessMapTransform = m.sheenRoughnessMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 20u );
+		m.iridescenceMapTransform = m.iridescenceMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 22u );
+		m.iridescenceThicknessMapTransform = m.iridescenceThicknessMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 24u );
+		m.specularColorMapTransform = m.specularColorMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 26u );
+		m.specularIntensityMapTransform = m.specularIntensityMap == - 1 ? mat3( 1.0 ) : readTextureTransform( tex, firstTextureTransformIdx + 28u );
+
+		return m;
+
+	}
+
+`,hr=`
+
+	struct SurfaceRecord {
+
+		// surface type
+		bool volumeParticle;
+
+		// geometry
+		vec3 faceNormal;
+		bool frontFace;
+		vec3 normal;
+		mat3 normalBasis;
+		mat3 normalInvBasis;
+
+		// cached properties
+		float eta;
+		float f0;
+
+		// material
+		float roughness;
+		float filteredRoughness;
+		float metalness;
+		vec3 color;
+		vec3 emission;
+
+		// transmission
+		float ior;
+		float transmission;
+		bool thinFilm;
+		vec3 attenuationColor;
+		float attenuationDistance;
+
+		// clearcoat
+		vec3 clearcoatNormal;
+		mat3 clearcoatBasis;
+		mat3 clearcoatInvBasis;
+		float clearcoat;
+		float clearcoatRoughness;
+		float filteredClearcoatRoughness;
+
+		// sheen
+		float sheen;
+		vec3 sheenColor;
+		float sheenRoughness;
+
+		// iridescence
+		float iridescence;
+		float iridescenceIor;
+		float iridescenceThickness;
+
+		// specular
+		vec3 specularColor;
+		float specularIntensity;
+	};
+
+	struct ScatterRecord {
+		float specularPdf;
+		float pdf;
+		vec3 direction;
+		vec3 color;
+	};
+
+`,dr=`
+
+	// samples the the given environment map in the given direction
+	vec3 sampleEquirectColor( sampler2D envMap, vec3 direction ) {
+
+		return texture2D( envMap, equirectDirectionToUv( direction ) ).rgb;
+
+	}
+
+	// gets the pdf of the given direction to sample
+	float equirectDirectionPdf( vec3 direction ) {
+
+		vec2 uv = equirectDirectionToUv( direction );
+		float theta = uv.y * PI;
+		float sinTheta = sin( theta );
+		if ( sinTheta == 0.0 ) {
+
+			return 0.0;
+
+		}
+
+		return 1.0 / ( 2.0 * PI * PI * sinTheta );
+
+	}
+
+	// samples the color given env map with CDF and returns the pdf of the direction
+	float sampleEquirect( vec3 direction, inout vec3 color ) {
+
+		float totalSum = envMapInfo.totalSum;
+		if ( totalSum == 0.0 ) {
+
+			color = vec3( 0.0 );
+			return 1.0;
+
+		}
+
+		vec2 uv = equirectDirectionToUv( direction );
+		color = texture2D( envMapInfo.map, uv ).rgb;
+
+		float lum = luminance( color );
+		ivec2 resolution = textureSize( envMapInfo.map, 0 );
+		float pdf = lum / totalSum;
+
+		return float( resolution.x * resolution.y ) * pdf * equirectDirectionPdf( direction );
+
+	}
+
+	// samples a direction of the envmap with color and retrieves pdf
+	float sampleEquirectProbability( vec2 r, inout vec3 color, inout vec3 direction ) {
+
+		// sample env map cdf
+		float v = texture2D( envMapInfo.marginalWeights, vec2( r.x, 0.0 ) ).x;
+		float u = texture2D( envMapInfo.conditionalWeights, vec2( r.y, v ) ).x;
+		vec2 uv = vec2( u, v );
+
+		vec3 derivedDirection = equirectUvToDirection( uv );
+		direction = derivedDirection;
+		color = texture2D( envMapInfo.map, uv ).rgb;
+
+		float totalSum = envMapInfo.totalSum;
+		float lum = luminance( color );
+		ivec2 resolution = textureSize( envMapInfo.map, 0 );
+		float pdf = lum / totalSum;
+
+		return float( resolution.x * resolution.y ) * pdf * equirectDirectionPdf( direction );
+
+	}
+`,mr=`
+
+	float getSpotAttenuation( const in float coneCosine, const in float penumbraCosine, const in float angleCosine ) {
+
+		return smoothstep( coneCosine, penumbraCosine, angleCosine );
+
+	}
+
+	float getDistanceAttenuation( const in float lightDistance, const in float cutoffDistance, const in float decayExponent ) {
+
+		// based upon Frostbite 3 Moving to Physically-based Rendering
+		// page 32, equation 26: E[window1]
+		// https://seblagarde.files.wordpress.com/2015/07/course_notes_moving_frostbite_to_pbr_v32.pdf
+		float distanceFalloff = 1.0 / max( pow( lightDistance, decayExponent ), EPSILON );
+
+		if ( cutoffDistance > 0.0 ) {
+
+			distanceFalloff *= pow2( saturate( 1.0 - pow4( lightDistance / cutoffDistance ) ) );
+
+		}
+
+		return distanceFalloff;
+
+	}
+
+	float getPhotometricAttenuation( sampler2DArray iesProfiles, int iesProfile, vec3 posToLight, vec3 lightDir, vec3 u, vec3 v ) {
+
+		float cosTheta = dot( posToLight, lightDir );
+		float angle = acos( cosTheta ) / PI;
+
+		return texture2D( iesProfiles, vec3( angle, 0.0, iesProfile ) ).r;
+
+	}
+
+	struct LightRecord {
+
+		float dist;
+		vec3 direction;
+		float pdf;
+		vec3 emission;
+		int type;
+
+	};
+
+	bool intersectLightAtIndex( sampler2D lights, vec3 rayOrigin, vec3 rayDirection, uint l, inout LightRecord lightRec ) {
+
+		bool didHit = false;
+		Light light = readLightInfo( lights, l );
+
+		vec3 u = light.u;
+		vec3 v = light.v;
+
+		// check for backface
+		vec3 normal = normalize( cross( u, v ) );
+		if ( dot( normal, rayDirection ) > 0.0 ) {
+
+			u *= 1.0 / dot( u, u );
+			v *= 1.0 / dot( v, v );
+
+			float dist;
+
+			// MIS / light intersection is not supported for punctual lights.
+			if(
+				( light.type == RECT_AREA_LIGHT_TYPE && intersectsRectangle( light.position, normal, u, v, rayOrigin, rayDirection, dist ) ) ||
+				( light.type == CIRC_AREA_LIGHT_TYPE && intersectsCircle( light.position, normal, u, v, rayOrigin, rayDirection, dist ) )
+			) {
+
+				float cosTheta = dot( rayDirection, normal );
+				didHit = true;
+				lightRec.dist = dist;
+				lightRec.pdf = ( dist * dist ) / ( light.area * cosTheta );
+				lightRec.emission = light.color * light.intensity;
+				lightRec.direction = rayDirection;
+				lightRec.type = light.type;
+
+			}
+
+		}
+
+		return didHit;
+
+	}
+
+	LightRecord randomAreaLightSample( Light light, vec3 rayOrigin, vec2 ruv ) {
+
+		vec3 randomPos;
+		if( light.type == RECT_AREA_LIGHT_TYPE ) {
+
+			// rectangular area light
+			randomPos = light.position + light.u * ( ruv.x - 0.5 ) + light.v * ( ruv.y - 0.5 );
+
+		} else if( light.type == CIRC_AREA_LIGHT_TYPE ) {
+
+			// circular area light
+			float r = 0.5 * sqrt( ruv.x );
+			float theta = ruv.y * 2.0 * PI;
+			float x = r * cos( theta );
+			float y = r * sin( theta );
+
+			randomPos = light.position + light.u * x + light.v * y;
+
+		}
+
+		vec3 toLight = randomPos - rayOrigin;
+		float lightDistSq = dot( toLight, toLight );
+		float dist = sqrt( lightDistSq );
+		vec3 direction = toLight / dist;
+		vec3 lightNormal = normalize( cross( light.u, light.v ) );
+
+		LightRecord lightRec;
+		lightRec.type = light.type;
+		lightRec.emission = light.color * light.intensity;
+		lightRec.dist = dist;
+		lightRec.direction = direction;
+
+		// TODO: the denominator is potentially zero
+		lightRec.pdf = lightDistSq / ( light.area * dot( direction, lightNormal ) );
+
+		return lightRec;
+
+	}
+
+	LightRecord randomSpotLightSample( Light light, sampler2DArray iesProfiles, vec3 rayOrigin, vec2 ruv ) {
+
+		float radius = light.radius * sqrt( ruv.x );
+		float theta = ruv.y * 2.0 * PI;
+		float x = radius * cos( theta );
+		float y = radius * sin( theta );
+
+		vec3 u = light.u;
+		vec3 v = light.v;
+		vec3 normal = normalize( cross( u, v ) );
+
+		float angle = acos( light.coneCos );
+		float angleTan = tan( angle );
+		float startDistance = light.radius / max( angleTan, EPSILON );
+
+		vec3 randomPos = light.position - normal * startDistance + u * x + v * y;
+		vec3 toLight = randomPos - rayOrigin;
+		float lightDistSq = dot( toLight, toLight );
+		float dist = sqrt( lightDistSq );
+
+		vec3 direction = toLight / max( dist, EPSILON );
+		float cosTheta = dot( direction, normal );
+
+		float spotAttenuation = light.iesProfile != - 1 ?
+			getPhotometricAttenuation( iesProfiles, light.iesProfile, direction, normal, u, v ) :
+			getSpotAttenuation( light.coneCos, light.penumbraCos, cosTheta );
+
+		float distanceAttenuation = getDistanceAttenuation( dist, light.distance, light.decay );
+		LightRecord lightRec;
+		lightRec.type = light.type;
+		lightRec.dist = dist;
+		lightRec.direction = direction;
+		lightRec.emission = light.color * light.intensity * distanceAttenuation * spotAttenuation;
+		lightRec.pdf = 1.0;
+
+		return lightRec;
+
+	}
+
+	LightRecord randomLightSample( sampler2D lights, sampler2DArray iesProfiles, uint lightCount, vec3 rayOrigin, vec3 ruv ) {
+
+		LightRecord result;
+
+		// pick a random light
+		uint l = uint( ruv.x * float( lightCount ) );
+		Light light = readLightInfo( lights, l );
+
+		if ( light.type == SPOT_LIGHT_TYPE ) {
+
+			result = randomSpotLightSample( light, iesProfiles, rayOrigin, ruv.yz );
+
+		} else if ( light.type == POINT_LIGHT_TYPE ) {
+
+			vec3 lightRay = light.u - rayOrigin;
+			float lightDist = length( lightRay );
+			float cutoffDistance = light.distance;
+			float distanceFalloff = 1.0 / max( pow( lightDist, light.decay ), 0.01 );
+			if ( cutoffDistance > 0.0 ) {
+
+				distanceFalloff *= pow2( saturate( 1.0 - pow4( lightDist / cutoffDistance ) ) );
+
+			}
+
+			LightRecord rec;
+			rec.direction = normalize( lightRay );
+			rec.dist = length( lightRay );
+			rec.pdf = 1.0;
+			rec.emission = light.color * light.intensity * distanceFalloff;
+			rec.type = light.type;
+			result = rec;
+
+		} else if ( light.type == DIR_LIGHT_TYPE ) {
+
+			LightRecord rec;
+			rec.dist = 1e10;
+			rec.direction = light.u;
+			rec.pdf = 1.0;
+			rec.emission = light.color * light.intensity;
+			rec.type = light.type;
+
+			result = rec;
+
+		} else {
+
+			// sample the light
+			result = randomAreaLightSample( light, rayOrigin, ruv.yz );
+
+		}
+
+		return result;
+
+	}
+
+`,pr=`
+
+	vec3 sampleHemisphere( vec3 n, vec2 uv ) {
+
+		// https://www.rorydriscoll.com/2009/01/07/better-sampling/
+		// https://graphics.pixar.com/library/OrthonormalB/paper.pdf
+		float sign = n.z == 0.0 ? 1.0 : sign( n.z );
+		float a = - 1.0 / ( sign + n.z );
+		float b = n.x * n.y * a;
+		vec3 b1 = vec3( 1.0 + sign * n.x * n.x * a, sign * b, - sign * n.x );
+		vec3 b2 = vec3( b, sign + n.y * n.y * a, - n.y );
+
+		float r = sqrt( uv.x );
+		float theta = 2.0 * PI * uv.y;
+		float x = r * cos( theta );
+		float y = r * sin( theta );
+		return x * b1 + y * b2 + sqrt( 1.0 - uv.x ) * n;
+
+	}
+
+	vec2 sampleTriangle( vec2 a, vec2 b, vec2 c, vec2 r ) {
+
+		// get the edges of the triangle and the diagonal across the
+		// center of the parallelogram
+		vec2 e1 = a - b;
+		vec2 e2 = c - b;
+		vec2 diag = normalize( e1 + e2 );
+
+		// pick the point in the parallelogram
+		if ( r.x + r.y > 1.0 ) {
+
+			r = vec2( 1.0 ) - r;
+
+		}
+
+		return e1 * r.x + e2 * r.y;
+
+	}
+
+	vec2 sampleCircle( vec2 uv ) {
+
+		float angle = 2.0 * PI * uv.x;
+		float radius = sqrt( uv.y );
+		return vec2( cos( angle ), sin( angle ) ) * radius;
+
+	}
+
+	vec3 sampleSphere( vec2 uv ) {
+
+		float u = ( uv.x - 0.5 ) * 2.0;
+		float t = uv.y * PI * 2.0;
+		float f = sqrt( 1.0 - u * u );
+
+		return vec3( f * cos( t ), f * sin( t ), u );
+
+	}
+
+	vec2 sampleRegularPolygon( int sides, vec3 uvw ) {
+
+		sides = max( sides, 3 );
+
+		vec3 r = uvw;
+		float anglePerSegment = 2.0 * PI / float( sides );
+		float segment = floor( float( sides ) * r.x );
+
+		float angle1 = anglePerSegment * segment;
+		float angle2 = angle1 + anglePerSegment;
+		vec2 a = vec2( sin( angle1 ), cos( angle1 ) );
+		vec2 b = vec2( 0.0, 0.0 );
+		vec2 c = vec2( sin( angle2 ), cos( angle2 ) );
+
+		return sampleTriangle( a, b, c, r.yz );
+
+	}
+
+	// samples an aperture shape with the given number of sides. 0 means circle
+	vec2 sampleAperture( int blades, vec3 uvw ) {
+
+		return blades == 0 ?
+			sampleCircle( uvw.xy ) :
+			sampleRegularPolygon( blades, uvw );
+
+	}
+
+
+`,gr=`
+
+	bool totalInternalReflection( float cosTheta, float eta ) {
+
+		float sinTheta = sqrt( 1.0 - cosTheta * cosTheta );
+		return eta * sinTheta > 1.0;
+
+	}
+
+	// https://google.github.io/filament/Filament.md.html#materialsystem/diffusebrdf
+	float schlickFresnel( float cosine, float f0 ) {
+
+		return f0 + ( 1.0 - f0 ) * pow( 1.0 - cosine, 5.0 );
+
+	}
+
+	vec3 schlickFresnel( float cosine, vec3 f0 ) {
+
+		return f0 + ( 1.0 - f0 ) * pow( 1.0 - cosine, 5.0 );
+
+	}
+
+	vec3 schlickFresnel( float cosine, vec3 f0, vec3 f90 ) {
+
+		return f0 + ( f90 - f0 ) * pow( 1.0 - cosine, 5.0 );
+
+	}
+
+	float dielectricFresnel( float cosThetaI, float eta ) {
+
+		// https://schuttejoe.github.io/post/disneybsdf/
+		float ni = eta;
+		float nt = 1.0;
+
+		// Check for total internal reflection
+		float sinThetaISq = 1.0f - cosThetaI * cosThetaI;
+		float sinThetaTSq = eta * eta * sinThetaISq;
+		if( sinThetaTSq >= 1.0 ) {
+
+			return 1.0;
+
+		}
+
+		float sinThetaT = sqrt( sinThetaTSq );
+
+		float cosThetaT = sqrt( max( 0.0, 1.0f - sinThetaT * sinThetaT ) );
+		float rParallel = ( ( nt * cosThetaI ) - ( ni * cosThetaT ) ) / ( ( nt * cosThetaI ) + ( ni * cosThetaT ) );
+		float rPerpendicular = ( ( ni * cosThetaI ) - ( nt * cosThetaT ) ) / ( ( ni * cosThetaI ) + ( nt * cosThetaT ) );
+		return ( rParallel * rParallel + rPerpendicular * rPerpendicular ) / 2.0;
+
+	}
+
+	// https://raytracing.github.io/books/RayTracingInOneWeekend.html#dielectrics/schlickapproximation
+	float iorRatioToF0( float eta ) {
+
+		return pow( ( 1.0 - eta ) / ( 1.0 + eta ), 2.0 );
+
+	}
+
+	vec3 evaluateFresnel( float cosTheta, float eta, vec3 f0, vec3 f90 ) {
+
+		if ( totalInternalReflection( cosTheta, eta ) ) {
+
+			return f90;
+
+		}
+
+		return schlickFresnel( cosTheta, f0, f90 );
+
+	}
+
+	// TODO: disney fresnel was removed and replaced with this fresnel function to better align with
+	// the glTF but is causing blown out pixels. Should be revisited
+	// float evaluateFresnelWeight( float cosTheta, float eta, float f0 ) {
+
+	// 	if ( totalInternalReflection( cosTheta, eta ) ) {
+
+	// 		return 1.0;
+
+	// 	}
+
+	// 	return schlickFresnel( cosTheta, f0 );
+
+	// }
+
+	// https://schuttejoe.github.io/post/disneybsdf/
+	float disneyFresnel( vec3 wo, vec3 wi, vec3 wh, float f0, float eta, float metalness ) {
+
+		float dotHV = dot( wo, wh );
+		if ( totalInternalReflection( dotHV, eta ) ) {
+
+			return 1.0;
+
+		}
+
+		float dotHL = dot( wi, wh );
+		float dielectricFresnel = dielectricFresnel( abs( dotHV ), eta );
+		float metallicFresnel = schlickFresnel( dotHL, f0 );
+
+		return mix( dielectricFresnel, metallicFresnel, metalness );
+
+	}
+
+`,vr=`
+
+	// Fast arccos approximation used to remove banding artifacts caused by numerical errors in acos.
+	// This is a cubic Lagrange interpolating polynomial for x = [-1, -1/2, 0, 1/2, 1].
+	// For more information see: https://github.com/gkjohnson/three-gpu-pathtracer/pull/171#issuecomment-1152275248
+	float acosApprox( float x ) {
+
+		x = clamp( x, -1.0, 1.0 );
+		return ( - 0.69813170079773212 * x * x - 0.87266462599716477 ) * x + 1.5707963267948966;
+
+	}
+
+	// An acos with input values bound to the range [-1, 1].
+	float acosSafe( float x ) {
+
+		return acos( clamp( x, -1.0, 1.0 ) );
+
+	}
+
+	float saturateCos( float val ) {
+
+		return clamp( val, 0.001, 1.0 );
+
+	}
+
+	float square( float t ) {
+
+		return t * t;
+
+	}
+
+	vec2 square( vec2 t ) {
+
+		return t * t;
+
+	}
+
+	vec3 square( vec3 t ) {
+
+		return t * t;
+
+	}
+
+	vec4 square( vec4 t ) {
+
+		return t * t;
+
+	}
+
+	vec2 rotateVector( vec2 v, float t ) {
+
+		float ac = cos( t );
+		float as = sin( t );
+		return vec2(
+			v.x * ac - v.y * as,
+			v.x * as + v.y * ac
+		);
+
+	}
+
+	// forms a basis with the normal vector as Z
+	mat3 getBasisFromNormal( vec3 normal ) {
+
+		vec3 other;
+		if ( abs( normal.x ) > 0.5 ) {
+
+			other = vec3( 0.0, 1.0, 0.0 );
+
+		} else {
+
+			other = vec3( 1.0, 0.0, 0.0 );
+
+		}
+
+		vec3 ortho = normalize( cross( normal, other ) );
+		vec3 ortho2 = normalize( cross( normal, ortho ) );
+		return mat3( ortho2, ortho, normal );
+
+	}
+
+`,xr=`
+
+	// Finds the point where the ray intersects the plane defined by u and v and checks if this point
+	// falls in the bounds of the rectangle on that same plane.
+	// Plane intersection: https://lousodrome.net/blog/light/2020/07/03/intersection-of-a-ray-and-a-plane/
+	bool intersectsRectangle( vec3 center, vec3 normal, vec3 u, vec3 v, vec3 rayOrigin, vec3 rayDirection, inout float dist ) {
+
+		float t = dot( center - rayOrigin, normal ) / dot( rayDirection, normal );
+
+		if ( t > EPSILON ) {
+
+			vec3 p = rayOrigin + rayDirection * t;
+			vec3 vi = p - center;
+
+			// check if p falls inside the rectangle
+			float a1 = dot( u, vi );
+			if ( abs( a1 ) <= 0.5 ) {
+
+				float a2 = dot( v, vi );
+				if ( abs( a2 ) <= 0.5 ) {
+
+					dist = t;
+					return true;
+
+				}
+
+			}
+
+		}
+
+		return false;
+
+	}
+
+	// Finds the point where the ray intersects the plane defined by u and v and checks if this point
+	// falls in the bounds of the circle on that same plane. See above URL for a description of the plane intersection algorithm.
+	bool intersectsCircle( vec3 position, vec3 normal, vec3 u, vec3 v, vec3 rayOrigin, vec3 rayDirection, inout float dist ) {
+
+		float t = dot( position - rayOrigin, normal ) / dot( rayDirection, normal );
+
+		if ( t > EPSILON ) {
+
+			vec3 hit = rayOrigin + rayDirection * t;
+			vec3 vi = hit - position;
+
+			float a1 = dot( u, vi );
+			float a2 = dot( v, vi );
+
+			if( length( vec2( a1, a2 ) ) <= 0.5 ) {
+
+				dist = t;
+				return true;
+
+			}
+
+		}
+
+		return false;
+
+	}
+
+`,yr=`
+
+	// add texel fetch functions for texture arrays
+	vec4 texelFetch1D( sampler2DArray tex, int layer, uint index ) {
+
+		uint width = uint( textureSize( tex, 0 ).x );
+		uvec2 uv;
+		uv.x = index % width;
+		uv.y = index / width;
+
+		return texelFetch( tex, ivec3( uv, layer ), 0 );
+
+	}
+
+	vec4 textureSampleBarycoord( sampler2DArray tex, int layer, vec3 barycoord, uvec3 faceIndices ) {
+
+		return
+			barycoord.x * texelFetch1D( tex, layer, faceIndices.x ) +
+			barycoord.y * texelFetch1D( tex, layer, faceIndices.y ) +
+			barycoord.z * texelFetch1D( tex, layer, faceIndices.z );
+
+	}
+
+`,Ge=`
+
+	// TODO: possibly this should be renamed something related to material or path tracing logic
+
+	#ifndef RAY_OFFSET
+	#define RAY_OFFSET 1e-4
+	#endif
+
+	// adjust the hit point by the surface normal by a factor of some offset and the
+	// maximum component-wise value of the current point to accommodate floating point
+	// error as values increase.
+	vec3 stepRayOrigin( vec3 rayOrigin, vec3 rayDirection, vec3 offset, float dist ) {
+
+		vec3 point = rayOrigin + rayDirection * dist;
+		vec3 absPoint = abs( point );
+		float maxPoint = max( absPoint.x, max( absPoint.y, absPoint.z ) );
+		return point + offset * ( maxPoint + 1.0 ) * RAY_OFFSET;
+
+	}
+
+	// https://github.com/KhronosGroup/glTF/blob/main/extensions/2.0/Khronos/KHR_materials_volume/README.md#attenuation
+	vec3 transmissionAttenuation( float dist, vec3 attColor, float attDist ) {
+
+		vec3 ot = - log( attColor ) / attDist;
+		return exp( - ot * dist );
+
+	}
+
+	vec3 getHalfVector( vec3 wi, vec3 wo, float eta ) {
+
+		// get the half vector - assuming if the light incident vector is on the other side
+		// of the that it's transmissive.
+		vec3 h;
+		if ( wi.z > 0.0 ) {
+
+			h = normalize( wi + wo );
+
+		} else {
+
+			// Scale by the ior ratio to retrieve the appropriate half vector
+			// From Section 2.2 on computing the transmission half vector:
+			// https://blog.selfshadow.com/publications/s2015-shading-course/burley/s2015_pbs_disney_bsdf_notes.pdf
+			h = normalize( wi + wo * eta );
+
+		}
+
+		h *= sign( h.z );
+		return h;
+
+	}
+
+	vec3 getHalfVector( vec3 a, vec3 b ) {
+
+		return normalize( a + b );
+
+	}
+
+	// The discrepancy between interpolated surface normal and geometry normal can cause issues when a ray
+	// is cast that is on the top side of the geometry normal plane but below the surface normal plane. If
+	// we find a ray like that we ignore it to avoid artifacts.
+	// This function returns if the direction is on the same side of both planes.
+	bool isDirectionValid( vec3 direction, vec3 surfaceNormal, vec3 geometryNormal ) {
+
+		bool aboveSurfaceNormal = dot( direction, surfaceNormal ) > 0.0;
+		bool aboveGeometryNormal = dot( direction, geometryNormal ) > 0.0;
+		return aboveSurfaceNormal == aboveGeometryNormal;
+
+	}
+
+	// ray sampling x and z are swapped to align with expected background view
+	vec2 equirectDirectionToUv( vec3 direction ) {
+
+		// from Spherical.setFromCartesianCoords
+		vec2 uv = vec2( atan( direction.z, direction.x ), acos( direction.y ) );
+		uv /= vec2( 2.0 * PI, PI );
+
+		// apply adjustments to get values in range [0, 1] and y right side up
+		uv.x += 0.5;
+		uv.y = 1.0 - uv.y;
+		return uv;
+
+	}
+
+	vec3 equirectUvToDirection( vec2 uv ) {
+
+		// undo above adjustments
+		uv.x -= 0.5;
+		uv.y = 1.0 - uv.y;
+
+		// from Vector3.setFromSphericalCoords
+		float theta = uv.x * 2.0 * PI;
+		float phi = uv.y * PI;
+
+		float sinPhi = sin( phi );
+
+		return vec3( sinPhi * cos( theta ), cos( phi ), sinPhi * sin( theta ) );
+
+	}
+
+	// power heuristic for multiple importance sampling
+	float misHeuristic( float a, float b ) {
+
+		float aa = a * a;
+		float bb = b * b;
+		return aa / ( aa + bb );
+
+	}
+
+	// tentFilter from Peter Shirley's 'Realistic Ray Tracing (2nd Edition)' book, pg. 60
+	// erichlof/THREE.js-PathTracing-Renderer/
+	float tentFilter( float x ) {
+
+		return x < 0.5 ? sqrt( 2.0 * x ) - 1.0 : 1.0 - sqrt( 2.0 - ( 2.0 * x ) );
+
+	}
+`,Me=`
+
+	// https://www.shadertoy.com/view/wltcRS
+	uvec4 WHITE_NOISE_SEED;
+
+	void rng_initialize( vec2 p, int frame ) {
+
+		// white noise seed
+		WHITE_NOISE_SEED = uvec4( p, uint( frame ), uint( p.x ) + uint( p.y ) );
+
+	}
+
+	// https://www.pcg-random.org/
+	void pcg4d( inout uvec4 v ) {
+
+		v = v * 1664525u + 1013904223u;
+		v.x += v.y * v.w;
+		v.y += v.z * v.x;
+		v.z += v.x * v.y;
+		v.w += v.y * v.z;
+		v = v ^ ( v >> 16u );
+		v.x += v.y*v.w;
+		v.y += v.z*v.x;
+		v.z += v.x*v.y;
+		v.w += v.y*v.z;
+
+	}
+
+	// returns [ 0, 1 ]
+	float pcgRand() {
+
+		pcg4d( WHITE_NOISE_SEED );
+		return float( WHITE_NOISE_SEED.x ) / float( 0xffffffffu );
+
+	}
+
+	vec2 pcgRand2() {
+
+		pcg4d( WHITE_NOISE_SEED );
+		return vec2( WHITE_NOISE_SEED.xy ) / float(0xffffffffu);
+
+	}
+
+	vec3 pcgRand3() {
+
+		pcg4d( WHITE_NOISE_SEED );
+		return vec3( WHITE_NOISE_SEED.xyz ) / float( 0xffffffffu );
+
+	}
+
+	vec4 pcgRand4() {
+
+		pcg4d( WHITE_NOISE_SEED );
+		return vec4( WHITE_NOISE_SEED ) / float( 0xffffffffu );
+
+	}
+`,br=`
+
+	uniform sampler2D stratifiedTexture;
+	uniform sampler2D stratifiedOffsetTexture;
+
+	uint sobolPixelIndex = 0u;
+	uint sobolPathIndex = 0u;
+	uint sobolBounceIndex = 0u;
+	vec4 pixelSeed = vec4( 0 );
+
+	vec4 rand4( int v ) {
+
+		ivec2 uv = ivec2( v, sobolBounceIndex );
+		vec4 stratifiedSample = texelFetch( stratifiedTexture, uv, 0 );
+		return fract( stratifiedSample + pixelSeed.r ); // blue noise + stratified samples
+
+	}
+
+	vec3 rand3( int v ) {
+
+		return rand4( v ).xyz;
+
+	}
+
+	vec2 rand2( int v ) {
+
+		return rand4( v ).xy;
+
+	}
+
+	float rand( int v ) {
+
+		return rand4( v ).x;
+
+	}
+
+	void rng_initialize( vec2 screenCoord, int frame ) {
+
+		// tile the small noise texture across the entire screen
+		ivec2 noiseSize = ivec2( textureSize( stratifiedOffsetTexture, 0 ) );
+		ivec2 pixel = ivec2( screenCoord.xy ) % noiseSize;
+		vec2 pixelWidth = 1.0 / vec2( noiseSize );
+		vec2 uv = vec2( pixel ) * pixelWidth + pixelWidth * 0.5;
+
+		// note that using "texelFetch" here seems to break Android for some reason
+		pixelSeed = texture( stratifiedOffsetTexture, uv );
+
+	}
+
+`,Tr=`
+
+	// diffuse
+	float diffuseEval( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf, inout vec3 color ) {
+
+		// https://schuttejoe.github.io/post/disneybsdf/
+		float fl = schlickFresnel( wi.z, 0.0 );
+		float fv = schlickFresnel( wo.z, 0.0 );
+
+		float metalFactor = ( 1.0 - surf.metalness );
+		float transFactor = ( 1.0 - surf.transmission );
+		float rr = 0.5 + 2.0 * surf.roughness * fl * fl;
+		float retro = rr * ( fl + fv + fl * fv * ( rr - 1.0f ) );
+		float lambert = ( 1.0f - 0.5f * fl ) * ( 1.0f - 0.5f * fv );
+
+		// TODO: subsurface approx?
+
+		// float F = evaluateFresnelWeight( dot( wo, wh ), surf.eta, surf.f0 );
+		float F = disneyFresnel( wo, wi, wh, surf.f0, surf.eta, surf.metalness );
+		color = ( 1.0 - F ) * transFactor * metalFactor * wi.z * surf.color * ( retro + lambert ) / PI;
+
+		return wi.z / PI;
+
+	}
+
+	vec3 diffuseDirection( vec3 wo, SurfaceRecord surf ) {
+
+		vec3 lightDirection = sampleSphere( rand2( 11 ) );
+		lightDirection.z += 1.0;
+		lightDirection = normalize( lightDirection );
+
+		return lightDirection;
+
+	}
+
+	// specular
+	float specularEval( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf, inout vec3 color ) {
+
+		// if roughness is set to 0 then D === NaN which results in black pixels
+		float metalness = surf.metalness;
+		float roughness = surf.filteredRoughness;
+
+		float eta = surf.eta;
+		float f0 = surf.f0;
+
+		vec3 f0Color = mix( f0 * surf.specularColor * surf.specularIntensity, surf.color, surf.metalness );
+		vec3 f90Color = vec3( mix( surf.specularIntensity, 1.0, surf.metalness ) );
+		vec3 F = evaluateFresnel( dot( wo, wh ), eta, f0Color, f90Color );
+
+		vec3 iridescenceF = evalIridescence( 1.0, surf.iridescenceIor, dot( wi, wh ), surf.iridescenceThickness, f0Color );
+		F = mix( F, iridescenceF,  surf.iridescence );
+
+		// PDF
+		// See 14.1.1 Microfacet BxDFs in https://www.pbr-book.org/
+		float incidentTheta = acos( wo.z );
+		float G = ggxShadowMaskG2( wi, wo, roughness );
+		float D = ggxDistribution( wh, roughness );
+		float G1 = ggxShadowMaskG1( incidentTheta, roughness );
+		float ggxPdf = D * G1 * max( 0.0, abs( dot( wo, wh ) ) ) / abs ( wo.z );
+
+		color = wi.z * F * G * D / ( 4.0 * abs( wi.z * wo.z ) );
+		return ggxPdf / ( 4.0 * dot( wo, wh ) );
+
+	}
+
+	vec3 specularDirection( vec3 wo, SurfaceRecord surf ) {
+
+		// sample ggx vndf distribution which gives a new normal
+		float roughness = surf.filteredRoughness;
+		vec3 halfVector = ggxDirection(
+			wo,
+			vec2( roughness ),
+			rand2( 12 )
+		);
+
+		// apply to new ray by reflecting off the new normal
+		return - reflect( wo, halfVector );
+
+	}
+
+
+	// transmission
+	/*
+	float transmissionEval( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf, inout vec3 color ) {
+
+		// See section 4.2 in https://www.cs.cornell.edu/~srm/publications/EGSR07-btdf.pdf
+
+		float filteredRoughness = surf.filteredRoughness;
+		float eta = surf.eta;
+		bool frontFace = surf.frontFace;
+		bool thinFilm = surf.thinFilm;
+
+		color = surf.transmission * surf.color;
+
+		float denom = pow( eta * dot( wi, wh ) + dot( wo, wh ), 2.0 );
+		return ggxPDF( wo, wh, filteredRoughness ) / denom;
+
+	}
+
+	vec3 transmissionDirection( vec3 wo, SurfaceRecord surf ) {
+
+		float filteredRoughness = surf.filteredRoughness;
+		float eta = surf.eta;
+		bool frontFace = surf.frontFace;
+
+		// sample ggx vndf distribution which gives a new normal
+		vec3 halfVector = ggxDirection(
+			wo,
+			vec2( filteredRoughness ),
+			rand2( 13 )
+		);
+
+		vec3 lightDirection = refract( normalize( - wo ), halfVector, eta );
+		if ( surf.thinFilm ) {
+
+			lightDirection = - refract( normalize( - lightDirection ), - vec3( 0.0, 0.0, 1.0 ), 1.0 / eta );
+
+		}
+
+		return normalize( lightDirection );
+
+	}
+	*/
+
+	// TODO: This is just using a basic cosine-weighted specular distribution with an
+	// incorrect PDF value at the moment. Update it to correctly use a GGX distribution
+	float transmissionEval( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf, inout vec3 color ) {
+
+		color = surf.transmission * surf.color;
+
+		// PDF
+		// float F = evaluateFresnelWeight( dot( wo, wh ), surf.eta, surf.f0 );
+		// float F = disneyFresnel( wo, wi, wh, surf.f0, surf.eta, surf.metalness );
+		// if ( F >= 1.0 ) {
+
+		// 	return 0.0;
+
+		// }
+
+		// return 1.0 / ( 1.0 - F );
+
+		// reverted to previous to transmission. The above was causing black pixels
+		float eta = surf.eta;
+		float f0 = surf.f0;
+		float cosTheta = min( wo.z, 1.0 );
+		float sinTheta = sqrt( 1.0 - cosTheta * cosTheta );
+		float reflectance = schlickFresnel( cosTheta, f0 );
+		bool cannotRefract = eta * sinTheta > 1.0;
+		if ( cannotRefract ) {
+
+			return 0.0;
+
+		}
+
+		return 1.0 / ( 1.0 - reflectance );
+
+	}
+
+	vec3 transmissionDirection( vec3 wo, SurfaceRecord surf ) {
+
+		float roughness = surf.filteredRoughness;
+		float eta = surf.eta;
+		vec3 halfVector = normalize( vec3( 0.0, 0.0, 1.0 ) + sampleSphere( rand2( 13 ) ) * roughness );
+		vec3 lightDirection = refract( normalize( - wo ), halfVector, eta );
+
+		if ( surf.thinFilm ) {
+
+			lightDirection = - refract( normalize( - lightDirection ), - vec3( 0.0, 0.0, 1.0 ), 1.0 / eta );
+
+		}
+		return normalize( lightDirection );
+
+	}
+
+	// clearcoat
+	float clearcoatEval( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf, inout vec3 color ) {
+
+		float ior = 1.5;
+		float f0 = iorRatioToF0( ior );
+		bool frontFace = surf.frontFace;
+		float roughness = surf.filteredClearcoatRoughness;
+
+		float eta = frontFace ? 1.0 / ior : ior;
+		float G = ggxShadowMaskG2( wi, wo, roughness );
+		float D = ggxDistribution( wh, roughness );
+		float F = schlickFresnel( dot( wi, wh ), f0 );
+
+		float fClearcoat = F * D * G / ( 4.0 * abs( wi.z * wo.z ) );
+		color = color * ( 1.0 - surf.clearcoat * F ) + fClearcoat * surf.clearcoat * wi.z;
+
+		// PDF
+		// See equation (27) in http://jcgt.org/published/0003/02/03/
+		return ggxPDF( wo, wh, roughness ) / ( 4.0 * dot( wi, wh ) );
+
+	}
+
+	vec3 clearcoatDirection( vec3 wo, SurfaceRecord surf ) {
+
+		// sample ggx vndf distribution which gives a new normal
+		float roughness = surf.filteredClearcoatRoughness;
+		vec3 halfVector = ggxDirection(
+			wo,
+			vec2( roughness ),
+			rand2( 14 )
+		);
+
+		// apply to new ray by reflecting off the new normal
+		return - reflect( wo, halfVector );
+
+	}
+
+	// sheen
+	vec3 sheenColor( vec3 wo, vec3 wi, vec3 wh, SurfaceRecord surf ) {
+
+		float cosThetaO = saturateCos( wo.z );
+		float cosThetaI = saturateCos( wi.z );
+		float cosThetaH = wh.z;
+
+		float D = velvetD( cosThetaH, surf.sheenRoughness );
+		float G = velvetG( cosThetaO, cosThetaI, surf.sheenRoughness );
+
+		// See equation (1) in http://www.aconty.com/pdf/s2017_pbs_imageworks_sheen.pdf
+		vec3 color = surf.sheenColor;
+		color *= D * G / ( 4.0 * abs( cosThetaO * cosThetaI ) );
+		color *= wi.z;
+
+		return color;
+
+	}
+
+	// bsdf
+	void getLobeWeights(
+		vec3 wo, vec3 wi, vec3 wh, vec3 clearcoatWo, SurfaceRecord surf,
+		inout float diffuseWeight, inout float specularWeight, inout float transmissionWeight, inout float clearcoatWeight
+	) {
+
+		float metalness = surf.metalness;
+		float transmission = surf.transmission;
+		// float fEstimate = evaluateFresnelWeight( dot( wo, wh ), surf.eta, surf.f0 );
+		float fEstimate = disneyFresnel( wo, wi, wh, surf.f0, surf.eta, surf.metalness );
+
+		float transSpecularProb = mix( max( 0.25, fEstimate ), 1.0, metalness );
+		float diffSpecularProb = 0.5 + 0.5 * metalness;
+
+		diffuseWeight = ( 1.0 - transmission ) * ( 1.0 - diffSpecularProb );
+		specularWeight = transmission * transSpecularProb + ( 1.0 - transmission ) * diffSpecularProb;
+		transmissionWeight = transmission * ( 1.0 - transSpecularProb );
+		clearcoatWeight = surf.clearcoat * schlickFresnel( clearcoatWo.z, 0.04 );
+
+		float totalWeight = diffuseWeight + specularWeight + transmissionWeight + clearcoatWeight;
+		diffuseWeight /= totalWeight;
+		specularWeight /= totalWeight;
+		transmissionWeight /= totalWeight;
+		clearcoatWeight /= totalWeight;
+	}
+
+	float bsdfEval(
+		vec3 wo, vec3 clearcoatWo, vec3 wi, vec3 clearcoatWi, SurfaceRecord surf,
+		float diffuseWeight, float specularWeight, float transmissionWeight, float clearcoatWeight, inout float specularPdf, inout vec3 color
+	) {
+
+		float metalness = surf.metalness;
+		float transmission = surf.transmission;
+
+		float spdf = 0.0;
+		float dpdf = 0.0;
+		float tpdf = 0.0;
+		float cpdf = 0.0;
+		color = vec3( 0.0 );
+
+		vec3 halfVector = getHalfVector( wi, wo, surf.eta );
+
+		// diffuse
+		if ( diffuseWeight > 0.0 && wi.z > 0.0 ) {
+
+			dpdf = diffuseEval( wo, wi, halfVector, surf, color );
+			color *= 1.0 - surf.transmission;
+
+		}
+
+		// ggx specular
+		if ( specularWeight > 0.0 && wi.z > 0.0 ) {
+
+			vec3 outColor;
+			spdf = specularEval( wo, wi, getHalfVector( wi, wo ), surf, outColor );
+			color += outColor;
+
+		}
+
+		// transmission
+		if ( transmissionWeight > 0.0 && wi.z < 0.0 ) {
+
+			tpdf = transmissionEval( wo, wi, halfVector, surf, color );
+
+		}
+
+		// sheen
+		color *= mix( 1.0, sheenAlbedoScaling( wo, wi, surf ), surf.sheen );
+		color += sheenColor( wo, wi, halfVector, surf ) * surf.sheen;
+
+		// clearcoat
+		if ( clearcoatWi.z >= 0.0 && clearcoatWeight > 0.0 ) {
+
+			vec3 clearcoatHalfVector = getHalfVector( clearcoatWo, clearcoatWi );
+			cpdf = clearcoatEval( clearcoatWo, clearcoatWi, clearcoatHalfVector, surf, color );
+
+		}
+
+		float pdf =
+			dpdf * diffuseWeight
+			+ spdf * specularWeight
+			+ tpdf * transmissionWeight
+			+ cpdf * clearcoatWeight;
+
+		// retrieve specular rays for the shadows flag
+		specularPdf = spdf * specularWeight + cpdf * clearcoatWeight;
+
+		return pdf;
+
+	}
+
+	float bsdfResult( vec3 worldWo, vec3 worldWi, SurfaceRecord surf, inout vec3 color ) {
+
+		if ( surf.volumeParticle ) {
+
+			color = surf.color / ( 4.0 * PI );
+			return 1.0 / ( 4.0 * PI );
+
+		}
+
+		vec3 wo = normalize( surf.normalInvBasis * worldWo );
+		vec3 wi = normalize( surf.normalInvBasis * worldWi );
+
+		vec3 clearcoatWo = normalize( surf.clearcoatInvBasis * worldWo );
+		vec3 clearcoatWi = normalize( surf.clearcoatInvBasis * worldWi );
+
+		vec3 wh = getHalfVector( wo, wi, surf.eta );
+		float diffuseWeight;
+		float specularWeight;
+		float transmissionWeight;
+		float clearcoatWeight;
+		getLobeWeights( wo, wi, wh, clearcoatWo, surf, diffuseWeight, specularWeight, transmissionWeight, clearcoatWeight );
+
+		float specularPdf;
+		return bsdfEval( wo, clearcoatWo, wi, clearcoatWi, surf, diffuseWeight, specularWeight, transmissionWeight, clearcoatWeight, specularPdf, color );
+
+	}
+
+	ScatterRecord bsdfSample( vec3 worldWo, SurfaceRecord surf ) {
+
+		if ( surf.volumeParticle ) {
+
+			ScatterRecord sampleRec;
+			sampleRec.specularPdf = 0.0;
+			sampleRec.pdf = 1.0 / ( 4.0 * PI );
+			sampleRec.direction = sampleSphere( rand2( 16 ) );
+			sampleRec.color = surf.color / ( 4.0 * PI );
+			return sampleRec;
+
+		}
+
+		vec3 wo = normalize( surf.normalInvBasis * worldWo );
+		vec3 clearcoatWo = normalize( surf.clearcoatInvBasis * worldWo );
+		mat3 normalBasis = surf.normalBasis;
+		mat3 invBasis = surf.normalInvBasis;
+		mat3 clearcoatNormalBasis = surf.clearcoatBasis;
+		mat3 clearcoatInvBasis = surf.clearcoatInvBasis;
+
+		float diffuseWeight;
+		float specularWeight;
+		float transmissionWeight;
+		float clearcoatWeight;
+		// using normal and basically-reflected ray since we don't have proper half vector here
+		getLobeWeights( wo, wo, vec3( 0, 0, 1 ), clearcoatWo, surf, diffuseWeight, specularWeight, transmissionWeight, clearcoatWeight );
+
+		float pdf[4];
+		pdf[0] = diffuseWeight;
+		pdf[1] = specularWeight;
+		pdf[2] = transmissionWeight;
+		pdf[3] = clearcoatWeight;
+
+		float cdf[4];
+		cdf[0] = pdf[0];
+		cdf[1] = pdf[1] + cdf[0];
+		cdf[2] = pdf[2] + cdf[1];
+		cdf[3] = pdf[3] + cdf[2];
+
+		if( cdf[3] != 0.0 ) {
+
+			float invMaxCdf = 1.0 / cdf[3];
+			cdf[0] *= invMaxCdf;
+			cdf[1] *= invMaxCdf;
+			cdf[2] *= invMaxCdf;
+			cdf[3] *= invMaxCdf;
+
+		} else {
+
+			cdf[0] = 1.0;
+			cdf[1] = 0.0;
+			cdf[2] = 0.0;
+			cdf[3] = 0.0;
+
+		}
+
+		vec3 wi;
+		vec3 clearcoatWi;
+
+		float r = rand( 15 );
+		if ( r <= cdf[0] ) { // diffuse
+
+			wi = diffuseDirection( wo, surf );
+			clearcoatWi = normalize( clearcoatInvBasis * normalize( normalBasis * wi ) );
+
+		} else if ( r <= cdf[1] ) { // specular
+
+			wi = specularDirection( wo, surf );
+			clearcoatWi = normalize( clearcoatInvBasis * normalize( normalBasis * wi ) );
+
+		} else if ( r <= cdf[2] ) { // transmission / refraction
+
+			wi = transmissionDirection( wo, surf );
+			clearcoatWi = normalize( clearcoatInvBasis * normalize( normalBasis * wi ) );
+
+		} else if ( r <= cdf[3] ) { // clearcoat
+
+			clearcoatWi = clearcoatDirection( clearcoatWo, surf );
+			wi = normalize( invBasis * normalize( clearcoatNormalBasis * clearcoatWi ) );
+
+		}
+
+		ScatterRecord result;
+		result.pdf = bsdfEval( wo, clearcoatWo, wi, clearcoatWi, surf, diffuseWeight, specularWeight, transmissionWeight, clearcoatWeight, result.specularPdf, result.color );
+		result.direction = normalize( surf.normalBasis * wi );
+
+		return result;
+
+	}
+
+`,wr=`
+
+	// returns the hit distance given the material density
+	float intersectFogVolume( Material material, float u ) {
+
+		// https://raytracing.github.io/books/RayTracingTheNextWeek.html#volumes/constantdensitymediums
+		return material.opacity == 0.0 ? INFINITY : ( - 1.0 / material.opacity ) * log( u );
+
+	}
+
+	ScatterRecord sampleFogVolume( SurfaceRecord surf, vec2 uv ) {
+
+		ScatterRecord sampleRec;
+		sampleRec.specularPdf = 0.0;
+		sampleRec.pdf = 1.0 / ( 2.0 * PI );
+		sampleRec.direction = sampleSphere( uv );
+		sampleRec.color = surf.color;
+		return sampleRec;
+
+	}
+
+`,Rr=`
+
+	// The GGX functions provide sampling and distribution information for normals as output so
+	// in order to get probability of scatter direction the half vector must be computed and provided.
+	// [0] https://www.cs.cornell.edu/~srm/publications/EGSR07-btdf.pdf
+	// [1] https://hal.archives-ouvertes.fr/hal-01509746/document
+	// [2] http://jcgt.org/published/0007/04/01/
+	// [4] http://jcgt.org/published/0003/02/03/
+
+	// trowbridge-reitz === GGX === GTR
+
+	vec3 ggxDirection( vec3 incidentDir, vec2 roughness, vec2 uv ) {
+
+		// TODO: try GGXVNDF implementation from reference [2], here. Needs to update ggxDistribution
+		// function below, as well
+
+		// Implementation from reference [1]
+		// stretch view
+		vec3 V = normalize( vec3( roughness * incidentDir.xy, incidentDir.z ) );
+
+		// orthonormal basis
+		vec3 T1 = ( V.z < 0.9999 ) ? normalize( cross( V, vec3( 0.0, 0.0, 1.0 ) ) ) : vec3( 1.0, 0.0, 0.0 );
+		vec3 T2 = cross( T1, V );
+
+		// sample point with polar coordinates (r, phi)
+		float a = 1.0 / ( 1.0 + V.z );
+		float r = sqrt( uv.x );
+		float phi = ( uv.y < a ) ? uv.y / a * PI : PI + ( uv.y - a ) / ( 1.0 - a ) * PI;
+		float P1 = r * cos( phi );
+		float P2 = r * sin( phi ) * ( ( uv.y < a ) ? 1.0 : V.z );
+
+		// compute normal
+		vec3 N = P1 * T1 + P2 * T2 + V * sqrt( max( 0.0, 1.0 - P1 * P1 - P2 * P2 ) );
+
+		// unstretch
+		N = normalize( vec3( roughness * N.xy, max( 0.0, N.z ) ) );
+
+		return N;
+
+	}
+
+	// Below are PDF and related functions for use in a Monte Carlo path tracer
+	// as specified in Appendix B of the following paper
+	// See equation (34) from reference [0]
+	float ggxLamda( float theta, float roughness ) {
+
+		float tanTheta = tan( theta );
+		float tanTheta2 = tanTheta * tanTheta;
+		float alpha2 = roughness * roughness;
+
+		float numerator = - 1.0 + sqrt( 1.0 + alpha2 * tanTheta2 );
+		return numerator / 2.0;
+
+	}
+
+	// See equation (34) from reference [0]
+	float ggxShadowMaskG1( float theta, float roughness ) {
+
+		return 1.0 / ( 1.0 + ggxLamda( theta, roughness ) );
+
+	}
+
+	// See equation (125) from reference [4]
+	float ggxShadowMaskG2( vec3 wi, vec3 wo, float roughness ) {
+
+		float incidentTheta = acos( wi.z );
+		float scatterTheta = acos( wo.z );
+		return 1.0 / ( 1.0 + ggxLamda( incidentTheta, roughness ) + ggxLamda( scatterTheta, roughness ) );
+
+	}
+
+	// See equation (33) from reference [0]
+	float ggxDistribution( vec3 halfVector, float roughness ) {
+
+		float a2 = roughness * roughness;
+		a2 = max( EPSILON, a2 );
+		float cosTheta = halfVector.z;
+		float cosTheta4 = pow( cosTheta, 4.0 );
+
+		if ( cosTheta == 0.0 ) return 0.0;
+
+		float theta = acosSafe( halfVector.z );
+		float tanTheta = tan( theta );
+		float tanTheta2 = pow( tanTheta, 2.0 );
+
+		float denom = PI * cosTheta4 * pow( a2 + tanTheta2, 2.0 );
+		return ( a2 / denom );
+
+	}
+
+	// See equation (3) from reference [2]
+	float ggxPDF( vec3 wi, vec3 halfVector, float roughness ) {
+
+		float incidentTheta = acos( wi.z );
+		float D = ggxDistribution( halfVector, roughness );
+		float G1 = ggxShadowMaskG1( incidentTheta, roughness );
+
+		return D * G1 * max( 0.0, dot( wi, halfVector ) ) / wi.z;
+
+	}
+
+`,Ir=`
+
+	// XYZ to sRGB color space
+	const mat3 XYZ_TO_REC709 = mat3(
+		3.2404542, -0.9692660,  0.0556434,
+		-1.5371385,  1.8760108, -0.2040259,
+		-0.4985314,  0.0415560,  1.0572252
+	);
+
+	vec3 fresnel0ToIor( vec3 fresnel0 ) {
+
+		vec3 sqrtF0 = sqrt( fresnel0 );
+		return ( vec3( 1.0 ) + sqrtF0 ) / ( vec3( 1.0 ) - sqrtF0 );
+
+	}
+
+	// Conversion FO/IOR
+	vec3 iorToFresnel0( vec3 transmittedIor, float incidentIor ) {
+
+		return square( ( transmittedIor - vec3( incidentIor ) ) / ( transmittedIor + vec3( incidentIor ) ) );
+
+	}
+
+	// ior is a value between 1.0 and 3.0. 1.0 is air interface
+	float iorToFresnel0( float transmittedIor, float incidentIor ) {
+
+		return square( ( transmittedIor - incidentIor ) / ( transmittedIor + incidentIor ) );
+
+	}
+
+	// Fresnel equations for dielectric/dielectric interfaces. See https://belcour.github.io/blog/research/2017/05/01/brdf-thin-film.html
+	vec3 evalSensitivity( float OPD, vec3 shift ) {
+
+		float phase = 2.0 * PI * OPD * 1.0e-9;
+
+		vec3 val = vec3( 5.4856e-13, 4.4201e-13, 5.2481e-13 );
+		vec3 pos = vec3( 1.6810e+06, 1.7953e+06, 2.2084e+06 );
+		vec3 var = vec3( 4.3278e+09, 9.3046e+09, 6.6121e+09 );
+
+		vec3 xyz = val * sqrt( 2.0 * PI * var ) * cos( pos * phase + shift ) * exp( - square( phase ) * var );
+		xyz.x += 9.7470e-14 * sqrt( 2.0 * PI * 4.5282e+09 ) * cos( 2.2399e+06 * phase + shift[ 0 ] ) * exp( - 4.5282e+09 * square( phase ) );
+		xyz /= 1.0685e-7;
+
+		vec3 srgb = XYZ_TO_REC709 * xyz;
+		return srgb;
+
+	}
+
+	// See Section 4. Analytic Spectral Integration, A Practical Extension to Microfacet Theory for the Modeling of Varying Iridescence, https://hal.archives-ouvertes.fr/hal-01518344/document
+	vec3 evalIridescence( float outsideIOR, float eta2, float cosTheta1, float thinFilmThickness, vec3 baseF0 ) {
+
+		vec3 I;
+
+		// Force iridescenceIor -> outsideIOR when thinFilmThickness -> 0.0
+		float iridescenceIor = mix( outsideIOR, eta2, smoothstep( 0.0, 0.03, thinFilmThickness ) );
+
+		// Evaluate the cosTheta on the base layer (Snell law)
+		float sinTheta2Sq = square( outsideIOR / iridescenceIor ) * ( 1.0 - square( cosTheta1 ) );
+
+		// Handle TIR:
+		float cosTheta2Sq = 1.0 - sinTheta2Sq;
+		if ( cosTheta2Sq < 0.0 ) {
+
+			return vec3( 1.0 );
+
+		}
+
+		float cosTheta2 = sqrt( cosTheta2Sq );
+
+		// First interface
+		float R0 = iorToFresnel0( iridescenceIor, outsideIOR );
+		float R12 = schlickFresnel( cosTheta1, R0 );
+		float R21 = R12;
+		float T121 = 1.0 - R12;
+		float phi12 = 0.0;
+		if ( iridescenceIor < outsideIOR ) {
+
+			phi12 = PI;
+
+		}
+
+		float phi21 = PI - phi12;
+
+		// Second interface
+		vec3 baseIOR = fresnel0ToIor( clamp( baseF0, 0.0, 0.9999 ) ); // guard against 1.0
+		vec3 R1 = iorToFresnel0( baseIOR, iridescenceIor );
+		vec3 R23 = schlickFresnel( cosTheta2, R1 );
+		vec3 phi23 = vec3( 0.0 );
+		if ( baseIOR[0] < iridescenceIor ) {
+
+			phi23[ 0 ] = PI;
+
+		}
+
+		if ( baseIOR[1] < iridescenceIor ) {
+
+			phi23[ 1 ] = PI;
+
+		}
+
+		if ( baseIOR[2] < iridescenceIor ) {
+
+			phi23[ 2 ] = PI;
+
+		}
+
+		// Phase shift
+		float OPD = 2.0 * iridescenceIor * thinFilmThickness * cosTheta2;
+		vec3 phi = vec3( phi21 ) + phi23;
+
+		// Compound terms
+		vec3 R123 = clamp( R12 * R23, 1e-5, 0.9999 );
+		vec3 r123 = sqrt( R123 );
+		vec3 Rs = square( T121 ) * R23 / ( vec3( 1.0 ) - R123 );
+
+		// Reflectance term for m = 0 (DC term amplitude)
+		vec3 C0 = R12 + Rs;
+		I = C0;
+
+		// Reflectance term for m > 0 (pairs of diracs)
+		vec3 Cm = Rs - T121;
+		for ( int m = 1; m <= 2; ++ m ) {
+
+			Cm *= r123;
+			vec3 Sm = 2.0 * evalSensitivity( float( m ) * OPD, float( m ) * phi );
+			I += Cm * Sm;
+
+		}
+
+		// Since out of gamut colors might be produced, negative color values are clamped to 0.
+		return max( I, vec3( 0.0 ) );
+
+	}
+
+`,Sr=`
+
+	// See equation (2) in http://www.aconty.com/pdf/s2017_pbs_imageworks_sheen.pdf
+	float velvetD( float cosThetaH, float roughness ) {
+
+		float alpha = max( roughness, 0.07 );
+		alpha = alpha * alpha;
+
+		float invAlpha = 1.0 / alpha;
+
+		float sqrCosThetaH = cosThetaH * cosThetaH;
+		float sinThetaH = max( 1.0 - sqrCosThetaH, 0.001 );
+
+		return ( 2.0 + invAlpha ) * pow( sinThetaH, 0.5 * invAlpha ) / ( 2.0 * PI );
+
+	}
+
+	float velvetParamsInterpolate( int i, float oneMinusAlphaSquared ) {
+
+		const float p0[5] = float[5]( 25.3245, 3.32435, 0.16801, -1.27393, -4.85967 );
+		const float p1[5] = float[5]( 21.5473, 3.82987, 0.19823, -1.97760, -4.32054 );
+
+		return mix( p1[i], p0[i], oneMinusAlphaSquared );
+
+	}
+
+	float velvetL( float x, float alpha ) {
+
+		float oneMinusAlpha = 1.0 - alpha;
+		float oneMinusAlphaSquared = oneMinusAlpha * oneMinusAlpha;
+
+		float a = velvetParamsInterpolate( 0, oneMinusAlphaSquared );
+		float b = velvetParamsInterpolate( 1, oneMinusAlphaSquared );
+		float c = velvetParamsInterpolate( 2, oneMinusAlphaSquared );
+		float d = velvetParamsInterpolate( 3, oneMinusAlphaSquared );
+		float e = velvetParamsInterpolate( 4, oneMinusAlphaSquared );
+
+		return a / ( 1.0 + b * pow( abs( x ), c ) ) + d * x + e;
+
+	}
+
+	// See equation (3) in http://www.aconty.com/pdf/s2017_pbs_imageworks_sheen.pdf
+	float velvetLambda( float cosTheta, float alpha ) {
+
+		return abs( cosTheta ) < 0.5 ? exp( velvetL( cosTheta, alpha ) ) : exp( 2.0 * velvetL( 0.5, alpha ) - velvetL( 1.0 - cosTheta, alpha ) );
+
+	}
+
+	// See Section 3, Shadowing Term, in http://www.aconty.com/pdf/s2017_pbs_imageworks_sheen.pdf
+	float velvetG( float cosThetaO, float cosThetaI, float roughness ) {
+
+		float alpha = max( roughness, 0.07 );
+		alpha = alpha * alpha;
+
+		return 1.0 / ( 1.0 + velvetLambda( cosThetaO, alpha ) + velvetLambda( cosThetaI, alpha ) );
+
+	}
+
+	float directionalAlbedoSheen( float cosTheta, float alpha ) {
+
+		cosTheta = saturate( cosTheta );
+
+		float c = 1.0 - cosTheta;
+		float c3 = c * c * c;
+
+		return 0.65584461 * c3 + 1.0 / ( 4.16526551 + exp( -7.97291361 * sqrt( alpha ) + 6.33516894 ) );
+
+	}
+
+	float sheenAlbedoScaling( vec3 wo, vec3 wi, SurfaceRecord surf ) {
+
+		float alpha = max( surf.sheenRoughness, 0.07 );
+		alpha = alpha * alpha;
+
+		float maxSheenColor = max( max( surf.sheenColor.r, surf.sheenColor.g ), surf.sheenColor.b );
+
+		float eWo = directionalAlbedoSheen( saturateCos( wo.z ), alpha );
+		float eWi = directionalAlbedoSheen( saturateCos( wi.z ), alpha );
+
+		return min( 1.0 - maxSheenColor * eWo, 1.0 - maxSheenColor * eWi );
+
+	}
+
+	// See Section 5, Layering, in http://www.aconty.com/pdf/s2017_pbs_imageworks_sheen.pdf
+	float sheenAlbedoScaling( vec3 wo, SurfaceRecord surf ) {
+
+		float alpha = max( surf.sheenRoughness, 0.07 );
+		alpha = alpha * alpha;
+
+		float maxSheenColor = max( max( surf.sheenColor.r, surf.sheenColor.g ), surf.sheenColor.b );
+
+		float eWo = directionalAlbedoSheen( saturateCos( wo.z ), alpha );
+
+		return 1.0 - maxSheenColor * eWo;
+
+	}
+
+`,_r=`
+
+#ifndef FOG_CHECK_ITERATIONS
+#define FOG_CHECK_ITERATIONS 30
+#endif
+
+// returns whether the given material is a fog material or not
+bool isMaterialFogVolume( sampler2D materials, uint materialIndex ) {
+
+	uint i = materialIndex * 45u;
+	vec4 s14 = texelFetch1D( materials, i + 14u );
+	return bool( int( s14.b ) & 4 );
+
+}
+
+// returns true if we're within the first fog volume we hit
+bool bvhIntersectFogVolumeHit(
+	vec3 rayOrigin, vec3 rayDirection,
+	usampler2D materialIndexAttribute, sampler2D materials,
+	inout Material material
+) {
+
+	material.fogVolume = false;
+
+	for ( int i = 0; i < FOG_CHECK_ITERATIONS; i ++ ) {
+
+		// find nearest hit
+		uvec4 faceIndices = uvec4( 0u );
+		vec3 faceNormal = vec3( 0.0, 0.0, 1.0 );
+		vec3 barycoord = vec3( 0.0 );
+		float side = 1.0;
+		float dist = 0.0;
+		bool hit = bvhIntersectFirstHit( bvh, rayOrigin, rayDirection, faceIndices, faceNormal, barycoord, side, dist );
+		if ( hit ) {
+
+			// if it's a fog volume return whether we hit the front or back face
+			uint materialIndex = uTexelFetch1D( materialIndexAttribute, faceIndices.x ).r;
+			if ( isMaterialFogVolume( materials, materialIndex ) ) {
+
+				material = readMaterialInfo( materials, materialIndex );
+				return side == - 1.0;
+
+			} else {
+
+				// move the ray forward
+				rayOrigin = stepRayOrigin( rayOrigin, rayDirection, - faceNormal, dist );
+
+			}
+
+		} else {
+
+			return false;
+
+		}
+
+	}
+
+	return false;
+
+}
+
+`,Fr=`
+
+	// step through multiple surface hits and accumulate color attenuation based on transmissive surfaces
+	// returns true if a solid surface was hit
+	bool attenuateHit(
+		RenderState state,
+		Ray ray, float rayDist,
+		out vec3 color
+	) {
+
+		// store the original bounce index so we can reset it after
+		uint originalBounceIndex = sobolBounceIndex;
+
+		int traversals = state.traversals;
+		int transmissiveTraversals = state.transmissiveTraversals;
+		bool isShadowRay = state.isShadowRay;
+		Material fogMaterial = state.fogMaterial;
+
+		vec3 startPoint = ray.origin;
+
+		// hit results
+		SurfaceHit surfaceHit;
+
+		color = vec3( 1.0 );
+
+		bool result = true;
+		for ( int i = 0; i < traversals; i ++ ) {
+
+			sobolBounceIndex ++;
+
+			int hitType = traceScene( ray, fogMaterial, surfaceHit );
+
+			if ( hitType == FOG_HIT ) {
+
+				result = true;
+				break;
+
+			} else if ( hitType == SURFACE_HIT ) {
+
+				float totalDist = distance( startPoint, ray.origin + ray.direction * surfaceHit.dist );
+				if ( totalDist > rayDist ) {
+
+					result = false;
+					break;
+
+				}
+
+				// TODO: attenuate the contribution based on the PDF of the resulting ray including refraction values
+				// Should be able to work using the material BSDF functions which will take into account specularity, etc.
+				// TODO: should we account for emissive surfaces here?
+
+				uint materialIndex = uTexelFetch1D( materialIndexAttribute, surfaceHit.faceIndices.x ).r;
+				Material material = readMaterialInfo( materials, materialIndex );
+
+				// adjust the ray to the new surface
+				bool isEntering = surfaceHit.side == 1.0;
+				ray.origin = stepRayOrigin( ray.origin, ray.direction, - surfaceHit.faceNormal, surfaceHit.dist );
+
+				#if FEATURE_FOG
+
+				if ( material.fogVolume ) {
+
+					fogMaterial = material;
+					fogMaterial.fogVolume = surfaceHit.side == 1.0;
+					i -= sign( transmissiveTraversals );
+					transmissiveTraversals --;
+					continue;
+
+				}
+
+				#endif
+
+				if ( ! material.castShadow && isShadowRay ) {
+
+					continue;
+
+				}
+
+				vec2 uv = textureSampleBarycoord( attributesArray, ATTR_UV, surfaceHit.barycoord, surfaceHit.faceIndices.xyz ).xy;
+				vec4 vertexColor = textureSampleBarycoord( attributesArray, ATTR_COLOR, surfaceHit.barycoord, surfaceHit.faceIndices.xyz );
+
+				// albedo
+				vec4 albedo = vec4( material.color, material.opacity );
+				if ( material.map != - 1 ) {
+
+					vec3 uvPrime = material.mapTransform * vec3( uv, 1 );
+					albedo *= texture2D( textures, vec3( uvPrime.xy, material.map ) );
+
+				}
+
+				if ( material.vertexColors ) {
+
+					albedo *= vertexColor;
+
+				}
+
+				// alphaMap
+				if ( material.alphaMap != - 1 ) {
+
+					albedo.a *= texture2D( textures, vec3( uv, material.alphaMap ) ).x;
+
+				}
+
+				// transmission
+				float transmission = material.transmission;
+				if ( material.transmissionMap != - 1 ) {
+
+					vec3 uvPrime = material.transmissionMapTransform * vec3( uv, 1 );
+					transmission *= texture2D( textures, vec3( uvPrime.xy, material.transmissionMap ) ).r;
+
+				}
+
+				// metalness
+				float metalness = material.metalness;
+				if ( material.metalnessMap != - 1 ) {
+
+					vec3 uvPrime = material.metalnessMapTransform * vec3( uv, 1 );
+					metalness *= texture2D( textures, vec3( uvPrime.xy, material.metalnessMap ) ).b;
+
+				}
+
+				float alphaTest = material.alphaTest;
+				bool useAlphaTest = alphaTest != 0.0;
+				float transmissionFactor = ( 1.0 - metalness ) * transmission;
+				if (
+					transmissionFactor < rand( 9 ) && ! (
+						// material sidedness
+						material.side != 0.0 && surfaceHit.side == material.side
+
+						// alpha test
+						|| useAlphaTest && albedo.a < alphaTest
+
+						// opacity
+						|| material.transparent && ! useAlphaTest && albedo.a < rand( 10 )
+					)
+				) {
+
+					result = true;
+					break;
+
+				}
+
+				if ( surfaceHit.side == 1.0 && isEntering ) {
+
+					// only attenuate by surface color on the way in
+					color *= mix( vec3( 1.0 ), albedo.rgb, transmissionFactor );
+
+				} else if ( surfaceHit.side == - 1.0 ) {
+
+					// attenuate by medium once we hit the opposite side of the model
+					color *= transmissionAttenuation( surfaceHit.dist, material.attenuationColor, material.attenuationDistance );
+
+				}
+
+				bool isTransmissiveRay = dot( ray.direction, surfaceHit.faceNormal * surfaceHit.side ) < 0.0;
+				if ( ( isTransmissiveRay || isEntering ) && transmissiveTraversals > 0 ) {
+
+					i -= sign( transmissiveTraversals );
+					transmissiveTraversals --;
+
+				}
+
+			} else {
+
+				result = false;
+				break;
+
+			}
+
+		}
+
+		// reset the bounce index
+		sobolBounceIndex = originalBounceIndex;
+		return result;
+
+	}
+
+`,Ar=`
+
+	vec3 ndcToRayOrigin( vec2 coord ) {
+
+		vec4 rayOrigin4 = cameraWorldMatrix * invProjectionMatrix * vec4( coord, - 1.0, 1.0 );
+		return rayOrigin4.xyz / rayOrigin4.w;
+	}
+
+	Ray getCameraRay() {
+
+		vec2 ssd = vec2( 1.0 ) / resolution;
+
+		// Jitter the camera ray by finding a uv coordinate at a random sample
+		// around this pixel's UV coordinate for AA
+		vec2 ruv = rand2( 0 );
+		vec2 jitteredUv = vUv + vec2( tentFilter( ruv.x ) * ssd.x, tentFilter( ruv.y ) * ssd.y );
+		Ray ray;
+
+		#if CAMERA_TYPE == 2
+
+			// Equirectangular projection
+			vec4 rayDirection4 = vec4( equirectUvToDirection( jitteredUv ), 0.0 );
+			vec4 rayOrigin4 = vec4( 0.0, 0.0, 0.0, 1.0 );
+
+			rayDirection4 = cameraWorldMatrix * rayDirection4;
+			rayOrigin4 = cameraWorldMatrix * rayOrigin4;
+
+			ray.direction = normalize( rayDirection4.xyz );
+			ray.origin = rayOrigin4.xyz / rayOrigin4.w;
+
+		#else
+
+			// get [- 1, 1] normalized device coordinates
+			vec2 ndc = 2.0 * jitteredUv - vec2( 1.0 );
+			ray.origin = ndcToRayOrigin( ndc );
+
+			#if CAMERA_TYPE == 1
+
+				// Orthographic projection
+				ray.direction = ( cameraWorldMatrix * vec4( 0.0, 0.0, - 1.0, 0.0 ) ).xyz;
+				ray.direction = normalize( ray.direction );
+
+			#else
+
+				// Perspective projection
+				ray.direction = normalize( mat3( cameraWorldMatrix ) * ( invProjectionMatrix * vec4( ndc, 0.0, 1.0 ) ).xyz );
+
+			#endif
+
+		#endif
+
+		#if FEATURE_DOF
+		{
+
+			// depth of field
+			vec3 focalPoint = ray.origin + normalize( ray.direction ) * physicalCamera.focusDistance;
+
+			// get the aperture sample
+			// if blades === 0 then we assume a circle
+			vec3 shapeUVW= rand3( 1 );
+			int blades = physicalCamera.apertureBlades;
+			float anamorphicRatio = physicalCamera.anamorphicRatio;
+			vec2 apertureSample = blades == 0 ? sampleCircle( shapeUVW.xy ) : sampleRegularPolygon( blades, shapeUVW );
+			apertureSample *= physicalCamera.bokehSize * 0.5 * 1e-3;
+
+			// rotate the aperture shape
+			apertureSample =
+				rotateVector( apertureSample, physicalCamera.apertureRotation ) *
+				saturate( vec2( anamorphicRatio, 1.0 / anamorphicRatio ) );
+
+			// create the new ray
+			ray.origin += ( cameraWorldMatrix * vec4( apertureSample, 0.0, 0.0 ) ).xyz;
+			ray.direction = focalPoint - ray.origin;
+
+		}
+		#endif
+
+		ray.direction = normalize( ray.direction );
+
+		return ray;
+
+	}
+
+`,Mr=`
+
+	vec3 directLightContribution( vec3 worldWo, SurfaceRecord surf, RenderState state, vec3 rayOrigin ) {
+
+		vec3 result = vec3( 0.0 );
+
+		// uniformly pick a light or environment map
+		if( lightsDenom != 0.0 && rand( 5 ) < float( lights.count ) / lightsDenom ) {
+
+			// sample a light or environment
+			LightRecord lightRec = randomLightSample( lights.tex, iesProfiles, lights.count, rayOrigin, rand3( 6 ) );
+
+			bool isSampleBelowSurface = ! surf.volumeParticle && dot( surf.faceNormal, lightRec.direction ) < 0.0;
+			if ( isSampleBelowSurface ) {
+
+				lightRec.pdf = 0.0;
+
+			}
+
+			// check if a ray could even reach the light area
+			Ray lightRay;
+			lightRay.origin = rayOrigin;
+			lightRay.direction = lightRec.direction;
+			vec3 attenuatedColor;
+			if (
+				lightRec.pdf > 0.0 &&
+				isDirectionValid( lightRec.direction, surf.normal, surf.faceNormal ) &&
+				! attenuateHit( state, lightRay, lightRec.dist, attenuatedColor )
+			) {
+
+				// get the material pdf
+				vec3 sampleColor;
+				float lightMaterialPdf = bsdfResult( worldWo, lightRec.direction, surf, sampleColor );
+				bool isValidSampleColor = all( greaterThanEqual( sampleColor, vec3( 0.0 ) ) );
+				if ( lightMaterialPdf > 0.0 && isValidSampleColor ) {
+
+					// weight the direct light contribution
+					float lightPdf = lightRec.pdf / lightsDenom;
+					float misWeight = lightRec.type == SPOT_LIGHT_TYPE || lightRec.type == DIR_LIGHT_TYPE || lightRec.type == POINT_LIGHT_TYPE ? 1.0 : misHeuristic( lightPdf, lightMaterialPdf );
+					result = attenuatedColor * lightRec.emission * state.throughputColor * sampleColor * misWeight / lightPdf;
+
+				}
+
+			}
+
+		} else if ( envMapInfo.totalSum != 0.0 && environmentIntensity != 0.0 ) {
+
+			// find a sample in the environment map to include in the contribution
+			vec3 envColor, envDirection;
+			float envPdf = sampleEquirectProbability( rand2( 7 ), envColor, envDirection );
+			envDirection = invEnvRotation3x3 * envDirection;
+
+			// this env sampling is not set up for transmissive sampling and yields overly bright
+			// results so we ignore the sample in this case.
+			// TODO: this should be improved but how? The env samples could traverse a few layers?
+			bool isSampleBelowSurface = ! surf.volumeParticle && dot( surf.faceNormal, envDirection ) < 0.0;
+			if ( isSampleBelowSurface ) {
+
+				envPdf = 0.0;
+
+			}
+
+			// check if a ray could even reach the surface
+			Ray envRay;
+			envRay.origin = rayOrigin;
+			envRay.direction = envDirection;
+			vec3 attenuatedColor;
+			if (
+				envPdf > 0.0 &&
+				isDirectionValid( envDirection, surf.normal, surf.faceNormal ) &&
+				! attenuateHit( state, envRay, INFINITY, attenuatedColor )
+			) {
+
+				// get the material pdf
+				vec3 sampleColor;
+				float envMaterialPdf = bsdfResult( worldWo, envDirection, surf, sampleColor );
+				bool isValidSampleColor = all( greaterThanEqual( sampleColor, vec3( 0.0 ) ) );
+				if ( envMaterialPdf > 0.0 && isValidSampleColor ) {
+
+					// weight the direct light contribution
+					envPdf /= lightsDenom;
+					float misWeight = misHeuristic( envPdf, envMaterialPdf );
+					result = attenuatedColor * environmentIntensity * envColor * state.throughputColor * sampleColor * misWeight / envPdf;
+
+				}
+
+			}
+
+		}
+
+		// Function changed to have a single return statement to potentially help with crashes on Mac OS.
+		// See issue #470
+		return result;
+
+	}
+
+`,Cr=`
+
+	#define SKIP_SURFACE 0
+	#define HIT_SURFACE 1
+	int getSurfaceRecord(
+		Material material, SurfaceHit surfaceHit, sampler2DArray attributesArray,
+		float accumulatedRoughness,
+		inout SurfaceRecord surf
+	) {
+
+		if ( material.fogVolume ) {
+
+			vec3 normal = vec3( 0, 0, 1 );
+
+			SurfaceRecord fogSurface;
+			fogSurface.volumeParticle = true;
+			fogSurface.color = material.color;
+			fogSurface.emission = material.emissiveIntensity * material.emissive;
+			fogSurface.normal = normal;
+			fogSurface.faceNormal = normal;
+			fogSurface.clearcoatNormal = normal;
+
+			surf = fogSurface;
+			return HIT_SURFACE;
+
+		}
+
+		// uv coord for textures
+		vec2 uv = textureSampleBarycoord( attributesArray, ATTR_UV, surfaceHit.barycoord, surfaceHit.faceIndices.xyz ).xy;
+		vec4 vertexColor = textureSampleBarycoord( attributesArray, ATTR_COLOR, surfaceHit.barycoord, surfaceHit.faceIndices.xyz );
+
+		// albedo
+		vec4 albedo = vec4( material.color, material.opacity );
+		if ( material.map != - 1 ) {
+
+			vec3 uvPrime = material.mapTransform * vec3( uv, 1 );
+			albedo *= texture2D( textures, vec3( uvPrime.xy, material.map ) );
+
+		}
+
+		if ( material.vertexColors ) {
+
+			albedo *= vertexColor;
+
+		}
+
+		// alphaMap
+		if ( material.alphaMap != - 1 ) {
+
+			albedo.a *= texture2D( textures, vec3( uv, material.alphaMap ) ).x;
+
+		}
+
+		// possibly skip this sample if it's transparent, alpha test is enabled, or we hit the wrong material side
+		// and it's single sided.
+		// - alpha test is disabled when it === 0
+		// - the material sidedness test is complicated because we want light to pass through the back side but still
+		// be able to see the front side. This boolean checks if the side we hit is the front side on the first ray
+		// and we're rendering the other then we skip it. Do the opposite on subsequent bounces to get incoming light.
+		float alphaTest = material.alphaTest;
+		bool useAlphaTest = alphaTest != 0.0;
+		if (
+			// material sidedness
+			material.side != 0.0 && surfaceHit.side != material.side
+
+			// alpha test
+			|| useAlphaTest && albedo.a < alphaTest
+
+			// opacity
+			|| material.transparent && ! useAlphaTest && albedo.a < rand( 3 )
+		) {
+
+			return SKIP_SURFACE;
+
+		}
+
+		// fetch the interpolated smooth normal
+		vec3 normal = normalize( textureSampleBarycoord(
+			attributesArray,
+			ATTR_NORMAL,
+			surfaceHit.barycoord,
+			surfaceHit.faceIndices.xyz
+		).xyz );
+
+		// roughness
+		float roughness = material.roughness;
+		if ( material.roughnessMap != - 1 ) {
+
+			vec3 uvPrime = material.roughnessMapTransform * vec3( uv, 1 );
+			roughness *= texture2D( textures, vec3( uvPrime.xy, material.roughnessMap ) ).g;
+
+		}
+
+		// metalness
+		float metalness = material.metalness;
+		if ( material.metalnessMap != - 1 ) {
+
+			vec3 uvPrime = material.metalnessMapTransform * vec3( uv, 1 );
+			metalness *= texture2D( textures, vec3( uvPrime.xy, material.metalnessMap ) ).b;
+
+		}
+
+		// emission
+		vec3 emission = material.emissiveIntensity * material.emissive;
+		if ( material.emissiveMap != - 1 ) {
+
+			vec3 uvPrime = material.emissiveMapTransform * vec3( uv, 1 );
+			emission *= texture2D( textures, vec3( uvPrime.xy, material.emissiveMap ) ).xyz;
+
+		}
+
+		// transmission
+		float transmission = material.transmission;
+		if ( material.transmissionMap != - 1 ) {
+
+			vec3 uvPrime = material.transmissionMapTransform * vec3( uv, 1 );
+			transmission *= texture2D( textures, vec3( uvPrime.xy, material.transmissionMap ) ).r;
+
+		}
+
+		// normal
+		if ( material.flatShading ) {
+
+			// if we're rendering a flat shaded object then use the face normals - the face normal
+			// is provided based on the side the ray hits the mesh so flip it to align with the
+			// interpolated vertex normals.
+			normal = surfaceHit.faceNormal * surfaceHit.side;
+
+		}
+
+		vec3 baseNormal = normal;
+		if ( material.normalMap != - 1 ) {
+
+			vec4 tangentSample = textureSampleBarycoord(
+				attributesArray,
+				ATTR_TANGENT,
+				surfaceHit.barycoord,
+				surfaceHit.faceIndices.xyz
+			);
+
+			// some provided tangents can be malformed (0, 0, 0) causing the normal to be degenerate
+			// resulting in NaNs and slow path tracing.
+			if ( length( tangentSample.xyz ) > 0.0 ) {
+
+				vec3 tangent = normalize( tangentSample.xyz );
+				vec3 bitangent = normalize( cross( normal, tangent ) * tangentSample.w );
+				mat3 vTBN = mat3( tangent, bitangent, normal );
+
+				vec3 uvPrime = material.normalMapTransform * vec3( uv, 1 );
+				vec3 texNormal = texture2D( textures, vec3( uvPrime.xy, material.normalMap ) ).xyz * 2.0 - 1.0;
+				texNormal.xy *= material.normalScale;
+				normal = vTBN * texNormal;
+
+			}
+
+		}
+
+		normal *= surfaceHit.side;
+
+		// clearcoat
+		float clearcoat = material.clearcoat;
+		if ( material.clearcoatMap != - 1 ) {
+
+			vec3 uvPrime = material.clearcoatMapTransform * vec3( uv, 1 );
+			clearcoat *= texture2D( textures, vec3( uvPrime.xy, material.clearcoatMap ) ).r;
+
+		}
+
+		// clearcoatRoughness
+		float clearcoatRoughness = material.clearcoatRoughness;
+		if ( material.clearcoatRoughnessMap != - 1 ) {
+
+			vec3 uvPrime = material.clearcoatRoughnessMapTransform * vec3( uv, 1 );
+			clearcoatRoughness *= texture2D( textures, vec3( uvPrime.xy, material.clearcoatRoughnessMap ) ).g;
+
+		}
+
+		// clearcoatNormal
+		vec3 clearcoatNormal = baseNormal;
+		if ( material.clearcoatNormalMap != - 1 ) {
+
+			vec4 tangentSample = textureSampleBarycoord(
+				attributesArray,
+				ATTR_TANGENT,
+				surfaceHit.barycoord,
+				surfaceHit.faceIndices.xyz
+			);
+
+			// some provided tangents can be malformed (0, 0, 0) causing the normal to be degenerate
+			// resulting in NaNs and slow path tracing.
+			if ( length( tangentSample.xyz ) > 0.0 ) {
+
+				vec3 tangent = normalize( tangentSample.xyz );
+				vec3 bitangent = normalize( cross( clearcoatNormal, tangent ) * tangentSample.w );
+				mat3 vTBN = mat3( tangent, bitangent, clearcoatNormal );
+
+				vec3 uvPrime = material.clearcoatNormalMapTransform * vec3( uv, 1 );
+				vec3 texNormal = texture2D( textures, vec3( uvPrime.xy, material.clearcoatNormalMap ) ).xyz * 2.0 - 1.0;
+				texNormal.xy *= material.clearcoatNormalScale;
+				clearcoatNormal = vTBN * texNormal;
+
+			}
+
+		}
+
+		clearcoatNormal *= surfaceHit.side;
+
+		// sheenColor
+		vec3 sheenColor = material.sheenColor;
+		if ( material.sheenColorMap != - 1 ) {
+
+			vec3 uvPrime = material.sheenColorMapTransform * vec3( uv, 1 );
+			sheenColor *= texture2D( textures, vec3( uvPrime.xy, material.sheenColorMap ) ).rgb;
+
+		}
+
+		// sheenRoughness
+		float sheenRoughness = material.sheenRoughness;
+		if ( material.sheenRoughnessMap != - 1 ) {
+
+			vec3 uvPrime = material.sheenRoughnessMapTransform * vec3( uv, 1 );
+			sheenRoughness *= texture2D( textures, vec3( uvPrime.xy, material.sheenRoughnessMap ) ).a;
+
+		}
+
+		// iridescence
+		float iridescence = material.iridescence;
+		if ( material.iridescenceMap != - 1 ) {
+
+			vec3 uvPrime = material.iridescenceMapTransform * vec3( uv, 1 );
+			iridescence *= texture2D( textures, vec3( uvPrime.xy, material.iridescenceMap ) ).r;
+
+		}
+
+		// iridescence thickness
+		float iridescenceThickness = material.iridescenceThicknessMaximum;
+		if ( material.iridescenceThicknessMap != - 1 ) {
+
+			vec3 uvPrime = material.iridescenceThicknessMapTransform * vec3( uv, 1 );
+			float iridescenceThicknessSampled = texture2D( textures, vec3( uvPrime.xy, material.iridescenceThicknessMap ) ).g;
+			iridescenceThickness = mix( material.iridescenceThicknessMinimum, material.iridescenceThicknessMaximum, iridescenceThicknessSampled );
+
+		}
+
+		iridescence = iridescenceThickness == 0.0 ? 0.0 : iridescence;
+
+		// specular color
+		vec3 specularColor = material.specularColor;
+		if ( material.specularColorMap != - 1 ) {
+
+			vec3 uvPrime = material.specularColorMapTransform * vec3( uv, 1 );
+			specularColor *= texture2D( textures, vec3( uvPrime.xy, material.specularColorMap ) ).rgb;
+
+		}
+
+		// specular intensity
+		float specularIntensity = material.specularIntensity;
+		if ( material.specularIntensityMap != - 1 ) {
+
+			vec3 uvPrime = material.specularIntensityMapTransform * vec3( uv, 1 );
+			specularIntensity *= texture2D( textures, vec3( uvPrime.xy, material.specularIntensityMap ) ).a;
+
+		}
+
+		surf.volumeParticle = false;
+
+		surf.faceNormal = surfaceHit.faceNormal;
+		surf.normal = normal;
+
+		surf.metalness = metalness;
+		surf.color = albedo.rgb;
+		surf.emission = emission;
+
+		surf.ior = material.ior;
+		surf.transmission = transmission;
+		surf.thinFilm = material.thinFilm;
+		surf.attenuationColor = material.attenuationColor;
+		surf.attenuationDistance = material.attenuationDistance;
+
+		surf.clearcoatNormal = clearcoatNormal;
+		surf.clearcoat = clearcoat;
+
+		surf.sheen = material.sheen;
+		surf.sheenColor = sheenColor;
+
+		surf.iridescence = iridescence;
+		surf.iridescenceIor = material.iridescenceIor;
+		surf.iridescenceThickness = iridescenceThickness;
+
+		surf.specularColor = specularColor;
+		surf.specularIntensity = specularIntensity;
+
+		// apply perceptual roughness factor from gltf. sheen perceptual roughness is
+		// applied by its brdf function
+		// https://registry.khronos.org/glTF/specs/2.0/glTF-2.0.html#microfacet-surfaces
+		surf.roughness = roughness * roughness;
+		surf.clearcoatRoughness = clearcoatRoughness * clearcoatRoughness;
+		surf.sheenRoughness = sheenRoughness;
+
+		// frontFace is used to determine transmissive properties and PDF. If no transmission is used
+		// then we can just always assume this is a front face.
+		surf.frontFace = surfaceHit.side == 1.0 || transmission == 0.0;
+		surf.eta = material.thinFilm || surf.frontFace ? 1.0 / material.ior : material.ior;
+		surf.f0 = iorRatioToF0( surf.eta );
+
+		// Compute the filtered roughness value to use during specular reflection computations.
+		// The accumulated roughness value is scaled by a user setting and a "magic value" of 5.0.
+		// If we're exiting something transmissive then scale the factor down significantly so we can retain
+		// sharp internal reflections
+		surf.filteredRoughness = applyFilteredGlossy( surf.roughness, accumulatedRoughness );
+		surf.filteredClearcoatRoughness = applyFilteredGlossy( surf.clearcoatRoughness, accumulatedRoughness );
+
+		// get the normal frames
+		surf.normalBasis = getBasisFromNormal( surf.normal );
+		surf.normalInvBasis = inverse( surf.normalBasis );
+
+		surf.clearcoatBasis = getBasisFromNormal( surf.clearcoatNormal );
+		surf.clearcoatInvBasis = inverse( surf.clearcoatBasis );
+
+		return HIT_SURFACE;
+
+	}
+`,Dr=`
+
+	struct Ray {
+
+		vec3 origin;
+		vec3 direction;
+
+	};
+
+	struct SurfaceHit {
+
+		uvec4 faceIndices;
+		vec3 barycoord;
+		vec3 faceNormal;
+		float side;
+		float dist;
+
+	};
+
+	struct RenderState {
+
+		bool firstRay;
+		bool transmissiveRay;
+		bool isShadowRay;
+		float accumulatedRoughness;
+		int transmissiveTraversals;
+		int traversals;
+		uint depth;
+		vec3 throughputColor;
+		Material fogMaterial;
+
+	};
+
+	RenderState initRenderState() {
+
+		RenderState result;
+		result.firstRay = true;
+		result.transmissiveRay = true;
+		result.isShadowRay = false;
+		result.accumulatedRoughness = 0.0;
+		result.transmissiveTraversals = 0;
+		result.traversals = 0;
+		result.throughputColor = vec3( 1.0 );
+		result.depth = 0u;
+		result.fogMaterial.fogVolume = false;
+		return result;
+
+	}
+
+`,Pr=`
+
+	#define NO_HIT 0
+	#define SURFACE_HIT 1
+	#define LIGHT_HIT 2
+	#define FOG_HIT 3
+
+	// Passing the global variable 'lights' into this function caused shader program errors.
+	// So global variables like 'lights' and 'bvh' were moved out of the function parameters.
+	// For more information, refer to: https://github.com/gkjohnson/three-gpu-pathtracer/pull/457
+	int traceScene(
+		Ray ray, Material fogMaterial, inout SurfaceHit surfaceHit
+	) {
+
+		int result = NO_HIT;
+		bool hit = bvhIntersectFirstHit( bvh, ray.origin, ray.direction, surfaceHit.faceIndices, surfaceHit.faceNormal, surfaceHit.barycoord, surfaceHit.side, surfaceHit.dist );
+
+		#if FEATURE_FOG
+
+		if ( fogMaterial.fogVolume ) {
+
+			// offset the distance so we don't run into issues with particles on the same surface
+			// as other objects
+			float particleDist = intersectFogVolume( fogMaterial, rand( 1 ) );
+			if ( particleDist + RAY_OFFSET < surfaceHit.dist ) {
+
+				surfaceHit.side = 1.0;
+				surfaceHit.faceNormal = normalize( - ray.direction );
+				surfaceHit.dist = particleDist;
+				return FOG_HIT;
+
+			}
+
+		}
+
+		#endif
+
+		if ( hit ) {
+
+			result = SURFACE_HIT;
+
+		}
+
+		return result;
+
+	}
+
+`;class Er extends oe{onBeforeRender(){this.setDefine("FEATURE_DOF",this.physicalCamera.bokehSize===0?0:1),this.setDefine("FEATURE_BACKGROUND_MAP",this.backgroundMap?1:0),this.setDefine("FEATURE_FOG",this.materials.features.isUsed("FOG")?1:0)}constructor(e){super({transparent:!0,depthWrite:!1,defines:{FEATURE_MIS:1,FEATURE_RUSSIAN_ROULETTE:1,FEATURE_DOF:1,FEATURE_BACKGROUND_MAP:0,FEATURE_FOG:1,RANDOM_TYPE:2,CAMERA_TYPE:0,DEBUG_MODE:0,ATTR_NORMAL:0,ATTR_TANGENT:1,ATTR_UV:2,ATTR_COLOR:3},uniforms:{resolution:{value:new l.Vector2},opacity:{value:1},bounces:{value:10},transmissiveBounces:{value:10},filterGlossyFactor:{value:0},physicalCamera:{value:new Mt},cameraWorldMatrix:{value:new l.Matrix4},invProjectionMatrix:{value:new l.Matrix4},bvh:{value:new Xe},attributesArray:{value:new Lt},materialIndexAttribute:{value:new Be},materials:{value:new jt},textures:{value:new Fe().texture},lights:{value:new Ht},iesProfiles:{value:new Fe(360,180,{type:l.HalfFloatType,wrapS:l.ClampToEdgeWrapping,wrapT:l.ClampToEdgeWrapping}).texture},environmentIntensity:{value:1},environmentRotation:{value:new l.Matrix4},envMapInfo:{value:new Pt},backgroundBlur:{value:0},backgroundMap:{value:null},backgroundAlpha:{value:1},backgroundIntensity:{value:1},backgroundRotation:{value:new l.Matrix4},seed:{value:0},sobolTexture:{value:null},stratifiedTexture:{value:new tr},stratifiedOffsetTexture:{value:new nr(64,1)}},vertexShader:`
+
+				varying vec2 vUv;
+				void main() {
+
+					vec4 mvPosition = vec4( position, 1.0 );
+					mvPosition = modelViewMatrix * mvPosition;
+					gl_Position = projectionMatrix * mvPosition;
+
+					vUv = uv;
+
+				}
+
+			`,fragmentShader:`
+				#define RAY_OFFSET 1e-4
+				#define INFINITY 1e20
+
+				precision highp isampler2D;
+				precision highp usampler2D;
+				precision highp sampler2DArray;
+				vec4 envMapTexelToLinear( vec4 a ) { return a; }
+				#include <common>
+
+				// bvh intersection
+				${Ze}
+				${et}
+				${Je}
+
+				// uniform structs
+				${lr}
+				${ur}
+				${cr}
+				${fr}
+				${hr}
+
+				// random
+				#if RANDOM_TYPE == 2 	// Stratified List
+
+					${br}
+
+				#elif RANDOM_TYPE == 1 	// Sobol
+
+					${Me}
+					${Le}
+					${St}
+
+					#define rand(v) sobol(v)
+					#define rand2(v) sobol2(v)
+					#define rand3(v) sobol3(v)
+					#define rand4(v) sobol4(v)
+
+				#else 					// PCG
+
+				${Me}
+
+					// Using the sobol functions seems to break the the compiler on MacOS
+					// - specifically the "sobolReverseBits" function.
+					uint sobolPixelIndex = 0u;
+					uint sobolPathIndex = 0u;
+					uint sobolBounceIndex = 0u;
+
+					#define rand(v) pcgRand()
+					#define rand2(v) pcgRand2()
+					#define rand3(v) pcgRand3()
+					#define rand4(v) pcgRand4()
+
+				#endif
+
+				// common
+				${yr}
+				${gr}
+				${Ge}
+				${vr}
+				${xr}
+
+				// environment
+				uniform EquirectHdrInfo envMapInfo;
+				uniform mat4 environmentRotation;
+				uniform float environmentIntensity;
+
+				// lighting
+				uniform sampler2DArray iesProfiles;
+				uniform LightsInfo lights;
+
+				// background
+				uniform float backgroundBlur;
+				uniform float backgroundAlpha;
+				#if FEATURE_BACKGROUND_MAP
+
+				uniform sampler2D backgroundMap;
+				uniform mat4 backgroundRotation;
+				uniform float backgroundIntensity;
+
+				#endif
+
+				// camera
+				uniform mat4 cameraWorldMatrix;
+				uniform mat4 invProjectionMatrix;
+				#if FEATURE_DOF
+
+				uniform PhysicalCamera physicalCamera;
+
+				#endif
+
+				// geometry
+				uniform sampler2DArray attributesArray;
+				uniform usampler2D materialIndexAttribute;
+				uniform sampler2D materials;
+				uniform sampler2DArray textures;
+				uniform BVH bvh;
+
+				// path tracer
+				uniform int bounces;
+				uniform int transmissiveBounces;
+				uniform float filterGlossyFactor;
+				uniform int seed;
+
+				// image
+				uniform vec2 resolution;
+				uniform float opacity;
+
+				varying vec2 vUv;
+
+				// globals
+				mat3 envRotation3x3;
+				mat3 invEnvRotation3x3;
+				float lightsDenom;
+
+				// sampling
+				${pr}
+				${dr}
+				${mr}
+
+				${_r}
+				${Rr}
+				${Sr}
+				${Ir}
+				${wr}
+				${Tr}
+
+				float applyFilteredGlossy( float roughness, float accumulatedRoughness ) {
+
+					return clamp(
+						max(
+							roughness,
+							accumulatedRoughness * filterGlossyFactor * 5.0 ),
+						0.0,
+						1.0
+					);
+
+				}
+
+				vec3 sampleBackground( vec3 direction, vec2 uv ) {
+
+					vec3 sampleDir = sampleHemisphere( direction, uv ) * 0.5 * backgroundBlur;
+
+					#if FEATURE_BACKGROUND_MAP
+
+					sampleDir = normalize( mat3( backgroundRotation ) * direction + sampleDir );
+					return backgroundIntensity * sampleEquirectColor( backgroundMap, sampleDir );
+
+					#else
+
+					sampleDir = normalize( envRotation3x3 * direction + sampleDir );
+					return environmentIntensity * sampleEquirectColor( envMapInfo.map, sampleDir );
+
+					#endif
+
+				}
+
+				${Dr}
+				${Ar}
+				${Pr}
+				${Fr}
+				${Mr}
+				${Cr}
+
+				void main() {
+
+					// init
+					rng_initialize( gl_FragCoord.xy, seed );
+					sobolPixelIndex = ( uint( gl_FragCoord.x ) << 16 ) | uint( gl_FragCoord.y );
+					sobolPathIndex = uint( seed );
+
+					// get camera ray
+					Ray ray = getCameraRay();
+
+					// inverse environment rotation
+					envRotation3x3 = mat3( environmentRotation );
+					invEnvRotation3x3 = inverse( envRotation3x3 );
+					lightsDenom =
+						( environmentIntensity == 0.0 || envMapInfo.totalSum == 0.0 ) && lights.count != 0u ?
+							float( lights.count ) :
+							float( lights.count + 1u );
+
+					// final color
+					gl_FragColor = vec4( 0, 0, 0, 1 );
+
+					// surface results
+					SurfaceHit surfaceHit;
+					ScatterRecord scatterRec;
+
+					// path tracing state
+					RenderState state = initRenderState();
+					state.transmissiveTraversals = transmissiveBounces;
+					#if FEATURE_FOG
+
+					state.fogMaterial.fogVolume = bvhIntersectFogVolumeHit(
+						ray.origin, - ray.direction,
+						materialIndexAttribute, materials,
+						state.fogMaterial
+					);
+
+					#endif
+
+					for ( int i = 0; i < bounces; i ++ ) {
+
+						sobolBounceIndex ++;
+
+						state.depth ++;
+						state.traversals = bounces - i;
+						state.firstRay = i == 0 && state.transmissiveTraversals == transmissiveBounces;
+
+						int hitType = traceScene( ray, state.fogMaterial, surfaceHit );
+
+						// check if we intersect any lights and accumulate the light contribution
+						// TODO: we can add support for light surface rendering in the else condition if we
+						// add the ability to toggle visibility of the the light
+						if ( ! state.firstRay && ! state.transmissiveRay ) {
+
+							LightRecord lightRec;
+							float lightDist = hitType == NO_HIT ? INFINITY : surfaceHit.dist;
+							for ( uint i = 0u; i < lights.count; i ++ ) {
+
+								if (
+									intersectLightAtIndex( lights.tex, ray.origin, ray.direction, i, lightRec ) &&
+									lightRec.dist < lightDist
+								) {
+
+									#if FEATURE_MIS
+
+									// weight the contribution
+									// NOTE: Only area lights are supported for forward sampling and can be hit
+									float misWeight = misHeuristic( scatterRec.pdf, lightRec.pdf / lightsDenom );
+									gl_FragColor.rgb += lightRec.emission * state.throughputColor * misWeight;
+
+									#else
+
+									gl_FragColor.rgb += lightRec.emission * state.throughputColor;
+
+									#endif
+
+								}
+
+							}
+
+						}
+
+						if ( hitType == NO_HIT ) {
+
+							if ( state.firstRay || state.transmissiveRay ) {
+
+								gl_FragColor.rgb += sampleBackground( ray.direction, rand2( 2 ) ) * state.throughputColor;
+								gl_FragColor.a = backgroundAlpha;
+
+							} else {
+
+								#if FEATURE_MIS
+
+								// get the PDF of the hit envmap point
+								vec3 envColor;
+								float envPdf = sampleEquirect( envRotation3x3 * ray.direction, envColor );
+								envPdf /= lightsDenom;
+
+								// and weight the contribution
+								float misWeight = misHeuristic( scatterRec.pdf, envPdf );
+								gl_FragColor.rgb += environmentIntensity * envColor * state.throughputColor * misWeight;
+
+								#else
+
+								gl_FragColor.rgb +=
+									environmentIntensity *
+									sampleEquirectColor( envMapInfo.map, envRotation3x3 * ray.direction ) *
+									state.throughputColor;
+
+								#endif
+
+							}
+							break;
+
+						}
+
+						uint materialIndex = uTexelFetch1D( materialIndexAttribute, surfaceHit.faceIndices.x ).r;
+						Material material = readMaterialInfo( materials, materialIndex );
+
+						#if FEATURE_FOG
+
+						if ( hitType == FOG_HIT ) {
+
+							material = state.fogMaterial;
+							state.accumulatedRoughness += 0.2;
+
+						} else if ( material.fogVolume ) {
+
+							state.fogMaterial = material;
+							state.fogMaterial.fogVolume = surfaceHit.side == 1.0;
+
+							ray.origin = stepRayOrigin( ray.origin, ray.direction, - surfaceHit.faceNormal, surfaceHit.dist );
+
+							i -= sign( state.transmissiveTraversals );
+							state.transmissiveTraversals -= sign( state.transmissiveTraversals );
+							continue;
+
+						}
+
+						#endif
+
+						// early out if this is a matte material
+						if ( material.matte && state.firstRay ) {
+
+							gl_FragColor = vec4( 0.0 );
+							break;
+
+						}
+
+						// if we've determined that this is a shadow ray and we've hit an item with no shadow casting
+						// then skip it
+						if ( ! material.castShadow && state.isShadowRay ) {
+
+							ray.origin = stepRayOrigin( ray.origin, ray.direction, - surfaceHit.faceNormal, surfaceHit.dist );
+							continue;
+
+						}
+
+						SurfaceRecord surf;
+						if (
+							getSurfaceRecord(
+								material, surfaceHit, attributesArray, state.accumulatedRoughness,
+								surf
+							) == SKIP_SURFACE
+						) {
+
+							// only allow a limited number of transparency discards otherwise we could
+							// crash the context with too long a loop.
+							i -= sign( state.transmissiveTraversals );
+							state.transmissiveTraversals -= sign( state.transmissiveTraversals );
+
+							ray.origin = stepRayOrigin( ray.origin, ray.direction, - surfaceHit.faceNormal, surfaceHit.dist );
+							continue;
+
+						}
+
+						scatterRec = bsdfSample( - ray.direction, surf );
+						state.isShadowRay = scatterRec.specularPdf < rand( 4 );
+
+						bool isBelowSurface = ! surf.volumeParticle && dot( scatterRec.direction, surf.faceNormal ) < 0.0;
+						vec3 hitPoint = stepRayOrigin( ray.origin, ray.direction, isBelowSurface ? - surf.faceNormal : surf.faceNormal, surfaceHit.dist );
+
+						// next event estimation
+						#if FEATURE_MIS
+
+						gl_FragColor.rgb += directLightContribution( - ray.direction, surf, state, hitPoint );
+
+						#endif
+
+						// accumulate a roughness value to offset diffuse, specular, diffuse rays that have high contribution
+						// to a single pixel resulting in fireflies
+						// TODO: handle transmissive surfaces
+						if ( ! surf.volumeParticle && ! isBelowSurface ) {
+
+							// determine if this is a rough normal or not by checking how far off straight up it is
+							vec3 halfVector = normalize( - ray.direction + scatterRec.direction );
+							state.accumulatedRoughness += max(
+								sin( acosApprox( dot( halfVector, surf.normal ) ) ),
+								sin( acosApprox( dot( halfVector, surf.clearcoatNormal ) ) )
+							);
+
+							state.transmissiveRay = false;
+
+						}
+
+						// accumulate emissive color
+						gl_FragColor.rgb += ( surf.emission * state.throughputColor );
+
+						// skip the sample if our PDF or ray is impossible
+						if ( scatterRec.pdf <= 0.0 || ! isDirectionValid( scatterRec.direction, surf.normal, surf.faceNormal ) ) {
+
+							break;
+
+						}
+
+						// if we're bouncing around the inside a transmissive material then decrement
+						// perform this separate from a bounce
+						bool isTransmissiveRay = ! surf.volumeParticle && dot( scatterRec.direction, surf.faceNormal * surfaceHit.side ) < 0.0;
+						if ( ( isTransmissiveRay || isBelowSurface ) && state.transmissiveTraversals > 0 ) {
+
+							state.transmissiveTraversals --;
+							i --;
+
+						}
+
+						//
+
+						// handle throughput color transformation
+						// attenuate the throughput color by the medium color
+						if ( ! surf.frontFace ) {
+
+							state.throughputColor *= transmissionAttenuation( surfaceHit.dist, surf.attenuationColor, surf.attenuationDistance );
+
+						}
+
+						#if FEATURE_RUSSIAN_ROULETTE
+
+						// russian roulette path termination
+						// https://www.arnoldrenderer.com/research/physically_based_shader_design_in_arnold.pdf
+						uint minBounces = 3u;
+						float depthProb = float( state.depth < minBounces );
+
+						float rrProb = luminance( state.throughputColor * scatterRec.color / scatterRec.pdf );
+						rrProb /= luminance( state.throughputColor );
+						rrProb = sqrt( rrProb );
+						rrProb = max( rrProb, depthProb );
+						rrProb = min( rrProb, 1.0 );
+						if ( rand( 8 ) > rrProb ) {
+
+							break;
+
+						}
+
+						// perform sample clamping here to avoid bright pixels
+						state.throughputColor *= min( 1.0 / rrProb, 20.0 );
+
+						#endif
+
+						// adjust the throughput and discard and exit if we find discard the sample if there are any NaNs
+						state.throughputColor *= scatterRec.color / scatterRec.pdf;
+						if ( any( isnan( state.throughputColor ) ) || any( isinf( state.throughputColor ) ) ) {
+
+							break;
+
+						}
+
+						//
+
+						// prepare for next ray
+						ray.direction = scatterRec.direction;
+						ray.origin = hitPoint;
+
+					}
+
+					gl_FragColor.a *= opacity;
+
+					#if DEBUG_MODE == 1
+
+					// output the number of rays checked in the path and number of
+					// transmissive rays encountered.
+					gl_FragColor.rgb = vec3(
+						float( state.depth ),
+						transmissiveBounces - state.transmissiveTraversals,
+						0.0
+					);
+					gl_FragColor.a = 1.0;
+
+					#endif
+
+				}
+
+			`}),this.setValues(e)}}function*kr(){const{_renderer:s,_fsQuad:e,_blendQuad:t,_primaryTarget:r,_blendTargets:i,_sobolTarget:o,_subframe:a,alpha:c,material:h}=this,d=new l.Vector4,m=new l.Vector4,f=t.material;let[n,x]=i;for(;;){c?(f.opacity=this._opacityFactor/(this.samples+1),h.blending=l.NoBlending,h.opacity=1):(h.opacity=this._opacityFactor/(this.samples+1),h.blending=l.NormalBlending);const[p,v,u,g]=a,y=r.width,b=r.height;h.resolution.set(y*u,b*g),h.sobolTexture=o.texture,h.stratifiedTexture.init(20,h.bounces+h.transmissiveBounces+5),h.stratifiedTexture.next(),h.seed++;const w=this.tiles.x||1,R=this.tiles.y||1,_=w*R,F=Math.ceil(y*u),Y=Math.ceil(b*g),le=Math.floor(p*y),ce=Math.floor(v*b),X=Math.ceil(F/w),Q=Math.ceil(Y/R);for(let W=0;W<R;W++)for(let L=0;L<w;L++){const ue=s.getRenderTarget(),Ve=s.autoClear,qe=s.getScissorTest();s.getScissor(d),s.getViewport(m);let K=L,fe=W;if(!this.stableTiles){const Z=this._currentTile%(w*R);K=Z%w,fe=~~(Z/w),this._currentTile=Z+1}const he=R-fe-1;r.scissor.set(le+K*X,ce+he*Q,Math.min(X,F-K*X),Math.min(Q,Y-he*Q)),r.viewport.set(le,ce,F,Y),s.setRenderTarget(r),s.setScissorTest(!0),s.autoClear=!1,e.render(s),s.setViewport(m),s.setScissor(d),s.setScissorTest(qe),s.setRenderTarget(ue),s.autoClear=Ve,c&&(f.target1=n.texture,f.target2=r.texture,s.setRenderTarget(x),t.render(s),s.setRenderTarget(ue)),this.samples+=1/_,L===w-1&&W===R-1&&(this.samples=Math.round(this.samples)),yield}[n,x]=[x,n]}}const Ce=new l.Color;class De{get material(){return this._fsQuad.material}set material(e){this._fsQuad.material.removeEventListener("recompilation",this._compileFunction),e.addEventListener("recompilation",this._compileFunction),this._fsQuad.material=e}get target(){return this._alpha?this._blendTargets[1]:this._primaryTarget}set alpha(e){this._alpha!==e&&(e||(this._blendTargets[0].dispose(),this._blendTargets[1].dispose()),this._alpha=e,this.reset())}get alpha(){return this._alpha}get isCompiling(){return!!this._compilePromise}constructor(e){this.camera=null,this.tiles=new l.Vector2(3,3),this.stableNoise=!1,this.stableTiles=!0,this.samples=0,this._subframe=new l.Vector4(0,0,1,1),this._opacityFactor=1,this._renderer=e,this._alpha=!1,this._fsQuad=new z(new Er),this._blendQuad=new z(new Rt),this._task=null,this._currentTile=0,this._compilePromise=null,this._sobolTarget=new Ft().generate(e),this._primaryTarget=new l.WebGLRenderTarget(1,1,{format:l.RGBAFormat,type:l.FloatType,magFilter:l.NearestFilter,minFilter:l.NearestFilter}),this._blendTargets=[new l.WebGLRenderTarget(1,1,{format:l.RGBAFormat,type:l.FloatType,magFilter:l.NearestFilter,minFilter:l.NearestFilter}),new l.WebGLRenderTarget(1,1,{format:l.RGBAFormat,type:l.FloatType,magFilter:l.NearestFilter,minFilter:l.NearestFilter})],this._compileFunction=()=>{const t=this.compileMaterial(this._fsQuad._mesh);t.then(()=>{this._compilePromise===t&&(this._compilePromise=null)}),this._compilePromise=t},this.material.addEventListener("recompilation",this._compileFunction)}compileMaterial(){return this._renderer.compileAsync(this._fsQuad._mesh)}setCamera(e){const{material:t}=this;t.cameraWorldMatrix.copy(e.matrixWorld),t.invProjectionMatrix.copy(e.projectionMatrixInverse),t.physicalCamera.updateFrom(e);let r=0;e.projectionMatrix.elements[15]>0&&(r=1),e.isEquirectCamera&&(r=2),t.setDefine("CAMERA_TYPE",r),this.camera=e}setSize(e,t){e=Math.ceil(e),t=Math.ceil(t),!(this._primaryTarget.width===e&&this._primaryTarget.height===t)&&(this._primaryTarget.setSize(e,t),this._blendTargets[0].setSize(e,t),this._blendTargets[1].setSize(e,t),this.reset())}getSize(e){e.x=this._primaryTarget.width,e.y=this._primaryTarget.height}dispose(){this._primaryTarget.dispose(),this._blendTargets[0].dispose(),this._blendTargets[1].dispose(),this._sobolTarget.dispose(),this._fsQuad.dispose(),this._blendQuad.dispose(),this._task=null}reset(){const{_renderer:e,_primaryTarget:t,_blendTargets:r}=this,i=e.getRenderTarget(),o=e.getClearAlpha();e.getClearColor(Ce),e.setRenderTarget(t),e.setClearColor(0,0),e.clearColor(),e.setRenderTarget(r[0]),e.setClearColor(0,0),e.clearColor(),e.setRenderTarget(r[1]),e.setClearColor(0,0),e.clearColor(),e.setClearColor(Ce,o),e.setRenderTarget(i),this.samples=0,this._task=null,this.material.stratifiedTexture.stableNoise=this.stableNoise,this.stableNoise&&(this.material.seed=0,this.material.stratifiedTexture.reset())}update(){this.material.onBeforeRender(),!this.isCompiling&&(this._task||(this._task=kr.call(this)),this._task.next())}}const P=new l.Vector2,Pe=new l.Vector2,$=new l.Spherical,j=new l.Color;class Or extends l.DataTexture{constructor(e=512,t=512){super(new Float32Array(e*t*4),e,t,l.RGBAFormat,l.FloatType,l.EquirectangularReflectionMapping,l.RepeatWrapping,l.ClampToEdgeWrapping,l.LinearFilter,l.LinearFilter),this.generationCallback=null}update(){this.dispose(),this.needsUpdate=!0;const{data:e,width:t,height:r}=this.image;for(let i=0;i<t;i++)for(let o=0;o<r;o++){Pe.set(t,r),P.set(i/t,o/r),P.x-=.5,P.y=1-P.y,$.theta=P.x*2*Math.PI,$.phi=P.y*Math.PI,$.radius=1,this.generationCallback($,P,Pe,j);const c=4*(o*t+i);e[c+0]=j.r,e[c+1]=j.g,e[c+2]=j.b,e[c+3]=1}}copy(e){return super.copy(e),this.generationCallback=e.generationCallback,this}}const Ee=new l.Vector3;class Br extends Or{constructor(e=512){super(e,e),this.topColor=new l.Color().set(16777215),this.bottomColor=new l.Color().set(0),this.exponent=2,this.generationCallback=(t,r,i,o)=>{Ee.setFromSpherical(t);const a=Ee.y*.5+.5;o.lerpColors(this.bottomColor,this.topColor,a**this.exponent)}}copy(e){return super.copy(e),this.topColor.copy(e.topColor),this.bottomColor.copy(e.bottomColor),this}}class Nr extends l.ShaderMaterial{get map(){return this.uniforms.map.value}set map(e){this.uniforms.map.value=e}get opacity(){return this.uniforms.opacity.value}set opacity(e){this.uniforms&&(this.uniforms.opacity.value=e)}constructor(e){super({uniforms:{map:{value:null},opacity:{value:1}},vertexShader:`
+				varying vec2 vUv;
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}
+			`,fragmentShader:`
+				uniform sampler2D map;
+				uniform float opacity;
+				varying vec2 vUv;
+
+				vec4 clampedTexelFatch( sampler2D map, ivec2 px, int lod ) {
+
+					vec4 res = texelFetch( map, ivec2( px.x, px.y ), 0 );
+
+					#if defined( TONE_MAPPING )
+
+					res.xyz = toneMapping( res.xyz );
+
+					#endif
+
+			  		return linearToOutputTexel( res );
+
+				}
+
+				void main() {
+
+					vec2 size = vec2( textureSize( map, 0 ) );
+					vec2 pxUv = vUv * size;
+					vec2 pxCurr = floor( pxUv );
+					vec2 pxFrac = fract( pxUv ) - 0.5;
+					vec2 pxOffset;
+					pxOffset.x = pxFrac.x > 0.0 ? 1.0 : - 1.0;
+					pxOffset.y = pxFrac.y > 0.0 ? 1.0 : - 1.0;
+
+					vec2 pxNext = clamp( pxOffset + pxCurr, vec2( 0.0 ), size - 1.0 );
+					vec2 alpha = abs( pxFrac );
+
+					vec4 p1 = mix(
+						clampedTexelFatch( map, ivec2( pxCurr.x, pxCurr.y ), 0 ),
+						clampedTexelFatch( map, ivec2( pxNext.x, pxCurr.y ), 0 ),
+						alpha.x
+					);
+
+					vec4 p2 = mix(
+						clampedTexelFatch( map, ivec2( pxCurr.x, pxNext.y ), 0 ),
+						clampedTexelFatch( map, ivec2( pxNext.x, pxNext.y ), 0 ),
+						alpha.x
+					);
+
+					gl_FragColor = mix( p1, p2, alpha.y );
+					gl_FragColor.a *= opacity;
+					#include <premultiplied_alpha_fragment>
+
+				}
+			`}),this.setValues(e)}}class zr extends l.ShaderMaterial{constructor(){super({uniforms:{envMap:{value:null},flipEnvMap:{value:-1}},vertexShader:`
+				varying vec2 vUv;
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}`,fragmentShader:`
+				#define ENVMAP_TYPE_CUBE_UV
+
+				uniform samplerCube envMap;
+				uniform float flipEnvMap;
+				varying vec2 vUv;
+
+				#include <common>
+				#include <cube_uv_reflection_fragment>
+
+				${Ge}
+
+				void main() {
+
+					vec3 rayDirection = equirectUvToDirection( vUv );
+					rayDirection.x *= flipEnvMap;
+					gl_FragColor = textureCube( envMap, rayDirection );
+
+				}`}),this.depthWrite=!1,this.depthTest=!1}}class ke{constructor(e){this._renderer=e,this._quad=new z(new zr)}generate(e,t=null,r=null){if(!e.isCubeTexture)throw new Error("CubeToEquirectMaterial: Source can only be cube textures.");const i=e.images[0],o=this._renderer,a=this._quad;t===null&&(t=4*i.height),r===null&&(r=2*i.height);const c=new l.WebGLRenderTarget(t,r,{type:l.FloatType,colorSpace:i.colorSpace}),h=i.height,d=Math.log2(h)-2,m=1/h,f=1/(3*Math.max(Math.pow(2,d),112));a.material.defines.CUBEUV_MAX_MIP=`${d}.0`,a.material.defines.CUBEUV_TEXEL_WIDTH=f,a.material.defines.CUBEUV_TEXEL_HEIGHT=m,a.material.uniforms.envMap.value=e,a.material.uniforms.flipEnvMap.value=e.isRenderTargetTexture?1:-1,a.material.needsUpdate=!0;const n=o.getRenderTarget(),x=o.autoClear;o.autoClear=!0,o.setRenderTarget(c),a.render(o),o.setRenderTarget(n),o.autoClear=x;const p=new Uint16Array(t*r*4),v=new Float32Array(t*r*4);o.readRenderTargetPixels(c,0,0,t,r,v),c.dispose();for(let g=0,y=v.length;g<y;g++)p[g]=l.DataUtils.toHalfFloat(v[g]);const u=new l.DataTexture(p,t,r,l.RGBAFormat,l.HalfFloatType);return u.minFilter=l.LinearMipMapLinearFilter,u.magFilter=l.LinearFilter,u.wrapS=l.RepeatWrapping,u.wrapT=l.RepeatWrapping,u.mapping=l.EquirectangularReflectionMapping,u.needsUpdate=!0,u}dispose(){this._quad.dispose()}}function Hr(s){return s.extensions.get("EXT_float_blend")}const O=new l.Vector2;class Wr{get multipleImportanceSampling(){return!!this._pathTracer.material.defines.FEATURE_MIS}set multipleImportanceSampling(e){this._pathTracer.material.setDefine("FEATURE_MIS",e?1:0)}get transmissiveBounces(){return this._pathTracer.material.transmissiveBounces}set transmissiveBounces(e){this._pathTracer.material.transmissiveBounces=e}get bounces(){return this._pathTracer.material.bounces}set bounces(e){this._pathTracer.material.bounces=e}get filterGlossyFactor(){return this._pathTracer.material.filterGlossyFactor}set filterGlossyFactor(e){this._pathTracer.material.filterGlossyFactor=e}get samples(){return this._pathTracer.samples}get target(){return this._pathTracer.target}get tiles(){return this._pathTracer.tiles}get stableNoise(){return this._pathTracer.stableNoise}set stableNoise(e){this._pathTracer.stableNoise=e}get isCompiling(){return!!this._pathTracer.isCompiling}constructor(e){this._renderer=e,this._generator=new yt,this._pathTracer=new De(e),this._queueReset=!1,this._clock=new l.Clock,this._compilePromise=null,this._lowResPathTracer=new De(e),this._lowResPathTracer.tiles.set(1,1),this._quad=new z(new Nr({map:null,transparent:!0,blending:l.NoBlending,premultipliedAlpha:e.getContextAttributes().premultipliedAlpha})),this._materials=null,this._previousEnvironment=null,this._previousBackground=null,this._internalBackground=null,this.renderDelay=100,this.minSamples=5,this.fadeDuration=500,this.enablePathTracing=!0,this.pausePathTracing=!1,this.dynamicLowRes=!1,this.lowResScale=.25,this.renderScale=1,this.synchronizeRenderSize=!0,this.rasterizeScene=!0,this.renderToCanvas=!0,this.textureSize=new l.Vector2(1024,1024),this.rasterizeSceneCallback=(t,r)=>{this._renderer.render(t,r)},this.renderToCanvasCallback=(t,r,i)=>{const o=r.autoClear;r.autoClear=!1,i.render(r),r.autoClear=o},this.setScene(new l.Scene,new l.PerspectiveCamera)}setBVHWorker(e){this._generator.setBVHWorker(e)}setScene(e,t,r={}){e.updateMatrixWorld(!0),t.updateMatrixWorld();const i=this._generator;if(i.setObjects(e),this._buildAsync)return i.generateAsync(r.onProgress).then(o=>this._updateFromResults(e,t,o));{const o=i.generate();return this._updateFromResults(e,t,o)}}setSceneAsync(...e){this._buildAsync=!0;const t=this.setScene(...e);return this._buildAsync=!1,t}setCamera(e){this.camera=e,this.updateCamera()}updateCamera(){const e=this.camera;e.updateMatrixWorld(),this._pathTracer.setCamera(e),this._lowResPathTracer.setCamera(e),this.reset()}updateMaterials(){const e=this._pathTracer.material,t=this._renderer,r=this._materials,i=this.textureSize,o=Vt(r);e.textures.setTextures(t,o,i.x,i.y),e.materials.updateFrom(r,o),this.reset()}updateLights(){const e=this.scene,t=this._renderer,r=this._pathTracer.material,i=qt(e),o=Gt(i);r.lights.updateFrom(i,o),r.iesProfiles.setTextures(t,o),this.reset()}updateEnvironment(){const e=this.scene,t=this._pathTracer.material;if(this._internalBackground&&(this._internalBackground.dispose(),this._internalBackground=null),t.backgroundBlur=e.backgroundBlurriness,t.backgroundIntensity=e.backgroundIntensity??1,t.backgroundRotation.makeRotationFromEuler(e.backgroundRotation).invert(),e.background===null)t.backgroundMap=null,t.backgroundAlpha=0;else if(e.background.isColor){this._colorBackground=this._colorBackground||new Br(16);const r=this._colorBackground;r.topColor.equals(e.background)||(r.topColor.set(e.background),r.bottomColor.set(e.background),r.update()),t.backgroundMap=r,t.backgroundAlpha=1}else if(e.background.isCubeTexture){if(e.background!==this._previousBackground){const r=new ke(this._renderer).generate(e.background);this._internalBackground=r,t.backgroundMap=r,t.backgroundAlpha=1}}else t.backgroundMap=e.background,t.backgroundAlpha=1;if(t.environmentIntensity=e.environmentIntensity??1,t.environmentRotation.makeRotationFromEuler(e.environmentRotation).invert(),this._previousEnvironment!==e.environment)if(e.environment!==null)if(e.environment.isCubeTexture){const r=new ke(this._renderer).generate(e.environment);t.envMapInfo.updateFrom(r)}else t.envMapInfo.updateFrom(e.environment);else t.environmentIntensity=0;this._previousEnvironment=e.environment,this._previousBackground=e.background,this.reset()}_updateFromResults(e,t,r){const{materials:i,geometry:o,bvh:a,bvhChanged:c}=r;this._materials=i;const d=this._pathTracer.material;return c&&(d.bvh.updateFrom(a),d.attributesArray.updateFrom(o.attributes.normal,o.attributes.tangent,o.attributes.uv,o.attributes.color),d.materialIndexAttribute.updateFrom(o.attributes.materialIndex)),this._previousScene=e,this.scene=e,this.camera=t,this.updateCamera(),this.updateMaterials(),this.updateEnvironment(),this.updateLights(),r}renderSample(){const e=this._lowResPathTracer,t=this._pathTracer,r=this._renderer,i=this._clock,o=this._quad;this._updateScale(),this._queueReset&&(t.reset(),e.reset(),this._queueReset=!1,o.material.opacity=0,i.start());const a=i.getDelta()*1e3,c=i.getElapsedTime()*1e3;if(!this.pausePathTracing&&this.enablePathTracing&&this.renderDelay<=c&&!this.isCompiling&&t.update(),t.alpha=t.material.backgroundAlpha!==1||!Hr(r),e.alpha=t.alpha,this.renderToCanvas){const h=this._renderer,d=this.minSamples;if(c>=this.renderDelay&&this.samples>=this.minSamples&&(this.fadeDuration!==0?o.material.opacity=Math.min(o.material.opacity+a/this.fadeDuration,1):o.material.opacity=1),!this.enablePathTracing||this.samples<d||o.material.opacity<1){if(this.dynamicLowRes&&!this.isCompiling){e.samples<1&&(e.material=t.material,e.update());const m=o.material.opacity;o.material.opacity=1-o.material.opacity,o.material.map=e.target.texture,o.render(h),o.material.opacity=m}(!this.dynamicLowRes&&this.rasterizeScene||this.dynamicLowRes&&this.isCompiling)&&this.rasterizeSceneCallback(this.scene,this.camera)}this.enablePathTracing&&o.material.opacity>0&&(o.material.opacity<1&&(o.material.blending=this.dynamicLowRes?l.AdditiveBlending:l.NormalBlending),o.material.map=t.target.texture,this.renderToCanvasCallback(t.target,h,o),o.material.blending=l.NoBlending)}}reset(){this._queueReset=!0,this._pathTracer.samples=0}dispose(){this._renderQuad.dispose(),this._renderQuad.material.dispose(),this._pathTracer.dispose()}_updateScale(){if(this.synchronizeRenderSize){this._renderer.getDrawingBufferSize(O);const e=Math.floor(this.renderScale*O.x),t=Math.floor(this.renderScale*O.y);if(this._pathTracer.getSize(O),O.x!==e||O.y!==t){const r=this.lowResScale;this._pathTracer.setSize(e,t),this._lowResPathTracer.setSize(Math.floor(e*r),Math.floor(t*r))}}}}function B(s){return s!==null&&typeof s=="object"&&"getInternalRenderer"in s&&typeof s.getInternalRenderer=="function"}class Lr{constructor(){this.pathTracer=null,this.enabled=!1,this.renderer=null,this.sampleCount=0,this.createAttempts=0,this.maxCreateAttempts=10,this.sceneInitialized=!1,this.environmentWaitFrames=0,this.maxEnvironmentWaitFrames=300,this.disposed=!1,this.disposeTimer=null,this.convertedEnvTexture=null,this.lastResetTime=0,this.events=new T.TypedEventEmitter,this.settings={samples:T.DEFAULT_PATH_TRACING_SAMPLES,bounces:4,transmissiveBounces:2,renderScale:.5,lowResScale:.5,dynamicLowRes:!0,enablePathTracing:!0}}convertToDataTexture(e){if(!e.image||!(e.image instanceof HTMLImageElement))return null;const t=e.image,r=document.createElement("canvas");r.width=t.width,r.height=t.height;const i=r.getContext("2d");if(!i)return null;i.drawImage(t,0,0);const o=i.getImageData(0,0,r.width,r.height),a=new Float32Array(o.data.length);for(let h=0;h<o.data.length;h++)a[h]=o.data[h]/255;const c=new A.DataTexture(a,r.width,r.height,A.RGBAFormat,A.FloatType);return c.mapping=e.mapping,c.wrapS=e.wrapS,c.wrapT=e.wrapT,c.magFilter=e.magFilter,c.minFilter=e.minFilter,c.anisotropy=e.anisotropy,c.needsUpdate=!0,c}async initialize(e){try{return this.renderer=e.renderer,this.enabled=e.enabled,this.isSupported()?T.Result.ok(void 0):T.Result.err(new T.ThreeViewerError("Path tracing is not supported in this environment",T.ErrorCode.PATH_TRACING_INIT_FAILED,{reason:"WebGL2 or required extensions not available"}))}catch(t){return T.Result.err(new T.ThreeViewerError("Failed to initialize path tracing",T.ErrorCode.PATH_TRACING_INIT_FAILED,{originalError:t}))}}createPathTracer(){if(this.disposed)return T.Result.err(new T.ThreeViewerError("Instance is disposed",T.ErrorCode.INVALID_STATE));try{const e=B(this.renderer)?this.renderer.getInternalRenderer():null;return e?(this.pathTracer=new Wr(e),this.pathTracer.tiles&&"set"in this.pathTracer.tiles&&this.pathTracer.tiles.set(1,1),this.pathTracer.bounces=this.settings.bounces,this.settings.transmissiveBounces!==void 0&&(this.pathTracer.transmissiveBounces=this.settings.transmissiveBounces),this.pathTracer.renderScale=this.settings.renderScale,this.pathTracer.dynamicLowRes=this.settings.dynamicLowRes,this.pathTracer.lowResScale=this.settings.lowResScale,this.pathTracer.environmentIntensity!==void 0&&(this.pathTracer.environmentIntensity=2),e.toneMapping!==A.ACESFilmicToneMapping&&(e.toneMapping=A.ACESFilmicToneMapping,e.toneMappingExposure=1.5),e.autoClear=!1,e.setRenderTarget(null),T.Result.ok(void 0)):T.Result.err(new T.ThreeViewerError("Renderer not ready, will retry",T.ErrorCode.RENDERER_NOT_INITIALIZED))}catch(e){return T.Result.err(new T.ThreeViewerError("Failed to create path tracer",T.ErrorCode.PATH_TRACING_INIT_FAILED,{originalError:e}))}}setEnabled(e){if(this.enabled!==e&&(this.enabled=e,!e)){this.reset(),this.sceneInitialized=!1;const t=B(this.renderer)?this.renderer.getInternalRenderer():null;t&&(t.autoClear=!0)}}disableAfterSelfPause(e){this.enabled=!1,this.events.emit("pathtracing:paused",{samples:this.sampleCount,reason:e})}updateSettings(e){this.settings={...this.settings,...e},e.enablePathTracing!==void 0&&this.setEnabled(e.enablePathTracing),this.pathTracer&&(e.bounces!==void 0&&(this.pathTracer.bounces=e.bounces),e.transmissiveBounces!==void 0&&(this.pathTracer.transmissiveBounces=e.transmissiveBounces),e.renderScale!==void 0&&(this.pathTracer.renderScale=e.renderScale),e.lowResScale!==void 0&&(this.pathTracer.lowResScale=e.lowResScale),e.dynamicLowRes!==void 0&&(this.pathTracer.dynamicLowRes=e.dynamicLowRes))}async render(e,t){if(this.disposed)return T.Result.ok(void 0);const r=this.renderer;if(!r)return T.Result.ok(void 0);if(!this.enabled)return this.renderWhileDisabled(e,t);if(!this.pathTracer){const i=this.ensurePathTracerCreated(r,e,t);if(i)return i}try{const i=this.extractThreeObjects(e,t);if(!i)return T.Result.err(new T.ThreeViewerError("Could not extract Three.js scene or camera",T.ErrorCode.INVALID_PARAMETER));const{threeScene:o,threeCamera:a}=i;if(!this.sceneInitialized){const c=await this.initializeSceneForPathTracing(r,e,t,o,a);if(c)return c}return this.sceneInitialized?this.accumulateSample(r,e,t):r.render(e,t)}catch(i){return T.Result.err(new T.ThreeViewerError("Failed to render with path tracing",T.ErrorCode.RENDER_FAILED,{originalError:i}))}}ensurePathTracerCreated(e,t,r){var o;this.createAttempts++;const i=this.createPathTracer();return i.ok?(this.createAttempts=0,null):(this.createAttempts<this.maxCreateAttempts&&((o=i.error)==null?void 0:o.message)==="Renderer not ready, will retry"||(this.disableAfterSelfPause("gave-up"),this.createAttempts=0),e.render(t,r))}extractThreeObjects(e,t){let r=null,i=null;return T.hasInternalRenderer(e)&&(r=e.getInternalRenderer()),T.hasInternalRenderer(t)&&(i=t.getInternalRenderer()),!r||!i?null:{threeScene:r,threeCamera:i}}async initializeSceneForPathTracing(e,t,r,i,o){var a,c,h;try{const d=i.__originalEnvironmentTexture;if(!i.environment&&!d)return this.environmentWaitFrames++,this.environmentWaitFrames>=this.maxEnvironmentWaitFrames&&(this.disableAfterSelfPause("gave-up"),this.environmentWaitFrames=0),e.render(t,r);if(this.environmentWaitFrames=0,d&&d.mapping===A.EquirectangularReflectionMapping){let m=d;if(d.image instanceof HTMLImageElement)if(this.convertedEnvTexture)m=this.convertedEnvTexture;else{const n=this.convertToDataTexture(d);if(!n)throw new Error("Failed to convert environment texture to DataTexture");this.convertedEnvTexture=n,m=n}const f=i.environment;i.environment=m;try{const n=d.image;n&&!("data"in n&&n.data)&&n instanceof HTMLImageElement&&!n.complete&&await new Promise(x=>{n.onload=()=>{d.needsUpdate=!0,x()},n.onerror=()=>{x()},n.complete&&x()}),(a=this.pathTracer)==null||a.setScene(i,o),this.sceneInitialized=!0,i.environment=f}catch(n){throw i.environment=f,n}}else if(((c=i.environment)==null?void 0:c.mapping)===A.EquirectangularReflectionMapping)(h=this.pathTracer)==null||h.setScene(i,o),this.sceneInitialized=!0;else return this.disableAfterSelfPause("gave-up"),e.render(t,r)}catch(d){return console.warn("Scene initialization error:",d),e.render(t,r)}return null}accumulateSample(e,t,r){var o;if(this.sampleCount===0)try{(o=this.pathTracer)==null||o.updateLights()}catch(a){console.warn("Failed to update lights for path tracing:",a)}const i=e.render(t,r);return i.ok?(this.accumulateOneSample(),this.sampleCount++,this.sampleCount===this.settings.samples?this.captureCompletedFrame():T.Result.ok(void 0)):i}renderWhileDisabled(e,t){return this.renderer?this.sampleCount>=this.settings.samples&&this.sampleCount>0?(this.pathTracer&&B(this.renderer)&&this.renderer.getInternalRenderer().setRenderTarget(null),T.Result.ok(void 0)):this.renderer.render(e,t):T.Result.ok(void 0)}accumulateOneSample(){const e=B(this.renderer)?this.renderer.getInternalRenderer():null;if(!e)throw new Error("Three.js renderer not available");const t=e.autoClear,r=e.getRenderTarget();e.autoClear=!1,this.pathTracer&&this.pathTracer.renderSample(),e.autoClear=t,e.setRenderTarget(r)}captureCompletedFrame(){const e=B(this.renderer)?this.renderer.getInternalRenderer():null;if(e&&this.pathTracer)try{const t=e.autoClear;e.autoClear=!0,e.setRenderTarget(null),e.clear(!0,!0,!0);const r=this.pathTracer.copyQuad;r&&typeof r.render=="function"?r.render(e):(e.autoClear=!1,this.pathTracer.renderSample()),e.autoClear=t}catch(t){console.warn("Failed to present path traced result:",t)}return this.disableAfterSelfPause("completed"),this.disposeTimer=setTimeout(()=>{this.disposeTimer=null,!this.disposed&&this.disposePathTracingResources()},100),T.Result.ok(void 0)}getSampleCount(){return this.sampleCount}isEnabled(){return this.enabled}isPathTracerDisposed(){return this.disposed}reset(){const e=performance.now();e-this.lastResetTime<50||(this.lastResetTime=e,this.sampleCount=0,this.pathTracer&&(this.pathTracer.reset(),this.sceneInitialized=!1))}disposePathTracingResources(){if(this.pathTracer){try{this.pathTracer.dispose()}catch(e){console.warn("Failed to dispose path tracer:",e)}this.pathTracer=null}if(this.convertedEnvTexture&&(this.convertedEnvTexture.dispose(),this.convertedEnvTexture=null),this.renderer){const e=B(this.renderer)?this.renderer.getInternalRenderer():null;e&&(e.autoClear=!0)}this.sceneInitialized=!1,this.createAttempts=0}dispose(){if(this.disposed=!0,this.events.removeAllListeners(),this.disposeTimer!==null&&(clearTimeout(this.disposeTimer),this.disposeTimer=null),this.pathTracer){try{this.pathTracer.dispose()}catch(e){console.warn("Failed to dispose path tracer during service disposal:",e)}this.pathTracer=null}this.convertedEnvTexture&&(this.convertedEnvTexture.dispose(),this.convertedEnvTexture=null),this.sampleCount=0,this.sceneInitialized=!1,this.createAttempts=0,this.renderer=null}isSupported(){return document.createElement("canvas").getContext("webgl2")!==null}}exports.ThreePathTracingService=Lr;