// Copyright © 2021 IWBTShyGuy // Attribution 4.0 International (CC BY 4.0) precision highp float; uniform vec2 resolution; uniform float time; uniform vec2 mouse; uniform sampler2D backbuffer; #define iTime time #define iResolution resolution const float PI = 3.141592653; struct Camera { vec3 position; vec3 direction; vec3 up_direction; // not require dot(direction, up_direction) == 0 float fov; float aspect; // x / y }; struct Ray { vec3 origin; vec3 direction; }; Camera newCamera(vec3 position, vec3 direction, vec3 up_direction, float fov, float aspect) { Camera camera; camera.position = position; camera.direction = direction; camera.up_direction = up_direction; camera.fov = fov; camera.aspect = aspect; return camera; } // perspective camera ray, uv = fragCoord / iResolution.xy // cf: https://qiita.com/aa_debdeb/items/301dfc54788f1219b554 Ray cameraRay(in Camera camera, in vec2 uv) { uv = uv * 2.0 - 1.0; float radian = camera.fov; float h = tan(radian * 0.5); float w = h * camera.aspect; vec3 right = normalize(cross(camera.direction, camera.up_direction)); vec3 up = normalize(cross(right, camera.direction)); vec3 direction = normalize(right * w * uv.x + up * h * uv.y + camera.direction); Ray ray; ray.origin = camera.position; ray.direction = direction; return ray; } float microfacet(vec3 normal, vec3 cameraDir, vec3 lightDir, float metal, float roughness) { vec3 middle = normalize(cameraDir + lightDir); float dotCN = clamp(-dot(cameraDir, normal), 0.0, 1.0); float dotLN = clamp(-dot(lightDir, normal), 0.0, 1.0); float dotNM = -dot(normal, middle); float dotCM = dot(cameraDir, middle); float alpha = roughness * roughness; // diffuse BRDF float diffuse = 1.0 - metal; // microfacet distribution float alpha2 = alpha * alpha; float tmp = 1.0 - dotNM * dotNM * (1.0 - alpha2); float distribution = alpha2 / (tmp * tmp); // schlick approxy & geometric decay float alp = alpha * 0.5; float sal = dotLN / (dotLN * (1.0 - alp) + alp); float sac = dotCN / (dotCN * (1.0 - alp) + alp); float decay = sal * sac; // fresnel float c = 1.0 - dotCM; c = c * c * c * c * c; float fresnel = metal + (1.0 - metal) * c; // specular BRDF tmp = 4.0 * dotCN * dotLN; float specular = distribution * decay / tmp * fresnel; if (tmp < 0.0001) specular = 0.0; return (diffuse + specular) * dotLN; } // Hash without Sine https://www.shadertoy.com/view/4djSRW float hash13(vec3 p3) { p3 = fract(p3 * .1031); p3 += dot(p3, p3.yzx + 33.33); return fract((p3.x + p3.y) * p3.z); } mat3 transpose(mat3 a) { float b = a[0][1]; a[0][1] = a[1][0]; a[1][0] = b; b = a[0][2]; a[0][2] = a[2][0]; a[2][0] = b; b = a[1][2]; a[1][2] = a[2][1]; a[2][1] = b; return a; } // ----------------------- Config ----------------------- // // Increasing these values will increase the visual impact, // and the load on the GPU will increase dramatically. // Farthest distance to draw const float FAR = 4.0; // the number of coins in one line const int LINE_COIN = 3; // --------------------- Basic SDF --------------------- // // https://www.iquilezles.org/www/articles/distfunctions/distfunctions.htm float sdCylinder(in vec3 p, in float r, in float h) { vec2 d = abs(vec2(length(p.xz),p.y)) - vec2(r,h); return min(max(d.x,d.y),0.0) + length(max(d,0.0)); } // https://www.iquilezles.org/www/articles/distfunctions2d/distfunctions2d.htm float sdStar5(in vec2 p, in float r, in float rf) { const vec2 k1 = vec2(0.809016994375, -0.587785252292); const vec2 k2 = vec2(-k1.x,k1.y); p.x = abs(p.x); p -= 2.0*max(dot(k1,p),0.0)*k1; p -= 2.0*max(dot(k2,p),0.0)*k2; p.x = abs(p.x); p.y -= r; vec2 ba = rf*vec2(-k1.y,k1.x) - vec2(0,1); float h = clamp( dot(p,ba)/dot(ba,ba), 0.0, r ); return length(p-ba*h) * sign(p.y*ba.x-p.x*ba.y); } const vec3 SUN_SHINE = normalize(vec3(1, -1, -1)); // ----------------------- PIPE ----------------------- // const float PIPE_SIZE = 0.075; // sdf float sdPipe(in vec3 p) { float dist0 = sdCylinder(p, PIPE_SIZE * 0.5, PIPE_SIZE * 0.51); float dist1 = sdCylinder( p - vec3(0, 0.3 * PIPE_SIZE, 0), 0.55 * PIPE_SIZE, 0.15 * PIPE_SIZE ) - 0.06 * PIPE_SIZE; return min(dist0, dist1); } // normal vec3 normalPipe(in vec3 p) { const vec2 h = vec2(0.0001,0); return normalize(vec3( sdPipe(p+h.xyy) - sdPipe(p-h.xyy), sdPipe(p+h.yxy) - sdPipe(p-h.yxy), sdPipe(p+h.yyx) - sdPipe(p-h.yyx) )); } // albedo vec3 albPipe(in vec3 p) { float c = smoothstep( 0.0, 1.0, (length(p.zx) - 0.42 * PIPE_SIZE) / (0.01 * PIPE_SIZE) ); return vec3(0, c, 0); } // radiance without distance decay float radPipe(in vec3 normal, in vec3 dir) { float c0 = microfacet( normal, dir, dir, 0.0, 0.5 ); float c1 = microfacet( normal, dir, SUN_SHINE, 0.9, 0.5 ); return c0 * 0.2 + c1 * 0.8; } // ----------------------- COIN ----------------------- // const float COIN_SIZE = 0.03; // sdf float sdCoin(in vec3 p) { return sdCylinder( p, 0.98 * COIN_SIZE, 0.08 * COIN_SIZE ) - 0.02 * COIN_SIZE; } const vec3 COIN_ALBEDO = vec3(241, 208, 12) / 255.0; // bumpmap vec3 normalCoin(in vec3 p) { float r = 0.98 * COIN_SIZE; float h = 0.08 * COIN_SIZE; vec2 e = vec2(0.0001,0); vec3 base = normalize(vec3( sdCoin(p + e.xyy) - sdCoin(p - e.xyy), sdCoin(p + e.yxy) - sdCoin(p - e.yxy), sdCoin(p + e.yyx) - sdCoin(p - e.yyx) )); if (abs(p.y) < h - 0.001) { return base; } else { vec2 uv = p.zx / r; float r0 = length(uv) - 0.7; float c0 = min( max(0.0, r0 - 0.05), max(0.0, 0.15 - r0) ) * 20.0; vec2 tmp = normalize(p.zx); vec3 n0 = vec3(tmp.y, 0, tmp.x); vec2 e = vec2(0.00001, 0); float r1 = sdStar5(uv, 0.6, 0.4); float c1 = -min( max(0.0, r1 + 0.025), max(0.0, 0.025 - r1) ) * 20.0; vec3 n1 = normalize(vec3( sdStar5(uv+e.yx, 0.6, 0.4) - sdStar5(uv-e.yx, 0.6, 0.4), 0, sdStar5(uv+e.xy, 0.6, 0.4) - sdStar5(uv-e.xy, 0.6, 0.4) )); return normalize(base - c0 * n0 - c1 * n1); } } // radiance without distance decay float radCoin(in vec3 normal, in vec3 dir) { float c0 = microfacet( normal, dir, dir, 0.0, 0.5 ); float c1 = microfacet( normal, dir, SUN_SHINE, 0.9, 0.5 ); return c0 * 0.2 + c1 * 0.8; } // ----------------------- SDF ----------------------- // #define rotate3D(n, t) mat3( \ cos(t) + n.x * n.x * (1.0 - cos(t)), n.x * n.y * (1.0 - cos(t)) + n.z * sin(t), n.z * n.x * (1.0 - cos(t)) - n.y * sin(t), \ n.x * n.y * (1.0 - cos(t)) - n.z * sin(t), cos(t) + n.y * n.y * (1.0 - cos(t)), n.y * n.z * (1.0 - cos(t)) + n.x * sin(t), \ n.z * n.x * (1.0 - cos(t)) + n.y * sin(t), n.y * n.z * (1.0 - cos(t)) - n.x * sin(t), cos(t) + n.z * n.z * (1.0 - cos(t)) \ ) const mat3 T0 = rotate3D(normalize(vec3(0, 1, -1)), PI); const mat3 T1 = rotate3D(normalize(vec3(1, 0, 1)), PI); const mat3 T2 = rotate3D(normalize(vec3(0, 1, 1)), PI); // hexahedron group mat3 S(int index) { if (index == 0) return mat3(1); else if (index == 1) return T2; else if(index == 2) return T1; else if (index == 3) return T2 * T1; else if (index == 4) return T1 * T2; else if (index == 5) return T2 * T1 * T2; else if (index == 6) return T0; else if (index == 7) return T2 * T0; else if (index == 8) return T1 * T0; else if (index == 9) return T2 * T1 * T0; else if (index == 10) return T1 * T2 * T0; else if (index == 11) return T2 * T1 * T2 * T0; else if (index == 12) return T0 * T1; else if (index == 13) return T2 * T0 * T1; else if (index == 14) return T1 * T0 * T1; else if (index == 15) return T2 * T1 * T0 * T1; else if (index == 16) return T1 * T2 * T0 * T1; else if (index == 17) return T2 * T1 * T2 * T0 * T1; else if (index == 18) return T0 * T1 * T2; else if (index == 19) return T2 * T0 * T1 * T2; else if (index == 20) return T1 * T0 * T1 * T2; else if (index == 21) return T2 * T1 * T0 * T1 * T2; else if (index == 22) return T1 * T2 * T0 * T1 * T2; else return T2 * T1 * T2 * T0 * T1 * T2; } struct MarchingOutput { float dist; int shape; mat3 mat; vec3 p; }; // cell-wise sdf MarchingOutput sDist(in vec3 p, in float t) { int index = int(hash13(floor(p) + 0.23) * 24.0); p = fract(p); MarchingOutput op, op0; // ----------------- PIPE ----------------- // op0.shape = 0; op0.mat = mat3(1); op0.p = p - vec3(0.5, 0.5 * PIPE_SIZE, 0.5); op0.dist = sdPipe(op0.p); op = op0; op0.mat = rotate3D(vec3(1, 0, 0), PI); op0.p = op0.mat * (p - vec3(0.5, 1.0 - 0.5 * PIPE_SIZE, 0.5)); op0.dist = sdPipe(op0.p); if (op0.dist < op.dist) op = op0; op0.mat = rotate3D(vec3(0, 0, 1), PI/2.0); op0.p = op0.mat * (p - vec3(0.5 * PIPE_SIZE, 0.5, 0.5)); op0.dist = sdPipe(op0.p); if (op0.dist < op.dist) op = op0; op0.mat = rotate3D(vec3(0, 0, 1), -PI/2.0); op0.p = op0.mat * (p - vec3(1.0 - 0.5 * PIPE_SIZE, 0.5, 0.5)); op0.dist = sdPipe(op0.p); if (op0.dist < op.dist) op = op0; op0.mat = rotate3D(vec3(1, 0, 0), -PI/2.0); op0.p = op0.mat * (p - vec3(0.5, 0.5, 0.5 * PIPE_SIZE)); op0.dist = sdPipe(op0.p); if (op0.dist < op.dist) op = op0; op0.mat = rotate3D(vec3(1, 0, 0), PI/2.0); op0.p = op0.mat * (p - vec3(0.5, 0.5, 1.0 - 0.5 * PIPE_SIZE)); op0.dist = sdPipe(op0.p); if (op0.dist < op.dist) op = op0; // ----------------- COIN ----------------- // op0.shape = 1; op0.mat = rotate3D(vec3(0, 0, 1), -PI / 2.0) * rotate3D(vec3(0, 1, 0), -t); mat3 rot = S(index); for (int i = 0; i < LINE_COIN; i++) { float k = float(i) / float(LINE_COIN); vec3 disp = vec3(0.5, fract(t * 0.1 + k), 0.5); disp = fract(rot * (disp - 0.5) + 0.5); op0.p = op0.mat * (p - disp); op0.dist = sdCoin(op0.p); if (op0.dist < op.dist) op = op0; } for (int i = 0; i < LINE_COIN; i++) { float k = float(i) / float(LINE_COIN); float theta = (fract(t * 0.1 + k) - 0.05) * PI * 0.55; vec2 disp2d = vec2(cos(theta), sin(theta)) * (0.5 - PIPE_SIZE) + PIPE_SIZE; vec3 disp = vec3(disp2d.y, 0.5, disp2d.x); disp = fract(rot * (disp - 0.5) + 0.5); op0.p = op0.mat * (p - disp); op0.dist = sdCoin(op0.p); if (op0.dist < op.dist) op = op0; } for (int i = 0; i < LINE_COIN; i++) { float k = float(i) / float(LINE_COIN); float theta = (fract(t * 0.1 + k) - 0.05) * PI * 0.55; vec2 disp2d = vec2(cos(theta), sin(theta)) * (0.5 - PIPE_SIZE) + PIPE_SIZE; disp2d = 1.0 - disp2d; vec3 disp = vec3(disp2d.y, 0.5, disp2d.x); disp = fract(rot * (disp - 0.5) + 0.5); op0.p = op0.mat * (p - disp); op0.dist = sdCoin(op0.p); if (op0.dist < op.dist) op = op0; } return op; } void mainImage(out vec4 fragColor, in vec2 fragCoord) { vec2 theta = vec2(0.15, 0.25) * iTime; theta.y = sin(theta.y) * 0.5; Camera camera = newCamera( vec3(sin(iTime * 0.1) * 0.1, -iTime * 0.1, cos(iTime * 0.1) * 0.1), vec3(cos(theta.x) * cos(theta.y), sin(theta.y), sin(theta.x) * cos(theta.y)), //vec3(0), //vec3(0, 0, 1), vec3(0, 1, 0), PI / 4.0, iResolution.x / iResolution.y ); Ray ray = cameraRay(camera, fragCoord / iResolution.xy); vec3 p = ray.origin; float dist = 0.0; MarchingOutput op; for (int _i = 0; _i < 100; _i++) { op = sDist(p, iTime); if (op.dist < 0.00001 || dist >= FAR) break; dist += op.dist; p = ray.origin + dist * ray.direction; } const vec3 BLUE0 = vec3(160, 216, 234) / 255.0; const vec3 BLUE1 = vec3(108, 155, 210) / 255.0; float y = (1.0 + ray.direction.y) * 0.5; vec3 col = pow((1.0 - y) * BLUE0 + y * BLUE1, vec3(2.2)); if (dist < FAR) { float decay = smoothstep(0.0, 1.0, dist / FAR); decay *= decay; if (op.shape == 0) { // pipe vec3 normal = transpose(op.mat) * normalPipe(op.p); float rad = radPipe(normal, ray.direction); vec3 albedo = albPipe(op.p); col = (1.0 - decay) * rad * albedo + decay * col; } else { // coin vec3 normal = transpose(op.mat) * normalCoin(op.p); float rad = radCoin(normal, ray.direction); col = (1.0 - decay) * rad * COIN_ALBEDO + decay * col; } } fragColor = vec4(pow(col, vec3(0.4545)), 1); } void main(void) { mainImage(gl_FragColor, gl_FragCoord.xy); }

3D truchet - Bonus Stage