symmetry_axis 0.1.2 → 0.1.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 7ac00b6ca3f12a6cf9ff5f8f34ef6e440bbc8e6d
-  data.tar.gz: 0da763b5764c2e3098aabce450ce456f1074f31d
+  metadata.gz: ec1f5b69d4fc36f4c0c526330060b996ed5c5718
+  data.tar.gz: 91850dc12093fc4b9fb4adb0174807fda74e7658
 SHA512:
-  metadata.gz: 4bad29fbf34ede0252c1db5ec83af74531e3aa1ea9c8c32ab723c287095ea0f0faf26215443cd88480e861e17ab73a4a401ca60db1e2e924b9bb5655771383fd
-  data.tar.gz: cabbcf950fb5891d326775907f34c86c1193dea00f94f0da442175a6420c77b8327af5d05861aa975a718b746f0832d867478b92615f43e895c9e845bcb434e9
+  metadata.gz: 838a5ebaca4fbe9137b3dd4f91bae665ca0385a088a35b2e1199b224570f9d86e430be2940d45981909cc169dfa9982b6fb876370c6c31da2e3a5888aa7f10e8
+  data.tar.gz: b0255602df7641adb7eeb17e4ffe2e8edf2f9fd812efce39e868dcb07cdc8677c603cd7cd83cfcbca6174879e3946a5e99e1e813eadad084d451efd8a12f1803

data/lib/symmetry_axis/finder.rb

@@ -0,0 +1,85 @@
+module SymmetryAxis
+  # The algorithm, based on principle, every point have a pair.
+  # It uses the first got point and iterates over +(N - 1)+
+  # other points to find a line,
+  # that another points lie by twos on lines parallel to.
+  # The perpendicular line to this line is the symmetry axis.
+  class Finder
+    attr_reader :line
+
+    def initialize(point_array)
+      @input = point_array
+      return if @input.size.odd?
+      @impossible_angles = []
+      @two_points_cursor = 1
+      loop do
+        return unless take_two_points # step 1
+        break if points_are_on_parallel_lines # step 2
+      end
+      find_axis # step 3
+    end
+
+    private
+
+    def take_two_points
+      @two_points = [@input.first]
+      (@two_points_cursor...(@input.size)).each do |i|
+        rad = Line.for_points(@two_points[0], @input[i]).angle_rad
+        next if impossible?(rad)
+        choose_point i
+        return true
+      end
+      false
+    end
+
+    def choose_point(i)
+      @two_points_cursor = i + 1
+      @two_points[1] = @input[i]
+      @two_points_line = Line.for_points(@two_points[0], @input[i])
+    end
+
+    def impossible?(rad)
+      return false if @impossible_angles.empty?
+      impossible_rad = @impossible_angles.min_by { |e| (e - rad).abs }
+      (rad - impossible_rad).abs < (2 * Line.EPS)
+    end
+
+    def points_are_on_parallel_lines
+      other_points = @input.select { |p| !@two_points.include?(p) }
+      @axis_perpendicular_points = @two_points.dup
+      check_perpendiculars(other_points)
+    end
+
+    def check_perpendiculars(points)
+      loop do
+        break if points.empty?
+        line = @two_points_line.parallel(*points.first)
+        size0 = points.size
+        on_line = points.select { |p2| line.on_me?(*p2) }
+        points -= on_line
+        return false if size0 - points.size < 2
+        save_perpendicular_points(on_line, line)
+      end
+      true
+    end
+
+    def save_perpendicular_points(points, line)
+      return if line != @two_points_line
+      points += @two_points
+      @axis_perpendicular_points = pppair(points, line)
+    end
+
+    def pppair(points, line)
+      if line.vertical?
+        [points.min_by { |p| p[1] }, points.max_by { |p| p[1] }]
+      else
+        [points.min_by { |p| p[0] }, points.max_by { |p| p[0] }]
+      end
+    end
+
+    def find_axis
+      p1, p2 = *@axis_perpendicular_points
+      @line = Line.symmetryof(p1[0], p1[1], p2[0], p2[1])
+    end
+  end
+end

data/lib/symmetry_axis/line.rb

@@ -1,126 +1,128 @@
-# Every line here defined by the coefficients of the equation:
-#     y = a + bx
-# But this equation can not describe vertical lines (infinite +b+).
-# So the vertical lines are defined by +x+.
-class Line
-  attr_reader :a, :b, :x
-  EPS, TOO_BIG_B = 0.000001, 100_000_000
-
-  def initialize(x1, y1, x2, y2)
-    @vertical = nil
-    if x1 < x2
-      resolve_coefs(x1, y1, x2, y2)
-    else
-      resolve_coefs(x2, y2, x1, y1)
+module SymmetryAxis
+  # The lines here are defined by the coefficients of the equation:
+  #     y = a + bx
+  # But this equation can not describe vertical lines (infinite +b+).
+  # So the vertical lines are defined by +x+.
+  class Line
+    attr_reader :a, :b, :x
+    EPS, TOO_BIG_B = 0.000001, 100_000_000
+
+    def initialize(x1, y1, x2, y2)
+      @vertical = nil
+      if x1 < x2
+        resolve_coefs(x1, y1, x2, y2)
+      else
+        resolve_coefs(x2, y2, x1, y1)
+      end
     end
