mass refactor

Author: monkstone

contributed/circles.rb

@@ -0,0 +1,46 @@
+load_library :circle
+
+def settings
+  size(800, 600, P2D)
+end
+
+def setup
+  sketch_title 'Circles'
+  color_mode(HSB, 360, 100, 100, 100)
+  reset
+  ellipse_mode(RADIUS)
+end
+
+def draw
+  fill(0, 0, 0)
+  no_stroke
+  reset if (frame_count % 8_000).zero?
+  @points.each do |point|
+    # change direction sometimes
+    point.direction Vec2D.random if rand > 0.96
+    point.update
+  end
+
+  # set the style of the circle
+  @dc = map1d(millis, 0..150_000, 0..360) # slowly changes hue
+  stroke((@c + @dc) % 360, 50, 100, 5)
+  no_fill
+
+  ## verifies if there is a circle and draw it
+  draw_circle @points unless @points.collinear?
+end
+
+def draw_circle(pts)
+  circumcircle = Circumcircle.new(@points.positions)
+  circumcircle.calculate
+  center_point = circumcircle.center
+  radius = circumcircle.radius
+  ellipse(center_point.x, center_point.y, radius, radius)
+end
+
+def reset
+  @c = rand(360)
+  @points = TrianglePoints.new
+  3.times { @points << TPoint.new(Vec2D.new(rand(5..width - 5), rand(5..height - 5))) }
+  background 0
+end

contributed/library/circle/circle.rb

@@ -0,0 +1,3 @@
+require_relative 'lib/circumcircle'
+require_relative 'lib/triangle_point'
+require_relative 'lib/t_points'

contributed/library/circle/lib/circumcircle.rb

@@ -0,0 +1,38 @@
+# frozen_string_literal: true
+require 'matrix'
+
+# Circumcircle from 3 points
+class Circumcircle
+  attr_reader :center, :radius, :points
+  def initialize(points)
+    @points = points
+  end
+
+  def calculate
+    @center = Vec2D.new(-(bx / am), -(by / am))
+    @radius = center.dist(points[2]) # points[2] = c
+  end
+
+  private
+
+  # Matrix math see matrix_math.md and in detail
+  # http://mathworld.wolfram.com/Circumcircle.html
+
+  def am
+    2 * Matrix[
+      *points.map { |pt| [pt.x, pt.y, 1] }
+    ].determinant
+  end
+
+  def bx
+    -Matrix[
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.y, 1] }
+    ].determinant
+  end
+
+  def by
+    Matrix[
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.x, 1] }
+    ].determinant
+  end
+end

contributed/library/circle/lib/t_points.rb

@@ -0,0 +1,34 @@
+# frozen_string_literal: true
+require 'forwardable'
+MAX_POINT = 3
+# A collection of a maximum of 3 points in the processing world
+# includes a collinearity test using Vec2D
+class TrianglePoints
+  extend Forwardable
+  def_delegators(:@points, :each, :map, :size, :shift, :clear, :[])
+  include Enumerable
+
+  attr_reader :points
+
+  def initialize
+    @points = []
+  end
+
+  def <<(pt)
+    points << pt
+    shift if size > MAX_POINT
+  end
+
+  def collinear?
+    full? ? (positions[0] - positions[1]).cross(positions[1] - positions[2]).zero? : false
+  end
+
+  # returns positions as an array of Vec2D
+  def positions
+    points.map(&:pos)
+  end
+
+  def full?
+    points.length == MAX_POINT
+  end
+end

contributed/library/circle/lib/triangle_point.rb

