scikit build

CMakeLists.txt

@@ -0,0 +1,11 @@
+cmake_minimum_required(VERSION 3.15...3.26)
+project(${SKBUILD_PROJECT_NAME} LANGUAGES CXX)
+
+set(PYBIND11_NEWPYTHON ON)
+find_package(pybind11 CONFIG REQUIRED)
+set(TARGET pixel_math_ext)
+
+pybind11_add_module(${TARGET} src/cpp/pixel_math.cpp )
+
+target_include_directories(${TARGET} PUBLIC src/cpp)
+install(TARGETS ${TARGET} LIBRARY DESTINATION ${SKBUILD_PROJECT_NAME})

pyproject.toml

@@ -52,10 +52,14 @@ dev = [
 
 [build-system]
 requires = [
-    "setuptools>=62", # Used to build and package the Python project
     "setuptools_scm>=6.2", # Gets release version from git. Makes it available programmatically
+    "scikit-build-core", "pybind11"
 ]
-build-backend = "setuptools.build_meta"
+build-backend = "scikit_build_core.build"
+
+[tool.scikit-build]
+metadata.version.provider = "scikit_build_core.metadata.setuptools_scm"
+sdist.include = ["src/hipscat/_version.py"]
 
 [tool.setuptools_scm]
 write_to = "src/hipscat/_version.py"

src/cpp/pixel_math.cpp

@@ -0,0 +1,16 @@
+#include <pybind11/numpy.h>
+#include <pybind11/pybind11.h>
+#include <pybind11/stl.h>
+
+#include <vector>
+#include "pixel_math.h"
+
+std::vector<std::vector<unsigned long long int> > generate_alignment(
+    const std::vector<unsigned long long int> &histogram,
+    unsigned short highest_order,
+    unsigned short lowest_order,
+    unsigned long long int threshold
+){
+    std::vector<std::vector<unsigned long long int> > result (order2npix(highest_order));
+    return result;
+}

src/cpp/pixel_math.h

@@ -0,0 +1,30 @@
+#ifndef PIXEL_MATH_H_
+#define  PIXEL_MATH_H_
+
+#include <pybind11/pybind11.h>
+
+#include <vector>
+
+namespace py = pybind11;
+
+
+float square(float x) { return x * x; }
+
+unsigned long order2npix(unsigned short pix) {
+    return 12 * 1 << (2 * pix);
+}
+
+std::vector<std::vector<unsigned long long int> > generate_alignment(
+    const std::vector<unsigned long long int> &histogram,
+    unsigned short highest_order,
+    unsigned short lowest_order,
+    unsigned long long int threshold
+);
+
+PYBIND11_MODULE(pixel_math_ext, m) {
+    m.def("square", &square);
+    m.def("order2npix", &order2npix);
+    m.def("generate_alignment", &generate_alignment);
+}
+
+#endif

tests/hipscat/pixel_math/test_the_cpp.py

@@ -0,0 +1,9 @@
+from hipscat.pixel_math_ext import order2npix, generate_alignment
+import numpy as np
+
+def test_order2npix():
+    assert order2npix(0) == 12
+    assert order2npix(4) == 3072
+
+def test_generate_alignment():
+    assert len(generate_alignment(np.arange(12), 0, 0, 4)) == 12