precision highp float; uniform vec2 resolution; uniform float time; uniform vec2 mouse; uniform sampler2D backbuffer; // Created by Danil (2021+) https://twitter.com/AruGL // License Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. // self https://www.shadertoy.com/view/Nt2Szm // this is my port to WebGL1 on neort neort.io/art/c4g6ors3p9ffolj068ag vec4 m4x(mat4 mx, int idx){ if(idx==0)return mx[0]; if(idx==1)return mx[1]; if(idx==2)return mx[2]; if(idx==3)return mx[3]; return mx[0]; } vec3 m3x(mat3 mx, int idx){ if(idx==0)return mx[0]; if(idx==1)return mx[1]; if(idx==2)return mx[2]; return mx[0]; } float v4x(vec4 v, int idx){ if(idx==0)return v[0]; if(idx==1)return v[1]; if(idx==2)return v[2]; if(idx==3)return v[3]; return v[0]; } int iv4x(ivec4 v, int idx){ if(idx==0)return v[0]; if(idx==1)return v[1]; if(idx==2)return v[2]; if(idx==3)return v[3]; return v[0]; } float v3x(vec3 v, int idx){ if(idx==0)return v[0]; if(idx==1)return v[1]; if(idx==2)return v[2]; return v[0]; } void setm3x(inout mat3 mx, int idx, vec3 v){ if(idx==0)mx[0]=v; else if(idx==1)mx[1]=v; else if(idx==2)mx[2]=v; } void setm4x(inout mat4 mx, int idx, vec4 v){ if(idx==0)mx[0]=v; else if(idx==1)mx[1]=v; else if(idx==2)mx[2]=v; else if(idx==3)mx[3]=v; } void setv4x(inout vec4 v, int idx, float f){ if(idx==0)v[0]=f; else if(idx==1)v[1]=f; else if(idx==2)v[2]=f; else if(idx==3)v[3]=f; } void setiv4x(inout ivec4 v, int idx, int f){ if(idx==0)v[0]=f; else if(idx==1)v[1]=f; else if(idx==2)v[2]=f; else if(idx==3)v[3]=f; } float fresnel(vec3 d, vec3 n) { float a = clamp(1.0-dot(n,-d), 0.0, 1.0); return clamp(exp((5.0*a)-5.0), 0.0, 1.0); } // ANGLE bug https://bugs.chromium.org/p/chromium/issues/detail?id=1238461 // ANGLE compile time slowdown 20x times because of single call of (0.5 * t * t * t) below float easeInOutCubic(float t) { if ((t *= 2.0) < 1.0) { //return 0.5 * t * t * t; // uncomment for bug return 0.5 * pow(abs(t),3.); //fix } else { return 0.5 * ((t -= 2.0) * t * t + 2.0); } } float hash21(vec2 p) { vec3 p3 = fract(vec3(p.xyx) * .1031); p3 += dot(p3, p3.yzx + 31.433); return fract((p3.x + p3.y) * p3.z); } float hash31(vec3 p3) { p3 = fract(p3 * .1031); p3 += dot(p3, p3.zyx + 13.32); return fract((p3.x + p3.y) * p3.z); } float get_rand(vec2 pos,vec2 box){ float tps = hash21(pos); tps = tps*0.1+0.9*smoothstep(0.2,.9,tps); tps *= min((abs(pos.x-box.x-3.)/box.x)*((pos.y+3.)/(box.y*2.)),1.); return 0.05+0.95*tps; } // https://www.shadertoy.com/view/XsX3zB vec3 random3(vec3 c) { float j = 4096.0*sin(dot(c,vec3(17.0, 59.4, 15.0))); vec3 r;r.z = fract(512.0*j);j *= .125;r.x = fract(512.0*j);j *= .125;r.y = fract(512.0*j); return r-0.5; } const float F3 = 0.3333333; const float G3 = 0.1666667; /* 3d simplex noise */ float simplex3d(vec3 p) { vec3 s = floor(p + dot(p, vec3(F3))); vec3 x = p - s + dot(s, vec3(G3)); vec3 e = step(vec3(0.0), x - x.yzx); vec3 i1 = e*(1.0 - e.zxy); vec3 i2 = 1.0 - e.zxy*(1.0 - e); vec3 x1 = x - i1 + G3; vec3 x2 = x - i2 + 2.0*G3; vec3 x3 = x - 1.0 + 3.0*G3; vec4 w, d; w.x = dot(x, x); w.y = dot(x1, x1); w.z = dot(x2, x2); w.w = dot(x3, x3); w = max(0.6 - w, 0.0); d.x = dot(random3(s), x); d.y = dot(random3(s + i1), x1); d.z = dot(random3(s + i2), x2); d.w = dot(random3(s + 1.0), x3); w *= w; w *= w; d *= w; return dot(d, vec4(52.0)); } const mat3 rot1 = mat3(-0.37, 0.36, 0.85,-0.14,-0.93, 0.34,0.92, 0.01,0.4); const mat3 rot2 = mat3(-0.55,-0.39, 0.74, 0.33,-0.91,-0.24,0.77, 0.12,0.63); const mat3 rot3 = mat3(-0.71, 0.52,-0.47,-0.08,-0.72,-0.68,-0.7,-0.45,0.56); float simplex3d_fractal(vec3 m) { return 0.5333333*simplex3d(m*rot1) +0.2666667*simplex3d(2.0*m*rot2) +0.1333333*simplex3d(4.0*m*rot3) +0.0666667*simplex3d(8.0*m); } #define MAX_DIST 1000. #define MIN_DIST .001 // ffix float precision fix, cy height shift animation, stimer shift timer 0-1 // ffix also used in reflections to control voxel size to remove white borders on glow reflections bool VoxelsIntersect_30(vec2 box, float max_h, float ffix, in vec3 ro, in vec3 rd, out vec3 normal , out float tnew, out vec2 idx, inout float cy, float stimer) { vec2 pos = floor(ro.xz); vec3 rdi = 1./rd; vec3 rda = abs(rdi); vec3 rds = sign(rd); vec2 dis = (pos - ro.xz + .5 + rds.xz*.5) * rdi.xz; vec3 roi = rdi*(ro-vec3(.5,0.,.5)); vec2 mm = vec2(0.0); for( int i=0; i<30; i++ ) { float tps = get_rand(pos, box); float tm=1.; if (any(greaterThan(pos,box*2.))||any(lessThan(pos,vec2(1.))))tm=0.; cy = mix(-(tm*max_h*tps-1.)+max_h*2.-2.,tm*max_h*tps-1.,stimer); vec3 n = roi - rdi * vec3(pos.x, cy, pos.y); vec3 k = rda*vec3(.5, -0.00099+tm*max_h*tps, .5); vec3 t1 = -n - k; vec3 t2 = -n + k; float tN = max( max( t1.x, t1.y ), t1.z )-ffix; float tF = min( min( t2.x, t2.y ), t2.z ); if ( tN < tF && tN>MIN_DIST && tNMIN_DIST && tN0.) m+=(mousex.xy/resx.y-0.5*resx.xy/resx.y)*2.5; m.y = -m.y; float fov=70.; float aspect = resx.x / resx.y; float screenSize = (1.0 / (tan(((180.-fov)* (3.1415926 / 180.0)) / 2.0))); vec3 rd = normalize(vec3(uv*screenSize, 1./aspect)); mat3 rotX = mat3(1.0, 0.0, 0.0, 0.0, cos(m.y), sin(m.y), 0.0, -sin(m.y), cos(m.y)); mat3 rotY = mat3(cos(m.x), 0.0, -sin(m.x), 0.0, 1.0, 0.0, sin(m.x), 0.0, cos(m.x)); rd = (rotY * rotX) * rd; ray_ro=ro; ray_rd=rd; } bool boxAABB(in vec3 dims, vec3 ro, vec3 rd) { vec3 n = ro / rd; vec3 k = dims / abs(rd); vec3 t1 = -k - n, t2 = k - n; float tN = max(max(t1.x, t1.y), t1.z); float tF = min(min(t2.x, t2.y), t2.z); return tN < tF && tF > 0.0; } void VoxelsIntersectMin_ref(vec2 box, float max_h , vec3 posl,vec3 ro, vec3 rd, float rfnl, float sky_pow, float col_pow, inout float refc, inout float tHitInfo_t, inout vec3 tHitInfo_norm, inout vec4 tHitInfo_color, inout int tHitInfo_obj_type, inout vec2 tHitInfo_exdata, float atimer) { float tnew; vec2 idx; vec3 normnew; float cy; if(VoxelsIntersect_6(box, max_h, -0.001, ro, rd, normnew, tnew, idx, cy, atimer)){ bool in_box=(idx.x<=box.x*2.)&&(idx.y<=box.y*2.); in_box=in_box&&all(greaterThan(idx,vec2(0.))); float fnl = fresnel(rd,normnew); if(in_box){ float alb = hash21(idx+1.); float trnd = get_rand(idx+20., box*0.7); float emi = step(0.65,trnd-0.75*trnd*hash21(idx+5.)); alb*=0.25; alb *= rfnl * (1.-fnl) * (.8 + .2 * dot(normnew, vec3(-.25916,.8639,-.4319))); emi *= col_pow; vec3 tuv = fract(ro+tnew*rd) * 2.0 - 1.0; float vig=max(tuv.x*tuv.x, tuv.z*tuv.z)*0.5; vec3 tc = white+blue; tc=mix(tc,tc*.65,vig); vec3 col=max(alb*white,emi*1.175*tc); tHitInfo_color.rgb += col; tHitInfo_color.a += emi*0.25; refc*=0.12; } else { vec3 bg=clamp(vec3(.5,.8,1)*0.+1.*blurred_background(rd),0.,1.); tHitInfo_color.rgb += sky_pow*clamp(bg * (rfnl*(5.-2.5*rd.y)),vec3(0.),bg); } } else { vec3 bg=clamp(vec3(.5,.8,1)*0.+1.*blurred_background(rd),0.,1.); tHitInfo_color.rgb += sky_pow*clamp(bg * (rfnl*(5.-2.5*rd.y)),vec3(0.),bg); } } void VoxelsIntersectMin(vec2 box, float max_h , vec3 posl,vec3 ro, vec3 rd, inout bool result, vec3 sky, inout float tHitInfo_t, inout vec3 tHitInfo_norm, inout vec4 tHitInfo_color, inout int tHitInfo_obj_type, inout vec2 tHitInfo_exdata, float atimer) { float tnew; vec2 idx; vec3 normnew; posl.xz+=1.; ro+=posl; if(!boxAABB(vec3(box.x,max_h,box.y), ro+vec3(-1.,-max_h+1.,-1.), rd))return; ro.xz+=box; float vros = length(ro.xz-(box+1.)); if(VoxelsIntersect_30(box, max_h, 0.0001, ro, rd, normnew, tnew, idx, tHitInfo_exdata.y, atimer)){ bool in_box=(idx.x<=box.x*2.)&&(idx.y<=box.y*2.); in_box=in_box&&all(greaterThan(idx,vec2(0.))); if (in_box&&(tnew>0.)) { float tmx = 0.; if((tHitInfo_obj_type==OBJ_BALL)&&(tnew>tHitInfo_t))tmx=clamp(dot(tHitInfo_color.rgb,vec3(1.))*0.25,0.,1.); tHitInfo_color = vec4(0.,0.,0.,0.); float fnl = fresnel(rd,normnew); float alb = hash21(idx+1.); float trnd = get_rand(idx+20., box*0.7); float emi = step(0.65,trnd-0.75*trnd*hash21(idx+5.)); vec3 col = vec3(0.); alb*=0.1*alb; alb *= (1.-fnl) * (.8 + .2 * dot(normnew, vec3(-.25916,.8639,-.4319))); vec3 tuv = fract(ro+tnew*rd) * 2.0 - 1.0; float vig=max(tuv.x*tuv.x, tuv.z*tuv.z)*0.5; vec3 tc = white+blue; tc=mix(tc,tc*.85,vig); col=max(alb*white,emi*2.35*(tc)); vec3 skyo = background_sky(rd); tHitInfo_color.rgb += mix(col,skyo,fnl*fnl); tHitInfo_color.a += emi; tHitInfo_obj_type = OBJ_FLOOR; result = true; tHitInfo_t = tnew; tHitInfo_norm = normnew; float vrds = length(tnew * rd.xz); if((emi<0.5)&&(fnl > .001)){ tHitInfo_obj_type = OBJ_REFLOOR; tHitInfo_norm = normnew; tHitInfo_exdata.x = fnl; } tHitInfo_color.rgb = mix(tHitInfo_color.rgb,sky,smoothstep(24.,27.5,vrds)); tHitInfo_color.rgb = mix(tHitInfo_color.rgb,sky,tmx); } } } float glow_sphere(vec3 rd, vec3 n){ float a = pow(0.7122 + max(dot(n, reflect(rd,n)),0.), 10.85); float b = pow(0.7122 + max(1.-dot(n, reflect(rd,n)),0.), 3.85); return clamp(min(a,b),0.,10.); } void SphereIntersectMin(vec3 SpPos, float SpRad, vec3 ro, vec3 rd, inout bool result, bool refl, bool sph2, float ot, inout float tHitInfo_t, inout vec3 tHitInfo_norm, inout vec4 tHitInfo_color, inout int tHitInfo_obj_type, inout vec2 tHitInfo_exdata) { float tnew; vec3 normnew; if (SphereIntersectInv(SpPos, SpRad, ro, rd, tnew, normnew)) { if(!refl){ if (tnew < tHitInfo_t||sph2&&(tnew < ot)) { tHitInfo_t = min(tnew,tHitInfo_t); tHitInfo_norm = normnew; if(sph2){ vec3 spotLightOrigin = (SpPos - vec3(0.,0.,1.)*22.9); vec3 Ln = normalize(spotLightOrigin-(ro+rd*tnew)); tHitInfo_color.rgb += color_spiral(Ln,20.)*4.; } else tHitInfo_color.rgb += glow_sphere(rd,-normnew)*(sqrt(blue))*3.; tHitInfo_obj_type = OBJ_BALL; result = true; } }else{ if(sph2) { vec3 spotLightOrigin = (SpPos - vec3(0.,0.,1.)*19.9); vec3 Ln = normalize(spotLightOrigin-(ro+rd*tnew)); tHitInfo_color.a += color_spiral(Ln,10.).b*.85; } else tHitInfo_color.a = glow_sphere(rd,-normnew); tHitInfo_obj_type = OBJ_BALL; result = true; } } } void render(vec3 ray_ro, vec3 ray_rd, inout vec3 sky, inout float tHitInfo_t, inout vec3 tHitInfo_norm, inout vec4 tHitInfo_color, inout int tHitInfo_obj_type, inout vec2 tHitInfo_exdata, float atimer) { tHitInfo_t = MAX_DIST; tHitInfo_obj_type = OBJ_SKY; tHitInfo_color = vec4(sky,0.); tHitInfo_norm = vec3(0.,1.,0.); tHitInfo_exdata = vec2(0.); vec4 col = vec4(0.); vec3 ro = ray_ro; vec3 rd = ray_rd; bool result = false; float angle=0.; vec3 spherex=vec3(0.,8.15,8.5+voxbox.y+15.); vec3 spherex2=vec3(0.,8.15,8.5+voxbox.y+15.-22.25); rd += 0.0001 * (1.0 - abs(sign(rd))); float ot=tHitInfo_t; SphereIntersectMin(spherex, 13.5, ro, rd, result, false, false,0., tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata); SphereIntersectMin(spherex2, 23.5, ro, rd, result, false, true,ot, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata); sky=tHitInfo_color.rgb; VoxelsIntersectMin(voxbox, voxmax_h , voxposl,ro, rd, result, sky, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata, atimer); } // reflections void render_refl(vec3 ray_ro, vec3 ray_rd, float fbmx, inout vec3 sky, inout float tHitInfo_t, inout vec3 tHitInfo_norm, inout vec4 tHitInfo_color, inout int tHitInfo_obj_type, inout vec2 tHitInfo_exdata, float atimer) { vec3 skyo = background_sky(ray_rd); if(tHitInfo_obj_type == OBJ_REFLOOR){ vec3 ro = ray_ro; ro = ro+ray_rd*tHitInfo_t; float vrds = length(tHitInfo_t * ray_rd.xz); fbmx*=1.-smoothstep(10.,14.5,vrds); float gfnl = fresnel((ray_rd),tHitInfo_norm); fbmx*=max(1.-gfnl*2.,0.); vec3 ord = reflect(ray_rd-(.0295*fbmx*tHitInfo_norm), tHitInfo_norm); vec3 rd = reflect(ray_rd-(.295*fbmx*tHitInfo_norm), tHitInfo_norm); rd += 0.0001 * (1.0 - abs(sign(rd))); bool result = false; float angle=0.1; vec3 spherex=vec3(0.,10.15,8.5+voxbox.y+15.); vec3 spherex2=vec3(0.,10.15-2.,8.5+voxbox.y+15.-22.25); float oemi=tHitInfo_color.a; float ot=tHitInfo_t; SphereIntersectMin(spherex, 13.5, ro, rd, result, true, false,0., tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata); SphereIntersectMin(spherex2, 23.5, ro, ord, result, true, true,ot, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata); float refc=0.; float fnl = 0.; if(tHitInfo_obj_type == OBJ_BALL){ fnl=1.;refc=tHitInfo_color.a;tHitInfo_color.a=oemi; } ro+=voxposl;ro.xz+=voxbox+1.; rd = reflect(ray_rd+(.695*fbmx*tHitInfo_norm), tHitInfo_norm); rd += 0.0001 * (1.0 - abs(sign(rd))); VoxelsIntersectMin_ref(voxbox, voxmax_h , voxposl+1.,ro, rd, tHitInfo_exdata.x, 01.5, 01.+0.25*fbmx,refc, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata, atimer); tHitInfo_color.rgb = 0.985*tHitInfo_color.rgb+0.015*tHitInfo_color.rgb*fbmx; fnl = max(fresnel(rd,tHitInfo_norm),fnl); tHitInfo_color.rgb = mix(tHitInfo_color.rgb,skyo,1.-sqrt(fnl))+refc*(sqrt(blue))*3.; tHitInfo_color.rgb = mix(tHitInfo_color.rgb,sky,smoothstep(24.,27.5,vrds)); } } void mainImage(out vec4 fragColor, in vec2 fragCoord, vec4 iMouse, vec2 iResolution, float iTime) { float atimer=1.; vec3 ret_col = vec3(0.0); float emi = 0.; vec3 tro = camera_anim(iTime, atimer); vec2 m=get_look(iTime); #if AA>1 const int ta=AA*AA; #else const int ta = 1; #endif vec4 HitInfo_t; mat4 HitInfo_norm; mat4 HitInfo_color; ivec4 HitInfo_obj_type; mat4 HitInfo_exdata; mat4 cam_ro; mat4 cam_rd; mat4 sky; { #if AA>1 for( int mx=0; mx1 int idx=mx*AA+nx; #else int idx=0; #endif vec3 tr;vec3 td; SetCamera(uv,tro,m,tr,td, iMouse, iResolution.xy); setm4x(cam_ro,idx,vec4(tr,0.)); setm4x(cam_rd,idx,vec4(td,0.)); setm4x(sky,idx,vec4(background_sky(m4x(cam_rd,idx).xyz),0.)); vec3 tsky=m4x(sky,idx).rgb; float tHitInfo_t=v4x(HitInfo_t,idx);vec3 tHitInfo_norm=m4x(HitInfo_norm,idx).xyz; vec4 tHitInfo_color=m4x(HitInfo_color,idx); int tHitInfo_obj_type=iv4x(HitInfo_obj_type,idx); vec2 tHitInfo_exdata=m4x(HitInfo_exdata,idx).xy; render(m4x(cam_ro,idx).xyz,m4x(cam_rd,idx).xyz,tsky, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata, atimer); setv4x(HitInfo_t,idx,tHitInfo_t); setm4x(HitInfo_norm,idx,vec4(tHitInfo_norm,0.)); setm4x(HitInfo_color,idx,tHitInfo_color); setiv4x(HitInfo_obj_type,idx,tHitInfo_obj_type); setm4x(HitInfo_exdata,idx,vec4(tHitInfo_exdata,0.,0.)); setm4x(sky,idx,vec4(tsky,0.)); #if AA>1 } #endif } vec3 ro = cam_ro[0].xyz; ro+=voxposl;ro.xz+=voxbox+1.;ro = ro+cam_rd[0].xyz*HitInfo_t.t; float fbmx = (simplex3d_fractal(5.*(ro-vec3(0.,HitInfo_exdata[0].y,0.))+vec3(0.,.5,1.)*iTime)); { #if AA>1 for( int mx=0; mx1 int idx=mx*AA+nx; #else int idx=0; #endif vec3 tsky=m4x(sky,idx).rgb; float tHitInfo_t=v4x(HitInfo_t,idx);vec3 tHitInfo_norm=m4x(HitInfo_norm,idx).xyz; vec4 tHitInfo_color=m4x(HitInfo_color,idx); int tHitInfo_obj_type=iv4x(HitInfo_obj_type,idx); vec2 tHitInfo_exdata=m4x(HitInfo_exdata,idx).xy; render_refl(m4x(cam_ro,idx).xyz,m4x(cam_rd,idx).xyz,fbmx,tsky, tHitInfo_t, tHitInfo_norm, tHitInfo_color, tHitInfo_obj_type, tHitInfo_exdata, atimer); setv4x(HitInfo_t,idx,tHitInfo_t); setm4x(HitInfo_norm,idx,vec4(tHitInfo_norm,0.)); setm4x(HitInfo_color,idx,tHitInfo_color); setiv4x(HitInfo_obj_type,idx,tHitInfo_obj_type); setm4x(HitInfo_exdata,idx,vec4(tHitInfo_exdata,0.,0.)); setm4x(sky,idx,vec4(tsky,0.)); #if AA>1 } #endif } { #if AA>1 for( int mx=0; mx1 int idx=mx*AA+nx; #else int idx=0; #endif ret_col += m4x(HitInfo_color,idx).rgb; emi += m4x(HitInfo_color,idx).a; #if AA>1 } ret_col /= float(AA*AA); emi /= float(AA*AA); #endif } vec3 glow = vec3(0.7, 0.8, 1.6)*max((cam_rd[0].y*0.49 + 0.51),0.001)*0.8; ret_col += glow*glow*glow / 4.0; ret_col *= 0.425; ret_col = clamp(ret_col,0.,1.); ret_col = sqrt(ret_col); vec2 vig = (fragCoord.xy / iResolution.xy) * 2.0 - 1.0; ret_col = mix(ret_col, ret_col*ret_col, dot(vig, vig)*0.5); ret_col = mix(ret_col, ret_col*ret_col, 0.5); ret_col = clamp(ret_col,0.,1.); ret_col = pow(ret_col, vec3(01.4545)); ret_col=mix(ret_col,ret_col.bgr,atimer); fragColor = vec4(ret_col,clamp(emi,0.,1.)); fragColor.a=1.; //bloom not used } void main(void) { vec4 col; mainImage(col, gl_FragCoord.xy,vec4(-1.),resolution.xy,time); gl_FragColor = vec4(col.rgb, 1.0); }

Bright cubes