-  end
 
-  # Shorthand for constructor.
-  # @return [Line]
-  def self.for_points(p1, p2)
-    Line.new(p1[0], p1[1], p2[0], p2[1])
-  end
+    # Shorthand for constructor.
+    # @return [Line]
+    def self.for_points(p1, p2)
+      Line.new(p1[0], p1[1], p2[0], p2[1])
+    end
 
-  def ==(other)
-    self.class.aeq(self, other) && self.class.beq(self, other)
-  end
+    def ==(other)
+      self.class.aeq(self, other) && self.class.beq(self, other)
+    end
 
-  # @return [String]
-  def to_s
-    if !vertical?
-      "(a=#{@a}, b=#{@b})"
-    else
-      "(x=#{@x})"
+    # @return [String]
+    def to_s
+      if !vertical?
+        "(a=#{@a}, b=#{@b})"
+      else
+        "(x=#{@x})"
+      end
     end
-  end
 
-  # Is the line strictly vertical. This line could be created from
-  # almost vertical pair of points (+b+ coefficient was just Too Big).
-  #
-  # If +true+, you cannot use +a+ and +b+ fields
-  # (you can read +x+ field instead).
-  def vertical?
-    @vertical = (nil == @b) || (TOO_BIG_B <= @b.abs) if nil == @vertical
-    @vertical
-  end
+    # Is the line strictly vertical. This line could be created from
+    # almost vertical pair of points (+b+ coefficient was just Too Big).
+    #
+    # If +true+, you cannot use +a+ and +b+ fields
+    # (you can read +x+ field instead).
+    def vertical?
+      @vertical = (nil == @b) || (TOO_BIG_B <= @b.abs) if nil == @vertical
+      @vertical
+    end
 
-  # Zero means three o'clock, +pi/2+ means noon.
-  # @return [Float]
-  def angle_rad
-    return (Math::PI / 2) if vertical?
-    Math.atan(@b)
-  end
+    # Zero means three o'clock, +pi/2+ means noon.
+    # @return [Float]
+    def angle_rad
+      return (Math::PI / 2) if vertical?
+      Math.atan(@b)
+    end
 
-  # Zero means three o'clock, 90 means noon.
-  # @return [Float]
-  def angle_degrees
-    angle_rad / Math::PI * 180
-  end
+    # Zero means three o'clock, 90 means noon.
+    # @return [Float]
+    def angle_degrees
+      angle_rad / Math::PI * 180
+    end
 
-  # Does the point lie on the line.
-  def on_me?(x, y)
-    if vertical?
-      (@x - x).abs < EPS
-    else
-      f = @a + @b * x
-      (f - y).abs < on_me_eps(x)
+    # Does the point lie on the line.
+    def on_me?(x, y)
+      if vertical?
+        (@x - x).abs < EPS
+      else
+        f = @a + @b * x
+        (f - y).abs < on_me_eps(x)
+      end
     end
-  end
 
-  def on_me_eps(x)
-    e = EPS
-    e += (EPS / @b).abs if @b != 0
-    e += (EPS / x).abs if x != 0
-    e
-  end
+    def on_me_eps(x)
+      e = EPS
+      e += (EPS / @b).abs if @b != 0
+      e += (EPS / x).abs if x != 0
+      e
+    end
 
-  # Symmetry line of two points.
-  # The perpendicular line drawn through the center
-  # of the line segment between the points.
-  # @return [Line]
-  def self.symmetryof(x1, y1, x2, y2)
-    line = Line.new(x1, y1, x2, y2)
-    x, y = (x1 + x2) / 2, (y1 + y2) / 2
-    if line.vertical?
-      Line.new(0, y, 1, y)
-    else
-      r = line.angle_rad + (Math::PI / 2)
-      return Line.new(x, 0, x, 1) if Math::PI / 2 == r.abs
-      return Line.new(x, y, x + 1, y + Math.tan(r))
+    # Symmetry line of two points.
+    # The perpendicular line drawn through the center
+    # of the line segment between the points.
+    # @return [Line]
+    def self.symmetryof(x1, y1, x2, y2)
+      line = Line.new(x1, y1, x2, y2)
+      x, y = (x1 + x2) / 2, (y1 + y2) / 2
+      if line.vertical?
+        Line.new(0, y, 1, y)
+      else
+        r = line.angle_rad + (Math::PI / 2)
+        return Line.new(x, 0, x, 1) if Math::PI / 2 == r.abs
+        return Line.new(x, y, x + 1, y + Math.tan(r))
+      end
     end
-  end
 
-  # Draw a parallel line through the point.
-  # @return [Line]
-  def parallel(x, y)
-    return Line.new(x, y, x, y + 1) if vertical?
-    Line.new(x, y, x + 1, y + @b)
-  end
+    # Draw a parallel line through the point.
+    # @return [Line]
+    def parallel(x, y)
+      return Line.new(x, y, x, y + 1) if vertical?
+      Line.new(x, y, x + 1, y + @b)
+    end
 
