mageo 0.0.1 → 0.0.2

data/lib/mageo/vector.rb

@@ -1,18 +1,16 @@
 #! /usr/bin/ruby
 
-require "matrix"
-
 class Vector
   #include Math
-  class ZeroOperation < Exception; end
+  class ZeroOperationError < Exception; end
   class SizeError < Exception; end
 
   # Get a unit vector.
   def unit_vector
     len = self.r
-    raise Vector::ZeroOperation if (len == 0)
+    raise Vector::ZeroOperationError if (len == 0)
     self * (1/len)
-    # Vector3D.new(@x*(1.0/len), @y*(1.0/len), @z*(1.0/len))
+    # Mageo::Vector3D.new(@x*(1.0/len), @y*(1.0/len), @z*(1.0/len))
   end
 
   def floor

data/lib/mageo/vector3d.rb

@@ -1,12 +1,5 @@
 #! /usr/bin/ruby
 
-# written by Ippei KISHIDA
-#
-require "matrix"
-require "mageo/polar2d.rb"
-require "mageo/polar3d.rb"
-require "mageo/axes.rb"
-require "mageo/vector.rb"
 require "rubygems"
 gem "malge"
 require "malge.rb"
@@ -15,18 +8,18 @@ require "malge.rb"
 
 # Open class to add "to_v3d" method.
 class Array
-  # Convert Array to Vector3D
+  # Convert Array to Mageo::Vector3D
   def to_v3d
-    Vector3D[*self]
-    #要素数チェックは Vector3D.[] 側でやっている。
+    Mageo::Vector3D[*self]
+    #要素数チェックは Mageo::Vector3D.[] 側でやっている。
   end
 end
 
 class Vector
-  # Return a new instance converted to Vector3D class.
+  # Return a new instance converted to Mageo::Vector3D class.
   def to_v3d
-    Vector3D[*self]
-    #要素数チェックは Vector3D.[] 側でやっている。
+    Mageo::Vector3D[*self]
+    #要素数チェックは Mageo::Vector3D.[] 側でやっている。
   end
 end
 
@@ -35,16 +28,16 @@ end
 # This class provide exterior_product method and others, which is not included
 # in native Vector class.
 # This class is constructed under the assumption in the Cartesian coordinate.
-# If you want to be in an internal coordinate, you can use Math/Vector3DInternal.rb .
+# If you want to be in an internal coordinate, you can use Math/Mageo::Vector3DInternal.rb .
 # 
 # Memo:
-#   Vector3DInternal との対比として、Vector3DCartesian という名前にすることも考えたが、
+#   Mageo::Vector3DInternal との対比として、Vector3DCartesian という名前にすることも考えたが、
 #   長くなるし、普通直交座標で考えるよね、と。
 # 
 #   インスタンス生成の時点で要素数をチェックし、要素の追加削除を禁止しているので
 #   要素数は常に3であることが保証されている。
 # 
-class Vector3D < Vector
+class Mageo::Vector3D < Vector
 
   class TypeError < Exception; end
   class ZeroOperation < Exception; end
@@ -60,7 +53,7 @@ class Vector3D < Vector
   # Get the exterior product.
   def self.exterior_product(vec0, vec1)
     [vec0, vec1].each_with_index do |vec, index|
-      unless (vec.class == Vector3D)
+      unless (vec.class == Mageo::Vector3D)
         raise TypeError, "Vector #{index}, #{vec.inspect}."
       end
     end
@@ -83,7 +76,7 @@ class Vector3D < Vector
   # Get the scalar triple product.
   def self.scalar_triple_product(vec0, vec1, vec2)
     [vec0, vec1, vec2].each_with_index do |vec, index|
-      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Vector3D)
+      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Mageo::Vector3D)
     end
 
     vec0.inner_product(vec1.exterior_product(vec2))
@@ -92,7 +85,7 @@ class Vector3D < Vector
   # Get the angle with radian between self and other vectors.
   def self.angle_radian(vec0, vec1)
     [vec0, vec1].each_with_index do |vec, index|
-      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Vector3D)
+      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Mageo::Vector3D)
       raise ZeroOperation, "#{index}th vector: #{vec.inspect}" if (vec.r == 0.0)
     end
 
@@ -102,7 +95,7 @@ class Vector3D < Vector
   # Get the angle with degree between self and other vectors.
   def self.angle_degree(vec0, vec1)
     [vec0, vec1].each_with_index do |vec, index|
