Class: Vector

Inherits:
Object
  • Object
show all
Extended by:
Matrix::ConversionHelper
Includes:
Enumerable, ExceptionForMatrix, Matrix::CoercionHelper
Defined in:
lib/matrix.rb

Overview

The Vector class represents a mathematical vector, which is useful in its own right, and also constitutes a row or column of a Matrix.

Method Catalogue

To create a Vector:

  • Vector.[](*array)

  • Vector.elements(array, copy = true)

To access elements:

  • [](i)

To enumerate the elements:

  • #each2(v)

  • #collect2(v)

Vector arithmetic:

  • *(x) "is matrix or number"

  • +(v)

  • -(v)

Vector functions:

  • #inner_product(v)

  • #collect

  • #magnitude

  • #map

  • #map2(v)

  • #norm

  • #normalize

  • #r

  • #size

Conversion to other data types:

  • #covector

  • #to_a

  • #coerce(other)

String representations:

  • #to_s

  • #inspect

Defined Under Namespace

Classes: ZeroVectorError

Class Method Summary collapse

Instance Method Summary collapse

Methods included from Matrix::CoercionHelper

coerce_to, coerce_to_int

Constructor Details

#initialize(array) ⇒ Vector

Vector.new is private; use Vector[] or Vector.elements to create.



1567
1568
1569
1570
 
# File 'lib/matrix.rb', line 1567

def initialize(array)
  # No checking is done at this point.
  @elements = array
end

Class Method Details

.[](*array) ⇒ Object

Creates a Vector from a list of elements.

Vector[7, 4, ...]



1552
1553
1554
 
# File 'lib/matrix.rb', line 1552

def Vector.[](*array)
  new convert_to_array(array, false)
end

.elements(array, copy = true) ⇒ Object

Creates a vector from an Array. The optional second argument specifies whether the array itself or a copy is used internally.



1560
1561
1562
 
# File 'lib/matrix.rb', line 1560

def Vector.elements(array, copy = true)
  new convert_to_array(array, copy)
end

Instance Method Details

#*(x) ⇒ Object

Multiplies the vector by x, where x is a number or another vector.



1674
1675
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
 
# File 'lib/matrix.rb', line 1674

def *(x)
  case x
  when Numeric
    els = @elements.collect{|e| e * x}
    Vector.elements(els, false)
  when Matrix
    Matrix.column_vector(self) * x
  when Vector
    Vector.Raise ErrOperationNotDefined, "*", self.class, x.class
  else
    apply_through_coercion(x, __method__)
  end
end

#+(v) ⇒ Object

Vector addition.



1691
1692
1693
1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
 
# File 'lib/matrix.rb', line 1691

def +(v)
  case v
  when Vector
    Vector.Raise ErrDimensionMismatch if size != v.size
    els = collect2(v) {|v1, v2|
      v1 + v2
    }
    Vector.elements(els, false)
  when Matrix
    Matrix.column_vector(self) + v
  else
    apply_through_coercion(v, __method__)
  end
end

#-(v) ⇒ Object

Vector subtraction.



1709
1710
1711
1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
 
# File 'lib/matrix.rb', line 1709

def -(v)
  case v
  when Vector
    Vector.Raise ErrDimensionMismatch if size != v.size
    els = collect2(v) {|v1, v2|
      v1 - v2
    }
    Vector.elements(els, false)
  when Matrix
    Matrix.column_vector(self) - v
  else
    apply_through_coercion(v, __method__)
  end
end

#/(x) ⇒ Object

Vector division.



1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
 
# File 'lib/matrix.rb', line 1727

def /(x)
  case x
  when Numeric
    els = @elements.collect{|e| e / x}
    Vector.elements(els, false)
  when Matrix, Vector
    Vector.Raise ErrOperationNotDefined, "/", self.class, x.class
  else
    apply_through_coercion(x, __method__)
  end
end

#==(other) ⇒ Object

Returns true iff the two vectors have the same elements in the same order.



1643
1644
1645
1646
 
# File 'lib/matrix.rb', line 1643

def ==(other)
  return false unless Vector === other
  @elements == other.elements
end

#[](i) ⇒ Object Also known as: element, component

Returns element number i (starting at zero) of the vector.



1577
1578
1579
 
# File 'lib/matrix.rb', line 1577

def [](i)
  @elements[i]
end

#clone ⇒ Object

Return a copy of the vector.



1656
1657
1658
 
# File 'lib/matrix.rb', line 1656

def clone
  Vector.elements(@elements)
end

#coerce(other) ⇒ Object

The coerce method provides support for Ruby type coercion. This coercion mechanism is used by Ruby to handle mixed-type numeric operations: it is intended to find a compatible common type between the two operands of the operator. See also Numeric#coerce.



1840
1841
1842
1843
1844
1845
1846
1847
 
# File 'lib/matrix.rb', line 1840

def coerce(other)
  case other
  when Numeric
    return Matrix::Scalar.new(other), self
  else
    raise TypeError, "#{self.class} can't be coerced into #{other.class}"
  end
end

#collect(&block) ⇒ Object Also known as: map

Like Array#collect.



1760
1761
1762
1763
1764
 
# File 'lib/matrix.rb', line 1760