-  private :on_me_eps
+    private :on_me_eps
 
-  private
+    private
 
-  def resolve_coefs(xmin, yxmin, xmax, yxmax)
-    if 0 == (xmax - xmin).abs
-      @b, @a = nil, yxmin
-    else
-      @b = (yxmax - yxmin) / Float(xmax - xmin)
-      @a = yxmin - @b * xmin
+    def resolve_coefs(xmin, yxmin, xmax, yxmax)
+      if 0 == (xmax - xmin).abs
+        @b, @a = nil, yxmin
+      else
+        @b = (yxmax - yxmin) / Float(xmax - xmin)
+        @a = yxmin - @b * xmin
+      end
+      @x, @a = xmin, 0 if vertical?
     end
-    @x, @a = xmin, 0 if vertical?
-  end
 
-  def self.aeq(o1, o2)
-    (o1.a - o2.a).abs < (2 * EPS)
-  end
+    def self.aeq(o1, o2)
+      (o1.a - o2.a).abs < (2 * EPS)
+    end
 
-  def self.beq(o1, o2)
-    if o1.vertical? || o2.vertical?
-      return false unless o1.vertical? && o2.vertical?
-      (o1.x - o2.x).abs < (2 * EPS)
-    else
-      (o1.b - o2.b).abs < (2 * EPS)
+    def self.beq(o1, o2)
+      if o1.vertical? || o2.vertical?
+        return false unless o1.vertical? && o2.vertical?
+        (o1.x - o2.x).abs < (2 * EPS)
+      else
+        (o1.b - o2.b).abs < (2 * EPS)
+      end
     end
   end
 end

data/lib/symmetry_axis.rb

@@ -1,85 +1,11 @@
-# The algorithm, based on principle, every point have a pair.
-# It uses the first got point and iterates over +(N - 1)+ 
-# other points to find a line,
-# all points lie by twos on lines parallel to.
-# The perpendicular line to this line is the symmetry axis.
-class SymmetryAxis
-  attr_reader :line
-
-  def initialize(point_array)
-    @input = point_array
-    return if @input.size.odd?
-    @impossible_angles = []
-    @two_points_cursor = 1
-    loop do
-      return unless take_two_points # step 1
-      break if points_are_on_parallel_lines # step 2
-    end
-    find_axis # step 3
-  end
-
-  private
-
-  def take_two_points
-    @two_points = [@input.first]
-    (@two_points_cursor...(@input.size)).each do |i|
-      rad = Line.for_points(@two_points[0], @input[i]).angle_rad
-      next if impossible?(rad)
-      choose_point i
-      return true
-    end
-    false
-  end
-
-  def choose_point(i)
-    @two_points_cursor = i + 1
-    @two_points[1] = @input[i]
-    @two_points_line = Line.for_points(@two_points[0], @input[i])
-  end
-
-  def impossible?(rad)
-    return false if @impossible_angles.empty?
-    impossible_rad = @impossible_angles.min_by { |e| (e - rad).abs }
-    (rad - impossible_rad).abs < (2 * Line.EPS)
-  end
-
-  def points_are_on_parallel_lines
-    other_points = @input.select { |p| !@two_points.include?(p) }
-    @axis_perpendicular_points = @two_points.dup
-    check_perpendiculars(other_points)
-  end
-
-  def check_perpendiculars(points)
-    loop do
-      break if points.empty?
-      line = @two_points_line.parallel(*points.first)
-      size0 = points.size
-      on_line = points.select { |p2| line.on_me?(*p2) }
-      points -= on_line
-      return false if size0 - points.size < 2
-      save_perpendicular_points(on_line, line)
-    end
-    true
-  end
-
-  def save_perpendicular_points(points, line)
-    return if line != @two_points_line
-    points += @two_points
-    @axis_perpendicular_points = pppair(points, line)
-  end
-
-  def pppair(points, line)
-    if line.vertical?
-      [points.min_by { |p| p[1] }, points.max_by { |p| p[1] }]
-    else
-      [points.min_by { |p| p[0] }, points.max_by { |p| p[0] }]
-    end
-  end
-
-  def find_axis
-    p1, p2 = *@axis_perpendicular_points
-    @line = Line.symmetryof(p1[0], p1[1], p2[0], p2[1])
+module SymmetryAxis
+  # Shorthand method for
+  #     Finder.new(points).line
+  # @return [Line]
+  def self.of(points)
+    Finder.new(points).line
   end
 end
 
 require 'symmetry_axis/line'
+require 'symmetry_axis/finder'

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: symmetry_axis
 version: !ruby/object:Gem::Version
-  version: 0.1.2
+  version: 0.1.3
 platform: ruby
 authors:
 - "Георгий Устинов"
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2015-08-29 00:00:00.000000000 Z
+date: 2015-09-02 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: rake
@@ -62,6 +62,7 @@ extensions: []
 extra_rdoc_files: []
 files:
 - lib/symmetry_axis.rb
+- lib/symmetry_axis/finder.rb
 - lib/symmetry_axis/line.rb
 homepage: https://github.com/georgy7/symmetry_axis
 licenses: