test_numpy_vectorize.cpp revision 12037
111986Sandreas.sandberg@arm.com/* 211986Sandreas.sandberg@arm.com tests/test_numpy_vectorize.cpp -- auto-vectorize functions over NumPy array 311986Sandreas.sandberg@arm.com arguments 411986Sandreas.sandberg@arm.com 511986Sandreas.sandberg@arm.com Copyright (c) 2016 Wenzel Jakob <wenzel.jakob@epfl.ch> 611986Sandreas.sandberg@arm.com 711986Sandreas.sandberg@arm.com All rights reserved. Use of this source code is governed by a 811986Sandreas.sandberg@arm.com BSD-style license that can be found in the LICENSE file. 911986Sandreas.sandberg@arm.com*/ 1011986Sandreas.sandberg@arm.com 1111986Sandreas.sandberg@arm.com#include "pybind11_tests.h" 1211986Sandreas.sandberg@arm.com#include <pybind11/numpy.h> 1311986Sandreas.sandberg@arm.com 1411986Sandreas.sandberg@arm.comdouble my_func(int x, float y, double z) { 1511986Sandreas.sandberg@arm.com py::print("my_func(x:int={}, y:float={:.0f}, z:float={:.0f})"_s.format(x, y, z)); 1611986Sandreas.sandberg@arm.com return (float) x*y*z; 1711986Sandreas.sandberg@arm.com} 1811986Sandreas.sandberg@arm.com 1911986Sandreas.sandberg@arm.comstd::complex<double> my_func3(std::complex<double> c) { 2011986Sandreas.sandberg@arm.com return c * std::complex<double>(2.f); 2111986Sandreas.sandberg@arm.com} 2211986Sandreas.sandberg@arm.com 2311986Sandreas.sandberg@arm.comtest_initializer numpy_vectorize([](py::module &m) { 2411986Sandreas.sandberg@arm.com // Vectorize all arguments of a function (though non-vector arguments are also allowed) 2511986Sandreas.sandberg@arm.com m.def("vectorized_func", py::vectorize(my_func)); 2611986Sandreas.sandberg@arm.com 2711986Sandreas.sandberg@arm.com // Vectorize a lambda function with a capture object (e.g. to exclude some arguments from the vectorization) 2811986Sandreas.sandberg@arm.com m.def("vectorized_func2", 2911986Sandreas.sandberg@arm.com [](py::array_t<int> x, py::array_t<float> y, float z) { 3011986Sandreas.sandberg@arm.com return py::vectorize([z](int x, float y) { return my_func(x, y, z); })(x, y); 3111986Sandreas.sandberg@arm.com } 3211986Sandreas.sandberg@arm.com ); 3311986Sandreas.sandberg@arm.com 3411986Sandreas.sandberg@arm.com // Vectorize a complex-valued function 3511986Sandreas.sandberg@arm.com m.def("vectorized_func3", py::vectorize(my_func3)); 3611986Sandreas.sandberg@arm.com 3711986Sandreas.sandberg@arm.com /// Numpy function which only accepts specific data types 3811986Sandreas.sandberg@arm.com m.def("selective_func", [](py::array_t<int, py::array::c_style>) { return "Int branch taken."; }); 3911986Sandreas.sandberg@arm.com m.def("selective_func", [](py::array_t<float, py::array::c_style>) { return "Float branch taken."; }); 4011986Sandreas.sandberg@arm.com m.def("selective_func", [](py::array_t<std::complex<float>, py::array::c_style>) { return "Complex float branch taken."; }); 4112037Sandreas.sandberg@arm.com 4212037Sandreas.sandberg@arm.com 4312037Sandreas.sandberg@arm.com // Internal optimization test for whether the input is trivially broadcastable: 4412037Sandreas.sandberg@arm.com py::enum_<py::detail::broadcast_trivial>(m, "trivial") 4512037Sandreas.sandberg@arm.com .value("f_trivial", py::detail::broadcast_trivial::f_trivial) 4612037Sandreas.sandberg@arm.com .value("c_trivial", py::detail::broadcast_trivial::c_trivial) 4712037Sandreas.sandberg@arm.com .value("non_trivial", py::detail::broadcast_trivial::non_trivial); 4812037Sandreas.sandberg@arm.com m.def("vectorized_is_trivial", []( 4912037Sandreas.sandberg@arm.com py::array_t<int, py::array::forcecast> arg1, 5012037Sandreas.sandberg@arm.com py::array_t<float, py::array::forcecast> arg2, 5112037Sandreas.sandberg@arm.com py::array_t<double, py::array::forcecast> arg3 5212037Sandreas.sandberg@arm.com ) { 5312037Sandreas.sandberg@arm.com size_t ndim; 5412037Sandreas.sandberg@arm.com std::vector<size_t> shape; 5512037Sandreas.sandberg@arm.com std::array<py::buffer_info, 3> buffers {{ arg1.request(), arg2.request(), arg3.request() }}; 5612037Sandreas.sandberg@arm.com return py::detail::broadcast(buffers, ndim, shape); 5712037Sandreas.sandberg@arm.com }); 5811986Sandreas.sandberg@arm.com}); 59