@@ -0,0 +1,46 @@
+# frozen_string_literal: true
+# particle and triangle point
+class TPoint
+  include Processing::Proxy
+  attr_reader :pos, :vel, :accel, :xbound, :ybound
+  # attr_reader :width, :height # uncomment for testing
+
+  def initialize(position)
+    @pos = position
+    @vel = Vec2D.new
+    @accel = Vec2D.random
+    @xbound = Boundary.new(0, width)
+    @ybound = Boundary.new(0, height)
+  end
+
+  def direction(acc)
+    # direction of the acceleration is defined by the new angle
+    @accel = acc
+    # magnitude of the acceleration is proportional to the angle between
+    # acceleration and velocity
+    dif = acc.angle_between(vel)
+    dif = map1d(dif, 0..PI, 0.1..0.001)
+    @accel = acc * dif
+  end
+
+  def update
+    @vel += accel
+    @vel.set_mag(1.5) { vel.mag > 1.5 }
+    @pos += vel
+    check_bounds
+  end
+
+  private
+
+  def check_bounds
+    @vel.x *= -1 if xbound.exclude? pos.x
+    @vel.y *= -1 if ybound.exclude? pos.y
+  end
+end
+
+# we are looking for excluded values
+Boundary = Struct.new(:lower, :upper) do
+  def exclude?(val)
+    true unless (lower...upper).cover? val
+  end
+end

processing_app/library/vecmath/vec2d/library/circle/lib/circumcircle.rb

@@ -1,6 +1,7 @@
-# Circumcircle from 3 points
+# frozen_string_literal: true
 require 'matrix'
 
+# Circumcircle from 3 points
 class Circumcircle
   attr_reader :center, :radius, :points
   def initialize(points)
@@ -19,25 +20,19 @@ def calculate
 
   def am
     2 * Matrix[
-      [points[0].x, points[0].y, 1],
-      [points[1].x, points[1].y, 1],
-      [points[2].x, points[2].y, 1]
+      *points.map { |pt| [pt.x, pt.y, 1] }
     ].determinant
   end
 
   def bx
     -Matrix[
-      [points[0].x * points[0].x + points[0].y * points[0].y, points[0].y, 1],
-      [points[1].x * points[1].x + points[1].y * points[1].y, points[1].y, 1],
-      [points[2].x * points[2].x + points[2].y * points[2].y, points[2].y, 1]
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.y, 1] }
     ].determinant
   end
 
   def by
     Matrix[
-      [points[0].x * points[0].x + points[0].y * points[0].y, points[0].x, 1],
-      [points[1].x * points[1].x + points[1].y * points[1].y, points[1].x, 1],
-      [points[2].x * points[2].x + points[2].y * points[2].y, points[2].x, 1]
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.x, 1] }
     ].determinant
   end
 end

processing_app/library/vecmath/vec2d/library/simple_circle/lib/circumcircle.rb

@@ -1,6 +1,7 @@
-# Circumcircle from 3 points
+# frozen_string_literal: true
 require 'matrix'
 
+# Circumcircle from 3 points
 class Circumcircle
   attr_reader :center, :radius, :points
   def initialize(points)
@@ -19,25 +20,19 @@ def calculate
 
   def am
     2 * Matrix[
-      [points[0].x, points[0].y, 1],
-      [points[1].x, points[1].y, 1],
-      [points[2].x, points[2].y, 1]
+      *points.map { |pt| [pt.x, pt.y, 1] }
     ].determinant
   end
 
   def bx
     -Matrix[
-      [points[0].x * points[0].x + points[0].y * points[0].y, points[0].y, 1],
-      [points[1].x * points[1].x + points[1].y * points[1].y, points[1].y, 1],
-      [points[2].x * points[2].x + points[2].y * points[2].y, points[2].y, 1]
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.y, 1] }
     ].determinant
   end
 
   def by
     Matrix[
-      [points[0].x * points[0].x + points[0].y * points[0].y, points[0].x, 1],
-      [points[1].x * points[1].x + points[1].y * points[1].y, points[1].x, 1],
-      [points[2].x * points[2].x + points[2].y * points[2].y, points[2].x, 1]
+      *points.map { |pt| [pt.x * pt.x + pt.y * pt.y, pt.x, 1] }
     ].determinant
   end
 end