test_numpy_vectorize.cpp revision 12037
113558Snikos.nikoleris@arm.com/* 28839Sandreas.hansson@arm.com tests/test_numpy_vectorize.cpp -- auto-vectorize functions over NumPy array 38839Sandreas.hansson@arm.com arguments 48839Sandreas.hansson@arm.com 58839Sandreas.hansson@arm.com Copyright (c) 2016 Wenzel Jakob <wenzel.jakob@epfl.ch> 68839Sandreas.hansson@arm.com 78839Sandreas.hansson@arm.com All rights reserved. Use of this source code is governed by a 88839Sandreas.hansson@arm.com BSD-style license that can be found in the LICENSE file. 98839Sandreas.hansson@arm.com*/ 108839Sandreas.hansson@arm.com 118839Sandreas.hansson@arm.com#include "pybind11_tests.h" 128839Sandreas.hansson@arm.com#include <pybind11/numpy.h> 133101Sstever@eecs.umich.edu 148579Ssteve.reinhardt@amd.comdouble my_func(int x, float y, double z) { 153101Sstever@eecs.umich.edu py::print("my_func(x:int={}, y:float={:.0f}, z:float={:.0f})"_s.format(x, y, z)); 163101Sstever@eecs.umich.edu return (float) x*y*z; 173101Sstever@eecs.umich.edu} 183101Sstever@eecs.umich.edu 193101Sstever@eecs.umich.edustd::complex<double> my_func3(std::complex<double> c) { 203101Sstever@eecs.umich.edu return c * std::complex<double>(2.f); 213101Sstever@eecs.umich.edu} 223101Sstever@eecs.umich.edu 233101Sstever@eecs.umich.edutest_initializer numpy_vectorize([](py::module &m) { 243101Sstever@eecs.umich.edu // Vectorize all arguments of a function (though non-vector arguments are also allowed) 253101Sstever@eecs.umich.edu m.def("vectorized_func", py::vectorize(my_func)); 263101Sstever@eecs.umich.edu 273101Sstever@eecs.umich.edu // Vectorize a lambda function with a capture object (e.g. to exclude some arguments from the vectorization) 283101Sstever@eecs.umich.edu m.def("vectorized_func2", 293101Sstever@eecs.umich.edu [](py::array_t<int> x, py::array_t<float> y, float z) { 303101Sstever@eecs.umich.edu return py::vectorize([z](int x, float y) { return my_func(x, y, z); })(x, y); 313101Sstever@eecs.umich.edu } 323101Sstever@eecs.umich.edu ); 333101Sstever@eecs.umich.edu 343101Sstever@eecs.umich.edu // Vectorize a complex-valued function 353101Sstever@eecs.umich.edu m.def("vectorized_func3", py::vectorize(my_func3)); 363101Sstever@eecs.umich.edu 373101Sstever@eecs.umich.edu /// Numpy function which only accepts specific data types 383101Sstever@eecs.umich.edu m.def("selective_func", [](py::array_t<int, py::array::c_style>) { return "Int branch taken."; }); 393101Sstever@eecs.umich.edu m.def("selective_func", [](py::array_t<float, py::array::c_style>) { return "Float branch taken."; }); 403101Sstever@eecs.umich.edu m.def("selective_func", [](py::array_t<std::complex<float>, py::array::c_style>) { return "Complex float branch taken."; }); 413101Sstever@eecs.umich.edu 427778Sgblack@eecs.umich.edu 438839Sandreas.hansson@arm.com // Internal optimization test for whether the input is trivially broadcastable: 443101Sstever@eecs.umich.edu py::enum_<py::detail::broadcast_trivial>(m, "trivial") 453101Sstever@eecs.umich.edu .value("f_trivial", py::detail::broadcast_trivial::f_trivial) 463101Sstever@eecs.umich.edu .value("c_trivial", py::detail::broadcast_trivial::c_trivial) 473101Sstever@eecs.umich.edu .value("non_trivial", py::detail::broadcast_trivial::non_trivial); 483101Sstever@eecs.umich.edu m.def("vectorized_is_trivial", []( 493101Sstever@eecs.umich.edu py::array_t<int, py::array::forcecast> arg1, 503101Sstever@eecs.umich.edu py::array_t<float, py::array::forcecast> arg2, 513101Sstever@eecs.umich.edu py::array_t<double, py::array::forcecast> arg3 523101Sstever@eecs.umich.edu ) { 533101Sstever@eecs.umich.edu size_t ndim; 543101Sstever@eecs.umich.edu std::vector<size_t> shape; 553101Sstever@eecs.umich.edu std::array<py::buffer_info, 3> buffers {{ arg1.request(), arg2.request(), arg3.request() }}; 563101Sstever@eecs.umich.edu return py::detail::broadcast(buffers, ndim, shape); 573101Sstever@eecs.umich.edu }); 583101Sstever@eecs.umich.edu}); 593101Sstever@eecs.umich.edu