def collect(&block) # :yield: e
  return to_enum(:collect) unless block_given?
  els = @elements.collect(&block)
  Vector.elements(els, false)
end

#collect2(v) ⇒ Object

Collects (as in Enumerable#collect) over the elements of this vector and v in conjunction.

Raises:

  • (TypeError) 


1627
1628
1629
1630
1631
1632
1633
1634
 
# File 'lib/matrix.rb', line 1627

def collect2(v) # :yield: e1, e2
  raise TypeError, "Integer is not like Vector" if v.kind_of?(Integer)
  Vector.Raise ErrDimensionMismatch if size != v.size
  return to_enum(:collect2, v) unless block_given?
  Array.new(size) do |i|
    yield @elements[i], v[i]
  end
end

#covector ⇒ Object

Creates a single-row matrix from this vector.



1807
1808
1809
 
# File 'lib/matrix.rb', line 1807

def covector
  Matrix.row_vector(self)
end

#each(&block) ⇒ Object

Iterate over the elements of this vector



1604
1605
1606
1607
1608
 
# File 'lib/matrix.rb', line 1604

def each(&block)
  return to_enum(:each) unless block_given?
  @elements.each(&block)
  self
end

#each2(v) ⇒ Object

Iterate over the elements of this vector and v in conjunction.

Raises:

  • (TypeError) 


1613
1614
1615
1616
1617
1618
1619
1620
1621
 
# File 'lib/matrix.rb', line 1613

def each2(v) # :yield: e1, e2
  raise TypeError, "Integer is not like Vector" if v.kind_of?(Integer)
  Vector.Raise ErrDimensionMismatch if size != v.size
  return to_enum(:each2, v) unless block_given?
  size.times do |i|
    yield @elements[i], v[i]
  end
  self
end

#elements_to_f ⇒ Object 



1818
1819
1820
1821
 
# File 'lib/matrix.rb', line 1818

def elements_to_f
  warn "#{caller(1)[0]}: warning: Vector#elements_to_f is deprecated"
  map(&:to_f)
end

#elements_to_i ⇒ Object 



1823
1824
1825
1826
 
# File 'lib/matrix.rb', line 1823

def elements_to_i
  warn "#{caller(1)[0]}: warning: Vector#elements_to_i is deprecated"
  map(&:to_i)
end

#elements_to_r ⇒ Object 



1828
1829
1830
1831
 
# File 'lib/matrix.rb', line 1828

def elements_to_r
  warn "#{caller(1)[0]}: warning: Vector#elements_to_r is deprecated"
  map(&:to_r)
end

#eql?(other) ⇒ Boolean

Returns:

  • (Boolean) 


1648
1649
1650
1651
 
# File 'lib/matrix.rb', line 1648

def eql?(other)
  return false unless Vector === other
  @elements.eql? other.elements
end

#hash ⇒ Object

Return a hash-code for the vector.



1663
1664
1665
 
# File 'lib/matrix.rb', line 1663

def hash
  @elements.hash
end

#inner_product(v) ⇒ Object

Returns the inner product of this vector with the other.

Vector[4,7].inner_product Vector[10,1]  => 47



1747
1748
1749
1750
1751
1752
1753
1754
1755
 
# File 'lib/matrix.rb', line 1747

def inner_product(v)
  Vector.Raise ErrDimensionMismatch if size != v.size

  p = 0
  each2(v) {|v1, v2|
    p += v1 * v2
  }
  p
end

#inspect ⇒ Object

Overrides Object#inspect



1863
1864
1865
 
# File 'lib/matrix.rb', line 1863

def inspect
  "Vector" + @elements.inspect
end

#magnitude ⇒ Object Also known as: r, norm

Returns the modulus (Pythagorean distance) of the vector.

Vector[5,8,2].r => 9.643650761



1771
1772
1773
 
# File 'lib/matrix.rb', line 1771

def magnitude
  Math.sqrt(@elements.inject(0) {|v, e| v + e*e})
end

#map2(v, &block) ⇒ Object

Like Vector#collect2, but returns a Vector instead of an Array.



1780
1781
1782
1783
1784
 
# File 'lib/matrix.rb', line 1780

def map2(v, &block) # :yield: e1, e2
  return to_enum(:map2, v) unless block_given?
  els = collect2(v, &block)
  Vector.elements(els, false)
end

#normalize ⇒ Object

Returns a new vector with the same direction but with norm 1.

v = Vector[5,8,2].normalize
# => Vector[0.5184758473652127, 0.8295613557843402, 0.20739033894608505]
v.norm => 1.0

Raises:



1794
1795
1796
1797
1798
 
# File 'lib/matrix.rb', line 1794

def normalize
  n = magnitude
  raise ZeroVectorError, "Zero vectors can not be normalized" if n == 0
  self / n
end

#size ⇒ Object

Returns the number of elements in the vector.



1593
1594
1595
 
# File 'lib/matrix.rb', line 1593

def size
  @elements.size
end

#to_a ⇒ Object

Returns the elements of the vector in an array.



1814
1815
1816
 
# File 'lib/matrix.rb', line 1814

def to_a
  @elements.dup
end

#to_s ⇒ Object

Overrides Object#to_s



1856
1857
1858
 
# File 'lib/matrix.rb', line 1856

def to_s
  "Vector[" + @elements.join(", ") + "]"
end