111986Sandreas.sandberg@arm.comimport pytest 212391Sjason@lowepower.comfrom pybind11_tests import numpy_vectorize as m 311986Sandreas.sandberg@arm.com 412037Sandreas.sandberg@arm.compytestmark = pytest.requires_numpy 512037Sandreas.sandberg@arm.com 611986Sandreas.sandberg@arm.comwith pytest.suppress(ImportError): 711986Sandreas.sandberg@arm.com import numpy as np 811986Sandreas.sandberg@arm.com 911986Sandreas.sandberg@arm.com 1011986Sandreas.sandberg@arm.comdef test_vectorize(capture): 1112391Sjason@lowepower.com assert np.isclose(m.vectorized_func3(np.array(3 + 7j)), [6 + 14j]) 1211986Sandreas.sandberg@arm.com 1312391Sjason@lowepower.com for f in [m.vectorized_func, m.vectorized_func2]: 1411986Sandreas.sandberg@arm.com with capture: 1511986Sandreas.sandberg@arm.com assert np.isclose(f(1, 2, 3), 6) 1611986Sandreas.sandberg@arm.com assert capture == "my_func(x:int=1, y:float=2, z:float=3)" 1711986Sandreas.sandberg@arm.com with capture: 1811986Sandreas.sandberg@arm.com assert np.isclose(f(np.array(1), np.array(2), 3), 6) 1911986Sandreas.sandberg@arm.com assert capture == "my_func(x:int=1, y:float=2, z:float=3)" 2011986Sandreas.sandberg@arm.com with capture: 2111986Sandreas.sandberg@arm.com assert np.allclose(f(np.array([1, 3]), np.array([2, 4]), 3), [6, 36]) 2211986Sandreas.sandberg@arm.com assert capture == """ 2311986Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=3) 2411986Sandreas.sandberg@arm.com my_func(x:int=3, y:float=4, z:float=3) 2511986Sandreas.sandberg@arm.com """ 2611986Sandreas.sandberg@arm.com with capture: 2712037Sandreas.sandberg@arm.com a = np.array([[1, 2], [3, 4]], order='F') 2812037Sandreas.sandberg@arm.com b = np.array([[10, 20], [30, 40]], order='F') 2912037Sandreas.sandberg@arm.com c = 3 3012037Sandreas.sandberg@arm.com result = f(a, b, c) 3112037Sandreas.sandberg@arm.com assert np.allclose(result, a * b * c) 3212037Sandreas.sandberg@arm.com assert result.flags.f_contiguous 3312037Sandreas.sandberg@arm.com # All inputs are F order and full or singletons, so we the result is in col-major order: 3412037Sandreas.sandberg@arm.com assert capture == """ 3512037Sandreas.sandberg@arm.com my_func(x:int=1, y:float=10, z:float=3) 3612037Sandreas.sandberg@arm.com my_func(x:int=3, y:float=30, z:float=3) 3712037Sandreas.sandberg@arm.com my_func(x:int=2, y:float=20, z:float=3) 3812037Sandreas.sandberg@arm.com my_func(x:int=4, y:float=40, z:float=3) 3912037Sandreas.sandberg@arm.com """ 4012037Sandreas.sandberg@arm.com with capture: 4111986Sandreas.sandberg@arm.com a, b, c = np.array([[1, 3, 5], [7, 9, 11]]), np.array([[2, 4, 6], [8, 10, 12]]), 3 4211986Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 4311986Sandreas.sandberg@arm.com assert capture == """ 4411986Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=3) 4511986Sandreas.sandberg@arm.com my_func(x:int=3, y:float=4, z:float=3) 4611986Sandreas.sandberg@arm.com my_func(x:int=5, y:float=6, z:float=3) 4711986Sandreas.sandberg@arm.com my_func(x:int=7, y:float=8, z:float=3) 4811986Sandreas.sandberg@arm.com my_func(x:int=9, y:float=10, z:float=3) 4911986Sandreas.sandberg@arm.com my_func(x:int=11, y:float=12, z:float=3) 5011986Sandreas.sandberg@arm.com """ 5111986Sandreas.sandberg@arm.com with capture: 5211986Sandreas.sandberg@arm.com a, b, c = np.array([[1, 2, 3], [4, 5, 6]]), np.array([2, 3, 4]), 2 5311986Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 5411986Sandreas.sandberg@arm.com assert capture == """ 5511986Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=2) 5611986Sandreas.sandberg@arm.com my_func(x:int=2, y:float=3, z:float=2) 5711986Sandreas.sandberg@arm.com my_func(x:int=3, y:float=4, z:float=2) 5811986Sandreas.sandberg@arm.com my_func(x:int=4, y:float=2, z:float=2) 5911986Sandreas.sandberg@arm.com my_func(x:int=5, y:float=3, z:float=2) 6011986Sandreas.sandberg@arm.com my_func(x:int=6, y:float=4, z:float=2) 6111986Sandreas.sandberg@arm.com """ 6211986Sandreas.sandberg@arm.com with capture: 6311986Sandreas.sandberg@arm.com a, b, c = np.array([[1, 2, 3], [4, 5, 6]]), np.array([[2], [3]]), 2 6411986Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 6511986Sandreas.sandberg@arm.com assert capture == """ 6611986Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=2) 6711986Sandreas.sandberg@arm.com my_func(x:int=2, y:float=2, z:float=2) 6811986Sandreas.sandberg@arm.com my_func(x:int=3, y:float=2, z:float=2) 6911986Sandreas.sandberg@arm.com my_func(x:int=4, y:float=3, z:float=2) 7011986Sandreas.sandberg@arm.com my_func(x:int=5, y:float=3, z:float=2) 7111986Sandreas.sandberg@arm.com my_func(x:int=6, y:float=3, z:float=2) 7211986Sandreas.sandberg@arm.com """ 7312037Sandreas.sandberg@arm.com with capture: 7412037Sandreas.sandberg@arm.com a, b, c = np.array([[1, 2, 3], [4, 5, 6]], order='F'), np.array([[2], [3]]), 2 7512037Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 7612037Sandreas.sandberg@arm.com assert capture == """ 7712037Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=2) 7812037Sandreas.sandberg@arm.com my_func(x:int=2, y:float=2, z:float=2) 7912037Sandreas.sandberg@arm.com my_func(x:int=3, y:float=2, z:float=2) 8012037Sandreas.sandberg@arm.com my_func(x:int=4, y:float=3, z:float=2) 8112037Sandreas.sandberg@arm.com my_func(x:int=5, y:float=3, z:float=2) 8212037Sandreas.sandberg@arm.com my_func(x:int=6, y:float=3, z:float=2) 8312037Sandreas.sandberg@arm.com """ 8412037Sandreas.sandberg@arm.com with capture: 8512037Sandreas.sandberg@arm.com a, b, c = np.array([[1, 2, 3], [4, 5, 6]])[::, ::2], np.array([[2], [3]]), 2 8612037Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 8712037Sandreas.sandberg@arm.com assert capture == """ 8812037Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=2) 8912037Sandreas.sandberg@arm.com my_func(x:int=3, y:float=2, z:float=2) 9012037Sandreas.sandberg@arm.com my_func(x:int=4, y:float=3, z:float=2) 9112037Sandreas.sandberg@arm.com my_func(x:int=6, y:float=3, z:float=2) 9212037Sandreas.sandberg@arm.com """ 9312037Sandreas.sandberg@arm.com with capture: 9412037Sandreas.sandberg@arm.com a, b, c = np.array([[1, 2, 3], [4, 5, 6]], order='F')[::, ::2], np.array([[2], [3]]), 2 9512037Sandreas.sandberg@arm.com assert np.allclose(f(a, b, c), a * b * c) 9612037Sandreas.sandberg@arm.com assert capture == """ 9712037Sandreas.sandberg@arm.com my_func(x:int=1, y:float=2, z:float=2) 9812037Sandreas.sandberg@arm.com my_func(x:int=3, y:float=2, z:float=2) 9912037Sandreas.sandberg@arm.com my_func(x:int=4, y:float=3, z:float=2) 10012037Sandreas.sandberg@arm.com my_func(x:int=6, y:float=3, z:float=2) 10112037Sandreas.sandberg@arm.com """ 10211986Sandreas.sandberg@arm.com 10311986Sandreas.sandberg@arm.com 10411986Sandreas.sandberg@arm.comdef test_type_selection(): 10512391Sjason@lowepower.com assert m.selective_func(np.array([1], dtype=np.int32)) == "Int branch taken." 10612391Sjason@lowepower.com assert m.selective_func(np.array([1.0], dtype=np.float32)) == "Float branch taken." 10712391Sjason@lowepower.com assert m.selective_func(np.array([1.0j], dtype=np.complex64)) == "Complex float branch taken." 10811986Sandreas.sandberg@arm.com 10911986Sandreas.sandberg@arm.com 11011986Sandreas.sandberg@arm.comdef test_docs(doc): 11112391Sjason@lowepower.com assert doc(m.vectorized_func) == """ 11212037Sandreas.sandberg@arm.com vectorized_func(arg0: numpy.ndarray[int32], arg1: numpy.ndarray[float32], arg2: numpy.ndarray[float64]) -> object 11311986Sandreas.sandberg@arm.com """ # noqa: E501 line too long 11412037Sandreas.sandberg@arm.com 11512037Sandreas.sandberg@arm.com 11612037Sandreas.sandberg@arm.comdef test_trivial_broadcasting(): 11712391Sjason@lowepower.com trivial, vectorized_is_trivial = m.trivial, m.vectorized_is_trivial 11812037Sandreas.sandberg@arm.com 11912037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, 2, 3) == trivial.c_trivial 12012037Sandreas.sandberg@arm.com assert vectorized_is_trivial(np.array(1), np.array(2), 3) == trivial.c_trivial 12112037Sandreas.sandberg@arm.com assert vectorized_is_trivial(np.array([1, 3]), np.array([2, 4]), 3) == trivial.c_trivial 12212037Sandreas.sandberg@arm.com assert trivial.c_trivial == vectorized_is_trivial( 12312037Sandreas.sandberg@arm.com np.array([[1, 3, 5], [7, 9, 11]]), np.array([[2, 4, 6], [8, 10, 12]]), 3) 12412037Sandreas.sandberg@arm.com assert vectorized_is_trivial( 12512037Sandreas.sandberg@arm.com np.array([[1, 2, 3], [4, 5, 6]]), np.array([2, 3, 4]), 2) == trivial.non_trivial 12612037Sandreas.sandberg@arm.com assert vectorized_is_trivial( 12712037Sandreas.sandberg@arm.com np.array([[1, 2, 3], [4, 5, 6]]), np.array([[2], [3]]), 2) == trivial.non_trivial 12812037Sandreas.sandberg@arm.com z1 = np.array([[1, 2, 3, 4], [5, 6, 7, 8]], dtype='int32') 12912037Sandreas.sandberg@arm.com z2 = np.array(z1, dtype='float32') 13012037Sandreas.sandberg@arm.com z3 = np.array(z1, dtype='float64') 13112037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1, z2, z3) == trivial.c_trivial 13212037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, z2, z3) == trivial.c_trivial 13312037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1, 1, z3) == trivial.c_trivial 13412037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1, z2, 1) == trivial.c_trivial 13512037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1[::2, ::2], 1, 1) == trivial.non_trivial 13612037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, 1, z1[::2, ::2]) == trivial.c_trivial 13712037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, 1, z3[::2, ::2]) == trivial.non_trivial 13812037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1, 1, z3[1::4, 1::4]) == trivial.c_trivial 13912037Sandreas.sandberg@arm.com 14012037Sandreas.sandberg@arm.com y1 = np.array(z1, order='F') 14112037Sandreas.sandberg@arm.com y2 = np.array(y1) 14212037Sandreas.sandberg@arm.com y3 = np.array(y1) 14312037Sandreas.sandberg@arm.com assert vectorized_is_trivial(y1, y2, y3) == trivial.f_trivial 14412037Sandreas.sandberg@arm.com assert vectorized_is_trivial(y1, 1, 1) == trivial.f_trivial 14512037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, y2, 1) == trivial.f_trivial 14612037Sandreas.sandberg@arm.com assert vectorized_is_trivial(1, 1, y3) == trivial.f_trivial 14712037Sandreas.sandberg@arm.com assert vectorized_is_trivial(y1, z2, 1) == trivial.non_trivial 14812037Sandreas.sandberg@arm.com assert vectorized_is_trivial(z1[1::4, 1::4], y2, 1) == trivial.f_trivial 14912037Sandreas.sandberg@arm.com assert vectorized_is_trivial(y1[1::4, 1::4], z2, 1) == trivial.c_trivial 15012037Sandreas.sandberg@arm.com 15112391Sjason@lowepower.com assert m.vectorized_func(z1, z2, z3).flags.c_contiguous 15212391Sjason@lowepower.com assert m.vectorized_func(y1, y2, y3).flags.f_contiguous 15312391Sjason@lowepower.com assert m.vectorized_func(z1, 1, 1).flags.c_contiguous 15412391Sjason@lowepower.com assert m.vectorized_func(1, y2, 1).flags.f_contiguous 15512391Sjason@lowepower.com assert m.vectorized_func(z1[1::4, 1::4], y2, 1).flags.f_contiguous 15612391Sjason@lowepower.com assert m.vectorized_func(y1[1::4, 1::4], z2, 1).flags.c_contiguous 15712391Sjason@lowepower.com 15812391Sjason@lowepower.com 15912391Sjason@lowepower.comdef test_passthrough_arguments(doc): 16012391Sjason@lowepower.com assert doc(m.vec_passthrough) == ( 16112391Sjason@lowepower.com "vec_passthrough(" + ", ".join([ 16212391Sjason@lowepower.com "arg0: float", 16312391Sjason@lowepower.com "arg1: numpy.ndarray[float64]", 16412391Sjason@lowepower.com "arg2: numpy.ndarray[float64]", 16512391Sjason@lowepower.com "arg3: numpy.ndarray[int32]", 16612391Sjason@lowepower.com "arg4: int", 16712391Sjason@lowepower.com "arg5: m.numpy_vectorize.NonPODClass", 16812391Sjason@lowepower.com "arg6: numpy.ndarray[float64]"]) + ") -> object") 16912391Sjason@lowepower.com 17012391Sjason@lowepower.com b = np.array([[10, 20, 30]], dtype='float64') 17112391Sjason@lowepower.com c = np.array([100, 200]) # NOT a vectorized argument 17212391Sjason@lowepower.com d = np.array([[1000], [2000], [3000]], dtype='int') 17312391Sjason@lowepower.com g = np.array([[1000000, 2000000, 3000000]], dtype='int') # requires casting 17412391Sjason@lowepower.com assert np.all( 17512391Sjason@lowepower.com m.vec_passthrough(1, b, c, d, 10000, m.NonPODClass(100000), g) == 17612391Sjason@lowepower.com np.array([[1111111, 2111121, 3111131], 17712391Sjason@lowepower.com [1112111, 2112121, 3112131], 17812391Sjason@lowepower.com [1113111, 2113121, 3113131]])) 17912391Sjason@lowepower.com 18012391Sjason@lowepower.com 18112391Sjason@lowepower.comdef test_method_vectorization(): 18212391Sjason@lowepower.com o = m.VectorizeTestClass(3) 18312391Sjason@lowepower.com x = np.array([1, 2], dtype='int') 18412391Sjason@lowepower.com y = np.array([[10], [20]], dtype='float32') 18512391Sjason@lowepower.com assert np.all(o.method(x, y) == [[14, 15], [24, 25]]) 18612391Sjason@lowepower.com 18712391Sjason@lowepower.com 18812391Sjason@lowepower.comdef test_array_collapse(): 18912391Sjason@lowepower.com assert not isinstance(m.vectorized_func(1, 2, 3), np.ndarray) 19012391Sjason@lowepower.com assert not isinstance(m.vectorized_func(np.array(1), 2, 3), np.ndarray) 19112391Sjason@lowepower.com z = m.vectorized_func([1], 2, 3) 19212391Sjason@lowepower.com assert isinstance(z, np.ndarray) 19312391Sjason@lowepower.com assert z.shape == (1, ) 19412391Sjason@lowepower.com z = m.vectorized_func(1, [[[2]]], 3) 19512391Sjason@lowepower.com assert isinstance(z, np.ndarray) 19612391Sjason@lowepower.com assert z.shape == (1, 1, 1) 197