geometry_3d 0.0.0 → 0.1.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA1:
-  metadata.gz: 2ce961e96e2afd1aa40e1ed69b026f01fe6431f3
-  data.tar.gz: f7694694f3a44e170e00fb11176d259501dd9cab
+  metadata.gz: b6be9663a382709ce88d58a88343d355c8c11889
+  data.tar.gz: 434c0d1cad5c1474523ecdf2e3d025d06cf25a6a
 SHA512:
-  metadata.gz: fd0c38048554b2716e5b8308703e3a894d913cb38e4746b13533ed552f4cc3deb5a68b91ea26ab1ef23afbccdd7f5bdef5d7421888e0118cdca076397d805f56
-  data.tar.gz: 0859769cf7a415e81fa902df66271f584ed86ae1d4c0a3ddb9f092ac010d4a19f66a7cf16d520a5df7867077ad680acbe6c2445d95421975ba5b9c027d18389d
+  metadata.gz: f45cfb2296046a337a1628970c8e7b371791b593d1399088f5a7897c7a9468914d61cce976b455fa74866702de54d45a56669e03dcc49dffab943879c7ec4bf4
+  data.tar.gz: 728cd85b9e419cde91e666a1ba269ef355207eef0fe158a318629b9059f3facebf0917df4af8b52103415c563654421a0a9df6ecfac64eade17c560279051d10

data/bin/geometry_3d

@@ -0,0 +1,2 @@
+#!/usr/bin/env ruby
+require 'geometry_3d'

data/lib/geometry_3d.rb

@@ -1,5 +1,7 @@
-class Geometry
-  def self.hi
-    puts "djani"
-  end
-end
+require_relative './geometry_3d/point.rb'
+require_relative './geometry_3d/vector'
+require_relative './geometry_3d/line'
+require_relative './geometry_3d/plane'
+require_relative './geometry_3d/triangle'
+require_relative './geometry_3d/circle'
+require_relative './geometry_3d/equation.rb'

data/lib/geometry_3d/circle.rb

@@ -0,0 +1,39 @@
+class Circle
+  attr_accessor :point
+  attr_accessor :radius
+  attr_accessor :vector
+
+  def initialize(point, vector, radius)
+    @point = point
+    @radius = radius
+    @vector = vector
+  end
+
+  # Find the area of a circle.
+  #
+  # Example:
+  #   >> Circle.new(Point.new(1, 2, 3), Vector.new(1, 1, 1), 3)
+  #   => 28.274333882308138
+  #
+  # Arguments:
+  #   no
+
+  def area
+    Math::PI * radius ** 2
+  end
+
+  # Find the perimeter of a circle.
+  #
+  # Example:
+  #   >> Circle.new(Point.new(1, 2, 3), Vector.new(1, 1, 1), 3)
+  #   => 18.84955592153876
+  #
+  # Arguments:
+  #   no
+
+  def perimeter
+    2 * Math::PI * radius
+  end
+
+end
+

data/lib/geometry_3d/equation.rb

@@ -0,0 +1,17 @@
+class LinearEquation
+  attr_accessor :a
+  attr_accessor :b
+  attr_accessor :c
+
+  def initialize(a, b, c)
+    @a = a
+    @b = b
+    @c = c
+  end
+
+  def solve_system(equation)
+    y = (equation.a * c - a * equation.c) / (a * equation.b - equation.a * b).to_f
+    x = (- c - b * y) / a.to_f
+    [x, y]
+  end
+end

data/lib/geometry_3d/line.rb

@@ -0,0 +1,49 @@
+class Line
+  attr_accessor :vector
+  attr_accessor :point
+
+  def initialize(point, vector)
+    @point = point
+    @vector = vector
+  end
+
+  # Construct a line with two points.
+  #
+  # Example:
+  #   >> Line.construct_with_two_points(Point.new(1, 2, 3), Point.new(3, 2, 1))
+  #   => Line.new(Point.new(1, 2, 3), Vector.new(2, 0, -2))
+  #
+  # Arguments:
+  #   first_point: (Point)
+  #   second_point: (Point)
+
+  def self.construct_with_two_points(first_point, second_point)
+    v = Vector.construct_with_two_points(first_point, second_point)
+    Line.new(first_point, v)
+  end
+
+  # Find intersecting point of two lines.
+  #
+  # Example:
+  #   >> Line.new(Point.new(1,-2,3), Vector.new(5,2,-1)).find_intersecting_point(Line.new(Point.new(-4,-4,1), Vector.new(-1,3,-1)))
+  #   => Point.new(-4, -4, 4)
+  #
+  # Arguments:
+  #   line: (Line)
+
+  def find_intersecting_point(line)
+    eq1 = LinearEquation.new(vector.x, -line.vector.x, point.x - line.point.x)
+    eq2 = LinearEquation.new(vector.y, -line.vector.y, point.y - line.point.y)
+    find_point (eq1.solve_system(eq2)[0])
+  end
+
+  def ==(line)
+    point == line.point and vector == line.vector
+  end
+
+  private
+  def find_point(parameter)
+    Point.new(*point.to_a.zip(vector.to_a).map { |a, b| a + b * parameter})
+  end
+
+end