-      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Vector3D)
+      raise TypeError, "#{index}th vector: #{vec.inspect}" unless (vec.class == Mageo::Vector3D)
       raise ZeroOperation, "#{index}th vector: #{vec.inspect}" if (vec.r == 0.0)
     end
 
@@ -123,33 +116,33 @@ class Vector3D < Vector
   
   
   # ベクトルが等しいかチェック。
-  # other として Vector3D クラス以外のインスタンス渡すと Vector3D::TypeError。
+  # other として Mageo::Vector3D クラス以外のインスタンス渡すと Vector3D::TypeError。
   # 両者の差分ベクトルの長さが tolerance 以下という判定になる。
   def equal_in_delta?(other, tolerance = 0.0)
-    raise TypeError if (other.class != Vector3D)
+    raise TypeError if (other.class != Mageo::Vector3D)
     return (other - self).r <= tolerance
   end
 
   # Vectorクラスで用意されているメソッドは Vectorクラスインスタンスを返すようになっているので、
-  # Vector3D クラスインスタンスを返すようにした + メソッド。
+  # Mageo::Vector3D クラスインスタンスを返すようにした + メソッド。
   def +(vec)
-    unless (vec.class == Vector3D)
+    unless (vec.class == Mageo::Vector3D)
       raise TypeError, "#{vec.inspect}."
     end
     super(vec).to_v3d
   end
 
   # Vectorクラスで用意されているメソッドは Vectorクラスインスタンスを返すようになっているので、
-  # Vector3D クラスインスタンスを返すようにした - メソッド。
+  # Mageo::Vector3D クラスインスタンスを返すようにした - メソッド。
   def -(vec)
-    unless (vec.class == Vector3D)
+    unless (vec.class == Mageo::Vector3D)
       raise TypeError, "#{vec.inspect}."
     end
     super(vec).to_v3d
   end
 
   # Vectorクラスで用意されているメソッドは Vectorクラスインスタンスを返すようになっているので、
-  # Vector3D クラスインスタンスを返すようにした * メソッド。
+  # Mageo::Vector3D クラスインスタンスを返すようにした * メソッド。
   # Argument 'val' must have :to_f method.
   def *(val)
     #raise TypeError if (val.class != Float)
@@ -158,18 +151,18 @@ class Vector3D < Vector
   end
 
   # Vectorクラスで用意されているメソッドは Vectorクラスインスタンスを返すようになっているので、
-  # Vector3D クラスインスタンスを返すようにした clone メソッド。
+  # Mageo::Vector3D クラスインスタンスを返すようにした clone メソッド。
   def clone
     super().to_v3d
   end
 
-  # Convert to Vector3DInternal. Non-destructive.
+  # Convert to Mageo::Vector3DInternal. Non-destructive.
   def to_v3di(axes)
-    #pp axes.is_a?(Axes)
-    raise TypeError unless axes.is_a?(Axes)
+    #pp axes.is_a?(Mageo::Axes)
+    raise TypeError unless axes.is_a?(Mageo::Axes)
 
     axes = axes.to_a
-    Vector3DInternal[ *(Malge::SimultaneousEquations.cramer(axes.transpose, self)) ]
+    Mageo::Vector3DInternal[ *(Malge::SimultaneousEquations.cramer(axes.transpose, self)) ]
   end
 
   #Return size, always 3.
@@ -196,17 +189,17 @@ class Vector3D < Vector
     self.class.angle_degree(self, vec)
   end
 
-  #3次元極座標への変換した Polar3D インスタンスを返す。
+  #3次元極座標への変換した Mageo::Polar3D インスタンスを返す。
   def to_p3d
     r = self.r
     if r == 0.0
       theta = 0.0
       phi = 0.0
     else
-      theta = Polar2D.minimum_radian(Math::acos(self[2] / r))
+      theta = Mageo::Polar2D.minimum_radian(Math::acos(self[2] / r))
       phi = Vector[ self[0], self[1] ].to_p2d.theta
     end
-    Polar3D.new(r, theta, phi)
+    Mageo::Polar3D.new(r, theta, phi)
   end
 
   #x, y, z 軸のいずれかで self を回転する。破壊的。

data/lib/mageo/vector3dinternal.rb

@@ -1,8 +1,8 @@
-require 'mageo/vector3d.rb'
+#require 'mageo/vector3d.rb'
 
 class Array
   def to_v3di
