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 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); 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);} vec3 random33(vec3 p) { p=vec3(dot(p,vec3(127.1,311.7,74.7)),dot(p,vec3(269.5,183.3,246.1)),dot(p,vec3(113.5,271.9,124.6))); return fract(sin(p)*43758.5453123); } vec2 min2(vec2 a,vec2 b){ if(a.x < b.x)return a; return b; } vec3 glp; vec3 gap; vec3 gp ; #define TT(x,y) floor(x)+pow(fract(x),y) #define HEAD 1. #define LEG 2. #define ARM 3. vec2 map(vec3 p) { float t = time * 5.; p.y+=1. + sin(t*2.)/6.; p.xz*= rot(time); gp = p; vec2 o=vec2(10.,0.); o.x=length(p)-1.; o=min2(o,vec2(length(p-vec3(0.,2.,0.))-1.2,HEAD)); vec3 dp=p; dp.y-=.2; dp.z+=.4; dp.yz*=rot(.3); vec3 ds=vec3(1.,.01,1.)/2.; // o=min2(o,vec2(length(dp-clamp(dp,-ds,ds))-.1,4.) ); vec3 ep=p; vec3 es=vec3(1.,1.,.01)/7.; ep.x=abs(ep.x)-.4; ep.z+=1.; ep.y+=.3; // o=min2(o,vec2(length(ep-clamp(ep,-es,es))-.01,4.)); vec3 ap=p; ap.y-=1.5; ap.yz *= rot(sin(t)/5. * sign(-ap.x)); gap = ap; ap.x=abs(ap.x)-1.; ap.xy*=rot(.3); ap.xz*=rot(-.6); o=min2(o,vec2(length(ap-clamp(ap,-vec3(.5,.1,.1),vec3(1.,.1,.1)))-.1,ARM)); vec3 lp=p; lp.y-=3.2; lp.yz *= rot(sin(t)/5. * sign(lp.x)); glp = lp; lp.x=abs(lp.x)-1.; lp.xz*=rot(-1.); o=min2(o,vec2(length(lp-clamp(lp,-vec3(1.,.03,.6),vec3(1.,.03,.1)))-.1,LEG)); return o; } vec2 march(vec3 cp,vec3 rd) { float depth=0.; for(int i=0;i<99;i++){ vec3 rp=cp+rd*depth; vec2 d=map(rp); if(abs(d.x)<.01)return vec2(depth,d.y); depth+=d.x; } return vec2(-depth,0.); } vec3 getNormal(vec3 pos) { float eps=.05; vec3 n=vec3( map(pos+vec3(eps,0.,0.)).x-map(pos-vec3(eps,0.,0.)).x, map(pos+vec3(0.,eps,0.)).x-map(pos-vec3(0.,eps,0.)).x, map(pos+vec3(0.,0.,eps)).x-map(pos-vec3(0.,0.,eps)).x ); return normalize(n); } void cellular3D(vec3 p,float scale,out vec3 cellId,out float edgeDist,out vec3 edgeN) { vec3 q=p*scale; vec3 fp=floor(q); vec3 sp=fract(q); float md=100.; vec3 mr=vec3(0.); vec3 mg=vec3(0.); for(int z=-1;z<=1;z++){ for(int y=-1;y<=1;y++){ for(int x=-1;x<=1;x++){ vec3 g=vec3(float(x),float(y),float(z)); vec3 r=g+random33(fp+g)-sp; float d=dot(r,r); if(d.0001){ vec3 n=normalize(dv); float d=dot((mr+r)*.5,n); if(d.25)col=c1; if(h.x>.50)col=c2; if(h.x>.75)col=c3; vec3 a=normalize(h-.5+vec3(.001)); vec3 axis=abs(a.y)<.8?vec3(0.,1.,0.):vec3(1.,0.,0.); vec3 u=normalize(cross(a,axis)); vec3 v=normalize(cross(a,u)); float weave=sin(dot(pos,u)*130.+h.y*20.)*sin(dot(pos,v)*130.+h.z*20.); col*=.92+weave*.08; return col; } //ちゃんと理解できてないゾーン // 縫い目 float stitch3D(vec3 pos,vec3 N,vec3 cellId,float edgeDist,vec3 edgeN,float aa) { vec3 tangent=cross(N,edgeN); float tl=length(tangent); if(tl<.001)return 0.; tangent/=tl; if(dot(tangent,vec3(.71,.43,.56))<0.)tangent=-tangent; float inset=.035; float spacing=.08; float stitchLength=.02; float stitchWidth=.004; float phase=random33(cellId).x*spacing; vec3 axis=abs(edgeN.y)<.9?vec3(0.,1.,0.):vec3(1.,0.,0.); vec3 alongAxis=normalize(cross(edgeN,axis)); if(dot(alongAxis,vec3(.71,.43,.56))<0.)alongAxis=-alongAxis; float along=dot(pos,alongAxis)+phase; float stitchPos=mod(along+spacing*.5,spacing)-spacing*.5; float across=abs(edgeDist-inset); vec2 stitchP=vec2(across-stitchWidth,abs(stitchPos)-stitchLength*.5); float stitchDist=max(stitchP.x,stitchP.y); return 1.-smoothstep(0.,aa,stitchDist); } //ここまで void main(void) { vec2 p=(gl_FragCoord.xy*2.-resolution.xy)/min(resolution.x,resolution.y); vec3 cp=vec3(0.,0.,-8.); vec3 target=vec3(0.); vec3 col=vec3(0.); vec3 cd=normalize(target-cp); vec3 cs=normalize(cross(cd,vec3(0.,1.,0.))); vec3 cu=normalize(cross(cd,cs)); float fov=2.5; vec3 rd=normalize(cd*fov+cs*p.x+cu*p.y); vec2 d=march(cp,rd); if(d.x>0.){ vec3 pos=cp+rd*d.x; vec3 pos1=gp; if(d.y==ARM){ pos1=gap; }else if(d.y==LEG){ pos1=glp; } vec3 worldN=getNormal(pos); vec3 stitchN=worldN; float t=time*5.; stitchN.xz*=rot(time); if(d.y==ARM){ stitchN.yz*=rot(sin(t)/5.*sign(-pos1.x)); }else if(d.y==LEG){ stitchN.yz*=rot(sin(t)/5.*sign(pos1.x)); } vec3 axis=abs(worldN.y)<.999?vec3(0.,1.,0.):vec3(1.,0.,0.); vec3 side=normalize(cross(axis,worldN)); vec3 side2=normalize(cross(worldN,side)); float fine=150.; float y1=noise(pos1*fine)-.5; float y2=noise(pos1.yxz*fine)-.5; vec3 fuwafuwaN=normalize(worldN+y1*side/1.3+y2*side2/1.3); float cellScale=1.2; vec3 cellId; vec3 edgeN; float edgeDist; cellular3D(pos1,cellScale,cellId,edgeDist,edgeN); vec3 base=patchColor(cellId,pos1); base=mix(base,vec3(.2),step(3.6,d.y)); vec3 sun=normalize(vec3(2.,4.,8.)); float diff=max(0.,dot(-sun,fuwafuwaN)); diff=mix(diff,1.,.1); float sp=max(0.,dot(rd,reflect(fuwafuwaN,sun))); sp=pow(sp,56.6); col=base*diff+sp*vec3(.4); float aa=2./min(resolution.x,resolution.y); float seam=1.-smoothstep(.01,.03+aa*3.,edgeDist); col*=1.-seam*.45; float stitch=stitch3D(pos1,stitchN,cellId,edgeDist,edgeN,aa); stitch=mix(stitch,0.,step(3.5,d.y)); vec3 threadColor=vec3(.5); col=mix(col,threadColor,stitch); col+=stitch*diff*.08; } gl_FragColor=vec4(col,1.); }

#stitch