eymiha_math3 0.1.0

data/test/tc_point3c.rb

@@ -0,0 +1,98 @@
+require 'test/unit'
+
+require 'point3c'
+
+class TC_Point3c < Test::Unit::TestCase
+  
+  understands ThreeDimensions
+  
+  def test_initialize_and_equality
+    p3c = Point3c.new
+    assert((p3c.c_radius == 0)&&(p3c.theta == 0)&&(p3c.z == 0))
+    p3c = Point3c.new 1,2,3
+    assert((p3c.c_radius == 1)&&(p3c.theta == 2)&&(p3c.z == 3))
+    p3a = Point3c.new 1,2,3
+    assert(p3c == p3a)
+    assert(p3c =~ p3a)
+    assert(p3c.approximately_equals?((point3c 1.005, 1.991, 3),0.01))
+    p3a = Point3c.new 3,2,1
+    assert(p3c != p3a)
+  end
+  
+  def test_Object_point3c
+    p3c = point3c
+    assert(p3c == origin)
+    p3c = point3c 1,2,3
+    assert((p3c.c_radius == 1)&&(p3c.theta == 2)&&(p3c.z == 3))
+    p3a = Point3c.new 1,2,3
+    assert((p3a.c_radius == 1)&&(p3a.theta == 2)&&(p3a.z == 3))
+    assert(p3c == p3a)
+  end
+  
+  def test_accessors
+    p3c = point3c
+    assert(p3c.point3c_like?)
+    assert((p3c.c_radius == 0)&&(p3c.theta == 0)&&(p3c.z == 0))
+    p3c.c_radius = 1
+    p3c.theta = 2
+    p3c.z = 3
+    assert((p3c.c_radius == 1)&&(p3c.theta == 2)&&(p3c.z == 3))
+    p3c.c_radius, p3c.theta, p3c.z =  2, 4, 6
+    assert((p3c.c_radius == 2)&&(p3c.theta == 4)&&(p3c.z == 6))
+  end
+  
+   def test_origin
+     p3c = point3c 0, 0, 0
+     assert(p3c == origin)
+     assert(origin == p3c)
+     p3c = point3c 0, 1, 0
+     assert(p3c == origin)
+     assert(origin == p3c)
+  end
+  
+  def test_set
+    p3c = point3c
+    p3c.set 1, 2, 3
+    assert((p3c.c_radius == 1)&&(p3c.theta == 2)&&(p3c.z == 3))
+    p3a = point3c p3c
+    assert((p3a.c_radius == 1)&&(p3a.theta == 2)&&(p3a.z == 3))
+    p3b = point3c
+    assert_raise(TypeError) { p3b.set "bad" }
+    assert_raise(TypeError) { point3c "bad" }
+  end
+
+  def test_copy
+    p3c = point3c 1, 2, 3
+    p3a = p3c.point3c
+    assert(p3c == p3a)
+  end
+  
+  def test_cartesian_and_spherical
+    p3c = point3c 1, 2, 3
+    assert(p3c.x =~ -0.416146836547142)
+    assert(p3c.y =~ 0.909297426825682)
+    assert(p3c.s_radius =~ 3.16227766016838)
+    assert(p3c.phi =~ 0.321750554396642)
+  end
+ 
+  def test_point3
+    p3c = point3c 1, 2, 3
+    p3 = p3c.point3
+    assert(p3.x =~ -0.416146836547142)
+    assert(p3.y =~ 0.909297426825682)
+    assert(p3.z =~ 3)
+    p3a = p3.point3c
+    assert(p3c =~ p3a)
+  end
+
+  def test_point3s
+    p3c = point3c 1, 2, 3
+    p3s = p3c.point3s
+    assert(p3s.s_radius =~ 3.16227766016838)
+    assert(p3s.theta =~ p3c.theta)
+    assert(p3s.phi =~ 0.321750554396642)
+    p3a = p3s.point3c
+    assert(p3c =~ p3a)
+  end
+
+end

data/test/tc_point3s.rb

