Feature #10074 ยป generalize_Vector#cross_product.patch
| lib/matrix.rb | ||
|---|---|---|
|
end
|
||
|
#
|
||
|
# Returns the cross product of this vector with the other.
|
||
|
# Returns the cross product of this vector with the others.
|
||
|
# Vector[1, 0, 0].cross_product Vector[0, 1, 0] => Vector[0, 0, 1]
|
||
|
#
|
||
|
def cross_product(v)
|
||
|
Vector.Raise ErrDimensionMismatch unless size == v.size && v.size == 3
|
||
|
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] ]
|
||
|
def cross_product(*vs)
|
||
|
Vector.Raise ErrDimensionMismatch unless vs.all?{ |v| v.size == size } && size == vs.count + 2
|
||
|
rows = Array.new(size) {|i| Vector.basis(size, i) }, self, *vs
|
||
|
Matrix.rows(rows).laplace_expansion(0, :row)
|
||
|
end
|
||
|
#
|
||
| test/matrix/test_vector.rb | ||
|---|---|---|
|
end
|
||
|
def test_cross_product
|
||
|
v = Vector[1, 0, 0].cross_product Vector[0, 1, 0]
|
||
|
assert_equal(Vector[0, 0, 1], v)
|
||
|
v2 = Vector[1, 2].cross_product
|
||
|
assert_equal(Vector[2, -1], v2)
|
||
|
v3 = Vector[1, 0, 0].cross_product Vector[0, 1, 0]
|
||
|
assert_equal(Vector[0, 0, 1], v3)
|
||
|
v4 = Vector[3, 5, 2, 1].cross_product(Vector[4, 3, 1, 8], Vector[2, 9, 4, 3])
|
||
|
assert_equal(Vector[-16, 65, -139, 1], v4)
|
||
|
assert_raise(Vector::ErrDimensionMismatch) { Vector[1, 2].cross_product(Vector[1, 4, 9]) }
|
||
|
assert_raise(Vector::ErrDimensionMismatch) { Vector[1, 2].cross_product(Vector[2, -1]) }
|
||
|
end
|
||
|
def test_r
|
||