pybind11_example.cpp

#include <CL/sycl.hpp>
#include <cstdint>
#include <pybind11/numpy.h>
#include <pybind11/pybind11.h>

#include "../_sycl_queue.h"
#include "../_sycl_queue_api.h"
#include "dpctl_sycl_types.h"

namespace py = pybind11;

size_t get_max_compute_units(py::object queue)
{
    PyObject *queue_ptr = queue.ptr();
    if (PyObject_TypeCheck(queue_ptr, &PySyclQueueType)) {
        DPCTLSyclQueueRef QRef =
            get_queue_ref(reinterpret_cast<PySyclQueueObject *>(queue_ptr));
        sycl::queue *q = reinterpret_cast<sycl::queue *>(QRef);

        return q->get_device()
            .get_info<sycl::info::device::max_compute_units>();
    }
    else {
        throw std::runtime_error("expected dpctl.SyclQueue as argument");
    }
}

py::array_t<int64_t>
offloaded_array_mod(py::object queue,
                    py::array_t<int64_t, py::array::c_style> array,
                    int64_t mod)
{
    sycl::queue *q_ptr;

    PyObject *queue_pycapi = queue.ptr();
    if (PyObject_TypeCheck(queue_pycapi, &PySyclQueueType)) {
        DPCTLSyclQueueRef QRef =
            get_queue_ref(reinterpret_cast<PySyclQueueObject *>(queue_pycapi));
        q_ptr = reinterpret_cast<sycl::queue *>(QRef);
    }
    else {
        throw std::runtime_error("expected dpctl.SyclQueue as argument");
    }

    py::buffer_info arg_pybuf = array.request();
    if (arg_pybuf.ndim != 1) {
        throw std::runtime_error("Expecting a vector");
    }
    if (mod <= 0) {
        throw std::runtime_error("Modulus must be non-negative");
    }

    size_t n = arg_pybuf.size;

    auto res = py::array_t<int64_t>(n);
    py::buffer_info res_pybuf = res.request();

    int64_t *a = static_cast<int64_t *>(arg_pybuf.ptr);
    int64_t *r = static_cast<int64_t *>(res_pybuf.ptr);

    {
        const sycl::property_list props = {
            sycl::property::buffer::use_host_ptr()};
        sycl::buffer<int64_t, 1> a_buf(a, sycl::range<1>(n), props);
        sycl::buffer<int64_t, 1> r_buf(r, sycl::range<1>(n), props);

        q_ptr
            ->submit([&](sycl::handler &cgh) {
                sycl::accessor a_acc(a_buf, cgh, sycl::read_only);
                sycl::accessor r_acc(r_buf, cgh, sycl::write_only,
                                     sycl::noinit);

                cgh.parallel_for(sycl::range<1>(n), [=](sycl::id<1> idx) {
                    r_acc[idx] = a_acc[idx] % mod;
                });
            })
            .wait_and_throw();
    }

    return res;
}

PYBIND11_MODULE(pybind11_example, m)
{
    // Import the dpctl._sycl_queue extension
    import_dpctl___sycl_queue();
    m.def("get_max_compute_units", &get_max_compute_units,
          "Computes max_compute_units property of the device underlying given "
          "dpctl.SyclQueue");
    m.def("offloaded_array_mod", &offloaded_array_mod,
          "Compute offloaded modular reduction of integer-valued NumPy array");
}