precision mediump float; uniform float time; uniform vec2 mouse; uniform vec2 resolution; #define NEAR 3.0 #define FAR 50.0 #define EPS 1e-3 #define EPS_N 1e-4 #define RAYSTEP 32 #define AOSTEP 2 #define MENGER_ITER 3 #define REFLECTION_ITER 2 #define PI 3.1415926535 #define DIRECTIONAL_LIGHT normalize(vec3(-1.0, 1.0, -1.0)) * 1.0 // util #define saturate(x) clamp(x, 0.0, 1.0) #define calcNormal(p, dFunc) normalize(vec2(EPS_N, -EPS_N).xyy * dFunc(p + vec2(EPS_N, -EPS_N).xyy)+ vec2(EPS_N, -EPS_N).yyx * dFunc(p + vec2(EPS_N, -EPS_N).yyx ) + vec2(EPS_N, -EPS_N).yxy * dFunc(p + vec2(EPS_N, -EPS_N).yxy) + vec2(EPS_N, -EPS_N).xxx * dFunc(p + vec2(EPS_N, -EPS_N).xxx)) vec3 hsv(float h, float s, float v){ vec4 t = vec4(1.0, 2.0 / 3.0, 1.0 / 3.0, 3.0); vec3 p = abs(fract(vec3(h) + t.xyz) * 6.0 - vec3(t.w)); return v * mix(vec3(t.x), clamp(p - vec3(t.x), 0.0, 1.0), s); } float noise(float _a){ float _noise = sin(-time*3.45+_a); return _noise; } mat2 rot(float a) { float c = cos(a), s = sin(a); return mat2(c,s,-s,c); } vec2 pmod(vec2 p, float r) { float a = atan(p.x, p.y) + PI/r; float n = PI*2.0 / r; a = floor(a/n)*n; return p*rot(-a); } float sdPlane(vec3 p, float y){ return p.y - y; } float sdBox(vec3 p, vec3 s){ vec3 d = abs(p) - s; return max(max(d.x,d.y), d.z); } float sdSphere(vec3 p, float r){ return length(p) - r; } float sdCylinder( vec3 p, vec3 c ) { return length(p.xz-c.xy)-c.z; } float sdCappedCylinder( vec3 p, vec2 h ) { vec2 d = abs(vec2(length(p.xz),p.y)) - h; return min(max(d.x,d.y),0.0) + length(max(d,0.0)); } float dBar(vec2 p, float width) { vec2 d = abs(p) - width; return min(max(d.x, d.y), 0.0) + length(max(d, 0.0)) + 0.01 * width; } float dCrossBar(vec3 p, float x) { float bar_x = dBar(p.yz, x); float bar_y = dBar(p.zx, x); float bar_z = dBar(p.xy, x); return min(bar_z, min(bar_x, bar_y)); } float PolarCrossBar(vec3 p, float s){ float d = FAR; for(int i = 0; i < 6; i++){ p = abs(p); p.xy *= rot(time*0.1); p.xz *= rot(time*0.07); p.yz *= rot(time*0.05); d = min(d, dCrossBar(mod(p, 0.5) - 0.25, s)); s *= 0.8; } return d; } float sdMenger(vec3 p, float s){ float d = sdBox(p, vec3(s)); for(int i = 0; i < MENGER_ITER; i++){ d = max(d, -dCrossBar(mod(p + s / 2.0, s * 2.0) - s / 2.0, s / 3.0)); s /= 3.0; } return d; } float PolarMenger(vec3 p, float s){ float d = FAR; p = abs(p); p.xy *= rot(time*0.1); p.xz *= rot(time*0.13); p.yz *= rot(time*0.18); d = min(d, sdMenger(mod(p, s * 1.5) - s * 0.75, s)); return d; } float sdMemex(vec3 p, vec2 s){ p.yz *= rot(PI / 2.0); float d = sdCappedCylinder(p - vec3(0.0, 0.0, s.x / 2.0), vec2( s.x * 1.0, s.y)); d = max(d, sdCappedCylinder(p - vec3(0.0, 0.0, -s.x / 2.0), vec2( s.x * 1.25, s.y))); d = max(d, -sdCappedCylinder(p - vec3(0.0, 0.0, s.x * 0.125), vec2( s.x * 0.63, s.y * 2.0))); s.x /= 6.0; p.z += 0.5; for(int i = 0; i < MENGER_ITER; i++){ d = max(d, -dCrossBar(mod(p + s.x / 2.0, s.x * 2.0) - s.x / 2.0, s.x / 3.0)); s.x /= 3.0; } return d; } float map(vec3 p){ float d = FAR; d = min(d, sdPlane(p, sin(-time*8.0+length(p)*2.0) * 0.1 - 1.5)); d = min(d, PolarMenger(vec3(pmod(p.xy, 6.0), p.z), 26.0)); p.xy *= rot(time*0.125); p.xz *= rot(time*0.07); p.yz *= rot(time*0.03); p -= vec3(0.0, 0.5, 0.0); d = min(d, sdMemex(p, vec2(2.5, 0.5))); return d; } float calcAO(vec3 p, vec3 n, float d, float k){ float occ = 0.0; for(int i = 0; i < AOSTEP; i++){ float len = float(i+1) / float(AOSTEP) * d; float d = map(n * len + p); occ += (len - d) * k; k *= 0.5; } return clamp(1.0 - occ, 0.0, 1.0);//saturate(1.0 - occ); } float sdPointLight(vec3 p, vec3 lightPos){ float d = FAR; p.xy *= rot(time*0.25); d = min(d, sdSphere(p - lightPos, 0.1)); p.xz *= rot(120.0*PI/180.0); d = min(d, sdSphere(p - lightPos, 0.1)); p.xz *= rot(120.0*PI/180.0); d = min(d, sdSphere(p - lightPos, 0.1)); return d; } float calcBloom(vec3 lightPos, vec3 cam, vec3 ray, float endDist){ float depth = 0.0; float bloom = 0.0; for(int i=0; i endDist) { break; } vec3 pos = cam + ray * depth; float d = sdPointLight(pos, lightPos); bloom += exp(-d*0.25); depth+=d; } return bloom / float(RAYSTEP); } vec3 frame(vec3 cameraPos, vec3 rayDirection, out vec3 rayPos, out vec3 normal, out float depth, out bool hit){ vec3 col = vec3(0.0);// = vec3(0.25 - abs(0.22 - uv.y), 0.2 - abs(0.25 - uv.y), 0.5 - abs(0.28 - uv.y)); float dist = 0.0; for(int i = 0; i < RAYSTEP; i++){ rayPos = cameraPos + rayDirection * depth; dist = map(rayPos); if(dist < EPS || depth > FAR || i == RAYSTEP - 1){ hit = dist < EPS; break; } depth += dist; } if(!hit){ col = vec3(0.0); return col; } // color float dAmount = smoothstep(FAR, NEAR, depth); // menger if(abs(dist - PolarMenger(vec3(pmod(rayPos.xy, 6.0), rayPos.z), 26.0)) < EPS){ col = vec3(dAmount) * vec3(0.1, 0.15, 0.6) * hsv(dAmount + time * 0.25, 0.2, 0.1); col += vec3(0.2, 0.5, 1.0) * pow(clamp(1.0 - abs(mod(depth - time * 5.0, 6.0) - 3.0), 0.0, 1.0), 5.0); // water }else if(abs(dist - sdPlane(rayPos, sin(-time*8.0+length(rayPos)*2.0) * 0.1 - 1.5)) < EPS){ col = vec3(dAmount) * vec3(0.0, 0.15, 0.12); // other }else{ col = vec3(dAmount) * hsv(depth*0.3-time*0.5, 0.5, 0.5); } // normal normal = calcNormal(rayPos, map); // Ambient Occlusion float ao = calcAO(rayPos, normal, 0.25, 2.5); col *= ao; // diffuse float diff = max(dot(normal, DIRECTIONAL_LIGHT) * 1.5, 0.5); col *= diff; // emission //near = min(near, sdPointLight(rayPos, pointLight)); return col; } vec3 canvas(vec2 uv){ vec3 cameraPos = vec3(0.0, sin(time*0.25) * 0.2, -6.0 + sin(time*0.333) * 1.0); cameraPos.xz *= rot(time*0.25); float sceneZ = 2.0; vec3 rayDirection = normalize(vec3(uv.x, uv.y, sceneZ)); rayDirection.xz *= rot(time*0.25); vec3 pointLight = vec3(sin(time) * 3.0 + exp(cos(time * 1.67)), sin(time*2.0) * 0.5, cos(time) * 3.0 + exp(sin(time * 1.34))); float near = FAR; float bloom = 0.0; vec3 light = normalize(vec3(-1.0, 1.0, -1.0)) * 1.0; vec3 rayPos = cameraPos; vec3 col = vec3(0.0); float alpha = 1.0; // reflection for(int j = 0; j < REFLECTION_ITER; j++){ vec3 normal = vec3(0.0); bool hit = false; float depth = 0.0; // color col += alpha * frame(cameraPos, rayDirection, rayPos, normal, depth, hit); // emission //col += vec3(1.0, 0.5, 0.7) * exp(-near * 0.25); // pointlight col += vec3(1.0, 0.5, 0.7) * calcBloom(pointLight, cameraPos, rayDirection, depth) * 5.0; // reflection alpha *= 0.5; rayDirection = normalize(reflect(rayDirection, normal)); cameraPos = rayPos + normal * EPS; } return col; } void main( void ) { vec2 uv = (gl_FragCoord.xy * 2.0 - resolution.xy) / min(resolution.x, resolution.y); float impact = max(smoothstep(1.5, 2.0, sin(time * 3.0)+1.0), 0.5); vec3 col = vec3(canvas(uv - impact * 0.04 * uv).x, canvas(uv - impact * 0.08 * uv).y, canvas(uv - impact * 0.12 * uv + step(mod(uv.y + time * 0.4, 0.2), sin(time+uv.y*5.0)*0.05+0.06) * vec2(0.05,0.0)).z); //col = pow(col, vec3(1.2)); col *= 1.0 - smoothstep(0.2, 1.5, length((gl_FragCoord.xy * 2.0 - resolution.xy) / max(resolution.x, resolution.y))) * 1.0; col *= max(step(sin(uv.y * 500.0 + mod(time*1.5, 2.78)), 0.0), 0.7); gl_FragColor = vec4( col, 1.0 ); }