def testCrossProduct(self):
"""Test the cross_product function"""
a = Numeric.array([1,2,3])
b = Numeric.array([4,5,6])
assert_eq(Numeric.cross_product(a,b), [-3, 6, -3])
a = Numeric.array([1,2])
b = Numeric.array([4,5])
assert_eq(Numeric.cross_product(a,b), -3)
a = Numeric.array([[1,2,3], [4,5,6]])
b = Numeric.array([7, 8, 9])
assert_eq(Numeric.cross_product(a,b), [[-6,12,-6],[-3,6,-3]])
a = Numeric.array([[1,2,3], [4,5,6]])
b = Numeric.array([[10,11,12], [7,8,9]])
assert_eq(Numeric.cross_product(a,b,axis1=0,axis2=0), [-33,-39,-45])
assert_eq(Numeric.cross_product(a,b), [[-9,18,-9], [-3,6,-3]])