data/lib/geometry_3d/plane.rb

@@ -0,0 +1,109 @@
+class Plane
+  attr_accessor :a
+  attr_accessor :b
+  attr_accessor :c
+  attr_accessor :d
+
+  def initialize(a, b, c, d)
+    @a = a
+    @b = b
+    @c = c
+    @d = d
+  end
+
+  # Construct a plane passing through three points.
+  #
+  # Example:
+  #   >> Plane.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9))
+  #   => Plane.new(30, -48, 17, 15)
+  #
+  # Arguments:
+  #   point1: (Point)
+  #   point2: (Point)
+  #   point3: (Point)
+
+  def self.construct_with_three_points(point1, point2, point3)
+    vector1 = Vector.construct_with_two_points(point1, point3)
+    vector2 = Vector.construct_with_two_points(point1, point2)
+    construct_with_two_vectors_and_point(vector1, vector2, point1)
+  end
+
+
+
+  # Construct plane with normal vector and passing through a point.
+  #
+  # Example:
+  #   >> Plane.construct_with_normal_vector_and_point(Vector.new(-2, 3, 1), Point.new(1, 1, 1))
+  #   => Plane.new(-2, 3, 1, -2)
+  #
+  # Arguments:
+  #   vector: (Vector)
+  #   point: (Point)
+
+  def self.construct_with_normal_vector_and_point(vector, point)
+    Plane.new(*vector.to_a, find_d(vector.to_a, point))
+  end
+
+  # Construct a plane parallel to two vectors and passing through a point.
+  #
+  # Example:
+  #   >> Plane.construct_with_two_vectors_and_point(Vector.new(1, 2, 3), Vector.new(-2, 1, 2), Point.new(1, 1, 1))
+  #   => Plane.new(30, -48, 17, 15)
+  #
+  # Arguments:
+  #   vector1: (Vector)
+  #   vector2: (Vector)
+  #   point: (Point)
+
+  def self.construct_with_two_vectors_and_point(vector1, vector2, point)
+    vector = vector1.cross_product(vector2)
+    construct_with_normal_vector_and_point(vector, point)
+  end
+
+  # Construct a plane parallel other plane and passing through a point.
+  #
+  # Example:
+  #   >> Plane.construct_with_plane_and_point(Plane.new(2, -3, 4, 5), Point.new(1, 1, 1))
+  #   => Plane.new(2, -3, 4, -3)
+  #
+  # Arguments:
+  #   vector1: (Vector)
+  #   vector2: (Vector)
+  #   point: (Point)
+
+  def self.construct_with_plane_and_point(plane, point)
+    Plane.new(*plane.to_a.take(3), find_d(plane.to_a.take(3), point))
+  end
+
+  def to_a
+    [a, b, c, d]
+  end
+
+  # Find normal vector of the plane
+  #
+  # Example:
+  #   >> Plane.new(2, -5, 1.5, 3)
+  #   => Vector.new(-2, -5, 1.5)
+  #
+  # Arguments:
+  #   no
+
+  def find_normal_vector()
+    Vector.construct_with_plane(self)
+  end
+
+
+
+  def ==(plane)
+    #a == plane.a and b == plane.b and c == plane.c and d == plane.d
+    to_a.zip(plane.to_a).map { |a, b| a == b}.reduce(:&)
+  end
+
+
+  private
+  def self.find_d(array, point)
+    - array.zip(point.to_a).map { |a, b| a * b }.reduce(:+)
+  end
+
+
+end

data/lib/geometry_3d/point.rb

