Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 1 | /* |
| 2 | example/example20.cpp -- Usage of structured numpy dtypes |
| 3 | |
| 4 | Copyright (c) 2016 Ivan Smirnov |
| 5 | |
| 6 | All rights reserved. Use of this source code is governed by a |
| 7 | BSD-style license that can be found in the LICENSE file. |
| 8 | */ |
| 9 | |
| 10 | #include "example.h" |
| 11 | |
| 12 | #include <pybind11/numpy.h> |
| 13 | #include <cstdint> |
| 14 | #include <iostream> |
| 15 | |
| 16 | namespace py = pybind11; |
| 17 | |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 18 | struct SimpleStruct { |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 19 | bool x; |
| 20 | uint32_t y; |
| 21 | float z; |
| 22 | }; |
| 23 | |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 24 | std::ostream& operator<<(std::ostream& os, const SimpleStruct& v) { |
| 25 | return os << "s:" << v.x << "," << v.y << "," << v.z; |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 26 | } |
| 27 | |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 28 | struct PackedStruct { |
| 29 | bool x; |
| 30 | uint32_t y; |
| 31 | float z; |
| 32 | } __attribute__((packed)); |
| 33 | |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 34 | std::ostream& operator<<(std::ostream& os, const PackedStruct& v) { |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 35 | return os << "p:" << v.x << "," << v.y << "," << v.z; |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 36 | } |
| 37 | |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 38 | struct NestedStruct { |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 39 | SimpleStruct a; |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 40 | PackedStruct b; |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 41 | } __attribute__((packed)); |
| 42 | |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 43 | std::ostream& operator<<(std::ostream& os, const NestedStruct& v) { |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 44 | return os << "n:a=" << v.a << ";b=" << v.b; |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 45 | } |
| 46 | |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 47 | struct PartialStruct { |
| 48 | bool x; |
| 49 | uint32_t y; |
| 50 | float z; |
| 51 | long dummy2; |
| 52 | }; |
| 53 | |
| 54 | struct PartialNestedStruct { |
| 55 | long dummy1; |
| 56 | PartialStruct a; |
| 57 | long dummy2; |
| 58 | }; |
| 59 | |
Ivan Smirnov | d0bafd9 | 2016-06-26 16:35:28 +0100 | [diff] [blame] | 60 | struct UnboundStruct { }; |
| 61 | |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 62 | template <typename T> |
| 63 | py::array mkarray_via_buffer(size_t n) { |
| 64 | return py::array(py::buffer_info(nullptr, sizeof(T), |
Ivan Smirnov | 5e71e17 | 2016-06-26 12:42:34 +0100 | [diff] [blame] | 65 | py::format_descriptor<T>::format(), |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 66 | 1, { n }, { sizeof(T) })); |
| 67 | } |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 68 | |
| 69 | template <typename S> |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 70 | py::array_t<S, 0> create_recarray(size_t n) { |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 71 | auto arr = mkarray_via_buffer<S>(n); |
| 72 | auto ptr = static_cast<S*>(arr.request().ptr); |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 73 | for (size_t i = 0; i < n; i++) { |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 74 | ptr[i].x = i % 2; ptr[i].y = (uint32_t) i; ptr[i].z = (float) i * 1.5f; |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 75 | } |
| 76 | return arr; |
| 77 | } |
| 78 | |
Ivan Smirnov | d0bafd9 | 2016-06-26 16:35:28 +0100 | [diff] [blame] | 79 | std::string get_format_unbound() { |
| 80 | return py::format_descriptor<UnboundStruct>::format(); |
| 81 | } |
| 82 | |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 83 | py::array_t<NestedStruct, 0> create_nested(size_t n) { |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 84 | auto arr = mkarray_via_buffer<NestedStruct>(n); |
| 85 | auto ptr = static_cast<NestedStruct*>(arr.request().ptr); |
| 86 | for (size_t i = 0; i < n; i++) { |
| 87 | ptr[i].a.x = i % 2; ptr[i].a.y = (uint32_t) i; ptr[i].a.z = (float) i * 1.5f; |
| 88 | ptr[i].b.x = (i + 1) % 2; ptr[i].b.y = (uint32_t) (i + 1); ptr[i].b.z = (float) (i + 1) * 1.5f; |
| 89 | } |
| 90 | return arr; |
Ivan Smirnov | bdc9902 | 2016-06-19 16:54:07 +0100 | [diff] [blame] | 91 | } |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 92 | |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 93 | py::array_t<PartialNestedStruct, 0> create_partial_nested(size_t n) { |
| 94 | auto arr = mkarray_via_buffer<PartialNestedStruct>(n); |
| 95 | auto ptr = static_cast<PartialNestedStruct*>(arr.request().ptr); |
| 96 | for (size_t i = 0; i < n; i++) { |
| 97 | ptr[i].a.x = i % 2; ptr[i].a.y = (uint32_t) i; ptr[i].a.z = (float) i * 1.5f; |
| 98 | } |
| 99 | return arr; |
| 100 | } |
| 101 | |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 102 | template <typename S> |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 103 | void print_recarray(py::array_t<S, 0> arr) { |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 104 | auto buf = arr.request(); |
| 105 | auto ptr = static_cast<S*>(buf.ptr); |
| 106 | for (size_t i = 0; i < buf.size; i++) |
| 107 | std::cout << ptr[i] << std::endl; |
| 108 | } |
Ivan Smirnov | bdc9902 | 2016-06-19 16:54:07 +0100 | [diff] [blame] | 109 | |
| 110 | void print_format_descriptors() { |
Ivan Smirnov | 5e71e17 | 2016-06-26 12:42:34 +0100 | [diff] [blame] | 111 | std::cout << py::format_descriptor<SimpleStruct>::format() << std::endl; |
| 112 | std::cout << py::format_descriptor<PackedStruct>::format() << std::endl; |
| 113 | std::cout << py::format_descriptor<NestedStruct>::format() << std::endl; |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 114 | std::cout << py::format_descriptor<PartialStruct>::format() << std::endl; |
| 115 | std::cout << py::format_descriptor<PartialNestedStruct>::format() << std::endl; |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 116 | } |
| 117 | |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 118 | void print_dtypes() { |
| 119 | auto to_str = [](py::object obj) { |
| 120 | return (std::string) (py::str) ((py::object) obj.attr("__str__"))(); |
| 121 | }; |
| 122 | std::cout << to_str(py::dtype_of<SimpleStruct>()) << std::endl; |
| 123 | std::cout << to_str(py::dtype_of<PackedStruct>()) << std::endl; |
| 124 | std::cout << to_str(py::dtype_of<NestedStruct>()) << std::endl; |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 125 | std::cout << to_str(py::dtype_of<PartialStruct>()) << std::endl; |
| 126 | std::cout << to_str(py::dtype_of<PartialNestedStruct>()) << std::endl; |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 127 | } |
| 128 | |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 129 | void init_ex20(py::module &m) { |
Ivan Smirnov | 5412a05 | 2016-07-02 16:18:42 +0100 | [diff] [blame] | 130 | PYBIND11_NUMPY_DTYPE(SimpleStruct, x, y, z); |
| 131 | PYBIND11_NUMPY_DTYPE(PackedStruct, x, y, z); |
| 132 | PYBIND11_NUMPY_DTYPE(NestedStruct, a, b); |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 133 | PYBIND11_NUMPY_DTYPE(PartialStruct, x, y, z); |
| 134 | PYBIND11_NUMPY_DTYPE(PartialNestedStruct, a); |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 135 | |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 136 | m.def("create_rec_simple", &create_recarray<SimpleStruct>); |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 137 | m.def("create_rec_packed", &create_recarray<PackedStruct>); |
Ivan Smirnov | 7f913ae | 2016-06-19 16:41:15 +0100 | [diff] [blame] | 138 | m.def("create_rec_nested", &create_nested); |
Ivan Smirnov | 8fa09cb | 2016-07-06 00:28:12 +0100 | [diff] [blame] | 139 | m.def("create_rec_partial", &create_recarray<PartialStruct>); |
| 140 | m.def("create_rec_partial_nested", &create_partial_nested); |
Ivan Smirnov | bdc9902 | 2016-06-19 16:54:07 +0100 | [diff] [blame] | 141 | m.def("print_format_descriptors", &print_format_descriptors); |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 142 | m.def("print_rec_simple", &print_recarray<SimpleStruct>); |
Ivan Smirnov | 669e142 | 2016-06-21 21:05:29 +0100 | [diff] [blame] | 143 | m.def("print_rec_packed", &print_recarray<PackedStruct>); |
| 144 | m.def("print_rec_nested", &print_recarray<NestedStruct>); |
Ivan Smirnov | 4f16421 | 2016-06-22 01:07:20 +0100 | [diff] [blame] | 145 | m.def("print_dtypes", &print_dtypes); |
Ivan Smirnov | d0bafd9 | 2016-06-26 16:35:28 +0100 | [diff] [blame] | 146 | m.def("get_format_unbound", &get_format_unbound); |
Ivan Smirnov | bb4015d | 2016-06-19 15:50:31 +0100 | [diff] [blame] | 147 | } |