diff options
Diffstat (limited to '3rdparty/pybind11/tests/test_numpy_array.cpp')
-rw-r--r-- | 3rdparty/pybind11/tests/test_numpy_array.cpp | 470 |
1 files changed, 278 insertions, 192 deletions
diff --git a/3rdparty/pybind11/tests/test_numpy_array.cpp b/3rdparty/pybind11/tests/test_numpy_array.cpp index a84de77f..69ddbe1e 100644 --- a/3rdparty/pybind11/tests/test_numpy_array.cpp +++ b/3rdparty/pybind11/tests/test_numpy_array.cpp @@ -7,12 +7,13 @@ BSD-style license that can be found in the LICENSE file. */ -#include "pybind11_tests.h" - #include <pybind11/numpy.h> #include <pybind11/stl.h> +#include "pybind11_tests.h" + #include <cstdint> +#include <utility> // Size / dtype checks. struct DtypeCheck { @@ -21,7 +22,7 @@ struct DtypeCheck { }; template <typename T> -DtypeCheck get_dtype_check(const char* name) { +DtypeCheck get_dtype_check(const char *name) { py::module_ np = py::module_::import("numpy"); DtypeCheck check{}; check.numpy = np.attr("dtype")(np.attr(name)); @@ -30,17 +31,15 @@ DtypeCheck get_dtype_check(const char* name) { } std::vector<DtypeCheck> get_concrete_dtype_checks() { - return { - // Normalization - get_dtype_check<std::int8_t>("int8"), - get_dtype_check<std::uint8_t>("uint8"), - get_dtype_check<std::int16_t>("int16"), - get_dtype_check<std::uint16_t>("uint16"), - get_dtype_check<std::int32_t>("int32"), - get_dtype_check<std::uint32_t>("uint32"), - get_dtype_check<std::int64_t>("int64"), - get_dtype_check<std::uint64_t>("uint64") - }; + return {// Normalization + get_dtype_check<std::int8_t>("int8"), + get_dtype_check<std::uint8_t>("uint8"), + get_dtype_check<std::int16_t>("int16"), + get_dtype_check<std::uint16_t>("uint16"), + get_dtype_check<std::int32_t>("int32"), + get_dtype_check<std::uint32_t>("uint32"), + get_dtype_check<std::int64_t>("int64"), + get_dtype_check<std::uint64_t>("uint64")}; } struct DtypeSizeCheck { @@ -79,43 +78,71 @@ using arr = py::array; using arr_t = py::array_t<uint16_t, 0>; static_assert(std::is_same<arr_t::value_type, uint16_t>::value, ""); -template<typename... Ix> arr data(const arr& a, Ix... index) { +template <typename... Ix> +arr data(const arr &a, Ix... index) { return arr(a.nbytes() - a.offset_at(index...), (const uint8_t *) a.data(index...)); } -template<typename... Ix> arr data_t(const arr_t& a, Ix... index) { +template <typename... Ix> +arr data_t(const arr_t &a, Ix... index) { return arr(a.size() - a.index_at(index...), a.data(index...)); } -template<typename... Ix> arr& mutate_data(arr& a, Ix... index) { - auto ptr = (uint8_t *) a.mutable_data(index...); - for (py::ssize_t i = 0; i < a.nbytes() - a.offset_at(index...); i++) +template <typename... Ix> +arr &mutate_data(arr &a, Ix... index) { + auto *ptr = (uint8_t *) a.mutable_data(index...); + for (py::ssize_t i = 0; i < a.nbytes() - a.offset_at(index...); i++) { ptr[i] = (uint8_t) (ptr[i] * 2); + } return a; } -template<typename... Ix> arr_t& mutate_data_t(arr_t& a, Ix... index) { +template <typename... Ix> +arr_t &mutate_data_t(arr_t &a, Ix... index) { auto ptr = a.mutable_data(index...); - for (py::ssize_t i = 0; i < a.size() - a.index_at(index...); i++) + for (py::ssize_t i = 0; i < a.size() - a.index_at(index...); i++) { ptr[i]++; + } + return a; +} + +template <typename... Ix> +py::ssize_t index_at(const arr &a, Ix... idx) { + return a.index_at(idx...); +} +template <typename... Ix> +py::ssize_t index_at_t(const arr_t &a, Ix... idx) { + return a.index_at(idx...); +} +template <typename... Ix> +py::ssize_t offset_at(const arr &a, Ix... idx) { + return a.offset_at(idx...); +} +template <typename... Ix> +py::ssize_t offset_at_t(const arr_t &a, Ix... idx) { + return a.offset_at(idx...); +} +template <typename... Ix> +py::ssize_t at_t(const arr_t &a, Ix... idx) { + return a.at(idx...); +} +template <typename... Ix> +arr_t &mutate_at_t(arr_t &a, Ix... idx) { + a.mutable_at(idx...)++; return a; } -template<typename... Ix> py::ssize_t index_at(const arr& a, Ix... idx) { return a.index_at(idx...); } -template<typename... Ix> py::ssize_t index_at_t(const arr_t& a, Ix... idx) { return a.index_at(idx...); } -template<typename... Ix> py::ssize_t offset_at(const arr& a, Ix... idx) { return a.offset_at(idx...); } -template<typename... Ix> py::ssize_t offset_at_t(const arr_t& a, Ix... idx) { return a.offset_at(idx...); } -template<typename... Ix> py::ssize_t at_t(const arr_t& a, Ix... idx) { return a.at(idx...); } -template<typename... Ix> arr_t& mutate_at_t(arr_t& a, Ix... idx) { a.mutable_at(idx...)++; return a; } - -#define def_index_fn(name, type) \ - sm.def(#name, [](type a) { return name(a); }); \ - sm.def(#name, [](type a, int i) { return name(a, i); }); \ - sm.def(#name, [](type a, int i, int j) { return name(a, i, j); }); \ +#define def_index_fn(name, type) \ + sm.def(#name, [](type a) { return name(a); }); \ + sm.def(#name, [](type a, int i) { return name(a, i); }); \ + sm.def(#name, [](type a, int i, int j) { return name(a, i, j); }); \ sm.def(#name, [](type a, int i, int j, int k) { return name(a, i, j, k); }); -template <typename T, typename T2> py::handle auxiliaries(T &&r, T2 &&r2) { - if (r.ndim() != 2) throw std::domain_error("error: ndim != 2"); +template <typename T, typename T2> +py::handle auxiliaries(T &&r, T2 &&r2) { + if (r.ndim() != 2) { + throw std::domain_error("error: ndim != 2"); + } py::list l; l.append(*r.data(0, 0)); l.append(*r2.mutable_data(0, 0)); @@ -133,16 +160,18 @@ template <typename T, typename T2> py::handle auxiliaries(T &&r, T2 &&r2) { static int data_i = 42; TEST_SUBMODULE(numpy_array, sm) { - try { py::module_::import("numpy"); } - catch (...) { return; } + try { + py::module_::import("numpy"); + } catch (const py::error_already_set &) { + return; + } // test_dtypes py::class_<DtypeCheck>(sm, "DtypeCheck") .def_readonly("numpy", &DtypeCheck::numpy) .def_readonly("pybind11", &DtypeCheck::pybind11) - .def("__repr__", [](const DtypeCheck& self) { - return py::str("<DtypeCheck numpy={} pybind11={}>").format( - self.numpy, self.pybind11); + .def("__repr__", [](const DtypeCheck &self) { + return py::str("<DtypeCheck numpy={} pybind11={}>").format(self.numpy, self.pybind11); }); sm.def("get_concrete_dtype_checks", &get_concrete_dtype_checks); @@ -150,41 +179,41 @@ TEST_SUBMODULE(numpy_array, sm) { .def_readonly("name", &DtypeSizeCheck::name) .def_readonly("size_cpp", &DtypeSizeCheck::size_cpp) .def_readonly("size_numpy", &DtypeSizeCheck::size_numpy) - .def("__repr__", [](const DtypeSizeCheck& self) { - return py::str("<DtypeSizeCheck name='{}' size_cpp={} size_numpy={} dtype={}>").format( - self.name, self.size_cpp, self.size_numpy, self.dtype); + .def("__repr__", [](const DtypeSizeCheck &self) { + return py::str("<DtypeSizeCheck name='{}' size_cpp={} size_numpy={} dtype={}>") + .format(self.name, self.size_cpp, self.size_numpy, self.dtype); }); sm.def("get_platform_dtype_size_checks", &get_platform_dtype_size_checks); // test_array_attributes - sm.def("ndim", [](const arr& a) { return a.ndim(); }); - sm.def("shape", [](const arr& a) { return arr(a.ndim(), a.shape()); }); - sm.def("shape", [](const arr& a, py::ssize_t dim) { return a.shape(dim); }); - sm.def("strides", [](const arr& a) { return arr(a.ndim(), a.strides()); }); - sm.def("strides", [](const arr& a, py::ssize_t dim) { return a.strides(dim); }); - sm.def("writeable", [](const arr& a) { return a.writeable(); }); - sm.def("size", [](const arr& a) { return a.size(); }); - sm.def("itemsize", [](const arr& a) { return a.itemsize(); }); - sm.def("nbytes", [](const arr& a) { return a.nbytes(); }); - sm.def("owndata", [](const arr& a) { return a.owndata(); }); + sm.def("ndim", [](const arr &a) { return a.ndim(); }); + sm.def("shape", [](const arr &a) { return arr(a.ndim(), a.shape()); }); + sm.def("shape", [](const arr &a, py::ssize_t dim) { return a.shape(dim); }); + sm.def("strides", [](const arr &a) { return arr(a.ndim(), a.strides()); }); + sm.def("strides", [](const arr &a, py::ssize_t dim) { return a.strides(dim); }); + sm.def("writeable", [](const arr &a) { return a.writeable(); }); + sm.def("size", [](const arr &a) { return a.size(); }); + sm.def("itemsize", [](const arr &a) { return a.itemsize(); }); + sm.def("nbytes", [](const arr &a) { return a.nbytes(); }); + sm.def("owndata", [](const arr &a) { return a.owndata(); }); // test_index_offset - def_index_fn(index_at, const arr&); - def_index_fn(index_at_t, const arr_t&); - def_index_fn(offset_at, const arr&); - def_index_fn(offset_at_t, const arr_t&); + def_index_fn(index_at, const arr &); + def_index_fn(index_at_t, const arr_t &); + def_index_fn(offset_at, const arr &); + def_index_fn(offset_at_t, const arr_t &); // test_data - def_index_fn(data, const arr&); - def_index_fn(data_t, const arr_t&); + def_index_fn(data, const arr &); + def_index_fn(data_t, const arr_t &); // test_mutate_data, test_mutate_readonly - def_index_fn(mutate_data, arr&); - def_index_fn(mutate_data_t, arr_t&); - def_index_fn(at_t, const arr_t&); - def_index_fn(mutate_at_t, arr_t&); + def_index_fn(mutate_data, arr &); + def_index_fn(mutate_data_t, arr_t &); + def_index_fn(at_t, const arr_t &); + def_index_fn(mutate_at_t, arr_t &); // test_make_c_f_array - sm.def("make_f_array", [] { return py::array_t<float>({ 2, 2 }, { 4, 8 }); }); - sm.def("make_c_array", [] { return py::array_t<float>({ 2, 2 }, { 8, 4 }); }); + sm.def("make_f_array", [] { return py::array_t<float>({2, 2}, {4, 8}); }); + sm.def("make_c_array", [] { return py::array_t<float>({2, 2}, {8, 4}); }); // test_empty_shaped_array sm.def("make_empty_shaped_array", [] { return py::array(py::dtype("f"), {}, {}); }); @@ -192,19 +221,17 @@ TEST_SUBMODULE(numpy_array, sm) { sm.def("scalar_int", []() { return py::array(py::dtype("i"), {}, {}, &data_i); }); // test_wrap - sm.def("wrap", [](py::array a) { - return py::array( - a.dtype(), - {a.shape(), a.shape() + a.ndim()}, - {a.strides(), a.strides() + a.ndim()}, - a.data(), - a - ); + sm.def("wrap", [](const py::array &a) { + return py::array(a.dtype(), + {a.shape(), a.shape() + a.ndim()}, + {a.strides(), a.strides() + a.ndim()}, + a.data(), + a); }); // test_numpy_view struct ArrayClass { - int data[2] = { 1, 2 }; + int data[2] = {1, 2}; ArrayClass() { py::print("ArrayClass()"); } ~ArrayClass() { py::print("~ArrayClass()"); } }; @@ -212,103 +239,121 @@ TEST_SUBMODULE(numpy_array, sm) { .def(py::init<>()) .def("numpy_view", [](py::object &obj) { py::print("ArrayClass::numpy_view()"); - auto &a = obj.cast<ArrayClass&>(); + auto &a = obj.cast<ArrayClass &>(); return py::array_t<int>({2}, {4}, a.data, obj); - } - ); + }); // test_cast_numpy_int64_to_uint64 - sm.def("function_taking_uint64", [](uint64_t) { }); + sm.def("function_taking_uint64", [](uint64_t) {}); // test_isinstance sm.def("isinstance_untyped", [](py::object yes, py::object no) { - return py::isinstance<py::array>(yes) && !py::isinstance<py::array>(no); + return py::isinstance<py::array>(std::move(yes)) + && !py::isinstance<py::array>(std::move(no)); }); - sm.def("isinstance_typed", [](py::object o) { + sm.def("isinstance_typed", [](const py::object &o) { return py::isinstance<py::array_t<double>>(o) && !py::isinstance<py::array_t<int>>(o); }); // test_constructors sm.def("default_constructors", []() { - return py::dict( - "array"_a=py::array(), - "array_t<int32>"_a=py::array_t<std::int32_t>(), - "array_t<double>"_a=py::array_t<double>() - ); + return py::dict("array"_a = py::array(), + "array_t<int32>"_a = py::array_t<std::int32_t>(), + "array_t<double>"_a = py::array_t<double>()); }); - sm.def("converting_constructors", [](py::object o) { - return py::dict( - "array"_a=py::array(o), - "array_t<int32>"_a=py::array_t<std::int32_t>(o), - "array_t<double>"_a=py::array_t<double>(o) - ); + sm.def("converting_constructors", [](const py::object &o) { + return py::dict("array"_a = py::array(o), + "array_t<int32>"_a = py::array_t<std::int32_t>(o), + "array_t<double>"_a = py::array_t<double>(o)); }); // test_overload_resolution - sm.def("overloaded", [](py::array_t<double>) { return "double"; }); - sm.def("overloaded", [](py::array_t<float>) { return "float"; }); - sm.def("overloaded", [](py::array_t<int>) { return "int"; }); - sm.def("overloaded", [](py::array_t<unsigned short>) { return "unsigned short"; }); - sm.def("overloaded", [](py::array_t<long long>) { return "long long"; }); - sm.def("overloaded", [](py::array_t<std::complex<double>>) { return "double complex"; }); - sm.def("overloaded", [](py::array_t<std::complex<float>>) { return "float complex"; }); - - sm.def("overloaded2", [](py::array_t<std::complex<double>>) { return "double complex"; }); - sm.def("overloaded2", [](py::array_t<double>) { return "double"; }); - sm.def("overloaded2", [](py::array_t<std::complex<float>>) { return "float complex"; }); - sm.def("overloaded2", [](py::array_t<float>) { return "float"; }); + sm.def("overloaded", [](const py::array_t<double> &) { return "double"; }); + sm.def("overloaded", [](const py::array_t<float> &) { return "float"; }); + sm.def("overloaded", [](const py::array_t<int> &) { return "int"; }); + sm.def("overloaded", [](const py::array_t<unsigned short> &) { return "unsigned short"; }); + sm.def("overloaded", [](const py::array_t<long long> &) { return "long long"; }); + sm.def("overloaded", + [](const py::array_t<std::complex<double>> &) { return "double complex"; }); + sm.def("overloaded", [](const py::array_t<std::complex<float>> &) { return "float complex"; }); + + sm.def("overloaded2", + [](const py::array_t<std::complex<double>> &) { return "double complex"; }); + sm.def("overloaded2", [](const py::array_t<double> &) { return "double"; }); + sm.def("overloaded2", + [](const py::array_t<std::complex<float>> &) { return "float complex"; }); + sm.def("overloaded2", [](const py::array_t<float> &) { return "float"; }); + + // [workaround(intel)] ICC 20/21 breaks with py::arg().stuff, using py::arg{}.stuff works. // Only accept the exact types: - sm.def("overloaded3", [](py::array_t<int>) { return "int"; }, py::arg().noconvert()); - sm.def("overloaded3", [](py::array_t<double>) { return "double"; }, py::arg().noconvert()); + sm.def( + "overloaded3", [](const py::array_t<int> &) { return "int"; }, py::arg{}.noconvert()); + sm.def( + "overloaded3", + [](const py::array_t<double> &) { return "double"; }, + py::arg{}.noconvert()); // Make sure we don't do unsafe coercion (e.g. float to int) when not using forcecast, but // rather that float gets converted via the safe (conversion to double) overload: - sm.def("overloaded4", [](py::array_t<long long, 0>) { return "long long"; }); - sm.def("overloaded4", [](py::array_t<double, 0>) { return "double"; }); + sm.def("overloaded4", [](const py::array_t<long long, 0> &) { return "long long"; }); + sm.def("overloaded4", [](const py::array_t<double, 0> &) { return "double"; }); // But we do allow conversion to int if forcecast is enabled (but only if no overload matches // without conversion) - sm.def("overloaded5", [](py::array_t<unsigned int>) { return "unsigned int"; }); - sm.def("overloaded5", [](py::array_t<double>) { return "double"; }); + sm.def("overloaded5", [](const py::array_t<unsigned int> &) { return "unsigned int"; }); + sm.def("overloaded5", [](const py::array_t<double> &) { return "double"; }); // test_greedy_string_overload // Issue 685: ndarray shouldn't go to std::string overload - sm.def("issue685", [](std::string) { return "string"; }); - sm.def("issue685", [](py::array) { return "array"; }); - sm.def("issue685", [](py::object) { return "other"; }); + sm.def("issue685", [](const std::string &) { return "string"; }); + sm.def("issue685", [](const py::array &) { return "array"; }); + sm.def("issue685", [](const py::object &) { return "other"; }); // test_array_unchecked_fixed_dims - sm.def("proxy_add2", [](py::array_t<double> a, double v) { - auto r = a.mutable_unchecked<2>(); - for (py::ssize_t i = 0; i < r.shape(0); i++) - for (py::ssize_t j = 0; j < r.shape(1); j++) - r(i, j) += v; - }, py::arg().noconvert(), py::arg()); + sm.def( + "proxy_add2", + [](py::array_t<double> a, double v) { + auto r = a.mutable_unchecked<2>(); + for (py::ssize_t i = 0; i < r.shape(0); i++) { + for (py::ssize_t j = 0; j < r.shape(1); j++) { + r(i, j) += v; + } + } + }, + py::arg{}.noconvert(), + py::arg()); sm.def("proxy_init3", [](double start) { - py::array_t<double, py::array::c_style> a({ 3, 3, 3 }); + py::array_t<double, py::array::c_style> a({3, 3, 3}); auto r = a.mutable_unchecked<3>(); - for (py::ssize_t i = 0; i < r.shape(0); i++) - for (py::ssize_t j = 0; j < r.shape(1); j++) - for (py::ssize_t k = 0; k < r.shape(2); k++) - r(i, j, k) = start++; + for (py::ssize_t i = 0; i < r.shape(0); i++) { + for (py::ssize_t j = 0; j < r.shape(1); j++) { + for (py::ssize_t k = 0; k < r.shape(2); k++) { + r(i, j, k) = start++; + } + } + } return a; }); sm.def("proxy_init3F", [](double start) { - py::array_t<double, py::array::f_style> a({ 3, 3, 3 }); + py::array_t<double, py::array::f_style> a({3, 3, 3}); auto r = a.mutable_unchecked<3>(); - for (py::ssize_t k = 0; k < r.shape(2); k++) - for (py::ssize_t j = 0; j < r.shape(1); j++) - for (py::ssize_t i = 0; i < r.shape(0); i++) - r(i, j, k) = start++; + for (py::ssize_t k = 0; k < r.shape(2); k++) { + for (py::ssize_t j = 0; j < r.shape(1); j++) { + for (py::ssize_t i = 0; i < r.shape(0); i++) { + r(i, j, k) = start++; + } + } + } return a; }); - sm.def("proxy_squared_L2_norm", [](py::array_t<double> a) { + sm.def("proxy_squared_L2_norm", [](const py::array_t<double> &a) { auto r = a.unchecked<1>(); double sumsq = 0; - for (py::ssize_t i = 0; i < r.shape(0); i++) + for (py::ssize_t i = 0; i < r.shape(0); i++) { sumsq += r[i] * r(i); // Either notation works for a 1D array + } return sumsq; }); @@ -332,51 +377,69 @@ TEST_SUBMODULE(numpy_array, sm) { // test_array_unchecked_dyn_dims // Same as the above, but without a compile-time dimensions specification: - sm.def("proxy_add2_dyn", [](py::array_t<double> a, double v) { - auto r = a.mutable_unchecked(); - if (r.ndim() != 2) throw std::domain_error("error: ndim != 2"); - for (py::ssize_t i = 0; i < r.shape(0); i++) - for (py::ssize_t j = 0; j < r.shape(1); j++) - r(i, j) += v; - }, py::arg().noconvert(), py::arg()); + sm.def( + "proxy_add2_dyn", + [](py::array_t<double> a, double v) { + auto r = a.mutable_unchecked(); + if (r.ndim() != 2) { + throw std::domain_error("error: ndim != 2"); + } + for (py::ssize_t i = 0; i < r.shape(0); i++) { + for (py::ssize_t j = 0; j < r.shape(1); j++) { + r(i, j) += v; + } + } + }, + py::arg{}.noconvert(), + py::arg()); sm.def("proxy_init3_dyn", [](double start) { - py::array_t<double, py::array::c_style> a({ 3, 3, 3 }); + py::array_t<double, py::array::c_style> a({3, 3, 3}); auto r = a.mutable_unchecked(); - if (r.ndim() != 3) throw std::domain_error("error: ndim != 3"); - for (py::ssize_t i = 0; i < r.shape(0); i++) - for (py::ssize_t j = 0; j < r.shape(1); j++) - for (py::ssize_t k = 0; k < r.shape(2); k++) - r(i, j, k) = start++; + if (r.ndim() != 3) { + throw std::domain_error("error: ndim != 3"); + } + for (py::ssize_t i = 0; i < r.shape(0); i++) { + for (py::ssize_t j = 0; j < r.shape(1); j++) { + for (py::ssize_t k = 0; k < r.shape(2); k++) { + r(i, j, k) = start++; + } + } + } return a; }); sm.def("proxy_auxiliaries2_dyn", [](py::array_t<double> a) { return auxiliaries(a.unchecked(), a.mutable_unchecked()); }); - sm.def("array_auxiliaries2", [](py::array_t<double> a) { - return auxiliaries(a, a); - }); + sm.def("array_auxiliaries2", [](py::array_t<double> a) { return auxiliaries(a, a); }); // test_array_failures - // Issue #785: Uninformative "Unknown internal error" exception when constructing array from empty object: + // Issue #785: Uninformative "Unknown internal error" exception when constructing array from + // empty object: sm.def("array_fail_test", []() { return py::array(py::object()); }); sm.def("array_t_fail_test", []() { return py::array_t<double>(py::object()); }); // Make sure the error from numpy is being passed through: - sm.def("array_fail_test_negative_size", []() { int c = 0; return py::array(-1, &c); }); + sm.def("array_fail_test_negative_size", []() { + int c = 0; + return py::array(-1, &c); + }); // test_initializer_list // Issue (unnumbered; reported in #788): regression: initializer lists can be ambiguous - sm.def("array_initializer_list1", []() { return py::array_t<float>(1); }); // { 1 } also works, but clang warns about it - sm.def("array_initializer_list2", []() { return py::array_t<float>({ 1, 2 }); }); - sm.def("array_initializer_list3", []() { return py::array_t<float>({ 1, 2, 3 }); }); - sm.def("array_initializer_list4", []() { return py::array_t<float>({ 1, 2, 3, 4 }); }); + sm.def("array_initializer_list1", []() { return py::array_t<float>(1); }); + // { 1 } also works for the above, but clang warns about it + sm.def("array_initializer_list2", []() { return py::array_t<float>({1, 2}); }); + sm.def("array_initializer_list3", []() { return py::array_t<float>({1, 2, 3}); }); + sm.def("array_initializer_list4", []() { return py::array_t<float>({1, 2, 3, 4}); }); // test_array_resize // reshape array to 2D without changing size sm.def("array_reshape2", [](py::array_t<double> a) { - const auto dim_sz = (py::ssize_t)std::sqrt(a.size()); - if (dim_sz * dim_sz != a.size()) - throw std::domain_error("array_reshape2: input array total size is not a squared integer"); + const auto dim_sz = (py::ssize_t) std::sqrt(a.size()); + if (dim_sz * dim_sz != a.size()) { + throw std::domain_error( + "array_reshape2: input array total size is not a squared integer"); + } a.resize({dim_sz, dim_sz}); }); @@ -394,45 +457,68 @@ TEST_SUBMODULE(numpy_array, sm) { return a; }); - sm.def("index_using_ellipsis", [](py::array a) { - return a[py::make_tuple(0, py::ellipsis(), 0)]; + sm.def("array_view", + [](py::array_t<uint8_t> a, const std::string &dtype) { return a.view(dtype); }); + + sm.def("reshape_initializer_list", [](py::array_t<int> a, size_t N, size_t M, size_t O) { + return a.reshape({N, M, O}); + }); + sm.def("reshape_tuple", [](py::array_t<int> a, const std::vector<int> &new_shape) { + return a.reshape(new_shape); }); + sm.def("index_using_ellipsis", + [](const py::array &a) { return a[py::make_tuple(0, py::ellipsis(), 0)]; }); + // test_argument_conversions - sm.def("accept_double", - [](py::array_t<double, 0>) {}, - py::arg("a")); - sm.def("accept_double_forcecast", - [](py::array_t<double, py::array::forcecast>) {}, - py::arg("a")); - sm.def("accept_double_c_style", - [](py::array_t<double, py::array::c_style>) {}, - py::arg("a")); - sm.def("accept_double_c_style_forcecast", - [](py::array_t<double, py::array::forcecast | py::array::c_style>) {}, - py::arg("a")); - sm.def("accept_double_f_style", - [](py::array_t<double, py::array::f_style>) {}, - py::arg("a")); - sm.def("accept_double_f_style_forcecast", - [](py::array_t<double, py::array::forcecast | py::array::f_style>) {}, - py::arg("a")); - sm.def("accept_double_noconvert", - [](py::array_t<double, 0>) {}, - py::arg("a").noconvert()); - sm.def("accept_double_forcecast_noconvert", - [](py::array_t<double, py::array::forcecast>) {}, - py::arg("a").noconvert()); - sm.def("accept_double_c_style_noconvert", - [](py::array_t<double, py::array::c_style>) {}, - py::arg("a").noconvert()); - sm.def("accept_double_c_style_forcecast_noconvert", - [](py::array_t<double, py::array::forcecast | py::array::c_style>) {}, - py::arg("a").noconvert()); - sm.def("accept_double_f_style_noconvert", - [](py::array_t<double, py::array::f_style>) {}, - py::arg("a").noconvert()); - sm.def("accept_double_f_style_forcecast_noconvert", - [](py::array_t<double, py::array::forcecast | py::array::f_style>) {}, - py::arg("a").noconvert()); + sm.def( + "accept_double", [](const py::array_t<double, 0> &) {}, py::arg("a")); + sm.def( + "accept_double_forcecast", + [](const py::array_t<double, py::array::forcecast> &) {}, + py::arg("a")); + sm.def( + "accept_double_c_style", + [](const py::array_t<double, py::array::c_style> &) {}, + py::arg("a")); + sm.def( + "accept_double_c_style_forcecast", + [](const py::array_t<double, py::array::forcecast | py::array::c_style> &) {}, + py::arg("a")); + sm.def( + "accept_double_f_style", + [](const py::array_t<double, py::array::f_style> &) {}, + py::arg("a")); + sm.def( + "accept_double_f_style_forcecast", + [](const py::array_t<double, py::array::forcecast | py::array::f_style> &) {}, + py::arg("a")); + sm.def( + "accept_double_noconvert", [](const py::array_t<double, 0> &) {}, "a"_a.noconvert()); + sm.def( + "accept_double_forcecast_noconvert", + [](const py::array_t<double, py::array::forcecast> &) {}, + "a"_a.noconvert()); + sm.def( + "accept_double_c_style_noconvert", + [](const py::array_t<double, py::array::c_style> &) {}, + "a"_a.noconvert()); + sm.def( + "accept_double_c_style_forcecast_noconvert", + [](const py::array_t<double, py::array::forcecast | py::array::c_style> &) {}, + "a"_a.noconvert()); + sm.def( + "accept_double_f_style_noconvert", + [](const py::array_t<double, py::array::f_style> &) {}, + "a"_a.noconvert()); + sm.def( + "accept_double_f_style_forcecast_noconvert", + [](const py::array_t<double, py::array::forcecast | py::array::f_style> &) {}, + "a"_a.noconvert()); + + // Check that types returns correct npy format descriptor + sm.def("test_fmt_desc_float", [](const py::array_t<float> &) {}); + sm.def("test_fmt_desc_double", [](const py::array_t<double> &) {}); + sm.def("test_fmt_desc_const_float", [](const py::array_t<const float> &) {}); + sm.def("test_fmt_desc_const_double", [](const py::array_t<const double> &) {}); } |