3rb 0.1.0

data/examples/14_utah_teapot.rb

@@ -0,0 +1,291 @@
+#!/usr/bin/env ruby
+# frozen_string_literal: true
+
+# Utah Teapot
+# Created by Martin Newell at the University of Utah in 1975.
+# Public domain icon in computer graphics.
+
+require_relative "../lib/3rb"
+
+class UtahTeapot
+  PATCHES = [
+    # Rim (top of body)
+    [[1.4, 2.25, 0.0], [1.3375, 2.38125, 0.0], [1.4375, 2.38125, 0.0], [1.5, 2.25, 0.0],
+     [1.4, 2.25, 0.784], [1.3375, 2.38125, 0.749], [1.4375, 2.38125, 0.805], [1.5, 2.25, 0.84],
+     [0.784, 2.25, 1.4], [0.749, 2.38125, 1.3375], [0.805, 2.38125, 1.4375], [0.84, 2.25, 1.5],
+     [0.0, 2.25, 1.4], [0.0, 2.38125, 1.3375], [0.0, 2.38125, 1.4375], [0.0, 2.25, 1.5]],
+    [[0.0, 2.25, 1.4], [0.0, 2.38125, 1.3375], [0.0, 2.38125, 1.4375], [0.0, 2.25, 1.5],
+     [-0.784, 2.25, 1.4], [-0.749, 2.38125, 1.3375], [-0.805, 2.38125, 1.4375], [-0.84, 2.25, 1.5],
+     [-1.4, 2.25, 0.784], [-1.3375, 2.38125, 0.749], [-1.4375, 2.38125, 0.805], [-1.5, 2.25, 0.84],
+     [-1.4, 2.25, 0.0], [-1.3375, 2.38125, 0.0], [-1.4375, 2.38125, 0.0], [-1.5, 2.25, 0.0]],
+    [[-1.4, 2.25, 0.0], [-1.3375, 2.38125, 0.0], [-1.4375, 2.38125, 0.0], [-1.5, 2.25, 0.0],
+     [-1.4, 2.25, -0.784], [-1.3375, 2.38125, -0.749], [-1.4375, 2.38125, -0.805], [-1.5, 2.25, -0.84],
+     [-0.784, 2.25, -1.4], [-0.749, 2.38125, -1.3375], [-0.805, 2.38125, -1.4375], [-0.84, 2.25, -1.5],
+     [0.0, 2.25, -1.4], [0.0, 2.38125, -1.3375], [0.0, 2.38125, -1.4375], [0.0, 2.25, -1.5]],
+    [[0.0, 2.25, -1.4], [0.0, 2.38125, -1.3375], [0.0, 2.38125, -1.4375], [0.0, 2.25, -1.5],
+     [0.784, 2.25, -1.4], [0.749, 2.38125, -1.3375], [0.805, 2.38125, -1.4375], [0.84, 2.25, -1.5],
+     [1.4, 2.25, -0.784], [1.3375, 2.38125, -0.749], [1.4375, 2.38125, -0.805], [1.5, 2.25, -0.84],
+     [1.4, 2.25, 0.0], [1.3375, 2.38125, 0.0], [1.4375, 2.38125, 0.0], [1.5, 2.25, 0.0]],
+    # Body upper
+    [[1.5, 2.25, 0.0], [1.75, 1.725, 0.0], [2.0, 1.2, 0.0], [2.0, 0.75, 0.0],
+     [1.5, 2.25, 0.84], [1.75, 1.725, 0.98], [2.0, 1.2, 1.12], [2.0, 0.75, 1.12],
+     [0.84, 2.25, 1.5], [0.98, 1.725, 1.75], [1.12, 1.2, 2.0], [1.12, 0.75, 2.0],
+     [0.0, 2.25, 1.5], [0.0, 1.725, 1.75], [0.0, 1.2, 2.0], [0.0, 0.75, 2.0]],
+    [[0.0, 2.25, 1.5], [0.0, 1.725, 1.75], [0.0, 1.2, 2.0], [0.0, 0.75, 2.0],
+     [-0.84, 2.25, 1.5], [-0.98, 1.725, 1.75], [-1.12, 1.2, 2.0], [-1.12, 0.75, 2.0],
+     [-1.5, 2.25, 0.84], [-1.75, 1.725, 0.98], [-2.0, 1.2, 1.12], [-2.0, 0.75, 1.12],
+     [-1.5, 2.25, 0.0], [-1.75, 1.725, 0.0], [-2.0, 1.2, 0.0], [-2.0, 0.75, 0.0]],
+    [[-1.5, 2.25, 0.0], [-1.75, 1.725, 0.0], [-2.0, 1.2, 0.0], [-2.0, 0.75, 0.0],
+     [-1.5, 2.25, -0.84], [-1.75, 1.725, -0.98], [-2.0, 1.2, -1.12], [-2.0, 0.75, -1.12],
+     [-0.84, 2.25, -1.5], [-0.98, 1.725, -1.75], [-1.12, 1.2, -2.0], [-1.12, 0.75, -2.0],
+     [0.0, 2.25, -1.5], [0.0, 1.725, -1.75], [0.0, 1.2, -2.0], [0.0, 0.75, -2.0]],
+    [[0.0, 2.25, -1.5], [0.0, 1.725, -1.75], [0.0, 1.2, -2.0], [0.0, 0.75, -2.0],
+     [0.84, 2.25, -1.5], [0.98, 1.725, -1.75], [1.12, 1.2, -2.0], [1.12, 0.75, -2.0],
+     [1.5, 2.25, -0.84], [1.75, 1.725, -0.98], [2.0, 1.2, -1.12], [2.0, 0.75, -1.12],
+     [1.5, 2.25, 0.0], [1.75, 1.725, 0.0], [2.0, 1.2, 0.0], [2.0, 0.75, 0.0]],
+    # Body lower
+    [[2.0, 0.75, 0.0], [2.0, 0.3, 0.0], [1.5, 0.075, 0.0], [1.5, 0.0, 0.0],
+     [2.0, 0.75, 1.12], [2.0, 0.3, 1.12], [1.5, 0.075, 0.84], [1.5, 0.0, 0.84],
+     [1.12, 0.75, 2.0], [1.12, 0.3, 2.0], [0.84, 0.075, 1.5], [0.84, 0.0, 1.5],
+     [0.0, 0.75, 2.0], [0.0, 0.3, 2.0], [0.0, 0.075, 1.5], [0.0, 0.0, 1.5]],
+    [[0.0, 0.75, 2.0], [0.0, 0.3, 2.0], [0.0, 0.075, 1.5], [0.0, 0.0, 1.5],
+     [-1.12, 0.75, 2.0], [-1.12, 0.3, 2.0], [-0.84, 0.075, 1.5], [-0.84, 0.0, 1.5],
+     [-2.0, 0.75, 1.12], [-2.0, 0.3, 1.12], [-1.5, 0.075, 0.84], [-1.5, 0.0, 0.84],
+     [-2.0, 0.75, 0.0], [-2.0, 0.3, 0.0], [-1.5, 0.075, 0.0], [-1.5, 0.0, 0.0]],
+    [[-2.0, 0.75, 0.0], [-2.0, 0.3, 0.0], [-1.5, 0.075, 0.0], [-1.5, 0.0, 0.0],
+     [-2.0, 0.75, -1.12], [-2.0, 0.3, -1.12], [-1.5, 0.075, -0.84], [-1.5, 0.0, -0.84],
+     [-1.12, 0.75, -2.0], [-1.12, 0.3, -2.0], [-0.84, 0.075, -1.5], [-0.84, 0.0, -1.5],
+     [0.0, 0.75, -2.0], [0.0, 0.3, -2.0], [0.0, 0.075, -1.5], [0.0, 0.0, -1.5]],
+    [[0.0, 0.75, -2.0], [0.0, 0.3, -2.0], [0.0, 0.075, -1.5], [0.0, 0.0, -1.5],
+     [1.12, 0.75, -2.0], [1.12, 0.3, -2.0], [0.84, 0.075, -1.5], [0.84, 0.0, -1.5],
+     [2.0, 0.75, -1.12], [2.0, 0.3, -1.12], [1.5, 0.075, -0.84], [1.5, 0.0, -0.84],
+     [2.0, 0.75, 0.0], [2.0, 0.3, 0.0], [1.5, 0.075, 0.0], [1.5, 0.0, 0.0]],
+    # Handle
+    [[-1.6, 1.875, 0.0], [-2.3, 1.875, 0.0], [-2.7, 1.875, 0.0], [-2.7, 1.65, 0.0],
+     [-1.6, 1.875, 0.3], [-2.3, 1.875, 0.3], [-2.7, 1.875, 0.3], [-2.7, 1.65, 0.3],
+     [-1.5, 2.1, 0.3], [-2.5, 2.1, 0.3], [-3.0, 2.1, 0.3], [-3.0, 1.65, 0.3],
+     [-1.5, 2.1, 0.0], [-2.5, 2.1, 0.0], [-3.0, 2.1, 0.0], [-3.0, 1.65, 0.0]],
+    [[-1.5, 2.1, 0.0], [-2.5, 2.1, 0.0], [-3.0, 2.1, 0.0], [-3.0, 1.65, 0.0],
+     [-1.5, 2.1, -0.3], [-2.5, 2.1, -0.3], [-3.0, 2.1, -0.3], [-3.0, 1.65, -0.3],
+     [-1.6, 1.875, -0.3], [-2.3, 1.875, -0.3], [-2.7, 1.875, -0.3], [-2.7, 1.65, -0.3],
+     [-1.6, 1.875, 0.0], [-2.3, 1.875, 0.0], [-2.7, 1.875, 0.0], [-2.7, 1.65, 0.0]],
+    [[-2.7, 1.65, 0.0], [-2.7, 1.425, 0.0], [-2.5, 0.975, 0.0], [-2.0, 0.75, 0.0],
+     [-2.7, 1.65, 0.3], [-2.7, 1.425, 0.3], [-2.5, 0.975, 0.3], [-2.0, 0.75, 0.3],
+     [-3.0, 1.65, 0.3], [-3.0, 1.2, 0.3], [-2.65, 0.7875, 0.3], [-1.9, 0.45, 0.3],
+     [-3.0, 1.65, 0.0], [-3.0, 1.2, 0.0], [-2.65, 0.7875, 0.0], [-1.9, 0.45, 0.0]],
+    [[-3.0, 1.65, 0.0], [-3.0, 1.2, 0.0], [-2.65, 0.7875, 0.0], [-1.9, 0.45, 0.0],
+     [-3.0, 1.65, -0.3], [-3.0, 1.2, -0.3], [-2.65, 0.7875, -0.3], [-1.9, 0.45, -0.3],
+     [-2.7, 1.65, -0.3], [-2.7, 1.425, -0.3], [-2.5, 0.975, -0.3], [-2.0, 0.75, -0.3],
+     [-2.7, 1.65, 0.0], [-2.7, 1.425, 0.0], [-2.5, 0.975, 0.0], [-2.0, 0.75, 0.0]],
+    # Spout
+    [[1.7, 1.275, 0.0], [2.6, 1.275, 0.0], [2.3, 1.95, 0.0], [2.7, 2.25, 0.0],
+     [1.7, 1.275, 0.66], [2.6, 1.275, 0.66], [2.3, 1.95, 0.25], [2.7, 2.25, 0.25],
+     [1.7, 0.45, 0.66], [3.1, 0.675, 0.66], [2.4, 1.875, 0.25], [3.3, 2.25, 0.25],
+     [1.7, 0.45, 0.0], [3.1, 0.675, 0.0], [2.4, 1.875, 0.0], [3.3, 2.25, 0.0]],
+    [[1.7, 0.45, 0.0], [3.1, 0.675, 0.0], [2.4, 1.875, 0.0], [3.3, 2.25, 0.0],
+     [1.7, 0.45, -0.66], [3.1, 0.675, -0.66], [2.4, 1.875, -0.25], [3.3, 2.25, -0.25],
+     [1.7, 1.275, -0.66], [2.6, 1.275, -0.66], [2.3, 1.95, -0.25], [2.7, 2.25, -0.25],
+     [1.7, 1.275, 0.0], [2.6, 1.275, 0.0], [2.3, 1.95, 0.0], [2.7, 2.25, 0.0]],
+    [[2.7, 2.25, 0.0], [2.8, 2.325, 0.0], [2.9, 2.325, 0.0], [2.8, 2.25, 0.0],
+     [2.7, 2.25, 0.25], [2.8, 2.325, 0.25], [2.9, 2.325, 0.15], [2.8, 2.25, 0.15],
+     [3.3, 2.25, 0.25], [3.525, 2.34375, 0.25], [3.45, 2.3625, 0.15], [3.2, 2.25, 0.15],
+     [3.3, 2.25, 0.0], [3.525, 2.34375, 0.0], [3.45, 2.3625, 0.0], [3.2, 2.25, 0.0]],
+    [[3.3, 2.25, 0.0], [3.525, 2.34375, 0.0], [3.45, 2.3625, 0.0], [3.2, 2.25, 0.0],
+     [3.3, 2.25, -0.25], [3.525, 2.34375, -0.25], [3.45, 2.3625, -0.15], [3.2, 2.25, -0.15],
+     [2.7, 2.25, -0.25], [2.8, 2.325, -0.25], [2.9, 2.325, -0.15], [2.8, 2.25, -0.15],
+     [2.7, 2.25, 0.0], [2.8, 2.325, 0.0], [2.9, 2.325, 0.0], [2.8, 2.25, 0.0]],
+    # Lid knob
+    [[0.0, 3.0, 0.0], [0.8, 3.0, 0.0], [0.0, 2.7, 0.0], [0.2, 2.55, 0.0],
+     [0.0, 3.0, 0.0], [0.8, 3.0, 0.45], [0.0, 2.7, 0.0], [0.2, 2.55, 0.112],
+     [0.0, 3.0, 0.0], [0.45, 3.0, 0.8], [0.0, 2.7, 0.0], [0.112, 2.55, 0.2],
+     [0.0, 3.0, 0.0], [0.0, 3.0, 0.8], [0.0, 2.7, 0.0], [0.0, 2.55, 0.2]],
+    [[0.0, 3.0, 0.0], [0.0, 3.0, 0.8], [0.0, 2.7, 0.0], [0.0, 2.55, 0.2],
+     [0.0, 3.0, 0.0], [-0.45, 3.0, 0.8], [0.0, 2.7, 0.0], [-0.112, 2.55, 0.2],
+     [0.0, 3.0, 0.0], [-0.8, 3.0, 0.45], [0.0, 2.7, 0.0], [-0.2, 2.55, 0.112],
+     [0.0, 3.0, 0.0], [-0.8, 3.0, 0.0], [0.0, 2.7, 0.0], [-0.2, 2.55, 0.0]],
+    [[0.0, 3.0, 0.0], [-0.8, 3.0, 0.0], [0.0, 2.7, 0.0], [-0.2, 2.55, 0.0],
+     [0.0, 3.0, 0.0], [-0.8, 3.0, -0.45], [0.0, 2.7, 0.0], [-0.2, 2.55, -0.112],
+     [0.0, 3.0, 0.0], [-0.45, 3.0, -0.8], [0.0, 2.7, 0.0], [-0.112, 2.55, -0.2],
+     [0.0, 3.0, 0.0], [0.0, 3.0, -0.8], [0.0, 2.7, 0.0], [0.0, 2.55, -0.2]],
+    [[0.0, 3.0, 0.0], [0.0, 3.0, -0.8], [0.0, 2.7, 0.0], [0.0, 2.55, -0.2],
+     [0.0, 3.0, 0.0], [0.45, 3.0, -0.8], [0.0, 2.7, 0.0], [0.112, 2.55, -0.2],
+     [0.0, 3.0, 0.0], [0.8, 3.0, -0.45], [0.0, 2.7, 0.0], [0.2, 2.55, -0.112],
+     [0.0, 3.0, 0.0], [0.8, 3.0, 0.0], [0.0, 2.7, 0.0], [0.2, 2.55, 0.0]],
+    # Lid
+    [[0.2, 2.55, 0.0], [0.4, 2.4, 0.0], [1.3, 2.4, 0.0], [1.3, 2.25, 0.0],
+     [0.2, 2.55, 0.112], [0.4, 2.4, 0.224], [1.3, 2.4, 0.728], [1.3, 2.25, 0.728],
+     [0.112, 2.55, 0.2], [0.224, 2.4, 0.4], [0.728, 2.4, 1.3], [0.728, 2.25, 1.3],
+     [0.0, 2.55, 0.2], [0.0, 2.4, 0.4], [0.0, 2.4, 1.3], [0.0, 2.25, 1.3]],
+    [[0.0, 2.55, 0.2], [0.0, 2.4, 0.4], [0.0, 2.4, 1.3], [0.0, 2.25, 1.3],
+     [-0.112, 2.55, 0.2], [-0.224, 2.4, 0.4], [-0.728, 2.4, 1.3], [-0.728, 2.25, 1.3],
+     [-0.2, 2.55, 0.112], [-0.4, 2.4, 0.224], [-1.3, 2.4, 0.728], [-1.3, 2.25, 0.728],
+     [-0.2, 2.55, 0.0], [-0.4, 2.4, 0.0], [-1.3, 2.4, 0.0], [-1.3, 2.25, 0.0]],
+    [[-0.2, 2.55, 0.0], [-0.4, 2.4, 0.0], [-1.3, 2.4, 0.0], [-1.3, 2.25, 0.0],
+     [-0.2, 2.55, -0.112], [-0.4, 2.4, -0.224], [-1.3, 2.4, -0.728], [-1.3, 2.25, -0.728],
+     [-0.112, 2.55, -0.2], [-0.224, 2.4, -0.4], [-0.728, 2.4, -1.3], [-0.728, 2.25, -1.3],
+     [0.0, 2.55, -0.2], [0.0, 2.4, -0.4], [0.0, 2.4, -1.3], [0.0, 2.25, -1.3]],
+    [[0.0, 2.55, -0.2], [0.0, 2.4, -0.4], [0.0, 2.4, -1.3], [0.0, 2.25, -1.3],
+     [0.112, 2.55, -0.2], [0.224, 2.4, -0.4], [0.728, 2.4, -1.3], [0.728, 2.25, -1.3],
+     [0.2, 2.55, -0.112], [0.4, 2.4, -0.224], [1.3, 2.4, -0.728], [1.3, 2.25, -0.728],
+     [0.2, 2.55, 0.0], [0.4, 2.4, 0.0], [1.3, 2.4, 0.0], [1.3, 2.25, 0.0]]
+  ].freeze
+
+  def self.create(tessellation: 8, wireframe: false)
+    positions = []
+    normals = []
+    indices = []
+    vertex_count = 0
+
+    PATCHES.each do |patch|
+      tessellate_patch(patch, tessellation, positions, normals, indices, vertex_count)
+      vertex_count = positions.length / 3
+    end
+
+    geometry = Sunrb::BufferGeometry.new
+    geometry.set_attribute(:position, Sunrb::Float32BufferAttribute.new(positions, 3))
+    geometry.set_attribute(:normal, Sunrb::Float32BufferAttribute.new(normals, 3))
+    geometry.set_index(indices)
+
+    material = Sunrb::MeshStandardMaterial.new(
+      color: 0xc0c0c0,
+      metalness: 0.3,
+      roughness: 0.4,
+      wireframe: wireframe
+    )
+
+    mesh = Sunrb::Mesh.new(geometry, material)
+    mesh.name = "utah_teapot"
+    mesh.scale.set(0.5, 0.5, 0.5)
+    mesh.position.y = -0.75
+    mesh
+  end
+
+  def self.tessellate_patch(patch, divisions, positions, normals, indices, base_vertex)
+    points = []
+    normal_vectors = []
+
+    (0..divisions).each do |i|
+      u = i.to_f / divisions
+      (0..divisions).each do |j|
+        v = j.to_f / divisions
+        point = evaluate_bezier(patch, u, v)
+        normal = calculate_normal(patch, u, v)
+        points << point
+        normal_vectors << normal
+      end
+    end
+
+    points.each { |p| positions.push(p[0], p[1], p[2]) }
+    normal_vectors.each { |n| normals.push(n[0], n[1], n[2]) }
+
+    (0...divisions).each do |i|
+      (0...divisions).each do |j|
+        a = base_vertex + i * (divisions + 1) + j
+        b = a + 1
+        c = a + (divisions + 1)
+        d = c + 1
+
+        indices.push(a, c, b)
+        indices.push(b, c, d)
+      end
+    end
+  end
+
+  def self.evaluate_bezier(patch, u, v)
+    bu = bezier_basis(u)
+    bv = bezier_basis(v)
+
+    x = y = z = 0.0
+    4.times do |i|
+      4.times do |j|
+        idx = i * 4 + j
+        weight = bu[i] * bv[j]
+        x += patch[idx][0] * weight
+        y += patch[idx][1] * weight
+        z += patch[idx][2] * weight
+      end
+    end
+
+    [x, y, z]
+  end
+
+  def self.bezier_basis(t)
+    t2 = t * t
+    t3 = t2 * t
+    mt = 1.0 - t
+    mt2 = mt * mt
+    mt3 = mt2 * mt
+
+    [mt3, 3.0 * mt2 * t, 3.0 * mt * t2, t3]
+  end
+
+  def self.calculate_normal(patch, u, v)
+    epsilon = 0.001
+    u0 = [u - epsilon, 0.0].max
+    u1 = [u + epsilon, 1.0].min
+    v0 = [v - epsilon, 0.0].max
+    v1 = [v + epsilon, 1.0].min
+
+    p0 = evaluate_bezier(patch, u0, v)
+    p1 = evaluate_bezier(patch, u1, v)
+    p2 = evaluate_bezier(patch, u, v0)
+    p3 = evaluate_bezier(patch, u, v1)
+
+    du = [p1[0] - p0[0], p1[1] - p0[1], p1[2] - p0[2]]
+    dv = [p3[0] - p2[0], p3[1] - p2[1], p3[2] - p2[2]]
+
+    nx = du[1] * dv[2] - du[2] * dv[1]
+    ny = du[2] * dv[0] - du[0] * dv[2]
+    nz = du[0] * dv[1] - du[1] * dv[0]
+
+    len = Math.sqrt(nx * nx + ny * ny + nz * nz)
+    return [0, 1, 0] if len < 0.0001
+
+    [nx / len, ny / len, nz / len]
+  end
+end
+
+scene = Sunrb::Scene.new
+scene.background = Sunrb::Color.new(0x2c3e50)
+
+camera = Sunrb::PerspectiveCamera.new(fov: 45, aspect: 800.0 / 600.0, near: 0.1, far: 1000)
+camera.position.set(4, 3, 4)
+camera.look_at(Sunrb::Vector3.new(0, 0, 0))
+
+ambient = Sunrb::AmbientLight.new(color: 0x404040)
+scene.add(ambient)
+
+key_light = Sunrb::DirectionalLight.new(color: 0xffffff, intensity: 0.8)
+key_light.position.set(5, 5, 5)
+scene.add(key_light)
+
+fill_light = Sunrb::DirectionalLight.new(color: 0x4488ff, intensity: 0.3)
+fill_light.position.set(-5, 3, -5)
+scene.add(fill_light)
+
+teapot_solid = UtahTeapot.create(tessellation: 16, wireframe: false)
+teapot_solid.position.set(-1.5, 0, 0)
+scene.add(teapot_solid)
+
+teapot_wire = UtahTeapot.create(tessellation: 8, wireframe: true)
+teapot_wire.position.set(1.5, 0, 0)
+scene.add(teapot_wire)
+
+floor_geo = Sunrb::PlaneGeometry.new(width: 20, height: 20)
+floor_mat = Sunrb::MeshStandardMaterial.new(color: 0x1a1a2e, roughness: 0.9)
+floor = Sunrb::Mesh.new(floor_geo, floor_mat)
+floor.rotation.x = -Math::PI / 2
+floor.position.y = -1.5
+scene.add(floor)
+
+puts "=== Utah Teapot (Bezier Patches) ==="
+puts "Left: Solid  |  Right: Wireframe"
+puts "Press ESC to exit"
+
+renderer = Sunrb::OpenGLRenderer.new(width: 800, height: 600, title: "Utah Teapot - Bezier Patches")
+
+renderer.run do |delta|
+  teapot_solid.rotation.y += delta * 0.3
+  teapot_wire.rotation.y += delta * 0.3
+  renderer.render(scene, camera)
+end

