All project tests now really passing

Author: tskisner

include/Projection.h

@@ -40,7 +40,7 @@ class SignalSpace {
     py::object ret_val;
 
 private:
-    bool _Validate(py::object input, std::string var_name,
+    constexpr bool _Validate(py::object input, std::string var_name,
                    int dtype);
 };
 

include/numpy_assist.h

@@ -133,7 +133,7 @@ class BufferWrapper {
     // Constructor with no shape or type checking.
     BufferWrapper(std::string name, const py::object &src, bool optional)
         : BufferWrapper() {
-        if (optional && (src.ptr() == Py_None))
+        if (optional && (src.is_none()))
             return;
         //std::cerr << "BufferWrapper ctor " << name << " src = " << src.ptr() << ")" << std::endl;
         if (PyObject_GetBuffer(src.ptr(), view.get(),
@@ -156,27 +156,27 @@ class BufferWrapper {
                   std::vector<int> shape)
         : BufferWrapper(name, src, optional) {
 
-        //std::cerr << "BufferWrapper start 4-arg ctor" << std::endl;
+        std::cerr << "BufferWrapper start 4-arg ctor" << std::endl;
 
         // "optional" items will cause the parent constructor to
         // succeed, but will leave buffer pointer unset.
         if (view->buf == NULL) {
-            //std::cerr << "BufferWrapper: view->buf == NULL, return" << std::endl;
+            std::cerr << "BufferWrapper: view->buf == NULL, return" << std::endl;
             return;
         }
 
         if (!check_buffer_type<T>(*(view.get()))) {
-            //std::cerr << "BufferWrapper: view invalid type = " << type_name<T>() << std::endl;
+            std::cerr << "BufferWrapper: view invalid type = " << type_name<T>() << std::endl;
             throw dtype_exception(name, type_name<T>());
         }
 
         std::vector<int> vshape;
-        //std::cerr << "BufferWrapper: vshape = ";
+        std::cerr << "BufferWrapper: vshape = ";
         for (int i=0; i<view->ndim; i++) {
-            //std::cerr << view->shape[i] << ",";
+            std::cerr << view->shape[i] << ",";
             vshape.push_back(view->shape[i]);
         }
-        //std::cerr << std::endl;
+        std::cerr << std::endl;
 
         // Note special value -1 is as in numpy -- matches a single
         // axis.  Special value -2 is treated as an ellipsis -- can be
@@ -205,6 +205,7 @@ class BufferWrapper {
             std::ostringstream s;
             s << "Expected " << shape_string(shape) << " but got " <<
                 shape_string(vshape) << ".";
+            std::cerr << s.str() << std::endl;
             throw shape_exception(name, s.str());
         }
     }
@@ -229,3 +230,7 @@ class BufferWrapper {
 private:
     std::shared_ptr<Py_buffer> view;
 };
+
+
+// Convert an n-dimensional array into a list of array slices.
+py::list list_of_arrays(py::object input);

pyproject.toml

@@ -57,7 +57,8 @@ so-hk-tool = "so3g.hk.cli:main"
 
 [tool.scikit-build]
 # A list of args to pass to CMake when configuring the project.
-cmake.args = []
+#cmake.args = []
+cmake.args = ["-DCMAKE_C_FLAGS='-O0 -g'", "-DCMAKE_CXX_FLAGS='-O0 -g'"]
 # A table of defines to pass to CMake when configuring the project. Additive.
 cmake.define = {}
 # The build type to use when building the project.
@@ -110,7 +111,7 @@ build.requires = []
 # The components to install.
 install.components = []
 # Whether to strip the binaries.
-install.strip = true
+install.strip = false
 # Add the python build environment site_packages folder to the CMake prefix paths.
 search.site-packages = true
 # Strictly check all config options.

python/so3g/proj/wcs.py

@@ -606,6 +606,8 @@ def for_geom(cls, shape, wcs, interpol=None):
             interpol = INTERPOLS[0]
         self.interpol = interpol
 
+        print(f"for_geom: naxis = {self.naxis}, cdelt = {self.cdelt}, crpix = {self.crpix}", flush=True)
+
         return self
 
     @classmethod