precision highp float; uniform vec2 resolution; uniform float time; uniform vec2 mouse; uniform sampler2D backbuffer; #define FAR 195. #define M(x,y) mod(x,y) - y/2. #define RR(x) (x - .5) *2. #define ITR 128 float dep,tt; float pi =acos(-1.); float pi2 = acos(-1.) * 2.; float owner = 0.; vec3 spos; float rsp(vec3 p, vec3 s,float r){return ( length(p / s) - r) * min(s.x,min(s.y,s.z));} vec2 rot(vec2 p,float a){return p * mat2(cos(a),sin(a),-sin(a),cos(a));} float bo( vec3 p, vec3 s ){p = abs(p) - s;return length(max(p,0.0)) + min(max(p.x,max(p.y,p.z)),0.0);} //https://www.iquilezles.org/www/articles/distfunctions/distfunctions.htm float sdCappedCylinder( vec3 p, float h, float r ) { vec2 d = abs(vec2(length(p.xz),p.y)) - vec2(h,r); return min(max(d.x,d.y),0.0) + length(max(d,0.0)); } vec2 random2(vec2 st){ st = vec2( dot(st,vec2(127.1,311.7)), dot(st,vec2(269.5,183.3)) ); return -1.0 + 2.0*fract(sin(st)*43758.5453123); } float noise(vec2 st) { vec2 p = floor(st); vec2 f = fract(st); vec2 u = f*f*(3.0-2.0*f); vec2 v00 = random2(p+vec2(0,0)); vec2 v10 = random2(p+vec2(1,0)); vec2 v01 = random2(p+vec2(0,1)); vec2 v11 = random2(p+vec2(1,1)); return mix( mix( dot( v00, f - vec2(0,0) ), dot( v10, f - vec2(1,0) ), u.x ), mix( dot( v01, f - vec2(0,1) ), dot( v11, f - vec2(1,1) ), u.x ), u.y)+0.5; } float fbm (vec2 st) { float f = 0.; vec2 q = st/2.; for(int i = 0 ;i < 4;i++){ f += noise(q)/pow(2.,float(i)); q = q * (2.00+float(i)/100.); } return f; } 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; } //https://www.shadertoy.com/view/XsX3zB /* skew constants for 3d simplex functions */ const float F3 = 0.3333333; const float G3 = 0.1666667; /* 3d simplex noise */ float simplex3d(vec3 p) { /* 1. find current tetrahedron T and it's four vertices */ /* s, s+i1, s+i2, s+1.0 - absolute skewed (integer) coordinates of T vertices */ /* x, x1, x2, x3 - unskewed coordinates of p relative to each of T vertices*/ /* calculate s and x */ vec3 s = floor(p + dot(p, vec3(F3,F3,F3))); vec3 x = p - s + dot(s, vec3(G3,G3,G3)); /* calculate i1 and i2 */ vec3 e = step(vec3(0.,0.,0.), x - x.yzx); vec3 i1 = e*(1.0 - e.zxy); vec3 i2 = 1.0 - e.zxy*(1.0 - e); /* x1, x2, x3 */ vec3 x1 = x - i1 + G3; vec3 x2 = x - i2 + 2.0*G3; vec3 x3 = x - 1.0 + 3.0*G3; /* 2. find four surflets and store them in d */ vec4 w, d; /* calculate surflet weights */ w.x = dot(x, x); w.y = dot(x1, x1); w.z = dot(x2, x2); w.w = dot(x3, x3); /* w fades from 0.6 at the center of the surflet to 0.0 at the margin */ w = max(0.6 - w, 0.0); /* calculate surflet components */ 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); /* multiply d by w^4 */ w *= w; w *= w; d *= w; /* 3. return the sum of the four surflets */ return dot(d, vec4(52.0,52.0,52.0,52.0)); } float smin(float d1, float d2, float k){ float h = exp(-k * d1) + exp(-k * d2); return -log(h) / k; } vec3 rotmat(vec3 p,vec3 axis, float angle) { axis = normalize(axis); float s = sin(angle); float c = cos(angle); float oc = 1.0 - c; return p * 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 ); } float replim( float p, float s, float l ) { return p-s*clamp(floor(p/s + .5),-l,l); } float c(float x, float f) { return x - (x - x * x) * -f; } float humanF(vec3 p,float s) { vec3 fp = p; vec3 size = vec3(1.,7.5,1.7); fp.xz = rot(fp.xz,sign(s) * fp.z/9.); float b = bo(p + vec3(0.,0.,.3),vec3(1.,10.5,.05)); float h = min(rsp(fp + vec3(0.,0.,-0.4),size,.5) ,rsp(p + vec3(0.,0.,.6),vec3(.7,7.5,.9),.5)); h = max(h,-b); p.y -= abs( sin(clamp( p.z,-.2,1. ) * 16.)/30. ); b = bo(p + vec3(0.,-0.6,0.),vec3(1.8,0.5,1.8)); h = max(h ,-b); return h; } float humanLL(vec3 p) { p.x += 0.5; float ttl = tt; ttl += p.y/13.; p.z -= (floor(ttl / pi2) + pow( fract(ttl / pi2),7.) )* pi2; p.y += 2.5 + clamp(cos(tt + pi/7.),0.,1.); float s = rsp(p - vec3(0.2,1.9,1.),vec3(.6,.5,1.5),.8); p.z += 0.5; float l = sdCappedCylinder(p,.4,1.9); l = min(s,l); return l; } float humanRL(vec3 p) { p.x -= 0.5; float ttl = tt; ttl += p.y/13.; p.z -= (floor(ttl / pi2 + .5) + pow( fract(ttl / pi2 + .5),7.) )* pi2 - pi; p.y += 2.5 + clamp(sin(tt - pi/7.),0.,1.); float s = rsp(p - vec3(0.2,1.9,1.),vec3(.6,.5,1.5),.8); p.z += 0.5; float l = sdCappedCylinder(p,.4,1.9); l = min(s,l); return l; } float human(vec3 p) { vec3 pp = p ; float s = sign(sin(p.z + pi/2.)); pp.z = replim(pp.z,pi,floor(tt/pi)); pp.x += sin(p.z)/5. + s/2.; float f = humanF(pp,-s); return f; } float hlegs(vec3 p) { float leg = humanLL(p); leg = min(humanRL(p),leg); return leg; } float rf(vec3 p) { p.y += 1.3; vec3 size = vec3(.4,2.8,.5); p.xz = rot(p.xz,pi/6.); float s = 100.; for(int i = 0; i< 3 ; i++) { s = smin(s, rsp(p,size,.5),16.); s = smin(s,rsp(p - vec3(0.,.7,.4),vec3(1.8,3.8,7.9),.2 ),8.3); p.z += 0.5; p.xz = rot(p.xz,-pi/5.); p.z -= 0.5; } return s; } float rr(vec3 p) { vec3 pp = p ; float s = sign(sin(p.z + pi/2.)); pp.z = replim(pp.z,pi,floor(tt/pi)); pp.x += sin(p.z)/5. + s/.8; float f = rf(pp); return f; } float rlegs(vec3 p) { float leg = humanLL(p/vec3(1.9 - p.y/9.,1.5,1.) + vec3(.1,.0,-1.)); leg = min(humanRL(p/vec3(1.9 - p.y/9.,1.5,1.)+ vec3(-.38,.0,-1.)),leg); return leg; } float map(vec3 p) { float s = -p.y - ( fbm(p.xz/1.) - .8)/3.; float hs = 100.; float h = human(p); float w = max( max(s,h) ,-h); hs = max(s,-h); hs = min(hs,w*.55); float rs = 100.; float r = rr(p); rs = max(s,-r); s = mix(hs,rs,owner); return s * .5; } float m(vec3 cp,vec3 rd) { float dd = 0.; for(int i = 0; i< ITR;i++) { vec3 rp = rd * dd + cp; float d = map(rp); if(abs(d) < 0.00001 || dd > FAR){break;} dd += d; } dep = dd; if(dd > FAR){dd *= -1.;} return dd ; } float shmap(vec3 p) { p.y *= -1.; float s = 100.; float hs = 100.; vec3 pp = p; s = hlegs(pp); pp -= vec3(0.,-7.,tt); pp.xz = rot(pp.xz,pp.y + tt + sin(tt /2.413)); float b = bo(pp,vec3(2.5,2.,0.7)); hs = min(s,b); float rs = 100.; s = rlegs(p); p -= vec3(0.,-7.,tt); p.x += sin(p.z/7. + tt/3.21); float c = rsp(p + vec3(0.,-4.,19.9),vec3(1.5,.7,9.),1.8); // p.xz = rot(p.xz,p.y + tt + sin(tt /2.413)); c = min(c,rsp(p,vec3(1.9,2.5,1.7),1.8)); rs = min(s,c); s = mix(hs,rs,owner); return s * .7; } float shm(vec3 cp,vec3 rd) { float dd = 0.; for(int i = 0; i< ITR;i++) { vec3 rp = rd * dd + cp; float d = shmap(rp); if(abs(d) < 0.001 || dd > FAR){break;} dd += d; } if(dd > FAR){dd *= -1.;} return dd; } float moya(vec3 cp,vec3 rd) { float o = 0.; for(int i = 1; i < 8; i++) { vec3 p = cp + rd * float(i) * .9 ; if(map(p) < 0.01){break;} o += simplex3d(p + tt/11.675) / ( float(i)+1. ); } return o; } void main(void) { vec2 p = (gl_FragCoord.xy * 2.0 - resolution.xy) / min(resolution.x, resolution.y); vec3 cp = vec3(.5,-18.9,-7.) ; // cp.xz = rot(cp.xz,(mouse.x - .5) * 4.); // cp.yz = rot(cp.yz,(mouse.y) * 2.5); vec3 cd; vec3 target ; float scene = mod(time,59.) + 30. ; tt = scene * 4.34; vec3 fw = vec3(0.,0.,tt); if(scene < 36.) {//1 cp.xz = rot(cp.xz,4.5 ); cp.yz = rot(cp.yz,.8 ); cp.z += scene * 1.6; target = vec3(.1,-0.5,6.6) ; cd = normalize(target - cp); cp += cd * -10.5; }else if(scene < 48.) {//2 cp = vec3(-.2,.9,-.5) ; target = vec3(-0.13,.8,6.1);// + fw; //cp += fw; cd = normalize(target - cp); cp += cd * 160.7; } else if(scene < 68.) {//walk main cp.xz = rot(cp.xz,4.5); cp.yz = rot(cp.yz,.8); cp = rotmat(cp,vec3(0.,1.,0.),scene/9.9 - .4); target = vec3(.1,-.5,0.) + fw; cp += fw; cd = normalize(target - cp); scene/=1.3; owner = clamp( abs( floor(2.* sin(sin(scene * 1.35)+scene) /2.) ) ,0.,1.); } else if( scene < 73.) { cp.xz = rot(cp.xz,4.5); cp.yz = rot(cp.yz,.8); cp = rotmat(cp,vec3(0.,1.,0.),-12.); cp.z -= 0.; target = vec3(.1,-.5,1.) ; target.z += 68. * 4.34; cp.z += 68. * 4.34; cd = normalize(target - cp); scene/=1.3; owner = clamp( abs( floor(sin(scene/2.) + sin(scene/2.)) + sin(floor(scene/2.)))/2. ,0.,1.); } else if(scene < 89.) { cp.xz = rot(cp.xz,4.5); cp.yz = rot(cp.yz,.8); cp = rotmat(cp,vec3(0.,1.,0.),-12.); cp.z -= 68. ; target = vec3(.1,-.5,1.) ; target.z += 68. * 4.34; cp.z += 68. * 4.34; cd = normalize(target - cp); scene/=1.3; owner = clamp( abs( floor(sin(scene/2.) + sin(scene/2.)) + sin(floor(scene/2.)))/2. ,0.,1.); } vec3 cs = normalize(cross(vec3(0.,1.,0.),cd)); vec3 cu = normalize(cross(cs,cd)); float fov = 2.5; vec3 rd = normalize(p.x * cs + p.y * cu + cd * fov); float d = m(cp,rd); vec3 sun = normalize(vec3(0.2,.4,.8)); vec3 skycolor = vec3(.9,.9,.9) ; vec3 c = skycolor; if(d > 0.) { spos = rd * dep + cp; vec2 e = vec2(0.,0.01); vec3 N = normalize(map(spos) - vec3(map(spos + e.xyy),map(spos + e.yxy) , map(spos + e.yyx))); c = vec3(1.) ; float sray = shm( (spos + N * .001) ,sun); float shadow = step(sray,0.) * .4 + .6; float diff = .1 + .8 * max( dot(sun,N),0.); float skyc = max( dot(vec3(0.,1.,0.),N),0.); float bouc = max( dot(vec3(0.,-1.,0.),N),0.); // c = N; c *= diff; c *= shadow; } c = mix(c,skycolor,1. - exp(-0.000003 * dep * dep * dep )); c = mix(c,vec3(moya(cp,rd)*.2+.8),.3); c = clamp(c,0.,1.); gl_FragColor = vec4(c, 1.0); }