#extension GL_OES_standard_derivatives : enable
precision highp float;
uniform float time;
uniform vec2 mouse;
uniform vec2 resolution;
#define iTime mod(time,360.)
//東京
mat2 rot(float a){float c=cos(a),s=sin(a);return mat2(c,-s,s,c);}
float box(vec2 p,vec2 b){
vec2 q = abs(p)-b;
return length(max(vec2(0.),q))+min(0.,max(q.x,q.y));
}
float box3(vec3 p,vec3 b){
vec3 q =abs(p)-b;
return length(max(vec3(0.),q))+min(0.,max(q.x,max(q.y,q.z)));
}
float tokyo(vec3 p){
p*=2.;
vec3 op = p;
p.x -=6.;
float d = box3(p-vec3(0.,4.,0.),vec3(5.,1.,1.));
float q = min(d,box3(p-vec3(0.,0.,0.),vec3(4.,3.,1.)));
q = max(q,-min(d,box3(p-vec3(0.,0.,0.),vec3(3.,2.,2.))));
d = min(d,q);
d =min(d,box3(p,vec3(4.,1.,1.)));
d = min(d,box3(p,vec3(1.,7.,1.)));
p.y-=-5.;
p.x = abs(p.x)-2.;
p.xy*=rot(-.885);
d = min(d,box3(p,vec3(1.,2.5,1.)));
p =op;
p.x+=6.;
float d2 = box3(p-vec3(0.,4.,0.),vec3(5.,1.,1.));
float q2 = min(d2,box3(p-vec3(0.,0.,0.),vec3(4.,3.,1.)));
q2 = max(q2,-min(d2,box3(p-vec3(0.,0.,0.),vec3(3.,2.,2.))));
d2 = min(d2,q2);
d2 = min(d2,box3(p-vec3(0.,5.,0.),vec3(1.,2.,1.)));
d2 = min(d2,box3(p-vec3(0.,-5.,0.),vec3(1.,2.,1.)));
d2 = min(d2,box3(p-vec3(1,-6.5,0.),vec3(2.,.5,1.)));
p.y-=-4.;
p.x = abs(p.x)-2.;
p.xy*=rot(-.885);
d2 = min(d2,box3(p,vec3(.5,2.5,1.)));
d= min(d2,d);
return d/2.;
}
float diam(vec3 p,float s){
p = abs(p);
return (p.x+p.y+p.z-s)*inversesqrt(3.);
}
vec3 op;
float di(vec3 p){
float q = 0.;
vec4 pp = vec4(p,1.1);
const float cc = 2.;
for(float i=1.; i<=cc ;i++){
q += (asin(sin(pp.x))+acos(cos(pp.z)))/pp.a;
pp*=1.1;
pp.xyz +=vec3(1.1,2.,3.);
pp.xyz = abs(pp.xyz);
pp.xz *=rot(-.785*i+i);
}
return q/(cc+pp.a);
}
vec2 sdf(vec3 p){
vec3 tp = p;
op = p;
vec2 h;
p.xy *=rot(p.z*.1);
p.y= -abs(p.y)+(3.+asin(sin(iTime+p.z*.1))*.5);
h.x = length(p+1.)-.1;
p.y +=2.+di(p);
h.x = min(h.x, box(p.xy,vec2(6.1,.5) ) );
h.x *=.65;
h.y = 1.;
vec2 t;
tp.xz*=rot(pow(fract(iTime*.3),.125)*3.14);
tp.xy*=rot(pow(fract(iTime*.3+.5),.125)*3.14);
t.x = mix(diam(tp,2.5),box3(tp,vec3(1.1)),asin(sin(floor(iTime)+pow(fract(iTime),.5)))*.3+.2);
t.x =mix(t.x , tokyo(-tp),clamp(sin(time),-.5,.5)+.5);
tp = op;
vec4 ttp = vec4(tp,1.);
ttp.yz*=rot(iTime*.01);
for(float i=0.;i<=6.;i++){
ttp*=1.3;
ttp.xyz = abs(ttp.xzy)-11.;
ttp.xz*=rot(-.785+i);
ttp.xy*=rot(i);
t.x = min(t.x,box3(ttp.xyz,vec3(.5))/ttp.a);
}
t.y = 2.;
h = t.x < h.x ? t:h;
return h;
}
#define q(s) s*sdf(p+s).x
#define ao(rp,n,k) clamp(sdf(rp+n*k).x/k,.0,1.)
#define gao(rp,n) (ao(rp,n,.1)+ao(rp,n,.5)+ao(rp,n,.7))
vec3 norm(vec3 p,float e){vec2 nv=vec2(-e,e);return normalize(q(nv.xyy)+q(nv.yxy)+q(nv.yyx)+q(nv.xxx));}
vec3 pal(float t){return .5+vec3(.5,.2,.1)*cos(6.28*(vec3(1.,.5,.6)*t+vec3(.0,.1,.2)));}
void main( void ) {
vec2 position = ( gl_FragCoord.xy / resolution.xy ) + mouse / 4.0;
vec2 uv = (2.*gl_FragCoord.xy-resolution.xy)/resolution.y;
vec3 ro = vec3(12.,0.+mouse.y,-6.+mouse.x);
vec3 rt = vec3(0.);
vec3 z = normalize(rt-ro);
vec3 x = (cross(z,vec3(0.,-1.,.0)));
vec3 y = (cross(z,x));
vec3 rp = ro;
vec3 rd = normalize(mat3(x,y,z)*vec3(uv,.9));
vec3 light = mat3(x,y,z)*vec3(-2.,-1.,-4.);
vec3 col = vec3(.1);
for(float i=0.;i<=128.;i++){
vec2 d = sdf(rp);
if(length(rp) > 50.) break;
if(d.x <.001){
vec3 n = norm(rp,.019);
vec3 nn = norm(rp,.02+abs(sin(iTime+rp.z*.2))*.01);
float dif = dot(normalize(light-rp),n);
float spc = max(0.,dot(normalize(rp-ro),reflect(normalize(-light),n)));
spc = pow(spc,32.);
float fr = pow(1.-dot(-rd,n),4.);
col = vec3(.1)+vec3(.15,.5,.7)*dif+spc*.75;
float l = length(n-nn);
float s = step(.01,l);
col = max(vec3(0.),mix(col-s,col*s, floor(rp.z*.1-.3+di(rp.zzz*.1)*.2)*(gao(rp,n)/3.)*.75+.5));
col = pow(col,vec3(.8));
if(d.y ==1.) {col *=smoothstep(10.5,.0,op.z);col += -fr+s*pal(op.z*.1+iTime)*smoothstep(-.1,.1,op.z);}
break;
}
rp +=rd*d.x;
}
gl_FragColor = vec4(col,1.0 );
}
