shop-components 0.0.24 → 0.0.26

package/dist/assets/lightroom_14b_512-e_XJNpR6.js

@@ -0,0 +1,3981 @@
+import{n as Ss}from"./property-EQZXZ3dt.js";import{an as Ms,M as zr,B as De,aV as Qi,aU as Rs,aQ as Ii,j as Is,V as C,aW as ze,aX as Oe,Q as hi,aY as di,k as L,l as Ps,e as re,aR as xe,aS as ft,n as Cs,i as Y,m as he,aT as Or,u as Nr,az as mi,aZ as mt,a_ as ae,as as U,a$ as Xt,F as q,b0 as pi,b1 as Fs,b2 as Ji,b3 as Ds,aN as V,b4 as gi,b5 as Lr,b6 as Wt,b7 as Hr,b8 as Es,b9 as Bs,ba as Je,a7 as ks,y as qt,N as vt,v as pt,P as Ur,a as ve,H as fe,b as le,af as _e,aw as Me,aL as zs,bb as Os,C as xt,bc as Ns,bd as Ls,be as Wr,E as Vr,bf as Hs,_ as Us,f as Ws,z as Vs}from"./three.module-9MfFqitO.js";/**
+ * @license
+ * Copyright 2017 Google LLC
+ * SPDX-License-Identifier: BSD-3-Clause
+ */function Vo(s){return Ss({...s,state:!0,attribute:!1})}class qo{constructor(){this.isPass=!0,this.enabled=!0,this.needsSwap=!0,this.clear=!1,this.renderToScreen=!1}setSize(){}render(){console.error("THREE.Pass: .render() must be implemented in derived pass.")}dispose(){}}const qs=new Ms(-1,1,1,-1,0,1);class js extends De{constructor(){super(),this.setAttribute("position",new Qi([-1,3,0,-1,-1,0,3,-1,0],3)),this.setAttribute("uv",new Qi([0,2,0,0,2,0],2))}}const Gs=new js;class et{constructor(e){this._mesh=new zr(Gs,e)}dispose(){this._mesh.geometry.dispose()}render(e){e.render(this._mesh,qs)}get material(){return this._mesh.material}set material(e){this._mesh.material=e}}const er={type:"change"},Kt={type:"start"},tr={type:"end"},bt=new Rs,ir=new Ii,Ys=Math.cos(70*Ps.DEG2RAD);class jo extends Is{constructor(e,i){super(),this.object=e,this.domElement=i,this.domElement.style.touchAction="none",this.enabled=!0,this.target=new C,this.cursor=new C,this.minDistance=0,this.maxDistance=1/0,this.minZoom=0,this.maxZoom=1/0,this.minTargetRadius=0,this.maxTargetRadius=1/0,this.minPolarAngle=0,this.maxPolarAngle=Math.PI,this.minAzimuthAngle=-1/0,this.maxAzimuthAngle=1/0,this.enableDamping=!1,this.dampingFactor=.05,this.enableZoom=!0,this.zoomSpeed=1,this.enableRotate=!0,this.rotateSpeed=1,this.enablePan=!0,this.panSpeed=1,this.screenSpacePanning=!0,this.keyPanSpeed=7,this.zoomToCursor=!1,this.autoRotate=!1,this.autoRotateSpeed=2,this.keys={LEFT:"ArrowLeft",UP:"ArrowUp",RIGHT:"ArrowRight",BOTTOM:"ArrowDown"},this.mouseButtons={LEFT:ze.ROTATE,MIDDLE:ze.DOLLY,RIGHT:ze.PAN},this.touches={ONE:Oe.ROTATE,TWO:Oe.DOLLY_PAN},this.target0=this.target.clone(),this.position0=this.object.position.clone(),this.zoom0=this.object.zoom,this._domElementKeyEvents=null,this.getPolarAngle=function(){return c.phi},this.getAzimuthalAngle=function(){return c.theta},this.getDistance=function(){return this.object.position.distanceTo(this.target)},this.listenToKeyEvents=function(x){x.addEventListener("keydown",$t),this._domElementKeyEvents=x},this.stopListenToKeyEvents=function(){this._domElementKeyEvents.removeEventListener("keydown",$t),this._domElementKeyEvents=null},this.saveState=function(){t.target0.copy(t.target),t.position0.copy(t.object.position),t.zoom0=t.object.zoom},this.reset=function(){t.target.copy(t.target0),t.object.position.copy(t.position0),t.object.zoom=t.zoom0,t.object.updateProjectionMatrix(),t.dispatchEvent(er),t.update(),a=n.NONE},this.update=function(){const x=new C,D=new hi().setFromUnitVectors(e.up,new C(0,1,0)),k=D.clone().invert(),O=new C,j=new hi,be=new C,ie=2*Math.PI;return function(_s=null){const Ki=t.object.position;x.copy(Ki).sub(t.target),x.applyQuaternion(D),c.setFromVector3(x),t.autoRotate&&a===n.NONE&&M(I(_s)),t.enableDamping?(c.theta+=l.theta*t.dampingFactor,c.phi+=l.phi*t.dampingFactor):(c.theta+=l.theta,c.phi+=l.phi);let de=t.minAzimuthAngle,me=t.maxAzimuthAngle;isFinite(de)&&isFinite(me)&&(de<-Math.PI?de+=ie:de>Math.PI&&(de-=ie),me<-Math.PI?me+=ie:me>Math.PI&&(me-=ie),de<=me?c.theta=Math.max(de,Math.min(me,c.theta)):c.theta=c.theta>(de+me)/2?Math.max(de,c.theta):Math.min(me,c.theta)),c.phi=Math.max(t.minPolarAngle,Math.min(t.maxPolarAngle,c.phi)),c.makeSafe(),t.enableDamping===!0?t.target.addScaledVector(f,t.dampingFactor):t.target.add(f),t.target.sub(t.cursor),t.target.clampLength(t.minTargetRadius,t.maxTargetRadius),t.target.add(t.cursor);let rt=!1;if(t.zoomToCursor&&_||t.object.isOrthographicCamera)c.radius=Ee(c.radius);else{const pe=c.radius;c.radius=Ee(c.radius*h),rt=pe!=c.radius}if(x.setFromSpherical(c),x.applyQuaternion(k),Ki.copy(t.target).add(x),t.object.lookAt(t.target),t.enableDamping===!0?(l.theta*=1-t.dampingFactor,l.phi*=1-t.dampingFactor,f.multiplyScalar(1-t.dampingFactor)):(l.set(0,0,0),f.set(0,0,0)),t.zoomToCursor&&_){let pe=null;if(t.object.isPerspectiveCamera){const st=x.length();pe=Ee(st*h);const yt=st-pe;t.object.position.addScaledVector(w,yt),t.object.updateMatrixWorld(),rt=!!yt}else if(t.object.isOrthographicCamera){const st=new C(T.x,T.y,0);st.unproject(t.object);const yt=t.object.zoom;t.object.zoom=Math.max(t.minZoom,Math.min(t.maxZoom,t.object.zoom/h)),t.object.updateProjectionMatrix(),rt=yt!==t.object.zoom;const Zi=new C(T.x,T.y,0);Zi.unproject(t.object),t.object.position.sub(Zi).add(st),t.object.updateMatrixWorld(),pe=x.length()}else console.warn("WARNING: OrbitControls.js encountered an unknown camera type - zoom to cursor disabled."),t.zoomToCursor=!1;pe!==null&&(this.screenSpacePanning?t.target.set(0,0,-1).transformDirection(t.object.matrix).multiplyScalar(pe).add(t.object.position):(bt.origin.copy(t.object.position),bt.direction.set(0,0,-1).transformDirection(t.object.matrix),Math.abs(t.object.up.dot(bt.direction))<Ys?e.lookAt(t.target):(ir.setFromNormalAndCoplanarPoint(t.object.up,t.target),bt.intersectPlane(ir,t.target))))}else if(t.object.isOrthographicCamera){const pe=t.object.zoom;t.object.zoom=Math.max(t.minZoom,Math.min(t.maxZoom,t.object.zoom/h)),pe!==t.object.zoom&&(t.object.updateProjectionMatrix(),rt=!0)}return h=1,_=!1,rt||O.distanceToSquared(t.object.position)>r||8*(1-j.dot(t.object.quaternion))>r||be.distanceToSquared(t.target)>r?(t.dispatchEvent(er),O.copy(t.object.position),j.copy(t.object.quaternion),be.copy(t.target),!0):!1}}(),this.dispose=function(){t.domElement.removeEventListener("contextmenu",$i),t.domElement.removeEventListener("pointerdown",Vi),t.domElement.removeEventListener("pointercancel",it),t.domElement.removeEventListener("wheel",qi),t.domElement.removeEventListener("pointermove",Yt),t.domElement.removeEventListener("pointerup",it),t.domElement.getRootNode().removeEventListener("keydown",ji,{capture:!0}),t._domElementKeyEvents!==null&&(t._domElementKeyEvents.removeEventListener("keydown",$t),t._domElementKeyEvents=null)};const t=this,n={NONE:-1,ROTATE:0,DOLLY:1,PAN:2,TOUCH_ROTATE:3,TOUCH_PAN:4,TOUCH_DOLLY_PAN:5,TOUCH_DOLLY_ROTATE:6};let a=n.NONE;const r=1e-6,c=new di,l=new di;let h=1;const f=new C,u=new L,o=new L,m=new L,v=new L,b=new L,d=new L,p=new L,g=new L,y=new L,w=new C,T=new L;let _=!1;const A=[],S={};let R=!1;function I(x){return x!==null?2*Math.PI/60*t.autoRotateSpeed*x:2*Math.PI/60/60*t.autoRotateSpeed}function P(x){const D=Math.abs(x*.01);return Math.pow(.95,t.zoomSpeed*D)}function M(x){l.theta-=x}function F(x){l.phi-=x}const E=function(){const x=new C;return function(k,O){x.setFromMatrixColumn(O,0),x.multiplyScalar(-k),f.add(x)}}(),B=function(){const x=new C;return function(k,O){t.screenSpacePanning===!0?x.setFromMatrixColumn(O,1):(x.setFromMatrixColumn(O,0),x.crossVectors(t.object.up,x)),x.multiplyScalar(k),f.add(x)}}(),z=function(){const x=new C;return function(k,O){const j=t.domElement;if(t.object.isPerspectiveCamera){const be=t.object.position;x.copy(be).sub(t.target);let ie=x.length();ie*=Math.tan(t.object.fov/2*Math.PI/180),E(2*k*ie/j.clientHeight,t.object.matrix),B(2*O*ie/j.clientHeight,t.object.matrix)}else t.object.isOrthographicCamera?(E(k*(t.object.right-t.object.left)/t.object.zoom/j.clientWidth,t.object.matrix),B(O*(t.object.top-t.object.bottom)/t.object.zoom/j.clientHeight,t.object.matrix)):(console.warn("WARNING: OrbitControls.js encountered an unknown camera type - pan disabled."),t.enablePan=!1)}}();function X(x){t.object.isPerspectiveCamera||t.object.isOrthographicCamera?h/=x:(console.warn("WARNING: OrbitControls.js encountered an unknown camera type - dolly/zoom disabled."),t.enableZoom=!1)}function ee(x){t.object.isPerspectiveCamera||t.object.isOrthographicCamera?h*=x:(console.warn("WARNING: OrbitControls.js encountered an unknown camera type - dolly/zoom disabled."),t.enableZoom=!1)}function te(x,D){if(!t.zoomToCursor)return;_=!0;const k=t.domElement.getBoundingClientRect(),O=x-k.left,j=D-k.top,be=k.width,ie=k.height;T.x=O/be*2-1,T.y=-(j/ie)*2+1,w.set(T.x,T.y,1).unproject(t.object).sub(t.object.position).normalize()}function Ee(x){return Math.max(t.minDistance,Math.min(t.maxDistance,x))}function Be(x){u.set(x.clientX,x.clientY)}function os(x){te(x.clientX,x.clientX),p.set(x.clientX,x.clientY)}function zi(x){v.set(x.clientX,x.clientY)}function cs(x){o.set(x.clientX,x.clientY),m.subVectors(o,u).multiplyScalar(t.rotateSpeed);const D=t.domElement;M(2*Math.PI*m.x/D.clientHeight),F(2*Math.PI*m.y/D.clientHeight),u.copy(o),t.update()}function ls(x){g.set(x.clientX,x.clientY),y.subVectors(g,p),y.y>0?X(P(y.y)):y.y<0&&ee(P(y.y)),p.copy(g),t.update()}function us(x){b.set(x.clientX,x.clientY),d.subVectors(b,v).multiplyScalar(t.panSpeed),z(d.x,d.y),v.copy(b),t.update()}function fs(x){te(x.clientX,x.clientY),x.deltaY<0?ee(P(x.deltaY)):x.deltaY>0&&X(P(x.deltaY)),t.update()}function hs(x){let D=!1;switch(x.code){case t.keys.UP:x.ctrlKey||x.metaKey||x.shiftKey?F(2*Math.PI*t.rotateSpeed/t.domElement.clientHeight):z(0,t.keyPanSpeed),D=!0;break;case t.keys.BOTTOM:x.ctrlKey||x.metaKey||x.shiftKey?F(-2*Math.PI*t.rotateSpeed/t.domElement.clientHeight):z(0,-t.keyPanSpeed),D=!0;break;case t.keys.LEFT:x.ctrlKey||x.metaKey||x.shiftKey?M(2*Math.PI*t.rotateSpeed/t.domElement.clientHeight):z(t.keyPanSpeed,0),D=!0;break;case t.keys.RIGHT:x.ctrlKey||x.metaKey||x.shiftKey?M(-2*Math.PI*t.rotateSpeed/t.domElement.clientHeight):z(-t.keyPanSpeed,0),D=!0;break}D&&(x.preventDefault(),t.update())}function Oi(x){if(A.length===1)u.set(x.pageX,x.pageY);else{const D=ke(x),k=.5*(x.pageX+D.x),O=.5*(x.pageY+D.y);u.set(k,O)}}function Ni(x){if(A.length===1)v.set(x.pageX,x.pageY);else{const D=ke(x),k=.5*(x.pageX+D.x),O=.5*(x.pageY+D.y);v.set(k,O)}}function Li(x){const D=ke(x),k=x.pageX-D.x,O=x.pageY-D.y,j=Math.sqrt(k*k+O*O);p.set(0,j)}function ds(x){t.enableZoom&&Li(x),t.enablePan&&Ni(x)}function ms(x){t.enableZoom&&Li(x),t.enableRotate&&Oi(x)}function Hi(x){if(A.length==1)o.set(x.pageX,x.pageY);else{const k=ke(x),O=.5*(x.pageX+k.x),j=.5*(x.pageY+k.y);o.set(O,j)}m.subVectors(o,u).multiplyScalar(t.rotateSpeed);const D=t.domElement;M(2*Math.PI*m.x/D.clientHeight),F(2*Math.PI*m.y/D.clientHeight),u.copy(o)}function Ui(x){if(A.length===1)b.set(x.pageX,x.pageY);else{const D=ke(x),k=.5*(x.pageX+D.x),O=.5*(x.pageY+D.y);b.set(k,O)}d.subVectors(b,v).multiplyScalar(t.panSpeed),z(d.x,d.y),v.copy(b)}function Wi(x){const D=ke(x),k=x.pageX-D.x,O=x.pageY-D.y,j=Math.sqrt(k*k+O*O);g.set(0,j),y.set(0,Math.pow(g.y/p.y,t.zoomSpeed)),X(y.y),p.copy(g);const be=(x.pageX+D.x)*.5,ie=(x.pageY+D.y)*.5;te(be,ie)}function ps(x){t.enableZoom&&Wi(x),t.enablePan&&Ui(x)}function gs(x){t.enableZoom&&Wi(x),t.enableRotate&&Hi(x)}function Vi(x){t.enabled!==!1&&(A.length===0&&(t.domElement.setPointerCapture(x.pointerId),t.domElement.addEventListener("pointermove",Yt),t.domElement.addEventListener("pointerup",it)),!As(x)&&(Ts(x),x.pointerType==="touch"?Yi(x):vs(x)))}function Yt(x){t.enabled!==!1&&(x.pointerType==="touch"?bs(x):xs(x))}function it(x){switch(ws(x),A.length){case 0:t.domElement.releasePointerCapture(x.pointerId),t.domElement.removeEventListener("pointermove",Yt),t.domElement.removeEventListener("pointerup",it),t.dispatchEvent(tr),a=n.NONE;break;case 1:const D=A[0],k=S[D];Yi({pointerId:D,pageX:k.x,pageY:k.y});break}}function vs(x){let D;switch(x.button){case 0:D=t.mouseButtons.LEFT;break;case 1:D=t.mouseButtons.MIDDLE;break;case 2:D=t.mouseButtons.RIGHT;break;default:D=-1}switch(D){case ze.DOLLY:if(t.enableZoom===!1)return;os(x),a=n.DOLLY;break;case ze.ROTATE:if(x.ctrlKey||x.metaKey||x.shiftKey){if(t.enablePan===!1)return;zi(x),a=n.PAN}else{if(t.enableRotate===!1)return;Be(x),a=n.ROTATE}break;case ze.PAN:if(x.ctrlKey||x.metaKey||x.shiftKey){if(t.enableRotate===!1)return;Be(x),a=n.ROTATE}else{if(t.enablePan===!1)return;zi(x),a=n.PAN}break;default:a=n.NONE}a!==n.NONE&&t.dispatchEvent(Kt)}function xs(x){switch(a){case n.ROTATE:if(t.enableRotate===!1)return;cs(x);break;case n.DOLLY:if(t.enableZoom===!1)return;ls(x);break;case n.PAN:if(t.enablePan===!1)return;us(x);break}}function qi(x){t.enabled===!1||t.enableZoom===!1||a!==n.NONE||(x.preventDefault(),t.dispatchEvent(Kt),fs(ys(x)),t.dispatchEvent(tr))}function ys(x){const D=x.deltaMode,k={clientX:x.clientX,clientY:x.clientY,deltaY:x.deltaY};switch(D){case 1:k.deltaY*=16;break;case 2:k.deltaY*=100;break}return x.ctrlKey&&!R&&(k.deltaY*=10),k}function ji(x){x.key==="Control"&&(R=!0,t.domElement.getRootNode().addEventListener("keyup",Gi,{passive:!0,capture:!0}))}function Gi(x){x.key==="Control"&&(R=!1,t.domElement.getRootNode().removeEventListener("keyup",Gi,{passive:!0,capture:!0}))}function $t(x){t.enabled===!1||t.enablePan===!1||hs(x)}function Yi(x){switch(Xi(x),A.length){case 1:switch(t.touches.ONE){case Oe.ROTATE:if(t.enableRotate===!1)return;Oi(x),a=n.TOUCH_ROTATE;break;case Oe.PAN:if(t.enablePan===!1)return;Ni(x),a=n.TOUCH_PAN;break;default:a=n.NONE}break;case 2:switch(t.touches.TWO){case Oe.DOLLY_PAN:if(t.enableZoom===!1&&t.enablePan===!1)return;ds(x),a=n.TOUCH_DOLLY_PAN;break;case Oe.DOLLY_ROTATE:if(t.enableZoom===!1&&t.enableRotate===!1)return;ms(x),a=n.TOUCH_DOLLY_ROTATE;break;default:a=n.NONE}break;default:a=n.NONE}a!==n.NONE&&t.dispatchEvent(Kt)}function bs(x){switch(Xi(x),a){case n.TOUCH_ROTATE:if(t.enableRotate===!1)return;Hi(x),t.update();break;case n.TOUCH_PAN:if(t.enablePan===!1)return;Ui(x),t.update();break;case n.TOUCH_DOLLY_PAN:if(t.enableZoom===!1&&t.enablePan===!1)return;ps(x),t.update();break;case n.TOUCH_DOLLY_ROTATE:if(t.enableZoom===!1&&t.enableRotate===!1)return;gs(x),t.update();break;default:a=n.NONE}}function $i(x){t.enabled!==!1&&x.preventDefault()}function Ts(x){A.push(x.pointerId)}function ws(x){delete S[x.pointerId];for(let D=0;D<A.length;D++)if(A[D]==x.pointerId){A.splice(D,1);return}}function As(x){for(let D=0;D<A.length;D++)if(A[D]==x.pointerId)return!0;return!1}function Xi(x){let D=S[x.pointerId];D===void 0&&(D=new L,S[x.pointerId]=D),D.set(x.pageX,x.pageY)}function ke(x){const D=x.pointerId===A[0]?A[1]:A[0];return S[D]}t.domElement.addEventListener("contextmenu",$i),t.domElement.addEventListener("pointerdown",Vi),t.domElement.addEventListener("pointercancel",it),t.domElement.addEventListener("wheel",qi,{passive:!1}),t.domElement.getRootNode().addEventListener("keydown",ji,{passive:!0,capture:!0}),this.update()}}const qr=0,$s=1,jr=2,rr=2,Zt=1.25,sr=1,Se=6*4+4+4,jt=65535,Xs=Math.pow(2,-24),Qt=Symbol("SKIP_GENERATION");function Gr(s){return s.index?s.index.count:s.attributes.position.count}function tt(s){return Gr(s)/3}function Yr(s,e=ArrayBuffer){return s>65535?new Uint32Array(new e(4*s)):new Uint16Array(new e(2*s))}function Ks(s,e){if(!s.index){const i=s.attributes.position.count,t=e.useSharedArrayBuffer?SharedArrayBuffer:ArrayBuffer,n=Yr(i,t);s.setIndex(new re(n,1));for(let a=0;a<i;a++)n[a]=a}}function $r(s){const e=tt(s),i=s.drawRange,t=i.start/3,n=(i.start+i.count)/3,a=Math.max(0,t),r=Math.min(e,n)-a;return[{offset:Math.floor(a),count:Math.floor(r)}]}function Xr(s){if(!s.groups||!s.groups.length)return $r(s);const e=[],i=new Set,t=s.drawRange,n=t.start/3,a=(t.start+t.count)/3;for(const c of s.groups){const l=c.start/3,h=(c.start+c.count)/3;i.add(Math.max(n,l)),i.add(Math.min(a,h))}const r=Array.from(i.values()).sort((c,l)=>c-l);for(let c=0;c<r.length-1;c++){const l=r[c],h=r[c+1];e.push({offset:Math.floor(l),count:Math.floor(h-l)})}return e}function Zs(s){if(s.groups.length===0)return!1;const e=tt(s),i=Xr(s).sort((a,r)=>a.offset-r.offset),t=i[i.length-1];t.count=Math.min(e-t.offset,t.count);let n=0;return i.forEach(({count:a})=>n+=a),e!==n}function Jt(s,e,i,t,n){let a=1/0,r=1/0,c=1/0,l=-1/0,h=-1/0,f=-1/0,u=1/0,o=1/0,m=1/0,v=-1/0,b=-1/0,d=-1/0;for(let p=e*6,g=(e+i)*6;p<g;p+=6){const y=s[p+0],w=s[p+1],T=y-w,_=y+w;T<a&&(a=T),_>l&&(l=_),y<u&&(u=y),y>v&&(v=y);const A=s[p+2],S=s[p+3],R=A-S,I=A+S;R<r&&(r=R),I>h&&(h=I),A<o&&(o=A),A>b&&(b=A);const P=s[p+4],M=s[p+5],F=P-M,E=P+M;F<c&&(c=F),E>f&&(f=E),P<m&&(m=P),P>d&&(d=P)}t[0]=a,t[1]=r,t[2]=c,t[3]=l,t[4]=h,t[5]=f,n[0]=u,n[1]=o,n[2]=m,n[3]=v,n[4]=b,n[5]=d}function Qs(s,e=null,i=null,t=null){const n=s.attributes.position,a=s.index?s.index.array:null,r=tt(s),c=n.normalized;let l;e===null?(l=new Float32Array(r*6*4),i=0,t=r):(l=e,i=i||0,t=t||r);const h=n.array,f=n.offset||0;let u=3;n.isInterleavedBufferAttribute&&(u=n.data.stride);const o=["getX","getY","getZ"];for(let m=i;m<i+t;m++){const v=m*3,b=m*6;let d=v+0,p=v+1,g=v+2;a&&(d=a[d],p=a[p],g=a[g]),c||(d=d*u+f,p=p*u+f,g=g*u+f);for(let y=0;y<3;y++){let w,T,_;c?(w=n[o[y]](d),T=n[o[y]](p),_=n[o[y]](g)):(w=h[d+y],T=h[p+y],_=h[g+y]);let A=w;T<A&&(A=T),_<A&&(A=_);let S=w;T>S&&(S=T),_>S&&(S=_);const R=(S-A)/2,I=y*2;l[b+I+0]=A+R,l[b+I+1]=R+(Math.abs(A)+R)*Xs}}return l}function H(s,e,i){return i.min.x=e[s],i.min.y=e[s+1],i.min.z=e[s+2],i.max.x=e[s+3],i.max.y=e[s+4],i.max.z=e[s+5],i}function nr(s){let e=-1,i=-1/0;for(let t=0;t<3;t++){const n=s[t+3]-s[t];n>i&&(i=n,e=t)}return e}function ar(s,e){e.set(s)}function or(s,e,i){let t,n;for(let a=0;a<3;a++){const r=a+3;t=s[a],n=e[a],i[a]=t<n?t:n,t=s[r],n=e[r],i[r]=t>n?t:n}}function Tt(s,e,i){for(let t=0;t<3;t++){const n=e[s+2*t],a=e[s+2*t+1],r=n-a,c=n+a;r<i[t]&&(i[t]=r),c>i[t+3]&&(i[t+3]=c)}}function nt(s){const e=s[3]-s[0],i=s[4]-s[1],t=s[5]-s[2];return 2*(e*i+i*t+t*e)}const ge=32,Js=(s,e)=>s.candidate-e.candidate,Te=new Array(ge).fill().map(()=>({count:0,bounds:new Float32Array(6),rightCacheBounds:new Float32Array(6),leftCacheBounds:new Float32Array(6),candidate:0})),wt=new Float32Array(6);function en(s,e,i,t,n,a){let r=-1,c=0;if(a===qr)r=nr(e),r!==-1&&(c=(e[r]+e[r+3])/2);else if(a===$s)r=nr(s),r!==-1&&(c=tn(i,t,n,r));else if(a===jr){const l=nt(s);let h=Zt*n;const f=t*6,u=(t+n)*6;for(let o=0;o<3;o++){const m=e[o],d=(e[o+3]-m)/ge;if(n<ge/4){const p=[...Te];p.length=n;let g=0;for(let w=f;w<u;w+=6,g++){const T=p[g];T.candidate=i[w+2*o],T.count=0;const{bounds:_,leftCacheBounds:A,rightCacheBounds:S}=T;for(let R=0;R<3;R++)S[R]=1/0,S[R+3]=-1/0,A[R]=1/0,A[R+3]=-1/0,_[R]=1/0,_[R+3]=-1/0;Tt(w,i,_)}p.sort(Js);let y=n;for(let w=0;w<y;w++){const T=p[w];for(;w+1<y&&p[w+1].candidate===T.candidate;)p.splice(w+1,1),y--}for(let w=f;w<u;w+=6){const T=i[w+2*o];for(let _=0;_<y;_++){const A=p[_];T>=A.candidate?Tt(w,i,A.rightCacheBounds):(Tt(w,i,A.leftCacheBounds),A.count++)}}for(let w=0;w<y;w++){const T=p[w],_=T.count,A=n-T.count,S=T.leftCacheBounds,R=T.rightCacheBounds;let I=0;_!==0&&(I=nt(S)/l);let P=0;A!==0&&(P=nt(R)/l);const M=sr+Zt*(I*_+P*A);M<h&&(r=o,h=M,c=T.candidate)}}else{for(let y=0;y<ge;y++){const w=Te[y];w.count=0,w.candidate=m+d+y*d;const T=w.bounds;for(let _=0;_<3;_++)T[_]=1/0,T[_+3]=-1/0}for(let y=f;y<u;y+=6){let _=~~((i[y+2*o]-m)/d);_>=ge&&(_=ge-1);const A=Te[_];A.count++,Tt(y,i,A.bounds)}const p=Te[ge-1];ar(p.bounds,p.rightCacheBounds);for(let y=ge-2;y>=0;y--){const w=Te[y],T=Te[y+1];or(w.bounds,T.rightCacheBounds,w.rightCacheBounds)}let g=0;for(let y=0;y<ge-1;y++){const w=Te[y],T=w.count,_=w.bounds,S=Te[y+1].rightCacheBounds;T!==0&&(g===0?ar(_,wt):or(_,wt,wt)),g+=T;let R=0,I=0;g!==0&&(R=nt(wt)/l);const P=n-g;P!==0&&(I=nt(S)/l);const M=sr+Zt*(R*g+I*P);M<h&&(r=o,h=M,c=w.candidate)}}}}else console.warn(`MeshBVH: Invalid build strategy value ${a} used.`);return{axis:r,pos:c}}function tn(s,e,i,t){let n=0;for(let a=e,r=e+i;a<r;a++)n+=s[a*6+t*2];return n/i}class ei{constructor(){this.boundingData=new Float32Array(6)}}function rn(s,e,i,t,n,a){let r=t,c=t+n-1;const l=a.pos,h=a.axis*2;for(;;){for(;r<=c&&i[r*6+h]<l;)r++;for(;r<=c&&i[c*6+h]>=l;)c--;if(r<c){for(let f=0;f<3;f++){let u=e[r*3+f];e[r*3+f]=e[c*3+f],e[c*3+f]=u}for(let f=0;f<6;f++){let u=i[r*6+f];i[r*6+f]=i[c*6+f],i[c*6+f]=u}r++,c--}else return r}}function sn(s,e,i,t,n,a){let r=t,c=t+n-1;const l=a.pos,h=a.axis*2;for(;;){for(;r<=c&&i[r*6+h]<l;)r++;for(;r<=c&&i[c*6+h]>=l;)c--;if(r<c){let f=s[r];s[r]=s[c],s[c]=f;for(let u=0;u<6;u++){let o=i[r*6+u];i[r*6+u]=i[c*6+u],i[c*6+u]=o}r++,c--}else return r}}function G(s,e){return e[s+15]===65535}function Z(s,e){return e[s+6]}function Q(s,e){return e[s+14]}function se(s){return s+8}function J(s,e){return e[s+6]}function Pi(s,e){return e[s+7]}let Kr,ht,Ut,Zr;const nn=Math.pow(2,32);function vi(s){return"count"in s?1:1+vi(s.left)+vi(s.right)}function an(s,e,i){return Kr=new Float32Array(i),ht=new Uint32Array(i),Ut=new Uint16Array(i),Zr=new Uint8Array(i),xi(s,e)}function xi(s,e){const i=s/4,t=s/2,n="count"in e,a=e.boundingData;for(let r=0;r<6;r++)Kr[i+r]=a[r];if(n)if(e.buffer){const r=e.buffer;Zr.set(new Uint8Array(r),s);for(let c=s,l=s+r.byteLength;c<l;c+=Se){const h=c/2;G(h,Ut)||(ht[c/4+6]+=i)}return s+r.byteLength}else{const r=e.offset,c=e.count;return ht[i+6]=r,Ut[t+14]=c,Ut[t+15]=jt,s+Se}else{const r=e.left,c=e.right,l=e.splitAxis;let h;if(h=xi(s+Se,r),h/4>nn)throw new Error("MeshBVH: Cannot store child pointer greater than 32 bits.");return ht[i+6]=h/4,h=xi(h,c),ht[i+7]=l,h}}function on(s,e){const i=(s.index?s.index.count:s.attributes.position.count)/3,t=i>2**16,n=t?4:2,a=e?new SharedArrayBuffer(i*n):new ArrayBuffer(i*n),r=t?new Uint32Array(a):new Uint16Array(a);for(let c=0,l=r.length;c<l;c++)r[c]=c;return r}function cn(s,e,i,t,n){const{maxDepth:a,verbose:r,maxLeafTris:c,strategy:l,onProgress:h,indirect:f}=n,u=s._indirectBuffer,o=s.geometry,m=o.index?o.index.array:null,v=f?sn:rn,b=tt(o),d=new Float32Array(6);let p=!1;const g=new ei;return Jt(e,i,t,g.boundingData,d),w(g,i,t,d),g;function y(T){h&&h(T/b)}function w(T,_,A,S=null,R=0){if(!p&&R>=a&&(p=!0,r&&(console.warn(`MeshBVH: Max depth of ${a} reached when generating BVH. Consider increasing maxDepth.`),console.warn(o))),A<=c||R>=a)return y(_+A),T.offset=_,T.count=A,T;const I=en(T.boundingData,S,e,_,A,l);if(I.axis===-1)return y(_+A),T.offset=_,T.count=A,T;const P=v(u,m,e,_,A,I);if(P===_||P===_+A)y(_+A),T.offset=_,T.count=A;else{T.splitAxis=I.axis;const M=new ei,F=_,E=P-_;T.left=M,Jt(e,F,E,M.boundingData,d),w(M,F,E,d,R+1);const B=new ei,z=P,X=A-E;T.right=B,Jt(e,z,X,B.boundingData,d),w(B,z,X,d,R+1)}return T}}function ln(s,e){const i=s.geometry;e.indirect&&(s._indirectBuffer=on(i,e.useSharedArrayBuffer),Zs(i)&&!e.verbose&&console.warn('MeshBVH: Provided geometry contains groups that do not fully span the vertex contents while using the "indirect" option. BVH may incorrectly report intersections on unrendered portions of the geometry.')),s._indirectBuffer||Ks(i,e);const t=e.useSharedArrayBuffer?SharedArrayBuffer:ArrayBuffer,n=Qs(i),a=e.indirect?$r(i):Xr(i);s._roots=a.map(r=>{const c=cn(s,n,r.offset,r.count,e),l=vi(c),h=new t(Se*l);return an(0,c,h),h})}class ye{constructor(){this.min=1/0,this.max=-1/0}setFromPointsField(e,i){let t=1/0,n=-1/0;for(let a=0,r=e.length;a<r;a++){const l=e[a][i];t=l<t?l:t,n=l>n?l:n}this.min=t,this.max=n}setFromPoints(e,i){let t=1/0,n=-1/0;for(let a=0,r=i.length;a<r;a++){const c=i[a],l=e.dot(c);t=l<t?l:t,n=l>n?l:n}this.min=t,this.max=n}isSeparated(e){return this.min>e.max||e.min>this.max}}ye.prototype.setFromBox=function(){const s=new C;return function(i,t){const n=t.min,a=t.max;let r=1/0,c=-1/0;for(let l=0;l<=1;l++)for(let h=0;h<=1;h++)for(let f=0;f<=1;f++){s.x=n.x*l+a.x*(1-l),s.y=n.y*h+a.y*(1-h),s.z=n.z*f+a.z*(1-f);const u=i.dot(s);r=Math.min(u,r),c=Math.max(u,c)}this.min=r,this.max=c}}();const un=function(){const s=new C,e=new C,i=new C;return function(n,a,r){const c=n.start,l=s,h=a.start,f=e;i.subVectors(c,h),s.subVectors(n.end,n.start),e.subVectors(a.end,a.start);const u=i.dot(f),o=f.dot(l),m=f.dot(f),v=i.dot(l),d=l.dot(l)*m-o*o;let p,g;d!==0?p=(u*o-v*m)/d:p=0,g=(u+p*o)/m,r.x=p,r.y=g}}(),Ci=function(){const s=new L,e=new C,i=new C;return function(n,a,r,c){un(n,a,s);let l=s.x,h=s.y;if(l>=0&&l<=1&&h>=0&&h<=1){n.at(l,r),a.at(h,c);return}else if(l>=0&&l<=1){h<0?a.at(0,c):a.at(1,c),n.closestPointToPoint(c,!0,r);return}else if(h>=0&&h<=1){l<0?n.at(0,r):n.at(1,r),a.closestPointToPoint(r,!0,c);return}else{let f;l<0?f=n.start:f=n.end;let u;h<0?u=a.start:u=a.end;const o=e,m=i;if(n.closestPointToPoint(u,!0,e),a.closestPointToPoint(f,!0,i),o.distanceToSquared(u)<=m.distanceToSquared(f)){r.copy(o),c.copy(u);return}else{r.copy(f),c.copy(m);return}}}}(),fn=function(){const s=new C,e=new C,i=new Ii,t=new xe;return function(a,r){const{radius:c,center:l}=a,{a:h,b:f,c:u}=r;if(t.start=h,t.end=f,t.closestPointToPoint(l,!0,s).distanceTo(l)<=c||(t.start=h,t.end=u,t.closestPointToPoint(l,!0,s).distanceTo(l)<=c)||(t.start=f,t.end=u,t.closestPointToPoint(l,!0,s).distanceTo(l)<=c))return!0;const b=r.getPlane(i);if(Math.abs(b.distanceToPoint(l))<=c){const p=b.projectPoint(l,e);if(r.containsPoint(p))return!0}return!1}}(),hn=1e-15;function ti(s){return Math.abs(s)<hn}class ue extends ft{constructor(...e){super(...e),this.isExtendedTriangle=!0,this.satAxes=new Array(4).fill().map(()=>new C),this.satBounds=new Array(4).fill().map(()=>new ye),this.points=[this.a,this.b,this.c],this.sphere=new Cs,this.plane=new Ii,this.needsUpdate=!0}intersectsSphere(e){return fn(e,this)}update(){const e=this.a,i=this.b,t=this.c,n=this.points,a=this.satAxes,r=this.satBounds,c=a[0],l=r[0];this.getNormal(c),l.setFromPoints(c,n);const h=a[1],f=r[1];h.subVectors(e,i),f.setFromPoints(h,n);const u=a[2],o=r[2];u.subVectors(i,t),o.setFromPoints(u,n);const m=a[3],v=r[3];m.subVectors(t,e),v.setFromPoints(m,n),this.sphere.setFromPoints(this.points),this.plane.setFromNormalAndCoplanarPoint(c,e),this.needsUpdate=!1}}ue.prototype.closestPointToSegment=function(){const s=new C,e=new C,i=new xe;return function(n,a=null,r=null){const{start:c,end:l}=n,h=this.points;let f,u=1/0;for(let o=0;o<3;o++){const m=(o+1)%3;i.start.copy(h[o]),i.end.copy(h[m]),Ci(i,n,s,e),f=s.distanceToSquared(e),f<u&&(u=f,a&&a.copy(s),r&&r.copy(e))}return this.closestPointToPoint(c,s),f=c.distanceToSquared(s),f<u&&(u=f,a&&a.copy(s),r&&r.copy(c)),this.closestPointToPoint(l,s),f=l.distanceToSquared(s),f<u&&(u=f,a&&a.copy(s),r&&r.copy(l)),Math.sqrt(u)}}();ue.prototype.intersectsTriangle=function(){const s=new ue,e=new Array(3),i=new Array(3),t=new ye,n=new ye,a=new C,r=new C,c=new C,l=new C,h=new C,f=new xe,u=new xe,o=new xe,m=new C;function v(b,d,p){const g=b.points;let y=0,w=-1;for(let T=0;T<3;T++){const{start:_,end:A}=f;_.copy(g[T]),A.copy(g[(T+1)%3]),f.delta(r);const S=ti(d.distanceToPoint(_));if(ti(d.normal.dot(r))&&S){p.copy(f),y=2;break}const R=d.intersectLine(f,m);if(!R&&S&&m.copy(_),(R||S)&&!ti(m.distanceTo(A))){if(y<=1)(y===1?p.start:p.end).copy(m),S&&(w=y);else if(y>=2){(w===1?p.start:p.end).copy(m),y=2;break}if(y++,y===2&&w===-1)break}}return y}return function(d,p=null,g=!1){this.needsUpdate&&this.update(),d.isExtendedTriangle?d.needsUpdate&&d.update():(s.copy(d),s.update(),d=s);const y=this.plane,w=d.plane;if(Math.abs(y.normal.dot(w.normal))>1-1e-10){const T=this.satBounds,_=this.satAxes;i[0]=d.a,i[1]=d.b,i[2]=d.c;for(let R=0;R<4;R++){const I=T[R],P=_[R];if(t.setFromPoints(P,i),I.isSeparated(t))return!1}const A=d.satBounds,S=d.satAxes;e[0]=this.a,e[1]=this.b,e[2]=this.c;for(let R=0;R<4;R++){const I=A[R],P=S[R];if(t.setFromPoints(P,e),I.isSeparated(t))return!1}for(let R=0;R<4;R++){const I=_[R];for(let P=0;P<4;P++){const M=S[P];if(a.crossVectors(I,M),t.setFromPoints(a,e),n.setFromPoints(a,i),t.isSeparated(n))return!1}}return p&&(g||console.warn("ExtendedTriangle.intersectsTriangle: Triangles are coplanar which does not support an output edge. Setting edge to 0, 0, 0."),p.start.set(0,0,0),p.end.set(0,0,0)),!0}else{const T=v(this,w,u);if(T===1&&d.containsPoint(u.end))return p&&(p.start.copy(u.end),p.end.copy(u.end)),!0;if(T!==2)return!1;const _=v(d,y,o);if(_===1&&this.containsPoint(o.end))return p&&(p.start.copy(o.end),p.end.copy(o.end)),!0;if(_!==2)return!1;if(u.delta(c),o.delta(l),c.dot(l)<0){let F=o.start;o.start=o.end,o.end=F}const A=u.start.dot(c),S=u.end.dot(c),R=o.start.dot(c),I=o.end.dot(c),P=S<R,M=A<I;return A!==I&&R!==S&&P===M?!1:(p&&(h.subVectors(u.start,o.start),h.dot(c)>0?p.start.copy(u.start):p.start.copy(o.start),h.subVectors(u.end,o.end),h.dot(c)<0?p.end.copy(u.end):p.end.copy(o.end)),!0)}}}();ue.prototype.distanceToPoint=function(){const s=new C;return function(i){return this.closestPointToPoint(i,s),i.distanceTo(s)}}();ue.prototype.distanceToTriangle=function(){const s=new C,e=new C,i=["a","b","c"],t=new xe,n=new xe;return function(r,c=null,l=null){const h=c||l?t:null;if(this.intersectsTriangle(r,h))return(c||l)&&(c&&h.getCenter(c),l&&h.getCenter(l)),0;let f=1/0;for(let u=0;u<3;u++){let o;const m=i[u],v=r[m];this.closestPointToPoint(v,s),o=v.distanceToSquared(s),o<f&&(f=o,c&&c.copy(s),l&&l.copy(v));const b=this[m];r.closestPointToPoint(b,s),o=b.distanceToSquared(s),o<f&&(f=o,c&&c.copy(b),l&&l.copy(s))}for(let u=0;u<3;u++){const o=i[u],m=i[(u+1)%3];t.set(this[o],this[m]);for(let v=0;v<3;v++){const b=i[v],d=i[(v+1)%3];n.set(r[b],r[d]),Ci(t,n,s,e);const p=s.distanceToSquared(e);p<f&&(f=p,c&&c.copy(s),l&&l.copy(e))}}return Math.sqrt(f)}}();class ${constructor(e,i,t){this.isOrientedBox=!0,this.min=new C,this.max=new C,this.matrix=new Y,this.invMatrix=new Y,this.points=new Array(8).fill().map(()=>new C),this.satAxes=new Array(3).fill().map(()=>new C),this.satBounds=new Array(3).fill().map(()=>new ye),this.alignedSatBounds=new Array(3).fill().map(()=>new ye),this.needsUpdate=!1,e&&this.min.copy(e),i&&this.max.copy(i),t&&this.matrix.copy(t)}set(e,i,t){this.min.copy(e),this.max.copy(i),this.matrix.copy(t),this.needsUpdate=!0}copy(e){this.min.copy(e.min),this.max.copy(e.max),this.matrix.copy(e.matrix),this.needsUpdate=!0}}$.prototype.update=function(){return function(){const e=this.matrix,i=this.min,t=this.max,n=this.points;for(let h=0;h<=1;h++)for(let f=0;f<=1;f++)for(let u=0;u<=1;u++){const o=1*h|2*f|4*u,m=n[o];m.x=h?t.x:i.x,m.y=f?t.y:i.y,m.z=u?t.z:i.z,m.applyMatrix4(e)}const a=this.satBounds,r=this.satAxes,c=n[0];for(let h=0;h<3;h++){const f=r[h],u=a[h],o=1<<h,m=n[o];f.subVectors(c,m),u.setFromPoints(f,n)}const l=this.alignedSatBounds;l[0].setFromPointsField(n,"x"),l[1].setFromPointsField(n,"y"),l[2].setFromPointsField(n,"z"),this.invMatrix.copy(this.matrix).invert(),this.needsUpdate=!1}}();$.prototype.intersectsBox=function(){const s=new ye;return function(i){this.needsUpdate&&this.update();const t=i.min,n=i.max,a=this.satBounds,r=this.satAxes,c=this.alignedSatBounds;if(s.min=t.x,s.max=n.x,c[0].isSeparated(s)||(s.min=t.y,s.max=n.y,c[1].isSeparated(s))||(s.min=t.z,s.max=n.z,c[2].isSeparated(s)))return!1;for(let l=0;l<3;l++){const h=r[l],f=a[l];if(s.setFromBox(h,i),f.isSeparated(s))return!1}return!0}}();$.prototype.intersectsTriangle=function(){const s=new ue,e=new Array(3),i=new ye,t=new ye,n=new C;return function(r){this.needsUpdate&&this.update(),r.isExtendedTriangle?r.needsUpdate&&r.update():(s.copy(r),s.update(),r=s);const c=this.satBounds,l=this.satAxes;e[0]=r.a,e[1]=r.b,e[2]=r.c;for(let o=0;o<3;o++){const m=c[o],v=l[o];if(i.setFromPoints(v,e),m.isSeparated(i))return!1}const h=r.satBounds,f=r.satAxes,u=this.points;for(let o=0;o<3;o++){const m=h[o],v=f[o];if(i.setFromPoints(v,u),m.isSeparated(i))return!1}for(let o=0;o<3;o++){const m=l[o];for(let v=0;v<4;v++){const b=f[v];if(n.crossVectors(m,b),i.setFromPoints(n,e),t.setFromPoints(n,u),i.isSeparated(t))return!1}}return!0}}();$.prototype.closestPointToPoint=function(){return function(e,i){return this.needsUpdate&&this.update(),i.copy(e).applyMatrix4(this.invMatrix).clamp(this.min,this.max).applyMatrix4(this.matrix),i}}();$.prototype.distanceToPoint=function(){const s=new C;return function(i){return this.closestPointToPoint(i,s),i.distanceTo(s)}}();$.prototype.distanceToBox=function(){const s=["x","y","z"],e=new Array(12).fill().map(()=>new xe),i=new Array(12).fill().map(()=>new xe),t=new C,n=new C;return function(r,c=0,l=null,h=null){if(this.needsUpdate&&this.update(),this.intersectsBox(r))return(l||h)&&(r.getCenter(n),this.closestPointToPoint(n,t),r.closestPointToPoint(t,n),l&&l.copy(t),h&&h.copy(n)),0;const f=c*c,u=r.min,o=r.max,m=this.points;let v=1/0;for(let d=0;d<8;d++){const p=m[d];n.copy(p).clamp(u,o);const g=p.distanceToSquared(n);if(g<v&&(v=g,l&&l.copy(p),h&&h.copy(n),g<f))return Math.sqrt(g)}let b=0;for(let d=0;d<3;d++)for(let p=0;p<=1;p++)for(let g=0;g<=1;g++){const y=(d+1)%3,w=(d+2)%3,T=p<<y|g<<w,_=1<<d|p<<y|g<<w,A=m[T],S=m[_];e[b].set(A,S);const I=s[d],P=s[y],M=s[w],F=i[b],E=F.start,B=F.end;E[I]=u[I],E[P]=p?u[P]:o[P],E[M]=g?u[M]:o[P],B[I]=o[I],B[P]=p?u[P]:o[P],B[M]=g?u[M]:o[P],b++}for(let d=0;d<=1;d++)for(let p=0;p<=1;p++)for(let g=0;g<=1;g++){n.x=d?o.x:u.x,n.y=p?o.y:u.y,n.z=g?o.z:u.z,this.closestPointToPoint(n,t);const y=n.distanceToSquared(t);if(y<v&&(v=y,l&&l.copy(t),h&&h.copy(n),y<f))return Math.sqrt(y)}for(let d=0;d<12;d++){const p=e[d];for(let g=0;g<12;g++){const y=i[g];Ci(p,y,t,n);const w=t.distanceToSquared(n);if(w<v&&(v=w,l&&l.copy(t),h&&h.copy(n),w<f))return Math.sqrt(w)}}return Math.sqrt(v)}}();class Fi{constructor(e){this._getNewPrimitive=e,this._primitives=[]}getPrimitive(){const e=this._primitives;return e.length===0?this._getNewPrimitive():e.pop()}releasePrimitive(e){this._primitives.push(e)}}class dn extends Fi{constructor(){super(()=>new ue)}}const ne=new dn;class mn{constructor(){this.float32Array=null,this.uint16Array=null,this.uint32Array=null;const e=[];let i=null;this.setBuffer=t=>{i&&e.push(i),i=t,this.float32Array=new Float32Array(t),this.uint16Array=new Uint16Array(t),this.uint32Array=new Uint32Array(t)},this.clearBuffer=()=>{i=null,this.float32Array=null,this.uint16Array=null,this.uint32Array=null,e.length!==0&&this.setBuffer(e.pop())}}}const N=new mn;let Ae,Qe;const Ne=[],At=new Fi(()=>new he);function pn(s,e,i,t,n,a){Ae=At.getPrimitive(),Qe=At.getPrimitive(),Ne.push(Ae,Qe),N.setBuffer(s._roots[e]);const r=yi(0,s.geometry,i,t,n,a);N.clearBuffer(),At.releasePrimitive(Ae),At.releasePrimitive(Qe),Ne.pop(),Ne.pop();const c=Ne.length;return c>0&&(Qe=Ne[c-1],Ae=Ne[c-2]),r}function yi(s,e,i,t,n=null,a=0,r=0){const{float32Array:c,uint16Array:l,uint32Array:h}=N;let f=s*2;if(G(f,l)){const o=Z(s,h),m=Q(f,l);return H(s,c,Ae),t(o,m,!1,r,a+s,Ae)}else{let I=function(M){const{uint16Array:F,uint32Array:E}=N;let B=M*2;for(;!G(B,F);)M=se(M),B=M*2;return Z(M,E)},P=function(M){const{uint16Array:F,uint32Array:E}=N;let B=M*2;for(;!G(B,F);)M=J(M,E),B=M*2;return Z(M,E)+Q(B,F)};const o=se(s),m=J(s,h);let v=o,b=m,d,p,g,y;if(n&&(g=Ae,y=Qe,H(v,c,g),H(b,c,y),d=n(g),p=n(y),p<d)){v=m,b=o;const M=d;d=p,p=M,g=y}g||(g=Ae,H(v,c,g));const w=G(v*2,l),T=i(g,w,d,r+1,a+v);let _;if(T===rr){const M=I(v),E=P(v)-M;_=t(M,E,!0,r+1,a+v,g)}else _=T&&yi(v,e,i,t,n,a,r+1);if(_)return!0;y=Qe,H(b,c,y);const A=G(b*2,l),S=i(y,A,p,r+1,a+b);let R;if(S===rr){const M=I(b),E=P(b)-M;R=t(M,E,!0,r+1,a+b,y)}else R=S&&yi(b,e,i,t,n,a,r+1);return!!R}}const at=new C,ii=new C;function gn(s,e,i={},t=0,n=1/0){const a=t*t,r=n*n;let c=1/0,l=null;if(s.shapecast({boundsTraverseOrder:f=>(at.copy(e).clamp(f.min,f.max),at.distanceToSquared(e)),intersectsBounds:(f,u,o)=>o<c&&o<r,intersectsTriangle:(f,u)=>{f.closestPointToPoint(e,at);const o=e.distanceToSquared(at);return o<c&&(ii.copy(at),c=o,l=u),o<a}}),c===1/0)return null;const h=Math.sqrt(c);return i.point?i.point.copy(ii):i.point=ii.clone(),i.distance=h,i.faceIndex=l,i}const Le=new C,He=new C,Ue=new C,_t=new L,St=new L,Mt=new L,cr=new C,lr=new C,ur=new C,Rt=new C;function vn(s,e,i,t,n,a,r,c){let l;if(a===Or?l=s.intersectTriangle(t,i,e,!0,n):l=s.intersectTriangle(e,i,t,a!==Nr,n),l===null)return null;const h=s.origin.distanceTo(n);return h<r||h>c?null:{distance:h,point:n.clone()}}function xn(s,e,i,t,n,a,r,c,l,h,f){Le.fromBufferAttribute(e,a),He.fromBufferAttribute(e,r),Ue.fromBufferAttribute(e,c);const u=vn(s,Le,He,Ue,Rt,l,h,f);if(u){t&&(_t.fromBufferAttribute(t,a),St.fromBufferAttribute(t,r),Mt.fromBufferAttribute(t,c),u.uv=ft.getInterpolation(Rt,Le,He,Ue,_t,St,Mt,new L)),n&&(_t.fromBufferAttribute(n,a),St.fromBufferAttribute(n,r),Mt.fromBufferAttribute(n,c),u.uv1=ft.getInterpolation(Rt,Le,He,Ue,_t,St,Mt,new L)),i&&(cr.fromBufferAttribute(i,a),lr.fromBufferAttribute(i,r),ur.fromBufferAttribute(i,c),u.normal=ft.getInterpolation(Rt,Le,He,Ue,cr,lr,ur,new C),u.normal.dot(s.direction)>0&&u.normal.multiplyScalar(-1));const o={a,b:r,c,normal:new C,materialIndex:0};ft.getNormal(Le,He,Ue,o.normal),u.face=o,u.faceIndex=a}return u}function Gt(s,e,i,t,n,a,r){const c=t*3;let l=c+0,h=c+1,f=c+2;const u=s.index;s.index&&(l=u.getX(l),h=u.getX(h),f=u.getX(f));const{position:o,normal:m,uv:v,uv1:b}=s.attributes,d=xn(i,o,m,v,b,l,h,f,e,a,r);return d?(d.faceIndex=t,n&&n.push(d),d):null}function W(s,e,i,t){const n=s.a,a=s.b,r=s.c;let c=e,l=e+1,h=e+2;i&&(c=i.getX(c),l=i.getX(l),h=i.getX(h)),n.x=t.getX(c),n.y=t.getY(c),n.z=t.getZ(c),a.x=t.getX(l),a.y=t.getY(l),a.z=t.getZ(l),r.x=t.getX(h),r.y=t.getY(h),r.z=t.getZ(h)}function yn(s,e,i,t,n,a,r,c){const{geometry:l,_indirectBuffer:h}=s;for(let f=t,u=t+n;f<u;f++)Gt(l,e,i,f,a,r,c)}function bn(s,e,i,t,n,a,r){const{geometry:c,_indirectBuffer:l}=s;let h=1/0,f=null;for(let u=t,o=t+n;u<o;u++){let m;m=Gt(c,e,i,u,null,a,r),m&&m.distance<h&&(f=m,h=m.distance)}return f}function Tn(s,e,i,t,n,a,r){const{geometry:c}=i,{index:l}=c,h=c.attributes.position;for(let f=s,u=e+s;f<u;f++){let o;if(o=f,W(r,o*3,l,h),r.needsUpdate=!0,t(r,o,n,a))return!0}return!1}function wn(s,e=null){e&&Array.isArray(e)&&(e=new Set(e));const i=s.geometry,t=i.index?i.index.array:null,n=i.attributes.position;let a,r,c,l,h=0;const f=s._roots;for(let o=0,m=f.length;o<m;o++)a=f[o],r=new Uint32Array(a),c=new Uint16Array(a),l=new Float32Array(a),u(0,h),h+=a.byteLength;function u(o,m,v=!1){const b=o*2;if(c[b+15]===jt){const p=r[o+6],g=c[b+14];let y=1/0,w=1/0,T=1/0,_=-1/0,A=-1/0,S=-1/0;for(let R=3*p,I=3*(p+g);R<I;R++){let P=t[R];const M=n.getX(P),F=n.getY(P),E=n.getZ(P);M<y&&(y=M),M>_&&(_=M),F<w&&(w=F),F>A&&(A=F),E<T&&(T=E),E>S&&(S=E)}return l[o+0]!==y||l[o+1]!==w||l[o+2]!==T||l[o+3]!==_||l[o+4]!==A||l[o+5]!==S?(l[o+0]=y,l[o+1]=w,l[o+2]=T,l[o+3]=_,l[o+4]=A,l[o+5]=S,!0):!1}else{const p=o+8,g=r[o+6],y=p+m,w=g+m;let T=v,_=!1,A=!1;e?T||(_=e.has(y),A=e.has(w),T=!_&&!A):(_=!0,A=!0);const S=T||_,R=T||A;let I=!1;S&&(I=u(p,m,T));let P=!1;R&&(P=u(g,m,T));const M=I||P;if(M)for(let F=0;F<3;F++){const E=p+F,B=g+F,z=l[E],X=l[E+3],ee=l[B],te=l[B+3];l[o+F]=z<ee?z:ee,l[o+F+3]=X>te?X:te}return M}}}function Re(s,e,i,t,n){let a,r,c,l,h,f;const u=1/i.direction.x,o=1/i.direction.y,m=1/i.direction.z,v=i.origin.x,b=i.origin.y,d=i.origin.z;let p=e[s],g=e[s+3],y=e[s+1],w=e[s+3+1],T=e[s+2],_=e[s+3+2];return u>=0?(a=(p-v)*u,r=(g-v)*u):(a=(g-v)*u,r=(p-v)*u),o>=0?(c=(y-b)*o,l=(w-b)*o):(c=(w-b)*o,l=(y-b)*o),a>l||c>r||((c>a||isNaN(a))&&(a=c),(l<r||isNaN(r))&&(r=l),m>=0?(h=(T-d)*m,f=(_-d)*m):(h=(_-d)*m,f=(T-d)*m),a>f||h>r)?!1:((h>a||a!==a)&&(a=h),(f<r||r!==r)&&(r=f),a<=n&&r>=t)}function An(s,e,i,t,n,a,r,c){const{geometry:l,_indirectBuffer:h}=s;for(let f=t,u=t+n;f<u;f++){let o=h?h[f]:f;Gt(l,e,i,o,a,r,c)}}function _n(s,e,i,t,n,a,r){const{geometry:c,_indirectBuffer:l}=s;let h=1/0,f=null;for(let u=t,o=t+n;u<o;u++){let m;m=Gt(c,e,i,l?l[u]:u,null,a,r),m&&m.distance<h&&(f=m,h=m.distance)}return f}function Sn(s,e,i,t,n,a,r){const{geometry:c}=i,{index:l}=c,h=c.attributes.position;for(let f=s,u=e+s;f<u;f++){let o;if(o=i.resolveTriangleIndex(f),W(r,o*3,l,h),r.needsUpdate=!0,t(r,o,n,a))return!0}return!1}function Mn(s,e,i,t,n,a,r){N.setBuffer(s._roots[e]),bi(0,s,i,t,n,a,r),N.clearBuffer()}function bi(s,e,i,t,n,a,r){const{float32Array:c,uint16Array:l,uint32Array:h}=N,f=s*2;if(G(f,l)){const o=Z(s,h),m=Q(f,l);yn(e,i,t,o,m,n,a,r)}else{const o=se(s);Re(o,c,t,a,r)&&bi(o,e,i,t,n,a,r);const m=J(s,h);Re(m,c,t,a,r)&&bi(m,e,i,t,n,a,r)}}const Rn=["x","y","z"];function In(s,e,i,t,n,a){N.setBuffer(s._roots[e]);const r=Ti(0,s,i,t,n,a);return N.clearBuffer(),r}function Ti(s,e,i,t,n,a){const{float32Array:r,uint16Array:c,uint32Array:l}=N;let h=s*2;if(G(h,c)){const u=Z(s,l),o=Q(h,c);return bn(e,i,t,u,o,n,a)}else{const u=Pi(s,l),o=Rn[u],v=t.direction[o]>=0;let b,d;v?(b=se(s),d=J(s,l)):(b=J(s,l),d=se(s));const g=Re(b,r,t,n,a)?Ti(b,e,i,t,n,a):null;if(g){const T=g.point[o];if(v?T<=r[d+u]:T>=r[d+u+3])return g}const w=Re(d,r,t,n,a)?Ti(d,e,i,t,n,a):null;return g&&w?g.distance<=w.distance?g:w:g||w||null}}const It=new he,We=new ue,Ve=new ue,ot=new Y,fr=new $,Pt=new $;function Pn(s,e,i,t){N.setBuffer(s._roots[e]);const n=wi(0,s,i,t);return N.clearBuffer(),n}function wi(s,e,i,t,n=null){const{float32Array:a,uint16Array:r,uint32Array:c}=N;let l=s*2;if(n===null&&(i.boundingBox||i.computeBoundingBox(),fr.set(i.boundingBox.min,i.boundingBox.max,t),n=fr),G(l,r)){const f=e.geometry,u=f.index,o=f.attributes.position,m=i.index,v=i.attributes.position,b=Z(s,c),d=Q(l,r);if(ot.copy(t).invert(),i.boundsTree)return H(s,a,Pt),Pt.matrix.copy(ot),Pt.needsUpdate=!0,i.boundsTree.shapecast({intersectsBounds:g=>Pt.intersectsBox(g),intersectsTriangle:g=>{g.a.applyMatrix4(t),g.b.applyMatrix4(t),g.c.applyMatrix4(t),g.needsUpdate=!0;for(let y=b*3,w=(d+b)*3;y<w;y+=3)if(W(Ve,y,u,o),Ve.needsUpdate=!0,g.intersectsTriangle(Ve))return!0;return!1}});for(let p=b*3,g=(d+b)*3;p<g;p+=3){W(We,p,u,o),We.a.applyMatrix4(ot),We.b.applyMatrix4(ot),We.c.applyMatrix4(ot),We.needsUpdate=!0;for(let y=0,w=m.count;y<w;y+=3)if(W(Ve,y,m,v),Ve.needsUpdate=!0,We.intersectsTriangle(Ve))return!0}}else{const f=s+8,u=c[s+6];return H(f,a,It),!!(n.intersectsBox(It)&&wi(f,e,i,t,n)||(H(u,a,It),n.intersectsBox(It)&&wi(u,e,i,t,n)))}}const Ct=new Y,ri=new $,ct=new $,Cn=new C,Fn=new C,Dn=new C,En=new C;function Bn(s,e,i,t={},n={},a=0,r=1/0){e.boundingBox||e.computeBoundingBox(),ri.set(e.boundingBox.min,e.boundingBox.max,i),ri.needsUpdate=!0;const c=s.geometry,l=c.attributes.position,h=c.index,f=e.attributes.position,u=e.index,o=ne.getPrimitive(),m=ne.getPrimitive();let v=Cn,b=Fn,d=null,p=null;n&&(d=Dn,p=En);let g=1/0,y=null,w=null;return Ct.copy(i).invert(),ct.matrix.copy(Ct),s.shapecast({boundsTraverseOrder:T=>ri.distanceToBox(T),intersectsBounds:(T,_,A)=>A<g&&A<r?(_&&(ct.min.copy(T.min),ct.max.copy(T.max),ct.needsUpdate=!0),!0):!1,intersectsRange:(T,_)=>{if(e.boundsTree)return e.boundsTree.shapecast({boundsTraverseOrder:S=>ct.distanceToBox(S),intersectsBounds:(S,R,I)=>I<g&&I<r,intersectsRange:(S,R)=>{for(let I=S,P=S+R;I<P;I++){W(m,3*I,u,f),m.a.applyMatrix4(i),m.b.applyMatrix4(i),m.c.applyMatrix4(i),m.needsUpdate=!0;for(let M=T,F=T+_;M<F;M++){W(o,3*M,h,l),o.needsUpdate=!0;const E=o.distanceToTriangle(m,v,d);if(E<g&&(b.copy(v),p&&p.copy(d),g=E,y=M,w=I),E<a)return!0}}}});{const A=tt(e);for(let S=0,R=A;S<R;S++){W(m,3*S,u,f),m.a.applyMatrix4(i),m.b.applyMatrix4(i),m.c.applyMatrix4(i),m.needsUpdate=!0;for(let I=T,P=T+_;I<P;I++){W(o,3*I,h,l),o.needsUpdate=!0;const M=o.distanceToTriangle(m,v,d);if(M<g&&(b.copy(v),p&&p.copy(d),g=M,y=I,w=S),M<a)return!0}}}}}),ne.releasePrimitive(o),ne.releasePrimitive(m),g===1/0?null:(t.point?t.point.copy(b):t.point=b.clone(),t.distance=g,t.faceIndex=y,n&&(n.point?n.point.copy(p):n.point=p.clone(),n.point.applyMatrix4(Ct),b.applyMatrix4(Ct),n.distance=b.sub(n.point).length(),n.faceIndex=w),t)}function kn(s,e=null){e&&Array.isArray(e)&&(e=new Set(e));const i=s.geometry,t=i.index?i.index.array:null,n=i.attributes.position;let a,r,c,l,h=0;const f=s._roots;for(let o=0,m=f.length;o<m;o++)a=f[o],r=new Uint32Array(a),c=new Uint16Array(a),l=new Float32Array(a),u(0,h),h+=a.byteLength;function u(o,m,v=!1){const b=o*2;if(c[b+15]===jt){const p=r[o+6],g=c[b+14];let y=1/0,w=1/0,T=1/0,_=-1/0,A=-1/0,S=-1/0;for(let R=p,I=p+g;R<I;R++){const P=3*s.resolveTriangleIndex(R);for(let M=0;M<3;M++){let F=P+M;F=t?t[F]:F;const E=n.getX(F),B=n.getY(F),z=n.getZ(F);E<y&&(y=E),E>_&&(_=E),B<w&&(w=B),B>A&&(A=B),z<T&&(T=z),z>S&&(S=z)}}return l[o+0]!==y||l[o+1]!==w||l[o+2]!==T||l[o+3]!==_||l[o+4]!==A||l[o+5]!==S?(l[o+0]=y,l[o+1]=w,l[o+2]=T,l[o+3]=_,l[o+4]=A,l[o+5]=S,!0):!1}else{const p=o+8,g=r[o+6],y=p+m,w=g+m;let T=v,_=!1,A=!1;e?T||(_=e.has(y),A=e.has(w),T=!_&&!A):(_=!0,A=!0);const S=T||_,R=T||A;let I=!1;S&&(I=u(p,m,T));let P=!1;R&&(P=u(g,m,T));const M=I||P;if(M)for(let F=0;F<3;F++){const E=p+F,B=g+F,z=l[E],X=l[E+3],ee=l[B],te=l[B+3];l[o+F]=z<ee?z:ee,l[o+F+3]=X>te?X:te}return M}}}function zn(s,e,i,t,n,a,r){N.setBuffer(s._roots[e]),Ai(0,s,i,t,n,a,r),N.clearBuffer()}function Ai(s,e,i,t,n,a,r){const{float32Array:c,uint16Array:l,uint32Array:h}=N,f=s*2;if(G(f,l)){const o=Z(s,h),m=Q(f,l);An(e,i,t,o,m,n,a,r)}else{const o=se(s);Re(o,c,t,a,r)&&Ai(o,e,i,t,n,a,r);const m=J(s,h);Re(m,c,t,a,r)&&Ai(m,e,i,t,n,a,r)}}const On=["x","y","z"];function Nn(s,e,i,t,n,a){N.setBuffer(s._roots[e]);const r=_i(0,s,i,t,n,a);return N.clearBuffer(),r}function _i(s,e,i,t,n,a){const{float32Array:r,uint16Array:c,uint32Array:l}=N;let h=s*2;if(G(h,c)){const u=Z(s,l),o=Q(h,c);return _n(e,i,t,u,o,n,a)}else{const u=Pi(s,l),o=On[u],v=t.direction[o]>=0;let b,d;v?(b=se(s),d=J(s,l)):(b=J(s,l),d=se(s));const g=Re(b,r,t,n,a)?_i(b,e,i,t,n,a):null;if(g){const T=g.point[o];if(v?T<=r[d+u]:T>=r[d+u+3])return g}const w=Re(d,r,t,n,a)?_i(d,e,i,t,n,a):null;return g&&w?g.distance<=w.distance?g:w:g||w||null}}const Ft=new he,qe=new ue,je=new ue,lt=new Y,hr=new $,Dt=new $;function Ln(s,e,i,t){N.setBuffer(s._roots[e]);const n=Si(0,s,i,t);return N.clearBuffer(),n}function Si(s,e,i,t,n=null){const{float32Array:a,uint16Array:r,uint32Array:c}=N;let l=s*2;if(n===null&&(i.boundingBox||i.computeBoundingBox(),hr.set(i.boundingBox.min,i.boundingBox.max,t),n=hr),G(l,r)){const f=e.geometry,u=f.index,o=f.attributes.position,m=i.index,v=i.attributes.position,b=Z(s,c),d=Q(l,r);if(lt.copy(t).invert(),i.boundsTree)return H(s,a,Dt),Dt.matrix.copy(lt),Dt.needsUpdate=!0,i.boundsTree.shapecast({intersectsBounds:g=>Dt.intersectsBox(g),intersectsTriangle:g=>{g.a.applyMatrix4(t),g.b.applyMatrix4(t),g.c.applyMatrix4(t),g.needsUpdate=!0;for(let y=b,w=d+b;y<w;y++)if(W(je,3*e.resolveTriangleIndex(y),u,o),je.needsUpdate=!0,g.intersectsTriangle(je))return!0;return!1}});for(let p=b,g=d+b;p<g;p++){const y=e.resolveTriangleIndex(p);W(qe,3*y,u,o),qe.a.applyMatrix4(lt),qe.b.applyMatrix4(lt),qe.c.applyMatrix4(lt),qe.needsUpdate=!0;for(let w=0,T=m.count;w<T;w+=3)if(W(je,w,m,v),je.needsUpdate=!0,qe.intersectsTriangle(je))return!0}}else{const f=s+8,u=c[s+6];return H(f,a,Ft),!!(n.intersectsBox(Ft)&&Si(f,e,i,t,n)||(H(u,a,Ft),n.intersectsBox(Ft)&&Si(u,e,i,t,n)))}}const Et=new Y,si=new $,ut=new $,Hn=new C,Un=new C,Wn=new C,Vn=new C;function qn(s,e,i,t={},n={},a=0,r=1/0){e.boundingBox||e.computeBoundingBox(),si.set(e.boundingBox.min,e.boundingBox.max,i),si.needsUpdate=!0;const c=s.geometry,l=c.attributes.position,h=c.index,f=e.attributes.position,u=e.index,o=ne.getPrimitive(),m=ne.getPrimitive();let v=Hn,b=Un,d=null,p=null;n&&(d=Wn,p=Vn);let g=1/0,y=null,w=null;return Et.copy(i).invert(),ut.matrix.copy(Et),s.shapecast({boundsTraverseOrder:T=>si.distanceToBox(T),intersectsBounds:(T,_,A)=>A<g&&A<r?(_&&(ut.min.copy(T.min),ut.max.copy(T.max),ut.needsUpdate=!0),!0):!1,intersectsRange:(T,_)=>{if(e.boundsTree){const A=e.boundsTree;return A.shapecast({boundsTraverseOrder:S=>ut.distanceToBox(S),intersectsBounds:(S,R,I)=>I<g&&I<r,intersectsRange:(S,R)=>{for(let I=S,P=S+R;I<P;I++){const M=A.resolveTriangleIndex(I);W(m,3*M,u,f),m.a.applyMatrix4(i),m.b.applyMatrix4(i),m.c.applyMatrix4(i),m.needsUpdate=!0;for(let F=T,E=T+_;F<E;F++){const B=s.resolveTriangleIndex(F);W(o,3*B,h,l),o.needsUpdate=!0;const z=o.distanceToTriangle(m,v,d);if(z<g&&(b.copy(v),p&&p.copy(d),g=z,y=F,w=I),z<a)return!0}}}})}else{const A=tt(e);for(let S=0,R=A;S<R;S++){W(m,3*S,u,f),m.a.applyMatrix4(i),m.b.applyMatrix4(i),m.c.applyMatrix4(i),m.needsUpdate=!0;for(let I=T,P=T+_;I<P;I++){const M=s.resolveTriangleIndex(I);W(o,3*M,h,l),o.needsUpdate=!0;const F=o.distanceToTriangle(m,v,d);if(F<g&&(b.copy(v),p&&p.copy(d),g=F,y=I,w=S),F<a)return!0}}}}}),ne.releasePrimitive(o),ne.releasePrimitive(m),g===1/0?null:(t.point?t.point.copy(b):t.point=b.clone(),t.distance=g,t.faceIndex=y,n&&(n.point?n.point.copy(p):n.point=p.clone(),n.point.applyMatrix4(Et),b.applyMatrix4(Et),n.distance=b.sub(n.point).length(),n.faceIndex=w),t)}function jn(){return typeof SharedArrayBuffer<"u"}const gt=new N.constructor,Vt=new N.constructor,we=new Fi(()=>new he),Ge=new he,Ye=new he,ni=new he,ai=new he;let oi=!1;function Gn(s,e,i,t){if(oi)throw new Error("MeshBVH: Recursive calls to bvhcast not supported.");oi=!0;const n=s._roots,a=e._roots;let r,c=0,l=0;const h=new Y().copy(i).invert();for(let f=0,u=n.length;f<u;f++){gt.setBuffer(n[f]),l=0;const o=we.getPrimitive();H(0,gt.float32Array,o),o.applyMatrix4(h);for(let m=0,v=a.length;m<v&&(Vt.setBuffer(a[f]),r=ce(0,0,i,h,t,c,l,0,0,o),Vt.clearBuffer(),l+=a[m].length,!r);m++);if(we.releasePrimitive(o),gt.clearBuffer(),c+=n[f].length,r)break}return oi=!1,r}function ce(s,e,i,t,n,a=0,r=0,c=0,l=0,h=null,f=!1){let u,o;f?(u=Vt,o=gt):(u=gt,o=Vt);const m=u.float32Array,v=u.uint32Array,b=u.uint16Array,d=o.float32Array,p=o.uint32Array,g=o.uint16Array,y=s*2,w=e*2,T=G(y,b),_=G(w,g);let A=!1;if(_&&T)f?A=n(Z(e,p),Q(e*2,g),Z(s,v),Q(s*2,b),l,r+e,c,a+s):A=n(Z(s,v),Q(s*2,b),Z(e,p),Q(e*2,g),c,a+s,l,r+e);else if(_){const S=we.getPrimitive();H(e,d,S),S.applyMatrix4(i);const R=se(s),I=J(s,v);H(R,m,Ge),H(I,m,Ye);const P=S.intersectsBox(Ge),M=S.intersectsBox(Ye);A=P&&ce(e,R,t,i,n,r,a,l,c+1,S,!f)||M&&ce(e,I,t,i,n,r,a,l,c+1,S,!f),we.releasePrimitive(S)}else{const S=se(e),R=J(e,p);H(S,d,ni),H(R,d,ai);const I=h.intersectsBox(ni),P=h.intersectsBox(ai);if(I&&P)A=ce(s,S,i,t,n,a,r,c,l+1,h,f)||ce(s,R,i,t,n,a,r,c,l+1,h,f);else if(I)if(T)A=ce(s,S,i,t,n,a,r,c,l+1,h,f);else{const M=we.getPrimitive();M.copy(ni).applyMatrix4(i);const F=se(s),E=J(s,v);H(F,m,Ge),H(E,m,Ye);const B=M.intersectsBox(Ge),z=M.intersectsBox(Ye);A=B&&ce(S,F,t,i,n,r,a,l,c+1,M,!f)||z&&ce(S,E,t,i,n,r,a,l,c+1,M,!f),we.releasePrimitive(M)}else if(P)if(T)A=ce(s,R,i,t,n,a,r,c,l+1,h,f);else{const M=we.getPrimitive();M.copy(ai).applyMatrix4(i);const F=se(s),E=J(s,v);H(F,m,Ge),H(E,m,Ye);const B=M.intersectsBox(Ge),z=M.intersectsBox(Ye);A=B&&ce(R,F,t,i,n,r,a,l,c+1,M,!f)||z&&ce(R,E,t,i,n,r,a,l,c+1,M,!f),we.releasePrimitive(M)}}return A}const Bt=new $,dr=new he,Yn={strategy:qr,maxDepth:40,maxLeafTris:10,useSharedArrayBuffer:!1,setBoundingBox:!0,onProgress:null,indirect:!1,verbose:!0};class Di{static serialize(e,i={}){i={cloneBuffers:!0,...i};const t=e.geometry,n=e._roots,a=e._indirectBuffer,r=t.getIndex();let c;return i.cloneBuffers?c={roots:n.map(l=>l.slice()),index:r?r.array.slice():null,indirectBuffer:a?a.slice():null}:c={roots:n,index:r?r.array:null,indirectBuffer:a},c}static deserialize(e,i,t={}){t={setIndex:!0,indirect:!!e.indirectBuffer,...t};const{index:n,roots:a,indirectBuffer:r}=e,c=new Di(i,{...t,[Qt]:!0});if(c._roots=a,c._indirectBuffer=r||null,t.setIndex){const l=i.getIndex();if(l===null){const h=new re(e.index,1,!1);i.setIndex(h)}else l.array!==n&&(l.array.set(n),l.needsUpdate=!0)}return c}get indirect(){return!!this._indirectBuffer}constructor(e,i={}){if(e.isBufferGeometry){if(e.index&&e.index.isInterleavedBufferAttribute)throw new Error("MeshBVH: InterleavedBufferAttribute is not supported for the index attribute.")}else throw new Error("MeshBVH: Only BufferGeometries are supported.");if(i=Object.assign({...Yn,[Qt]:!1},i),i.useSharedArrayBuffer&&!jn())throw new Error("MeshBVH: SharedArrayBuffer is not available.");this.geometry=e,this._roots=null,this._indirectBuffer=null,i[Qt]||(ln(this,i),!e.boundingBox&&i.setBoundingBox&&(e.boundingBox=this.getBoundingBox(new he))),this.resolveTriangleIndex=i.indirect?t=>this._indirectBuffer[t]:t=>t}refit(e=null){return(this.indirect?kn:wn)(this,e)}traverse(e,i=0){const t=this._roots[i],n=new Uint32Array(t),a=new Uint16Array(t);r(0);function r(c,l=0){const h=c*2,f=a[h+15]===jt;if(f){const u=n[c+6],o=a[h+14];e(l,f,new Float32Array(t,c*4,6),u,o)}else{const u=c+Se/4,o=n[c+6],m=n[c+7];e(l,f,new Float32Array(t,c*4,6),m)||(r(u,l+1),r(o,l+1))}}}raycast(e,i=mi,t=0,n=1/0){const a=this._roots,r=this.geometry,c=[],l=i.isMaterial,h=Array.isArray(i),f=r.groups,u=l?i.side:i,o=this.indirect?zn:Mn;for(let m=0,v=a.length;m<v;m++){const b=h?i[f[m].materialIndex].side:u,d=c.length;if(o(this,m,b,e,c,t,n),h){const p=f[m].materialIndex;for(let g=d,y=c.length;g<y;g++)c[g].face.materialIndex=p}}return c}raycastFirst(e,i=mi,t=0,n=1/0){const a=this._roots,r=this.geometry,c=i.isMaterial,l=Array.isArray(i);let h=null;const f=r.groups,u=c?i.side:i,o=this.indirect?Nn:In;for(let m=0,v=a.length;m<v;m++){const b=l?i[f[m].materialIndex].side:u,d=o(this,m,b,e,t,n);d!=null&&(h==null||d.distance<h.distance)&&(h=d,l&&(d.face.materialIndex=f[m].materialIndex))}return h}intersectsGeometry(e,i){let t=!1;const n=this._roots,a=this.indirect?Ln:Pn;for(let r=0,c=n.length;r<c&&(t=a(this,r,e,i),!t);r++);return t}shapecast(e){const i=ne.getPrimitive(),t=this.indirect?Sn:Tn;let{boundsTraverseOrder:n,intersectsBounds:a,intersectsRange:r,intersectsTriangle:c}=e;if(r&&c){const u=r;r=(o,m,v,b,d)=>u(o,m,v,b,d)?!0:t(o,m,this,c,v,b,i)}else r||(c?r=(u,o,m,v)=>t(u,o,this,c,m,v,i):r=(u,o,m)=>m);let l=!1,h=0;const f=this._roots;for(let u=0,o=f.length;u<o;u++){const m=f[u];if(l=pn(this,u,a,r,n,h),l)break;h+=m.byteLength}return ne.releasePrimitive(i),l}bvhcast(e,i,t){let{intersectsRanges:n,intersectsTriangles:a}=t;const r=ne.getPrimitive(),c=this.geometry.index,l=this.geometry.attributes.position,h=this.indirect?v=>{const b=this.resolveTriangleIndex(v);W(r,b*3,c,l)}:v=>{W(r,v*3,c,l)},f=ne.getPrimitive(),u=e.geometry.index,o=e.geometry.attributes.position,m=e.indirect?v=>{const b=e.resolveTriangleIndex(v);W(f,b*3,u,o)}:v=>{W(f,v*3,u,o)};if(a){const v=(b,d,p,g,y,w,T,_)=>{for(let A=p,S=p+g;A<S;A++){m(A),f.a.applyMatrix4(i),f.b.applyMatrix4(i),f.c.applyMatrix4(i),f.needsUpdate=!0;for(let R=b,I=b+d;R<I;R++)if(h(R),r.needsUpdate=!0,a(r,f,R,A,y,w,T,_))return!0}return!1};if(n){const b=n;n=function(d,p,g,y,w,T,_,A){return b(d,p,g,y,w,T,_,A)?!0:v(d,p,g,y,w,T,_,A)}}else n=v}return Gn(this,e,i,n)}intersectsBox(e,i){return Bt.set(e.min,e.max,i),Bt.needsUpdate=!0,this.shapecast({intersectsBounds:t=>Bt.intersectsBox(t),intersectsTriangle:t=>Bt.intersectsTriangle(t)})}intersectsSphere(e){return this.shapecast({intersectsBounds:i=>e.intersectsBox(i),intersectsTriangle:i=>i.intersectsSphere(e)})}closestPointToGeometry(e,i,t={},n={},a=0,r=1/0){return(this.indirect?qn:Bn)(this,e,i,t,n,a,r)}closestPointToPoint(e,i={},t=0,n=1/0){return gn(this,e,i,t,n)}getBoundingBox(e){return e.makeEmpty(),this._roots.forEach(t=>{H(0,new Float32Array(t),dr),e.union(dr)}),e}}function $n(s){switch(s){case 1:return"R";case 2:return"RG";case 3:return"RGBA";case 4:return"RGBA"}throw new Error}function Xn(s){switch(s){case 1:return Wt;case 2:return Lr;case 3:return V;case 4:return V}}function mr(s){switch(s){case 1:return Es;case 2:return Hr;case 3:return gi;case 4:return gi}}class Qr extends ae{constructor(){super(),this.minFilter=U,this.magFilter=U,this.generateMipmaps=!1,this.overrideItemSize=null,this._forcedType=null}updateFrom(e){const i=this.overrideItemSize,t=e.itemSize,n=e.count;if(i!==null){if(t*n%i!==0)throw new Error("VertexAttributeTexture: overrideItemSize must divide evenly into buffer length.");e.itemSize=i,e.count=n*t/i}const a=e.itemSize,r=e.count,c=e.normalized,l=e.array.constructor,h=l.BYTES_PER_ELEMENT;let f=this._forcedType,u=a;if(f===null)switch(l){case Float32Array:f=q;break;case Uint8Array:case Uint16Array:case Uint32Array:f=mt;break;case Int8Array:case Int16Array:case Int32Array:f=Xt;break}let o,m,v,b,d=$n(a);switch(f){case q:v=1,m=Xn(a),c&&h===1?(b=l,d+="8",l===Uint8Array?o=pi:(o=Ji,d+="_SNORM")):(b=Float32Array,d+="32F",o=q);break;case Xt:d+=h*8+"I",v=c?Math.pow(2,l.BYTES_PER_ELEMENT*8-1):1,m=mr(a),h===1?(b=Int8Array,o=Ji):h===2?(b=Int16Array,o=Ds):(b=Int32Array,o=Xt);break;case mt:d+=h*8+"UI",v=c?Math.pow(2,l.BYTES_PER_ELEMENT*8-1):1,m=mr(a),h===1?(b=Uint8Array,o=pi):h===2?(b=Uint16Array,o=Fs):(b=Uint32Array,o=mt);break}u===3&&(m===V||m===gi)&&(u=4);const p=Math.ceil(Math.sqrt(r))||1,g=u*p*p,y=new b(g),w=e.normalized;e.normalized=!1;for(let T=0;T<r;T++){const _=u*T;y[_]=e.getX(T)/v,a>=2&&(y[_+1]=e.getY(T)/v),a>=3&&(y[_+2]=e.getZ(T)/v,u===4&&(y[_+3]=1)),a>=4&&(y[_+3]=e.getW(T)/v)}e.normalized=w,this.internalFormat=d,this.format=m,this.type=o,this.image.width=p,this.image.height=p,this.image.data=y,this.needsUpdate=!0,this.dispose(),e.itemSize=t,e.count=n}}class Jr extends Qr{constructor(){super(),this._forcedType=mt}}class es extends Qr{constructor(){super(),this._forcedType=q}}class Kn{constructor(){this.index=new Jr,this.position=new es,this.bvhBounds=new ae,this.bvhContents=new ae,this._cachedIndexAttr=null,this.index.overrideItemSize=3}updateFrom(e){const{geometry:i}=e;if(Qn(e,this.bvhBounds,this.bvhContents),this.position.updateFrom(i.attributes.position),e.indirect){const t=e._indirectBuffer;if(this._cachedIndexAttr===null||this._cachedIndexAttr.count!==t.length)if(i.index)this._cachedIndexAttr=i.index.clone();else{const n=Yr(Gr(i));this._cachedIndexAttr=new re(n,1,!1)}Zn(i,t,this._cachedIndexAttr),this.index.updateFrom(this._cachedIndexAttr)}else this.index.updateFrom(i.index)}dispose(){const{index:e,position:i,bvhBounds:t,bvhContents:n}=this;e&&e.dispose(),i&&i.dispose(),t&&t.dispose(),n&&n.dispose()}}function Zn(s,e,i){const t=i.array,n=s.index?s.index.array:null;for(let a=0,r=e.length;a<r;a++){const c=3*a,l=3*e[a];for(let h=0;h<3;h++)t[c+h]=n?n[l+h]:l+h}}function Qn(s,e,i){const t=s._roots;if(t.length!==1)throw new Error("MeshBVHUniformStruct: Multi-root BVHs not supported.");const n=t[0],a=new Uint16Array(n),r=new Uint32Array(n),c=new Float32Array(n),l=n.byteLength/Se,h=2*Math.ceil(Math.sqrt(l/2)),f=new Float32Array(4*h*h),u=Math.ceil(Math.sqrt(l)),o=new Uint32Array(2*u*u);for(let m=0;m<l;m++){const v=m*Se/4,b=v*2,d=v;for(let p=0;p<3;p++)f[8*m+0+p]=c[d+0+p],f[8*m+4+p]=c[d+3+p];if(G(b,a)){const p=Q(b,a),g=Z(v,r),y=4294901760|p;o[m*2+0]=y,o[m*2+1]=g}else{const p=4*J(v,r)/Se,g=Pi(v,r);o[m*2+0]=g,o[m*2+1]=p}}e.image.data=f,e.image.width=h,e.image.height=h,e.format=V,e.type=q,e.internalFormat="RGBA32F",e.minFilter=U,e.magFilter=U,e.generateMipmaps=!1,e.needsUpdate=!0,e.dispose(),i.image.data=o,i.image.width=u,i.image.height=u,i.format=Hr,i.type=mt,i.internalFormat="RG32UI",i.minFilter=U,i.magFilter=U,i.generateMipmaps=!1,i.needsUpdate=!0,i.dispose()}const Jn=`
+
+// 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;
+
+}
+`,ea=`
+
+#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 ptr = 0;
+	uint stack[ BVH_STACK_DEPTH ];
+	stack[ 0 ] = 0u;
+
+	float triangleDistance = INFINITY;
+	bool found = false;
+	while ( ptr > - 1 && ptr < BVH_STACK_DEPTH ) {
+
+		uint currNodeIndex = stack[ ptr ];
+		ptr --;
+
+		// 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 = 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
+			ptr ++;
+			stack[ ptr ] = c2;
+
+			ptr ++;
+			stack[ ptr ] = c1;
+
+		}
+
+	}
+
+	return found;
+
+}
+`,ta=`
+struct BVH {
+
+	usampler2D index;
+	sampler2D position;
+
+	sampler2D bvhBounds;
+	usampler2D bvhContents;
+
+};
+`;function ts(s,e,i=0){if(s.isInterleavedBufferAttribute){const t=s.itemSize;for(let n=0,a=s.count;n<a;n++){const r=n+i;e.setX(r,s.getX(n)),t>=2&&e.setY(r,s.getY(n)),t>=3&&e.setZ(r,s.getZ(n)),t>=4&&e.setW(r,s.getW(n))}}else{const t=e.array,n=t.constructor,a=t.BYTES_PER_ELEMENT*s.itemSize*i;new n(t.buffer,a,s.array.length).set(s.array)}}function dt(s,e=null){const i=s.array.constructor,t=s.normalized,n=s.itemSize,a=e===null?s.count:e;return new re(new i(n*a),n,t)}function Ze(s,e){if(!s&&!e)return!0;if(!!s!=!!e)return!1;const i=s.count===e.count,t=s.normalized===e.normalized,n=s.array.constructor===e.array.constructor,a=s.itemSize===e.itemSize;return!(!i||!t||!n||!a)}function ia(s){const e=s[0].index!==null,i=new Set(Object.keys(s[0].attributes));if(!s[0].getAttribute("position"))throw new Error("StaticGeometryGenerator: position attribute is required.");for(let t=0;t<s.length;++t){const n=s[t];let a=0;if(e!==(n.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 r in n.attributes){if(!i.has(r))throw new Error('StaticGeometryGenerator: All geometries must have compatible attributes; make sure "'+r+'" attribute exists among all geometries, or in none of them.');a++}if(a!==i.size)throw new Error("StaticGeometryGenerator: All geometries must have the same number of attributes.")}}function ra(s){let e=0;for(let i=0,t=s.length;i<t;i++)e+=s[i].getIndex().count;return e}function sa(s){let e=0;for(let i=0,t=s.length;i<t;i++)e+=s[i].getAttribute("position").count;return e}function na(s,e,i){s.index&&s.index.count!==e&&s.setIndex(null);const t=s.attributes;for(const n in t)t[n].count!==i&&s.deleteAttribute(n)}function aa(s,e={},i=new De){const{useGroups:t=!1,forceUpdate:n=!1,skipAssigningAttributes:a=[],overwriteIndex:r=!0}=e;ia(s);const c=s[0].index!==null,l=c?ra(s):-1,h=sa(s);if(na(i,l,h),t){let u=0;for(let o=0,m=s.length;o<m;o++){const v=s[o];let b;c?b=v.getIndex().count:b=v.getAttribute("position").count,i.addGroup(u,b,o),u+=b}}if(c){let u=!1;if(i.index||(i.setIndex(new re(new Uint32Array(l),1,!1)),u=!0),u||r){let o=0,m=0;const v=i.getIndex();for(let b=0,d=s.length;b<d;b++){const p=s[b],g=p.getIndex();if(!(!n&&!u&&a[b]))for(let w=0;w<g.count;++w)v.setX(o+w,g.getX(w)+m);o+=g.count,m+=p.getAttribute("position").count}}}const f=Object.keys(s[0].attributes);for(let u=0,o=f.length;u<o;u++){let m=!1;const v=f[u];if(!i.getAttribute(v)){const p=s[0].getAttribute(v);i.setAttribute(v,dt(p,h)),m=!0}let b=0;const d=i.getAttribute(v);for(let p=0,g=s.length;p<g;p++){const y=s[p],w=!n&&!m&&a[p],T=y.getAttribute(v);w||ts(T,d,b),b+=T.count}}}function oa(s,e,i){const t=s.index,a=s.attributes.position.count,r=t?t.count:a;let c=s.groups;c.length===0&&(c=[{count:r,start:0,materialIndex:0}]);let l=s.getAttribute("materialIndex");if(!l||l.count!==a){let f;i.length<=255?f=new Uint8Array(a):f=new Uint16Array(a),l=new re(f,1,!1),s.deleteAttribute("materialIndex"),s.setAttribute("materialIndex",l)}const h=l.array;for(let f=0;f<c.length;f++){const u=c[f],o=u.start,m=u.count,v=Math.min(m,r-o),b=Array.isArray(e)?e[u.materialIndex]:e,d=i.indexOf(b);for(let p=0;p<v;p++){let g=o+p;t&&(g=t.getX(g)),h[g]=d}}}function ca(s,e){if(!s.index){const i=s.attributes.position.count,t=new Array(i);for(let n=0;n<i;n++)t[n]=n;s.setIndex(t)}if(!s.attributes.normal&&e&&e.includes("normal")&&s.computeVertexNormals(),!s.attributes.uv&&e&&e.includes("uv")){const i=s.attributes.position.count;s.setAttribute("uv",new re(new Float32Array(i*2),2,!1))}if(!s.attributes.uv2&&e&&e.includes("uv2")){const i=s.attributes.position.count;s.setAttribute("uv2",new re(new Float32Array(i*2),2,!1))}if(!s.attributes.tangent&&e&&e.includes("tangent"))if(s.attributes.uv&&s.attributes.normal)s.computeTangents();else{const i=s.attributes.position.count;s.setAttribute("tangent",new re(new Float32Array(i*4),4,!1))}if(!s.attributes.color&&e&&e.includes("color")){const i=s.attributes.position.count,t=new Float32Array(i*4);t.fill(1),s.setAttribute("color",new re(t,4))}}function Ei(s){let e=0;if(s.byteLength!==0){const i=new Uint8Array(s);for(let t=0;t<s.byteLength;t++){const n=i[t];e=(e<<5)-e+n,e|=0}}return e}function pr(s){let e=s.uuid;const i=Object.values(s.attributes);s.index&&(i.push(s.index),e+=`index|${s.index.version}`);const t=Object.keys(i).sort();for(const n of t){const a=i[n];e+=`${n}_${a.version}|`}return e}function gr(s){const e=s.skeleton;return e?(e.boneTexture||e.computeBoneTexture(),`${Ei(e.boneTexture.image.data.buffer)}_${e.boneTexture.uuid}`):null}class la{constructor(e=null){this.matrixWorld=new Y,this.geometryHash=null,this.skeletonHash=null,this.primitiveCount=-1,e!==null&&this.updateFrom(e)}updateFrom(e){const i=e.geometry,t=(i.index?i.index.count:i.attributes.position.count)/3;this.matrixWorld.copy(e.matrixWorld),this.geometryHash=pr(i),this.primitiveCount=t,this.skeletonHash=gr(e)}didChange(e){const i=e.geometry,t=(i.index?i.index.count:i.attributes.position.count)/3;return!(this.matrixWorld.equals(e.matrixWorld)&&this.geometryHash===pr(i)&&this.skeletonHash===gr(e)&&this.primitiveCount===t)}}const Ie=new C,Pe=new C,Ce=new C,vr=new Je,kt=new C,ci=new C,xr=new Je,yr=new Je,zt=new Y,br=new Y;function Tr(s,e,i){const t=s.skeleton,n=s.geometry,a=t.bones,r=t.boneInverses;xr.fromBufferAttribute(n.attributes.skinIndex,e),yr.fromBufferAttribute(n.attributes.skinWeight,e),zt.elements.fill(0);for(let c=0;c<4;c++){const l=yr.getComponent(c);if(l!==0){const h=xr.getComponent(c);br.multiplyMatrices(a[h].matrixWorld,r[h]),ua(zt,br,l)}}return zt.multiply(s.bindMatrix).premultiply(s.bindMatrixInverse),i.transformDirection(zt),i}function li(s,e,i,t,n){kt.set(0,0,0);for(let a=0,r=s.length;a<r;a++){const c=e[a],l=s[a];c!==0&&(ci.fromBufferAttribute(l,t),i?kt.addScaledVector(ci,c):kt.addScaledVector(ci.sub(n),c))}n.add(kt)}function ua(s,e,i){const t=s.elements,n=e.elements;for(let a=0,r=n.length;a<r;a++)t[a]+=n[a]*i}function fa(s){const{index:e,attributes:i}=s;if(e)for(let t=0,n=e.count;t<n;t+=3){const a=e.getX(t),r=e.getX(t+2);e.setX(t,r),e.setX(t+2,a)}else for(const t in i){const n=i[t],a=n.itemSize;for(let r=0,c=n.count;r<c;r+=3)for(let l=0;l<a;l++){const h=n.getComponent(r,l),f=n.getComponent(r+2,l);n.setComponent(r,l,f),n.setComponent(r+2,l,h)}}return s}function ha(s,e={},i=new De){e={applyWorldTransforms:!0,attributes:[],...e};const t=s.geometry,n=e.applyWorldTransforms,a=e.attributes.includes("normal"),r=e.attributes.includes("tangent"),c=t.attributes,l=i.attributes;for(const g in i.attributes)(!e.attributes.includes(g)||!(g in t.attributes))&&i.deleteAttribute(g);!i.index&&t.index&&(i.index=t.index.clone()),l.position||i.setAttribute("position",dt(c.position)),a&&!l.normal&&c.normal&&i.setAttribute("normal",dt(c.normal)),r&&!l.tangent&&c.tangent&&i.setAttribute("tangent",dt(c.tangent)),Ze(t.index,i.index),Ze(c.position,l.position),a&&Ze(c.normal,l.normal),r&&Ze(c.tangent,l.tangent);const h=c.position,f=a?c.normal:null,u=r?c.tangent:null,o=t.morphAttributes.position,m=t.morphAttributes.normal,v=t.morphAttributes.tangent,b=t.morphTargetsRelative,d=s.morphTargetInfluences,p=new Bs;p.getNormalMatrix(s.matrixWorld),t.index&&i.index.array.set(t.index.array);for(let g=0,y=c.position.count;g<y;g++)Ie.fromBufferAttribute(h,g),f&&Pe.fromBufferAttribute(f,g),u&&(vr.fromBufferAttribute(u,g),Ce.fromBufferAttribute(u,g)),d&&(o&&li(o,d,b,g,Ie),m&&li(m,d,b,g,Pe),v&&li(v,d,b,g,Ce)),s.isSkinnedMesh&&(s.applyBoneTransform(g,Ie),f&&Tr(s,g,Pe),u&&Tr(s,g,Ce)),n&&Ie.applyMatrix4(s.matrixWorld),l.position.setXYZ(g,Ie.x,Ie.y,Ie.z),f&&(n&&Pe.applyNormalMatrix(p),l.normal.setXYZ(g,Pe.x,Pe.y,Pe.z)),u&&(n&&Ce.transformDirection(s.matrixWorld),l.tangent.setXYZW(g,Ce.x,Ce.y,Ce.z,vr.w));for(const g in e.attributes){const y=e.attributes[g];y==="position"||y==="tangent"||y==="normal"||!(y in c)||(l[y]||i.setAttribute(y,dt(c[y])),Ze(c[y],l[y]),ts(c[y],l[y]))}return s.matrixWorld.determinant()<0&&fa(i),i}class da extends De{constructor(){super(),this.version=0,this.hash=null,this._diff=new la}isCompatible(e,i){const t=e.geometry;for(let n=0;n<i.length;n++){const a=i[n],r=t.attributes[a],c=this.attributes[a];if(r&&!Ze(r,c))return!1}return!0}updateFrom(e,i){const t=this._diff;return t.didChange(e)?(ha(e,i,this),t.updateFrom(e),this.version++,this.hash=`${this.uuid}_${this.version}`,!0):!1}}const Mi=0,is=1,rs=2;function ma(s,e){for(let i=0,t=s.length;i<t;i++)s[i].traverseVisible(a=>{a.isMesh&&e(a)})}function pa(s){const e=[];for(let i=0,t=s.length;i<t;i++){const n=s[i];Array.isArray(n.material)?e.push(...n.material):e.push(n.material)}return e}function ga(s,e,i){if(s.length===0){e.setIndex(null);const t=e.attributes;for(const n in t)e.deleteAttribute(n);for(const n in i.attributes)e.setAttribute(i.attributes[n],new re(new Float32Array(0),4,!1))}else aa(s,i,e);for(const t in e.attributes)e.attributes[t].needsUpdate=!0}class va{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 ks,i=new De;i.setAttribute("position",new re(new Float32Array(9),3)),this._dummyMesh=new zr(i,e)}return this._dummyMesh}_getMeshes(){const e=[];return ma(this.objects,i=>{e.push(i)}),e.sort((i,t)=>i.uuid>t.uuid?1:i.uuid<t.uuid?-1:0),e.length===0&&e.push(this._getDummyMesh()),e}_updateIntermediateGeometries(){const{_intermediateGeometry:e}=this,i=this._getMeshes(),t=new Set(e.keys()),n={attributes:this.attributes,applyWorldTransforms:this.applyWorldTransforms};for(let a=0,r=i.length;a<r;a++){const c=i[a],l=c.uuid;t.delete(l);let h=e.get(l);(!h||!h.isCompatible(c,this.attributes))&&(h&&h.dispose(),h=new da,e.set(l,h)),h.updateFrom(c,n)&&this.generateMissingAttributes&&ca(h,this.attributes)}t.forEach(a=>{e.delete(a)})}setObjects(e){Array.isArray(e)?this.objects=[...e]:this.objects=[e]}generate(e=new De){const{useGroups:i,overwriteIndex:t,_intermediateGeometry:n,_geometryMergeSets:a}=this,r=this._getMeshes(),c=[],l=[],h=a.get(e)||[];this._updateIntermediateGeometries();let f=!1;r.length!==h.length&&(f=!0);for(let o=0,m=r.length;o<m;o++){const v=r[o],b=n.get(v.uuid);l.push(b);const d=h[o];!d||d.uuid!==b.uuid?(c.push(!1),f=!0):d.version!==b.version?c.push(!1):c.push(!0)}ga(l,e,{useGroups:i,forceUpdate:f,skipAssigningAttributes:c,overwriteIndex:t}),f&&e.dispose(),a.set(e,l.map(o=>({version:o.version,uuid:o.uuid})));let u=Mi;return f?u=rs:c.includes(!1)&&(u=is),{changeType:u,materials:pa(r),geometry:e}}}function xa(s){const e=new Set;for(let i=0,t=s.length;i<t;i++){const n=s[i];for(const a in n){const r=n[a];r&&r.isTexture&&e.add(r)}}return Array.from(e)}function ya(s){const e=[],i=new Set;for(let n=0,a=s.length;n<a;n++)s[n].traverse(r=>{r.visible&&(r.isRectAreaLight||r.isSpotLight||r.isPointLight||r.isDirectionalLight)&&(e.push(r),r.iesMap&&i.add(r.iesMap))});const t=Array.from(i).sort((n,a)=>n.uuid<a.uuid?1:n.uuid>a.uuid?-1:0);return{lights:e,iesTextures:t}}class ba{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 De,this.staticGeometryGenerator=new va(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 i=this.generate(e);return this._buildAsync=!1,i.bvh=this.bvh=await i.bvh,i}}generate(e=null){const{staticGeometryGenerator:i,geometry:t,attributes:n}=this,a=i.objects;i.attributes=n,a.forEach(u=>{u.traverse(o=>{o.isSkinnedMesh&&o.skeleton&&o.skeleton.update()})});const r=i.generate(t),c=r.materials,l=xa(c),{lights:h,iesTextures:f}=ya(a);if(r.changeType!==Mi&&oa(t,c,c),this.generateBVH){if(this.bvh instanceof Promise)throw new Error("PathTracingSceneGenerator: BVH is already building asynchronously.");if(r.changeType===rs){const u={strategy:jr,maxLeafTris:1,indirect:!0,onProgress:e,...this.bvhOptions};this._buildAsync?this.bvh=this._bvhWorker.generate(t,u):this.bvh=new Di(t,u)}else r.changeType===is&&this.bvh.refit()}return{bvhChanged:r.changeType!==Mi,bvh:this.bvh,lights:h,iesTextures:f,geometry:t,materials:c,textures:l,objects:a}}}class Bi extends qt{set needsUpdate(e){super.needsUpdate=!0,this.dispatchEvent({type:"recompilation"})}constructor(e){super(e);for(const i in this.uniforms)Object.defineProperty(this,i,{get(){return this.uniforms[i].value},set(t){this.uniforms[i].value=t}})}setDefine(e,i=void 0){if(i==null){if(e in this.defines)return delete this.defines[e],this.needsUpdate=!0,!0}else if(this.defines[e]!==i)return this.defines[e]=i,this.needsUpdate=!0,!0;return!1}}class Ta extends Bi{constructor(e){super({blending:vt,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 Ot(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 Nt(s=1){let e="uint",i="float",t="",n=".r",a="1u";return s>1&&(e="uvec"+s,i="vec"+s,t=s+"",s===2?(n=".rg",a="uvec2( 1u, 2u )"):s===3?(n=".rgb",a="uvec3( 1u, 2u, 3u )"):(n="",a="uvec4( 1u, 2u, 3u, 4u )")),`
+
+		${i} sobol${t}( 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 );
+			${i} sobol_pt = sobolGetTexturePoint( shuffled_index )${n};
+			${e} result = ${e}( sobol_pt * 16777216.0 );
+
+			${e} seed2 = sobolHashCombine( seed, ${a} );
+			result = nestedUniformScrambleBase2( result, seed2 );
+
+			return SOBOL_FACTOR * ${i}( result >> 8 );
+
+		}
+	`}const ss=`
+
+	// Utils
+	const float SOBOL_FACTOR = 1.0 / 16777216.0;
+	const uint SOBOL_MAX_POINTS = 256u * 256u;
+
+	${Ot(1)}
+	${Ot(2)}
+	${Ot(3)}
+	${Ot(4)}
+
+	uint sobolHash( uint x ) {
+
+		// finalizer from murmurhash3
+		x ^= x >> 16;
+		x *= 0x85ebca6bu;
+		x ^= x >> 13;
+		x *= 0xc2b2ae35u;
+		x ^= x >> 16;
+		return x;
+
+	}
+
+`,wa=`
+
+	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;
+
+	}
+
+`,Aa=`
+
+	// 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 );
+
+	}
+
+	${Nt(1)}
+	${Nt(2)}
+	${Nt(3)}
+	${Nt(4)}
+
+`;class _a extends Bi{constructor(){super({blending:vt,uniforms:{resolution:{value:new L}},vertexShader:`
+
+				varying vec2 vUv;
+				void main() {
+
+					vUv = uv;
+					gl_Position = projectionMatrix * modelViewMatrix * vec4( position, 1.0 );
+
+				}
+			`,fragmentShader:`
+
+				${ss}
+				${wa}
+
+				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 Sa{generate(e,i=256){const t=new pt(i,i,{type:q,format:V,minFilter:U,magFilter:U,generateMipmaps:!1}),n=e.getRenderTarget();e.setRenderTarget(t);const a=new et(new _a);return a.material.resolution.set(i,i),a.render(e),e.setRenderTarget(n),a.dispose(),t}}class Ma extends Ur{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,i){return super.copy(e,i),this.fStop=e.fStop,this.apertureBlades=e.apertureBlades,this.apertureRotation=e.apertureRotation,this.focusDistance=e.focusDistance,this.anamorphicRatio=e.anamorphicRatio,this}}class Ra{constructor(){this.bokehSize=0,this.apertureBlades=0,this.apertureRotation=0,this.focusDistance=10,this.anamorphicRatio=1}updateFrom(e){e instanceof Ma?(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 ui(s){const e=new Uint16Array(s.length);for(let i=0,t=s.length;i<t;++i)e[i]=ve.toHalfFloat(s[i]);return e}function wr(s,e,i=0,t=s.length){let n=i,a=i+t-1;for(;n<a;){const r=n+a>>1;s[r]<e?n=r+1:a=r}return n-i}function Ia(s,e,i){return .2126*s+.7152*e+.0722*i}function Pa(s,e=fe){const i=s.clone();i.source=new zs({...i.image});const{width:t,height:n,data:a}=i.image;let r=a;if(i.type!==e){e===fe?r=new Uint16Array(a.length):r=new Float32Array(a.length);let c;a instanceof Int8Array||a instanceof Int16Array||a instanceof Int32Array?c=2**(8*a.BYTES_PER_ELEMENT-1)-1:c=2**(8*a.BYTES_PER_ELEMENT)-1;for(let l=0,h=a.length;l<h;l++){let f=a[l];i.type===fe&&(f=ve.fromHalfFloat(a[l])),i.type!==q&&i.type!==fe&&(f/=c),e===fe&&(r[l]=ve.toHalfFloat(f))}i.image.data=r,i.type=e}if(i.flipY){const c=r;r=r.slice();for(let l=0;l<n;l++)for(let h=0;h<t;h++){const f=n-l-1,u=4*(l*t+h),o=4*(f*t+h);r[o+0]=c[u+0],r[o+1]=c[u+1],r[o+2]=c[u+2],r[o+3]=c[u+3]}i.flipY=!1,i.image.data=r}return i}class Ca{constructor(){const e=new ae(ui(new Float32Array([0,0,0,0])),1,1);e.type=fe,e.format=V,e.minFilter=le,e.magFilter=le,e.wrapS=_e,e.wrapT=_e,e.generateMipmaps=!1,e.needsUpdate=!0;const i=new ae(ui(new Float32Array([0,1])),1,2);i.type=fe,i.format=Wt,i.minFilter=le,i.magFilter=le,i.generateMipmaps=!1,i.needsUpdate=!0;const t=new ae(ui(new Float32Array([0,0,1,1])),2,2);t.type=fe,t.format=Wt,t.minFilter=le,t.magFilter=le,t.generateMipmaps=!1,t.needsUpdate=!0,this.map=e,this.marginalWeights=i,this.conditionalWeights=t,this.totalSum=0}dispose(){this.marginalWeights.dispose(),this.conditionalWeights.dispose(),this.map.dispose()}updateFrom(e){const i=Pa(e);i.wrapS=_e,i.wrapT=Me;const{width:t,height:n,data:a}=i.image,r=new Float32Array(t*n),c=new Float32Array(t*n),l=new Float32Array(n),h=new Float32Array(n);let f=0,u=0;for(let d=0;d<n;d++){let p=0;for(let g=0;g<t;g++){const y=d*t+g,w=ve.fromHalfFloat(a[4*y+0]),T=ve.fromHalfFloat(a[4*y+1]),_=ve.fromHalfFloat(a[4*y+2]),A=Ia(w,T,_);p+=A,f+=A,r[y]=A,c[y]=p}if(p!==0)for(let g=d*t,y=d*t+t;g<y;g++)r[g]/=p,c[g]/=p;u+=p,l[d]=p,h[d]=u}if(u!==0)for(let d=0,p=l.length;d<p;d++)l[d]/=u,h[d]/=u;const o=new Uint16Array(n),m=new Uint16Array(t*n);for(let d=0;d<n;d++){const p=(d+1)/n,g=wr(h,p);o[d]=ve.toHalfFloat((g+.5)/n)}for(let d=0;d<n;d++)for(let p=0;p<t;p++){const g=d*t+p,y=(p+1)/t,w=wr(c,y,d*t,t);m[g]=ve.toHalfFloat((w+.5)/t)}this.dispose();const{marginalWeights:v,conditionalWeights:b}=this;v.image={width:n,height:1,data:o},v.needsUpdate=!0,b.image={width:t,height:n,data:m},b.needsUpdate=!0,this.totalSum=f,this.map=i}}const fi=6,Fa=0,Da=1,Ea=2,Ba=3,ka=4,oe=new C,K=new C,Ar=new Y,$e=new hi,_r=new C,Xe=new C,za=new C(0,1,0);class Oa{constructor(){const e=new ae(new Float32Array(4),1,1);e.format=V,e.type=q,e.wrapS=Me,e.wrapT=Me,e.generateMipmaps=!1,e.minFilter=U,e.magFilter=U,this.tex=e,this.count=0}updateFrom(e,i=[]){const t=this.tex,n=Math.max(e.length*fi,1),a=Math.ceil(Math.sqrt(n));t.image.width!==a&&(t.dispose(),t.image.data=new Float32Array(a*a*4),t.image.width=a,t.image.height=a);const r=t.image.data;for(let l=0,h=e.length;l<h;l++){const f=e[l],u=l*fi*4;let o=0;for(let v=0;v<fi*4;v++)r[u+v]=0;f.getWorldPosition(K),r[u+o++]=K.x,r[u+o++]=K.y,r[u+o++]=K.z;let m=Fa;if(f.isRectAreaLight&&f.isCircular?m=Da:f.isSpotLight?m=Ea:f.isDirectionalLight?m=Ba:f.isPointLight&&(m=ka),r[u+o++]=m,r[u+o++]=f.color.r,r[u+o++]=f.color.g,r[u+o++]=f.color.b,r[u+o++]=f.intensity,f.getWorldQuaternion($e),f.isRectAreaLight)oe.set(f.width,0,0).applyQuaternion($e),r[u+o++]=oe.x,r[u+o++]=oe.y,r[u+o++]=oe.z,o++,K.set(0,f.height,0).applyQuaternion($e),r[u+o++]=K.x,r[u+o++]=K.y,r[u+o++]=K.z,r[u+o++]=oe.cross(K).length()*(f.isCircular?Math.PI/4:1);else if(f.isSpotLight){const v=f.radius||0;_r.setFromMatrixPosition(f.matrixWorld),Xe.setFromMatrixPosition(f.target.matrixWorld),Ar.lookAt(_r,Xe,za),$e.setFromRotationMatrix(Ar),oe.set(1,0,0).applyQuaternion($e),r[u+o++]=oe.x,r[u+o++]=oe.y,r[u+o++]=oe.z,o++,K.set(0,1,0).applyQuaternion($e),r[u+o++]=K.x,r[u+o++]=K.y,r[u+o++]=K.z,r[u+o++]=Math.PI*v*v,r[u+o++]=v,r[u+o++]=f.decay,r[u+o++]=f.distance,r[u+o++]=Math.cos(f.angle),r[u+o++]=Math.cos(f.angle*(1-f.penumbra)),r[u+o++]=f.iesMap?i.indexOf(f.iesMap):-1}else if(f.isPointLight){const v=oe.setFromMatrixPosition(f.matrixWorld);r[u+o++]=v.x,r[u+o++]=v.y,r[u+o++]=v.z,o++,o+=4,o+=1,r[u+o++]=f.decay,r[u+o++]=f.distance}else if(f.isDirectionalLight){const v=oe.setFromMatrixPosition(f.matrixWorld),b=K.setFromMatrixPosition(f.target.matrixWorld);Xe.subVectors(v,b).normalize(),r[u+o++]=Xe.x,r[u+o++]=Xe.y,r[u+o++]=Xe.z}}this.count=e.length;const c=Ei(r.buffer);return this.hash!==c?(this.hash=c,t.needsUpdate=!0,!0):!1}}function Sr(s,e,i,t,n){if(e>t)throw new Error;const a=s.length/e,r=s.constructor.BYTES_PER_ELEMENT*8;let c=1;switch(s.constructor){case Uint8Array:case Uint16Array:case Uint32Array:c=2**r-1;break;case Int8Array:case Int16Array:case Int32Array:c=2**(r-1)-1;break}for(let l=0;l<a;l++){const h=4*l,f=e*l;for(let u=0;u<t;u++)i[n+h+u]=e>=u+1?s[f+u]/c:0}}class Na extends Os{constructor(){super(),this._textures=[],this.type=q,this.format=V,this.internalFormat="RGBA32F"}updateAttribute(e,i){const t=this._textures[e];t.updateFrom(i);const n=t.image,a=this.image;if(n.width!==a.width||n.height!==a.height)throw new Error("FloatAttributeTextureArray: Attribute must be the same dimensions when updating single layer.");const{width:r,height:c,data:l}=a,f=r*c*4*e;let u=i.itemSize;u===3&&(u=4),Sr(t.image.data,u,l,4,f),this.dispose(),this.needsUpdate=!0}setAttributes(e){const i=e[0].count,t=e.length;for(let u=0,o=t;u<o;u++)if(e[u].count!==i)throw new Error("FloatAttributeTextureArray: All attributes must have the same item count.");const n=this._textures;for(;n.length<t;){const u=new es;n.push(u)}for(;n.length>t;)n.pop();for(let u=0,o=t;u<o;u++)n[u].updateFrom(e[u]);const r=n[0].image,c=this.image;(r.width!==c.width||r.height!==c.height||r.depth!==t)&&(c.width=r.width,c.height=r.height,c.depth=t,c.data=new Float32Array(c.width*c.height*c.depth*4));const{data:l,width:h,height:f}=c;for(let u=0,o=t;u<o;u++){const m=n[u],b=h*f*4*u;let d=e[u].itemSize;d===3&&(d=4),Sr(m.image.data,d,l,4,b)}this.dispose(),this.needsUpdate=!0}}class La extends Na{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,i,t,n){this.setAttributes([e,i,t,n])}}function ki(s,e){return s.uuid<e.uuid?1:s.uuid>e.uuid?-1:0}function Ri(s){return`${s.source.uuid}:${s.colorSpace}`}function Ha(s){const e=new Set,i=[];for(let t=0,n=s.length;t<n;t++){const a=s[t],r=Ri(a);e.has(r)||(e.add(r),i.push(a))}return i}function Ua(s){const e=s.map(t=>t.iesMap||null).filter(t=>t),i=new Set(e);return Array.from(i).sort(ki)}function Wa(s){const e=new Set;for(let t=0,n=s.length;t<n;t++){const a=s[t];for(const r in a){const c=a[r];c&&c.isTexture&&e.add(c)}}const i=Array.from(e);return Ha(i).sort(ki)}function Va(s){const e=[];return s.traverse(i=>{i.visible&&(i.isRectAreaLight||i.isSpotLight||i.isPointLight||i.isDirectionalLight)&&e.push(i)}),e.sort(ki)}const ns=45,Mr=ns*4;class qa{constructor(){this._features={}}isUsed(e){return e in this._features}setUsed(e,i=!0){i===!1?delete this._features[e]:this._features[e]=!0}reset(){this._features={}}}class ja extends ae{constructor(){super(new Float32Array(4),1,1),this.format=V,this.type=q,this.wrapS=Me,this.wrapT=Me,this.minFilter=U,this.magFilter=U,this.generateMipmaps=!1,this.features=new qa}updateFrom(e,i){function t(v,b,d=-1){if(b in v&&v[b]){const p=Ri(v[b]);return u[p]}else return d}function n(v,b,d){return b in v?v[b]:d}function a(v,b,d,p){const g=v[b]&&v[b].isTexture?v[b]:null;if(g){g.matrixAutoUpdate&&g.updateMatrix();const y=g.matrix.elements;let w=0;d[p+w++]=y[0],d[p+w++]=y[3],d[p+w++]=y[6],w++,d[p+w++]=y[1],d[p+w++]=y[4],d[p+w++]=y[7],w++}return 8}let r=0;const c=e.length*ns,l=Math.ceil(Math.sqrt(c))||1,{image:h,features:f}=this,u={};for(let v=0,b=i.length;v<b;v++)u[Ri(i[v])]=v;h.width!==l&&(this.dispose(),h.data=new Float32Array(l*l*4),h.width=l,h.height=l);const o=h.data;f.reset();for(let v=0,b=e.length;v<b;v++){const d=e[v];if(d.isFogVolumeMaterial){f.setUsed("FOG");for(let y=0;y<Mr;y++)o[r+y]=0;o[r+0*4+0]=d.color.r,o[r+0*4+1]=d.color.g,o[r+0*4+2]=d.color.b,o[r+2*4+3]=n(d,"emissiveIntensity",0),o[r+3*4+0]=d.emissive.r,o[r+3*4+1]=d.emissive.g,o[r+3*4+2]=d.emissive.b,o[r+13*4+1]=d.density,o[r+13*4+3]=0,o[r+14*4+2]=4,r+=Mr;continue}o[r++]=d.color.r,o[r++]=d.color.g,o[r++]=d.color.b,o[r++]=t(d,"map"),o[r++]=n(d,"metalness",0),o[r++]=t(d,"metalnessMap"),o[r++]=n(d,"roughness",0),o[r++]=t(d,"roughnessMap"),o[r++]=n(d,"ior",1.5),o[r++]=n(d,"transmission",0),o[r++]=t(d,"transmissionMap"),o[r++]=n(d,"emissiveIntensity",0),"emissive"in d?(o[r++]=d.emissive.r,o[r++]=d.emissive.g,o[r++]=d.emissive.b):(o[r++]=0,o[r++]=0,o[r++]=0),o[r++]=t(d,"emissiveMap"),o[r++]=t(d,"normalMap"),"normalScale"in d?(o[r++]=d.normalScale.x,o[r++]=d.normalScale.y):(o[r++]=1,o[r++]=1),o[r++]=n(d,"clearcoat",0),o[r++]=t(d,"clearcoatMap"),o[r++]=n(d,"clearcoatRoughness",0),o[r++]=t(d,"clearcoatRoughnessMap"),o[r++]=t(d,"clearcoatNormalMap"),"clearcoatNormalScale"in d?(o[r++]=d.clearcoatNormalScale.x,o[r++]=d.clearcoatNormalScale.y):(o[r++]=1,o[r++]=1),r++,o[r++]=n(d,"sheen",0),"sheenColor"in d?(o[r++]=d.sheenColor.r,o[r++]=d.sheenColor.g,o[r++]=d.sheenColor.b):(o[r++]=0,o[r++]=0,o[r++]=0),o[r++]=t(d,"sheenColorMap"),o[r++]=n(d,"sheenRoughness",0),o[r++]=t(d,"sheenRoughnessMap"),o[r++]=t(d,"iridescenceMap"),o[r++]=t(d,"iridescenceThicknessMap"),o[r++]=n(d,"iridescence",0),o[r++]=n(d,"iridescenceIOR",1.3);const p=n(d,"iridescenceThicknessRange",[100,400]);o[r++]=p[0],o[r++]=p[1],"specularColor"in d?(o[r++]=d.specularColor.r,o[r++]=d.specularColor.g,o[r++]=d.specularColor.b):(o[r++]=1,o[r++]=1,o[r++]=1),o[r++]=t(d,"specularColorMap"),o[r++]=n(d,"specularIntensity",1),o[r++]=t(d,"specularIntensityMap");const g=n(d,"thickness",0)===0&&n(d,"attenuationDistance",1/0)===1/0;if(o[r++]=Number(g),r++,"attenuationColor"in d?(o[r++]=d.attenuationColor.r,o[r++]=d.attenuationColor.g,o[r++]=d.attenuationColor.b):(o[r++]=1,o[r++]=1,o[r++]=1),o[r++]=n(d,"attenuationDistance",1/0),o[r++]=t(d,"alphaMap"),o[r++]=d.opacity,o[r++]=d.alphaTest,!g&&d.transmission>0)o[r++]=0;else switch(d.side){case mi:o[r++]=1;break;case Or:o[r++]=-1;break;case Nr:o[r++]=0;break}o[r++]=Number(n(d,"matte",!1)),o[r++]=Number(n(d,"castShadow",!0)),o[r++]=Number(d.vertexColors)|Number(d.flatShading)<<1,o[r++]=Number(d.transparent),r+=a(d,"map",o,r),r+=a(d,"metalnessMap",o,r),r+=a(d,"roughnessMap",o,r),r+=a(d,"transmissionMap",o,r),r+=a(d,"emissiveMap",o,r),r+=a(d,"normalMap",o,r),r+=a(d,"clearcoatMap",o,r),r+=a(d,"clearcoatNormalMap",o,r),r+=a(d,"clearcoatRoughnessMap",o,r),r+=a(d,"sheenColorMap",o,r),r+=a(d,"sheenRoughnessMap",o,r),r+=a(d,"iridescenceMap",o,r),r+=a(d,"iridescenceThicknessMap",o,r),r+=a(d,"specularColorMap",o,r),r+=a(d,"specularIntensityMap",o,r)}const m=Ei(o.buffer);return this.hash!==m?(this.hash=m,this.needsUpdate=!0,!0):!1}}const Rr=new xt;function Ga(s){return s?`${s.uuid}:${s.version}`:null}function Ya(s,e){for(const i in e)i in s&&(s[i]=e[i])}class Ir extends Ns{constructor(e,i,t){const n={format:V,type:pi,minFilter:le,magFilter:le,wrapS:_e,wrapT:_e,generateMipmaps:!1,...t};super(e,i,1,n),Ya(this.texture,n),this.texture.setTextures=(...r)=>{this.setTextures(...r)},this.hashes=[null];const a=new et(new $a);this.fsQuad=a}setTextures(e,i,t=this.width,n=this.height){const a=e.getRenderTarget(),r=e.toneMapping,c=e.getClearAlpha();e.getClearColor(Rr);const l=i.length||1;(t!==this.width||n!==this.height||this.depth!==l)&&(this.setSize(t,n,l),this.hashes=new Array(l).fill(null)),e.setClearColor(0,0),e.toneMapping=Ls;const h=this.fsQuad,f=this.hashes;let u=!1;for(let o=0,m=l;o<m;o++){const v=i[o],b=Ga(v);v&&(f[o]!==b||v.isWebGLRenderTarget)&&(v.matrixAutoUpdate=!1,v.matrix.identity(),h.material.map=v,e.setRenderTarget(this,o),h.render(e),v.updateMatrix(),v.matrixAutoUpdate=!0,f[o]=b,u=!0)}return h.material.map=null,e.setClearColor(Rr,c),e.setRenderTarget(a),e.toneMapping=r,u}dispose(){super.dispose(),this.fsQuad.dispose()}}class $a extends qt{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 Xa(s,e=Math.random()){for(let i=s.length-1;i>0;i--){const t=Math.floor(e()*(i+1)),n=s[i];s[i]=s[t],s[t]=n}return s}class Ka{constructor(e,i,t=Math.random){const n=e**i,a=new Uint16Array(n);let r=n;for(let c=0;c<n;c++)a[c]=c;this.samples=new Float32Array(i),this.strataCount=e,this.reset=function(){for(let c=0;c<n;c++)a[c]=c;r=0},this.reshuffle=function(){r=0},this.next=function(){const{samples:c}=this;r>=a.length&&(Xa(a,t),this.reshuffle());let l=a[r++];for(let h=0;h<i;h++)c[h]=(l%e+t())/e,l=Math.floor(l/e);return c}}}class Za{constructor(e,i,t=Math.random){let n=0;for(const l of i)n+=l;const a=new Float32Array(n),r=[];let c=0;for(const l of i){const h=new Ka(e,l,t);h.samples=new Float32Array(a.buffer,c,h.samples.length),c+=h.samples.length*4,r.push(h)}this.samples=a,this.strataCount=e,this.next=function(){for(const l of r)l.next();return a},this.reshuffle=function(){for(const l of r)l.reshuffle()},this.reset=function(){for(const l of r)l.reset()}}}class Qa{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 Ja extends ae{constructor(e=1,i=1,t=8){super(new Float32Array(1),1,1,V,q),this.minFilter=U,this.magFilter=U,this.strata=t,this.sampler=null,this.generator=new Qa,this.stableNoise=!1,this.random=()=>this.stableNoise?this.generator.nextFloat():Math.random(),this.init(e,i,t)}init(e=this.image.height,i=this.image.width,t=this.strata){const{image:n}=this;if(n.width===i&&n.height===e&&this.sampler!==null)return;const a=new Array(e*i).fill(4),r=new Za(t,a,this.random);n.width=i,n.height=e,n.data=r.samples,this.sampler=r,this.dispose(),this.next()}next(){this.sampler.next(),this.needsUpdate=!0}reset(){this.sampler.reset(),this.generator.seed=0}}function eo(s,e=Math.random){for(let i=s.length-1;i>0;i--){const t=~~((e()-1e-6)*i),n=s[i];s[i]=s[t],s[t]=n}}function to(s,e){s.fill(0);for(let i=0;i<e;i++)s[i]=1}class Pr{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:i}=this;let t=1/0,n=-1;for(let a=0,r=i.length;a<r;a++){if(i[a]!==0)continue;const c=e[a];c<t&&(t=c,n=a)}return n}findCluster(){const{score:e,binaryPattern:i}=this;let t=-1/0,n=-1;for(let a=0,r=i.length;a<r;a++){if(i[a]!==1)continue;const c=e[a];c>t&&(t=c,n=a)}return n}setSigma(e){if(e===this.sigma)return;const i=~~(Math.sqrt(10*2*e**2)+1),t=2*i+1,n=new Float32Array(t*t),a=e*e;for(let r=-i;r<=i;r++)for(let c=-i;c<=i;c++){const l=(i+c)*t+r+i,h=r*r+c*c;n[l]=Math.E**(-h/(2*a))}this.lookupTable=n,this.sigma=e,this.radius=i}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:i,size:t}=this;i.fill(0);for(let n=0,a=e.length;n<a;n++)if(e[n]===0){const r=~~(n/t),c=n-r*t;this.updateScore(c,r,1),e[n]=1}else e[n]=0}updateScore(e,i,t){const{size:n,score:a,lookupTable:r}=this,c=this.radius,l=2*c+1;for(let h=-c;h<=c;h++)for(let f=-c;f<=c;f++){const u=(c+f)*l+h+c,o=r[u];let m=e+h;m=m<0?n+m:m%n;let v=i+f;v=v<0?n+v:v%n;const b=v*n+m;a[b]+=t*o}}addPointIndex(e){this.binaryPattern[e]=1;const i=this.size,t=~~(e/i),n=e-t*i;this.updateScore(n,t,1),this.count++}removePointIndex(e){this.binaryPattern[e]=0;const i=this.size,t=~~(e/i),n=e-t*i;this.updateScore(n,t,-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 io{constructor(){this.random=Math.random,this.sigma=1.5,this.size=64,this.majorityPointsRatio=.1,this.samples=new Pr(1),this.savedSamples=new Pr(1)}generate(){const{samples:e,savedSamples:i,sigma:t,majorityPointsRatio:n,size:a}=this;e.resize(a),e.setSigma(t);const r=Math.floor(a*a*n),c=e.binaryPattern;to(c,r),eo(c,this.random);for(let u=0,o=c.length;u<o;u++)c[u]===1&&e.addPointIndex(u);for(;;){const u=e.findCluster();e.removePointIndex(u);const o=e.findVoid();if(u===o){e.addPointIndex(u);break}e.addPointIndex(o)}const l=new Uint32Array(a*a);i.copy(e);let h;for(h=e.count-1;h>=0;){const u=e.findCluster();e.removePointIndex(u),l[u]=h,h--}const f=a*a;for(h=i.count;h<f/2;){const u=i.findVoid();i.addPointIndex(u),l[u]=h,h++}for(i.invert();h<f;){const u=i.findCluster();i.removePointIndex(u),l[u]=h,h++}return{data:l,maxValue:f}}}function ro(s){return s>=3?4:s}function so(s){switch(s){case 1:return Wt;case 2:return Lr;default:return V}}class no extends ae{constructor(e=64,i=1){super(new Float32Array(4),1,1,V,q),this.minFilter=U,this.magFilter=U,this.size=e,this.channels=i,this.update()}update(){const e=this.channels,i=this.size,t=new io;t.channels=e,t.size=i;const n=ro(e),a=so(n);(this.image.width!==i||a!==this.format)&&(this.image.width=i,this.image.height=i,this.image.data=new Float32Array(i**2*n),this.format=a,this.dispose());const r=this.image.data;for(let c=0,l=e;c<l;c++){const h=t.generate(),f=h.data,u=h.maxValue;for(let o=0,m=f.length;o<m;o++){const v=f[o]/u;r[o*n+c]=v}}this.needsUpdate=!0}}const ao=`
+
+	struct PhysicalCamera {
+
+		float focusDistance;
+		float anamorphicRatio;
+		float bokehSize;
+		int apertureBlades;
+		float apertureRotation;
+
+	};
+
+`,oo=`
+
+	struct EquirectHdrInfo {
+
+		sampler2D marginalWeights;
+		sampler2D conditionalWeights;
+		sampler2D map;
+
+		float totalSum;
+
+	};
+
+`,co=`
+
+	#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;
+
+	}
+
+`,lo=`
+
+	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;
+
+	}
+
+`,uo=`
+
+	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;
+	};
+
+`,fo=`
+
+	// 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 );
+
+	}
+`,ho=`
+
+	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;
+
+	}
+
+`,mo=`
+
+	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 );
+
+	}
+
+
+`,po=`
+
+	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 );
+
+	}
+
+`,go=`
+
+	// 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 );
+
+	}
+
+`,vo=`
+
+	// 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;
+
+	}
+
+`,xo=`
+
+	// 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 );
+
+	}
+
+`,as=`
+
+	// 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 ) );
+
+	}
+`,Cr=`
+
+	// 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 );
+
+	}
+`,yo=`
+
+	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 );
+
+	}
+
+`,bo=`
+
+	// 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;
+
+	}
+
+`,To=`
+
+	// 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;
+
+	}
+
+`,wo=`
+
+	// 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;
+
+	}
+
+`,Ao=`
+
+	// 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 ) );
+
+	}
+
+`,_o=`
+
+	// 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;
+
+	}
+
+`,So=`
+
+#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;
+
+}
+
+`,Mo=`
+
+	// 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;
+
+	}
+
+`,Ro=`
+
+	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;
+
+	}
+
+`,Io=`
+
+	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;
+
+	}
+
+`,Po=`
+
+	#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;
+
+	}
+`,Co=`
+
+	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;
+
+	}
+
+`,Fo=`
+
+	#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 Do extends Bi{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},opacity:{value:1},bounces:{value:10},transmissiveBounces:{value:10},filterGlossyFactor:{value:0},physicalCamera:{value:new Ra},cameraWorldMatrix:{value:new Y},invProjectionMatrix:{value:new Y},bvh:{value:new Kn},attributesArray:{value:new La},materialIndexAttribute:{value:new Jr},materials:{value:new ja},textures:{value:new Ir().texture},lights:{value:new Oa},iesProfiles:{value:new Ir(360,180,{type:fe,wrapS:Me,wrapT:Me}).texture},environmentIntensity:{value:1},environmentRotation:{value:new Y},envMapInfo:{value:new Ca},backgroundBlur:{value:0},backgroundMap:{value:null},backgroundAlpha:{value:1},backgroundIntensity:{value:1},backgroundRotation:{value:new Y},seed:{value:0},sobolTexture:{value:null},stratifiedTexture:{value:new Ja},stratifiedOffsetTexture:{value:new no(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
+				${Jn}
+				${ta}
+				${ea}
+
+				// uniform structs
+				${ao}
+				${co}
+				${oo}
+				${lo}
+				${uo}
+
+				// random
+				#if RANDOM_TYPE == 2 	// Stratified List
+
+					${yo}
+
+				#elif RANDOM_TYPE == 1 	// Sobol
+
+					${Cr}
+					${ss}
+					${Aa}
+
+					#define rand(v) sobol(v)
+					#define rand2(v) sobol2(v)
+					#define rand3(v) sobol3(v)
+					#define rand4(v) sobol4(v)
+
+				#else 					// PCG
+
+				${Cr}
+
+					// 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
+				${xo}
+				${po}
+				${as}
+				${go}
+				${vo}
+
+				// 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
+				${mo}
+				${fo}
+				${ho}
+
+				${So}
+				${wo}
+				${_o}
+				${Ao}
+				${To}
+				${bo}
+
+				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
+
+				}
+
+				${Co}
+				${Ro}
+				${Fo}
+				${Mo}
+				${Io}
+				${Po}
+
+				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*Eo(){const{_renderer:s,_fsQuad:e,_blendQuad:i,_primaryTarget:t,_blendTargets:n,_sobolTarget:a,_subframe:r,alpha:c,material:l}=this,h=new Je,f=new Je,u=i.material;let[o,m]=n;for(;;){c?(u.opacity=this._opacityFactor/(this.samples+1),l.blending=vt,l.opacity=1):(l.opacity=this._opacityFactor/(this.samples+1),l.blending=Wr);const[v,b,d,p]=r,g=t.width,y=t.height;l.resolution.set(g*d,y*p),l.sobolTexture=a.texture,l.stratifiedTexture.init(20,l.bounces+l.transmissiveBounces+5),l.stratifiedTexture.next(),l.seed++;const w=this.tiles.x||1,T=this.tiles.y||1,_=w*T,A=Math.ceil(g*d),S=Math.ceil(y*p),R=Math.floor(v*g),I=Math.floor(b*y),P=Math.ceil(A/w),M=Math.ceil(S/T);for(let F=0;F<T;F++)for(let E=0;E<w;E++){const B=s.getRenderTarget(),z=s.autoClear,X=s.getScissorTest();s.getScissor(h),s.getViewport(f);let ee=E,te=F;if(!this.stableTiles){const Be=this._currentTile%(w*T);ee=Be%w,te=~~(Be/w),this._currentTile=Be+1}const Ee=T-te-1;t.scissor.set(R+ee*P,I+Ee*M,Math.min(P,A-ee*P),Math.min(M,S-Ee*M)),t.viewport.set(R,I,A,S),s.setRenderTarget(t),s.setScissorTest(!0),s.autoClear=!1,e.render(s),s.setViewport(f),s.setScissor(h),s.setScissorTest(X),s.setRenderTarget(B),s.autoClear=z,c&&(u.target1=o.texture,u.target2=t.texture,s.setRenderTarget(m),i.render(s),s.setRenderTarget(B)),this.samples+=1/_,E===w-1&&F===T-1&&(this.samples=Math.round(this.samples)),yield}[o,m]=[m,o]}}const Fr=new xt;class Dr{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(3,3),this.stableNoise=!1,this.stableTiles=!0,this.samples=0,this._subframe=new Je(0,0,1,1),this._opacityFactor=1,this._renderer=e,this._alpha=!1,this._fsQuad=new et(new Do),this._blendQuad=new et(new Ta),this._task=null,this._currentTile=0,this._compilePromise=null,this._sobolTarget=new Sa().generate(e),this._primaryTarget=new pt(1,1,{format:V,type:q,magFilter:U,minFilter:U}),this._blendTargets=[new pt(1,1,{format:V,type:q,magFilter:U,minFilter:U}),new pt(1,1,{format:V,type:q,magFilter:U,minFilter:U})],this._compileFunction=()=>{const i=this.compileMaterial(this._fsQuad._mesh);i.then(()=>{this._compilePromise===i&&(this._compilePromise=null)}),this._compilePromise=i},this.material.addEventListener("recompilation",this._compileFunction)}compileMaterial(){return this._renderer.compileAsync(this._fsQuad._mesh)}setCamera(e){const{material:i}=this;i.cameraWorldMatrix.copy(e.matrixWorld),i.invProjectionMatrix.copy(e.projectionMatrixInverse),i.physicalCamera.updateFrom(e);let t=0;e.projectionMatrix.elements[15]>0&&(t=1),e.isEquirectCamera&&(t=2),i.setDefine("CAMERA_TYPE",t),this.camera=e}setSize(e,i){e=Math.ceil(e),i=Math.ceil(i),!(this._primaryTarget.width===e&&this._primaryTarget.height===i)&&(this._primaryTarget.setSize(e,i),this._blendTargets[0].setSize(e,i),this._blendTargets[1].setSize(e,i),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:i,_blendTargets:t}=this,n=e.getRenderTarget(),a=e.getClearAlpha();e.getClearColor(Fr),e.setRenderTarget(i),e.setClearColor(0,0),e.clearColor(),e.setRenderTarget(t[0]),e.setClearColor(0,0),e.clearColor(),e.setRenderTarget(t[1]),e.setClearColor(0,0),e.clearColor(),e.setClearColor(Fr,a),e.setRenderTarget(n),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=Eo.call(this)),this._task.next())}}const Fe=new L,Er=new L,Lt=new di,Ht=new xt;class Bo extends ae{constructor(e=512,i=512){super(new Float32Array(e*i*4),e,i,V,q,Vr,_e,Me,le,le),this.generationCallback=null}update(){this.dispose(),this.needsUpdate=!0;const{data:e,width:i,height:t}=this.image;for(let n=0;n<i;n++)for(let a=0;a<t;a++){Er.set(i,t),Fe.set(n/i,a/t),Fe.x-=.5,Fe.y=1-Fe.y,Lt.theta=Fe.x*2*Math.PI,Lt.phi=Fe.y*Math.PI,Lt.radius=1,this.generationCallback(Lt,Fe,Er,Ht);const c=4*(a*i+n);e[c+0]=Ht.r,e[c+1]=Ht.g,e[c+2]=Ht.b,e[c+3]=1}}copy(e){return super.copy(e),this.generationCallback=e.generationCallback,this}}const Br=new C;class ko extends Bo{constructor(e=512){super(e,e),this.topColor=new xt().set(16777215),this.bottomColor=new xt().set(0),this.exponent=2,this.generationCallback=(i,t,n,a)=>{Br.setFromSpherical(i);const r=Br.y*.5+.5;a.lerpColors(this.bottomColor,this.topColor,r**this.exponent)}}copy(e){return super.copy(e),this.topColor.copy(e.topColor),this.bottomColor.copy(e.bottomColor),this}}class zo extends qt{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 Oo extends qt{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>
+
+				${as}
+
+				void main() {
+
+					vec3 rayDirection = equirectUvToDirection( vUv );
+					rayDirection.x *= flipEnvMap;
+					gl_FragColor = textureCube( envMap, rayDirection );
+
+				}`}),this.depthWrite=!1,this.depthTest=!1}}class kr{constructor(e){this._renderer=e,this._quad=new et(new Oo)}generate(e,i=null,t=null){if(!e.isCubeTexture)throw new Error("CubeToEquirectMaterial: Source can only be cube textures.");const n=e.images[0],a=this._renderer,r=this._quad;i===null&&(i=4*n.height),t===null&&(t=2*n.height);const c=new pt(i,t,{type:q,colorSpace:n.colorSpace}),l=n.height,h=Math.log2(l)-2,f=1/l,u=1/(3*Math.max(Math.pow(2,h),7*16));r.material.defines.CUBEUV_MAX_MIP=`${h}.0`,r.material.defines.CUBEUV_TEXEL_WIDTH=u,r.material.defines.CUBEUV_TEXEL_HEIGHT=f,r.material.uniforms.envMap.value=e,r.material.uniforms.flipEnvMap.value=e.isRenderTargetTexture?1:-1,r.material.needsUpdate=!0;const o=a.getRenderTarget(),m=a.autoClear;a.autoClear=!0,a.setRenderTarget(c),r.render(a),a.setRenderTarget(o),a.autoClear=m;const v=new Uint16Array(i*t*4),b=new Float32Array(i*t*4);a.readRenderTargetPixels(c,0,0,i,t,b),c.dispose();for(let p=0,g=b.length;p<g;p++)v[p]=ve.toHalfFloat(b[p]);const d=new ae(v,i,t,V,fe);return d.minFilter=Hs,d.magFilter=le,d.wrapS=_e,d.wrapT=_e,d.mapping=Vr,d.needsUpdate=!0,d}dispose(){this._quad.dispose()}}function No(s){return s.extensions.get("EXT_float_blend")}const Ke=new L;class Go{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 ba,this._pathTracer=new Dr(e),this._queueReset=!1,this._clock=new Us,this._compilePromise=null,this._lowResPathTracer=new Dr(e),this._lowResPathTracer.tiles.set(1,1),this._quad=new et(new zo({map:null,transparent:!0,blending:vt,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(1024,1024),this.rasterizeSceneCallback=(i,t)=>{this._renderer.render(i,t)},this.renderToCanvasCallback=(i,t,n)=>{const a=t.autoClear;t.autoClear=!1,n.render(t),t.autoClear=a},this.setScene(new Ws,new Ur)}setBVHWorker(e){this._generator.setBVHWorker(e)}setScene(e,i,t={}){e.updateMatrixWorld(!0),i.updateMatrixWorld();const n=this._generator;if(n.setObjects(e),this._buildAsync)return n.generateAsync(t.onProgress).then(a=>this._updateFromResults(e,i,a));{const a=n.generate();return this._updateFromResults(e,i,a)}}setSceneAsync(...e){this._buildAsync=!0;const i=this.setScene(...e);return this._buildAsync=!1,i}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,i=this._renderer,t=this._materials,n=this.textureSize,a=Wa(t);e.textures.setTextures(i,a,n.x,n.y),e.materials.updateFrom(t,a),this.reset()}updateLights(){const e=this.scene,i=this._renderer,t=this._pathTracer.material,n=Va(e),a=Ua(n);t.lights.updateFrom(n,a),t.iesProfiles.setTextures(i,a),this.reset()}updateEnvironment(){const e=this.scene,i=this._pathTracer.material;if(this._internalBackground&&(this._internalBackground.dispose(),this._internalBackground=null),i.backgroundBlur=e.backgroundBlurriness,i.backgroundIntensity=e.backgroundIntensity??1,i.backgroundRotation.makeRotationFromEuler(e.backgroundRotation).invert(),e.background===null)i.backgroundMap=null,i.backgroundAlpha=0;else if(e.background.isColor){this._colorBackground=this._colorBackground||new ko(16);const t=this._colorBackground;t.topColor.equals(e.background)||(t.topColor.set(e.background),t.bottomColor.set(e.background),t.update()),i.backgroundMap=t,i.backgroundAlpha=1}else if(e.background.isCubeTexture){if(e.background!==this._previousBackground){const t=new kr(this._renderer).generate(e.background);this._internalBackground=t,i.backgroundMap=t,i.backgroundAlpha=1}}else i.backgroundMap=e.background,i.backgroundAlpha=1;if(i.environmentIntensity=e.environmentIntensity??1,i.environmentRotation.makeRotationFromEuler(e.environmentRotation).invert(),this._previousEnvironment!==e.environment)if(e.environment!==null)if(e.environment.isCubeTexture){const t=new kr(this._renderer).generate(e.environment);i.envMapInfo.updateFrom(t)}else i.envMapInfo.updateFrom(e.environment);else i.environmentIntensity=0;this._previousEnvironment=e.environment,this._previousBackground=e.background,this.reset()}_updateFromResults(e,i,t){const{materials:n,geometry:a,bvh:r,bvhChanged:c}=t;this._materials=n;const h=this._pathTracer.material;return c&&(h.bvh.updateFrom(r),h.attributesArray.updateFrom(a.attributes.normal,a.attributes.tangent,a.attributes.uv,a.attributes.color),h.materialIndexAttribute.updateFrom(a.attributes.materialIndex)),this._previousScene=e,this.scene=e,this.camera=i,this.updateCamera(),this.updateMaterials(),this.updateEnvironment(),this.updateLights(),t}renderSample(){const e=this._lowResPathTracer,i=this._pathTracer,t=this._renderer,n=this._clock,a=this._quad;this._updateScale(),this._queueReset&&(i.reset(),e.reset(),this._queueReset=!1,a.material.opacity=0,n.start());const r=n.getDelta()*1e3,c=n.getElapsedTime()*1e3;if(!this.pausePathTracing&&this.enablePathTracing&&this.renderDelay<=c&&!this.isCompiling&&i.update(),i.alpha=i.material.backgroundAlpha!==1||!No(t),e.alpha=i.alpha,this.renderToCanvas){const l=this._renderer,h=this.minSamples;if(c>=this.renderDelay&&this.samples>=this.minSamples&&(this.fadeDuration!==0?a.material.opacity=Math.min(a.material.opacity+r/this.fadeDuration,1):a.material.opacity=1),!this.enablePathTracing||this.samples<h||a.material.opacity<1){if(this.dynamicLowRes&&!this.isCompiling){e.samples<1&&(e.material=i.material,e.update());const f=a.material.opacity;a.material.opacity=1-a.material.opacity,a.material.map=e.target.texture,a.render(l),a.material.opacity=f}(!this.dynamicLowRes&&this.rasterizeScene||this.dynamicLowRes&&this.isCompiling)&&this.rasterizeSceneCallback(this.scene,this.camera)}this.enablePathTracing&&a.material.opacity>0&&(a.material.opacity<1&&(a.material.blending=this.dynamicLowRes?Vs:Wr),a.material.map=i.target.texture,this.renderToCanvasCallback(i.target,l,a),a.material.blending=vt)}}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(Ke);const e=Math.floor(this.renderScale*Ke.x),i=Math.floor(this.renderScale*Ke.y);if(this._pathTracer.getSize(Ke),Ke.x!==e||Ke.y!==i){const t=this.lowResScale;this._pathTracer.setSize(e,i),this._lowResPathTracer.setSize(Math.floor(e*t),Math.floor(i*t))}}}}const Yo=""+new URL("lightroom_14b_512-nMVn659Q.hdr",import.meta.url).href;export{ue as E,et as F,Yo as H,Di as M,jo as O,qo as P,Go as W,Vo as r};