symmetry_axis 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA1:
+  metadata.gz: 6c13dde182042eff38651139bade52ae10525937
+  data.tar.gz: 7de66b90f6e33fdb5647cfd2f7442c507e7a6f5a
+SHA512:
+  metadata.gz: 251442207e32769923a89800e3b2d0a8c71518b53a3750146d771c10a537987d25e34209300ae6a30fc67d2aefe684764d564cad7ffd3136381b7d96d2cc7592
+  data.tar.gz: 5897c6c376f89a0c5bb97b66aa52a0285031b4c5c21b0d041b98abcf6b2eb3abb89470603c516a30d5685b3cba8327bc7613ae2ce54205810ca393030c1f6e38

data/lib/symmetry_axis/line.rb

@@ -0,0 +1,101 @@
+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
+
+  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 to_s
+    if !vertical
+      "(a=#{@a}, b=#{@b})"
+    else
+      "(x=#{@x})"
+    end
+  end
+
+  def vertical
+    @vertical = (nil == @b) || (TOO_BIG_B <= @b.abs) if nil == @vertical
+    @vertical
+  end
+
+  def angle_rad
+    return (Math::PI / 2) if vertical
+    Math.atan(@b)
+  end
+
+  def angle_degrees
+    angle_rad / Math::PI * 180
+  end
+
+  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
+
+  def on_me_eps(x)
+    e = EPS
+    e += (EPS / @b).abs if @b != 0
+    e += (EPS / x).abs if x != 0
+    e
+  end
+
+  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
+
+  def parallel(x, y)
+    return Line.new(x, y, x, y + 1) if vertical
+    Line.new(x, y, x + 1, y + @b)
+  end
+
+  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
+    end
+    @x, @a = xmin, 0 if vertical
+  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)
+    end
+  end
+end

data/lib/symmetry_axis.rb

@@ -0,0 +1,81 @@
+class SymmetryAxis
+  attr_reader :line
+
+  def initialize(point_array)
+    @input = point_array
+    @impossible_angles = []
+    @two_points_cursor = 1
+    @s1counter, @s2counter, @s2loopcounter = 0, 0, 0
+    loop do
+      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
+    @s1counter += 1
+    @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
+    end
+    fail 'There is no suitable pair of points.'
+  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
+    @s2counter += 1
+    points = @input.select { |p| !@two_points.include?(p) }
+    @perpendicular_points_pairs_list = [@two_points.dup]
+    loop do
+      break if points.empty?
+      @s2loopcounter += 1
+      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
+    @perpendicular_points_pairs_list[0] = 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
+    # puts "\n   Iterations:\n    Step1 = #{@s1counter}\n    Step2 = #{@s2counter}\n    Step2Loop = #{@s2loopcounter}"
+    points = @perpendicular_points_pairs_list[0]
+    p1, p2 = points[0], points[1]
+    @line = Line.symmetryof(p1[0], p1[1], p2[0], p2[1])
+  end
+end
+
+require 'symmetry_axis/line'

metadata

@@ -0,0 +1,75 @@
+--- !ruby/object:Gem::Specification
+name: symmetry_axis
+version: !ruby/object:Gem::Version
+  version: 0.1.0
+platform: ruby
+authors:
+- "Георгий Устинов"
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2015-08-28 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '10'
+- !ruby/object:Gem::Dependency
+  name: rspec
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '3'
+description: 'It has the restriction: no points must be on axis. Thus every point
+  must have a pair. Also you can separately use the Line class from the gem. It allows
+  you to compare for equality the lines, defined by different points on itself.'
+email: httpcatharsis@gmail.com
+executables: []
+extensions: []
+extra_rdoc_files: []
+files:
+- lib/symmetry_axis.rb
+- lib/symmetry_axis/line.rb
+homepage: https://github.com/georgy7/symmetry_axis
+licenses:
+- MIT
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.4.5
+signing_key: 
+specification_version: 4
+summary: Get axis from a point list.
+test_files: []