-    Vector3DInternal[ *self ]
+    Mageo::Vector3DInternal[ *self ]
   end
 end
 
@@ -11,16 +11,16 @@ end
 # (直交座標と間違えて)内部座標で内積を求めるなど誤った演算を例外で止めるのが目的。
 # 座標軸の情報は自身では持たない。
 # Actually, this class is very similar to Array class, except for +, -, and * methods.
-# Vector, Vector3D との混在計算は例外を発生。
+# Vector, Mageo::Vector3D との混在計算は例外を発生。
 #
 # 軸の情報は持たない。
 # 軸に関係ない抽象的な内部座標について議論することもありうるし、
 # 軸情報が必要なのは to_v3d メソッドくらいなので。
-class Vector3DInternal < Vector3D
+class Mageo::Vector3DInternal < Mageo::Vector3D
   class RangeError < Exception; end
   class TypeError < Exception; end
 
-  #要素数3以外では例外 Vector3DInternalSizeError を発生。
+  #要素数3以外では例外 Mageo::Vector3DInternalSizeError を発生。
   def self.[]( *args )
     raise RangeError unless args.size == 3
     super *args
@@ -36,7 +36,7 @@ class Vector3DInternal < Vector3D
     return self
   end
 
-  # Convert to Vector3D. Non-destructive.
+  # Convert to Mageo::Vector3D. Non-destructive.
   def to_v3d( axes )
     result = [ 0.0, 0.0, 0.0 ]
     3.times do |i|
@@ -44,7 +44,7 @@ class Vector3DInternal < Vector3D
         result[i] += ( self[j] * axes[j][i] )
       end
     end
-    Vector3D[ *result ]
+    Mageo::Vector3D[ *result ]
   end
 
   ## Return an array converted from self.
@@ -62,13 +62,13 @@ class Vector3DInternal < Vector3D
   # return true
   #end
 
-  #0〜2 以外の要素にアクセスしようとすると例外 Vector3DInternal::RangeError を発生。
+  #0〜2 以外の要素にアクセスしようとすると例外 Mageo::Vector3DInternal::RangeError を発生。
   def []( index )
     raise RangeError if ( index < 0 || 2 < index )
     super index
   end
 
-  #0〜2 以外の要素にアクセスしようとすると例外 Vector3DInternal::RangeError を発生。
+  #0〜2 以外の要素にアクセスしようとすると例外 Mageo::Vector3DInternal::RangeError を発生。
   def []=( index, val )
     raise RangeError if ( index < 0 || 2 < index )
     super index, val
@@ -76,11 +76,11 @@ class Vector3DInternal < Vector3D
 
   #ベクトルとしての加算
   def +(vec)
-    unless vec.class == Vector3DInternal
+    unless vec.class == Mageo::Vector3DInternal
       raise TypeError,
         "Argument: #{vec.inspect}, #{vec.class}."
     end
-    result = Vector3DInternal[0.0, 0.0, 0.0]
+    result = Mageo::Vector3DInternal[0.0, 0.0, 0.0]
     3.times do |i|
       result[ i ] = (self[ i ] + vec[ i ])
     end
@@ -89,11 +89,11 @@ class Vector3DInternal < Vector3D
 
   #ベクトルとしての減算
   def -( vec )
-    unless vec.class == Vector3DInternal
+    unless vec.class == Mageo::Vector3DInternal
       raise TypeError,
         "Argument: #{vec.inspect}, #{vec.class}."
     end
-    result = Vector3DInternal[0.0, 0.0, 0.0]
+    result = Mageo::Vector3DInternal[0.0, 0.0, 0.0]
     3.times do |i|
       result[ i ] = (self[ i ] - vec[ i ])
     end
@@ -102,7 +102,7 @@ class Vector3DInternal < Vector3D
 
   #ベクトルとしての乗算
   def *( val )
-    result = Vector3DInternal[0.0, 0.0, 0.0]
+    result = Mageo::Vector3DInternal[0.0, 0.0, 0.0]
     3.times do |i|
       result[ i ] = ( self[ i ] * val )
     end

data/mageo.gemspec

@@ -5,12 +5,12 @@
 
 Gem::Specification.new do |s|
   s.name = "mageo"
-  s.version = "0.0.1"
+  s.version = "0.0.2"
 
   s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
   s.authors = ["ippei94da"]
