// This shader enhances the rendering of a raymarched hexagonal Truchet pattern by incorporating subsurface scattering effects. // Height-based displacement and elaborate shading including Voronoi noise, ambient occlusion, and fake environment mapping are used for intricate shading. precision mediump float; varying vec2 uv; uniform float time, width, height; float heightMap(in vec2 p); float map(vec3 p); vec3 getNormal(vec3 p, inout float edge, inout float crv); float calculateAO(in vec3 p, in vec3 n); float n3D(vec3 p); vec3 envMap(vec3 rd, vec3 sn); vec2 hash22(vec2 p); float Voronoi(in vec2 p); vec3 subsurfaceScattering(in vec3 pos, in vec3 normal, in vec3 lightDir); void mainImage(out vec4 fragColor, in vec2 fragCoord); float heightMap(in vec2 p) { p *= 3.0; vec2 h = vec2(p.x + p.y * 0.57735, p.y * 1.1547); vec2 fh = floor(h); vec2 f = h - fh; h = fh; float c = fract((h.x + h.y) / 3.0); h = c < 0.666 ? c < 0.333 ? h : h + 1.0 : h + step(f.yx, f); p -= vec2(h.x - h.y * 0.5, h.y * 0.8660254); c = fract(cos(dot(h, vec2(41, 289))) * 43758.5453); p -= p * step(c, 0.5) * 2.0; p -= vec2(-1, 0); c = dot(p, p); p -= vec2(1.5, 0.8660254); c = min(c, dot(p, p)); p -= vec2(0, -1.73205); c = min(c, dot(p, p)); return sqrt(c); } float map(vec3 p) { float c = heightMap(p.xy); c = cos(c * 6.2831589) + cos(c * 6.2831589 * 2.0); c = clamp(c * 0.6 + 0.5, 0.0, 1.0); return 1.0 - p.z - c * 0.025; } vec3 getNormal(vec3 p, inout float edge, inout float crv) { vec2 e = vec2(0.01, 0.0); float d1 = map(p + e.xyy), d2 = map(p - e.xyy); float d3 = map(p + e.yxy), d4 = map(p - e.yxy); float d5 = map(p + e.yyx), d6 = map(p - e.yyx); float d = map(p) * 2.0; edge = abs(d1 + d2 - d) + abs(d3 + d4 - d) + abs(d5 + d6 - d); edge = smoothstep(0.0, 1.0, sqrt(edge / e.x * 2.0)); crv = clamp((d1 + d2 + d3 + d4 + d5 + d6 - d * 3.0) * 32.0 + 0.6, 0.0, 1.0); e = vec2(0.0025, 0.0); d1 = map(p + e.xyy), d2 = map(p - e.xyy); d3 = map(p + e.yxy), d4 = map(p - e.yxy); d5 = map(p + e.yyx), d6 = map(p - e.yyx); return normalize(vec3(d1 - d2, d3 - d4, d5 - d6)); } float calculateAO(in vec3 p, in vec3 n) { float sca = 2.0, occ = 0.0; for(float i = 0.0; i < 5.0; i++) { float hr = 0.01 + i * 0.5 / 4.0; float dd = map(n * hr + p); occ += (hr - dd) * sca; sca *= 0.7; } return clamp(1.0 - occ, 0.0, 1.0); } float n3D(vec3 p) { const vec3 s = vec3(7, 157, 113); vec3 ip = floor(p); p -= ip; vec4 h = vec4(0.0, s.yz, s.y + s.z) + dot(ip, s); p = p * p * (3.0 - 2.0 * p); h = mix(fract(sin(mod(h, 6.2831589)) * 43758.5453), fract(sin(mod(h + s.x, 6.2831589)) * 43758.5453), p.x); h.xy = mix(h.xz, h.yw, p.y); return mix(h.x, h.y, p.z); } vec3 envMap(vec3 rd, vec3 sn) { vec3 sRd = rd; rd.xy -= time * 0.25; rd *= 3.0; float c = n3D(rd) * 0.57 + n3D(rd * 2.0) * 0.28 + n3D(rd * 4.0) * 0.15; c = smoothstep(0.4, 1.0, c); vec3 col = vec3(c, c * c, c * c * c * c); return mix(col, col.yzx, sRd * 0.25 + 0.25); } vec2 hash22(vec2 p) { float n = sin(mod(dot(p, vec2(41, 289)), 6.2831589)); return fract(vec2(262144, 32768) * n) * 0.75 + 0.25; } float Voronoi(in vec2 p) { vec2 g = floor(p), o; p -= g; vec3 d = vec3(1.0); for(int y = -1; y <= 1; y++) { for(int x = -1; x <= 1; x++) { o = vec2(float(x), float(y)); o += hash22(g + o) - p; d.z = dot(o, o); d.y = max(d.x, min(d.y, d.z)); d.x = min(d.x, d.z); } } return max(d.y / 1.2 - d.x * 1.0, 0.0) / 1.2; } vec3 subsurfaceScattering(in vec3 pos, in vec3 normal, in vec3 lightDir) { float thickness = 0.03; // Control thickness of subsurface effect vec3 lightOffset = pos + normal * thickness; float lightAmount = dot(lightDir, getNormal(lightOffset, thickness, thickness)); return lightAmount * vec3(1.0, 0.5, 0.5); // Reddish scattering color } void mainImage(out vec4 fragColor, in vec2 fragCoord) { vec2 resolution = vec2(width, height); vec3 rd = normalize(vec3(2.0 * fragCoord - resolution.xy, resolution.y)); float tm = time / 2.0; vec2 a = sin(vec2(1.570796, 0.0) + sin(tm / 4.0) * 0.3); rd.xy = mat2(a, -a.y, a.x) * rd.xy; vec3 ro = vec3(tm, cos(tm / 4.0), 0.0); vec3 lp = ro + vec3(cos(tm / 2.0) * 0.5, sin(tm / 2.0) * 0.5, -0.5); float d, t = 0.0; for(int j = 0; j < 32; j++) { d = map(ro + rd * t); t += d * 0.7; if(d < 0.001) break; } float edge, crv; vec3 sp = ro + rd * t; vec3 sn = getNormal(sp, edge, crv); vec3 ld = lp - sp; float c = heightMap(sp.xy); vec3 fold = cos(vec3(1.0, 2.0, 4.0) * c * 6.2831589); float c2 = heightMap((sp.xy + sp.z * 0.025) * 6.0); c2 = cos(c2 * 6.2831589 * 3.0); c2 = clamp(c2 + 0.5, 0.0, 1.0); vec3 oC = vec3(1.0); if(fold.x > 0.0) oC = vec3(1.0, 0.05, 0.1) * c2; if(fold.x < 0.05 && (fold.y) < 0.0) oC = vec3(1.0, 0.7, 0.45) * (c2 * 0.25 + 0.75); else if(fold.x < 0.0) oC = vec3(1.0, 0.8, 0.4) * c2; float p1 = 1.0 - smoothstep(0.0, 0.1, fold.x * 0.5 + 0.5); float p2 = 1.0 - smoothstep(0.0, 0.1, Voronoi(sp.xy * 4.0 + vec2(tm, cos(tm / 4.0)))); p1 = (p2 + 0.25) * p1; oC += oC.yzx * p1 * p1; float lDist = max(length(ld), 0.001); float atten = 1.0 / (1.0 + lDist * 0.125); ld /= lDist; float diff = max(dot(ld, sn), 0.0); float spec = pow(max(dot(reflect(-ld, sn), -rd), 0.0), 16.0); float fre = pow(clamp(dot(sn, rd) + 1.0, 0.0, 1.0), 3.0); crv = crv * 0.9 + 0.1; float ao = calculateAO(sp, sn); vec3 sss = subsurfaceScattering(sp, sn, ld); vec3 col = oC * (diff + 0.5) + vec3(1.0, 0.7, 0.4) * spec * 2.0 + vec3(0.4, 0.7, 1.0) * fre; col += (oC * 0.5 + 0.5) * envMap(reflect(rd, sn), sn) * 6.0; col *= 1.0 - edge * 0.85; col *= (atten * crv * ao); col += sss * 0.3; // Blend subsurface scattering effect fragColor = vec4(sqrt(clamp(col, 0.0, 1.0)), 1.0); } void main() { vec2 fragCoord = vec2((uv.x - 0.5) * max(width / height, 1.0) + 0.5, (uv.y - 0.5) * max(height / width, 1.0) + 0.5); vec4 color; mainImage(color, fragCoord * vec2(width, height)); gl_FragColor = color; }

teal Teal Tint

edit

edit

tesseract effect Hypercube shadows

edit

tap to change color palette

edit

Subsurface Scattering

edit

edit

edit

like the auto bahn Speed Lines

edit

Gradient Overlay

edit

click to explore, then edit