xtensor-python/test/main.cpp at c72ec9b67edd0d0e2fd183079f670fd6264e7182 · wolfv/xtensor-python · GitHub

Repository navigation

History

#define NPY_NO_DEPRECATED_API NPY_1_7_API_VERSION

#include "numpy/arrayobject.h"

#include "xtensor/xmath.hpp"

#include "xtensor/xarray.hpp"

#include "xtensor-python/pyarray.hpp"

#include "xtensor-python/pyvectorize.hpp"

#include <numeric>

namespace py = pybind11;

using complex_t = std::complex<double>;

// Examples

double example1(xt::pyarray<double>& m)

{

return m(0);

}

xt::pyarray<double> example2(xt::pyarray<double>& m)

{

return m + 2;

}

// Readme Examples

double readme_example1(xt::pyarray<double>& m)

{

auto sines = xt::sin(m);

return std::accumulate(sines.begin(), sines.end(), 0.0);

}

double readme_example2(double i, double j)

{

return std::sin(i) - std::cos(j);

}

// Vectorize Examples

int add(int i, int j)

{

return i + j;

}

PYBIND11_PLUGIN(xtensor_python_test)

{

if(_import_array() < 0)

{

PyErr_SetString(PyExc_ImportError, "numpy.core.multiarray failed to import");

return nullptr;

}

py::module m("xtensor_python_test", "Test module for xtensor python bindings");

m.def("example1", example1, "");

m.def("example2", example2, "");

m.def("readme_example1", readme_example1, "");

m.def("readme_example2", xt::pyvectorize(readme_example2), "");

m.def("vectorize_example1", xt::pyvectorize(add), "");

m.def("rect_to_polar", [](xt::pyarray<complex_t> const& a) {

return py::make_tuple(xt::pyvectorize([](complex_t x) { return std::abs(x); })(a),

xt::pyvectorize([](complex_t x) { return std::arg(x); })(a));

});

return m.ptr();

}

xtensor-python/test/main.cpp at c72ec9b67edd0d0e2fd183079f670fd6264e7182 · wolfv/xtensor-python · GitHub

Repository navigation

Expand file tree

Breadcrumbs

main.cpp

Latest commit

History

Breadcrumbs

main.cpp

File metadata and controls