-  s.date = "2013-02-02"
-  s.description = "MAth GEOmetry library to deal with 2 and 3 dimension space.\n    Cartesian and internal coordinate systems can be used.\n    This includes besic objects in 3 dimensional space.\n  "
+  s.date = "2013-04-18"
+  s.description = "MAth GEOmetry library to deal with 2 and 3 dimensional spaces.\n    Cartesian and internal coordinate systems can be used.\n    This includes besic objects in 3 dimensional space.\n  "
   s.email = "ippei94da@gmail.com"
   s.extra_rdoc_files = [
     "LICENSE.txt",
@@ -58,34 +58,31 @@ Gem::Specification.new do |s|
   s.licenses = ["MIT"]
   s.require_paths = ["lib"]
   s.rubygems_version = "1.8.11"
-  s.summary = "MAth GEOmetry library to deal with 2 and 3 dimension space."
+  s.summary = "MAth GEOmetry library to deal with 2 and 3 dimensional spaces."
 
   if s.respond_to? :specification_version then
     s.specification_version = 3
 
     if Gem::Version.new(Gem::VERSION) >= Gem::Version.new('1.2.0') then
-      s.add_development_dependency(%q<builtinextension>, [">= 0"])
-      s.add_development_dependency(%q<bundler>, ["~> 1.2.2"])
+      s.add_development_dependency(%q<builtinextension>, [">= 0.1.0"])
+      s.add_development_dependency(%q<bundler>, ["~> 1.3.5"])
       s.add_development_dependency(%q<jeweler>, ["~> 1.8.3"])
-      s.add_development_dependency(%q<malge>, [">= 0.0.2"])
-      s.add_development_dependency(%q<psych>, [">= 0"])
+      s.add_development_dependency(%q<malge>, [">= 0.0.8"])
       s.add_development_dependency(%q<rdoc>, ["~> 3.12"])
       s.add_development_dependency(%q<simplecov>, [">= 0"])
     else
-      s.add_dependency(%q<builtinextension>, [">= 0"])
-      s.add_dependency(%q<bundler>, ["~> 1.2.2"])
+      s.add_dependency(%q<builtinextension>, [">= 0.1.0"])
+      s.add_dependency(%q<bundler>, ["~> 1.3.5"])
       s.add_dependency(%q<jeweler>, ["~> 1.8.3"])
-      s.add_dependency(%q<malge>, [">= 0.0.2"])
-      s.add_dependency(%q<psych>, [">= 0"])
+      s.add_dependency(%q<malge>, [">= 0.0.8"])
       s.add_dependency(%q<rdoc>, ["~> 3.12"])
       s.add_dependency(%q<simplecov>, [">= 0"])
     end
   else
-    s.add_dependency(%q<builtinextension>, [">= 0"])
-    s.add_dependency(%q<bundler>, ["~> 1.2.2"])
+    s.add_dependency(%q<builtinextension>, [">= 0.1.0"])
+    s.add_dependency(%q<bundler>, ["~> 1.3.5"])
     s.add_dependency(%q<jeweler>, ["~> 1.8.3"])
-    s.add_dependency(%q<malge>, [">= 0.0.2"])
-    s.add_dependency(%q<psych>, [">= 0"])
+    s.add_dependency(%q<malge>, [">= 0.0.8"])
     s.add_dependency(%q<rdoc>, ["~> 3.12"])
     s.add_dependency(%q<simplecov>, [">= 0"])
   end

data/test/test_axes.rb

@@ -2,17 +2,17 @@
 # coding: utf-8
 
 require "helper"
-require "test/unit"
-require "mageo/axes.rb"
+#require "test/unit"
+#require "mageo/axes.rb"
 
 class TC_Axes < Test::Unit::TestCase
   $tolerance = 1.0 * 10.0**(-10)
 
   def setup
-    @a10 = Axes.new([[1.0]])
-    @a20 = Axes.new([[1.0, 0.0], [0.0, 1.0]])
-    @a30 = Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]])
-    @a40 = Axes.new([
+    @a10 = Mageo::Axes.new([[1.0]])
+    @a20 = Mageo::Axes.new([[1.0, 0.0], [0.0, 1.0]])
+    @a30 = Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]])
+    @a40 = Mageo::Axes.new([
         [1.0, 0.0, 0.0, 0.0],
         [0.0, 1.0, 0.0, 0.0],
         [0.0, 0.0, 1.0, 0.0],
@@ -20,12 +20,12 @@ class TC_Axes < Test::Unit::TestCase
       ]
     )
 
