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.



1570
1571
1572
1573
 
# File 'lib/matrix.rb', line 1570

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, ...]



1555
1556
1557
 
# File 'lib/matrix.rb', line 1555

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.



1563
1564
1565
 
# File 'lib/matrix.rb', line 1563

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.



1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
 
# File 'lib/matrix.rb', line 1677

def *(x)
  case x
  when Numeric
    els = @elements.collect{|e| e * x}
    self.class.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.



1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
1705
1706
1707
 
# File 'lib/matrix.rb', line 1694

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

#-(v) ⇒ Object

Vector subtraction.



1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
1723
1724
1725
 
# File 'lib/matrix.rb', line 1712

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

#/(x) ⇒ Object

Vector division.



1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
 
# File 'lib/matrix.rb', line 1730

def /(x)
  case x
  when Numeric
    els = @elements.collect{|e| e / x}
    self.class.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.



1646
1647
1648
1649
 
# File 'lib/matrix.rb', line 1646

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.



1580
1581
1582
 
# File 'lib/matrix.rb', line 1580

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

#clone ⇒ Object

Return a copy of the vector.



1659
1660
1661
 
# File 'lib/matrix.rb', line 1659

def clone
  self.class.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.



1843
1844
1845
1846
1847
1848
1849
1850
 
# File 'lib/matrix.rb', line 1843

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.



1763
1764
1765
1766
1767
 
# File 'lib/matrix.rb', line 1763

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

#collect2(v) ⇒ Object

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

Raises:

  • (TypeError) 


1630
1631
1632
1633
1634
1635
1636
1637
 
# File 'lib/matrix.rb', line 1630

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.



1810
1811
1812
 
# File 'lib/matrix.rb', line 1810

def covector
  Matrix.row_vector(self)
end

#each(&block) ⇒ Object

Iterate over the elements of this vector



1607
1608
1609
1610
1611
 
# File 'lib/matrix.rb', line 1607

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) 


1616
1617
1618
1619
1620
1621
1622
1623
1624
 
# File 'lib/matrix.rb', line 1616

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 



1821
1822
1823
1824
 
# File 'lib/matrix.rb', line 1821

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

#elements_to_i ⇒ Object 



1826
1827
1828
1829
 
# File 'lib/matrix.rb', line 1826

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

#elements_to_r ⇒ Object 



1831
1832
1833
1834
 
# File 'lib/matrix.rb', line 1831

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

#eql?(other) ⇒ Boolean

Returns:

  • (Boolean) 


1651
1652
1653
1654
 
# File 'lib/matrix.rb', line 1651

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

#hash ⇒ Object

Return a hash-code for the vector.



1666
1667
1668
 
# File 'lib/matrix.rb', line 1666

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



1750
1751
1752
1753
1754
1755
1756
1757
1758
 
# File 'lib/matrix.rb', line 1750

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

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

#inspect ⇒ Object

Overrides Object#inspect



1866
1867
1868
 
# File 'lib/matrix.rb', line 1866

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



1774
1775
1776
 
# File 'lib/matrix.rb', line 1774

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

#map2(v, &block) ⇒ Object

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



1783
1784
1785
1786
1787
 
# File 'lib/matrix.rb', line 1783

def map2(v, &block) # :yield: e1, e2
  return to_enum(:map2, v) unless block_given?
  els = collect2(v, &block)
  self.class.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:



1797
1798
1799
1800
1801
 
# File 'lib/matrix.rb', line 1797

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.



1596
1597
1598
 
# File 'lib/matrix.rb', line 1596

def size
  @elements.size
end

#to_a ⇒ Object

Returns the elements of the vector in an array.



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

def to_a
  @elements.dup
end

#to_s ⇒ Object

Overrides Object#to_s



1859
1860
1861
 
# File 'lib/matrix.rb', line 1859

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