precision highp float; uniform vec2 resolution; uniform float time; uniform vec2 mouse; uniform sampler2D backbuffer; mat2 rot(float a){return mat2(cos(a),sin(a),-sin(a),cos(a));} float bo(vec3 p,vec3 s){p = abs(p)-s;return max(p.x,max(p.y,p.z));} float bar(vec2 p,vec2 s){p = abs(p)-s;return max(p.x,p.y);} float rbo(vec3 p,vec3 s,float r){p = abs(p)-s;return length(max(p,0.0)) + min(max(p.x,max(p.y,p.z)),0.0) - r;} float pi = acos(-1.); vec2 pmod(vec2 p, float r){float a = atan(p.x,p.y)+pi/r;float n = (pi*2.)/r;return p*rot(-(floor(a/n)*n));} vec2 smoothpmod(vec2 p, float r,float e) { float a = atan(p.x,p.y)+pi/r; float n = (pi*2.)/r; float angle = floor(a/n)*n; float x = a/n * n; x = a/n; angle = floor(x) + pow(fract(x),e); vec2 rp = p * rot(-( angle * n )); return rp; } float rand(vec3 p) { p = fract( p*0.3183099+.1 ); p *= 17.0; return fract( p.x*p.y*p.z*(p.x+p.y+p.z) ); } float noise(vec3 p) { vec3 fp = floor(p); vec3 fs = fract(p); vec3 sf = vec3(1.,0.,0.); float a = rand(fp); float b = rand(fp + sf); float c = rand(fp + sf.yxz); float d = rand(fp + sf.xxz); float e = rand(fp + sf.yyx); float f = rand(fp + sf.xyx); float g = rand(fp + sf.yxx); float h = rand(fp + sf.xxx); vec3 u = smoothstep(0.,1.,fs); // u = pow(fs,vec3(3.)); return mix(mix( mix(a,b,u.x), mix(c,d,u.x),u.y), mix( mix(e,f,u.x), mix(g,h,u.x),u.y),u.z ); } float fbm(vec3 p) { float v = 0.; float a = .5; float f = 0.; for(int i=0;i<3;i++) { v+=a*noise(p); p*=2.; a*=.5; } return v; } float menger(vec3 p,vec3 offset) { float scale = 1.65; vec4 z = vec4(p,1.); for(int i = 0;i < 3;i++) { //z.yz *= rot(float(i) * 2.); z = abs(z); if(z.x < z.y)z.xy = z.yx; if(z.x < z.z)z.xz = z.zx; //if(z.y < z.z)z.yz = z.zy; z *= scale; z.xyz -= offset * (scale - 1.); if(z.z < -.5 * offset.z * (scale - 1.) ) z.z += offset.z * (scale - 1.); } return (length(max(abs(z.xyz) - vec3(1.0, 1.0, 1.0), 0.0))) / z.w; } vec3 deMengerSponge3(vec3 p, vec3 offset, float scale) { vec4 z = vec4(p, 1.0); for (int i = 0; i < 2; i++) { z = abs(z); z.xz *= rot(2.); if (z.x < z.y) z.xy = z.yx; //if (z.x < z.z) z.xz = z.zx; if (z.y < z.z) z.yz = z.zy; z *= scale; z.xyz -= offset * (scale - 1.0); if (z.z < -0.5 * offset.z * (scale - 1.0)) z.z += offset.z * (scale - 1.0); } return z.xyz / z.w; } vec3 pp; mat3 RotMat(vec3 axis, float angle) { // http://www.neilmendoza.com/glsl-rotation-about-an-arbitrary-axis/ axis = normalize(axis); float s = sin(angle); float c = cos(angle); float oc = 1.0 - c; return mat3(oc * axis.x * axis.x + c, oc * axis.x * axis.y - axis.z * s, oc * axis.z * axis.x + axis.y * s, oc * axis.x * axis.y + axis.z * s, oc * axis.y * axis.y + c, oc * axis.y * axis.z - axis.x * s, oc * axis.z * axis.x - axis.y * s, oc * axis.y * axis.z + axis.x * s, oc * axis.z * axis.z + c ); } vec3 random33(vec3 st) { st = vec3(dot(st, vec3(127.1, 311.7,811.5)), dot(st, vec3(269.5, 183.3,211.91)), dot(st, vec3(511.3, 631.19,431.81)) ); return -1.0 + 2.0 * fract(sin(st) * 43758.5453123); } vec4 celler3D(vec3 i,vec3 sepc) { vec3 sep = i * sepc; vec3 fp = floor(sep); vec3 sp = fract(sep); float dist = 5.; vec3 mp = vec3(0.); for (int z = -1; z <= 1; z++) { for (int y = -1; y <= 1; y++) { for (int x = -1; x <= 1; x++) { vec3 neighbor = vec3(x, y ,z); vec3 pos = vec3(random33(fp+neighbor)); pos = sin( (pos*6. +time/2.) )* 0.5 + 0.5; vec3 pp = neighbor + pos - sp; float divs = length(neighbor + pos - sp); //divs = max(pp.x,pp.y); mp = (dist >divs)?pos:mp; dist = (dist > divs)?divs:dist; } } } return vec4(mp,dist); } vec3 gp; #define TT(x,y) floor(x) + pow(fract(x),y) vec2 map(vec3 p) { float o = 10.; float id = 0.; //p.y -= 10.5; vec3 op = p; //p.x += time * 30.; vec3 mp = p; p.x += time * 50.; //p -= deMengerSponge3(p,vec3(3.),3.); p.y += sin(time + p.x/10.) * 10.; p.xz = sin(p.xz/20.) * 20.; //p = abs(p) - 7.; p.xz *= rot(sin(time + op.x/100. + op.z/100.)); p.yz *= rot(p.x/40.); p.xz = smoothpmod(p.xz,2.,5.); p.yz = smoothpmod(p.yz,6.,2.); //p.xz = pmod(p.xz,8.); //p.yz = pmod(p.yz,8.); p.xz *= rot(time/32.); //p.xy = smoothpmod(p.xy,3.,.1); p.xz = sin(p.xz) * 5.; //p.x += pi/3. * step(abs(sin(p.x/2.)),0.); p -= clamp(p,-4.,4.); o = length(p) - 1.3; gp = p; mp.y += 45.; float t2 = TT(time/pi/2.,2.); mp.x += t2; mp.xy *= rot(t2); mp.xz *= rot(time/3.); mp.yz *= rot(time/4.); float sp = length(mp) - 62.3; //sp = menger(mp/13.,vec3(3.))*13.; //mp.xz *= rot(pi/4.); //mp.yz *= rot(pi/4.); mp -= clamp(mp,-36.,36.); sp = max(sp,(length(mp) - 3.)); if(o > sp) { o = sp; id = 1.; } return vec2(o,id); } vec2 march(vec3 cp , vec3 rd) { float depth = 0.; for(int i = 0 ; i< 128 ; i++) { vec3 rp = cp + rd * depth; vec2 d = map(rp); if(abs(d.x) < 0.01) { return vec2(depth,d.y); } if(depth > 500.)break; depth += d.x; } return vec2(-depth , 0.); } vec2 skymap(vec3 p) { float o = 10.; float id = 0.; //p.y -= 10.5; vec3 op = p; //p.x += time * 30.; vec3 mp = p; p.y += sin(time + p.x/10.) * 10.; p.xz = sin(p.xz/20.) * 20.; //p = abs(p) - 7.; p.xz *= rot(sin(time + op.x/100. + op.z/100.)); p.yz *= rot(p.x/40.); p.xz = smoothpmod(p.xz,2.,5.); p.yz = smoothpmod(p.yz,6.,2.); //p.xz = pmod(p.xz,8.); //p.yz = pmod(p.yz,8.); p.xz *= rot(time/32.); //p.xy = smoothpmod(p.xy,3.,.1); p.xz = sin(p.xz) * 5.; //p.x += pi/3. * step(abs(sin(p.x/2.)),0.); p -= clamp(p,-4.,4.); o = length(p) - 1.3; gp = p; // mp.y += 20.; // float sp = length(mp) - 50.; // if(o > sp) // { // o = sp; // id = 1.; // } return vec2(o,id); } vec2 skymarch(vec3 cp , vec3 rd) { float depth = 0.; for(int i = 0 ; i< 128 ; i++) { vec3 rp = cp + rd * depth; vec2 d = skymap(rp); if(abs(d.x) < 0.01) { return vec2(depth,d.y); } if(depth > 200.)break; depth += d.x; } return vec2(-depth , 0.); } vec2 shadowmarch(vec3 cp , vec3 rd) { float depth = 0.; for(int i = 0 ; i< 44 ; i++) { vec3 rp = cp + rd * depth; vec2 d = map(rp); if(abs(d.x) < 0.01) { return vec2(depth,d.y); } if(depth > 300.)break; depth += d.x; } return vec2(-depth , 0.); } //https://www.shadertoy.com/view/tlXBRl float calculateObstruction(vec3 pos, vec3 lpos, float lrad) { // A homemade algorithm to compute obstruction // Raymarch to the light source, and // record the largest obstruction. // We assume that if the ray passes through an object at depth // d (negative distance), then the object obstructs light // proportional to the relative size of d projected on the light // as given by Thales's theorem. vec3 toLight = normalize(lpos-pos); float distToLight = length(lpos-pos); float d, t=lrad*.1; float obstruction=0.; for(int j=0; j<128; j++) { d = map(pos + t*toLight).x; obstruction = max(.5+(-d)*distToLight/(2.*lrad*t), obstruction); if(obstruction >= 1.){break;} // If we're stuck, advance by the characteristic // size of an obstructing object t += max(d, lrad*t/distToLight); if(t >= distToLight) break; } return clamp(obstruction, 0.,1.); } vec3 getStar(vec3 rd) { rd.xz *= rot(time/10.); return vec3(0.,0.5,1.) + sin(rand(rd )*10.); } vec3 getColor(vec3 cp,vec3 rd) { vec3 col = vec3(0.); vec2 d = skymarch(cp,rd); if(d.x > 0.) { vec2 e = vec2(1.0, -1.0) * 0.005; vec3 pos = cp + rd * d.x; vec3 N = normalize( e.xyy * map(pos + e.xyy).x + e.yyx * map(pos + e.yyx).x + e.yxy * map(pos + e.yxy).x + e.xxx * map(pos + e.xxx).x); float rim = pow(clamp(1. - dot(N, -rd), 0., 1.), 3.)/6.; vec3 mat = vec3(1.); pp.y = abs(pp.y)-6.; vec3 pointlight = normalize(vec3(.1,0.1,0.)); pointlight.xz *= rot(time); vec3 ldir = normalize(pos - pointlight); float lrad = 0.2; float distToLight = length(pointlight + lrad - pos); float lightStrength = 3.1; float diff = max(0.,dot(ldir,N))/(distToLight*distToLight)*lightStrength; float sp = max(0.,dot(reflect(ldir,N) , -rd)); //diff = mix(diff, 1. , 0.); sp = pow(sp , 31.); diff = pow(diff , 3.); col = log2(cos(length(gp))*vec3(2.,3.,5.) ); col = mix(vec3(1.),col ,d.y ); float obstruction = calculateObstruction(pos,pointlight,lrad); float level = diff * (1. - obstruction); col = level * col + sp * col; col = 1.-exp(-2.*col); col += rim * .3; } return col; } void main(void) { vec2 p = (gl_FragCoord.xy * 2.0 - resolution.xy) / min(resolution.x, resolution.y); //p.y *= mix(1.,p.y,clamp( (time - timeline[7] + 4.),0.,1.) ); vec3 f = vec3(time/pi/2.,0.,0.) * 0.; vec3 cp = vec3(30.1,-193.,-270.) + f; cp.xz *= rot(time/14.); cp.yz *= rot(time/10.); vec3 target = vec3(0.,-40.,0.) + f; vec3 col = vec3(0.); vec3 cd = normalize(vec3(target - cp)); vec3 cs = normalize(cross(cd , vec3(0.,1.,0.))); vec3 cu = normalize(cross(cd,cs)); float fov =3.5-dot(p,p); vec3 rd = normalize(cd * fov + cs * p.x + cu * p.y); //rd = deMengerSponge3(rd,vec3(3.),3.); vec2 d = march(cp,rd); vec3 rd2 = rd; rd2.xz *= rot(time/10.); vec3 star = vec3(0.,0.5,1.) + sin(rand(rd2 )*10.); if( d.x > 0.) { vec2 e = vec2(1.0, -1.0) * 0.005; vec3 pos = cp + rd * d.x; vec3 N = normalize( e.xyy * map(pos + e.xyy).x + e.yyx * map(pos + e.yyx).x + e.yxy * map(pos + e.yxy).x + e.xxx * map(pos + e.xxx).x); vec3 sun = normalize(vec3(0.2,1.,.5)); float rim = pow(clamp(1. - dot(N, -rd), 0., 1.), 3.)/6.; vec3 mat = vec3(1.); pp.y = abs(pp.y)-6.; vec3 pointlight = normalize(vec3(0.1,50.1,0.)); //pointlight.xz *= rot(time); vec3 ldir = normalize(pos - pointlight); float lrad = 0.9; float distToLight = length(pointlight + lrad - pos); float lightStrength = 3000.; float diff = max(0.,dot(ldir,N))/(distToLight*distToLight)*lightStrength; float sp = max(0.,dot(reflect(ldir,N) , -rd)); //diff = mix(diff, 1. , 0.); sp = pow(sp , 31.); diff = pow(diff , 3.); col = log2(cos(length(gp))*vec3(2.,3.,5.) ); col = mix(vec3(1.),col ,d.y ); float obstruction = calculateObstruction(pos,pointlight,lrad); float level = diff * (1. - obstruction); col = level * col + sp * col; col = 1.-exp(-2.*col); col += rim * 2.; if(d.y > 0.) { col = getColor(pos - N*0.1 , refract(rd,N, 0.018)); col = col + sp * vec3(1.); col = 1.-exp(-2.*col); col += rim * 2.; } } col *= vec3(1.,1.4,1.2); //col *= mix(1.,0.,clamp( (time - length(p) + sin(atan(p.x,p.y) * pi ) * 0.3 ),0.,1.) ); gl_FragColor = vec4(col, 1.0); }