-    @a31 = Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] )
-    @a32 = Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.5], [0.0, 0.5, 0.5] ] )
-    @a33 = Axes.new( [ [1.0, 0.0, 0.0], [1.0, 0.0, 0.0], [0.0, 0.0, 1.0] ] )
-    @a34 = Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.0], [0.0, 0.5, 0.0] ] )
-    @a35 = Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] )
-    @a36 = Axes.new( [ [1.0, 1.0, 1.0], [0.0,-1.0,-1.0], [0.0, 0.0, 1.0] ] )
+    @a31 = Mageo::Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] )
+    @a32 = Mageo::Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.5], [0.0, 0.5, 0.5] ] )
+    @a33 = Mageo::Axes.new( [ [1.0, 0.0, 0.0], [1.0, 0.0, 0.0], [0.0, 0.0, 1.0] ] )
+    @a34 = Mageo::Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.0], [0.0, 0.5, 0.0] ] )
+    @a35 = Mageo::Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] )
+    @a36 = Mageo::Axes.new( [ [1.0, 1.0, 1.0], [0.0,-1.0,-1.0], [0.0, 0.0, 1.0] ] )
 
     @vec_x = Vector[ 1.0, 0.0, 0.0 ]
     @vec_y = Vector[ 0.0, 1.0, 0.0 ]
@@ -35,56 +35,56 @@ class TC_Axes < Test::Unit::TestCase
   end
 
   def test_initialize
