#version 300 es #define PI 3.14159265359 #define PI2 6.28318530718 precision highp float; uniform float time; uniform float iGlobalTime; in vec2 UV; out vec4 out_color; uniform vec3 iResolution; #define FACE_MATERIAL 0 #define HAT_MATERIAL 1 #define EYE_MATERIAL 2 // Mat4 data // distance, material, object index, extra info for material // x offset, y offset, z offset, available // available, available, available, available // available, available, available, available mat4 traceVec2ToMat4(vec2 positionAndMaterial) { return mat4(positionAndMaterial, 0.0, 0.0, vec4(0.0), vec4(0.0), vec4(0.0)); } vec2 traceMat4ToVec2(mat4 mat) { return mat[0].xy; } mat4 hat1(vec3 pos){ pos.xyz -= vec3(0.0,0.15,0.0); if (pos.y < -0.7) { return traceVec2ToMat4(vec2(1e10, float(HAT_MATERIAL))); } pos.y += 0.02 * cos(pos.x * 10.0); pos.y += 0.04 * cos(pos.x * 2.0); pos.xz *= 1.0 - 0.7 * (-pos.y - 0.0); float f = length(pos.xz) - 0.3; float verticalConstraint = abs(pos.y + 0.10) - 0.18; f = max(f, verticalConstraint); return traceVec2ToMat4(vec2(f, float(HAT_MATERIAL))); } mat4 hat2(vec3 pos){ pos.xyz -= vec3(0.0,-0.3,0.0); pos.y += 0.01 * cos(pos.x * 10.0); pos.y += 0.5 * cos(pos.x * 2.4) * pos.y; float r = length(pos.xz); float f = max(r - 0.7, 0.0); pos.y -= 0.5 * pow(1.1 * r, 2.0) + 0.4 * (2.0 * pow(1.0 - clamp(r, 0.3, 1.0), 5.0)); float verticalConstraint = max(abs(pos.y) - 0.003, 0.0); f = max(f, verticalConstraint); return traceVec2ToMat4(vec2(f, HAT_MATERIAL)); } void hat_Material(inout vec4 col, vec3 p, float AOLight, vec3 r, float t, mat4 trace_data) { vec4 hat = vec4(0.0); // Fog function hat.rgb = 0.7 * vec3(0.25,0.13,0.04); hat.rgb *= 0.6; col = hat; col *= 0.2 + AOLight * 1.0; col.a += 2.0/t; col = clamp(col, vec4(0.0), vec4(1.0)); } mat4 wire(vec3 pos){ // This part has to match hat2 pos.xyz -= vec3(0.0,-0.3,0.0); pos.y += 0.01 * cos(pos.x * 10.0); pos.y += 0.5 * cos(pos.x * 2.4) * pos.y; float r = length(pos.xz); pos.y -= 0.5 * pow(1.1 * r, 2.0) + 0.4 * (2.0 * pow(1.0 - clamp(r, 0.3, 1.0), 5.0)); float a = atan(pos.z, pos.x); float offsetNormal = 0.001 - 0.013 * sin(a * 30.0); float offsetUp = 0.001 + 0.014 * cos(a * 30.0); float wireDistNormal = abs(r - 0.7 + offsetNormal); float wireDistUp = abs(pos.y + offsetUp); float wireDist = max(length(vec2(wireDistNormal, wireDistUp)) - 0.004,0.0); return traceVec2ToMat4(vec2(wireDist, HAT_MATERIAL)); } // https://iquilezles.org/articles/distfunctions/ // Thanks I.Quilez for this function! float opSmoothUnion( float d1, float d2, float k ) { float h = clamp( 0.5 + 0.5*(d2-d1)/k, 0.0, 1.0 ); return mix(d2, d1, h ) - k*h*(1.0-h); } // https://iquilezles.org/articles/distfunctions/ // Thanks I.Quilez for this function! float sdOctahedron( vec3 p, float s) { p = abs(p); float m = p.x+p.y+p.z-s; vec3 q; if( 3.0*p.x < m ) q = p.xyz; else if( 3.0*p.y < m ) q = p.yzx; else if( 3.0*p.z < m ) q = p.zxy; else return m*0.57735027; float k = clamp(0.5*(q.z-q.y+s),0.0,s); return length(vec3(q.x,q.y-s+k,q.z-k)); } // https://iquilezles.org/articles/distfunctions/ // Thanks I.Quilez for this function! float opSmoothSubtraction( float d1, float d2, float k ) { float h = clamp( 0.5 - 0.5*(d2+d1)/k, 0.0, 1.0 ); return mix( d2, -d1, h ) + k*h*(1.0-h); } // https://iquilezles.org/articles/distfunctions/ // Thanks I.Quilez for this function! 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)); } float dEar(vec3 p) { p.x *= 10.0; p.x -= 0.4 * p.y; p.x -= 2.0 * p.z; float d = length(p) - 0.1; d += - 0.02 - 0.02 * cos(31.0 * length(p.zy) + 4.0); return d; } vec2 rotateXY(vec2 p, float a){ mat2 rotation_xy = mat2( cos(a), -sin(a), sin(a), cos(a) ); return rotation_xy * p; } float dLips(vec3 p) { float d = 0.0; vec3 lipsP = p + vec3(0.0,0.3,0.25); lipsP.z += 0.02 * cos(p.x * 10.0); lipsP.y += 0.01 * cos(p.x * 10.0); lipsP.xy = rotateXY(lipsP.xy, 1.58); float r = 0.06 * ( 0.08 - 2.4 * pow(abs(p.x), 2.0)); float lipTop = sdCappedCylinder(lipsP, r, 0.14); float lipBottom = sdCappedCylinder(lipsP + vec3(0.016 - 0.2 * (1.0 - cos(3.0*p.x)),0.0,0.0), r, 0.135); d = lipTop; d = opSmoothUnion(d, lipBottom, 0.001); return d; } mat4 face(vec3 p){ p.y += 0.4; // Distored position vec3 dp = p * vec3(1.3,0.8,1.5); float dSphere = length(dp) - 0.41; float d = dSphere; // p.z *= 3.0 + 3.0 * abs(cos(10.0 * p.y)); // p.x *= 0.4 + 0.9 * cos(p.y * 10.0); vec3 q = abs(dp * vec3(1.0,0.9,1.0)) - 0.17; float dBox = length(max(q, 0.0)) - 0.18; //d = dSphere; d = opSmoothUnion(d, dBox, 0.07); // For the nose, I'm using: // Octahedron - bound (not exact) // Clipped by a box // https://iquilezles.org/articles/distfunctions/ vec3 noseP = p - vec3(0.0,-0.19,-0.25); float dNose; noseP.y *= 0.3; noseP.y += 0.1 * noseP.z; dNose = sdOctahedron(noseP, 0.09); q = abs(noseP + vec3(0.0,0.3,0.0)) - 0.3; float noseBox = length(max(q, 0.0)) + min(max(q.x, max(q.y,q.z)), 0.0); dNose = opSmoothSubtraction(noseBox, dNose, 0.02) - 0.003; // Smooth union nose & face d = opSmoothUnion(d, dNose, 0.01); // Ears float earR = dEar(p - vec3(0.33,0.0,0.0)); float earL = dEar((p + vec3(0.33,0.0,0.0)) * vec3(-1.0,1.0,1.0)); float ear = min(earR, earL); d = opSmoothUnion(d, ear, 0.01); // Mouth float lipsDistance = dLips(p); d = opSmoothUnion(d, lipsDistance, 0.01); // Chin vec3 pChin = p + vec3(0.0,0.41,0.14); pChin.x *= 0.7; float dChin = length(pChin) - 0.12; d = opSmoothUnion(d, dChin, 0.036); return mat4(d, FACE_MATERIAL, 0.0, lipsDistance, vec4(0.0), vec4(0.0), vec4(0.0)); } void face_Material(inout vec4 col, vec3 p, float AOLight, vec3 r, float t, mat4 trace_data) { float d = trace_data[0][0]; vec4 faceCol = vec4(0.8, 0.65, 0.46, 1.0); faceCol.rgb *= 0.6 + 0.3 * pow(AOLight, 0.7); // SubSurface Estimation 1: // high res texture + low res texture mixed together, but // Low Res has more red channel impact float lf = 200.0; float hf = lf * 4.0; float a = atan(p.x, p.z); float h = p.y; float highRes = 0.2 * abs(cos(h * hf) * cos(a * hf)); float lowRes = abs(cos(h * lf)); float blendFactor = 0.1; faceCol.rgb *= (1.0 - blendFactor) + blendFactor/2.0 * highRes; faceCol.rgb += blendFactor/2.0 * lowRes * vec3(1.0, 0.1, 0.1); // SubSurface Estimation 2: // Thin parts let more lights come through, and mostly red // Depth estimation // Using a point far, hopefully at the other side // (considering object size << 300.0 * length from p to o) vec3 p1 = p; vec3 p2 = 0.3 * normalize(r) + p1; float sdfP2 = face(p2)[0].x; float dp2p1 = length(p2 - p1); float depth = dp2p1 - sdfP2; faceCol.rgb = clamp(faceCol.rgb, 0.0, 1.0); float ao = 1.0 * pow(clamp(1.0 - depth/0.3, 0.0, 1.0), 2.0); faceCol.r += ao * 1.0; faceCol.g += ao * 0.5; faceCol.b += ao * 0.5; float lipsDistance = trace_data[0][3]; float f = 1.0 - clamp(lipsDistance, 0.0, 0.01) / 0.01; faceCol.rgb += f * vec3(-0.0, -0.1, 0.1); faceCol.rgb *= 0.9; faceCol.a = 2.0/t; col = faceCol * faceCol.a; } vec3 eye_Offset = vec3(0.11, 0.37, 0.27); void eye_Material(inout vec4 col, vec3 p, float AOLight, vec3 r, float t, mat4 trace_data) { vec3 offset = trace_data[1].xyz; p += offset; float pupil = clamp(1.0 - (length(p.xy) - 0.018) / 0.002, 0.0, 1.0); float iris = clamp(1.0 - (length(p.xy) - 0.035) / 0.002, 0.0, 1.0); col.rgb = vec3(0.5,0.2,0.1) * iris + (1.0 - iris) * vec3(1.0); col.rgb *= 1.0 - pupil; col.rgb *= col.a = 1.0; } mat4 eyeLeft(vec3 p){ p += eye_Offset; p.z *= 2.5; float d = length(p) - 0.07; return mat4(d, float(EYE_MATERIAL), 0.0, 0.0, eye_Offset, 0.0, vec4(0.0), vec4(0.0)); } mat4 eyeRight(vec3 p){ vec3 eyeRightOffset = eye_Offset; eyeRightOffset.x *= -1.0; p += eyeRightOffset; p.z *= 2.5; float d = length(p) - 0.07; return mat4(d, float(EYE_MATERIAL), 0.0, 0.0, eyeRightOffset, 0.0, vec4(0.0), vec4(0.0)); } mat4 DISTANCE_FROM_ALL_OBJECTS(vec3 p) { mat4 trace_data[] = mat4[]( hat1(p), hat2(p), wire(p), face(p), eyeLeft(p), eyeRight(p) ); float curr_min = 10e3; float d; float curr_min_index = -1.0; float curr_min_material = -1.0; mat4 current; mat4 min_data; for(int i = 0; i < trace_data.length(); i++) { current = trace_data[i]; d = current[0][0]; if (d < curr_min) { min_data = current; curr_min = d; curr_min_material = current[0][1]; curr_min_index = float(i); } } min_data[0][2] = curr_min_index; return min_data; } vec3 randomize(vec3 p) { float s = 1.0 * mod((cos(p.x * 1120.0 + p.y * 1100.0 + p.z * 1022.0)), 1.0); p.x += 0.2 * mod(s + 1.0, 1.0) - 0.5; p.y += 0.2 * mod(s + 2.0, 1.0) - 0.5; p.z += 0.2 * mod(s + 3.0, 1.0) - 0.5; return p; } #define MAX_ITERATIONS 130 #define MARCH_SPEED 0.27 vec3 worldTransform(vec3 p) { p.y += 0.0; vec2 offset = vec2(0.0, 0.0); p.xz -= vec2(offset); float a = atan(p.z, p.x); float l = length(p.xz); a -= 0.0; a += time * 6.2832; p.x = l * cos(a); p.z = l * sin(a); p.xz += vec2(offset); return p; } /* Walk the ray through the scene */ vec4 trace(vec3 o, vec3 r){ vec3 p; float t = 0.0; float d; int i = 0; mat4 trace_data; for(; i < MAX_ITERATIONS; i++){ p = o + r * t; vec3 noise = vec3( cos(p.x * 300333.0 + float(i)) + cos(p.y * 10000.0), cos(p.y * 33020.0 + float(i)) * cos(p.y * 2e3), cos(p.z * 33000.0 + float(i)) + cos(p.z * 2e5) ); // Noise reduces banding artifacts p.xyz += noise * 0.0002; trace_data = DISTANCE_FROM_ALL_OBJECTS(p); d = trace_data[0][0]; t += MARCH_SPEED * d; if(d < 0.002){ break; } } int material = int(trace_data[0][1]); vec4 col = vec4(0.0); float AOLight = 1.0 / (float(i)/float(MAX_ITERATIONS)); switch (material) { case HAT_MATERIAL: hat_Material(col, p, AOLight, r, t, trace_data); break; case FACE_MATERIAL: face_Material(col, p, AOLight, r, t, trace_data); break; case EYE_MATERIAL: eye_Material(col, p, AOLight, r, t, trace_data); break; } col.a = clamp(col.a, 0.0, 1.0); vec3 background = vec3(0.6, 0.5, 0.45); col.rgb = col.rgb * col.a + (1.0 - col.a) * background; col.a = 1.0; return col; } void mainImage( out vec4 fragColor, in vec2 fragCoord ) { vec2 uv = fragCoord.xy / iResolution.xy; float ratio = iResolution.x / iResolution.y; uv.x *= ratio; vec4 col = vec4(1.0); vec2 pos = uv - vec2(0.5 * ratio, 0.5); // Create ray (with position and origin) vec3 r,o; float a = 0.7 * cos(time * 6.2832); mat2 rotation_xz = mat2( cos(a), -sin(a), sin(a), cos(a) ); r = normalize(vec3(pos, 1.0)); o = vec3(0.0, -0.4, -2.0); r.xz = rotation_xz * r.xz; o.xz = rotation_xz * o.xz; col = trace(o, r); fragColor = col; } void main(){ vec2 uv = UV.xy * iResolution.xy; vec4 col; mainImage(col, uv); col.rgba += 1.0 * (1.0 - col.a); out_color = col; }