класс Vector

# File lib/matrix.rb, line 2136
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
    raise ErrOperationNotDefined, ["*", self.class, x.class]
  else
    apply_through_coercion(x, __method__)
  end
end

# File lib/matrix.rb, line 2153
def +(v)
  case v
  when 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

# File lib/matrix.rb, line 2201
def +@
  self
end

# File lib/matrix.rb, line 2171
def -(v)
  case v
  when 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

# File lib/matrix.rb, line 2205
def -@
  collect {|e| -e }
end

# File lib/matrix.rb, line 2189
def /(x)
  case x
  when Numeric
    els = @elements.collect{|e| e / x}
    self.class.elements(els, false)
  when Matrix, Vector
    raise ErrOperationNotDefined, ["/", self.class, x.class]
  else
    apply_through_coercion(x, __method__)
  end
end

# File lib/matrix.rb, line 2112
def ==(other)
  return false unless Vector === other
  @elements == other.elements
end

# File lib/matrix.rb, line 1951
def [](i)
  @elements[i]
end

# File lib/matrix.rb, line 1966
def []=(i, v)
  raise FrozenError, "can't modify frozen Vector" if frozen?
  if i.is_a?(Range)
    range = Matrix::CoercionHelper.check_range(i, size, :vector)
    set_range(range, v)
  else
    index = Matrix::CoercionHelper.check_int(i, size, :index)
    set_value(index, v)
  end
end

# File lib/matrix.rb, line 2321
def angle_with(v)
  raise TypeError, "Expected a Vector, got a #{v.class}" unless v.is_a?(Vector)
  raise ErrDimensionMismatch if size != v.size
  prod = magnitude * v.magnitude
  raise ZeroVectorError, "Can't get angle of zero vector" if prod == 0
  dot = inner_product(v)
  if dot.abs >= prod
    dot.positive? ? 0 : Math::PI
  else
    Math.acos(dot / prod)
  end
end

# File lib/matrix.rb, line 2381
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

# File lib/matrix.rb, line 2265
def collect(&block) # :yield: e
  return to_enum(:collect) unless block_given?
  els = @elements.collect(&block)
  self.class.elements(els, false)
end

# File lib/matrix.rb, line 2275
def collect!(&block)
  return to_enum(:collect!) unless block_given?
  raise FrozenError, "can't modify frozen Vector" if frozen?
  @elements.collect!(&block)
  self
end

# File lib/matrix.rb, line 2040
def collect2(v) # :yield: e1, e2
  raise TypeError, "Integer is not like Vector" if v.kind_of?(Integer)
  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

# File lib/matrix.rb, line 2341
def covector
  Matrix.row_vector(self)
end

# File lib/matrix.rb, line 2240
def cross_product(*vs)
  raise ErrOperationNotDefined, "cross product is not defined on vectors of dimension #{size}" unless size >= 2
  raise ArgumentError, "wrong number of arguments (#{vs.size} for #{size - 2})" unless vs.size == size - 2
  vs.each do |v|
    raise TypeError, "expected Vector, got #{v.class}" unless v.is_a? Vector
    raise ErrDimensionMismatch unless v.size == size
  end
  case size
  when 2
    Vector[-@elements[1], @elements[0]]
  when 3
    v = vs[0]
    Vector[ v[2]*@elements[1] - v[1]*@elements[2],
      v[0]*@elements[2] - v[2]*@elements[0],
      v[1]*@elements[0] - v[0]*@elements[1] ]
  else
    rows = self, *vs, Array.new(size) {|i| Vector.basis(size: size, index: i) }
    Matrix.rows(rows).laplace_expansion(row: size - 1)
  end
end

# File lib/matrix.rb, line 2017
def each(&block)
  return to_enum(:each) unless block_given?
  @elements.each(&block)
  self
end

# File lib/matrix.rb, line 2026
def each2(v) # :yield: e1, e2
  raise TypeError, "Integer is not like Vector" if v.kind_of?(Integer)
  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

# File lib/matrix.rb, line 2359
def elements_to_f
  warn "Vector#elements_to_f is deprecated", uplevel: 1
  map(&:to_f)
end

# File lib/matrix.rb, line 2364
def elements_to_i
  warn "Vector#elements_to_i is deprecated", uplevel: 1
  map(&:to_i)
end

# File lib/matrix.rb, line 2369
def elements_to_r
  warn "Vector#elements_to_r is deprecated", uplevel: 1
  map(&:to_r)
end

# File lib/matrix.rb, line 2117
def eql?(other)
  return false unless Vector === other
  @elements.eql? other.elements
end

# File lib/matrix.rb, line 2091
def freeze
  @elements.freeze
  super
end

# File lib/matrix.rb, line 2125
def hash
  @elements.hash
end

# File lib/matrix.rb, line 2080
def independent?(*vs)
  self.class.independent?(self, *vs)
end

# File lib/matrix.rb, line 2217
def inner_product(v)
  raise ErrDimensionMismatch if size != v.size

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

# File lib/matrix.rb, line 2404
def inspect
  "Vector" + @elements.inspect
end

# File lib/matrix.rb, line 2287
def magnitude
  Math.sqrt(@elements.inject(0) {|v, e| v + e.abs2})
end

# File lib/matrix.rb, line 2296
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

# File lib/matrix.rb, line 2310
def normalize
  n = magnitude
  raise ZeroVectorError, "Zero vectors can not be normalized" if n == 0
  self / n
end

# File lib/matrix.rb, line 1999
def round(ndigits=0)
  map{|e| e.round(ndigits)}
end

# File lib/matrix.rb, line 2006
def size
  @elements.size
end

# File lib/matrix.rb, line 2348
def to_a
  @elements.dup
end

# File lib/matrix.rb, line 2355
def to_matrix
  Matrix.column_vector(self)
end

# File lib/matrix.rb, line 2397
def to_s
  "Vector[" + @elements.join(", ") + "]"
end

# File lib/matrix.rb, line 2087
def zero?
  all?(&:zero?)
end