-    assert_raise(Axes::InitializeError) { Axes.new([[]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [1.0, 0.0, 0.0], [0.0, 1.0, 1.0]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0], [0.0, 0.0, 1.0]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0]]) }
-    assert_raise(Axes::InitializeError) { Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0, 0.0]]) }
-
-    assert_nothing_raised{ Axes.new([["1.0", "0.0", "0.0"], [0.0, 1.0, 0.0], [0.0, 0.0, 0.0]]) }
-    assert_nothing_raised{ Axes.new([[1, 0, 0], [0, 1, 0], [0, 0, 0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [1.0, 0.0, 0.0], [0.0, 1.0, 1.0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0], [0.0, 0.0, 1.0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0]]) }
+    assert_raise(Mageo::Axes::InitializeError) { Mageo::Axes.new([[1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0, 0.0]]) }
+
+    assert_nothing_raised{ Mageo::Axes.new([["1.0", "0.0", "0.0"], [0.0, 1.0, 0.0], [0.0, 0.0, 0.0]]) }
+    assert_nothing_raised{ Mageo::Axes.new([[1, 0, 0], [0, 1, 0], [0, 0, 0]]) }
   end
 
   def test_self_dependent?
     # Righthand
-    assert_equal( false, Axes.dependent?( [ @vec_x, @vec_y, @vec_z ] ) )
-    assert_equal( false, Axes.dependent?( [ @vec_y, @vec_z, @vec_x ] ) )
-    assert_equal( false, Axes.dependent?( [ @vec_z, @vec_x, @vec_y ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_x, @vec_y, @vec_z ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_y, @vec_z, @vec_x ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_z, @vec_x, @vec_y ] ) )
     
     # Lefthand
-    assert_equal( false, Axes.dependent?( [ @vec_z, @vec_y, @vec_x ] ) )
-    assert_equal( false, Axes.dependent?( [ @vec_x, @vec_z, @vec_y ] ) )
-    assert_equal( false, Axes.dependent?( [ @vec_y, @vec_x, @vec_z ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_z, @vec_y, @vec_x ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_x, @vec_z, @vec_y ] ) )
+    assert_equal( false, Mageo::Axes.dependent?( [ @vec_y, @vec_x, @vec_z ] ) )
     
     # Including zero vector.
-    assert_equal( true , Axes.dependent?( [ @vec_0, @vec_y, @vec_z ] ) )
-    assert_equal( true , Axes.dependent?( [ @vec_0, @vec_z, @vec_x ] ) )
-    assert_equal( true , Axes.dependent?( [ @vec_0, @vec_x, @vec_y ] ) )
+    assert_equal( true , Mageo::Axes.dependent?( [ @vec_0, @vec_y, @vec_z ] ) )
+    assert_equal( true , Mageo::Axes.dependent?( [ @vec_0, @vec_z, @vec_x ] ) )
+    assert_equal( true , Mageo::Axes.dependent?( [ @vec_0, @vec_x, @vec_y ] ) )
 
     # One vector is on the plane of residual two vectors.
-    assert_equal( true , Axes.dependent?( [ @vec_x, @vec_y, @vec_1 ] ) )
+    assert_equal( true , Mageo::Axes.dependent?( [ @vec_x, @vec_y, @vec_1 ] ) )
   end
 
   def test_self_independent?
     # Righthand
-    assert_equal( true , Axes.independent?( [ @vec_x, @vec_y, @vec_z ] ) )
-    assert_equal( true , Axes.independent?( [ @vec_y, @vec_z, @vec_x ] ) )
-    assert_equal( true , Axes.independent?( [ @vec_z, @vec_x, @vec_y ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_x, @vec_y, @vec_z ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_y, @vec_z, @vec_x ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_z, @vec_x, @vec_y ] ) )
     
     # Lefthand
-    assert_equal( true , Axes.independent?( [ @vec_z, @vec_y, @vec_x ] ) )
-    assert_equal( true , Axes.independent?( [ @vec_x, @vec_z, @vec_y ] ) )
-    assert_equal( true , Axes.independent?( [ @vec_y, @vec_x, @vec_z ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_z, @vec_y, @vec_x ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_x, @vec_z, @vec_y ] ) )
+    assert_equal( true , Mageo::Axes.independent?( [ @vec_y, @vec_x, @vec_z ] ) )
     
     # Including zero vector.
-    assert_equal( false, Axes.independent?( [ @vec_0, @vec_y, @vec_z ] ) )
-    assert_equal( false, Axes.independent?( [ @vec_0, @vec_z, @vec_x ] ) )
-    assert_equal( false, Axes.independent?( [ @vec_0, @vec_x, @vec_y ] ) )
+    assert_equal( false, Mageo::Axes.independent?( [ @vec_0, @vec_y, @vec_z ] ) )
+    assert_equal( false, Mageo::Axes.independent?( [ @vec_0, @vec_z, @vec_x ] ) )
+    assert_equal( false, Mageo::Axes.independent?( [ @vec_0, @vec_x, @vec_y ] ) )
 
     # One vector is on the plane of residual two vectors.
-    assert_equal( false, Axes.independent?( [ @vec_x, @vec_y, @vec_1 ] ) )
+    assert_equal( false, Mageo::Axes.independent?( [ @vec_x, @vec_y, @vec_1 ] ) )
   end
 
   def test_size
@@ -95,17 +95,17 @@ class TC_Axes < Test::Unit::TestCase
   end
 
   def test_equal
-    assert_equal( true , @a31 == Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
-    assert_equal( false, @a31 == Axes.new( [ [0.5, 0.5, 0.0], [0.0, 0.5, 0.5], [0.5, 0.0, 0.5] ] ))
-    assert_equal( false, @a31 == Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( true , @a31 == Mageo::Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( false, @a31 == Mageo::Axes.new( [ [0.5, 0.5, 0.0], [0.0, 0.5, 0.5], [0.5, 0.0, 0.5] ] ))
+    assert_equal( false, @a31 == Mageo::Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
 
-    assert_equal( false, @a32 == Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
-    assert_equal( true , @a32 == Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.5], [0.0, 0.5, 0.5] ] ))
-    assert_equal( false, @a32 == Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( false, @a32 == Mageo::Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( true , @a32 == Mageo::Axes.new( [ [0.5, 0.5, 0.0], [0.5, 0.0, 0.5], [0.0, 0.5, 0.5] ] ))
+    assert_equal( false, @a32 == Mageo::Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
 
-    assert_equal( false, @a35 == Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
-    assert_equal( false, @a35 == Axes.new( [ [0.5, 0.5, 0.0], [0.0, 0.5, 0.5], [0.5, 0.0, 0.5] ] ))
-    assert_equal( true , @a35 == Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( false, @a35 == Mageo::Axes.new( [ [1.0, 0.0, 0.0], [0.0, 1.0, 0.0], [0.0, 0.0, 1.0] ] ))
+    assert_equal( false, @a35 == Mageo::Axes.new( [ [0.5, 0.5, 0.0], [0.0, 0.5, 0.5], [0.5, 0.0, 0.5] ] ))
+    assert_equal( true , @a35 == Mageo::Axes.new( [ [1.0, 1.0, 1.0], [0.0, 1.0, 1.0], [0.0, 0.0, 1.0] ] ))
   end
 
   def test_dependent?