// "Ocean Structure" by dr2 - 2017 (Optimized) // License: Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License // Optimizations: Reduced raymarching steps, simplified reflections, removed smoke effects #define REAL_WAVE 0 // real (=1) or fake (=0, less work) water waves #define N_REFL 1 // number of reflections (reduced from 3 to 1) float PrBoxDf (vec3 p, vec3 b); float PrSphAnDf (vec3 p, float r, float w); float PrCylDf (vec3 p, float r, float h); float PrTorusDf (vec3 p, float ri, float rc); float PrCapsDf (vec3 p, float r, float h); float PrCapsAnDf (vec3 p, float r, float w, float h); float PrFlatCylDf (vec3 p, float rhi, float rlo, float h); float Noisefv2 (vec2 p); float Noisefv3 (vec3 p); float Fbm2 (vec2 p); vec3 HsvToRgb (vec3 c); vec2 Rot2D (vec2 q, float a); vec3 qHit, sunDir; float dstFar, tCur, tWav; int idObj; bool isNite; const int idBase = 1, idPlat = 2, isShel = 3, idFrm = 4, idDway = 5, idTwr = 6, idBrg = 7, idBrCab = 8, idRdw = 9; const float pi = 3.14159; float BridgeDf (vec3 p, float dMin) { float cbRad = 0.06; vec3 q, qq; float d, wRd; wRd = 1.; q = p; q.x = abs (q.x) - wRd; q.y -= -0.3; d = PrBoxDf (q, vec3 (0.15, 5., 0.15)); q = p; q.y -= 4.85; d = min (d, PrBoxDf (q, vec3 (wRd + 0.15, 0.15, 0.15))); if (d < dMin) { dMin = d; qHit = q; idObj = idBrg; } qq = p; qq.x = abs (qq.x) - wRd + 0.07; qq.z = abs (qq.z); q = qq; q.y -= 4.92; q.yz = Rot2D (q.yz, -0.1 * pi); q.z -= 9.5; d = PrCylDf (q, cbRad, 9.5); if (d < dMin) { dMin = d; qHit = q; idObj = idBrCab; } return dMin; } float CentStrucDf (vec3 p, float dMin, float ar) { float cRad = 6., cLen = 8., wThk = 0.05, doorHt = 1.6, doorWd = 1.4; vec3 q; vec2 qo; float d, dd; q = p; q.y -= -1.05; d = PrCylDf (q.xzy, 8.5, 0.15); if (d < dMin) { dMin = d; qHit = q; idObj = idPlat; } q = p; qo = Rot2D (q.xz, 2. * pi * (floor (4. * ar) + 0.5) / 4.); q.xz = qo; q.y -= cLen + 1.2 * doorHt - 9.; dd = PrFlatCylDf (q.yzx, doorHt, doorWd, 0.); q = p; q.y -= cLen - 9.; d = max (max (max (PrCapsAnDf (q.xzy, cRad, wThk, cLen), - q.y), q.y - 1.635 * cLen), - dd); if (d < dMin) { dMin = d; qHit = q; idObj = isShel; } return dMin; } float CornStrucDf (vec3 p, float dMin) { vec3 q; float d, a; q = p; q.y -= -1.2; d = PrCylDf (q.xzy, 3.2, 0.15); if (d < dMin) { dMin = d; qHit = q; idObj = idPlat; } q = p; q.y -= 1.; d = max (PrCapsAnDf (q.xzy, 2.5, 0.1, 3.), -2.2 - q.y); q = p; q.y -= 7.; d = min (d, max (PrCapsDf (q.xzy, 0.7, 2.), -1. - q.y)); q = p; a = (length (q.xz) > 0.) ? atan (q.z, - q.x) / (2. * pi) : 0.; q.xz = Rot2D (q.xz, 2. * pi * (floor (4. * a + 0.5) / 4.)); d = max (d, - PrFlatCylDf (q.yzx, 1.5, 1., 0.)); if (d < dMin) { dMin = d; qHit = q; idObj = idTwr; } return dMin; } float ObjDf (vec3 p) { vec3 q; float d, dMin, hs, wRd, ar; dMin = dstFar; hs = 5.; q = p; q.xz = abs (q.xz) - 4.; d = max (PrSphAnDf (q, 4.85, 0.15), - min (3.9 - q.y, q.y)); q.y -= 0.5; d = max (d, 2.2 - min (length (q.yz), length (q.xy))); if (d < dMin) { dMin = d; qHit = q; idObj = idBase; } q = p; q.xz = abs (q.xz) - 20.; d = max (PrSphAnDf (q, 4.85, 0.15), - min (3.9 - q.y, q.y)); q.y -= 0.5; d = max (d, 2.2 - min (length (q.yz), length (q.xy))); if (d < dMin) { dMin = d; qHit = q; idObj = idBase; } wRd = 1.; q = p; q.xz = abs (q.xz) - 20.7 + wRd; d = max (max (length (max (q.xz - wRd, 0.)) - 0.3, - (length (max (q.xz + wRd, 0.)) - 0.3)), abs (q.y - hs + 1.) - 0.1); if (d < dMin) { dMin = d; qHit = q; idObj = idRdw; } q = p; ar = (length (p.xz) > 0.) ? atan (p.z, - p.x) / (2. * pi) : 0.; q.y -= hs; dMin = CentStrucDf (q, dMin, ar); q = p; q.y -= hs; q.xz = abs (q.xz) - vec2 (20.); dMin = CornStrucDf (q, dMin); q = p; q.y -= hs; q.xz = Rot2D (q.xz, 2. * pi * (floor (4. * ar + 0.5) / 4.)); q.x += 20.; dMin = BridgeDf (q, dMin); return dMin; } float ObjRay (vec3 ro, vec3 rd) { float dHit, d; dHit = 0.; for (int j = 0; j < 60; j ++) { // Reduced from 120 to 60 d = ObjDf (ro + dHit * rd); dHit += d; if (d < 0.002 || dHit > dstFar) break; // Increased threshold } return dHit; } vec3 ObjNf (vec3 p) { vec4 v; vec3 e = vec3 (0.002, -0.002, 0.); // Larger epsilon v = vec4 (ObjDf (p + e.xxx), ObjDf (p + e.xyy), ObjDf (p + e.yxy), ObjDf (p + e.yyx)); return normalize (vec3 (v.x - v.y - v.z - v.w) + 2. * v.yzw); } float WaveHt (vec2 p) { mat2 qRot = mat2 (0.8, -0.6, 0.6, 0.8); vec4 t4, v4; vec2 t; float wFreq, wAmp, ht; wFreq = 1.; wAmp = 1.; ht = 0.; for (int j = 0; j < 2; j ++) { // Reduced from 3 to 2 p *= qRot; t = tWav * vec2 (1., -1.); t4 = (p.xyxy + t.xxyy) * wFreq; t = vec2 (Noisefv2 (t4.xy), Noisefv2 (t4.zw)); t4 += 2. * t.xxyy - 1.; v4 = (1. - abs (sin (t4))) * (abs (sin (t4)) + abs (cos (t4))); ht += wAmp * dot (pow (1. - sqrt (v4.xz * v4.yw), vec2 (8.)), vec2 (1.)); wFreq *= 2.; wAmp *= 0.5; } return ht; } vec3 WaveNf (vec3 p, float d) { vec3 vn; vec2 e; e = vec2 (max (0.02, 0.01 * d * d), 0.); // Larger epsilon p *= 0.3; vn.xz = 0.4 * (WaveHt (p.xz) - vec2 (WaveHt (p.xz + e.xy), WaveHt (p.xz + e.yx))); vn.y = e.x; return normalize (vn); } vec4 ObjCol (vec3 n) { vec4 col; if (idObj == idBase) col = vec4 (0.3, 0.4, 0.1, 0.1); else if (idObj == idPlat) col = vec4 (0.4, 0.4, 0.3, 0.1); else if (idObj == isShel) col = vec4 (0.5, 0.5, 0.5, 0.3); else if (idObj == idFrm) col = vec4 (0.8, 0.8, 0.9, 0.5); else if (idObj == idDway) col = vec4 (0.7, 0.8, 0.7, 0.3); else if (idObj == idTwr) col = vec4 (0.7, 0.7, 0.6, 0.3); else if (idObj == idBrg) col = vec4 (1., 0.3, 0.3, 0.3); else if (idObj == idBrCab) col = vec4 (0.9, 0.9, 1., 0.5); else if (idObj == idRdw) col = vec4 (0.4, 0.3, 0.3, 0.1); return col; } vec4 AurCol (vec3 ro, vec3 rd) { vec4 col, mCol; vec3 p, dp; float ar; dp = rd / rd.y; p = ro + (40. - ro.y) * dp; col = vec4 (0.); mCol = vec4 (0.); tWav = 0.05 * tCur; for (float ns = 0.; ns < 25.; ns ++) { // Reduced from 50 to 25 p += dp; ar = 0.05 - clamp (0.06 * WaveHt (0.01 * p.xz), 0., 0.04); mCol = mix (mCol, ar * vec4 (HsvToRgb (vec3 (0.34 + 0.007 * ns, 1., 1. - 0.02 * ns)), 1.), 0.5); col += mCol; } return col; } vec3 NtSkyCol (vec3 rd) { vec3 rds; rds = floor (2000. * rd); rds = 0.00015 * rds + 0.1 * Noisefv3 (0.0005 * rds.yzx); for (int j = 0; j < 10; j ++) rds = abs (rds) / dot (rds, rds) - 0.9; // Reduced from 19 to 10 return 0.3 * vec3 (1., 1., 0.9) * min (1., 0.5e-3 * pow (min (6., length (rds)), 5.)); } vec3 BgCol (vec3 ro, vec3 rd) { vec3 col; float f, a; if (rd.y > 0.) { a = atan (rd.x, - rd.z); if (rd.y < 0.03 * 0.5 + 0.005) // Simplified Fbm1 col = (isNite ? vec3 (0.07, 0.1, 0.07) : vec3 (0.4, 0.5, 0.7)); else { if (isNite) { vec4 aCol = AurCol (ro, rd); col = (1. - 0.5 * aCol.a) * NtSkyCol (rd) + 0.6 * aCol.rgb; } else { ro.xz += 2. * tCur; col = vec3 (0.2, 0.3, 0.6) + 0.2 * (1. - max (rd.y, 0.)) + 0.1 * pow (max (dot (rd, sunDir), 0.), 16.); f = Fbm2 (0.02 * (ro.xz + rd.xz * (100. - ro.y) / max (rd.y, 0.01))); col = mix (col, vec3 (1.), 0.2 + 0.8 * f * rd.y); } } } else { col = vec3 (0.6, 0.5, 0.3); if (- ro.y / rd.y < dstFar) { ro += - (ro.y / rd.y) * rd; col *= 1.1 - 0.2 * Noisefv2 (30. * ro.xz); } col = mix (col, 0.9 * (vec3 (0.4, 0.2, 0.1) + 0.2) + 0.1, pow (1. + rd.y, 5.)); } return col; } vec4 GlowCol (vec3 ro, vec3 rd, float dstObj) { vec3 gloDir; float gloDist, wGlow, s; wGlow = 0.; for (float j = 0.; j < 4.; j ++) { gloDir = vec3 (20., 9.3, 20.) * (1. - 2. * vec3 (floor (j / 2.), 0., mod (j, 2.))) - ro; gloDist = length (gloDir); s = dot (rd, normalize (gloDir)); if (s > 0. && gloDist < dstObj) wGlow += 1. - smoothstep (1., 2., sqrt (1. - s * s) * gloDist); } gloDir = vec3 (0., 15.5, 0.) - ro; gloDist = length (gloDir); s = dot (rd, normalize (gloDir)); if (s > 0. && gloDist < dstObj) wGlow += 1. - smoothstep (2., 3., sqrt (1. - s * s) * gloDist); return (0.6 + 0.4 * sin (0.3 * 2. * pi * tCur)) * clamp (wGlow, 0., 1.) * vec4 (1., 0.5, 0.3, 1.); } vec3 ShowScene (vec3 ro, vec3 rd) { vec4 objCol, glwCol, glwColW; vec3 col, vn; float dstObj, dstWat, reflCol; bool wRefl; col = vec3 (0.2, 0.2, 0.); wRefl = false; dstObj = ObjRay (ro, rd); glwCol = GlowCol (ro, rd, dstObj); glwColW = vec4 (0.); tWav = 0.4 * tCur; reflCol = 1.; if (N_REFL >= 1 && rd.y < 0.) { dstWat = - ro.y / rd.y; if (dstWat < min (dstObj, dstFar)) { wRefl = true; ro += dstWat * rd; vn = WaveNf (ro, dstWat); rd = reflect (rd, vn); ro += 0.01 * rd; dstObj = ObjRay (ro, rd); glwColW = GlowCol (ro, rd, dstObj); reflCol *= 0.7; } } if (dstObj < dstFar) { ro += dstObj * rd; vn = ObjNf (ro); objCol = ObjCol (vn); if (isNite) col = objCol.rgb * vec3 (0.3, 0.35, 0.3) * (0.2 + 0.8 * max (0.,vn.y)); else col = objCol.rgb * (0.2 + 0.8 * max (0., max (dot (vn, sunDir), 0.))) + objCol.a * pow (max (dot (normalize (sunDir - rd), vn), 0.), 64.); } else if (rd.y > 0.) { col = BgCol (ro, rd); } else { dstWat = - ro.y / rd.y; col = BgCol (ro + dstWat * rd, reflect (rd, vec3 (0., 1., 0.))); reflCol = 0.7; } col = clamp (reflCol * col, 0., 1.); col = mix (col, glwCol.rgb, glwCol.a); if (wRefl) { col = mix (col, reflCol * glwColW.rgb, glwColW.a); col = mix (mix (vec3 (0., 0.1, 0.), vec3 (0., 0.05, 0.05), smoothstep (0.4, 0.6, Fbm2 (0.5 * ro.xz))), col, 1. - pow (abs (rd.y), 4.)); } return clamp (col, 0., 1.); } void mainImage (out vec4 fragColor, in vec2 fragCoord) { mat3 vuMat; vec4 mPtr; vec3 ro, rd; vec2 uv, ori, ca, sa; float el, az; uv = 2. * fragCoord.xy / iResolution.xy - 1.; uv.x *= iResolution.x / iResolution.y; tCur = iTime; mPtr = iMouse; mPtr.xy = mPtr.xy / iResolution.xy - 0.5; dstFar = 120.; // Reduced from 140 isNite = true; az = 0.33 * pi; el = -0.016 * pi; if (mPtr.z > 0.) { if (mPtr.x > 0.45 && mPtr.y < -0.45) isNite = false; else { az += pi * mPtr.x; el += 0.05 * pi * mPtr.y; } } else { az += 0.002 * pi * tCur; el += 0.01 * pi * sin (0.01 * pi * tCur); } el = clamp (el, -0.4 * pi, -0.01 * pi); ori = vec2 (el, az); ca = cos (ori); sa = sin (ori); vuMat = mat3 (ca.y, 0., - sa.y, 0., 1., 0., sa.y, 0., ca.y) * mat3 (1., 0., 0., 0., ca.x, - sa.x, 0., sa.x, ca.x); ro = vuMat * vec3 (0., 10., -100.); rd = vuMat * normalize (vec3 (uv, 4.2)); sunDir = vuMat * normalize (vec3 (1., 1., -1.)); fragColor = vec4 (ShowScene (ro, rd), 1.); } float PrBoxDf (vec3 p, vec3 b) { vec3 d; d = abs (p) - b; return min (max (d.x, max (d.y, d.z)), 0.) + length (max (d, 0.)); } float PrSphAnDf (vec3 p, float r, float w) { return abs (length (p) - r) - w; } float PrCylDf (vec3 p, float r, float h) { return max (length (p.xy) - r, abs (p.z) - h); } float PrCapsDf (vec3 p, float r, float h) { return length (p - vec3 (0., 0., h * clamp (p.z / h, -1., 1.))) - r; } float PrCapsAnDf (vec3 p, float r, float w, float h) { p.z = abs (p.z); return max (length (p - vec3 (0., 0., min (p.z, h + w))) - r, - length (p - vec3 (0., 0., min (p.z, h - w))) + r) - w; } float PrFlatCylDf (vec3 p, float rhi, float rlo, float h) { float d; d = length (p.xy - vec2 (rhi * clamp (p.x / rhi, -1., 1.), 0.)) - rlo; if (h > 0.) d = max (d, abs (p.z) - h); return d; } float PrTorusDf (vec3 p, float ri, float rc) { return length (vec2 (length (p.xy) - rc, p.z)) - ri; } const float cHashM = 43758.54; vec2 Hashv2v2 (vec2 p) { vec2 cHashVA2 = vec2 (37., 39.); return fract (sin (vec2 (dot (p, cHashVA2), dot (p + vec2 (1., 0.), cHashVA2))) * cHashM); } vec4 Hashv4v3 (vec3 p) { vec3 cHashVA3 = vec3 (37., 39., 41.); vec2 e = vec2 (1., 0.); return fract (sin (vec4 (dot (p + e.yyy, cHashVA3), dot (p + e.xyy, cHashVA3), dot (p + e.yxy, cHashVA3), dot (p + e.xxy, cHashVA3))) * cHashM); } float Noisefv2 (vec2 p) { vec2 t, ip, fp; ip = floor (p); fp = fract (p); fp = fp * fp * (3. - 2. * fp); t = mix (Hashv2v2 (ip), Hashv2v2 (ip + vec2 (0., 1.)), fp.y); return mix (t.x, t.y, fp.x); } float Noisefv3 (vec3 p) { vec4 t; vec3 ip, fp; ip = floor (p); fp = fract (p); fp *= fp * (3. - 2. * fp); t = mix (Hashv4v3 (ip), Hashv4v3 (ip + vec3 (0., 0., 1.)), fp.z); return mix (mix (t.x, t.y, fp.x), mix (t.z, t.w, fp.x), fp.y); } float Fbm2 (vec2 p) { float f, a; f = 0.; a = 1.; for (int i = 0; i < 3; i ++) { // Reduced from 5 to 3 f += a * Noisefv2 (p); a *= 0.5; p *= 2.; } return f * (1. / 1.75); } vec3 HsvToRgb (vec3 c) { vec3 p; p = abs (fract (c.xxx + vec3 (1., 2./3., 1./3.)) * 6. - 3.); return c.z * mix (vec3 (1.), clamp (p - 1., 0., 1.), c.y); } vec2 Rot2D (vec2 q, float a) { return q * cos (a) + q.yx * sin (a) * vec2 (-1., 1.); }

Use it