data/examples/15_stanford_bunny.rb

@@ -0,0 +1,200 @@
+#!/usr/bin/env ruby
+# frozen_string_literal: true
+
+# Stanford Bunny - Neural Network SDF using SIREN
+# Based on https://www.shadertoy.com/view/wtVyWK
+
+require_relative "../lib/3rb"
+
+VERTEX_SHADER = <<~GLSL
+  #version 330 core
+
+  layout(location = 0) in vec3 position;
+  layout(location = 2) in vec2 uv;
+
+  out vec2 vUv;
+
+  void main() {
+    vUv = uv;
+    gl_Position = vec4(position.xy * 2.0, 0.0, 1.0);
+  }
+GLSL
+
+FRAGMENT_SHADER = <<~GLSL
+  #version 330 core
+
+  in vec2 vUv;
+
+  uniform float iTime;
+  uniform vec2 iResolution;
+  uniform vec3 iMouse;
+
+  out vec4 fragColor;
+
+  vec3 erot(vec3 p, vec3 ax, float ro) {
+    return mix(dot(p, ax) * ax, p, cos(ro)) + sin(ro) * cross(ax, p);
+  }
+
+  float scene(vec3 p) {
+    if (length(p) > 1.0) {
+      return length(p) - 0.8;
+    }
+
+    vec4 f00 = sin(p.y * vec4(-3.02, 1.95, -3.42, -0.60) + p.z * vec4(3.08, 0.85, -2.25, -0.24) - p.x * vec4(-0.29, 1.16, -3.74, 2.89) + vec4(-0.71, 4.50, -3.24, -3.50));
+    vec4 f01 = sin(p.y * vec4(-0.40, -3.61, 3.23, -0.14) + p.z * vec4(-0.36, 3.64, -3.91, 2.66) - p.x * vec4(2.90, -0.54, -2.75, 2.71) + vec4(7.02, -5.41, -1.12, -7.41));
+    vec4 f02 = sin(p.y * vec4(-1.77, -1.28, -4.29, -3.20) + p.z * vec4(-3.49, -2.81, -0.64, 2.79) - p.x * vec4(3.15, 2.14, -3.85, 1.83) + vec4(-2.07, 4.49, 5.33, -2.17));
+    vec4 f03 = sin(p.y * vec4(-0.49, 0.68, 3.05, 0.42) + p.z * vec4(-2.87, 0.78, 3.78, -3.41) - p.x * vec4(-2.65, 0.33, 0.07, -0.64) + vec4(-3.24, -5.90, 1.14, -4.71));
+
+    vec4 f10 = sin(mat4(-0.34, 0.06, -0.59, -0.76, 0.10, -0.19, -0.12, 0.44, 0.64, -0.02, -0.26, 0.15, -0.16, 0.21, 0.91, 0.15) * f00 +
+        mat4(0.01, 0.54, -0.77, 0.11, 0.06, -0.14, 0.43, 0.51, -0.18, 0.08, 0.39, 0.20, 0.33, -0.49, -0.10, 0.19) * f01 +
+        mat4(0.27, 0.22, 0.43, 0.53, 0.18, -0.17, 0.23, -0.64, -0.14, 0.02, -0.10, 0.16, -0.13, -0.06, -0.04, -0.36) * f02 +
+        mat4(-0.13, 0.29, -0.29, 0.08, 1.13, 0.02, -0.83, 0.32, -0.32, 0.04, -0.31, -0.16, 0.14, -0.03, -0.20, 0.39) * f03 +
+        vec4(0.73, -4.28, -1.56, -1.80)) / 1.0 + f00;
+
+    vec4 f11 = sin(mat4(-1.11, 0.55, -0.12, -1.00, 0.16, 0.15, -0.30, 0.31, -0.01, 0.01, 0.31, -0.42, -0.29, 0.38, -0.04, 0.71) * f00 +
+        mat4(0.96, -0.02, 0.86, 0.52, -0.14, 0.60, 0.44, 0.43, 0.02, -0.15, -0.49, -0.05, -0.06, -0.25, -0.03, -0.22) * f01 +
+        mat4(0.52, 0.44, -0.05, -0.11, -0.56, -0.10, -0.61, -0.40, -0.04, 0.55, 0.32, -0.07, -0.02, 0.28, 0.26, -0.49) * f02 +
+        mat4(0.02, -0.32, 0.06, -0.17, -0.59, 0.00, -0.24, 0.60, -0.06, 0.13, -0.21, -0.27, -0.12, -0.14, 0.58, -0.55) * f03 +
+        vec4(-2.24, -3.48, -0.80, 1.41)) / 1.0 + f01;
+
+    vec4 f12 = sin(mat4(0.44, -0.06, -0.79, -0.46, 0.05, -0.60, 0.30, 0.36, 0.35, 0.12, 0.02, 0.12, 0.40, -0.26, 0.63, -0.21) * f00 +
+        mat4(-0.48, 0.43, -0.73, -0.40, 0.11, -0.01, 0.71, 0.05, -0.25, 0.25, -0.28, -0.20, 0.32, -0.02, -0.84, 0.16) * f01 +
+        mat4(0.39, -0.07, 0.90, 0.36, -0.38, -0.27, -1.86, -0.39, 0.48, -0.20, -0.05, 0.10, -0.00, -0.21, 0.29, 0.63) * f02 +
+        mat4(0.46, -0.32, 0.06, 0.09, 0.72, -0.47, 0.81, 0.78, 0.90, 0.02, -0.21, 0.08, -0.16, 0.22, 0.32, -0.13) * f03 +
+        vec4(3.38, 1.20, 0.84, 1.41)) / 1.0 + f02;
+
+    vec4 f13 = sin(mat4(-0.41, -0.24, -0.71, -0.25, -0.24, -0.75, -0.09, 0.02, -0.27, -0.42, 0.02, 0.03, -0.01, 0.51, -0.12, -1.24) * f00 +
+        mat4(0.64, 0.31, -1.36, 0.61, -0.34, 0.11, 0.14, 0.79, 0.22, -0.16, -0.29, -0.70, 0.02, -0.37, 0.49, 0.39) * f01 +
+        mat4(0.79, 0.47, 0.54, -0.47, -1.13, -0.35, -1.03, -0.22, -0.67, -0.26, 0.10, 0.21, -0.07, -0.73, -0.11, 0.72) * f02 +
+        mat4(0.43, -0.23, 0.13, 0.09, 1.38, -0.63, 1.57, -0.20, 0.39, -0.14, 0.42, 0.13, -0.57, -0.08, -0.21, 0.21) * f03 +
+        vec4(-0.34, -3.28, 0.43, -0.52)) / 1.0 + f03;
+
+    f00 = sin(mat4(-0.72, 0.23, -0.89, 0.52, 0.38, 0.19, -0.16, -0.88, 0.26, -0.37, 0.09, 0.63, 0.29, -0.72, 0.30, -0.95) * f10 +
+        mat4(-0.22, -0.51, -0.42, -0.73, -0.32, 0.00, -1.03, 1.17, -0.20, -0.03, -0.13, -0.16, -0.41, 0.09, 0.36, -0.84) * f11 +
+        mat4(-0.21, 0.01, 0.33, 0.47, 0.05, 0.20, -0.44, -1.04, 0.13, 0.12, -0.13, 0.31, 0.01, -0.34, 0.41, -0.34) * f12 +
+        mat4(-0.13, -0.06, -0.39, -0.22, 0.48, 0.25, 0.24, -0.97, -0.34, 0.14, 0.42, -0.00, -0.44, 0.05, 0.09, -0.95) * f13 +
+        vec4(0.48, 0.87, -0.87, -2.06)) / 1.4 + f10;
+
+    f01 = sin(mat4(-0.27, 0.29, -0.21, 0.15, 0.34, -0.23, 0.85, -0.09, -1.15, -0.24, -0.05, -0.25, -0.12, -0.73, -0.17, -0.37) * f10 +
+        mat4(-1.11, 0.35, -0.93, -0.06, -0.79, -0.03, -0.46, -0.37, 0.60, -0.37, -0.14, 0.45, -0.03, -0.21, 0.02, 0.59) * f11 +
+        mat4(-0.92, -0.17, -0.58, -0.18, 0.58, 0.60, 0.83, -1.04, -0.80, -0.16, 0.23, -0.11, 0.08, 0.16, 0.76, 0.61) * f12 +
+        mat4(0.29, 0.45, 0.30, 0.39, -0.91, 0.66, -0.35, -0.35, 0.21, 0.16, -0.54, -0.63, 1.10, -0.38, 0.20, 0.15) * f13 +
+        vec4(-1.72, -0.14, 1.92, 2.08)) / 1.4 + f11;
+
+    f02 = sin(mat4(1.00, 0.66, 1.30, -0.51, 0.88, 0.25, -0.67, 0.03, -0.68, -0.08, -0.12, -0.14, 0.46, 1.15, 0.38, -0.10) * f10 +
+        mat4(0.51, -0.57, 0.41, -0.09, 0.68, -0.50, -0.04, -1.01, 0.20, 0.44, -0.60, 0.46, -0.09, -0.37, -1.30, 0.04) * f11 +
+        mat4(0.14, 0.29, -0.45, -0.06, -0.65, 0.33, -0.37, -0.95, 0.71, -0.07, 1.00, -0.60, -1.68, -0.20, -0.00, -0.70) * f12 +
+        mat4(-0.31, 0.69, 0.56, 0.13, 0.95, 0.36, 0.56, 0.59, -0.63, 0.52, -0.30, 0.17, 1.23, 0.72, 0.95, 0.75) * f13 +
+        vec4(-0.90, -3.26, -0.44, -3.11)) / 1.4 + f12;
+
+    f03 = sin(mat4(0.51, -0.98, -0.28, 0.16, -0.22, -0.17, -1.03, 0.22, 0.70, -0.15, 0.12, 0.43, 0.78, 0.67, -0.85, -0.25) * f10 +
+        mat4(0.81, 0.60, -0.89, 0.61, -1.03, -0.33, 0.60, -0.11, -0.06, 0.01, -0.02, -0.44, 0.73, 0.69, 1.02, 0.62) * f11 +
+        mat4(-0.10, 0.52, 0.80, -0.65, 0.40, -0.75, 0.47, 1.56, 0.03, 0.05, 0.08, 0.31, -0.03, 0.22, -1.63, 0.07) * f12 +
+        mat4(-0.18, -0.07, -1.22, 0.48, -0.01, 0.56, 0.07, 0.15, 0.24, 0.25, -0.09, -0.54, 0.23, -0.08, 0.20, 0.36) * f13 +
+        vec4(-1.11, -4.28, 1.02, -0.23)) / 1.4 + f13;
+
+    return dot(f00, vec4(0.09, 0.12, -0.07, -0.03)) + dot(f01, vec4(-0.04, 0.07, -0.08, 0.05)) +
+        dot(f02, vec4(-0.01, 0.06, -0.02, 0.07)) + dot(f03, vec4(-0.05, 0.07, 0.03, 0.04)) - 0.16;
+  }
+
+  vec3 norm(vec3 p) {
+    mat3 k = mat3(p, p, p) - mat3(0.001);
+    return normalize(scene(p) - vec3(scene(k[0]), scene(k[1]), scene(k[2])));
+  }
+
+  void main() {
+    vec2 fragCoord = vUv * iResolution;
+    vec2 uv = (fragCoord - 0.5 * iResolution.xy) / iResolution.y;
+    vec2 mouse = (iMouse.xy - 0.5 * iResolution.xy) / iResolution.y;
+
+    vec3 cam = normalize(vec3(1.5, uv));
+    vec3 init = vec3(-5.0, 0.0, 0.0);
+
+    float yrot = 0.5;
+    float zrot = iTime * 0.2;
+    if (iMouse.z > 0.0) {
+      yrot += -4.0 * mouse.y;
+      zrot = 4.0 * mouse.x;
+    }
+    cam = erot(cam, vec3(0.0, 1.0, 0.0), yrot);
+    init = erot(init, vec3(0.0, 1.0, 0.0), yrot);
+    cam = erot(cam, vec3(0.0, 0.0, 1.0), zrot);
+    init = erot(init, vec3(0.0, 0.0, 1.0), zrot);
+
+    vec3 p = init;
+    bool hit = false;
+    for (int i = 0; i < 150 && !hit; i++) {
+      float dist = scene(p);
+      hit = dist * dist < 1e-6;
+      p += dist * cam;
+      if (distance(p, init) > 5.0) break;
+    }
+
+    vec3 n = norm(p);
+    vec3 r = reflect(cam, n);
+    vec3 nz = p - erot(p, vec3(1.0), 2.0) + erot(p, vec3(1.0), 4.0);
+    float spec = length(sin(r * 3.5 + sin(nz * 120.0) * 0.15) * 0.4 + 0.6) / sqrt(3.0);
+    spec *= smoothstep(-0.3, 0.2, scene(p + r * 0.2));
+    vec3 col = vec3(0.1, 0.1, 0.12) * spec + pow(spec, 8.0);
+    float bgdot = length(sin(cam * 8.0) * 0.4 + 0.6) / 2.0;
+    vec3 bg = vec3(0.1, 0.1, 0.11) * bgdot + pow(bgdot, 10.0);
+
+    vec4 result = hit ? vec4(col, 1.0) : vec4(bg, 1.0);
+    fragColor = smoothstep(-0.02, 1.05, sqrt(result)) * (1.0 - dot(uv, uv) * 0.5);
+    fragColor.a = 1.0;
+  }
+GLSL
+
+scene = Sunrb::Scene.new
+scene.background = Sunrb::Color.new(0x000000)
+
+camera = Sunrb::OrthographicCamera.new(
+  left: -5, right: 5, top: 5, bottom: -5, near: 0.1, far: 20
+)
+camera.position.z = 8
+
+geometry = Sunrb::PlaneGeometry.new(width: 2, height: 2)
+
+shader_material = Sunrb::ShaderMaterial.new(
+  vertex_shader: VERTEX_SHADER,
+  fragment_shader: FRAGMENT_SHADER,
+  uniforms: {
+    iTime: 0.0,
+    iResolution: Sunrb::Vector2.new(800, 600),
+    iMouse: Sunrb::Vector3.new(0, 0, 0)
+  }
+)
+
+quad = Sunrb::Mesh.new(geometry, shader_material)
+scene.add(quad)
+
+puts "=== Stanford Bunny - Neural Network SDF ==="
+puts "Based on SIREN networks: https://vsitzmann.github.io/siren/"
+puts "Drag mouse to rotate, press ESC to exit"
+
+renderer = Sunrb::OpenGLRenderer.new(width: 800, height: 600, title: "Stanford Bunny - Neural Network SDF")
+
+time = 0.0
+mouse_x = 0.0
+mouse_y = 0.0
+mouse_pressed = false
+
+renderer.on_mouse_move do |x, y|
+  mouse_x = x
+  mouse_y = 600 - y
+end
+
+renderer.on_mouse_button do |button, action, _mods|
+  mouse_pressed = (button == 0 && action == 1)
+end
+
+renderer.run do |delta|
+  time += delta
+  shader_material.uniforms[:iTime] = time
+  shader_material.uniforms[:iMouse] = Sunrb::Vector3.new(
+    mouse_x,
+    mouse_y,
+    mouse_pressed ? 1.0 : 0.0
+  )
+  renderer.render(scene, camera)
+end