@@ -0,0 +1,101 @@
+require 'test/unit'
+
+require 'point3s'
+
+class TC_Point3s < Test::Unit::TestCase
+  
+  understands ThreeDimensions
+  
+  def test_initialize_and_equality
+    p3s = Point3s.new
+    assert((p3s.s_radius == 0)&&(p3s.theta == 0)&&(p3s.phi == 0))
+    p3s = Point3s.new 1,2,3
+    assert((p3s.s_radius == 1)&&(p3s.theta == 2)&&(p3s.phi == 3))
+    p3a = Point3s.new 1,2,3
+    assert(p3s == p3a)
+    assert(p3s =~ p3a)
+    assert(p3s.approximately_equals?((point3s 1.005, 1.991, 3),0.01))
+    p3a = Point3s.new 3,2,1
+    assert(p3s != p3a)
+  end
+  
+  def test_Object_point3s
+    p3s = point3s
+    assert(p3s == origin)
+    p3s = point3s 1,2,3
+    assert((p3s.s_radius == 1)&&(p3s.theta == 2)&&(p3s.phi == 3))
+    p3a = Point3s.new 1,2,3
+    assert((p3a.s_radius == 1)&&(p3a.theta == 2)&&(p3a.phi == 3))
+    assert(p3s == p3a)
+  end
+  
+  def test_accessors
+    p3s = point3s
+    assert(p3s.point3s_like?)
+    assert((p3s.s_radius == 0)&&(p3s.theta == 0)&&(p3s.phi == 0))
+    p3s.s_radius = 1
+    p3s.theta = 2
+    p3s.phi = 3
+    assert((p3s.s_radius == 1)&&(p3s.theta == 2)&&(p3s.phi == 3))
+    p3s.s_radius, p3s.theta, p3s.phi =  2, 4, 6
+    assert((p3s.s_radius == 2)&&(p3s.theta == 4)&&(p3s.phi == pi))
+    p3s.s_radius, p3s.theta, p3s.phi =  -2, -4, -6
+    assert((p3s.s_radius == -2)&&(p3s.theta =~ 2.28318530717959)&&
+           (p3s.phi == 0))
+  end
+  
+   def test_origin
+     p3s = point3s 0, 0, 0
+     assert(p3s == origin)
+     assert(origin == p3s)
+     p3s = point3s 0, 1, 2
+     assert(p3s == origin)
+     assert(origin == p3s)
+  end
+  
+  def test_set
+    p3s = point3s
+    p3s.set 1, 2, 3
+    assert((p3s.s_radius == 1)&&(p3s.theta == 2)&&(p3s.phi == 3))
+    p3a = point3s p3s
+    assert((p3a.s_radius == 1)&&(p3a.theta == 2)&&(p3a.phi == 3))
+    p3b = point3s
+    assert_raise(TypeError) { p3b.set "bad" }
+    assert_raise(TypeError) { point3s "bad" }
+  end
+
+  def test_copy
+    p3s = point3s 1, 2, 3
+    p3a = p3s.point3s
+    assert(p3s == p3a)
+  end
+  
+  def test_cartesian_and_cylindrical
+    p3s = point3s 1, 2, 3
+    assert(p3s.x =~ -0.058726644927621)
+    assert(p3s.y =~ 0.128320060202457)
+    assert(p3s.z =~ -0.989992496600445)
+    assert(p3s.c_radius =~ 0.141120008059867)
+  end
+ 
+  def test_point3
+    p3s = point3s 1, 2, 3
+    p3 = p3s.point3
+    assert(p3.x =~ -0.058726644927621)
+    assert(p3.y =~ 0.128320060202457)
+    assert(p3.z =~ -0.989992496600445)
+    p3a = p3.point3s
+    assert(p3s =~ p3a)
+  end
+
+  def test_point3c
+    p3s = point3s 1, 2, 3
+    p3c = p3s.point3c
+    assert(p3c.c_radius =~ 0.141120008059867)
+    assert(p3c.theta =~ p3s.theta)
+    assert(p3c.z =~ -0.989992496600445)
+    p3a = p3c.point3s
+    assert(p3s =~ p3a)
+  end
+
+end

data/test/tc_quaternion.rb