@@ -0,0 +1,92 @@
+class Point
+  attr_accessor :x
+  attr_accessor :y
+  attr_accessor :z
+
+  def initialize(x, y, z)
+    @x = x
+    @y = y
+    @z = z
+  end
+
+  # Translate a point with a vector
+  #
+  # Example:
+  #   >> Point.new(2, -1, 3).translate(Vector.new(1, 2, 3))
+  #   => Point.new(3, 1, 6)
+  #
+  # Arguments:
+  #   vector: (Vector)
+
+  def translate(vector)
+    Point.new(*to_a.zip(vector.to_a).map { |a, b| a + b})
+  end
+
+  # Check if a point is on a plane
+  #
+  # Example:
+  #   >> Point.new(1, 1, 1).is_on_plane?(Plane.new(1, 2, 3, -6))
+  #   => true
+  #
+  # Arguments:
+  #   plane: (Plane)
+
+  def is_on_plane?(plane)
+    (plane.a * x + plane.b * y + plane.c * z + plane.d) == 0
+  end
+
+  # Check if a point is on a line
+  #
+  # Example:
+  #   >> Point.new(1, 1, 4).is_on_line?(Line.new(Point.new(2, -1, 3), Vector.new(-1, 2, 1)))
+  #   => true
+  #
+  # Arguments:
+  #   line: (Line)
+
+  def is_on_line?(line)
+    to_a.zip(line.point.to_a, line.vector.to_a).map { |a, b, c| (a - b) / c }.uniq.length == 1
+  end
+
+  # Find the distance between two points
+  #
+  # Example:
+  #   >> Point.new(-7, -4, 3).distance_to(Point.new(17, 6, 2.5))
+  #   => 26.004807247891687
+  #
+  # Arguments:
+  #   point: (Point)
+
+  def distance_to(point)
+    to_a.zip(point.to_a).map{|a, b| (b - a) ** 2}.reduce(:+) ** 0.5
+  end
+
+  # Find the midpoint between two points
+  #
+  # Example:
+  #   >> Point.new(-7, -4, 3).get_midpoint(Point.new(17, 6, 2.5))
+  #   => Point.new(5, 1, 2.75)
+  #
+  # Arguments:
+  #   point: (Point)
+
+  def get_midpoint(point)
+    Point.new(*to_a.zip(point.to_a).map { |a, b| (a + b) / 2.0 })
+  end
+
+  def to_a
+    [x, y, z]
+  end
+
+
+  def ==(point)
+    @x == point.x and @y == point.y and @z == point.z
+  end
+
+  def round
+    Point.new(*to_a.map{|x| x.round(4)})
+  end
+
+
+end
+

data/lib/geometry_3d/triangle.rb

@@ -0,0 +1,115 @@
+class Triangle
+  attr_accessor :point_a
+  attr_accessor :point_b
+  attr_accessor :point_c
+
+  def initialize(point_a, point_b, point_c)
+    @point_a = point_a
+    @point_b = point_b
+    @point_c = point_c
+  end
+
+  # Find the plane in which the triangle lies.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_plane
+  #   => Plane.new(30, -48, 17, 15)
+  #
+  # Arguments:
+  #   no
+
+  def get_plane
+    Plane.construct_with_three_points(point_a, point_b, point_c)
+  end
+
+  # Find the line in which the AB side of the triangle lies.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_line_ab
+  #   => Line.new(Point.new(1, 2, 3), Vector.new(3, 4, 6))
+  #
+  # Arguments:
+  #   no
+
+  def get_line_ab
+    Line.construct_with_two_points(point_a, point_b)
+  end
+
+  # Find the line in which the BC side of the triangle lies.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_line_bc
+  #   => Line.new(Point.new(4, 6, 9), Vector.new(8, 5, 0))
+  #
+  # Arguments:
+  #   no
+
+  def get_line_bc
+    Line.construct_with_two_points(point_b, point_c)
+  end
+
+  # Find the line in which the AC side of the triangle lies.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_line_ac
+  #   => Line.new(Point.new(1, 2, 3), Vector.new(11, 9, 6))
+  #
+  # Arguments:
+  #   no
+
+  def get_line_ac
+    Line.construct_with_two_points(point_a, point_c)
+  end
+
+  # Find the median of the triangle from the point A
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_median_a
+  #   => Line.new(Point.new(1, 2, 3), Vector.new(7, 6.5, 6))
+  #
+  # Arguments:
+  #   no
+
+  def get_median_a
+    Line.construct_with_two_points(point_a, point_b.get_midpoint(point_c))
+  end
+
+  # Find the median of the triangle from the point B.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_median_b
+  #   => Line.new(Point.new(4, 6, 9), Vector.new(2.5, 0.5, -3))
+  #
+  # Arguments:
+  #   no
+#Line.new(Point.new(4, 6, 9), Vector.new(2, 0, -3)).find_intersecting_point(Line.new(Point.new(12, 11, 9), Vector.new(-10, -7, -3)))
+  def get_median_b
+    Line.construct_with_two_points(point_b, point_a.get_midpoint(point_c))
+  end
+
+  # Find the median of the triangle from the point C.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_median_c
+  #   => Line.new(Point.new(12, 11, 9), Vector.new(-9.5, -7, -3))
+  #
+  # Arguments:
+  #   no
+
+  def get_median_c
+    Line.construct_with_two_points(point_c, point_a.get_midpoint(point_b))
+  end
+
+  # Find the centroid of the triangle.
+  #
+  # Example:
+  #   >> Triangle.new(Point.new(1, 2, 3), Point.new(4, 6, 9), Point.new(12, 11, 9)).get_centroid
+  #   => Point.new(5.666666666, 6.333333333, 7)
+  #
+  # Arguments:
+  #   no
+
+  def get_centroid
+    get_median_a.find_intersecting_point(get_median_b)
+  end
+end

data/lib/geometry_3d/vector.rb

@@ -0,0 +1,76 @@
+class Vector
+  attr_accessor :x
+  attr_accessor :y
+  attr_accessor :z
+
+  def initialize(x, y, z)
+    @x = x
+    @y = y
+    @z = z
+  end
+
+  # Construct a vector with two points.
+  #
+  # Example:
+  #   >> Vector.construct_with_two_points(Point.new(1, 2, 3), Point.new(3, 2, 1))
+  #   => Vector.new(2, 0, -2)
+  #
+  # Arguments:
+  #   first_point: (Point)
+  #   second_point: (Point)
+
+  def self.construct_with_two_points(first_point, second_point)
+    Vector.new(*second_point.to_a.zip(first_point.to_a).map { |a, b| a - b} )
+  end
+
+  # Construct a vector as normal vector to a plane.
+  #
+  # Example:
+  #   >> Vector.construct_with_plane(Plane.new(-2, 5, 2.5, 7))
+  #   => Vector.new(-2, 5, 2.5)
+  #
+  # Arguments:
+  #   first_point: (Point)
+  #   second_point: (Point)
+
+  def self.construct_with_plane(plane)
+    Vector.new(*plane.to_a.take(3))
+  end
+
+  # Finds the sum vector of two vectors.
+  #
+  # Example:
+  #   >> Vector.new(1, 2, 3) + (Vector.new(-2, 1, 2))
+  #   => Vector.new(-1, 3, 5)
+  #
+  # Arguments:
+  #   vector: (Vector)
+
+  def +(vector)
+    Vector.new(*to_a.zip(vector.to_a).map { |a, b| a + b } )
+  end
+
+  # Finds a vector perpendicular to both vectors.
+  #
+  # Example:
+  #   >> Vector.new(1, 2, 3).cross_product(Vector.new(-2, 1, 2))
+  #   => Vector.new(1, -8, 5)
+  #
+  # Arguments:
+  #   vector: (Vector)
+
+  def cross_product(vector)
+    coordinate_x = y * vector.z - z * vector.y
+    coordinate_y = z * vector.x - x * vector.z
+    coordinate_z = x * vector.y - y * vector.x
+    Vector.new(coordinate_x, coordinate_y, coordinate_z)
+  end
+
+  def to_a()
+    [x, y, z]
+  end
+
+  def ==(vector)
+    x == vector.x and y == vector.y and z == vector.z
+  end
+end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: geometry_3d
 version: !ruby/object:Gem::Version
-  version: 0.0.0
+  version: 0.1.1
 platform: ruby
 authors:
 - Guldjan Kupen
@@ -12,11 +12,20 @@ date: 2015-01-26 00:00:00.000000000 Z
 dependencies: []
 description: A gem implementing basic 3D algorithms in Ruby.
 email: guldjan.kupen@gmail.com
-executables: []
+executables:
+- geometry_3d
 extensions: []
 extra_rdoc_files: []
 files:
+- bin/geometry_3d
 - lib/geometry_3d.rb
+- lib/geometry_3d/circle.rb
+- lib/geometry_3d/equation.rb
+- lib/geometry_3d/line.rb
+- lib/geometry_3d/plane.rb
+- lib/geometry_3d/point.rb
+- lib/geometry_3d/triangle.rb
+- lib/geometry_3d/vector.rb
 homepage: http://rubygems.org/gems/geometry_3d
 licenses:
 - MIT