@@ -0,0 +1,148 @@
+require 'test/unit'
+
+require 'quaternion'
+
+class TC_Quaternion < Test::Unit::TestCase
+  
+  understands ThreeDimensions
+  
+  def test_initialize_and_equality
+    q = Quaternion.new
+    assert((q.axis == origin)&&(q.real == 0))
+    qa = Quaternion.new(origin,0)
+    assert(q == qa)
+    p3 = point3 1,2,3
+    qb = Quaternion.new p3,4
+    assert(qb.axis == p3)
+    assert(qb.real == 4)
+    assert(q != qb)
+    qc = quaternion p3,4
+    assert(qc.axis == p3)
+    assert(qc.real == 4)
+    assert(q != qb)
+    assert(qb == qc)
+  end
+
+  def test_accessors
+    p3 = point3 1,2,3
+    q = quaternion p3,4
+    q.axis.scale! 2
+    q.real -= 3
+    assert(q.axis == point3(2,4,6))
+    assert(q.real == 1)
+    qb = quaternion q
+    assert(q == qb)
+  end
+
+  def test_axis_angle
+    p3 = point3 1,2,3
+    q = Quaternion.from_axis_angle(p3,pi/8)
+    assert(q.axis.x =~ 0.195090322016128)
+    assert(q.axis.y =~ 0.390180644032256)
+    assert(q.axis.z =~ 0.585270966048385)
+    assert(q.real =~ 0.98078528040323)
+    qb = q.to_axis_angle
+    assert(qb.axis =~ p3)
+    assert(qb.real =~ pi/8)
+  end
+
+  def test_add
+    p3a,ra = point3(1,2,3),4
+    p3b,rb = point3(4,-6,0),-3
+    qa = quaternion p3a,ra
+    qb = quaternion p3b,rb
+    p3c,rc = p3a+p3b,ra+rb
+    qc = qa+qb
+    assert(qc.axis =~ p3c)
+    assert(qc.real =~ rc)
+    qa.add! qb
+    assert(qa =~ qc)
+  end
+
+  def test_subtract
+    p3a,ra = point3(1,2,3),4
+    p3b,rb = point3(4,-6,0),-3
+    qa = quaternion p3a,ra
+    qb = quaternion p3b,rb
+    p3c,rc = p3a-p3b,ra-rb
+    qc = qa-qb
+    assert(qc.axis =~ p3c)
+    assert(qc.real =~ rc)
+    qa.subtract! qb
+    assert(qa =~ qc)
+  end
+
+  def test_multiply
+    p3a,ra = point3(1,2,3),4
+    p3b,rb = point3(4,-6,0),-3
+    qa = quaternion p3a,ra
+    qb = quaternion p3b,rb
+    p3c,rc = point3(-1,30,-23),-4
+    qc = qa*qb
+    assert(qc.axis =~ p3c)
+    assert(qc.real =~ rc)
+    qa.multiply! qb
+    assert(qa =~ qc)
+  end
+
+  def test_scale
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    p3c,rc = p3a.scale(5.5),ra*5.5
+    qc = qa.scale(5.5)
+    assert(qc.axis =~ p3c)
+    assert(qc.real =~ rc)
+    qa.scale! 5.5
+    assert(qa =~ qc)
+  end
+
+  def test_conjugate
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    p3c,rc = point3(-1,-2,-3),ra
+    qc = qa.conjugate
+    assert(qc.axis =~ p3c)
+    assert(qc.real =~ rc)
+    qa.conjugate!
+    assert(qa =~ qc)
+  end
+
+  def test_norm
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    assert(qa.norm =~ 19.7416573867739)
+  end
+
+  def test_abs
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    assert(qa.abs =~ 4.44315849219606)
+  end
+
+  def test_inverse
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    q = qa.inverse
+    assert(q.axis.x =~ (-0.0506543083191159))
+    assert(q.axis.y =~ (-0.101308616638232))
+    assert(q.axis.z =~ (-0.151962924957348))
+    assert(q.real =~ 0.202617233276464)
+    qa.inverse!
+    assert(q =~ qa)
+  end
+
+  def test_divide
+    p3a,ra = point3(1,2,3),4
+    qa = quaternion p3a,ra
+    p3b,rb = point3(10,0,-10),0
+    qb = quaternion p3b,rb
+    qc = qa.divide qb
+    assert(qc.axis.x =~ 2.41421356237309)
+    assert(qc.axis.y =~ (-2.82842712474619))
+    assert(qc.axis.z =~ 0.414213562373095)
+    assert(qc.real =~ -1.41421356237309)
+    qa.divide! qb
+    assert(qa =~ qc)
+  end
+
+end

metadata

@@ -0,0 +1,75 @@
+--- !ruby/object:Gem::Specification 
+rubygems_version: 0.9.0
+specification_version: 1
+name: eymiha_math3
+version: !ruby/object:Gem::Version 
+  version: 0.1.0
+date: 2007-05-22 00:00:00 -04:00
+summary: Emiyha - basic 3D math extensions
+require_paths: 
+- lib
+email: dave@eymiha.com
+homepage: http://www.eymiha.com
+rubyforge_project: cori
+description: 
+autorequire: eymiha_math3
+default_executable: 
+bindir: bin
+has_rdoc: true
+required_ruby_version: !ruby/object:Gem::Version::Requirement 
+  requirements: 
+  - - ">"
+    - !ruby/object:Gem::Version 
+      version: 0.0.0
+  version: 
+platform: ruby
+signing_key: 
+cert_chain: 
+post_install_message: 
+authors: 
+- Dave Anderson
+files: 
+- gem_package.rb
+- rakefile.rb
+- lib/envelope3.rb
+- lib/eymiha_math3.rb
+- lib/point3.rb
+- lib/point3c.rb
+- lib/point3s.rb
+- lib/quaternion.rb
+- test/tc_envelope3.rb
+- test/tc_point3.rb
+- test/tc_point3c.rb
+- test/tc_point3s.rb
+- test/tc_quaternion.rb
+test_files: []
+
+rdoc_options: 
+- --all
+extra_rdoc_files: []
+
+executables: []
+
+extensions: []
+
+requirements: []
+
+dependencies: 
+- !ruby/object:Gem::Dependency 
+  name: eymiha
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Version::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: 0.1.0
+    version: 
+- !ruby/object:Gem::Dependency 
+  name: eymiha_math
+  version_requirement: 
+  version_requirements: !ruby/object:Gem::Version::Requirement 
+    requirements: 
+    - - ">="
+      - !ruby/object:Gem::Version 
+        version: 0.1.0
+    version: